Assortative mating and parental genetic relatedness contribute to the pathogenicity of variably expressive variants.

We examined more than 97,000 families from four neurodevelopmental disease cohorts and the UK Biobank to identify phenotypic and genetic patterns in parents contributing to neurodevelopmental disease risk in children. We identified within- and cross-disorder correlations between six phenotypes in parents and children, such as obsessive-compulsive disorder (R = 0.32-0.38, p < 10-126). We also found that measures of sub-clinical autism features in parents are associated with several autism severity measures in children, including biparental mean Social Responsiveness Scale scores and proband Repetitive Behaviors Scale scores (regression coefficient = 0.14, p = 3.38 × 10-4). We further describe patterns of phenotypic similarity between spouses, where spouses show correlations for six neurological and psychiatric phenotypes, including a within-disorder correlation for depression (R = 0.24-0.68, p < 0.001) and a cross-disorder correlation between anxiety and bipolar disorder (R = 0.09-0.22, p < 10-92). Using a simulated population, we also found that assortative mating can lead to increases in disease liability over generations and the appearance of "genetic anticipation" in families carrying rare variants. We identified several families in a neurodevelopmental disease cohort where the proband inherited multiple rare variants in disease-associated genes from each of their affected parents. We further identified parental relatedness as a risk factor for neurodevelopmental disorders through its inverse relationship with variant pathogenicity and propose that parental relatedness modulates disease risk by increasing genome-wide homozygosity in children (R = 0.05-0.26, p < 0.05). Our results highlight the utility of assessing parent phenotypes and genotypes toward predicting features in children who carry rare variably expressive variants and implicate assortative mating as a risk factor for increased disease severity in these families.

A Missense Variant Affecting the N-Terminal Domain of the Laminin-332 ß3 Chain Results in a Distinct Form of Junctional Epidermolysis Bullosa With Altered Granulation Tissue Response and No New Blistering: A Second Family Report.

Junctional epidermolysis bullosa (JEB) is a rare genodermatosis characterized by fragility of the skin and mucous membranes due to alterations in the dermal epidermal junction. This condition manifests as mechanically induced bullous lesions that heal with hypertrophic granulation tissue and/or atrophic scars. Here, we report two brothers carrying a homozygous LAMB3 missense variant, p.Gly254Asp, which affects the N-terminal end of the laminin-332 (LM332) ß3 chain, previously described in another JEB family sharing a common ethnic origin and LAMB3 haplotype with the siblings reported here. Moreover, all affected patients with p.Gly254Asp mutation from both families exhibits a distinct phenotype consisting of a few localized long-standing skin lesions characterized by excessive granulation tissue formation or keloid scars, without new blistering, and associated with amelogenesis imperfecta. Our patients also showed nail dystrophy, expanding the phenotypic spectrum and confirming the peculiar role of the N-terminal end of the ß3 chain in regulating the granulation tissue response associated with the wound healing process.

Exome sequencing efficacy and phenotypic expansions involving esophageal atresia/tracheoesophageal fistula plus.

Esophageal atresia/tracheoesophageal fistula (EA/TEF) is a life-threatening birth defect that often occurs with other major birth defects (EA/TEF+). Despite advances in genetic testing, a molecular diagnosis can only be made in a minority of EA/TEF+ cases. Here, we analyzed clinical exome sequencing data and data from the DECIPHER database to determine the efficacy of exome sequencing in cases of EA/TEF+ and to identify phenotypic expansions involving EA/TEF. Among 67 individuals with EA/TEF+ referred for clinical exome sequencing, a definitive or probable diagnosis was made in 11 cases for an efficacy rate of 16% (11/67). This efficacy rate is significantly lower than that reported for other major birth defects, suggesting that polygenic, multifactorial, epigenetic, and/or environmental factors may play a particularly important role in EA/TEF pathogenesis. Our cohort included individuals with pathogenic or likely pathogenic variants that affect TCF4 and its downstream target NRXN1, and FANCA, FANCB, and FANCC, which are associated with Fanconi anemia. These cases, previously published case reports, and comparisons to other EA/TEF genes made using a machine learning algorithm, provide evidence in support of a potential pathogenic role for these genes in the development of EA/TEF.

Mutational spectrum and clinical signatures in 114 families with hereditary multiple osteochondromas: insights into molecular properties of selected exostosin variants.

