Elephant Genomes Reveal Accelerated Evolution in Mechanisms Underlying Disease Defenses.

Disease susceptibility and resistance are important factors for the conservation of endangered species, including elephants. We analyzed pathology data from 26 zoos and report that Asian elephants have increased neoplasia and malignancy prevalence compared with African bush elephants. This is consistent with observed higher susceptibility to tuberculosis and elephant endotheliotropic herpesvirus (EEHV) in Asian elephants. To investigate genetic mechanisms underlying disease resistance, including differential responses between species, among other elephant traits, we sequenced multiple elephant genomes. We report a draft assembly for an Asian elephant, and defined 862 and 1,017 conserved potential regulatory elements in Asian and African bush elephants, respectively. In the genomes of both elephant species, conserved elements were significantly enriched with genes differentially expressed between the species. In Asian elephants, these putative regulatory regions were involved in immunity pathways including tumor-necrosis factor, which plays an important role in EEHV response. Genomic sequences of African bush, forest, and Asian elephant genomes revealed extensive sequence conservation at TP53 retrogene loci across three species, which may be related to TP53 functionality in elephant cancer resistance. Positive selection scans revealed outlier genes related to additional elephant traits. Our study suggests that gene regulation plays an important role in the differential inflammatory response of Asian and African elephants, leading to increased infectious disease and cancer susceptibility in Asian elephants. These genomic discoveries can inform future functional and translational studies aimed at identifying effective treatment approaches for ill elephants, which may improve conservation.

Haplotyping the Vitis collinear core genome with rhAmpSeq improves marker transferability in a diverse genus.

Transferable DNA markers are essential for breeding and genetics. Grapevine (Vitis) breeders utilize disease resistance alleles from congeneric species ~20 million years divergent, but existing Vitis marker platforms have cross-species transfer rates as low as 2%. Here, we apply a marker strategy targeting the inferred Vitis core genome. Incorporating seven linked-read de novo assemblies and three existing assemblies, the Vitis collinear core genome is estimated to converge at 39.8 Mb (8.67% of the genome). Adding shotgun genome sequences from 40 accessions enables identification of conserved core PCR primer binding sites flanking polymorphic haplotypes with high information content. From these target regions, we develop 2,000 rhAmpSeq markers as a PCR multiplex and validate the panel in four biparental populations spanning the diversity of the Vitis genus, showing transferability increases to 91.9%. This marker development strategy should be widely applicable for genetic studies in many taxa, particularly those ~20 million years divergent.

The maize W22 genome provides a foundation for functional genomics and transposon biology.

The maize W22 inbred has served as a platform for maize genetics since the mid twentieth century. To streamline maize genome analyses, we have sequenced and de novo assembled a W22 reference genome using short-read sequencing technologies. We show that significant structural heterogeneity exists in comparison to the B73 reference genome at multiple scales, from transposon composition and copy number variation to single-nucleotide polymorphisms. The generation of this reference genome enables accurate placement of thousands of Mutator (Mu) and Dissociation (Ds) transposable element insertions for reverse and forward genetics studies. Annotation of the genome has been achieved using RNA-seq analysis, differential nuclease sensitivity profiling and bisulfite sequencing to map open reading frames, open chromatin sites and DNA methylation profiles, respectively. Collectively, the resources developed here integrate W22 as a community reference genome for functional genomics and provide a foundation for the maize pan-genome.

Regulatory Divergence in Wound-Responsive Gene Expression between Domesticated and Wild Tomato.

The evolution of transcriptional regulatory mechanisms is central to how stress response and tolerance differ between species. However, it remains largely unknown how divergence in cis-regulatory sites and, subsequently, transcription factor (TF) binding specificity contribute to stress-responsive expression divergence, particularly between wild and domesticated species. By profiling wound-responsive gene transcriptomes in wild Solanum pennellii and domesticated S. lycopersicum, we found extensive wound response divergence and identified 493 S. lycopersicum and 278 S. pennellii putative cis-regulatory elements (pCREs) that were predictive of wound-responsive gene expression. Only 24-52% of these wound response pCREs (depending on wound response patterns) were consistently enriched in the putative promoter regions of wound-responsive genes across species. In addition, between these two species, their differences in pCRE site sequences were significantly and positively correlated with differences in wound-responsive gene expression. Furthermore, ∼11-39% of pCREs were specific to only one of the species and likely bound by TFs from different families. These findings indicate substantial regulatory divergence in these two plant species that diverged ∼3-7 million years ago. Our study provides insights into the mechanistic basis of how the transcriptional response to wounding is regulated and, importantly, the contribution of cis-regulatory components to variation in wound-responsive gene expression between a wild and a domesticated plant species.

