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.

Chromatin conformation capture in the clinic: 4C-seq/HiC distinguishes pathogenic from neutral duplications at the GPR101 locus.

BACKGROUND: X-linked acrogigantism (X-LAG; MIM: 300942) is a severe form of pituitary gigantism caused by chromosome Xq26.3 duplications involving GPR101. X-LAG-associated duplications disrupt the integrity of the topologically associating domain (TAD) containing GPR101 and lead to the formation of a neo-TAD that drives pituitary GPR101 misexpression and gigantism. As X-LAG is fully penetrant and heritable, duplications involving GPR101 identified on prenatal screening studies, like amniocentesis, can pose an interpretation challenge for medical geneticists and raise important concerns for patients and families. Therefore, providing robust information on the functional genomic impact of such duplications has important research and clinical value with respect to gene regulation and triplosensitivity traits. METHODS: We employed 4C/HiC-seq as a clinical tool to determine the functional impact of incidentally discovered GPR101 duplications on TAD integrity in three families. After defining duplications and breakpoints around GPR101 by clinical-grade and high-density aCGH, we constructed 4C/HiC chromatin contact maps for our study population and compared them with normal and active (X-LAG) controls. RESULTS: We showed that duplications involving GPR101 that preserved the centromeric invariant TAD boundary did not generate a pathogenic neo-TAD and that ectopic enhancers were not adopted. This allowed us to discount presumptive/suspected X-LAG diagnoses and GPR101 misexpression, obviating the need for intensive clinical follow-up. CONCLUSIONS: This study highlights the importance of TAD boundaries and chromatin interactions in determining the functional impact of copy number variants and provides proof-of-concept for using 4C/HiC-seq as a clinical tool to acquire crucial information for genetic counseling and to support clinical decision-making in cases of suspected TADopathies.

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.

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

We examined more than 38,000 spouse pairs from four neurodevelopmental disease cohorts and the UK Biobank to identify phenotypic and genetic patterns in parents associated with neurodevelopmental disease risk in children. We identified correlations between six phenotypes in parents and children, including correlations of clinical diagnoses such as obsessive-compulsive disorder (R=0.31-0.49, p<0.001), and two measures of sub-clinical autism features in parents affecting several autism severity measures in children, such as bi-parental mean Social Responsiveness Scale (SRS) scores affecting proband SRS scores (regression coefficient=0.11, p=0.003). We further describe patterns of phenotypic and genetic similarity between spouses, where spouses show both within- and cross-disorder correlations for seven neurological and psychiatric phenotypes, including a within-disorder correlation for depression (R=0.25-0.72, p<0.001) and a cross-disorder correlation between schizophrenia and personality disorder (R=0.20-0.57, p<0.001). Further, these spouses with similar phenotypes were significantly correlated for rare variant burden (R=0.07-0.57, p<0.0001). We propose that assortative mating on these features may drive the increases in genetic risk over generations and the appearance of "genetic anticipation" associated with many variably expressive variants. We further identified parental relatedness as a risk factor for neurodevelopmental disorders through its inverse correlations with burden and pathogenicity of rare variants and propose that parental relatedness drives disease risk by increasing genome-wide homozygosity in children (R=0.09-0.30, p<0.001). Our results highlight the utility of assessing parent phenotypes and genotypes in predicting features in children carrying variably expressive variants and counseling families carrying these variants.

Understanding the new BRD4-related syndrome: Clinical and genomic delineation with an international cohort study.

BRD4 is part of a multiprotein complex involved in loading the cohesin complex onto DNA, a fundamental process required for cohesin-mediated loop extrusion and formation of Topologically Associating Domains. Pathogenic variations in this complex have been associated with a growing number of syndromes, collectively known as cohesinopathies, the most classic being Cornelia de Lange syndrome. However, no cohort study has been conducted to delineate the clinical and molecular spectrum of BRD4-related disorder. We formed an international collaborative study, and collected 14 new patients, including two fetuses. We performed phenotype and genotype analysis, integrated prenatal findings from fetopathological examinations, phenotypes of pediatric patients and adults. We report the first cohort of patients with BRD4-related disorder and delineate the dysmorphic features at different ages. This work extends the phenotypic spectrum of cohesinopathies and characterize a new clinically relevant and recognizable pattern, distinguishable from the other cohesinopathies.

