Genome sequencing is a sensitive first-line test to diagnose individuals with intellectual disability.

PURPOSE: Individuals with intellectual disability (ID) and/or neurodevelopment disorders (NDDs) are currently investigated with several different approaches in clinical genetic diagnostics. METHODS: We compared the results from 3 diagnostic pipelines in patients with ID/NDD: genome sequencing (GS) first (N = 100), GS as a secondary test (N = 129), or chromosomal microarray (CMA) with or without FMR1 analysis (N = 421). RESULTS: The diagnostic yield was 35% (GS-first), 26% (GS as a secondary test), and 11% (CMA/FMR1). Notably, the age of diagnosis was delayed by 1 year when GS was performed as a secondary test and the cost per diagnosed individual was 36% lower with GS first than with CMA/FMR1. Furthermore, 91% of those with a negative result after CMA/FMR1 analysis (338 individuals) have not yet been referred for additional genetic testing and remain undiagnosed. CONCLUSION: Our findings strongly suggest that genome analysis outperforms other testing strategies and should replace traditional CMA and FMR1 analysis as a first-line genetic test in individuals with ID/NDD. GS is a sensitive, time- and cost-effective method that results in a confirmed molecular diagnosis in 35% of all referred patients.

DNA Methylation Signature for JARID2-Neurodevelopmental Syndrome.

JARID2 (Jumonji, AT Rich Interactive Domain 2) pathogenic variants cause a neurodevelopmental syndrome, that is characterized by developmental delay, cognitive impairment, hypotonia, autistic features, behavior abnormalities and dysmorphic facial features. JARID2 encodes a transcriptional repressor protein that regulates the activity of various histone methyltransferase complexes. However, the molecular etiology is not fully understood, and JARID2-neurodevelopmental syndrome may vary in its typical clinical phenotype. In addition, the detection of variants of uncertain significance (VUSs) often results in a delay of final diagnosis which could hamper the appropriate care. In this study we aim to detect a specific and sensitive DNA methylation signature for JARID2-neurodevelopmental syndrome. Peripheral blood DNA methylation profiles from 56 control subjects, 8 patients with (likely) pathogenic JARID2 variants and 3 patients with JARID2 VUSs were analyzed. DNA methylation analysis indicated a clear and robust separation between patients with (likely) pathogenic variants and controls. A binary model capable of classifying patients with the JARID2-neurodevelopmental syndrome was constructed on the basis of the identified episignature. Patients carrying VUSs clustered with the control group. We identified a distinct DNA methylation signature associated with JARID2-neurodevelopmental syndrome, establishing its utility as a biomarker for this syndrome and expanding the EpiSign diagnostic test.

JARID2 haploinsufficiency is associated with a clinically distinct neurodevelopmental syndrome.

PURPOSE: JARID2, located on chromosome 6p22.3, is a regulator of histone methyltransferase complexes that is expressed in human neurons. So far, 13 individuals sharing clinical features including intellectual disability (ID) were reported with de novo heterozygous deletions in 6p22-p24 encompassing the full length JARID2 gene (OMIM 601594). However, all published individuals to date have a deletion of at least one other adjoining gene, making it difficult to determine if JARID2 is the critical gene responsible for the shared features. We aim to confirm JARID2 as a human disease gene and further elucidate the associated clinical phenotype. METHODS: Chromosome microarray analysis, exome sequencing, and an online matching platform (GeneMatcher) were used to identify individuals with single-nucleotide variants or deletions involving JARID2. RESULTS: We report 16 individuals in 15 families with a deletion or single-nucleotide variant in JARID2. Several of these variants are likely to result in haploinsufficiency due to nonsense-mediated messenger RNA (mRNA) decay. All individuals have developmental delay and/or ID and share some overlapping clinical characteristics such as facial features with those who have larger deletions involving JARID2. CONCLUSION: We report that JARID2 haploinsufficiency leads to a clinically distinct neurodevelopmental syndrome, thus establishing gene-disease validity for the purpose of diagnostic reporting.

Replicative and non-replicative mechanisms in the formation of clustered CNVs are indicated by whole genome characterization.

