Retro-miRs: novel and functional miRNAs originating from mRNA retrotransposition.

BACKGROUND: Reverse-transcribed gene copies (retrocopies) have emerged as major sources of evolutionary novelty. MicroRNAs (miRNAs) are small and highly conserved RNA molecules that serve as key post-transcriptional regulators of gene expression. The origin and subsequent evolution of miRNAs have been addressed but not fully elucidated. RESULTS: In this study, we performed a comprehensive investigation of miRNA origination through retroduplicated mRNA sequences (retro-miRs). We identified 17 retro-miRs that emerged from the mRNA retrocopies. Four of these retro-miRs had de novo origins within retrocopied sequences, while 13 retro-miRNAs were located within exon regions and duplicated along with their host mRNAs. We found that retro-miRs were primate-specific, including five retro-miRs conserved among all primates and two human-specific retro-miRs. All retro-miRs were expressed, with predicted and experimentally validated target genes except miR-10527. Notably, the target genes of retro-miRs are involved in key biological processes such as metabolic processes, cell signaling, and regulation of neurotransmitters in the central nervous system. Additionally, we found that these retro-miRs play a potential oncogenic role in cancer by targeting key cancer genes and are overexpressed in several cancer types, including liver hepatocellular carcinoma and stomach adenocarcinoma. CONCLUSIONS: Our findings demonstrated that mRNA retrotransposition is a key mechanism for the generation of novel miRNAs (retro-miRs) in primates. These retro-miRs are expressed, conserved, have target genes with important cellular functions, and play important roles in cancer.

Genomic analyses of new genes and their phenotypic effects reveal rapid evolution of essential functions in Drosophila development.

It is a conventionally held dogma that the genetic basis underlying development is conserved in a long evolutionary time scale. Ample experiments based on mutational, biochemical, functional, and complementary knockdown/knockout approaches have revealed the unexpectedly important role of recently evolved new genes in the development of Drosophila. The recent progress in the genome-wide experimental testing of gene effects and improvements in the computational identification of new genes (< 40 million years ago, Mya) open the door to investigate the evolution of gene essentiality with a phylogenetically high resolution. These advancements also raised interesting issues in techniques and concepts related to phenotypic effect analyses of genes, particularly of those that recently originated. Here we reported our analyses of these issues, including reproducibility and efficiency of knockdown experiment and difference between RNAi libraries in the knockdown efficiency and testing of phenotypic effects. We further analyzed a large data from knockdowns of 11,354 genes (~75% of the Drosophila melanogaster total genes), including 702 new genes (~66% of the species total new genes that aged < 40 Mya), revealing a similarly high proportion (~32.2%) of essential genes that originated in various Sophophora subgenus lineages and distant ancestors beyond the Drosophila genus. The transcriptional compensation effect from CRISPR knockout were detected for highly similar duplicate copies. Knockout of a few young genes detected analogous essentiality in various functions in development. Taken together, our experimental and computational analyses provide valuable data for detection of phenotypic effects of genes in general and further strong evidence for the concept that new genes in Drosophila quickly evolved essential functions in viability during development.

Clinal and seasonal changes are correlated in Drosophila melanogaster natural populations.

Spatial and seasonal variations in the environment are ubiquitous. Environmental heterogeneity can affect natural populations and lead to covariation between environment and allele frequencies. Drosophila melanogaster is known to harbor polymorphisms that change both with latitude and seasons. Identifying the role of selection in driving these changes is not trivial, because nonadaptive processes can cause similar patterns. Given the environment changes in similar ways across seasons and along the latitudinal gradient, one promising approach may be to look for parallelism between clinal and seasonal changes. Here, we test whether there is a genome-wide correlation between clinal and seasonal changes, and whether the pattern is consistent with selection. Allele frequency estimates were obtained from pooled samples from seven different locations along the east coast of the United States, and across seasons within Pennsylvania. We show that there is a genome-wide correlation between clinal and seasonal variations, which cannot be explained by linked selection alone. This pattern is stronger in genomic regions with higher functional content, consistent with natural selection. We derive a way to biologically interpret these correlations and show that around 3.7% of the common, autosomal variants could be under parallel seasonal and spatial selection. Our results highlight the contribution of natural selection in driving fluctuations in allele frequencies in natural fly populations and point to a shared genomic basis to climate adaptation that happens over space and time in D. melanogaster.

