Whole-exome sequencing of 2,000 Danish individuals and the role of rare coding variants in type 2 diabetes

Dissecting heritability, environmental risk, and air pollution causal effects using > 50 million individuals in MarketScan

Large national-level electronic health record (EHR) datasets offer new opportunities for disentangling the role of genes and environment through deep phenotype information and approximate pedigree structures. Here we use the approximate geographical locations of patients as a proxy for spatially correlated community-level environmental risk factors. We develop a spatial mixed linear effect (SMILE) model that incorporates both genetics and environmental contribution. We extract EHR and geographical locations from 257,620 nuclear families and compile 1083 disease outcome measurements from the MarketScan dataset. We augment the EHR with publicly available environmental data, including levels of particulate matter 2.5 (PM2.5), nitrogen dioxide (NO2), climate, and sociodemographic data. We refine the estimates of genetic heritability and quantify community-level environmental contributions. We also use wind speed and direction as instrumental variables to assess the causal effects of air pollution. In total, we find PM2.5 or NO2 have statistically significant causal effects on 135 diseases, including respiratory, musculoskeletal, digestive, metabolic, and sleep disorders, where PM2.5 and NO2 tend to affect biologically distinct disease categories. These analyses showcase several robust strategies for jointly modeling genetic and environmental effects on disease risk using large EHR datasets and will benefit upcoming biobank studies in the era of precision medicine.

TINF2 is a major susceptibility gene in Danish patients with multiple primary melanoma

TINF2 encodes the TINF2 protein, which is a subunit in the shelterin complex critical for telomere regulation. Three recent studies have associated six truncating germline variants in TINF2 that have previously been associated with a cancer predisposition syndrome (CPS) caused by elongation of the telomeres. This has added TINF2 to the long telomere syndrome genes, together with other telomere maintenance genes such as ACD, POT1, TERF2IP, and TERT. We report a clinical study of 102 Danish patients with multiple primary melanoma (MPM) in which a germline truncating variant in TINF2 (p.(Arg265Ter)) was identified in four unrelated participants. The telomere lengths of three variant carriers were >90% percentile. In a routine diagnostic setting, the variant was identified in two more families, including an additional MPM patient and monozygotic twins with thyroid cancer and other cancer types. A total of 10 individuals from six independent families were confirmed carriers, all with cancer history, predominantly melanoma. Our findings suggest a major role of TINF2 in Danish patients with MPM. In addition to melanoma, other cancers in the six families include thyroid, renal, breast, and sarcoma, supporting a CPS in which melanoma, thyroid cancer, and sarcoma predominate. Further studies are needed to establish the full spectrum of associated cancer types and characterize lifetime cancer risk in carriers.

Whole exome sequencing of 28 families of Danish descent reveals novel candidate genes and pathways in developmental dysplasia of the hip

Developmental dysplasia of the hip (DDH) is a common condition involving instability of the hip with multifactorial etiology. Early diagnosis and treatment are critical as undetected DDH is an important cause of long-term hip complications. Better diagnostics may be achieved through genetic methods, especially for patients with positive family history. Several candidate genes have been reported but the exact molecular etiology of the disease is yet unknown. In the present study, we performed whole exome sequencing of DDH patients from 28 families with at least two affected first-degree relatives. Four genes previously not associated with DDH (METTL21B, DIS3L2, PPP6R2, and TM4SF19) were identified with the same variants shared among affected family members, in more than two families. Among known association genes, we found damaging variants in DACH1, MYH10, NOTCH2, TBX4, EVC2, OTOG, and SHC3. Mutational burden analysis across the families identified 322 candidate genes, and enriched pathways include the extracellular matrix, cytoskeleton, ion-binding, and detection of mechanical stimulus. Taken altogether, our data suggest a polygenic mode of inheritance for DDH, and we propose that an impaired transduction of the mechanical stimulus is involved in the etiopathological mechanism. Our findings refine our current understanding of candidate causal genes in DDH, and provide a foundation for downstream functional studies.

Disease Status at Diagnosis in Danish Children with α 1 -antitrypsin Deficiency

OBJECTIVES

The aim of this cross-sectional study was to assess the state of disease at the time of diagnosis in Danish children with α 1 -antitrypsin deficiency as Denmark has a high prevalence of ZZ-homozygosity.

METHODS

Children either heterozygous, compound heterozygous, or homozygous for Z- and S-variants in the SERPINA1 -gene were included. Clinical characteristics, SERPINA1 -genotype, and blood serum (S) concentrations were recorded concurrently with genetic testing. Serum liver marker concentrations were compared using T tests and Wilcoxon-Mann-Whitney tests. Generalized estimating equation (GEE) linear regression models, both univariable and multivariable adjusted for age and sex, were applied to identify correlations with serum α 1 -antitrypsin (S-AAT). The relationship between S-AAT concentration and genotype was assessed using logistic regression with GEE.

RESULTS

The study included 183 of 225 children genetically tested for alpha-1-antitrypsin deficiency (AATD). Of these, 36.6% were homozygous for the Z-variant. Of the heterozygotes, 89.7% had a ZM genotype and the remaining had either an MS genotype or were compound heterozygous. At diagnosis, ZZ-homozygous children had higher serum concentrations of liver enzymes and conjugated bilirubin, but lower concentrations of S-AAT compared with heterozygotes. Serum concentrations of conjugated bilirubin and liver enzymes were negatively associated with S-AAT. Children under 6 months of age had higher total S-bilirubin concentrations than children over 6 months of age.

CONCLUSIONS

A low S-AAT concentration is a strong indicator of homozygosity, and homozygous children have higher enzymatic and cholestatic parameters compared with heterozygous children at diagnosis. This underlines the importance of measuring the S-AAT concentration in children with prolonged neonatal jaundice.

A combined polygenic score of 21,293 rare and 22 common variants improves diabetes diagnosis based on hemoglobin A1C levels

Polygenic scores (PGSs) combine the effects of common genetic variants1,2 to predict risk or treatment strategies for complex diseases3-7. Adding rare variation to PGSs has largely unknown benefits and is methodically challenging. Here, we developed a method for constructing rare variant PGSs and applied it to calculate genetically modified hemoglobin A1C thresholds for type 2 diabetes (T2D) diagnosis7-10. The resultant rare variant PGS is highly polygenic (21,293 variants across 154 genes), depends on ultra-rare variants (72.7% observed in fewer than three people) and identifies significantly more undiagnosed T2D cases than expected by chance (odds ratio = 2.71; P = 1.51 × 10-6). A PGS combining common and rare variants is expected to identify 4.9 million misdiagnosed T2D cases in the United States-nearly 1.5-fold more than the common variant PGS alone. These results provide a method for constructing complex trait PGSs from rare variants and suggest that rare variants will augment common variants in precision medicine approaches for common disease.

Rare loss of function variants in the hepatokine gene INHBE protect from abdominal obesity

Identifying genetic variants associated with lower waist-to-hip ratio can reveal new therapeutic targets for abdominal obesity. We use exome sequences from 362,679 individuals to identify genes associated with waist-to-hip ratio adjusted for BMI (WHRadjBMI), a surrogate for abdominal fat that is causally linked to type 2 diabetes and coronary heart disease. Predicted loss of function (pLOF) variants in INHBE associate with lower WHRadjBMI and this association replicates in data from AMP-T2D-GENES. INHBE encodes a secreted protein, the hepatokine activin E. In vitro characterization of the most common INHBE pLOF variant in our study, indicates an in-frame deletion resulting in a 90% reduction in secreted protein levels. We detect associations with lower WHRadjBMI for variants in ACVR1C, encoding an activin receptor, further highlighting the involvement of activins in regulating fat distribution. These findings highlight activin E as a potential therapeutic target for abdominal obesity, a phenotype linked to cardiometabolic disease.

The Scaffold Protein PICK1 as a Target in Chronic Pain