Hereditary multiple osteochondromas (HMO) is a rare autosomal dominant skeletal disorder, caused by heterozygous variants in either EXT1 or EXT2, which encode proteins involved in the biogenesis of heparan sulphate. Pathogenesis and genotype-phenotype correlations remain poorly understood. We studied 114 HMO families (158 affected individuals) with causative EXT1 or EXT2 variants identified by Sanger sequencing, or multiplex ligation-dependent probe amplification and qPCR. Eighty-seven disease-causative variants (55 novel and 32 known) were identified including frameshift (42%), nonsense (32%), missense (11%), splicing (10%) variants and genomic rearrangements (5%). Informative clinical features were available for 42 EXT1 and 27 EXT2 subjects. Osteochondromas were more frequent in EXT1 as compared to EXT2 patients. Anatomical distribution of lesions showed significant differences based on causative gene. Microscopy analysis for selected EXT1 and EXT2 variants verified that EXT1 and EXT2 mutants failed to co-localize each other and loss Golgi localization by surrounding the nucleus and/or assuming a diffuse intracellular distribution. In a cell viability study, cells expressing EXT1 and EXT2 mutants proliferated more slowly than cells expressing wild-type proteins. This confirms the physiological relevance of EXT1 and EXT2 Golgi co-localization and the key role of these proteins in the cell cycle. Taken together, our data expand genotype-phenotype correlations, offer further insights in the pathogenesis of HMO and open the path to future therapies.

When to test fetuses for RASopathies? Proposition from a systematic analysis of 352 multicenter cases and a postnatal cohort.

PURPOSE: Recent studies have identified suggestive prenatal features of RASopathies (e.g., increased nuchal translucency [NT], cystic hygroma [CH], hydrops, effusions, congenital heart diseases [CHD], polyhydramnios, renal anomalies). Our objective is to clarify indications for RASopathy prenatal testing. We compare genotype distributions between pre- and postnatal populations and propose genotype-phenotype correlations. METHODS: Three hundred fifty-two chromosomal microarray-negative cases sent for prenatal RASopathy testing between 2012 and 2019 were collected. For most, 11 RASopathy genes were tested. Postnatal cohorts (25 patients with available prenatal information and 108 institutional database genotypes) and the NSeuroNet database were used for genotypic comparisons. RESULTS: The overall diagnostic yield was 14% (50/352), with rates >20% for effusions, hydrops, and CHD. Diagnostic yield was significantly improved in presence of hypertrophic cardiomyopathy (HCM), persistent or associated CH, any suggestive finding combined with renal anomaly or polyhydramnios, or ≥2 ultrasound findings. Largest prenatal contributors of pathogenic variants were PTPN11 (30%), RIT1 (16%), RAF1 (14%), and HRAS (12%), which considerably differ from their prevalence in postnatal populations. HRAS, LZTR1, and RAF1 variants correlated with hydrops/effusions, and RIT1 with prenatal onset HCM. CONCLUSION: After normal chromosomal microarray, RASopathies should be considered when any ultrasound finding of lymphatic dysplasia or suggestive CHD is found alone or in association.

Copy number variation analysis implicates novel pathways in patients with oculo-auriculo-vertebral-spectrum and congenital heart defects.

Oculo-auriculo-vertebral spectrum (OAVS) is a developmental disorder of craniofacial morphogenesis. Its etiology is unclear, but assumed to be complex and heterogeneous, with contribution of both genetic and environmental factors. We assessed the occurrence of copy number variants (CNVs) in a cohort of 19 unrelated OAVS individuals with congenital heart defect. Chromosomal microarray analysis identified pathogenic CNVs in 2/19 (10.5%) individuals, and CNVs classified as variants of uncertain significance in 7/19 (36.9%) individuals. Remarkably, two subjects had small intragenic CNVs involving DACH1 and DACH2, two paralogs coding for key components of the PAX-SIX-EYA-DACH network, a transcriptional regulatory pathway controlling developmental processes relevant to OAVS and causally associated with syndromes characterized by craniofacial involvement. Moreover, a third patient showed a large duplication encompassing DMBX1/OTX3, encoding a transcriptional repressor of OTX2, another transcription factor functionally connected to the DACH-EYA-PAX network. Among the other relevant CNVs, a deletion encompassing HSD17B6, a gene connected with the retinoic acid signaling pathway, whose dysregulation has been implicated in craniofacial malformations, was also identified. Our findings suggest that CNVs affecting gene dosage likely contribute to the genetic heterogeneity of OAVS, and implicate the PAX-SIX-EYA-DACH network as novel pathway involved in the etiology of this developmental trait.

