Search

The chromatin environment plays a central role in regulating developmental gene expression in metazoans. Yet, the basal regulatory landscape of metazoan embryogenesis is unknown. Here, we generate chromatin accessibility profiles for six embryonic, plus larval and adult stages in the sponge Amphimedon queenslandica. These profiles are reproducible within stages, reflect histone modifications, and identify transcription factor (TF) binding sequence motifs predictive of cis-regulatory elements during embryogenesis in other metazoans, but not the unicellular relative Capsaspora. Motif analysis of chromatin accessibility profiles across Amphimedon embryogenesis identifies three major developmental periods. As in bilaterian embryogenesis, early development in Amphimedon involves activating and repressive chromatin in regions both proximal and distal to transcription start sites. Transcriptionally repressive elements (‘silencers’) are prominent during late embryogenesis. They coincide with an increase in cis-regulatory regions harbouring metazoan TF binding motifs, and an increase in the expression of metazoan-specific genes. Changes in chromatin state and gene expression in Amphimedon suggest the conservation of distal enhancers, dynamically silenced chromatin, and TF-DNA binding specificity in animal embryogenesis.

Precise patterns of gene expression are driven by interactions between transcription factors, regulatory DNA sequences, and chromatin. How DNA mutations affecting any one of these regulatory “layers” are buffered or propagated to gene expression remains unclear. To address this, we quantified allele-specific changes in chromatin accessibility, histone modifications, and gene expression in F1 embryos generated from eight Drosophila crosses at three embryonic stages, yielding a comprehensive data set of 240 samples spanning multiple regulatory layers. Genetic variation (allelic imbalance) impacts gene expression more frequently than chromatin features, with metabolic and environmental response genes being most often affected. Allelic imbalance in cis-regulatory elements (enhancers) is common and highly heritable, yet its functional impact does not generally propagate to gene expression. When it does, genetic variation impacts RNA levels through two alternative mechanisms involving either H3K4me3 or chromatin accessibility and H3K27ac. Changes in RNA are more predictive of variation in H3K4me3 than vice versa, suggesting a role for H3K4me3 downstream from transcription. The impact of a substantial proportion of genetic variation is consistent across embryonic stages, with 50% of allelic imbalanced features at one stage being also imbalanced at subsequent developmental stages. Crucially, buffering, as well as the magnitude and evolutionary impact of genetic variants, is influenced by regulatory complexity (i.e., number of enhancers regulating a gene), with transcription factors being most robust to cis-acting, but most influenced by trans-acting, variation.

Precise patterns of gene expression are driven by interactions between transcription factors, regulatory DNA sequence, and chromatin. How DNA mutations affecting any one of these regulatory ‘layers’ is buffered or propagated to gene expression remains unclear. To address this, we quantified allele-specific changes in chromatin accessibility, histone modifications, and gene expression in F1 embryos generated from eightDrosophilacrosses, at three embryonic stages, yielding a comprehensive dataset of 240 samples spanning multiple regulatory layers. Genetic variation incis-regulatory elements is common, highly heritable, and surprisingly consistent in its effects across embryonic stages. Much of this variation does not propagate to gene expression. When it does, it acts through H3K4me3 or alternatively through chromatin accessibility and H3K27ac. The magnitude and evolutionary impact of mutations is influenced by a genes’ regulatory complexity (i.e. enhancer number), with transcription factors being most robust tocis-acting, and most influenced bytrans-acting, variation. Overall, the impact of genetic variation on regulatory phenotypes appears context-dependent even within the constraints of embryogenesis.

Narirutin is one of the most common flavanones found in citrus fruits. The vascular effects of its analogues naringenin and naringin have been reported but its effects on the cardiovascular system are largely unknown. In this study, relaxation effect of narirutin and its mechanisms of action were investigated by measuring isometric tension in rat mesenteric arteries. Patch-clamping was also used to study the effect of narirutin on potassium channels in vascular smooth muscle cells. Moreover, its effects on phosphorylation of endothelial nitric oxide synthase, cAMP level and phosphodiesterase activity in rat mesenteric arteries were studied by Western blot and biochemical assays. The results showed that pre-incubation of rat mesenteric arteries with narirutin had no influence on acetylcholine-induced endothelial-dependent relaxation. However, narirutin caused a direct concentration-dependent relaxation in rat mesenteric arteries. This relaxation effect was comparable to that of narirutin's structural analogue naringenin. Narirutin-induced relaxation was reduced by the removal of endothelium, NG-nitro-L-arginine methyl ester (a nitric oxide synthase inhibitor), and 4-aminopyridine (a voltage-gated potassium channel blocker). In addition, narirutin increased the phosphorylation of endothelial nitric oxide synthase and increased the voltage-dependent potassium current in mesenteric arterial smooth muscle cells. These effects were abolished by protein kinase A inhibitor. Furthermore, narirutin could increase cAMP level and inhibit phosphodiesterase activity in rat mesenteric arteries. In conclusion, narirutin has vasorelaxing effect and the mechanism involves the inhibition of phosphodiesterase, which increases intracellular cAMP, thereby stimulating the endothelial nitric oxide synthase and activating the voltage-gated potassium channels in vascular smooth muscle cells.

