Conservation of Glutathione Transferase mRNA and Protein Sequences Similar to Human and Horse Alpha Class GST A3-3 across Dog, Goat, and Opossum Species

The glutathione transferase A3-3 (GST A3-3) homodimeric enzyme is the most efficient enzyme that catalyzes isomerization of the precursors of testosterone, estradiol, and progesterone in the gonads of humans and horses. However, the presence of GST A3-3 orthologs with equally high ketosteroid isomerase activity has not been verified in other mammalian species, even though pig and cattle homologs have been cloned and studied. Identifying GSTA3 genes is a challenge because of multiple GSTA gene duplications (e.g., 12 in the human genome); consequently, the GSTA3 gene is not annotated in most genomes. To improve our understanding of GSTA3 gene products and their functions across diverse mammalian species, we cloned homologs of the horse and human GSTA3 mRNAs from the testes of a dog, goat, and gray short-tailed opossum, the genomes of which all currently lack GSTA3 gene annotations. The resultant novel GSTA3 mRNA and inferred protein sequences had a high level of conservation with human GSTA3 mRNA and protein sequences (≥70% and ≥64% identities, respectively). Sequence conservation was also apparent for the 12 residues of the "H-site" in the 222 amino acid GSTA3 protein that is known to interact with the steroid substrates. Modeling predicted that the dog GSTA3-3 may be a more active ketosteroid isomerase than the corresponding goat or opossum enzymes. However, expression of the GSTA3 gene was higher in liver than in other dog tissue. Our results improve understanding of the active sites of mammalian GST A3-3 enzymes, inhibitors of which might be useful for reducing steroidogenesis for medical purposes, such as fertility control or treatment of steroid-dependent diseases.

Origin and Evolution of Marsupial-specific Imprinting Clusters Through Lineage-specific Gene Duplications and Acquisition of Promoter Differential Methylation

Genomic imprinting is a parent-of-origin-specific expression phenomenon that plays fundamental roles in many biological processes. In animals, imprinting is only observed in therian mammals, with ∼200 imprinted genes known in humans and mice. The imprinting pattern in marsupials has been minimally investigated by examining orthologs to known eutherian imprinted genes. To identify marsupial-specific imprinting in an unbiased way, we performed RNA-seq studies on samples of fetal brain and placenta from the reciprocal cross progeny of two laboratory opossum stocks. We inferred allele-specific expression for >3,000 expressed genes and discovered/validated 13 imprinted genes, including three previously known imprinted genes, Igf2r, Peg10, and H19. We estimate that marsupials imprint ∼60 autosomal genes, which is a much smaller set compared with eutherians. Among the nine novel imprinted genes, three noncoding RNAs have no known homologs in eutherian mammals, while the remaining genes have important functions in pluripotency, transcription regulation, nucleolar homeostasis, and neural differentiation. Methylation analyses at promoter CpG islands revealed differentially methylated regions in five of these marsupial-specific imprinted genes, suggesting that differential methylation is a common mechanism in the epigenetic regulation of marsupial imprinting. Clustering and co-regulation were observed at marsupial imprinting loci Pou5f3-Npdc1 and Nkrfl-Ipncr2, but eutherian-type multi-gene imprinting clusters were not detected. Also differing from eutherian mammals, the brain and placenta imprinting profiles are remarkably similar in opossums, presumably due to the shared origin of these organs from the trophectoderm. Our results contribute to a fuller understanding of the origin, evolution, and mechanisms of genomic imprinting in therian mammals.

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.

The generation of haploid gametes is a hallmark of sexual reproduction. And this is accomplished by a complex, albeit tightly regulated, reductional cell division called meiosis. Although meiosis has been extensively studied in eutherian mammal model species, our understanding of the mechanisms regulating chromosome synapsis, recombination and segregation during meiosis progression is still incomplete especially in non-eutherian mammals. To fill this gap and capture the diversity of meiotic strategies among mammals, we study previously uncharacterised representative marsupial species, an evolutionary assemblage that last shared a common ancestry more than 80 million years ago. We uncover novel, hence non-canonical, strategies for sex chromosome pairing, DNA repair, recombination and transcription. Most importantly, we reveal the uniqueness of marsupial meiosis, which includes the unprecedented detection of alternative mechanism (ALT) for the paternal control of telomere length during prophase I. Our findings suggest that ALT (previously only associated to cancer cells) could play a role in telomere homeostasis in mammalian germ cells.

Evolutionary Characterization of the Short Protein SPAAR

Microproteins (<100 amino acids) are receiving increasing recognition as important participants in numerous biological processes, but their evolutionary dynamics are poorly understood. SPAAR is a recently discovered microprotein that regulates muscle regeneration and angiogenesis through interactions with conserved signaling pathways. Interestingly, SPAAR does not belong to any known protein family and has known homologs exclusively among placental mammals. This lack of distant homology could be caused by challenges in homology detection of short sequences, or it could indicate a recent de novo emergence from a noncoding sequence. By integrating syntenic alignments and homology searches, we identify SPAAR orthologs in marsupials and monotremes, establishing that SPAAR has existed at least since the emergence of mammals. SPAAR shows substantial primary sequence divergence but retains a conserved protein structure. In primates, we infer two independent evolutionary events leading to the de novo origination of 5' elongated isoforms of SPAAR from a noncoding sequence and find evidence of adaptive evolution in this extended region. Thus, SPAAR may be of ancient origin, but it appears to be experiencing continual evolutionary innovation in mammals.

Developmental and evolutionary dynamics of cis-regulatory elements in mouse cerebellar cells

Organ development is orchestrated by cell- and time-specific gene regulatory networks. In this study, we investigated the regulatory basis of mouse cerebellum development from early neurogenesis to adulthood. By acquiring snATAC-seq (single-nucleus assay for transposase accessible chromatin using sequencing) profiles for ~90,000 cells spanning 11 stages, we mapped cerebellar cell types and identified candidate cis-regulatory elements (CREs). We detected extensive spatiotemporal heterogeneity among progenitor cells and a gradual divergence in the regulatory programs of cerebellar neurons during differentiation. Comparisons to vertebrate genomes and snATAC-seq profiles for ∼20,000 cerebellar cells from the marsupial opossum revealed a shared decrease in CRE conservation during development and differentiation as well as differences in constraint between cell types. Our work delineates the developmental and evolutionary dynamics of gene regulation in cerebellar cells and provides insights into mammalian organ development.

A Chromosome-Level Genome of the Agile Gracile Mouse Opossum (Gracilinanus agilis)

There are more than 100 species of American didelphid marsupials (opossums and mouse opossums). Limited genomic resources for didelphids exists, with only two publicly available genome assemblies compared with dozens in the case of their Australasian counterparts. This discrepancy impedes evolutionary and ecological research. To address this gap, we assembled a high-quality chromosome-level genome of the agile gracile mouse opossum (Gracilinanus agilis) using a combination of stLFR sequencing, polishing with mate-pair data, and anchoring onto pseudochromosomes using Hi-C. This species employs a rare life-history strategy, semelparity, and all G. agilis males and most females die at the end of their first breeding season after succumbing to stress and exhaustion. The 3.7-Gb chromosome-level assembly, with 92.6% anchored onto pseudochromosomes, has a scaffold N50 of 683.5 Mb and a contig N50 of 56.9 kb. The genome assembly shows high completeness, with a mammalian BUSCO score of 88.1%. Around 49.7% of the genome contains repetitive elements. Gene annotation yielded 24,425 genes, of which 83.9% were functionally annotated. The G. agilis genome is an important resource for future studies of marsupial biology, evolution, and conservation.

Molecular Evolution of CatSper in Mammals and Function of Sperm Hyperactivation in Gray Short-Tailed Opossum

