Three dimensions of thermolabile sex determination.

The molecular mechanism of temperature-dependent sex determination (TSD) is a long-standing mystery. How is the thermal signal sensed, captured and transduced to regulate key sex genes? Although there is compelling evidence for pathways via which cells capture the temperature signal, there is no known mechanism by which cells transduce those thermal signals to affect gene expression. Here we propose a novel hypothesis we call 3D-TSD (the three dimensions of thermolabile sex determination). We postulate that the genome has capacity to remodel in response to temperature by changing 3D chromatin conformation, perhaps via temperature-sensitive transcriptional condensates. This could rewire enhancer-promoter interactions to alter the expression of key sex-determining genes. This hypothesis can accommodate monogenic or multigenic thermolabile sex-determining systems, and could be combined with upstream thermal sensing and transduction to the epigenome to commit gonadal fate.

Strategies for meiotic sex chromosome dynamics and telomeric elongation in Marsupials.

During meiotic prophase I, homologous chromosomes pair, synapse and recombine in a tightly regulated process that ensures the generation of genetically variable haploid gametes. Although the mechanisms underlying meiotic cell division have been well studied in model species, our understanding of the dynamics of meiotic prophase I in non-traditional model mammals remains in its infancy. Here, we reveal key meiotic features in previously uncharacterised marsupial species (the tammar wallaby and the fat-tailed dunnart), plus the fat-tailed mouse opossum, with a focus on sex chromosome pairing strategies, recombination and meiotic telomere homeostasis. We uncovered differences between phylogroups with important functional and evolutionary implications. First, sex chromosomes, which lack a pseudo-autosomal region in marsupials, had species specific pairing and silencing strategies, with implications for sex chromosome evolution. Second, we detected two waves of γH2AX accumulation during prophase I. The first wave was accompanied by low γH2AX levels on autosomes, which correlated with the low recombination rates that distinguish marsupials from eutherian mammals. In the second wave, γH2AX was restricted to sex chromosomes in all three species, which correlated with transcription from the X in tammar wallaby. This suggests non-canonical functions of γH2AX on meiotic sex chromosomes. Finally, we uncover evidence for telomere elongation in primary spermatocytes of the fat-tailed dunnart, a unique strategy within mammals. Our results provide new insights into meiotic progression and telomere homeostasis in marsupials, highlighting the importance of capturing the diversity of meiotic strategies within mammals.

Microchromosomes are building blocks of bird, reptile, and mammal chromosomes.

Microchromosomes, once considered unimportant shreds of the chicken genome, are gene-rich elements with a high GC content and few transposable elements. Their origin has been debated for decades. We used cytological and whole-genome sequence comparisons, and chromosome conformation capture, to trace their origin and fate in genomes of reptiles, birds, and mammals. We find that microchromosomes as well as macrochromosomes are highly conserved across birds and share synteny with single small chromosomes of the chordate amphioxus, attesting to their origin as elements of an ancient animal genome. Turtles and squamates (snakes and lizards) share different subsets of ancestral microchromosomes, having independently lost microchromosomes by fusion with other microchromosomes or macrochromosomes. Patterns of fusions were quite different in different lineages. Cytological observations show that microchromosomes in all lineages are spatially separated into a central compartment at interphase and during mitosis and meiosis. This reflects higher interaction between microchromosomes than with macrochromosomes, as observed by chromosome conformation capture, and suggests some functional coherence. In highly rearranged genomes fused microchromosomes retain most ancestral characteristics, but these may erode over evolutionary time; surprisingly, de novo microchromosomes have rapidly adopted high interaction. Some chromosomes of early-branching monotreme mammals align to several bird microchromosomes, suggesting multiple microchromosome fusions in a mammalian ancestor. Subsequently, multiple rearrangements fueled the extraordinary karyotypic diversity of therian mammals. Thus, microchromosomes, far from being aberrant genetic elements, represent fundamental building blocks of amniote chromosomes, and it is mammals, rather than reptiles and birds, that are atypical.

Meiotic Executioner Genes Protect the Y from Extinction.