By analyzing multi-tissue gene expression and genome-wide genetic variation data in samples from a vervet monkey pedigree, we generated a transcriptome resource and produced the first catalogue of expression quantitative trait loci (eQTLs) in a non-human primate model. This catalogue contains more genome-wide significant eQTLs, per sample, than comparable human resources, and reveals sex and age-related expression patterns. Findings include a master regulatory locus that likely plays a role in immune function, and a locus regulating hippocampal long non-coding RNAs (lncRNAs), whose expression correlates with hippocampal volume. This resource will facilitate genetic investigation of quantitative traits, including brain and behavioral phenotypes relevant to neuropsychiatric disorders.

Transcription factors (TFs) bind DNA enhancer sequences to regulate gene transcription in animals. Unlike TFs, the evolution of enhancers has been difficult to trace because of their rapid evolution. Here, we show enhancers from the sponge Amphimedon queenslandica can drive cell type-specific reporter gene expression in zebrafish and mouse, despite sponge and vertebrate lineages diverging over 700 million years ago. Although sponge enhancers, which are present in both highly conserved syntenic gene regions (Islet–Scaper, Ccne1–Uri and Tdrd3–Diaph3) and sponge-specific intergenic regions, have no significant sequence identity with vertebrate genomic sequences, the type and frequency of TF binding motifs in the sponge enhancer allow for the identification of homologous enhancers in bilaterians. Islet enhancers identified in human and mouse Scaper genes drive zebrafish reporter expression patterns that are almost identical to the sponge Islet enhancer. The existence of homologous enhancers in these disparate metazoans suggests animal development is controlled by TF-enhancer DNA interactions that were present in the first multicellular animals.One-sentence summaryEnhancer activity is conserved across 700 million years of trans-phyletic divergence.

Enhancer function, from sponges to humans Identifying the function of enhancers, DNA regions that help to regulate gene expression and evolve rapidly, has been difficult. This area of research has been hampered by the difficultly in identifying functional conservation. Wong et al. now show that despite low sequence conservation, enhancer function is strongly conserved through the animal kingdom (see the Perspective by Harmston). Transgenic expression of sponge enhancers in zebrafish and mice demonstrates that these sequences can drive cell type–specific gene expression across species. These results suggest an unexpectedly deep level of conservation of gene regulation across the animal kingdom maintained over the course of metazoan evolution. Science , this issue p. eaax8137 ; see also p. 657

Differentiation into diverse cell lineages requires orchestration of gene regulatory networks guiding cell fate choices. Here, we present the dissection of cellular composition and gene networks from transcriptomic data of 43,168 cells across five discrete time points during cardiacdirected differentiation. We utilize unsupervised clustering and implement a lineage trajectory prediction algorithm that integrates transcription factor networks to predict cell fate progression of 15 subpopulations that correlate with germ layer and cardiovascular differentiation in vivo. These data reveal transcriptional networks underlying lineage derivation of mesoderm, definitive endoderm, vascular endothelium, cardiac precursors, and definitive cell types that comprise cardiomyocytes and a previously unrecognized cardiac outflow tract population. Single cell analysis of genetic regulators governing cardiac fate diversification identified the non-DNA binding homeodomain protein, HOPX, as functionally necessary for endothelial specification. Our findings further implicate dysregulation of HOPX during in vitro cardiac-directed differentiation underlying the molecular and physiological immaturity of stem cell-derived cardiomyocytes.

Background The α2-adrenergic sedative/anesthetic agent dexmedetomidine exerts biphasic effects on isolated arteries, causing endothelium-dependent relaxations at concentrations at or below 30 nM, followed by contractions at higher concentrations. l-arginine is a common substrate of endothelial nitric oxide synthase and arginases. This study was designed to investigate the role of l-arginine in modulating the overall vascular response to dexmedetomidine. Methods Isometric tension was measured in isolated aortic rings of Sprague Dawley rats. Cumulative concentrations of dexmedetomidine (10 nM to 10 μM) were added to quiescent rings (with and without endothelium) after previous incubation with vehicle, Nω-nitro-l-arginine methyl ester hydrochloride (l-NAME; nitric oxide synthase inhibitor), prazosin (α1-adrenergic antagonist), rauwolscine (α2-adrenergic antagonist), l-arginine, (S)-(2-boronethyl)-l-cysteine hydrochloride (arginase inhibitor), NG-hydroxy-l-arginine (arginase inhibitor), urea and/or ornithine. In some preparations, immunofluorescent staining, immunoblotting, or measurement of urea content were performed. Results Dexmedetomidine did not contract control rings with endothelium but evoked concentration-dependent increases in tension in such rings treated with l-NAME (Emax 50 ± 4%) or after endothelium-removal (Emax 74 ± 5%; N = 7 to 12). Exogenous l-arginine augmented the dexmedetomidine-induced contractions in the presence of l-NAME (Emax 75 ± 3%). This potentiation was abolished by (S)-(2-boronethyl)-l-cysteine hydrochloride (Emax 16 ± 4%) and NG-hydroxy-l-arginine (Emax 18 ± 4%). Either urea or ornithine, the downstream arginase products, had a similar potentiating effect as l-arginine. Immunoassay measurements demonstrated an upregulation of arginase I by l-arginine treatment in the presence of l-NAME (N = 4). Conclusions These results suggest that when vascular nitric oxide homeostasis is impaired, the potentiation of the vasoconstrictor effect of dexmedetomidine by l-arginine depends on arginase activity and the production of urea and ornithine.