Predicting gene structure changes resulting from genetic variants via exon definition features.

Motivation: Genetic variation that disrupts gene function by altering gene splicing between individuals can substantially influence traits and disease. In those cases, accurately predicting the effects of genetic variation on splicing can be highly valuable for investigating the mechanisms underlying those traits and diseases. While methods have been developed to generate high quality computational predictions of gene structures in reference genomes, the same methods perform poorly when used to predict the potentially deleterious effects of genetic changes that alter gene splicing between individuals. Underlying that discrepancy in predictive ability are the common assumptions by reference gene finding algorithms that genes are conserved, well-formed and produce functional proteins. Results: We describe a probabilistic approach for predicting recent changes to gene structure that may or may not conserve function. The model is applicable to both coding and non-coding genes, and can be trained on existing gene annotations without requiring curated examples of aberrant splicing. We apply this model to the problem of predicting altered splicing patterns in the genomes of individual humans, and we demonstrate that performing gene-structure prediction without relying on conserved coding features is feasible. The model predicts an unexpected abundance of variants that create de novo splice sites, an observation supported by both simulations and empirical data from RNA-seq experiments. While these de novo splice variants are commonly misinterpreted by other tools as coding or non-coding variants of little or no effect, we find that in some cases they can have large effects on splicing activity and protein products and we propose that they may commonly act as cryptic factors in disease. Availability and implementation: The software is available from geneprediction.org/SGRF. Supplementary information: Supplementary information is available at Bioinformatics online.

Improved maize reference genome with single-molecule technologies.

Complete and accurate reference genomes and annotations provide fundamental tools for characterization of genetic and functional variation. These resources facilitate the determination of biological processes and support translation of research findings into improved and sustainable agricultural technologies. Many reference genomes for crop plants have been generated over the past decade, but these genomes are often fragmented and missing complex repeat regions. Here we report the assembly and annotation of a reference genome of maize, a genetic and agricultural model species, using single-molecule real-time sequencing and high-resolution optical mapping. Relative to the previous reference genome, our assembly features a 52-fold increase in contig length and notable improvements in the assembly of intergenic spaces and centromeres. Characterization of the repetitive portion of the genome revealed more than 130,000 intact transposable elements, allowing us to identify transposable element lineage expansions that are unique to maize. Gene annotations were updated using 111,000 full-length transcripts obtained by single-molecule real-time sequencing. In addition, comparative optical mapping of two other inbred maize lines revealed a prevalence of deletions in regions of low gene density and maize lineage-specific genes.

Correction: Genome, Functional Gene Annotation, and Nuclear Transformation of the Heterokont Oleaginous Alga Nannochloropsis oceanica CCMP1779.

[This corrects the article DOI: 10.1371/journal.pgen.1003064.].

Cellular and ultrastructural characterization of the grey-morph phenotype in southern right whales (Eubalaena australis).

Southern right whales (SRWs, Eubalena australis) are polymorphic for an X-linked pigmentation pattern known as grey morphism. Most SRWs have completely black skin with white patches on their bellies and occasionally on their backs; these patches remain white as the whale ages. Grey morphs (previously referred to as partial albinos) appear mostly white at birth, with a splattering of rounded black marks; but as the whales age, the white skin gradually changes to a brownish grey color. The cellular and developmental bases of grey morphism are not understood. Here we describe cellular and ultrastructural features of grey-morph skin in relation to that of normal, wild-type skin. Melanocytes were identified histologically and counted, and melanosomes were measured using transmission electron microscopy. Grey-morph skin had fewer melanocytes when compared to wild-type skin, suggesting reduced melanocyte survival, migration, or proliferation in these whales. Grey-morph melanocytes had smaller melanosomes relative to wild-type skin, normal transport of melanosomes to surrounding keratinocytes, and normal localization of melanin granules above the keratinocyte nuclei. These findings indicate that SRW grey-morph pigmentation patterns are caused by reduced numbers of melanocytes in the skin, as well as by reduced amounts of melanin production and/or reduced sizes of mature melanosomes. Grey morphism is distinct from piebaldism and albinism found in other species, which are genetic pigmentation conditions resulting from the local absence of melanocytes, or the inability to synthesize melanin, respectively.

High-throughput interpretation of gene structure changes in human and nonhuman resequencing data, using ACE.

