Multiplex Amplicon Quantification (MAQ), a fast and efficient method for the simultaneous detection of copy number alterations in neuroblastoma

On the core segmentation algorithms of copy number variation detection tools

Shotgun sequencing is a high-throughput method used to detect copy number variants (CNVs). Although there are numerous CNV detection tools based on shotgun sequencing, their quality varies significantly, leading to performance discrepancies. Therefore, we conducted a comprehensive analysis of next-generation sequencing-based CNV detection tools over the past decade. Our findings revealed that the majority of mainstream tools employ similar detection rationale: calculates the so-called read depth signal from aligned sequencing reads and then segments the signal by utilizing either circular binary segmentation (CBS) or hidden Markov model (HMM). Hence, we compared the performance of those two core segmentation algorithms in CNV detection, considering varying sequencing depths, segment lengths and complex types of CNVs. To ensure a fair comparison, we designed a parametrical model using mainstream statistical distributions, which allows for pre-excluding bias correction such as guanine-cytosine (GC) content during the preprocessing step. The results indicate the following key points: (1) Under ideal conditions, CBS demonstrates high precision, while HMM exhibits a high recall rate. (2) For practical conditions, HMM is advantageous at lower sequencing depths, while CBS is more competitive in detecting small variant segments compared to HMM. (3) In case involving complex CNVs resembling real sequencing, HMM demonstrates more robustness compared with CBS. (4) When facing large-scale sequencing data, HMM costs less time compared with the CBS, while their memory usage is approximately equal. This can provide an important guidance and reference for researchers to develop new tools for CNV detection.

Copy number variant analysis and expression profiling of the olfactory receptor-rich 11q11 region in obesity predisposition

Genome-wide copy number surveys associated chromosome 11q11 with obesity. As this is an olfactory receptor-rich region, we hypothesize that genetic variation in olfactory receptor genes might be implicated in the pathogenesis of obesity. Multiplex Amplicon Quantification analysis was applied to screen for copy number variants at chromosome 11q11 in 627 patients with obesity and 330 healthy-weight individuals. A ± 80 kb deletion with an internally 1.3 kb retained segment was identified, covering the three olfactory receptor genes OR4C11, OR4P4, and OR4S2. A significant increase in copy number loss(es) was perceived in our patient cohort (MAF = 27%; p = 0.02). Gene expression profiling in metabolic relevant tissues was performed to evaluate the functional impact of the obesity susceptible locus. All three 11q11 genes were present in visceral and subcutaneous adipose tissue while no expression was perceived in the liver. These results support the 'metabolic system' hypothesis and imply that gene disruption of OR4C11, OR4P4, and OR4S2 will negatively influence energy metabolism, ultimately leading to fat accumulation and obesity. Our study thus demonstrates a role for structural variation within olfactory receptor-rich regions in complex diseases and defines the 11q11 deletion as a risk factor for obesity.

Novel compound heterozygous mutations in OCA2 gene associated with non-syndromic oculocutaneous albinism in a Chinese Han patient: a case report

Background

Oculocutaneous albinism (OCA) is a group of rare genetically heterogeneous disorders. The present study aimed to identify the genetic cause of a Chinese Han family with non-syndromic oculocutaneous albinism (OCA).

Case presentation

Here, we report an 11-month-old male proband from a Chinese Han non-consanguineous family, who presented with milky skin, yellow white hair, nystagmus, astigmatism, and hypermetropia. We performed the targeted next-generation sequencing (NGS) on the proband and identified two novel compound heterozygous variants (c.1865 T > C (p.Leu622Pro) and exons 17–21 deletion) in OCA2 gene associated with OCA type 2 (OCA2, OMIM 203200). Meanwhile, a previously reported heterozygous mutation (c.4805G > A) in MYO7 gene related with Usher syndrome type 1B was found. The online tools SIFT, PolyPhen-2, and Mutation Taster predicted variant c.1865 T > C was probably damaging. The residue p.Leu622 was in a highly conserved region among species by CLUSTALW. Three-dimensional homology model with I-TASSER indicated that p.Leu622Pro variant disturbed the formation of the α-helix, resulting in a random coil structure. The gross deletion (exons 17–21) in OCA2 gene has was not been reported previously. These two novel variants in OCA2 gene were inherited from each parent respectively, after verification by Sanger sequencing and quantitative PCR (qPCR) in the family.

