High-density rhesus macaque oligonucleotide microarray design using early-stage rhesus genome sequence information and human genome annotations

Genomic Tools for the Use of Nonhuman Primates in Translational Research

Nonhuman primates (NHPs) are important preclinical models for understanding the etiology of human diseases and for developing therapies and vaccines to cure or eliminate disease. Most human diseases have genetic components. Therefore, to be of maximal utility, the NHP species used for translational science should be as well characterized in regard to their genome and transcriptome as possible. This article reviews the current status of genomic information for the five NHP species used most often in translational research: rhesus macaque, cynomolgus macaque, vervet (African green) monkey, baboon, and marmoset NHP. These species have published whole genome sequences (with the exception of the baboon) and relatively well-characterized transcriptomes. Some have also been characterized in regard to specific genetic loci that are particularly related to translational concerns, such as the major histocompatability complex and the cytochrome P40 genes. Genomic resources to aid in stratifying captive populations in regard to genetic and phenotypic characteristics have been developed as an aid to enhancing reproducibility and facilitating more efficient use of animals. Taken together, the current genomic resources and numerous studies currently underway to improve them should enhance the value of NHPs as preclinical models of human disease.

Plant responses to environmental stresses-from gene to biotechnology

Increasing global population, urbanization and industrialization are increasing the rate of conversion of arable land into wasteland. Supplying food to an ever-increasing population is one of the biggest challenges that agriculturalists and plant scientists are currently confronting. Environmental stresses make this situation even graver. Despite the induction of several tolerance mechanisms, sensitive plants often fail to survive under environmental extremes. New technological approaches are imperative. Conventional breeding methods have a limited potential to improve plant genomes against environmental stress. Recently, genetic engineering has contributed enormously to the development of genetically modified varieties of different crops such as cotton, maize, rice, canola and soybean. The identification of stress-responsive genes and their subsequent introgression or overexpression within sensitive crop species are now being widely carried out by plant scientists. Engineering of important tolerance pathways, like antioxidant enzymes, osmolyte accumulation, membrane-localized transporters for efficient compartmentation of deleterious ions and accumulation of essential elements and resistance against pests or pathogens is also an area that has been intensively researched. In this review, the role of biotechnology and its successes, prospects and challenges in developing stress-tolerant crop cultivars are discussed.

Global transcriptional expression in ovarian follicles from Tsaiya ducks (Anas platyrhynchos) with a high-fertilization rate

Novel candidates for biomarkers of a high-fertilization rate were identified here through global transcriptional profiling of ovarian follicles. Some other differentially expressed candidate genes were first noted to influence animal reproduction in our previous cDNA microarray analysis and are now recognized as markers for marker-assisted selection. In the present study, we compared gene expression in ovarian follicles from animals with high- and low-fertilization rates using an oligonucleotide array. On the basis of a fold change of greater than 1.2 and less than -1.2, a difference of >100 Affymetrix arbitrary units between the two groups, and a P value of less than 0.05, 47 genes were found to be associated with fertilization rate. GOEAST and MetaCore software were further used to identify the functional categories of genes that were differentially expressed. Then, we focused on three interesting genes associated with a high-fertilization rate: one of these genes was discovered to participate in signaling pathways of fertilization, and two genes take roles in lipid metabolism. An oligonucleotide array showed that the levels of orthodenticle homeobox 2 (OTX2) and lecithin:cholesterol acyltransferase (LCAT) gene expression were 1.62-fold and 1.95-fold higher in the high-fertilization rate group than in the low-fertilization rate group, respectively (P < 0.05). The level of apolipoprotein A-I (APOA1) gene expression was also higher in the high-fertilization rate group, with a difference of 2.31-fold (P < 0.05). The data were validated through quantitative polymerase chain reaction analysis. These results confirm the usefulness of the array technique and data mining methods in the discovery of new biomarkers and add knowledge to our understanding of the factors affecting fertilization rates in ovarian follicles.

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.

Exosomes secreted by human cells transport largely mRNA fragments that are enriched in the 3'-untranslated regions

