Tissue-specific transcriptome sequencing analysis expands the non-human primate reference transcriptome resource (NHPRTR)

Compared to other NHEJ factors, DNA-PK protein and RNA levels are markedly increased in all higher primates, but not in prosimians or other mammals

The DNA dependent protein kinase (DNA-PK) initiates non-homologous recombination (NHEJ), the predominate DNA double-strand break (DSBR) pathway in higher vertebrates. It has been known for decades that the enzymatic activity of DNA-PK [that requires its three component polypeptides, Ku70, Ku80 (that comprise the DNA-end binding Ku heterodimer), and the catalytic subunit (DNA-PKcs)] is present in humans at 10-50 times the level observed in other mammals. Here, we show that the high level of DNA-PKcs protein expression appears evolutionarily in mammals between prosimians and higher primates. Moreover, the RNAs encoding the three component polypeptides of DNA-PK are present at similarly high levels in hominids, new-, and old-world monkeys, but expression of these RNAs in prosimians is ∼5-50 fold less, analogous to the levels observed in other non-primate species. This is reminiscent of the appearance of Alu repeats in primate genomes -- abundant in higher primates, but present at much lower density in prosimians. Alu repeats are well-known for their capacity to promote non-allelic homologous recombination (NAHR) a process known to be inhibited by DNA-PK. Nanopore sequence analyses of cultured cells proficient or deficient in DNA-PK revealed an increase of inter-chromosomal translocations caused by NAHR. Although the high levels of DNA-PK in primates may have many functions, we posit that high levels of DNA-PK may function to restrain deleterious NAHR events between Alu elements.

TERT Expression and Clinical Outcome in Pulmonary Carcinoids

PURPOSE

The clinical course of pulmonary carcinoids ranges from indolent to fatal disease, suggesting that specific molecular alterations drive progression toward the fully malignant state. A similar spectrum of clinical phenotypes occurs in pediatric neuroblastoma, in which activation of telomerase reverse transcriptase (TERT) is decisive in determining the course of disease. We therefore investigated whether TERT expression defines the clinical fate of patients with pulmonary carcinoid.

METHODS

TERT expression was examined by RNA sequencing in a test cohort and a validation cohort of pulmonary carcinoids (n = 88 and n = 105, respectively). A natural TERT expression cutoff was determined in the test cohort on the basis of the distribution of TERT expression, and its prognostic value was assessed by Kaplan-Meier survival estimates and multivariable analyses. Telomerase activity was validated by telomere repeat amplification protocol assay.

RESULTS

Similar to neuroblastoma, TERT expression exhibited a bimodal distribution in pulmonary carcinoids, separating tumors into TERT-high and TERT-low subgroups. A natural TERT cutoff discriminated unfavorable from favorable clinical courses with high accuracy both in the test cohort (5-year overall survival [OS], 0.547 ± 0.132 v 1.0; P < .001) and the validation cohort (5-year OS, 0.788 ± 0.063 v 0.913 ± 0.048; P < .001). In line with these findings, telomerase activity was largely absent in TERT-low tumors, whereas it was readily detectable in TERT-high carcinoids. In multivariable analysis considering TERT expression, histology (typical v atypical carcinoid), and stage (≤IIA v ≥IIB), high TERT expression was an independent prognostic marker for poor survival, with a hazard ratio of 5.243 (95% CI, 1.943 to 14.148; P = .001).

CONCLUSION

Our data demonstrate that high TERT expression defines clinically aggressive pulmonary carcinoids with fatal outcome, similar to neuroblastoma, indicating that activation of TERT may be a defining feature of lethal cancers.

Lineage-specific splicing regulation of MAPT gene in the primate brain

Divergence of precursor messenger RNA (pre-mRNA) alternative splicing (AS) is widespread in mammals, including primates, but the underlying mechanisms and functional impact are poorly understood. Here, we modeled cassette exon inclusion in primate brains as a quantitative trait and identified 1,170 (∼3%) exons with lineage-specific splicing shifts under stabilizing selection. Among them, microtubule-associated protein tau (MAPT) exons 2 and 10 underwent anticorrelated, two-step evolutionary shifts in the catarrhine and hominoid lineages, leading to their present inclusion levels in humans. The developmental-stage-specific divergence of exon 10 splicing, whose dysregulation can cause frontotemporal lobar degeneration (FTLD), is mediated by divergent distal intronic MBNL-binding sites. Competitive binding of these sites by CRISPR-dCas13d/gRNAs effectively reduces exon 10 inclusion, potentially providing a therapeutically compatible approach to modulate tau isoform expression. Our data suggest adaptation of MAPT function and, more generally, a role for AS in the evolutionary expansion of the primate brain.

Correlated evolution of social organization and lifespan in mammals

Discerning the relationship between sociality and longevity would permit a deeper understanding of how animal life history evolved. Here, we perform a phylogenetic comparative analysis of ~1000 mammalian species on three states of social organization (solitary, pair-living, and group-living) and longevity. We show that group-living species generally live longer than solitary species, and that the transition rate from a short-lived state to a long-lived state is higher in group-living than non-group-living species, altogether supporting the correlated evolution of social organization and longevity. The comparative brain transcriptomes of 94 mammalian species identify 31 genes, hormones and immunity-related pathways broadly involved in the association between social organization and longevity. Further selection features reveal twenty overlapping pathways under selection for both social organization and longevity. These results underscore a molecular basis for the influence of the social organization on longevity.

Distinct dosage compensations of ploidy-sensitive and -insensitive X chromosome genes during development and in diseases

The active X chromosome in mammals is upregulated to balance its dosage to autosomes during evolution. However, it is elusive why the known dosage compensation machinery showed uneven and small influence on X genes. Here, based on >20,000 transcriptomes, we identified two X gene groups (ploidy-sensitive [PSX] and ploidy-insensitive [PIX]), showing distinct but evolutionarily conserved dosage compensations (termed XAR). We demonstrated that XAR-PIX was downregulated whereas XAR-PSX upregulated at both RNA and protein levels across cancer types, in contrast with their trends during stem cell differentiation. XAR-PIX, but not XAR-PSX, was lower and correlated with autoantibodies and inflammation in patients of lupus, suggesting that insufficient dosage of PIX genes contribute to lupus pathogenesis. We further identified and experimentally validated two XAR regulators, TP53 and ATRX. Collectively, we provided insights into X dosage compensation in mammals and demonstrated different regulation of PSX and PIX and their pathophysiological roles in human diseases.

Current advances in primate genomics: novel approaches for understanding evolution and disease

Primate genomics holds the key to understanding fundamental aspects of human evolution and disease. However, genetic diversity and functional genomics data sets are currently available for only a few of the more than 500 extant primate species. Concerted efforts are under way to characterize primate genomes, genetic polymorphism and divergence, and functional landscapes across the primate phylogeny. The resulting data sets will enable the connection of genotypes to phenotypes and provide new insight into aspects of the genetics of primate traits, including human diseases. In this Review, we describe the existing genome assemblies as well as genetic variation and functional genomic data sets. We highlight some of the challenges with sample acquisition. Finally, we explore how technological advances in single-cell functional genomics and induced pluripotent stem cell-derived organoids will facilitate our understanding of the molecular foundations of primate biology.

Convergent evolution of antiviral machinery derived from endogenous retrovirus truncated envelope genes in multiple species

Host genetic resistance to viral infection controls the pathogenicity and epidemic dynamics of infectious diseases. Refrex-1 is a restriction factor against feline leukemia virus subgroup D (FeLV-D) and an endogenous retrovirus (ERV) in domestic cats (ERV-DC). Refrex-1 is encoded by a subset of ERV-DC loci with truncated envelope genes and secreted from cells as a soluble protein. Here, we identified the copper transporter CTR1 as the entry receptor for FeLV-D and genotype I ERV-DCs. We also identified CTR1 as a receptor for primate ERVs from crab-eating macaques and rhesus macaques, which were found in a search of intact envelope genes capable of forming infectious viruses. Refrex-1 counteracted infection by FeLV-D and ERV-DCs via competition for the entry receptor CTR1; the antiviral effects extended to primate ERVs with CTR1-dependent entry. Furthermore, truncated ERV envelope genes found in chimpanzee, bonobo, gorilla, crab-eating macaque, and rhesus macaque genomes could also block infection by feline and primate retroviruses. Genetic analyses showed that these ERV envelope genes were acquired in a species- or genus-specific manner during host evolution. These results indicated that soluble envelope proteins could suppress retroviral infection across species boundaries, suggesting that they function to control retroviral spread. Our findings revealed that several mammalian species acquired antiviral machinery from various ancient retroviruses, leading to convergent evolution for host defense.

