Copy Number Variation of CCL3-like Genes Affects Rate of Progression to Simian-AIDS in Rhesus Macaques (Macaca mulatta)

Variation in genes underlying host immunity can lead to marked differences in susceptibility to HIV infection among humans. Despite heavy reliance on non-human primates as models for HIV/AIDS, little is known about which host factors are shared and which are unique to a given primate lineage. Here, we investigate whether copy number variation (CNV) at CCL3-like genes (CCL3L), a key genetic host factor for HIV/AIDS susceptibility and cell-mediated immune response in humans, is also a determinant of time until onset of simian-AIDS in rhesus macaques. Using a retrospective study of 57 rhesus macaques experimentally infected with SIVmac, we find that CCL3L CNV explains approximately 18% of the variance in time to simian-AIDS (p<0.001) with lower CCL3L copy number associating with more rapid disease course. We also find that CCL3L copy number varies significantly (p<10⁻⁶) among rhesus subpopulations, with Indian-origin macaques having, on average, half as many CCL3L gene copies as Chinese-origin macaques. Lastly, we confirm that CCL3L shows variable copy number in humans and chimpanzees and report on CCL3L CNV within and among three additional primate species. On the basis of our findings we suggest that (1) the difference in population level copy number may explain previously reported observations of longer post-infection survivorship of Chinese-origin rhesus macaques, (2) stratification by CCL3L copy number in rhesus SIV vaccine trials will increase power and reduce noise due to non-vaccine-related differences in survival, and (3) CCL3L CNV is an ancestral component of the primate immune response and, therefore, copy number variation has not been driven by HIV or SIV per se.     

Contributions of Mamu-A*01 Status and TRIM5 Allele Expression,But Not CCL3L Copy Number Variation,to the Control of SIVmac251 Replication in Indian-Origin Rhesus Monkeys

CCL3 is a ligand for the HIV-1 co-receptor CCR5. There have recently been conflicting reports in the literature concerning whether CCL3-like gene (CCL3L) copy number variation (CNV) is associated with resistance to HIV-1 acquisition and with both viral load and disease progression following infection with HIV-1. An association has also been reported between CCL3L CNV and clinical sequelae of the simian immunodeficiency virus (SIV) infection in vivo in rhesus monkeys. The present study was initiated to explore the possibility of an association of CCL3L CNV with the control of virus replication and AIDS progression in a carefully defined cohort of SIVmac251-infected, Indian-origin rhesus monkeys. Although we demonstrated extensive variation in copy number of CCL3L in this cohort of monkeys, CCL3L CNV was not significantly associated with either peak or set-point plasma SIV RNA levels in these monkeys when MHC class I allele Mamu-A*01 was included in the models or progression to AIDS in these monkeys. With 66 monkeys in the study, there was adequate power for these tests if the correlation of CCL3L and either peak or set-point plasma SIV RNA levels was 0.34 or 0.36, respectively. These findings call into question the premise that CCL3L CNV is important in HIV/SIV pathogenesis.     

CCL3L1 copy number variation and susceptibility to HIV-1 infection: a meta-analysis

Background

Although several studies have investigated whether CCL3L1 copy number variation (CNV) influences the risk of HIV-1 infection, there are still no clear conclusions. Therefore, we performed a meta-analysis using two models to generate a more robust estimate of the association between CCL3L1 CNV and susceptibility to HIV-1 infection.

Methods

We divided the cases and controls into two parts as individuals with CCL3L1 gene copy number (GCN) above the population specific median copy number (PMN) and individuals with CCL3L1 GCN below PMN, respectively. Odds ratios (ORs) with 95% confidence intervals (95% CIs) were given for the main analysis. We also conducted stratified analyses by ethnicity, age group and sample size. Relevant literatures were searched through PubMed and ISI Web of Knowledge up to March 2010.

Results

In total, 9 studies with 2434 cases and 4029 controls were included. ORs for the main analysis were 1.35 (95% CI, 1.02–1.78, model: GCN ≤ PMN Vs. GCN > PMN) and 1.70 (95% CI, 1.30–2.23, model: GCN < PMN Vs. GCN ≥ PMN), respectively. Either in stratified analysis, statistically significant results can be detected in some subgroups.