The Y has been described as a wimpy degraded relic of the X, with imminent demise should it lose sex-determining function. Why then has it persisted in almost all mammals? Here we present a novel mechanistic explanation for its evolutionary perseverance: the persistent Y hypothesis. The Y chromosome bears genes that act as their own judge, jury, and executioner in the tightly regulated meiotic surveillance pathways. These executioners are crucial for successful meiosis, yet need to be silenced during the meiotic sex chromosome inactivation window, otherwise germ cells die. Only rare transposition events to the X, where they remain subject to obligate meiotic silencing, are heritable, posing strong evolutionary constraint for the Y chromosome to persist.

Fragile, unfaithful and persistent Ys-on how meiosis can shape sex chromosome evolution.

Sex-linked inheritance is a stark exception to Mendel's Laws of Heredity. Here we discuss how the evolution of heteromorphic sex chromosomes (mainly the Y) has been shaped by the intricacies of the meiotic programme. We propose that persistence of Y chromosomes in distantly related mammalian phylogroups can be explained in the context of pseudoautosomal region (PAR) size, meiotic pairing strategies, and the presence of Y-borne executioner genes that regulate meiotic sex chromosome inactivation. We hypothesise that variation in PAR size can be an important driver for the evolution of recombination frequencies genome wide, imposing constraints on Y fate. If small PAR size compromises XY segregation during male meiosis, the stress of producing aneuploid gametes could drive function away from the Y (i.e., a fragile Y). The Y chromosome can avoid fragility either by acquiring an achiasmatic meiotic XY pairing strategy to reduce aneuploid gamete production, or gain meiotic executioner protection (a persistent Y). Persistent Ys will then be under strong pressure to maintain high recombination rates in the PAR (and subsequently genome wide), as improper segregation has fatal consequences for germ cells. In the event that executioner protection is lost, the Y chromosome can be maintained in the population by either PAR rejuvenation (extension by addition of autosome material) or gaining achiasmatic meiotic pairing, the alternative is Y loss. Under this dynamic cyclic evolutionary scenario, understanding the meiotic programme in vertebrate and invertebrate species will be crucial to further understand the plasticity of the rise and fall of heteromorphic sex chromosomes.

Nonlinear sequence similarity between the Xist and Rsx long noncoding RNAs suggests shared functions of tandem repeat domains.

The marsupial inactive X chromosome expresses a long noncoding RNA (lncRNA) called Rsx that has been proposed to be the functional analog of eutherian Xist Despite the possibility that Xist and Rsx encode related functions, the two lncRNAs harbor no linear sequence similarity. However, both lncRNAs harbor domains of tandemly repeated sequence. In Xist, these repeat domains are known to be critical for function. Using k-mer based comparison, we show that the repeat domains of Xist and Rsx unexpectedly partition into two major clusters that each harbor substantial levels of nonlinear sequence similarity. Xist Repeats B, C, and D were most similar to each other and to Rsx Repeat 1, whereas Xist Repeats A and E were most similar to each other and to Rsx Repeats 2, 3, and 4. Similarities at the level of k-mers corresponded to domain-specific enrichment of protein-binding motifs. Within individual domains, protein-binding motifs were often enriched to extreme levels. Our data support the hypothesis that Xist and Rsx encode similar functions through different spatial arrangements of functionally analogous protein-binding domains. We propose that the two clusters of repeat domains in Xist and Rsx function in part to cooperatively recruit PRC1 and PRC2 to chromatin. The physical manner in which these domains engage with protein cofactors may be just as critical to the function of the domains as the protein cofactors themselves. The general approaches we outline in this report should prove useful in the study of any set of RNAs.

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.

Waking the sleeping dragon: gene expression profiling reveals adaptive strategies of the hibernating reptile Pogona vitticeps.

