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

Donkey genome and insight into the imprinting of fast karyotype evolution

The donkey, like the horse, is a promising model for exploring karyotypic instability. We report the de novo whole-genome assemblies of the donkey and the Asiatic wild ass. Our results reflect the distinct characteristics of donkeys, including more effective energy metabolism and better immunity than horses. The donkey shows a steady demographic trajectory. We detected abundant satellite sequences in some inactive centromere regions but not in neocentromere regions, while ribosomal RNAs frequently emerged in neocentromere regions but not in the obsolete centromere regions. Expanded miRNA families and five newly discovered miRNA target genes involved in meiosis may be associated with fast karyotype evolution. APC/C, controlling sister chromatid segregation, cytokinesis, and the establishment of the G1 cell cycle phase were identified by analysis of miRNA targets and rapidly evolving genes.

Evolution of Vertebrate Adam Genes; Duplication of Testicular Adams from Ancient Adam9/9-like Loci

Members of the disintegrin metalloproteinase (ADAM) family have important functions in regulating cell-cell and cell-matrix interactions as well as cell signaling. There are two major types of ADAMs: the somatic ADAMs (sADAMs) that have a significant presence in somatic tissues, and the testicular ADAMs (tADAMs) that are expressed predominantly in the testis. Genes encoding tADAMs can be further divided into two groups: group I (intronless) and group II (intron-containing). To date, tAdams have only been reported in placental mammals, and their evolutionary origin and relationship to sAdams remain largely unknown. Using phylogenetic and syntenic tools, we analyzed the Adam genes in various vertebrates ranging from fishes to placental mammals. Our analyses reveal duplication and loss of some sAdams in certain vertebrate species. In particular, there exists an Adam9-like gene in non-mammalian vertebrates but not mammals. We also identified putative group I and group II tAdams in all amniote species that have been examined. These tAdam homologues are more closely related to Adams 9 and 9-like than to other sAdams. In all amniote species examined, group II tAdams lie in close vicinity to Adam9 and hence likely arose from tandem duplication, whereas group I tAdams likely originated through retroposition because of their lack of introns. Clusters of multiple group I tAdams are also common, suggesting tandem duplication after retroposition. Therefore, Adam9/9-like and some of the derived tAdam loci are likely preferred targets for tandem duplication and/or retroposition. Consistent with this hypothesis, we identified a young retroposed gene that duplicated recently from Adam9 in the opossum. As a result of gene duplication, some tAdams were pseudogenized in certain species, whereas others acquired new expression patterns and functions. The rapid duplication of Adam genes has a major contribution to the diversity of ADAMs in various vertebrate species.

Draft De Novo Transcriptome of the Rat Kangaroo Potorous tridactylus as a Tool for Cell Biology

The rat kangaroo (long-nosed potoroo, Potorous tridactylus) is a marsupial native to Australia. Cultured rat kangaroo kidney epithelial cells (PtK) are commonly used to study cell biological processes. These mammalian cells are large, adherent, and flat, and contain large and few chromosomes—and are thus ideal for imaging intra-cellular dynamics such as those of mitosis. Despite this, neither the rat kangaroo genome nor transcriptome have been sequenced, creating a challenge for probing the molecular basis of these cellular dynamics. Here, we present the sequencing, assembly and annotation of the draft rat kangaroo de novo transcriptome. We sequenced 679 million reads that mapped to 347,323 Trinity transcripts and 20,079 Unigenes. We present statistics emerging from transcriptome-wide analyses, and analyses suggesting that the transcriptome covers full-length sequences of most genes, many with multiple isoforms. We also validate our findings with a proof-of-concept gene knockdown experiment. We expect that this high quality transcriptome will make rat kangaroo cells a more tractable system for linking molecular-scale function and cellular-scale dynamics.

Marsupials and monotremes possess a novel family of MHC class I genes that is lost from the eutherian lineage

Background

Major histocompatibility complex (MHC) class I genes are found in the genomes of all jawed vertebrates. The evolution of this gene family is closely tied to the evolution of the vertebrate genome. Family members are frequently found in four paralogous regions, which were formed in two rounds of genome duplication in the early vertebrates, but in some species class Is have been subject to additional duplication or translocation, creating additional clusters. The gene family is traditionally grouped into two subtypes: classical MHC class I genes that are usually MHC-linked, highly polymorphic, expressed in a broad range of tissues and present endogenously-derived peptides to cytotoxic T-cells; and non-classical MHC class I genes generally have lower polymorphism, may have tissue-specific expression and have evolved to perform immune-related or non-immune functions. As immune genes can evolve rapidly and are subject to different selection pressure, we hypothesised that there may be divergent, as yet unannotated or uncharacterised class I genes.

Results

Application of a novel method of sensitive genome searching of available vertebrate genome sequences revealed a new, extensive sub-family of divergent MHC class I genes, denoted as UT, which has not previously been characterized. These class I genes are found in both American and Australian marsupials, and in monotremes, at an evolutionary chromosomal breakpoint, but are not present in non-mammalian genomes and have been lost from the eutherian lineage. We show that UT family members are expressed in the thymus of the gray short-tailed opossum and in other immune tissues of several Australian marsupials. Structural homology modelling shows that the proteins encoded by this family are predicted to have an open, though short, antigen-binding groove.

Conclusions

We have identified a novel sub-family of putatively non-classical MHC class I genes that are specific to marsupials and monotremes. This family was present in the ancestral mammal and is found in extant marsupials and monotremes, but has been lost from the eutherian lineage. The function of this family is as yet unknown, however, their predicted structure may be consistent with presentation of antigens to T-cells.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1745-4) contains supplementary material, which is available to authorized users.

Physical Mapping and Refinement of the Painted Turtle Genome (Chrysemys picta) Inform Amniote Genome Evolution and Challenge Turtle-Bird Chromosomal Conservation

Comparative genomics continues illuminating amniote genome evolution, but for many lineages our understanding remains incomplete. Here, we refine the assembly (CPI 3.0.3 NCBI AHGY00000000.2) and develop a cytogenetic map of the painted turtle (Chrysemys picta-CPI) genome, the first in turtles and in vertebrates with temperature-dependent sex determination. A comparison of turtle genomes with those of chicken, selected nonavian reptiles, and human revealed shared and novel genomic features, such as numerous chromosomal rearrangements. The largest conserved syntenic blocks between birds and turtles exist in four macrochromosomes, whereas rearrangements were evident in these and other chromosomes, disproving that turtles and birds retain fully conserved macrochromosomes for greater than 300 Myr. C-banding revealed large heterochromatic blocks in the centromeric region of only few chromosomes. The nucleolar-organizing region (NOR) mapped to a single CPI microchromosome, whereas in some turtles and lizards the NOR maps to nonhomologous sex-chromosomes, thus revealing independent translocations of the NOR in various reptilian lineages. There was no evidence for recent chromosomal fusions as interstitial telomeric-DNA was absent. Some repeat elements (CR1-like, Gypsy) were enriched in the centromeres of five chromosomes, whereas others were widespread in the CPI genome. Bacterial artificial chromosome (BAC) clones were hybridized to 18 of the 25 CPI chromosomes and anchored to a G-banded ideogram. Several CPI sex-determining genes mapped to five chromosomes, and homology was detected between yet other CPI autosomes and the globally nonhomologous sex chromosomes of chicken, other turtles, and squamates, underscoring the independent evolution of vertebrate sex-determining mechanisms.

