The non-human primate reference transcriptome resource (NHPRTR) for comparative functional genomics

Discovery of NRG1-VII: the myeloid-derived class of NRG1

The growth factor Neuregulin-1 (NRG1) has pleiotropic roles in proliferation and differentiation of the stem cell niche in different tissues. It has been implicated in gut, brain and muscle development and repair. Six isoform classes of NRG1 and over 28 protein isoforms have been previously described. Here we report a new class of NRG1, designated NRG1-VII to denote that these NRG1 isoforms arise from a myeloid-specific transcriptional start site (TSS) previously uncharacterized. Long-read sequencing was used to identify eight high-confidence NRG1-VII transcripts. These transcripts presented major structural differences from one another, through the use of cassette exons and alternative stop codons. Expression of NRG1-VII was confirmed in primary human monocytes and tissue resident macrophages and induced pluripotent stem cell-derived macrophages (iPSC-derived macrophages). Isoform switching via cassette exon usage and alternate polyadenylation was apparent during monocyte maturation and macrophage differentiation. NRG1-VII is the major class expressed by the myeloid lineage, including tissue-resident macrophages. Analysis of public gene expression data indicates that monocytes and macrophages are a primary source of NRG1. The size and structure of class VII isoforms suggests that they may be more diffusible through tissues than other NRG1 classes. However, the specific roles of class VII variants in tissue homeostasis and repair have not yet been determined.

The identification of retro-DNAs in primate genomes as DNA transposons mobilizing via retrotransposition

Background

Mobile elements (MEs) constitute a major portion of the genome in primates and other higher eukaryotes, and they play important role in genome evolution and gene function. MEs can be divided into two fundamentally different classes: DNA transposons which transpose in the genome in a "cut-and-paste" style, and retrotransposons which propagate in a "copy-and-paste" fashion via a process involving transcription and reverse-transcription. In primate genomes, DNA transposons are mostly dead, while many retrotransposons are still highly active. We report here the identification of a unique group of MEs, which we call "retro-DNAs", for their combined characteristics of these two fundamentally different ME classes.

Methods

A comparative computational genomic approach was used to analyze the reference genome sequences of 10 primate species consisting of five apes, four monkeys, and marmoset.

Results

From our analysis, we identified a total of 1,750 retro-DNAs, representing 748 unique insertion events in the genomes of ten primate species including human. These retro-DNAs contain sequences of DNA transposons but lack the terminal inverted repeats (TIRs), the hallmark of DNA transposons. Instead, they show characteristics of retrotransposons, such as polyA tails, longer target-site duplications (TSDs), and the "TT/AAAA" insertion site motif, suggesting the use of the L1-based target- primed reverse transcription (TPRT) mechanism. At least 40% of these retro-DNAs locate into genic regions, presenting potentials for impacting gene function. More interestingly, some retro-DNAs, as well as their parent sites, show certain levels of expression, suggesting that they have the potential to create more retro-DNA copies in the present primate genomes.

Conclusions

Although small in number, the identification of these retro-DNAs reveals a new mean for propagating DNA transposons in primate genomes without active canonical DNA transposon activity. Our data also suggest that the TPRT machinery may transpose a wider variety of DNA sequences in the genomes.

Binding and sequestration of poison frog alkaloids by a plasma globulin

Alkaloids are important bioactive molecules throughout the natural world, and in many animals they serve as a source of chemical defense against predation. Dendrobatid poison frogs bioaccumulate alkaloids from their diet to make themselves toxic or unpalatable to predators. Despite the proposed roles of plasma proteins as mediators of alkaloid trafficking and bioavailability, the responsible proteins have not been identified. We use chemical approaches to show that a ~50 kDa plasma protein is the principal alkaloid-binding molecule in blood of poison frogs. Proteomic and biochemical studies establish this plasma protein to be a liver-derived alkaloid-binding globulin (ABG) that is a member of the serine-protease inhibitor (serpin) family. In addition to alkaloid-binding activity, ABG sequesters and regulates the bioavailability of 'free' plasma alkaloids in vitro. Unexpectedly, ABG is not related to saxiphilin, albumin, or other known vitamin carriers, but instead exhibits sequence and structural homology to mammalian hormone carriers and amphibian biliverdin-binding proteins. ABG represents a new small molecule binding functionality in serpin proteins, a novel mechanism of plasma alkaloid transport in poison frogs, and more broadly points toward serpins acting as tunable scaffolds for small molecule binding and transport across different organisms.

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.

Loop-Mediated Isothermal Amplification Detection of SARS-CoV-2 and Myriad Other Applications

As the second year of the COVID-19 pandemic begins, it remains clear that a massive increase in the ability to test for SARS-CoV-2 infections in a myriad of settings is critical to controlling the pandemic and to preparing for future outbreaks. The current gold standard for molecular diagnostics is the polymerase chain reaction (PCR), but the extraordinary and unmet demand for testing in a variety of environments means that both complementary and supplementary testing solutions are still needed. This review highlights the role that loop-mediated isothermal amplification (LAMP) has had in filling this global testing need, providing a faster and easier means of testing, and what it can do for future applications, pathogens, and the preparation for future outbreaks. This review describes the current state of the art for research of LAMP-based SARS-CoV-2 testing, as well as its implications for other pathogens and testing. The authors represent the global LAMP (gLAMP) Consortium, an international research collective, which has regularly met to share their experiences on LAMP deployment and best practices; sections are devoted to all aspects of LAMP testing, including preanalytic sample processing, target amplification, and amplicon detection, then the hardware and software required for deployment are discussed, and finally, a summary of the current regulatory landscape is provided. Included as well are a series of first-person accounts of LAMP method development and deployment. The final discussion section provides the reader with a distillation of the most validated testing methods and their paths to implementation. This review also aims to provide practical information and insight for a range of audiences: for a research audience, to help accelerate research through sharing of best practices; for an implementation audience, to help get testing up and running quickly; and for a public health, clinical, and policy audience, to help convey the breadth of the effect that LAMP methods have to offer.

