Skip to main content

Advertisement

A clinical survey of mosaic single nucleotide variants in disease-causing genes detected by exome sequencing

Abstract

Background

Although mosaic variation has been known to cause disease for decades, high-throughput sequencing technologies with the analytical sensitivity to consistently detect variants at reduced allelic fractions have only recently emerged as routine clinical diagnostic tests. To date, few systematic analyses of mosaic variants detected by diagnostic exome sequencing for diverse clinical indications have been performed.

Methods

To investigate the frequency, type, allelic fraction, and phenotypic consequences of clinically relevant somatic mosaic single nucleotide variants (SNVs) and characteristics of the corresponding genes, we retrospectively queried reported mosaic variants from a cohort of ~ 12,000 samples submitted for clinical exome sequencing (ES) at Baylor Genetics.

Results

We found 120 mosaic variants involving 107 genes, including 80 mosaic SNVs in proband samples and 40 in parental/grandparental samples. Average mosaic alternate allele fraction (AAF) detected in autosomes and in X-linked disease genes in females was 18.2% compared with 34.8% in X-linked disease genes in males. Of these mosaic variants, 74 variants (61.7%) were classified as pathogenic or likely pathogenic and 46 (38.3%) as variants of uncertain significance. Mosaic variants occurred in disease genes associated with autosomal dominant (AD) or AD/autosomal recessive (AR) (67/120, 55.8%), X-linked (33/120, 27.5%), AD/somatic (10/120, 8.3%), and AR (8/120, 6.7%) inheritance. Of note, 1.7% (2/120) of variants were found in genes in which only somatic events have been described. Nine genes had recurrent mosaic events in unrelated individuals which accounted for 18.3% (22/120) of all detected mosaic variants in this study. The proband group was enriched for mosaicism affecting Ras signaling pathway genes.

Conclusions

In sum, an estimated 1.5% of all molecular diagnoses made in this cohort could be attributed to a mosaic variant detected in the proband, while parental mosaicism was identified in 0.3% of families analyzed. As ES design favors breadth over depth of coverage, this estimate of the prevalence of mosaic variants likely represents an underestimate of the total number of clinically relevant mosaic variants in our cohort.

Background

Mosaicism is defined by the presence of different genotypic variants among cells of an individual that are derived from the same zygote [1]. Depending on the timing of mutation acquisition, mosaicism may be restricted to the germline (gonadal mosaicism) or non-germ cell tissues (somatic mosaicism) or may involve both (gonosomal mosaicism) [2]. It is estimated that three base substitution mutations arise per cell division in early human embryogenesis [3]. Postzygotic mutations dynamically accumulate and/or are negatively selected during the developmental process [4, 5], rendering each individual a complex mosaic of multiple genetically unique cell lines [1, 4].

Somatic mutations have been well known for their critical role in tumorigenesis [6] and overgrowth syndromes [5]. Mosaic variation has been reported also in asymptomatic individuals. In healthy donors, mutant allele fractions within organ samples ranged from 1.0 to 29.7% [7]. Mosaic variants may be clinically silent for several possible reasons: (1) the mutation is functionally inconsequential, (2) it is restricted to tissues not pertinent to the gene in which the mutation has arisen, (3) it may have occurred after a critical time frame for gene function, or (4) the mutation may be so disadvantageous that selective pressures favor survival and proliferation of cells carrying the reference allele.

Clinically relevant mosaicism is easily recognizable when cutaneous manifestations are present as with segmental neurofibromatosis or McCune-Albright syndrome [8]. However, in the absence of overt skin findings, recognizing underlying mosaicism may present a clinical challenge, particularly when the expressed phenotype deviates substantially from what has been reported in patients with non-mosaic variation. As patients with atypical phenotypes are often referred for exome sequencing (ES), an assessment of the performance of ES for detecting mosaic variation is warranted. Previous studies have evaluated the frequency and type of mosaic variation detectable by ES in specific disease populations, including neurodevelopmental disorders [9], autism [10, 11], and congenital heart disease [12]. However, few systematic analyses of mosaic variants detected by diagnostic ES for diverse clinical indications have been performed [13].

To address this gap in the literature and to lay a framework for additional studies of mosaicism in clinically relevant genes, we present a retrospective review of all reported mosaic variants detected in nearly 12,000 consecutive patients referred for diagnostic ES at Baylor Genetics (BG).

Methods

Study cohort

Laboratory reports for 11,992 consecutive unrelated patients referred for ES were queried to ascertain all clinically relevant mosaic variants reported between Nov 2011 and Aug 2018. Exome analyses were performed as trio ES in 19.8% (n = 2373) and proband-only ES in 80.2% (n = 9619) of cases. One hundred twenty clinical reports with mosaic variants were analyzed for this study; this included 30 cases (25%) analyzed by trio ES and 90 cases (75%) by proband-only ES. Only mosaic variants detected in DNA samples from peripheral blood were analyzed.