Death-associated protein kinases and intestinal epithelial homeostasis

The family of death-associated protein kinases (DAPKs) and DAPK-related apoptosis-inducing protein kinases (DRAKs) act as molecular switches for a multitude of cellular processes, including apoptotic and autophagic cell death events. This review summarizes the mechanisms for kinase activity regulation and discusses recent molecular investigations of DAPK and DRAK family members in the intestinal epithelium. In general, recent literature convincingly supports the importance of this family of protein kinases in the homeostatic processes that govern the proper function of the intestinal epithelium. Each of the DAPK family of proteins possesses distinct biochemical properties, and we compare similarities in the information available as well as those cases where functional distinctions are apparent. As the prototypical member of the family, DAPK1 is noteworthy for its tumor suppressor function and association with colorectal cancer. In the intestinal epithelium, DAPK2 is associated with programmed cell death, potential tumor-suppressive functions, and a unique influence on granulocyte biology. The impact of the DRAKs in the epithelium is understudied, but recent studies support a role for DRAK1 in inflammation-mediated tumor growth and metastasis. A commentary is provided on the potential importance of DAPK3 in facilitating epithelial restitution and wound healing during the resolution of colitis. An update on efforts to develop selective pharmacologic effectors of individual DAPK members is also supplied.

Low complexity regions in mammalian proteins are associated with low protein abundance and high transcript abundance

Low Complexity Regions (LCRs) are present in a surprisingly large number of eukaryotic proteins. These highly repetitive and compositionally biased sequences are often structurally disordered, bind promiscuously, and evolve rapidly. Frequently studied in terms of evolutionary dynamics, little is known about how LCRs affect the expression of the proteins which contain them. It would be expected that rapidly evolving LCRs are unlikely to be tolerated in strongly conserved, highly abundant proteins, leading to lower overall abundance in proteins which contain LCRs. To test this hypothesis and examine the associations of protein abundance and transcript abundance with the presence of LCRs, we have integrated high-throughput data from across mammals. We have found that LCRs are indeed associated with reduced protein abundance, but are also associated with elevated transcript abundance. These associations are qualitatively consistent across 12 human tissues and nine mammalian species. The differential impacts of LCRs on abundance at the protein and transcript level are not explained by differences in either protein degradation rates or the inefficiency of translation for LCR containing proteins. We suggest that rapidly evolving LCRs are a source of selective pressure on the regulatory mechanisms which maintain steady-state protein abundance levels.

Molecular characterization of matrix metalloproteinase gene family across primates

Deregulation of matrix metalloproteinases (MMPs) contributes considerably to cancers, psychiatric disorders, macular degeneration and bone diseases. The use of humans in the development of MMPs as prognostic biomarkers and therapeutic targets is complicated by many factors, while primate models can be useful alternatives for this purpose. Here, we performed genome-enabled identification of putative MMPs across primate species, and comprehensively investigated the genes. Phylogenetic topology of the MMP family showed each type formulates a distinct clade, and was further clustered to classes, largely agreeing with classification based on biochemical properties and domain organization. Across primates, the excess of candidate sites of positive selection was detected for MMP-19, in addition to 1-3 sites in MMP-8, MMP-10 and MMP-26. MMP-26 showed Ka/Ks value above 1 between human and chimpanzee copies. We observed two copies of MMP-19 in the old-world monkey genomes, suggesting gene duplication at the early stage of or prior to the emergence of the lineage. Furin-activatable MMPs demonstrate the most variable properties regarding Domain organization and gene structure. During human aging, MMP-11 showed gradually decreased expression in testis, so as MMP-2, MMP-14, MMP15 and MMP-28 in ovary, while MMP-7 and MMP-21 showed elevated expression, implying their distinct roles in different reproductive organs. Co-expression clusters were formed among human MMPs both within and across classes, and expression correlation was observed in MMP genes across primates. Our results illuminate the utilization of MMPs for the discovery of prognostic biomarkers and therapeutic targets for aging-related diseases and carry new messages on MMP classification.

Comparative molecular genomic analyses of a spontaneous rhesus macaque model of mismatch repair-deficient colorectal cancer