Conclusions

Our analyses indicate that CCL3L1 CNV is associated with susceptibility to HIV-1 infection. A lower copy number is associated with an increased risk of HIV-1 infection, while a higher copy number is associated with reduced risk for acquiring HIV-1.     

Copy number variations of <Emphasis Type="Italic">CCL3L1</Emphasis> and long-term prognosis of HIV-1 infection in asymptomatic HIV-infected Japanese with hemophilia

We set up a cohort of HIV-infected, asymptomatic Japanese patients with hemophilia for follow-up study in 1995. All subjects who had been infected with HIV-1 for more than 10 years met the criteria for long-term nonprogressors (LTNPs) at the time of entry; however, some of them later developed lymphopenia and required antiretroviral treatment during five more years of observation. In this study, we investigated the impacts of the CCL3L1 dose on the long-term prognosis in the subjects with chronic HIV-1 infection. We collected genomic DNA from 95 long-term survivors including 48 nonprogressors and 47 subjects receiving antiretroviral treatment. The distributions of CCL3L1 copy number significantly differed between the 95 HIV-1-infected subjects with hemophilia and 205 controls. Average copy number of CCL3L1 in the HIV-1-infected subjects was significantly lower than in control (5.00 ± 0.22 vs 3.35 ± 0.24, p < 0.001). Moreover, the subjects possessing two or less copies of CCL3L1 had significantly higher risk of acquiring HIV-1. However, CCL3L1 copy number variations had no significant effect on the disease progression among the LTNP subjects who had been afflicted with chronic HIV-1 infection for more than 15 years, when compared between nonprogressors and patients under treatment (3.68 ± 0.37 vs 3.02 ± 0.29, ns). Furthermore, variations in the CCL3L1 copy number had little effect on the levels of HIV-1 load among them. We conclude that variation in the CCL3L1 copy number is apparently not a factor that determines the prognosis of chronic HIV-1 infection, even though it is linked to HIV-1 susceptibility.     

Serotonin Transporter Genotype Modulates Social Reward and Punishment in Rhesus Macaques

Background

Serotonin signaling influences social behavior in both human and nonhuman primates. In humans, variation upstream of the promoter region of the serotonin transporter gene (5-HTTLPR) has recently been shown to influence both behavioral measures of social anxiety and amygdala response to social threats. Here we show that length polymorphisms in 5-HTTLPR predict social reward and punishment in rhesus macaques, a species in which 5-HTTLPR variation is analogous to that of humans.

Methodology/Principal Findings

In contrast to monkeys with two copies of the long allele (L/L), monkeys with one copy of the short allele of this gene (S/L) spent less time gazing at face than non-face images, less time looking in the eye region of faces, and had larger pupil diameters when gazing at photos of a high versus low status male macaques. Moreover, in a novel primed gambling task, presentation of photos of high status male macaques promoted risk-aversion in S/L monkeys but promoted risk-seeking in L/L monkeys. Finally, as measured by a “pay-per-view” task, S/L monkeys required juice payment to view photos of high status males, whereas L/L monkeys sacrificed fluid to see the same photos.

Conclusions/Significance

These data indicate that genetic variation in serotonin function contributes to social reward and punishment in rhesus macaques, and thus shapes social behavior in humans and rhesus macaques alike.     

Multiplex Paralogue Ratio Tests for accurate measurement of multiallelic CNVs

Demonstrating an association between a polymorphism and a disease phenotype through case-control studies requires reliable large-scale genotyping, but accurate measurement of copy number variation has proven to be technically challenging. Here we build on our previous experience with Paralogue Ratio Tests (PRT) to develop PRT copy number determination at the CCL3L1/CCL4L1 copy number variant. A multiplex PRT assay based on four independent comparative PCRs results in a convenient, accurate and robust method of multiallelic copy number measurement suitable for use in large-scale case-control studies, which can unambiguously assign virtually all samples tested to discrete copy number classes.     

Variation of CNV distribution in five different ethnic populations