Genome-Wide DNA Methylation Analysis of a Cohort of 41 Patients Affected by Oculo-Auriculo-Vertebral Spectrum (OAVS).

Oculo-auriculo-vertebral-spectrum (OAVS; OMIM 164210) is a rare disorder originating from abnormal development of the first and second branchial arch. The clinical phenotype is extremely heterogeneous with ear anomalies, hemifacial microsomia, ocular defects, and vertebral malformations being the main features. MYT1, AMIGO2, and ZYG11B gene variants were reported in a few OAVS patients, but the etiology remains largely unknown. A multifactorial origin has been proposed, including the involvement of environmental and epigenetic mechanisms. To identify the epigenetic mechanisms contributing to OAVS, we evaluated the DNA-methylation profiles of 41 OAVS unrelated affected individuals by using a genome-wide microarray-based methylation approach. The analysis was first carried out comparing OAVS patients with controls at the group level. It revealed a moderate epigenetic variation in a large number of genes implicated in basic chromatin dynamics such as DNA packaging and protein-DNA organization. The alternative analysis in individual profiles based on the searching for Stochastic Epigenetic Variants (SEV) identified an increased number of SEVs in OAVS patients compared to controls. Although no recurrent deregulated enriched regions were found, isolated patients harboring suggestive epigenetic deregulations were identified. The recognition of a different DNA methylation pattern in the OAVS cohort and the identification of isolated patients with suggestive epigenetic variations provide consistent evidence for the contribution of epigenetic mechanisms to the etiology of this complex and heterogeneous disorder.

Expanding the phenotype of Wiedemann-Steiner syndrome: Craniovertebral junction anomalies.

Wiedemann-Steiner syndrome (WDSTS) is a rare autosomal dominant condition caused by heterozygous loss of function variants in the KMT2A (MLL) gene, encoding a lysine N-methyltransferase that mediates a histone methylation pattern specific for epigenetic transcriptional activation. WDSTS is characterized by a distinctive facial phenotype, hypertrichosis, short stature, developmental delay, intellectual disability, congenital malformations, and skeletal anomalies. Recently, a few patients have been reported having abnormal skeletal development of the cervical spine. Here we describe 11 such individuals, all with KMT2A de novo loss-of-function variants: 10 showed craniovertebral junction anomalies, while an 11th patient had a cervical abnormality in C7. By evaluating clinical and diagnostic imaging data we characterized these anomalies, which consist primarily of fused cervical vertebrae, C1 and C2 abnormalities, small foramen magnum and Chiari malformation type I. Craniovertebral anomalies in WDSTS patients have been largely disregarded so far, but the increasing number of reports suggests that they may be an intrinsic feature of this syndrome. Specific investigation strategies should be considered for early identification and prevention of craniovertebral junction complications in WDSTS patients.

A Customized Next-Generation Sequencing-Based Panel to Identify Novel Genetic Variants in Dementing Disorders: A Pilot Study.

Purpose: The advancements in the next-generation sequencing (NGS) techniques have allowed for rapid, efficient, and cost-time-effective genetic variant detection. However, in both clinical practice and research setting, sequencing is still often limited to the use of gene panels clinically targeted on the genes underlying the disease of interest. Methods: We performed a neurogenetic study through an ad hoc NGS-based custom sequencing gene panel in order to screen 16 genes in 8 patients with different types of degenerative cognitive disorders (Alzheimer's disease, mild cognitive impairment, frontotemporal dementia, and dementia associated with Parkinson's disease). The study protocol was based on previous evidence showing a high sensitivity and specificity of the technique even when the panel is limited to some hotspot exons. Results: We found variants of the TREM2 and APP genes in three patients; these have been previously identified as pathogenic or likely pathogenic and, therefore, considered "disease causing." In the remaining subjects, the pathogenicity was evaluated according to the guidelines of the American College of Medical Genetics (ACMG). In one patient, the p.R205W variant in the CHMP2B gene was found to be likely pathogenic of the disease. A variant in the CSF1R and SERPINI1 genes found in two patients was classified as benign, whereas the other two (in the GRN and APP genes) were classified as likely pathogenic according to the ACMG. Conclusions: Notwithstanding the preliminary value of this study, some rare genetic variants with a probable disease association were detected. Although future application of NGS-based sensors and further replication of these experimental data are needed, this approach seems to offer promising translational perspectives in the diagnosis and management of a wide range of neurodegenerative disorders.