Exome sequencing and analysis

ES was performed at BG laboratories as previously described [14, 15] (Additional file 1: Supplementary Methods). The validated ES protocol achieves a mean coverage of 130× with over 95% of targeted regions, including coding and untranslated exons, reaching a minimum coverage of 20×. All samples were concurrently analyzed by the HumanOmni1-Quad or HumanExome-12 v1 array (Illumina) for sample identity confirmation and to screen for copy-number variants and regions of homozygosity. Variant classification was performed in accordance with the American College of Medical Genetics and Genomics (ACMG) and Association for Molecular Pathology (AMP) guidelines for variant interpretation [16]. Mosaic variants of uncertain significance in our cohort that were reported prior to the publication of the ACMG/AMP guidelines were reassessed and classified according to the updated criteria. Common SNPs were filtered out from the analysis.

Mosaic variants reporting/selection criteria

  1. 1.

    Alternate allele fraction (AAF) (mosaic variant reads/total reads) was calculated for each mosaic variant using the data generated by exome sequencing or PCR amplicon-based next-generation sequencing (NGS). For autosomal variants and X-linked variants in females, a variant was considered possibly mosaic if the AAF was less than 36% or greater than 64% by NGS analysis (Additional file 1: Supplementary Methods), while AAF higher than 10% was used as a threshold to identify mosaic variants in X-linked genes in males.

  2. 2.

    Mosaic variants detected by ES were orthogonally confirmed by Sanger sequencing. For mosaic variants ascertained by Sanger sequencing, a substantial and consistent reduction in the electropherogram peak height for the variant allele generated by the Mutation Quantifier function of the Mutation Surveyor software (SoftGenetics, State College, PA, USA) was deemed consistent with mosaicism. Mosaicism detected by Sanger sequencing was also confirmed by subsequent PCR amplicon-based NGS.

  3. 3.

    Only clinically reported mosaic variants were included in the analysis. Mosaic variants detected in disease genes not related to patient phenotype or in candidate disease genes and/or genes of uncertain significance were excluded from the analysis.

  4. 4.

    Mosaic variants detected in non-blood tissues were excluded from the study.

NGS amplicon sequencing

PCR primers targeting mosaic variants were designed using “Primer 3” and synthesized by Sigma Genosys, Woodlands, TX, USA. For each sample, 40 ng of genomic DNA was amplified using Roche’s FastStart kit and/or GC-Rich PCR System for PCR. For SLC6A8 and TUBB (genes with significantly homology to other regions of the genome), long-range PCR (TaKaRa long range PCR kit) followed by nested PCR was used. Amplicon size was checked by gel electrophoresis. PCR products were treated with Exonuclease-Shrimp Alkaline Phosphatase (New England’s BioLabs), and the SPRI bead purified products (Beckman and Coulter Inc. Brea, CA, USA) were used for bar-coding using Illumina compatible index adapters (Sigma Genosys, Woodlands, TX, USA). Barcoded samples were quantified by Qubit (Invitrogen, Life Technologies Corporation, Eugene, OR, USA) and sequenced using the Illumina HiSeq 2500 sequencing system with 100-bp paired-end reads (Illumina, San Diego, CA, USA).

Computational analyses

To better assess the somatic mosaicism burden in ES data, we performed additional computational analyses of AAF distribution for heterozygous single nucleotide variants (SNVs) in 900 ES trios and simulation experiments for evaluating the effect of potential alignment biases.

Results

A total of 120 reported mosaic variants in 107 disease genes were detected in this cohort. Eighty-seven variants were detected by ES and 82 were confirmed by Sanger sequencing (Tables 1 and 2, Fig. 1), whereas 33 mosaic variants (in parental samples) were initially detected by Sanger sequencing. Thirty-two of 33 mosaic variants detected by Sanger sequencing were further validated using PCR amplicon-based NGS analysis (Table 2). For the 87 variants detected by ES, the average coverage at the site of the variant was approximately 202× (range 24–854×) while the average coverage of 32 variants assessed by amplicon-based NGS exceeded 10,000×. Average AAF of variants detected on autosomal chromosomes and in X-linked disease genes in females was 18.2% ± 9.5% (range 3.1–79.7%) compared with 34.8% ± 25.1% (range 10.0–85.0%) for X-linked disease gene variants detected in males. The AAF calculated based on the NGS data was significantly correlated (Spearman rho = 0.93, p = 0) with that quantified by Sanger sequencing (Additional file 2: Figure S1).

Table 1 The 80 mosaic variants detected in the probands
Table 2 The 40 mosaic variants detected in the parental or grandparental samples
Fig. 1
figure1

Overview of the SNV selection strategy

