- Research article
- Open Access
Diversity and regulatory impact of copy number variation in the primate Macaca fascicularis
© The Author(s). 2017
Received: 30 April 2016
Accepted: 1 February 2017
Published: 10 February 2017
Copy number variations (CNVs) are a significant source of genetic diversity and commonly found in mammalian genomes. We have generated a genome-wide CNV map for Cynomolgus monkeys (Macaca fascicularis). This crab-eating macaque is the closest animal model to humans that is used in biomedical research.
We show that Cynomolgus monkey CNVs are in general much smaller in size than gene loci and are specific to the population of origin. Genome-wide expression data from five vitally important organs demonstrates that CNVs in close proximity to transcription start sites associate strongly with expression changes. Among these eQTL genes we find an overrepresentation of genes involved in metabolism, receptor activity, and transcription.
These results provide evidence that CNVs shape tissue transcriptomes in monkey populations, potentially offering an adaptive advantage. We suggest that this genetic diversity should be taken into account when using Cynomolgus macaques as models.
Copy number variations (CNVs) are genetic differences in the normal population displayed as microscopically invisible deletions or amplifications of stretches of genomic DNA ranging from 1 kilobase up to the megabase scale . CNVs are commonly found in the genomes of humans , primates , rodents , or arthropods like Drosophila melanogaster . In humans, more than 6.2 million different CNVs mapping to ~500’000 genomic regions have so far been identified . They significantly contribute to genetic variation, covering more nucleotide content per genome than single nucleotide polymorphisms (e.g. approximately 0.8% of the length of the human genome differs between two human individuals) . Furthermore CNVs exhibit a higher per-locus mutation rate than single nucleotide polymorphisms (SNPs) . Since CNVs can reside in genomic regions harboring genes they can alter gene dosage, disrupt coding sequences or modify the level and timing of gene expression for genes within the CNV [9, 10] and on its flanks [4, 11–15]. These effects of CNVs are difficult to understand and not necessarily predictable, but relevant for many diseases [16–19] and pharmacological responses like in the case of CYP2D6 CNVs .
Cynomolgus monkeys (Macaca fascicularis) are well-established translational models for biomedical research and drug testing. These non-human primates are one of the closest animal model to humans with high genetic similarity (~93% in nucleotide sequence identity), similar anatomies, and very similar physiologies [21–23]. These animals offer great promise as models for many aspects of human health and disease. Cynomolgus monkeys are outbred species, caught in the wild in many different places of peninsular Southeast Asia, the Philippines, and Mauritius, and used to found and continuously refresh breeding programs [23, 24]. They exhibit substantial levels of genetic variation which can affect the outcome and interpretation of biomedical studies [25–28]. Understanding of the contribution of this variation to phenotypes is lagging behind in Cynomolgus monkeys compared to the knowledge about human genetic and genomic variation . Genome-wide catalogs of SNPs start to emerge for Cynomolgus monkeys with more and more genome sequencing projects published [21, 29–32]. However, information on structural variants, such as CNVs, is not available for Cynomolgus monkeys despite their prominent role in phenotypic variation. In this study, we comprehensively assess for the first time genome-wide copy number variation among Cynomolgus monkeys from cohorts used in pharmaceutical studies using a custom 4.2 million probes comparative genomic hybridization (CGH) array. To investigate the potential functional implications of the detected copy number variation, we used a Cynomolgus monkey specific gene expression microarray to associate CNV genotypes with expression changes of proximal genes using a cis expression quantitative trait loci (cis-eQTL) mapping approach.
CNVs in diverse Cynomolgus monkey populations
To minimize genotyping errors we then calculated “quantitative genotypes”, defined as an individual’s median log2-ratio of all CGH probes within a given CNVR, and used these for all further analyses (Additional file 2: Table S3). To shed light on the relationship between individual monkeys and their CNV profiles, we used unsupervised hierarchical clustering and principal component analysis (PCA). Hierarchical clustering of the CNV signals showed a clear grouping of the samples by geographical origin (Fig. 1e). PCA revealed also a clear separation of the individuals by their geographical origin, particularly in monkeys originating from Mauritius (Additional file 1: Figure S1). Interestingly, this finding was mainly driven by the second principal component, whilst the first principal component correlated with the duplication/deletion ratio of individual monkeys (Fig. 1b). We assumed different population genetic properties of duplications and deletions, similar to previous observations in great apes and humans [37, 38]. Therefore we performed PCA on deletions and duplications separately to assess which structural variant serves as more reliable population genetic marker. In this analysis deletion genotypes segregate most accurately Cynomolgus monkeys according to the different island and mainland populations (Fig. 1f and Additional file 1: Figure S1).
Gene expression and CNVs
Closer analysis of the CNV-eQTL associations revealed increasing significance and frequency for CNVs in close proximity to the transcription start site (TSS) (Fig. 2c). The highest density of eQTL CNVs was found around 200 kb upstream of the TSS (Fig. 2d), which usually marks DNaseI hypersensitive regulatory DNA regions . However, eQTL CNVs were generally not found to be located more often upstream of the TSS than downstream (Binomial test, p = 0.56).
Previous research suggests that CNVs might affect gene expression changes of multiple proximal genes . We therefore sought to discover additional, weaker associations within the detected eQTL regions, which we might have missed in the genome-wide eQTL mapping. For each eQTL region, all associations between genes within 1 Mb of an eQTL gene’s TSS and the eQTL CNV were assessed and corrected for multiple testing per region by Bonferroni correction. A total of 12 additional associations were detected in ten out of 31 non-overlapping eQTL regions. Within these regions, on average 11.3% of genes were also associated with the eQTL CNV and all associations showed the same directionality as the eQTL.
Taking advantage of the recently sequenced Cynomolgus monkey genome, we developed species-specific microarrays for CNV detection and transcriptome profiling to explore CNV diversity, genome variation, and transcriptional consequences in natural monkey populations. For this approach we developed a computational pipeline combining careful data normalization and 3 different, non-redundant CNV calling methods enabling us to assess copy number variation in Cynomolgus monkeys at very high resolution (>2 kb). According to our knowledge, this is the first study to investigate natural copy number variation in this species.
We predominately find small CNVs among our individuals with a median length of about 8 kb. Technically, this finding highlights the importance of a meticulous CNV calling approach, since we operate close to the resolution limit of the array with many CNVs encompassing only 5 probes. Biologically, this finding is in line with current human genetic research, which suggests that individuals from normal and healthy populations carry mostly short CNVs (median length of 7.4 kb) . Genome-wide mapping of CNVs to copy number variable regions revealed that more than half of all CNVs were population specific with a majority detected in only one individual, suggesting that many of the population specific CNVs are probably found at low frequency at the respective loci. Genotyping of Cynomolgus monkeys based on CNVs revealed a clustering by origin, clearly separating the Mauritian and Philippine island populations from mainland populations, an effect probably mediated by genetic bottlenecks and geographical isolation. For example Mauritian monkeys are descendants from a small number of individuals brought to the island in the 16th century . Separation between the Southeast Asian and Vietnamese mainland populations is less pronounced, in line with the fact that these two populations share geographically adjacent biotops . In this context we also discovered that duplications are less informative than deletions as population genetic markers in agreement with CNV profiling studies in great apes and humans [37, 38]. These studies have shown that deletions are less likely than duplications to be subjected to recurrent mutation events. Hence deletions cannot change their copy number state so dynamically over a short time period and are more likely to exhibit identity by descent as a result of an ancestral mutation event.
Since the median gene locus size in Cynomolgus monkeys is 47 kb, the majority of CNVs are too small to delete or alter the copy number of an entire gene making them less likely to become an evolutionary constraint. Furthermore, short CNVs are more frequently generated de novo than large CNVs (>500 kb) , which indicates that they are not under strong purifying selection in contrast to deleterious large CNVs . Therefore the genomic regions where most CNVs are found might be those where copy number is less important or where CNVs affect regulatory regions (enhancers, insulators, silencers) that have modulatory impact on gene expression. We find that a relatively small number of CNV loci are associated with gene expression changes, suggesting that most copy number variation have no link to gene expression regulation, similar to single nucleotide variants. When looking at the association strength of our eQTLs, it became evident that the strongest associations are detected with variants close to their targets’ TSS. This relationship between genetic distance to TSS and stronger association effects appears to be similar between CNVs and SNPs. However eQTL-SNPs distribute more symmetrically in very close proximity to the TSS (within +/− 100 kb) [13, 43, 44], whilst our eQTL-CNV distributions peak ~200 kb upstream of the TSS. Therefore it is possible that eQTL-CNVs indicate more tissue-specific regulatory DNA regions or enhancer elements, whilst the majority of eQTL-SNPs mark general, tissue-independent regulatory elements which tend to be found in very close proximity to the TSS [45, 46]. Indeed, almost all of the CNV-eQTLs identified in our study exhibit tissue-specific effects. Furthermore, the eQTL-genes show generally lower expression levels than the respective tissue average. This observation might indicate, that highly expressed genes are either depleted of regulatory variants by purifying selection or that gene regulatory networks can buffer the effects of regulatory variants. However, lower expressed genes with insufficient network mediated buffering might be more responsive to CNVs in regulatory DNA regions. The resulting change in expression level might offer an adaptive advantage in a tissue-specific manner. In this respect, we highlight two particularly interesting eQTL regions on Cynomolgus monkey chromosome 7 and 3. A duplication event on chromosome 7 changes the expression levels of a group of OR genes in lung and kidney from a very low level to a level close to the average global gene expression. Although olfactory receptors are typically not expected to be expressed in internal organs, they belong to genes relevant to the immediate environment and have been found overrepresented amongst copy number variable genes . Recent studies highlighted a role of such receptors in kidney where their chemosensory function plays a role in tracking of the chemical composition of blood and tubular fluids and the modulation of renin secretion, blood pressure, and glomerular filtration rate [48, 49]. On chromosome 3 we detected a deletion event associated with increased expression of ABCB4 in lung. This gene product is known to act as tumor suppressor once overexpressed in lung cancer . Since ABCB4 expression was shown to be regulated by epigenetic silencing, a deletion 480 kb upstream of this gene might abolish epigenetic control of transcription, possibly resulting in constitutive overexpression. Carrier of such copy number variations modulating expression levels of a tumor suppressor gene might benefit from protective advantages.
We describe for the first time substantial copy number variation in natural Macaca fascicularis populations as an additional genetic source of diversity and interindividual variation. We report several tissue specific associations between CNVs and expression levels of proximal genes providing a molecular link for variable transcriptional programs between individuals. As an example a genomic region on chromosome 7, which harbors several OR genes, shows an association between gene expression and a close-by duplication event in both kidney and lung. Our data suggest that CNVs shape tissue transcriptomes of vitally important organs, possibly offering an adaptive advantage.
All peripheral blood and tissue samples (heart, kidney, liver, lung, spleen) used in this study were supplied by AAALAC accredited contract research organizations. Animal samples were taken from healthy, untreated animals of GLP drug-safety studies in accordance with local, national and international regulations. All procedures were approved by the Institutional Animal Care and Use Committee (IACUC) and governmental agencies responsible for animal welfare in compliance with local laws and regulations.
Tissue samples for CGH and expression analysis were taken from Cynomolus monkeys obtained from breeding centers located in the Philippines (3 females and 3 males), in Vietnam (2 males, 2 females), in China for animals from Mainland Southeast Asia (3 females), or in Mauritius (4 females and 4 males) (Fig. 1a). Blood samples for CGH analysis were taken from Cynomolgus monkeys with Mauritian origin (25 males). Details (gender, weight, age, origin) of all animals and their suppliers are on record and were part of the data submitted to public databases.
NimbleGen gene expression analysis
Cynomolgus monkey tissues were homogenized in tubes prefilled with 1.4 mm ceramic beads and QiaGen’s lysis reagent RLT using a FastPrep-24 instrument (MP Biomedicals, Solon, OH, USA). Total RNA from lysates was extracted using the RNeasy Mini kit combined with DNase treatment on a solid support (Qiagen Inc., Valencia, CA, USA). RNA quality assessment and quantification was performed using microfluidic chip analysis on an Agilent 2100 bioanalyzer (Agilent Technologies Inc., Santa Clara, CA, USA). On a Biomek FXp workstation (Beckman Coulter Inc., Brea, CA, USA), 10 ng of total RNA was used to prepare cDNA with the NuGen Ovation Pico WTA System V2 (NuGEN Technologies, Inc., SanCarlos, CA, USA), followed Cy3 labeling of cDNA with the Roche NimbleGen One Color DNA Labeling Kit. NimbleGen 12x135K gene expression microarrays (design: 120419_Cynomolgus_v5_TH_exp_HX12) were hybridized with 4 μg of Cy3-labeled cDNA for 16 h at 42 °C and were washed and dried according to the manufacturer’s instruction. Microarray data were collected by confocal scanning using the Roche NimbleGen MS200 Microarray scanner at 2 μm pixel resolution (Roche NimbleGen, Inc., Madison, WI, USA). NimbleGen probe intensities were subjected to Robust Multi-Array Analysis (RMA) with background correction and quantile normalization as implemented in the NimbleScan Software, version 2.6 (Roche NimbleGen, Inc., Madison, WI). Averaged gene-level signal intensities were summarized into gene calls and log2 transformed. To summarize the major biology of genes, functional annotation was performed using the DAVID knowledgebase version 6.8 (https://david.ncifcrf.gov/).
Comparative genomic hybridization arrays
Cynomolgus monkey spleen tissues were homogenized in tubes prefilled with 1.4 mm ceramic beads and QiaGen’s lysis reagent ALT using a FastPrep-24 instrument (MP Biomedicals, Solon, OH, USA) and then incubated with Proteinase K at 55 °C for 1 h followed by RNAse A treatment at 25 °C for 2 min (Qiagen Inc., Valencia, CA). Cynomolgus monkey blood specimens (200 μl) were incubated at 70 °C for 10 min in QiaGen’s lysis reagent ALT with Proteinase K and RNAse A. Genomic DNA from lysates was extracted using the QIAamp Mini kit (Qiagen Inc., Valencia, CA, USA). Assessment of unfragmented, high molecular weight DNA and quantification was performed using microfluidic chip analysis on an Agilent 2100 bioanalyzer (Agilent Technologies Inc., Santa Clara, CA, USA). 0.5 μg of DNA from one animal tissue at a time and 0.5 μg of reference DNA - pooled DNA from blood specimens of 25 male Cynomolgus monkeys — were used for labeling by an isothermal Klenow fill-in reaction with either Cy3 or Cy5 random nonamer primer using the Roche NimbleGen Dual color labeling kit (Roche NimbleGen, Inc., Madison, WI). Labeling Hybridization Controls were spiked-in as quality controls for copy number variation detection (Roche NimbleGen, Inc., Madison, WI).
NimbleGen 4.2 M CGH microarrays (design: 120405_Cynomolgus5_CGH_UX1) were hybridized with 34 μg of Cy3- and 34 μg of Cy5-labeled DNA for 72 h at 42 °C. After hybridization, microarrays were washed and dried according to the manufacturer’s instruction, whereat 150 mM 1,3,5-Triaza-7 phospha-adamantane was included in the last washing step to avoid interference of ozone with the Cy5 dye during drying and scanning. Microarray data were collected by confocal scanning using the Roche NimbleGen MS200 Microarray Scanner at 2 μm pixel resolution (Roche NimbleGen, Inc., Madison, WI, USA).
SNR i denotes the signal-to-noise ratio of a given probe i, and σ Ci the standard deviation of all probes on the same chromosome as probe i. The average SNR of all CNV probes per CGH array was used as metric to evaluate the normalization performance. Using this approach, a fraction of 4000 probes per model fitting step resulted in the largest median SNR improvement (1.1%) and was therefore chosen (Additional file 1: Figure S4, Table S5).
CNV calling was performed with 3 inherently different approaches to mitigate method specific errors: R-GADA  was used with the following parameters: alpha = 0.2, T = 4.5, minseglen = 5. DNAcopy  was used with minseglen = 5, undosd = 3, undoprune = 0.05 and data smoothing was applied prior to CNV calling. CopyMap  was used with r = 20, T = 4, m = 5, a = 2.1, P = 0.001. Furthermore, for R-GADA and DNAcopy z-scores were calculated for all CNV calls based on the mean log2 − ratio of the CNV, and only CNVs with z-scores >1.5 or <−1.5 retained. For CNVs called by CopyMap a carrier probability of at least 0.8 was required.
The three obtained CNV calling profiles per individual were then merged and only CNVs called by at least two methods were kept, and loci with conflicting copy number states were removed. These resulting profiles were then further merged between individuals to obtain CNV regions that could be genotyped across individuals. In cases where an individual carried more than one CNV in a CNV region, the locus was marked as a complex locus and removed from subsequent steps. Additionally CNV loci located on the X chromosome or within array probe gaps larger than 500 kb + − 250 kb (e.g. centromeres) were removed. To avoid potential calling mistakes the median log2 − ratio of each CNV was used as genotype rather than the discrete copy number state provided by the CNV calling methods. PCA was performed with the R-function princomp to visualize and evaluate the inferred CNV genotypes together with hierarchical clustering of the Pearson correlation between log2 − ratios of samples.
Where E it is the measured expression level of a gene E in tissue t for the i th individual, C i is the genotype of a proximal CNV C for the i th individual and ε it a random error. To account for systematic confounding variation between samples, the loadings of the first principal component for individual i for the expression levels of all genes of the same tissue (G it ) was added as covariate. This was done to account for global changes of the transcriptome resulting from batch effects, population stratification or metabolic state of an individual, which could lead to spurious associations. Multiple combinations of principal component covariates were tested and inspection of the results (data not shown) revealed the most credible eQTLs when using the loadings for the first principal component for gene expression levels as covariate. An adapted version of the fastQTL software  was used to test all possible associations using this model. Correction for multiple testing was carried out in two steps, where first local permutations were applied to correct for multiple variants per gene  and then the false discovery rate (FDR) (q-value R-package, Storey J., 2015) was calculated per tissue to account for multiple tested genes. Only eQTLs below an FDR of 10% were considered as significant.
To further investigate the impact of CNVs on the gene expression landscape, genes within the regions of detected eQTLs were investigated for further associations with the eQTL CNV. The expression levels of all genes within 1 Mb from the TSS of an eQTL gene were tested for an association with the eQTL CNV with the same linear model as used for eQTL mapping. Bonferroni correction was calculated for all tested associations per region and association with a corrected p-value <0.05 were considered significant.
This work was supported by Swiss National Science Foundation grant 31003A_129835 (A.R.). All authors affiliated with F. Hoffmann-La Roche Ltd did not receive any external support or funding. We thank Tom Albert, Benjamin Gifford, and Ezra Lyon from NimbleGen for their support in the construction of customized Cynomolgus monkey microarrays. Further we thank Dr. Eugenia Migliavacca for her advice in the data analysis. Finally we are grateful to Dr. Thomas Singer from Roche for scientific and managerial support during the entire study.
Availability of data and materials
The microarray data from this study have been deposited at the NCBI Gene Expression Omnibus (GEO) (http://www.ncbi.nlm.nih.gov/geo/) under accession numbers GSE76560 for gene expression data and GSE76707 for array CGH data.
TH designed the Cynomolgus monkey microarrays, AS performed most of the experimental work, and AG performed the data analysis. TH, UC, and AR managed the entire project. AG, TH, UC, AR wrote the manuscript. All authors read and approved the final manuscript.
The authors declare that they have no competing interests.
Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
- Scherer SW, Lee C, Birney E, Altshuler DM, Eichler EE, Carter NP, et al. Challenges and standards in integrating surveys of structural variation. Nat Genet. 2007;39:S7–S15.View ArticlePubMedPubMed CentralGoogle Scholar
- Redon R, Ishikawa S, Fitch KR, Feuk L, Perry GH, Andrews TD, et al. Global variation in copy number in the human genome. Nature. 2006;444:444–54.View ArticlePubMedPubMed CentralGoogle Scholar
- Gokcumen O, Lee C. Copy number variants (CNVs) in primate species using array-based comparative genomic hybridization. Methods. 2009;49:18–25.View ArticlePubMedPubMed CentralGoogle Scholar
- Henrichsen CN, Vinckenbosch N, Zöllner S, Chaignat E, Pradervand S, Schütz F, et al. Segmental copy number variation shapes tissue transcriptomes. Nat Genet. 2009;41:424–9.View ArticlePubMedGoogle Scholar
- Cardoso-Moreira M, Arguello JR, Clark AG. Mutation spectrum of Drosophila CNVs revealed by breakpoint sequencing. Genome Biol. 2012;13:R119.View ArticlePubMedPubMed CentralGoogle Scholar
- MacDonald JR, Ziman R, Yuen RKC, Feuk L, Scherer SW. The Database of Genomic Variants: a curated collection of structural variation in the human genome. Nucleic Acids Res. 2014;42:D986–92.View ArticlePubMedGoogle Scholar
- Conrad DF, Pinto D, Redon R, Feuk L, Gokcumen O, Zhang Y, et al. Origins and functional impact of copy number variation in the human genome. Nature. 2010;464:704–12.View ArticlePubMedGoogle Scholar
- Lupski JR. Genomic rearrangements and sporadic disease. Nat Genet. 2007;39:S43–7.View ArticlePubMedGoogle Scholar
- Henrichsen CN, Chaignat E, Reymond A. Copy number variants, diseases and gene expression. Hum Mol Genet. 2009;18:R1–8.View ArticlePubMedGoogle Scholar
- Chaignat E, Yahya-Graison EA, Henrichsen CN, Chrast J, Schütz F, Pradervand S, et al. Copy number variation modifies expression time courses. Genome Res. 2011;21:106–13.View ArticlePubMedPubMed CentralGoogle Scholar
- Cahan P, Li Y, Izumi M, Graubert TA. The impact of copy number variation on local gene expression in mouse hematopoietic stem and progenitor cells. Nat Genet. 2009;41:430–7.View ArticlePubMedPubMed CentralGoogle Scholar
- Merla G, Howald C, Henrichsen CN, Lyle R, Wyss C, Zabot M-T, et al. Submicroscopic deletion in patients with Williams-Beuren syndrome influences expression levels of the nonhemizygous flanking genes. Am J Hum Genet. 2006;79:332–41.View ArticlePubMedPubMed CentralGoogle Scholar
- Stranger BE, Forrest MS, Dunning M, Ingle CE, Beazley C, Thorne N, et al. Relative impact of nucleotide and copy number variation on gene expression phenotypes. Science. 2007;315:848–53.View ArticlePubMedPubMed CentralGoogle Scholar
- Orozco LD, Cokus SJ, Ghazalpour A, Ingram-Drake L, Wang S, van Nas A, et al. Copy number variation influences gene expression and metabolic traits in mice. Hum Mol Genet. 2009;18:4118–29.View ArticlePubMedPubMed CentralGoogle Scholar
- Ricard G, Molina J, Chrast J, Gu W, Gheldof N, Pradervand S, et al. Phenotypic consequences of copy number variation: insights from Smith-Magenis and Potocki-Lupski syndrome mouse models. Plos Biol. 2010;8:e1000543.View ArticlePubMedPubMed CentralGoogle Scholar
- Usher CL, McCarroll SA. Complex and multi-allelic copy number variation in human disease. Brief Funct Genomics. 2015;14:329–38.View ArticlePubMedPubMed CentralGoogle Scholar
- Gonzaga-Jauregui C, Lupski JR, Gibbs RA. Human genome sequencing in health and disease. Annu Rev Med. 2012;63:35–61.View ArticlePubMedPubMed CentralGoogle Scholar
- Girirajan S, Rosenfeld JA, Coe BP, Parikh S, Friedman N, Goldstein A, et al. Phenotypic heterogeneity of genomic disorders and rare copy-number variants. N Engl J Med. 2012;367:1321–31.View ArticlePubMedPubMed CentralGoogle Scholar
- Männik K, Mägi R, Macé A, Cole B, Guyatt AL, Shihab HA, et al. Copy number variations and cognitive phenotypes in unselected populations. JAMA. 2015;313:2044–54.View ArticlePubMedPubMed CentralGoogle Scholar
- Johansson I, Ingelman-Sundberg M. Genetic polymorphism and toxicology—with emphasis on cytochrome p450. Toxicol Sci. 2011;120:1–13.View ArticlePubMedGoogle Scholar
- Ebeling M, Küng E, See A, Broger C, Steiner G, Berrera M, et al. Genome-based analysis of the nonhuman primate Macaca fascicularis as a model for drug safety assessment. Genome Res. 2011;21:1746–56.View ArticlePubMedPubMed CentralGoogle Scholar
- Shively CA, Clarkson TB. The unique value of primate models in translational research. Nonhuman primate models of women’s health: introduction and overview. Shively CA, Clarkson TB, editors. Am J Primatol. 2009;71:715–21.View ArticlePubMedGoogle Scholar
- Haus T, Ferguson B, Rogers J, Doxiadis G, Certa U, Rose NJ, et al. Genome typing of nonhuman primate models: implications for biomedical research. Trends Genet. 2014;30:482–7.View ArticlePubMedGoogle Scholar
- Stevison LS, Kohn MH. Determining genetic background in captive stocks of cynomolgus macaques (Macaca fascicularis). J Med Primatol. 2008;37:311–7.PubMedGoogle Scholar
- Liu YY, Wang Z, Thomas J, Goodwin KJ, Stavrou S, Neville DM. Polymorphisms of CD3epsilon in cynomolgus and rhesus monkeys and their relevance to anti-CD3 antibodies and immunotoxins. Immunol Cell Biol. 2007;85:357–62.View ArticlePubMedGoogle Scholar
- Menninger K, Wieczorek G, Riesen S, Kunkler A, Audet M, Blancher A, et al. The origin of cynomolgus monkey affects the outcome of kidney allografts under Neoral immunosuppression. Transplant Proc. 2002;34:2887–8.View ArticlePubMedGoogle Scholar
- Drevon-Gaillot E, Perron-Lepage M-F, Clément C, Burnett R. A review of background findings in cynomolgus monkeys (Macaca fascicularis) from three different geographical origins. Exp Toxicol Pathol. 2006;58:77–88.View ArticlePubMedGoogle Scholar
- Heckel T, Singh A, Gschwind A, Reymond A, Certa U. Genetic Variations in the Macaca fascicularis Genome Related to Biomedical Research. In: Bluemel J, Korte S, Schenck E, Weinbauer G, editors. The Nonhuman Primate in Nonclinical Drug Development and Safety Assessment. San Diego: Academic Press; 2015. p. 53–64.Google Scholar
- Yan G, Zhang G, Fang X, Zhang Y, Li C, Ling F, et al. Genome sequencing and comparison of two nonhuman primate animal models, the cynomolgus and Chinese rhesus macaques. Nat Biotechnol. 2011;29:1019–89.View ArticlePubMedGoogle Scholar
- Ericsen AJ, Starrett GJ, Greene JM, Lauck M, Raveendran M, Deiros DR, et al. Whole genome sequencing of SIV-infected macaques identifies candidate loci that may contribute to host control of virus replication. Genome Biol. 2014;15:478.View ArticlePubMedPubMed CentralGoogle Scholar
- Osada N, Hettiarachchi N, Adeyemi Babarinde I, Saitou N, Blancher A. Whole-genome sequencing of six Mauritian Cynomolgus macaques (Macaca fascicularis) reveals a genome-wide pattern of polymorphisms under extreme population bottleneck. Genome Biol Evol. 2015;7:821–30.View ArticlePubMedGoogle Scholar
- Higashino A, Sakate R, Kameoka Y, Takahashi I, Hirata M, Tanuma R, et al. Whole-genome sequencing and analysis of the Malaysian cynomolgus macaque (Macaca fascicularis) genome. Genome Biol. 2012;13:R58.View ArticlePubMedPubMed CentralGoogle Scholar
- Marioni JC, Thorne NP, Valsesia A, Fitzgerald T, Redon R, Fiegler H, et al. Breaking the waves: improved detection of copy number variation from microarray-based comparative genomic hybridization. Genome Biol. 2007;8:R228.View ArticlePubMedPubMed CentralGoogle Scholar
- Pique-Regi R, Cáceres A, González JR. R-Gada: a fast and flexible pipeline for copy number analysis in association studies. BMC Bioinformatics. 2010;11:380.View ArticlePubMedPubMed CentralGoogle Scholar
- Venkatraman ES, Olshen AB. A faster circular binary segmentation algorithm for the analysis of array CGH data. Bioinformatics. 2007;23:657–63.View ArticlePubMedGoogle Scholar
- Zöllner S. CopyMap: localization and calling of copy number variation by joint analysis of hybridization data from multiple individuals. Bioinformatics. 2010;26:2776–7.View ArticlePubMedPubMed CentralGoogle Scholar
- Sudmant PH, Huddleston J, Catacchio CR, Malig M, Hillier LW, Baker C, et al. Evolution and diversity of copy number variation in the great ape lineage. Genome Res. 2013;23:1373–82.View ArticlePubMedPubMed CentralGoogle Scholar
- Sudmant PH, Mallick S, Nelson BJ, Hormozdiari F, Krumm N, Huddleston J, et al. Global diversity, population stratification, and selection of human copy-number variation. Science. 2015. doi:10.1126/science.aab3761.
- Ongen H, Buil A, Brown A, Dermitzakis E, Delaneau O. Fast and efficient QTL mapper for thousands of molecular phenotypes. Bioinformatics. 2015;32:1479–85.Google Scholar
- Huang DW, Sherman BT, Lempicki RA. Systematic and integrative analysis of large gene lists using DAVID bioinformatics resources. Nat Protoc. 2009;4:44–57.View ArticleGoogle Scholar
- Thurman RE, Rynes E, Humbert R, Vierstra J, Maurano MT, Haugen E, et al. The accessible chromatin landscape of the human genome. Nature. 2012;489:75–82.View ArticlePubMedPubMed CentralGoogle Scholar
- Itsara A, Wu H, Smith JD, Nickerson DA, Romieu I, London SJ, et al. De novo rates and selection of large copy number variation. Genome Res. 2010;20:1469–81.View ArticlePubMedPubMed CentralGoogle Scholar
- Stranger BE, Nica AC, Forrest MS, Dimas A, Bird CP, Beazley C, et al. Population genomics of human gene expression. Nat Genet. 2007;39:1217–24.View ArticlePubMedPubMed CentralGoogle Scholar
- Stranger BE, Montgomery SB, Dimas AS, Parts L, Stegle O, Ingle CE, et al. Patterns of cis regulatory variation in diverse human populations. Barsh GS, editor. PLoS Genet. 2012. doi:10.1371/journal.pgen.1002639.
- Dimas AS, Deutsch S, Stranger BE, Montgomery SB, Borel C, Attar-Cohen H, et al. Common regulatory variation impacts gene expression in a cell type-dependent manner. Science. 2009;325:1246–50.View ArticlePubMedPubMed CentralGoogle Scholar
- ENCODE Project Consortium, Birney E, Stamatoyannopoulos JA, Dutta A, Guigo R, Gingeras TR, et al. Identification and analysis of functional elements in 1% of the human genome by the ENCODE pilot project. Nature. 2007;447:799–816.View ArticleGoogle Scholar
- Nguyen D-Q, Webber C, Ponting CP. Bias of selection on human copy-number variants. PLoS Genet. 2006. doi:10.1371/journal.pgen.0020020.
- Pluznick JL, Zou D-J, Zhang X, Yan Q, Rodriguez-Gil DJ, Eisner C, et al. Functional expression of the olfactory signaling system in the kidney. Proc Natl Acad Sci U S A. 2009;106:2059–64.View ArticlePubMedPubMed CentralGoogle Scholar
- Pluznick JL, Protzko RJ, Gevorgyan H, Peterlin Z, Sipos A, Han J, et al. Olfactory receptor responding to gut microbiota-derived signals plays a role in renin secretion and blood pressure regulation. Proc Natl Acad Sci U S A. 2013;110:4410–5.View ArticlePubMedPubMed CentralGoogle Scholar
- Kiehl S, Herkt SC, Richter AM, Fuhrmann L, El-Nikhely N, Seeger W, et al. ABCB4 is frequently epigenetically silenced in human cancers and inhibits tumor growth. Sci Rep. 2014;4:6899.View ArticlePubMedPubMed CentralGoogle Scholar