Males have evolved species-specifical sperm morphology and swimming patterns to adapt to different fertilization environments. In eutherians, only a small fraction of the sperm overcome the diverse obstacles in the female reproductive tract and successfully migrate to the fertilizing site. Sperm arriving at the fertilizing site show hyperactivated motility, a unique motility pattern displaying asymmetric beating of sperm flagella with increased amplitude. This motility change is triggered by Ca²⁺ influx through the sperm-specific ion channel, CatSper. However, the current understanding of the CatSper function and its molecular regulation is limited in eutherians. Here, we report molecular evolution and conservation of the CatSper channel in the genome throughout eutherians and marsupials. Sequence analyses reveal that CatSper proteins are slowly evolved in marsupials. Using an American marsupial, gray short-tailed opossum (Monodelphis domestica), we demonstrate the expression of CatSper in testes and its function in hyperactivation and unpairing of sperm. We demonstrate that a conserved IQ-like motif in CatSperζ is required for CatSperζ interaction with the pH-tuned Ca²⁺ sensor, EFCAB9, for regulating CatSper activity. Recombinant opossum EFCAB9 can interact with mouse CatSperζ despite high sequence divergence of CatSperζ among CatSper subunits in therians. Our finding suggests that molecular characteristics and functions of CatSper are evolutionarily conserved in gray short-tailed opossum, unraveling the significance of sperm hyperactivation and fertilization in marsupials for the first time.

Didelphis albiventris: an overview of unprecedented transcriptome sequencing of the white-eared opossum

BACKGROUND

The white-eared opossum (Didelphis albiventris) is widely distributed throughout Brazil and South America. It has been used as an animal model for studying different scientific questions ranging from the restoration of degraded green areas to medical aspects of Chagas disease, leishmaniasis and resistance against snake venom. As a marsupial, D. albiventris can also contribute to the understanding of the molecular mechanisms that govern the different stages of organogenesis. Opossum joeys are born after only 13 days, and the final stages of organogenesis occur when the neonates are inside the pouch, depending on lactation. As neither the genome of this opossum species nor its transcriptome has been completely sequenced, the use of D. albiventris as an animal model is limited. In this work, we sequenced the D. albiventris transcriptome by RNA-seq to obtain the first catalogue of differentially expressed (DE) genes and gene ontology (GO) annotations during the neonatal stages of marsupial development.

RESULTS

The D. albiventris transcriptome was obtained from whole neonates harvested at birth (P0), at 5 days of age (P5) and at 10 days of age (P10). The de novo assembly of these transcripts generated 85,338 transcripts. Approximately 30% of these transcripts could be mapped against the amino acid sequences of M. domestica, the evolutionarily closest relative of D. albiventris to be sequenced thus far. Among the expressed transcripts, 2077 were found to be DE between P0 and P5, 13,780 between P0 and P10, and 1453 between P5 and P10. The enriched GO terms were mainly related to the immune system, blood tissue development and differentiation, vision, hearing, digestion, the CNS and limb development.

CONCLUSIONS

The elucidation of opossum transcriptomes provides an out-group for better understanding the distinct characteristics associated with the evolution of mammalian species. This study provides the first transcriptome sequences and catalogue of genes for a marsupial species at different neonatal stages, allowing the study of the mechanisms involved in organogenesis.

The SCJ Small Parsimony Problem for Weighted Gene Adjacencies

Reconstructing ancestral gene orders in a given phylogeny is a classical problem in comparative genomics. Most existing methods compare conserved features in extant genomes in the phylogeny to define potential ancestral gene adjacencies, and either try to reconstruct all ancestral genomes under a global evolutionary parsimony criterion, or, focusing on a single ancestral genome, use a scaffolding approach to select a subset of ancestral gene adjacencies, generally aiming at reducing the fragmentation of the reconstructed ancestral genome. In this paper, we describe an exact algorithm for the Small Parsimony Problem that combines both approaches. We consider that gene adjacencies at internal nodes of the species phylogeny are weighted, and we introduce an objective function defined as a convex combination of these weights and the evolutionary cost under the Single-Cut-or-Join (SCJ) model. The weights of ancestral gene adjacencies can, e.g., be obtained through the recent availability of ancient DNA sequencing data, which provide a direct hint at the genome structure of the considered ancestor, or through probabilistic analysis of gene adjacencies evolution. We show the NP-hardness of our problem variant and propose a Fixed-Parameter Tractable algorithm based on the Sankoff-Rousseau dynamic programming algorithm that also allows to sample co-optimal solutions. We apply our approach to mammalian and bacterial data providing different degrees of complexity. We show that including adjacency weights in the objective has a significant impact in reducing the fragmentation of the reconstructed ancestral gene orders. An implementation is available at http://github.com/nluhmann/PhySca.

Molecular evolution of DNMT1 in vertebrates: Duplications in marsupials followed by positive selection

DNA methylation is mediated by a conserved family of DNA methyltransferases (Dnmts). The human genome encodes three active Dnmts (Dnmt1, Dnmt3a and Dnmt3b), the tRNA methyltransferase Dnmt2, and the regulatory protein Dnmt3L. Despite their high degree of conservation among different species, genes encoding Dnmts have been duplicated and/or lost in multiple lineages throughout evolution, indicating that the DNA methylation machinery has some potential to undergo evolutionary change. However, little is known about the extent to which this machinery, or the methylome, varies among vertebrates. Here, we study the molecular evolution of Dnmt1, the enzyme responsible for maintenance of DNA methylation patterns after replication, in 79 vertebrate species. Our analyses show that all studied species exhibit a single copy of the DNMT1 gene, with the exception of tilapia and marsupials (tammar wallaby, koala, Tasmanian devil and opossum), each of which displays two apparently functional DNMT1 copies. Our phylogenetic analyses indicate that DNMT1 duplicated before the radiation of major marsupial groups (i.e., at least ~75 million years ago), thus giving rise to two DNMT1 copies in marsupials (copy 1 and copy 2). In the opossum lineage, copy 2 was lost, and copy 1 recently duplicated again, generating three DNMT1 copies: two putatively functional genes (copy 1a and 1b) and one pseudogene (copy 1ψ). Both marsupial copies (DNMT1 copies 1 and 2) are under purifying selection, and copy 2 exhibits elevated rates of evolution and signatures of positive selection, suggesting a scenario of neofunctionalization. This gene duplication might have resulted in modifications in marsupial methylomes and their dynamics.

Superior ab initio identification, annotation and characterisation of TEs and segmental duplications from genome assemblies

Transposable Elements (TEs) are mobile DNA sequences that make up significant fractions of amniote genomes. However, they are difficult to detect and annotate ab initio because of their variable features, lengths and clade-specific variants. We have addressed this problem by refining and developing a Comprehensive ab initio Repeat Pipeline (CARP) to identify and cluster TEs and other repetitive sequences in genome assemblies. The pipeline begins with a pairwise alignment using krishna, a custom aligner. Single linkage clustering is then carried out to produce families of repetitive elements. Consensus sequences are then filtered for protein coding genes and then annotated using Repbase and a custom library of retrovirus and reverse transcriptase sequences. This process yields three types of family: fully annotated, partially annotated and unannotated. Fully annotated families reflect recently diverged/young known TEs present in Repbase. The remaining two types of families contain a mixture of novel TEs and segmental duplications. These can be resolved by aligning these consensus sequences back to the genome to assess copy number vs. length distribution. Our pipeline has three significant advantages compared to other methods for ab initio repeat identification: 1) we generate not only consensus sequences, but keep the genomic intervals for the original aligned sequences, allowing straightforward analysis of evolutionary dynamics, 2) consensus sequences represent low-divergence, recently/currently active TE families, 3) segmental duplications are annotated as a useful by-product. We have compared our ab initio repeat annotations for 7 genome assemblies to other methods and demonstrate that CARP compares favourably with RepeatModeler, the most widely used repeat annotation package.

The transcriptome of the Didelphis virginiana opossum kidney OK proximal tubule cell line