Colorectal cancer (CRC) remains the third most common cancer in the US with 15% of cases displaying Microsatellite Instability (MSI) secondary to Lynch Syndrome (LS) or somatic hypermethylation of the MLH1 promoter. A cohort of rhesus macaques from our institution developed spontaneous mismatch repair deficient (MMRd) CRC with a notable fraction harboring a pathogenic germline mutation in MLH1 (c.1029C Collapse

Key Words Collapse

MESH Headings

• Animals
• Colonic Neoplasms/genetics
• Colorectal Neoplasms/genetics
• Colorectal Neoplasms/pathology
• Colorectal Neoplasms, Hereditary Nonpolyposis/genetics
• Colorectal Neoplasms, Hereditary Nonpolyposis/pathology
• DNA Methylation/genetics
• DNA Mismatch Repair/genetics
• Genomics
• Humans
• Macaca mulatta/genetics
• Mice
• Microsatellite Instability
• MutL Protein Homolog 1/genetics

Collapse

Grants

• P30 CA016672 NCI NIH HHS
• R24 OD011173 NIH HHS
• Feinberg Family Foundation
• MDACC Institutional Research Grant Program
• Cattleman for Cancer Research
• National Institutes of Health
• National Cancer Institute

Collapse

Affiliation(s)

Nejla Ozirmak Lermi Department of Clinical Cancer Prevention, The University of Texas MD Anderson Cancer Center, Houston, Texas, United States of America School of Health Professions, The University of Texas MD Anderson Cancer Center, Houston, Texas, United States of America
Stanton B. Gray Michale E. Keeling Center for Comparative Medicine and Research, The University of Texas MD Anderson Cancer Center, Houston, Texas, United States of America
Charles M. Bowen Department of Clinical Cancer Prevention, The University of Texas MD Anderson Cancer Center, Houston, Texas, United States of America
Laura Reyes-Uribe Department of Clinical Cancer Prevention, The University of Texas MD Anderson Cancer Center, Houston, Texas, United States of America
Beth K. Dray Charles River Laboratories, Ashland, Ohio, United States of America
Nan Deng Department of Clinical Cancer Prevention, The University of Texas MD Anderson Cancer Center, Houston, Texas, United States of America
R. Alan Harris Human Genome Sequencing Center and Department of Molecular and Human Genetics, Baylor College of Medicine, Houston, Texas, United States of America
Muthuswamy Raveendran Human Genome Sequencing Center and Department of Molecular and Human Genetics, Baylor College of Medicine, Houston, Texas, United States of America
Fernando Benavides Department of Epigenetics and Molecular Carcinogenesis, The University of Texas MD Anderson Cancer Center, Houston, Texas, United States of America
Carolyn L. Hodo Michale E. Keeling Center for Comparative Medicine and Research, The University of Texas MD Anderson Cancer Center, Houston, Texas, United States of America
Melissa W. Taggart Department of Pathology, The University of Texas MD Anderson Cancer Center, Houston, Texas, United States of America
Karen Colbert Maresso Department of Clinical Cancer Prevention, The University of Texas MD Anderson Cancer Center, Houston, Texas, United States of America
Krishna M. Sinha Department of Clinical Cancer Prevention, The University of Texas MD Anderson Cancer Center, Houston, Texas, United States of America
Jeffrey Rogers Human Genome Sequencing Center and Department of Molecular and Human Genetics, Baylor College of Medicine, Houston, Texas, United States of America
Eduardo Vilar Department of Clinical Cancer Prevention, The University of Texas MD Anderson Cancer Center, Houston, Texas, United States of America Clinical Cancer Genetics Program, The University of Texas MD Anderson Cancer Center, Houston, Texas, United States of America

Collapse

Landscape of adenosine-to-inosine RNA recoding across human tissues

RNA editing by adenosine deaminases changes the information encoded in the mRNA from its genomic blueprint. Editing of protein-coding sequences can introduce novel, functionally distinct, protein isoforms and diversify the proteome. The functional importance of a few recoding sites has been appreciated for decades. However, systematic methods to uncover these sites perform poorly, and the full repertoire of recoding in human and other mammals is unknown. Here we present a new detection approach, and analyze 9125 GTEx RNA-seq samples, to produce a highly-accurate atlas of 1517 editing sites within the coding region and their editing levels across human tissues. Single-cell RNA-seq data shows protein recoding contributes to the variability across cell subpopulations. Most highly edited sites are evolutionary conserved in non-primate mammals, attesting for adaptation. This comprehensive set can facilitate understanding of the role of recoding in human physiology and diseases.

Gabay et al. provide a highly-accurate atlas of recoding by A-to-I RNA editing in human, profiled across tissues and cell subpopulations. Most highly edited sites are evolutionary conserved in non-primate mammals, attesting for adaptation.

Genomic resources for rhesus macaques (Macaca mulatta)

Rhesus macaques (Macaca mulatta) are among the most extensively studied of nonhuman primates. This species has been the subject of many investigations concerning basic primate biology and behavior, including studies of social organization, developmental psychology, physiology, endocrinology, and neurodevelopment. Rhesus macaques are also critically important as a nonhuman primate model of human health and disease, including use in studies of infectious diseases, metabolic diseases, aging, and drug or alcohol abuse. Current research addressing fundamental biological and/or applied biomedical questions benefits from various genetic and genomic analyses. As a result, the genome of rhesus macaques has been the subject of more study than most nonhuman primates. This paper briefly discusses a number of information resources that can provide interested researchers with access to genetic and genomic data describing the content of the rhesus macaque genome, available information regarding genetic variation within the species, results from studies of gene expression, and other aspects of genomic analysis. Specific online databases are discussed, including the US National Center for Biotechnology Information, the University of California Santa Cruz genome browser, Ensembl genome browser, the Macaque Genotype and Phenotype database (mGAP), Rhesusbase, and others.

Characterization of the Testis-specific LINC01016 Gene Reveals Isoform-specific Roles in Controlling Biological Processes

Long noncoding RNAs (lncRNAs) have emerged as critical regulators of biological processes. However, the aberrant expression of an isoform from the same lncRNA gene could lead to RNA with altered functions due to changes in their conformations, leading to diseases. Here, we describe a detailed characterization of the gene that encodes long intergenic non-protein-coding RNA 01016 (LINC01016, also known as LncRNA1195) with a focus on its structure, exon usage, and expression in human and macaque tissues. In this study we show that it is among the highly expressed lncRNAs in the testis, exclusively conserved among nonhuman primates, suggesting its recent evolution and is processed into 12 distinct RNAs in testis, cervix, and uterus tissues. Further, we integrate de novo annotation of expressed LINC01016 transcripts and isoform-dependent gene expression analyses to show that human LINC01016 is a multiexon gene, processed through differential exon usage with isoform-specific roles. Furthermore, in cervical, testicular, and uterine cancers, LINC01016 isoforms are differentially expressed, and their expression is predictive of survival in these cancers. This study has revealed an essential aspect of lncRNA biology, rarely associated with coding RNAs, that lncRNA genes are precisely processed to generate isoforms with distinct biological roles in specific tissues.

Evolutionary dynamics of circular RNAs in primates

Many primate genes produce circular RNAs (circRNAs). However, the extent of circRNA conservation between closely related species remains unclear. By comparing tissue-specific transcriptomes across over 70 million years of primate evolution, we identify that within 3 million years circRNA expression profiles diverged such that they are more related to species identity than organ type. However, our analysis also revealed a subset of circRNAs with conserved neural expression across tens of millions of years of evolution. By comparing to species-specific circRNAs, we identified that the downstream intron of the conserved circRNAs display a dramatic lengthening during evolution due to the insertion of novel retrotransposons. Our work provides comparative analyses of the mechanisms promoting circRNAs to generate increased transcriptomic complexity in primates.

Elevated metallothionein expression in long-lived species mediates the influence of cadmium accumulation on aging

Cadmium (Cd) accumulates with aging and is elevated in long-lived species. Metallothioneins (MTs), small cysteine-rich proteins involved in metal homeostasis and Cd detoxification, are known to be related to longevity. However, the relationship between Cd accumulation, the role of MTs, and aging is currently unclear. Specifically, we do not know if long-lived species evolved an efficient metal stress response by upregulating their MT levels to reduce the toxic effects of environmental pollutants, such as Cd, that accumulate over their longer life span. It is also unknown if the number of MT genes, their expression, or both protect the organisms from potentially damaging effects during aging. To address these questions, we reanalyzed several cross-species studies and obtained data on MT expression and Cd accumulation in long-lived mouse models. We confirmed a relationship between species maximum life span in captive mammals and their Cd content in liver and kidney. We found that although the number of MT genes does not affect longevity, gene expression and protein amount of specific MT paralogs are strongly related to life span in mammals. MT expression rather than gene number may influence the high Cd levels and longevity of some species. In support of this, we found that overexpression of MT-1 accelerated Cd accumulation in mice and that tissue Cd was higher in long-lived mouse strains with high MT expression. We conclude that long-lived species have evolved a more efficient stress response by upregulating the expression of MT genes in presence of Cd, which contributes to elevated tissue Cd levels.

REDIportal: millions of novel A-to-I RNA editing events from thousands of RNAseq experiments

RNA editing is a relevant epitranscriptome phenomenon able to increase the transcriptome and proteome diversity of eukaryotic organisms. ADAR mediated RNA editing is widespread in humans in which millions of A-to-I changes modify thousands of primary transcripts. RNA editing has pivotal roles in the regulation of gene expression or modulation of the innate immune response or functioning of several neurotransmitter receptors. Massive transcriptome sequencing has fostered the research in this field. Nonetheless, different aspects of the RNA editing biology are still unknown and need to be elucidated. To support the study of A-to-I RNA editing we have updated our REDIportal catalogue raising its content to about 16 millions of events detected in 9642 human RNAseq samples from the GTEx project by using a dedicated pipeline based on the HPC version of the REDItools software. REDIportal now allows searches at sample level, provides overviews of RNA editing profiles per each RNAseq experiment, implements a Gene View module to look at individual events in their genic context and hosts the CLAIRE database. Starting from this novel version, REDIportal will start collecting non-human RNA editing changes for comparative genomics investigations. The database is freely available at http://srv00.recas.ba.infn.it/atlas/index.html.

LRF/ZBTB7A conservation accentuates its potential as a therapeutic target for the hematopoietic disorders

Conserved sequences across species have always provided valuable insights to improve our understanding on the human genome's entity and the interplay among different loci. Lymphoma/leukemia related factor (LRF) is encoded by ZBTB7A gene and belongs to an evolutionarily conserved family of transcription factors, implicated in vital cellular functions. The present data, demonstrating the wide-spread and the high overlap of the LRF/ZBTB7A recognition sites with genomic segments identified as CpG islands in the human genome, suggest that its binding capacity strongly depends on a specific sequence-encoded feature within CpGs. We have previously shown that de-methylation of the CpG island 326 lying in the ZBTB7A gene promoter is associated with impaired pharmacological induction of fetal hemoglobin in β-type hemoglobinopathies patients. Within this context we aimed to investigate the extent of the LRF/ZBTB7A conservation among primates and mouse genome, focusing our interest also on the CpG island flanking the gene's promoter region, in an effort to further establish its epigenetic regulatory role in human hematopoiesis and pharmacological involvement in hematopoietic disorders. Comparative analysis of the human ZBTB7A nucleotide and amino acid sequences and orthologous sequences among non-human primates and mouse, exhibited high conservation scores. Pathway analysis, clearly indicated that LRF/ZBTB7A influences conserved cellular processes. These data in conjunction with the high levels of expression foremost in hematopoietic tissues, highlighted LRF/ZBTB7A as an essential factor operating indisputably during hematopoiesis.

An improved de novo genome assembly of the common marmoset genome yields improved contiguity and increased mapping rates of sequence data

BACKGROUND

The common marmoset (Callithrix jacchus) is one of the most studied primate model organisms. However, the marmoset genomes available in the public databases are highly fragmented and filled with sequence gaps, hindering research advances related to marmoset genomics and transcriptomics.

RESULTS

Here we utilize single-molecule, long-read sequence data to improve and update the existing genome assembly and report a near-complete genome of the common marmoset. The assembly is of 2.79 Gb size, with a contig N50 length of 6.37 Mb and a chromosomal scaffold N50 length of 143.91 Mb, representing the most contiguous and high-quality marmoset genome up to date. Approximately 90% of the assembled genome was represented in contigs longer than 1 Mb, with approximately 104-fold improvement in contiguity over the previously published marmoset genome. More than 98% of the gaps from the previously published genomes were filled successfully, which improved the mapping rates of genomic and transcriptomic data on to the assembled genome.

CONCLUSIONS

Altogether the updated, high-quality common marmoset genome assembly provide improvements at various levels over the previous versions of the marmoset genome assemblies. This will allow researchers working on primate genomics to apply the genome more efficiently for their genomic and transcriptomic sequence data.

Resources for functional genomic studies of health and development in nonhuman primates

Primates display a wide range of phenotypic variation underlaid by complex genetically regulated mechanisms. The links among DNA sequence, gene function, and phenotype have been of interest from an evolutionary perspective, to understand functional genome evolution and its phenotypic consequences, and from a biomedical perspective to understand the shared and human-specific roots of health and disease. Progress in methods for characterizing genetic, transcriptomic, and DNA methylation (DNAm) variation is driving the rapid development of extensive omics resources, which are now increasingly available from humans as well as a growing number of nonhuman primates (NHPs). The fast growth of large-scale genomic data is driving the emergence of integrated tools and databases, thus facilitating studies of gene functionality across primates. This review describes NHP genomic resources that can aid in exploration of how genes shape primate phenotypes. It focuses on the gene expression trajectories across development in different tissues, the identification of functional genetic variation (including variants deleterious for protein function and regulatory variants modulating gene expression), and DNAm profiles as an emerging tool to understand the process of aging. These resources enable comparative functional genomics approaches to identify species-specific and primate-shared gene functionalities associated with health and development.

Transcriptomic Analysis of Cell-free Fetal RNA in the Amniotic Fluid of Vervet Monkeys (Chlorocebus sabaeus)

NHP are important translational models for understanding the genomic underpinnings of growth, development, fetal programming, and predisposition to disease, with potential for the development of early health biomarkers. Understanding how prenatal gene expression is linked to pre- and postnatal health and development requires methods for assessing the fetal transcriptome. Here we used RNAseq methodology to analyze the expression of cell-free fetal RNA in the amniotic fluid supernatant (AFS) of vervet monkeys. Despite the naturally high level of degradation of free-floating RNA, we detected more than 10,000 gene transcripts in vervet AFS. The most highly expressed genes were H19, IGF2, and TPT1, which are involved in embryonic growth and glycemic health. We noted global similarities in expression profiles between vervets and humans, with genes involved in embryonic growth and glycemic health among the genes most highly expressed in AFS. Our study demonstrates both the feasibility and usefulness of prenatal transcriptomic profiles, by using amniocentesis procedures to obtain AFS and cell-free fetal RNA from pregnant vervets.

NHP-immunome: A translational research-oriented database of non-human primate immune system proteins

We are currently living the advent of a new age for medicine in which basic research is being quickly translated into marketable drugs, and the widespread access to genomics data is allowing the design and implementation of personalized solutions to medical conditions. Non-human primates (NHP) have gained an essential role in drug discovery and safety testing due to their close phylogenetic relationship to humans. In this study, a collection of well characterized genes of the human immune system was used to define the orthology-based immunome in four NHP species, with carefully curated annotations available based on multi-tissue RNA-seq datasets. A broad variation in the frequency of expressed protein isoforms was observed between species. Finally, this analysis also revealed the lack of expression of at least four different chemokines in new-world primates. In addition, transcripts corresponding to four genes including interleukin 12 subunit alpha were expressed in humans but no other primate species analyzed. Access to the non-human primate immunome is available in http://www.fidic.org.co:90/proyecto/.

Conservation of cell-intrinsic immune responses in diverse nonhuman primate species

The transcriptomic response of diverse nonhuman primate (NHP) species to poly(I:C) is highly conserved, and this novel RNA sequencing dataset will help improve NHP genome annotations.

Differences in immune responses across species can contribute to the varying permissivity of species to the same viral pathogen. Understanding how our closest evolutionary relatives, nonhuman primates (NHPs), confront pathogens and how these responses have evolved over time could shed light on host range barriers, especially for zoonotic infections. Here, we analyzed cell-intrinsic immunity of primary cells from the broadest panel of NHP species interrogated to date, including humans, great apes, and Old and New World monkeys. Our analysis of their transcriptomes after poly(I:C) transfection revealed conservation in the functional consequences of their response. In mapping reads to either the human or the species-specific genomes, we observed that with the current state of NHP annotations, the percent of reads assigned to a genetic feature was largely similar regardless of the method. Together, these data provide a baseline for the cell-intrinsic responses elicited by a potent immune stimulus across multiple NHP donors, including endangered species, and serve as a resource for refining and furthering the existing annotations of NHP genomes.

Absence of a human ortholog of rodent Kupffer cell galactose-binding receptor encoded by the CLEC4f gene

The murine CLEC4f gene encodes the Kupffer cell receptor, a galactose-binding receptor containing a C-type carbohydrate-recognition domain. Orthologs have been identified in nearly 100 species. The receptors from rat and mouse have previously been characterized and data presented here show that functional CLEC4f protein is expressed in domestic cattle (Bos taurus). However, the human CLEC4f gene does not encode a functional receptor because a mutation in the splice acceptor site of the final exon prevents appropriate splicing and a missense mutation disrupts the sugar-binding site. Transcriptomic and PCR analysis of transcripts confirms the absence of a spliced transcript containing the final exon and only background levels of transcripts are detected in human tissues. These mutations are also present in the CLEC4f gene in Neanderthals. In contrast to humans, closely related species, including chimpanzees, do have CLEC4f genes that encode full-length receptors. Affinity chromatography and glycan array results demonstrate that the chimpanzee, bovine and murine proteins all bind to galactose, but they show preferences for different subsets of galactose-containing glycans. In non-human primates, the receptor is expressed in spleen rather than in liver. The results indicate that the CLEC4f protein probably has distinct functions in different species. Absence of the receptor precludes using it for targeting of glycoconjugates to cells in human liver. The fact that CLEC4f protein is expressed in spleen in non-human primates and the close evolutionary relationship of the CLEC4f protein to langerin (CD207) suggest that it may function in the immune system, possibly as a pathogen receptor.

Epigenetic Modifications in Acute Myeloid Leukemia: Prognosis, Treatment, and Heterogeneity

Leukemia, specifically acute myeloid leukemia (AML), is a common malignancy that can be differentiated into multiple subtypes based on leukemogenic history and etiology. Although genetic aberrations, particularly cytogenetic abnormalities and mutations in known oncogenes, play an integral role in AML development, epigenetic processes have been shown as a significant and sometimes independent dynamic in AML pathophysiology. Here, we summarize how tumors evolve and describe AML through an epigenetic lens, including discussions on recent discoveries that include prognostics from epialleles, changes in RNA function for hematopoietic stem cells and the epitranscriptome, and novel epigenetic treatment options. We further describe the limitations of treatment in the context of the high degree of heterogeneity that characterizes acute myeloid leukemia.

Combinatorial Genetics Reveals a Scaling Law for the Effects of Mutations on Splicing

Despite a wealth of molecular knowledge, quantitative laws for accurate prediction of biological phenomena remain rare. Alternative pre-mRNA splicing is an important regulated step in gene expression frequently perturbed in human disease. To understand the combined effects of mutations during evolution, we quantified the effects of all possible combinations of exonic mutations accumulated during the emergence of an alternatively spliced human exon. This revealed that mutation effects scale non-monotonically with the inclusion level of an exon, with each mutation having maximum effect at a predictable intermediate inclusion level. This scaling is observed genome-wide for cis and trans perturbations of splicing, including for natural and disease-associated variants. Mathematical modeling suggests that competition between alternative splice sites is sufficient to cause this non-linearity in the genotype-phenotype map. Combining the global scaling law with specific pairwise interactions between neighboring mutations allows accurate prediction of the effects of complex genotype changes involving >10 mutations.

The conservation landscape of the human ribosomal RNA gene repeats

Ribosomal RNA gene repeats (rDNA) encode ribosomal RNA, a major component of ribosomes. Ribosome biogenesis is central to cellular metabolic regulation, and several diseases are associated with rDNA dysfunction, notably cancer, However, its highly repetitive nature has severely limited characterization of the elements responsible for rDNA function. Here we make use of phylogenetic footprinting to provide a comprehensive list of novel, potentially functional elements in the human rDNA. Complete rDNA sequences for six non-human primate species were constructed using de novo whole genome assemblies. These new sequences were used to determine the conservation profile of the human rDNA, revealing 49 conserved regions in the rDNA intergenic spacer (IGS). To provide insights into the potential roles of these conserved regions, the conservation profile was integrated with functional genomics datasets. We find two major zones that contain conserved elements characterised by enrichment of transcription-associated chromatin factors, and transcription. Conservation of some IGS transcripts in the apes underpins the potential functional significance of these transcripts and the elements controlling their expression. Our results characterize the conservation landscape of the human IGS and suggest that noncoding transcription and chromatin elements are conserved and important features of this unique genomic region.

Improved Macaca fascicularis gene annotation reveals evolution of gene expression profiles in multiple tissues

Backgrounds

Macaca fascicularis (M. fascicularis) is a primate model organism that played important role in studying human health. It is vital to better understand the similarity and differences of gene regulation between M. fascicularis and human. Current comparative study of gene regulation between the two species are limited by low quality of gene annotation and lack of regulatory element data on M. fascicularis genome.

Results

In this study, we improved the M. fascicularis gene annotation with 57 gene expression data from multiple tissues and, more importantly, a manual curation procedure. The new annotation enabled us to map gene expression and identify gene location more accurately.

Conclusions

Comparing with human gene expression data from the same cell types, we characterized the evolution of expression patterns of homologous genes.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-5183-y) contains supplementary material, which is available to authorized users.