Dynamic sex chromosome expression in Drosophila male germ cells.

Given their copy number differences and unique modes of inheritance, the evolved gene content and expression of sex chromosomes is unusual. In many organisms the X and Y chromosomes are inactivated in spermatocytes, possibly as a defense mechanism against insertions into unpaired chromatin. In addition to current sex chromosomes, Drosophila has a small gene-poor X-chromosome relic (4th) that re-acquired autosomal status. Here we use single cell RNA-Seq on fly larvae to demonstrate that the single X and pair of 4th chromosomes are specifically inactivated in primary spermatocytes, based on measuring all genes or a set of broadly expressed genes in testis we identified. In contrast, genes on the single Y chromosome become maximally active in primary spermatocytes. Reduced X transcript levels are due to failed activation of RNA-Polymerase-II by phosphorylation of Serine 2 and 5.

Haploid selection drives new gene male germline expression.

New genes are a major source of novelties, and a disproportionate amount of them are known to show testis expression in later phases of male gametogenesis in different groups such as mammals and plants. Here, we propose that this enhanced expression is a consequence of haploid selection during the latter stages of male gametogenesis. Because emerging adaptive mutations will be fixed faster if their phenotypes are expressed by haploid rather than diploid genotypes, new genes with advantageous functions arising during this unique stage of development have a better chance to become fixed. To test this hypothesis, expression levels of genes of differing evolutionary age were examined at various stages of Drosophila spermatogenesis. We found, consistent with a model based on haploid selection, that new Drosophila genes are both expressed in later haploid phases of spermatogenesis and harbor a significant enrichment of adaptive mutations. Additionally, the observed overexpression of new genes in the latter phases of spermatogenesis was limited to the autosomes. Because all male cells exhibit hemizygous expression for X-linked genes (and therefore effectively haploid), there is no expectation that selection acting on late spermatogenesis will have a different effect on X-linked genes in comparison to initial diploid phases. Together, our proposed hypothesis and the analyzed data suggest that natural selection in haploid cells elucidates several aspects of the origin of new genes by explaining the general prevalence of their testis expression, and a parsimonious solution for new alleles to avoid being lost by genetic drift or pseudogenization.

Early X chromosome inactivation during human preimplantation development revealed by single-cell RNA-sequencing.

In female mammals, one X chromosome is transcriptionally inactivated (XCI), leading to dosage compensation between sexes, fundamental for embryo viability. A previous study using single-cell RNA-sequencing (scRNA-seq) data proposed that female human preimplantation embryos achieve dosage compensation by downregulating both Xs, a phenomenon named dampening of X expression. Using a novel pipeline on those data, we identified a decrease in the proportion of biallelically expressed X-linked genes during development, consistent with XCI. Moreover, we show that while the expression sum of biallelically expressed X-linked genes decreases with embryonic development, their median expression remains constant, rejecting the hypothesis of X dampening. In addition, analyses of a different dataset of scRNA-seq suggest the appearance of X-linked monoallelic expression by the late blastocyst stage in females, another hallmark of initiation of XCI. Finally, we addressed the issue of dosage compensation between the single active X and autosomes in males and females for the first time during human preimplantation development, showing emergence of X to autosome dosage compensation by the upregulation of the active X chromosome in both male and female embryonic stem cells. Our results show compelling evidence of an early process of X chromosome inactivation during human preimplantation development.

Unveiling the Impact of the Genomic Architecture on the Evolution of Vertebrate microRNAs.

Eukaryotic genomes frequently exhibit interdependency between transcriptional units, as evidenced by regions of high gene density. It is well recognized that vertebrate microRNAs (miRNAs) are usually embedded in those regions. Recent work has shown that the genomic context is of utmost importance to determine miRNA expression in time and space, thus affecting their evolutionary fates over long and short terms. Consequently, understanding the inter- and intraspecific changes on miRNA genomic architecture may bring novel insights on the basic cellular processes regulated by miRNAs, as well as phenotypic evolution and disease-related mechanisms.

