Copy Number Variations Breakpoints Research Articles

Three-dimensional genome landscape comprehensively reveals patterns of spatial gene regulation in papillary and anaplastic thyroid cancers: a study using representative cell lines for each cancer type

BackgroundSpatial chromatin structure is intricately linked with somatic aberrations, and somatic mutations of various cancer-related genes, termed co-mutations (CoMuts), occur in certain patterns during cancer initiation and progression. The functional mechanisms underlying these genetic events remain largely unclear in thyroid cancer (TC). With discrepant differentiation, papillary thyroid cancer (PTC) and anaplastic thyroid cancer (ATC) differ greatly in characteristics and prognosis. We aimed to reveal the spatial gene alterations and regulations between the two TC subtypes.MethodsWe systematically investigated and compared the spatial co-mutations between ATC (8305C), PTC (BCPAP and TPC-1), and normal thyroid cells (Nthy-ori-3–1). We constructed a framework integrating whole-genome sequencing (WGS), high-throughput chromosome conformation capture (Hi-C), and transcriptome sequencing, to systematically detect the associations between the somatic co-mutations of cancer-related genes, structural variations (SVs), copy number variations (CNVs), and high-order chromatin conformation.ResultsSpatial co-mutation hotspots were enriched around topologically associating domains (TADs) in TC. A common set of 227 boundaries were identified in both ATC and PTC, with significant overlaps between them. The spatial proximities of the co-mutated gene pairs in the two TC types were significantly greater than in the gene-level and overall backgrounds, and ATC cells had higher TAD contact frequency with CoMuts > 10 compared with PTC cells. Compared with normal thyroid cells, in ATC the number of the created novel three-dimensional chromatin structural domains increased by 10%, and the number of shifted TADs decreased by 7%. We found five TAD blocks with CoMut genes/events specific to ATC with certain mutations in genes including MAST-NSUN4, AM129B/TRUB2, COL5A1/PPP1R26, PPP1R26/GPSM1/CCDC183, and PRAC2/DLX4. For the majority of ATC and PTC cells, the HOXA10 and HIF2α signals close to the transcription start sites of CoMut genes within TADs were significantly stronger than those at the background. CNV breakpoints significantly overlapped with TAD boundaries in both TC subtypes. ATCs had more CNV losses overlapping with TAD boundaries, and noncoding SVs involved in intrachromosomal SVs, amplified inversions, and tandem duplication differed between ATC and PTC. TADs with short range were more abundant in ATC than PTC. More switches of A/B compartment types existed in ATC cells compared with PTC. Gene expression was significantly synchronized, and orchestrated by complex epigenetics and regulatory elements.ConclusionChromatin interactions and gene alterations and regulations are largely heterogeneous in TC. CNVs and complex SVs may function in the TC genome by interplaying with TADs, and are largely different between ATC and PTC. Complexity of TC genomes, which are highly organized by 3D genome-wide interactions mediating mutational and structural variations and gene activation, may have been largely underappreciated. Our comprehensive analysis may provide key evidence and targets for more customized diagnosis and treatment of TC.

Cas9-Mediated Nanopore Sequencing Enables Precise Characterization of Structural Variants in CCM Genes.

Deletions in the CCM1, CCM2, and CCM3 genes are a common cause of familial cerebral cavernous malformations (CCMs). In current molecular genetic laboratories, targeted next-generation sequencing or multiplex ligation-dependent probe amplification are mostly used to identify copy number variants (CNVs). However, both techniques are limited in their ability to specify the breakpoints of CNVs and identify complex structural variants (SVs). To overcome these constraints, we established a targeted Cas9-mediated nanopore sequencing approach for CNV detection with single nucleotide resolution. Using a MinION device, we achieved complete coverage for the CCM genes and determined the exact size of CNVs in positive controls. Long-read sequencing for a CCM1 and CCM2 CNV revealed that the adjacent ANKIB1 and NACAD genes were also partially or completely deleted. In addition, an interchromosomal insertion and an inversion in CCM2 were reliably re-identified by long-read sequencing. The refinement of CNV breakpoints by long-read sequencing enabled fast and inexpensive PCR-based variant confirmation, which is highly desirable to reduce costs in subsequent family analyses. In conclusion, Cas9-mediated nanopore sequencing is a cost-effective and flexible tool for molecular genetic diagnostics which can be easily adapted to various target regions.