A Comparative Assessment of Human and Chimpanzee iPSC-derived Cardiomyocytes with Primary Heart Tissues

Comparative genomic studies in primates have the potential to reveal the genetic and mechanistic basis for human specific traits. These studies may also help us better understand inter-species phenotypic differences that are clinically relevant. Unfortunately, the obvious limitation on sample collection and experimentation in humans and non-human apes severely restrict our ability to perform dynamic comparative studies in primates. Induced pluripotent stem cells (iPSCs), and their corresponding differentiated cells, may provide a suitable alternative system for dynamic comparative studies. Yet, to effectively use iPSCs and differentiated cells for comparative studies, one must characterize the extent to which these systems faithfully represent biological processes in primary tissues. To do so, we compared gene expression data from primary adult heart tissue and iPSC-derived cardiomyocytes from multiple human and chimpanzee individuals. We determined that gene expression in cultured cardiomyocytes from both human and chimpanzee is most similar to that of adult hearts compared to other adult tissues. Using a comparative framework, we found that 50% of gene regulatory differences between human and chimpanzee hearts are also observed between species in cultured cardiomyocytes; conversely, inter-species regulatory differences seen in cardiomyocytes are found significantly more often in hearts than in other primary tissues. Our work provides a detailed description of the utility and limitation of differentiated cardiomyocytes as a system for comparative functional genomic studies in primates.