Clustered copy number variants (CNVs) as detected by chromosomal microarray analysis (CMA) are often reported as germline chromothripsis. However, such cases might need further investigations by massive parallel whole genome sequencing (WGS) in order to accurately define the underlying complex rearrangement, predict the occurrence mechanisms and identify additional complexities. Here, we utilized WGS to delineate the rearrangement structure of 21 clustered CNV carriers first investigated by CMA and identified a total of 83 breakpoint junctions (BPJs). The rearrangements were further sub-classified depending on the patterns observed: I) Cases with only deletions (n = 8) often had additional structural rearrangements, such as insertions and inversions typical to chromothripsis; II) cases with only duplications (n = 7) or III) combinations of deletions and duplications (n = 6) demonstrated mostly interspersed duplications and BPJs enriched with microhomology. In two cases the rearrangement mutational signatures indicated both a breakage-fusion-bridge cycle process and haltered formation of a ring chromosome. Finally, we observed two cases with Alu- and LINE-mediated rearrangements as well as two unrelated individuals with seemingly identical clustered CNVs on 2p25.3, possibly a rare European founder rearrangement. In conclusion, through detailed characterization of the derivative chromosomes we show that multiple mechanisms are likely involved in the formation of clustered CNVs and add further evidence for chromoanagenesis mechanisms in both "simple" and highly complex chromosomal rearrangements. Finally, WGS characterization adds positional information, important for a correct clinical interpretation and deciphering mechanisms involved in the formation of these rearrangements.

Small 6q16.1 Deletions Encompassing POU3F2 Cause Susceptibility to Obesity and Variable Developmental Delay with Intellectual Disability.

Genetic studies of intellectual disability and identification of monogenic causes of obesity in humans have made immense contribution toward the understanding of the brain and control of body mass. The leptin > melanocortin > SIM1 pathway is dysregulated in multiple monogenic human obesity syndromes but its downstream targets are still unknown. In ten individuals from six families, with overlapping 6q16.1 deletions, we describe a disorder of variable developmental delay, intellectual disability, and susceptibility to obesity and hyperphagia. The 6q16.1 deletions segregated with the phenotype in multiplex families and were shown to be de novo in four families, and there was dramatic phenotypic overlap among affected individuals who were independently ascertained without bias from clinical features. Analysis of the deletions revealed a ∼350 kb critical region on chromosome 6q16.1 that encompasses a gene for proneuronal transcription factor POU3F2, which is important for hypothalamic development and function. Using morpholino and mutant zebrafish models, we show that POU3F2 lies downstream of SIM1 and controls oxytocin expression in the hypothalamic neuroendocrine preoptic area. We show that this finding is consistent with the expression patterns of POU3F2 and related genes in the human brain. Our work helps to further delineate the neuro-endocrine control of energy balance/body mass and demonstrates that this molecular pathway is conserved across multiple species.

Prognosis of Preschool Eczema and Factors of Importance for Remission.

Information on factors of importance for remission of eczema is scarce. This study explored factors related to the remission and course of preschool eczema (PSE) (eczema at 1, 2 and/or 4 years of age) to 16 years of age (n = 889) in a Swedish cohort. Half of the children were in complete remission by school age (at age 8, 12, and 16 years). In multivariate prognostic models, persistent PSE (eczema at 1, 2 and 4 years of age) (odds ratio 0.27 (95% confidence interval 0.18-0.41)), PSE with sleep disturbance (due to itch at least once a week at 1, 2 and/or 4 years of age) (0.59 (0.43-0.81)), parental allergy (0.73 (0.55-0.96)), parental smoking at child's birth (0.70 (0.50-0.99)) and filaggrin mutation (R501X, R2447X, 2282del4) (0.47 (0.26-0.85)) were inversely associated with complete remission by school age. Male sex (1.37 (1.03-1.82)) and exclusive breastfeeding ≥4 months (1.44 (1.01-2.05)) were positively associated with complete remission by school age. In conclusion, half of the children with PSE were in complete remission by school age. The most important prognostic factors were persistent PSE and PSE with sleep disturbance due to itch.

A genome-wide association study of atopic dermatitis identifies loci with overlapping effects on asthma and psoriasis.