Cardiac differentiation of human pluripotent stem cells (hPSCs) requires orchestration of dynamic gene regulatory networks during stepwise fate transitions, but often generates immature cell types that do not fully recapitulate properties of their adult counterparts, suggesting incomplete activation of key transcriptional networks. We performed extensive single-cell transcriptomic analyses to map fate choices and gene expression programs during cardiac differentiation of hPSCs, and identified strategies to improve in vitro cardiomyocyte differentiation. Utilizing genetic gain- and loss-of-function approaches, we found that hypertrophic signaling is not effectively activated during monolayer-based cardiac differentiation, thereby preventing expression of HOPX and its activation of downstream genes that govern late stages of cardiomyocyte maturation. This study therefore provides a key transcriptional roadmap of in vitro cardiac differentiation at single-cell resolution, revealing fundamental mechanisms underlying heart development and differentiation of hPSC-derived cardiomyocytes.

How genetic variation contributes to phenotypic variation is a central question in genetics. Association signals for a complex trait are found throughout the majority of the genome suggesting much of the genome is under some degree of genetic constraint. Here, we develop a intraspecific population genetics approach to define a measure of population structure for each single nucleotide polymorphism (SNP). Using this approach, we test for evidence of stabilizing selection at complex traits and pleiotropic loci arising from the evolutionary history of 47 complex traits and common diseases. Our approach allowed us to identify traits and regions under stabilizing selection towards both global and subpopulation optima. Strongest depletion of allelic diversity was found at disease loci, indicating stabilizing selection has acted on these phenotypes in all subpopulations. Pleiotropic loci predominantly displayed evidence of stabilizing selection, often contributed to multiple disease risks, and sometimes also affected non-disease traits such as height. Risk alleles at pleiotropic disease loci displayed a more consistent direction of effect than expected by chance suggesting that stabilizing selection acting on pleiotropic loci is amplified through multiple disease phenotypes.

To understand the evolutionary dynamics between transcription factor (TF) binding and gene expression in mammals, we compared transcriptional output and the binding intensities for three tissue-specific TFs in livers from four closely related mouse species. For each transcription factor, TF-dependent genes and the TF binding sites most likely to influence mRNA expression were identified by comparing mRNA expression levels between wild-type and TF knockout mice. Independent evolution was observed genome-wide between the rate of change in TF binding and the rate of change in mRNA expression across taxa, with the exception of a small number of TF-dependent genes. We also found that binding intensities are preferentially conserved near genes whose expression is dependent on the TF, and the conservation is shared among binding peaks in close proximity to each other near the TSS. Expression of TF-dependent genes typically showed an increased sensitivity to changes in binding levels as measured by mRNA abundance. Taken together, these results highlight a significant tolerance to evolutionary changes in TF binding intensity in mammalian transcriptional networks and suggest that some TF-dependent genes may be largely regulated by a single TF across evolution.

Noncoding regulatory variants play a central role in the genetics of human diseases and in evolution. Here we measure allele-specific transcription factor binding occupancy of three liver-specific transcription factors between crosses of two inbred mouse strains to elucidate the regulatory mechanisms underlying transcription factor binding variations in mammals. Our results highlight the pre-eminence of cis-acting variants on transcription factor occupancy divergence. Transcription factor binding differences linked to cis-acting variants generally exhibit additive inheritance, while those linked to trans-acting variants are most often dominantly inherited. Cis-acting variants lead to local coordination of transcription factor occupancies that decay with distance; distal coordination is also observed and may be modulated by long-range chromatin contacts. Our results reveal the regulatory mechanisms that interplay to drive transcription factor occupancy, chromatin state, and gene expression in complex mammalian cell states., “Variation in the noncoding regulatory sequences in the genome plays important roles in human disease and evolution. Here, the authors use F1 mouse hybrids to shed light on the regulatory mechanisms mediating transcription factor binding, chromatin state and gene expression in mammalian cells.”

The platypus is a venomous monotreme. Male platypuses possess a spur on their hind legs that is connected to glands in the pelvic region. They produce venom only during the breeding season, presumably to fight off conspecifics. We have taken advantage of this unique seasonal production of venom to compare the transcriptomes of in- and out-of-season venom glands, in conjunction with proteomic analysis, to identify previously undiscovered venom genes. Comparison of the venom glands revealed distinct gene expression profiles that are consistent with changes in venom gland morphology and venom volumes in and out of the breeding season. Venom proteins were identified through shot-gun sequenced venom proteomes of three animals using RNA-seq-derived transcripts for peptide-spectral matching. 5,157 genes were expressed in the venom glands, 1,821 genes were up-regulated in the in-season gland, and 10 proteins were identified in the venom. New classes of platypus-venom proteins identified included antimicrobials, amide oxidase, serpin protease inhibitor, proteins associated with the mammalian stress response pathway, cytokines, and other immune molecules. Five putative toxins have only been identified in platypus venom: growth differentiation factor 15, nucleobindin-2, CD55, a CXC-chemokine, and corticotropin-releasing factor-binding protein. These novel venom proteins have potential biomedical and therapeutic applications and provide insights into venom evolution.