Mosaic variants occurred in genes associated with all types of inheritance, including autosomal dominant (AD) or AD/autosomal recessive (AR) (67/120, 55.8%), X-linked (33/120, 27.5%), AD/somatic (10/120, 8.3%), and AR (8/120, 6.7%) inheritance (Additional file 3: Table S1). Two of the 120 identified mosaic variants involved the IDH1 (MIM 137800) and TET2 (MIM 614286) genes in which only somatic events have been described. Nine genes, including CACNA1A, CREBBP, MTOR, and PIK3CA (n = 3 each), and DDX3X, DNM1, DYRK1A, GRIA3, and KMT2D (n = 2 each) harbored recurrent mosaic events in unrelated individuals. The observed mosaic variants included missense 67.5% (81/120), nonsense 14.1% (17/120), frameshift or in-frame del/dup 13.3% (16/120), and splice 5.0% (6/120) changes (Additional file 3: Table S2). Simulation experiments did not show potential alignment bias of different types of mutations (Additional file 2: Figure S2-S4). Of all single nucleotide substitution variants, 33.7% (35/104) involved CpG sites (Additional file 3: Table S2), and nucleotide C/G>T/A was the most common substitution change (Additional file 3: Table S3).

Mosaic variants in probands

In proband samples, 80 mosaic variants were found in 72 genes in 33 female patients, 45 male patients, and two fetuses. The vast majority were reported in genes associated with AD (47.5%) and X-linked (30.0%) disorders. Mean AAF in proband samples was 32.6% ± 24.4% (n = 15) for X-linked variants in males and 20.2% ± 9.8% (n = 65) for autosomal variants and variants in X-linked disease genes in females (Table 1, Additional file 3: Table S4). For 65 of the 80 probands with mosaic variants, both parental samples were available for inheritance determination. Eight probands had only one parental sample available, and 7 probands had no parental samples available for analysis. The majority of mosaic variants detected in probands (63/65) were deemed de novo due to the absence of the variant in parental DNA by Sanger sequencing. Parental chromosome of origin could not be determined due to a lack of informative SNPs flanking the mosaic variants. In patient 55F, a c.1077dupT (p.L362fs) change in ZMPSTE24 (an autosomal recessive disease gene) was found at an AAF of 80% due to suspected uniparental disomy (UPD) involving chromosome 1. In patient 52F, an inherited c.1129A>T (p.K377*) change in COX15 (also an autosomal recessive disease gene) was found at an AAF of 12% due to suspected segmental UPD involving chromosome 10.

Of the mosaic variants detected in the proband samples, 58.8% (n = 47) were classified as pathogenic (P) or likely pathogenic (LP), and 41.3% (n = 33) as variants of uncertain significance (VOUS). For probands with a mosaic VOUS, 36.4% (12/33) were reported together with one or more non-mosaic P/LP mutations, including de novo or biallelic changes that could explain the core phenotype in four cases, and a heterozygous P/LP variant in an autosomal recessive disease gene in eight cases.

Genotype-phenotype analysis was performed for 47 patients with mosaic P/LP variants (Additional file 4) [17]. Eighty-three percent of the patients had core phenotypes that were consistent with what had been previously reported in association with heterozygous variants, with no evidence of disease attenuation related to the mosaic status of the variant. However, patient 43F carrying a c.38G>A (p.G13D) variant with an AAF 20.8% in HRAS had an apparently attenuated Costello syndrome phenotype, mirroring but less severe than typical for patients with germline mutations in this gene. Three patients had mosaic variants that, even if fully penetrant, would not have explained the full scope of the clinical presentation, including patient 12U with a c.67+2T>G variant in ENG; patient 69M with a c.583C>T (p.R195*) in DMD; and patient 79M with a c.87881T>C (p.V29294A) variant in TTN. We also found three patients with dual molecular diagnoses in whom a second non-mosaic pathogenic variant was considered contributory to the patient’s phenotype (patients 12U, 27F, and 35M). Two patients had multiple mosaic variants detected, including patient 3M who had 17 mosaic variants, only two of which were clinically reported and included in this analysis (see “Discussion”). Patient 12U had eight mosaic variants detected, but only one was found in a known disease-associated gene; the remaining mosaic variants were excluded from this analysis. In both cases, it was unclear whether the mosaic variants had contributed to the patient’s phenotype or if they were a consequence of an underlying predisposition to somatic mutation in the context of a pre-cancerous or cancerous state.

Mosaic variants in parental samples