Well-tolerated and effective drugs for treating chronic pain conditions are urgently needed. Most chronic pain patients are not effectively relieved from their pain and suffer from debilitating drug side effects. This has not only drastic negative consequences for the patients’ quality of life, but also constitute an enormous burden on society. It is therefore of great interest to explore new potent targets for effective pain treatment with fewer side effects and without addiction liability. A critical component of chronic pain conditions is central sensitization, which involves the reorganization and strengthening of synaptic transmission within nociceptive pathways. Such changes are considered as maladaptive and depend on changes in the surface expression and signaling of AMPA-type glutamate receptors (AMPARs). The PDZ-domain scaffold protein PICK1 binds the AMPARs and has been suggested to play a key role in these maladaptive changes. In the present paper, we review the regulation of AMPARs by PICK1 and its relation to pain pathology. Moreover, we highlight other pain-relevant PICK1 interactions, and we evaluate various compounds that target PICK1 and have been successfully tested in pain models. Finally, we evaluate the potential on-target side effects of interfering with the action of PICK1 action in CNS and beyond. We conclude that PICK1 constitutes a valid drug target for the treatment of inflammatory and neuropathic pain conditions without the side effects and abuse liability associated with current pain medication.

Whole-Exome Sequencing Implicates Neuronal Calcium Channel with Familial Atrial Fibrillation

Background: Atrial Fibrillation (AF) is the most prevalent sustained cardiac arrhythmia, responsible for considerable morbidity and mortality. The heterogenic and complex pathogenesis of AF remains poorly understood, which contributes to the current limitation in effective treatments. We aimed to identify rare genetic variants associated with AF in patients with familial AF.

Methods and results: We performed whole exome sequencing in a large family with familial AF and identified a rare variant in the gene CACNA1A c.5053G > A which co-segregated with AF. The gene encodes for the protein variants Ca_V2.1-V1686M, and is important in neuronal function. Functional characterization of the CACNA1A, using patch-clamp recordings on transiently transfected mammalian cells, revealed a modest loss-of-function of Ca_V2.1-V1686M.

Conclusion: We identified a rare loss-of-function variant associated with AF in a gene previously linked with neuronal function. The results allude to a novel link between dysfunction of an ion channel previously associated with neuronal functions and increased risk of developing AF.

Coding variants identified in diabetic patients alter PICK1 BAR domain function in insulin granule biogenesis

Bin/amphiphysin/Rvs (BAR) domains are positively charged crescent-shaped modules that mediate curvature of negatively charged lipid membranes during remodeling processes. The BAR domain proteins PICK1, ICA69, and the arfaptins have recently been demonstrated to coordinate the budding and formation of immature secretory granules (ISGs) at the trans-Golgi network. Here, we identify 4 coding variants in the PICK1 gene from a whole-exome screening of Danish patients with diabetes that each involve a change in positively charged residues in the PICK1 BAR domain. All 4 coding variants failed to rescue insulin content in INS-1E cells upon knock down of endogenous PICK1. Moreover, 2 variants showed dominant-negative properties. In vitro assays addressing BAR domain function suggested that the coding variants compromised BAR domain function but increased the capacity to cause fission of liposomes. Live confocal microscopy and super-resolution microscopy further revealed that PICK1 resides transiently on ISGs before egress via vesicular budding events. Interestingly, this egress of PICK1 was accelerated in the coding variants. We propose that PICK1 assists in or complements the removal of excess membrane and generic membrane trafficking proteins, and possibly also insulin, from ISGs during the maturation process; and that the coding variants may cause premature budding, possibly explaining their dominant-negative function.

Lipolysis drives expression of the constitutively active receptor GPR3 to induce adipose thermogenesis

Thermogenic adipocytes possess a therapeutically appealing, energy-expending capacity, which is canonically cold-induced by ligand-dependent activation of β-adrenergic G protein-coupled receptors (GPCRs). Here, we uncover an alternate paradigm of GPCR-mediated adipose thermogenesis through the constitutively active receptor, GPR3. We show that the N terminus of GPR3 confers intrinsic signaling activity, resulting in continuous Gs-coupling and cAMP production without an exogenous ligand. Thus, transcriptional induction of Gpr3 represents the regulatory parallel to ligand-binding of conventional GPCRs. Consequently, increasing Gpr3 expression in thermogenic adipocytes is alone sufficient to drive energy expenditure and counteract metabolic disease in mice. Gpr3 transcription is cold-stimulated by a lipolytic signal, and dietary fat potentiates GPR3-dependent thermogenesis to amplify the response to caloric excess. Moreover, we find GPR3 to be an essential, adrenergic-independent regulator of human brown adipocytes. Taken together, our findings reveal a noncanonical mechanism of GPCR control and thermogenic activation through the lipolysis-induced expression of constitutively active GPR3.

KATK: Fast genotyping of rare variants directly from unmapped sequencing reads

KATK is a fast and accurate software tool for calling variants directly from raw next-generation sequencing reads. It uses predefined k-mers to retrieve only the reads of interest from the FASTQ file and calls genotypes by aligning retrieved reads locally. KATK does not use data about known polymorphisms and has NC (no call) as the default genotype. The reference or variant allele is called only if there is sufficient evidence for their presence in data. Thus it is not biased against rare variants or de-novo mutations. With simulated datasets, we achieved a false-negative rate of 0.23% (sensitivity 99.77%) and a false discovery rate of 0.19%. Calling all human exonic regions with KATK requires 1-2 h, depending on sequencing coverage.

PHACTR1 genetic variability is not critical in small vessel ischemic disease patients and PcomA recruitment in C57BL/6J mice

Recently, several genome-wide association studies identified PHACTR1 as key locus for five diverse vascular disorders: coronary artery disease, migraine, fibromuscular dysplasia, cervical artery dissection and hypertension. Although these represent significant risk factors or comorbidities for ischemic stroke, PHACTR1 role in brain small vessel ischemic disease and ischemic stroke most important survival mechanism, such as the recruitment of brain collateral arteries like posterior communicating arteries (PcomAs), remains unknown. Therefore, we applied exome and genome sequencing in a multi-ethnic cohort of 180 early-onset independent familial and apparently sporadic brain small vessel ischemic disease and CADASIL-like Caucasian patients from US, Portugal, Finland, Serbia and Turkey and in 2 C57BL/6J stroke mouse models (bilateral common carotid artery stenosis [BCCAS] and middle cerebral artery occlusion [MCAO]), characterized by different degrees of PcomAs patency. We report 3 very rare coding variants in the small vessel ischemic disease-CADASIL-like cohort (p.Glu198Gln, p.Arg204Gly, p.Val251Leu) and a stop-gain mutation (p.Gln273*) in one MCAO mouse. These coding variants do not cluster in PHACTR1 known pathogenic domains and are not likely to play a critical role in small vessel ischemic disease or brain collateral circulation. We also exclude the possibility that copy number variants (CNVs) or a variant enrichment in Phactr1 may be associated with PcomA recruitment in BCCAS mice or linked to diverse vascular traits (cerebral blood flow pre-surgery, PcomA size, leptomeningeal microcollateral length and junction density during brain hypoperfusion) in C57BL/6J mice, respectively. Genetic variability in PHACTR1 is not likely to be a common susceptibility factor influencing small vessel ischemic disease in patients and PcomA recruitment in C57BL/6J mice. Nonetheless, rare variants in PHACTR1 RPEL domains may influence the stroke outcome and are worth investigating in a larger cohort of small vessel ischemic disease patients, different ischemic stroke subtypes and with functional studies.

Clinical genetic diagnostics in Danish autosomal dominant polycystic kidney disease patients reveal possible founder variants

BACKGROUND

Autosomal dominant polycystic kidney disease (ADPKD) is the most common heritable kidney disease. ADPKD leads to cysts, kidney enlargement and end-stage renal disease. ADPKD is mainly caused by variants in PKD1 and PKD2, with truncating PKD1 variants causing the most severe phenotype. This study aimed to characterize variants in Danish patients referred for screening of genes related to cystic kidney disease.

METHODS

147 families were analysed for variants in PKD1, PKD2 and GANAB using next generation sequencing and multiplex ligation-dependent probe amplification. If a variant was identified, relatives were analysed for the specific variant using Sanger sequencing.

RESULTS