Recent studies have revealed a new type of variation in the human genome encompassing relatively large genomic segments ( approximately 100 kb-2.5 Mb), commonly referred to as copy number variation (CNV). The full nature and extent of CNV and its frequency in different ethnic populations is still largely unknown. In this study we surveyed a set of 12 CNVs previously detected by array-CGH. More than 300 individuals from five different ethnic populations, including three distinct European, one Asian and one African population, were tested for the occurrence of CNV using multiplex ligation-dependent probe amplification (MLPA). Seven of these loci indeed showed CNV, i.e., showed copy numbers that deviated from the population median. More precise estimations of the actual genomic copy numbers for (part of) the NSF gene locus, revealed copy numbers ranging from two to at least seven. Additionally, significant inter-population differences in the distribution of these copy numbers were observed. These data suggest that insight into absolute DNA copy numbers for loci exhibiting CNV is required to determine their potential contribution to normal phenotypic variation and, in addition, disease susceptibility.     

CCL3L1 and CCL4L1 chemokine genes are located in a segmental duplication at chromosome 17q12

Sixteen CC chemokine genes localize to a 2.06-Mb interval at 17q11.2-q12 on genomic contig NT_010799.13. Four of these genes comprise two closely related paralogous pairs: CCL3-CCL3L1 and CCL4-CCL4L1. Members within each pair share 95% sequence identity at both the genomic and the amino acid levels. One BAC clone (AC131056.5) on the contig with substantial internal sequence duplication contains two complete copies of CCL3L1 and CCL4L1 and one truncated copy of CCL3L1, while a partially overlapping clone (AC003976.1) contains one copy each of CCL3 and CCL4. Dot-matrix comparison of the regions of AC131056.5 with those of AC003976.1 containing the four genes reveals 90% sequence similarity over 37 kb. These observations support the idea that the multiple copies of CCL3L1 and CCL4L1 present in a single diploid genome are the result of segmental duplication.     

A Large‐Scale SNP‐Based Genomic Admixture Analysis of the Captive Rhesus Macaque Colony at the California National Primate Research Center

Some breeding facilities in the United States have crossbred Chinese and Indian rhesus macaque (Macaca mulatta) founders either purposefully or inadvertently. Genetic variation that reflects geographic origins among research subjects has the potential to influence experimental outcomes. The use of animals from different geographic regions, their hybrids, and animals of varying degrees of kinship in an experiment can obscure treatment effects under study because high interanimal genetic variance can increase phenotypic variance among the research subjects. The intent of this study, based on a broad genomic analysis of 2,808 single nucleotide polymorphisms (SNPs), is to ensure that only animals estimated to be of pure Indian or Chinese ancestry, based on both demographic and genetic information, are used as sources of infants for derivation and expansion of the California National Primate Research Center's (CNPRC) super‐Specific Pathogen Free (SSPF) rhesus macaque colony. Studies of short tandem repeats (STRs) in Indian and Chinese rhesus macaques have reported that heterozygosity of STRs is higher in Chinese rhesus macaques than in Indian rhesus macaques. The present study shows that heterozygosity of SNPs is actually higher in Indian than in Chinese rhesus macaques and that the Chinese SSPF rhesus macaque colony is far less differentiated from their founders compared to the Indian‐origin animals. The results also reveal no evidence of recent gene flow from long‐tailed and pig‐tailed macaques into the source populations of the SSPF rhesus macaques. This study indicates that many of the long‐tailed macaques held in the CNPRC are closely related individuals. Most polymorphisms shared among the captive rhesus, long‐tailed, and pig‐tailed macaques likely predate the divergence among these groups. Am. J. Primatol. 74:747‐757, 2012. © 2012 Wiley Periodicals, Inc.     

HaplotypeCN: Copy Number Haplotype Inference with Hidden Markov Model and Localized Haplotype Clustering