MOTIVATION: The accurate interpretation of genetic variants is critical for characterizing genotype-phenotype associations. Because the effects of genetic variants can depend strongly on their local genomic context, accurate genome annotations are essential. Furthermore, as some variants have the potential to disrupt or alter gene structure, variant interpretation efforts stand to gain from the use of individualized annotations that account for differences in gene structure between individuals or strains. RESULTS: We describe a suite of software tools for identifying possible functional changes in gene structure that may result from sequence variants. ACE ('Assessing Changes to Exons') converts phased genotype calls to a collection of explicit haplotype sequences, maps transcript annotations onto them, detects gene-structure changes and their possible repercussions, and identifies several classes of possible loss of function. Novel transcripts predicted by ACE are commonly supported by spliced RNA-seq reads, and can be used to improve read alignment and transcript quantification when an individual-specific genome sequence is available. Using publicly available RNA-seq data, we show that ACE predictions confirm earlier results regarding the quantitative effects of nonsense-mediated decay, and we show that predicted loss-of-function events are highly concordant with patterns of intolerance to mutations across the human population. ACE can be readily applied to diverse species including animals and plants, making it a broadly useful tool for use in eukaryotic population-based resequencing projects, particularly for assessing the joint impact of all variants at a locus. AVAILABILITY AND IMPLEMENTATION: ACE is written in open-source C ++ and Perl and is available from geneprediction.org/ACE. CONTACT: myandell@genetics.utah.edu or tim.reddy@duke.edu. SUPPLEMENTARY INFORMATION: Supplementary information is available at Bioinformatics online.

The spotted gar genome illuminates vertebrate evolution and facilitates human-teleost comparisons.

To connect human biology to fish biomedical models, we sequenced the genome of spotted gar (Lepisosteus oculatus), whose lineage diverged from teleosts before teleost genome duplication (TGD). The slowly evolving gar genome has conserved in content and size many entire chromosomes from bony vertebrate ancestors. Gar bridges teleosts to tetrapods by illuminating the evolution of immunity, mineralization and development (mediated, for example, by Hox, ParaHox and microRNA genes). Numerous conserved noncoding elements (CNEs; often cis regulatory) undetectable in direct human-teleost comparisons become apparent using gar: functional studies uncovered conserved roles for such cryptic CNEs, facilitating annotation of sequences identified in human genome-wide association studies. Transcriptomic analyses showed that the sums of expression domains and expression levels for duplicated teleost genes often approximate the patterns and levels of expression for gar genes, consistent with subfunctionalization. The gar genome provides a resource for understanding evolution after genome duplication, the origin of vertebrate genomes and the function of human regulatory sequences.

An Introduction to Genome Annotation.

Genome projects have evolved from large international undertakings to tractable endeavors for a single lab. Accurate genome annotation is critical for successful genomic, genetic, and molecular biology experiments. These annotations can be generated using a number of approaches and available software tools. This unit describes methods for genome annotation and a number of software tools commonly used in gene annotation.

Potential Mechanisms for Cancer Resistance in Elephants and Comparative Cellular Response to DNA Damage in Humans.