Small secreted membrane vesicles called exosomes have recently attracted a great interest after the discovery that they transfer mRNA that can be translated into protein in recipient cells. Surprisingly, we found that for the majority of exosomal mRNAs only a fraction of their corresponding probes is detectable on the expression microarrays. Exosomal mRNA fragmentation is characterized with a specific structural pattern. The closer to the 3'-end of the transcript the fragments are localized, the larger fraction among the secreted RNAs they constitute. Since the 3'-ends of transcripts contain elements conferring subcellular localization of mRNA and are rich in miRNA-binding sites, exosomal RNA may act as competing RNA to regulate stability, localization and translation activity of mRNAs in recipient cells.

REVIEWERS

This article was reviewed by Neil Smalheiser and Sandor Pongor.

'Omics investigations of HIV and SIV pathogenesis and innate immunity

In the 30 years since the advent of the AIDS epidemic, the biomedical community has put forward a battery of molecular therapies that are based on the accumulated knowledge of a limited number of viral targets. Despite these accomplishments, the community still confronts unanswered foundational questions about HIV infection. What are the cellular or biomolecular processes behind HIV pathogenesis? Can we elucidate the characteristics that distinguish those individuals who are naturally resistant to either infection or disease progression? The discovery of simian immunodeficiency viruses (SIVs) and the ensuing development of in vivo, nonhuman primate (NHP) infection models was a tremendous advance, especially in abetting the exploration of vaccine strategies. And while there have been numerous NHP infection models and vaccine trials performed, fundamental questions remain regarding host-virus interactions and immune correlates of protection. These issues are, perhaps, most starkly illustrated with the appreciation that many species of African nonhuman primates are naturally infected with strains of SIV that do not cause any appreciable disease while replicating to viral loads that match or exceed those seen with pathogenic SIV infections in Asian species of nonhuman primates. The last decade has seen the establishment of high-throughput molecular profiling tools, such as microarrays for transcriptomics, SNP arrays for genome features, and LC-MS techniques for proteins or metabolites. These provide the capacity to interrogate a biological model at a comprehensive, systems level, in contrast to historical approaches that characterized a few genes or proteins in an experiment. These methods have already had revolutionary impacts in understanding human diseases originating within the host genome such as genetic disorders and cancer, and the methods are finding increasing application in the context of infectious disease. We will provide a review of the use of such 'omics investigations as applied to understanding of HIV pathogenesis and innate immunity, drawing from our own research as well as the literature examples that utilized in vitro cell-based models or studies in nonhuman primates. We will also discuss the potential for systems biology to help guide strategies for HIV vaccines that offer significant protection by either preventing acquisition or strongly suppressing viral replication levels post-infection.

Characterization of single-nucleotide variation in Indian-origin rhesus macaques (Macaca mulatta)

BACKGROUND

Rhesus macaques are the most widely utilized nonhuman primate model in biomedical research. Previous efforts have validated fewer than 900 single nucleotide polymorphisms (SNPs) in this species, which limits opportunities for genetic studies related to health and disease. Extensive information about SNPs and other genetic variation in rhesus macaques would facilitate valuable genetic analyses, as well as provide markers for genome-wide linkage analysis and the genetic management of captive breeding colonies.

RESULTS

We used the available rhesus macaque draft genome sequence, new sequence data from unrelated individuals and existing published sequence data to create a genome-wide SNP resource for Indian-origin rhesus monkeys. The original reference animal and two additional Indian-origin individuals were resequenced to low coverage using SOLiD™ sequencing. We then used three strategies to validate SNPs: comparison of potential SNPs found in the same individual using two different sequencing chemistries, and comparison of potential SNPs in different individuals identified with either the same or different sequencing chemistries. Our approach validated approximately 3 million SNPs distributed across the genome. Preliminary analysis of SNP annotations suggests that a substantial number of these macaque SNPs may have functional effects. More than 700 non-synonymous SNPs were scored by Polyphen-2 as either possibly or probably damaging to protein function and these variants now constitute potential models for studying functional genetic variation relevant to human physiology and disease.

CONCLUSIONS

Resequencing of a small number of animals identified greater than 3 million SNPs. This provides a significant new information resource for rhesus macaques, an important research animal. The data also suggests that overall genetic variation is high in this species. We identified many potentially damaging non-synonymous coding SNPs, providing new opportunities to identify rhesus models for human disease.

Design, validation and annotation of transcriptome-wide oligonucleotide probes for the oligochaete annelid Eisenia fetida

