Scripps Genome ADVISER: Annotation and Distributed Variant Interpretation SERver

Interrogating the Human Diplome: Computational Methods, Emerging Applications, and Challenges

Human DNA sequencing protocols have revolutionized human biology, biomedical science, and clinical practice, but still have very important limitations. One limitation is that most protocols do not separate or assemble (i.e., "phase") the nucleotide content of each of the maternally and paternally derived chromosomal homologs making up the 22 autosomal pairs and the chromosomal pair making up the pseudo-autosomal region of the sex chromosomes. This has led to a dearth of studies and a consequent underappreciation of many phenomena of fundamental importance to basic and clinical genomic science. We discuss a few protocols for obtaining phase information as well as their limitations, including those that could be used in tumor phasing settings. We then describe a number of biological and clinical phenomena that require phase information. These include phenomena that require precise knowledge of the nucleotide sequence in a chromosomal segment from germline or somatic cells, such as DNA binding events, and insight into unique cis vs. trans-acting functionally impactful variant combinations-for example, variants implicated in a phenotype governed by compound heterozygosity. In addition, we also comment on the need for reliable and consensus-based diploid-context computational workflows for variant identification as well as the need for laboratory-based functional verification strategies for validating cis vs. trans effects of variant combinations. We also briefly describe available resources, example studies, as well as areas of further research, and ultimately argue that the science behind the study of human diploidy, referred to as "diplomics," which will be enabled by nucleotide-level resolution of phased genomes, is a logical next step in the analysis of human genome biology.

Gain of function of a metalloproteinase associated with multiple myeloma, bicuspid aortic valve, and Von Hippel-Lindau syndrome

A patient diagnosed with multiple myeloma, bicuspid aortic valve, and Von Hippel-Lindau syndrome underwent whole-exome sequencing seeking a unified genetic cause for these three pathologies. The patient possessed a single-point mutation of arginine to cysteine (R24C) in the N-terminal region(pro-domain) of matrix metalloproteinase 9 (MMP-9). The pro-domain interacts with the catalytic site of this enzyme rendering it inactive. MMP-9 has previously been associated with all three pathologies suffered by the patient. We hypothesized that the observed mutation in the pro-domain would influence the activity of this enzyme. We expressed recombinant versions of MMP-9 and an investigation of their biochemical properties revealed that MMP-9 R24C is a constitutively active zymogen. To our knowledge, this is the first example of a mutation that discloses catalytic activity in the pro-form in any of the 24 human MMPs.

Genetic loci for alcohol-related life events and substance-induced affective symptoms: indexing the "dark side" of addiction

A limited number of genetic variants have been identified in traditional GWAS as risk or protective factors for alcohol use disorders (AUD) and related phenotypes. We herein report whole-genome association and rare-variant analyses on AUD traits in American Indians (AI) and European Americans (EA). We evaluated 742 AIs and 1711 EAs using low-coverage whole-genome sequencing. Phenotypes included: (1) a metric based on the occurrence of 36 alcohol-related life events that reflect AUD severity; (2) two alcohol-induced affective symptoms that accompany severe AUDs. We identified two new loci for alcohol-related life events with converging evidence from both cohorts: rare variants of K2P channel gene KCNK2, and rare missense and splice-site variants in pro-inflammatory mediator gene PDE4C. A NAF1-FSTL5 intergenic variant and an FSTL5 variant were respectively associated with alcohol-related life events in AI and EA. PRKG2 of serine/threonine protein kinase family, and rare variants in interleukin subunit gene EBI3 (IL-27B) were uniquely associated with alcohol-induced affective symptoms in AI. LncRNA LINC02347 on 12q24.32 was uniquely associated with alcohol-induced depression in EA. The top GWAS findings were primarily rare/low-frequency variants in AI, and common variants in EA. Adrenal gland was the most enriched in tissue-specific gene expression analysis for alcohol-related life events, and nucleus accumbens was the most enriched for alcohol-induced affective states in AI. Prefrontal cortex was the most enriched in EA for both traits. These studies suggest that whole-genome sequencing can identify novel, especially uncommon, variants associated with severe AUD phenotypes although the findings may be population specific.