BACKGROUND: Hibernation is a physiological state exploited by many animals exposed to prolonged adverse environmental conditions associated with winter. Large changes in metabolism and cellular function occur, with many stress response pathways modulated to tolerate physiological challenges that might otherwise be lethal. Many studies have sought to elucidate the molecular mechanisms of mammalian hibernation, but detailed analyses are lacking in reptiles. Here we examine gene expression in the Australian central bearded dragon (Pogona vitticeps) using mRNA-seq and label-free quantitative mass spectrometry in matched brain, heart and skeletal muscle samples from animals at late hibernation, 2 days post-arousal and 2 months post-arousal. RESULTS: We identified differentially expressed genes in all tissues between hibernation and post-arousal time points; with 4264 differentially expressed genes in brain, 5340 differentially expressed genes in heart, and 5587 differentially expressed genes in skeletal muscle. Furthermore, we identified 2482 differentially expressed genes across all tissues. Proteomic analysis identified 743 proteins (58 differentially expressed) in brain, 535 (57 differentially expressed) in heart, and 337 (36 differentially expressed) in skeletal muscle. Tissue-specific analyses revealed enrichment of protective mechanisms in all tissues, including neuroprotective pathways in brain, cardiac hypertrophic processes in heart, and atrophy protective pathways in skeletal muscle. In all tissues stress response pathways were induced during hibernation, as well as evidence for gene expression regulation at transcription, translation and post-translation. CONCLUSIONS: These results reveal critical stress response pathways and protective mechanisms that allow for maintenance of both tissue-specific function, and survival during hibernation in the central bearded dragon. Furthermore, we provide evidence for multiple levels of gene expression regulation during hibernation, particularly enrichment of miRNA-mediated translational repression machinery; a process that would allow for rapid and energy efficient reactivation of translation from mature mRNA molecules at arousal. This study is the first molecular investigation of its kind in a hibernating reptile, and identifies strategies not yet observed in other hibernators to cope stress associated with this remarkable state of metabolic depression.

Landscape of DNA Methylation on the Marsupial X.

DNA methylation plays a key role in maintaining transcriptional silence on the inactive X chromosome of eutherian mammals. Beyond eutherians, there are limited genome wide data on DNA methylation from other vertebrates. Previous studies of X borne genes in various marsupial models revealed no differential DNA methylation of promoters between the sexes, leading to the conclusion that CpG methylation plays no role in marsupial X-inactivation. Using reduced representation bisulfite sequencing, we generated male and female CpG methylation profiles in four representative vertebrates (mouse, gray short-tailed opossum, platypus, and chicken). A variety of DNA methylation patterns were observed. Platypus and chicken displayed no large-scale differential DNA methylation between the sexes on the autosomes or the sex chromosomes. As expected, a metagene analysis revealed hypermethylation at transcription start sites (TSS) of genes subject to X-inactivation in female mice. This contrasted with the opossum, in which metagene analysis did not detect differential DNA methylation between the sexes at TSSs of genes subject to X-inactivation. However, regions flanking TSSs of these genes were hypomethylated. Our data are the first to demonstrate that, for genes subject to X-inactivation in both eutherian and marsupial mammals, there is a consistent difference between DNA methylation levels at TSSs and immediate flanking regions, which we propose has a silencing effect in both groups.

Comparative analysis of mammalian Y chromosomes illuminates ancestral structure and lineage-specific evolution.

Although more than thirty mammalian genomes have been sequenced to draft quality, very few of these include the Y chromosome. This has limited our understanding of the evolutionary dynamics of gene persistence and loss, our ability to identify conserved regulatory elements, as well our knowledge of the extent to which different types of selection act to maintain genes within this unique genomic environment. Here, we present the first MSY (male-specific region of the Y chromosome) sequences from two carnivores, the domestic dog and cat. By combining these with other available MSY data, our multiordinal comparison allows for the first accounting of levels of selection constraining the evolution of eutherian Y chromosomes. Despite gene gain and loss across the phylogeny, we show the eutherian ancestor retained a core set of 17 MSY genes, most being constrained by negative selection for nearly 100 million years. The X-degenerate and ampliconic gene classes are partitioned into distinct chromosomal domains in most mammals, but were radically restructured on the human lineage. We identified multiple conserved noncoding elements that potentially regulate eutherian MSY genes. The acquisition of novel ampliconic gene families was accompanied by signatures of positive selection and has differentially impacted the degeneration and expansion of MSY gene repertoires in different species.

