Whole exome sequencing of families with 1q21.1 microdeletion or microduplication.

Recurrent microduplications/microdeletions of 1q21.1 are characterized by variable phenotypes ranging from normal development to developmental delay (DD) and congenital anomalies. Their interpretation is challenging especially in families with affected and unaffected carriers. We used whole exome sequencing (WES) to look for sequence variants in two male probands with inherited 1q21.1 CNVs that could explain their more severe phenotypes. One proband had a 1q21.1 deletion transmitted from maternal grandmother, while the other had a paternal duplication. We found mutations in five genes (SMPD1, WNK3, NOS1, ATF6, and EFHC1) that could contribute to the more severe phenotype in the probands in comparison to their mildly affected or unaffected 1q21.1 CNV carrying relatives. Interestingly, all genes have roles in stress responses (oxidative/Endoplasmic Reticulum (ER)/osmotic). One of the variants was in an X-linked gene WNK3 and segregated with the developmental features and X inactivation pattern in the family with 1q21.1 deletion transmitted from maternal grandmother. In silico analysis of all rare deleterious variants in both probands identified enrichment in nervous system diseases, metabolic pathways, protein processing in the ER and protein export. Our studies suggest that rare deleterious variants outside of the 1q21.1 CNV, individually or as a pool, could contribute to phenotypic variability in carriers of this CNV. Rare deleterious variants in stress response genes are of interest and raise the possibility of susceptibility of carriers to variable environmental influences. Next generation sequencing of additional familial cases with 1q21.1 CNV could further help determine the possible causes of phenotypic variability in carriers of this CNV.

Whole exome sequencing in recurrent early pregnancy loss.