Novel Loss of Function Variant in BCKDK Causes a Treatable Developmental and Epileptic Encephalopathy.

Branched-chain amino acids (BCAA) are essential amino acids playing crucial roles in protein synthesis and brain neurotransmission. Branched-chain ketoacid dehydrogenase (BCKDH), the flux-generating step of BCAA catabolism, is tightly regulated by reversible phosphorylation of its E1α-subunit. BCKDK is the kinase responsible for the phosphorylation-mediated inactivation of BCKDH. In three siblings with severe developmental delays, microcephaly, autism spectrum disorder and epileptic encephalopathy, we identified a new homozygous in-frame deletion (c.999_1001delCAC; p.Thr334del) of BCKDK. Plasma and cerebrospinal fluid concentrations of BCAA were markedly reduced. Hyperactivity of BCKDH and over-consumption of BCAA were demonstrated by functional tests in cells transfected with the mutant BCKDK. Treatment with pharmacological doses of BCAA allowed the restoring of BCAA concentrations and greatly improved seizure control. Behavioral and developmental skills of the patients improved to a lesser extent. Importantly, a retrospective review of the newborn screening results allowed the identification of a strong decrease in BCAA concentrations on dried blood spots, suggesting that BCKDK is a new treatable metabolic disorder probably amenable to newborn screening programs.

Three years pilot of spinal muscular atrophy newborn screening turned into official program in Southern Belgium.

Three new therapies for spinal muscular atrophy (SMA) have been approved by the United States Food and Drug Administration and the European Medicines Agency since 2016. Although these new therapies improve the quality of life of patients who are symptomatic at first treatment, administration before the onset of symptoms is significantly more effective. As a consequence, newborn screening programs have been initiated in several countries. In 2018, we launched a 3-year pilot program to screen newborns for SMA in the Belgian region of Liège. This program was rapidly expanding to all of Southern Belgium, a region of approximately 55,000 births annually. During the pilot program, 136,339 neonates were tested for deletion of exon 7 of SMN1, the most common cause of SMA. Nine SMA cases with homozygous deletion were identified through this screen. Another patient was identified after presenting with symptoms and was shown to be heterozygous for the SMN1 exon 7 deletion and a point mutation on the opposite allele. These ten patients were treated. The pilot program has now successfully transitioned into the official neonatal screening program in Southern Belgium. The lessons learned during implementation of this pilot program are reported.

Newborn screening of duchenne muscular dystrophy specifically targeting deletions amenable to exon-skipping therapy.

Duchenne Muscular Dystrophy (DMD) is a lethal progressive muscle-wasting disease. New treatment strategies relying on DMD gene exon-skipping therapy have recently been approved and about 30% of patients could be amenable to exon 51, 53 or 45 skipping. We evaluated the spectrum of deletions reported in DMD registries, and designed a method to screen newborns and identify DMD deletions amenable to exon 51, 53 and 45 skipping. We developed a multiplex qPCR assay identifying hemi(homo)-zygotic deletions of the flanking exons of these therapeutic targets in DMD exons (i.e. exons 44, 46, 50, 52 and 54). We conducted an evaluation of our new method in 51 male patients with a DMD phenotype, 50 female carriers of a DMD deletion and 19 controls. Studies were performed on dried blood spots with patient's consent. We analyzed qPCR amplification curves of controls, carriers, and DMD patients to discern the presence or the absence of the target exons. Analysis of the exons flanking the exon-skipping targets permitted the identification of patients that could benefit from exon-skipping. All samples were correctly genotyped, with either presence or absence of amplification of the target exon. This proof-of-concept study demonstrates that this new assay is a highly sensitive method to identify DMD patients carrying deletions that are rescuable by exon-skipping treatment. The method is easily scalable to population-based screening. This targeted screening approach could address the new management paradigm in DMD, and could help to optimize the beneficial therapeutic effect of DMD therapies by permitting pre-symptomatic care.