The OK cell line derived from the kidney of a female opossum Didelphis virginiana has proven to be a useful model in which to investigate the unique regulation of ion transport and membrane trafficking mechanisms in the proximal tubule (PT). Sequence data and comparison of the transcriptome of this cell line to eutherian mammal PTs would further broaden the utility of this culture model. However, the genomic sequence for D. virginiana is not available and although a draft genome sequence for the opossum Monodelphis domestica (sequenced in 2012 by the Broad Institute) exists, transcripts sequenced from both species show significant divergence. The M. domestica sequence is not highly annotated, and the majority of transcripts are predicted rather than experimentally validated. Using deep RNA sequencing of the D. virginiana OK cell line, we characterized its transcriptome via de novo transcriptome assembly and alignment to the M. domestica genome. The quality of the de novo assembled transcriptome was assessed by the extent of homology to sequences in nucleotide and protein databases. Gene expression levels in the OK cell line, from both the de novo transcriptome and genes aligned to the M. domestica genome, were compared with publicly available rat kidney nephron segment expression data. Our studies demonstrate the expression in OK cells of numerous PT-specific ion transporters and other key proteins relevant for rodent and human PT function. Additionally, the sequence and expression data reported here provide an important resource for genetic manipulation and other studies on PT cell function using these cells.

Transcriptomic Changes Associated with Pregnancy in a Marsupial, the Gray Short-Tailed Opossum Monodelphis domestica

Live birth has emerged as a reproductive strategy many times across vertebrate evolution; however, mammals account for the majority of viviparous vertebrates. Marsupials are a mammalian lineage that last shared a common ancestor with eutherians (placental mammals) over 148 million years ago. Marsupials are noted for giving birth to highly altricial young after a short gestation, and represent humans’ most distant viviparous mammalian relatives. Here we ask what insight can be gained into the evolution of viviparity in mammals specifically and vertebrates in general by analyzing the global uterine transcriptome in a marsupial. Transcriptome analyses were performed using NextGen sequencing of uterine RNA samples from the gray short-tailed opossum, Monodelphis domestica. Samples were collected from late stage pregnant, virgin, and non-pregnant experienced breeders. Three different algorithms were used to determine differential expression, and results were confirmed by quantitative PCR. Over 900 opossum gene transcripts were found to be significantly more abundant in the pregnant uterus than non-pregnant, and over 1400 less so. Most with increased abundance were genes related to metabolism, immune systems processes, and transport. This is the first study to characterize the transcriptomic differences between pregnant, non-pregnant breeders, and virgin marsupial uteruses and helps to establish a set of pregnancy-associated genes in the opossum. These observations allowed for comparative analyses of the differentially transcribed genes with other mammalian and non-mammalian viviparous species, revealing similarities in pregnancy related gene expression over 300 million years of amniote evolution.

Euarchontan Opsin Variation Brings New Focus to Primate Origins

Debate on the adaptive origins of primates has long focused on the functional ecology of the primate visual system. For example, it is hypothesized that variable expression of short- (SWS1) and middle-to-long-wavelength sensitive (M/LWS) opsins, which confer color vision, can be used to infer ancestral activity patterns and therefore selective ecological pressures. A problem with this approach is that opsin gene variation is incompletely known in the grandorder Euarchonta, that is, the orders Scandentia (treeshrews), Dermoptera (colugos), and Primates. The ancestral state of primate color vision is therefore uncertain. Here, we report on the genes (OPN1SW and OPN1LW) that encode SWS1 and M/LWS opsins in seven species of treeshrew, including the sole nocturnal scandentian Ptilocercus lowii. In addition, we examined the opsin genes of the Central American woolly opossum (Caluromys derbianus), an enduring ecological analogue in the debate on primate origins. Our results indicate: 1) retention of ultraviolet (UV) visual sensitivity in C. derbianus and a shift from UV to blue spectral sensitivities at the base of Euarchonta; 2) ancient pseudogenization of OPN1SW in the ancestors of P. lowii, but a signature of purifying selection in those of C. derbianus; and, 3) the absence of OPN1LW polymorphism among diurnal treeshrews. These findings suggest functional variation in the color vision of nocturnal mammals and a distinctive visual ecology of early primates, perhaps one that demanded greater spatial resolution under light levels that could support cone-mediated color discrimination.

The Relationship between Gene Network Structure and Expression Variation among Individuals and Species

Variation among individuals is a prerequisite of evolution by natural selection. As such, identifying the origins of variation is a fundamental goal of biology. We investigated the link between gene interactions and variation in gene expression among individuals and species using the mammalian limb as a model system. We first built interaction networks for key genes regulating early (outgrowth; E9.5-11) and late (expansion and elongation; E11-13) limb development in mouse. This resulted in an Early (ESN) and Late (LSN) Stage Network. Computational perturbations of these networks suggest that the ESN is more robust. We then quantified levels of the same key genes among mouse individuals and found that they vary less at earlier limb stages and that variation in gene expression is heritable. Finally, we quantified variation in gene expression levels among four mammals with divergent limbs (bat, opossum, mouse and pig) and found that levels vary less among species at earlier limb stages. We also found that variation in gene expression levels among individuals and species are correlated for earlier and later limb development. In conclusion, results are consistent with the robustness of the ESN buffering among-individual variation in gene expression levels early in mammalian limb development, and constraining the evolution of early limb development among mammalian species.

PhyLAT: a phylogenetic local alignment tool

MOTIVATION

The expansion of DNA sequencing capacity has enabled the sequencing of whole genomes from a number of related species. These genomes can be combined in a multiple alignment that provides useful information about the evolutionary history at each genomic locus. One area in which evolutionary information can productively be exploited is in aligning a new sequence to a database of existing, aligned genomes. However, existing high-throughput alignment tools are not designed to work effectively with multiple genome alignments.

RESULTS

We introduce PhyLAT, the phylogenetic local alignment tool, to compute local alignments of a query sequence against a fixed multiple-genome alignment of closely related species. PhyLAT uses a known phylogenetic tree on the species in the multiple alignment to improve the quality of its computed alignments while also estimating the placement of the query on this tree. It combines a probabilistic approach to alignment with seeding and expansion heuristics to accelerate discovery of significant alignments. We provide evidence, using alignments of human chromosome 22 against a five-species alignment from the UCSC Genome Browser database, that PhyLAT's alignments are more accurate than those of other commonly used programs, including BLAST, POY, MAFFT, MUSCLE and CLUSTAL. PhyLAT also identifies more alignments in coding DNA than does pairwise alignment alone. Finally, our tool determines the evolutionary relationship of query sequences to the database more accurately than do POY, RAxML, EPA or pplacer.

A first-generation integrated tammar wallaby map and its use in creating a tammar wallaby first-generation virtual genome map

BACKGROUND

The limited (2X) coverage of the tammar wallaby (Macropus eugenii) genome sequence dataset currently presents a challenge for assembly and anchoring onto chromosomes. To provide a framework for this assembly, it would be a great advantage to have a dense map of the tammar wallaby genome. However, only limited mapping data are available for this non-model species, comprising a physical map and a linkage map.

RESULTS

We combined all available tammar wallaby mapping data to create a tammar wallaby integrated map, using the Location DataBase (LDB) strategy. This first-generation integrated map combines all available information from the second-generation tammar wallaby linkage map with 148 loci, and extensive FISH mapping data for 492 loci, especially for genes likely to be located at the ends of wallaby chromosomes or at evolutionary breakpoints inferred from comparative information. For loci whose positions are only approximately known, their location in the integrated map was refined on the basis of comparative information from opossum (Monodelphis domestica) and human. Interpolation of segments from the opossum and human assemblies into the integrated map enabled the subsequent construction of a tammar wallaby first-generation virtual genome map, which comprises 14336 markers, including 13783 genes recruited from opossum and human assemblies. Both maps are freely available at http://compldb.angis.org.au.

CONCLUSIONS

The first-generation integrated map and the first-generation virtual genome map provide a backbone for the chromosome assembly of the tammar wallaby genome sequence. For example, 78% of the 10257 gene-scaffolds in the Ensembl annotation of the tammar wallaby genome sequence (including 10522 protein-coding genes) can now be given a chromosome location in the tammar wallaby virtual genome map.