Mutations of the mitochondrial carrier translocase channel subunit TIM22 cause early-onset mitochondrial myopathy

Protein import into mitochondria is facilitated by translocases within the outer and the inner mitochondrial membranes that are dedicated to a highly specific subset of client proteins. The mitochondrial carrier translocase (TIM22 complex) inserts multispanning proteins, such as mitochondrial metabolite carriers and translocase subunits (TIM23, TIM17A/B and TIM22), into the inner mitochondrial membrane. Both types of substrates are essential for mitochondrial metabolic function and biogenesis. Here, we report on a subject, diagnosed at 1.5 years, with a neuromuscular presentation, comprising hypotonia, gastroesophageal reflux disease and persistently elevated serum and Cerebrospinal fluid lactate (CSF). Patient fibroblasts displayed reduced oxidative capacity and altered mitochondrial morphology. Using trans-mitochondrial cybrid cell lines, we excluded a candidate variant in mitochondrial DNA as causative of these effects. Whole-exome sequencing identified compound heterozygous variants in the TIM22 gene (NM_013337), resulting in premature truncation in one allele (p.Tyr25Ter) and a point mutation in a conserved residue (p.Val33Leu), within the intermembrane space region, of the TIM22 protein in the second allele. Although mRNA transcripts of TIM22 were elevated, biochemical analyses revealed lower levels of TIM22 protein and an even greater deficiency of TIM22 complex formation. In agreement with a defect in carrier translocase function, carrier protein amounts in the inner membrane were found to be reduced. This is the first report of pathogenic variants in the TIM22 pore-forming subunit of the carrier translocase affecting the biogenesis of inner mitochondrial membrane proteins critical for metabolite exchange.

Dysregulation of NRXN1 by mutant MIR8485 leads to calcium overload in pre-synapses inducing neurodegeneration in Multiple sclerosis

OBJECTIVES

To identify Damaging mutations in microRNAs (miRNAs) and 3' untranslated regions (UTRs) of target genes to establish Multiple sclerosis (MS) disease pathway.

METHODS

Female aged 16, with Relapsing Remitting Multiple sclerosis (RRMS) was reported with initial symptoms of blurred vision, severe immobility, upper and lower limb numbness and backache. Whole Exome Sequencing (WES) and disease pathway analysis was performed to identify mutations in miRNAs and UTRs.

RESULTS

We identified Deleterious/Damaging multibase mutations in MIR8485 and NRXN1. miR-8485 was found carrying frameshift homozygous deletion of bases CA, while NRXN1 was found carrying nonframeshift homozygous substitution of bases CT to TC in exon 8 replacing Serine with Leucine.

CONCLUSIONS

Mutations in miR-8485 and NRXN1 was found to alter calcium homeostasis and NRXN1/NLGN1 cell adhesion molecule binding affinities. The miR-8485 mutation leads to overexpression of NRXN1 altering pre-synaptic Ca²⁺ homeostasis, inducing neurodegeneration.

Molecular Autopsy for Sudden Death in the Young: Is Data Aggregation the Key?

The Scripps molecular autopsy study seeks to incorporate genetic testing into the postmortem examination of cases of sudden death in the young (<45 years old). Here, we describe the results from the first 2 years of the study, which consisted of whole exome sequencing (WES) of a cohort of 50 cases predominantly from San Diego County. Apart from the individual description of cases, we analyzed the data at the cohort-level, which brought new perspectives on the genetic causes of sudden death. We investigated the advantages and disadvantages of using WES compared to a gene panel for cardiac disease (usually the first genetic test used by medical examiners). In an attempt to connect complex clinical phenotypes with genotypes, we classified samples by their genetic fingerprint. Finally, we studied the benefits of analyzing the mitochondrial DNA genome. In this regard, we found that half of the cases clinically diagnosed as sudden infant death syndrome had an increased ratio of heteroplasmic variants, and that the variants were also present in the mothers. We believe that community-based data aggregation and sharing will eventually lead to an improved classification of variants. Allele frequencies for the all cases can be accessed via our genomics browser at https://genomics.scripps.edu/browser.