Conclusions

This study indicates the two novel compound heterozygous mutations in OCA2 gene may be responsible for clinical manifestations of OCA2. It expands the mutation spectrum of OCA2 gene and is helpful to screen for large deletions with targeted NGS protocol in monogenic disease. It also assists the genetic counselling, carrier screening and personalized healthcare of the disease.

Electronic supplementary material

The online version of this article (10.1186/s12881-019-0850-7) contains supplementary material, which is available to authorized users.

A Whole Germline BRCA2 Gene Deletion: How to Learn from CNV In Silico Analysis

BRCA1/2 screening in Hereditary Breast and Ovarian Syndrome (HBOC) is an essential step for effective patients’ management. Next-Generation Sequencing (NGS) can rapidly provide high throughput and reliable information about the qualitative and quantitative status of tumor-associated genes. Straightforwardly, bioinformatics methods play a key role in molecular diagnostics pipelines. BRCA1/2 genes were evaluated with our NGS workflow, coupled with Multiplex Amplicon Quantification (MAQ) and Multiplex Ligation-dependent Probe Amplification (MLPA) assays. Variant calling was performed on Amplicon Suite, while Copy Number Variant (CNV) prediction by in house and commercial CNV tools, before confirmatory MAQ/MLPA testing. The germline profile of BRCA genes revealed a unique HBOC pattern. Although variant calling analysis pinpointed heterozygote and homozygote polymorphisms on BRCA1 and BRCA2, respectively, the CNV predicted by our script suggested two conflicting interpretations: BRCA1 duplication and/or BRCA2 deletion. Our commercial software reported a BRCA1 duplication, in contrast with variant calling results. Finally, the MAQ/MLPA assays assessed a whole BRCA2 copy loss. In silico CNV analysis is a time and cost-saving procedure to powerfully identify possible Large Rearrangements using robust and efficient NGS pipelines. Our layout shows as bioinformatics algorithms alone cannot completely and correctly identify whole BRCA1/2 deletions/duplications. In particular, the complete deletion of an entire gene, like in our case, cannot be solved without alternative strategies as MLPA/MAQ. These findings support the crucial role of bioinformatics in deciphering pitfalls within NGS data analysis.

Detection of clinically relevant copy-number variants by exome sequencing in a large cohort of genetic disorders

Purpose:

Copy-number variation is a common source of genomic variation and an important genetic cause of disease. Microarray-based analysis of copy-number variants (CNVs) has become a first-tier diagnostic test for patients with neurodevelopmental disorders, with a diagnostic yield of 10–20%. However, for most other genetic disorders, the role of CNVs is less clear and most diagnostic genetic studies are generally limited to the study of single-nucleotide variants (SNVs) and other small variants. With the introduction of exome and genome sequencing, it is now possible to detect both SNVs and CNVs using an exome- or genome-wide approach with a single test.

Methods:

We performed exome-based read-depth CNV screening on data from 2,603 patients affected by a range of genetic disorders for which exome sequencing was performed in a diagnostic setting.

Results:

In total, 123 clinically relevant CNVs ranging in size from 727 bp to 15.3 Mb were detected, which resulted in 51 conclusive diagnoses and an overall increase in diagnostic yield of ~2% (ranging from 0 to –5.8% per disorder).

Conclusions:

This study shows that CNVs play an important role in a broad range of genetic disorders and that detection via exome-based CNV profiling results in an increase in the diagnostic yield without additional testing, bringing us closer to single-test genomics.

Genet Med advance online publication 27 October 2016