A pathogenic or possibly pathogenic variant was identified in 87% (103/118) of patients suspected to suffer from ADPKD, according to the requisition form. In total, 112 pathogenic or possibly pathogenic variants were observed, of which 94 were unique; 74 (79%) in PKD1 and 20 (21%) in PKD2, while 41 variants were novel. No variants in GANAB were observed. Ten recurrent variants were observed in 26 (26%) families. These were either PKD2 variants (N = 6) or non-truncating PKD1 variants (N = 4). Five of these were likely founder variants.

CONCLUSIONS

The distribution of pathogenic or possibly pathogenic variants in the Danish ADPKD population is similar to that in other populations, except that recurrent truncating PKD1 variants appear to be rare, i.e. founder variants tend to be variant types associated with a mild phenotype. Patients with a mild phenotype may remain undiagnosed, consequently the frequency of founder variants and prevalence of ADPKD may be underestimated.

Germline RBBP8 variants associated with early-onset breast cancer compromise replication fork stability

Haploinsufficiency of factors governing genome stability underlies hereditary breast and ovarian cancer. One significant pathway that is disabled as a result is homologous recombination repair (HRR). With the aim of identifying new candidate genes, we examined early-onset breast cancer patients negative for BRCA1 and BRCA2 pathogenic variants. Here, we focused on CtIP (RBBP8 gene), which mediates HRR through the end resection of DNA double-strand breaks (DSBs). Notably, these patients exhibited a number of rare germline RBBP8 variants. Functional analysis revealed that these variants did not affect DNA DSB end resection efficiency. However, expression of a subset of variants led to deleterious nucleolytic degradation of stalled DNA replication forks in a manner similar to that of cells lacking BRCA1 or BRCA2. In contrast to BRCA1 and BRCA2, CtIP deficiency promoted the helicase-driven destabilization of RAD51 nucleofilaments at damaged DNA replication forks. Taken together, our work identifies CtIP as a critical regulator of DNA replication fork integrity, which, when compromised, may predispose to the development of early-onset breast cancer.

Systems genetics analysis identifies calcium-signaling defects as novel cause of congenital heart disease

BACKGROUND

Congenital heart disease (CHD) occurs in almost 1% of newborn children and is considered a multifactorial disorder. CHD may segregate in families due to significant contribution of genetic factors in the disease etiology. The aim of the study was to identify pathophysiological mechanisms in families segregating CHD.

METHODS

We used whole exome sequencing to identify rare genetic variants in ninety consenting participants from 32 Danish families with recurrent CHD. We applied a systems biology approach to identify developmental mechanisms influenced by accumulation of rare variants. We used an independent cohort of 714 CHD cases and 4922 controls for replication and performed functional investigations using zebrafish as in vivo model.

RESULTS

We identified 1785 genes, in which rare alleles were shared between affected individuals within a family. These genes were enriched for known cardiac developmental genes, and 218 of these genes were mutated in more than one family. Our analysis revealed a functional cluster, enriched for proteins with a known participation in calcium signaling. Replication in an independent cohort confirmed increased mutation burden of calcium-signaling genes in CHD patients. Functional investigation of zebrafish orthologues of ITPR1, PLCB2, and ADCY2 verified a role in cardiac development and suggests a combinatorial effect of inactivation of these genes.

CONCLUSIONS

The study identifies abnormal calcium signaling as a novel pathophysiological mechanism in human CHD and confirms the complex genetic architecture underlying CHD.

A rare missense variant in APC interrupts splicing and causes AFAP in two Danish families

Background

We report the first case of a missense variant in the APC gene that interrupts splicing by creating a new cryptic acceptor site. The variant, c.289G>A, p.(Gly97Arg), is located in exon 3, and qualitative and semi-quantitative RNA splicing analysis reveal that the variant results in skipping of the last 70 nucleotides of the exon, which leads to the introduction of a frameshift and a premature stop codon.

Case presentation

The variant was detected in two, apparently unrelated, Danish families with an accumulation of colorectal cancers, colonic adenomas and other cancers. The families both have an attenuated familial adenomatous polyposis phenotype, which is consistent with the association of pathogenic variants in the 5′ end of the gene.

One variant-carrier also had Caroli Disease and a Caroli Disease associated hepatic mucinous cystadenocarcinoma. This is the first description of a person with both Caroli Disease and a pathogenic APC variant, and although the APC variant is not known to be connected to the development of the hepatic malformations in Caroli Disease, it remains unclear whether the variant could have contributed to the carcinogenesis of the liver tumour.

Conclusions

Based on functional and co-segregation data we classify the APC c.289G>A, p.(Gly97Arg) variant as pathogenic (class 5). Our findings emphasize the importance of a functional evaluation of missense variants although located far from the exon-intron boundaries.

A mutation map for human glycoside hydrolase genes

Glycoside hydrolases (GHs) are found in all domains of life, and at least 87 distinct genes encoding proteins related to GHs are found in the human genome. GHs serve diverse functions from digestion of dietary polysaccharides to breakdown of intracellular oligosaccharides, glycoproteins, proteoglycans and glycolipids. Congenital disorders of GHs (CDGHs) represent more than 30 rare diseases caused by mutations in one of the GH genes. We previously used whole-exome sequencing of a homogenous Danish population of almost 2000 individuals to probe the incidence of deleterious mutations in the human glycosyltransferases (GTs) and developed a mutation map of human GT genes (GlyMAP-I). While deleterious disease-causing mutations in the GT genes were very rare, and in many cases lethal, we predicted deleterious mutations in GH genes to be less rare and less severe given the higher incidence of CDGHs reported worldwide. To probe the incidence of GH mutations, we constructed a mutation map of human GH-related genes (GlyMAP-II) using the Danish WES data, and correlating this with reported disease-causing mutations confirmed the higher prevalence of disease-causing mutations in several GH genes compared to GT genes. We identified 76 novel nonsynonymous single-nucleotide variations (nsSNVs) in 32 GH genes that have not been associated with a CDGH phenotype, and we experimentally validated two novel potentially damaging nsSNVs in the congenital sucrase-isomaltase deficiency gene, SI. Our study provides a global view of human GH genes and disease-causing mutations and serves as a discovery tool for novel damaging nsSNVs in CDGHs.

A Novel Loss-of-Function Variant in the Chloride Ion Channel Gene Clcn2 Associates with Atrial Fibrillation

Atrial Fibrillation (AF) is the most common cardiac arrhythmia. Its pathogenesis is complex and poorly understood. Whole exome sequencing of Danish families with AF revealed a novel four nucleotide deletion c.1041_1044del in CLCN2 shared by affected individuals. We aimed to investigate the role of genetic variation of CLCN2 encoding the inwardly rectifying chloride channel ClC-2 as a risk factor for the development of familiar AF. The effect of the CLCN2 variant was evaluated by electrophysiological recordings on transiently transfected cells. We used quantitative PCR to assess CLCN2 mRNA expression levels in human atrial and ventricular tissue samples. The nucleotide deletion CLCN2 c.1041_1044del results in a frame-shift and premature stop codon. The truncated ClC-2 p.V347fs channel does not conduct current. Co-expression with wild-type ClC-2, imitating the heterozygote state of the patients, resulted in a 50% reduction in macroscopic current, suggesting an inability of truncated ClC-2 protein to form channel complexes with wild type channel subunits. Quantitative PCR experiments using human heart tissue from healthy donors demonstrated that CLCN2 is expressed across all four heart chambers. Our genetic and functional data points to a possible link between loss of ClC-2 function and an increased risk of developing AF.

Targeted exon sequencing in deceased schizophrenia patients in Denmark

