Modeling The Correlation Structure Of Rna Sequencing Data Using A Multivariate Poisson-Lognormal Model