Rich annotation of DNA sequencing variants by leveraging the Ensembl Variant Effect Predictor with plugins

Grau M, Sánchez-Guixé M, López-Elorduy C, Reyes-Salazar I, Deu-Pons J, López-Bigas N, González-Pérez A. OpenVariant: a toolkit to parse and operate multiple input file formats

SUMMARY

Advances in high-throughput DNA sequencing technologies and decreasing costs have fueled the identification of small genetic variants (such as single nucleotide variants and indels) across tumors. Despite efforts to standardize variant formats and vocabularies, many sources of variability persist across databases and computational tools that annotate variants, hindering their integration within cancer genomic analyses. In this context, we present OpenVariant, an easily extendable Python package that facilitates seamless reading, parsing and refinement of diverse input file formats in a customizable structure, all within a single process.

AVAILABILITY AND IMPLEMENTATION

OpenVariant is an open-source package available at https://github.com/bbglab/openvariant. Documentation may be found at https://openvariant.readthedocs.io.

Strategies in Aggregation Tests for Rare Variants

Genome-wide association studies (GWAS) successfully identified numerous common variants involved in complex diseases, but only limited heritability was explained by these findings. Advances in high-throughput sequencing technology made it possible to assess the contribution of rare variants in common diseases. However, study of rare variants introduces challenges due to low frequency of rare variants. Well-established common variant methods were underpowered to identify the rare variants in GWAS. To address this challenge, several new methods have been developed to examine the role of rare variants in complex diseases. These approaches are based on testing the aggregate effect of multiple rare variants in a predefined genetic region. Provided here is an overview of statistical approaches and the protocols explaining step-by-step analysis of aggregations tests with the hands-on experience using R scripts in four categories: burden tests, adaptive burden tests, variance-component tests, and combined tests. Also explained are the concepts of rare variants, permutation tests, kernel methods, and genetic variant annotation. At the end we discuss relevant topics of bioinformatics tools for annotation, family-based design of rare-variant analysis, population stratification adjustment, and meta-analysis. © 2023 The Authors. Current Protocols published by Wiley Periodicals LLC.

Large-scale exome array summary statistics resources for glycemic traits to aid effector gene prioritization

Background

Genome-wide association studies for glycemic traits have identified hundreds of loci associated with these biomarkers of glucose homeostasis. Despite this success, the challenge remains to link variant associations to genes, and underlying biological pathways.

Methods

To identify coding variant associations which may pinpoint effector genes at both novel and previously established genome-wide association loci, we performed meta-analyses of exome-array studies for four glycemic traits: glycated hemoglobin (HbA1c, up to 144,060 participants), fasting glucose (FG, up to 129,665 participants), fasting insulin (FI, up to 104,140) and 2hr glucose post-oral glucose challenge (2hGlu, up to 57,878). In addition, we performed network and pathway analyses.

Results

Single-variant and gene-based association analyses identified coding variant associations at more than 60 genes, which when combined with other datasets may be useful to nominate effector genes. Network and pathway analyses identified pathways related to insulin secretion, zinc transport and fatty acid metabolism. HbA1c associations were strongly enriched in pathways related to blood cell biology.

Conclusions

Our results provided novel glycemic trait associations and highlighted pathways implicated in glycemic regulation. Exome-array summary statistic results are being made available to the scientific community to enable further discoveries.

A report of 5 Indian families with multicentric carpotarsal osteolysis syndrome

Multicentric carpotarsal osteolysis syndrome (MCTO) is a rare autosomal dominant skeletal dysplasia characterised by swelling and restriction of movement in the wrist and ankle joints, as well as osteolysis of the carpal and tarsal bones, that can be misdiagnosed as juvenile idiopathic arthritis. We describe five Indian families with heterozygous nonrecurrent missense pathogenic variants in exon 1 of MAF bZIP transcription factor B (MAFB).

Integration of Protein Structure and Population-Scale DNA Sequence Data for Disease Gene Discovery and Variant Interpretation

The experimental and computational techniques for capturing information about protein structures and genetic variation within the human genome have advanced dramatically in the past 20 years, generating extensive new data resources. In this review, we discuss these advances, along with new approaches for determining the impact a genetic variant has on protein function. We focus on the potential of new methods that integrate human genetic variation into protein structures to discover relationships to disease, including the discovery of mutational hotspots in cancer-related proteins, the localization of protein-altering variants within protein regions for common complex diseases, and the assessment of variants of unknown significance for Mendelian traits. We expect that approaches that integrate these data sources will play increasingly important roles in disease gene discovery and variant interpretation. Expected final online publication date for the Annual Review of Biomedical Data Science, Volume 5 is August 2022. Please see http://www.annualreviews.org/page/journal/pubdates for revised estimates.

Three chromosome-level duck genome assemblies provide insights into genomic variation during domestication

Domestic ducks are raised for meat, eggs and feather down, and almost all varieties are descended from the Mallard (Anas platyrhynchos). Here, we report chromosome-level high-quality genome assemblies for meat and laying duck breeds, and the Mallard. Our new genomic databases contain annotations for thousands of new protein-coding genes and recover a major percentage of the presumed "missing genes" in birds. We obtain the entire genomic sequences for the C-type lectin (CTL) family members that regulate eggshell biomineralization. Our population and comparative genomics analyses provide more than 36 million sequence variants between duck populations. Furthermore, a mutant cell line allows confirmation of the predicted anti-adipogenic function of NR2F2 in the duck, and uncovered mutations specific to Pekin duck that potentially affect adipose deposition. Our study provides insights into avian evolution and the genetics of oviparity, and will be a rich resource for the future genetic improvement of commercial traits in the duck.

