CCL3L1 copy number, HIV load, and immune reconstitution in sub-Saharan Africans

Copy number variation of CCL3L1 among three major ethnic groups in Malaysia

Background

C-C motif Chemokine Ligand 3 Like 1 (CCL3L1) is a multiallelic copy number variable, which plays a crucial role in immunoregulatory and hosts defense through the production of macrophage inflammatory protein (MIP)-1α. Variable range of the CCL3L1 copies from 0 to 14 copies have been documented in several different populations. However, there is still lack of report on the range of CCL3L1 copy number exclusively among Malaysians who are a multi-ethnic population. Thus, this study aims to extensively examine the distribution of CCL3L1 copy number in the three major populations from Malaysia namely Malay, Chinese and Indian. A diploid copy number of CCL3L1 for 393 Malaysians (Malay = 178, Indian = 90, and Chinese = 125) was quantified using Paralogue Ratio Tests (PRTs) and then validated with microsatellites analysis.

Results

To our knowledge, this is the first report on the CCL3L1 copy number that has been attempted among Malaysians and the Chinese ethnic group exhibits a diverse pattern of CCL3L1 distribution copy number from the Malay and Indian (p < 0.0001). The CCL3L1 ranged from 0 to 8 copies for both the Malay and Indian ethnic groups while 0 to 10 copies for the Chinese ethnic. Consequently, the CCL3L1 copy number among major ethnic groups in the Malaysian population is found to be significantly varied when compared to the European population (p < 0.0001). The mean/median reported for the Malay, Chinese, Indian, and European are 2.759/2.869, 3.453/3.290, 2.437/1.970 and 2.001/1.940 respectively.

Conclusion

This study reveals the existence of genetic variation of CCL3L1 in the Malaysian population, and suggests by examining genetic diversity on the ethnicity, and specific geographical region could help in reconstructing human evolutionary history and for the prediction of disease risk related to the CCL3L1 copy number.

Population-wide copy number variation calling using variant call format files from 6,898 individuals

Copy number variants (CNVs) play an important role in a number of human diseases, but the accurate calling of CNVs remains challenging. Most current approaches to CNV detection use raw read alignments, which are computationally intensive to process. We use a regression tree-based approach to call germline CNVs from whole-genome sequencing (WGS, >18x) variant call sets in 6,898 samples across four European cohorts, and describe a rich large variation landscape comprising 1,320 CNVs. Eighty-one percent of detected events have been previously reported in the Database of Genomic Variants. Twenty-three percent of high-quality deletions affect entire genes, and we recapitulate known events such as the GSTM1 and RHD gene deletions. We test for association between the detected deletions and 275 protein levels in 1,457 individuals to assess the potential clinical impact of the detected CNVs. We describe complex CNV patterns underlying an association with levels of the CCL3 protein (MAF = 0.15, p = 3.6x10^-12 ) at the CCL3L3 locus, and a novel cis-association between a low-frequency NOMO1 deletion and NOMO1 protein levels (MAF = 0.02, p = 2.2x10^-7 ). This study demonstrates that existing population-wide WGS call sets can be mined for germline CNVs with minimal computational overhead, delivering insight into a less well-studied, yet potentially impactful class of genetic variant.

Human CCL3L1 copy number variation, gene expression, and the role of the CCL3L1-CCR5 axis in lung function

Background: The CCL3L1-CCR5 signaling axis is important in a number of inflammatory responses, including macrophage function, and T-cell-dependent immune responses. Small molecule CCR5 antagonists exist, including the approved antiretroviral drug maraviroc, and therapeutic monoclonal antibodies are in development. Repositioning of drugs and targets into new disease areas can accelerate the availability of new therapies and substantially reduce costs. As it has been shown that drug targets with genetic evidence supporting their involvement in the disease are more likely to be successful in clinical development, using genetic association studies to identify new target repurposing opportunities could be fruitful. Here we investigate the potential of perturbation of the CCL3L1-CCR5 axis as treatment for respiratory disease. Europeans typically carry between 0 and 5 copies of CCL3L1 and this multi-allelic variation is not detected by widely used genome-wide single nucleotide polymorphism studies.  Methods: We directly measured the complex structural variation of CCL3L1 using the Paralogue Ratio Test and imputed (with validation) CCR5del32 genotypes in 5,000 individuals from UK Biobank, selected from the extremes of the lung function distribution, and analysed DNA and RNAseq data for CCL3L1 from the 1000 Genomes Project. Results: We confirmed the gene dosage effect of CCL3L1 copy number on CCL3L1 mRNA expression levels.  We found no evidence for association of CCL3L1 copy number or CCR5del32 genotype with lung function. Conclusions: These results suggest that repositioning CCR5 antagonists is unlikely to be successful for the treatment of airflow obstruction.