Atopic dermatitis (AD) is the most common dermatological disease of childhood. Many children with AD have asthma and AD shares regions of genetic linkage with psoriasis, another chronic inflammatory skin disease. We present here a genome-wide association study (GWAS) of childhood-onset AD in 1563 European cases with known asthma status and 4054 European controls. Using Illumina genotyping followed by imputation, we generated 268 034 consensus genotypes and in excess of 2 million single nucleotide polymorphisms (SNPs) for analysis. Association signals were assessed for replication in a second panel of 2286 European cases and 3160 European controls. Four loci achieved genome-wide significance for AD and replicated consistently across all cohorts. These included the epidermal differentiation complex (EDC) on chromosome 1, the genomic region proximal to LRRC32 on chromosome 11, the RAD50/IL13 locus on chromosome 5 and the major histocompatibility complex (MHC) on chromosome 6; reflecting action of classical HLA alleles. We observed variation in the contribution towards co-morbid asthma for these regions of association. We further explored the genetic relationship between AD, asthma and psoriasis by examining previously identified susceptibility SNPs for these diseases. We found considerable overlap between AD and psoriasis together with variable coincidence between allergic rhinitis (AR) and asthma. Our results indicate that the pathogenesis of AD incorporates immune and epidermal barrier defects with combinations of specific and overlapping effects at individual loci.

Rare copy number variants are common in young children with autism spectrum disorder.

AIM: Several studies have suggested that rare copy number variants (CNVs) are an important genetic contributor to autism spectrum disorders. The aims of the study were to use chromosomal microarray to investigate the presence of rare copy number variants in a population-based cohort of well-characterised young children with autism spectrum disorders and to relate the genetic results to neurodevelopmental profiles and medical conditions. METHODS: We performed chromosomal microarray on samples from 162 children who had been referred to the Stockholm Autism Centre for Young Children in Sweden after being diagnosed with autism spectrum disorder between 20 and 54 months of age. RESULTS: Pathogenic aberrations were detected in 8.6% of the children and variants of uncertain significance were present in another 8.6%. CNVs were more frequent in children with congenital malformations or dysmorphic features as well as in the subgroup with intellectual disability. CONCLUSION: Our results support the use of chromosomal microarray methods for the first tier genetic analysis of autism spectrum disorder. However, it is likely in the near future that chromosomal microarray methods will probably be replaced by whole-exome and whole-genome sequencing technologies in clinical genetic testing.

Intragenic duplication--a novel causative mechanism for SATB2-associated syndrome.

Previous studies have shown that genetic aberrations involving the special AT-rich sequence-binding protein 2 (SATB2) gene result in a variable phenotype of syndromic intellectual disability. Although only a small number of patients have been described, there is already considerable variation in regard to the underlying molecular mechanism spanning from structural variation to point mutations. We here describe a male patient with intellectual disability, speech and language impairment, cleft palate, malformed teeth, and oligodontia. Array CGH analysis identified a small intragenic duplication in the SATB2 gene that included three coding exons. The result was confirmed by multiplex ligation-dependent probe amplification and low coverage whole genome mate pair sequencing. WGS breakpoint analysis directly confirmed the duplication as intragenic. This is the first reported patient with an intragenic duplication in SATB2 in combination with a phenotype that is highly similar to previously described patients with small deletions or point mutations of the same gene. Our findings expand the spectra of SATB2 mutations and confirm the presence of a distinct SATB2-phenotype with severe ID and speech impairment, cleft palate and/or high arched palate, and abnormalities of the teeth. For patients that present with this clinical picture, a high-resolution exon targeted array CGH and/or WGS, in addition to sequencing of SATB2, should be considered.

Molecular and cytogenetic analysis in stillbirth: results from 481 consecutive cases.

INTRODUCTION: The underlying causes of stillbirth are heterogeneous and in many cases unexplained. Our aim was to conclude clinical results from karyotype and quantitative fluorescence-polymerase chain reaction (QF-PCR) analysis of all stillbirths occurring in Stockholm County between 2008 and 2012. By screening a subset of cases, we aimed to study the possible benefits of chromosomal microarray (CMA) in the analysis of the etiology of stillbirth. METHODS: During 2008-2012, 481 stillbirths in Stockholm County were investigated according to a clinical protocol including karyotype or QF-PCR analysis. CMA screening was performed on a subset of 90 cases, corresponding to all stillbirths from 2010 without a genetic diagnosis. RESULTS: Chromosomal aberrations were detected by karyotype or QF-PCR analysis in 7.5% of the stillbirths. CMA analysis additionally identified two known syndromes, one aberration disrupting a known disease gene, and 26 variants of unknown significance. Furthermore, CMA had a significantly higher success rate than karyotyping (100 vs. 80%, p < 0.001). DISCUSSION: In the analysis of stillbirth, conventional karyotyping is prone to failure, and QF-PCR is a useful complement. We show that CMA has a higher success rate and aberration detection frequency than these methods, and conclude that CMA is a valuable tool for identification of chromosomal aberrations in stillbirth.