Mitchell-Riley syndrome iPSCs exhibit reduced pancreatic endoderm differentiation due to a mutation in RFX6

Mitchell-Riley syndrome (MRS) is caused by recessive mutations in the regulatory factor X6 gene (RFX6) and is characterised by pancreatic hypoplasia and neonatal diabetes. To determine why individuals with MRS specifically lack pancreatic endocrine cells, we micro-CT imaged a 12-week-old foetus homozygous for the nonsense mutation RFX6 c.1129C>T, which revealed loss of the pancreas body and tail. From this foetus, we derived iPSCs and show that differentiation of these cells in vitro proceeds normally until generation of pancreatic endoderm, which is significantly reduced. We additionally generated an RFX6HA reporter allele by gene targeting in wild-type H9 cells to precisely define RFX6 expression and in parallel performed in situ hybridisation for RFX6 in the dorsal pancreatic bud of a Carnegie stage 14 human embryo. Both in vitro and in vivo, we find that RFX6 specifically labels a subset of PDX1-expressing pancreatic endoderm. In summary, RFX6 is essential for efficient differentiation of pancreatic endoderm, and its absence in individuals with MRS specifically impairs formation of endocrine cells of the pancreas head and tail.

A weighted empirical Bayes risk prediction model using multiple traits

With rapid advances in high-throughput sequencing technology, millions of single-nucleotide variants (SNVs) can be simultaneously genotyped in a sequencing study. These SNVs residing in functional genomic regions such as exons may play a crucial role in biological process of the body. In particular, non-synonymous SNVs are closely related to the protein sequence and its function, which are important in understanding the biological mechanism of sequence evolution. Although statistically challenging, models incorporating such SNV annotation information can improve the estimation of genetic effects, and multiple responses may further strengthen the signals of these variants on the assessment of disease risk. In this work, we develop a new weighted empirical Bayes method to integrate SNV annotation information in a multi-trait design. The performance of this proposed model is evaluated in simulation as well as a real sequencing data; thus, the proposed method shows improved prediction accuracy compared to other approaches.

Hadoop and PySpark for reproducibility and scalability of genomic sequencing studies

Modern genomic studies are rapidly growing in scale, and the analytical approaches used to analyze genomic data are increasing in complexity. Genomic data management poses logistic and computational challenges, and analyses are increasingly reliant on genomic annotation resources that create their own data management and versioning issues. As a result, genomic datasets are increasingly handled in ways that limit the rigor and reproducibility of many analyses. In this work, we examine the use of the Spark infrastructure for the management, access, and analysis of genomic data in comparison to traditional genomic workflows on typical cluster environments. We validate the framework by reproducing previously published results from the Alzheimer's Disease Sequencing Project. Using the framework and analyses designed using Jupyter notebooks, Spark provides improved workflows, reduces user-driven data partitioning, and enhances the portability and reproducibility of distributed analyses required for large-scale genomic studies.

PBX-WNT-P63-IRF6 pathway in nonsyndromic cleft lip and palate

Nonsyndromic cleft lip and palate (NSCLP) is one of the most common craniofacial anomalies in humans, affecting more than 135,000 newborns worldwide. NSCLP has a multifactorial etiology with more than 50 genes postulated to play an etiologic role. The genetic pathway comprised of Pbx-Wnt-p63-Irf6 genes was shown to control facial morphogenesis in mice and proposed as a regulatory pathway for NSCLP. Based on these findings, we investigated whether variation in PBX1, PBX2, and TP63, and their proposed interactions were associated with NSCLP. Fourteen single nucleotide variants (SNVs) in/nearby PBX1, PBX2, and TP63 were genotyped in 780 NSCLP families of nonHispanic white (NHW) and Hispanic ethnicities. Family-based association tests were performed for individual SNVs stratified by ethnicity and family history of NSCLP. Gene-gene interactions were also tested. A significant association was found for PBX2 rs3131300 and NSCLP in combined Hispanic families (p = .003) while nominal association was found for TP63 rs9332461 in multiplex Hispanic families (p = .005). Significant haplotype associations were observed for PBX2 in NHW (p = .0002) and Hispanic families (p = .003), and for TP63 in multiplex Hispanic families (.003). An independent case-control group was used to validate findings, and significant associations were found with PBX1 rs6426870 (p = .007) and TP63 rs9332461 (p = .03). Gene-gene interactions were detected between PBX1/PBX2/TP63 with IRF6 in NHW families, and between PBX1 with WNT9B in both NHW and Hispanic families (p < .0018). This study provides the first evidence for a role of PBX1 and PBX2, additional evidence for the role of TP63, and support for the proposed PBX-WNT-TP63-IRF6 regulatory pathway in the etiology of NSCLP.

Genome-wide discovery of somatic regulatory variants in diffuse large B-cell lymphoma

Diffuse large B-cell lymphoma (DLBCL) is an aggressive cancer originating from mature B-cells. Prognosis is strongly associated with molecular subgroup, although the driver mutations that distinguish the two main subgroups remain poorly defined. Through an integrative analysis of whole genomes, exomes, and transcriptomes, we have uncovered genes and non-coding loci that are commonly mutated in DLBCL. Our analysis has identified novel cis-regulatory sites, and implicates recurrent mutations in the 3′ UTR of NFKBIZ as a novel mechanism of oncogene deregulation and NF-κB pathway activation in the activated B-cell (ABC) subgroup. Small amplifications associated with over-expression of FCGR2B (the Fcγ receptor protein IIB), primarily in the germinal centre B-cell (GCB) subgroup, correlate with poor patient outcomes suggestive of a novel oncogene. These results expand the list of subgroup driver mutations that may facilitate implementation of improved diagnostic assays and could offer new avenues for the development of targeted therapeutics.