Human CCL3L1 copy number variation, gene expression, and the role of the CCL3L1-CCR5 axis in lung function

Background: The CCL3L1-CCR5 signaling axis is important in a number of inflammatory responses, including macrophage function, and T-cell-dependent immune responses. Small molecule CCR5 antagonists exist, including the approved antiretroviral drug maraviroc, and therapeutic monoclonal antibodies are in development. Repositioning of drugs and targets into new disease areas can accelerate the availability of new therapies and substantially reduce costs. As it has been shown that drug targets with genetic evidence supporting their involvement in the disease are more likely to be successful in clinical development, using genetic association studies to identify new target repurposing opportunities could be fruitful. Here we investigate the potential of perturbation of the CCL3L1-CCR5 axis as treatment for respiratory disease. Europeans typically carry between 0 and 5 copies of CCL3L1 and this multi-allelic variation is not detected by widely used genome-wide single nucleotide polymorphism studies.  Methods: We directly measured the complex structural variation of CCL3L1 using the Paralogue Ratio Test and imputed (with validation) CCR5del32 genotypes in 5,000 individuals from UK Biobank, selected from the extremes of the lung function distribution, and analysed DNA and RNAseq data for CCL3L1 from the 1000 Genomes Project. Results: We confirmed the gene dosage effect of CCL3L1 copy number on CCL3L1 mRNA expression levels.  We found no evidence for association of CCL3L1 copy number or CCR5del32 genotype with lung function. Conclusions: These results suggest that repositioning CCR5 antagonists is unlikely to be successful for the treatment of airflow obstruction.

Population genetics of immune-related multilocus copy number variation in Native Americans

While multiallelic copy number variation (mCNV) loci are a major component of genomic variation, quantifying the individual copy number of a locus and defining genotypes is challenging. Few methods exist to study how mCNV genetic diversity is apportioned within and between populations (i.e. to define the population genetic structure of mCNV). These inferences are critical in populations with a small effective size, such as Amerindians, that may not fit the Hardy-Weinberg model due to inbreeding, assortative mating, population subdivision, natural selection or a combination of these evolutionary factors. We propose a likelihood-based method that simultaneously infers mCNV allele frequencies and the population structure parameter f, which quantifies the departure of homozygosity from the Hardy-Weinberg expectation. This method is implemented in the freely available software CNVice, which also infers individual genotypes using information from both the population and from trios, if available. We studied the population genetics of five immune-related mCNV loci associated with complex diseases (beta-defensins, CCL3L1/CCL4L1, FCGR3A, FCGR3B and FCGR2C) in 12 traditional Native American populations and found that the population structure parameters inferred for these mCNVs are comparable to but lower than those for single nucleotide polymorphisms studied in the same populations.

No Evidence for Association of β-Defensin Genomic Copy Number with HIV Susceptibility, HIV Load during Clinical Latency, or Progression to AIDS

Common single-nucleotide variation in the host accounts for 25% of the variability in the plasma levels of HIV during the clinical latency stage (viral load set point). However, the role of rare variants and copy number variants remains relatively unexplored. Previous work has suggested copy number variation of a cluster of β-defensin genes affects HIV load in treatment-naïve sub-Saharan Africans and rate of response to antiretroviral treatment. Here we analyse a total of 1827 individuals from two cohorts of HIV-infected individuals from Europe and sub-Saharan Africa to investigate the role of β-defensin copy number variation on HIV load at set point. We find no evidence for association of copy number with viral load. We also compare distribution of β-defensin copy number between European cases and controls and find no differences, arguing against a role of β-defensin copy number in HIV acquisition. Taken together, our data argue against an effect of copy number variation of the β-defensin region in the spontaneous control of HIV infection.

Analysis of Copy Number Variation Using the Paralogue Ratio Test (PRT)