Adaptive evolution of the venom-targeted vWF protein in opossums that eat pitvipers

The rapid evolution of venom toxin genes is often explained as the result of a biochemical arms race between venomous animals and their prey. However, it is not clear that an arms race analogy is appropriate in this context because there is no published evidence for rapid evolution in genes that might confer toxin resistance among routinely envenomed species. Here we report such evidence from an unusual predator-prey relationship between opossums (Marsupialia: Didelphidae) and pitvipers (Serpentes: Crotalinae). In particular, we found high ratios of replacement to silent substitutions in the gene encoding von Willebrand Factor (vWF), a venom-targeted hemostatic blood protein, in a clade of opossums known to eat pitvipers and to be resistant to their hemorrhagic venom. Observed amino-acid substitutions in venom-resistant opossums include changes in net charge and hydrophobicity that are hypothesized to weaken the bond between vWF and one of its toxic snake-venom ligands, the C-type lectin-like protein botrocetin. Our results provide the first example of rapid adaptive evolution in any venom-targeted molecule, and they support the notion that an evolutionary arms race might be driving the rapid evolution of snake venoms. However, in the arms race implied by our results, venomous snakes are prey, and their venom has a correspondingly defensive function in addition to its usual trophic role.

Five-vertebrate ChIP-seq reveals the evolutionary dynamics of transcription factor binding

Transcription factors (TFs) direct gene expression by binding to DNA regulatory regions. To explore the evolution of gene regulation, we used chromatin immunoprecipitation with high-throughput sequencing (ChIP-seq) to determine experimentally the genome-wide occupancy of two TFs, CCAAT/enhancer-binding protein alpha and hepatocyte nuclear factor 4 alpha, in the livers of five vertebrates. Although each TF displays highly conserved DNA binding preferences, most binding is species-specific, and aligned binding events present in all five species are rare. Regions near genes with expression levels that are dependent on a TF are often bound by the TF in multiple species yet show no enhanced DNA sequence constraint. Binding divergence between species can be largely explained by sequence changes to the bound motifs. Among the binding events lost in one lineage, only half are recovered by another binding event within 10 kilobases. Our results reveal large interspecies differences in transcriptional regulation and provide insight into regulatory evolution.

Comparative genomic analysis of mammalian NKG2D ligand family genes provides insights into their origin and evolution

NKG2D is a major activating receptor of natural killer cells. Its ligands are major histocompatibility complex (MHC) class I-like molecules whose expression is induced by cellular stresses such as infections and tumorigenesis. Humans have two families of NKG2D ligands (NKG2DL): MHC class I-related chains (MIC) encoded in the MHC and UL16-binding proteins (ULBP) encoded outside the MHC. By contrast, mice have only the latter family of ligands; instead, they have non-MHC-encoded MILL molecules that are closely related to MIC, but do not function as NKG2DL. To gain insights into the origin and evolution of MIC, ULBP, and MILL gene families, we conducted comparative genomic analysis of NKG2DL family genes in five mammalian species. In the opossum MHC, we identified a ULBP-like gene adjacent to a previously described MIC-like gene, suggesting that ULBP genes were originally encoded in the MHC. The opossum genome also contained a transcribed MILL-like gene in a region syntenic to the rodent regions encoding MILL molecules. These observations indicate that MIC-, ULBP-, and MILL-like genes emerged before the divergence of placental and marsupial mammals. Comparison of the human, cattle, rat, mouse, and opossum genomes indicates that after emigration from the MHC, ULBP genes underwent extensive duplications in each species. In mice, some of the ULBP genes appear to have been translocated telomerically on the same chromosome, forming a major cluster of existent NKG2DL genes.

Long-range regulation is a major driving force in maintaining genome integrity

BACKGROUND

The availability of newly sequenced vertebrate genomes, along with more efficient and accurate alignment algorithms, have enabled the expansion of the field of comparative genomics. Large-scale genome rearrangement events modify the order of genes and non-coding conserved regions on chromosomes. While certain large genomic regions have remained intact over much of vertebrate evolution, others appear to be hotspots for genomic breakpoints. The cause of the non-uniformity of breakpoints that occurred during vertebrate evolution is poorly understood.

RESULTS

We describe a machine learning method to distinguish genomic regions where breakpoints would be expected to have deleterious effects (called breakpoint-refractory regions) from those where they are expected to be neutral (called breakpoint-susceptible regions). Our predictor is trained using breakpoints that took place along the human lineage since amniote divergence. Based on our predictions, refractory and susceptible regions have very distinctive features. Refractory regions are significantly enriched for conserved non-coding elements as well as for genes involved in development, whereas susceptible regions are enriched for housekeeping genes, likely to have simpler transcriptional regulation.

CONCLUSION

We postulate that long-range transcriptional regulation strongly influences chromosome break fixation. In many regions, the fitness cost of altering the spatial association between long-range regulatory regions and their target genes may be so high that rearrangements are not allowed. Consequently, only a limited, identifiable fraction of the genome is susceptible to genome rearrangements.

Evolutionary conservation of orthoretroviral long terminal repeats (LTRs) and ab initio detection of single LTRs in genomic data

Background

Retroviral LTRs, paired or single, influence the transcription of both retroviral and non-retroviral genomic sequences. Vertebrate genomes contain many thousand endogenous retroviruses (ERVs) and their LTRs. Single LTRs are difficult to detect from genomic sequences without recourse to repetitiveness or presence in a proviral structure. Understanding of LTR structure increases understanding of LTR function, and of functional genomics. Here we develop models of orthoretroviral LTRs useful for detection in genomes and for structural analysis.

Principal Findings

Although mutated, ERV LTRs are more numerous and diverse than exogenous retroviral (XRV) LTRs. Hidden Markov models (HMMs), and alignments based on them, were created for HML- (human MMTV-like), general-beta-, gamma- and lentiretroviruslike LTRs, plus a general-vertebrate LTR model. Training sets were XRV LTRs and RepBase LTR consensuses. The HML HMM was most sensitive and detected 87% of the HML LTRs in human chromosome 19 at 96% specificity. By combining all HMMs with a low cutoff, for screening, 71% of all LTRs found by RepeatMasker in chromosome 19 were found. HMM consensus sequences had a conserved modular LTR structure. Target site duplications (TG-CA), TATA (occasionally absent), an AATAAA box and a T-rich region were prominent features. Most of the conservation was located in, or adjacent to, R and U5, with evidence for stem loops. Several of the long HML LTRs contained long ORFs inserted after the second A rich module. HMM consensus alignment allowed comparison of functional features like transcriptional start sites (sense and antisense) between XRVs and ERVs.

Conclusion

The modular conserved and redundant orthoretroviral LTR structure with three A-rich regions is reminiscent of structurally relaxed Giardia promoters. The five HMMs provided a novel broad range, repeat-independent, ab initio LTR detection, with prospects for greater generalisation, and insight into LTR structure, which may aid development of LTR-targeted pharmaceuticals.

Exon 3 of the growth hormone receptor (GH-R) is specific to eutherian mammals

Growth hormone receptor (GH-R) plays a critical role in the control of growth and metabolism in all vertebrates. GH-R consists of 9 coding exons (2-10) in all eutherian mammals, while the chicken only has 8 coding exons, and does not have an orthologous region to exon 3 in eutherians. To further understand the evolutionary origins of exon 3 of the GH-R in eutherians we cloned the full-length GH-R sequence in a marsupial, the tammar wallaby to determine whether exon 3 was present or absent in marsupial liver cDNA. There was no evidence for the presence of an exon 3 containing mRNA in sequence of tammar pouch young and adult livers. We next examined the genomes of the platypus (a monotreme mammal) and the grey short-tailed opossum (another marsupial). Like the tammar, the GH-R gene of neither species contained an exon 3. GH receptor can obviously function in the absence of this exon, raising speculation about the function of this domain, if any, in eutherians. A comparison of exon 3 protein sequences within 16 species of eutherian mammals showed that there was approximately 75% homology in the domain but only 3 of the 21 amino acids were identical (Leu12, Gln13 and Pro17). Interestingly, we detected greater evolutionary divergence in exon 3 sequences from species that have variants of GH or prolactin (PRL) in their placentas. These data show that exon 3 was inserted into the GH-R after the divergence of the marsupial and eutherian lineages at least 130 million years ago.