Venoms contain highly complex mixtures that typically include hundreds of different components and have evolved independently in a diverse range of animals including platypuses, shrews, snakes, lizards, fishes, echinoderms, spiders, wasps, centipedes, sea snails, cephalopods, jellyfish and sea anemones. Many venom genes evolved through gene duplication. Gene duplication occurs in all domains of life and provides the raw substrate from which novel function arise. In this review, we focus on the role that gene duplication has played in the origin and diversification of venom genes. We outline the selective advantages of venom gene duplicates and the role that selection has played in the retention of these duplicates. We use toxin gene intermediates to help trace the evolution of toxin innovation. We also focus on other genomic processes, such as exon and domain duplications, in venom evolution. Finally, we conclude by focusing on the use of high throughput sequencing technology in understanding venom evolution.

Gene duplication followed by adaptive selection is believed to be the primary driver of venom evolution. However, to date, no studies have evaluated the importance of gene duplications for venom evolution using a genomic approach. The availability of a sequenced genome and a venom gland transcriptome for the enigmatic platypus provides a unique opportunity to explore the role that gene duplication plays in venom evolution. Here, we identify gene duplication events and correlate them with expressed transcripts in an in-season venom gland. Gene duplicates (1,508) were identified. These duplicated pairs (421), including genes that have undergone multiple rounds of gene duplications, were expressed in the venom gland. The majority of these genes are involved in metabolism and protein synthesis not toxin functions. Twelve secretory genes including serine proteases, metalloproteinases, and protease inhibitors likely to produce symptoms of envenomation such as vasodilation and pain were detected. Only 16 of 107 platypus genes with high similarity to known toxins evolved through gene duplication. Platypus venom C-type natriuretic peptides and nerve growth factor do not possess lineage-specific gene duplicates. Extensive duplications, believed to increase the potency of toxic content and promote toxin diversification, were not found. This is the first study to take a genome-wide approach in order to examine the impact of gene duplication on venom evolution. Our findings support the idea that adaptive selection acts on gene duplicates to drive the independent evolution and functional diversification of similar venom genes in venomous species. However, gene duplications alone do not explain the “venome” of the platypus. Other mechanisms, such as alternative splicing and mutation, may be important in venom innovation.

Dexmedetomidine is an α(2)-adrenoceptor agonist and anesthetic. The present study was designed to characterize the receptor subtypes and the downstream mechanisms of the vascular effects of dexmedetomidine in small (mesenteric artery) and large (aorta) arteries ex vivo. Isometric tension was measured in Sprague-Dawley rat mesenteric and aortic rings (with or without endothelium). To study relaxations, cumulative concentrations of dexmedetomidine, 5-bromo-N-(2-imidazolin-2-yl)-6-quinoxalinamine, (UK14304), or clonidine were added to rings contracted with 9,11-dideoxy-9α,11α-methanoepoxy prostaglandin F(2α) (U46619) in the presence or absence of indomethacin; N(ω)-nitro-l-arginine methyl ester hydrochloride (l-NAME); 2-[2H-(1-methyl-1,3-dihydroisoindole)methyl]-4,5-dihydroimidazole maleate (BRL44408); 2-[2-(4-(2-methoxyphenyl)piperazin-1-yl)ethyl]-4,4-dimethyl-1,3-(2H,4H)-isoquinolindione dihydrochloride (ARC239); l-657,743, (2S-trans)-1,3,4,5',6,6',7,12b-octahydro-1',3'-dimethyl-spiro[2H-benzofuro[2,3-a]quinolizine-2,4'(1'H)-pyrimidin]-2'(3'H)-one hydrochloride hydrate (MK912); rauwolscine; prazosin; or pertussis toxin. To study contractions, dexmedetomidine was added to quiescent rings without endothelium or after incubation with l-NAME, rauwolscine, prazosin, indomethacin, or 3-[(6-amino-(4-chlorobenzensulfonyl)-2-methyl-5,6,7,8-tetrahydronaphth)-1-yl]propionic acid (S18886). Dexmedetomidine evoked relaxation at low concentrations (10 pM-30 nM) followed by contraction at higher concentrations (>30 nM) in the mesenteric artery. In the aorta, the relaxation was significantly smaller. The relaxation to dexmedetomidine depended on nitric oxide, endothelium, and G(i) protein, and it was mediated by α(2A/D)-adrenoceptors and possibly α(2B)-adrenoceptors. The contraction was mediated mainly by α(2B)- and α(1)-adrenoceptors and involved the action of prostanoids. UK14304 and clonidine induced greater and smaller relaxations, respectively, than dexmedetomidine. In conclusion, depending on the concentration used and the presence of functional endothelium, dexmedetomidine may evoke both relaxation and contraction in isolated arteries. The vascular effects also vary depending on the blood vessel studied. Its vascular effect is different from that of clonidine and UK14304.

