FOXP1 Haploinsufficiency Contributes to the Development of Congenital Diaphragmatic Hernia.

FOXP1 encodes a transcription factor involved in tissue regulation and cell-type-specific functions. Haploinsufficiency of FOXP1 is associated with a neurodevelopmental disorder: autosomal dominant mental retardation with language impairment with or without autistic features. More recently, heterozygous FOXP1 variants have also been shown to cause a variety of structural birth defects including central nervous system (CNS) anomalies, congenital heart defects, congenital anomalies of the kidney and urinary tract, cryptorchidism, and hypospadias. In this report, we present a previously unpublished case of an individual with congenital diaphragmatic hernia (CDH) who carries an approximately 3.8 Mb deletion. Based on this deletion, and deletions previously reported in two other individuals with CDH, we define a CDH critical region on chromosome 3p13 that includes FOXP1 and four other protein-coding genes. We also provide detailed clinical descriptions of two previously reported individuals with CDH who carry de novo, pathogenic variants in FOXP1 that are predicted to trigger nonsense-mediated mRNA decay. A subset of individuals with putatively deleterious FOXP4 variants has also been shown to develop CDH. Since FOXP proteins function as homo- or heterodimers and the homologs of FOXP1 and FOXP4 are expressed at the same time points in the embryonic mouse diaphragm, they may function together as a dimer, or in parallel as homodimers, to regulate gene expression during diaphragm development. Not all individuals with heterozygous, loss-of-function changes in FOXP1 develop CDH. Hence, we conclude that FOXP1 acts as a susceptibility factor that contributes to the development of CDH in conjunction with other genetic, epigenetic, environmental, and/or stochastic factors.

Targeted accurate RNA consensus sequencing (tARC-seq) reveals mechanisms of replication error affecting SARS-CoV-2 divergence.

RNA viruses, like SARS-CoV-2, depend on their RNA-dependent RNA polymerases (RdRp) for replication, which is error prone. Monitoring replication errors is crucial for understanding the virus's evolution. Current methods lack the precision to detect rare de novo RNA mutations, particularly in low-input samples such as those from patients. Here we introduce a targeted accurate RNA consensus sequencing method (tARC-seq) to accurately determine the mutation frequency and types in SARS-CoV-2, both in cell culture and clinical samples. Our findings show an average of 2.68 × 10-5 de novo errors per cycle with a C > T bias that cannot be solely attributed to APOBEC editing. We identified hotspots and cold spots throughout the genome, correlating with high or low GC content, and pinpointed transcription regulatory sites as regions more susceptible to errors. tARC-seq captured template switching events including insertions, deletions and complex mutations. These insights shed light on the genetic diversity generation and evolutionary dynamics of SARS-CoV-2.

High Clinical Exome Sequencing Diagnostic Rates and Novel Phenotypic Expansions for Nonisolated Microphthalmia, Anophthalmia, and Coloboma.

Purpose: A molecular diagnosis is only made in a subset of individuals with nonisolated microphthalmia, anophthalmia, and coloboma (MAC). This may be due to underutilization of clinical (whole) exome sequencing (cES) and an incomplete understanding of the genes that cause MAC. The purpose of this study is to determine the efficacy of cES in cases of nonisolated MAC and to identify new MAC phenotypic expansions. Methods: We determined the efficacy of cES in 189 individuals with nonisolated MAC. We then used cES data, a validated machine learning algorithm, and previously published expression data, case reports, and animal models to determine which candidate genes were most likely to contribute to the development of MAC. Results: We found the efficacy of cES in nonisolated MAC to be between 32.3% (61/189) and 48.1% (91/189). Most genes affected in our cohort were not among genes currently screened in clinically available ophthalmologic gene panels. A subset of the genes implicated in our cohort had not been clearly associated with MAC. Our analyses revealed sufficient evidence to support low-penetrance MAC phenotypic expansions involving nine of these human disease genes. Conclusions: We conclude that cES is an effective means of identifying a molecular diagnosis in individuals with nonisolated MAC and may identify putatively damaging variants that would be missed if only a clinically available ophthalmologic gene panel was obtained. Our data also suggest that deleterious variants in BRCA2, BRIP1, KAT6A, KAT6B, NSF, RAC1, SMARCA4, SMC1A, and TUBA1A can contribute to the development of MAC.

