Dystonia-causing mutations in the transcription factor THAP1 disrupt HCFC1 cofactor recruitment and alter gene expression.

Thanatos-associated protein domain containing, apoptosis-associated protein 1 (THAP1), the gene mutated in DYT6 dystonia, encodes a transcription factor. While the N-terminal THAP domain allows for specific DNA-binding, the functional relevance of the other regions is largely unknown. The C-terminus contains a 4-amino-acid-spanning host cell factor 1 (HCFC1)-binding domain (HBM) that mediates the interaction with HCFC1. Interestingly, three mutations affecting the HBM (p.N136S, p.N136K, p.Y137C) have been reported in dystonia patients. We investigated the consequences of these mutations on the interaction of THAP1 with HCFC1 and demonstrated that all three mutations abolished HCFC1-THAP1 complex formation. Notably, HCFC1 co-localization was found in >90% of the almost 3,500 chromatin regions loaded with THAP1 in publicly available genome-wide ChIP data. By siRNA-mediated depletion of HCFC1, we detected an increase of THAP1 expression, indicating a co-repressor activity of HCFC1 for THAP1. Quantitative ChIP on selected promoters revealed that none of the mutations significantly decreased the DNA-binding ability of THAP1 while HCFC1 binding was highly reduced. Our findings indicate a THAP1-mediated recruitment of HCFC1 to THAP1 target sites. Of note, dystonia-causing mutations within the HBM in THAP1 abolished this interaction. Thus, we demonstrate disrupted THAP1-HCFC1 complex formation as another mechanism of dystonia-causing mutations leading to transcriptional dysregulation.

Mutations in chromatin regulators functionally link Cornelia de Lange syndrome and clinically overlapping phenotypes.

The coordinated tissue-specific regulation of gene expression is essential for the proper development of all organisms. Mutations in multiple transcriptional regulators cause a group of neurodevelopmental disorders termed "transcriptomopathies" that share core phenotypical features including growth retardation, developmental delay, intellectual disability and facial dysmorphism. Cornelia de Lange syndrome (CdLS) belongs to this class of disorders and is caused by mutations in different subunits or regulators of the cohesin complex. Herein, we report on the clinical and molecular characterization of seven patients with features overlapping with CdLS who were found to carry mutations in chromatin regulators previously associated to other neurodevelopmental disorders that are frequently considered in the differential diagnosis of CdLS. The identified mutations affect the methyltransferase-encoding genes KMT2A and SETD5 and different subunits of the SWI/SNF chromatin-remodeling complex. Complementary to this, a patient with Coffin-Siris syndrome was found to carry a missense substitution in NIPBL. Our findings indicate that mutations in a variety of chromatin-associated factors result in overlapping clinical phenotypes, underscoring the genetic heterogeneity that should be considered when assessing the clinical and molecular diagnosis of neurodevelopmental syndromes. It is clear that emerging molecular mechanisms of chromatin dysregulation are central to understanding the pathogenesis of these clinically overlapping genetic disorders.

Prediction of Causal Candidate Genes in Coronary Artery Disease Loci.

OBJECTIVE: Genome-wide association studies have to date identified 159 significant and suggestive loci for coronary artery disease (CAD). We now report comprehensive bioinformatics analyses of sequence variation in these loci to predict candidate causal genes. APPROACH AND RESULTS: All annotated genes in the loci were evaluated with respect to protein-coding single-nucleotide polymorphism and gene expression parameters. The latter included expression quantitative trait loci, tissue specificity, and miRNA binding. High priority candidate genes were further identified based on literature searches and our experimental data. We conclude that the great majority of causal variations affecting CAD risk occur in noncoding regions, with 41% affecting gene expression robustly versus 6% leading to amino acid changes. Many of these genes differed from the traditionally annotated genes, which was usually based on proximity to the lead single-nucleotide polymorphism. Indeed, we obtained evidence that genetic variants at CAD loci affect 98 genes which had not been linked to CAD previously. CONCLUSIONS: Our results substantially revise the list of likely candidates for CAD and suggest that genome-wide association studies efforts in other diseases may benefit from similar bioinformatics analyses.

Hidden mutations in Cornelia de Lange syndrome limitations of sanger sequencing in molecular diagnostics.