DNA Diagnostics of Hereditary Hearing Loss: A Targeted Resequencing Approach Combined with a Mutation Classification System

Although there are nearly 100 different causative genes identified for nonsyndromic hearing loss (NSHL), Sanger sequencing-based DNA diagnostics usually only analyses three, namely, GJB2, SLC26A4, and OTOF. As this is seen as inadequate, there is a need for high-throughput diagnostic methods to detect disease-causing variations, including single-nucleotide variations (SNVs), insertions/deletions (Indels), and copy-number variations (CNVs). In this study, a targeted resequencing panel for hearing loss was developed including 79 genes for NSHL and selected forms of syndromic hearing loss. One-hundred thirty one presumed autosomal-recessive NSHL (arNSHL) patients of Western-European ethnicity were analyzed for SNVs, Indels, and CNVs. In addition, we established a straightforward variant classification system to deal with the large number of variants encountered. We estimate that combining prescreening of GJB2 with our panel leads to a diagnosis in 25%-30% of patients. Our data show that after GJB2, the most commonly mutated genes in a Western-European population are TMC1, MYO15A, and MYO7A (3.1%). CNV analysis resulted in the identification of causative variants in two patients in OTOA and STRC. One of the major challenges for diagnostic gene panels is assigning pathogenicity for variants. A collaborative database collecting all identified variants from multiple centers could be a valuable resource for hearing loss diagnostics.

Targeted next generation sequencing reveals a novel intragenic deletion of the LAMA2 gene in a patient with congenital muscular dystrophy

Mutations in the LAMA2 gene cause laminin α‑2 (merosin)‑deficient congenital muscular dystrophies, which are autosomal recessive muscle disorders. Laminin α‑2 is widely expressed in the basement membrane of skeletal muscle, the myotendinous junctions and extra‑synaptically at neuromuscular synapses. In the present study, target next‑generation sequencing was used for mutation detection, and polymerase chain reaction (PCR) analysis and Sanger sequencing were used in the identification of small deletions. Subsequently, quantitative PCR (qPCR) was performed to characterize the identified deletion encompassing exon five of the LAMA2 gene. Two causative mutations were identified using target region sequencing which provided the additional information required to facilitate clinical diagnosis. One heterozygous mutation (p. Lys682LysfsX22) was identified and confirmed by Sanger sequencing, and another heterozygous mutation (Exon5del) was found and validated by qPCR. Co‑segregation analysis indicated that the Exon5del mutation originated from the proband's mother and the previously reported frameshift mutation (p. Lys682LysfsX22) was inherited from the proband's father. To the best of our knowledge, the present study was the first to report an entire exon five deletion in the LAMA2 gene.

A pan-European study of the C9orf72 repeat associated with FTLD: geographic prevalence, genomic instability, and intermediate repeats

We assessed the geographical distribution of C9orf72 G(4) C(2) expansions in a pan-European frontotemporal lobar degeneration (FTLD) cohort (n = 1,205), ascertained by the European Early-Onset Dementia (EOD) consortium. Next, we performed a meta-analysis of our data and that of other European studies, together 2,668 patients from 15 Western European countries. The frequency of the C9orf72 expansions in Western Europe was 9.98% in overall FTLD, with 18.52% in familial, and 6.26% in sporadic FTLD patients. Outliers were Finland and Sweden with overall frequencies of respectively 29.33% and 20.73%, but also Spain with 25.49%. In contrast, prevalence in Germany was limited to 4.82%. In addition, we studied the role of intermediate repeats (7-24 repeat units), which are strongly correlated with the risk haplotype, on disease and C9orf72 expression. In vitro reporter gene expression studies demonstrated significantly decreased transcriptional activity of C9orf72 with increasing number of normal repeat units, indicating that intermediate repeats might act as predisposing alleles and in favor of the loss-of-function disease mechanism. Further, we observed a significantly increased frequency of short indels in the GC-rich low complexity sequence adjacent to the G(4) C(2) repeat in C9orf72 expansion carriers (P < 0.001) with the most common indel creating one long contiguous imperfect G(4) C(2) repeat, which is likely more prone to replication slippage and pathological expansion.