Similarity and variation in the insulin-like growth factor 2 - H19 locus in primates

Insulin-like growth factor 2 (IGF2), a small, secreted protein, is critical for fetal and prenatal growth in humans and other mammals. The IGF2 gene and its mouse homolog comprise part of a conserved linkage group that is regulated by parental imprinting, with IGF2/ Igf2 being expressed from the paternal chromosome, and the adjacent H19 gene from the maternal chromosome. By using information extracted from public genomic and gene expression databases, I have now analyzed this locus in nine nonhuman primate species representing over 60 million years of evolutionary divergence from a common progenitor. Both IGF2 and H19 genes and the entire locus have been conserved among these primates. Each primate IGF2 gene except for gibbon and marmoset is composed of 10 exons and contains five potential promoters, each with distinctive 5'-untranslated exons. Similarly, except for marmoset and mouse lemur, H19 consists of six exons and has two promoters. DNA sequence conservation is high, not only in orthologous exons and promoters, but also in a putative imprinting control region located 5' to H19 and in multiple potential distal enhancer elements found 3' to H19. Collectively, these results support the hypothesis that common regulatory processes shaped the IGF2 - H19 locus before the onset of primate speciation more than 85 million years ago. This study also leads to the conclusion that inaccuracies in data presentation in genetic repositories could limit our ability to develop novel insights about roles of individual genes and multigene loci in mammalian physiology and disease.

Sooty mangabey genome sequence provides insight into AIDS resistance in a natural SIV host

In contrast to infections with human immunodeficiency virus (HIV) in humans and simian immunodeficiency virus (SIV) in macaques, SIV infection of a natural host, sooty mangabeys (Cercocebus atys), is non-pathogenic despite high viraemia¹. Here we sequenced and assembled the genome of a captive sooty mangabey. We conducted genome-wide comparative analyses of transcript assemblies from C. atys and AIDS-susceptible species, such as humans and macaques, to identify candidates for host genetic factors that influence susceptibility. We identified several immune-related genes in the genome of C. atys that show substantial sequence divergence from macaques or humans. One of these sequence divergences, a C-terminal frameshift in the toll-like receptor-4 (TLR4) gene of C. atys, is associated with a blunted in vitro response to TLR-4 ligands. In addition, we found a major structural change in exons 3–4 of the immune-regulatory protein intercellular adhesion molecule 2 (ICAM-2); expression of this variant leads to reduced cell surface expression of ICAM-2. These data provide a resource for comparative genomic studies of HIV and/or SIV pathogenesis and may help to elucidate the mechanisms by which SIV-infected sooty mangabeys avoid AIDS.

Caloric Restriction Engages Hepatic RNA Processing Mechanisms in Rhesus Monkeys

Caloric restriction (CR) extends lifespan and delays the onset of age-related disorders in diverse species. Metabolic regulatory pathways have been implicated in the mechanisms of CR, but the molecular details have not been elucidated. Here, we show that CR engages RNA processing of genes associated with a highly integrated reprogramming of hepatic metabolism. We conducted molecular profiling of liver biopsies collected from adult male rhesus monkeys (Macaca mulatta) at baseline and after 2 years on control or CR (30% restricted) diet. Quantitation of over 20,000 molecules from the hepatic transcriptome, proteome, and metabolome indicated that metabolism and RNA processing are major features of the response to CR. Predictive models identified lipid, branched-chain amino acid, and short-chain carbon metabolic pathways, with alternate transcript use for over half of the genes in the CR network. We conclude that RNA-based mechanisms are central to the CR response and integral in metabolic reprogramming.

Identify Down syndrome transcriptome associations using integrative analysis of microarray database and correlation-interaction network

BACKGROUND

Long non-coding RNAs (lncRNAs) have previously been emerged as key players in a series of biological processes. Dysregulation of lncRNA is correlated to human diseases including neurological disorders. Here, we developed a multi-step bioinformatics analysis to study the functions of a particular Down syndrome-associated gene DSCR9 including the lncRNAs. The method is named correlation-interaction-network (COIN), based on which a pipeline is implemented. Co-expression gene network analysis and biological network analysis results are presented.

METHODS

We identified the regulation function of DSCR9, a lncRNA transcribed from the Down syndrome critical region (DSCR) of chromosome 21, by analyzing its co-expression genes from over 1700 sets and nearly 60,000 public Affymetrix human U133-Plus 2 transcriptional profiling microarrays. After proper evaluations, a threshold is chosen to filter the data and get satisfactory results. Microarray data resource is from EBI database and protein-protein interaction (PPI) network information is incorporated from the most complete network databases. PPI integration strategy guarantees complete information regarding DSCR9. Enrichment analysis is performed to identify significantly correlated pathways.

RESULTS

We found that the most significant pathways associated with the top DSCR9 co-expressed genes were shown to be involved in neuro-active ligand-receptor interaction (GLP1R, HTR4, P2RX2, UCN3, and UTS2R), calcium signaling pathway (CACNA1F, CACNG4, HTR4, P2RX2, and SLC8A3), neuronal system (KCNJ5 and SYN1) by the KEGG, and GO analysis. The A549 and U251 cell lines with stable DSCR9 overexpression were constructed. We validated 10 DSCR9 co-expression genes by qPCR in both cell lines with over 70% accuracy.