Forty mosaic variants in 37 genes were detected in 40 parental samples, including one variant detected in a grandparental sample (Table 2). Seven mosaic variants were identified by trio ES analysis whereas the remaining 33 variants were found by Sanger sequencing. Thirty-two of 33 mosaic variants detected by Sanger sequencing were confirmed by PCR-based amplicon NGS. The average AAF of variants detected in autosomal chromosomes and in X-linked disease genes in maternal samples was 14.6 ± 8.0% (Additional file 3: Table S4). One father (120F-Fa) had a mosaic variant with an AAF of 67.8% in the X-linked disease gene, COL4A5, which was detected as a heterozygous change in his daughter. 67.5% (27/40) of mosaic variants detected in parental samples were classified as P/LP in the proband. However, the majority of parents harboring mosaic variants were reported to be clinically unaffected. Only two parents with mosaic variants exhibited phenotypes related to the mosaic change. The father of patient 120F (120F-Fa) with a c.2365A>C (p.T789P) variant in COL4A5 associated with X-linked Alport syndrome (MIM:301050), was reported to have a renal defect. The mother of patient 82M (82M-Mo) was reported to have seizures, muscle weakness, leg weakness, and a clumsy gait; she was found to have a mosaic c.410C>A (p.S137Y) variant in ATP1A3 with an AAF of 14.9%. ATP1A3 is associated with the autosomal dominant disorders, Dystonia 12 (DYT12) [MIM:128235] and cerebellar ataxia, areflexia, pes cavus, optic atrophy, and sensorineural hearing loss (CAPOS) [MIM:601338]. Interestingly, mosaic variants in the CACNA1A gene with AAFs ranging from 15.7 to 29.5% were exclusively detected in parental samples (n = 3). In contrast, mosaic variants in MTOR with comparable AAFs ranging from 16.0 to 32.0% were exclusively detected in proband samples.

Discussion

Each cell division brings with it a risk of a new mutation. Mutations that occur after fertilization lead to the formation of distinct cell lineages or a state of genetic mosaicism. Depending on the functional consequence of the mutation, the timing of its acquisition, and its tissue distribution, the effect of a mosaic variant on patient phenotype can range from negligible to catastrophic. Although mosaic variation has been known to cause disease for decades, high-throughput sequencing technologies with the analytical sensitivity to consistently detect variants at reduced allelic fractions have only recently emerged as routine clinical diagnostic tests. Therefore, empirical studies of the frequency of mosaicism in large patient populations are only now being performed and published. The incidence of mosaic CNVs and aneuploidy found in patients referred for microarray testing has been estimated at 0.55–1% [18, 19]. Without additional verification studies, it is challenging in routine ES analyses to distinguish real somatic variants from apparently de novo heterozygous variants with highly skewed (lower than 0.36) AAF. Therefore, we have focused here only on clinically relevant SNVs. A systematic assessment of the rate of clinically relevant mosaic variant detection in large cohorts of individuals referred for ES with heterogeneous clinical presentations needs more investigations [13].

We endeavored to study the frequency, type, allelic fraction, and phenotypic consequences of reportable mosaic SNVs in a cohort of nearly 12,000 consecutive unrelated patients referred for clinical ES. A total of 120 mosaic variants in 107 established disease genes were detected and reported in either proband (n = 80) or parental (n = 39)/grandparental (n = 1) samples. Mosaic variation was considered definitely or possibly contributory to disease in approximately 1% of 11,992 subjects in this study. Assuming a molecular diagnosis was ascertained in 25% of patients in this cohort [14], an estimated 1.5% of all molecular diagnoses could be attributed to a mosaic variant detected in the proband samples. The fact that these estimates are low relative to other published cohorts was anticipated, as existing reports have studied mosaicism in specific genes [9, 20] or phenotypes [10, 11, 21], and/or have assessed the frequency of rare mosaic variants [11] but not specifically clinically reportable variants.

To assess the phenotypic effects of mosaicism in our cohort, we analyzed the provided clinical information and compared the phenotype of each patient to descriptions in the literature and/or in Online Mendelian Inheritance in Man (OMIM) of individuals with predominantly non-mosaic mutations. In the vast majority of probands with mosaic P/LP variants in AD/X-linked/somatic genes and no confounding factors (e.g., presence of multiple mosaic variants, underlying structural variation), the clinical presentation was not appreciably diminished in severity. In contrast, among parents with mosaic variants, only two (82M-Mo, 120F-Fa) were reported to have a phenotype that could be attributed to the identified mosaic mutation. Excluding mosaic variants detected in X-linked genes in males, a comparison of the AAF of mosaic variants in parental samples (14.6% ± 8.0%) relative to proband samples (20.0% ± 9.8%) showed that unaffected parents with mosaic variants have a significantly lower AAF (p = 0.004, t-test). It is intriguing that mosaic variants with ~ 5% lower AAFs can result in mild or absent phenotypes or can cause clinically significant manifestations. One explanation would be that the impact of any given postzygotic variant is likely to be dependent on the biological function of the gene and the distribution of the mutation in critical tissues. This notion is supported by the mosaic variants found in MTOR, PIK3CA, and CACNA1A in our study. Mosaic variants in MTOR and PIK3CA with AAFs ranging from 12.7 to 24.4% were detected in affected probands with Smith-Kingsmore syndrome [MIM: 616638], Cowden syndrome 5 [MIM: 615108], and/or megalencephaly-capillary malformation-polymicrogyria syndrome [MIM: 602501]. Conversely, mosaic variants in CACNA1A with similar AAFs ranging from 15.7 to 29.5% were all detected in asymptomatic parents. The contrasting severity of phenotypes seen in probands versus clinically unaffected parents highlights the challenge of predicting phenotypic outcomes based on genetic testing alone. It also raises the question of how variant mosaicism should be weighed in the course of variant classification given that both pathogenic and benign effects are possible depending on the clinical context in which the variant is detected.