Hematopoietic stem and progenitor cells confer cross-protective trained immunity in mouse models.

Recent studies suggest that infection reprograms hematopoietic stem and progenitor cells (HSPCs) to enhance innate immune responses upon secondary infectious challenge, a process called "trained immunity." However, the specificity and cell types responsible for this response remain poorly defined. We established a model of trained immunity in mice in response to Mycobacterium avium infection. scRNA-seq analysis revealed that HSPCs activate interferon gamma-response genes heterogeneously upon primary challenge, while rare cell populations expand. Macrophages derived from trained HSPCs demonstrated enhanced bacterial killing and metabolism, and a single dose of recombinant interferon gamma exposure was sufficient to induce similar training. Mice transplanted with influenza-trained HSPCs displayed enhanced immunity against M. avium challenge and vice versa, demonstrating cross protection against antigenically distinct pathogens. Together, these results indicate that heterogeneous responses to infection by HSPCs can lead to long-term production of bone marrow derived macrophages with enhanced function and confer cross-protection against alternative pathogens.

Clinical exome sequencing efficacy and phenotypic expansions involving anomalous pulmonary venous return.

Anomalous pulmonary venous return (APVR) frequently occurs with other congenital heart defects (CHDs) or extra-cardiac anomalies. While some genetic causes have been identified, the optimal approach to genetic testing in individuals with APVR remains uncertain, and the etiology of most cases of APVR is unclear. Here, we analyzed molecular data from 49 individuals to determine the diagnostic yield of clinical exome sequencing (ES) for non-isolated APVR. A definitive or probable diagnosis was made for 8 of those individuals yielding a diagnostic efficacy rate of 16.3%. We then analyzed molecular data from 62 individuals with APVR accrued from three databases to identify novel APVR genes. Based on data from this analysis, published case reports, mouse models, and/or similarity to known APVR genes as revealed by a machine learning algorithm, we identified 3 genes-EFTUD2, NAA15, and NKX2-1-for which there is sufficient evidence to support phenotypic expansion to include APVR. We also provide evidence that 3 recurrent copy number variants contribute to the development of APVR: proximal 1q21.1 microdeletions involving RBM8A and PDZK1, recurrent BP1-BP2 15q11.2 deletions, and central 22q11.2 deletions involving CRKL. Our results suggest that ES and chromosomal microarray analysis (or genome sequencing) should be considered for individuals with non-isolated APVR for whom a genetic etiology has not been identified, and that genetic testing to identify an independent genetic etiology of APVR is not warranted in individuals with EFTUD2-, NAA15-, and NKX2-1-related disorders.

Sox7-positive endothelial progenitors establish coronary arteries and govern ventricular compaction.

The cardiac endothelium influences ventricular chamber development by coordinating trabeculation and compaction. However, the endothelial-specific molecular mechanisms mediating this coordination are not fully understood. Here, we identify the Sox7 transcription factor as a critical cue instructing cardiac endothelium identity during ventricular chamber development. Endothelial-specific loss of Sox7 function in mice results in cardiac ventricular defects similar to non-compaction cardiomyopathy, with a change in the proportions of trabecular and compact cardiomyocytes in the mutant hearts. This phenotype is paralleled by abnormal coronary artery formation. Loss of Sox7 function disrupts the transcriptional regulation of the Notch pathway and connexins 37 and 40, which govern coronary arterial specification. Upon Sox7 endothelial-specific deletion, single-nuclei transcriptomics analysis identifies the depletion of a subset of Sox9/Gpc3-positive endocardial progenitor cells and an increase in erythro-myeloid cell lineages. Fate mapping analysis reveals that a subset of Sox7-null endothelial cells transdifferentiate into hematopoietic but not cardiomyocyte lineages. Our findings determine that Sox7 maintains cardiac endothelial cell identity, which is crucial to the cellular cross-talk that drives ventricular compaction and coronary artery development.