Schizophrenia patients have higher mortality rates and lower life expectancy than the general population. However, forensic investigations of their deaths often fail to determine the cause of death, hindering prevention. As schizophrenia is a highly heritable condition and given recent advances in our understanding of the genetics of schizophrenia, it is now possible to investigate how genetic factors may contribute to mortality. We made use of findings from genome-wide association studies (GWAS) to design a targeted panel (PsychPlex) for sequencing of exons of 451 genes near index single nucleotide polymorphisms (SNPs) identified with GWAS. We sequenced the DNA of 95 deceased schizophrenia patients included in SURVIVE, a prospective, autopsy-based study of mentally ill persons in Denmark. We compared the allele frequencies of 1039 SNPs in these cases with the frequencies of 2000 Danes without psychiatric diseases and calculated their deleteriousness (CADD) scores. For 81 SNPs highly associated with schizophrenia and CADD scores above 15, expression profiles in the Genotype-Tissue Expression (GTEx) Project indicated that these variants were in exons, whose expressions are increased in several types of brain tissues, particularly in the cerebellum. Molecular pathway analysis indicated the involvement of 163 different pathways. As for rare SNP variants, most variants were scored as either benign or likely benign with an average of 17 variants of unknown significance per individual and no pathogenic variant. Our results highlight the potential of DNA sequencing of an exon panel to discover genetic factors that may be involved in the development of schizophrenia.

Reappraisal of variants previously linked with sudden infant death syndrome: results from three population-based cohorts

We aimed to investigate the pathogenicity of cardiac ion channel variants previously associated with SIDS. We reviewed SIDS-associated variants previously reported in databases and the literature in three large population-based cohorts; The ExAC database, the Inter99 study, and the UK Biobank (UKBB). Variants were classified according to the American College of Medical Genetics and Genomics (ACMG) guidelines. Of the 92 SIDS-associated variants, 59 (64%) were present in ExAC, 18 (20%) in Inter99, and 24 (26%) in UKBB. Using the Inter99 cohort, we found no difference in J-point amplitude and QTc-interval between carriers and non-carriers for 14/18 variants. There was no difference in the risk of syncope (P = 0.32), malignant ventricular arrhythmia (P = 0.96), and all-cause mortality (P = 0.59) between carriers and non-carriers. The ACMG guidelines reclassified 75% of all variants as variant-of-uncertain significance, likely benign, and benign. We identified ~2/3 of variants previously associated with SIDS and found no significant associations with electrocardiographic traits, syncope, malignant ventricular arrhythmia, or all-cause mortality. These data indicate that many of these variants are not highly penetrant, monogenic causes of SIDS and underline the importance of frequent reappraisal of genetic variants to avoid future misdiagnosis.

T2D Risk Genes: Exome Sequencing Goes Straight to the Source

Genome-wide association studies have identified hundreds of genomic variants associated with human T2D risk, but translating such findings to clinically useful information has proved challenging. A new study in Nature (Flannick et al., 2019) breaks this gridlock, using direct exome sequencing to identify functional coding variants, providing critical complementary gene-level information.

Exome sequencing of 20,791 cases of type 2 diabetes and 24,440 controls

Protein-coding genetic variants that strongly affect disease risk can yield relevant clues to disease pathogenesis. Here we report exome-sequencing analyses of 20,791 individuals with type 2 diabetes (T2D) and 24,440 non-diabetic control participants from 5 ancestries. We identify gene-level associations of rare variants (with minor allele frequencies of less than 0.5%) in 4 genes at exome-wide significance, including a series of more than 30 SLC30A8 alleles that conveys protection against T2D, and in 12 gene sets, including those corresponding to T2D drug targets (P = 6.1 × 10-3) and candidate genes from knockout mice (P = 5.2 × 10-3). Within our study, the strongest T2D gene-level signals for rare variants explain at most 25% of the heritability of the strongest common single-variant signals, and the gene-level effect sizes of the rare variants that we observed in established T2D drug targets will require 75,000-185,000 sequenced cases to achieve exome-wide significance. We propose a method to interpret these modest rare-variant associations and to incorporate these associations into future target or gene prioritization efforts.

The Contribution of Low-Frequency and Rare Coding Variation to Susceptibility to Type 2 Diabetes

PURPOSE OF REVIEW

Soon after the first genome-wide association study (GWAS) for type 2 diabetes (T2D) was published, it was hypothesized that rare and low-frequency variants might explain a substantial proportion of disease risk. Rare coding variants in particular were emphasized given their large expected role in disease. This review summarizes the extent to which recent T2D genetic studies provide evidence for or against this hypothesis.

RECENT FINDINGS

Following a comprehensive study of T2D genetic architecture using three sequencing and genotyping technologies, four even larger studies have provided a yet higher resolution view of the role of rare and low-frequency coding variation in T2D susceptibility. Empirical evidence strongly suggests that common regulatory variants are the dominant contributor to T2D heritability. However, rare coding variants may nonetheless be pervasive across T2D-relevant genes. A strategy using common variants to map disease genes, and rare coding variants to link molecular gene perturbations to cellular and phenotypic effects, may be an effective means to investigate T2D pathogenesis and potential new therapies.

Whole Genome Sequence, Variant Discovery and Annotation in Mapuche-Huilliche Native South Americans

Whole human genome sequencing initiatives help us understand population history and the basis of genetic diseases. Current data mostly focuses on Old World populations, and the information of the genomic structure of Native Americans, especially those from the Southern Cone is scant. Here we present annotation and variant discovery from high-quality complete genome sequences of a cohort of 11 Mapuche-Huilliche individuals (HUI) from Southern Chile. We found approximately 3.1 × 10⁶ single nucleotide variants (SNVs) per individual and identified 403,383 (6.9%) of novel SNVs events. Analyses of large-scale genomic events detected 680 copy number variants (CNVs) and 4,514 structural variants (SVs), including 398 and 1,910 novel events, respectively. Global ancestry composition of HUI genomes revealed that the cohort represents a sample from a marginally admixed population from the Southern Cone, whose main genetic component derives from Native American ancestors. Additionally, we found that HUI genomes contain variants in genes associated with 5 of the 6 leading causes of noncommunicable diseases in Chile, which may have an impact on the risk of prevalent diseases in Chilean and Amerindian populations. Our data represents a useful resource that can contribute to population-based studies and for the design of early diagnostics or prevention tools for Native and admixed Latin American populations.

A study of familial Char syndrome involving the TFAP2B gene with a focus on facial shape characteristics

In this case study, we investigate a child presenting with patent ductus arteriosus, short philtrum, duck-bill lips, strabismus, a flat nasal bridge, a broad forehead, low-set ears, hypertelorism, up-slanting palpebral fissures, almond-shaped eyes, and hypodontia, all leading to the clinical diagnosis of Char syndrome. Genetic analysis showed heterozygosity for the novel variant c.851T>C, p. Leu284Ser in the TFAP2B gene. Family analysis suggested that at least 20 members, extending six generations back, were affected. All 10 members available for genetic testing were heterozygous for the novel pathogenic variant. Qualitative analysis of the facial dysmorphology in the proband and three of the affected family members using three-dimensional surface scanning showed that the major deviations were observed in the forehead/eyebrow, nose, upper lip, and chin regions with, for example, a flattened nose and reduced height of the upper lip and the face. Furthermore, it is suggested that Char syndrome is associated with disturbances of tooth formation and eruption.

Rare truncating variants in the sarcomeric protein titin associate with familial and early-onset atrial fibrillation

A family history of atrial fibrillation constitutes a substantial risk of developing the disease, however, the pathogenesis of this complex disease is poorly understood. We perform whole-exome sequencing on 24 families with at least three family members diagnosed with atrial fibrillation (AF) and find that titin-truncating variants (TTNtv) are significantly enriched in these patients (P = 1.76 × 10⁻⁶). This finding is replicated in an independent cohort of early-onset lone AF patients (n = 399; odds ratio = 36.8; P = 4.13 × 10⁻⁶). A CRISPR/Cas9 modified zebrafish carrying a truncating variant of titin is used to investigate TTNtv effect in atrial development. We observe compromised assembly of the sarcomere in both atria and ventricle, longer PR interval, and heterozygous adult zebrafish have a higher degree of fibrosis in the atria, indicating that TTNtv are important risk factors for AF. This aligns with the early onset of the disease and adds an important dimension to the understanding of the molecular predisposition for AF.