The driver mutations for the two main molecular subgroups of diffuse large B-cell lymphoma (DLBCL) are poorly defined. Here, an integrative genomics analysis identifies 3′ UTR NFKBIZ mutations within the activated B-cell DLBCL subgroup and small FCGR2B amplifications in the germinal centre B-cell DLBCL subgroup.

Functional annotation of genomic variants in studies of late-onset Alzheimer's disease

Motivation

Annotation of genomic variants is an increasingly important and complex part of the analysis of sequence-based genomic analyses. Computational predictions of variant function are routinely incorporated into gene-based analyses of rare-variants, though to date most studies use limited information for assessing variant function that is often agnostic of the disease being studied.

Results

In this work, we outline an annotation process motivated by the Alzheimer's Disease Sequencing Project, illustrate the impact of including tissue-specific transcript sets and sources of gene regulatory information and assess the potential impact of changing genomic builds on the annotation process. While these factors only impact a small proportion of total variant annotations (∼5%), they influence the potential analysis of a large fraction of genes (∼25%).

Availability and implementation

Individual variant annotations are available via the NIAGADS GenomicsDB, at https://www.niagads.org/genomics/ tools-and-software/databases/genomics-database. Annotations are also available for bulk download at https://www.niagads.org/datasets. Annotation processing software is available at http://www.icompbio.net/resources/software-and-downloads/.

Supplementary information

Supplementary data are available at Bioinformatics online.

Rare genetic variants in the endocannabinoid system genes CNR1 and DAGLA are associated with neurological phenotypes in humans

Rare genetic variants in the core endocannabinoid system genes CNR1, CNR2, DAGLA, MGLL and FAAH were identified in molecular testing data from 6,032 patients with a broad spectrum of neurological disorders. The variants were evaluated for association with phenotypes similar to those observed in the orthologous gene knockouts in mice. Heterozygous rare coding variants in CNR1, which encodes the type 1 cannabinoid receptor (CB1), were found to be significantly associated with pain sensitivity (especially migraine), sleep and memory disorders—alone or in combination with anxiety—compared to a set of controls without such CNR1 variants. Similarly, heterozygous rare variants in DAGLA, which encodes diacylglycerol lipase alpha, were found to be significantly associated with seizures and neurodevelopmental disorders, including autism and abnormalities of brain morphology, compared to controls. Rare variants in MGLL, FAAH and CNR2 were not associated with any neurological phenotypes in the patients tested. Diacylglycerol lipase alpha synthesizes the endocannabinoid 2-AG in the brain, which interacts with CB1 receptors. The phenotypes associated with rare CNR1 variants are reminiscent of those implicated in the theory of clinical endocannabinoid deficiency syndrome. The severe phenotypes associated with rare DAGLA variants underscore the critical role of rapid 2-AG synthesis and the endocannabinoid system in regulating neurological function and development. Mapping of the variants to the 3D structure of the type 1 cannabinoid receptor, or primary structure of diacylglycerol lipase alpha, reveals clustering of variants in certain structural regions and is consistent with impacts to function.

Knowledge Discovery in Biological Databases for Revealing Candidate Genes Linked to Complex Phenotypes

Genetics and “omics” studies designed to uncover genotype to phenotype relationships often identify large numbers of potential candidate genes, among which the causal genes are hidden. Scientists generally lack the time and technical expertise to review all relevant information available from the literature, from key model species and from a potentially wide range of related biological databases in a variety of data formats with variable quality and coverage. Computational tools are needed for the integration and evaluation of heterogeneous information in order to prioritise candidate genes and components of interaction networks that, if perturbed through potential interventions, have a positive impact on the biological outcome in the whole organism without producing negative side effects. Here we review several bioinformatics tools and databases that play an important role in biological knowledge discovery and candidate gene prioritization. We conclude with several key challenges that need to be addressed in order to facilitate biological knowledge discovery in the future.

Mutations in DYNC2H1, the cytoplasmic dynein 2, heavy chain 1 motor protein gene, cause short-rib polydactyly type I, Saldino-Noonan type

The short-rib polydactyly syndromes (SRPS) are autosomal recessively inherited, genetically heterogeneous skeletal ciliopathies. SRPS phenotypes were historically categorized as types I-IV, with type I first delineated by Saldino and Noonan in 1972. Characteristic findings among all forms of SRP include short horizontal ribs, short limbs and polydactyly. The SRP type I phenotype is characterized by a very small thorax, extreme micromelia, very short, poorly mineralized long bones, and multiple organ system anomalies. To date, the molecular basis of this most severe type of SRP, also known as Saldino-Noonan syndrome, has not been determined. We identified three SRP cases that fit the original phenotypic description of SRP type I. In all three cases, exome sequence analysis revealed compound heterozygosity for mutations in DYNC2H1, which encodes the main component of the retrograde IFT A motor, cytoplasmic dynein 2 heavy chain 1. Thus SRP type I, II, III and asphyxiating thoracic dystrophy (ATD), which also result from DYNC2H1 mutations. Herein we describe the phenotypic features, radiographic findings, and molecular basis of SRP type I.