Improving Jitter Distribution in the Breakpoints of Genome CNVs

The copy number variations (CNVs) are a form of structural genetic changes which are recognized to have an importance for diagnosing human disease. Therefore, accurate estimation of the CNVs using high resolution technologies has been under peer attention in both research and clinical applications during last decades. We propose a more accurate approximation for jitter distribution in the CNVs breakpoints based on the modif ed Bessel function of the second kind and zeroth order. We show that the modif ed distribution allows improving the estimates of the CNVs when the segmental signal-to-noise ratio is small and extremely small.

Identification of primary copy number variations reveal enrichment of Calcium, and MAPK pathways sensitizing secondary sites for autism

BackgroundAutism is a neurodevelopmental condition with genetic heterogeneity. It is characterized by difficulties in reciprocal social interactions with strong repetitive behaviors and stereotyped interests. Copy number variations (CNVs) are genomic structural variations altering the genomic structure either by duplication or deletion. De novo or inherited CNVs are found in 5–10% of autistic subjects with a size range of few kilobases to several megabases. CNVs predispose humans to various diseases by altering gene regulation, generation of chimeric genes, and disruption of the coding region or through position effect. Although, CNVs are not the initiating event in pathogenesis; additional preceding mutations might be essential for disease manifestation. The present study is aimed to identify the primary CNVs responsible for autism susceptibility in healthy cohorts to sensitize secondary-hits. In the current investigation, primary-hit autism gene CNVs are characterized in 1715 healthy cohorts of varying ethnicities across 12 populations using Affymetrix high-resolution array study. Thirty-eight individuals from twelve families residing in Karnataka, India, with the age group of 13–73 years are included for the comparative CNV analysis. The findings are validated against global 179 autism whole-exome sequence datasets derived from Simons Simplex Collection. These datasets are deposited at the Simons Foundation Autism Research Initiative (SFARI) database.ResultsThe study revealed that 34.8% of the subjects carried 2% primary-hit CNV burden with 73 singleton-autism genes in different clusters. Of these, three conserved CNV breakpoints were identified with ARHGAP11B, DUSP22, and CHRNA7 as the target genes across 12 populations. Enrichment analysis of the population-specific autism genes revealed two signaling pathways—calcium and mitogen-activated protein kinases (MAPK) in the CNV identified regions. These impaired pathways affected the downstream cascades of neuronal function and physiology, leading to autism behavior. The pathway analysis of enriched genes unravelled complex protein interaction networks, which sensitized secondary sites for autism. Further, the identification of miRNA targets associated with autism gene CNVs added severity to the condition.ConclusionThese findings contribute to an atlas of primary-hit genes to detect autism susceptibility in healthy cohorts, indicating their impact on secondary sites for manifestation.

Epilepsy subtype-specific copy number burden observed in a genome-wide study of 17 458 subjects.

Cytogenic testing is routinely applied in most neurological centres for severe paediatric epilepsies. However, which characteristics of copy number variants (CNVs) confer most epilepsy risk and which epilepsy subtypes carry the most CNV burden, have not been explored on a genome-wide scale. Here, we present the largest CNV investigation in epilepsy to date with 10 712 European epilepsy cases and 6746 ancestry-matched controls. Patients with genetic generalized epilepsy, lesional focal epilepsy, non-acquired focal epilepsy, and developmental and epileptic encephalopathy were included. All samples were processed with the same technology and analysis pipeline. All investigated epilepsy types, including lesional focal epilepsy patients, showed an increase in CNV burden in at least one tested category compared to controls. However, we observed striking differences in CNV burden across epilepsy types and investigated CNV categories. Genetic generalized epilepsy patients have the highest CNV burden in all categories tested, followed by developmental and epileptic encephalopathy patients. Both epilepsy types also show association for deletions covering genes intolerant for truncating variants. Genome-wide CNV breakpoint association showed not only significant loci for genetic generalized and developmental and epileptic encephalopathy patients but also for lesional focal epilepsy patients. With a 34-fold risk for developing genetic generalized epilepsy, we show for the first time that the established epilepsy-associated 15q13.3 deletion represents the strongest risk CNV for genetic generalized epilepsy across the whole genome. Using the human interactome, we examined the largest connected component of the genes overlapped by CNVs in the four epilepsy types. We observed that genetic generalized epilepsy and non-acquired focal epilepsy formed disease modules. In summary, we show that in all common epilepsy types, 1.5-3% of patients carry epilepsy-associated CNVs. The characteristics of risk CNVs vary tremendously across and within epilepsy types. Thus, we advocate genome-wide genomic testing to identify all disease-associated types of CNVs.