Interestingly, recurrent mosaic variants in a subset of 9 genes: MTOR, CREBBP, CACNA1A, DDX3X, DNM1, DYRK1A, GRIA3, KMT2D, and PIK3CA accounted for 18.3% (22/120) of all detected mosaic variants in the analyzed cohort. Mosaic variants in several of these genes have been reported previously in the literature: MTOR [11], CREBBP [22], CACNA1A [23], DNM1 [24], KMT2D [25], and PIK3CA [26]. In some cases, e.g., the MTOR and PIK3CA genes, somatic variants are the predominant or the only form of disease-causing mutation described in affected individuals. We have also noted that 10 (12.5%) of the 80 de novo mosaic variants detected in the proband samples were found in a gene associated with the Ras or PI3K-AKT-mTOR pathway, including one variant each in BRAF, NF1, HRAS, and KRAS, and three variants in PIK3CA and MTOR. Heterozygous variants in the same six genes were reported in less than 1% of the entire cohort, indicating that mosaic variation is disproportionately likely to affect this pathway. In fact, mosaic events in this pathway have been commonly observed [27]. The reason for enrichment of mosaicism in the Ras or PI3K-AKT-mTOR signaling pathway is unclear; possible explanations include (1) preferential expansion of hematologic clones with variants in these genes increasing the likelihood of mosaic variant detection, (2) high penetrance of mosaic variants in Ras pathway genes relative to other genes, and (3) a preponderance of intragenic mutation-prone residues.

The recognition that certain genes are more prone to pathogenic postzygotic mutation critically informs recurrence risk counseling and enables optimization of test development and data interpretation in the diagnostic lab setting. Panel-based tests targeting genes with recurrent mosaic variants should have sufficient depth of coverage and, to account for the risk of parental mosaicism, should include recommendations for parental testing. AAF filters are often utilized for comprehensive genomic assays such as exome and whole genome sequencing to exclude variants that are likely to represent sequencing artifact, a practice that can preclude detection of low-level mosaicism. Even with an average ES read depth of 130×, mosaic variants with AAF of less than 10% may be filtered out and excluded from review. For these methodologies, relaxing AAF filters for a defined subset of phenotypically relevant genes in which recurrent mosaic events are known to occur may help to optimize mosaic variant detection. Additionally, testing of tissues distant from the hematopoietic lineage (e.g., urine or hair follicles) could be performed to confirm mosaic status [7].

Adding to the complexity of mosaic variant interpretation, several patients in our cohort were found to harbor more than one mosaic variant. One patient (12U) with multiple congenital malformations was found to have compound heterozygous variants in RAD51C, a gene associated with Fanconi anemia [28], a mosaic VOUS in ENG, and seven additional mosaic variants in genes with no definitive disease association. Genomic instability resulting from spontaneous chromosome breakage is a hallmark of FA [29] and previous studies have shown an increased risk of mosaic copy-number and structural variants in affected individuals [30]. However, the impact of underlying FA on acquisition of somatic single nucleotide and small insertion/deletion variants has not been clearly elucidated. Therefore, although likely, the mosaic variants detected in this patient cannot be unequivocally attributed to the FA diagnosis. Multiple mosaic variants (n = 17) were also detected in patient 3M referred for ES with a history of malignant astrocytoma, myelodysplasia, and dysmorphic features. The mosaic mutations detected in this individual were likely related to the patient’s recent history of myelodysplastic syndrome. Although the phenomenon of mutation acquisition in pre-cancerous and cancerous states is not novel [31], multiple mosaic events stemming from malignancy can be an unexpected finding on assays like ES that are generally performed for the detection of germline, rather than somatic mutations. These findings are also challenging from the standpoint of clinical follow-up, as guidelines do not exist to direct management of incidentally ascertained cancer variants in individuals without a known malignancy.

Finally, we have noted that SNV mosaicism can also be explained by chromosomal abnormalities. Patient 52F with developmental delay and microcephaly was found to have a pathogenic variant in the COX15 gene detected at an AAF of 12%. Analysis of the parental samples for the pathogenic change indicated that the father was heterozygous and the mother was negative for the variant. Due to the unexpectedly low AAF in the proband of the purportedly inherited COX15 variant, review of the SNP array data was performed and the mosaic maternal uniparental disomy of distal chromosome 10q encompassing the COX15 gene was found. In a second case, patient 55F with macrocephaly, dysmorphic features, and digital anomalies was found to have a mosaic pathogenic variant in ZMPSTE24 at an AAF of 80%. The pathogenic variant was found to be heterozygous in the mother and negative in the father. Analysis of the SNP array data again revealed mosaic copy neutral AOH suspicious for UPD involving chromosome 1 and encompassing the ZMPSTE24 gene, which presumably served as the “second hit” for the autosomal recessive disorder.