CONCLUSIONS

DSCR9 was highly correlated with genes that were known as important factors in the developments and functions of nervous system, indicating that DSCR9 may regulate neurological proteins regarding Down syndrome and other neurological-related diseases. The pipeline can be properly adjusted to other applications.

Transcription start site profiling of 15 anatomical regions of the Macaca mulatta central nervous system

Rhesus macaque was the second non-human primate whose genome has been fully sequenced and is one of the most used model organisms to study human biology and disease, thanks to the close evolutionary relationship between the two species. But compared to human, where several previously unknown RNAs have been uncovered, the macaque transcriptome is less studied. Publicly available RNA expression resources for macaque are limited, even for brain, which is highly relevant to study human cognitive abilities. In an effort to complement those resources, FANTOM5 profiled 15 distinct anatomical regions of the aged macaque central nervous system using Cap Analysis of Gene Expression, a high-resolution, annotation-independent technology that allows monitoring of transcription initiation events with high accuracy. We identified 25,869 CAGE peaks, representing bona fide promoters. For each peak we provide detailed annotation, expanding the landscape of ‘known’ macaque genes, and we show concrete examples on how to use the resulting data. We believe this data represents a useful resource to understand the central nervous system in macaque.

Identification and characterization of a subtelomeric satellite DNA in Callitrichini monkeys

Repetitive DNAs are abundant fast-evolving components of eukaryotic genomes, which often possess important structural and functional roles. Despite their ubiquity, repetitive DNAs are poorly studied when compared with the genic fraction of genomes. Here, we took advantage of the availability of the sequenced genome of the common marmoset Callithrix jacchus to assess its satellite DNAs (satDNAs) and their distribution in Callitrichini. After clustering analysis of all reads and comparisons by similarity, we identified a satDNA composed by 171 bp motifs, named MarmoSAT, which composes 1.09% of the C. jacchus genome. Fluorescent in situ hybridization on chromosomes of species from the genera Callithrix, Mico and Callimico showed that MarmoSAT had a subtelomeric location. In addition to the common monomeric, we found that MarmoSAT was also organized in higher-order repeats of 338 bp in Callimico goeldii. Our phylogenetic analyses showed that MarmoSAT repeats from C. jacchus lack chromosome-specific features, suggesting exchange events among subterminal regions of non-homologous chromosomes. MarmoSAT is transcribed in several tissues of C. jacchus, with the highest transcription levels in spleen, thymus and heart. The transcription profile and subtelomeric location suggest that MarmoSAT may be involved in the regulation of telomerase and modulation of telomeric chromatin.

Elucidating the Role of Host Long Non-Coding RNA during Viral Infection: Challenges and Paths Forward

Research over the past decade has clearly shown that long non-coding RNAs (lncRNAs) are functional. Many lncRNAs can be related to immunity and the host response to viral infection, but their specific functions remain largely elusive. The vast majority of lncRNAs are annotated with extremely limited knowledge and tend to be expressed at low levels, making ad hoc experimentation difficult. Changes to lncRNA expression during infection can be systematically profiled using deep sequencing; however, this often produces an intractable number of candidate lncRNAs, leaving no clear path forward. For these reasons, it is especially important to prioritize lncRNAs into high-confidence “hits” by utilizing multiple methodologies. Large scale perturbation studies may be used to screen lncRNAs involved in phenotypes of interest, such as resistance to viral infection. Single cell transcriptome sequencing quantifies cell-type specific lncRNAs that are less abundant in a mixture. When coupled with iterative experimental validations, new computational strategies for efficiently integrating orthogonal high-throughput data will likely be the driver for elucidating the functional role of lncRNAs during viral infection. This review highlights new high-throughput technologies and discusses the potential for integrative computational analysis to streamline the identification of infection-related lncRNAs and unveil novel targets for antiviral therapeutics.

Host Transcriptional Response to Ebola Virus Infection

Ebola virus disease (EVD) is a serious illness that causes severe disease in humans and non-human primates (NHPs) and has mortality rates up to 90%. EVD is caused by the Ebolavirus and currently there are no licensed therapeutics or vaccines to treat EVD. Due to its high mortality rates and potential as a bioterrorist weapon, a better understanding of the disease is of high priority. Multiparametric analysis techniques allow for a more complete understanding of a disease and the host response. Analysis of RNA species present in a sample can lead to a greater understanding of activation or suppression of different states of the immune response. Transcriptomic analyses such as microarrays and RNA-Sequencing (RNA-Seq) have been important tools to better understand the global gene expression response to EVD. In this review, we outline the current knowledge gained by transcriptomic analysis of EVD.

RNA sequencing (RNA-Seq) of lymph node, spleen, and thymus transcriptome from wild Peninsular Malaysian cynomolgus macaque (Macaca fascicularis)

The cynomolgus macaque (Macaca fascicularis) is an extensively utilised nonhuman primate model for biomedical research due to its biological, behavioural, and genetic similarities to humans. Genomic information of cynomolgus macaque is vital for research in various fields; however, there is presently a shortage of genomic information on the Malaysian cynomolgus macaque. This study aimed to sequence, assemble, annotate, and profile the Peninsular Malaysian cynomolgus macaque transcriptome derived from three tissues (lymph node, spleen, and thymus) using RNA sequencing (RNA-Seq) technology. A total of 174,208,078 paired end 70 base pair sequencing reads were obtained from the Illumina Hi-Seq 2500 sequencer. The overall mapping percentage of the sequencing reads to the M. fascicularis reference genome ranged from 53–63%. Categorisation of expressed genes to Gene Ontology (GO) and KEGG pathway categories revealed that GO terms with the highest number of associated expressed genes include Cellular process, Catalytic activity, and Cell part, while for pathway categorisation, the majority of expressed genes in lymph node, spleen, and thymus fall under the Global overview and maps pathway category, while 266, 221, and 138 genes from lymph node, spleen, and thymus were respectively enriched in the Immune system category. Enriched Immune system pathways include Platelet activation pathway, Antigen processing and presentation, B cell receptor signalling pathway, and Intestinal immune network for IgA production. Differential gene expression analysis among the three tissues revealed 574 differentially expressed genes (DEG) between lymph and spleen, 5402 DEGs between lymph and thymus, and 7008 DEGs between spleen and thymus. Venn diagram analysis of expressed genes revealed a total of 2,630, 253, and 279 tissue-specific genes respectively for lymph node, spleen, and thymus tissues. This is the first time the lymph node, spleen, and thymus transcriptome of the Peninsular Malaysian cynomolgus macaque have been sequenced via RNA-Seq. Novel transcriptomic data will further enrich the present M. fascicularis genomic database and provide future research potentials, including novel transcript discovery, comparative studies, and molecular markers development.

Mitochondrial thioredoxin reductase 2 is elevated in long-lived primate as well as rodent species and extends fly mean lifespan

In a survey of enzymes related to protein oxidation and cellular redox state, we found activity of the redox enzyme thioredoxin reductase (TXNRD) to be elevated in cells from long‐lived species of rodents, primates, and birds. Elevated TXNRD activity in long‐lived species reflected increases in the mitochondrial form, TXNRD2, rather than the cytosolic forms TXNRD1 and TXNRD3. Analysis of published RNA‐Seq data showed elevated TXNRD2 mRNA in multiple organs of longer‐lived primates, suggesting that the phenomenon is not limited to skin‐derived fibroblasts. Elevation of TXNRD2 activity and protein levels was also noted in liver of three different long‐lived mutant mice, and in normal male mice treated with a drug that extends lifespan in males. Overexpression of mitochondrial TXNRD2 in Drosophila melanogaster extended median (but not maximum) lifespan in female flies with a small lifespan extension in males; in contrast, overexpression of the cytosolic form, TXNRD1, did not produce a lifespan extension.

Evolutionary acquisition of promoter-associated non-coding RNA (pancRNA) repertoires diversifies species-dependent gene activation mechanisms in mammals

Background

Recent transcriptome analyses have shown that long non-coding RNAs (ncRNAs) play extensive roles in transcriptional regulation. In particular, we have reported that promoter-associated ncRNAs (pancRNAs) activate the partner gene expression via local epigenetic changes.

Results