Copy number variation (CNV), where a segment of DNA differs in copy number between different individuals, is an extensive and often underappreciated source of genetic variation within species. However, reliably determining copy number of a particular DNA sequence for a large number of samples can be challenging. Here, I describe and review the paralogue ratio test (PRT) in detail. PRT was developed to robustly type the CNV of the beta-defensin locus using small amounts of genomic DNA in a high-throughput manner, and has been applied successfully at many other loci. I discuss the strategies for designing successful PRT assays using both manual and bioinformatics methods, how to optimize experimental conditions, and approaches for analyzing the data. I discuss strengths and weaknesses of the approach, and how to troubleshoot results, as well as the range of problems to which PRT can be a potential solution.

SRBreak: A Read-Depth and Split-Read Framework to Identify Breakpoints of Different Events Inside Simple Copy-Number Variable Regions

Copy-number variation (CNV) has been associated with increased risk of complex diseases. High-throughput sequencing (HTS) technologies facilitate the detection of copy-number variable regions (CNVRs) and their breakpoints. This helps in understanding genome structure as well as their evolution process. Various approaches have been proposed for detecting CNV breakpoints, but currently it is still challenging for tools based on a single analysis method to identify breakpoints of CNVs. It has been shown, however, that pipelines which integrate multiple approaches are able to report more reliable breakpoints. Here, based on HTS data, we have developed a pipeline to identify approximate breakpoints (±10 bp) relating to different ancestral events within a specific CNVR. The pipeline combines read-depth and split-read information to infer breakpoints, using information from multiple samples to allow an imputation approach to be taken. The main steps involve using a normal mixture model to cluster samples into different groups, followed by simple kernel-based approaches to maximize information obtained from read-depth and split-read approaches, after which common breakpoints of groups are inferred. The pipeline uses split-read information directly from CIGAR strings of BAM files, without using a re-alignment step. On simulated data sets, it was able to report breakpoints for very low-coverage samples including those for which only single-end reads were available. When applied to three loci from existing human resequencing data sets (NEGR1, LCE3, IRGM) the pipeline obtained good concordance with results from the 1000 Genomes Project (92, 100, and 82%, respectively). The package is available at https://github.com/hoangtn/SRBreak, and also as a docker-based application at https://registry.hub.docker.com/u/hoangtn/srbreak/.

Copy number variation of genes involved in the hepatitis C virus-human interactome

Copy number variation (CNV) is a newly discovered form of intra-species genetic polymorphism that is defined as deletions or duplications of genome segments ranging from 1 kbp to several Mbp. CNV accounts for the majority of the genetic variation observed in humans (CNV regions cover more than 10% of the human genome); therefore, it may significantly influence both the phenotype and susceptibility to various diseases. Unfortunately, the impact of CNV on a number of diseases, including hepatitis C virus (HCV) infection, remains largely unexplored. Here, we analyzed 421 human genes encoding proteins that have been shown to interact with HCV proteins or genomic RNA (proteins from the HCV-human interactome). We found that 19 of the 421 candidate genes are located in putative CNV regions. For all of these genes, copy numbers were determined for European, Asiatic and African populations using the multiplex ligation-dependent amplification (MLPA) method. As a result, we identified 4 genes, IGLL1, MLLT4, PDPK1, PPP1R13L, for which the CN-genotype ranged from 1 to 6. All of these genes are involved in host-virus interaction; thus, their polymorphism has a potential impact on the development of HCV infection and/or therapy outcome.

Defensin gene variation and HIV/AIDS: a comprehensive perspective needed

Both α- and β-defensins have anti-human immunodeficiency virus activity. These defensins achieve human immunodeficiency virus inhibition through a variety of mechanisms, including direct binding with virions, binding to and modulation of host cell-surface receptors with disruption of intracellular signaling, and functioning as chemokines or cytokines to augment and alter adaptive immune responses. Polymorphisms in the defensin genes have been associated with susceptibility to human immunodeficiency virus infection and disease progression. However, the roles that these defensins and their genetic polymorphisms have in influencing human immunodeficiency virus/acquired immunodeficiency syndrome outcomes are not straightforward and, at times, appear contradictory. Differences in populations, study designs, and techniques for genotyping defensin gene polymorphisms may have contributed to this lack of clarity. In addition, a comprehensive approach, where both subfamilies of defensins and their all-inclusive genetic polymorphism profiles are analyzed, is lacking. Such an approach may reveal whether the human immunodeficiency virus inhibitory activities of α- and β-defensins are based on parallel or divergent mechanisms and may provide further insights into how the genetic predisposition for susceptibility or resistance to human immunodeficiency virus/acquired immunodeficiency syndrome is orchestrated between these molecules.