The many variables that complicate mosaic variant interpretation can also be leveraged in research studies to make inferences about variant pathogenicity and to provide insights into gene function. For example, from the observation that activating mutations in GNAS (associated with McCune-Albright syndrome, OMIM 174800) are detected only in the mosaic state, one can infer that constitutional activating mutations in this gene are incompatible with life [8, 32]. It is plausible that studies of affected individuals, including analyses of AAF by tissue type, would help to define key aspects of gene function, including after what critical developmental period the mutation must occur to ensure viability. For example, conditional PIK3CA activation in mouse cortex showed that abnormal mTOR activation in excitatory neurons and glia, but not interneurons, is sufficient for abnormal cortical overgrowth [33].

Although our cohort is comprised of nearly 12,000 families and we have detected and reported 120 mosaic mutations, only a minority of individuals were found to have mosaic variants in the same gene, which limits our ability to draw conclusions about gene function from analysis of mosaic variation in this cohort specifically. Moreover, causative mutations may be restricted to brain or other tissues that are not commonly studied sources of DNA [34]. As such, additional studies dedicated to assessing mosaicism including larger cohorts of affected and unaffected individuals will be necessary to accumulate the evidence needed to make broad conclusions about gene function based on mosaic variation in the population. Such studies may also allow the use of quantitative information, such as AAF, to predict clinical phenotype, particularly if multiple tissues can be analyzed. Finally, single-cell sequencing will permit a more accurate evaluation of the role of somatic mutations in neurodevelopmental disorders and during normal brain development [35].

Conclusions

In summary, in our cohort of nearly 12,000 patients/families referred for clinical diagnostic ES, mosaic variants considered likely or definitively contributory to phenotype were detected in approximately 1.5% of probands in whom a molecular diagnosis was ascertained. Parental mosaicism was identified in 0.3% of families analyzed. We observed that certain genes, pathways, and even individuals were prone to mosaic variation and that SNV mosaicism can be an indication of underlying structural variation. Since clinical ES by design favors breadth over depth of coverage and only blood samples were analyzed in this study, this analysis likely underestimates the true frequency of clinically relevant mosaicism in our cohort. As sequencing strategies evolve and directed efforts to detect mosaicism are implemented, an increased contribution of mosaic variants to genetic disease will undoubtedly be uncovered.

Availability of data and materials

The datasets supporting the conclusions of this article are included within the article and its additional files. Our raw data cannot be submitted to publicly available databases because the patient families were not consented for sharing their raw data, which can potentially identify the individuals.

Abbreviations

AAF:

Alternate allele fraction

AD:

Autosomal dominant

AOH:

Absence of heterozygosity

AR:

Autosomal recessive

ES:

Exome sequencing

NGS:

Next-generation sequencing

OMIM:

Online Mendelian Inheritance in Man

P:

Pathogenic

SNV:

Single nucleotide variant

UPD:

Uniparental disomy

VOUS:

Variants of uncertain significance

XL:

X-linked