Evolution of V genes from the TRV loci of mammals

Information concerning the evolution of T lymphocyte receptors (TCR) can be deciphered from that part of the molecule that recognizes antigen presented by major histocompatibility complex (MHC), namely the variable (V) regions. The genes that code for these variable regions are found within the TCR loci. Here, we describe a study of the evolutionary origin of V genes that code for the α and β chains of the TCR loci of mammals. In particular, we demonstrate that most of the 35 TRAV and 25 TRBV conserved genes found in Primates are also found in other Eutheria, while in Marsupials, Monotremes, and Reptiles, these genes diversified in a different manner. We also show that in mammals, all TRAV genes are derived from five ancestral genes, while all TRBV genes originate from four such genes. In Reptiles, the five TRAV and three out of the four TRBV ancestral genes exist, as well as other V genes not found in mammals. We also studied the TRGV and TRDV loci from all mammals, and we show a relationship of the TRDV to the TRAV locus throughout evolutionary time.

Comparative genomics for biodiversity conservation

Genomic approaches are gathering momentum in biology and emerging opportunities lie in the creative use of comparative molecular methods for revealing the processes that influence diversity of wildlife. However, few comparative genomic studies are performed with explicit and specific objectives to aid conservation of wild populations. Here I provide a brief overview of comparative genomic approaches that offer specific benefits to biodiversity conservation. Because conservation examples are few, I draw on research from other areas to demonstrate how comparing genomic data across taxa may be used to inform the characterisation of conservation units and studies of hybridisation, as well as studies that provide conservation outcomes from a better understanding of the drivers of divergence. A comparative approach can also provide valuable insight into the threatening processes that impact rare species, such as emerging diseases and their management in conservation. In addition to these opportunities, I note areas where additional research is warranted. Overall, comparing and contrasting the genomic composition of threatened and other species provide several useful tools for helping to preserve the molecular biodiversity of the global ecosystem.

LINEs in mice: features, families, and potential roles in early development

Approximately half of the mammalian genome is composed of repetitive elements, including LINE-1 (L1) elements. Because of their potential ability to transpose and integrate into other regions of the genome, their activation represents a threat to genome stability. Molecular pathways have emerged to tightly regulate and repress their transcriptional activity, including DNA methylation, histone modifications, and RNA pathways. It has become evident that Line-L1 elements are evolutionary diverse and dedicated repression pathways have been recently uncovered that discriminate between evolutionary old and young elements, with RNA-directed silencing mechanisms playing a prominent role. During periods of epigenetic reprogramming in development, specific classes of repetitive elements are upregulated, presumably due to the loss of most heterochromatic marks in this process. While we have learnt a lot on the molecular mechanisms that regulate Line-L1 expression over the last years, it is still unclear whether reactivation of Line-L1 after fertilization serves a functional purpose or it is a simple side effect of reprogramming.

Differential SINE evolution in vesper and non-vesper bats

Background

Short interspersed elements (SINEs) have a powerful influence on genome evolution and can be useful markers for phylogenetic inference and population genetic analyses. In this study, we examined survey sequence and whole genome data to determine the evolutionary dynamics of Ves SINEs in the genomes of 11 bats, nine from Vespertilionidae.

Results

We identified 41 subfamilies of Ves and linked several to specific lineages. We also revealed substantial differences among lineages including the observation that Ves accumulation and Ves subfamily diversity is significantly higher in vesper as opposed to non-vesper bats. This is especially interesting when one considers the increased transposable element diversity of vesper bats in general.

Conclusions

Our data suggest that survey sequencing and genome mining are valuable tools to investigate SINE evolution among related lineages and can provide substantial information about the ability of SINEs to proliferate in diverse genomes. This method would also be a useful first step in determining which subfamilies would be the best to target when developing SINEs as markers for phylogenetic and population genetic analyses.

Electronic supplementary material

The online version of this article (doi:10.1186/s13100-015-0038-4) contains supplementary material, which is available to authorized users.

Principles of long noncoding RNA evolution derived from direct comparison of transcriptomes in 17 species

The inability to predict long noncoding RNAs from genomic sequence has impeded the use of comparative genomics for studying their biology. Here, we develop methods that use RNA sequencing (RNA-seq) data to annotate the transcriptomes of 16 vertebrates and the echinoid sea urchin, uncovering thousands of previously unannotated genes, most of which produce long intervening noncoding RNAs (lincRNAs). Although in each species, >70% of lincRNAs cannot be traced to homologs in species that diverged >50 million years ago, thousands of human lincRNAs have homologs with similar expression patterns in other species. These homologs share short, 5'-biased patches of sequence conservation nested in exonic architectures that have been extensively rewired, in part by transposable element exonization. Thus, over a thousand human lincRNAs are likely to have conserved functions in mammals, and hundreds beyond mammals, but those functions require only short patches of specific sequences and can tolerate major changes in gene architecture.

Genetic changes shaping the human brain

The development and function of our brain are governed by a genetic blueprint, which reflects dynamic changes over the history of evolution. Recent progress in genetics and genomics, facilitated by next-generation sequencing and single-cell sorting, has identified numerous genomic loci that are associated with a neuroanatomical or neurobehavioral phenotype. Here, we review some of the genetic changes in both protein-coding and noncoding regions that affect brain development and evolution, as well as recent progress in brain transcriptomics. Understanding these genetic changes may provide novel insights into neurological and neuropsychiatric disorders, such as autism and schizophrenia.

Interspersed repeats in the horse (Equus caballus); spatial correlations highlight conserved chromosomal domains

The interspersed repeat content of mammalian genomes has been best characterized in human, mouse and cow. In this study, we carried out de novo identification of repeated elements in the equine genome and identified previously unknown elements present at low copy number. The equine genome contains typical eutherian mammal repeats, but also has a significant number of hybrid repeats in addition to clade-specific Long Interspersed Nuclear Elements (LINE). Equus caballus clade specific LINE 1 (L1) repeats can be classified into approximately five subfamilies, three of which have undergone significant expansion. There are 1115 full-length copies of these equine L1, but of the 103 presumptive active copies, 93 fall within a single subfamily, indicating a rapid recent expansion of this subfamily. We also analysed both interspersed and simple sequence repeats (SSR) genome-wide, finding that some repeat classes are spatially correlated with each other as well as with G+C content and gene density. Based on these spatial correlations, we have confirmed that recently-described ancestral vs. clade-specific genome territories can be defined by their repeat content. The clade-specific Short Interspersed Nuclear Element correlations were scattered over the genome and appear to have been extensively remodelled. In contrast, territories enriched for ancestral repeats tended to be contiguous domains. To determine if the latter territories were evolutionarily conserved, we compared these results with a similar analysis of the human genome, and observed similar ancestral repeat enriched domains. These results indicate that ancestral, evolutionarily conserved mammalian genome territories can be identified on the basis of repeat content alone. Interspersed repeats of different ages appear to be analogous to geologic strata, allowing identification of ancient vs. newly remodelled regions of mammalian genomes.