Rare copy number variants in neuropsychiatric disorders: Specific phenotype or not?

From a number of genome-wide association studies it was shown that de novo and/or rare copy number variants (CNVs) are found at an increased frequency in neuropsychiatric diseases. In this study we examined the prevalence of CNVs in six genomic regions (1q21.1, 2p16.3, 3q29, 15q11.2, 15q13.3, and 16p11.2) previously implicated in neuropsychiatric diseases. Hereto, a cohort of four neuropsychiatric disorders (schizophrenia, bipolar disorder, major depressive disorder, and intellectual disability) and control individuals from three different populations was used in combination with Multilpex Amplicon Quantifiaction (MAQ) assays, capable of high resolution (kb range) and custom-tailored CNV detection. Our results confirm the etiological candidacy of the six selected CNV regions for neuropsychiatric diseases. It is possible that CNVs in these regions can result in disturbed brain development and in this way lead to an increased susceptibility for different neuropsychiatric disorders, dependent on additional genetic and environmental factors. Our results also suggest that the neurodevelopmental component is larger in the etiology of schizophrenia and intellectual disability than in mood disorders. Finally, our data suggest that deletions are in general more pathogenic than duplications. Given the high frequency of the examined CNVs (1-2%) in patients of different neuropsychiatric disorders, screening of large cohorts with an affordable and feasible method like the MAQ assays used in this study is likely to result in important progress in unraveling the genetic factors leading to an increased susceptibility for several psychiatric disorders.

Applying microsatellite multiplex PCR analysis (MMPA) for determining allele copy-number status and percentage of normal cells within tumors

The study of somatic genetic alterations in tumors contributes to the understanding and management of cancer. Genetic alterations, such us copy number or copy neutral changes, generate allelic imbalances (AIs) that can be determined using polymorphic markers. Here we report the development of a simple set of calculations for analyzing microsatellite multiplex PCR data from control-tumor pairs that allows us to obtain accurate information not only regarding the AI status of tumors, but also the percentage of tumor-infiltrating normal cells, the locus copy-number status and the mechanism involved in AI. We validated this new approach by re-analyzing a set of Neurofibromatosis type 1-associated dermal neurofibromas and comparing newly generated data with results obtained for the same tumors in a previous study using MLPA, Paralog Ratio Analysis and SNP-array techniques.Microsatellite multiplex PCR analysis (MMPA) should be particularly useful for analyzing specific regions of the genome containing tumor suppressor genes and also for determining the percentage of infiltrating normal cells within tumors allowing them to be sorted before they are analyzed by more expensive techniques.

Alzheimer risk associated with a copy number variation in the complement receptor 1 increasing C3b/C4b binding sites

Two multicentre genome-wide association (GWA) studies provided substantial evidence, implicating the complement receptor 1 gene (CR1) in Alzheimer disease (AD) genetic etiology. CR1 encodes a large transmembrane receptor with a crucial role in the immune complement cascade. We performed a genetic follow-up of the GWA CR1 association in a Flanders-Belgian cohort (n=1883), and investigated the effect of single-nucleotide polymorphisms (SNPs) located in the CR1 locus on AD risk and cerebrospinal fluid (CSF) biomarker levels. We obtained significant association (P(adj)<0.03; odds ratio (OR)=1.24 (95% confidence interval (CI): 1.02-1.51)) for one CR1 risk haplotype, and haplotype association was strongest in individuals carrying apolipoprotein E (APOE) ɛ4 alleles (P(adj)<0.006; OR=1.50 (95% CI: 1.08-2.09)). Also, four SNPs correlated with increased CSF amyloid Aβ₁₋₄₂ levels, suggesting a role for the CR1 protein in Aβ metabolism. Moreover, we quantified a low-copy repeat (LCR)-associated copy number variation (CNV) in CR1, producing different CR1 isoforms, CR1-F and CR1-S, and obtained significant association in carriers of CR1-S. We replicated the CR1 CNV association finding in a French cohort (n=2003) and calculated in the combined cohorts, an OR of 1.32; 95% CI: 1.10-1.59 (P=0.0025). Our data showed that the common AD risk association may well be explained by the presence of CR1-S increasing the number of C3b/C4b and cofactor activity sites and AD risk with 30% in CR1-S carriers. How precisely the different functional role of CR1-S in the immune complement cascade contributes to AD pathogenesis will need additional functional studies.