Clinical genome sequencing: Three years' experience at a tertiary children's hospital.

PURPOSE: Genome sequencing (GS) may shorten the diagnostic odyssey for patients, but clinical experience with this assay in nonresearch settings remains limited. Texas Children's Hospital began offering GS as a clinical test to admitted patients in 2020, providing an opportunity to study GS utilization, possibilities for test optimization, and testing outcomes. METHODS: We retrospectively reviewed GS orders for admitted patients for a nearly 3-year period from March 2020 through December 2022. We gathered anonymized clinical data from the electronic health record to answer the study questions. RESULTS: The diagnostic yield over 97 admitted patients was 35%. The majority of GS clinical indications were neurologic or metabolic (61%) and most patients were in intensive care (58%). Tests were often characterized as candidates for intervention/improvement (56%), frequently because of redundancy with prior testing. Patients receiving GS without prior exome sequencing (ES) had higher diagnostic rates (45%) than the cohort as a whole. In 2 cases, GS revealed a molecular diagnosis that is unlikely to be detected by ES. CONCLUSION: The performance of GS in clinical settings likely justifies its use as a first-line diagnostic test, but the incremental benefit for patients with prior ES may be limited.

High molecular diagnostic yields and novel phenotypic expansions involving syndromic anorectal malformations.

Evidence suggests that genetic factors contribute to the development of anorectal malformations (ARMs). However, the etiology of the majority of ARMs cases remains unclear. Exome sequencing (ES) may be underutilized in the diagnostic workup of ARMs due to uncertainty regarding its diagnostic yield. In a clinical database of ~17,000 individuals referred for ES, we identified 130 individuals with syndromic ARMs. A definitive or probable diagnosis was made in 45 of these individuals for a diagnostic yield of 34.6% (45/130). The molecular diagnostic yield of individuals who initially met criteria for VACTERL association was lower than those who did not (26.8% vs 44.1%; p = 0.0437), suggesting that non-genetic factors may play an important role in this subset of syndromic ARM cases. Within this cohort, we identified two individuals who carried de novo pathogenic frameshift variants in ADNP, two individuals who were homozygous for pathogenic variants in BBS1, and single individuals who carried pathogenic or likely pathogenic variants in CREBBP, EP300, FANCC, KDM6A, SETD2, and SMARCA4. The association of these genes with ARMs was supported by previously published cases, and their similarity to known ARM genes as demonstrated using a machine learning algorithm. These data suggest that ES should be considered for all individuals with syndromic ARMs in whom a molecular diagnosis has not been made, and that ARMs represent a low penetrance phenotype associated with Helsmoortel-van der Aa syndrome, Bardet-Biedl syndrome 1, Rubinstein-Taybi syndromes 1 and 2, Fanconi anemia group C, Kabuki syndrome 2, SETD2-related disorders, and Coffin-Siris syndrome 4.

Identifying phenotypic expansions for congenital diaphragmatic hernia plus (CDH+) using DECIPHER data.

Congenital diaphragmatic hernia (CDH) can occur in isolation or in conjunction with other birth defects (CDH+). A molecular etiology can only be identified in a subset of CDH cases. This is due, in part, to an incomplete understanding of the genes that contribute to diaphragm development. Here, we used clinical and molecular data from 36 individuals with CDH+ who are cataloged in the DECIPHER database to identify genes that may play a role in diaphragm development and to discover new phenotypic expansions. Among this group, we identified individuals who carried putatively deleterious sequence or copy number variants affecting CREBBP, SMARCA4, UBA2, and USP9X. The role of these genes in diaphragm development was supported by their expression in the developing mouse diaphragm, their similarity to known CDH genes using data from a previously published and validated machine learning algorithm, and/or the presence of CDH in other individuals with their associated genetic disorders. Our results demonstrate how data from DECIPHER, and other public databases, can be used to identify new phenotypic expansions and suggest that CREBBP, SMARCA4, UBA2, and USP9X play a role in diaphragm development.

Clinical exome sequencing data reveal high diagnostic yields for congenital diaphragmatic hernia plus (CDH+) and new phenotypic expansions involving CDH.