Partial tetrasomy 14 associated with multiple malformations.

We report on an 8-year-old female patient with multiple malformations including bilateral cleft lip and palate, coloboma, and craniosynostosis. She presented with severe intellectual disability, seizures, and gastrointestinal dysfunction. Mitochondrial investigations in a muscle biopsy revealed reduced activity in complex I of the mitochondrial respiratory chain. Chromosome analysis and fluorescent in situ hybridization (FISH) studies showed an isodicentric marker chromosome 14 that was identified in all cells analyzed in peripheral blood lymphocytes and cultured fibroblasts. Parental chromosome studies were normal. To further characterize the marker chromosome and determine its origin, we performed array-based comparative genomic hybridization (CGH) and polymorphic marker analysis with quantitative fluorescent PCR (QF-PCR). The combined results from cytogenetic and array-CGH analyses showed tetrasomy 14p13q13.1 and results from the QF-PCR point to formation of the marker chromosome in the maternal meiosis. Isodicentric chromosomes involving partial 14q have previously been reported in four cases; however, this is the first patient with tetrasomy 14p13q13.1 in non-mosaic form surviving beyond infancy.

Cell-free tumour DNA analysis detects copy number alterations in gastro-oesophageal cancer patients.

BACKGROUND: Analysis of cell-free tumour DNA, a liquid biopsy, is a promising biomarker for cancer. We have performed a proof-of principle study to test the applicability in the clinical setting, analysing copy number alterations (CNAs) in plasma and tumour tissue from 44 patients with gastro-oesophageal cancer. METHODS: DNA was isolated from blood plasma and a tissue sample from each patient. Array-CGH was applied to the tissue DNA. The cell-free plasma DNA was sequenced by low-coverage whole-genome sequencing using a clinical pipeline for non-invasive prenatal testing. WISECONDOR and ichorCNA, two bioinformatic tools, were used to process the output data and were compared to each other. RESULTS: Cancer-associated CNAs could be seen in 59% (26/44) of the tissue biopsies. In the plasma samples, a targeted approach analysing 61 regions of special interest in gastro-oesophageal cancer detected cancer-associated CNAs with a z-score >5 in 11 patients. Broadening the analysis to a whole-genome view, 17/44 patients (39%) had cancer-associated CNAs using WISECONDOR and 13 (30%) using ichorCNA. Of the 26 patients with tissue-verified cancer-associated CNAs, 14 (54%) had corresponding CNAs in plasma. Potentially clinically actionable amplifications overlapping the genes VEGFA, EGFR and FGFR2 were detected in the plasma from three patients. CONCLUSIONS: We conclude that low-coverage whole-genome sequencing without prior knowledge of the tumour alterations could become a useful tool for cell-free tumour DNA analysis of total CNAs in plasma from patients with gastro-oesophageal cancer.

Identification of putative pathogenic single nucleotide variants (SNVs) in genes associated with heart disease in 290 cases of stillbirth.

The incidence of stillbirth in Sweden has essentially remained constant since the 1980's, and despite thorough investigation, many cases remain unexplained. It has been suggested that a proportion of stillbirth cases is caused by heart disease, mainly channelopathies. The aim of this study was to analyze DNA from 290 stillbirth cases without chromosomal abnormalities for pathogenic single nucleotide variants (SNVs) in 70 genes associated with cardiac channelopathies and cardiomyopathies. The HaloPlex Target Enrichment System (Agilent Technologies) was utilized to prepare sequencing libraries which were sequenced on the Illumina NextSeq platform. We found that 12.1% of the 290 investigated stillbirth cases had one (n = 31) or two (n = 4) variants with evidence supporting pathogenicity, i.e. loss-of-function variants (nonsense, frameshift, splice site substitutions), evidence from functional studies, or previous identification of the variants in affected individuals. Regarding identified putative pathogenic variants in genes associated with channelopathies, the prevalence was significantly higher in the stillbirth cohort (n = 23, 7.93%) than the corresponding prevalence of the same variants in the non-Finnish European population of the Exome Aggregation Consortium (2.70%, p<0.001) and SweGen, (2.30%, p<0.001). Our results give further support to the hypothesis that cardiac channelopathies might contribute to stillbirth. Screening for pathogenic SNVs in genes associated with heart disease might be a valuable complement for stillbirth cases where today's conventional investigation does not reveal the underlying cause of fetal demise.