The tweety Gene Family: From Embryo to Disease

The tweety genes encode gated chloride channels that are found in animals, plants, and even simple eukaryotes, signifying their deep evolutionary origin. In vertebrates, the tweety gene family is highly conserved and consists of three members—ttyh1, ttyh2, and ttyh3—that are important for the regulation of cell volume. While research has elucidated potential physiological functions of ttyh1 in neural stem cell maintenance, proliferation, and filopodia formation during neural development, the roles of ttyh2 and ttyh3 are less characterized, though their expression patterns during embryonic and fetal development suggest potential roles in the development of a wide range of tissues including a role in the immune system in response to pathogen-associated molecules. Additionally, members of the tweety gene family have been implicated in various pathologies including cancers, particularly pediatric brain tumors, and neurodegenerative diseases such as Alzheimer’s and Parkinson’s disease. Here, we review the current state of research using information from published articles and open-source databases on the tweety gene family with regard to its structure, evolution, expression during development and adulthood, biochemical and cellular functions, and role in human disease. We also identify promising areas for further research to advance our understanding of this important, yet still understudied, family of genes.

Functional and Adaptive Significance of Promoter Mutations That Affect Divergent Myocardial Expressions of TRIM72 in Primates

Cis-regulatory elements play important roles in tissue-specific gene expression and in the evolution of various phenotypes, and mutations in promoters and enhancers may be responsible for adaptations of species to environments. TRIM72 is a highly conserved protein that is involved in energy metabolism. Its expression in the heart varies considerably in primates, with high levels of expression in Old World monkeys and near absence in hominids. Here, we combine phylogenetic hypothesis testing and experimentation to demonstrate that mutations in promoter are responsible for the differences among primate species in the heart-specific expression of TRIM72. Maximum likelihood estimates of lineage-specific substitution rates under local-clock models show that relative to the evolutionary rate of introns, the rate of promoter was accelerated by 78% in the common ancestor of Old World monkeys, suggesting a role for positive selection in the evolution of the TRIM72 promoter, possibly driven by selective pressure due to changes in cardiac physiology after species divergence. We demonstrate that mutations in the TRIM72 promoter account for the differential myocardial TRIM72 expression of the human and the rhesus macaque. Furthermore, changes in TRIM72 expression alter the expression of genes involved in oxidative phosphorylation, which in turn affects mitochondrial respiration and cardiac energy capacity. On a broader timescale, phylogenetic regression analyses of data from 29 mammalian species show that mammals with high cardiac expression of TRIM72 have high heart rate, suggesting that the expression changes of TRIM72 may be related to differences in the heart physiology of those species.

Human-chimpanzee fused cells reveal cis-regulatory divergence underlying skeletal evolution

Gene regulatory divergence is thought to play a central role in determining human-specific traits. However, our ability to link divergent regulation to divergent phenotypes is limited. Here, we utilized human-chimpanzee hybrid induced pluripotent stem cells to study gene expression separating these species. The tetraploid hybrid cells allowed us to separate cis- from trans-regulatory effects, and to control for non-genetic confounding factors. We differentiated these cells into cranial neural crest cells (CNCCs), the primary cell type giving rise to the face. We discovered evidence of lineage-specific selection on the hedgehog signaling pathway, including a human-specific 6-fold down-regulation of EVC2 (LIMBIN), a key hedgehog gene. Inducing a similar down-regulation of EVC2 substantially reduced hedgehog signaling output. Mice and humans lacking functional EVC2 show striking phenotypic parallels to human-chimpanzee craniofacial differences, suggesting that the regulatory divergence of hedgehog signaling may have contributed to the unique craniofacial morphology of humans.

Multiple Alu Exonization in 3'UTR of a Primate-Specific Isoform of CYP20A1 Creates a Potential miRNA Sponge

