Accurate gene model annotation of reference genomes is critical for making them useful. The modENCODE project has improved the D. melanogaster genome annotation by using deep and diverse high-throughput data. Since transcriptional activity that has been evolutionarily conserved is likely to have an advantageous function, we have performed large-scale interspecific comparisons to increase confidence in predicted annotations. To support comparative genomics, we filled in divergence gaps in the Drosophila phylogeny by generating draft genomes for eight new species. For comparative transcriptome analysis, we generated mRNA expression profiles on 81 samples from multiple tissues and developmental stages of 15 Drosophila species, and we performed cap analysis of gene expression in D. melanogaster and D. pseudoobscura. We also describe conservation of four distinct core promoter structures composed of combinations of elements at three positions. Overall, each type of genomic feature shows a characteristic divergence rate relative to neutral models, highlighting the value of multispecies alignment in annotating a target genome that should prove useful in the annotation of other high priority genomes, especially human and other mammalian genomes that are rich in noncoding sequences. We report that the vast majority of elements in the annotation are evolutionarily conserved, indicating that the annotation will be an important springboard for functional genetic testing by the Drosophila community.

A generally accepted approach to the analysis of RNA-Seq read count data does not yet exist. We sequenced the mRNA of 726 individuals from the Drosophila Genetic Reference Panel in order to quantify differences in gene expression among single flies. One of our experimental goals was to identify the optimal analysis approach for the detection of differential gene expression among the factors we varied in the experiment: genotype, environment, sex, and their interactions. Here we evaluate three different filtering strategies, eight normalization methods, and two statistical approaches using our data set. We assessed differential gene expression among factors and performed a statistical power analysis using the eight biological replicates per genotype, environment, and sex in our data set. We found that the most critical considerations for the analysis of RNA-Seq read count data were the normalization method, underlying data distribution assumption, and numbers of biological replicates, an observation consistent with previous RNA-Seq and microarray analysis comparisons. Some common normalization methods, such as Total Count, Quantile, and RPKM normalization, did not align the data across samples. Furthermore, analyses using the Median, Quantile, and Trimmed Mean of M-values normalization methods were sensitive to the removal of low-expressed genes from the data set. Although it is robust in many types of analysis, the normal data distribution assumption produced results vastly different than the negative binomial distribution. In addition, at least three biological replicates per condition were required in order to have sufficient statistical power to detect expression differences among the three-way interaction of genotype, environment, and sex. The best analysis approach to our data was to normalize the read counts using the DESeq method and apply a generalized linear model assuming a negative binomial distribution using either edgeR or DESeq software. Genes having very low read counts were removed after normalizing the data and fitting it to the negative binomial distribution. We describe the results of this evaluation and include recommended analysis strategies for RNA-Seq read count data.

Men demonstrate higher incidence and mortality rates of colorectal cancer (CRC) than women. This study aims to explain the potential causes of such sexual dimorphism in CRC from the perspective of sex‐biased gut microbiota and metabolites. The results show that sexual dimorphism in colorectal tumorigenesis is observed in both Apc^Min/+ mice and azoxymethane (AOM)/dextran sulfate sodium (DSS)‐treated mice with male mice have significantly larger and more tumors, accompanied by more impaired gut barrier function. Moreover, pseudo‐germ mice receiving fecal samples from male mice or patients show more severe intestinal barrier damage and higher level of inflammation. A significant change in gut microbiota composition is found with increased pathogenic bacteria Akkermansia muciniphila and deplets probiotic Parabacteroides goldsteinii in both male mice and pseudo‐germ mice receiving fecal sample from male mice. Sex‐biased gut metabolites in pseudo‐germ mice receiving fecal sample from CRC patients or CRC mice contribute to sex dimorphism in CRC tumorigenesis through glycerophospholipids metabolism pathway. Sexual dimorphism in tumorigenesis of CRC mouse models. In conclusion, the sex‐biased gut microbiome and metabolites contribute to sexual dimorphism in CRC. Modulating sex‐biased gut microbiota and metabolites could be a potential sex‐targeting therapeutic strategy of CRC.

Background: Pseudogenes are excellent markers for genome evolution, which are emerging as crucial regulators of development and disease, especially cancer. However, systematic functional characterization and evolution of pseudogenes remain largely unexplored. Results: To systematically characterize pseudogenes, we date the origin of human and mouse pseudogenes across vertebrates and observe a burst of pseudogene gain in these two lineages. Based on a hybrid sequencing dataset combining full-length PacBio sequencing, sample-matched Illumina sequencing, and public time-course transcriptome data, we observe that abundant mammalian pseudogenes could be transcribed, which contribute to the establishment of organ identity. Our analyses reveal that developmentally dynamic pseudogenes are evolutionarily conserved and show an increasing weight during development. Besides, they are involved in complex transcriptional and post-transcriptional modulation, exhibiting the signatures of functional enrichment. Coding potential evaluation suggests that 19% of human pseudogenes could be translated, thus serving as a new way for protein innovation. Moreover, pseudogenes carry disease-associated SNPs and conduce to cancer transcriptome perturbation. Conclusions: Our discovery reveals an unexpectedly high abundance of mammalian pseudogenes that can be transcribed and translated, and these pseudogenes represent a novel regulatory layer. Our study also prioritizes developmentally dynamic pseudogenes with signatures of functional enrichment and provides a hybrid sequencing dataset for further unraveling their biological mechanisms in organ development and carcinogenesis in the future.

Differences in phenotype among genetically identical individuals exposed to the same environmental condition are often noted in genetic studies. Despite this commonplace observation, little is known about the causes of this variability, which has been termed microenvironmental plasticity. One possibility is that stochastic or technical sources of variance produce these differences. A second possibility is that this variation has a genetic component. We have explored gene expression robustness in the transcriptomes of 730 individual Drosophila melanogaster of 16 fixed genotypes, 9 of which are infected with Wolbachia. Three replicates of flies were grown controlling for food, day/night cycles, humidity, temperature, sex, mating status, social exposure, and circadian timing of RNA extraction. Despite the use of inbred genotypes and carefully controlled experimental conditions, thousands of genes were differentially expressed, revealing a unique and dynamic transcriptional signature for each individual fly. We found that 23% of the transcriptome was differentially expressed among individuals, and that the variability in gene expression among individuals is influenced by genotype. This transcriptional variation originated from specific gene pathways, suggesting a plastic response to the micro-environment; but there was also evidence of gene expression differences due to stochastic fluctuations. These observations reveal previously unappreciated genetic sources of variability in gene expression among individuals, which has implications for complex trait genetics and precision medicine.