Small supernumerary marker chromosomes: A legacy of trisomy rescue?

We studied by a whole genomic approach and trios genotyping, 12 de novo, nonrecurrent small supernumerary marker chromosomes (sSMC), detected as mosaics during pre- or postnatal diagnosis and associated with increased maternal age. Four sSMCs contained pericentromeric portions only, whereas eight had additional non-contiguous portions of the same chromosome, assembled together in a disordered fashion by repair-based mechanisms in a chromothriptic event. Maternal hetero/isodisomy was detected with a paternal origin of the sSMC in some cases, whereas in others two maternal alleles in the sSMC region and biparental haplotypes of the homologs were detected. In other cases, the homologs were biparental while the sSMC had the same haplotype of the maternally inherited chromosome. These findings strongly suggest that most sSMCs are the result of a multiple-step mechanism, initiated by maternal meiotic nondisjunction followed by postzygotic anaphase lagging of the supernumerary chromosome and its subsequent chromothripsis.

Rare variants in the genetic background modulate cognitive and developmental phenotypes in individuals carrying disease-associated variants.

PURPOSE: To assess the contribution of rare variants in the genetic background toward variability of neurodevelopmental phenotypes in individuals with rare copy-number variants (CNVs) and gene-disruptive variants. METHODS: We analyzed quantitative clinical information, exome sequencing, and microarray data from 757 probands and 233 parents and siblings who carry disease-associated variants. RESULTS: The number of rare likely deleterious variants in functionally intolerant genes ("other hits") correlated with expression of neurodevelopmental phenotypes in probands with 16p12.1 deletion (n=23, p=0.004) and in autism probands carrying gene-disruptive variants (n=184, p=0.03) compared with their carrier family members. Probands with 16p12.1 deletion and a strong family history presented more severe clinical features (p=0.04) and higher burden of other hits compared with those with mild/no family history (p=0.001). The number of other hits also correlated with severity of cognitive impairment in probands carrying pathogenic CNVs (n=53) or de novo pathogenic variants in disease genes (n=290), and negatively correlated with head size among 80 probands with 16p11.2 deletion. These co-occurring hits involved known disease-associated genes such as SETD5, AUTS2, and NRXN1, and were enriched for cellular and developmental processes. CONCLUSION: Accurate genetic diagnosis of complex disorders will require complete evaluation of the genetic background even after a candidate disease-associated variant is identified.

Melorheostosis and Osteopoikilosis Clinical and Molecular Description of an Italian Case Series.

Melorheostosis (MEL) is an uncommon, sclerosing disease, characterised by hyperostosis of long bones, resembling the flowing of candle wax. The disease is sporadic and the pathogenesis is still poorly understood. Occasionally, the same family can include individuals with MEL and Osteopoikilosis (OPK), a disease characterised by multiple round foci of increased bone density. LEMD3 gene mutations are related to OPK and Buschke-Ollendorff Syndrome, a genetic condition in which an association between MEL, OPK and skin lesions is observed. In rare cases, LEMD3 mutations and recently mosaic MAP2K1 gene mutations have been correlated to MEL suggesting that somatic mosaicism could be causative of the disease. In this study, we described the clinical, radiological and molecular findings of 19 individuals with MEL and 8 with OPK and compared the results to the medical literature. The molecular analyses of this case series corroborate the available data in the medical literature, indicating that LEMD3 germline mutations are not a major cause of isolated MEL and reporting five further cases of OPK caused by LEMD3 germline mutations.

Gene-targeted deletion in mice of the Ets-1 transcription factor, a candidate gene in the Jacobsen syndrome kidney "critical region," causes abnormal kidney development.