Complex and multi-allelic copy number variation in human disease

Hundreds of copy number variants are complex and multi-allelic, in that they have many structural alleles and have rearranged multiple times in the ancestors who contributed chromosomes to current humans. Not only are the relationships of these multi-allelic CNVs (mCNVs) to phenotypes generally unknown, but many mCNVs have not yet been described at the basic levels—alleles, allele frequencies, structural features—that support genetic investigation. To date, most reported disease associations to these variants have been ascertained through candidate gene studies. However, only a few associations have reached the level of acceptance defined by durable replications in many cohorts. This likely stems from longstanding challenges in making precise molecular measurements of the alleles individuals have at these loci. However, approaches for mCNV analysis are improving quickly, and some of the unique characteristics of mCNVs may assist future association studies. Their various structural alleles are likely to have different magnitudes of effect, creating a natural allelic series of growing phenotypic impact and giving investigators a set of natural predictions and testable hypotheses about the extent to which each allele of an mCNV predisposes to a phenotype. Also, mCNVs’ low-to-modest correlation to individual single-nucleotide polymorphisms (SNPs) may make it easier to distinguish between mCNVs and nearby SNPs as the drivers of an association signal, and perhaps, make it possible to preliminarily screen candidate loci, or the entire genome, for the many mCNV–disease relationships that remain to be discovered.

An evolutionary history of defensins: a role for copy number variation in maximizing host innate and adaptive immune responses

Defensins represent an evolutionary ancient family of antimicrobial peptides that play diverse roles in human health and disease. Defensins are cationic cysteine-containing multifunctional peptides predominantly expressed by epithelial cells or neutrophils. Defensins play a key role in host innate immune responses to infection and, in addition to their classically described role as antimicrobial peptides, have also been implicated in immune modulation, fertility, development, and wound healing. Aberrant expression of defensins is important in a number of inflammatory diseases as well as modulating host immune responses to bacteria, unicellular pathogens, and viruses. In parallel with their role in immunity, in other species, defensins have evolved alternative functions, including the control of coat color in dogs. Defensin genes reside in complex genomic regions that are prone to structural variations and some defensin family members exhibit copy number variation (CNV). Structural variations have mediated, and continue to influence, the diversification and expression of defensin family members. This review highlights the work currently being done to better understand the genomic architecture of the β-defensin locus. It evaluates current evidence linking defensin CNV to autoimmune disease (i.e., Crohn’s disease and psoriasis) as well as the contribution CNV has in influencing immune responses to HIV infection.

Genetics of Sub-Saharan African Human Population: Implications for HIV/AIDS, Tuberculosis, and Malaria

Sub-Saharan Africa has continued leading in prevalence and incidence of major infectious disease killers such as HIV/AIDS, tuberculosis, and malaria. Epidemiological triad of infectious diseases includes susceptible host, pathogen, and environment. It is imperative that all aspects of vertices of the infectious disease triad are analysed to better understand why this is so. Studies done to address this intriguing reality though have mainly addressed pathogen and environmental components of the triad. Africa is the most genetically diverse region of the world as well as being the origin of modern humans. Malaria is relatively an ancient infection in this region as compared to TB and HIV/AIDS; from the evolutionary perspective, we would draw lessons that this ancestrally unique population now under three important infectious diseases both ancient and exotic will be skewed into increased genetic diversity; moreover, other evolutionary forces are also still at play. Host genetic diversity resulting from many years of malaria infection has been well documented in this population; we are yet to account for genetic diversity from the trio of these infections. Effect of host genetics on treatment outcome has been documented. Host genetics of sub-Saharan African population and its implication to infectious diseases are an important aspect that this review seeks to address.

Inferring mechanisms of copy number change from haplotype structures at the human DEFA1A3 locus

BACKGROUND