Pattern of the divergence of olfactory receptor genes during tetrapod evolution

The olfactory receptor (OR) multigene family is responsible for the sense of smell in vertebrate species. OR genes are scattered widely in our chromosomes and constitute one of the largest gene families in eutherian genomes. Some previous studies revealed that eutherian OR genes diverged mainly during early mammalian evolution. However, the exact period when, and the ecological reason why eutherian ORs strongly diverged has remained unclear. In this study, I performed a strict data mining effort for marsupial opossum OR sequences and bootstrap analyses to estimate the periods of chromosomal migrations and gene duplications of OR genes during tetrapod evolution. The results indicate that chromosomal migrations occurred mainly during early vertebrate evolution before the monotreme-placental split, and that gene duplications occurred mainly during early mammalian evolution between the bird-mammal split and marsupial-placental split, coinciding with the reduction of opsin genes in primitive mammals. It could be thought that the previous chromosomal dispersal allowed the OR genes to subsequently expand easily, and the nocturnal adaptation of early mammals might have triggered the OR gene expansion.

Mammalian genome projects reveal new growth hormone (GH) sequences. Characterization of the GH-encoding genes of armadillo (Dasypus novemcinctus), hedgehog (Erinaceus europaeus), bat (Myotis lucifugus), hyrax (Procavia capensis), shrew (Sorex araneus), ground squirrel (Spermophilus tridecemlineatus), elephant (Loxodonta africana), cat (Felis catus) and opossum (Monodelphis domestica)

Mammalian growth hormone (GH) sequences have been shown previously to display episodic evolution: the sequence is generally strongly conserved but on at least two occasions during mammalian evolution (on lineages leading to higher primates and ruminants) bursts of rapid evolution occurred. However, the number of mammalian orders studied previously has been relatively limited, and the availability of sequence data via mammalian genome projects provides the potential for extending the range of GH gene sequences examined. Complete or nearly complete GH gene sequences for six mammalian species for which no data were previously available have been extracted from the genome databases-Dasypus novemcinctus (nine-banded armadillo), Erinaceus europaeus (western European hedgehog), Myotis lucifugus (little brown bat), Procavia capensis (cape rock hyrax), Sorex araneus (European shrew), Spermophilus tridecemlineatus (13-lined ground squirrel). In addition incomplete data for several other species have been extended. Examination of the data in detail and comparison with previously available sequences has allowed assessment of the reliability of deduced sequences. Several of the new sequences differ substantially from the consensus sequence previously determined for eutherian GHs, indicating greater variability than previously recognised, and confirming the episodic pattern of evolution. The episodic pattern is not seen for signal sequences, 5' upstream sequence or synonymous substitutions-it is specific to the mature protein sequence, suggesting that it relates to the hormonal function. The substitutions accumulated during the course of GH evolution have occurred mainly on the side of the hormone facing away from the receptor, in a non-random fashion, and it is suggested that this may reflect interaction of the receptor-bound hormone with other proteins or small ligands.

Ancient fossil specimens of extinct species are genetically more distant to an outgroup than extant sister species are

There exists a remarkable correlation between genetic distance as measured by protein or DNA dissimilarity and time of species divergence as inferred from fossil records. This observation has provoked the molecular clock hypothesis. However, data inconsistent with the hypothesis have steadily accumulated in recent years from studies of extant organisms. Here the published DNA and protein sequences from ancient fossil specimens were examined to see if they would support the molecular clock hypothesis. The hypothesis predicts that ancient specimens cannot be genetically more distant to an outgroup than extant sister species are. Also, two distinct ancient specimens cannot be genetically more distant than their extant sister species are. The findings here do not conform to these predictions. Neanderthals are more distant to chimpanzees and gorillas than modern humans are. Dinosaurs are more distant to frogs than extant birds are. Mastodons are more distant to opossums than other placental mammals are. The genetic distance between dinosaurs and mastodons is greater than that between extant birds and mammals. Therefore, while the molecular clock hypothesis is consistent with some data from extant organisms, it has yet to find support from ancient fossils. Far more damaging to the hypothesis than data from extant organisms, which merely question the constancy of mutation rate, the study of ancient fossil organisms here challenges for the first time the fundamental premise of modern evolution theory that genetic distances had always increased with time in the past history of life on Earth.

Ancient exaptation of a CORE-SINE retroposon into a highly conserved mammalian neuronal enhancer of the proopiomelanocortin gene

The proopiomelanocortin gene (POMC) is expressed in the pituitary gland and the ventral hypothalamus of all jawed vertebrates, producing several bioactive peptides that function as peripheral hormones or central neuropeptides, respectively. We have recently determined that mouse and human POMC expression in the hypothalamus is conferred by the action of two 5′ distal and unrelated enhancers, nPE1 and nPE2. To investigate the evolutionary origin of the neuronal enhancer nPE2, we searched available vertebrate genome databases and determined that nPE2 is a highly conserved element in placentals, marsupials, and monotremes, whereas it is absent in nonmammalian vertebrates. Following an in silico paleogenomic strategy based on genome-wide searches for paralog sequences, we discovered that opossum and wallaby nPE2 sequences are highly similar to members of the superfamily of CORE-short interspersed nucleotide element (SINE) retroposons, in particular to MAR1 retroposons that are widely present in marsupial genomes. Thus, the neuronal enhancer nPE2 originated from the exaptation of a CORE-SINE retroposon in the lineage leading to mammals and remained under purifying selection in all mammalian orders for the last 170 million years. Expression studies performed in transgenic mice showed that two nonadjacent nPE2 subregions are essential to drive reporter gene expression into POMC hypothalamic neurons, providing the first functional example of an exapted enhancer derived from an ancient CORE-SINE retroposon. In addition, we found that this CORE-SINE family of retroposons is likely to still be active in American and Australian marsupial genomes and that several highly conserved exonic, intronic and intergenic sequences in the human genome originated from the exaptation of CORE-SINE retroposons. Together, our results provide clear evidence of the functional novelties that transposed elements contributed to their host genomes throughout evolution.

One of the most striking observations derived from the genomic era is the overwhelming contribution of transposed elements to mammalian genomes. For example, 45% of the human genome is derived from mobile element fragments. Although historically viewed as “junk DNA,” transposed elements could also contribute to novel advantageous functional elements in their host genomes, a process called exaptation. Functionally proven examples of exaptation derived from ancient retroposition events are rare. Using an in silico paleogenomic strategy, we unraveled the evolutionary origin of nPE2, a neuronal enhancer of the proopiomelancortin gene that participates in the production of hypothalamic peptides involved in feeding behavior and stress-induced analgesia. We demonstrate that nPE2 originated from the exaptation of a SINE retroposon in the lineage leading to mammals and remained under purifying selection for the last 170 million years. The difficulty in detecting nPE2 origin as an exapted retroposon illustrates the underestimation of this phenomenon and encourages the finding of the many thousands of retroposon-derived functional elements still hidden within the genomes. Their discovery will contribute to a better understanding of the dynamics of gene evolution and, at a larger scale, the origin of macroevolutionary novelties that lead to the appearance of new species, orders, or classes.

SINEs, evolution and genome structure in the opossum