Common genetic variants in structural proteins contribute to risk of atrial fibrillation (AF). Here, using whole-exome sequencing, the authors identify rare truncating variants in TTN that associate with familial and early-onset AF and show defects in cardiac sarcomere assembly in ttn.2-mutant zebrafish.

Imputation-Aware Tag SNP Selection To Improve Power for Large-Scale, Multi-ethnic Association Studies

The emergence of very large cohorts in genomic research has facilitated a focus on genotype-imputation strategies to power rare variant association. These strategies have benefited from improvements in imputation methods and association tests, however little attention has been paid to ways in which array design can increase rare variant association power. Therefore, we developed a novel framework to select tag SNPs using the reference panel of 26 populations from Phase 3 of the 1000 Genomes Project. We evaluate tag SNP performance via mean imputed r2 at untyped sites using leave-one-out internal validation and standard imputation methods, rather than pairwise linkage disequilibrium. Moving beyond pairwise metrics allows us to account for haplotype diversity across the genome for improve imputation accuracy and demonstrates population-specific biases from pairwise estimates. We also examine array design strategies that contrast multi-ethnic cohorts vs. single populations, and show a boost in performance for the former can be obtained by prioritizing tag SNPs that contribute information across multiple populations simultaneously. Using our framework, we demonstrate increased imputation accuracy for rare variants (frequency < 1%) by 0.5-3.1% for an array of one million sites and 0.7-7.1% for an array of 500,000 sites, depending on the population. Finally, we show how recent explosive growth in non-African populations means tag SNPs capture on average 30% fewer other variants than in African populations. The unified framework presented here will enable investigators to make informed decisions for the design of new arrays, and help empower the next phase of rare variant association for global health.

Genetic and Phenotypic Characterization of Community Hospital Patients With QT Prolongation

Background Congenital long- QT syndrome ( LQTS ) is a genetic disorder characterized by prolongation of the corrected QT interval ( QT c) on an ECG . The aim of the present study was to estimate the prevalence of pathogenic and likely pathogenic sequence variants in patients who had at least 1 ECG with a QT c ≥500 ms. Methods and Results Telemark Hospital Trust is a community hospital within the Norwegian national health system, serving ≈173 000 inhabitants. We searched the ECG database at Telemark Hospital Trust, Norway, from January 2004 to December 2014, and identified 1531 patients with at least 1 ECG with a QT c ≥500 ms. At the time of inclusion in this study (2015), 766 patients were alive. A total of 733 patients were invited to participate, and 475 accepted. The 17 genes that have been reported to cause monogenic LQTS were sequenced among the patients. Pro- QT c score was calculated for each patient. A molecular genetic cause of LQTS was detected in 31 (6.5%) of 475 patients. These patients had a lower pro- QT c score than those without pathogenic or likely pathogenic variants (1.7±1.0 versus 2.8±1.6; P<0.001). Conclusions Compared with the general population, hospitalized patients with a QT c ≥500 ms in at least 1 ECG recording had an increased likelihood for pathogenic and likely pathogenic variants in LQTS genes. We recommend increased awareness of the possibility of LQTS in patients with at least 1 ECG with a QT c ≥500 ms.

Identification of Novel Candidate Markers of Type 2 Diabetes and Obesity in Russia by Exome Sequencing with a Limited Sample Size

Type 2 diabetes (T2D) and obesity are common chronic disorders with multifactorial etiology. In our study, we performed an exome sequencing analysis of 110 patients of Russian ethnicity together with a multi-perspective approach based on biologically meaningful filtering criteria to detect novel candidate variants and loci for T2D and obesity. We have identified several known single nucleotide polymorphisms (SNPs) as markers for obesity (rs11960429), T2D (rs9379084, rs1126930), and body mass index (BMI) (rs11553746, rs1956549 and rs7195386) (p < 0.05). We show that a method based on scoring of case-specific variants together with selection of protein-altering variants can allow for the interrogation of novel and known candidate markers of T2D and obesity in small samples. Using this method, we identified rs328 in LPL (p = 0.023), rs11863726 in HBQ1 (p = 8 × 10⁻⁵), rs112984085 in VAV3 (p = 4.8 × 10⁻⁴) for T2D and obesity, rs6271 in DBH (p = 0.043), rs62618693 in QSER1 (p = 0.021), rs61758785 in RAD51B (p = 1.7 × 10⁻⁴), rs34042554 in PCDHA1 (p = 1 × 10⁻⁴), and rs144183813 in PLEKHA5 (p = 1.7 × 10⁻⁴) for obesity; and rs9379084 in RREB1 (p = 0.042), rs2233984 in C6orf15 (p = 0.030), rs61737764 in ITGB6 (p = 0.035), rs17801742 in COL2A1 (p = 8.5 × 10⁻⁵), and rs685523 in ADAMTS13 (p = 1 × 10⁻⁶) for T2D as important susceptibility loci in Russian population. Our results demonstrate the effectiveness of whole exome sequencing (WES) technologies for searching for novel markers of multifactorial diseases in cohorts of limited size in poorly studied populations.

DNAJC3 mutation in Thai familial type 2 diabetes mellitus

Type 2 diabetes mellitus (T2D) is a heterogeneous disease, with certain cases presenting an autosomal dominant type. The rare coding variants of disease‑causing genes in T2D remain mostly unclear. The present study aimed to identify the disease‑causing gene conducting whole exome sequencing in a Thai T2D family with an autosomal dominant transmission of T2D with no evidence of mutations in known maturity‑onset diabetes of the young (MODY) genes. Candidate variants were selected according to certain criteria of mutation prediction programs, followed by segregation analysis with diabetes in the family. The results demonstrated that, of the 68,817 variants obtained, 122 were considered as candidate variants subsequent to the filtering processes. Genotyping of these variants revealed that DnaJ homolog subfamily C member 3 (DNAJC3) p.H238N segregated with diabetes in the family. This mutation was also identified in another proband from the autosomal dominant T2D family without mutation in known MODY genes and was segregated with diabetes. This variant was also identified in 14/1,000 older‑onset T2D patients [minor allele frequency (MAF)=0.007], 2/500 non‑diabetic controls (MAF=0.002) and 3 prediabetic individuals who were previously classified as non‑diabetic controls. In silico mutagenesis and protein modeling of p.H238N revealed changes of the polar contacts across the tetratricopeptide repeat (TPR) motif and TPR subdomains, which may affect the protein tertiary structure. Furthermore, the expression of DNAJC3 H238N protein was 0.68±0.08 fold (P<0.05) lower when compared with that of the wild‑type, possibly due to protein instability. Thus, DNAJC3 p.H238N is likely to be a variant causing diabetes.

Whole-Exome Sequencing in Age-Related Macular Degeneration Identifies Rare Variants in COL8A1, a Component of Bruch's Membrane

PURPOSE

Genome-wide association studies and targeted sequencing studies of candidate genes have identified common and rare variants that are associated with age-related macular degeneration (AMD). Whole-exome sequencing (WES) studies allow a more comprehensive analysis of rare coding variants across all genes of the genome and will contribute to a better understanding of the underlying disease mechanisms. To date, the number of WES studies in AMD case-control cohorts remains scarce and sample sizes are limited. To scrutinize the role of rare protein-altering variants in AMD cause, we performed the largest WES study in AMD to date in a large European cohort consisting of 1125 AMD patients and 1361 control participants.

DESIGN

Genome-wide case-control association study of WES data.

PARTICIPANTS

One thousand one hundred twenty-five AMD patients and 1361 control participants.

METHODS

A single variant association test of WES data was performed to detect variants that are associated individually with AMD. The cumulative effect of multiple rare variants with 1 gene was analyzed using a gene-based CMC burden test. Immunohistochemistry was performed to determine the localization of the Col8a1 protein in mouse eyes.

MAIN OUTCOME MEASURES

Genetic variants associated with AMD.

RESULTS

We detected significantly more rare protein-altering variants in the COL8A1 gene in patients (22/2250 alleles [1.0%]) than in control participants (11/2722 alleles [0.4%]; P = 7.07×10^-5). The association of rare variants in the COL8A1 gene is independent of the common intergenic variant (rs140647181) near the COL8A1 gene previously associated with AMD. We demonstrated that the Col8a1 protein localizes at Bruch's membrane.