Ets-1 is a member of the Ets family of transcription factors and has critical roles in multiple biological functions. Structural kidney defects occur at an increased frequency in Jacobsen syndrome (OMIM #147791), a rare chromosomal disorder caused by deletions in distal 11q, implicating at least one causal gene in distal 11q. In this study, we define an 8.1 Mb "critical region" for kidney defects in Jacobsen syndrome, which spans ~50 genes. We demonstrate that gene-targeted deletion of Ets-1 in mice results in some of the most common congenital kidney defects occurring in Jacobsen syndrome, including: duplicated kidney, hypoplastic kidney, and dilated renal pelvis and calyces. Taken together, our results implicate Ets-1 in normal mammalian kidney development and, potentially, in the pathogenesis of some of the most common types of human structural kidney defects.

Small 4p16.3 deletions: Three additional patients and review of the literature.

Wolf-Hirschhorn syndrome is a well-defined disorder due to 4p16.3 deletion, characterized by distinct facial features, intellectual disability, prenatal and postnatal growth retardation, and seizures. Genotype-phenotype correlations based on differently sized deletions have been attempted, and some candidate genes have been suggested. We report on clinical characteristics of three patients with pure interstitial submicroscopic 4p16.3 deletions, ranging in size from 68 to 166 kb, involving WHSCR1 and/or part of WHSCR2, and review published cases with overlapping 4p16.3 losses. The present study highlights a major role of NSD2 gene in the pathogenesis of the WHS main features and predicts that loss-of-function mutations affecting NSD2 gene could result in microcephaly, prenatal and postnatal growth retardation, psychomotor and language delay, and craniofacial features. Absent seizures in all subjects corroborate the suggestion that this specific feature is causally linked with at least one additional causative gene. Finally, we suggest that mir-943 could play a role in the pathogenesis of CHD in some of these patients.

Familial 18q12.2 deletion supports the role of RNA-binding protein CELF4 in autism spectrum disorders.

Deletion of 18q12.2 is an increasingly recognized condition with a distinct neuropsychiatric phenotype. Twenty-two patients have been described with overlapping neurobehavioral disturbances including developmental delay, intellectual disability of variable degree, seizures, motor coordination disorder, behavioral/emotional disturbances, and autism spectrum disorders. The CUGBP Elav-like family member 4 (CELF4) gene at 18q12.2 encodes a RNA-binding protein that links to RNA subsets involved in pre- and postsynaptic neurotransmission including almost 30% of potential autism-related genes. Haploinsufficiency of CELF4 was associated with an autism or autistic behavior diagnosis in two adult patients with de novo 18q12.2 deletions. We report on a girl and her mildly affected mother with a 275 kb deletion at 18q12.2 involving CELF4 and KIAA1328 whose disruption is not associated with any known disease. The child was diagnosed with syndromic intellectual disability and autism at 6 years of age. Her mother had minor dysmorphisms, mild intellectual disability, and autistic behavior. The deleted region reported in this family is one of the smallest so far reported at 18q12.2. This is also the first full clinical description of maternally inherited CELF4 haploinsufficiency. The present study refines the molecular and neuropsychiatric phenotype associated with 18q12.2 deletion leading to CELF4 haploinsufficiency and provides evidence for a role for CELF4 in brain development and autism spectrum disorders.

Opposite effects on facial morphology due to gene dosage sensitivity.

Sequencing technology is increasingly demonstrating the impact of genomic copy number variation (CNV) on phenotypes. Opposing variation in growth, head size, cognition and behaviour is known to result from deletions and reciprocal duplications of some genomic regions. We propose normative inversion of face shape, opposing difference from a matched norm, as a basis for investigating the effects of gene dosage on craniofacial development. We use dense surface modelling techniques to match any face (or part of a face) to a facial norm of unaffected individuals of matched age, sex and ethnicity and then we reverse the individual's face shape differences from the matched norm to produce the normative inversion. We demonstrate for five genomic regions, 4p16.3, 7q11.23, 11p15, 16p13.3 and 17p11.2, that such inversion for individuals with a duplication or (epi)-mutation produces facial forms remarkably similar to those associated with a deletion or opposite (epi-)mutation of the same region, and vice versa. The ability to visualise and quantify face shape effects of gene dosage is of major benefit for determining whether a CNV is the cause of the phenotype of an individual and for predicting reciprocal consequences. It enables face shape to be used as a relatively simple and inexpensive functional analysis of the gene(s) involved.

TBR1 is the candidate gene for intellectual disability in patients with a 2q24.2 interstitial deletion.

Interstitial deletion of 2q24.2 is a rarely described cytogenetic aberration in patients with intellectual disability (ID). Previously reported genotype-phenotype correlation identified a minimum deleted region of 2.65 Mb including 15 genes. Recently, a patient with a de novo 2q24.2 microdeletion of 0.4 Mb encompassing only three genes was described. However, the precise relationship between most deleted genes and the clinical features remains unclear. Here we describe a 12-year-old male patient diagnosed with growth retardation and ID. He also showed microcephaly, right palpebral ptosis, scapular winging, and pectus excavatum. Single nucleotide polymorphisms (SNP) array analysis showed a de novo interstitial deletion of 0.122 Mb at 2q24.2 region harboring only TBR1 (T-box, brain, 1; OMIM: 604616), which encodes a T-box family transcription factor expressed in post-mitotic projection neurons and functionally significant in embryologic corticogenesis. This is the first case of a deletion at 2q24.2 involving only TBR1. This finding narrows the smallest region of overlap (SRO) for deletions in this region and strengthens the previously suggested hypothesis that this gene is a strong candidate for the ID phenotype. The identification of TBR1 as candidate for ID encourages further molecular studies to identify novel mutations to understand the pathogenic effects of its haploinsufficiency. Finally, this report provides a review on 10 2q24.2 microdeletion patients.

A corpus of GA4GH phenopackets: Case-level phenotyping for genomic diagnostics and discovery.

The Global Alliance for Genomics and Health (GA4GH) Phenopacket Schema was released in 2022 and approved by ISO as a standard for sharing clinical and genomic information about an individual, including phenotypic descriptions, numerical measurements, genetic information, diagnoses, and treatments. A phenopacket can be used as an input file for software that supports phenotype-driven genomic diagnostics and for algorithms that facilitate patient classification and stratification for identifying new diseases and treatments. There has been a great need for a collection of phenopackets to test software pipelines and algorithms. Here, we present Phenopacket Store. Phenopacket Store v.0.1.19 includes 6,668 phenopackets representing 475 Mendelian and chromosomal diseases associated with 423 genes and 3,834 unique pathogenic alleles curated from 959 different publications. This represents the first large-scale collection of case-level, standardized phenotypic information derived from case reports in the literature with detailed descriptions of the clinical data and will be useful for many purposes, including the development and testing of software for prioritizing genes and diseases in diagnostic genomics, machine learning analysis of clinical phenotype data, patient stratification, and genotype-phenotype correlations. This corpus also provides best-practice examples for curating literature-derived data using the GA4GH Phenopacket Schema.

A corpus of GA4GH Phenopackets: case-level phenotyping for genomic diagnostics and discovery.

The Global Alliance for Genomics and Health (GA4GH) Phenopacket Schema was released in 2022 and approved by ISO as a standard for sharing clinical and genomic information about an individual, including phenotypic descriptions, numerical measurements, genetic information, diagnoses, and treatments. A phenopacket can be used as an input file for software that supports phenotype-driven genomic diagnostics and for algorithms that facilitate patient classification and stratification for identifying new diseases and treatments. There has been a great need for a collection of phenopackets to test software pipelines and algorithms. Here, we present phenopacket-store. Version 0.1.12 of phenopacket-store includes 4916 phenopackets representing 277 Mendelian and chromosomal diseases associated with 236 genes, and 2872 unique pathogenic alleles curated from 605 different publications. This represents the first large-scale collection of case-level, standardized phenotypic information derived from case reports in the literature with detailed descriptions of the clinical data and will be useful for many purposes, including the development and testing of software for prioritizing genes and diseases in diagnostic genomics, machine learning analysis of clinical phenotype data, patient stratification, and genotype-phenotype correlations. This corpus also provides best-practice examples for curating literature-derived data using the GA4GH Phenopacket Schema.

Genetic modifiers and ascertainment drive variable expressivity of complex disorders.

Variable expressivity of disease-associated variants implies a role for secondary variants that modify clinical features. We assessed the effects of modifier variants towards clinical outcomes of 2,252 individuals with primary variants. Among 132 families with the 16p12.1 deletion, distinct rare and common variant classes conferred risk for specific developmental features, including short tandem repeats for neurological defects and SNVs for microcephaly, while additional disease-associated variants conferred multiple genetic diagnoses. Within disease and population cohorts of 773 individuals with the 16p12.1 deletion, we found opposing effects of secondary variants towards clinical features across ascertainments. Additional analysis of 1,479 probands with other primary variants, such as 16p11.2 deletion and CHD8 variants, and 1,084 without primary variants, showed that phenotypic associations differed by primary variant context and were influenced by synergistic interactions between primary and secondary variants. Our study provides a paradigm to dissect the genomic architecture of complex disorders towards personalized treatment.