Here, we identify thousands of genes under pancRNA-mediated transcriptional activation in five mammalian species in common. In the mouse, 1) pancRNA-partnered genes confined their expression pattern to certain tissues compared to pancRNA-lacking genes, 2) expression of pancRNAs was significantly correlated with the enrichment of active chromatin marks, H3K4 trimethylation and H3K27 acetylation, at the promoter regions of the partner genes, 3) H3K4me1 marked the pancRNA-partnered genes regardless of their expression level, and 4) C- or G-skewed motifs were exclusively overrepresented between−200 and−1 bp relative to the transcription start sites of the pancRNA-partnered genes. More importantly, the comparative transcriptome analysis among five different mammalian species using a total of 25 counterpart tissues showed that the overall pancRNA expression profile exhibited extremely high species-specificity compared to that of total mRNA, suggesting that interspecies difference in pancRNA repertoires might lead to the diversification of mRNA expression profiles.

Conclusions

The present study raises the interesting possibility that the gain and/or loss of gene-activation-associated pancRNA repertoires, caused by formation or disruption of the genomic GC-skewed structure in the course of evolution, finely shape the tissue-specific pattern of gene expression according to a given species.

Electronic supplementary material

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

The value of new genome references

Genomic information has become a ubiquitous and almost essential aspect of biological research. Over the last 10-15 years, the cost of generating sequence data from DNA or RNA samples has dramatically declined and our ability to interpret those data increased just as remarkably. Although it is still possible for biologists to conduct interesting and valuable research on species for which genomic data are not available, the impact of having access to a high quality whole genome reference assembly for a given species is nothing short of transformational. Research on a species for which we have no DNA or RNA sequence data is restricted in fundamental ways. In contrast, even access to an initial draft quality genome (see below for definitions) opens a wide range of opportunities that are simply not available without that reference genome assembly. Although a complete discussion of the impact of genome sequencing and assembly is beyond the scope of this short paper, the goal of this review is to summarize the most common and highest impact contributions that whole genome sequencing and assembly has had on comparative and evolutionary biology.

Inter-platform concordance of gene expression data for the prediction of chemical mode of action

BACKGROUND

It is interesting to study the consistency of outcomes arising from two genomic platforms: Microarray and RNAseq, which are established on fundamentally different technologies. This topic has been frequently discussed from the prospect of comparing differentially expressed genes (DEGs). In this study, we explore the inter-platform concordance between microarray and RNASeq in their ability to classify samples based on genomic information. We use a set of 7 standard multi-class classifiers and an adaptive ensemble classifier developed around them to predict Chemical Modes of Actions (MOA) of data profiled by microarray and RNASeq platforms from Rat Liver samples exposed to a variety of chemical compounds. We study the concordance between microarray and RNASeq data in various forms, based on classifier's performance between two platforms.

RESULTS

Using an ensemble classifier we observe improved prediction performance compared to a set of standard classifiers. We discover a clear concordance between each individual classifier's performances in two genomic platforms. Additionally, we identify a set of important genes those specifies MOAs, by focusing on their impact on the classification and later we find that some of these top genes have direct associations with the presence of toxic compounds in the liver.

CONCLUSION

Overall there appears to be fair amount of concordance between the two platforms as far as classification is concerned. We observe widely different classification performances among individual classifiers, which reflect the unreliability of restricting to a single classifier in the case of high dimensional classification problems.

REVIEWERS

An extended abstract of this research paper was selected for the CAMDA Satellite Meeting to ISMB 2015 by the CAMDA Programme Committee. The full research paper then underwent two rounds of Open Peer Review under a responsible CAMDA Programme Committee member, Lan Hu, PhD (Bio-Rad Laboratories, Digital Biology Center-Cambridge). Open Peer Review was provided by Yiyi Liu and Partha Dey. The Reviewer Comments section shows the full reviews and author responses.

A human microprotein that interacts with the mRNA decapping complex

Proteomic detection of non-annotated microproteins indicates the translation of hundreds of small open reading frames (smORFs) in human cells, but whether these microproteins are functional or not is unknown. Here, we report the discovery and characterization of a 7-kDa human microprotein we named non-annotated P-body dissociating polypeptide (NoBody). NoBody interacts with mRNA decapping proteins, which remove the 5' cap from mRNAs to promote 5'-to-3' decay. Decapping proteins participate in mRNA turnover and nonsense-mediated decay (NMD). NoBody localizes to mRNA-decay-associated RNA-protein granules called P-bodies. Modulation of NoBody levels reveals that its abundance is anticorrelated with cellular P-body numbers and alters the steady-state levels of a cellular NMD substrate. These results implicate NoBody as a novel component of the mRNA decapping complex and demonstrate potential functionality of a newly discovered microprotein.

The beagle dog MicroRNA tissue atlas: identifying translatable biomarkers of organ toxicity

Background

MicroRNAs (miRNA) are varied in length, under 25 nucleotides, single-stranded noncoding RNA that regulate post-transcriptional gene expression via translational repression or mRNA degradation. Elevated levels of miRNAs can be detected in systemic circulation after tissue injury, suggesting that miRNAs are released following cellular damage. Because of their remarkable stability, ease of detection in biofluids, and tissue specific expression patterns, miRNAs have the potential to be specific biomarkers of organ injury. The identification of miRNA biomarkers requires a systematic approach: 1) determine the miRNA tissue expression profiles within a mammalian species via next generation sequencing; 2) identify enriched and/or specific miRNA expression within organs of toxicologic interest, and 3) in vivo validation with tissue-specific toxicants. While miRNA tissue expression has been reported in rodents and humans, little data exists on miRNA tissue expression in the dog, a relevant toxicology species. The generation and evaluation of the first dog miRNA tissue atlas is described here.

Results

Analysis of 16 tissues from five male beagle dogs identified 106 tissue enriched miRNAs, 60 of which were highly enriched in a single organ, and thus may serve as biomarkers of organ injury. A proof of concept study in dogs dosed with hepatotoxicants evaluated a qPCR panel of 15 tissue enriched miRNAs specific to liver, heart, skeletal muscle, pancreas, testes, and brain. Dogs with elevated serum levels of miR-122 and miR-885 had a correlative increase of alanine aminotransferase, and microscopic analysis confirmed liver damage. Other non-liver enriched miRNAs included in the screening panel were unaffected. Eli Lilly authors created a complimentary Sprague Dawely rat miRNA tissue atlas and demonstrated increased pancreas enriched miRNA levels in circulation, following caerulein administration in rat and dog.

Conclusion

The dog miRNA tissue atlas provides a resource for biomarker discovery and can be further mined with refinement of dog genome annotation. The 60 highly enriched tissue miRNAs identified within the dog miRNA tissue atlas could serve as diagnostic biomarkers and will require further validation by in vivo correlation to histopathology. Once validated, these tissue enriched miRNAs could be combined into a powerful qPCR screening panel to identify organ toxicity during early drug development.

Electronic supplementary material

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

Identification of Tissue-Specific Protein-Coding and Noncoding Transcripts across 14 Human Tissues Using RNA-seq

Many diseases and adverse drug reactions exhibit tissue specificity. To better understand the tissue-specific expression characteristics of transcripts in different human tissues, we deeply sequenced RNA samples from 14 different human tissues. After filtering many lowly expressed transcripts, 24,729 protein-coding transcripts and 1,653 noncoding transcripts were identified. By analyzing highly expressed tissue-specific protein-coding transcripts (TSCTs) and noncoding transcripts (TSNTs), we found that testis expressed the highest numbers of TSCTs and TSNTs. Brain, monocytes, ovary, and heart expressed more TSCTs than the rest tissues, whereas brain, placenta, heart, and monocytes expressed more TSNTs than other tissues. Co-expression network constructed based on the TSCTs and TSNTs showed that each hub TSNT was co-expressed with several TSCTs, allowing functional annotation of TSNTs. Important biological processes and KEGG pathways highly related to the specific functions or diseases of each tissue were enriched with the corresponding TSCTs. These TSCTs and TSNTs may participate in the tissue-specific physiological or pathological processes. Our study provided a unique data set and systematic analysis of expression characteristics and functions of both TSCTs and TSNTs based on 14 distinct human tissues, and could facilitate future investigation of the mechanisms behind tissue-specific diseases and adverse drug reactions.