Independent evolution of transcriptional inactivation on sex chromosomes in birds and mammals.

X chromosome inactivation in eutherian mammals has been thought to be tightly controlled, as expected from a mechanism that compensates for the different dosage of X-borne genes in XX females and XY males. However, many X genes escape inactivation in humans, inactivation of the X in marsupials is partial, and the unrelated sex chromosomes of monotreme mammals have incomplete and gene-specific inactivation of X-linked genes. The bird ZW sex chromosome system represents a third independently evolved amniote sex chromosome system with dosage compensation, albeit partial and gene-specific, via an unknown mechanism (i.e. upregulation of the single Z in females, down regulation of one or both Zs in males, or a combination). We used RNA-fluorescent in situ hybridization (RNA-FISH) to demonstrate, on individual fibroblast cells, inactivation of 11 genes on the chicken Z and 28 genes on the X chromosomes of platypus. Each gene displayed a reproducible frequency of 1Z/1X-active and 2Z/2X-active cells in the homogametic sex. Our results indicate that the probability of inactivation is controlled on a gene-by-gene basis (or small domains) on the chicken Z and platypus X chromosomes. This regulatory mechanism must have been exapted independently to the non-homologous sex chromosomes in birds and mammals in response to an over-expressed Z or X in the homogametic sex, highlighting the universal importance that (at least partial) silencing plays in the evolution on amniote dosage compensation and, therefore, the differentiation of sex chromosomes.

Long intervening non-coding RNA 00320 is human brain-specific and highly expressed in the cortical white matter.

Pervasive transcription of the genome produces a diverse array of functional non-coding RNAs (ncRNAs). One particular class of ncRNAs, long intervening non-coding RNAs (lincRNAs) are thought to play a role in regulating gene expression and may be a major contributor to organism and tissue complexity. The human brain with its heterogeneous cellular make-up is a rich source of lincRNAs; however, the functions of the majority of lincRNAs are unknown. Recently, by completing RNA sequencing (RNA-Seq) of the human frontal cortex, we identified linc00320 as being highly expressed in the white matter compared to grey matter in multiple system atrophy (MSA) brain. Here, we further investigate the expression patterns of linc00320 and conclude that it is involved in specific brain regions rather than having involvement in the MSA disease process. We also show that the full-length linc00320 is only expressed in human brain tissue and not in other primates, suggesting that it may be involved in improved functional connectivity for higher human brain cognition.

Evolutionary history of novel genes on the tammar wallaby Y chromosome: Implications for sex chromosome evolution.

We report here the isolation and sequencing of 10 Y-specific tammar wallaby (Macropus eugenii) BAC clones, revealing five hitherto undescribed tammar wallaby Y genes (in addition to the five genes already described) and several pseudogenes. Some genes on the wallaby Y display testis-specific expression, but most have low widespread expression. All have partners on the tammar X, along with homologs on the human X. Nonsynonymous and synonymous substitution ratios for nine of the tammar XY gene pairs indicate that they are each under purifying selection. All 10 were also identified as being on the Y in Tasmanian devil (Sarcophilus harrisii; a distantly related Australian marsupial); however, seven have been lost from the human Y. Maximum likelihood phylogenetic analyses of the wallaby YX genes, with respective homologs from other vertebrate representatives, revealed that three marsupial Y genes (HCFC1X/Y, MECP2X/Y, and HUWE1X/Y) were members of the ancestral therian pseudoautosomal region (PAR) at the time of the marsupial/eutherian split; three XY pairs (SOX3/SRY, RBMX/Y, and ATRX/Y) were isolated from each other before the marsupial/eutherian split, and the remaining three (RPL10X/Y, PHF6X/Y, and UBA1/UBE1Y) have a more complex evolutionary history. Thus, the small marsupial Y chromosome is surprisingly rich in ancient genes that are retained in at least Australian marsupials and evolved from testis-brain expressed genes on the X.