BACKGROUND: Congenital diaphragmatic hernia (CDH) is a life-threatening birth defect that often co-occurs with non-hernia-related anomalies (CDH+). While copy number variant (CNV) analysis is often employed as a diagnostic test for CDH+, clinical exome sequencing (ES) has not been universally adopted. METHODS: We analysed a clinical database of ~12 000 test results to determine the diagnostic yields of ES in CDH+ and to identify new phenotypic expansions. RESULTS: Among the 76 cases with an indication of CDH+, a molecular diagnosis was made in 28 cases for a diagnostic yield of 37% (28/76). A provisional diagnosis was made in seven other cases (9%; 7/76). Four individuals had a diagnosis of Kabuki syndrome caused by frameshift variants in KMT2D. Putatively deleterious variants in ALG12 and EP300 were each found in two individuals, supporting their role in CDH development. We also identified individuals with de novo pathogenic variants in FOXP1 and SMARCA4, and compound heterozygous pathogenic variants in BRCA2. The role of these genes in CDH development is supported by the expression of their mouse homologs in the developing diaphragm, their high CDH-specific pathogenicity scores generated using a previously validated algorithm for genome-scale knowledge synthesis and previously published case reports. CONCLUSION: We conclude that ES should be ordered in cases of CDH+ when a specific diagnosis is not suspected and CNV analyses are negative. Our results also provide evidence in favour of phenotypic expansions involving CDH for genes associated with ALG12-congenital disorder of glycosylation, Rubinstein-Taybi syndrome, Fanconi anaemia, Coffin-Siris syndrome and FOXP1-related disorders.

Selection for or against escape from nonsense mediated decay is a novel signature for the detection of cancer genes.

Escape from nonsense mediated mRNA decay (NMD-) can produce activated or inactivated gene products, and bias in rates of escape can identify functionally important genes in germline disease. We hypothesized that the same would be true of cancer genes, and tested for NMD- bias within The Cancer Genome Atlas pan-cancer somatic mutation dataset. We identify 29 genes that show significantly elevated or suppressed rates of NMD-. This novel approach to cancer gene discovery reveals genes not previously cataloged as potentially tumorigenic, and identifies many potential driver mutations in known cancer genes for functional characterization.

A joint modeling approach for longitudinal microbiome data improves ability to detect microbiome associations with disease.

Changes in the composition of the microbiome over time are associated with myriad human illnesses. Unfortunately, the lack of analytic techniques has hindered researchers' ability to quantify the association between longitudinal microbial composition and time-to-event outcomes. Prior methodological work developed the joint model for longitudinal and time-to-event data to incorporate time-dependent biomarker covariates into the hazard regression approach to disease outcomes. The original implementation of this joint modeling approach employed a linear mixed effects model to represent the time-dependent covariates. However, when the distribution of the time-dependent covariate is non-Gaussian, as is the case with microbial abundances, researchers require different statistical methodology. We present a joint modeling framework that uses a negative binomial mixed effects model to determine longitudinal taxon abundances. We incorporate these modeled microbial abundances into a hazard function with a parameterization that not only accounts for the proportional nature of microbiome data, but also generates biologically interpretable results. Herein we demonstrate the performance improvements of our approach over existing alternatives via simulation as well as a previously published longitudinal dataset studying the microbiome during pregnancy. The results demonstrate that our joint modeling framework for longitudinal microbiome count data provides a powerful methodology to uncover associations between changes in microbial abundances over time and the onset of disease. This method offers the potential to equip researchers with a deeper understanding of the associations between longitudinal microbial composition changes and disease outcomes. This new approach could potentially lead to new diagnostic biomarkers or inform clinical interventions to help prevent or treat disease.

Bayesian modelling of high-throughput sequencing assays with malacoda.