High density oligonucleotide probe arrays have increasingly become an important tool in genomics studies. In organisms with incomplete genome sequence, one strategy for oligo probe design is to reduce the number of unique probes that target every non-redundant transcript through bioinformatic analysis and experimental testing. Here we adopted this strategy in making oligo probes for the earthworm Eisenia fetida, a species for which we have sequenced transcriptome-scale expressed sequence tags (ESTs). Our objectives were to identify unique transcripts as targets, to select an optimal and non-redundant oligo probe for each of these target ESTs, and to annotate the selected target sequences. We developed a streamlined and easy-to-follow approach to the design, validation and annotation of species-specific array probes. Four 244K-formatted oligo arrays were designed using eArray and were hybridized to a pooled E. fetida cRNA sample. We identified 63,541 probes with unsaturated signal intensities consistently above the background level. Target transcripts of these probes were annotated using several sequence alignment algorithms. Significant hits were obtained for 37,439 (59%) probed targets. We validated and made publicly available 63.5K oligo probes so the earthworm research community can use them to pursue ecological, toxicological, and other functional genomics questions. Our approach is efficient, cost-effective and robust because it (1) does not require a major genomics core facility; (2) allows new probes to be easily added and old probes modified or eliminated when new sequence information becomes available, (3) is not bioinformatics-intensive upfront but does provide opportunities for more in-depth annotation of biological functions for target genes; and (4) if desired, EST orthologs to the UniGene clusters of a reference genome can be identified and selected in order to improve the target gene specificity of designed probes. This approach is particularly applicable to organisms with a wealth of EST sequences but unfinished genome.

Genomic analysis reveals pre- and postchallenge differences in a rhesus macaque AIDS vaccine trial: insights into mechanisms of vaccine efficacy

We have employed global transcriptional profiling of whole blood to identify biologically relevant changes in cellular gene expression in response to alternative AIDS vaccine strategies with subsequent viral challenge in a rhesus macaque vaccine model. Samples were taken at day 0 (prechallenge), day 14 (peak viremia), and week 12 (set point) from animals immunized with replicating adenovirus type 5 host range (Ad5hr) recombinant viruses expressing human immunodeficiency virus HIV(env)(89.6P), simian immunodeficiency virus SIV(gag)(239), or SIV(nef)(239) alone or in combination with two intramuscular boosts with HIV(89.6P)gp140ΔCFI protein (L. J. Patterson et al., Virology 374:322-337, 2008), and each treatment resulted in significant control of viremia following simian-human immunodeficiency virus SHIV(89.6P) challenge (six animals per group plus six controls). At day 0, 8 weeks after the last treatment, the microarray profiles revealed significant prechallenge differences between treatment groups; data from the best-protected animals led to identification of a network of genes related to B cell development and lymphocyte survival. At peak viremia, expression profiles of the immunized groups were extremely similar, and comparisons to control animals reflected immunological differences other than effector T cell functions. Suggested protective mechanisms for vaccinated animals included upregulation of interleukin-27, a cytokine known to inhibit lentivirus replication, and increased expression of complement components, which may synergize with vaccine-induced antibodies. Divergent expression profiles at set point for the immunized groups implied distinct immunological responses despite phenotypic similarities in viral load and CD4(+) T cell levels. Data for the gp140-boosted group provided evidence for antibody-dependent, cell-mediated viral control, whereas animals immunized with only the replicating Ad5hr recombinants exhibited a different evolution of the B cell compartment even at 3 months postchallenge. This study demonstrates the sensitivity and discrimination of gene expression profiling of whole blood as an analytical tool in AIDS vaccine trials, providing unique insights into in vivo mechanisms and potential correlates of protection.

Collection of Macaca fascicularis cDNAs derived from bone marrow, kidney, liver, pancreas, spleen, and thymus

Background

Consolidating transcriptome data of non-human primates is essential to annotate primate genome sequences, and will facilitate research using non-human primates in the genomic era. Macaca fascicularis is a macaque monkey that is commonly used for biomedical and ecological research.

Findings

We constructed cDNA libraries of Macaca fascicularis, derived from tissues obtained from bone marrow, liver, pancreas, spleen, and thymus of a young male, and kidney of a young female. In total, 5'-end sequences of 56,856 clones were determined. Including the previously established cDNA libraries from brain and testis, we have isolated 112,587 cDNAs of Macaca fascicularis, which correspond to 56% of the curated human reference genes.