Germline NLRP1 Mutations Cause Skin Inflammatory and Cancer Susceptibility Syndromes via Inflammasome Activation

Inflammasome complexes function as key innate immune effectors that trigger inflammation in response to pathogen- and danger-associated signals. Here, we report that germline mutations in the inflammasome sensor NLRP1 cause two overlapping skin disorders: multiple self-healing palmoplantar carcinoma (MSPC) and familial keratosis lichenoides chronica (FKLC). We find that NLRP1 is the most prominent inflammasome sensor in human skin, and all pathogenic NLRP1 mutations are gain-of-function alleles that predispose to inflammasome activation. Mechanistically, NLRP1 mutations lead to increased self-oligomerization by disrupting the PYD and LRR domains, which are essential in maintaining NLRP1 as an inactive monomer. Primary keratinocytes from patients experience spontaneous inflammasome activation and paracrine IL-1 signaling, which is sufficient to cause skin inflammation and epidermal hyperplasia. Our findings establish a group of non-fever inflammasome disorders, uncover an unexpected auto-inhibitory function for the pyrin domain, and provide the first genetic evidence linking NLRP1 to skin inflammatory syndromes and skin cancer predisposition.

Phenotypic spectrum of POLR3B mutations: isolated hypogonadotropic hypogonadism without neurological or dental anomalies

BACKGROUND

A constellation of neurodegenerative disorders exists (Gordon Holmes syndrome, 4H leucodystrophy, Boucher-Neuhauser syndrome) in which patients suffer from both neurological disease (typically manifested by ataxia) and reproductive failure (idiopathic hypogonadotropic hypogonadism (IHH)). POLR3B, which encodes the second largest subunit of RNA polymerase III (pol III), and POLR3A, which forms the pol III catalytic centre, are associated with 4H leucodystrophy.

METHODS

Whole exome sequencing was performed on a large cohort of subjects with IHH (n=565). Detailed neuroendocrine studies were performed in some individuals within this cohort.

RESULTS

Four individuals (two of them siblings) were identified with two rare nucleotide variants in POLR3B. On initial evaluation, all subjects were free of neurological disease. One patient underwent treatment with exogenous pulsatile gonadotropin-releasing hormone for 8 weeks which failed to result in normalisation of his sex steroid milieu due to pituitary resistance.

CONCLUSIONS

These findings suggest that the spectrum of phenotypes resulting from POLR3B mutations is wider than previously believed and that POLR3B can be associated exclusively with disorders characterised by abnormal gonadotropin secretion.

A Genome-Wide Association Study for Regulators of Micronucleus Formation in Mice

In mammals the regulation of genomic instability plays a key role in tumor suppression and also controls genome plasticity, which is important for recombination during the processes of immunity and meiosis. Most studies to identify regulators of genomic instability have been performed in cells in culture or in systems that report on gross rearrangements of the genome, yet subtle differences in the level of genomic instability can contribute to whole organism phenotypes such as tumor predisposition. Here we performed a genome-wide association study in a population of 1379 outbred Crl:CFW(SW)-US_P08 mice to dissect the genetic landscape of micronucleus formation, a biomarker of chromosomal breaks, whole chromosome loss, and extranuclear DNA. Variation in micronucleus levels is a complex trait with a genome-wide heritability of 53.1%. We identify seven loci influencing micronucleus formation (false discovery rate <5%), and define candidate genes at each locus. Intriguingly at several loci we find evidence for sexual dimorphism in micronucleus formation, with a locus on chromosome 11 being specific to males.

Truncating mutations in APP cause a distinct neurological phenotype

Dominant missense mutations in the amyloid β (Aβ) precursor protein (APP) gene have been implicated in early onset Alzheimer disease. These mutations alter protein structure to favor the pathologic production of Aβ. We report that homozygous nonsense mutations in APP are associated with decreased somatic growth, microcephaly, hypotonia, developmental delay, thinning of the corpus callosum, and seizures. We compare the phenotype of this case to those reported in mouse models and demonstrate multiple similarities, strengthening the role of amyloid precursor protein in normal brain function and development. Ann Neurol 2016;80:456-460.

Common Genetic Polymorphisms Influence Blood Biomarker Measurements in COPD

Implementing precision medicine for complex diseases such as chronic obstructive lung disease (COPD) will require extensive use of biomarkers and an in-depth understanding of how genetic, epigenetic, and environmental variations contribute to phenotypic diversity and disease progression. A meta-analysis from two large cohorts of current and former smokers with and without COPD [SPIROMICS (N = 750); COPDGene (N = 590)] was used to identify single nucleotide polymorphisms (SNPs) associated with measurement of 88 blood proteins (protein quantitative trait loci; pQTLs). PQTLs consistently replicated between the two cohorts. Features of pQTLs were compared to previously reported expression QTLs (eQTLs). Inference of causal relations of pQTL genotypes, biomarker measurements, and four clinical COPD phenotypes (airflow obstruction, emphysema, exacerbation history, and chronic bronchitis) were explored using conditional independence tests. We identified 527 highly significant (p < 8 X 10-10) pQTLs in 38 (43%) of blood proteins tested. Most pQTL SNPs were novel with low overlap to eQTL SNPs. The pQTL SNPs explained >10% of measured variation in 13 protein biomarkers, with a single SNP (rs7041; p = 10-392) explaining 71%-75% of the measured variation in vitamin D binding protein (gene = GC). Some of these pQTLs [e.g., pQTLs for VDBP, sRAGE (gene = AGER), surfactant protein D (gene = SFTPD), and TNFRSF10C] have been previously associated with COPD phenotypes. Most pQTLs were local (cis), but distant (trans) pQTL SNPs in the ABO blood group locus were the top pQTL SNPs for five proteins. The inclusion of pQTL SNPs improved the clinical predictive value for the established association of sRAGE and emphysema, and the explanation of variance (R2) for emphysema improved from 0.3 to 0.4 when the pQTL SNP was included in the model along with clinical covariates. Causal modeling provided insight into specific pQTL-disease relationships for airflow obstruction and emphysema. In conclusion, given the frequency of highly significant local pQTLs, the large amount of variance potentially explained by pQTL, and the differences observed between pQTLs and eQTLs SNPs, we recommend that protein biomarker-disease association studies take into account the potential effect of common local SNPs and that pQTLs be integrated along with eQTLs to uncover disease mechanisms. Large-scale blood biomarker studies would also benefit from close attention to the ABO blood group.