Cornelia de Lange syndrome (CdLS) is a well-characterized developmental disorder. The genetic cause of CdLS is a mutation in one of five associated genes (NIPBL, SMC1A, SMC3, RAD21, and HDAC8) accounting for about 70% of cases. To improve our current molecular diagnostic and to analyze some of CdLS candidate genes, we developed and established a gene panel approach. Because recent data indicate a high frequency of mosaic NIPBL mutations that were not detected by conventional sequencing approaches of blood DNA, we started to collect buccal mucosa (BM) samples of our patients that were negative for mutations in the known CdLS genes. Here, we report the identification of three mosaic NIPBL mutations by our high-coverage gene panel sequencing approach that were undetected by classical Sanger sequencing analysis of BM DNA. All mutations were confirmed by the use of highly sensitive SNaPshot fragment analysis using DNA from BM, urine, and fibroblast samples. In blood samples, we could not detect the respective mutation. Finally, in fibroblast samples from all three patients, Sanger sequencing could identify all the mutations. Thus, our study highlights the need for highly sensitive technologies in molecular diagnostic of CdLS to improve genetic diagnosis and counseling of patients and their families.

Whole-exome sequencing in an extended family with myocardial infarction unmasks familial hypercholesterolemia.

BACKGROUND: Familial hypercholesterolemia (FH) is an autosomal-dominant disease leading to markedly elevated low-density lipoprotein (LDL) cholesterol levels and increased risk for premature myocardial infarction (MI). Mutation carriers display variable LDL cholesterol levels, which may obscure the diagnosis. We examined by whole-exome sequencing a family in which multiple myocardial infarctions occurred at a young age with unclear etiology. METHODS: Whole-exome sequencing of three affected family members, validation of the identified variant with Sanger-sequencing, and subsequent co-segregation analysis in the family. RESULTS: The index patient (LDL cholesterol 188 mg/dL) was referred for molecular-genetic investigations. He had coronary artery bypass graft (CABG) at the age of 59 years; 12 out of 15 1st, 2nd and 3rd degree relatives were affected with coronary artery disease (CAD) and/or premature myocardial infarction (MI). We sequenced the whole-exome of the patient and two cousins with premature MI. After filtering, we were left with a potentially disease causing variant in the LDL receptor (LDLR) gene, which we validated by Sanger-sequencing (nucleotide substitution in the acceptor splice-site of exon 10, c.1359-1G > A). Sequencing of all family members available for genetic analysis revealed co-segregation of the variant with CAD (LOD 3.0) and increased LDLC (>190 mg/dL), following correction for statin treatment (LOD 4.3). Interestingly, mutation carriers presented with highly variable corrected (183-354 mg/dL) and on-treatment LDL levels (116-274 mg/dL) such that the diagnosis of FH in this family was made only after the molecular-genetic analysis. CONCLUSION: Even in families with unusual clustering of CAD FH remains to be underdiagnosed, which underscores the need for implementation of systematic screening programs. Whole-exome sequencing may facilitate identification of disease-causing variants in families with unclear etiology of MI and enable preventive treatment of mutation carriers in a more timely fashion.

Analysis of candidate genes for cleft lip ± cleft palate using murine single-cell expression data.