Host gene constraints and genomic context impact the expression and evolution of human microRNAs.

Increasing evidence has shown that recent miRNAs tend to emerge within coding genes. Here we conjecture that human miRNA evolution is tightly influenced by the genomic context, especially by host genes. Our findings show a preferential emergence of intragenic miRNAs within old genes. We found that miRNAs within old host genes are significantly more broadly expressed than those within young ones. Young miRNAs within old genes are more broadly expressed than their intergenic counterparts, suggesting that young miRNAs have an initial advantage by residing in old genes, and benefit from their hosts' expression control and from the exposure to diverse cellular contexts and target genes. Our results demonstrate that host genes may provide stronger expression constraints to intragenic miRNAs in the long run. We also report associated functional implications, highlighting the genomic context and host genes as driving factors for the expression and evolution of human miRNAs.

Meiotic sex chromosome inactivation in Drosophila.

In several different taxa, there is indubitable evidence of transcriptional silencing of the X and Y chromosomes in male meiotic cells of spermatogenesis. However, the so called meiotic sex chromosome inactivation (MSCI) has been recently a hot bed for debate in Drosophila melanogaster. This review covers cytological and genetic observations, data from transgenic constructs with testis-specific promoters, global expression profiles obtained from mutant, wild-type, larvae and adult testes as well as from cells of different stages of spermatogenesis. There is no dispute on that D. melanogaster spermatogenesis presents a down-regulation of X chromosome that does not result from the lack of dosage compensation. However, the issue is currently focused on the level of reduction of X-linked expression, the precise time it occurs and how many genes are affected. The deep examination of data and experiments in this review exposes the limitations intrinsic to the methods of studying MSCI in D. melanogaster. The current methods do not allow us to affirm anything else than the X chromosome down-regulation in meiosis (MSCI). Therefore, conclusion about level, degree or precise timing is inadequate until new approaches are implemented to know the details of MSCI or other processes involved for D. melanogaster model.

A novel dataset for identifying sex-biased genes in Drosophila.

Phenotypic differences between males and females of sexually dimorphic species are caused in large part by differences in gene expression between the sexes, most of which occurs in the gonads. To accurately identify genes differentially expressed between males and females in Drosophila, we sequenced the testis and ovary transcriptomes of D. yakuba, D. pseudoobscura, and D. ananassae and used them to identify sex-biased genes in the latter two species. We highlight the increased sensitivity and improved power of sex-biased gene detection methods when using our testis/ovary data versus male and female whole body transcriptome data. We thus provide a resource specifically designed to accurately identify and characterize sex-biased genes across Drosophila. This dataset is available through NCBI GEO accession GSE52058.

A long-term demasculinization of X-linked intergenic noncoding RNAs in Drosophila melanogaster.

Recent studies have revealed key roles of noncoding RNAs in sex-related pathways, but little is known about the evolutionary forces acting on these noncoding RNAs. Profiling the transcriptome of Drosophila melanogaster with whole-genome tiling arrays found that 15% of male-biased transcribed fragments are intergenic noncoding RNAs (incRNAs), suggesting a potentially important role for incRNAs in sex-related biological processes. Statistical analysis revealed a paucity of male-biased incRNAs and coding genes on the X chromosome, suggesting that similar evolutionary forces could be affecting the genomic organization of both coding and noncoding genes. Expression profiling across germline and somatic tissues further suggested that both male meiotic sex chromosome inactivation (MSCI) and sexual antagonism could contribute to the chromosomal distribution of male-biased incRNAs. Comparative sequence analysis showed that the evolutionary age of male-biased incRNAs is a significant predictor of their chromosomal locations. In addition to identifying abundant sex-biased incRNAs in the fly genome, our work unveils a global picture of the complex interplay between noncoding RNAs and sexual chromosome evolution.

New gene evolution: little did we know.