An inactivating mutation in intestinal cell kinase, ICK, impairs hedgehog signalling and causes short rib-polydactyly syndrome

The short rib polydactyly syndromes (SRPS) are a group of recessively inherited, perinatal-lethal skeletal disorders primarily characterized by short ribs, shortened long bones, varying types of polydactyly and concomitant visceral abnormalities. Mutations in several genes affecting cilia function cause SRPS, revealing a role for cilia function in skeletal development. To identify additional SRPS genes and discover novel ciliary molecules required for normal skeletogenesis, we performed exome sequencing in a cohort of patients and identified homozygosity for a missense mutation, p.E80K, in Intestinal Cell Kinase, ICK, in one SRPS family. The p.E80K mutation abolished serine/threonine kinase activity, resulting in altered ICK subcellular and ciliary localization, increased cilia length, aberrant cartilage growth plate structure, defective Hedgehog and altered ERK signalling. These data identify ICK as an SRPS-associated gene and reveal that abnormalities in signalling pathways contribute to defective skeletogenesis.

Reproducible pharmacogenomic profiling of cancer cell line panels

The use of large-scale genomic and drug response screening of cancer cell lines depends crucially on the reproducibility of results. Here we consider two previously published screens, plus a later critique of these studies. Using independent data, we show that consistency is achievable, and provide a systematic description of the best laboratory and analysis practices for future studies.

The Ensembl Variant Effect Predictor

The Ensembl Variant Effect Predictor is a powerful toolset for the analysis, annotation, and prioritization of genomic variants in coding and non-coding regions. It provides access to an extensive collection of genomic annotation, with a variety of interfaces to suit different requirements, and simple options for configuring and extending analysis. It is open source, free to use, and supports full reproducibility of results. The Ensembl Variant Effect Predictor can simplify and accelerate variant interpretation in a wide range of study designs.

The Ensembl Variant Effect Predictor

The Ensembl Variant Effect Predictor is a powerful toolset for the analysis, annotation, and prioritization of genomic variants in coding and non-coding regions. It provides access to an extensive collection of genomic annotation, with a variety of interfaces to suit different requirements, and simple options for configuring and extending analysis. It is open source, free to use, and supports full reproducibility of results. The Ensembl Variant Effect Predictor can simplify and accelerate variant interpretation in a wide range of study designs.

Exome Sequencing Identified a Splice Site Mutation in FHL1 that Causes Uruguay Syndrome, an X-Linked Disorder With Skeletal Muscle Hypertrophy and Premature Cardiac Death

BACKGROUND

Previously, we reported a rare X-linked disorder, Uruguay syndrome in a single family. The main features are pugilistic facies, skeletal deformities, and muscular hypertrophy despite a lack of exercise and cardiac ventricular hypertrophy leading to premature death.

METHODS AND RESULTS

An ≈19 Mb critical region on X chromosome was identified through identity-by-descent analysis of 3 affected males. Exome sequencing was conducted on one affected male to identify the disease-causing gene and variant. A splice site variant (c.502-2A>G) in the FHL1 gene was highly suspicious among other candidate genes and variants. FHL1A is the predominant isoform of FHL1 in cardiac and skeletal muscle. Sequencing cDNA showed the splice site variant led to skipping of exons 6 of the FHL1A isoform, equivalent to the FHL1C isoform. Targeted analysis showed that this splice site variant cosegregated with disease in the family. Western blot and immunohistochemical analysis of muscle from the proband showed a significant decrease in protein expression of FHL1A. Real-time polymerase chain reaction analysis of different isoforms of FHL1 demonstrated that the FHL1C is markedly increased.

CONCLUSIONS

Mutations in the FHL1 gene have been reported in disorders with skeletal and cardiac myopathy but none has the skeletal or facial phenotype seen in patients with Uruguay syndrome. Our data suggest that a novel FHL1 splice site variant results in the absence of FHL1A and the abundance of FHL1C, which may contribute to the complex and severe phenotype. Mutation screening of the FHL1 gene should be considered for patients with uncharacterized myopathies and cardiomyopathies.

A novel follicle-stimulating hormone receptor mutation causing primary ovarian failure: a fertility application of whole exome sequencing

STUDY QUESTION

Can whole exome sequencing (WES) and in vitro validation studies be used to find the causative genetic etiology in a patient with primary ovarian failure and infertility?

SUMMARY ANSWER

A novel follicle-stimulating hormone receptor (FSHR) mutation was found by WES and shown, via in vitro flow cytometry studies, to affect membrane trafficking.