Introduction: Cleft lip ± cleft palate (CL/P) is one of the most common birth defects. Although research has identified multiple genetic risk loci for different types of CL/P (i.e., syndromic or non-syndromic forms), determining the respective causal genes and understanding the relevant functional networks remain challenging. The recent introduction of single-cell RNA sequencing (scRNA-seq) has provided novel opportunities to study gene expression patterns at cellular resolution. The aims of our study were to: (i) aggregate available scRNA-seq data from embryonic mice and provide this as a resource for the craniofacial community; and (ii) demonstrate the value of these data in terms of the investigation of the gene expression patterns of CL/P candidate genes. Methods and Results: First, two published scRNA-seq data sets from embryonic mice were re-processed, i.e., data representing the murine time period of craniofacial development: (i) facial data from embryonic day (E) E11.5; and (ii) whole embryo data from E9.5-E13.5 from the Mouse Organogenesis Cell Atlas (MOCA). Marker gene expression analyses demonstrated that at E11.5, the facial data were a high-resolution representation of the MOCA data. Using CL/P candidate gene lists, distinct groups of genes with specific expression patterns were identified. Among others we identified that a co-expression network including Irf6, Grhl3 and Tfap2a in the periderm, while it was limited to Irf6 and Tfap2a in palatal epithelia, cells of the ectodermal surface, and basal cells at the fusion zone. The analyses also demonstrated that additional CL/P candidate genes (e.g., Tpm1, Arid3b, Ctnnd1, and Wnt3) were exclusively expressed in Irf6+ facial epithelial cells (i.e., as opposed to Irf6- epithelial cells). The MOCA data set was finally used to investigate differences in expression profiles for candidate genes underlying different types of CL/P. These analyses showed that syndromic CL/P genes (syCL/P) were expressed in significantly more cell types than non-syndromic CL/P candidate genes (nsCL/P). Discussion: The present study illustrates how scRNA-seq data can empower research on craniofacial development and disease.

Human variation in population-wide gene expression data predicts gene perturbation phenotype.

Population-scale datasets of healthy individuals capture genetic and environmental factors influencing gene expression. The expression variance of a gene of interest (GOI) can be exploited to set up a quasi loss- or gain-of-function "in population" experiment. We describe here an approach, huva (human variation), taking advantage of population-scale multi-layered data to infer gene function and relationships between phenotypes and expression. Within a reference dataset, huva derives two experimental groups with LOW or HIGH expression of the GOI, enabling the subsequent comparison of their transcriptional profile and functional parameters. We demonstrate that this approach robustly identifies the phenotypic relevance of a GOI allowing the stratification of genes according to biological functions, and we generalize this concept to almost 16,000 genes in the human transcriptome. Additionally, we describe how huva predicts monocytes to be the major cell type in the pathophysiology of STAT1 mutations, evidence validated in a clinical cohort.

Alveolar macrophages in early stage COPD show functional deviations with properties of impaired immune activation.

Despite its high prevalence, the cellular and molecular mechanisms of chronic obstructive pulmonary disease (COPD) are far from being understood. Here, we determine disease-related changes in cellular and molecular compositions within the alveolar space and peripheral blood of a cohort of COPD patients and controls. Myeloid cells were the largest cellular compartment in the alveolar space with invading monocytes and proliferating macrophages elevated in COPD. Modeling cell-to-cell communication, signaling pathway usage, and transcription factor binding predicts TGF-ß1 to be a major upstream regulator of transcriptional changes in alveolar macrophages of COPD patients. Functionally, macrophages in COPD showed reduced antigen presentation capacity, accumulation of cholesteryl ester, reduced cellular chemotaxis, and mitochondrial dysfunction, reminiscent of impaired immune activation.

A familial congenital heart disease with a possible multigenic origin involving a mutation in BMPR1A.

The genetics of many congenital heart diseases (CHDs) can only unsatisfactorily be explained by known chromosomal or Mendelian syndromes. Here, we present sequencing data of a family with a potentially multigenic origin of CHD. Twelve of nineteen family members carry a familial mutation [NM_004329.2:c.1328 G > A (p.R443H)] which encodes a predicted deleterious variant of BMPR1A. This mutation co-segregates with a linkage region on chromosome 1 that associates with the emergence of severe CHDs including Ebstein's anomaly, atrioventricular septal defect, and others. We show that the continuous overexpression of the zebrafish homologous mutation bmpr1aap.R438H within endocardium causes a reduced AV valve area, a downregulation of Wnt/ß-catenin signalling at the AV canal, and growth of additional tissue mass in adult zebrafish hearts. This finding opens the possibility of testing genetic interactions between BMPR1A and other candidate genes within linkage region 1 which may provide a first step towards unravelling more complex genetic patterns in cardiovascular disease aetiology.

New NR5A1 mutations and phenotypic variations of gonadal dysgenesis.