16p13.11 microduplication in 45 new patients: refined clinical significance and genotype-phenotype correlations.

BACKGROUND: The clinical significance of 16p13.11 duplications remains controversial while frequently detected in patients with developmental delay (DD), intellectual deficiency (ID) or autism spectrum disorder (ASD). Previously reported patients were not or poorly characterised. The absence of consensual recommendations leads to interpretation discrepancy and makes genetic counselling challenging. This study aims to decipher the genotype-phenotype correlations to improve genetic counselling and patients' medical care. METHODS: We retrospectively analysed data from 16 013 patients referred to 12 genetic centers for DD, ID or ASD, and who had a chromosomal microarray analysis. The referring geneticists of patients for whom a 16p13.11 duplication was detected were asked to complete a questionnaire for detailed clinical and genetic data for the patients and their parents. RESULTS: Clinical features are mainly speech delay and learning disabilities followed by ASD. A significant risk of cardiovascular disease was noted. About 90% of the patients inherited the duplication from a parent. At least one out of four parents carrying the duplication displayed a similar phenotype to the propositus. Genotype-phenotype correlations show no impact of the size of the duplicated segment on the severity of the phenotype. However, NDE1 and miR-484 seem to have an essential role in the neurocognitive phenotype. CONCLUSION: Our study shows that 16p13.11 microduplications are likely pathogenic when detected in the context of DD/ID/ASD and supports an essential role of NDE1 and miR-484 in the neurocognitive phenotype. Moreover, it suggests the need for cardiac evaluation and follow-up and a large study to evaluate the aortic disease risk.

Estimating the effect size of the 15Q11.2 BP1-BP2 deletion and its contribution to neurodevelopmental symptoms: recommendations for practice.

BACKGROUND: The 15q11.2 deletion is frequently identified in the neurodevelopmental clinic. Case-control studies have associated the 15q11.2 deletion with neurodevelopmental disorders, and clinical case series have attempted to delineate a microdeletion syndrome with considerable phenotypic variability. The literature on this deletion is extensive and confusing, which is a challenge for genetic counselling. The aim of this study was to estimate the effect size of the 15q11.2 deletion and quantify its contribution to neurodevelopmental disorders. METHODS: We performed meta-analyses on new and previously published case-control studies and used statistical models trained in unselected populations with cognitive assessments. We used new (n=241) and previously published (n=150) data from a clinically referred group of deletion carriers. 15q11.2 duplications (new n=179 and previously published n=35) were used as a neutral control variant. RESULTS: The deletion decreases IQ by 4.3 points. The estimated ORs and respective frequencies in deletion carriers for intellectual disabilities, schizophrenia and epilepsy are 1.7 (3.4%), 1.5 (2%) and 3.1 (2.1%), respectively. There is no increased risk for heart malformations and autism. In the clinically referred group, the frequency and nature of symptoms in deletions are not different from those observed in carriers of the 15q11.2 duplication suggesting that most of the reported symptoms are due to ascertainment bias. CONCLUSIONS: We recommend that the deletion should be classified as 'pathogenic of mild effect size'. Since it explains only a small proportion of the phenotypic variance in carriers, it is not worth discussing in the developmental clinic or in a prenatal setting.

VPS51 biallelic variants cause microcephaly with brain malformations: A confirmatory report.

Whole exome sequencing undertaken in two siblings with delayed psychomotor development, absent speech, severe intellectual disability and postnatal microcephaly, with brain malformations consisting of cerebellar atrophy in the eldest affected and hypoplastic corpus callosum in the younger sister; revealed a homozygous intragenic deletion in VPS51, which encodes the vacuolar protein sorting-associated protein, one the four subunits of the Golgi-associated retrograde protein (GARP) and endosome-associated recycling protein (EARP) complexes that promotes the fusion of endosome-derived vesicles with the trans-Golgi network (GARP) and recycling endosomes (EARP). This observation supports a pathogenic effect of VPS51 variants, which has only been reported previously once, in a single child with microcephaly. It confirms the key role of membrane trafficking in normal brain development and homeostasis.