Alu repeats contribute to phylogenetic novelties in conserved regulatory networks in primates. Our study highlights how exonized Alus could nucleate large-scale mRNA-miRNA interactions. Using a functional genomics approach, we characterize a transcript isoform of an orphan gene, CYP20A1 (CYP20A1_Alu-LT) that has exonization of 23 Alus in its 3'UTR. CYP20A1_Alu-LT, confirmed by 3'RACE, is an outlier in length (9 kb 3'UTR) and widely expressed. Using publically available data sets, we demonstrate its expression in higher primates and presence in single nucleus RNA-seq of 15,928 human cortical neurons. miRanda predicts ∼4,700 miRNA recognition elements (MREs) for ∼1,000 miRNAs, primarily originated within these 3'UTR-Alus. CYP20A1_Alu-LT could be a potential multi-miRNA sponge as it harbors ≥10 MREs for 140 miRNAs and has cytosolic localization. We further tested whether expression of CYP20A1_Alu-LT correlates with mRNAs harboring similar MRE targets. RNA-seq with conjoint miRNA-seq analysis was done in primary human neurons where we observed CYP20A1_Alu-LT to be downregulated during heat shock response and upregulated in HIV1-Tat treatment. In total, 380 genes were positively correlated with its expression (significantly downregulated in heat shock and upregulated in Tat) and they harbored MREs for nine expressed miRNAs which were also enriched in CYP20A1_Alu-LT. MREs were significantly enriched in these 380 genes compared with random sets of differentially expressed genes (P = 8.134e-12). Gene ontology suggested involvement of these genes in neuronal development and hemostasis pathways thus proposing a novel component of Alu-miRNA-mediated transcriptional modulation that could govern specific physiological outcomes in higher primates.

Alu master copies serve as the drivers of differential SINE transposition in recent primate genomes

Alu elements, averaging ~300bp in length, are a family of primate-specific short intersperse nuclear elements (SINEs) with more than one million copies and contributing to ~11% of primate genomes. Despite mostly being shared among primates, our recent study revealed highly differential recent Alu transposition among the genomes of primates from Hominidae and Cercopithecidae families. To understand the underlying mechanism, we analyzed six primate genomes and revealed species- and lineage-specific Alu profile exclusively defined by AluY composition. Among all Alus from the 6 genomes, we identified 5401 Alu master copies with 99% being from the AluY subfamily. The numbers of Alu master copies are positively correlated to the number of AluY elements in the genomes with the baboon genome having the largest number of most recent Alu master copies at high activities, while the crab-eating macaque genome having a low number of Alu master copies with low activity. Furthermore, the expression level of Alu master copies is positively correlated with their transposition activity. Our results support the concept that Alu transposition in primate genomes is driven by a small number of master copies, the number and relative activity of which contribute to the differential Alu transposition in recent primate genomes.

A mouse tissue atlas of small noncoding RNA

We report a systematic unbiased analysis of small RNA molecule expression in 11 different tissues of the model organism mouse. We discovered uncharacterized noncoding RNA molecules and identified that ∼30% of total noncoding small RNA transcriptome are distributed across the body in a tissue-specific manner with some also being sexually dimorphic. Distinct distribution patterns of small RNA across the body suggest the existence of tissue-specific mechanisms involved in noncoding RNA processing.

Small noncoding RNAs (ncRNAs) play a vital role in a broad range of biological processes both in health and disease. A comprehensive quantitative reference of small ncRNA expression would significantly advance our understanding of ncRNA roles in shaping tissue functions. Here, we systematically profiled the levels of five ncRNA classes (microRNA [miRNA], small nucleolar RNA [snoRNA], small nuclear RNA [snRNA], small Cajal body-specific RNA [scaRNA], and transfer RNA [tRNA] fragments) across 11 mouse tissues by deep sequencing. Using 14 biological replicates spanning both sexes, we identified that ∼30% of small ncRNAs are distributed across the body in a tissue-specific manner with some also being sexually dimorphic. We found that some miRNAs are subject to “arm switching” between healthy tissues and that tRNA fragments are retained within tissues in both a gene- and a tissue-specific manner. Out of 11 profiled tissues, we confirmed that brain contains the largest number of unique small ncRNA transcripts, some of which were previously annotated while others are identified in this study. Furthermore, by combining these findings with single-cell chromatin accessibility (scATAC-seq) data, we were able to connect identified brain-specific ncRNAs with their cell types of origin. These results yield the most comprehensive characterization of specific and ubiquitous small RNAs in individual murine tissues to date, and we expect that these data will be a resource for the further identification of ncRNAs involved in tissue function in health and dysfunction in disease.

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.

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.

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.

Brain Banks Spur New Frontiers in Neuropsychiatric Research and Strategies for Analysis and Validation

Neuropsychiatric disorders affect hundreds of millions of patients and families worldwide. To decode the molecular framework of these diseases, many studies use human postmortem brain samples. These studies reveal brain-specific genetic and epigenetic patterns via high-throughput sequencing technologies. Identifying best practices for the collection of postmortem brain samples, analyzing such large amounts of sequencing data, and interpreting these results are critical to advance neuropsychiatry. We provide an overview of human brain banks worldwide, including progress in China, highlighting some well-known projects using human postmortem brain samples to understand molecular regulation in both normal brains and those with neuropsychiatric disorders. Finally, we discuss future research strategies, as well as state-of-the-art statistical and experimental methods that are drawn upon brain bank resources to improve our understanding of the agents of neuropsychiatric disorders.

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.

A quantitative framework for characterizing the evolutionary history of mammalian gene expression

The evolutionary history of a gene helps predict its function and relationship to phenotypic traits. Although sequence conservation is commonly used to decipher gene function and assess medical relevance, methods for functional inference from comparative expression data are lacking. Here, we use RNA-seq across seven tissues from 17 mammalian species to show that expression evolution across mammals is accurately modeled by the Ornstein–Uhlenbeck process, a commonly proposed model of continuous trait evolution. We apply this model to identify expression pathways under neutral, stabilizing, and directional selection. We further demonstrate novel applications of this model to quantify the extent of stabilizing selection on a gene's expression, parameterize the distribution of each gene's optimal expression level, and detect deleterious expression levels in expression data from individual patients. Our work provides a statistical framework for interpreting expression data across species and in disease.

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.

Copy number variants implicate cardiac function and development pathways in earthquake-induced stress cardiomyopathy

The pathophysiology of stress cardiomyopathy (SCM), also known as takotsubo syndrome, is poorly understood. SCM usually occurs sporadically, often in association with a stressful event, but clusters of cases are reported after major natural disasters. There is some evidence that this is a familial condition. We have examined three possible models for an underlying genetic predisposition to SCM. Our primary study cohort consists of 28 women who suffered SCM as a result of two devastating earthquakes that struck the city of Christchurch, New Zealand, in 2010 and 2011. To seek possible underlying genetic factors we carried out exome analysis, genotyping array analysis, and array comparative genomic hybridization on these subjects. The most striking finding was the observation of a markedly elevated rate of rare, heterogeneous copy number variants (CNV) of uncertain clinical significance (in 12/28 subjects). Several of these CNVs impacted on genes of cardiac relevance including RBFOX1, GPC5, KCNRG, CHODL, and GPBP1L1. There is no physical overlap between the CNVs, and the genes they impact do not appear to be functionally related. The recognition that SCM predisposition may be associated with a high rate of rare CNVs offers a novel perspective on this enigmatic condition.

The aquatic animals' transcriptome resource for comparative functional analysis

Background

Aquatic animals have great economic and ecological importance. Among them, non-model organisms have been studied regarding eco-toxicity, stress biology, and environmental adaptation. Due to recent advances in next-generation sequencing techniques, large amounts of RNA-seq data for aquatic animals are publicly available. However, currently there is no comprehensive resource exist for the analysis, unification, and integration of these datasets. This study utilizes computational approaches to build a new resource of transcriptomic maps for aquatic animals. This aquatic animal transcriptome map database dbATM provides de novo assembly of transcriptome, gene annotation and comparative analysis of more than twenty aquatic organisms without draft genome.

Results

To improve the assembly quality, three computational tools (Trinity, Oases and SOAPdenovo-Trans) were employed to enhance individual transcriptome assembly, and CAP3 and CD-HIT-EST software were then used to merge these three assembled transcriptomes. In addition, functional annotation analysis provides valuable clues to gene characteristics, including full-length transcript coding regions, conserved domains, gene ontology and KEGG pathways. Furthermore, all aquatic animal genes are essential for comparative genomics tasks such as constructing homologous gene groups and blast databases and phylogenetic analysis.

Conclusion

In conclusion, we establish a resource for non model organism aquatic animals, which is great economic and ecological importance and provide transcriptomic information including functional annotation and comparative transcriptome analysis. The database is now publically accessible through the URL http://dbATM.mbc.nctu.edu.tw/.

Electronic supplementary material

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

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.

Tracing the evolution of the heterotrimeric G protein α subunit in Metazoa

BACKGROUND

Heterotrimeric G proteins are fundamental signaling proteins composed of three subunits, Gα and a Gβγ dimer. The role of Gα as a molecular switch is critical for transmitting and amplifying intracellular signaling cascades initiated by an activated G protein Coupled Receptor (GPCR). Despite their biochemical and therapeutic importance, the study of G protein evolution has been limited to the scope of a few model organisms. Furthermore, of the five primary Gα subfamilies, the underlying gene structure of only two families has been thoroughly investigated outside of Mammalia evolution. Therefore our understanding of Gα emergence and evolution across phylogeny remains incomplete.

RESULTS

We have computationally identified the presence and absence of every Gα gene (GNA-) across all major branches of Deuterostomia and evaluated the conservation of the underlying exon-intron structures across these phylogenetic groups. We provide evidence of mutually exclusive exon inclusion through alternative splicing in specific lineages. Variations of splice site conservation and isoforms were found for several paralogs which coincide with conserved, putative motifs of DNA-/RNA-binding proteins. In addition to our curated gene annotations, within Primates, we identified 15 retrotranspositions, many of which have undergone pseudogenization. Most importantly, we find numerous deviations from previous findings regarding the presence and absence of individual GNA- genes, nuanced differences in phyla-specific gene copy numbers, novel paralog duplications and subsequent intron gain and loss events.

CONCLUSIONS

Our curated annotations allow us to draw more accurate inferences regarding the emergence of all Gα family members across Metazoa and to present a new, updated theory of Gα evolution. Leveraging this, our results are critical for gaining new insights into the co-evolution of the Gα subunit and its many protein binding partners, especially therapeutically relevant G protein - GPCR signaling pathways which radiated in Vertebrata evolution.

Generation and characterization of the blood transcriptome of Macaca thibetana and comparative analysis with M. mulatta

Transcriptome profiles provide a large transcript sequence data set for genomic study, particularly in organisms that have no accurate genome data published. The Tibetan macaque (Macaca thibetana) is commonly considered to be an endemic species to China and an important animal in biomedical research in the present day. In the present study, we report the de novo assembly and characterization of the blood transcriptome of the Tibetan macaque from three individuals, and we also sequenced the blood transcriptome of the rhesus macaque (Macaca mulatta) for comparison. Using RNA-seq technology, 138 million sequencing reads of the M. thibetana transcriptome were generated. The assembly included 327 871 transcripts with an N50 of 1571 bp. According to the sequence similarity search, 80 317 (24.5%) transcripts were annotated in the nr protein database. All transcripts from M. thibetana and M. mulatta were functionally classified and compared using GO and KEGG analyses. The two transcriptomes were different in the GO term of nutrient reservoir activity, and in the KEGG subcategories of signaling molecules and interaction, infectious diseases, cell growth and death, and immune system. The transcriptomes in this study would provide a valuable resource for future functional and comparative genomic studies, and even for biological studies of this non-human primate.

Diurnal transcriptome atlas of a primate across major neural and peripheral tissues

Diurnal gene expression patterns underlie time-of-the-day-specific functional specialization of tissues. However, available circadian gene expression atlases of a few organs are largely from nocturnal vertebrates. We report the diurnal transcriptome of 64 tissues, including 22 brain regions, sampled every 2 hours over 24 hours, from the primate Papio anubis (baboon). Genomic transcription was highly rhythmic, with up to 81.7% of protein-coding genes showing daily rhythms in expression. In addition to tissue-specific gene expression, the rhythmic transcriptome imparts another layer of functional specialization. Most ubiquitously expressed genes that participate in essential cellular functions exhibit rhythmic expression in a tissue-specific manner. The peak phases of rhythmic gene expression clustered around dawn and dusk, with a "quiescent period" during early night. Our findings also unveil a different temporal organization of central and peripheral tissues between diurnal and nocturnal animals.

Cytoplasmic Linker Protein CLIP170 Negatively Regulates TLR4 Signaling by Targeting the TLR Adaptor Protein TIRAP

Cytoplasmic linker protein 170 (CLIP170) is a CAP-Gly domain-containing protein that is associated with the plus end of growing microtubules and implicated in various cellular processes, including the regulation of microtubule dynamics, cell migration, and intracellular transport. Our studies revealed a previously unrecognized property and role of CLIP170. We identified CLIP170 as one of the interacting partners of Brucella effector protein TcpB that negatively regulates TLR2 and TLR4 signaling. In this study, we demonstrate that CLIP170 interacts with the TLR2 and TLR4 adaptor protein TIRAP. Furthermore, our studies revealed that CLIP170 induces ubiquitination and subsequent degradation of TIRAP to negatively regulate TLR4-mediated proinflammatory responses. Overexpression of CLIP170 in mouse macrophages suppressed the LPS-induced expression of IL-6 and TNF-α whereas silencing of endogenous CLIP170 potentiated the levels of proinflammatory cytokines. In vivo silencing of CLIP170 in C57BL/6 mice by CLIP170-specific small interfering RNA enhanced LPS-induced IL-6 and TNF-α expression. Furthermore, we found that LPS modulates the expression of CLIP170 in mouse macrophages. Overall, our experimental data suggest that CLIP170 serves as an intrinsic negative regulator of TLR4 signaling that targets TIRAP.

Brain gene expression signature on primate genomic sequence evolution

Considering the overwhelming changes that occurred during primate evolution in brain structure and function, one might expect corresponding changes at the molecular level. Surprisingly, a relatively constrained gene expression pattern is observed in brain compared with other tissues among primates, an observation that calls for reassessment of RNA expression influence on primate genome evolution. We built phylogenetic trees based on genomic sequences of functional genomic regions and tissue-specific RNA expression in eight tissue types for six primate species. Comparisons of the phylogenetic trees from brain tissues revealed that DNA- and RNA-based trees were significantly similar. The similarity was specific for promoter regions and cerebellum and frontal cortex expression, suggesting a major impact of gene regulation in the brain on genome shaping along the primate branch.

Dynamic landscape and regulation of RNA editing in mammals

Adenosine-to-inosine (A-to-I) RNA editing is a conserved post-transcriptional mechanism mediated by ADAR enzymes that diversifies the transcriptome by altering selected nucleotides in RNA molecules. Although many editing sites have recently been discovered, the extent to which most sites are edited and how the editing is regulated in different biological contexts are not fully understood. Here we report dynamic spatiotemporal patterns and new regulators of RNA editing, discovered through an extensive profiling of A-to-I RNA editing in 8,551 human samples (representing 53 body sites from 552 individuals) from the Genotype-Tissue Expression (GTEx) project and in hundreds of other primate and mouse samples. We show that editing levels in non-repetitive coding regions vary more between tissues than editing levels in repetitive regions. Globally, ADAR1 is the primary editor of repetitive sites and ADAR2 is the primary editor of non-repetitive coding sites, whereas the catalytically inactive ADAR3 predominantly acts as an inhibitor of editing. Cross-species analysis of RNA editing in several tissues revealed that species, rather than tissue type, is the primary determinant of editing levels, suggesting stronger cis-directed regulation of RNA editing for most sites, although the small set of conserved coding sites is under stronger trans-regulation. In addition, we curated an extensive set of ADAR1 and ADAR2 targets and showed that many editing sites display distinct tissue-specific regulation by the ADAR enzymes in vivo. Further analysis of the GTEx data revealed several potential regulators of editing, such as AIMP2, which reduces editing in muscles by enhancing the degradation of the ADAR proteins. Collectively, our work provides insights into the complex cis- and trans-regulation of A-to-I editing.

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.

Fetal brain lesions after subcutaneous inoculation of Zika virus in a pregnant nonhuman primate

We describe the development of fetal brain lesions after Zika virus (ZIKV) inoculation in a pregnant pigtail macaque. Periventricular lesions developed within 10 d and evolved asymmetrically in the occipital-parietal lobes. Fetal autopsy revealed ZIKV in the brain and significant cerebral white matter hypoplasia, periventricular white matter gliosis, and axonal and ependymal injury. Our observation of ZIKV-associated fetal brain lesions in a nonhuman primate provides a model for therapeutic evaluation.

The SLC2A14 gene: genomic locus, tissue expression, splice variants, and subcellular localization of the protein

The SLC2A14 gene encodes for GLUT14, an orphan member of the facilitated membrane glucose transporter family, which was originally described to be exclusively expressed in human testis. However, genetic variations in SLC2A14 are associated with chronic diseases such as Alzheimer's disease and Inflammatory Bowel Disease, which cannot be explained by a strictly testicular expression. Therefore we analyzed available information on the SLC2A14 gene to update knowledge of the locus and its encoded products. This report presents an expanded SLC2A14 gene locus and a more diverse tissue expression, concurring with the existing evidence for disease associations. The exon utilization is tissue specific, with major expression in testis. When the 2 major testicular protein isoforms were expressed in mammalian cells, they located to the plasmalemma membrane, providing early evidence that GLUT14 could function as a membrane transporter.

Resequencing of the common marmoset genome improves genome assemblies and gene-coding sequence analysis

The first draft of the common marmoset (Callithrix jacchus) genome was published by the Marmoset Genome Sequencing and Analysis Consortium. The draft was based on whole-genome shotgun sequencing, and the current assembly version is Callithrix_jacches-3.2.1, but there still exist 187,214 undetermined gap regions and supercontigs and relatively short contigs that are unmapped to chromosomes in the draft genome. We performed resequencing and assembly of the genome of common marmoset by deep sequencing with high-throughput sequencing technology. Several different sequence runs using Illumina sequencing platforms were executed, and 181 Gbp of high-quality bases including mate-pairs with long insert lengths of 3, 8, 20, and 40 Kbp were obtained, that is, approximately 60× coverage. The resequencing significantly improved the MGSAC draft genome sequence. The N50 of the contigs, which is a statistical measure used to evaluate assembly quality, doubled. As a result, 51% of the contigs (total length: 299 Mbp) that were unmapped to chromosomes in the MGSAC draft were merged with chromosomal contigs, and the improved genome sequence helped to detect 5,288 new genes that are homologous to human cDNAs and the gaps in 5,187 transcripts of the Ensembl gene annotations were completely filled.

High-throughput RNA sequencing reveals structural differences of orthologous brain-expressed genes between western lowland gorillas and humans

The human brain and human cognitive abilities are strikingly different from those of other great apes despite relatively modest genome sequence divergence. However, little is presently known about the interspecies divergence in gene structure and transcription that might contribute to these phenotypic differences. To date, most comparative studies of gene structure in the brain have examined humans, chimpanzees, and macaque monkeys. To add to this body of knowledge, we analyze here the brain transcriptome of the western lowland gorilla (Gorilla gorilla gorilla), an African great ape species that is phylogenetically closely related to humans, but with a brain that is approximately one-third the size. Manual transcriptome curation from a sample of the planum temporale region of the neocortex revealed 12 protein-coding genes and one noncoding-RNA gene with exons in the gorilla unmatched by public transcriptome data from the orthologous human loci. These interspecies gene structure differences accounted for a total of 134 amino acids in proteins found in the gorilla that were absent from protein products of the orthologous human genes. Proteins varying in structure between human and gorilla were involved in immunity and energy metabolism, suggesting their relevance to phenotypic differences. This gorilla neocortical transcriptome comprises an empirical, not homology- or prediction-driven, resource for orthologous gene comparisons between human and gorilla. These findings provide a unique repository of the sequences and structures of thousands of genes transcribed in the gorilla brain, pointing to candidate genes that may contribute to the traits distinguishing humans from other closely related great apes.

What to compare and how: Comparative transcriptomics for Evo-Devo

Evolutionary developmental biology has grown historically from the capacity to relate patterns of evolution in anatomy to patterns of evolution of expression of specific genes, whether between very distantly related species, or very closely related species or populations. Scaling up such studies by taking advantage of modern transcriptomics brings promising improvements, allowing us to estimate the overall impact and molecular mechanisms of convergence, constraint or innovation in anatomy and development. But it also presents major challenges, including the computational definitions of anatomical homology and of organ function, the criteria for the comparison of developmental stages, the annotation of transcriptomics data to proper anatomical and developmental terms, and the statistical methods to compare transcriptomic data between species to highlight significant conservation or changes. In this article, we review these challenges, and the ongoing efforts to address them, which are emerging from bioinformatics work on ontologies, evolutionary statistics, and data curation, with a focus on their implementation in the context of the development of our database Bgee (http://bgee.org). J. Exp. Zool. (Mol. Dev. Evol.) 324B: 372–382, 2015. © 2015 The Authors. J. Exp. Zool. (Mol. Dev. Evol.) published by Wiley Periodicals, Inc.

GASS: genome structural annotation for Eukaryotes based on species similarity

Background

With the development of high-throughput sequencing techniques, more and more genomes were sequenced and assembled. However, annotating a genome’s structure rapidly and expressly remains challenging. Current eukaryotic genome annotations require various, abundant supporting data, such as: species-specific and cross-species protein sequences, ESTs, cDNA and RNA-Seq data. Collecting those data and merging their analytical results to achieve a consistent complete annotation is a complex, time and cost consuming task.

Results

In our study, we proposed a fast and easy-to-use computational tool: GASS (Genome Annotation based on Species Similarity). It annotates a eukaryotic genome based on only the annotations from another similar species. With aligning the exons’ sequences of an annotated similar species to the un-annotated genome, GASS detects the optimal transcript annotations with a shortest-path model. In our study, GASS was used to achieve the rhesus annotations based on the human annotations. The produced annotations were evaluated by comparing them to the two existing rhesus annotation databases (RefSeq and Ensembl) directly and being aligned with three RNA-Seq data of rhesus. The experiment results showed that more than 65% RefSeq exons and splicing junctions were exactly found by GASS. GASS’s sensitivity was higher than RefSeq’s, and was close to Ensembl’s. GASS had higher specificities than Ensembl at gene, transcript, exon and splicing junction levels. We also found the mis-assemblies of rheMac3 genome, which led to the 2 bp shifts in annotating position on exons’ boundary and then the incomplete splicing canonical sites in Refseq annotations. These detections were further supported by various data sources.

Conclusions

GASS quickly produces structural genome annotations in sufficient abundance and accuracy. With simple and rapid running of GASS, small labs can create quick views of genome annotations for an un-annotated species, without the necessity to create, collect, analyze and synthesize extra various data sources, or wait several months for the annotations from professional organizations. GASS can be applied to many study occasions, such as the analysis of RNA-Seq datasets from the unannotated species whose genome drafts are available but the annotations are not.

Electronic supplementary material

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

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

The non-human primate reference transcriptome resource (NHPRTR, available online at http://nhprtr.org/) aims to generate comprehensive RNA-seq data from a wide variety of non-human primates (NHPs), from lemurs to hominids. In the 2012 Phase I of the NHPRTR project, 19 billion fragments or 3.8 terabases of transcriptome sequences were collected from pools of ∼20 tissues in 15 species and subspecies. Here we describe a major expansion of NHPRTR by adding 10.1 billion fragments of tissue-specific RNA-seq data. For this effort, we selected 11 of the original 15 NHP species and subspecies and constructed total RNA libraries for the same ∼15 tissues in each. The sequence quality is such that 88% of the reads align to human reference sequences, allowing us to compute the full list of expression abundance across all tissues for each species, using the reads mapped to human genes. This update also includes improved transcript annotations derived from RNA-seq data for rhesus and cynomolgus macaques, two of the most commonly used NHP models and additional RNA-seq data compiled from related projects. Together, these comprehensive reference transcriptomes from multiple primates serve as a valuable community resource for genome annotation, gene dynamics and comparative functional analysis.

Enabling large-scale next-generation sequence assembly with Blacklight

A variety of extremely challenging biological sequence analyses were conducted on the XSEDE large shared memory resource Blacklight, using current bioinformatics tools and encompassing a wide range of scientific applications. These include genomic sequence assembly, very large metagenomic sequence assembly, transcriptome assembly, and sequencing error correction. The data sets used in these analyses included uncategorized fungal species, reference microbial data, very large soil and human gut microbiome sequence data, and primate transcriptomes, composed of both short-read and long-read sequence data. A new parallel command execution program was developed on the Blacklight resource to handle some of these analyses. These results, initially reported previously at XSEDE13 and expanded here, represent significant advances for their respective scientific communities. The breadth and depth of the results achieved demonstrate the ease of use, versatility, and unique capabilities of the Blacklight XSEDE resource for scientific analysis of genomic and transcriptomic sequence data, and the power of these resources, together with XSEDE support, in meeting the most challenging scientific problems.

Detecting and correcting systematic variation in large-scale RNA sequencing data

High-throughput RNA sequencing (RNA-seq) enables comprehensive scans of entire transcriptomes, but best practices for analyzing RNA-seq data have not been fully defined, particularly for data collected with multiple sequencing platforms or at multiple sites. Here we used standardized RNA samples with built-in controls to examine sources of error in large-scale RNA-seq studies and their impact on the detection of differentially expressed genes (DEGs). Analysis of variations in guanine-cytosine content, gene coverage, sequencing error rate and insert size allowed identification of decreased reproducibility across sites. Moreover, commonly used methods for normalization (cqn, EDASeq, RUV2, sva, PEER) varied in their ability to remove these systematic biases, depending on sample complexity and initial data quality. Normalization methods that combine data from genes across sites are strongly recommended to identify and remove site-specific effects and can substantially improve RNA-seq studies.

Multi-platform assessment of transcriptome profiling using RNA-seq in the ABRF next-generation sequencing study

High-throughput RNA sequencing (RNA-seq) greatly expands the potential for genomics discoveries, but the wide variety of platforms, protocols and performance capabilitites has created the need for comprehensive reference data. Here we describe the Association of Biomolecular Resource Facilities next-generation sequencing (ABRF-NGS) study on RNA-seq. We carried out replicate experiments across 15 laboratory sites using reference RNA standards to test four protocols (poly-A-selected, ribo-depleted, size-selected and degraded) on five sequencing platforms (Illumina HiSeq, Life Technologies PGM and Proton, Pacific Biosciences RS and Roche 454). The results show high intraplatform (Spearman rank R > 0.86) and inter-platform (R > 0.83) concordance for expression measures across the deep-count platforms, but highly variable efficiency and cost for splice junction and variant detection between all platforms. For intact RNA, gene expression profiles from rRNA-depletion and poly-A enrichment are similar. In addition, rRNA depletion enables effective analysis of degraded RNA samples. This study provides a broad foundation for cross-platform standardization, evaluation and improvement of RNA-seq.

Infection with MERS-CoV causes lethal pneumonia in the common marmoset

The availability of a robust disease model is essential for the development of countermeasures for Middle East respiratory syndrome coronavirus (MERS-CoV). While a rhesus macaque model of MERS-CoV has been established, the lack of uniform, severe disease in this model complicates the analysis of countermeasure studies. Modeling of the interaction between the MERS-CoV spike glycoprotein and its receptor dipeptidyl peptidase 4 predicted comparable interaction energies in common marmosets and humans. The suitability of the marmoset as a MERS-CoV model was tested by inoculation via combined intratracheal, intranasal, oral and ocular routes. Most of the marmosets developed a progressive severe pneumonia leading to euthanasia of some animals. Extensive lesions were evident in the lungs of all animals necropsied at different time points post inoculation. Some animals were also viremic; high viral loads were detected in the lungs of all infected animals, and total RNAseq demonstrated the induction of immune and inflammatory pathways. This is the first description of a severe, partially lethal, disease model of MERS-CoV, and as such will have a major impact on the ability to assess the efficacy of vaccines and treatment strategies as well as allowing more detailed pathogenesis studies.

The development of vaccines and treatment strategies is aided by robust animal disease models that accurately depict the illness that is observed in humans. Here we describe a new, improved model for MERS-CoV using the common marmoset, whereby the severe, and even lethal, illness that has been observed in many human cases is recapitulated. Prior to the development of this model, the only available animal models for MERS-CoV infection were the rhesus macaque and a mouse model that requires adenovirus-transduced expression of the human version of the protein required for virus entry. The rhesus macaque model more closely mimics the mild to moderate disease observed in some patients—mainly those without significant comorbidities. The increased severity of illness in the common marmoset model is an important advance in the ability to evaluate potential therapeutic agents against MERS-CoV, as discrimination between successfully treated and control animals should be more apparent. In addition, the closer models recapitulate the disease observed in humans, the more likely findings can be eventually translated into use in humans.

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

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

In transition: primate genomics at a time of rapid change

The field of nonhuman primate genomics is undergoing rapid change and making impressive progress. Exploiting new technologies for DNA sequencing, researchers have generated new whole-genome sequence assemblies for multiple primate species over the past 6 years. In addition, investigations of within-species genetic variation, gene expression and RNA sequences, conservation of non-protein-coding regions of the genome, and other aspects of comparative genomics are moving at an accelerating speed. This progress is opening a wide array of new research opportunities in the analysis of comparative primate genome content and evolution. It also creates new possibilities for the use of nonhuman primates as model organisms in biomedical research. This transition, based on both new technology and the new information being generated in regard to human genetics, provides an important justification for reevaluating the research goals, strategies, and study designs used in primate genetics and genomics.

Old world monkeys and new age science: the evolution of nonhuman primate systems virology

Nonhuman primate (NHP) biomedical models are critical to our understanding of human health and disease, yet we are still in the early stages of developing sufficient tools to support primate genomic research that allow us to better understand the basis of phenotypic traits in NHP models of disease. A mere 7 years ago, the limited NHP transcriptome profiling that was being performed was done using complementary DNA arrays based on human genome sequences, and the lack of NHP genomic information and immunologic reagents precluded the use of NHPs in functional genomic studies. Since then, significant strides have been made in developing genomics capabilities for NHP research, from the rhesus macaque genome sequencing project to the construction of the first macaque-specific high-density oligonucleotide microarray, paving the way for further resource development and additional primate sequencing projects. Complete published draft genome sequences are now available for the chimpanzee ( Chimpanzee Sequencing Analysis Consortium 2005), bonobo ( Prufer et al. 2012), gorilla ( Scally et al. 2012), and baboon ( Ensembl.org 2013), along with the recently completed draft genomes for the cynomolgus macaque and Chinese rhesus macaque. Against this backdrop of both expanding sequence data and the early application of sequence-derived DNA microarrays tools, we will contextualize the development of these community resources and their application to infectious disease research through a literature review of NHP models of acquired immune deficiency syndrome and models of respiratory virus infection. In particular, we will review the use of -omics approaches in studies of simian immunodeficiency virus and respiratory virus pathogenesis and vaccine development, emphasizing the acute and innate responses and the relationship of these to the course of disease and to the evolution of adaptive immunity.

Deep transcriptional sequencing of mucosal challenge compartment from rhesus macaques acutely infected with simian immunodeficiency virus implicates loss of cell adhesion preceding immune activation

Pathology resulting from human immunodeficiency virus (HIV) infection is driven by protracted inflammation; the primary loss of CD4(+) T cells is caused by activation-driven apoptosis. Recent studies of nonhuman primates (NHPs) have suggested that during the acute phase of infection, antiviral mucosal immunity restricts viral replication in the primary infection compartment. These studies imply that HIV achieves systemic infection as a consequence of a failure in host antiviral immunity. Here, we used high-dose intrarectal inoculation of rhesus macaques with simian immunodeficiency virus (SIV) SIVmac251 to examine how the mucosal immune system is overcome by SIV during acute infection. The host response in rectal mucosa was characterized by deep mRNA sequencing (mRNA-seq) at 3 and 12 days postinoculation (dpi) in 4 animals for each time point. While we observed a strong host transcriptional response at 3 dpi, functions relating to antiviral immunity were absent. Instead, we observed a significant number of differentially expressed genes relating to cell adhesion and reorganization of the cytoskeleton. We also observed downregulation of genes encoding members of the claudin family of cell adhesion molecules, which are coexpressed with genes associated with pathology in the colorectal mucosa, and a large number of noncoding transcripts. In contrast, at 12 dpi the differentially expressed genes were enriched in those involved with immune system functions, in particular, functions relating to T cells, B cells, and NK cells. Our findings indicate that host responses that negatively affect mucosal integrity occur before inflammation. Consequently, when inflammation is activated at peak viremia, mucosal integrity is already compromised, potentially enabling rapid tissue damage, driving further inflammation. Importance: The HIV pandemic is one of the major threats to human health, causing over a million deaths per year. Recent studies have suggested that mucosal antiviral immune responses play an important role in preventing systemic infection after exposure to the virus. Yet, despite their potential role in decreasing transmission rates between individuals, these antiviral mechanisms are poorly understood. Here, we carried out the first deep mRNA sequencing analysis of mucosal host responses in the primary infection compartment during acute SIV infection. We found that during acute infection, a significant host response was mounted in the mucosa before inflammation was triggered. Our analysis indicated that the response has a detrimental effect on tissue integrity, causing increased permeability, tissue damage, and recruitment of SIV target cells. These results emphasize the importance of mucosal host responses preceding immune activation in preventing systemic SIV infection.