Mutations in NR5A1 have been reported as a frequent cause of 46,XY disorders of sex development (DSD) associated to a broad phenotypic spectrum ranging from infertility, ambiguous genitalia, anorchia to gonadal dygenesis and female genitalia. Here we present the clinical follow up of four 46,XY DSD patients with three novel heterozygous mutations in the NR5A1 gene leading to a p.T40P missense mutation and a p.18DKVSG22del nonframeshift deletion in the DNA-binding domain and a familiar p.Y211Tfs*83 frameshift mutation. Functional analysis of the missense and nonframeshift mutation revealed a deleterious character with loss of DNA-binding and transactivation capacity. Both, the mutations in the DNA-binding domain, as well as the familiar frameshift mutation are associated with highly variable endocrine values and phenotypic appearance. Phenotypes vary from males with spontaneous puberty, substantial testosterone production and possible fertility to females with and without Müllerian structures and primary amenorrhea. Exome sequencing of the sibling's family revealed TBX2 as a possible modifier of gonadal development in patients with NR5A1 mutations.

A genomic exploration identifies mechanisms that may explain adverse cardiovascular effects of COX-2 inhibitors.

Cyclooxygenase-2 inhibitors (coxibs) are characterized by multiple molecular off-target effects and increased coronary artery disease (CAD) risk. Here, we systematically explored common variants of genes representing molecular targets of coxibs for association with CAD. Given a broad spectrum of pleiotropic effects of coxibs, our intention was to narrow potential mechanisms affecting CAD risk as we hypothesized that the affected genes may also display genomic signals of coronary disease risk. A Drug Gene Interaction Database search identified 47 gene products to be affected by coxibs. We traced association signals in 200-kb regions surrounding these genes in 84,813 CAD cases and 202,543 controls. Based on a threshold of 1 × 10-5 (Bonferroni correction for 3131 haplotype blocks), four gene loci yielded significant associations. The lead SNPs were rs7270354 (MMP9), rs4888383 (BCAR1), rs6905288 (VEGFA1), and rs556321 (CACNA1E). By additional genotyping, rs7270354 at MMP9 and rs4888383 at BCAR1 also reached the established GWAS threshold for genome-wide significance. The findings demonstrate overlap of genes affected by coxibs and those mediating CAD risk and points to further mechanisms, which are potentially responsible for coxib-associated CAD risk. The novel approach furthermore suggests that genetic studies may be useful to explore the clinical relevance of off-target drug effects.

Systematic analysis of variants related to familial hypercholesterolemia in families with premature myocardial infarction.

Familial hypercholesterolemia (FH) is an oligogenic disorder characterized by markedly elevated low-density lipoprotein cholesterol (LDLC) levels. Variants in four genes have been reported to cause the classical autosomal-dominant form of the disease. FH is largely under-diagnosed in European countries. As FH increases the risk for coronary artery disease (CAD) and myocardial infarction (MI), it might be specifically overlooked in the large number of such patients. Here, we systematically examined the frequency of potential FH-causing variants by exome sequencing in 255 German patients with premature MI and a positive family history for CAD. We further performed co-segregation analyses in an average of 5.5 family members per MI patient. In total, we identified 11 potential disease-causing variants that co-segregate within the families, that is, 5% of patients with premature MI and positive CAD family history had FH. Eight variants were previously reported as disease-causing and three are novel (LDLR.c.811G>A p.(V271I)), PCSK9.c.610G>A (p.(D204N)) and STAP1.c.139A>G (p.(T47A))). Co-segregation analyses identified multiple additional family members carrying one of these FH variants and the clinical phenotype of either FH (n=2) or FH and premature CAD (n=15). However, exome sequencing also revealed that some variants in FH genes, which have been reported to cause FH, do not co-segregate with FH. The data reveal that a large proportion of FH patients escape the diagnosis, even when they have premature MI. Hence, systematic molecular-genetic screening for FH in such patients may reveal a substantial number of cases and thereby allow a timely LDLC-lowering in both FH/MI patients as well as their variant-carrying family members.

46,XY Gonadal Dysgenesis due to a Homozygous Mutation in Desert Hedgehog (DHH) Identified by Exome Sequencing.