Abstract P6-08-13: Personalized and cost-effective detection of copy number variants by molecular-barcode next-generation sequencing and long-read nanopore sequencing

Abstract Background: Germline copy number variants (CNV) of hereditary breast and ovarian cancer genes are a known class of clinically significant mutations. However, failure to detect CNV is a known limitation of Sanger sequencing and conventional amplicon next-generation sequencing (NGS). Multiplex ligation-dependent probe amplification (MLPA), while robust for CNV detection, is labor-intensive, costly to scale up and has limited breakpoint resolution at exon level only. We evaluated the clinical utility of molecular-barcode NGS and long-read nanopore sequencing in our genetic testing algorithm, specifically on CNV detection and breakpoint characterization of probands and family members. Methods: Both molecular-barcode NGS (BRCA1, BRCA2, TP53, PTEN, PALB2 and CDH1) and MLPA (BRCA1 and BRCA2) were performed in parallel on 200 high-risk breast and/or ovarian cancer probands and 11 CNV positive controls. CNV calling from MLPA and NGS data was performed by Coffalyser.Net and a novel in-house bioinformatics algorithm BAMClipper-SE, respectively. Nanopore amplicon or whole-genome sequencing was attempted to map the exact breakpoint of CNV controls and any newly identified CNV positive probands (up to May 2019). Personalized breakpoint PCR test for specific CNVs was designed and validated on probands and subsequently applied to corresponding family members for carrier testing. Results: Molecular-barcode NGS was 100% sensitive in detecting all BRCA1/2 CNVs from 11 positive controls. A total of 4 BRCA1/2 CNVs were detected by NGS from the 200 probands and matched the corresponding MLPA results. In addition, a PALB2 CNV (exons 4-6 deletion) was also detected by NGS during the parallel comparison. Up to May 2019, a total of 22 probands with detectable CNV in BRCA1, BRCA2, PALB2 or CHEK2 were accrued. Exact CNV breakpoint at nucleotide resolution was identified from 20 probands by nanopore and Sanger sequencing (BRCA1 n=12, BRCA2 n=5, PALB2 n=2 and CHEK2 n=1). High-resolution breakpoint mapping was applicable to inform variant pathogenicity and to enable recurrent CNV identification. Personalized breakpoint PCR designed for the 20 probands were applied to 42 family members, of whom 16 were CNV carriers and 26 were non-carriers. Breakpoint PCR results and any available MLPA results of these family members were 100% concordant. Conclusion: Molecular-barcode NGS allows a streamlined workflow to detect both sequence mutations and CNVs. It proves to be a reliable and cost-effective replacement of BRCA1/2 MLPA in genetic testing algorithm of new probands. Simultaneously the NGS workflow expands CNV detection coverage to additional genes (e.g. PALB2) without costly scale-up of MLPA gene panels. Nanopore sequencing is a practical alternative to MLPA in orthogonal validation of newly detected CNV. The high-resolution breakpoint mapping transforms CNV carrier testing to a personalized, simple and cost-effective breakpoint PCR. Citation Format: Ava Kwong, Chun H Au, Dona N Ho, Elaine YL Wong, Yvonne Chung, Fian BF Law, Cecilia YS Ho, Jiawei Chen, Isabella W Cheuk, Vivian Y Shin, Tsun L Chan, Hextan YS Ngan, Edmond SK Ma. Personalized and cost-effective detection of copy number variants by molecular-barcode next-generation sequencing and long-read nanopore sequencing [abstract]. In: Proceedings of the 2019 San Antonio Breast Cancer Symposium; 2019 Dec 10-14; San Antonio, TX. Philadelphia (PA): AACR; Cancer Res 2020;80(4 Suppl):Abstract nr P6-08-13.