Conclusion

These sequences were deposited in the public sequence database as well as in-house macaque genome database . These data will become valuable resources for identifying functional parts of the genome of macaque monkeys in future studies.

Application of chicken microarrays for gene expression analysis in other avian species

BACKGROUND

With the threat of emerging infectious diseases such as avian influenza, whose natural hosts are thought to be a variety of wild water birds including duck, we are armed with very few genomic resources to investigate large scale immunological gene expression studies in avian species. Multiple options exist for conducting large gene expression studies in chickens and in this study we explore the feasibility of using one of these tools to investigate gene expression in other avian species.

RESULTS

In this study we utilised a whole genome long oligonucleotide chicken microarray to assess the utility of cross species hybridisation (CSH). We successfully hybridised a number of different avian species to this array, obtaining reliable signals. We were able to distinguish ducks that were infected with avian influenza from uninfected ducks using this microarray platform. In addition, we were able to detect known chicken immunological genes in all of the hybridised avian species.

CONCLUSION

Cross species hybridisation using long oligonucleotide microarrays is a powerful tool to study the immune response in avian species with little available genomic information. The present study validated the use of the whole genome long oligonucleotide chicken microarray to investigate gene expression in a range of avian species.

Transcriptional profiling in pathogenic and non-pathogenic SIV infections reveals significant distinctions in kinetics and tissue compartmentalization

Simian immunodeficiency virus (SIV) infection leads to AIDS in experimentally infected macaques, whereas natural reservoir hosts exhibit limited disease and pathology. It is, however, unclear how natural hosts can sustain high viral loads, comparable to those observed in the pathogenic model, without developing severe disease. We performed transcriptional profiling on lymph node, blood, and colon samples from African green monkeys (natural host model) and Asian pigtailed macaques (pathogenic model) to directly compare gene expression patterns during acute pathogenic versus non-pathogenic SIV infection. The majority of gene expression changes that were unique to either model were detected in the lymph nodes at the time of peak viral load. Results suggest a shift toward cellular stress pathways and Th1 profiles during pathogenic infection, with strong and sustained type I and II interferon responses. In contrast, a strong type I interferon response was initially induced during non-pathogenic infection but resolved after peak viral load. The natural host also exhibited controlled Th1 profiles and better preservation of overall cell homeostasis. This study identified gene expression patterns that are specific to disease susceptibility, tissue compartmentalization, and infection duration. These patterns provide a unique view of how host responses differ depending upon lentiviral infection outcome.

Simian immunodeficiency virus (SIV) does not cause disease in African green monkeys (a natural host for the virus), whereas experimentally infected Asian macaques (a non-natural host) develop a progressive disease that is similar to that which occurs in HIV-infected humans. Insight into how HIV causes disease and leads to development of AIDS may therefore be gained by comparing the response of natural and non-natural hosts to SIV infection. To this end, we examined changes that occurred in gene expression levels over time and in multiple tissues derived from African green monkeys and Asian macaques experimentally infected with SIV. Infection leads to host-specific gene expression patterns in lymph nodes, blood, and colon. The natural and non-natural hosts differed with respect to the timing, intensity, and duration of infection-induced gene expression changes associated with inflammation and response to stress.

Large-scale analysis of Macaca fascicularis transcripts and inference of genetic divergence between M. fascicularis and M. mulatta

Background

Cynomolgus macaques (Macaca fascicularis) are widely used as experimental animals in biomedical research and are closely related to other laboratory macaques, such as rhesus macaques (M. mulatta). We isolated 85,721 clones and determined 9407 full-insert sequences from cynomolgus monkey brain, testis, and liver. These sequences were annotated based on homology to human genes and stored in a database, QFbase .

Results

We found that 1024 transcripts did not represent any public human cDNA sequence and examined their expression using M. fascicularis oligonucleotide microarrays. Significant expression was detected for 544 (51%) of the unidentified transcripts. Moreover, we identified 226 genes containing exon alterations in the untranslated regions of the macaque transcripts, despite the highly conserved structure of the coding regions. Considering the polymorphism in the common ancestor of cynomolgus and rhesus macaques and the rate of PCR errors, the divergence time between the two species was estimated to be around 0.9 million years ago.

Conclusion