Whole genome sequence association and ancestry-informed polygenic profile of EEG alpha in a Native American population

EEG alpha activity is the dominant oscillation in most adult humans, is highly heritable, and has been associated with a number of cognitive functions. Two EEG phenotypes, low- and high-voltage alpha (LVA & HVA), have been demonstrated to have high heritabilities. They have different prevalence depending on a population's ancestral origins. In the present study we assessed the influence of ancestry admixture on EEG alpha power, and conducted a whole genome sequencing association analysis and an ancestry-informed polygenic study on those phenotypes in a Native American (NA) population that has a high prevalence of LVA. Seven common variants, in LD with each other upstream from gene ASIC2, reached genome-wide significance (p = 2 × 10-8 ) having a positive association with alpha voltage. They had lower minor allele frequencies in the NAs than in a global population sample. Overall correlations between lower degrees of NA (higher degree European) ancestry and HVA, and higher degrees of NA and LVA were also found. Additionally a rare-variant gene-based study identified gene TIA1 being negatively associated with LVA. Approximately 3% of SNPs exhibited a 15-fold enrichment that explained nearly half of the total SNP-heritability for EEG alpha. These regions showed the most significant anti-correlations between NA ancestry and alpha voltage, and were enriched for genes and pathways mediating cognitive functions. Our findings suggested that these regions likely harbor causal variants for HVA, and lacking of such variants could explain the high prevalence of LVA in this NA population, possibly illuminating the ancestral origin and genetic basis for EEG alpha.

Mutation of WIF1: a potential novel cause of a Nail-Patella-like disorder

PURPOSE

Nail-Patella syndrome is a dominantly inherited genetic disorder characterized by abnormalities of the nails, knees, elbows, and pelvis. Nail abnormalities are the most constant feature of Nail-Patella syndrome. Pathogenic mutations in a single gene, LMX1B, a mesenchymal determinant of dorsal-ventral patterning, explain approximately 95% of Nail-Patella syndrome cases. However, 5% of cases remain unexplained.

METHODS

Here, we present exome sequencing and analysis of four generations of a family with a dominantly inherited Nail-Patella-like disorder (nail dysplasia with some features of Nail-Patella syndrome) who tested negative for LMX1B mutation.

RESULTS

We identify a loss-of-function mutation in WIF1 (NM_007191 p.W15*), which is involved in mesoderm segmentation, as the suspected cause of the Nail-Patella-like disorder observed in this family.

CONCLUSIONS

Mutation of WIF1 is a potential novel cause of a Nail-Patella-like disorder. Testing of additional patients negative for LMX1B mutation is needed to confirm this finding and further clarify the phenotype.Genet Med advance online publication 06 April 2017.

Influence of donor age on induced pluripotent stem cells

Induced pluripotent stem cells (iPSCs) are being pursued as a source of cells for autologous therapies, many of which will be aimed at aged patients. To explore the impact of age on iPSC quality, we produced iPSCs from blood cells of 16 donors aged 21-100. We find that iPSCs from older donors retain an epigenetic signature of age, which can be reduced through passaging. Clonal expansion via reprogramming also enables the discovery of somatic mutations present in individual donor cells, which are missed by bulk sequencing methods. We show that exomic mutations in iPSCs increase linearly with age, and all iPSC lines analyzed carry at least one gene-disrupting mutation, several of which have been associated with cancer or dysfunction. Unexpectedly, elderly donors (>90 yrs) harbor fewer mutations than predicted, likely due to a contracted blood progenitor pool. These studies establish that donor age is associated with an increased risk of abnormalities in iPSCs and will inform clinical development of reprogramming technology.

Whole-Genome Sequencing of a Healthy Aging Cohort