WWOX, the FRA16D gene: A target of and a contributor to genomic instability.

WWOX is one of the largest human genes spanning over 1.11 Mbp in length at chr16q23.1-q23.2 and containing FRA16D, the second most common chromosomal fragile site. FRA16D is a hot spot of genomic instability, prone to breakage and for causing germline and somatic copy number variations (CNVs). Consequentially WWOX is frequent target for deletions in cancer. Esophageal, stomach, colon, bladder, ovarian, and uterine cancers are those most commonly affected by WWOX deep focal deletions. WWOX deletions significantly correlate with various clinicopathological features in esophageal carcinoma. WWOX is also a common target for translocations in multiple myeloma. By mapping R-loop (RNA:DNA hybrid) forming sequences (RFLS) we observe this to be a consistent feature aligning with germline and somatic CNV break points at the edges and core of FRA16D spanning from introns 5 to 8 of WWOX. Germline CNV polymorphisms affecting WWOX are extremely common in humans across different ethnic groups. Importantly, structural variants datasets allowed us to identify a specific hot spot for germline duplications and deletions within intron 5 of WWOX coinciding with the 5' edge of the FRA16D core and various RFLS. Recently, multiple pathogenic CNVs spanning WWOX have been identified associated with neurological conditions such as autism spectrum disorder, infantile epileptic encephalopathies, and other developmental anomalies. Loss of WWOX function has recently been associated with DNA damage repair abnormalities, increased genomic instability, and resistance to chemoradiotherapy. The described observations place WWOX both as a target of and a contributor to genomic instability. Both of these aspects will be discussed in this review.

886: Detection of copy-number variants in expanded carrier screening maximizes identification of cystic fibrosis carriers

Cystic fibrosis is one of the most common autosomal recessive conditions. Medical-society guidelines recommend routine carrier screening for cystic fibrosis via targeted genotyping of 23 frequent single-nucleotide variants (SNVs) and short insertions or deletions (indels) in the CFTR gene. Screening for copy-number variants (CNVs) is recommended only when a reproductive partner is a known carrier. Here we assess the performance and clinical impact of routinely screening for SNVs, indels, and CNVs in a next-generation sequencing (NGS)-based expanded carrier screen (ECS). Pathogenic variants in CFTR from 103,718 patients were discovered via a validated NGS-based ECS. A custom algorithm identified CNVs via relative deviations in NGS read depth. Approximate CNV breakpoints were inferred from the NGS-depth profile. Positive CNVs were orthogonally assessed via multiplex ligation-dependent probe amplification (MLPA). CFTR CNV sensitivity was explored across a range of length scales via in silico simulations. Among carriers of cystic fibrosis, 98.7% had pathogenic SNVs or indels in the CFTR gene (79% had one of the 23 common variants), and the remaining 1.3% harbored a pathogenic CNV spanning at least one exon. In total, we observed 25 unique CNVs with a diversity of breakpoints, suggesting that algorithms must be configured to detect novel CNVs. Further, single-exon deletions were observed for seven different CFTR exons; analysis of confidence scores for these empirical deletions—coupled with extensive simulations—demonstrated that the bioinformatics pipeline was both accurate and robust, even for short CNVs. All CNVs identified via NGS were confirmed via MLPA. CNV detection maximizes identification of cystic fibrosis carriers and can be applied to all patients undergoing an NGS-based ECS test validated to detect these complicated variants. Clinical guidelines recommending screening of only the 23 most frequent variants miss critical identification of carriers and should be revisited.

Single-cell copy number variant detection reveals the dynamics and diversity of adaptation.

