Large scale assessment of regulatory evolution and transcriptome complexity in mammals

AbstractIn addition to genetic changes affecting the function of gene products, changes in gene expression have been suggested to underlie many or even most of the phenotypic differences among mammals. However, detailed gene expression comparisons were, until recently, restricted to closely related species, owing to technological limitations. Thus, we took advantage of the latest technologies (RNA-Seq) to generate extensive qualitative and quantitative transcriptome data for a unique collection of somatic and germline tissues from representatives of all major mammalian lineages (placental mammals, marsupials and monotremes) and birds, the evolutionary outgroup.In the first major project of my thesis, we performed global comparative analyses of gene expression levels based on these data. Our analyses provided fundamental insights into the dynamics of transcriptome change during mammalian evolution (e.g., the rate of expression change across species, tissues and chromosomes) and allowed the exploration of the functional relevance and phenotypic implications of transcription changes at a genome-wide scale (e.g., we identified numerous potentially selectively driven expression switches).In a second project of my thesis, which was also based on the unique transcriptome data generated in the context of the first project we focused on the evolution of alternative splicing in mammals. Alternative splicing contributes to transcriptome complexity by generating several transcript isoforms from a single gene, which can, thus, perform various functions. To complete the global comparative analysis of gene expression changes, we explored patterns of alternative splicing evolution. This work uncovered several general and unexpected patterns of alternative splicing evolution (e.g., we found that alternative splicing evolves extremely rapidly) as well as a large number of conserved alternative isoforms that may be crucial for the functioning of mammalian organs.Finally, the third and final project of my PhD consisted in analyzing in detail the unique functional and evolutionary properties of the testis by exploring the extent of its transcriptome complexity. This organ was previously shown to evolve rapidly both at the phenotypic and molecular level, apparently because of the specific pressures that act on this organ and are associated with its reproductive function. Moreover, my analyses of the amniote tissue transcriptome data described above, revealed strikingly widespread transcriptional activity of both functional and nonfunctional genomic elements in the testis compared to the other organs. To elucidate the cellular source and mechanisms underlying this promiscuous transcription in the testis, we generated deep coverage RNA-Seq data for all major testis cell types as well as epigenetic data (DNA and histone methylation) using the mouse as model system. The integration of these complete dataset revealed that meiotic and especially post-meiotic germ cells are the major contributors to the widespread functional and nonfunctional transcriptome complexity of the testis, and that this "promiscuous" spermatogenic transcription is resulting, at least partially, from an overall transcriptionally permissive chromatin state. We hypothesize that this particular open state of the chromatin results from the extensive chromatin remodeling that occurs during spermatogenesis which ultimately leads to the replacement of histones by protamines in the mature spermatozoa. Our results have important functional and evolutionary implications (e.g., regarding new gene birth and testicular gene expression evolution).Generally, these three large-scale projects of my thesis provide complete and massive datasets that constitute valuables resources for further functional and evolutionary analyses of mammalian genomes.

Conserved microRNA editing in mammalian evolution, development and disease.

BACKGROUND: Mammalian microRNAs (miRNAs) are sometimes subject to adenosine-to-inosine RNA editing, which can lead to dramatic changes in miRNA target specificity or expression levels. However, although a few miRNAs are known to be edited at identical positions in human and mouse, the evolution of miRNA editing has not been investigated in detail. In this study, we identify conserved miRNA editing events in a range of mammalian and non-mammalian species. RESULTS: We demonstrate deep conservation of several site-specific miRNA editing events, including two that date back to the common ancestor of mammals and bony fishes some 450 million years ago. We also find evidence of a recent expansion of an edited miRNA family in placental mammals and show that editing of these miRNAs is associated with changes in target mRNA expression during primate development and aging. While global patterns of miRNA editing tend to be conserved across species, we observe substantial variation in editing frequencies depending on tissue, age and disease state: editing is more frequent in neural tissues compared to heart, kidney and testis; in older compared to younger individuals; and in samples from healthy tissues compared to tumors, which together suggests that miRNA editing might be associated with a reduced rate of cell proliferation. CONCLUSIONS: Our results show that site-specific miRNA editing is an evolutionarily conserved mechanism, which increases the functional diversity of mammalian miRNA transcriptomes. Furthermore, we find that although miRNA editing is rare compared to editing of long RNAs, miRNAs are greatly overrepresented among conserved editing targets.

Mechanisms and evolutionary patterns of mammalian and avian dosage compensation.