WHAT IS KNOWN ALREADY

WES may diagnose up to 25-35% of patients with suspected disorders of sex development (DSD). FSHR mutations are an extremely rare cause of 46, XX gonadal dysgenesis with primary amenorrhea due to hypergonadotropic ovarian failure.

STUDY DESIGN, SIZE, DURATION

A WES study was followed by flow cytometry studies of mutant protein function.

PARTICIPANTS/MATERIALS, SETTING, METHODS

The study subjects were two Turkish sisters with hypergonadotropic primary amenorrhea, their parents and two unaffected sisters. The affected siblings and both parents were sequenced (trio-WES). Transient transfection of HEK 293T cells was performed with a vector containing wild-type FSHR as well as the novel FSHR variant that was discovered by WES. Cellular localization of FSHR protein as well as FSH-stimulated cyclic AMP (cAMP) production was evaluated using flow cytometry.

MAIN RESULTS AND THE ROLE OF CHANCE

Both affected sisters were homozygous for a previously unreported missense mutation (c.1222G>T, p.Asp408Tyr) in the second transmembrane domain of FSHR. Modeling predicted disrupted secondary structure. Flow cytometry demonstrated an average of 48% reduction in cell-surface signal detection (P < 0.01). The mean fluorescent signal for cAMP (second messenger of FSHR), stimulated by FSH, was reduced by 50% in the mutant-transfected cells (P < 0.01).

LIMITATIONS, REASONS FOR CAUTION

This is an in vitro validation. All novel purported genetic variants can be clinically reported only as 'variants of uncertain significance' until more patients with a similar phenotype are discovered with the same variant.

WIDER IMPLICATIONS OF THE FINDINGS

We report the first WES-discovered FSHR mutation, validated by quantitative flow cytometry. WES is a valuable tool for diagnosis of rare genetic diseases, and flow cytometry allows for quantitative characterization of purported variants. WES-assisted diagnosis allows for treatments aimed at the underlying molecular etiology of disease. Future studies should focus on pharmacological and assisted reproductive treatments aimed at the disrupted FSHR, so that patients with FSH resistance can be treated by personalized medicine.

STUDY FUNDING/COMPETING INTERESTS

E.V. is partially funded by the DSD Translational Research Network (NICHD 1R01HD068138). M.S.B. is funded by the Neuroendocrinology, Sex Differences and Reproduction training grant (NICHD 5T32HD007228). The authors have no competing interests to disclose.

Next-generation sequencing using a pre-designed gene panel for the molecular diagnosis of congenital disorders in pediatric patients

BACKGROUND

Next-generation sequencing (NGS) has revolutionized genetic research and offers enormous potential for clinical application. Sequencing the exome has the advantage of casting the net wide for all known coding regions while targeted gene panel sequencing provides enhanced sequencing depths and can be designed to avoid incidental findings in adult-onset conditions. A HaloPlex panel consisting of 180 genes within commonly altered chromosomal regions is available for use on both the Ion Personal Genome Machine (PGM) and MiSeq platforms to screen for causative mutations in these genes.

METHODS

We used this Haloplex ICCG panel for targeted sequencing of 15 patients with clinical presentations indicative of an abnormality in one of the 180 genes. Sequencing runs were done using the Ion 318 Chips on the Ion Torrent PGM. Variants were filtered for known polymorphisms and analysis was done to identify possible disease-causing variants before validation by Sanger sequencing. When possible, segregation of variants with phenotype in family members was performed to ascertain the pathogenicity of the variant.

RESULTS

More than 97% of the target bases were covered at >20×. There was an average of 9.6 novel variants per patient. Pathogenic mutations were identified in five genes for six patients, with two novel variants. There were another five likely pathogenic variants, some of which were unreported novel variants.

CONCLUSIONS

In a cohort of 15 patients, we were able to identify a likely genetic etiology in six patients (40%). Another five patients had candidate variants for which further evaluation and segregation analysis are ongoing. Our results indicate that the HaloPlex ICCG panel is useful as a rapid, high-throughput and cost-effective screening tool for 170 of the 180 genes. There is low coverage for some regions in several genes which might have to be supplemented by Sanger sequencing. However, comparing the cost, ease of analysis, and shorter turnaround time, it is a good alternative to exome sequencing for patients whose features are suggestive of a genetic etiology involving one of the genes in the panel.

Mutadelic: mutation analysis using description logic inferencing capabilities

MOTIVATION

As next generation sequencing gains a foothold in clinical genetics, there is a need for annotation tools to characterize increasing amounts of patient variant data for identifying clinically relevant mutations. While existing informatics tools provide efficient bulk variant annotations, they often generate excess information that may limit their scalability.

RESULTS

We propose an alternative solution based on description logic inferencing to generate workflows that produce only those annotations that will contribute to the interpretation of each variant. Workflows are dynamically generated using a novel abductive reasoning framework called a basic framework for abductive workflow generation (AbFab). Criteria for identifying disease-causing variants in Mendelian blood disorders were identified and implemented as AbFab services. A web application was built allowing users to run workflows generated from the criteria to analyze genomic variants. Significant variants are flagged and explanations provided for why they match or fail to match the criteria.

AVAILABILITY AND IMPLEMENTATION

The Mutadelic web application is available for use at http://krauthammerlab.med.yale.edu/mutadelic.

CONTACT

michael.krauthammer@yale.edu.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Mutations in DYNC2LI1 disrupt cilia function and cause short rib polydactyly syndrome