Genes are perpetually added to and deleted from genomes during evolution. Thus, it is important to understand how new genes are formed and how they evolve to be critical components of the genetic systems that determine the biological diversity of life. Two decades of effort have shed light on the process of new gene origination and have contributed to an emerging comprehensive picture of how new genes are added to genomes, ranging from the mechanisms that generate new gene structures to the presence of new genes in different organisms to the rates and patterns of new gene origination and the roles of new genes in phenotypic evolution. We review each of these aspects of new gene evolution, summarizing the main evidence for the origination and importance of new genes in evolution. We highlight findings showing that new genes rapidly change existing genetic systems that govern various molecular, cellular, and phenotypic functions.

Segmental dataset and whole body expression data do not support the hypothesis that non-random movement is an intrinsic property of Drosophila retrogenes.

BACKGROUND: Several studies in Drosophila have shown excessive movement of retrogenes from the X chromosome to autosomes, and that these genes are frequently expressed in the testis. This phenomenon has led to several hypotheses invoking natural selection as the process driving male-biased genes to the autosomes. Metta and Schlötterer (BMC Evol Biol 2010, 10:114) analyzed a set of retrogenes where the parental gene has been subsequently lost. They assumed that this class of retrogenes replaced the ancestral functions of the parental gene, and reported that these retrogenes, although mostly originating from movement out of the X chromosome, showed female-biased or unbiased expression. These observations led the authors to suggest that selective forces (such as meiotic sex chromosome inactivation and sexual antagonism) were not responsible for the observed pattern of retrogene movement out of the X chromosome. RESULTS: We reanalyzed the dataset published by Metta and Schlötterer and found several issues that led us to a different conclusion. In particular, Metta and Schlötterer used a dataset combined with expression data in which significant sex-biased expression is not detectable. First, the authors used a segmental dataset where the genes selected for analysis were less testis-biased in expression than those that were excluded from the study. Second, sex-biased expression was defined by comparing male and female whole-body data and not the expression of these genes in gonadal tissues. This approach significantly reduces the probability of detecting sex-biased expressed genes, which explains why the vast majority of the genes analyzed (parental and retrogenes) were equally expressed in both males and females. Third, the female-biased expression observed by Metta and Schlötterer is mostly found for parental genes located on the X chromosome, which is known to be enriched with genes with female-biased expression. Fourth, using additional gonad expression data, we found that autosomal genes analyzed by Metta and Schlötterer are less up regulated in ovaries and have higher chance to be expressed in meiotic cells of spermatogenesis when compared to X-linked genes. CONCLUSIONS: The criteria used to select retrogenes and the sex-biased expression data based on whole adult flies generated a segmental dataset of female-biased and unbiased expressed genes that was unable to detect the higher propensity of autosomal retrogenes to be expressed in males. Thus, there is no support for the authors' view that the movement of new retrogenes, which originated from X-linked parental genes, was not driven by selection. Therefore, selection-based genetic models remain the most parsimonious explanations for the observed chromosomal distribution of retrogenes.

Re-analysis of the larval testis data on meiotic sex chromosome inactivation revealed evidence for tissue-specific gene expression related to the drosophila X chromosome.

BACKGROUND: Meiotic sex chromosome inactivation (MSCI) during spermatogenesis has been proposed as one of the evolutionary driving forces behind both the under-representation of male-biased genes on, and the gene movement out of, the X chromosome in Drosophila. However, the relevance of MSCI in shaping sex chromosome evolution is controversial. Here we examine two aspects of a recent study on testis gene expression (Mikhaylova and Nurminsky, BMC Biol 2011, 9:29) that failed to support the MSCI in Drosophila. First, Mikhaylova and Nurminsky found no differences between X-linked and autosomal genes based on the transcriptional profiling of the early testis development, and thus concluded that MSCI does not occur in D. melanogaster. Second, they also analyzed expression data from several D. melanogaster tissues and concluded that under-representation on the X chromosome is not an exclusive property of testis-biased genes, but instead, a general property of tissue-specific genes. RESULTS: By re-analyzing the Mikhaylova and Nurminsky's testis data and the expression data on several D. melanogaster tissues, we made two major findings that refuted their original claims. First, the developmental testis data has generally greater experimental error than conventional analyses, which reduced significantly the power to detect chromosomal differences in expression. Nevertheless, our re-analysis observed significantly lower expression of the X chromosome in the genomic transcriptomes of later development stages of the testis, which is consistent with the MSCI hypothesis. Second, tissue-specific genes are also in general enriched with genes more expressed in testes than in ovaries, that is testis-biased genes. By completely excluding from the analyses the testis-biased genes, which are known to be under-represented in the X, we found that all the other tissue-specific genes are randomly distributed between the X chromosome and the autosomes. CONCLUSIONS: Our findings negate the original study of Mikhaylova and Nurminsky, which concluded a lack of MSCI and generalized the pattern of paucity in the X chromosome for tissue-specific genes in Drosophila. Therefore, MSCI and other selection-based models such as sexual antagonism, dosage compensation, and meiotic-drive continue to be viable models as driving forces shaping the genomic distribution of male-related genes in Drosophila.