Transcript data from Old World monkeys provide a means not only to determine the evolutionary difference between human and non-human primates but also to unveil hidden transcripts in the human genome. Increasing the genomic resources and information of macaque monkeys will greatly contribute to the development of evolutionary biology and biomedical sciences.

Use of functional genomics to understand influenza-host interactions

Infection with influenza typically results in mild‐to‐moderate illness in healthy individuals; however, it is responsible for 30,000–40,000 deaths each year in the United States. In extreme cases, such as the influenza pandemic of 1918, tens of millions of people have died from the infection. To prepare for future influenza outbreaks, it is necessary to understand how the virus interacts with the host and to determine what makes certain strains of influenza highly pathogenic. Functional genomics provides a unique approach to this effort by allowing researchers to examine the effect of influenza infection on global host mRNA levels. Researchers are making increasing use of this approach to study virus–host interactions using a variety of model systems. For example, data obtained using microarray technology, in combination with mouse and macaque infection models, is providing exciting new insights into the pathogenicity of the 1918 virus. These studies suggest that the lethality associated with this virus is in part due to an aberrant and unchecked immune response. Progress is also being made toward using functional genomics in the diagnosis and prognosis of acute lung infections and in the development of more effective influenza vaccines and antivirals.

Lost in the world of functional genomics, systems biology, and translational research: is there life after the Milstein award?

We have always wanted to save the world from the scourges of virus infection by developing better drugs and vaccines. But fully understanding the intricacies of virus-host interactions, the first step in achieving this goal, requires the ability to view the process on a grand scale. The advent of high-throughput technologies, such as DNA microarrays and mass spectrometry, provided the first opportunities to obtain such a view. Here, we describe our efforts to use these tools to focus on the changes in cellular gene expression and protein abundance that occur in response to virus infection. By examining these changes in a comprehensive manner, we have been able to discover exciting new insights into innate immunity, interferon and cytokine signaling, and the strategies used by viruses to overcome these cellular defenses. Functional genomics may yet save the world from killer viruses.

Development of the first marmoset-specific DNA microarray (EUMAMA): a new genetic tool for large-scale expression profiling in a non-human primate

Background

The common marmoset monkey (Callithrix jacchus), a small non-endangered New World primate native to eastern Brazil, is becoming increasingly used as a non-human primate model in biomedical research, drug development and safety assessment. In contrast to the growing interest for the marmoset as an animal model, the molecular tools for genetic analysis are extremely limited.

Results

Here we report the development of the first marmoset-specific oligonucleotide microarray (EUMAMA) containing probe sets targeting 1541 different marmoset transcripts expressed in hippocampus. These 1541 transcripts represent a wide variety of different functional gene classes. Hybridisation of the marmoset microarray with labelled RNA from hippocampus, cortex and a panel of 7 different peripheral tissues resulted in high detection rates of 85% in the neuronal tissues and on average 70% in the non-neuronal tissues. The expression profiles of the 2 neuronal tissues, hippocampus and cortex, were highly similar, as indicated by a correlation coefficient of 0.96. Several transcripts with a tissue-specific pattern of expression were identified. Besides the marmoset microarray we have generated 3215 ESTs derived from marmoset hippocampus, which have been annotated and submitted to GenBank [GenBank: EF214838 – EF215447, EH380242 – EH382846].

Conclusion

We have generated the first marmoset-specific DNA microarray and demonstrated its use to characterise large-scale gene expression profiles of hippocampus but also of other neuronal and non-neuronal tissues. In addition, we have generated a large collection of ESTs of marmoset origin, which are now available in the public domain. These new tools will facilitate molecular genetic research into this non-human primate animal model.

Evolutionary and biomedical insights from the rhesus macaque genome

The rhesus macaque (Macaca mulatta) is an abundant primate species that diverged from the ancestors of Homo sapiens about 25 million years ago. Because they are genetically and physiologically similar to humans, rhesus monkeys are the most widely used nonhuman primate in basic and applied biomedical research. We determined the genome sequence of an Indian-origin Macaca mulatta female and compared the data with chimpanzees and humans to reveal the structure of ancestral primate genomes and to identify evidence for positive selection and lineage-specific expansions and contractions of gene families. A comparison of sequences from individual animals was used to investigate their underlying genetic diversity. The complete description of the macaque genome blueprint enhances the utility of this animal model for biomedical research and improves our understanding of the basic biology of the species.