No evidence for GNAS copy number variants in patients with features of Albright's hereditary osteodystrophy and abnormal platelet Gs activity

Albright's hereditary osteodystrophy (AHO) is characterized by short stature, round face, calcifications, obesity, brachydactyly and intellectual disability. AHO without hormone resistance is called pseudopseudohypoparathyroidism (PPHP), a rare clinical condition difficult to diagnose with highly variable features. PPHP is caused by paternally inherited loss-of-function mutations in the GNAS. Patients with 2q37 microdeletions or HDAC4 mutations are also defined as having an AHO-like phenotype with normal stimulatory G (Gs) function. We have studied 256 patients with AHO features but no other diagnosis. Their platelet Gs activity was determined via the aggregation-inhibition test showing Gs hypo- or hyperfuncton in 24% and 15% of the patients, respectively. Before initiating with detailed (epi)genetic GNAS studies, we here wanted to excluded copy number variants (CNVs) in GNAS as cause of AHO with a novel large-scale screening technique. Multiplex amplicon quantification (MAQ) for CNVs screening was developed for the 20q13.3 region including GNAS and potential long-range imprinting control elements such as STX16. This is the first large-scale GNAS CNV study in patients with common AHO features but no CNVs were detected. In conclusion, CNVs in the GNAS region are not likely to cause an AHO-like phenotype with or without abnormal platelet Gs activity. Future studies will be undertaken to find out whether these AHO patients with abnormal Gs function are characterized by GNAS coding or methylation defects.

DLB and PDD: a role for mutations in dementia and Parkinson disease genes?

Based on the substantial overlap in clinical and pathological characteristics of dementia with Lewy bodies (DLB) and Parkinson disease with dementia (PDD) with Alzheimer disease (AD) and Parkinson disease (PD) we hypothesized that these disorders might share underlying genetic factors. The contribution of both sequence and copy number variants (CNVs) in known AD and PD genes to the genetic etiology of DLB and PDD however is currently unclear. Therefore, we performed a gene-based mutation analysis of all major AD and PD genes in 99 DLB and 75 PDD patients, including familial and sporadic forms, from Flanders, Belgium. Also, copy number variants in APP, SNCA, and PARK2 were determined. In the AD genes we detected proven pathogenic missense mutations in PSEN1 and PSEN2, and 2 novel missense variants in PSEN2 and MAPT. In the PD genes we identified 1 SNCA duplication, the LRRK2 R1441C founder mutation and 4 novel heterozygous missense variants with unknown pathogenicity. Our results suggest a contribution of established AD and PD genes to the genetic etiology of DLB and PDD though to a limited extent. They do support the hypothesis of a genetic overlap between members of the Lewy body disease spectrum, but additional genes still have to exist.

Genome structural variation discovery and genotyping

Comparisons of human genomes show that more base pairs are altered as a result of structural variation - including copy number variation - than as a result of point mutations. Here we review advances and challenges in the discovery and genotyping of structural variation. The recent application of massively parallel sequencing methods has complemented microarray-based methods and has led to an exponential increase in the discovery of smaller structural-variation events. Some global discovery biases remain, but the integration of experimental and computational approaches is proving fruitful for accurate characterization of the copy, content and structure of variable regions. We argue that the long-term goal should be routine, cost-effective and high quality de novo assembly of human genomes to comprehensively assess all classes of structural variation.