The short rib polydactyly syndromes (SRPSs) are a heterogeneous group of autosomal recessive, perinatal lethal skeletal disorders characterized primarily by short, horizontal ribs, short limbs and polydactyly. Mutations in several genes affecting intraflagellar transport (IFT) cause SRPS but they do not account for all cases. Here we identify an additional SRPS gene and further unravel the functional basis for IFT. We perform whole-exome sequencing and identify mutations in a new disease-producing gene, cytoplasmic dynein-2 light intermediate chain 1, DYNC2LI1, segregating with disease in three families. Using primary fibroblasts, we show that DYNC2LI1 is essential for dynein-2 complex stability and that mutations in DYNC2LI1 result in variable length, including hyperelongated, cilia, Hedgehog pathway impairment and ciliary IFT accumulations. The findings in this study expand our understanding of SRPS locus heterogeneity and demonstrate the importance of DYNC2LI1 in dynein-2 complex stability, cilium function, Hedgehog regulation and skeletogenesis.

Inherited CHST11/MIR3922 deletion is associated with a novel recessive syndrome presenting with skeletal malformation and malignant lymphoproliferative disease

Glycosaminoglycans (GAGs) such as chondroitin are ubiquitous disaccharide carbohydrate chains that contribute to the formation and function of proteoglycans at the cell membrane and in the extracellular matrix. Although GAG-modifying enzymes are required for diverse cellular functions, the role of these proteins in human development and disease is less well understood. Here, we describe two sisters out of seven siblings affected by congenital limb malformation and malignant lymphoproliferative disease. Using Whole-Genome Sequencing (WGS), we identified in the proband deletion of a 55 kb region within chromosome 12q23 that encompasses part of CHST11 (encoding chondroitin-4-sulfotransferase 1) and an embedded microRNA (MIR3922). The deletion was homozygous in the proband but not in each of three unaffected siblings. Genotyping data from the 1000 Genomes Project suggest that deletions inclusive of both CHST11 and MIR3922 are rare events. Given that CHST11 deficiency causes severe chondrodysplasia in mice that is similar to human limb malformation, these results underscore the importance of chondroitin modification in normal skeletal development. Our findings also potentially reveal an unexpected role for CHST11 and/or MIR3922 as tumor suppressors whose disruption may contribute to malignant lymphoproliferative disease.

SNPchiMp v.3: integrating and standardizing single nucleotide polymorphism data for livestock species

BACKGROUND

In recent years, the use of genomic information in livestock species for genetic improvement, association studies and many other fields has become routine. In order to accommodate different market requirements in terms of genotyping cost, manufacturers of single nucleotide polymorphism (SNP) arrays, private companies and international consortia have developed a large number of arrays with different content and different SNP density. The number of currently available SNP arrays differs among species: ranging from one for goats to more than ten for cattle, and the number of arrays available is increasing rapidly. However, there is limited or no effort to standardize and integrate array- specific (e.g. SNP IDs, allele coding) and species-specific (i.e. past and current assemblies) SNP information.

RESULTS

Here we present SNPchiMp v.3, a solution to these issues for the six major livestock species (cow, pig, horse, sheep, goat and chicken). Original data was collected directly from SNP array producers and specific international genome consortia, and stored in a MySQL database. The database was then linked to an open-access web tool and to public databases. SNPchiMp v.3 ensures fast access to the database (retrieving within/across SNP array data) and the possibility of annotating SNP array data in a user-friendly fashion.

CONCLUSIONS

This platform allows easy integration and standardization, and it is aimed at both industry and research. It also enables users to easily link the information available from the array producer with data in public databases, without the need of additional bioinformatics tools or pipelines. In recognition of the open-access use of Ensembl resources, SNPchiMp v.3 was officially credited as an Ensembl E!mpowered tool. Availability at http://bioinformatics.tecnoparco.org/SNPchimp.

Loss of ADAM17 is associated with severe multiorgan dysfunction

ADAM metallopeptidase domain 17 (ADAM17) is responsible for processing large numbers of proteins. Recently, 1 family involving 2 patients with a homozygous mutation in ADAM17 were described, presenting with skin lesions and diarrhea. In this report, we describe a second family confirming the existence of this syndrome. The proband presented with severe diarrhea, skin rash, and recurrent sepsis, eventually leading to her death at the age of 10 months. We performed exome sequencing and detailed pathological and immunological investigations. We identified a novel homozygous frameshift mutation in ADAM17 (NM_003183.4:c.308dupA) leading to a premature stop codon. CD4(+) and CD8(+) T-cell stimulation assays showed severely diminished tumor necrosis factor-α and interleukin-2 production. Skin biopsies indicated a focal neutrophilic infiltrate and spongiotic dermatitis. Interestingly, the patient developed unexplained systolic hypertension and nonspecific hepatitis with apoptosis. This report provides evidence for an important role of ADAM17 in human immunological response and underscores its multiorgan involvement.

Exome sequencing for the diagnosis of 46,XY disorders of sex development

CONTEXT

Disorders of sex development (DSD) are clinical conditions where there is a discrepancy between the chromosomal sex and the phenotypic (gonadal or genital) sex of an individual. Such conditions can be stressful for patients and their families and have historically been difficult to diagnose, especially at the genetic level. In particular, for cases of 46,XY gonadal dysgenesis, once variants in SRY and NR5A1 have been ruled out, there are few other single gene tests available.

OBJECTIVE

We used exome sequencing followed by analysis with a list of all known human DSD-associated genes to investigate the underlying genetic etiology of 46,XY DSD patients who had not previously received a genetic diagnosis.