Disentangling the relationship of the Australian marsupial orders using retrotransposon and evolutionary network analyses

The ancestors to the Australian marsupials entered Australia around 60 (54-72) Ma from Antarctica, and radiated into the four living orders Peramelemorphia, Dasyuromorphia, Diprotodontia, and Notoryctemorphia. The relationship between the four Australian marsupial orders has been a long-standing question, because different phylogenetic studies have not been able to consistently reconstruct the same topology. Initial in silico analysis of the Tasmanian devil genome and experimental screening in the seven marsupial orders revealed 20 informative transposable element insertions for resolving the inter- and intraordinal relationships of Australian and South American orders. However, the retrotransposon insertions support three conflicting topologies regarding Peramelemorphia, Dasyuromorphia, and Notoryctemorphia, indicating that the split between the three orders may be best understood as a network. This finding is supported by a phylogenetic reanalysis of nuclear gene sequences, using a consensus network approach that allows depicting hidden phylogenetic conflict, otherwise lost when forcing the data into a bifurcating tree. The consensus network analysis agrees with the transposable element analysis in that all possible topologies regarding Peramelemorphia, Dasyuromorphia, and Notoryctemorphia in a rooted four-taxon topology are equally well supported. In addition, retrotransposon insertion data support the South American order Didelphimorphia being the sistergroup to all other living marsupial orders. The four Australian orders originated within 3 Myr at the Cretaceous-Paleogene boundary. The rapid divergences left conflicting phylogenetic information in the genome possibly generated by incomplete lineage sorting or introgressive hybridization, leaving the relationship among Australian marsupial orders unresolvable as a bifurcating process millions of years later.

Immunomics of the koala (Phascolarctos cinereus)

The study of the koala transcriptome has the potential to advance our understanding of its immunome--immunological reaction of a given host to foreign antigens--and to help combat infectious diseases (e.g., chlamydiosis) that impede ongoing conservation efforts. We used Illumina sequencing of cDNA to characterize genes expressed in two different koala tissues of immunological importance, blood and spleen. We generated nearly 600 million raw sequence reads, and about 285 million of these were subsequently assembled and condensed into ~70,000 subcomponents that represent putative transcripts. We annotated ~16% of these subcomponents and identified those related to infection and the immune response, including Toll-like receptors (TLRs), RIG-I-like receptors (RLRs), major histocompatibility complex (MHC) genes, and koala retrovirus (KoRV). Using phylogenetic analyses, we identified 29 koala genes in these target categories and report their concordance with currently accepted gene groups. By mapping multiple sequencing reads to transcripts, we identified 56 putative SNPs in genes of interest. The distribution of these SNPs indicates that MHC genes (34 SNPs) are more diverse than KoRV (12 SNPs), TLRs (8 SNPs), or RLRs (2 SNPs). Our sequence data also indicate that KoRV sequences are highly expressed in the transcriptome. Our efforts have produced full-length sequences for potentially important immune genes in koala, which should serve as targets for future investigations that aim to conserve koala populations.

The opossum MHC genomic region revisited

The gray short-tailed opossum Monodelphis domestica is one of the few marsupial species for which a high quality whole genome sequence is available and the major histocompatibility complex (MHC) region has been annotated. Previous analyses revealed only a single locus within the opossum MHC region, designated Modo-UA1, with the features expected for encoding a functionally classical class I α-chain. Nine other class I genes found within the MHC are highly divergent and have features usually associated with non-classical roles. The original annotation, however, was based on an early version of the opossum genome assembly. More recent analyses of allelic variation in individual opossums revealed too many Modo-UA1 sequences per individual to be accounted for by a single MHC class I locus found in the genome assembly. A reanalysis of a later generation assembly, MonDom5, revealed the presence of two additional loci, now designated Modo-UA3 and UA4, in a region that was expanded and more complete than in the earlier assembly. Modo-UA1, UA3, and UA4 are all transcribed, although Modo-UA4 transcripts are rarer. Modo-UA4 is also relatively non-polymorphic. Evidence presented support the accuracy of the later assembly and the existence of three related class I genes in the opossum, making opossums more typical of mammals and most tetrapods by having multiple apparent classical MHC class I loci.

Exogenous retinoic acid induces digit reduction in opossums (Monodelphis domestica) by disrupting cell death and proliferation, and apical ectodermal ridge and zone of polarizing activity function

BACKGROUND

Retinoic acid (RA) is a vitamin A derivative. Exposure to exogenous RA generates congenital limb malformations (CLMs) in species from frogs to humans. These CLMs include but are not limited to oligodactyly and long-bone hypoplasia. The processes by which exogenous RA induces CLMs in mammals have been best studied in mouse, but as of yet remain unresolved.

METHODS

We investigated the impact of exogenous RA on the cellular and molecular development of the limbs of a nonrodent model mammal, the opossum Monodelphis domestica. Opossums exposed to exogenous retinoic acid display CLMs including oligodactly, and results are consistent with opossum development being more susceptible to RA-induced disruptions than mouse development.

RESULTS

Exposure of developing opossums to exogenous RA leads to an increase in cell death in the limb mesenchyme that is most pronounced in the zone of polarizing activity, and a reduction in cell proliferation throughout the limb mesenchyme. Exogenous RA also disrupts the expression of Shh in the zone of polarizing activity, and Fgf8 in the apical ectodermal ridge, and other genes with roles in the regulation of limb development and cell death.

CONCLUSION

Results are consistent with RA inducing CLMs in opossum limbs by disrupting the functions of the apical ectodermal ridge and zone of polarizing activity, and driving an increase in cell death and reduction of cell proliferation in the mesenchyme of the developing limb.

Ancient transposable elements transformed the uterine regulatory landscape and transcriptome during the evolution of mammalian pregnancy

A major challenge in biology is determining how evolutionarily novel characters originate; however, mechanistic explanations for the origin of new characters are almost completely unknown. The evolution of pregnancy is an excellent system in which to study the origin of novelties because mammals preserve stages in the transition from egg laying to live birth. To determine the molecular bases of this transition, we characterized the pregnant/gravid uterine transcriptome from tetrapods to trace the evolutionary history of uterine gene expression. We show that thousands of genes evolved endometrial expression during the origins of mammalian pregnancy, including genes that mediate maternal-fetal communication and immunotolerance. Furthermore, thousands of cis-regulatory elements that mediate decidualization and cell-type identity in decidualized stromal cells are derived from ancient mammalian transposable elements (TEs). Our results indicate that one of the defining mammalian novelties evolved from DNA sequences derived from ancient mammalian TEs coopted into hormone-responsive regulatory elements distributed throughout the genome.