Copy number variants (CNVs) are a pervasive source of genetic variation and evolutionary potential, but the dynamics and diversity of CNVs within evolving populations remain unclear. Long-term evolution experiments in chemostats provide an ideal system for studying the molecular processes underlying CNV formation and the temporal dynamics with which they are generated, selected, and maintained. Here, we developed a fluorescent CNV reporter to detect de novo gene amplifications and deletions in individual cells. We used the CNV reporter in Saccharomyces cerevisiae to study CNV formation at the GAP1 locus, which encodes the general amino acid permease, in different nutrient-limited chemostat conditions. We find that under strong selection, GAP1 CNVs are repeatedly generated and selected during the early stages of adaptive evolution, resulting in predictable dynamics. Molecular characterization of CNV-containing lineages shows that the CNV reporter detects different classes of CNVs, including aneuploidies, nonreciprocal translocations, tandem duplications, and complex CNVs. Despite GAP1’s proximity to repeat sequences that facilitate intrachromosomal recombination, breakpoint analysis revealed that short inverted repeat sequences mediate formation of at least 50% of GAP1 CNVs. Inverted repeat sequences are also found at breakpoints at the DUR3 locus, where CNVs are selected in urea-limited chemostats. Analysis of 28 CNV breakpoints indicates that inverted repeats are typically 8 nucleotides in length and separated by 40 bases. The features of these CNVs are consistent with origin-dependent inverted-repeat amplification (ODIRA), suggesting that replication-based mechanisms of CNV formation may be a common source of gene amplification. We combined the CNV reporter with barcode lineage tracking and found that 102–104 independent CNV-containing lineages initially compete within populations, resulting in extreme clonal interference. However, only a small number (18–21) of CNV lineages ever constitute more than 1% of the CNV subpopulation, and as selection progresses, the diversity of CNV lineages declines. Our study introduces a novel means of studying CNVs in heterogeneous cell populations and provides insight into their dynamics, diversity, and formation mechanisms in the context of adaptive evolution.

Genetic and epigenetic changes in clonal descendants of irradiated human fibroblasts

To study delayed genetic and epigenetic radiation effects, which may trigger radiation-induced carcinogenesis, we have established single-cell clones from irradiated and non-irradiated primary human fibroblasts. Stable clones were endowed with the same karyotype in all analyzed metaphases after 20 population doublings (PDs), whereas unstable clones displayed mosaics of normal and abnormal karyotypes. To account for variation in radiation sensitivity, all experiments were performed with two different fibroblast strains. After a single X-ray dose of 2 Gy more than half of the irradiated clones exhibited radiation-induced genome instability (RIGI). Irradiated clones displayed an increased rate of loss of chromosome Y (LOY) and copy number variations (CNVs), compared to controls. CNV breakpoints clustered in specific chromosome regions, in particular 3p14.2 and 7q11.21, coinciding with common fragile sites. CNVs affecting the FHIT gene in FRA3B were observed in independent unstable clones and may drive RIGI. Bisulfite pyrosequencing of control clones and the respective primary culture revealed global hypomethylation of ALU, LINE-1, and alpha-satellite repeats as well as rDNA hypermethylation during in vitro ageing. Irradiated clones showed further reduced ALU and alpha-satellite methylation and increased rDNA methylation, compared to controls. Methylation arrays identified several hundred differentially methylated genes and several enriched pathways associated with in vitro ageing. Methylation changes in 259 genes and the MAP kinase signaling pathway were associated with delayed radiation effects (after 20 PDs). Collectively, our results suggest that both genetic (LOY and CNVs) and epigenetic changes occur in the progeny of exposed cells that were not damaged directly by irradiation, likely contributing to radiation-induced carcinogenesis. We did not observe epigenetic differences between stable and unstable irradiated clones. The fact that the DNA methylation (DNAm) age of clones derived from the same primary culture varied greatly suggests that DNAm age of a single cell (represented by a clone) can be quite different from the DNAm age of a tissue. We propose that DNAm age reflects the emergent property of a large number of individual cells whose respective DNAm ages can be highly variable.

3D genome of multiple myeloma reveals spatial genome disorganization associated with copy number variations