Newborn screening for SMA in Southern Belgium.

Approval was recently granted for a new treatment for spinal muscular atrophy (SMA). Given that the treatment is effective when administered early and the societal burden of SMA-related disability, the implementation of a newborn screening program is warranted. We describe the stepwise process that led us to launch a newborn screening program for SMA in Southern Belgium. Different political, ethical, and clinical partners were informed about this project and were involved in its governance, as were genetic and screening labs. We developed and validated a newborn screening method to specifically recognize homozygous deletions of exon 7 in the SMN1 gene. Subsequently, a 3-year pilot study has been recently initiated in one Belgian neonatal screening laboratory to cover 17.000 neonates per year. Coverage extension to all of Southern Belgium to screen 55.000 babies each year is underway.

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.

Transcriptome-wide analysis of natural antisense transcripts shows their potential role in breast cancer.

Non-coding RNAs (ncRNA) represent 1/5 of the mammalian transcript number, and 90% of the genome length is transcribed. Many ncRNAs play a role in cancer. Among them, non-coding natural antisense transcripts (ncNAT) are RNA sequences that are complementary and overlapping to those of either protein-coding (PCT) or non-coding transcripts. Several ncNATs were described as regulating protein coding gene expression on the same loci, and they are expected to act more frequently in cis compared to other ncRNAs that commonly function in trans. In this work, 22 breast cancers expressing estrogen receptors and their paired adjacent non-malignant tissues were analyzed by strand-specific RNA sequencing. To highlight ncNATs potentially playing a role in protein coding gene regulations that occur in breast cancer, three different data analysis methods were used: differential expression analysis of ncNATs between tumor and non-malignant tissues, differential correlation analysis of paired ncNAT/PCT between tumor and non-malignant tissues, and ncNAT/PCT read count ratio variation between tumor and non-malignant tissues. Each of these methods yielded lists of ncNAT/PCT pairs that were enriched in survival-associated genes. This work highlights ncNAT lists that display potential to affect the expression of protein-coding genes involved in breast cancer pathology.

Genomic studies of multiple myeloma reveal an association between X chromosome alterations and genomic profile complexity.

The genomic profile of multiple myeloma (MM) has prognostic value by dividing patients into a good prognosis hyperdiploid group and a bad prognosis nonhyperdiploid group with a higher incidence of IGH translocations. This classification, however, is inadequate and many other parameters like mutations, epigenetic modifications, and genomic heterogeneity may influence the prognosis. We performed a genomic study by array-based comparative genomic hybridization on a cohort of 162 patients to evaluate the frequency of genomic gains and losses. We identified a high frequency of X chromosome alterations leading to partial Xq duplication, often associated with inactive X (Xi) deletion in female patients. This partial X duplication could be a cytogenetic marker of aneuploidy as it is correlated with a high number of chromosomal breakages. Patient with high level of chromosomal breakage had reduced survival regardless the region implicated. A higher transcriptional level was shown for genes with potential implication in cancer and located in this altered region. Among these genes, IKBKG and IRAK1 are members of the NFKB pathway which plays an important role in MM and is a target for specific treatments. © 2016 Wiley Periodicals, Inc.

Exome copy number variation detection: Use of a pool of unrelated healthy tissue as reference sample.