CONCLUSIONS

This study supported a role for protein-altering variants in the COL8A1 gene in AMD pathogenesis. We demonstrated the presence of Col8a1 in Bruch's membrane, further supporting the role of COL8A1 variants in AMD pathogenesis. Protein-altering variants in COL8A1 may alter the integrity of Bruch's membrane, contributing to the accumulation of drusen and the development of AMD.

Loss-of-function variants in ADCY3 increase risk of obesity and type 2 diabetes

We have identified a variant in ADCY3 (encoding adenylate cyclase 3) associated with markedly increased risk of obesity and type 2 diabetes in the Greenlandic population. The variant disrupts a splice acceptor site, and carriers have decreased ADCY3 RNA expression. Additionally, we observe an enrichment of rare ADCY3 loss-of-function variants among individuals with type 2 diabetes in trans-ancestry cohorts. These findings provide new information on disease etiology relevant for future treatment strategies.

Targeted next generation sequencing in a young population with suspected inherited malignant cardiac arrhythmias

Aborted sudden cardiac death in the young often is due to inherited heart disease. However, the clinical phenotype in these patients is not always evident. The aim of this study was to identify pathogenic molecular genetic variants in a population with suspected inherited cardiac arrhythmias. Eligible patients were admitted to Aarhus University Hospital, Denmark during the period 1999–2013 with arrhythmias assumed caused by a hereditary heart disease, and in whom no genotype had been established. We used the Danish national pacemaker and ICD registry to identify this cohort. One third (24/80) of the study population had first-line genetic testing with a targeted next-generation sequencing (NGS) panel, and two-third (56/80) of the study population had second-line genetic testing with NGS where prior Sanger sequencing did not reveal a causative variant. Variants were assessed according to the American College of Medical Genetics and Genomics (ACMG) guidelines. We included 80 patients. Median age (IQR) was 38 (28–43) years, 54 (68%) were males. First-line genetic testing identified a genetic variant in 33% (8/24) of the cases and second-line genetic testing revealed a variant in 20% (11/56) of the cases. Eleven variants were considered pathogenic, three likely pathogenic and 10 were variants of unknown significance (VUS). Seventeen variants were very rare with a minor allele frequency (MAF) ≤0.02% in all population databases used in the study. Molecular genetic testing of patients with suspected inherited cardiac arrhythmias with NGS identifies a molecular-genetic cause in a significant proportion of patients.

Hidden genetic variation shapes the structure of functional elements in Drosophila

Mutations that add, subtract, rearrange, or otherwise refashion genome structure often affect phenotypes, although the fragmented nature of most contemporary assemblies obscures them. To discover such mutations, we assembled the first new reference-quality genome of Drosophila melanogaster since its initial sequencing. By comparing this new genome to the existing D. melanogaster assembly, we created a structural variant map of unprecedented resolution and identified extensive genetic variation that has remained hidden until now. Many of these variants constitute candidates underlying phenotypic variation, including tandem duplications and a transposable element insertion that amplifies the expression of detoxification-related genes associated with nicotine resistance. The abundance of important genetic variation that still evades discovery highlights how crucial high-quality reference genomes are to deciphering phenotypes.

Causes of iron overload in blood donors - a clinical study

BACKGROUND AND OBJECTIVES

Despite the obligate iron loss from blood donation, some donors present with hyperferritinaemia that can result from a wide range of acute and chronic conditions including hereditary haemochromatosis (HH). The objective of our study was to investigate the causes of hyperferritinaemia in the blood donor population and explore the value of extensive HH mutational analyses.

MATERIALS AND METHODS

Forty-nine consecutive donors (f = 6, m = 43) were included prospectively from the Capital Regional Blood Center. Inclusion criteria were a single ferritin value >1000 μg/l or repeated hyperferritinaemia with at least one value >500 μg/l. All donors were questioned about their medical history and underwent a physical examination, biochemical investigations and next-generation sequencing of HH-related genes, including the HFE gene, the haemojuvelin gene (HFE2/HJV), the hepcidin gene (HAMP), the ferroportin 1 gene (SLC40A1) and the transferrin receptor 2 gene (TFR2).

RESULTS

Forty of 49 donors were mutation positive with a combined 69 mutations, 54 of which were located in the HFE gene. There were 11 mutations in the TFR2 gene, two mutations in the HFE2 gene and two mutations in the HAMP gene. Only four donors had apparent alternative causes of hyperferritinaemia.

CONCLUSION

HH-related mutations were the most frequent cause of hyperferritinaemia in a Danish blood donor population, and it appears that several different HH-genotypes can contribute to hyperferritinaemia. HH screening in blood donors with high ferritin levels could be warranted. HH-related iron overload should not in itself result in donor ineligibility.

Independent impacts of aging on mitochondrial DNA quantity and quality in humans

BACKGROUND

The accumulation of mitochondrial DNA (mtDNA) mutations, and the reduction of mtDNA copy number, both disrupt mitochondrial energetics, and may contribute to aging and age-associated phenotypes. However, there are few genetic and epidemiological studies on the spectra of blood mtDNA heteroplasmies, and the distribution of mtDNA copy numbers in different age groups and their impact on age-related phenotypes. In this work, we used whole-genome sequencing data of isolated peripheral blood mononuclear cells (PBMCs) from the UK10K project to investigate in parallel mtDNA heteroplasmy and copy number in 1511 women, between 17 and 85 years old, recruited in the TwinsUK cohorts.

RESULTS

We report a high prevalence of pathogenic mtDNA heteroplasmies in this population. We also find an increase in mtDNA heteroplasmies with age (β = 0.011, P = 5.77e-6), and showed that, on average, individuals aged 70-years or older had 58.5% more mtDNA heteroplasmies than those under 40-years old. Conversely, mtDNA copy number decreased by an average of 0.4 copies per year (β = -0.395, P = 0.0097). Multiple regression analyses also showed that age had independent effects on mtDNA copy number decrease and heteroplasmy accumulation. Finally, mtDNA copy number was positively associated with serum bicarbonate level (P = 4.46e-5), and inversely correlated with white blood cell count (P = 0.0006). Moreover, the aggregated heteroplasmy load was associated with blood apolipoprotein B level (P = 1.33e-5), linking the accumulation of mtDNA mutations to age-related physiological markers.

CONCLUSIONS

Our population-based study indicates that both mtDNA quality and quantity are influenced by age. An open question for the future is whether interventions that would contribute to maintain optimal mtDNA copy number and prevent the expansion of heteroplasmy could promote healthy aging.

SCA28: Novel Mutation in the AFG3L2 Proteolytic Domain Causes a Mild Cerebellar Syndrome with Selective Type-1 Muscle Fiber Atrophy

The spinocerebellar ataxias (SCA) are a group of rare inherited neurodegenerative diseases characterized by slowly progressive cerebellar ataxia, resulting in unsteady gait, clumsiness, and dysarthria. The disorders are predominantly inherited in an autosomal dominant manner. Mutations in the gene AFG3L2 that encodes a subunit of the mitochondrial m-AAA protease have previously been shown to cause spinocerebellar ataxia type 28 (SCA28). Here, we present the clinical phenotypes of three patients from a family with autosomal dominant cerebellar ataxia and show by molecular genetics and in silico modelling that this is caused by a novel missense mutation in the AFG3L2 gene. Furthermore, we show, for the first time, fluorodeoxyglucose-positron emission tomography (FDG-PET) scans of the brain and selective type I fiber atrophy of skeletal muscle of SCA28 patients indicating non-nervous-system involvement in SCA28 as well.

Prioritising Causal Genes at Type 2 Diabetes Risk Loci

PURPOSE OF REVIEW

Genome-wide association studies (GWAS) for type 2 diabetes (T2D) risk have identified a large number of genetic loci associated with disease susceptibility. However, progress moving from association signals through causal genes to functional understanding has so far been slow, hindering clinical translation. This review discusses the benefits and limitations of emerging, unbiased approaches for prioritising causal genes at T2D risk loci.