Paternal X inactivation does not correlate with X chromosome evolutionary strata in marsupials.

BACKGROUND: X chromosome inactivation is the transcriptional silencing of one X chromosome in the somatic cells of female mammals. In eutherian mammals (e.g. humans) one of the two X chromosomes is randomly chosen for silencing, with about 15% (usually in younger evolutionary strata of the X chromosome) of genes escaping this silencing. In contrast, in the distantly related marsupial mammals the paternally derived X is silenced, although not as completely as the eutherian X. A chromosome wide examination of X inactivation, using RNA-seq, was recently undertaken in grey short-tailed opossum (Monodelphis domestica) brain and extraembryonic tissues. However, no such study has been conduced in Australian marsupials, which diverged from their American cousins ~80 million years ago, leaving a large gap in our understanding of marsupial X inactivation. RESULTS: We used RNA-seq data from blood or liver of a family (mother, father and daughter) of tammar wallabies (Macropus eugenii), which in conjunction with available genome sequence from the mother and father, permitted genotyping of 42 expressed heterozygous SNPs on the daughter's X. These 42 SNPs represented 34 X loci, of which 68% (23 of the 34) were confirmed as inactivated on the paternally derived X in the daughter's liver; the remaining 11 X loci escaped inactivation. Seven of the wallaby loci sampled were part of the old X evolutionary stratum, of which three escaped inactivation. Three loci were classified as part of the newer X stratum, of which two escaped inactivation. A meta-analysis of previously published opossum X inactivation data revealed that 5 of 52 genes in the old X stratum escaped inactivation. CONCLUSIONS: We demonstrate that chromosome wide inactivation of the paternal X is common to an Australian marsupial representative, but that there is more escape from inactivation than reported for opossum (32% v 14%). We also provide evidence that, unlike the human X chromosome, the location of loci within the oldest evolutionary stratum on the marsupial X does not correlate with their probability of escape from inactivation.

Genome analysis of the platypus reveals unique signatures of evolution.

We present a draft genome sequence of the platypus, Ornithorhynchus anatinus. This monotreme exhibits a fascinating combination of reptilian and mammalian characters. For example, platypuses have a coat of fur adapted to an aquatic lifestyle; platypus females lactate, yet lay eggs; and males are equipped with venom similar to that of reptiles. Analysis of the first monotreme genome aligned these features with genetic innovations. We find that reptile and platypus venom proteins have been co-opted independently from the same gene families; milk protein genes are conserved despite platypuses laying eggs; and immune gene family expansions are directly related to platypus biology. Expansions of protein, non-protein-coding RNA and microRNA families, as well as repeat elements, are identified. Sequencing of this genome now provides a valuable resource for deep mammalian comparative analyses, as well as for monotreme biology and conservation.

Identification of the RSX interactome in a marsupial shows functional coherence with the Xist interactome during X inactivation.

The marsupial specific RSX lncRNA is the functional analogue of the eutherian specific XIST, which coordinates X chromosome inactivation. We characterized the RSX interactome in a marsupial representative (the opossum Monodelphis domestica), identifying 135 proteins, of which 54 had orthologues in the XIST interactome. Both interactomes were enriched for biological pathways related to RNA processing, regulation of translation, and epigenetic transcriptional silencing. This represents a remarkable example showcasing the functional coherence of independently evolved lncRNAs in distantly related mammalian lineages.

A cross-species comparison of escape from X inactivation in Eutheria: implications for evolution of X chromosome inactivation.