DESIGN

Samples were either submitted to the research laboratory or submitted as clinical samples to the UCLA Clinical Genomic Center. Sequencing data were filtered using a list of genes known to be involved in DSD.

RESULTS

We were able to identify a likely genetic diagnosis in more than a third of cases, including 22.5% with a pathogenic finding, an additional 12.5% with likely pathogenic findings, and 15% with variants of unknown clinical significance.

CONCLUSIONS

Early identification of the genetic cause of a DSD will in many cases streamline and direct the clinical management of the patient, with more focused endocrine and imaging studies and better-informed surgical decisions. Exome sequencing proved an efficient method toward such a goal in 46,XY DSD patients.

Clinical exome sequencing for genetic identification of rare Mendelian disorders

IMPORTANCE

Clinical exome sequencing (CES) is rapidly becoming a common molecular diagnostic test for individuals with rare genetic disorders.

OBJECTIVE

To report on initial clinical indications for CES referrals and molecular diagnostic rates for different indications and for different test types.

DESIGN, SETTING, AND PARTICIPANTS

Clinical exome sequencing was performed on 814 consecutive patients with undiagnosed, suspected genetic conditions at the University of California, Los Angeles, Clinical Genomics Center between January 2012 and August 2014. Clinical exome sequencing was conducted as trio-CES (both parents and their affected child sequenced simultaneously) to effectively detect de novo and compound heterozygous variants or as proband-CES (only the affected individual sequenced) when parental samples were not available.

MAIN OUTCOMES AND MEASURES

Clinical indications for CES requests, molecular diagnostic rates of CES overall and for phenotypic subgroups, and differences in molecular diagnostic rates between trio-CES and proband-CES.

RESULTS

Of the 814 cases, the overall molecular diagnosis rate was 26% (213 of 814; 95% CI, 23%-29%). The molecular diagnosis rate for trio-CES was 31% (127 of 410 cases; 95% CI, 27%-36%) and 22% (74 of 338 cases; 95% CI, 18%-27%) for proband-CES. In cases of developmental delay in children (<5 years, n = 138), the molecular diagnosis rate was 41% (45 of 109; 95% CI, 32%-51%) for trio-CES cases and 9% (2 of 23, 95% CI, 1%-28%) for proband-CES cases. The significantly higher diagnostic yield (P value = .002; odds ratio, 7.4 [95% CI, 1.6-33.1]) of trio-CES was due to the identification of de novo and compound heterozygous variants.

CONCLUSIONS AND RELEVANCE

In this sample of patients with undiagnosed, suspected genetic conditions, trio-CES was associated with higher molecular diagnostic yield than proband-CES or traditional molecular diagnostic methods. Additional studies designed to validate these findings and to explore the effect of this approach on clinical and economic outcomes are warranted.

Expanding the mutational spectrum of LZTR1 in schwannomatosis

Schwannomatosis is characterized by the development of multiple non-vestibular, non-intradermal schwannomas. Constitutional inactivating variants in two genes, SMARCB1 and, very recently, LZTR1, have been reported. We performed exome sequencing of 13 schwannomatosis patients from 11 families without SMARCB1 deleterious variants. We identified four individuals with heterozygous loss-of-function variants in LZTR1. Sequencing of the germline of 60 additional patients identified 18 additional heterozygous variants in LZTR1. We identified LZTR1 variants in 43% and 30% of familial (three of the seven families) and sporadic patients, respectively. In addition, we tested LZTR1 protein immunostaining in 22 tumors from nine unrelated patients with and without LZTR1 deleterious variants. Tumors from individuals with LZTR1 variants lost the protein expression in at least a subset of tumor cells, consistent with a tumor suppressor mechanism. In conclusion, our study demonstrates that molecular analysis of LZTR1 may contribute to the molecular characterization of schwannomatosis patients, in addition to NF2 mutational analysis and the detection of chromosome 22 losses in tumor tissue. It will be especially useful in differentiating schwannomatosis from mosaic Neurofibromatosis type 2 (NF2). However, the role of LZTR1 in the pathogenesis of schwannomatosis needs further elucidation.

De Novo variants in the KMT2A (MLL) gene causing atypical Wiedemann-Steiner syndrome in two unrelated individuals identified by clinical exome sequencing

Background

Wiedemann-Steiner Syndrome (WSS) is characterized by short stature, a variety of dysmorphic facial and skeletal features, characteristic hypertrichosis cubiti (excessive hair on the elbows), mild-to-moderate developmental delay and intellectual disability. [MIM#: 605130]. Here we report two unrelated children for whom clinical exome sequencing of parent-proband trios was performed at UCLA, resulting in a molecular diagnosis of WSS and atypical clinical presentation.

Case presentation

For patient 1, clinical features at 9 years of age included developmental delay, craniofacial abnormalities, and multiple minor anomalies. Patient 2 presented at 1 year of age with developmental delay, microphthalmia, partial 3–4 left hand syndactyly, and craniofacial abnormalities. A de novo missense c.4342T>C variant and a de novo splice site c.4086+G>A variant were identified in the KMT2A gene in patients 1 and 2, respectively.

Conclusions

Based on the clinical and molecular findings, both patients appear to have novel presentations of WSS. As the hallmark hypertrichosis cubiti was not initially appreciated in either case, this syndrome was not suspected during the clinical evaluation. This report expands the phenotypic spectrum of the clinical phenotypes and KMT2A variants associated with WSS.