As a result of sex chromosome differentiation from ancestral autosomes, male mammalian cells only contain one X chromosome. It has long been hypothesized that X-linked gene expression levels have become doubled in males to restore the original transcriptional output, and that the resulting X overexpression in females then drove the evolution of X inactivation (XCI). However, this model has never been directly tested and patterns and mechanisms of dosage compensation across different mammals and birds generally remain little understood. Here we trace the evolution of dosage compensation using extensive transcriptome data from males and females representing all major mammalian lineages and birds. Our analyses suggest that the X has become globally upregulated in marsupials, whereas we do not detect a global upregulation of this chromosome in placental mammals. However, we find that a subset of autosomal genes interacting with X-linked genes have become downregulated in placentals upon the emergence of sex chromosomes. Thus, different driving forces may underlie the evolution of XCI and the highly efficient equilibration of X expression levels between the sexes observed for both of these lineages. In the egg-laying monotremes and birds, which have partially homologous sex chromosome systems, partial upregulation of the X (Z in birds) evolved but is largely restricted to the heterogametic sex, which provides an explanation for the partially sex-biased X (Z) expression and lack of global inactivation mechanisms in these lineages. Our findings suggest that dosage reductions imposed by sex chromosome differentiation events in amniotes were resolved in strikingly different ways.

Evolution of mRNA expression levels of autosomal and sex-linked genes in amniotes

In addition to differences in protein-coding gene sequences, changes in expression resulting from mutations in regulatory sequences have long been hypothesized to be responsible for phenotypic differences between species. However, unlike comparison of genome sequences, few studies, generally restricted to pairwise comparisons of closely related mammalian species, have assessed between-species differences at the transcriptome level. They reported that gene expression evolves at different rates in various organs and in a pattern that is overall consistent with neutral models of evolution. In the first part of my thesis, I investigated the evolution of gene expression in therian mammals (i.e.7 placental and marsupials), based on microarray data from human, mouse and the gray short-tailed opossum (Monodelphis domestica). In addition to autosomal genes, a special focus was given to the evolution of X-linked genes. The therian X chromosome was recently shown to be younger than previously thought and to harbor a specific gene content (e.g., genes involved in brain or reproductive functions) that is thought to have been shaped by specific sex-related evolutionary forces. Sex chromosomes derive from ordinary autosomes and their differentiation led to the degeneration of the Y chromosome (in mammals) or W chromosome (in birds). Consequently, X- or Z-linked genes differ in gene dose between males and females such that the heterogametic sex has half the X/Z gene dose compared to the ancestral state. To cope with this dosage imbalance, mammals have been reported to have evolved mechanisms of dosage compensation.¦In the first project, I could first show that transcriptomes evolve at different rates in different organs. Out of the five tissues I investigated, the testis is the most rapidly evolving organ at the gene expression level while the brain has the most conserved transcriptome. Second, my analyses revealed that mammalian gene expression evolution is compatible with a neutral model, where the rates of change in gene expression levels is linked to the efficiency of purifying selection in a given lineage, which, in turn, is determined by the long-term effective population size in that lineage. Thus, the rate of DNA sequence evolution, which could be expected to determine the rate of regulatory sequence change, does not seem to be a major determinant of the rate of gene expression evolution. Thus, most gene expression changes seem to be (slightly) deleterious. Finally, X-linked genes seem to have experienced elevated rates of gene expression change during the early stage of X evolution. To further investigate the evolution of mammalian gene expression, we generated an extensive RNA-Seq gene expression dataset for nine mammalian species and a bird. The analyses of this dataset confirmed the patterns previously observed with microarrays and helped to significantly deepen our view on gene expression evolution.¦In a specific project based on these data, I sought to assess in detail patterns of evolution of dosage compensation in amniotes. My analyses revealed the absence of male to female dosage compensation in monotremes and its presence in marsupials and, in addition, confirmed patterns previously described for placental mammals and birds. I then assessed the global level of expression of X/Z chromosomes and contrasted this with its ancestral gene expression levels estimated from orthologous autosomal genes in species with non-homologous sex chromosomes. This analysis revealed a lack of up-regulation for placental mammals, the level of expression of X-linked genes being proportional to gene dose. Interestingly, the ancestral gene expression level was at least partially restored in marsupials as well as in the heterogametic sex of monotremes and birds. Finally, I investigated alternative mechanisms of dosage compensation and found that gene duplication did not seem to be a widespread mechanism to restore the ancestral gene dose. However, I could show that placental mammals have preferentially down-regulated autosomal genes interacting with X-linked genes which underwent gene expression decrease, and thus identified a novel alternative mechanism of dosage compensation.

The current molecular phylogeny of Eutherian mammals challenges previous interpretations of placental evolution.