Comparing human and macaque placental transcriptomes to disentangle preterm birth pathology from gestational age effects

INTRODUCTION

A major issue in the transcriptomic study of spontaneous preterm birth (sPTB) in humans is the inability to collect healthy control tissue at the same gestational age (GA) to compare with pathologic preterm tissue. Thus, gene expression differences identified after the standard comparison of sPTB and term tissues necessarily reflect differences in both sPTB pathology and GA. One potential solution is to use GA-matched controls from a closely related species to tease apart genes that are dysregulated during sPTB from genes that are expressed differently as a result of GA effects.

METHODS

To disentangle genes whose expression levels are associated with sPTB pathology from those linked to GA, we compared RNA sequencing data from human preterm placentas, human term placentas, and rhesus macaque placentas at 80% completed gestation (serving as healthy non-human primate GA-matched controls). We first compared sPTB and term human placental transcriptomes to identify significantly differentially expressed genes. We then overlaid the results of the comparison between human sPTB and macaque placental transcriptomes to identify sPTB-specific candidates. Finally, we overlaid the results of the comparison between human term and macaque placental transcriptomes to identify GA-specific candidates.

RESULTS

Examination of relative expression for all human genes with macaque orthologs identified 267 candidate genes that were significantly differentially expressed between preterm and term human placentas. 29 genes were identified as sPTB-specific candidates and 37 as GA-specific candidates. Altogether, the 267 differentially expressed genes were significantly enriched for a variety of developmental, metabolic, reproductive, immune, and inflammatory functions. Although there were no notable differences between the functions of the 29 sPTB-specific and 37 GA-specific candidate genes, many of these candidates have been previously shown to be dysregulated in diverse pregnancy-associated pathologies.

DISCUSSION

By comparing human sPTB and term transcriptomes with GA-matched control transcriptomes from a closely related species, this study disentangled the confounding effects of sPTB pathology and GA, leading to the identification of 29 promising sPTB-specific candidate genes and 37 genes potentially related to GA effects. The apparent similarity in functions of the sPTB and GA candidates may suggest that the effects of sPTB and GA do not correspond to biologically distinct processes. Alternatively, it may reflect the poor state of knowledge of the transcriptional landscape underlying placental development and disease.

NeuN+ neuronal nuclei in non-human primate prefrontal cortex and subcortical white matter after clozapine exposure

Increased neuronal densities in subcortical white matter have been reported for some cases with schizophrenia. The underlying cellular and molecular mechanisms remain unresolved. We exposed 26 young adult macaque monkeys for 6 months to either clozapine, haloperidol or placebo and measured by structural MRI frontal gray and white matter volumes before and after treatment, followed by observer-independent, flow-cytometry-based quantification of neuronal and non-neuronal nuclei and molecular fingerprinting of cell-type specific transcripts. After clozapine exposure, the proportion of nuclei expressing the neuronal marker NeuN increased by approximately 50% in subcortical white matter, in conjunction with a more subtle and non-significant increase in overlying gray matter. Numbers and proportions of nuclei expressing the oligodendrocyte lineage marker, OLIG2, and cell-type specific RNA expression patterns, were maintained after antipsychotic drug exposure. Frontal lobe gray and white matter volumes remained indistinguishable between antipsychotic-drug-exposed and control groups. Chronic clozapine exposure increases the proportion of NeuN+ nuclei in frontal subcortical white matter, without alterations in frontal lobe volumes or cell type-specific gene expression. Further exploration of neurochemical plasticity in non-human primate brain exposed to antipsychotic drugs is warranted.

MTD: a mammalian transcriptomic database to explore gene expression and regulation

A systematic transcriptome survey is essential for the characterization and comprehension of the molecular basis underlying phenotypic variations. Recently developed RNA-seq methodology has facilitated efficient data acquisition and information mining of transcriptomes in multiple tissues/cell lines. Current mammalian transcriptomic databases are either tissue-specific or species-specific, and they lack in-depth comparative features across tissues and species. Here, we present a mammalian transcriptomic database (MTD) that is focused on mammalian transcriptomes, and the current version contains data from humans, mice, rats and pigs. Regarding the core features, the MTD browses genes based on their neighboring genomic coordinates or joint KEGG pathway and provides expression information on exons, transcripts and genes by integrating them into a genome browser. We developed a novel nomenclature for each transcript that considers its genomic position and transcriptional features. The MTD allows a flexible search of genes or isoforms with user-defined transcriptional characteristics and provides both table-based descriptions and associated visualizations. To elucidate the dynamics of gene expression regulation, the MTD also enables comparative transcriptomic analysis in both intraspecies and interspecies manner. The MTD thus constitutes a valuable resource for transcriptomic and evolutionary studies. The MTD is freely accessible at http://mtd.cbi.ac.cn.

Telomerase activation by genomic rearrangements in high-risk neuroblastoma

Neuroblastoma is a malignant paediatric tumour of the sympathetic nervous system. Roughly half of these tumours regress spontaneously or are cured by limited therapy. By contrast, high-risk neuroblastomas have an unfavourable clinical course despite intensive multimodal treatment, and their molecular basis has remained largely elusive. Here we have performed whole-genome sequencing of 56 neuroblastomas (high-risk, n = 39; low-risk, n = 17) and discovered recurrent genomic rearrangements affecting a chromosomal region at 5p15.33 proximal of the telomerase reverse transcriptase gene (TERT). These rearrangements occurred only in high-risk neuroblastomas (12/39, 31%) in a mutually exclusive fashion with MYCN amplifications and ATRX mutations, which are known genetic events in this tumour type. In an extended case series (n = 217), TERT rearrangements defined a subgroup of high-risk tumours with particularly poor outcome. Despite a large structural diversity of these rearrangements, they all induced massive transcriptional upregulation of TERT. In the remaining high-risk tumours, TERT expression was also elevated in MYCN-amplified tumours, whereas alternative lengthening of telomeres was present in neuroblastomas without TERT or MYCN alterations, suggesting that telomere lengthening represents a central mechanism defining this subtype. The 5p15.33 rearrangements juxtapose the TERT coding sequence to strong enhancer elements, resulting in massive chromatin remodelling and DNA methylation of the affected region. Supporting a functional role of TERT, neuroblastoma cell lines bearing rearrangements or amplified MYCN exhibited both upregulated TERT expression and enzymatic telomerase activity. In summary, our findings show that remodelling of the genomic context abrogates transcriptional silencing of TERT in high-risk neuroblastoma and places telomerase activation in the centre of transformation in a large fraction of these tumours.

A tree-like Bayesian structure learning algorithm for small-sample datasets from complex biological model systems

Background

There are increasing efforts to bring high-throughput systems biology techniques to bear on complex animal model systems, often with a goal of learning about underlying regulatory network structures (e.g., gene regulatory networks). However, complex animal model systems typically have significant limitations on cohort sizes, number of samples, and the ability to perform follow-up and validation experiments. These constraints are particularly problematic for many current network learning approaches, which require large numbers of samples and may predict many more regulatory relationships than actually exist.

Results

Here, we test the idea that by leveraging the accuracy and efficiency of classifiers, we can construct high-quality networks that capture important interactions between variables in datasets with few samples. We start from a previously-developed tree-like Bayesian classifier and generalize its network learning approach to allow for arbitrary depth and complexity of tree-like networks. Using four diverse sample networks, we demonstrate that this approach performs consistently better at low sample sizes than the Sparse Candidate Algorithm, a representative approach for comparison because it is known to generate Bayesian networks with high positive predictive value. We develop and demonstrate a resampling-based approach to enable the identification of a viable root for the learned tree-like network, important for cases where the root of a network is not known a priori. We also develop and demonstrate an integrated resampling-based approach to the reduction of variable space for the learning of the network. Finally, we demonstrate the utility of this approach via the analysis of a transcriptional dataset of a malaria challenge in a non-human primate model system, Macaca mulatta, suggesting the potential to capture indicators of the earliest stages of cellular differentiation during leukopoiesis.

Conclusions

We demonstrate that by starting from effective and efficient approaches for creating classifiers, we can identify interesting tree-like network structures with significant ability to capture the relationships in the training data. This approach represents a promising strategy for inferring networks with high positive predictive value under the constraint of small numbers of samples, meeting a need that will only continue to grow as more high-throughput studies are applied to complex model systems.

Electronic supplementary material

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