An increasing number of bioinformatic tools designed to detect CNVs (copy number variants) in tumor samples based on paired exome data where a matched healthy tissue constitutes the reference have been published in the recent years. The idea of using a pool of unrelated healthy DNA as reference has previously been formulated but not thoroughly validated. As of today, the gold standard for CNV calling is still aCGH but there is an increasing interest in detecting CNVs by exome sequencing. We propose to design a metric allowing the comparison of two CNV profiles, independently of the technique used and assessed the validity of using a pool of unrelated healthy DNA instead of a matched healthy tissue as reference in exome-based CNV detection. We compared the CNV profiles obtained with three different approaches (aCGH, exome sequencing with a matched healthy tissue as reference, exome sequencing with a pool of eight unrelated healthy tissue as reference) on three multiple myeloma samples. We show that the usual analyses performed to compare CNV profiles (deletion/amplification ratios and CNV size distribution) lack in precision when confronted with low LRR values, as they only consider the binary status of each CNV. We show that the metric-based distance constitutes a more accurate comparison of two CNV profiles. Based on these analyses, we conclude that a reliable picture of CNV alterations in multiple myeloma samples can be obtained from whole-exome sequencing in the absence of a matched healthy sample.

Combined treatment with octreotide LAR and pegvisomant in patients with pituitary gigantism: clinical evaluation and genetic screening.

INTRODUCTION: Pituitary gigantism is a rare condition caused by growth hormone secreting hypersecretion, usually by a pituitary tumor. Acromegaly and gigantism cases that have a genetic cause are challenging to treat, due to large tumor size and poor responses to some medical therapies (e.g. AIP mutation affected cases and those with X-linked acrogigantism syndrome). MATERIALS AND METHODS: We performed a retrospective study to identify gigantism cases among 160 somatotropinoma patients treated between 1985 and 2015 at the University Hospital of Caracas, Venezuela. We studied clinical details at diagnosis, hormonal responses to therapy and undertook targeted genetic testing. Among the 160 cases, eight patients (six males; 75 %) were diagnosed with pituitary gigantism and underwent genetic analysis that included array comparative genome hybridization for Xq26.3 duplications. RESULTS: All patients had GH secreting pituitary macroadenomas that were difficult to control with conventional treatment options, such as surgery or primary somatostatin receptor ligand (SRL) therapy. Combined therapy (long-acting SRL and pegvisomant) as primary treatment or after pituitary surgery and radiotherapy permitted the normalization of IGF-1 levels and clinical improvement. Novel AIP mutations were the found in three patients. None of the patients had Xq26.3 microduplications. CONCLUSIONS: Treatment of pituitary gigantism is frequently challenging; delayed control increases the harmful effects of GH excess, such as, excessive stature and symptom burden, so early diagnosis and effective treatment are particularly important in these cases.

Somatic mosaicism underlies X-linked acrogigantism syndrome in sporadic male subjects.

Somatic mosaicism has been implicated as a causative mechanism in a number of genetic and genomic disorders. X-linked acrogigantism (XLAG) syndrome is a recently characterized genomic form of pediatric gigantism due to aggressive pituitary tumors that is caused by submicroscopic chromosome Xq26.3 duplications that include GPR101 We studied XLAG syndrome patients (n= 18) to determine if somatic mosaicism contributed to the genomic pathophysiology. Eighteen subjects with XLAG syndrome caused by Xq26.3 duplications were identified using high-definition array comparative genomic hybridization (HD-aCGH). We noted that males with XLAG had a decreased log2ratio (LR) compared with expected values, suggesting potential mosaicism, whereas females showed no such decrease. Compared with familial male XLAG cases, sporadic males had more marked evidence for mosaicism, with levels of Xq26.3 duplication between 16.1 and 53.8%. These characteristics were replicated using a novel, personalized breakpoint junction-specific quantification droplet digital polymerase chain reaction (ddPCR) technique. Using a separate ddPCR technique, we studied the feasibility of identifying XLAG syndrome cases in a distinct patient population of 64 unrelated subjects with acromegaly/gigantism, and identified one female gigantism patient who had had increased copy number variation (CNV) threshold for GPR101 that was subsequently diagnosed as having XLAG syndrome on HD-aCGH. Employing a combination of HD-aCGH and novel ddPCR approaches, we have demonstrated, for the first time, that XLAG syndrome can be caused by variable degrees of somatic mosaicism for duplications at chromosome Xq26.3. Somatic mosaicism was shown to occur in sporadic males but not in females with XLAG syndrome, although the clinical characteristics of the disease were similarly severe in both sexes.