Reshaping of global gene expression networks and sex-biased gene expression by integration of a young gene.

New genes originate frequently across diverse taxa. Given that genetic networks are typically comprised of robust, co-evolved interactions, the emergence of new genes raises an intriguing question: how do new genes interact with pre-existing genes? Here, we show that a recently originated gene rapidly evolved new gene networks and impacted sex-biased gene expression in Drosophila. This 4-6 million-year-old factor, named Zeus for its role in male fecundity, originated through retroposition of a highly conserved housekeeping gene, Caf40. Zeus acquired male reproductive organ expression patterns and phenotypes. Comparative expression profiling of mutants and closely related species revealed that Zeus has recruited a new set of downstream genes, and shaped the evolution of gene expression in germline. Comparative ChIP-chip revealed that the genomic binding profile of Zeus diverged rapidly from Caf40. These data demonstrate, for the first time, how a new gene quickly evolved novel networks governing essential biological processes at the genomic level.

Retrogenes moved out of the z chromosome in the silkworm.

Previous studies on organisms with well-differentiated X and Y chromosomes, such as Drosophila and mammals, consistently detected an excess of genes moving out of the X chromosome and gaining testis-biased expression. Several selective evolutionary mechanisms were shown to be associated with this nonrandom gene traffic, which contributed to the evolution of the X chromosome and autosomes. If selection drives gene traffic, such traffic should also exist in species with Z and W chromosomes, where the females are the heterogametic sex. However, no previous studies on gene traffic in species with female heterogamety have found any nonrandom chromosomal gene movement. Here, we report an excess of retrogenes moving out of the Z chromosome in an organism with the ZW sex determination system, Bombyx mori. In addition, we showed that those "out of Z" retrogenes tended to have ovary-biased expression, which is consistent with the pattern of non-retrogene traffic recently reported in birds and symmetrical to the retrogene movement in mammals and fruit flies out of the X chromosome evolving testis functions. These properties of gene traffic in the ZW system suggest a general role for the heterogamety of sex chromosomes in determining the chromosomal locations and the evolution of sex-biased genes.

Accelerated recruitment of new brain development genes into the human genome.

How the human brain evolved has attracted tremendous interests for decades. Motivated by case studies of primate-specific genes implicated in brain function, we examined whether or not the young genes, those emerging genome-wide in the lineages specific to the primates or rodents, showed distinct spatial and temporal patterns of transcription compared to old genes, which had existed before primate and rodent split. We found consistent patterns across different sources of expression data: there is a significantly larger proportion of young genes expressed in the fetal or infant brain of humans than in mouse, and more young genes in humans have expression biased toward early developing brains than old genes. Most of these young genes are expressed in the evolutionarily newest part of human brain, the neocortex. Remarkably, we also identified a number of human-specific genes which are expressed in the prefrontal cortex, which is implicated in complex cognitive behaviors. The young genes upregulated in the early developing human brain play diverse functional roles, with a significant enrichment of transcription factors. Genes originating from different mechanisms show a similar expression bias in the developing brain. Moreover, we found that the young genes upregulated in early brain development showed rapid protein evolution compared to old genes also expressed in the fetal brain. Strikingly, genes expressed in the neocortex arose soon after its morphological origin. These four lines of evidence suggest that positive selection for brain function may have contributed to the origination of young genes expressed in the developing brain. These data demonstrate a striking recruitment of new genes into the early development of the human brain.