Evolutionary histories of transposable elements in the genome of the largest living marsupial carnivore, the Tasmanian devil

The largest living carnivorous marsupial, the Tasmanian devil (Sarcophilus harrisii), is the sole survivor of a lineage originating about 12 Ma. We set out to investigate the spectrum of transposable elements found in the Tasmanian devil genome, the first high-coverage genome of an Australian marsupial. Marsupial genomes have been shown to have the highest amount of transposable elements among vertebrates. We analyzed the horizontally transmitted DNA transposons OC1 and hAT-1_MEu in the Tasmanian devil genome. OC1 is present in all carnivorous marsupials, while having a very limited distribution among the remaining Australian marsupial orders. In contrast, hAT-1_MEu is present in all Australian marsupial orders, and has so far only been identified in a few placental mammals. We screened 158 introns for phylogenetically informative retrotransposons in the order Dasyuromorphia, and found that the youngest SINE (Short INterspersed Element), WSINE1, is no longer active in the subfamily Dasyuridae. The lack of detectable WSINE1 activity in this group may be due to a retrotransposon inactivation event approximately 30 Ma. We found that the Tasmanian devil genome contains a relatively low number of continuous full-length LINE-1 (Long INterspersed Element 1, L1) retrotransposons compared with the opossum genome. Furthermore, all L1 elements in the Tasmanian devil appeared to be nonfunctional. Hidden Markov Model approaches suggested that other potential sources of functional reverse transcriptase are absent from the genome. We discuss the issues associated with assembling long, highly similar L1 copies from short read Illumina data and describe how assembly artifacts can potentially lead to erroneous conclusions.

The Genome 10K Project: a way forward

The Genome 10K Project was established in 2009 by a consortium of biologists and genome scientists determined to facilitate the sequencing and analysis of the complete genomes of 10,000 vertebrate species. Since then the number of selected and initiated species has risen from ∼26 to 277 sequenced or ongoing with funding, an approximately tenfold increase in five years. Here we summarize the advances and commitments that have occurred by mid-2014 and outline the achievements and present challenges of reaching the 10,000-species goal. We summarize the status of known vertebrate genome projects, recommend standards for pronouncing a genome as sequenced or completed, and provide our present and future vision of the landscape of Genome 10K. The endeavor is ambitious, bold, expensive, and uncertain, but together the Genome 10K Consortium of Scientists and the worldwide genomics community are moving toward their goal of delivering to the coming generation the gift of genome empowerment for many vertebrate species.

Divergent evolutionary rates in vertebrate and mammalian specific conserved non-coding elements (CNEs) in echolocating mammals

BACKGROUND

The majority of DNA contained within vertebrate genomes is non-coding, with a certain proportion of this thought to play regulatory roles during development. Conserved Non-coding Elements (CNEs) are an abundant group of putative regulatory sequences that are highly conserved across divergent groups and thus assumed to be under strong selective constraint. Many CNEs may contain regulatory factor binding sites, and their frequent spatial association with key developmental genes - such as those regulating sensory system development - suggests crucial roles in regulating gene expression and cellular patterning. Yet surprisingly little is known about the molecular evolution of CNEs across diverse mammalian taxa or their role in specific phenotypic adaptations. We examined 3,110 vertebrate-specific and ~82,000 mammalian-specific CNEs across 19 and 9 mammalian orders respectively, and tested for changes in the rate of evolution of CNEs located in the proximity of genes underlying the development or functioning of auditory systems. As we focused on CNEs putatively associated with genes underlying the development/functioning of auditory systems, we incorporated echolocating taxa in our dataset because of their highly specialised and derived auditory systems.

RESULTS

Phylogenetic reconstructions of concatenated CNEs broadly recovered accepted mammal relationships despite high levels of sequence conservation. We found that CNE substitution rates were highest in rodents and lowest in primates, consistent with previous findings. Comparisons of CNE substitution rates from several genomic regions containing genes linked to auditory system development and hearing revealed differences between echolocating and non-echolocating taxa. Wider taxonomic sampling of four CNEs associated with the homeobox genes Hmx2 and Hmx3 - which are required for inner ear development - revealed family-wise variation across diverse bat species. Specifically within one family of echolocating bats that utilise frequency-modulated echolocation calls varying widely in frequency and intensity high levels of sequence divergence were found.

CONCLUSIONS

Levels of selective constraint acting on CNEs differed both across genomic locations and taxa, with observed variation in substitution rates of CNEs among bat species. More work is needed to determine whether this variation can be linked to echolocation, and wider taxonomic sampling is necessary to fully document levels of conservation in CNEs across diverse taxa.

Comparative genomics reveals insights into avian genome evolution and adaptation

Birds are the most species-rich class of tetrapod vertebrates and have wide relevance across many research fields. We explored bird macroevolution using full genomes from 48 avian species representing all major extant clades. The avian genome is principally characterized by its constrained size, which predominantly arose because of lineage-specific erosion of repetitive elements, large segmental deletions, and gene loss. Avian genomes furthermore show a remarkably high degree of evolutionary stasis at the levels of nucleotide sequence, gene synteny, and chromosomal structure. Despite this pattern of conservation, we detected many non-neutral evolutionary changes in protein-coding genes and noncoding regions. These analyses reveal that pan-avian genomic diversity covaries with adaptations to different lifestyles and convergent evolution of traits.

Multiple lineages of ancient CR1 retroposons shaped the early genome evolution of amniotes

Chicken repeat 1 (CR1) retroposons are long interspersed elements (LINEs) that are ubiquitous within amniote genomes and constitute the most abundant family of transposed elements in birds, crocodilians, turtles, and snakes. They are also present in mammalian genomes, where they reside as numerous relics of ancient retroposition events. Yet, despite their relevance for understanding amniote genome evolution, the diversity and evolution of CR1 elements has never been studied on an amniote-wide level. We reconstruct the temporal and quantitative activity of CR1 subfamilies via presence/absence analyses across crocodilian phylogeny and comparative analyses of 12 crocodilian genomes, revealing relative genomic stasis of retroposition during genome evolution of extant Crocodylia. Our large-scale phylogenetic analysis of amniote CR1 subfamilies suggests the presence of at least seven ancient CR1 lineages in the amniote ancestor; and amniote-wide analyses of CR1 successions and quantities reveal differential retention (presence of ancient relics or recent activity) of these CR1 lineages across amniote genome evolution. Interestingly, birds and lepidosaurs retained the fewest ancient CR1 lineages among amniotes and also exhibit smaller genome sizes. Our study is the first to analyze CR1 evolution in a genome-wide and amniote-wide context and the data strongly suggest that the ancestral amniote genome contained myriad CR1 elements from multiple ancient lineages, and remnants of these are still detectable in the relatively stable genomes of crocodilians and turtles. Early mammalian genome evolution was thus characterized by a drastic shift from CR1 prevalence to dominance and hyperactivity of L2 LINEs in monotremes and L1 LINEs in therians.