References

  1. 1.

    Forsberg LA, Gisselsson D, Dumanski JP. Mosaicism in health and disease - clones picking up speed. Nat Rev Genet. 2017;18:128–42.

  2. 2.

    Biesecker LG, Spinner NB. A genomic view of mosaicism and human disease. Nat Rev Genet. 2013;14:307–20.

  3. 3.

    Ju YS, Martincorena I, Gerstung M, Petljak M, Alexandrov LB, Rahbari R, et al. Somatic mutations reveal asymmetric cellular dynamics in the early human embryo. Nature. 2017;543:714–8.

  4. 4.

    Campbell IM, Shaw CA, Stankiewicz P, Lupski JR. Somatic mosaicism: implications for disease and transmission genetics. Trends Genet. 2015;31:382–92.

  5. 5.

    Erickson RP. Recent advances in the study of somatic mosaicism and diseases other than cancer. Curr Opin Genet Dev. 2014;26:73–8.

  6. 6.

    Tate JG, Bamford S, Jubb HC, Sondka Z, Beare DM, Bindal N, et al. COSMIC: the catalogue of somatic mutations in Cancer. Nucleic Acids Res. 2018;47:D941–7.

  7. 7.

    Huang AY, Xu X, Ye AY, Wu Q, Yan L, Zhao B, et al. Postzygotic single-nucleotide mosaicisms in whole-genome sequences of clinically unremarkable individuals. Cell Res. 2014;24:1311–27.

  8. 8.

    Happle R. Mosaicism in human skin. Understanding the patterns and mechanisms. Arch Dermatol. 1993;129:1460–70.

  9. 9.

    Stosser MB, Lindy AS, Butler E, Retterer K, Piccirillo-Stosser CM, Richard G, et al. High frequency of mosaic pathogenic variants in genes causing epilepsy-related neurodevelopmental disorders. Genet Med. 2018;20:403–10.

  10. 10.

    Freed D, Pevsner J. The contribution of mosaic variants to autism spectrum disorder. PLoS Genet. 2016;12(9):e1006245.

  11. 11.

    Lim ET, Uddin M, De Rubeis S, Chan Y, Kamumbu AS, Zhang X, et al. Rates, distribution and implications of postzygotic mosaic mutations in autism spectrum disorder. Nat Neurosci. 2017;20:1217–24.

  12. 12.

    Manheimer KB, Richter F, Edelmann LJ, D'Souza SL, Shi L, Shen Y, et al. Robust identification of mosaic variants in congenital heart disease. Hum Genet. 2018;137:183–93.

  13. 13.

    Acuna-Hidalgo R, Bo T, Kwint MP, van de Vorst M, Pinelli M, Veltman JA, et al. Post-zygotic point mutations are an underrecognized source of de novo genomic variation. Am J Hum Genet. 2015;97:67–74.

  14. 14.

    Yang Y, Muzny DM, Xia F, Niu Z, Person R, Ding Y, et al. Molecular findings among patients referred for clinical whole-exome sequencing. JAMA. 2014;312:1870–9.

  15. 15.

    Yang Y, Muzny DM, Reid JG, Bainbridge MN, Willis A, Ward PA, et al. Clinical whole-exome sequencing for the diagnosis of mendelian disorders. N Engl J Med. 2013;369:1502–11.

  16. 16.

    Richards S, Aziz N, Bale S, Bick D, Das S, Gastier-Foster J, et al. Standards and guidelines for the interpretation of sequence variants: a joint consensus recommendation of the American College of Medical Genetics and Genomics and the Association for Molecular Pathology. Genet Med. 2015;17:405–24.

  17. 17.

    Pinero J, Bravo A, Queralt-Rosinach N, Gutierrez-Sacristan A, Deu-Pons J, Centeno E, et al. DisGeNET: a comprehensive platform integrating information on human disease-associated genes and variants. Nucleic Acids Res. 2017;45:D833–9.

  18. 18.

    Conlin LK, Thiel BD, Bonnemann CG, Medne L, Ernst LM, Zackai EH, et al. Mechanisms of mosaicism, chimerism and uniparental disomy identified by single nucleotide polymorphism array analysis. Hum Mol Genet. 2010;19:1263–75.

  19. 19.

    Pham J, Shaw C, Pursley A, Hixson P, Sampath S, Roney E, et al. Somatic mosaicism detected by exon-targeted, high-resolution aCGH in 10,362 consecutive cases. Eur J Hum Genet. 2014;22:969–78.

  20. 20.

    Huisman SA, Redeker EJ, Maas SM, Mannens MM, Hennekam RC. High rate of mosaicism in individuals with Cornelia de Lange syndrome. J Med Genet. 2013;50:339–44.

  21. 21.

    Krupp DR, Barnard RA, Duffourd Y, Evans SA, Mulqueen RM, Bernier R, et al. Exonic mosaic mutations contribute risk for autism spectrum disorder. Am J Hum Genet. 2017;101:369–90.

  22. 22.

    Bartsch O, Kress W, Kempf O, Lechno S, Haaf T, Zechner U. Inheritance and variable expression in Rubinstein-Taybi syndrome. Am J Med Genet Part A. 2010;152a:2254–61.

  23. 23.

    Consortium EK. De novo mutations in SLC1A2 and CACNA1A are important causes of epileptic encephalopathies. Am J Hum Genet. 2016;99:287–98.

  24. 24.

    von Spiczak S, Helbig KL, Shinde DN, Huether R, Pendziwiat M, Lourenco C, et al. DNM1 encephalopathy: a new disease of vesicle fission. Neurology. 2017;89:385–94.

  25. 25.

    Lepri FR, Cocciadiferro D, Augello B, Alfieri P, Pes V, Vancini A, et al. Clinical and neurobehavioral features of three novel Kabuki syndrome patients with mosaic KMT2D mutations and a review of literature. Int J Mol Sci. 2018;19(1):E82.

  26. 26.

    Mirzaa G, Timms AE, Conti V, Boyle EA, Girisha KM, Martin B, et al. PIK3CA-associated developmental disorders exhibit distinct classes of mutations with variable expression and tissue distribution. JCI Insight. 2016;1(9): e87623.

  27. 27.

    McConnell MJ, Moran JV, Abyzov A, Akbarian S, Bae T, Cortes-Ciriano I, et al. Intersection of diverse neuronal genomes and neuropsychiatric disease: the brain somatic mosaicism network. Science. 2017;356:eaaal1641.

  28. 28.

    Jacquinet A, Brown L, Sawkins J, Liu P, Pugash D, Van Allen MI, et al. Expanding the FANCO/RAD51C associated phenotype: cleft lip and palate and lobar holoprosencephaly, two rare findings in Fanconi anemia. Eur J Med Genet. 2018;61:257–61.

  29. 29.

    Kalb R, Neveling K, Nanda I, Schindler D, Hoehn H. Fanconi anemia: causes and consequences of genetic instability. Genome Dyn. 2006;1:218–42.

  30. 30.

    Reina-Castillon J, Pujol R, Lopez-Sanchez M, Rodriguez-Santiago B, Aza-Carmona M, Gonzalez JR, et al. Detectable clonal mosaicism in blood as a biomarker of cancer risk in Fanconi anemia. Blood Adv. 2017;1:319–29.

  31. 31.

    Feinberg AP, Ohlsson R, Henikoff S. The epigenetic progenitor origin of human cancer. Nat Rev Genet. 2006;7:21–33.

  32. 32.

    Happle R. The McCune-Albright syndrome: a lethal gene surviving by mosaicism. Clin Genet. 1986;29:321–4.

  33. 33.

    D'Gama AM, Woodworth MB, Hossain AA, Bizzotto S, Hatem NE, LaCoursiere CM, et al. Somatic mutations activating the mTOR pathway in dorsal telencephalic progenitors cause a continuum of cortical dysplasias. Cell Rep. 2017;21:3754–66.

  34. 34.

    Rodin RE, Walsh CA. Somatic mutation in pediatric neurological diseases. Pediatr Neurol. 2018;87:20–2.

  35. 35.

    Poduri A, Evrony GD, Cai X, Walsh CA. Somatic mutation, genomic variation, and neurological disease. Science. 2013;341:1237758.