NGS studies have uncovered an ever-growing catalog of human variation while leaving an enormous gap between observed variation and experimental characterization of variant function. High-throughput screens powered by NGS have greatly increased the rate of variant functionalization, but the development of comprehensive statistical methods to analyze screen data has lagged. In the massively parallel reporter assay (MPRA), short barcodes are counted by sequencing DNA libraries transfected into cells and the cell's output RNA in order to simultaneously measure the shifts in transcription induced by thousands of genetic variants. These counts present many statistical challenges, including overdispersion, depth dependence, and uncertain DNA concentrations. So far, the statistical methods used have been rudimentary, employing transformations on count level data and disregarding experimental and technical structure while failing to quantify uncertainty in the statistical model. We have developed an extensive framework for the analysis of NGS functionalization screens available as an R package called malacoda (available from github.com/andrewGhazi/malacoda). Our software implements a probabilistic, fully Bayesian model of screen data. The model uses the negative binomial distribution with gamma priors to model sequencing counts while accounting for effects from input library preparation and sequencing depth. The method leverages the high-throughput nature of the assay to estimate the priors empirically. External annotations such as ENCODE data or DeepSea predictions can also be incorporated to obtain more informative priors-a transformative capability for data integration. The package also includes quality control and utility functions, including automated barcode counting and visualization methods. To validate our method, we analyzed several datasets using malacoda and alternative MPRA analysis methods. These data include experiments from the literature, simulated assays, and primary MPRA data. We also used luciferase assays to experimentally validate several hits from our primary data, as well as variants for which the various methods disagree and variants detectable only with the aid of external annotations.

Increased Moraxella and Streptococcus species abundance after severe bronchiolitis is associated with recurrent wheezing.

BACKGROUND: The role of the airway microbiome in the development of recurrent wheezing and asthma remains uncertain, particularly in the high-risk group of infants hospitalized for bronchiolitis. OBJECTIVE: We sought to examine the relation of the nasal microbiota at bronchiolitis-related hospitalization and 3 later points to the risk of recurrent wheezing by age 3 years. METHODS: In 17 US centers researchers collected clinical data and nasal swabs from infants hospitalized for bronchiolitis. Trained parents collected nasal swabs 3 weeks after hospitalization and, when healthy, during the summer and 1 year after hospitalization. We applied 16S rRNA gene sequencing to all nasal swabs. We used joint modeling to examine the relation of longitudinal nasal microbiota abundances to the risk of recurrent wheezing. RESULTS: Among 842 infants hospitalized for bronchiolitis, there was 88% follow-up at 3 years, and 31% had recurrent wheezing. The median age at enrollment was 3.2 months (interquartile range, 1.7-5.8 months). In joint modeling analyses adjusting for 16 covariates, including viral cause, a 10% increase in relative abundance of Moraxella or Streptococcus species 3 weeks after day 1 of hospitalization was associated with an increased risk of recurrent wheezing (hazard ratio [HR] of 1.38 and 95% high-density interval [HDI] of 1.11-1.85 and HR of 1.76 and 95% HDI of 1.13-3.19, respectively). Increased Streptococcus species abundance the summer after hospitalization was also associated with a greater risk of recurrent wheezing (HR, 1.76; 95% HDI, 1.15-3.27). CONCLUSIONS: Enrichment of Moraxella or Streptococcus species after bronchiolitis hospitalization was associated with recurrent wheezing by age 3 years, possibly providing new avenues to ameliorate the long-term respiratory outcomes of infants with severe bronchiolitis.

Validation Studies for Single Circulating Trophoblast Genetic Testing as a Form of Noninvasive Prenatal Diagnosis.

It has long been appreciated that genetic analysis of fetal or trophoblast cells in maternal blood could revolutionize prenatal diagnosis. We implemented a protocol for single circulating trophoblast (SCT) testing using positive selection by magnetic-activated cell sorting and single-cell low-coverage whole-genome sequencing to detect fetal aneuploidies and copy-number variants (CNVs) at ∼1 Mb resolution. In 95 validation cases, we identified on average 0.20 putative trophoblasts/mL, of which 55% were of high quality and scorable for both aneuploidy and CNVs. We emphasize the importance of analyzing individual cells because some cells are apoptotic, in S-phase, or otherwise of poor quality. When two or more high-quality trophoblast cells were available for singleton pregnancies, there was complete concordance between all trophoblasts unless there was evidence of confined placental mosaicism. SCT results were highly concordant with available clinical data from chorionic villus sampling (CVS) or amniocentesis procedures. Although determining the exact sensitivity and specificity will require more data, this study further supports the potential for SCT testing to become a diagnostic prenatal test.