The Hi-C method is widely used to study the functional roles of the three-dimensional (3D) architecture of genomes. Here, we integrate Hi-C, whole-genome sequencing (WGS) and RNA-seq to study the 3D genome architecture of multiple myeloma (MM) and how it associates with genomic variation and gene expression. Our results show that Hi-C interaction matrices are biased by copy number variations (CNVs) and can be used to detect CNVs. Also, combining Hi-C and WGS data can improve the detection of translocations. We find that CNV breakpoints significantly overlap with topologically associating domain (TAD) boundaries. Compared to normal B cells, the numbers of TADs increases by 25% in MM, the average size of TADs is smaller, and about 20% of genomic regions switch their chromatin A/B compartment types. In summary, we report a 3D genome interaction map of aneuploid MM cells and reveal the relationship among CNVs, translocations, 3D genome reorganization, and gene expression regulation.

Analysis of copy number variations in Holstein-Friesian cow genomes based on whole-genome sequence data

Thirty-two whole genome DNA sequences of cows were analyzed to evaluate inter-individual variability in the distribution and length of copy number variations (CNV) and to functionally annotate CNV breakpoints. The total number of deletions per individual varied between 9,731 and 15,051, whereas the number of duplications was between 1,694 and 5,187. Most of the deletions (81%) and duplications (86%) were unique to a single cow. No relation between the pattern of variant sharing and a family relationship or disease status was found. The animal-averaged length of deletions was from 5,234 to 9,145 bp and the average length of duplications was between 7,254 and 8,843 bp. Highly significant inter-individual variation in length and number of CNV was detected for both deletions and duplications. The majority of deletion and duplication breakpoints were located in intergenic regions and introns, whereas fewer were identified in noncoding transcripts and splice regions. Only 1.35 and 0.79% of the deletion and duplication breakpoints were observed within coding regions. A gene with the highest number of deletion breakpoints codes for protein kinase cGMP-dependent type I, whereas the T-cell receptor α constant gene had the most duplication breakpoints. The functional annotation of genes with the largest incidence of deletion/duplication breakpoints identified 87/112 Kyoto Encyclopedia of Genes and Genomes pathways, but none of the pathways were significantly enriched or depleted with breakpoints. The analysis of Gene Ontology (GO) terms revealed that a cluster with the highest enrichment score among genes with many deletion breakpoints was represented by GO terms related to ion transport, whereas the GO term cluster mostly enriched among the genes with many duplication breakpoints was related to binding of macromolecules. Furthermore, when considering the number of deletion breakpoints per gene functional category, no significant differences were observed between the "housekeeping" and "strong selection" categories, but genes representing the "low selection pressure" group showed a significantly higher number of breakpoints.

Efficient CNV breakpoint analysis reveals unexpected structural complexity and correlation of dosage-sensitive genes with clinical severity in genomic disorders.

Genomic disorders are the clinical conditions manifested by submicroscopic genomic rearrangements including copy number variants (CNVs). The CNVs can be identified by array-based comparative genomic hybridization (aCGH), the most commonly used technology for molecular diagnostics of genomic disorders. However, clinical aCGH only informs CNVs in the probe-interrogated regions. Neither orientational information nor the resulting genomic rearrangement structure is provided, which is a key to uncovering mutational and pathogenic mechanisms underlying genomic disorders. Long-range polymerase chain reaction (PCR) is a traditional approach to obtain CNV breakpoint junction, but this method is inefficient when challenged by structural complexity such as often found at the PLP1 locus in association with Pelizaeus-Merzbacher disease (PMD). Here we introduced 'capture and single-molecule real-time sequencing' (cap-SMRT-seq) and newly developed 'asymmetry linker-mediated nested PCR walking' (ALN-walking) for CNV breakpoint sequencing in 49 subjects with PMD-associated CNVs. Remarkably, 29 (94%) of the 31 CNV breakpoint junctions unobtainable by conventional long-range PCR were resolved by cap-SMRT-seq and ALN-walking. Notably, unexpected CNV complexities, including inter-chromosomal rearrangements that cannot be resolved by aCGH, were revealed by efficient breakpoint sequencing. These sequence-based structures of PMD-associated CNVs further support the role of DNA replicative mechanisms in CNV mutagenesis, and facilitate genotype-phenotype correlation studies. Intriguingly, the lengths of gained segments by CNVs are strongly correlated with clinical severity in PMD, potentially reflecting the functional contribution of other dosage-sensitive genes besides PLP1. Our study provides new efficient experimental approaches (especially ALN-walking) for CNV breakpoint sequencing and highlights their importance in uncovering CNV mutagenesis and pathogenesis in genomic disorders.

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.