IMPORTANCE: Evolutionary medicine may provide insights into human physiology and pathophysiology, including tumor biology. OBJECTIVE: To identify mechanisms for cancer resistance in elephants and compare cellular response to DNA damage among elephants, healthy human controls, and cancer-prone patients with Li-Fraumeni syndrome (LFS). DESIGN, SETTING, AND PARTICIPANTS: A comprehensive survey of necropsy data was performed across 36 mammalian species to validate cancer resistance in large and long-lived organisms, including elephants (n = 644). The African and Asian elephant genomes were analyzed for potential mechanisms of cancer resistance. Peripheral blood lymphocytes from elephants, healthy human controls, and patients with LFS were tested in vitro in the laboratory for DNA damage response. The study included African and Asian elephants (n = 8), patients with LFS (n = 10), and age-matched human controls (n = 11). Human samples were collected at the University of Utah between June 2014 and July 2015. EXPOSURES: Ionizing radiation and doxorubicin. MAIN OUTCOMES AND MEASURES: Cancer mortality across species was calculated and compared by body size and life span. The elephant genome was investigated for alterations in cancer-related genes. DNA repair and apoptosis were compared in elephant vs human peripheral blood lymphocytes. RESULTS: Across mammals, cancer mortality did not increase with body size and/or maximum life span (eg, for rock hyrax, 1% [95% CI, 0%-5%]; African wild dog, 8% [95% CI, 0%-16%]; lion, 2% [95% CI, 0%-7%]). Despite their large body size and long life span, elephants remain cancer resistant, with an estimated cancer mortality of 4.81% (95% CI, 3.14%-6.49%), compared with humans, who have 11% to 25% cancer mortality. While humans have 1 copy (2 alleles) of TP53, African elephants have at least 20 copies (40 alleles), including 19 retrogenes (38 alleles) with evidence of transcriptional activity measured by reverse transcription polymerase chain reaction. In response to DNA damage, elephant lymphocytes underwent p53-mediated apoptosis at higher rates than human lymphocytes proportional to TP53 status (ionizing radiation exposure: patients with LFS, 2.71% [95% CI, 1.93%-3.48%] vs human controls, 7.17% [95% CI, 5.91%-8.44%] vs elephants, 14.64% [95% CI, 10.91%-18.37%]; P < .001; doxorubicin exposure: human controls, 8.10% [95% CI, 6.55%-9.66%] vs elephants, 24.77% [95% CI, 23.0%-26.53%]; P < .001). CONCLUSIONS AND RELEVANCE: Compared with other mammalian species, elephants appeared to have a lower-than-expected rate of cancer, potentially related to multiple copies of TP53. Compared with human cells, elephant cells demonstrated increased apoptotic response following DNA damage. These findings, if replicated, could represent an evolutionary-based approach for understanding mechanisms related to cancer suppression.

Automated update, revision, and quality control of the maize genome annotations using MAKER-P improves the B73 RefGen_v3 gene models and identifies new genes.

The large size and relative complexity of many plant genomes make creation, quality control, and dissemination of high-quality gene structure annotations challenging. In response, we have developed MAKER-P, a fast and easy-to-use genome annotation engine for plants. Here, we report the use of MAKER-P to update and revise the maize (Zea mays) B73 RefGen_v3 annotation build (5b+) in less than 3 h using the iPlant Cyberinfrastructure. MAKER-P identified and annotated 4,466 additional, well-supported protein-coding genes not present in the 5b+ annotation build, added additional untranslated regions to 1,393 5b+ gene models, identified 2,647 5b+ gene models that lack any supporting evidence (despite the use of large and diverse evidence data sets), identified 104,215 pseudogene fragments, and created an additional 2,522 noncoding gene annotations. We also describe a method for de novo training of MAKER-P for the annotation of newly sequenced grass genomes. Collectively, these results lead to the 6a maize genome annotation and demonstrate the utility of MAKER-P for rapid annotation, management, and quality control of grasses and other difficult-to-annotate plant genomes.

Genome Annotation and Curation Using MAKER and MAKER-P.

This unit describes how to use the genome annotation and curation tools MAKER and MAKER-P to annotate protein-coding and noncoding RNA genes in newly assembled genomes, update/combine legacy annotations in light of new evidence, add quality metrics to annotations from other pipelines, and map existing annotations to a new assembly. MAKER and MAKER-P can rapidly annotate genomes of any size, and scale to match available computational resources.

Gibbon genome and the fast karyotype evolution of small apes.

Gibbons are small arboreal apes that display an accelerated rate of evolutionary chromosomal rearrangement and occupy a key node in the primate phylogeny between Old World monkeys and great apes. Here we present the assembly and analysis of a northern white-cheeked gibbon (Nomascus leucogenys) genome. We describe the propensity for a gibbon-specific retrotransposon (LAVA) to insert into chromosome segregation genes and alter transcription by providing a premature termination site, suggesting a possible molecular mechanism for the genome plasticity of the gibbon lineage. We further show that the gibbon genera (Nomascus, Hylobates, Hoolock and Symphalangus) experienced a near-instantaneous radiation ∼5 million years ago, coincident with major geographical changes in southeast Asia that caused cycles of habitat compression and expansion. Finally, we identify signatures of positive selection in genes important for forelimb development (TBX5) and connective tissues (COL1A1) that may have been involved in the adaptation of gibbons to their arboreal habitat.

Intrauterine growth restriction increases TNF α and activates the unfolded protein response in male rat pups.