Copy number variation (CNV) has been reported to be associated with disease and various cancers. Hence, identifying the accurate position and the type of CNV is currently a critical issue. There are many tools targeting on detecting CNV regions, constructing haplotype phases on CNV regions, or estimating the numerical copy numbers. However, none of them can do all of the three tasks at the same time. This paper presents a method based on Hidden Markov Model to detect parent specific copy number change on both chromosomes with signals from SNP arrays. A haplotype tree is constructed with dynamic branch merging to model the transition of the copy number status of the two alleles assessed at each SNP locus. The emission models are constructed for the genotypes formed with the two haplotypes. The proposed method can provide the segmentation points of the CNV regions as well as the haplotype phasing for the allelic status on each chromosome. The estimated copy numbers are provided as fractional numbers, which can accommodate the somatic mutation in cancer specimens that usually consist of heterogeneous cell populations. The algorithm is evaluated on simulated data and the previously published regions of CNV of the 270 HapMap individuals. The results were compared with five popular methods: PennCNV, genoCN, COKGEN, QuantiSNP and cnvHap. The application on oral cancer samples demonstrates how the proposed method can facilitate clinical association studies. The proposed algorithm exhibits comparable sensitivity of the CNV regions to the best algorithm in our genome-wide study and demonstrates the highest detection rate in SNP dense regions. In addition, we provide better haplotype phasing accuracy than similar approaches. The clinical association carried out with our fractional estimate of copy numbers in the cancer samples provides better detection power than that with integer copy number states.     

Development of a Chinese-Indian hybrid (Chindian) rhesus macaque colony at the California National Primate Research Center by introgression

Background  Fullbred Chinese and Indian rhesus macaques represent genetically distinct populations. The California National Primate Research Center introduced Chinese founders into its Indian-derived rhesus colony in response to the 1978 Indian embargo on exportation of animals for research and the concern that loss of genetic variation in the closed colony would hamper research efforts. The resulting hybrid rhesus now number well over a thousand animals and represent a growing proportion of the animals in the colony.
Methods  We characterized the population genetic structure of the hybrid colony and compared it with that of their pure Indian and Chinese progenitors.
Results  The hybrid population contains higher genetic diversity and linkage disequilibrium than their full Indian progenitors and represents a resource with unique research applications.
Conclusions  The genetic diversity of the hybrids indicates that the strategy to introduce novel genes into the colony by hybridizing Chinese founders and their hybrid offspring with Indian-derived animals was successful.     

Relationship between complement membrane attack complex, chemokine (C-C motif) ligand 2 (CCL2) and vascular endothelial growth factor in mouse model of laser-induced choroidal neovascularization

The present study investigated the interactions among the complement membrane attack complex (MAC), CCL2, and VEGF that occur in vivo during the development of choroidal neovascularization (CNV). We first investigated the sequential expression of MAC, CCL2, and VEGF during laser-induced CNV in C57BL/6 mice. Increased MAC deposition was detected at 1 h, CCL2 increased at 3 h, and VEGF was up-regulated at day 3 post-laser treatment. These results suggested that during laser-induced CNV, MAC, CCL2 and VEGF are formed and/or expressed in the following order: MAC → CCL2 → VEGF. To determine the cross-talk between MAC, CCL2, and VEGF during laser-induced CNV, neutralizing antibodies were injected both systemically and locally to block the bioactivity of each molecule. Blocking MAC formation inhibited CCL2 and VEGF expression and also limited CNV formation, whereas neutralization of CCL2 bioactivity did not affect MAC deposition; however, it reduced VEGF expression and CNV formation. When bioactivity of VEGF was blocked, CNV formation was significantly inhibited, but MAC deposition was not affected. Together, our results demonstrate that MAC is an upstream mediator and effect of MAC on the development of laser-induced CNV can be attributed to its direct effect on VEGF as well as its effect on VEGF that is mediated by CCL2. Understanding the interplay between immune mediators is critical to gain insight into the pathogenesis of CNV.     

Human gene copy number variation and infectious disease