Early mesozoic coexistence of amniotes and hepadnaviridae

Hepadnaviridae are double-stranded DNA viruses that infect some species of birds and mammals. This includes humans, where hepatitis B viruses (HBVs) are prevalent pathogens in considerable parts of the global population. Recently, endogenized sequences of HBVs (eHBVs) have been discovered in bird genomes where they constitute direct evidence for the coexistence of these viruses and their hosts from the late Mesozoic until present. Nevertheless, virtually nothing is known about the ancient host range of this virus family in other animals. Here we report the first eHBVs from crocodilian, snake, and turtle genomes, including a turtle eHBV that endogenized >207 million years ago. This genomic “fossil” is >125 million years older than the oldest avian eHBV and provides the first direct evidence that Hepadnaviridae already existed during the Early Mesozoic. This implies that the Mesozoic fossil record of HBV infection spans three of the five major groups of land vertebrates, namely birds, crocodilians, and turtles. We show that the deep phylogenetic relationships of HBVs are largely congruent with the deep phylogeny of their amniote hosts, which suggests an ancient amniote–HBV coexistence and codivergence, at least since the Early Mesozoic. Notably, the organization of overlapping genes as well as the structure of elements involved in viral replication has remained highly conserved among HBVs along that time span, except for the presence of the X gene. We provide multiple lines of evidence that the tumor-promoting X protein of mammalian HBVs lacks a homolog in all other hepadnaviruses and propose a novel scenario for the emergence of X via segmental duplication and overprinting of pre-existing reading frames in the ancestor of mammalian HBVs. Our study reveals an unforeseen host range of prehistoric HBVs and provides novel insights into the genome evolution of hepadnaviruses throughout their long-lasting association with amniote hosts.

Viruses are not known to leave physical fossil traces, which makes our understanding of their evolutionary prehistory crucially dependent on the detection of endogenous viruses. Ancient endogenous viruses, also known as paleoviruses, are relics of viral genomes or fragments thereof that once infiltrated their host's germline and then remained as molecular “fossils” within the host genome. The massive genome sequencing of recent years has unearthed vast numbers of paleoviruses from various animal genomes, including the first endogenous hepatitis B viruses (eHBVs) in bird genomes. We screened genomes of land vertebrates (amniotes) for the presence of paleoviruses and identified ancient eHBVs in the recently sequenced genomes of crocodilians, snakes, and turtles. We report an eHBV that is >207 million years old, making it the oldest endogenous virus currently known. Furthermore, our results provide direct evidence that the Hepadnaviridae virus family infected birds, crocodilians and turtles during the Mesozoic Era, and suggest a long-lasting coexistence of these viruses and their amniote hosts at least since the Early Mesozoic. We challenge previous views on the origin of the oncogenic X gene and provide an evolutionary explanation as to why only mammalian hepatitis B infection leads to hepatocellular carcinoma.

Molecular analysis of the koala reproductive hormones and their receptors: gonadotrophin-releasing hormone (GnRH), follicle-stimulating hormone β and luteinising hormone β with localisation of GnRH

During evolution, reproductive hormones and their receptors in the brain-pituitary-gonadal axis have been altered by genetic mechanisms. To understand how the neuroendocrine control of reproduction evolved in mammals, it is important to examine marsupials, the closest group to placental mammals. We hypothesised that at least some of the hormones and receptors found in placental mammals would be present in koala, a marsupial. We examined the expression of koala mRNA for the reproductive molecules. Koala cDNAs were cloned from brain for gonadotrophin-releasing hormones (GnRH1 and GnRH2) or from pituitary for GnRH receptors, types I and II, follicle-stimulating hormone (FSH)β and luteinising hormone (LH)β, and from gonads for FSH and LH receptors. Deduced proteins were compared by sequence alignment and phylogenetic analysis with those of other vertebrates. In conclusion, the koala expressed mRNA for these eight putative reproductive molecules, whereas at least one of these molecules is missing in some species in the amniote lineage, including humans. In addition, GnRH1 and 2 are shown by immunohistochemistry to be expressed as proteins in the brain.

Tempo and mode of gene duplication in mammalian ribosomal protein evolution

Gene duplication has been widely recognized as a major driver of evolutionary change and organismal complexity through the generation of multi-gene families. Therefore, understanding the forces that govern the evolution of gene families through the retention or loss of duplicated genes is fundamentally important in our efforts to study genome evolution. Previous work from our lab has shown that ribosomal protein (RP) genes constitute one of the largest classes of conserved duplicated genes in mammals. This result was surprising due to the fact that ribosomal protein genes evolve slowly and transcript levels are very tightly regulated. In our present study, we identified and characterized all RP duplicates in eight mammalian genomes in order to investigate the tempo and mode of ribosomal protein family evolution. We show that a sizable number of duplicates are transcriptionally active and are very highly conserved. Furthermore, we conclude that existing gene duplication models do not readily account for the preservation of a very large number of intact retroduplicated ribosomal protein (RT-RP) genes observed in mammalian genomes. We suggest that selection against dominant-negative mutations may underlie the unexpected retention and conservation of duplicated RP genes, and may shape the fate of newly duplicated genes, regardless of duplication mechanism.

A transcriptome resource for the koala (Phascolarctos cinereus): insights into koala retrovirus transcription and sequence diversity

Background

The koala, Phascolarctos cinereus, is a biologically unique and evolutionarily distinct Australian arboreal marsupial. The goal of this study was to sequence the transcriptome from several tissues of two geographically separate koalas, and to create the first comprehensive catalog of annotated transcripts for this species, enabling detailed analysis of the unique attributes of this threatened native marsupial, including infection by the koala retrovirus.

Results

RNA-Seq data was generated from a range of tissues from one male and one female koala and assembled de novo into transcripts using Velvet-Oases. Transcript abundance in each tissue was estimated. Transcripts were searched for likely protein-coding regions and a non-redundant set of 117,563 putative protein sequences was produced. In similarity searches there were 84,907 (72%) sequences that aligned to at least one sequence in the NCBI nr protein database. The best alignments were to sequences from other marsupials. After applying a reciprocal best hit requirement of koala sequences to those from tammar wallaby, Tasmanian devil and the gray short-tailed opossum, we estimate that our transcriptome dataset represents approximately 15,000 koala genes. The marsupial alignment information was used to look for potential gene duplications and we report evidence for copy number expansion of the alpha amylase gene, and of an aldehyde reductase gene.

Koala retrovirus (KoRV) transcripts were detected in the transcriptomes. These were analysed in detail and the structure of the spliced envelope gene transcript was determined. There was appreciable sequence diversity within KoRV, with 233 sites in the KoRV genome showing small insertions/deletions or single nucleotide polymorphisms. Both koalas had sequences from the KoRV-A subtype, but the male koala transcriptome has, in addition, sequences more closely related to the KoRV-B subtype. This is the first report of a KoRV-B-like sequence in a wild population.

Conclusions