Download references

Acknowledgements

We are thankful to our colleagues who provided their expertise that greatly assisted this research work.

Funding

This study is supported by the Institutes of Health (Eunice Kennedy Shriver National Institute of Child Health & Human Development grant R01HD087292 to Dr. Stankiewicz).

Author information

YC, MT, PL, and PS conceived the project and designed the experiments. YC, MT, ES, RG, TC, YF, EG, FG, PW, AB, XW, LM, RX, WB, FX, CE, YY, TG, CS, PL, and PS performed the experiments and analyzed the data. YC, MT, PL, and PS wrote the manuscript. All authors read and approved the final manuscript.

Correspondence to Pengfei Liu or Pawel Stankiewicz.

Ethics declarations

Ethics approval and consent to participate

Our study, which is the review of aggregate clinical data, was approved by Baylor College of Medicine Institutional Review Board with the waiver of informed consent granted. De-identified reporting of demographic and molecular data from this laboratory was approved by the Institutional Review Board at Baylor College of Medicine (H-42680 and H-41191). For clinical testing, exome tests involving minor or fetal sample required informed consent, which was obtained from parents. This research conformed with the principles of the Declaration of Helsinki.

Consent for publication

Not applicable.

Competing interests

BCM and Miraca Holdings Inc. have formed a joint venture with shared ownership and governance of Baylor Genetics (BG), formerly the Baylor Miraca Genetics Laboratories (BMGL), which performs chromosomal microarray analysis and clinical exome sequencing. PW, LM, RX, WB, FX, CS, CE, PL, and PS are current employees of BCM and derive support through a professional service agreement with BG. TC, YF, EG, and FG are current employees of BG. The remaining authors declare that they have no competing interests.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Additional files

Additional file 1:

Description of mosaic alternate allele fraction cutoff, and exome sequencing analysis. (DOCX 14 kb)

Additional file 2:

Figure S1. Correlation of Sanger AAFs with the NGS AAFs of mosaic variants. Figure S2. The relationship of the AAF of heterozygous variants to total read depth. Figure S3. Estimated AAF for randomly selected 13 mosaic variants. Figure S4. Simulation of AAF distribution on SNVs and Indels. Figure S5. The distribution of AAF of all heterozygous variants detected in the 900 ES trios. Figure S6 The AAF distribution of the mosaic variants from Tables 1 and 2. (DOCX 392 kb)

Additional file 3:

Table S1. Summary of mosaic variants and genes according to the inheritance pattern. Table S2. Distribution of mosaic mutation types in probands and parents. Table S3. Spectrum of different single nucleotide substitutions in proband and parent samples. Table S4. Alternate allele fraction of the variants reported in this study. Table S5. Mutations spectrum of apparently de novo heterozygous and mosaic autosomal variants in 900 ES trios. (DOCX 38 kb)

Additional file 4

The list of HPO terms of 80 proband phenotypes. (XLSX 16 kb)

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Keywords

  • AOH
  • CACNA1A
  • CpG site
  • Somatic mosaicism
  • Genotype-phenotype correlation
  • PI3K-AKT-mTOR pathway
  • RASopathies
  • UPD