Variability in the susceptibility to infectious disease and its clinical manifestation can be determined by variation in the environment and by genetic variation in the pathogen and the host. Despite several successes based on candidate gene studies, defining the host variation affecting infectious disease has not been as successful as for other multifactorial diseases. Both single nucleotide variation and copy number variation (CNV) of the host contribute to the host’s susceptibility to infectious disease. In this review we focus on CNV, particularly on complex multiallelic CNV that is often not well characterised either directly by hybridisation methods or indirectly by analysis of genotypes and flanking single nucleotide variants. We summarise the well-known examples, such as α-globin deletion and susceptibility to severe malaria, as well as more recent controversies, such as the extensive CNV of the chemokine gene CCL3L1 and HIV infection. We discuss the potential biological mechanisms that could underly any genetic association and reflect on the extensive complexity and functional variation generated by a combination of CNV and sequence variation, as illustrated by the Fc gamma receptor genes FCGR3A, FCGR3B and FCGR2C. We also highlight some understudied areas that might prove fruitful areas for further research.     

CCL3L1 dose and HIV-1 susceptibility

Several genetic factors influence HIV-1 susceptibility or AIDS disease progression. A recent study reported on what could be a particularly important genetic determinant for HIV-1 susceptibility and disease progression: copy number of a chemokine gene termed CCL3L1. Individuals with low copy numbers of the gene, relative to their ethnic background, were associated with markedly enhanced HIV-1/AIDS susceptibility. These findings define an important new genetic determinant of HIV-1 susceptibility and further emphasize the importance of the chemokine system, either as elements that inhibit HIV-1 infection or that modulate antiviral immune responses.     

MOCSphaser: a haplotype inference tool from a mixture of copy number variation and single nucleotide polymorphism data

Detailed analyses of the population-genetic nature of copy number variations (CNVs) and the linkage disequilibrium between CNV and single nucleotide polymorphism (SNP) loci from high-throughput experimental data require a computational tool to accurately infer alleles of CNVs and haplotypes composed of both CNV alleles and SNP alleles. Here we developed a new tool to infer population frequencies of such alleles and haplotypes from observed copy numbers and SNP genotypes, using the expectation-maximization algorithm. This tool can also handle copy numbers ambiguously determined, such as 2 or 3 copies, due to experimental noise. AVAILABILITY: http://emu.src.riken.jp/MOCSphaser/MOCSphaser.zip.     

Determination of beta-defensin genomic copy number in different populations: a comparison of three methods

Background

There have been conflicting reports in the literature on association of gene copy number with disease, including CCL3L1 and HIV susceptibility, and β-defensins and Crohn''s disease. Quantification of precise gene copy numbers is important in order to define any association of gene copy number with disease. At present, real-time quantitative PCR (QPCR) is the most commonly used method to determine gene copy number, however the Paralogue Ratio Test (PRT) is being used in more and more laboratories.

Findings

In this study we compare a Pyrosequencing-based Paralogue Ratio Test (PPRT) for determining beta-defensin gene copy number with two currently used methods for gene copy number determination, QPCR and triplex PRT by typing five different cohorts (UK, Danish, Portuguese, Ghanaian and Czech) of DNA from a total of 576 healthy individuals. We found a systematic measurement bias between DNA cohorts revealed by QPCR, but not by the PRT-based methods. Using PRT, copy number ranged from 2 to 9 copies, with a modal copy number of 4 in all populations.

Conclusions

QPCR is very sensitive to quality of the template DNA, generating systematic biases that could produce false-positive or negative disease associations. Both triplex PRT and PPRT do not show this systematic bias, and type copy number within the correct range, although triplex PRT appears to be a more precise and accurate method to type beta-defensin copy number.     

Comparison of the frequencies of major histocompatibility (MHC) class-II DQA1 and DQB1 alleles in Indian and Chinese rhesus macaques (Macaca mulatta)

The major histocompatibility complex (MHC) comprises related gene families, some of which are highly polymorphic, whose protein products mediate immune response. Rhesus macaques (Macaca mulatta) are a vital animal model for research in human diseases and are native to regions extending from Afghanistan in the west to the Eastern Plains of China and from Peking to the north, southward through islands of Southeast Asia. The distributions of MHC class-II Mamu DQA1 and Mamu DQB1 alleles in two groups of domestically bred rhesus macaques of Indian and Chinese origin and the Mamu DQA1 genotypes of a small number of Burmese rhesus macaques were compared. Major allelic differences were observed between the Indian and Chinese rhesus macaques, and gene diversity decreased from east to west. These and other intra-specific genetic differences among regional populations of rhesus macaques might influence the outcome of biomedical research in which they are used as subjects, and illustrate the importance of completely genetically characterizing subjects used as animal models in biomedical research.     