The determination of structural haplotypes at copy number variable regions can indicate the mechanisms responsible for changes in copy number, as well as explain the relationship between gene copy number and expression. However, obtaining spatial information at regions displaying extensive copy number variation, such as the DEFA1A3 locus, is complex, because of the difficulty in the phasing and assembly of these regions. The DEFA1A3 locus is intriguing in that it falls within a region of high linkage disequilibrium, despite its high variability in copy number (n = 3-16); hence, the mechanisms responsible for changes in copy number at this locus are unclear.

RESULTS

In this study, a region flanking the DEFA1A3 locus was sequenced across 120 independent haplotypes with European ancestry, identifying five common classes of DEFA1A3 haplotype. Assigning DEFA1A3 class to haplotypes within the 1000 Genomes project highlights a significant difference in DEFA1A3 class frequencies between populations with different ancestry. The features of each DEFA1A3 class, for example, the associated DEFA1A3 copy numbers, were initially assessed in a European cohort (n = 599) and replicated in the 1000 Genomes samples, showing within-class similarity, but between-class and between-population differences in the features of the DEFA1A3 locus. Emulsion haplotype fusion-PCR was used to generate 61 structural haplotypes at the DEFA1A3 locus, showing a high within-class similarity in structure.

CONCLUSIONS

Structural haplotypes across the DEFA1A3 locus indicate that intra-allelic rearrangement is the predominant mechanism responsible for changes in DEFA1A3 copy number, explaining the conservation of linkage disequilibrium across the locus. The identification of common structural haplotypes at the DEFA1A3 locus could aid studies into how DEFA1A3 copy number influences expression, which is currently unclear.

Identification of sequence variants in the CCL3 chemokine gene family in the HapMap West African reference population

Gene copy number variation (CNV) of the CC Chemokine ligand 3-Like-1 (CCL3L1) gene located on chromosome 17q12 has been associated with many diseases, including viral infections and autoimmune diseases. High sequence homology between CCL3L1 and three other related genes within the same cluster, CCL3, CCL3L2, and CCL3L3, make it difficult to determine the copy number of each gene as well as distinguishing variants within each gene versus between genes. We identified a total of 50SNPs, 31 known and 19 novel SNPs, in a subset of West Africa Reference (Yoruba individuals from Ibadan, Nigeria (YRI)) samples from HapMap. One of these previously unidentified variations is a non-synonymous change while several other unreported variations are located near potential regulatory sites. The variations identified in these immune-related genes from this study will shed light in the understanding of both structural and nucleotide polymorphisms that can be used in association studies of diseases in populations.

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.

Comparison of a quantitative Real-Time PCR assay and droplet digital PCR for copy number analysis of the CCL4L genes

The controversy surrounding the findings that copy number variation, of the CCL3 encoding genes, influences HIV-1 infection and disease progression has been in part attributed to the variable results obtained from methods used for copy number evaluation. Like CCL3, the genes encoding the CC chemokine CCL4, also a natural ligand of the CCR5 receptor, are found to occur in population-specific multiple copy number and have been shown to play a protective role against HIV-1. This study evaluated the standard method of quantitative Real-Time PCR (qPCR) and droplet digital PCR (ddPCR) for CCL4L gene copy number determination. The CCL4 encoding genes are CCL4, occurring in two copies per diploid genome (pdg), and the non-allelic CCL4L genes, comprised of CCL4L1 and CCL4L2, which are both found in multiple copies pdg. Copy number of CCL4L, CCL4L1 and CCL4L2 was determined in a cohort of HIV-1-uninfected individuals from the South African Black (n=23) and Caucasian (n=32) population groups using qPCR and ddPCR. A stronger correlation between the number of CCL4L copies and the sum of CCL4L1 and CCL4L2 copies generated by ddPCR (r=0.99, p<0.0001) compared to qPCR (r=0.87, p<0.0001) was observed. Real-Time qPCR exhibited greater inaccuracy at higher copy numbers which is particularly relevant to our cohort of Black individuals who have a higher range of CCL4L copies (3-6) compared to Caucasians (0-4) and a higher population median (4 and 2, respectively). Medians and ranges of CCL4L1 (Black: 2, 0-4, Caucasian: 0, 0-2) and CCL4L2 (Black: 2, 1-5, Caucasian: 2, 0-3) were also higher in the Black population. Droplet digital PCR was shown to be a far superior method to qPCR for assessment of CCL4 gene copy number variation, the accuracy of which is essential for studies of the contribution of variable gene copy number to phenotypic outcomes of host infection and disease course.