Cathepsin L proteases secreted by the helminth pathogen Fasciola hepatica have functions in parasite virulence including tissue invasion and suppression of host immune responses. Using proteomics methods alongside phylogenetic studies we characterized the profile of cathepsin L proteases secreted by adult F. hepatica and hence identified those involved in host-pathogen interaction. Phylogenetic analyses showed that the Fasciola cathepsin L gene family expanded by a series of gene duplications followed by divergence that gave rise to three clades associated with mature adult worms (Clades 1, 2, and 5) and two clades specific to infective juvenile stages (Clades 3 and 4). Consistent with these observations our proteomics studies identified representatives from Clades 1, 2, and 5 but not from Clades 3 and 4 in adult F. hepatica secretory products. Clades 1 and 2 account for 67.39 and 27.63% of total secreted cathepsin Ls, respectively, suggesting that their expansion was positively driven and that these proteases are most critical for parasite survival and adaptation. Sequence comparison studies revealed that the expansion of cathepsin Ls by gene duplication was followed by residue changes in the S2 pocket of the active site. Our biochemical studies showed that these changes result in alterations in substrate binding and suggested that the divergence of the cathepsin L family produced a repertoire of enzymes with overlapping and complementary substrate specificities that could cleave host macromolecules more efficiently. Although the cathepsin Ls are produced as zymogens containing a prosegment and mature domain, all secreted enzymes identified by MS were processed to mature active enzymes. The prosegment region was highly conserved between the clades except at the boundary of prosegment and mature enzyme. Despite the lack of conservation at this section, sites for exogenous cleavage by asparaginyl endopeptidases and a Leu-Ser↓His motif for autocatalytic cleavage by cathepsin Ls were preserved.

When the platypus (Ornithorhynchus anatinus) was first discovered, it was thought to be a taxidermist’s hoax, as it has a blend of mammalian and reptilian features. It is a most remarkable mammal, not only because it lays eggs but also because it is venomous. Rather than delivering venom through a bite, as do snakes and shrews, male platypuses have venomous spurs on each hind leg. The platypus genome sequence provides a unique opportunity to unravel the evolutionary history of many of these interesting features. While searching the platypus genome for the sequences of antimicrobial defensin genes, we identified three Ornithorhynchus venom defensin-like peptide (OvDLP) genes, which produce the major components of platypus venom. We show that gene duplication and subsequent functional diversification of beta-defensins gave rise to these platypus OvDLPs. The OvDLP genes are located adjacent to the beta-defensins and share similar gene organization and peptide structures. Intriguingly, some species of snakes and lizards also produce venoms containing similar molecules called crotamines and crotamine-like peptides. This led us to trace the evolutionary origins of other components of platypus and reptile venom. Here we show that several venom components have evolved separately in the platypus and reptiles. Convergent evolution has repeatedly selected genes coding for proteins containing specific structural motifs as templates for venom molecules.

