Technical considerations for genotyping multi-allelic copy number variation (CNV), in regions of segmental duplication

Association of Structural Forms of 17q21.31 with the Risk of Progressive Supranuclear Palsy and MAPT Sub-haplotypes

Importance

The chromosome 17q21.31 region, containing a 900 Kb inversion that defines H1 and H2 haplotypes, represents the strongest genetic risk locus in progressive supranuclear palsy (PSP). In addition to H1 and H2, various structural forms of 17q21.31, characterized by the copy number of α, β, and γ duplications, have been identified. However, the specific effect of each structural form on the risk of PSP has never been evaluated in a large cohort study.

Objective

To assess the association of different structural forms of 17q.21.31, defined by the copy numbers of α, β, and γ duplications, with the risk of PSP and MAPT sub-haplotypes.

Design setting and participants

Utilizing whole genome sequencing data of 1,684 (1,386 autopsy confirmed) individuals with PSP and 2,392 control subjects, a case-control study was conducted to investigate the association of copy numbers of α, β, and γ duplications and structural forms of 17q21.31 with the risk of PSP. All study subjects were selected from the Alzheimer's Disease Sequencing Project (ADSP) Umbrella NG00067.v7. Data were analyzed between March 2022 and November 2023.

Main outcomes and measures

The main outcomes were the risk (odds ratios [ORs]) for PSP with 95% CIs. Risks for PSP were evaluated by logistic regression models.

Results

The copy numbers of α and β were associated with the risk of PSP only due to their correlation with H1 and H2, while the copy number of γ was independently associated with the increased risk of PSP. Each additional duplication of γ was associated with 1.10 (95% CI, 1.04-1.17; P = 0.0018) fold of increased risk of PSP when conditioning H1 and H2. For the H1 haplotype, addition γ duplications displayed a higher odds ratio for PSP: the odds ratio increases from 1.21 (95%CI 1.10-1.33, P = 5.47 × 10-5) for H1β1γ1 to 1.29 (95%CI 1.16-1.43, P = 1.35 × 10-6) for H1β1γ2, 1.45 (95%CI 1.27-1.65, P = 3.94 × 10-8) for H1β1γ3, and 1.57 (95%CI 1.10-2.26, P = 1.35 × 10-2) for H1β1γ4. Moreover, H1β1γ3 is in linkage disequilibrium with H1c (R2 = 0.31), a widely recognized MAPT sub-haplotype associated with increased risk of PSP. The proportion of MAPT sub-haplotypes associated with increased risk of PSP (i.e., H1c, H1d, H1g, H1o, and H1h) increased from 34% in H1β1γ1 to 77% in H1β1γ4.

Conclusions and relevance

This study revealed that the copy number of γ was associated with the risk of PSP independently from H1 and H2. The H1 haplotype with more γ duplications showed a higher odds ratio for PSP and were associated with MAPT sub-haplotypes with increased risk of PSP. These findings expand our understanding of how the complex structure at 17q21.31 affect the risk of PSP.

Quantitative PCR from human genomic DNA: The determination of gene copy numbers for congenital adrenal hyperplasia and RCCX copy number variation