Ataxia in Patients With Bi-Allelic NFASC Mutations and Absence of Full-Length NF186.

The etiology of hereditary ataxia syndromes is heterogeneous, and the mechanisms underlying these disorders are often unknown. Here, we utilized exome sequencing in two siblings with progressive ataxia and muscular weakness and identified a novel homozygous splice mutation (c.3020-1G > A) in neurofascin (NFASC). In RNA extracted from fibroblasts, we showed that the mutation resulted in inframe skipping of exon 26, with a deprived expression of the full-length transcript that corresponds to NFASC isoform NF186. To further investigate the disease mechanisms, we reprogrammed fibroblasts from one affected sibling to induced pluripotent stem cells, directed them to neuroepithelial stem cells and finally differentiated to neurons. In early neurogenesis, differentiating cells with selective depletion of the NF186 isoform showed significantly reduced neurite outgrowth as well as fewer emerging neurites. Furthermore, whole-cell patch-clamp recordings of patient-derived neuronal cells revealed a lower threshold for openings, indicating altered Na+ channel kinetics, suggesting a lower threshold for openings as compared to neuronal cells without the NFASC mutation. Taken together, our results suggest that loss of the full-length NFASC isoform NF186 causes perturbed neurogenesis and impaired neuronal biophysical properties resulting in a novel early-onset autosomal recessive ataxia syndrome.

From cytogenetics to cytogenomics: whole-genome sequencing as a first-line test comprehensively captures the diverse spectrum of disease-causing genetic variation underlying intellectual disability.

BACKGROUND: Since different types of genetic variants, from single nucleotide variants (SNVs) to large chromosomal rearrangements, underlie intellectual disability, we evaluated the use of whole-genome sequencing (WGS) rather than chromosomal microarray analysis (CMA) as a first-line genetic diagnostic test. METHODS: We analyzed three cohorts with short-read WGS: (i) a retrospective cohort with validated copy number variants (CNVs) (cohort 1, n = 68), (ii) individuals referred for monogenic multi-gene panels (cohort 2, n = 156), and (iii) 100 prospective, consecutive cases referred to our center for CMA (cohort 3). Bioinformatic tools developed include FindSV, SVDB, Rhocall, Rhoviz, and vcf2cytosure. RESULTS: First, we validated our structural variant (SV)-calling pipeline on cohort 1, consisting of three trisomies and 79 deletions and duplications with a median size of 850 kb (min 500 bp, max 155 Mb). All variants were detected. Second, we utilized the same pipeline in cohort 2 and analyzed with monogenic WGS panels, increasing the diagnostic yield to 8%. Next, cohort 3 was analyzed by both CMA and WGS. The WGS data was processed for large (> 10 kb) SVs genome-wide and for exonic SVs and SNVs in a panel of 887 genes linked to intellectual disability as well as genes matched to patient-specific Human Phenotype Ontology (HPO) phenotypes. This yielded a total of 25 pathogenic variants (SNVs or SVs), of which 12 were detected by CMA as well. We also applied short tandem repeat (STR) expansion detection and discovered one pathologic expansion in ATXN7. Finally, a case of Prader-Willi syndrome with uniparental disomy (UPD) was validated in the WGS data. Important positional information was obtained in all cohorts. Remarkably, 7% of the analyzed cases harbored complex structural variants, as exemplified by a ring chromosome and two duplications found to be an insertional translocation and part of a cryptic unbalanced translocation, respectively. CONCLUSION: The overall diagnostic rate of 27% was more than doubled compared to clinical microarray (12%). Using WGS, we detected a wide range of SVs with high accuracy. Since the WGS data also allowed for analysis of SNVs, UPD, and STRs, it represents a powerful comprehensive genetic test in a clinical diagnostic laboratory setting.

Solute carriers (SLC) in inflammatory bowel disease: a potential target of probiotics?