Spider neurotoxins are commonly used as pharmacological tools and are a popular source of novel compounds with therapeutic and agrochemical potential. Since venom peptides are inherently toxic, the host spider must employ strategies to avoid adverse effects prior to venom use. It is partly for this reason that most spider toxins encode a protective proregion that upon enzymatic cleavage is excised from the mature peptide. In order to identify the mature toxin sequence directly from toxin transcripts, without resorting to protein sequencing, the propeptide cleavage site in the toxin precursor must be predicted bioinformatically. We evaluated different machine learning strategies (support vector machines, hidden Markov model and decision tree) and developed an algorithm (SpiderP) for prediction of propeptide cleavage sites in spider toxins. Our strategy uses a support vector machine (SVM) framework that combines both local and global sequence information. Our method is superior or comparable to current tools for prediction of propeptide sequences in spider toxins. Evaluation of the SVM method on an independent test set of known toxin sequences yielded 96% sensitivity and 100% specificity. Furthermore, we sequenced five novel peptides (not used to train the final predictor) from the venom of the Australian tarantula Selenotypus plumipes to test the accuracy of the predictor and found 80% sensitivity and 99.6% 8-mer specificity. Finally, we used the predictor together with homology information to predict and characterize seven groups of novel toxins from the deeply sequenced venom gland transcriptome of S. plumipes, which revealed structural complexity and innovations in the evolution of the toxins. The precursor prediction tool (SpiderP) is freely available on ArachnoServer (http://www.arachnoserver.org/spiderP.html), a web portal to a comprehensive relational database of spider toxins. All training data, test data, and scripts used are available from the SpiderP website.

Within the mammalian immune system, natural killer (NK) cells contribute to the first line of defence against infectious agents and tumours. Their activity is regulated, in part, by cell surface NK cell receptors. NK receptors can be divided into two unrelated, but functionally analogous superfamilies based on the structure of their extracellular ligand-binding domains. Receptors belonging to the C-type lectin superfamily are predominantly encoded in the natural killer complex (NKC), while receptors belonging to the immunoglobulin superfamily are predominantly encoded in the leukocyte receptor complex (LRC). Natural killer cell receptors are emerging as a rapidly evolving gene family which can display significant intra- and interspecific variation. To date, most studies have focused on eutherian mammals, with significantly less known about the evolution of these receptors in marsupials. Here, we describe the identification of 43 immunoglobulin domain-containing LRC genes in the genome of the Tasmanian devil (Sarcophilus harrisii), the largest remaining marsupial carnivore and only the second marsupial species to be studied. We also identify orthologs of NKC genes KLRK1, CD69, CLEC4E, CLEC1B, CLEC1A and an ortholog of an opossum NKC receptor. Characterisation of these regions in a second, distantly related marsupial provides new insights into the dynamic evolutionary histories of these receptors in mammals. Understanding the functional role of these genes is also important for the development of therapeutic agents against Devil Facial Tumour Disease, a contagious cancer that threatens the Tasmanian devil with extinction.

BACKGROUND: We present the genome sequence of the tammar wallaby, Macropus eugenii, which is a member of the kangaroo family and the first representative of the iconic hopping mammals that symbolize Australia to be sequenced. The tammar has many unusual biological characteristics, including the longest period of embryonic diapause of any mammal, extremely synchronized seasonal breeding and prolonged and sophisticated lactation within a well-defined pouch. Like other marsupials, it gives birth to highly altricial young, and has a small number of very large chromosomes, making it a valuable model for genomics, reproduction and development. RESULTS: The genome has been sequenced to 2 × coverage using Sanger sequencing, enhanced with additional next generation sequencing and the integration of extensive physical and linkage maps to build the genome assembly. We also sequenced the tammar transcriptome across many tissues and developmental time points. Our analyses of these data shed light on mammalian reproduction, development and genome evolution: there is innovation in reproductive and lactational genes, rapid evolution of germ cell genes, and incomplete, locus-specific X inactivation. We also observe novel retrotransposons and a highly rearranged major histocompatibility complex, with many class I genes located outside the complex. Novel microRNAs in the tammar HOX clusters uncover new potential mammalian HOX regulatory elements. CONCLUSIONS: Analyses of these resources enhance our understanding of marsupial gene evolution, identify marsupial-specific conserved non-coding elements and critical genes across a range of biological systems, including reproduction, development and immunity, and provide new insight into marsupial and mammalian biology and genome evolution.

BackgroundTo overcome the increasing resistance of pathogens to existing antibiotics the 10×'20 Initiative declared the urgent need for a global commitment to develop 10 new antimicrobial drugs by the year 2020. Naturally occurring animal antibiotics are an obvious place to start. The recently sequenced genomes of mammals that are divergent from human and mouse, including the tammar wallaby and the platypus, provide an opportunity to discover novel antimicrobials. Marsupials and monotremes are ideal potential sources of new antimicrobials because they give birth to underdeveloped immunologically naïve young that develop outside the sterile confines of a uterus in harsh pathogen-laden environments. While their adaptive immune system develops innate immune factors produced either by the mother or by the young must play a key role in protecting the immune-compromised young. In this study we focus on the cathelicidins, a key family of antimicrobial peptide genes.Principal findingWe identified 14 cathelicidin genes in the tammar wallaby genome and 8 in the platypus genome. The tammar genes were expressed in the mammary gland during early lactation before the adaptive immune system of the young develops, as well as in the skin of the pouch young. Both platypus and tammar peptides were effective in killing a broad range of bacterial pathogens. One potent peptide, expressed in the early stages of tammar lactation, effectively killed multidrug-resistant clinical isolates of Pseudomonas aeruginosa, Klebsiella pneumoniae and Acinetobacter baumannii.Conclusions and significanceMarsupial and monotreme young are protected by antimicrobial peptides that are potent, broad spectrum and salt resistant. The genomes of our distant relatives may hold the key for the development of novel drugs to combat multidrug-resistant pathogens.

The cytokine repertoire of marsupials is largely unknown. The sequencing of the opossum genome has expedited the identification of many immune genes. However, many genes have not been identified using automated annotation pipelines because of high levels of sequence divergence. To fill gaps in our knowledge of the cytokine gene complement in marsupials, we searched the genome assembly of the gray short-tailed opossum for chemokine, interleukin, colony-stimulating factor, tumor necrosis factor, and transforming growth factor genes. In particular, we focused on genes that were not previously identified through Ensembl's automatic annotations. We report that the vast majority of cytokines are conserved, with direct orthologs between therian species. The major exceptions are chemokine genes, which show lineage-specific duplication/loss. Thirty-six chemokines were identified in opossum, including a lineage-specific expansion of macrophage inflammatory protein family genes. Divergent cytokines IL7, IL9, IL31, IL33, and CSF2 were identified. This is the first time IL31 and IL33 have been described outside of eutherian species. The high levels of similarities between the cytokine gene repertoires of therians suggest that the marsupial immune response is highly similar to eutherians.

Background The thymus plays a critical role in the development and maturation of T-cells. Humans have a single thoracic thymus and presence of a second thymus is considered an anomaly. However, many vertebrates have multiple thymuses. The tammar wallaby has two thymuses: a thoracic thymus (typically found in all mammals) and a dominant cervical thymus. Researchers have known about the presence of the two wallaby thymuses since the 1800s, but no genome-wide research has been carried out into possible functional differences between the two thymic tissues. Here, we used pyrosequencing to compare the transcriptomes of a cervical and thoracic thymus from a single 178 day old tammar wallaby. Results We show that both the tammar thoracic and the cervical thymuses displayed gene expression profiles consistent with roles in T-cell development. Both thymuses expressed genes that mediate distinct phases of T-cells differentiation, including the initial commitment of blood stem cells to the T-lineage, the generation of T-cell receptor diversity and development of thymic epithelial cells. Crucial immune genes, such as chemokines were also present. Comparable patterns of expression of non-coding RNAs were seen. 67 genes differentially expressed between the two thymuses were detected, and the possible significance of these results are discussed. Conclusion This is the first study comparing the transcriptomes of two thymuses from a single individual. Our finding supports that both thymuses are functionally equivalent and drive T-cell development. These results are an important first step in the understanding of the genetic processes that govern marsupial immunity, and also allow us to begin to trace the evolution of the mammalian immune system.

The sequencing of the platypus genome represents a significant milestone in the study of mammalian evolution. The platypus is an egg-laying mammal, a member of the order Monotremata, the most divergent mammalian clade. It is the first monotreme to have its genome sequenced. Its unique evolutionary position, as an offshoot between the early divergence of birds and the later emergence of therian (eutherian and marsupials) mammals, provides an unrivalled opportunity to understand the evolution of all mammals. Much like the physical characteristics of the platypus, the platypus genome shows traces of an amalgamation of mammalian and reptilian traits. Key findings to emerge from the genome project include the complex evolution of mammalian sex chromosomes, the transition from egg-laying to live birth in mammals, and the evolution of immune genes and venom molecules in this unique species.

Marsupial immune genes, including cytokines, Natural Killer cell receptors and antimicrobial peptides, have been difficult to isolate due to the wide divergence of these genes from their eutherian counterparts. Early studies using conventional laboratory techniques only isolated the most conserved genes. With the sequencing of opossum and tammar wallaby genomes, many of these divergent genes have now been identified. Knowledge of these gene sequences will enable us to study immune function in marsupials. We can develop immunological reagents for studying the immune system of healthy marsupials and look at the immune response of diseased marsupials. We can also investigate marsupial specific immune functions such as the immune protection of pouch young. Additionally, antimicrobial therapies for both humans and marsupials may be designed around some of these immune peptides.

Background MHC class I antigens are encoded by a rapidly evolving gene family comprising classical and non-classical genes that are found in all vertebrates and involved in diverse immune functions. However, there is a fundamental difference between the organization of class I genes in mammals and non-mammals. Non-mammals have a single classical gene responsible for antigen presentation, which is linked to the antigen processing genes, including TAP. This organization allows co-evolution of advantageous class Ia/TAP haplotypes. In contrast, mammals have multiple classical genes within the MHC, which are separated from the antigen processing genes by class III genes. It has been hypothesized that separation of classical class I genes from antigen processing genes in mammals allowed them to duplicate. We investigated this hypothesis by characterizing the class I genes of the tammar wallaby, a model marsupial that has a novel MHC organization, with class I genes located within the MHC and 10 other chromosomal locations. Results Sequence analysis of 14 BACs containing 15 class I genes revealed that nine class I genes, including one to three classical class I, are not linked to the MHC but are scattered throughout the genome. Kangaroo Endogenous Retroviruses (KERVs) were identified flanking the MHC un-linked class I. The wallaby MHC contains four non-classical class I, interspersed with antigen processing genes. Clear orthologs of non-classical class I are conserved in distant marsupial lineages. Conclusion We demonstrate that classical class I genes are not linked to antigen processing genes in the wallaby and provide evidence that retroviral elements were involved in their movement. The presence of retroviral elements most likely facilitated the formation of recombination hotspots and subsequent diversification of class I genes. The classical class I have moved away from antigen processing genes in eutherian mammals and the wallaby independently, but both lineages appear to have benefited from this loss of linkage by increasing the number of classical genes, perhaps enabling response to a wider range of pathogens. The discovery of non-classical orthologs between distantly related marsupial species is unusual for the rapidly evolving class I genes and may indicate an important marsupial specific function.

The CD33-related sialic acid binding Ig-like lectins (CD33rSiglecs) are predominantly inhibitory receptors expressed on leukocytes. They are distinguishable from conserved Siglecs, such as Sialoadhesin and MAG, by their rapid evolution. A comparison of the CD33rSiglec gene cluster in different mammalian species showed that it can be divided into subclusters, A and B. The two subclusters, inverted in relation to each other, each encode a set of CD33rSiglec genes arranged head-to-tail. Two regions of strong correspondence provided evidence for a large-scale inverse duplication, encompassing the framework CEACAM-18 (CE18) and ATPBD3 (ATB3) genes that seeded the mammalian CD33rSiglec cluster. Phylogenetic analysis was consistent with the predicted inversion. Rodents appear to have undergone wholesale loss of CD33rSiglec genes after the inverse duplication. In contrast, CD33rSiglecs expanded in primates and many are now pseudogenes with features consistent with activating receptors. In contrast to mammals, the fish CD33rSiglecs clusters show no evidence of an inverse duplication. They display greater variation in cluster size and structure than mammals. The close arrangement of other Siglecs and CD33rSiglecs in fish is consistent with a common ancestral region for Siglecs. Expansion of mammalian CD33rSiglecs appears to have followed a large inverse duplication of a smaller primordial cluster over 180 million years ago, prior to eutherian/marsupial divergence. Inverse duplications in general could potentially have a stabilizing effect in maintaining the size and structure of large gene clusters, facilitating the rapid evolution of immune gene families.

Natural killer (NK) cell receptors belong to two unrelated, but functionally analogous gene families: the immunoglobulin superfamily, situated in the leukocyte receptor complex (LRC) and the C-type lectin superfamily, located in the natural killer complex (NKC). Here, we describe the largest NK receptor gene expansion seen to date. We identified 213 putative C-type lectin NK receptor homologs in the genome of the platypus. Many have arisen as the result of a lineage-specific expansion. Orthologs of OLR1, CD69, KLRE, CLEC12B, and CLEC16p genes were also identified. The NKC is split into at least two regions of the genome: 34 genes map to chromosome 7, two map to a small autosome, and the remainder are unanchored in the current genome assembly. No NK receptor genes from the LRC were identified. The massive C-type lectin expansion and lack of Ig-domain-containing NK receptors represents the most extreme polarization of NK receptors found to date. We have used this new data from platypus to trace the possible evolutionary history of the NK receptor clusters.

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

The availability of the first marsupial genome sequence has allowed us to characterize the immunome of the gray short-tailed opossum (Monodelphis domestica). Here we report the identification of key immune genes, including the highly divergent chemokines, defensins, cathelicidins, and Natural Killer cell receptors. It appears that the increase in complexity of the mammalian immune system occurred prior to the divergence of the marsupial and eutherian lineages ∼180 million years ago. Genomes of ancestral mammals most likely contained all of the key mammalian immune gene families, with evolution on different continents, in the presence of different pathogens leading to lineage specific expansions and contractions, resulting in some minor differences in gene number and composition between different mammalian lineages. Gene expansion and extensive heterogeneity in opossum antimicrobial peptide genes may have evolved as a consequence of the newborn young needing to survive without an adaptive immune system in a pathogen laden environment. Given the similarities in the genomic architecture of the marsupial and eutherian immune systems, we propose that marsupials are ideal model organisms for the study of developmental immunology.

Background: To date, few peptides in the complex mixture of platypus venom have been identified and sequenced, in part due to the limited amounts of platypus venom available to study. We have constructed and sequenced a cDNA library from an active platypus venom gland to identify the remaining components. Results: We identified 83 novel putative platypus venom genes from 13 toxin families, which are homologous to known toxins from a wide range of vertebrates (fish, reptiles, insectivores) and invertebrates (spiders, sea anemones, starfish). A number of these are expressed in tissues other than the venom gland, and at least three of these families (those with homology to toxins from distant invertebrates) may play non-toxin roles. Thus, further functional testing is required to confirm venom activity. However, thepresence of similar putative toxins in such widely divergent species provides further evidence for the hypothesis that there are certain protein families that are selected preferentially during evolution to become venom peptides. We have also used homology with known proteins to speculate on the contributions of each venom component to the symptoms of platypus envenomation. Conclusions: This study represents a step towards fully characterizing the first mammal venom transcriptome. We have found similarities between putative platypus toxins and those of a number of unrelated species, providing insight into the evolution of mammalian venom.

The sequencing of the platypus genome has spurred investigations into the characterisation of the monotreme immune response. As the most divergent of extant mammals, the characterisation of the monotreme immune repertoire allows us to trace the evolutionary history of immunity in mammals and provide insights into the immune gene complement of ancestral mammals. The immune system of monotremes has remained largely uncharacterised due to the lack of specific immunological reagents and limited access to animals for experimentation. Early immunological studies focussed on the anatomy and physiology of the lymphoid system in the platypus. More recent molecular studies have focussed on characterisation of individual immunoglobulin, T-cell receptor and MHC genes in both the platypus and short-beaked echidna. Here, we review the published literature on the monotreme immune gene repertoire and provide new data generated from genome analysis on cytokines, Fc receptors and immunoglobulins. We present an overview of key gene families responsible for innate and adaptive immunity including the cathelicidins, defensins, T-cell receptors and the major histocompatibility complex (MHC) Class I and Class II antigens. We comment on the usefulness of these sequences for future studies into immunity, health and disease in monotremes.

Nature 453, 175–183 (2008) In this Article, Mikhail Nefedov and Pieter J. de Jong were omitted from the author list.

Background: Cytokines are small proteins that regulate immunity in vertebrate species. Marsupial and eutherian mammals last shared a common ancestor more than 180 million years ago, so it is not surprising that attempts to isolate many key marsupial cytokines using traditional laboratory techniques have been unsuccessful. This paucity of molecular data has led some authors to suggest that the marsupial immune system is \'primitive\' and not on par with the sophisticated immune system of eutherian (placental) mammals. Results: The sequencing of the first marsupial genome has allowed us to identify highly divergent immune genes. We used gene prediction methods that incorporate the identification of gene location using BLAST, SYNTENY + BLAST and HMMER to identify 23 key marsupial immune genes, including IFN-γ, IL-2, IL-4, IL-6, IL-12 and IL-13, in the genome of the grey short-tailed opossum (Monodelphis domestica). Many of these genes were not predicted in the publicly available automated annotations. Conclusion: The power of this approach was demonstrated by the identification of orthologous cytokines between marsupials and eutherians that share only 30% identity at the amino acid level. Furthermore, the presence of key immunological genes suggests that marsupials do indeed possess a sophisticated immune system, whose function may parallel that of eutherian mammals.