Functional redundancy of the human CCL4 and CCL4L1 chemokine genes

CCL4 and CCL4L1 are two CC chemokine genes located at chromosome 17q21 whose mature proteins differ at only a single amino acid. Abundant functional information exists for CCL4, however, CCL4L1 has only recently been recognized as a distinct gene, thus information describing it is wanting. The CCL4L1 protein was synthesized in Escherichia coli and compared with the CCL4 protein. Competitive binding studies using HEK-293/CCR5 cells produced comparable EC50 values for the two proteins. Similarly, chemotaxis assays with cells expressing CCR1, CCR3, or CCR5 revealed no substantial differences. CCL4L1 was somewhat more effective at inhibiting HIV-1 replication in PBMCs than was CCL4, however the difference was not statistically significant. These data combined with the observation of individual variation in CCL4L1 gene copy number [Eur. J. Immunol. 32 (2002) 3016, Genomics 83 (2004) 735] support the contention that the CCL4 and CCL4L1 proteins have redundant functions.     

Identification of Country of Origin and Admixture Between Indian and Chinese Rhesus Macaques

We describe a restriction analysis that distinguishes between rhesus macaques of unmixed Indian and Chinese ancestry and between western and eastern Chinese ancestry. We amplified a 254-bp fragment of mitochondrial DNA (mtDNA) that contains restriction sites hypothesized to be diagnostic of country of origin for samples from 534 and 567 individuals alleged to be of solely Indian or solely Chinese origin, respectively. After digestion with the MaeIII, SmlI, and BccI restriction enzymes, the polymerase chain reaction (PCR) products of only 3 of the 1101 samples exhibited restriction patterns uncharacteristic of their alleged country of origin. A sample comprising 392 of these rhesus macaques was genotyped for 24 nuclear microsatellite (STR) loci. Principal coordinates analysis confirmed marked genetic similarity of regional populations within each country but a substantial difference between Indian and Chinese rhesus macaques. Using STRUCTURE (Pritchard and Wen, 2003),we assigned probabilities of Chinese and Indian ancestry to each sample based on its STR genotypes. We assigned all the unmixed rhesus macaques to their correct countries of origin with probabilities >0.95. We constructed an artificial sample of 1st-generation hybrid Indian/Chinese rhesus macaques by randomly sampling from the genotypes of Indian and Chinese individuals. STRUCTURE assigned robabilities of Chinese and Indian ancestry to hybrids that closely corresponded with the proportions of alleles in that sample drawn from unmixed Chinese and Indian rhesus macaques.     

Global copy number analyses by next generation sequencing provide insight into pig genome variation

Background

Copy number variations (CNVs) confer significant effects on genetic innovation and phenotypic variation. Previous CNV studies in swine seldom focused on in-depth characterization of global CNVs.

Results

Using whole-genome assembly comparison (WGAC) and whole-genome shotgun sequence detection (WSSD) approaches by next generation sequencing (NGS), we probed formation signatures of both segmental duplications (SDs) and individualized CNVs in an integrated fashion, building the finest resolution CNV and SD maps of pigs so far. We obtained copy number estimates of all protein-coding genes with copy number variation carried by individuals, and further confirmed two genes with high copy numbers in Meishan pigs through an enlarged population. We determined genome-wide CNV hotspots, which were significantly enriched in SD regions, suggesting evolution of CNV hotspots may be affected by ancestral SDs. Through systematically enrichment analyses based on simulations and bioinformatics analyses, we revealed CNV-related genes undergo a different selective constraint from those CNV-unrelated regions, and CNVs may be associated with or affect pig health and production performance under recent selection.

Conclusions

Our studies lay out one way for characterization of CNVs in the pig genome, provide insight into the pig genome variation and prompt CNV mechanisms studies when using pigs as biomedical models for human diseases.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-593) contains supplementary material, which is available to authorized users.