Functionalization of CD36 cardiovascular disease and expression associated variants by interdisciplinary high throughput analysis.

CD36 is a platelet membrane glycoprotein whose engagement with oxidized low-density lipoprotein (oxLDL) results in platelet activation. The CD36 gene has been associated with platelet count, platelet volume, as well as lipid levels and CVD risk by genome-wide association studies. Platelet CD36 expression levels have been shown to be associated with both the platelet oxLDL response and an elevated risk of thrombo-embolism. Several genomic variants have been identified as associated with platelet CD36 levels, however none have been conclusively demonstrated to be causative. We screened 81 expression quantitative trait loci (eQTL) single nucleotide polymorphisms (SNPs) associated with platelet CD36 expression by a Massively Parallel Reporter Assay (MPRA) and analyzed the results with a novel Bayesian statistical method. Ten eQTLs located 13kb to 55kb upstream of the CD36 transcriptional start site of transcript ENST00000309881 and 49kb to 92kb upstream of transcript ENST00000447544, demonstrated significant transcription shifts between their minor and major allele in the MPRA assay. Of these, rs2366739 and rs1194196, separated by only 20bp, were confirmed by luciferase assay to alter transcriptional regulation. In addition, electromobility shift assays demonstrated differential DNA:protein complex formation between the two alleles of this locus. Furthermore, deletion of the genomic locus by CRISPR/Cas9 in K562 and Meg-01 cells results in upregulation of CD36 transcription. These data indicate that we have identified a variant that regulates expression of CD36, which in turn affects platelet function. To assess the clinical relevance of our findings we used the PhenoScanner tool, which aggregates large scale GWAS findings; the results reinforce the clinical relevance of our variants and the utility of the MPRA assay. The study demonstrates a generalizable paradigm for functional testing of genetic variants to inform mechanistic studies, support patient management and develop precision therapies.

Copy number variant and runs of homozygosity detection by microarrays enabled more precise molecular diagnoses in 11,020 clinical exome cases.

BACKGROUND: Exome sequencing (ES) has been successfully applied in clinical detection of single nucleotide variants (SNVs) and small indels. However, identification of copy number variants (CNVs) using ES data remains challenging. The purpose of this study is to understand the contribution of CNVs and copy neutral runs of homozygosity (ROH) in molecular diagnosis of patients referred for ES. METHODS: In a cohort of 11,020 consecutive ES patients, an Illumina SNP array analysis interrogating mostly coding SNPs was performed as a quality control (QC) measurement and for CNV/ROH detection. Among these patients, clinical chromosomal microarray analysis (CMA) was performed at Baylor Genetics (BG) on 3229 patients, either before, concurrently, or after ES. We retrospectively analyzed the findings from CMA and the QC array. RESULTS: The QC array can detect ~ 70% of pathogenic/likely pathogenic CNVs (PCNVs) detectable by CMA. Out of the 11,020 ES cases, the QC array identified PCNVs in 327 patients and uniparental disomy (UPD) disorder-related ROH in 10 patients. The overall PCNV/UPD detection rate was 5.9% in the 3229 ES patients who also had CMA at BG; PCNV/UPD detection rate was higher in concurrent ES and CMA than in ES with prior CMA (7.2% vs 4.6%). The PCNVs/UPD contributed to the molecular diagnoses in 17.4% (189/1089) of molecularly diagnosed ES cases with CMA and were estimated to contribute in 10.6% of all molecularly diagnosed ES cases. Dual diagnoses with both PCNVs and SNVs were detected in 38 patients. PCNVs affecting single recessive disorder genes in a compound heterozygous state with SNVs were detected in 4 patients, and homozygous deletions (mostly exonic deletions) were detected in 17 patients. A higher PCNV detection rate was observed for patients with syndromic phenotypes and/or cardiovascular abnormalities. CONCLUSIONS: Our clinical genomics study demonstrates that detection of PCNV/UPD through the QC array or CMA increases ES diagnostic rate, provides more precise molecular diagnosis for dominant as well as recessive traits, and enables more complete genetic diagnoses in patients with dual or multiple molecular diagnoses. Concurrent ES and CMA using an array with exonic coverage for disease genes enables most effective detection of both CNVs and SNVs and therefore is recommended especially in time-sensitive clinical situations.