Short INterspersed Elements (SINEs) are non-autonomous retrotransposons, usually between 100 and 500 base pairs (bp) in length, which are ubiquitous components of eukaryotic genomes. Their activity, distribution, and evolution can be highly informative on genomic structure and evolutionary processes. To determine recent activity, we amplified more than one hundred SINE1 loci in a panel of 43 M. domestica individuals derived from five diverse geographic locations. The SINE1 family has expanded recently enough that many loci were polymorphic, and the SINE1 insertion-based genetic distances among populations reflected geographic distance. Genome-wide comparisons of SINE1 densities and GC content revealed that high SINE1 density is associated with high GC content in a few long and many short spans. Young SINE1s, whether fixed or polymorphic, showed an unbiased GC content preference for insertion, indicating that the GC preference accumulates over long time periods, possibly in periodic bursts. SINE1 evolution is thus broadly similar to human Alu evolution, although it has an independent origin. High GC content adjacent to SINE1s is strongly correlated with bias towards higher AT to GC substitutions and lower GC to AT substitutions. This is consistent with biased gene conversion, and also indicates that like chickens, but unlike eutherian mammals, GC content heterogeneity (isochore structure) is reinforced by substitution processes in the M. domestica genome. Nevertheless, both high and low GC content regions are apparently headed towards lower GC content equilibria, possibly due to a relative shift to lower recombination rates in the recent Monodelphis ancestral lineage. Like eutherians, metatherian (marsupial) mammals have evolved high CpG substitution rates, but this is apparently a convergence in process rather than a shared ancestral state.

Sex chromosome silencing in the marsupial male germ line

In marsupials, dosage compensation involves silencing of the father's X-chromosome. Because no XIST orthologue has been found, how imprinted X-inactivation occurs is unknown. In eutherians, the X is subject to meiotic sex chromosome inactivation (MSCI) in the paternal germ line and persists thereafter as postmeiotic sex chromatin (PMSC). One hypothesis proposes that the paternal X is inherited by the eutherian zygote as a preinactive X and raises the possibility of a similar process in the marsupial germ line. Here we demonstrate that MSCI and PMSC occur in the opossum. Surprisingly, silencing occurs before X-Y association. After MSCI, the X and Y fuse through a dense plate without obvious synapsis. Significantly, sex chromosome silencing continues after meiosis, with the opossum PMSC sharing features of eutherian PMSC. These results reveal a common gametogenic program in two diverse clades of mammals and support the idea that male germ-line silencing may have provided an ancestral form of mammalian dosage compensation.

Genome of the marsupial Monodelphis domestica reveals innovation in non-coding sequences

We report a high-quality draft of the genome sequence of the grey, short-tailed opossum (Monodelphis domestica). As the first metatherian ('marsupial') species to be sequenced, the opossum provides a unique perspective on the organization and evolution of mammalian genomes. Distinctive features of the opossum chromosomes provide support for recent theories about genome evolution and function, including a strong influence of biased gene conversion on nucleotide sequence composition, and a relationship between chromosomal characteristics and X chromosome inactivation. Comparison of opossum and eutherian genomes also reveals a sharp difference in evolutionary innovation between protein-coding and non-coding functional elements. True innovation in protein-coding genes seems to be relatively rare, with lineage-specific differences being largely due to diversification and rapid turnover in gene families involved in environmental interactions. In contrast, about 20% of eutherian conserved non-coding elements (CNEs) are recent inventions that postdate the divergence of Eutheria and Metatheria. A substantial proportion of these eutherian-specific CNEs arose from sequence inserted by transposable elements, pointing to transposons as a major creative force in the evolution of mammalian gene regulation.

V2R gene families degenerated in primates, dog and cow, but expanded in opossum

The V2R genes are expressed in the mammalian vomeronasal organ, and their products are involved in detecting pheromones. Here, we describe the evolution of the V2R gene family. We have found that the human, chimpanzee, macaque, cow and dog V2R gene families have completely degenerated. Each now contains 9-20 pseudogenes but no intact V2R genes. By contrast, opossum has approximately 90 intact V2R genes that mostly arose by duplication after opossum and rodent lineages diverged. One V2R gene subfamily with unusual biology evolved atypically, showing limited expansion in rodents and persistence of a single, albeit sometimes dysfunctional, ortholog in all other species examined.

Base-compositional heterogeneity in the RAG1 locus among didelphid marsupials: implications for phylogenetic inference and the evolution of GC content

Although theoretical studies have suggested that base-compositional heterogeneity can adversely affect phylogenetic reconstruction, only a few empirical examples of this phenomenon, mostly among ancient lineages (with divergence dates > 100 Mya), have been reported. In the course of our phylogenetic research on the New World marsupial family Didelphidae, we sequenced 2790 bp of the RAG1 exon from exemplar species of most extant genera. Phylogenetic analysis of these sequences recovered an anomalous node consisting of two clades previously shown to be distantly related based on analyses of other molecular data. These two clades show significantly increased GC content at RAG1 third codon positions, and the resulting convergence in base composition is strong enough to overwhelm phylogenetic signal from other genes (and morphology) in most analyses of concatenated datasets. This base-compositional convergence occurred relatively recently (over tens rather than hundreds of millions of years), and the affected gene region is still in a state of evolutionary disequilibrium. Both mutation rate and substitution rate are higher in GC-rich didelphid taxa, observations consistent with RAG1 sequences having experienced a higher rate of recombination in the convergent lineages.

Phylogenetic appearance of neuropeptide S precursor proteins in tetrapods

Sleep and emotional behavior are two hallmarks of vertebrate animal behavior, implying that specialized neuronal circuits and dedicated neurochemical messengers may have been developed during evolution to regulate such complex behaviors. Neuropeptide S (NPS) is a newly identified peptide transmitter that activates a typical G protein-coupled receptor. Central administration of NPS produces profound arousal, enhances wakefulness and suppresses all stages of sleep. In addition, NPS can alleviate behavioral responses to stress by producing anxiolytic-like effects. A bioinformatic analysis of current genome databases revealed that the NPS peptide precursor gene is present in all vertebrates with the exception of fish. A high level of sequence conservation, especially of aminoterminal structures was detected, indicating stringent requirements for agonist-induced receptor activation. Duplication of the NPS precursor gene was only found in one out of two marsupial species with sufficient genome coverage (Monodelphis domestica; opossum), indicating that the duplicated opossum NPS sequence might have arisen as an isolated event. Pharmacological analysis of both Monodelphis NPS peptides revealed that only the closely related NPS peptide retained agonistic activity at NPS receptors. The duplicated precursor might be either a pseudogene or could have evolved different receptor selectivity. Together, these data show that NPS is a relatively recent gene in vertebrate evolution whose appearance might coincide with its specialized physiological functions in terrestrial vertebrates.

Species-specific evolution of the FcR family in endothermic vertebrates

In primates and rodents, the extended FcR family is comprised of three subsets: classical FcRs, structurally diverse cell surface receptors currently designated FCRL1-FCRL6, and intracellular proteins FCRLA and FCRLB. Using bioinformatic analysis, we revealed the FcR-like genes of the same three subsets in the genome of dog, another representative of placental mammals, and in the genome of short-tailed opossum, a representative of marsupials. In contrast, a single FcR-like gene was found in the current version of the chicken genome. This in silico finding was confirmed by the gene cloning and subsequent Southern blot hybridization. The chicken FCRL gene encodes a cell surface receptor with the extracellular region composed of four Ig-like domains of the D1-, D2-, D3-, and D4-subtypes. The gene is expressed in lymphoid and non-lymphoid tissues. Phylogenetic analysis of the mammalian and chicken genes suggested that classical FcRs, FCRLA, and FCRLB emerged after the mammalian-avian split but before the eutherian-marsupial radiation. The data obtained show that the repertoire of the classical FcRs and surface FcR-like proteins in mammalian species was shaped by an extensive recombination process, which resulted in domain shuffling and species-specific gain and loss of distinct exons or entire genes.

Cranial endocasts from a growth series ofMonodelphis domestica (Didelphidae, Marsupialia): A study of individual and ontogenetic variation