This transcriptomic dataset is a useful resource for molecular genetic studies of the koala, for evolutionary genetic studies of marsupials, for validation and annotation of the koala genome sequence, and for investigation of koala retrovirus. Annotated transcripts can be browsed and queried at http://koalagenome.org.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-786) contains supplementary material, which is available to authorized users.

Relating gene expression evolution with CpG content changes

Background

Previous studies have shown that CpG dinucleotides are enriched in a subset of promoters and the CpG content of promoters is positively correlated with gene expression levels. But the relationship between divergence of CpG content and gene expression evolution has not been investigated. Here we calculate the normalized CpG (nCpG) content in DNA regions around transcription start site (TSS) and transcription terminal site (TTS) of genes in nine organisms, and relate them with expression levels measured by RNA-seq.

Results

The nCpG content of TSS shows a bimodal distribution in all organisms except platypus, whereas the nCpG content of TTS only has a single peak. When the nCpG contents are compared between different organisms, we observe a different evolution pattern between TSS and TTS: compared with TTS, TSS exhibits a faster divergence rate between closely related species but are more conserved between distant species. More importantly, we demonstrate the link between gene expression evolution and nCpG content changes: up-/down- regulation of genes in an organism is accompanied by the nCpG content increase/decrease in their TSS and TTS proximal regions.

Conclusions

Our results suggest that gene expression changes between different organisms are correlated with the alterations in normalized CpG contents of promoters. Our analyses provide evidences for the impact of nCpG content on gene expression evolution.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-693) contains supplementary material, which is available to authorized users.

Frequent loss and alteration of the MOXD2 gene in catarrhines and whales: a possible connection with the evolution of olfaction

The MOXD2 gene encodes a membrane-bound monooxygenase similar to dopamine-β-hydroxylase, and has been proposed to be associated with olfaction. In this study, we analyzed MOXD2 genes from 64 mammalian species, and identified loss-of-function mutations in apes (humans, Sumatran and Bornean orangutans, and five gibbon species from the four major gibbon genera), toothed whales (killer whales, bottlenose dolphins, finless porpoises, baijis, and sperm whales), and baleen whales (minke whales and fin whales). We also identified a shared 13-nt deletion in the last exon of Old World cercopithecine monkeys that results in conversion of a membrane-bound protein to a soluble form. We hypothesize that the frequent inactivation and alteration of MOXD2 genes in catarrhines and whales may be associated with the evolution of olfaction in these clades.

Noncoding origins of anthropoid traits and a new null model of transposon functionalization

Little is known about novel genetic elements that drove the emergence of anthropoid primates. We exploited the sequencing of the marmoset genome to identify 23,849 anthropoid-specific constrained (ASC) regions and confirmed their robust functional signatures. Of the ASC base pairs, 99.7% were noncoding, suggesting that novel anthropoid functional elements were overwhelmingly cis-regulatory. ASCs were highly enriched in loci associated with fetal brain development, motor coordination, neurotransmission, and vision, thus providing a large set of candidate elements for exploring the molecular basis of hallmark primate traits. We validated ASC192 as a primate-specific enhancer in proliferative zones of the developing brain. Unexpectedly, transposable elements (TEs) contributed to >56% of ASCs, and almost all TE families showed functional potential similar to that of nonrepetitive DNA. Three L1PA repeat-derived ASCs displayed coherent eye-enhancer function, thus demonstrating that the "gene-battery" model of TE functionalization applies to enhancers in vivo. Our study provides fundamental insights into genome evolution and the origins of anthropoid phenotypes and supports an elegantly simple new null model of TE exaptation.

Noncoding RNAs and epigenetic mechanisms during X-chromosome inactivation

In mammals, the process of X-chromosome inactivation ensures equivalent levels of X-linked gene expression between males and females through the silencing of one of the two X chromosomes in female cells. The process is established early in development and is initiated by a unique locus, which produces a long noncoding RNA, Xist. The Xist transcript triggers gene silencing in cis by coating the future inactive X chromosome. It also induces a cascade of chromatin changes, including posttranslational histone modifications and DNA methylation, and leads to the stable repression of all X-linked genes throughout development and adult life. We review here recent progress in our understanding of the molecular mechanisms involved in the initiation of Xist expression, the propagation of the Xist RNA along the chromosome, and the cis-elements and trans-acting factors involved in the maintenance of the repressed state. We also describe the diverse strategies used by nonplacental mammals for X-chromosome dosage compensation and highlight the common features and differences between eutherians and metatherians, in particular regarding the involvement of long noncoding RNAs.

Age-dependent transcriptome and proteome following transection of neonatal spinal cord of Monodelphis domestica (South American grey short-tailed opossum)

This study describes a combined transcriptome and proteome analysis of Monodelphis domestica response to spinal cord injury at two different postnatal ages. Previously we showed that complete transection at postnatal day 7 (P7) is followed by profuse axon growth across the lesion with near-normal locomotion and swimming when adult. In contrast, at P28 there is no axon growth across the lesion, the animals exhibit weight-bearing locomotion, but cannot use hind limbs when swimming. Here we examined changes in gene and protein expression in the segment of spinal cord rostral to the lesion at 24 h after transection at P7 and at P28. Following injury at P7 only forty genes changed (all increased expression); most were immune/inflammatory genes. Following injury at P28 many more genes changed their expression and the magnitude of change for some genes was strikingly greater. Again many were associated with the immune/inflammation response. In functional groups known to be inhibitory to regeneration in adult cords the expression changes were generally muted, in some cases opposite to that required to account for neurite inhibition. For example myelin basic protein expression was reduced following injury at P28 both at the gene and protein levels. Only four genes from families with extracellular matrix functions thought to influence neurite outgrowth in adult injured cords showed substantial changes in expression following injury at P28: Olfactomedin 4 (Olfm4, 480 fold compared to controls), matrix metallopeptidase (Mmp1, 104 fold), papilin (Papln, 152 fold) and integrin α4 (Itga4, 57 fold). These data provide a resource for investigation of a priori hypotheses in future studies of mechanisms of spinal cord regeneration in immature animals compared to lack of regeneration at more mature stages.

Heterochromatin dynamics during the differentiation process revealed by the DNA methylation reporter mouse, MethylRO

In mammals, DNA is methylated at CpG sites, which play pivotal roles in gene silencing and chromatin organization. Furthermore, DNA methylation undergoes dynamic changes during development, differentiation, and in pathological processes. The conventional methods represent snapshots; therefore, the dynamics of this marker within living organisms remains unclear. To track this dynamics, we made a knockin mouse that expresses a red fluorescent protein (RFP)-fused methyl-CpG-binding domain (MBD) protein from the ROSA26 locus ubiquitously; we named it MethylRO (methylation probe in ROSA26 locus). Using this mouse, we performed RFP-mediated methylated DNA immunoprecipitation sequencing (MeDIP-seq), whole-body section analysis, and live-cell imaging. We discovered that mobility and pattern of heterochromatin as well as DNA methylation signal intensity inside the nuclei can be markers for cellular differentiation status. Thus, the MethylRO mouse represents a powerful bioresource and technique for DNA methylation dynamics studies in developmental biology, stem cell biology, as well as in disease states.