Based on histology, the placentae of eutherians are currently grouped in epitheliochorial, endotheliochorial and haemochorial placentae. In a haeckelian sense, the epitheliochorial contact with marked histiotrophic feeding by uterine milk is generally considered as primitive, especially since similar contacts exist in Marsupials. In contrast, the more intimate endotheliochorial and haemochorial contact, facilitating haemotrophic nutrition, is interpreted as a derived state. A cladistic analysis based on the phylogenetic relationships established by molecular analyses reveals that the basic clades are all characterized by an endotheliochorial or haemochorial placenta, and that the epitheliochorial placenta evolved at least three times in a convergent manner. This evolution may be explained by the fact that the epitheliochorial placenta in eutherians is more efficient in nutritional transfer (flow rate by exchange surface). Moreover, this arrangement may confer an advantage to the mother who can probably reduce the degree of manipulation by a genetically imprinted embryo.

Origins and functional evolution of Y chromosomes across mammals.

Y chromosomes underlie sex determination in mammals, but their repeat-rich nature has hampered sequencing and associated evolutionary studies. Here we trace Y evolution across 15 representative mammals on the basis of high-throughput genome and transcriptome sequencing. We uncover three independent sex chromosome originations in mammals and birds (the outgroup). The original placental and marsupial (therian) Y, containing the sex-determining gene SRY, emerged in the therian ancestor approximately 180 million years ago, in parallel with the first of five monotreme Y chromosomes, carrying the probable sex-determining gene AMH. The avian W chromosome arose approximately 140 million years ago in the bird ancestor. The small Y/W gene repertoires, enriched in regulatory functions, were rapidly defined following stratification (recombination arrest) and erosion events and have remained considerably stable. Despite expression decreases in therians, Y/W genes show notable conservation of proto-sex chromosome expression patterns, although various Y genes evolved testis-specificities through differential regulatory decay. Thus, although some genes evolved novel functions through spatial/temporal expression shifts, most Y genes probably endured, at least initially, because of dosage constraints.

Species identification in mammals from mixed biological samples based on mitochondrial DNA control region length polymorphism.

The ability to identify the species origin of an unknown biological sample is relevant in the fields of human and wildlife forensics. However, the detection of several species mixed in the same sample still remains a challenge. We developed and tested a new approach for mammal DNA identification in mixtures of two or three species, based on the analysis of mitochondrial DNA control region interspecific length polymorphism followed by direct sequencing. Contrary to other published methods dealing with species mixtures, our protocol requires a single universal primer pair and is not based on a pre-defined panel of species. Amplicons can be separated either on agarose gels or using CE. The advantages and limitations of the assay are discussed under different conditions, such as variable template concentration, amplicon sizes and size difference among the amplicons present in the mixture. For the first time, this protocol provides a simple, reliable and flexible method for simultaneous identification of multiple mammalian species from mixtures, without any prior knowledge of the species involved.

Immunoglobulin classes in aquatic mammals. Characterization by serologic cross-reactivity, molecular size and binding of human free secretory component.

On the basis of serologic cross-reactivity, three immunoglobulin classes homologous to human IgG, IgM and IgA were identified in two species of acquatic mammal representing the orders Cetacea (dolphin) and Pinnipedea (sea lion). Molecular size was estimated by sucrose density gradient ultracentrifugation and Sephadex G-200 chromatography, indicating a 7S IgG, 19S IgM and heterogeneous serum IgA. Human secretory component was readily bound to the IgM of both species and to an apparently lesser extent to the larger molecular size populations of IgA. No binding was observed with IgG. Several antisera specific for human γ-chains gave a single precipitin line with the sea lion IgG but when made to react with dolphin serum produced two lines, suggesting the presence of two different subclasses of IgG in this species.

Chromosomal gene movements reflect the recent origin and biology of therian sex chromosomes

Mammalian sex chromosomes stem from ancestral autosomes and have substantially differentiated. It was shown that X-linked genes have generated duplicate intronless gene copies (retrogenes) on autosomes due to this differentiation. However, the precise driving forces for this out-of-X gene "movement" and its evolutionary onset are not known. Based on expression analyses of male germ-cell populations, we here substantiate and extend the hypothesis that autosomal retrogenes functionally compensate for the silencing of their X-linked housekeeping parental genes during, but also after, male meiotic sex chromosome inactivation (MSCI). Thus, sexually antagonistic forces have not played a major role for the selective fixation of X-derived gene copies in mammals. Our dating analyses reveal that although retrogenes were produced ever since the common mammalian ancestor, selectively driven retrogene export from the X only started later, on the placental mammal (eutherian) and marsupial (metatherian) lineages, respectively. Together, these observations suggest that chromosome-wide MSCI emerged close to the eutherian-marsupial split approximately 180 million years ago. Given that MSCI probably reflects the spread of the recombination barrier between the X and Y, crucial for their differentiation, our data imply that these chromosomes became more widely differentiated only late in the therian ancestor, well after the divergence of the monotreme lineage. Thus, our study also provides strong independent support for the recent notion that our sex chromosomes emerged, not in the common ancestor of all mammals, but rather in the therian ancestor, and therefore are much younger than previously thought