RECENT FINDINGS

Candidate causal genes can be identified by a number of different strategies that rely on genetic data, genomic annotations, and functional screening of selected genes. To overcome the limitations of each particular method, integration of multiple data sets is proving essential for establishing confidence in the prioritised genes. Previous studies have also highlighted the need to support these efforts through identification of causal variants and disease-relevant tissues. Prioritisation of causal genes at T2D risk loci by integrating complementary lines of evidence promises to accelerate our understanding of disease pathology and promote translation into new therapeutics.

Association Analysis of the MHC in Lupus Nephritis

Lupus nephritis (LN) is one of the most prevalent and serious complications of SLE, with significant effects on patient and renal survival. Although a large number of genetic variants associated with SLE have been identified, biomarkers that correlate with LN are extremely limited. In this study, we performed a comprehensive sequencing analysis of the whole MHC region in 1331 patients with LN and 1296 healthy controls and validated the independent associations in another 950 patients with LN and 1000 controls. We discovered five independent risk variants for LN within the MHC region, including HLA-DRβ1 amino acid 11 (P_omnibus<0.001), HLA-DQβ1 amino acid 45 (P<0.001; odds ratio, 0.58; 95% confidence interval, 0.52 to 0.65), HLA-A amino acid 156 (P_omnibus<0.001), HLA-DPβ1 amino acid 76 (P_omnibus<0.001), and a missense variant in PRRC2A (rs114580964; P<0.001; odds ratio, 0.38; 95% confidence interval, 0.30 to 0.49) at genome-wide significance. These data implicate aberrant peptide presentation by MHC classes 1 and 2 molecules and sex hormone modulation in the development of LN.

Mutation of the Planar Cell Polarity Gene VANGL1 in Adolescent Idiopathic Scoliosis

STUDY DESIGN

Mutation analysis of a candidate disease gene in a cohort of patients with moderate to severe Adolescent idiopathic scoliosis (AIS).

OBJECTIVE

To investigate if damaging mutations in the planar cell polarity gene VANGL1 could be identified in AIS patients.

SUMMARY OF BACKGROUND DATA

AIS is a spinal deformity which occurs in 1% to 3% of the population. The cause of AIS is often unknown, but genetic factors are important in the etiology. Rare variants in genes encoding regulators of WNT/planar cell polarity (PCP) signaling were recently identified in AIS patients.

METHODS

We analyzed the coding region of the VANGL1 gene for mutations using Sanger sequencing in 157 unrelated patients with moderate to severe AIS. The frequency of mutations in the patient cohort was compared with their frequency in a large cohort of controls. Functional effect of mutations were predicted in silico and analyzed in vitro by transfection of normal and mutant recombinant VANGL1 protein in Madin-Darby Canine Kidney (MDCK) cells. Cellular localization of recombinant proteins was analyzed by immunofluorescence microscopy analysis.

RESULTS

In the patient cohort, we identified two rare missense mutations in VANGL1, encoding a receptor involved in WNT/PCP signaling. The mutations, p.I136N and p.F440 V, are very rare in the normal population. Both mutations are predicted to be damaging, and to affect evolutionary conserved amino acid residues of VANGL1. Functional analysis in MDCK cells showed that the mutations abolished the normal translocation of VANGL1 to the cell membrane.

CONCLUSION

Our data support that mutations in genes involved in WNT/PCP signaling may be associated with AIS, but replication in other patient cohorts and further analysis of the role of WNT/PCP signaling in AIS is needed.

LEVEL OF EVIDENCE

4.

Clinical Characteristics, Mutation Spectrum, and Prevalence of Åland Eye Disease/Incomplete Congenital Stationary Night Blindness in Denmark

Purpose

To assess clinical characteristics, foveal structure, mutation spectrum, and prevalence rate of Åland eye disease (AED)/incomplete congenital stationary night blindness (iCSNB).

Methods

A retrospective survey included individuals diagnosed with AED at a national low-vision center from 1980 to 2014. A subset of affected males underwent ophthalmologic examinations including psychophysical tests, full-field electroretinography, and spectral-domain optical coherence tomography.

Results

Over the 34-year period, 74 individuals from 35 families were diagnosed with AED. Sixty individuals from 29 families participated in a follow-up study of whom 59 harbored a CACNA1F mutation and 1 harbored a CABP4 mutation. Among the subjects with a CACNA1F mutation, subnormal visual acuity was present in all, nystagmus was present in 63%, and foveal hypoplasia was observed in 25/43 subjects. Foveal pit volume was significantly reduced as compared to normal (P < 0.0001). Additionally, outer segment length at the fovea was measured in 46 subjects and found to be significantly reduced as compared to normal (P < 0.001). Twenty-nine CACNA1F variations were detected among 34 families in the total cohort, and a novel CABP4 variation was identified in one family. The estimated mean birth prevalence rate was 1 per 22,000 live-born males.

Conclusions

Our data support the viewpoint that AED, iCSNB, and X-linked cone–rod dystrophy 3 are designations that refer to a broad, continuous spectrum of clinical appearances caused in the majority by a variety of mutations in CACNA1F. We argue that the original designation AED should be used for this entity.

Comparison of single-marker and multi-marker tests in rare variant association studies of quantitative traits

In genetic association studies of rare variants, low statistical power and potential violations of established estimator properties are among the main challenges of association tests. Multi-marker tests (MMTs) have been proposed to target these challenges, but any comparison with single-marker tests (SMTs) has to consider that their aim is to identify causal genomic regions instead of variants. Valid power comparisons have been performed for the analysis of binary traits indicating that MMTs have higher power, but there is a lack of conclusive studies for quantitative traits. The aim of our study was therefore to fairly compare SMTs and MMTs in their empirical power to identify the same causal loci associated with a quantitative trait. The results of extensive simulation studies indicate that previous results for binary traits cannot be generalized. First, we show that for the analysis of quantitative traits, conventional estimation methods and test statistics of single-marker approaches have valid properties yielding association tests with valid type I error, even when investigating singletons or doubletons. Furthermore, SMTs lead to more powerful association tests for identifying causal genes than MMTs when the effect sizes of causal variants are large, and less powerful tests when causal variants have small effect sizes. For moderate effect sizes, whether SMTs or MMTs have higher power depends on the sample size and percentage of causal SNVs. For a more complete picture, we also compare the power in studies of quantitative and binary traits, and the power to identify causal genes with the power to identify causal rare variants. In a genetic association analysis of systolic blood pressure in the Genetic Analysis Workshop 19 data, SMTs yielded smaller p-values compared to MMTs for most of the investigated blood pressure genes, and were least influenced by the definition of gene regions.

Deep sequencing of atrial fibrillation patients with mitral valve regurgitation shows no evidence of mosaicism but reveals novel rare germline variants

BACKGROUND

Atrial fibrillation (AF) is the most common cardiac arrhythmia. Valvular heart disease is a strong predictor, yet the underlying molecular mechanisms are unknown.

OBJECTIVE

The purpose of this study was to investigate the prevalence of somatic variants in AF candidate genes in an AF patient population undergoing surgery for mitral valve regurgitation (MVR) to determine whether these patients are genetically predisposed to AF.

METHODS

DNA was extracted from blood and left atrial tissue from 44 AF patients with MVR. Using next-generation sequencing, we investigated 110 genes using the HaloPlex Target Enrichment System. MuTect software was used for identification of somatic point variants. We functionally characterized selected variants using electrophysiologic techniques.

RESULTS

No somatic variants were identified in the cardiac tissue. Thirty-three patients (75%) had a rare germline variation in ≥1 candidate genes. Fourteen variants were novel. Fifteen variants were predicted damaging or likely damaging in ≥6 in silico predictions. We identified rare variants in genes never directly associated with AF: KCNE4, SCN4B, NEURL1, and CAND2. Interestingly, 7 patients (16%) had variants in genes involved in cellular potassium handling. The variants KCNQ1 (p.G272S) and KCNH2 (p.A913V) resulted in gain of function due to faster activation (KCNQ1) and slowed deactivation kinetics (KCNQ1, KCNH2).