BACKGROUND: 46,XY disorders of sex development (DSD) comprise a heterogeneous group of congenital conditions. Mutations in a variety of genes can affect gonadal development or androgen biosynthesis/action and thereby influence the development of the internal and external genital organs. OBJECTIVE: The objective of the study was to identify the genetic cause in two 46,XY sisters of a consanguineous family with DSD and gonadal tumor formation. METHODS: We used a next-generation sequencing approach by exome sequencing. Electrophysiological and high-resolution ultrasound examination of peripheral nerves as well as histopathological examination of the gonads were performed. RESULTS: We identified a novel homozygous R124Q mutation in the desert hedgehog gene (DHH), which alters a conserved residue among the three mammalian Hedgehog ligands sonic hedgehog, Indian hedgehog, and desert hedgehog. No other relevant mutations in DSD-related genes were encountered. The gonads of one patient showed partial gonadal dysgenesis with loss of Leydig cells in tubular areas with seminoma in situ and a hyperplasia of Leydig cell-like cells expressing CYP17A1 in more dysgenetic parts of the gonad. In addition, both patients suffer from a polyneuropathy. High-resolution ultrasound revealed a structural change of peripheral nerve structure that fits well to a minifascicle formation of peripheral nerves. CONCLUSION: Mutations in DHH play a role in 46,XY gonadal dysgenesis and are associated with seminoma formation and a neuropathy with minifascicle formation. Gonadal dysgenesis in these cases may be due to impairment of Sertoli cell-Leydig cell interaction during gonadal development.

A positive genotype-phenotype correlation in a large cohort of patients with Pseudohypoparathyroidism Type Ia and Pseudo-pseudohypoparathyroidism and 33 newly identified mutations in the GNAS gene.

Maternally inherited inactivating GNAS mutations are the most common cause of parathyroid hormone (PTH) resistance and Albright hereditary osteodystrophy (AHO) leading to pseudohypoparathyroidism type Ia (PHPIa) due to Gsα deficiency. Paternally inherited inactivating mutations lead to isolated AHO signs characterizing pseudo-pseudohypoparathyroidism (PPHP). Mutations are distributed throughout the Gsα coding exons of GNAS and there is a lack of genotype-phenotype correlation. In this study, we sequenced exon 1-13 of GNAS in a large cohort of PHPIa- and PPHP patients and identified 58 different mutations in 88 patients and 27 relatives. Thirty-three mutations including 15 missense mutations were newly discovered. Furthermore, we found three hot spots: a known hotspot (p.D190MfsX14), a second at codon 166 (p.R166C), and a third at the exon 5 acceptor splice site (c.435 + 1G>A), found in 15, 5, and 4 unrelated patients, respectively. Comparing the clinical features to the molecular genetic data, a significantly higher occurrence of subcutaneous calcifications in patients harboring truncating versus missense mutations was demonstrated. Thus, in the largest cohort of PHPIa patients described to date, we extend the spectrum of known GNAS mutations and hot spots and demonstrate for the first time a correlation between the genetic defects and the expression of a clinical AHO-feature.

Novel Insights into 46,XY Disorders of Sex Development due to NR5A1 Gene Mutation.

The differential diagnosis of 46,XY disorders of sex development (DSD) is based on the distinction between forms of gonadal dysgenesis and disorders of androgen biosynthesis and action. However, clinical and endocrine evaluations are often not conclusive. Here, we describe an adolescent female with hirsutism and hyperandrogenization at puberty. Her karyotype was 46,XY, and clinical investigation demonstrated clitoromegaly, but no uterine remnants were detected. Histology of the gonads revealed a testicular structure with a Sertoli-cell-only pattern. Endocrine evaluation showed hypergonadotropic hypogonadism, and the Sertoli cell markers inhibin B and anti-Müllerian hormone were also low. Several molecular genetic studies were initiated. While analyses of the androgen receptor gene, the SRD5A2 gene and HSD17B3 gene were uninformative, a novel p.L230R mutation was found in the NR5A1 gene. A mutant construct proved a severe dysfunction of this variant in functional analysis after recreation and transfection into HeLa cells. We conclude that the NR5A1 p.L230R mutation most likely leads to a spatial and time-dependent Leydig cell and Sertoli cell dysfunction during development not causing the classical gonadal dysgenesis phenotype. This case demonstrates that the current classification should be updated to encompass the overlapping phenotypes of some genetic conditions within 46,XY DSD.