BackgroundIntraspecies 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.ResultsWe 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.ConclusionsWe 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 materialThe online version of this article (doi:10.1186/s12864-016-3221-1) contains supplementary material, which is available to authorized users.

Ready to clone: CNV detection and breakpoint fine-mapping in breast and ovarian cancer susceptibility genes by high-resolution array CGH

Detection of predisposing copy number variants (CNV) in 330 families affected with hereditary breast and ovarian cancer (HBOC). In order to complement mutation detection with Illumina's TruSight Cancer panel, we designed a customized high-resolution 8×60k array for CGH (aCGH) that covers all 94 genes from the panel. Copy number variants with immediate clinical relevance were detected in 12 families (3.6%). Besides 3 known CNVs in CHEK2, RAD51C, and BRCA1, we identified 3 novel pathogenic CNVs in BRCA1 (deletion of exons 4-13, deletion of exons 12-18) and ATM (deletion exons 57-63) plus an intragenic duplication of BRCA2 (exons 3-11) and an intronic BRCA1 variant with unknown pathogenicity. The precision of high-resolution aCGH enabled straight forward breakpoint amplification of a BRCA1 deletion which subsequently allowed for fast and economic CNV verification in family members of the index patient. Furthermore, we used our aCGH data to validate an algorithm that was able to detect all identified copy number changes from next-generation sequencing (NGS) data. Copy number detection is a mandatory analysis in HBOC families at least if no predisposing mutations were found by sequencing. Currently, high-resolution array CGH is our first choice of method of analysis due to unmatched detection precision. Although it seems possible to detect CNV from sequencing data, there currently is no satisfying tool to do so in a routine diagnostic setting.

Fast detection of deletion breakpoints using quantitative PCR.

The routine detection of large and medium copy number variants (CNVs) is well established. Hemizygotic deletions or duplications in the large Duchenne muscular dystrophy DMD gene responsible for Duchenne and Becker muscular dystrophies are routinely identified using multiple ligation probe amplification and array-based comparative genomic hybridization. These methods only map deleted or duplicated exons, without providing the exact location of breakpoints. Commonly used methods for the detection of CNV breakpoints include long-range PCR and primer walking, their success being limited by the deletion size, GC content and presence of DNA repeats. Here, we present a strategy for detecting the breakpoints of medium and large CNVs regardless of their size. The hemizygous deletion of exons 45-50 in the DMD gene and the large autosomal heterozygous PARK2 deletion were used to demonstrate the workflow that relies on real-time quantitative PCR to narrow down the deletion region and Sanger sequencing for breakpoint confirmation. The strategy is fast, reliable and cost-efficient, making it amenable to widespread use in genetic laboratories.

Alu-mediated diverse and complex pathogenic copy-number variants within human chromosome 17 at p13.3.

Alu repetitive elements are known to be major contributors to genome instability by generating Alu-mediated copy-number variants (CNVs). Most of the reported Alu-mediated CNVs are simple deletions and duplications, and the mechanism underlying Alu-Alu-mediated rearrangement has been attributed to non-allelic homologous recombination (NAHR). Chromosome 17 at the p13.3 genomic region lacks extensive low-copy repeat architecture; however, it is highly enriched for Alu repetitive elements, with a fraction of 30% of total sequence annotated in the human reference genome, compared with the 10% genome-wide and 18% on chromosome 17. We conducted mechanistic studies of the 17p13.3 CNVs by performing high-density oligonucleotide array comparative genomic hybridization, specifically interrogating the 17p13.3 region with ∼150 bp per probe density; CNV breakpoint junctions were mapped to nucleotide resolution by polymerase chain reaction and Sanger sequencing. Studied rearrangements include 5 interstitial deletions, 14 tandem duplications, 7 terminal deletions and 13 complex genomic rearrangements (CGRs). Within the 17p13.3 region, Alu-Alu-mediated rearrangements were identified in 80% of the interstitial deletions, 46% of the tandem duplications and 50% of the CGRs, indicating that this mechanism was a major contributor for formation of breakpoint junctions. Our studies suggest that Alu repetitive elements facilitate formation of non-recurrent CNVs, CGRs and other structural aberrations of chromosome 17 at p13.3. The common observation of Alu-mediated rearrangement in CGRs and breakpoint junction sequences analysis further demonstrates that this type of mechanism is unlikely attributed to NAHR, but rather may be due to a recombination-coupled DNA replicative repair process.