Intrauterine growth restriction (IUGR) programs adult disease, including obesity and insulin resistance. Our group previously demonstrated that IUGR dysregulates adipose deposition in male, but not female, weanling rats. Dysregulated adipose deposition is often accompanied by the release of proinflammatory signaling molecules, such as tumor necrosis factor alpha (TNF α ). TNF α contributes to adipocyte inflammation and impaired insulin signaling. TNF α has also been implicated in the activation of the unfolded protein response (UPR), which impairs insulin signaling. We hypothesized that, in male rat pups, IUGR would increase TNF α , TNFR1, and components of the UPR (Hspa5, ATF6, p-eIF2 α , and Ddit3) prior to the onset of obesity. We further hypothesized that impaired glucose tolerance would occur after the onset of adipose dysfunction in male IUGR rats. To test this hypothesis, we used a well-characterized rat model of uteroplacental insufficiency-induced IUGR. Our primary findings are that, in male rats, IUGR (1) increased circulating and adipose TNF α , (2) increased mRNA levels of UPR components as well as p-eIF2a, and (3) impaired glucose tolerance after observed TNF α increased and after UPR activation. We speculate that programmed dysregulation of TNF α and UPR contributed to the development of glucose intolerance in male IUGR rats.

Maternal tobacco smoke increased visceral adiposity and serum corticosterone levels in adult male rat offspring.

BACKGROUND: Maternal tobacco smoke (MTS) predisposes human and rat offspring to visceral obesity in early adulthood. Glucocorticoid excess also causes visceral obesity. We hypothesized that in utero MTS would increase visceral adiposity and alter the glucocorticoid pathway in young adult rats. METHODS: We developed a novel model of in utero MTS exposure in pregnant rats by exposing them to cigarette smoke from E11.5 to term. Neonatal rats were cross-fostered to control dams and weaned to standard rat chow through young adulthood (postnatal day 60). RESULTS: We demonstrated increased visceral adiposity (193%)*, increased visceral adipose 11-ß hydroxysteroid dehydrogenase 1 mRNA (204%)*, increased serum corticosterone (147%)*, and no change in glucocorticoid receptor protein in adult male MTS rat offspring. Female rats exposed to MTS in utero demonstrated no change in visceral or subcutaneous adiposity, decreased serum corticosterone (60%)*, and decreased adipose glucocorticoid receptor protein (66%)*. *P < 0.05. CONCLUSION: We conclude that in utero MTS exposure increased visceral adiposity and altered in the glucocorticoid pathway in a sex-specific manner. We speculate that in utero MTS exposure programs adipose dysfunction in adult male rat offspring via alteration in the glucocorticoid pathway.

MAKER-P: a tool kit for the rapid creation, management, and quality control of plant genome annotations.

We have optimized and extended the widely used annotation engine MAKER in order to better support plant genome annotation efforts. New features include better parallelization for large repeat-rich plant genomes, noncoding RNA annotation capabilities, and support for pseudogene identification. We have benchmarked the resulting software tool kit, MAKER-P, using the Arabidopsis (Arabidopsis thaliana) and maize (Zea mays) genomes. Here, we demonstrate the ability of the MAKER-P tool kit to automatically update, extend, and revise the Arabidopsis annotations in light of newly available data and to annotate pseudogenes and noncoding RNAs absent from The Arabidopsis Informatics Resource 10 build. Our results demonstrate that MAKER-P can be used to manage and improve the annotations of even Arabidopsis, perhaps the best-annotated plant genome. We have also installed and benchmarked MAKER-P on the Texas Advanced Computing Center. We show that this public resource can de novo annotate the entire Arabidopsis and maize genomes in less than 3 h and produce annotations of comparable quality to those of the current The Arabidopsis Information Resource 10 and maize V2 annotation builds.

The African coelacanth genome provides insights into tetrapod evolution.

The discovery of a living coelacanth specimen in 1938 was remarkable, as this lineage of lobe-finned fish was thought to have become extinct 70 million years ago. The modern coelacanth looks remarkably similar to many of its ancient relatives, and its evolutionary proximity to our own fish ancestors provides a glimpse of the fish that first walked on land. Here we report the genome sequence of the African coelacanth, Latimeria chalumnae. Through a phylogenomic analysis, we conclude that the lungfish, and not the coelacanth, is the closest living relative of tetrapods. Coelacanth protein-coding genes are significantly more slowly evolving than those of tetrapods, unlike other genomic features. Analyses of changes in genes and regulatory elements during the vertebrate adaptation to land highlight genes involved in immunity, nitrogen excretion and the development of fins, tail, ear, eye, brain and olfaction. Functional assays of enhancers involved in the fin-to-limb transition and in the emergence of extra-embryonic tissues show the importance of the coelacanth genome as a blueprint for understanding tetrapod evolution.