STUDY HYPOTHESIS: Exome sequencing can identify genetic causes of idiopathic recurrent pregnancy loss (RPL). STUDY FINDING: We identified compound heterozygous deleterious mutations affecting DYNC2H1 and ALOX15 in two out of four families with RPL. Both genes have a role in early development. Bioinformatics analysis of all genes with rare and putatively pathogenic mutations in miscarriages and couples showed enrichment in pathways relevant to pregnancy loss, including the complement and coagulation cascades pathways. WHAT IS KNOWN ALREADY: Next generation sequencing (NGS) is increasingly being used to identify known and novel gene mutations in children with developmental delay and in fetuses with ultrasound-detected anomalies. In contrast, NGS is rarely used to study pregnancy loss. Chromosome microarray analysis detects putatively causative DNA copy number variants (CNVs) in ∼2% of miscarriages and CNVs of unknown significance (predominantly parental in origin) in up to 40% of miscarriages. Therefore, a large number of miscarriages still have an unknown cause. STUDY DESIGN, SAMPLES/MATERIALS, METHODS: Whole exome sequencing (WES) was performed using Illumina HiSeq 2000 platform on seven euploid miscarriages from four families with RPL. Golden Helix SVS v8.1.5 was used for data assessment and inheritance analysis for deleterious DNA variants predicted to severely disrupt protein-coding genes by introducing a frameshift, loss of the stop codon, gain of the stop codon, changes in splicing or the initial codon. Webgestalt (http://bioinfo.vanderbilt.edu/webgestalt/) was used for pathway and disease association enrichment analysis of a gene pool containing putatively pathogenic variants in miscarriages and couples in comparison to control gene pools. MAIN RESULTS AND THE ROLE OF CHANCE: Compound heterozygous mutations in DYNC2H1 and ALOX15 were identified in miscarriages from two families with RPL. DYNC2H1 is involved in cilia biogenesis and has been associated with fetal lethality in humans. ALOX15 is expressed in placenta and its dysregulation has been associated with inflammation, placental, dysfunction, abnormal oxidative stress response and angiogenesis. The pool of putatively pathogenic single nucleotide variants (SNVs) and small insertions and deletions (indels) detected in the miscarriages showed enrichment in 'complement and coagulation cascades pathway', and 'ciliary motility disorders'. We conclude that CNVs, individual SNVs and pool of deleterious gene mutations identified by exome sequencing could contribute to RPL. LIMITATIONS, REASONS FOR CAUTION: The size of our sample cohort is small. The functional effect of candidate mutations should be evaluated to determine whether the mutations are causative. WIDER IMPLICATIONS OF THE FINDINGS: This is the first study to assess whether SNVs may contribute to the pathogenesis of miscarriage. Furthermore, our findings suggest that collective effect of mutations in relevant biological pathways could be implicated in RPL. STUDY FUNDING AND COMPETING INTERESTS: The study was funded by Canadian Institutes of Health Research (grant MOP 106467) and Michael Smith Foundation of Health Research Career Scholar salary award to ERS.

Exome sequencing identifies pathogenic variants of VPS13B in a patient with familial 16p11.2 duplication.

BACKGROUND: The recurrent microduplication of 16p11.2 (dup16p11.2) is associated with a broad spectrum of neurodevelopmental disorders (NDD) confounded by incomplete penetrance and variable expressivity. This inter- and intra-familial clinical variability highlights the importance of personalized genetic counselling in individuals at-risk. CASE PRESENTATION: In this study, we performed whole exome sequencing (WES) to look for other genomic alterations that could explain the clinical variability in a family with a boy presenting with NDD who inherited the dup16p11.2 from his apparently healthy mother. We identified novel splicing variants of VPS13B (8q22.2) in the proband with compound heterozygous inheritance. Two VPS13B mutations abolished the canonical splice sites resulting in low RNA expression in transformed lymphoblasts of the proband. VPS13B mutation causes Cohen syndrome (CS) consistent with the proband's phenotype (intellectual disability (ID), microcephaly, facial gestalt, retinal dystrophy, joint hypermobility and neutropenia). The new diagnosis of CS has important health implication for the proband, provides the opportunity for more meaningful and accurate genetic counselling for the family; and underscores the importance of longitudinally following patients for evolving phenotypic features. CONCLUSIONS: This is the first report of a co-occurrence of pathogenic variants with familial dup16p11.2. Our finding suggests that the variable expressivity among carriers of rare putatively pathogenic CNVs such as dup16p11.2 warrants further study by WES and individualized genetic counselling of families with such CNVs.

Variant ATRX syndrome with dysfunction of ATRX and MAGT1 genes.

A 0.8 kb intronic duplication in MAGT1 and a single base pair deletion in the last exon of ATRX were identified using a chromosome X-specific microarray and exome sequencing in a family with five males demonstrating intellectual disability (ID) and unusual skin findings (e.g., generalized pruritus). MAGT1 is an Mg²âº transporter previously associated with primary immunodeficiency and ID, whereas mutations in ATRX cause ATRX-ID syndrome. In patient cells, the function of ATRX was demonstrated to be abnormal based on altered RNA/protein expression, hypomethylation of rDNA, and abnormal cytokinesis. Dysfunction of MAGT1 was reflected in reduced RNA/protein expression and Mg²âº influx. The mutation in ATRX most likely explains the ID, whereas MAGT1 disruption could be linked to abnormal skin findings, as normal magnesium homeostasis is necessary for skin health. This work supports observations that multiple mutations collectively contribute to the phenotypic variability of syndromic ID, and emphasizes the importance of correlating clinical phenotype with genomic and cell function analyses.

Complex Autism Spectrum Disorder with Epilepsy, Strabismus and Self-Injurious Behaviors in a Patient with a De Novo Heterozygous POLR2A Variant.

Autism spectrum disorder (ASD) describes a complex and heterogenous group of neurodevelopmental disorders. Whole genome sequencing continues to shed light on the multifactorial etiology of ASD. Dysregulated transcriptional pathways have been implicated in neurodevelopmental disorders. Emerging evidence suggests that de novo POLR2A variants cause a newly described phenotype called 'Neurodevelopmental Disorder with Hypotonia and Variable Intellectual and Behavioral Abnormalities' (NEDHIB). The variable phenotype manifests with a spectrum of features; primarily early onset hypotonia and delay in developmental milestones. In this study, we investigate a patient with complex ASD involving epilepsy and strabismus. Whole genome sequencing of the proband−parent trio uncovered a novel de novo POLR2A variant (c.1367T>C, p. Val456Ala) in the proband. The variant appears deleterious according to in silico tools. We describe the phenotype in our patient, who is now 31 years old, draw connections between the previously reported phenotypes and further delineate this emerging neurodevelopmental phenotype. This study sheds new insights into this neurodevelopmental disorder, and more broadly, the genetic etiology of ASD.

Contribution of Multiple Inherited Variants to Autism Spectrum Disorder (ASD) in a Family with 3 Affected Siblings.

Autism Spectrum Disorder (ASD) is the most common neurodevelopmental disorder in children and shows high heritability. However, how inherited variants contribute to ASD in multiplex families remains unclear. Using whole-genome sequencing (WGS) in a family with three affected children, we identified multiple inherited DNA variants in ASD-associated genes and pathways (RELN, SHANK2, DLG1, SCN10A, KMT2C and ASH1L). All are shared among the three children, except ASH1L, which is only present in the most severely affected child. The compound heterozygous variants in RELN, and the maternally inherited variant in SHANK2, are considered to be major risk factors for ASD in this family. Both genes are involved in neuron activities, including synaptic functions and the GABAergic neurotransmission system, which are highly associated with ASD pathogenesis. DLG1 is also involved in synapse functions, and KMT2C and ASH1L are involved in chromatin organization. Our data suggest that multiple inherited rare variants, each with a subthreshold and/or variable effect, may converge to certain pathways and contribute quantitatively and additively, or alternatively act via a 2nd-hit or multiple-hits to render pathogenicity of ASD in this family. Additionally, this multiple-hits model further supports the quantitative trait hypothesis of a complex genetic, multifactorial etiology for the development of ASDs.

Exome sequencing identified a de novo mutation of PURA gene in a patient with familial Xp22.31 microduplication.

The clinical significance of Xp22.31 microduplication is controversial as it is reported in subjects with developmental delay (DD), their unaffected relatives and unrelated controls. We performed multifaceted studies in a family of a boy with hypotonia, dysmorphic features and DD who carried a 600 Kb Xp22.31 microduplication (7515787-8123310bp, hg19) containing two genes, VCX and PNPLA4. The duplication was transmitted from his cognitively normal maternal grandfather. We found no evidence of the duplication causing the proband's DD and congenital anomalies based on unaltered expression of PNPLA4 in the proband and his mother in comparison to controls and preferential activation of the paternal chromosome X with Xp22.31 duplication in proband's mother. However, a de novo, previously reported deleterious, missense mutation in Pur-alpha gene (PURA) (5q31.2), with a role in neuronal differentiation was detected in the proband by exome sequencing. We propose that the variability in the phenotype in carriers of Xp22.31 microduplication can be due to a second and more deleterious genetic mutation in more severely affected carriers. Widespread use of whole genome next generation sequencing in families with Xp22.31 CNV could help identify such cases.

Identifying candidate genes for 2p15p16.1 microdeletion syndrome using clinical, genomic, and functional analysis.

The 2p15p16.1 microdeletion syndrome has a core phenotype consisting of intellectual disability, microcephaly, hypotonia, delayed growth, common craniofacial features, and digital anomalies. So far, more than 20 cases of 2p15p16.1 microdeletion syndrome have been reported in the literature; however, the size of the deletions and their breakpoints vary, making it difficult to identify the candidate genes. Recent reports pointed to 4 genes (XPO1, USP34, BCL11A, and REL) that were included, alone or in combination, in the smallest deletions causing the syndrome. Here, we describe 8 new patients with the 2p15p16.1 deletion and review all published cases to date. We demonstrate functional deficits for the above 4 candidate genes using patients' lymphoblast cell lines (LCLs) and knockdown of their orthologs in zebrafish. All genes were dosage sensitive on the basis of reduced protein expression in LCLs. In addition, deletion of XPO1, a nuclear exporter, cosegregated with nuclear accumulation of one of its cargo molecules (rpS5) in patients' LCLs. Other pathways associated with these genes (e.g., NF-κB and Wnt signaling as well as the DNA damage response) were not impaired in patients' LCLs. Knockdown of xpo1a, rel, bcl11aa, and bcl11ab resulted in abnormal zebrafish embryonic development including microcephaly, dysmorphic body, hindered growth, and small fins as well as structural brain abnormalities. Our multifaceted analysis strongly implicates XPO1, REL, and BCL11A as candidate genes for 2p15p16.1 microdeletion syndrome.

Nitroglycerin upregulates matrix metalloproteinase expression by human macrophages.

OBJECTIVES: This study aimed to determine whether nitroglycerin (NTG) treatment affects matrix metalloproteinase (MMP) gene expression and activities in human macrophages. BACKGROUND: Nitroglycerin is one of the most frequently used therapeutic agents for the symptomatic relief of stable or unstable coronary artery disease; however, its effects on vascular biology are poorly characterized. Despite its powerful vasodilator activity, NTG has not been shown to improve outcomes in coronary disease. We now describe evidence that NTG has potentially pro-inflammatory effects in human monocyte-derived macrophages (MDMs). METHODS: Human monocytes were isolated from whole blood by elutriation and allowed to differentiate into macrophages over eight to 10 days. The MDMs were then treated for 4 or 24 h with control media, pharmacologically relevant doses of NTG or other nitric oxide donors. Matrix metalloproteinase activity was measured by zymography, protein levels measured by enzyme-linked immunosorbent assay and messenger ribonucleic acid (mRNA) levels were quantified by competitive reverse transcription-polymerase chain reaction. RESULTS: The major MMP expressed by MDMs was MMP-9. Nitroglycerin treatment stimulated a dose-dependent increase in MMP-9 mRNA levels (NTG 200 pmol: 193 +/- 6% and NTG 2,000 pmol: 372 +/- 9% compared to controls, p < 0.005) and MMP-9 activity (NTG 200: 142 +/- 5.5% and NTG 2,000: 167 +/- 11% compared to controls, p < 0.005). Nitroglycerin 2,000 pmol also increased MMP-2 and MMP-7 mRNA levels to 187 +/- 8% and 183 +/- 21% of control values, respectively (p < 0.05). Furthermore, tissue inhibitor of metalloproteinase (TIMP)-1 (the major tissue inhibitor of MMPs) mRNA and protein levels were decreased in NTG 2,000 pmol-treated MDMs compared with control cells (mRNA: 67 +/- 7%, p < 0.005; protein: 45 +/- 5%, p < 0.005). CONCLUSIONS: Nitroglycerin in pharmacologically relevant concentrations activates MMP but represses TIMP expression in human macrophages. The subsequent imbalance in MMP/TIMP expression associated with NTG treatment could promote matrix degradation, with potentially adverse effects on plaque stability.

Functional consequences of copy number variants in miscarriage.

BACKGROUND: The presence of unique copy number variations (CNVs) in miscarriages suggests that their integral genes have a role in maintaining early pregnancy. In our previous work, we identified 19 unique CNVs in ~40% of studied euploid miscarriages, which were predominantly familial in origin. In our current work, we assessed their relevance to miscarriage by expression analysis of 14 genes integral to CNVs in available miscarriage chorionic villi. As familial CNVs could cause miscarriage due to imprinting effect, we investigated the allelic expression of one of the genes (TIMP2) previously suggested to be maternally expressed in placenta and involved in placental remodelling and embryo development. RESULTS: Six out of fourteen genes had detectable expression in villi and for three genes the RNA and protein expression was altered due to maternal CNVs. These genes were integral to duplication on Xp22.2 (TRAPPC2 and OFD1) or disrupted by a duplication mapping to 17q25.3 (TIMP2). RNA and protein expression was increased for TRAPPC2 and OFD1 and reduced for TIMP2 in carrier miscarriages. The three genes have roles in processes important for pregnancy development such as extracellular matrix homeostasis (TIMP2 and TRAPPC2) and cilia function (OFD1). TIMP2 allelic expression was not affected by the CNV in miscarriages in comparison to control elective terminations. CONCLUSION: We propose that functional studies of CNVs could help determine if and how the miscarriage CNVs affect the expression of integral genes. In case of parental CNVs, assessment of the function of their integral genes in parental reproductive tissues should be also considered in the future, especially if they affect processes relevant for pregnancy development and support.

Phenotypic and functional consequences of haploinsufficiency of genes from exocyst and retinoic acid pathway due to a recurrent microdeletion of 2p13.2.

BACKGROUND: Rare, recurrent genomic imbalances facilitate the association of genotype with abnormalities at the "whole body" level. However, at the cellular level, the functional consequences of recurrent genomic abnormalities and how they can be linked to the phenotype are much less investigated. METHOD AND RESULTS: We report an example of a functional analysis of two genes from a new, overlapping microdeletion of 2p13.2 region (from 72,140,702-72,924,626). The subjects shared intellectual disability (ID), language delay, hyperactivity, facial asymmetry, ear malformations, and vertebral and/or craniofacial abnormalities. The overlapping region included two genes, EXOC6B and CYP26B1, which are involved in exocytosis/Notch signaling and retinoic acid (RA) metabolism, respectively, and are of critical importance for early morphogenesis, symmetry as well as craniofacial, skeleton and brain development. The abnormal function of EXOC6B was documented in patient lymphoblasts by its reduced expression and with perturbed expression of Notch signaling pathway genes HES1 and RBPJ, previously noted to be the consequence of EXOC6B dysfunction in animal and cell line models. Similarly, the function of CYP26B1 was affected by the deletion since the retinoic acid induced expression of this gene in patient lymphoblasts was significantly lower compared to controls (8% of controls). CONCLUSION: Haploinsufficiency of CYP26B1 and EXOC6B genes involved in retinoic acid and exocyst/Notch signaling pathways, respectively, has not been reported previously in humans. The developmental anomalies and phenotypic features of our subjects are in keeping with the dysfunction of these genes, considering their known role. Documenting their dysfunction at the cellular level in patient cells enhanced our understanding of biological processes which contribute to the clinical phenotype.

Understanding the impact of 1q21.1 copy number variant.

BACKGROUND: 1q21.1 Copy Number Variant (CNV) is associated with a highly variable phenotype ranging from congenital anomalies, learning deficits/intellectual disability (ID), to a normal phenotype. Hence, the clinical significance of this CNV can be difficult to evaluate. Here we described the consequences of the 1q21.1 CNV on genome-wide gene expression and function of selected candidate genes within 1q21.1 using cell lines from clinically well described subjects. METHODS AND RESULTS: Eight subjects from 3 families were included in the study: six with a 1q21.1 deletion and two with a 1q21.1 duplication. High resolution Affymetrix 2.7M array was used to refine the 1q21.1 CNV breakpoints and exclude the presence of secondary CNVs of pathogenic relevance. Whole genome expression profiling, studied in lymphoblast cell lines (LBCs) from 5 subjects, showed enrichment of genes from 1q21.1 in the top 100 genes ranked based on correlation of expression with 1q21.1 copy number. The function of two top genes from 1q21.1, CHD1L/ALC1 and PRKAB2, was studied in detail in LBCs from a deletion and a duplication carrier. CHD1L/ALC1 is an enzyme with a role in chromatin modification and DNA damage response while PRKAB2 is a member of the AMP kinase complex, which senses and maintains systemic and cellular energy balance. The protein levels for CHD1L/ALC1 and PRKAB2 were changed in concordance with their copy number in both LBCs. A defect in chromatin remodeling was documented based on impaired decatenation (chromatid untangling) checkpoint (DCC) in both LBCs. This defect, reproduced by CHD1L/ALC1 siRNA, identifies a new role of CHD1L/ALC1 in DCC. Both LBCs also showed elevated levels of micronuclei following treatment with a Topoisomerase II inhibitor suggesting increased DNA breaks. AMP kinase function, specifically in the deletion containing LBCs, was attenuated. CONCLUSION: Our studies are unique as they show for the first time that the 1q21.1 CNV not only causes changes in the expression of its key integral genes, associated with changes at the protein level, but also results in changes in their known function, in the case of AMPK, and newly identified function such as DCC activation in the case of CHD1L/ALC1. Our results support the use of patient lymphoblasts for dissecting the functional sequelae of genes integral to CNVs in carrier cell lines, ultimately enhancing understanding of biological processes which may contribute to the clinical phenotype.