Sex chromosome dosage compensation in both eutherian and marsupial mammals is achieved by X chromosome inactivation (XCI)--transcriptional repression that silences one of the two X chromosomes in the somatic cells of females. We recently used RNA fluorescent in situ hybridization (FISH) to show, in individual nuclei, that marsupial X inactivation (in the absence of XIST) occurs on a gene-by-gene basis, and that escape from inactivation is stochastic and independent of gene location. In the absence of similar data from fibroblast cell lines of eutherian representatives, a meaningful comparison is lacking. We therefore used RNA-FISH to examine XCI in fibroblast cell lines obtained from three distantly related eutherian model species: African savannah elephant (Loxodonta africana), mouse (Mus musculus) and human (Homo sapiens). We show that, unlike the orthologous marsupial X, inactivation of the X conserved region (XCR) in eutherians generally is complete. Two-colour RNA-FISH on female human, mouse and elephant interphase nuclei showed that XCR loci have monoallelic expression in almost all nuclei. However, we found that many loci located in the evolutionarily distinct recently added region (XAR) displayed reproducible locus-specific frequencies of nuclei with either one or two active X alleles. We propose that marsupial XCI retains features of an ancient incomplete silencing mechanism that was augmented by the evolution of the XIST gene that progressively stabilized the eutherian XCR. In contrast, the recently added region of the eutherian X displays an incomplete inactivation profile similar to that observed on the evolutionarily distinct marsupial X and the independently evolved monotreme X chromosomes.

Divergent patterns of meiotic double strand breaks and synapsis initiation dynamics suggest an evolutionary shift in the meiosis program between American and Australian marsupials.

In eutherian mammals, hundreds of programmed DNA double-strand breaks (DSBs) are generated at the onset of meiosis. The DNA damage response is then triggered. Although the dynamics of this response is well studied in eutherian mammals, recent findings have revealed different patterns of DNA damage signaling and repair in marsupial mammals. To better characterize these differences, here we analyzed synapsis and the chromosomal distribution of meiotic DSBs markers in three different marsupial species (Thylamys elegans, Dromiciops gliorides, and Macropus eugenii) that represent South American and Australian Orders. Our results revealed inter-specific differences in the chromosomal distribution of DNA damage and repair proteins, which were associated with differing synapsis patterns. In the American species T. elegans and D. gliroides, chromosomal ends were conspicuously polarized in a bouquet configuration and synapsis progressed exclusively from the telomeres towards interstitial regions. This was accompanied by sparse H2AX phosphorylation, mainly accumulating at chromosomal ends. Accordingly, RAD51 and RPA were mainly localized at chromosomal ends throughout prophase I in both American marsupials, likely resulting in reduced recombination rates at interstitial positions. In sharp contrast, synapsis initiated at both interstitial and distal chromosomal regions in the Australian representative M. eugenii, the bouquet polarization was incomplete and ephemeral, γH2AX had a broad nuclear distribution, and RAD51 and RPA foci displayed an even chromosomal distribution. Given the basal evolutionary position of T. elegans, it is likely that the meiotic features reported in this species represent an ancestral pattern in marsupials and that a shift in the meiotic program occurred after the split of D. gliroides and the Australian marsupial clade. Our results open intriguing questions about the regulation and homeostasis of meiotic DSBs in marsupials. The low recombination rates observed at the interstitial chromosomal regions in American marsupials can result in the formation of large linkage groups, thus having an impact in the evolution of their genomes.

The admixed brushtail possum genome reveals invasion history in New Zealand and novel imprinted genes.

Combining genome assembly with population and functional genomics can provide valuable insights to development and evolution, as well as tools for species management. Here, we present a chromosome-level genome assembly of the common brushtail possum (Trichosurus vulpecula), a model marsupial threatened in parts of their native range in Australia, but also a major introduced pest in New Zealand. Functional genomics reveals post-natal activation of chemosensory and metabolic genes, reflecting unique adaptations to altricial birth and delayed weaning, a hallmark of marsupial development. Nuclear and mitochondrial analyses trace New Zealand possums to distinct Australian subspecies, which have subsequently hybridised. This admixture allowed phasing of parental alleles genome-wide, ultimately revealing at least four genes with imprinted, parent-specific expression not yet detected in other species (MLH1, EPM2AIP1, UBP1 and GPX7). We find that reprogramming of possum germline imprints, and the wider epigenome, is similar to eutherian mammals except onset occurs after birth. Together, this work is useful for genetic-based control and conservation of possums, and contributes to understanding of the evolution of novel mammalian epigenetic traits.