Studies of long-lived individuals have revealed few genetic mechanisms for protection against age-associated disease. Therefore, we pursued genome sequencing of a related phenotype-healthy aging-to understand the genetics of disease-free aging without medical intervention. In contrast with studies of exceptional longevity, usually focused on centenarians, healthy aging is not associated with known longevity variants, but is associated with reduced genetic susceptibility to Alzheimer and coronary artery disease. Additionally, healthy aging is not associated with a decreased rate of rare pathogenic variants, potentially indicating the presence of disease-resistance factors. In keeping with this possibility, we identify suggestive common and rare variant genetic associations implying that protection against cognitive decline is a genetic component of healthy aging. These findings, based on a relatively small cohort, require independent replication. Overall, our results suggest healthy aging is an overlapping but distinct phenotype from exceptional longevity that may be enriched with disease-protective genetic factors. VIDEO ABSTRACT.

Genetic variants associated with sleep disorders

OBJECTIVE

The diagnostic boundaries of sleep disorders are under considerable debate. The main sleep disorders are partly heritable; therefore, defining heritable pathophysiologic mechanisms could delineate diagnoses and suggest treatment. We collected clinical data and DNA from consenting patients scheduled to undergo clinical polysomnograms, to expand our understanding of the polymorphisms associated with the phenotypes of particular sleep disorders.

METHODS

Patients at least 21 years of age were recruited to contribute research questionnaires, and to provide access to their medical records, saliva for deoxyribonucleic acid (DNA), and polysomnographic data. From these complex data, 38 partly overlapping phenotypes were derived indicating complaints, subjective and objective sleep timing, and polysomnographic disturbances. A custom chip was used to genotype 768 single-nucleotide polymorphisms (SNPs). Additional assays derived ancestry-informative markers (eg, 751 participants of European ancestry). Linear regressions controlling for age, gender, and ancestry were used to assess the associations of each phenotype with each of the SNPs, highlighting those with Bonferroni-corrected significance.

RESULTS

In peroxisome proliferator-activated receptor gamma, coactivator 1 beta (PPARGC1B), rs6888451 was associated with several markers of obstructive sleep apnea. In aryl hydrocarbon receptor nuclear translocator-like (ARNTL), rs10766071 was associated with decreased polysomnographic sleep duration. The association of rs3923809 in BTBD9 with periodic limb movements in sleep was confirmed. SNPs in casein kinase 1 delta (CSNK1D rs11552085), cryptochrome 1 (CRY1 rs4964515), and retinoic acid receptor-related orphan receptor A (RORA rs11071547) were less persuasively associated with sleep latency and time of falling asleep.

CONCLUSIONS

SNPs associated with several sleep phenotypes were suggested, but due to risks of false discovery, independent replications are needed before the importance of these associations can be assessed, followed by investigation of molecular mechanisms.

Circadian polymorphisms in night owls, in bipolars, and in non-24-hour sleep cycles

People called night owls habitually have late bedtimes and late times of arising, sometimes suffering a heritable circadian disturbance called delayed sleep phase syndrome (DSPS). Those with DSPS, those with more severe progressively-late non-24-hour sleep-wake cycles, and those with bipolar disorder may share genetic tendencies for slowed or delayed circadian cycles. We searched for polymorphisms associated with DSPS in a case-control study of DSPS research participants and a separate study of Sleep Center patients undergoing polysomnography. In 45 participants, we resequenced portions of 15 circadian genes to identify unknown polymorphisms that might be associated with DSPS, non-24-hour rhythms, or bipolar comorbidities. We then genotyped single nucleotide polymorphisms (SNPs) in both larger samples, using Illumina Golden Gate assays. Associations of SNPs with the DSPS phenotype and with the morningness-eveningness parametric phenotype were computed for both samples, then combined for meta-analyses. Delayed sleep and "eveningness" were inversely associated with loci in circadian genes NFIL3 (rs2482705) and RORC (rs3828057). A group of haplotypes overlapping BHLHE40 was associated with non-24-hour sleep-wake cycles, and less robustly, with delayed sleep and bipolar disorder (e.g., rs34883305, rs34870629, rs74439275, and rs3750275 were associated with n=37, p=4.58E-09, Bonferroni p=2.95E-06). Bright light and melatonin can palliate circadian disorders, and genetics may clarify the underlying circadian photoperiodic mechanisms. After further replication and identification of the causal polymorphisms, these findings may point to future treatments for DSPS, non-24-hour rhythms, and possibly bipolar disorder or depression.