Interchromosomal template-switching as a novel molecular mechanism for imprinting perturbations associated with Temple syndrome.

BACKGROUND: Intrachromosomal triplications (TRP) can contribute to disease etiology via gene dosage effects, gene disruption, position effects, or fusion gene formation. Recently, post-zygotic de novo triplications adjacent to copy-number neutral genomic intervals with runs of homozygosity (ROH) have been shown to result in uniparental isodisomy (UPD). The genomic structure of these complex genomic rearrangements (CGRs) shows a consistent pattern of an inverted triplication flanked by duplications (DUP-TRP/INV-DUP) formed by an iterative DNA replisome template-switching mechanism during replicative repair of a single-ended, double-stranded DNA (seDNA), the ROH results from an interhomolog or nonsister chromatid template switch. It has been postulated that these CGRs may lead to genetic abnormalities in carriers due to dosage-sensitive genes mapping within the copy-number variant regions, homozygosity for alleles at a locus causing an autosomal recessive (AR) disease trait within the ROH region, or imprinting-associated diseases. METHODS: Here, we report a family wherein the affected subject carries a de novo 2.2-Mb TRP followed by 42.2 Mb of ROH and manifests clinical features overlapping with those observed in association with chromosome 14 maternal UPD (UPD(14)mat). UPD(14)mat can cause clinical phenotypic features enabling a diagnosis of Temple syndrome. This CGR was then molecularly characterized by high-density custom aCGH, genome-wide single-nucleotide polymorphism (SNP) and methylation arrays, exome sequencing (ES), and the Oxford Nanopore long-read sequencing technology. RESULTS: We confirmed the postulated DUP-TRP/INV-DUP structure by multiple orthogonal genomic technologies in the proband. The methylation status of known differentially methylated regions (DMRs) on chromosome 14 revealed that the subject shows the typical methylation pattern of UPD(14)mat. Consistent with these molecular findings, the clinical features overlap with those observed in Temple syndrome, including speech delay. CONCLUSIONS: These data provide experimental evidence that, in humans, triplication can lead to segmental UPD and imprinting disease. Importantly, genotype/phenotype analyses further reveal how a post-zygotically generated complex structural variant, resulting from a replication-based mutational mechanism, contributes to expanding the clinical phenotype of known genetic syndromes. Mechanistically, such events can distort transmission genetics resulting in homozygosity at a locus for which only one parent is a carrier as well as cause imprinting diseases.

Megabase Length Hypermutation Accompanies Human Structural Variation at 17p11.2.

DNA rearrangements resulting in human genome structural variants (SVs) are caused by diverse mutational mechanisms. We used long- and short-read sequencing technologies to investigate end products of de novo chromosome 17p11.2 rearrangements and query the molecular mechanisms underlying both recurrent and non-recurrent events. Evidence for an increased rate of clustered single-nucleotide variant (SNV) mutation in cis with non-recurrent rearrangements was found. Indel and SNV formation are associated with both copy-number gains and losses of 17p11.2, occur up to ∼1 Mb away from the breakpoint junctions, and favor C > G transversion substitutions; results suggest that single-stranded DNA is formed during the genesis of the SV and provide compelling support for a microhomology-mediated break-induced replication (MMBIR) mechanism for SV formation. Our data show an additional mutational burden of MMBIR consisting of hypermutation confined to the locus and manifesting as SNVs and indels predominantly within genes.

Publisher Correction: Non-invasive prenatal sequencing for multiple Mendelian monogenic disorders using circulating cell-free fetal DNA.

In the version of this article originally published, some cases that were presented in Fig. 3 should have been underlined but were not. The appropriate cases have now been underlined. The error has been corrected in the print, PDF and HTML versions of the article.