GROM-RD: resolving genomic biases to improve read depth detection of copy number variants.

Amplifications or deletions of genome segments, known as copy number variants (CNVs), have been associated with many diseases. Read depth analysis of next-generation sequencing (NGS) is an essential method of detecting CNVs. However, genome read coverage is frequently distorted by various biases of NGS platforms, which reduce predictive capabilities of existing approaches. Additionally, the use of read depth tools has been somewhat hindered by imprecise breakpoint identification. We developed GROM-RD, an algorithm that analyzes multiple biases in read coverage to detect CNVs in NGS data. We found non-uniform variance across distinct GC regions after using existing GC bias correction methods and developed a novel approach to normalize such variance. Although complex and repetitive genome segments complicate CNV detection, GROM-RD adjusts for repeat bias and uses a two-pipeline masking approach to detect CNVs in complex and repetitive segments while improving sensitivity in less complicated regions. To overcome a typical weakness of RD methods, GROM-RD employs a CNV search using size-varying overlapping windows to improve breakpoint resolution. We compared our method to two widely used programs based on read depth methods, CNVnator and RDXplorer, and observed improved CNV detection and breakpoint accuracy for GROM-RD. GROM-RD is available at http://grigoriev.rutgers.edu/software/.

Rearrangement structure-independent strategy of CNV breakpoint analysis.

Rare copy number variations (CNVs) generated by human genomic rearrangements have been shown to play an important role in pathogenesis of human diseases and cancers. CNV breakpoint analysis can help define genomic location, genetic content and sequence structure of pathogenic CNVs. This process is vital to elucidate CNV mutational mechanism and etiology of CNV-associated disorders. However, it is technically challenging to map CNV breakpoints at base-pair level, especially in the genomic regions with sequence complexity. In this study, we developed a new method of capture and breakpoint approaching sequencing (CBAS) to efficiently obtain CNV breakpoint sequences. This strategy is independent of CNV structures and applicable to various CNV types. As was demonstrated in CNV-associated patients with neurological disorders, CBAS achieved fine mapping of breakpoint sequences for compound deletion, complex duplication, and translocation. Intriguingly, CBAS also revealed unexpected CNV complexity involving long-range DNA rearrangement. Our observations showed that CBAS is an efficient method for obtaining CNV breakpoint sequence and mapping insertional events as well. This method can facilitate the researches on CNV-associated human diseases and cancers. CBAS is also applicable to mapping the integration sites of retrovirus (such as HIV) and transgenes in model organisms.

Comment on “Genomic Hypomethylation in the Human Germline Associates with Selective Structural Mutability in the Human Genome”

Copy number variants (CNVs) are a dynamic feature of the human genome that play important roles in human adaptation and susceptibility to both common and rare disease [1]. The distribution of CNVs in mammalian genomes is nonrandom, and several sequence features have been associated with CNV breakpoints and regions of high structural mutability [2]–[8]. Based on an analysis of DNA methylation patterns in human sperm, Li et al. recently reported a significant relationship between CNVs and hypomethylation in the male germline [9], leading to the suggestion that DNA hypomethylation plays a causative role in the generation of structural variation. Given the potentially profound implications of this report for the study of human disease, we read the findings of Li et al. with great interest. However, after systematically reanalyzing the relationship between CNVs and DNA methylation patterns in sperm, we have identified several cryptic confounders in the data that we believe seriously undermine the conclusions of Li et al. We outline and discuss each of these in detail below.