Intraspecific variation (e.g., ontogenetic, individual, sexual dimorphic) is rarely examined among cranial endocasts (infillings of the braincase cavity) because of the difficulty in obtaining multiple specimens of a species, particularly fossil taxa. We extracted digital cranial endocasts from CT scans of a growth series of skulls of Monodelphis domestica, the gray short-tailed opossum, as a preliminary assessment of the amount of intraspecific variation in mammalian endocranial morphology. The goals of this study were 1) to provide an anatomical description to document developmental changes in endocranial morphology of M. domestica and 2) to examine ontogenetic and individual variation with respect to phylogenetic characters of endocranial cavities that are known to be variable between different mammalian taxa. In this study, "ontogenetic variation" refers to variation between specimens of different ages whereas "individual variation" (i.e., polymorphism) is restricted to variation between specimens of comparable age. Aside from size, changes in shape account for the greatest amount of morphological variation between the endocasts of different ages. Endocast length, width, and volume increase with age for the growth series. Relative olfactory bulb cast size increases with age in the growth series, but the relative size of the parafloccular casts shows a slight negative allometric trend through ontogeny. More than one-third of the phylogenetic characters of the endocranial cavity we examined showed some sort of variation (ontogenetic, individual, or both). This suggests that although endocasts are potentially informative for systematics, both ontogenetic and individual variation affect how endocranial characters are scored for phylogenetic analysis. Further studies such as this are necessary to determine the taxonomic extent of significant intraspecific variation of these endocranial characters.

Evolution and comparative analysis of the MHC Class III inflammatory region

BACKGROUND

The Major Histocompatibility Complex (MHC) is essential for immune function. Historically, it has been subdivided into three regions (Class I, II, and III), but a cluster of functionally related genes within the Class III region has also been referred to as the Class IV region or "inflammatory region". This group of genes is involved in the inflammatory response, and includes members of the tumour necrosis family. Here we report the sequencing, annotation and comparative analysis of a tammar wallaby BAC containing the inflammatory region. We also discuss the extent of sequence conservation across the entire region and identify elements conserved in evolution.

RESULTS

Fourteen Class III genes from the tammar wallaby inflammatory region were characterised and compared to their orthologues in other vertebrates. The organisation and sequence of genes in the inflammatory region of both the wallaby and South American opossum are highly conserved compared to known genes from eutherian ("placental") mammals. Some minor differences separate the two marsupial species. Eight genes within the inflammatory region have remained tightly clustered for at least 360 million years, predating the divergence of the amphibian lineage. Analysis of sequence conservation identified 354 elements that are conserved. These range in size from 7 to 431 bases and cover 15.6% of the inflammatory region, representing approximately a 4-fold increase compared to the average for vertebrate genomes. About 5.5% of this conserved sequence is marsupial-specific, including three cases of marsupial-specific repeats. Highly Conserved Elements were also characterised.

CONCLUSION

Using comparative analysis, we show that a cluster of MHC genes involved in inflammation, including TNF, LTA (or its putative teleost homolog TNF-N), APOM, and BAT3 have remained together for over 450 million years, predating the divergence of mammals from fish. The observed enrichment in conserved sequences within the inflammatory region suggests conservation at the transcriptional regulatory level, in addition to the functional level.

Ancient genomic architecture for mammalian olfactory receptor clusters

A new tool for genome-wide definition of genomic gene clusters conserved in multiple species was applied to olfactory receptors in five mammals, demonstrating that most mammalian olfactory receptor clusters have a common ancestry.

Background

Mammalian olfactory receptor (OR) genes reside in numerous genomic clusters of up to several dozen genes. Whole-genome sequence alignment nets of five mammals allow their comprehensive comparison, aimed at reconstructing the ancestral olfactory subgenome.

Results

We developed a new and general tool for genome-wide definition of genomic gene clusters conserved in multiple species. Syntenic orthologs, defined as gene pairs showing conservation of both genomic location and coding sequence, were subjected to a graph theory algorithm for discovering CLICs (clusters in conservation). When applied to ORs in five mammals, including the marsupial opossum, more than 90% of the OR genes were found within a framework of 48 multi-species CLICs, invoking a general conservation of gene order and composition. A detailed analysis of individual CLICs revealed multiple differences among species, interpretable through species-specific genomic rearrangements and reflecting complex mammalian evolutionary dynamics. One significant instance involves CLIC #1, which lacks a human member, implying the human-specific deletion of an OR cluster, whose mouse counterpart has been tentatively associated with isovaleric acid odorant detection.

Conclusion

The identified multi-species CLICs demonstrate that most of the mammalian OR clusters have a common ancestry, preceding the split between marsupials and placental mammals. However, only two of these CLICs were capable of incorporating chicken OR genes, parsimoniously implying that all other CLICs emerged subsequent to the avian-mammalian divergence.

Comparative analysis of the paired immunoglobulin-like receptor (PILR) locus in six mammalian genomes: duplication, conversion, and the birth of new genes

Manyaspects of the immune system are controlled by homologous cell surface receptors that mediate inhibitory and activating pathways. The paired immunoglobulin-like receptor (PILR) locus at 7q22 encodes both PILRA, an inhibitory receptor, and PILRB, its activating counterpart. Mouse Pilrb1 is a novel immune system regulator, and its ligand Cd99 participates in the recruitment of T-cells to inflamed tissue. We characterized the PILR locus in six mammalian genomes and investigated the structure and mRNA expression of human PILRB. Synteny at the PILR locus is conserved in the human, chimpanzee, dog, mouse and rat genomes. The absence of the PILR locus in opossum and chicken genomes suggests it arose after the divergence of placental and nonplacental mammals. In humans, a Williams-Beuren syndrome-related segmental duplication has created a complex chimeric transcript representing the predominantly expressed form of PILRB. Unlike PILRA, PILRB transcripts were detected in a wide variety of tissues including cells of the lymphoid lineage. In the mouse genome, a second activating gene, Pilrb2, and six pseudogenes were found. Extensive gene duplications in the rat genome have resulted in at least 27 Pilrb genes and or pseudogenes. Abundant gene duplication events involving novel CD99-related genes were also detected in the rat genome. In addition to duplication, we show that gene conversion has played a persistent role in the evolution of the PILR genes. Overall, we demonstrate that the PILR locus is dynamically evolving via multiple evolutionary mechanisms in several mammalian genomes.

Germ cell restricted expression of chick Nanog

Nanog is a pluripotency-associated factor expressed in embryonic stem cells and in the epiblast and primordial germ cells of the mouse embryo. We have identified the chick orthologue of Nanog and found that its expression is limited to primordial germ cells in the early embryo and is not found throughout the epiblast. Genomic analysis has shown that Nanog is an amniote-specific gene and is absent from anamniotes and invertebrates. Furthermore, other pluripotency associated genes that are located in close proximity to Nanog in human and mouse are absent from the chick genome. Such observations lead to a scenario of sequential addition of novel genes to a genomic region associated to pluripotency. These results have profound implications for the study of the evolution of pluripotent lineages in the embryo and of vertebrate stem cells.

Using several pair-wise informant sequences for de novo prediction of alternatively spliced transcripts

Background

As part of the ENCODE Genome Annotation Assessment Project (EGASP), we developed the MARS extension to the Twinscan algorithm. MARS is designed to find human alternatively spliced transcripts that are conserved in only one or a limited number of extant species. MARS is able to use an arbitrary number of informant sequences and predicts a number of alternative transcripts at each gene locus.

Results

MARS uses the mouse, rat, dog, opossum, chicken, and frog genome sequences as pairwise informant sources for Twinscan and combines the resulting transcript predictions into genes based on coding (CDS) region overlap. Based on the EGASP assessment, MARS is one of the more accurate dual-genome prediction programs. Compared to the GENCODE annotation, we find that predictive sensitivity increases, while specificity decreases, as more informant species are used. MARS correctly predicts alternatively spliced transcripts for 11 of the 236 multi-exon GENCODE genes that are alternatively spliced in the coding region of their transcripts. For these genes a total of 24 correct transcripts are predicted.