CONCLUSION

We did not find any somatic variants in patients with AF and MVR. Surprisingly, we found that our cohort of non-lone AF patients might, like lone AF patients, be predisposed to AF by rare germline variants. Our findings emphasize the extent of still unknown factors in the pathogenesis of AF.

High Prevalence of Diabetes-Predisposing Variants in MODY Genes Among Danish Women With Gestational Diabetes Mellitus

Context:

Gestational diabetes mellitus (GDM), defined as any degree of glucose intolerance with first recognition during pregnancy, is a heterogeneous form of diabetes characterized by various degrees of β-cell dysfunction.

Objectives:

We aimed to estimate the prevalence of possibly pathogenic variants in the maturity-onset diabetes of the young genes GCK, HNF1A, HNF4A, HNF1B, and INS among women with GDM. Furthermore, we examined the glucose tolerance status in variant carriers vs noncarriers at follow-up.

Design, Setting, and Patients:

We sequenced the coding regions and intron/exon boundaries of GCK, HNF1A, HNF4A, HNF1B, and INS using targeted region capture and next-generation sequencing in 354 Danish women with diet-treated GDM. Glucose tolerance was examined at follow-up 10 years after the index pregnancy.

Main Outcome Measures:

The prevalence of possibly pathogenic variants in GCK, HNF1A, HNF4A, HNF1B, and INS was estimated, and differences in anthropometric traits, high-sensitivity C-Reactive Protein (CRP), and glucose metabolism were measured.

Results:

At baseline, 17 possibly disease-causing variants were found in 21 women, revealing a combined GCK, HNF1A, HNF4A, HNF1B, and INS variant prevalence of 5.9% (95% confidence interval: 3.5% to 8.4%). At follow-up, 15 out of 135 women with diabetes (11%) were carriers of variants in GCK, HNF1A, HNF4A, HNF1B, or INS.

Conclusions:

Almost 6% of Danish women with diet-treated GDM have possibly pathogenic variants in GCK, HNF1A, HNF4A, HNF1B, or INS. These women are at high risk of developing diabetes after pregnancy. Thus screening for variants in GCK, HNF1A, HNF4A, HNF1B, and INS should be considered among women with GDM.

Gene and Variant Annotation for Mendelian Disorders in the Era of Advanced Sequencing Technologies

Comprehensive annotations of genetic and noncoding regions and corresponding accurate variant classification for Mendelian diseases are the next big challenge in the new genomic era of personalized medicine. Progress in the development of faster and more accurate pipelines for genome annotation and variant classification will lead to the discovery of more novel disease associations and candidate therapeutic targets. This ultimately will facilitate better patient recruitment in clinical trials. In this review, we describe the trends in research at the intersection of basic and clinical genomics that aims to increase understanding of overall genomic complexity, complex inheritance patterns of disease, and patient-phenotype-specific genomic associations. We describe the emerging field of translational functional genomics, which integrates other functional "-omics" approaches that support next-generation sequencing genomic data in order to facilitate personalized diagnostics, disease management, biomarker discovery, and medicine. We also discuss the utility of this integrated approach for diagnostic clinics and medical databases and its role in the future of personalized medicine.

The Weighting is the Hardest Part: On the Behavior of the Likelihood Ratio Test and the Score Test Under a Data-Driven Weighting Scheme in Sequenced Samples

Sequence-based association studies are at a critical inflexion point with the increasing availability of exome-sequencing data. A popular test of association is the sequence kernel association test (SKAT). Weights are embedded within SKAT to reflect the hypothesized contribution of the variants to the trait variance. Because the true weights are generally unknown, and so are subject to misspecification, we examined the efficiency of a data-driven weighting scheme. We propose the use of a set of theoretically defensible weighting schemes, of which, we assume, the one that gives the largest test statistic is likely to capture best the allele frequency-functional effect relationship. We show that the use of alternative weights obviates the need to impose arbitrary frequency thresholds. As both the score test and the likelihood ratio test (LRT) may be used in this context, and may differ in power, we characterize the behavior of both tests. The two tests have equal power, if the weights in the set included weights resembling the correct ones. However, if the weights are badly specified, the LRT shows superior power (due to its robustness to misspecification). With this data-driven weighting procedure the LRT detected significant signal in genes located in regions already confirmed as associated with schizophrenia - the PRRC2A (p = 1.020e-06) and the VARS2 (p = 2.383e-06) - in the Swedish schizophrenia case-control cohort of 11,040 individuals with exome-sequencing data. The score test is currently preferred for its computational efficiency and power. Indeed, assuming correct specification, in some circumstances, the score test is the most powerful test. However, LRT has the advantageous properties of being generally more robust and more powerful under weight misspecification. This is an important result given that, arguably, misspecified models are likely to be the rule rather than the exception in weighting-based approaches.

Phenotypes, genotypes, and prevalence of congenital myopathies older than 5 years in Denmark

Objective:

Congenital myopathy as a nosologic entity has long been recognized, but knowledge of overall and subtype prevalence and phenotype-genotype relationship is scarce, especially in the adult population.

Methods:

A national cohort of 107 patients ≥5 years diagnosed with congenital myopathy were prospectively assessed clinically, histologically, and genetically.

Results:

Twenty-five patients were excluded because of atypical features or alternative etiologies. The remaining 82 were on average 28 years old. Histologic examination revealed 14 (17%) with core disease, 15 (18%) centronuclear myopathy, 12 (15%) nemaline rods, 27 (33%) congenital fiber-type disproportion or type I predominance, and 14 (17%) nonspecific myopathic changes. Genetic etiology was identified in 46 patients (56.1%); 22.0% were heterozygous or compound heterozygous for mutations in RYR1, 7.3% had DNM2 mutations, and 7.3% NEB mutations. Less than 5% had mutations in ACTA1, TPM2/3, MTM1, TTN, SEPN1, or SC4NA. A genetic cause was established in 83% with specific histology (cores/rods/centronuclear myopathy) vs 29% with unspecific histology. The detailed clinical examination found gene-dependent discrepancies in the pattern of muscle affection and walking ability. Although walking ability was delayed in patients with ACTA1, TPM2/3, and RYR1 mutations, it was within normal limits in patients with NEB and DNM2 mutations.

Conclusions:

We found that overall, genetic and histologic prevalence of congenital myopathy in Denmark differs from previous retrospective reports. Less RYR1 and more DNM2 and NEB mutations and less core histology were present in our cohort. These differences may be explained by our prospective design, the older cohort of patients, and by differences in genetic background.

On the association analysis of genome-sequencing data: A spatial clustering approach for partitioning the entire genome into nonoverlapping windows

For the association analysis of whole-genome sequencing (WGS) studies, we propose an efficient and fast spatial-clustering algorithm. Compared to existing analysis approaches for WGS data, that define the tested regions either by sliding or consecutive windows of fixed sizes along variants, a meaningful grouping of nearby variants into consecutive regions has the advantage that, compared to sliding window approaches, the number of tested regions is likely to be smaller. In comparison to consecutive, fixed-window approaches, our approach is likely to group nearby variants together. Given existing biological evidence that disease-associated mutations tend to physically cluster in specific regions along the chromosome, the identification of meaningful groups of nearby located variants could thus lead to a potential power gain for association analysis. Our algorithm defines consecutive genomic regions based on the physical positions of the variants, assuming an inhomogeneous Poisson process and groups together nearby variants. As parameters are estimated locally, the algorithm takes the differing variant density along the chromosome into account and provides locally optimal partitioning of variants into consecutive regions. An R-implementation of the algorithm is provided. We discuss the theoretical advances of our algorithm compared to existing, window-based approaches and show the performance and advantage of our introduced algorithm in a simulation study and by an application to Alzheimer's disease WGS data. Our analysis identifies a region in the ITGB3 gene that potentially harbors disease susceptibility loci for Alzheimer's disease. The region-based association signal of ITGB3 replicates in an independent data set and achieves formally genome-wide significance. Software Implementation: An implementation of the algorithm in R is available at: https://github.com/heidefier/cluster_wgs_data.