Benedikt Reiz

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.

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.

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.

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.

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.

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.

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.

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.

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.