Conclusion

The MARS algorithm is able to predict alternatively spliced transcripts without the use of expressed sequence information, although the number of loci in which multiple predicted transcripts match multiple alternatively spliced transcripts in the GENCODE annotation is relatively small.

The evolving doublecortin (DCX) superfamily

Background

Doublecortin (DCX) domains serve as protein-interaction platforms. Mutations in members of this protein superfamily are linked to several genetic diseases. Mutations in the human DCX gene result in abnormal neuronal migration, epilepsy, and mental retardation; mutations in RP1 are associated with a form of inherited blindness, and DCDC2 has been associated with dyslectic reading disabilities.

Results

The DCX-repeat gene family is composed of eleven paralogs in human and in mouse. Its evolution was followed across vertebrates, invertebrates, and was traced to unicellular organisms, thus enabling following evolutionary additions and losses of genes or domains. The N-terminal and C-terminal DCX domains have undergone sub-specialization and divergence. Developmental in situ hybridization data for nine genes was generated. In addition, a novel co-expression analysis for most human and mouse DCX superfamily-genes was performed using high-throughput expression data extracted from Unigene. We performed an in-depth study of a complete gene superfamily using several complimentary methods.

Conclusion

This study reveals the existence and conservation of multiple members of the DCX superfamily in different species. Sequence analysis combined with expression analysis is likely to be a useful tool to predict correlations between human disease and mouse models. The sub-specialization of some members due to restricted expression patterns and sequence divergence may explain the successful addition of genes to this family throughout evolution.

New members of the neurexin superfamily: multiple rodent homologues of the human CASPR5 gene

Proteins of the Caspr family are involved in cell contacts and communication in the nervous system. We identified and, by in silico reconstruction, compiled three orthologues of the human CASPR5 gene from the mouse genome, four from the rat genome, and one each from the chimpanzee, dog, opossum, and chicken genomes. Obviously, Caspr5 gene duplications have taken place during evolution of the rodent lineage. In the rat, the four paralogues are located in one chromosome arm, Chr 13p. In the mouse, however, the three Caspr5 genes are located in two chromosomes, Chr 1 and Chr 17. RT-PCR shows that all three mouse paralogues are being expressed. Common expression is found in brain tissue but different expression patterns are seen in other organs during fetal development and in the adult stage. Tissue specificity of expression has diverged during evolution of this young rodent gene family.

The Xist RNA gene evolved in eutherians by pseudogenization of a protein-coding gene

The Xist noncoding RNA is the key initiator of the process of X chromosome inactivation in eutherian mammals, but its precise function and origin remain unknown. Although Xist is well conserved among eutherians, until now, no homolog has been identified in other mammals. We show here that Xist evolved, at least partly, from a protein-coding gene and that the loss of protein-coding function of the proto-Xist coincides with the four flanking protein genes becoming pseudogenes. This event occurred after the divergence between eutherians and marsupials, which suggests that mechanisms of dosage compensation have evolved independently in both lineages.

Comparative phylogenetic analysis reveals multiple non-imprinted isoforms of opossum Dlk1

Imprinted genes are monoallelically expressed in a parent-of-origin manner and were previously identified in both marsupials and eutherians, but not in monotremes. Phylogenetic comparison of imprinted domains is a powerful tool for investigating the molecular and adaptive evolution of this unique gene regulatory mechanism. Herein, we report that multiple transcripts of Dlk1 (Delta, Drosophila, Homolog-like 1) are expressed in the opossum, but none are imprinted. Thus, we provide the first example of a reciprocally imprinted gene domain in which imprinting evolved in a common ancestor to eutherian rather than therian mammals. Moreover, the reciprocally imprinted Meg3 (Maternally Expressed Gene 3), found downstream of Dlk1 in eutherian mammals, is absent in the opossum. We propose that the Meg3 sequence integrated into the eutherian Dlk1 domain via a LINE-1 element and that Dlk1 became imprinted in eutherian mammals only after this downstream integration. These findings clearly demonstrate that imprinted genes did not all evolve before the divergence of marsupials and eutherians.

Mannose 6-phosphate receptors in an ancient vertebrate, zebrafish

The endosome/lysosome system plays key roles in embryonic development, but difficulties posed by inaccessible mammalian embryos have hampered detailed studies. The accessible, transparent embryos of Danio rerio, together with the genetic and experimental approaches possible with this organism, provide many advantages over rodents. In mammals, mannose 6-phosphate receptors (MPRs) target acid hydrolases to endosomes and lysosomes, but nothing is known of acid hydrolase targeting in zebrafish. Here, we describe the sequence of the zebrafish cation-dependent MPR (CD-MPR) and cation-independent MPR (CI-MPR), and compare them with their mammalian orthologs. We show that all residues critical for mannose 6-phosphate (M6P) recognition are present in the extracellular domains of the zebrafish receptors, and that trafficking signals in the cytoplasmic tails are also conserved. This suggests that the teleost receptors possess M6P binding sites with properties similar to those of mammalian MPRs, and that targeting of lysosomal enzymes by MPRs represents an ancient pathway in vertebrate cell biology. We also determined the expression patterns of the CD-MPR and CI-MPR during embryonic development in zebrafish. Both genes are expressed from the one-cell stage through to the hatching period. In early embryos, expression is ubiquitous, but in later stages, expression of both receptors is restricted to the anterior region of the embryo, covering the forebrain, midbrain and hindbrain. The expression patterns suggest time- and tissue-specific functions for the receptors, with particular evidence for roles in neural development. Our study establishes zebrafish as a novel, genetically tractable model for in vivo studies of MPR function and lysosome biogenesis.

Transposon-free regions in mammalian genomes

Despite the presence of over 3 million transposons separated on average by approximately 500 bp, the human and mouse genomes each contain almost 1000 transposon-free regions (TFRs) over 10 kb in length. The majority of human TFRs correlate with orthologous TFRs in the mouse, despite the fact that most transposons are lineage specific. Many human TFRs also overlap with orthologous TFRs in the marsupial opossum, indicating that these regions have remained refractory to transposon insertion for long evolutionary periods. Over 90% of the bases covered by TFRs are noncoding, much of which is not highly conserved. Most TFRs are not associated with unusual nucleotide composition, but are significantly associated with genes encoding developmental regulators, suggesting that they represent extended regions of regulatory information that are largely unable to tolerate insertions, a conclusion difficult to reconcile with current conceptions of gene regulation.

New DNA data from a transthyretin nuclear intron suggest an Oligocene to Miocene diversification of living South America opossums (Marsupialia: Didelphidae)

Phylogenetic relationships of 19 species of didelphid marsupials were studied using two nuclear markers, the non-coding transthyretin intron 1 (TTR) and the coding interphotoreceptor retinoid binding protein exon 1 (IRBP), and two mitochondrial genes, the protein-coding cytochrome b (cyt-b) and the structural 12S ribosomal DNA (12S rDNA). Evolutionary dynamics of these four markers were compared to each other, revealing the appropriate properties presented by TTR intron 1 together with its well supported and resolved phylogenetic signal. Nuclear markers supported the monophyly of medium and large-sized opossums Metachirus+(Chironectes, Lutreolina, Didelphis, Philander), and the paraphyly of mouse-sized opossums, with the genera Gracilinanus, Thylamys, and Marmosops as a sister group to medium and large-sized didelphids. Conflicting branching patterns between mitochondrial and nuclear data involved the phylogenetic position of Marmosa-Micoureus-Monodelphis relative to other mouse-sized opossums. Nuclear phylogenetic inferences among genera were confirmed by the presence of synapomorphic indels observed in TTR intron 1. A Bayesian relaxed molecular clock dating of didelphid evolution using nuclear markers estimated their origin in the Middle Eocene (39.8 million years ago), with subsequent diversification during the Oligocene (Deseadan) and Miocene.