Evolutionary patterns of RNA-based duplication in non-mammalian chordates.

The role of RNA-based duplication, or retroposition, in the evolution of new gene functions in mammals, plants, and Drosophila has been widely reported. However, little is known about RNA-based duplication in non-mammalian chordates. In this study, we screened ten non-mammalian chordate genomes for retrocopies and investigated their evolutionary patterns. We identified numerous retrocopies in these species. Examination of the age distribution of these retrocopies revealed no burst of young retrocopies in ancient chordate species. Upon comparing these non-mammalian chordate species to the mammalian species, we observed that a larger fraction of the non-mammalian retrocopies was under strong evolutionary constraints than mammalian retrocopies are, as evidenced by signals of purifying selection and expression profiles. For the Western clawed frog, Medaka, and Sea squirt, many retrogenes have evolved gonad and brain expression patterns, similar to what was observed in human. Testing of retrogene movement in the Medaka genome, where the nascent sex chrosomes have been well assembled, did not reveal any significant gene movement. Taken together, our analyses demonstrate that RNA-based duplication generates many functional genes and can make a significant contribution to the evolution of non-mammalian genomes.

A cautionary note for retrocopy identification: DNA-based duplication of intron-containing genes significantly contributes to the origination of single exon genes.

MOTIVATION: Retrocopies are important genes in the genomes of almost all higher eukaryotes. However, the annotation of such genes is a non-trivial task. Intronless genes have often been considered to be retroposed copies of intron-containing paralogs. Such categorization relies on the implicit premise that alignable regions of the duplicates should be long enough to cover exon-exon junctions of the intron-containing genes, and thus intron loss events can be inferred. Here, we examined the alternative possibility that intronless genes could be generated by partial DNA-based duplication of intron-containing genes in the fruitfly genome. RESULTS: By building pairwise protein-, transcript- and genome-level DNA alignments between intronless genes and their corresponding intron-containing paralogs, we found that alignments do not cover exon-exon junctions in 40% of cases and thus no intron loss could be inferred. For these cases, the candidate parental proteins tend to be partially duplicated, and intergenic sequences or neighboring genes are included in the intronless paralog. Moreover, we observed that it is significantly less likely for these paralogs to show inter-chromosomal duplication and testis-dominant transcription, compared to the remaining 60% of cases with evidence of clear intron loss (retrogenes). These lines of analysis reveal that DNA-based duplication contributes significantly to the 40% of cases of single exon gene duplication. Finally, we performed an analogous survey in the human genome and the result is similar, wherein 34% of the cases do not cover exon-exon junctions. Thus, genome annotation for retrogene identification should discard candidates without clear evidence of intron loss. CONTACT: mlong@uchicago.edu; zhangy@uchicago.edu

Chromosomal redistribution of male-biased genes in mammalian evolution with two bursts of gene gain on the X chromosome.

Mammalian X chromosomes evolved under various mechanisms including sexual antagonism, the faster-X process, and meiotic sex chromosome inactivation (MSCI). These forces may contribute to nonrandom chromosomal distribution of sex-biased genes. In order to understand the evolution of gene content on the X chromosome and autosome under these forces, we dated human and mouse protein-coding genes and miRNA genes on the vertebrate phylogenetic tree. We found that the X chromosome recently acquired a burst of young male-biased genes, which is consistent with fixation of recessive male-beneficial alleles by sexual antagonism. For genes originating earlier, however, this pattern diminishes and finally reverses with an overrepresentation of the oldest male-biased genes on autosomes. MSCI contributes to this dynamic since it silences X-linked old genes but not X-linked young genes. This demasculinization process seems to be associated with feminization of the X chromosome with more X-linked old genes expressed in ovaries. Moreover, we detected another burst of gene originations after the split of eutherian mammals and opossum, and these genes were quickly incorporated into transcriptional networks of multiple tissues. Preexisting X-linked genes also show significantly higher protein-level evolution during this period compared to autosomal genes, suggesting positive selection accompanied the early evolution of mammalian X chromosomes. These two findings cast new light on the evolutionary history of the mammalian X chromosome in terms of gene gain, sequence, and expressional evolution.