Quantitative PCR (qPCR) is used for the determination of gene copy number (GCN). GCNs contribute to human disorders, and characterize copy number variation (CNV). The single laboratory method validations of duplex qPCR assays with hydrolysis probes on CYP21A1P and CYP21A2 genes, residing a CNV (RCCX CNV) and related to congenital adrenal hyperplasia, were performed using 46 human genomic DNA samples. We also performed the verifications on 5 qPCR assays for the genetic elements of RCCX CNV; C4A, C4B, CNV breakpoint, HERV-K(C4) CNV deletion and insertion alleles. Precision of each qPCR assay was under 1.01 CV%. Accuracy (relative error) ranged from 4.96±4.08% to 9.91±8.93%. Accuracy was not tightly linked to precision, but was significantly correlated with the efficiency of normalization using the RPPH1 internal reference gene (Spearman's ρ: 0.793-0.940, p>0.0001), ambiguity (ρ = 0.671, p = 0.029) and misclassification (ρ = 0.769, p = 0.009). A strong genomic matrix effect was observed, and target-singleplex (one target gene in one assay) qPCR was able to appropriately differentiate 2 GCN from 3 GCN at best. The analysis of all GCNs from the 7 qPCR assays using a multiplex approach increased the resolution of differentiation, and produced 98% of GCNs unambiguously, and all of which were in 100% concordance with GCNs measured by Southern blot, MLPA and aCGH. We conclude that the use of an internal (in one assay with the target gene) reference gene, the use of allele-specific primers or probes, and the multiplex approach (in one assay or different assays) are crucial for GCN determination using qPCR or other methods.

Genome structural variation in human evolution

Structural variation (SV) is a large difference (typically >100 bp) in the genomic structure of two genomes and includes both copy number variation and variation that does not change copy number of a genomic region, such as an inversion. Improved reference genomes, combined with widespread genome sequencing using short-read sequencing technology, and increasingly using long-read sequencing, have reignited interest in SV. Recent large-scale studies and functional focused analyses have highlighted the role of SV in human evolution. In this review, we highlight human-specific SVs involved in changes in the brain, population-specific SVs that affect response to the environment, including adaptation to diet and infectious diseases, and summarise the contribution of archaic hominin admixture to present-day human SV.

Genotyping complex structural variation at the malaria-associated human glycophorin locus using a PCR-based strategy

Structural variation in the human genome can affect risk of disease. An example is a complex structural variant of the human glycophorin gene cluster, called DUP4, which is associated with a clinically significant level of protection against severe malaria. The human glycophorin gene cluster harbours at least 23 distinct structural variants, and accurate genotyping of this complex structural variation remains a challenge. Here, we use a polymerase chain reaction-based strategy to genotype structural variation at the human glycophorin gene cluster, including the alleles responsible for the U- blood group. We validate our approach, based on a triplex paralogue ratio test, on publically available samples from the 1000 Genomes project. We then genotype 574 individuals from a longitudinal birth cohort (Tori-Bossito cohort) using small amounts of DNA at low cost. Our approach readily identifies known deletions and duplications, and can potentially identify novel variants for further analysis. It will allow exploration of genetic variation at the glycophorin locus, and investigation of its relationship with malaria, in large sample sets at minimal cost, using standard molecular biology equipment.

Structural variation of the malaria-associated human glycophorin A-B-E region

BACKGROUND

Approximately 5% of the human genome shows common structural variation, which is enriched for genes involved in the immune response and cell-cell interactions. A well-established region of extensive structural variation is the glycophorin gene cluster, comprising three tandemly-repeated regions about 120 kb in length and carrying the highly homologous genes GYPA, GYPB and GYPE. Glycophorin A (encoded by GYPA) and glycophorin B (encoded by GYPB) are glycoproteins present at high levels on the surface of erythrocytes, and they have been suggested to act as decoy receptors for viral pathogens. They are receptors for the invasion of the protist parasite Plasmodium falciparum, a causative agent of malaria. A particular complex structural variant, called DUP4, creates a GYPB-GYPA fusion gene known to confer resistance to malaria. Many other structural variants exist across the glycophorin gene cluster, and they remain poorly characterised.

RESULTS

Here, we analyse sequences from 3234 diploid genomes from across the world for structural variation at the glycophorin locus, confirming 15 variants in the 1000 Genomes project cohort, discovering 9 new variants, and characterising a selection of these variants using fibre-FISH and breakpoint mapping at the sequence level. We identify variants predicted to create novel fusion genes and a common inversion duplication variant at appreciable frequencies in West Africans. We show that almost all variants can be explained by non-allelic homologous recombination and by comparing the structural variant breakpoints with recombination hotspot maps, confirm the importance of a particular meiotic recombination hotspot on structural variant formation in this region.

CONCLUSIONS

We identify and validate large structural variants in the human glycophorin A-B-E gene cluster which may be associated with different clinical aspects of malaria.

Sudden Cardiac Death and Copy Number Variants: What Do We Know after 10 Years of Genetic Analysis?

Over the last ten years, analysis of copy number variants has increasingly been applied to the study of arrhythmogenic pathologies associated with sudden death, mainly due to significant advances in the field of massive genetic sequencing. Nevertheless, few published reports have focused on the prevalence of copy number variants associated with sudden cardiac death. As a result, the frequency of these genetic alterations in arrhythmogenic diseases as well as their genetic interpretation and clinical translation has not been established. This review summarizes the current available data concerning copy number variants in sudden cardiac death-related diseases.

Spinal Muscular Atrophy in the Black South African Population: A Matter of Rearrangement?

Spinal muscular atrophy (SMA) is a neuromuscular disorder, characterized by muscle atrophy and impaired mobility. A homozygous deletion of survival motor neuron 1 (SMN1), exon 7 is the main cause of SMA in ~94% of patients worldwide, but only accounts for 51% of South African (SA) black patients. SMN1 and its highly homologous centromeric copy, survival motor neuron 2 (SMN2), are located in a complex duplicated region. Unusual copy number variations (CNVs) have been reported in black patients, suggesting the presence of complex pathogenic rearrangements. The aim of this study was to further investigate the genetic cause of SMA in the black SA population. Multiplex ligation-dependent probe amplification (MLPA) testing was performed on 197 unrelated black patients referred for SMA testing (75 with a homozygous deletion of SMN1, exon 7; 50 with a homozygous deletion of SMN2, exon 7; and 72 clinically suggestive patients with no homozygous deletions). Furthermore, 122 black negative controls were tested. For comparison, 68 white individuals (30 with a homozygous deletion of SMN1, exon 7; 8 with a homozygous deletion of SMN2, exon 7 and 30 negative controls) were tested. Multiple copies (>2) of SMN1, exon 7 were observed in 50.8% (62/122) of black negative controls which could mask heterozygous SMN1 deletions and potential pathogenic CNVs. MLPA is not a reliable technique for detecting carriers in the black SA population. Large deletions extending into the rest of SMN1 and neighboring genes were more frequently observed in black patients with homozygous SMN1, exon 7 deletions when compared to white patients. Homozygous SMN2, exon 7 deletions were commonly observed in black individuals. No clear pathogenic CNVs were identified in black patients but discordant copy numbers of exons suggest complex rearrangements, which may potentially interrupt the SMN1 gene. Only 8.3% (6/72) of clinically suggestive patients had heterozygous deletions of SMN1, exon 7 (1:0) which is lower than previous SA reports of 69.5%. This study emphasizes the lack of understanding of the architecture of the SMN region as well as the cause of SMA in the black SA population. These factors need to be taken into account when counseling and performing diagnostic testing in black populations.

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.

Divergent patterns of genic copy number variation in KCNIP1 gene reveal risk locus of type 2 diabetes in Chinese population

Copy number variation (CNV) has emerged as another important genetic marker in addition to SNP for understanding etiology of complex disease. Kv channel interacting protein 1 (KCNIP1) is a Ca²⁺-dependent transcriptional modulator that contributes to the regulation of insulin secretion. Previous genome-wide CNV assay identified the KCNIP1 gene encompassing a CNV region, however, its further effect and risk rate on type 2 diabetes (T2D) have rarely been addressed, especially in Chinese population. The current study aims to detect and excavate genetic distribution profile of KCNIP1 CNV in Chinese T2D and control populations, and further to investigate the associations with clinical characteristics. Divergent patterns of the KCNIP1 CNV were identified (p < 0.01), in which the copy number gain was predominant in T2D, while the copy number normal accounted for the most in control group. Consistently, the individuals with copy number gain showed significant risk on T2D (OR = 4.550, p < 0.01). The KCNIP1 copy numbers presented significantly positive correlations with fasting plasma glucose and glycated hemoglobin in T2D. For OGTT test, the T2D patients with copy number gain had remarkably elevated glucose contents (60, 120, 180-min, p < 0.05 or p < 0.01) and diminished insulin levels (60, 120-min, p < 0.05) than those with copy number loss and normal, which suggested that the KCNIP1 CNV was correlated with the glucose and insulin action. This is the first CNV association study of the KCNIP1 gene in Chinese population, and these data indicated that KCNIP1 might function as a T2D-susceptibility gene whose dysregulation alters insulin production.

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.

Complement receptor 1 gene (CR1) intragenic duplication and risk of Alzheimer's disease

Single nucleotide variants (SNVs) within and surrounding the complement receptor 1 (CR1) gene show some of the strongest genome-wide association signals with late-onset Alzheimer’s disease. Some studies have suggested that this association signal is due to a duplication allele (CR1-B) of a low copy repeat (LCR) within the CR1 gene, which increases the number of complement C3b/C4b-binding sites in the mature receptor. In this study, we develop a triplex paralogue ratio test assay for CR1 LCR copy number allowing large numbers of samples to be typed with a limited amount of DNA. We also develop a CR1-B allele-specific PCR based on the junction generated by an historical non-allelic homologous recombination event between CR1 LCRs. We use these methods to genotype CR1 and measure CR1-B allele frequency in both late-onset and early-onset cases and unaffected controls from the United Kingdom. Our data support an association of late-onset Alzheimer’s disease with the CR1-B allele, and confirm that this allele occurs most frequently on the risk haplotype defined by SNV alleles. Furthermore, regression models incorporating CR1-B genotype provide a better fit to our data compared to incorporating the SNV-defined risk haplotype, supporting the CR1-B allele as the variant underlying the increased risk of late-onset Alzheimer’s disease.

Droplet digital PCR is a powerful technique to demonstrate frequent FGFR1 duplication in dysembryoplastic neuroepithelial tumors

Dysembryoplastic neuroepithelial tumors (DNT) share V600E mutation in the BRAF gene with other low grade neuroepithelial tumors (LGNTs). FGFR1 internal tandem duplication of the tyrosine-kinase domain (FGFR1-ITD), another genetic alteration that also leads to MAP kinase pathway alteration, has been previously reported in LGNTs by whole-genome sequencing. In the present study we searched for FGFR1-ITD by droplet digital PCR (DDPCR™) and for FGFR1 point mutations by HRM-sequencing in a series of formalin-fixed paraffin-embedded (FFPE) LGNTs including 12 DNT, 2 oligodendrogliomas lacking IDH mutation and 1p/19q co- deletion (pediatric-type oligodendrogliomas; PTOs), 3 pediatric diffuse astrocytomas (PDAs), 14 gangliogliomas (GGs) and 5 pilocytic astrocytomas (PAs). We showed by DDPCR™ that 5/12 DNT, but none of the other LGNTs, demonstrated FGFR1-ITD. In addition, these cases also accumulated phosphorylated-FGFR1 protein as shown by immunohistochemistry. FGFR1^G539R point mutation was only recorded in one DNT that also showed FGFR1-ITD. Interestingly, these FGFR1 alterations were mutually exclusive from BRAF^V600E mutation that was recorded in 13 LGNTs (3 DNTs, 1 PTO, 2 PDAs, 5 GGs and 2 PAs). Therefore, FGFR1 alteration mainly represented by FGFR1-ITD is a frequent event in DNT. DDPCR™ is an easy and alternative method than whole-genome sequencing to detect FGFR1-ITD in FFPE brain tumors, in routine practice.

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.

Copy Number Variation Disorders

Purpose of Review

Copy number variation (CNV) disorders arise from the dosage imbalance of one or more gene(s), resulting from deletions, duplications or other genomic rearrangements that lead to the loss or gain of genetic material. Several disorders, characterized by multiple birth defects and neurodevelopmental abnormalities, have been associated with relatively large (>1 Mb) and often recurrent CNVs. CNVs have also been implicated in the etiology of neuropsychiatric disorders including autism and schizophrenia as well as other common complex diseases. Thus, CNVs have a significant impact on human health and disease.

Recent Findings

The use of increasingly higher resolution, genomewide analysis has greatly enhanced the detection of genetic variation, including CNVs. Furthermore, the availability of comprehensive genetic variation data from large cohorts of healthy controls has the potential to greatly improve the identification of disease associated genetic variants in patient samples.

Summary

This review discusses the current knowledge about CNV disorders, including the mechanisms underlying their formation and phenotypic outcomes, and the advantages and limitations of current methods of detection and disease association.

CBF2A-CBF4B genomic region copy numbers alongside the circadian clock play key regulatory mechanisms driving expression of FR-H2 CBFs

The C-Repeat Binding Factors (CBFs) are DNA-binding transcriptional activators that were identified using Arabidopsis thaliana. In barley, Hordeum vulgare, a cluster of CBF genes reside at FROST RESISTANCE-H2, one of two loci having major effects on winter-hardiness. FR-H2 was revealed in a population derived from the winter barley 'Nure' and the spring barley 'Trèmois'. 'Nure' harbors two to three copies of CBF2A and CBF4B as a consequence of tandem iteration of the genomic region encompassing these genes whereas 'Trèmois' harbors single copies, and these copy number differences are associated with their transcript level differences. Here we explore further the relationship between FR-H2 CBF gene copy number and transcript levels using 'Admire', a winter barley accumulating FR-H2 CBF gene transcripts to very high levels, and a group of lines related to 'Admire' through descent. DNA blot hybridization indicated the CBF2A-CBF4B genomic region is present in 7-8 copies in 'Admire' and is highly variable in copy number across the lines related to 'Admire'. At normal growth temperatures transcript levels of CBF12, CBF14, and CBF16 were higher in lines having greater CBF2A-CBF4B genomic region copy numbers than in lines having fewer copy numbers at peak expression level time points controlled by the circadian clock. Chromatin immunoprecipitation indicated CBF2 was at the CBF12 and CBF16 promoters at normal growth temperatures. These data support a scenario in which CBF2A-CBF4B genomic region copy numbers affect expression of other FR-H2 CBFs through a mechansim in which these other FR-H2 CBFs are activated by those in the copy number variable unit.

MLPA-Based Analysis of Copy Number Variation in Plant Populations

Copy number variants (CNVs) are intraspecies duplications/deletions of large DNA segments (>1 kb). A growing number of reports highlight the functional and evolutionary impact of CNV in plants, increasing the need for appropriate tools that enable locus-specific CNV genotyping on a population scale. Multiplex ligation-dependent probe amplification (MLPA) is considered a gold standard in genotyping CNV in humans. Consequently, numerous commercial MLPA assays for CNV-related human diseases have been created. We routinely genotype complex multiallelic CNVs in human and plant genomes using the modified MLPA procedure based on fully synthesized oligonucleotide probes (90-200 nt), which greatly simplifies the design process and allows for the development of custom assays. Here, we present a step-by-step protocol for gene-specific MLPA probe design, multiplexed assay setup and data analysis in a copy number genotyping experiment in plants. As a case study, we present the results of a custom assay designed to genotype the copy number status of 12 protein coding genes in a population of 80 Arabidopsis accessions. The genes were pre-selected based on whole genome sequencing data and are localized in the genomic regions that display different levels of population-scale variation (non-variable, biallelic, or multiallelic, as well as CNVs overlapping whole genes or their fragments). The presented approach is suitable for population-scale validation of the CNV regions inferred from whole genome sequencing data analysis and for focused analysis of selected genes of interest. It can also be very easily adopted for any plant species, following optimization of the template amount and design of the appropriate control probes, according to the general guidelines presented in this paper.

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.

Dosage of Sex Chromosomal Genes in Blood Deposited on Filter Paper for Neonatal Screening of Sex Chromosome Aneuploidy

AIMS

In this study, we examined the doses of the stature homeobox (SHOX), vesicle-associated membrane protein 7 (VAMP7), and SRY genes to establish a protocol for using peripheral blood samples deposited on filter paper for the screening of sex chromosome aneuploidy in neonates. We also measured correlations with karyotypes to assess this method as a neonatal screening strategy.

MATERIALS AND METHODS

This was an observational, descriptive, comparative blind study. Thirty-two healthy young adults (17 women, 15 men; age, ≥18 years), four patients with known sex chromosome aneuploidy (positive control group), and 1000 healthy newborns were included. Gene dosages were determined using quantitative real-time polymerase chain reaction (RT-PCR). Values with standard deviations (SDs) of three or more were considered abnormal.

RESULTS

Men and women differed in the gene dosage of the SRY gene. Cases with Turner syndrome showed values below 3 SDs for SHOX and VAMP7 genes, and cases with Klinefelter syndrome showed values above 3 SDs for SHOX and VAMP7 genes. Two suspected cases of sex chromosome aneuploidy were diagnosed using our neonatal screening strategy; these cases were confirmed as Turner syndrome and 47,XYY syndrome by karyotyping.

CONCLUSIONS

Our data establish a basis for the determination of chromosomal sex and neonatal screening of sex chromosome aneuploidy using RT-PCR.

Arabidopsis thaliana population analysis reveals high plasticity of the genomic region spanning MSH2, AT3G18530 and AT3G18535 genes and provides evidence for NAHR-driven recurrent CNV events occurring in this location

Background

Intraspecies copy number variations (CNVs), defined as unbalanced structural variations of specific genomic loci, ≥1 kb in size, are present in the genomes of animals and plants. A growing number of examples indicate that CNVs may have functional significance and contribute to phenotypic diversity. In the model plant Arabidopsis thaliana at least several hundred protein-coding genes might display CNV; however, locus-specific genotyping studies in this plant have not been conducted.

Results

We analyzed the natural CNVs in the region overlapping MSH2 gene that encodes the DNA mismatch repair protein, and AT3G18530 and AT3G18535 genes that encode poorly characterized proteins. By applying multiplex ligation-dependent probe amplification and droplet digital PCR we genotyped those genes in 189 A. thaliana accessions. We found that AT3G18530 and AT3G18535 were duplicated (2–14 times) in 20 and deleted in 101 accessions. MSH2 was duplicated in 12 accessions (up to 12-14 copies) but never deleted. In all but one case, the MSH2 duplications were associated with those of AT3G18530 and AT3G18535. Considering the structure of the CNVs, we distinguished 5 genotypes for this region, determined their frequency and geographical distribution. We defined the CNV breakpoints in 35 accessions with AT3G18530 and AT3G18535 deletions and tandem duplications and showed that they were reciprocal events, resulting from non-allelic homologous recombination between 99 %-identical sequences flanking these genes. The widespread geographical distribution of the deletions supported by the SNP and linkage disequilibrium analyses of the genomic sequence confirmed the recurrent nature of this CNV.

Conclusions

We characterized in detail for the first time the complex multiallelic CNV in Arabidopsis genome. The region encoding MSH2, AT3G18530 and AT3G18535 genes shows enormous variation of copy numbers among natural ecotypes, being a remarkable example of high Arabidopsis genome plasticity. We provided the molecular insight into the mechanism underlying the recurrent nature of AT3G18530-AT3G18535 duplications/deletions. We also performed the first direct comparison of the two leading experimental methods, suitable for assessing the DNA copy number status. Our comprehensive case study provides foundation information for further analyses of CNV evolution in Arabidopsis and other plants, and their possible use in plant breeding.

Electronic supplementary material

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

A Genome-Wide Association Study to Identify Potential Germline Copy Number Variants for Sporadic Breast Cancer Susceptibility

Breast cancer (BC) predisposition in populations arises from both genetic and nongenetic risk factors. Structural variations such as copy number variations (CNVs) are heritable determinants for disease susceptibility. The primary objectives of this study are (1) to identify CNVs associated with sporadic BC using a genome-wide association study (GWAS) design; (2) to utilize 2 distinct CNV calling algorithms to identify concordant CNVs as a strategy to reduce false positive associations in the hypothesis-generating GWAS discovery phase, and (3) to identify potential candidate CNVs for follow-up replication studies. We used Affymetrix SNP Array 6.0 data profiled on Caucasian subjects (422 cases/348 controls) to call CNVs using algorithms implemented in Nexus Copy Number and Partek Genomics Suite software. Nexus algorithm identified CNVs associated with BC (731 autosomal CNVs with >5% frequency in the total sample and Q < 0.05). Thirteen CNVs were identified when Partek algorithm-called CNVs were overlapped with Nexus-identified CNVs; these CNVs showed concordances for frequency, effect size, and direction. Coding genes present within BC-associated CNVs were known to play a role in disease etiology and prognosis. Long noncoding RNAs identified within CNVs showed tissue-specific expression, indicating potential functional relevance of the findings. The identified candidate CNVs warrant independent replication.

Analysıs of gene copy number changes ın head and neck cancer

OBJECTIVE

Chromosomal alterations and copy number changes are frequent events in tumors, leading to amplification of focal regions containing several oncogenes. Gains and losses of several regions have been reported in head and neck cancer (HNC) but the copy number changes of the individual genes located in these regions have not been analyzed so far. In this study we aimed to analyze the copy number variations in patients with HNC.

DESIGN

Prospective study SETTING: University hospital PARTICIPANTS: 50 patients with squamous cell carcinoma of the head and neck METHODS: Copy number changes and amplifications of 22 genes in tumors and matched tissue were analyzed by MLPA which allows simultaneous analysis of gene copy numbers in multiple genetic regions.

RESULTS

Amplifications were observed in 52% and losses were detected in 20% of the samples. Chromosome 8 was found to harbor the most frequent copy number alterations. The most frequently amplified genes were CCND1 and the MED1 genes followed by the MTDH and MYC genes on the long arm and ZNF703 on the short arm of chromosome 8. Amplification of the ZNF703, PRDM14 and MYC genes were highly correlated suggesting that the genes displaying high copy number changes on chromosome 8 collaborate during carcinogenesis.

CONCLUSIONS

The alterations found in our study supports the contribution of gene amplifications and indicate cooperation between certain oncogenes in the pathogenesis of HNSCC. Correlations between amplification of less familiar genes and known oncogenes warrant further investigation. This article is protected by copyright. All rights reserved.

Multiallelic copy number variation in the complement component 4A (C4A) gene is associated with late-stage age-related macular degeneration (AMD)

BACKGROUND

Age-related macular degeneration (AMD) is the leading cause of vision loss in Western societies with a strong genetic component. Candidate gene studies as well as genome-wide association studies strongly implicated genetic variations in complement genes to be involved in disease risk. So far, no association of AMD with complement component 4 (C4) was reported probably due to the complex nature of the C4 locus on chromosome 6.

METHODS

We used multiplex ligation-dependent probe amplification (MLPA) to determine the copy number of the C4 gene as well as of both relevant isoforms, C4A and C4B, and assessed their association with AMD using logistic regression models.

RESULTS

Here, we report on the analysis of 2645 individuals (1536 probands and 1109 unaffected controls), across three different centers, for multiallelic copy number variation (CNV) at the C4 locus. We find strong statistical significance for association of increased copy number of C4A (OR 0.81 (0.73; 0.89);P = 4.4 × 10(-5)), with the effect most pronounced in individuals over 78 years (OR 0.67 (0.55; 0.81)) and females (OR 0.77 (0.68; 0.87)). Furthermore, this association is independent of known AMD-associated risk variants in the nearby CFB/C2 locus, particularly in females and in individuals over 78 years.

CONCLUSIONS

Our data strengthen the notion that complement dysregulation plays a crucial role in AMD etiology, an important finding for early intervention strategies and future therapeutics. In addition, for the first time, we provide evidence that multiallelic CNVs are associated with AMD pathology.

Characterization of copy number variants for CCL3L1 gene in rheumatoid arthritis for French trio families and Tunisian cases and controls

Analyses of copy number variants (CNVs) for candidate genes in complex diseases are currently a promising research field. CNVs of C-C chemokine ligand 3-like 1 (CCL3L1) gene are candidate genomic factors in rheumatoid arthritis (RA). We investigated CCL3L1 CNVs association with a case-control study in Tunisians and a transmission analysis in French trio families. Relative copy number (rCN) of CCL3L1 gene was quantified by droplet digital PCR (ddPCR) in 100 French trio families (RA patients and their two parents) and in 166 RA cases and 102 healthy controls from Tunisia. We calculated odds ratio (OR) to investigate association risk for CCL3L1 CNVs in RA. rCN identified varied from 0 to 4 in the French population and from 0 to 7 in the Tunisian population. A significant difference was observed in the distribution of these rCNs between the two populations (p = 2.34 × 10(-10)), as when rCN from French and Tunisian RA patients were compared (p = 2.83 × 10(-5)). CNVs transmission in French RA trios allowed the characterization of genotypes with the presence of tandem duplication and triplication on the same chromosome. RA association tests highlighted a protective effect of rCN = 5 for CCL3L1 gene in the Tunisian population (OR = 0.056; CI 95 % [0.01-0.46]). Characterization of CCL3L1 CNVs with ddPCR methodology highlighted specific CN genotypes in a French family sample. A copy number polymorphism of a RA candidate gene was quantified, and its significant association with RA was revealed in a Tunisian sample.

Determining multiallelic complex copy number and sequence variation from high coverage exome sequencing data

BACKGROUND

Copy number variation (CNV) is a major component of genomic variation, yet methods to accurately type genomic CNV lag behind methods that type single nucleotide variation. High-throughput sequencing can contribute to these methods by using sequence read depth, which takes the number of reads that map to a given part of the reference genome as a proxy for copy number of that region, and compares across samples. Furthermore, high-throughput sequencing also provides information on the sequence differences between copies within and between individuals.

METHODS

In this study we use high-coverage phase 3 exome sequences of the 1000 Genomes project to infer diploid copy number of the beta-defensin genomic region, a well-studied CNV that carries several beta-defensin genes involved in the antimicrobial response, signalling, and fertility. We also use these data to call sequence variants, a particular challenge given the multicopy nature of the region.

RESULTS

We confidently call copy number and sequence variation of the beta-defensin genes on 1285 samples from 26 global populations, validate copy number using Nanostring nCounter and triplex paralogue ratio test data. We use the copy number calls to verify the genomic extent of the CNV and validate sequence calls using analysis of cloned PCR products. We identify novel variation, mostly individually rare, predicted to alter amino-acid sequence in the beta-defensin genes. Such novel variants may alter antimicrobial properties or have off-target receptor interactions, and may contribute to individuality in immunological response and fertility.

CONCLUSIONS

Given that 81% of identified sequence variants were not previously in dbSNP, we show that sequence variation in multiallelic CNVs represent an unappreciated source of genomic diversity.

Counting copy number and calories

Copy number variation (CNV) at several genomic loci has been associated with different human traits and diseases, but in many cases the findings could not be replicated. A new study provides insights into the degree of variation present at the amylase locus and calls into question a previous association between amylase copy number and body mass index.

The impact of human copy number variation on gene expression

Recent years have witnessed a flurry of important technological and methodological developments in the discovery and analysis of copy number variations (CNVs), which are increasingly enabling the systematic evaluation of their impact on a broad range of phenotypes from molecular-level (intermediate) traits to higher-order clinical phenotypes. Like single nucleotide variants in the human genome, CNVs have been linked to complex traits in humans, including disease and drug response. These recent developments underscore the importance of incorporating complex forms of genetic variation into disease mapping studies and promise to transform our understanding of genome function and the genetic basis of disease. Here we review some of the findings that have emerged from transcriptome studies of CNVs facilitated by the rapid advances in -omics technologies and corresponding methodologies.

Evolution of the rapidly mutating human salivary agglutinin gene (DMBT1) and population subsistence strategy

The dietary change resulting from the domestication of plant and animal species and development of agriculture at different locations across the world was one of the most significant changes in human evolution. An increase in dietary carbohydrates caused an increase in dental caries following the development of agriculture, mediated by the cariogenic oral bacterium Streptococcus mutans. Salivary agglutinin [SAG, encoded by the deleted in malignant brain tumors 1 (DMBT1) gene] is an innate immune receptor glycoprotein that binds a variety of bacteria and viruses, and mediates attachment of S. mutans to hydroxyapatite on the surface of the tooth. In this study we show that multiallelic copy number variation (CNV) within DMBT1 is extensive across all populations and is predicted to result in between 7-20 scavenger-receptor cysteine-rich (SRCR) domains within each SAG molecule. Direct observation of de novo mutation in multigeneration families suggests these CNVs have a very high mutation rate for a protein-coding locus, with a mutation rate of up to 5% per gamete. Given that the SRCR domains bind S. mutans and hydroxyapatite in the tooth, we investigated the association of sequence diversity at the SAG-binding gene of S. mutans, and DMBT1 CNV. Furthermore, we show that DMBT1 CNV is also associated with a history of agriculture across global populations, suggesting that dietary change as a result of agriculture has shaped the pattern of CNV at DMBT1, and that the DMBT1-S. mutans interaction is a promising model of host-pathogen-culture coevolution in humans.

A comparison of assays for accurate copy number measurement of the low-affinity Fc gamma receptor genes FCGR3A and FCGR3B

The FCGR3 locus encoding the low affinity activating receptor FcγRIII, plays a vital role in immunity triggered by cellular effector and regulatory functions. Copy number of the genes FCGR3A and FCGR3B has previously been reported to affect susceptibility to several autoimmune diseases and chronic inflammatory conditions. However, such genetic association studies often yield inconsistent results; hence require assays that are robust with low error rate. We investigated the accuracy and efficiency in estimating FCGR3 CNV by comparing Sequenom MassARRAY and paralogue ratio test-restriction enzyme digest variant ratio (PRT-REDVR). In addition, since many genetic association studies of FCGR3B CNV were carried out using real-time quantitative PCR, we have also included the evaluation of that method’s performance in estimating the multi-allelic CNV of FCGR3B. The qPCR assay exhibited a considerably broader distribution of signal intensity, potentially introducing error in estimation of copy number and higher false positive rates. Both Sequenom and PRT-REDVR showed lesser systematic bias, but Sequenom skewed towards copy number normal (CN = 2). The discrepancy between Sequenom and PRT-REDVR might be attributed either to batch effects noise in individual measurements. Our study suggests that PRT-REDVR is more robust and accurate in genotyping the CNV of FCGR3, but highlights the needs of multiple independent assays for extensive validation when performing a genetic association study with multi-allelic CNVs.