Plutonium from above-ground nuclear tests in milk teeth: investigation of placental transfer in children born between 1951 and 1995 in Switzerland.

BACKGROUND: Occupational risks, the present nuclear threat, and the potential danger associated with nuclear power have raised concerns regarding the metabolism of plutonium in pregnant women. OBJECTIVE: We measured plutonium levels in the milk teeth of children born between 1951 and 1995 to assess the potential risk that plutonium incorporated by pregnant women might pose to the radiosensitive tissues of the fetus through placenta transfer. METHODS: We used milk teeth, whose enamel is formed during pregnancy, to investigate the transfer of plutonium from the mother's blood plasma to the fetus. We measured plutonium using sensitive sector field inductively coupled plasma mass spectrometry techniques. We compared our results with those of a previous study on strontium-90 ((90)Sr) released into the atmosphere after nuclear bomb tests. RESULTS: Results show that plutonium activity peaks in the milk teeth of children born about 10 years before the highest recorded levels of plutonium fallout. By contrast, (90)Sr, which is known to cross the placenta barrier, manifests differently in milk teeth, in accordance with (90)Sr fallout deposition as a function of time. CONCLUSIONS: These findings demonstrate that plutonium found in milk teeth is caused by fallout that was inhaled around the time the milk teeth were shed and not from any accumulation during pregnancy through placenta transfer. Thus, plutonium may not represent a radiologic risk for the radiosensitive tissues of the fetus.

Detection of suspected placental invasion by MRI: Do the results depend on observer' experience?

PURPOSE: To evaluate the diagnostic value of previously described MR features used for detecting suspected placental invasion according to observers' experience. MATERIALS AND METHODS: Our population included 25 pregnant women (mean age 35.16) investigated by prenatal MRI (1.5T, T1- and T2-weighted MR-sequences without i.v. contrast), among them 12 with histopathologically proven placental invasion and 13 women (52%) without placental invasion used as control group. Two senior and two junior radiologists blindly and independently reviewed MR-examinations in view of 6 previously defined MR-features indicating presence and degree of placental invasion (placenta increta, accreta or percreta). For each reader the sensibility, specificity, and receiver operating curve (ROC) were calculated. Interobserver agreements between senior and junior readers were determined. Stepwise logistic regression was performed including the 6 MR-features predictive of placental invasion. RESULTS: Demographics between both groups were statistically equivalent. Overall sensitivity and specificity for placental invasion was 90.9% and 75.0% for seniors and 81.8% and 61.8% for juniors, respectively. The best single MR-feature indicating placental invasion was T2-hypointense placental bands (r(2)=0.28), followed by focally interrupted myometrial border, infiltration of pelvic organs and tenting of the bladder (r(2)=0.36). Interobserver agreement for detecting placental invasion was 0.64 for seniors and 0.41 for juniors, thus substantial and moderate, respectively. Seniors detected placental invasion and depth of infiltration with significantly higher diagnostic certitude than juniors (p=0.0002 and p=0.0282, respectively). CONCLUSION: MRI can be a reliable and reproducible tool for the detection of suspected placental invasion, but the diagnostic value significantly depends on observers' experience.

The emergence of new genes on the young therian X.

Recent studies have revealed that our sex chromosomes differentiated relatively recently from ancestral autosomes in the common ancestor of placental and marsupial mammals (therians). Here, we show that the therian X started to accumulate new retroduplicate genes with overall sex-biased expression upon therian sex chromosome differentiation. This process reached its peak within the first approximately 90 million years of sex chromosome evolution and then leveled off. Taken together, our observations suggest that the major sex-related functional remodeling of the X was completed relatively soon after the origination of therian sex chromosomes.

Role of Waddlia chondrophila placental infection in miscarriage.

Waddlia chondrophila is an intracellular bacterium suspected to cause human and bovine abortion. We confirmed an association between antibodies against W. chondrophila and human miscarriage and identified this organism in placenta or genital tract of women who had had miscarriages. These results suggest a possible role of W. chondrophila infection in miscarriage.