•

Changes in DNA methylation are tracked in living mice

•

Heterochromatin structure changes dynamically during development and differentiation

•

Heterochromatin of preimplantation embryonic cells is highly dynamic than ESCs

•

Heterochromatin pattern in nucleus can be a marker for cell differentiation states

Tracing the evolution of amniote chromosomes

A great deal of diversity in chromosome number and arrangement is observed across the amniote phylogeny. Understanding how this diversity is generated is important for determining the role of chromosomal rearrangements in generating phenotypic variation and speciation. Gaining this understanding is achieved by reconstructing the ancestral genome arrangement based on comparisons of genome organization of extant species. Ancestral karyotypes for several amniote lineages have been reconstructed, mainly from cross-species chromosome painting data. The availability of anchored whole genome sequences for amniote species has increased the evolutionary depth and confidence of ancestral reconstructions from those made solely from chromosome painting data. Nonetheless, there are still several key lineages where the appropriate data required for ancestral reconstructions is lacking. This review highlights the progress that has been made towards understanding the chromosomal changes that have occurred during amniote evolution and the reconstruction of ancestral karyotypes.

Analysis of horse genomes provides insight into the diversification and adaptive evolution of karyotype

Karyotypic diversification is more prominent in Equus species than in other mammals. Here, using next generation sequencing technology, we generated and de novo assembled quality genomes sequences for a male wild horse (Przewalski's horse) and a male domestic horse (Mongolian horse), with about 93-fold and 91-fold coverage, respectively. Portion of Y chromosome from wild horse assemblies (3 M bp) and Mongolian horse (2 M bp) were also sequenced and de novo assembled. We confirmed a Robertsonian translocation event through the wild horse's chromosomes 23 and 24, which contained sequences that were highly homologous with those on the domestic horse's chromosome 5. The four main types of rearrangement, insertion of unknown origin, inserted duplication, inversion, and relocation, are not evenly distributed on all the chromosomes, and some chromosomes, such as the X chromosome, contain more rearrangements than others, and the number of inversions is far less than the number of insertions and relocations in the horse genome. Furthermore, we discovered the percentages of LINE_L1 and LTR_ERV1 are significantly increased in rearrangement regions. The analysis results of the two representative Equus species genomes improved our knowledge of Equus chromosome rearrangement and karyotype evolution.

Drosophila 3' UTRs are more complex than protein-coding sequences

The 3′ UTRs of eukaryotic genes participate in a variety of post-transcriptional (and some transcriptional) regulatory interactions. Some of these interactions are well characterised, but an undetermined number remain to be discovered. While some regulatory sequences in 3′ UTRs may be conserved over long evolutionary time scales, others may have only ephemeral functional significance as regulatory profiles respond to changing selective pressures. Here we propose a sensitive segmentation methodology for investigating patterns of composition and conservation in 3′ UTRs based on comparison of closely related species. We describe encodings of pairwise and three-way alignments integrating information about conservation, GC content and transition/transversion ratios and apply the method to three closely related Drosophila species: D. melanogaster, D. simulans and D. yakuba. Incorporating multiple data types greatly increased the number of segment classes identified compared to similar methods based on conservation or GC content alone. We propose that the number of segments and number of types of segment identified by the method can be used as proxies for functional complexity. Our main finding is that the number of segments and segment classes identified in 3′ UTRs is greater than in the same length of protein-coding sequence, suggesting greater functional complexity in 3′ UTRs. There is thus a need for sustained and extensive efforts by bioinformaticians to delineate functional elements in this important genomic fraction. C code, data and results are available upon request.

Realistic artificial DNA sequences as negative controls for computational genomics

A common practice in computational genomic analysis is to use a set of ‘background’ sequences as negative controls for evaluating the false-positive rates of prediction tools, such as gene identification programs and algorithms for detection of cis-regulatory elements. Such ‘background’ sequences are generally taken from regions of the genome presumed to be intergenic, or generated synthetically by ‘shuffling’ real sequences. This last method can lead to underestimation of false-positive rates. We developed a new method for generating artificial sequences that are modeled after real intergenic sequences in terms of composition, complexity and interspersed repeat content. These artificial sequences can serve as an inexhaustible source of high-quality negative controls. We used artificial sequences to evaluate the false-positive rates of a set of programs for detecting interspersed repeats, ab initio prediction of coding genes, transcribed regions and non-coding genes. We found that RepeatMasker is more accurate than PClouds, Augustus has the lowest false-positive rate of the coding gene prediction programs tested, and Infernal has a low false-positive rate for non-coding gene detection. A web service, source code and the models for human and many other species are freely available at http://repeatmasker.org/garlic/.

Conserved noncoding elements follow power-law-like distributions in several genomes as a result of genome dynamics

Conserved, ultraconserved and other classes of constrained elements (collectively referred as CNEs here), identified by comparative genomics in a wide variety of genomes, are non-randomly distributed across chromosomes. These elements are defined using various degrees of conservation between organisms and several thresholds of minimal length. We here investigate the chromosomal distribution of CNEs by studying the statistical properties of distances between consecutive CNEs. We find widespread power-law-like distributions, i.e. linearity in double logarithmic scale, in the inter-CNE distances, a feature which is connected with fractality and self-similarity. Given that CNEs are often found to be spatially associated with genes, especially with those that regulate developmental processes, we verify by appropriate gene masking that a power-law-like pattern emerges irrespectively of whether elements found close or inside genes are excluded or not. An evolutionary model is put forward for the understanding of these findings that includes segmental or whole genome duplication events and eliminations (loss) of most of the duplicated CNEs. Simulations reproduce the main features of the observed size distributions. Power-law-like patterns in the genomic distributions of CNEs are in accordance with current knowledge about their evolutionary history in several genomes.

Why do a wide variety of animals retain multiple isoforms of cyclooxygenase?

Cyclooxygenase (COX) has been cloned from the phyla Cnidaria, Mollusca, Arthropoda, and Chordata of the animal kingdom. Many organisms have multiple COX isoforms that have arisen from gene duplication. It is not well understood why there are multiple COX isoforms in the same organism, or when duplication of the COX gene occurred. Here, we summarize the current knowledge of the evolutionary history of COX in the animal kingdom and discuss the reasons why the multiple COX system has been retained so widely. The phylogenetic analysis suggests that all COX genes in animals may descend from a common ancestor and that the duplication of an ancestral COX gene might occur within each lineage after the divergence of the animal. In most instances, the expressions of multiple COX isoforms are separately regulated and these isoforms play different and important pathophysiological roles in each organism. This may be the reason why multiple COX isoforms are widely retained.

Photic preference of the short-tailed opossum (Monodelphis domestica)

The gray short-tailed opossum (Monodelphis domestica) is a nocturnal South American marsupial that has been gaining popularity as a laboratory animal. However, compared to traditional laboratory animals like rats, very little is known about its behavior, either in the wild or in a laboratory setting. Here we investigated the photic preference of the short-tailed opossum. Opossums were placed in a circular testing arena and allowed to move freely between dark (0 lux) and light (∼1.4, 40, or 400 lux) sides of the arena. In each of these conditions opossums spent significantly more time in the dark than in the illuminated side and a greater proportion of time in the dark than would be expected by chance. In the high-contrast (∼400 lux) illumination condition, the mean bout length (i.e., duration of one trip on the light or dark side) was significantly longer on the dark side than on the light side. When we examined the number of bouts greater than 30 and 60s in duration, we found a significant difference between the light and dark sides in all light contrast conditions. These data indicate that the short-tailed opossum prefers the dark to the light, and can also detect very slight differences in light intensity. We conclude that although rats and opossums share many similar characteristics, including ecological niche, their divergent evolutionary heritage results in vastly different behavioral capabilities. Only by observing the behavioral capabilities and preferences of opossums will we be able to manipulate the experimental environment to best elicit and elucidate their behavior and alterations in behavior that can arise from experimental manipulations.

The phylogenetic distribution and evolution of enzymes within the thymidine kinase 2-like gene family in metazoa

Deoxyribonucleoside kinases (dNKs) carry out the rate-determining step in the nucleoside salvage pathway within all domains of life where the pathway is present, and, hence, are an indication on whether or not a species/genus retains the ability to salvage deoxyribonucleosides. Here, a phylogenetic tree is constructed for the thymidine kinase 2-like dNK gene family in metazoa. Each enzyme class (deoxycytidine, deoxyguanosine, and deoxythymidine kinases, as well as the multisubstrate dNKs) falls into a monophyletic clade. However, in vertebrates, dCK contains an apparent duplication with one paralog lost in mammals, and a number of crustacean genomes (like Caligus rogercresseyi and Lepeophtheirus salmonis) unexpectedly contain not only the multisubstrate dNKs, related to Drosophila multisubstrate dNK, but also a TK2-like kinase. Additionally, crustaceans (Daphnia, Caligus, and Lepeophtheirus) and some insects (Tribolium, Danaus, Pediculus, and Acyrthosiphon) contain several multisubstrate dNK-like enzymes which group paraphyletically within the arthropod clade. This might suggest that the multisubstrate dNKs underwent multiple rounds of duplications with differential retention of duplicate copies between insect families and more complete retention within some crustaceans and insects. Genomes of several basal animalia contain more than one dNK-like sequence, some of which group outside the remaining eukaryotes (both plants and animals) and/or with bacterial dNKs. Within the vertebrates, the mammalian genomes do not contain the second dCK, while birds, fish, and amphibians do retain it. Phasianidae (chicken and turkey) have lost dGK, while it has been retained in other bird lineages, like zebra finch. Reconstruction of the ancestral sequence between the multisubstrate arthropod dNKs and the TK2 clade of vertebrates followed by homology modeling and discrete molecular dynamics calculations on this sequence were performed to examine the evolutionary path which led to the two different enzyme classes. The structural models showed that the carboxyl terminus of the ancestral sequence is more helical than dNK, in common with TK2, although any implications of this for enzyme specificity will require biochemical validation. Finally, rate-shift and conservation-shift analysis between clades with different specificities uncovered candidate residues outside the active site pocket which may have contributed to differentiation in substrate specificity between enzyme clades.

Evolutionarily dynamic alternative splicing of GPR56 regulates regional cerebral cortical patterning

The human neocortex has numerous specialized functional areas whose formation is poorly understood. Here, we describe a 15-base pair deletion mutation in a regulatory element of GPR56 that selectively disrupts human cortex surrounding the Sylvian fissure bilaterally including "Broca's area," the primary language area, by disrupting regional GPR56 expression and blocking RFX transcription factor binding. GPR56 encodes a heterotrimeric guanine nucleotide-binding protein (G protein)-coupled receptor required for normal cortical development and is expressed in cortical progenitor cells. GPR56 expression levels regulate progenitor proliferation. GPR56 splice forms are highly variable between mice and humans, and the regulatory element of gyrencephalic mammals directs restricted lateral cortical expression. Our data reveal a mechanism by which control of GPR56 expression pattern by multiple alternative promoters can influence stem cell proliferation, gyral patterning, and, potentially, neocortex evolution.

Genome-wide histone state profiling of fibroblasts from the opossum, Monodelphis domestica, identifies the first marsupial-specific imprinted gene

Background

Imprinted genes have been extensively documented in eutherian mammals and found to exhibit significant interspecific variation in the suites of genes that are imprinted and in their regulation between tissues and developmental stages. Much less is known about imprinted loci in metatherian (marsupial) mammals, wherein studies have been limited to a small number of genes previously known to be imprinted in eutherians. We describe the first ab initio search for imprinted marsupial genes, in fibroblasts from the opossum, Monodelphis domestica, based on a genome-wide ChIP-seq strategy to identify promoters that are simultaneously marked by mutually exclusive, transcriptionally opposing histone modifications.

Results

We identified a novel imprinted gene (Meis1) and two additional monoallelically expressed genes, one of which (Cstb) showed allele-specific, but non-imprinted expression. Imprinted vs. allele-specific expression could not be resolved for the third monoallelically expressed gene (Rpl17). Transcriptionally opposing histone modifications H3K4me3, H3K9Ac, and H3K9me3 were found at the promoters of all three genes, but differential DNA methylation was not detected at CpG islands at any of these promoters.

Conclusions

In generating the first genome-wide histone modification profiles for a marsupial, we identified the first gene that is imprinted in a marsupial but not in eutherian mammals. This outcome demonstrates the practicality of an ab initio discovery strategy and implicates histone modification, but not differential DNA methylation, as a conserved mechanism for marking imprinted genes in all therian mammals. Our findings suggest that marsupials use multiple epigenetic mechanisms for imprinting and support the concept that lineage-specific selective forces can produce sets of imprinted genes that differ between metatherian and eutherian lines.

Recognition of vitamin B metabolites by mucosal-associated invariant T cells

The mucosal-associated invariant T-cell antigen receptor (MAIT TCR) recognizes MR1 presenting vitamin B metabolites. Here we describe the structures of a human MAIT TCR in complex with human MR1 presenting a non-stimulatory ligand derived from folic acid and an agonist ligand derived from a riboflavin metabolite. For both vitamin B antigens, the MAIT TCR docks in a conserved manner above MR1, thus acting as an innate-like pattern recognition receptor. The invariant MAIT TCR α-chain usage is attributable to MR1-mediated interactions that prise open the MR1 cleft to allow contact with the vitamin B metabolite. Although the non-stimulatory antigen does not contact the MAIT TCR, the stimulatory antigen does. This results in a higher affinity of the MAIT TCR for a stimulatory antigen in comparison with a non-stimulatory antigen. We formally demonstrate a structural basis for MAIT TCR recognition of vitamin B metabolites, while illuminating how TCRs recognize microbial metabolic signatures.