Transporter proteins of the solute carriers (SLCs) family play a role in epithelial permeability and barrier function in the intestine, and polymorphisms in SLC genes are associated with inflammatory bowel disease. Many SLCs also mediate the bioavailability of pharmaceutical compounds, and the modulation of such transport systems to increase drug efficacy is, therefore, of great interest. We have undertaken a large-scale project to evaluate whether bacteria can modulate the expression of SLCs in the intestine. Here we report the effect of VSL[sharp]3 (a high-potency probiotic preparation) on the expression of 3 large solute carrier families, SLC4, SLC21, and SLC22, which are involved in the transport of bicarbonates, organic anions and cations, and affect the bioavailability of several pharmaceutical compounds. Two groups of animals (VSL[sharp]3 and phosphate-buffered saline controls) were studied for SLC expression in the intestine by Real-Time PCR at the beginning (day 1) and at the end (day 20) of the treatment, and 7 days after the interruption of the treatment. An effect of VSL[sharp]3 administration was detected on the expression of 10% of the studied genes. This reached statistical significance (P=0.01) for the poorly characterized sodium-borate cotransporter SLC4A11, which showed a 5-times lower expression in VSL[sharp]3 than in control mice on day 1 of probiotic treatment. VSL[sharp]3-driven changes in the expression levels of SLC transporters might contribute to its reported effects on intestinal permeability. The elucidation of SLC4A11 function in the intestine will be the key to fully evaluate the relevance of specific findings.

TSHZ3 deletion causes an autism syndrome and defects in cortical projection neurons.

TSHZ3, which encodes a zinc-finger transcription factor, was recently positioned as a hub gene in a module of the genes with the highest expression in the developing human neocortex, but its functions remained unknown. Here we identify TSHZ3 as the critical region for a syndrome associated with heterozygous deletions at 19q12-q13.11, which includes autism spectrum disorder (ASD). In Tshz3-null mice, differentially expressed genes include layer-specific markers of cerebral cortical projection neurons (CPNs), and the human orthologs of these genes are strongly associated with ASD. Furthermore, mice heterozygous for Tshz3 show functional changes at synapses established by CPNs and exhibit core ASD-like behavioral abnormalities. These findings highlight essential roles for Tshz3 in CPN development and function, whose alterations can account for ASD in the newly defined TSHZ3 deletion syndrome.

Deletion of Wiskott-Aldrich syndrome protein triggers Rac2 activity and increased cross-presentation by dendritic cells.

Wiskott-Aldrich syndrome (WAS) is caused by loss-of-function mutations in the WASp gene. Decreased cellular responses in WASp-deficient cells have been interpreted to mean that WASp directly regulates these responses in WASp-sufficient cells. Here, we identify an exception to this concept and show that WASp-deficient dendritic cells have increased activation of Rac2 that support cross-presentation to CD8(+) T cells. Using two different skin pathology models, WASp-deficient mice show an accumulation of dendritic cells in the skin and increased expansion of IFNγ-producing CD8(+) T cells in the draining lymph node and spleen. Specific deletion of WASp in dendritic cells leads to marked expansion of CD8(+) T cells at the expense of CD4(+) T cells. WASp-deficient dendritic cells induce increased cross-presentation to CD8(+) T cells by activating Rac2 that maintains a near neutral pH of phagosomes. Our data reveals an intricate balance between activation of WASp and Rac2 signalling pathways in dendritic cells.

No evidence for mosaic pathogenic copy number variations in cardiac tissue from patients with congenital heart malformations.

The aim of this study was to investigate if pathogenic copy number variations (CNVs) are present in mosaic form in patients with congenital heart malformations. We have collected cardiac tissue and blood samples from 23 patients with congenital heart malformations that underwent cardiac surgery and screened for mosaic gene dose alterations restricted to cardiac tissue using array comparative genomic hybridization (array CGH). We did not find evidence of CNVs in mosaic form after array CGH analysis. Pathogenic CNVs that were present in both cardiac tissue and blood were detected in 2/23 patients (9%), and in addition we found several constitutional CNVs of unclear clinical significance. This is the first study investigating mosaicism for CNVs in heart tissue compared to peripheral blood and the results do not indicate that pathogenic mosaic copy number changes are common in patients with heart malformations. Importantly, in line with previous studies, our results show that constitutional pathogenic CNVs are important factors contributing to congenital heart malformations.