Search

Large-scale studies spanning diverse project sites, populations, languages, and measurements are increasingly important to relate psychological to biological variables. National and international consortia already are collecting and executing mega-analyses on aggregated data from individuals, with different measures on each person. In this research, we show that Asparouhov and Muthén's alignment method can be adapted to align data from disparate item sets and response formats. We argue that with these adaptations, the alignment method is well suited for combining data across multiple sites even when they use different measurement instruments. The approach is illustrated using data from the Whole Genome Sequencing in Psychiatric Disorders consortium and a real-data-based simulation is used to verify accurate parameter recovery. Factor alignment appears to increase precision of measurement and validity of scores with respect to external criteria. The resulting parameter estimates may further inform development of more effective and efficient methods to assess the same constructs in prospectively designed studies. [The GROUP Investigators and the WGSPD Consortium contributed to the writing of this article.]

The flipped classroom is a relatively new active learning pedagogical intervention, gaining popularity as a blended learning methodology. The flipped classroom comprises two distinct parts, directed learning carried out at the student's own pace away from the classroom and an interactive, class-based activity encouraging problem-solving and experiential learning. This research presents a 1-year study to measure student performance and perception toward a flipped classroom approach to teaching core biochemical calculations to first-year undergraduate biochemistry and genetics students. A post-task questionnaire showed an overall positive student perception with an associated significant improvement in the end of module summative assessment. These results suggest that this teaching approach offers some advantages over more traditional teaching pedagogies. [This article also appears in the virtual issue of "Biochemistry and Molecular Biology Education," "Teaching in the Time of COVID-19."]

Polymerization of nucleotides under prebiotic conditions simulating the early Earth has been extensively studied. Several independent methods have been used to verify that RNA-like polymers can be produced by hot wet-dry cycling of nucleotides. However, it has not been shown that these RNA-like polymers are similar to biological RNA with 3'-5' phosphodiester bonds. In the results described here, RNA-like polymers were generated from 5'-monophosphate nucleosides AMP and UMP. To confirm that the polymers resemble biological RNA, ribonuclease A should catalyze hydrolysis of the 3'-5' phosphodiester bonds between pyrimidine nucleotides to each other or to purine nucleotides, but not purine-purine nucleotide bonds. Here we show AFM images of specific polymers produced by hot wet-dry cycling of AMP, UMP and AMP/UMP (1:1) solutions on mica surfaces, before and after exposure to ribonuclease A. AMP polymers were unaffected by ribonuclease A but UMP polymers disappeared. This indicates that a major fraction of the bonds in the UMP polymers is indeed 3'-5' phosphodiester bonds. Some of the polymers generated from the AMP/UMP mixture also showed clear signs of cleavage. Because ribonuclease A recognizes the ester bonds in the polymers, we show for the first time that these prebiotically produced polymers are in fact similar to biological RNA but are likely to be linked by a mixture of 3'-5' and 2'-5' phosphodiester bonds.

Leukemia-initiating cells (LICs) are regarded as the origin of leukemia relapse and therapeutic resistance. Identifying direct stemness determinants that fuel LIC self-renewal is critical for developing targeted approaches. Here, we show that the RNA-editing enzyme ADAR1 is a crucial stemness factor that promotes LIC self-renewal by attenuating aberrant double-stranded RNA (dsRNA) sensing. Elevated adenosine-to-inosine editing is a common attribute of relapsed T cell acute lymphoblastic leukemia (T-ALL) regardless of molecular subtype. Consequently, knockdown of ADAR1 severely inhibits LIC self-renewal capacity and prolongs survival in T-ALL patient-derived xenograft models. Mechanistically, ADAR1 directs hyper-editing of immunogenic dsRNA to avoid detection by the innate immune sensor melanoma differentiation-associated protein 5 (MDA5). Moreover, we uncover that the cell-intrinsic level of MDA5 dictates the dependency on the ADAR1-MDA5 axis in T-ALL. Collectively, our results show that ADAR1 functions as a self-renewal factor that limits the sensing of endogenous dsRNA. Thus, targeting ADAR1 presents an effective therapeutic strategy for eliminating T-ALL LICs.

Abstract Background Pubertal growth patterns correlate with future health outcomes. However, the genetic mechanisms mediating growth trajectories remain largely unknown. Here, we modeled longitudinal height growth with Super-Imposition by Translation And Rotation (SITAR) growth curve analysis on ~ 56,000 trans-ancestry samples with repeated height measurements from age 5 years to adulthood. We performed genetic analysis on six phenotypes representing the magnitude, timing, and intensity of the pubertal growth spurt. To investigate the lifelong impact of genetic variants associated with pubertal growth trajectories, we performed genetic correlation analyses and phenome-wide association studies in the Penn Medicine BioBank and the UK Biobank. Results Large-scale growth modeling enables an unprecedented view of adolescent growth across contemporary and 20th-century pediatric cohorts. We identify 26 genome-wide significant loci and leverage trans-ancestry data to perform fine-mapping. Our data reveals genetic relationships between pediatric height growth and health across the life course, with different growth trajectories correlated with different outcomes. For instance, a faster tempo of pubertal growth correlates with higher bone mineral density, HOMA-IR, fasting insulin, type 2 diabetes, and lung cancer, whereas being taller at early puberty, taller across puberty, and having quicker pubertal growth were associated with higher risk for atrial fibrillation. Conclusion We report novel genetic associations with the tempo of pubertal growth and find that genetic determinants of growth are correlated with reproductive, glycemic, respiratory, and cardiac traits in adulthood. These results aid in identifying specific growth trajectories impacting lifelong health and show that there may not be a single “optimal” pubertal growth pattern.

Abstract Background The presence of coronary plaques with high-risk characteristics is strongly associated with adverse cardiac events beyond the identification of coronary stenosis. Testing by coronary computed tomography angiography (CCTA) enables the identification of high-risk plaques (HRP). Referral for CCTA is presently based on pre-test probability estimates including clinical risk factors (CRFs); however, proteomics and/or genetic information could potentially improve patient selection for CCTA and, hence, identification of HRP. We aimed to (1) identify proteomic and genetic features associated with HRP presence and (2) investigate the effect of combining CRFs, proteomics, and genetics to predict HRP presence. Methods Consecutive chest pain patients (n = 1462) undergoing CCTA to diagnose obstructive coronary artery disease (CAD) were included. Coronary plaques were assessed using a semi-automatic plaque analysis tool. Measurements of 368 circulating proteins were obtained with targeted Olink panels, and DNA genotyping was performed in all patients. Imputed genetic variants were used to compute a multi-trait multi-ancestry genome-wide polygenic score (GPSMult). HRP presence was defined as plaques with two or more high-risk characteristics (low attenuation, spotty calcification, positive remodeling, and napkin ring sign). Prediction of HRP presence was performed using the glmnet algorithm with repeated fivefold cross-validation, using CRFs, proteomics, and GPSMult as input features. Results HRPs were detected in 165 (11%) patients, and 15 input features were associated with HRP presence. Prediction of HRP presence based on CRFs yielded a mean area under the receiver operating curve (AUC) ± standard error of 73.2 ± 0.1, versus 69.0 ± 0.1 for proteomics and 60.1 ± 0.1 for GPSMult. Combining CRFs with GPSMult increased prediction accuracy (AUC 74.8 ± 0.1 (P = 0.004)), while the inclusion of proteomics provided no significant improvement to either the CRF (AUC 73.2 ± 0.1, P = 1.00) or the CRF + GPSMult (AUC 74.6 ± 0.1, P = 1.00) models, respectively. Conclusions In patients with suspected CAD, incorporating genetic data with either clinical or proteomic data improves the prediction of high-risk plaque presence. Trial registration https://clinicaltrials.gov/ct2/show/NCT02264717 (September 2014).

Abstract Background Despite a clear appreciation of the impact of human pathogens on community health, efforts to understand pathogen dynamics within populations often follow a narrow-targeted approach and rely on the deployment of specific molecular probes for quantitative detection or rely on clinical detection and reporting. Main text Genomic analysis of wastewater samples for the broad detection of viruses, bacteria, fungi, and antibiotic resistance genes of interest/concern is inherently difficult, and while deep sequencing of wastewater provides a wealth of information, a robust and cooperative foundation is needed to support healthier communities. In addition to furthering the capacity of high-throughput sequencing wastewater-based epidemiology to detect human pathogens in an unbiased and agnostic manner, it is critical that collaborative networks among public health agencies, researchers, and community stakeholders be fostered to prepare communities for future public health emergencies or for the next pandemic. A more inclusive public health infrastructure must be built for better data reporting where there is a global human health risk burden. Conclusions As wastewater platforms continue to be developed and refined, high-throughput sequencing of human pathogens in wastewater samples will emerge as a gold standard for understanding community health.

A major challenge of genome-wide association studies (GWASs) is to translate phenotypic associations into biological insights. Here, we integrate a large GWAS on blood lipids involving 1.6 million individuals from five ancestries with a wide array of functional genomic datasets to discover regulatory mechanisms underlying lipid associations. We first prioritize lipid-associated genes with expression quantitative trait locus (eQTL) colocalizations and then add chromatin interaction data to narrow the search for functional genes. Polygenic enrichment analysis across 697 annotations from a host of tissues and cell types confirms the central role of the liver in lipid levels and highlights the selective enrichment of adipose-specific chromatin marks in high-density lipoprotein cholesterol and triglycerides. Overlapping transcription factor (TF) binding sites with lipid-associated loci identifies TFs relevant in lipid biology. In addition, we present an integrative framework to prioritize causal variants at GWAS loci, producing a comprehensive list of candidate causal genes and variants with multiple layers of functional evidence. We highlight two of the prioritized genes, CREBRF and RRBP1, which show convergent evidence across functional datasets supporting their roles in lipid biology.

Abstract Background Whole genome sequencing is increasingly being used for the diagnosis of patients with rare diseases. However, the diagnostic yields of many studies, particularly those conducted in a healthcare setting, are often disappointingly low, at 25–30%. This is in part because although entire genomes are sequenced, analysis is often confined to in silico gene panels or coding regions of the genome. Methods We undertook WGS on a cohort of 122 unrelated rare disease patients and their relatives (300 genomes) who had been pre-screened by gene panels or arrays. Patients were recruited from a broad spectrum of clinical specialties. We applied a bioinformatics pipeline that would allow comprehensive analysis of all variant types. We combined established bioinformatics tools for phenotypic and genomic analysis with our novel algorithms (SVRare, ALTSPLICE and GREEN-DB) to detect and annotate structural, splice site and non-coding variants. Results Our diagnostic yield was 43/122 cases (35%), although 47/122 cases (39%) were considered solved when considering novel candidate genes with supporting functional data into account. Structural, splice site and deep intronic variants contributed to 20/47 (43%) of our solved cases. Five genes that are novel, or were novel at the time of discovery, were identified, whilst a further three genes are putative novel disease genes with evidence of causality. We identified variants of uncertain significance in a further fourteen candidate genes. The phenotypic spectrum associated with RMND1 was expanded to include polymicrogyria. Two patients with secondary findings in FBN1 and KCNQ1 were confirmed to have previously unidentified Marfan and long QT syndromes, respectively, and were referred for further clinical interventions. Clinical diagnoses were changed in six patients and treatment adjustments made for eight individuals, which for five patients was considered life-saving. Conclusions Genome sequencing is increasingly being considered as a first-line genetic test in routine clinical settings and can make a substantial contribution to rapidly identifying a causal aetiology for many patients, shortening their diagnostic odyssey. We have demonstrated that structural, splice site and intronic variants make a significant contribution to diagnostic yield and that comprehensive analysis of the entire genome is essential to maximise the value of clinical genome sequencing.

Increased blood lipid levels are heritable risk factors of cardiovascular disease with varied prevalence worldwide owing to different dietary patterns and medication use1. Despite advances in prevention and treatment, in particular through reducing low-density lipoprotein cholesterol levels2, heart disease remains the leading cause of death worldwide3. Genome-wideassociation studies (GWAS) of blood lipid levels have led to important biological and clinical insights, as well as new drug targets, for cardiovascular disease. However, most previous GWAS4-23 have been conducted in European ancestry populations and may have missed genetic variants that contribute to lipid-level variation in other ancestry groups. These include differences in allele frequencies, effect sizes and linkage-disequilibrium patterns24. Here we conduct a multi-ancestry, genome-wide genetic discovery meta-analysis of lipid levels in approximately 1.65 million individuals, including 350,000 of non-European ancestries. We quantify the gain in studying non-European ancestries and provide evidence to support the expansion of recruitment of additional ancestries, even with relatively small sample sizes. We find that increasing diversity rather than studying additional individuals of European ancestry results in substantial improvements in fine-mapping functional variants and portability of polygenic prediction (evaluated in approximately 295,000 individuals from 7 ancestry groupings). Modest gains in the number of discovered loci and ancestry-specific variants were also achieved. As GWAS expand emphasis beyond the identification of genes and fundamental biology towards the use of genetic variants for preventive and precision medicine25, we anticipate that increased diversity of participants will lead to more accurate and equitable26 application of polygenic scores in clinical practice.

IntroductionAs the clinical applications of medical genetics and genomics continue to expand, nongenetics professionals increasingly find themselves in the position of managing patients with genetic conditions. To prepare medical students to handle this future practice demand, it is imperative that they obtain skills and confidence in utilizing credible medical genetics resources to care for patients with genetic conditions. To this end, we developed active learning materials to introduce first-year medical students to these resources.MethodsThis approximately 2-hour session targeted first-year medical students (123 students) and combined flipped classroom and small-group collaborative case-based learning models. Students first completed a hands-on preclass exercise, which guided them in navigating the Online Mendelian Inheritance in Man website, and then attended an in-person small-group classroom activity, which provided the opportunity to apply information obtained from credible medical genetics resources to a patient case. At the conclusion of the classroom activity, students voluntarily completed an anonymous survey.ResultsResults of student postsession surveys showed that, regardless of previous exposure to medical genetics resources, this session increased both confidence in skills and future intention to use medical genetics resources.DiscussionSince the majority of students were unfamiliar with using specialized medical genetics resources prior to this educational intervention, the session functioned as a practical introduction to these essential resources. We propose that equipping medical students with skills that support inquiry-oriented learning, particularly in the early stage of training, can cultivate the practice of lifelong learning in medical genetics.

Inflammation-dependent base deaminases promote therapeutic resistance in many malignancies. However, their roles in human pre-leukemia stem cell (pre-LSC) evolution to acute myeloid leukemia stem cells (LSCs) had not been elucidated. Comparative whole-genome and whole-transcriptome sequencing analyses of FACS-purified pre-LSCs from myeloproliferative neoplasm (MPN) patients reveal APOBEC3C upregulation, an increased C-to-T mutational burden, and hematopoietic stem and progenitor cell (HSPC) proliferation during progression, which can be recapitulated by lentiviral APOBEC3C overexpression. In pre-LSCs, inflammatory splice isoform overexpression coincides with APOBEC3C upregulation and ADAR1p150-induced A-to-I RNA hyper-editing. Pre-LSC evolution to LSCs is marked by STAT3 editing, STAT3β isoform switching, elevated phospho-STAT3, and increased ADAR1p150 expression, which can be prevented by JAK2/STAT3 inhibition with ruxolitinib or fedratinib or lentiviral ADAR1 shRNA knockdown. Conversely, lentiviral ADAR1p150 expression enhances pre-LSC replating and STAT3 splice isoform switching. Thus, pre-LSC evolution to LSCs is fueled by primate-specific APOBEC3C-induced pre-LSC proliferation and ADAR1-mediated splicing deregulation.

Background: Fragile X-associated tremor/ataxia syndrome (FXTAS) is an adult-onset neurodegenerative disorder associated with premutation CGG-repeat expansions (55-200 repeats) in the 5' non-coding portion of the fragile X mental retardation 1 (FMR1) gene. Core features of FXTAS include progressive tremor/ataxia, cognitive decline, variable brain volume loss, and white matter disease. The principal histopathological feature of FXTAS is the presence of central nervous system (CNS) and non-CNS intranuclear inclusions. Objective: To further elucidate the molecular underpinnings of FXTAS through the proteomic characterization of human FXTAS cortexes. Results: Proteomic analysis of FXTAS brain cortical tissue (n = 8) identified minor differences in protein abundance compared to control brains (n = 6). Significant differences in FXTAS relative to control brain predominantly involved decreased abundance of proteins, with the greatest decreases observed for tenascin-C (TNC), cluster of differentiation 38 (CD38), and phosphoserine aminotransferase 1 (PSAT1); proteins typically increased in other neurodegenerative diseases. Proteins with the greatest increased abundance include potentially novel neurodegeneration-related proteins and small ubiquitin-like modifier 1/2 (SUMO1/2). The FMRpolyG peptide, proposed in models of FXTAS pathogenesis but only identified in trace amounts in the earlier study of FXTAS inclusions, was not identified in any of the FXTAS or control brains in the current study. Discussion: The observed proteomic shifts, while generally relatively modest, do show a bias toward decreased protein abundance with FXTAS. Such shifts in protein abundance also suggest altered RNA binding as well as loss of cell-cell adhesion/structural integrity. Unlike other neurodegenerative diseases, the proteome of end-stage FXTAS does not suggest a strong inflammation-mediated degenerative response.

Autosomal dominant hyper-IgE syndrome (AD-HIES) is typically caused by dominant-negative (DN) STAT3 mutations. Patients suffer from cold staphylococcal lesions and mucocutaneous candidiasis, severe allergy, and skeletal abnormalities. We report 12 patients from 8 unrelated kindreds with AD-HIES due to DN IL6ST mutations. We identified seven different truncating mutations, one of which was recurrent. The mutant alleles encode GP130 receptors bearing the transmembrane domain but lacking both the recycling motif and all four STAT3-recruiting tyrosine residues. Upon overexpression, the mutant proteins accumulate at the cell surface and are loss of function and DN for cellular responses to IL-6, IL-11, LIF, and OSM. Moreover, the patients' heterozygous leukocytes and fibroblasts respond poorly to IL-6 and IL-11. Consistently, patients with STAT3 and IL6ST mutations display infectious and allergic manifestations of IL-6R deficiency, and some of the skeletal abnormalities of IL-11R deficiency. DN STAT3 and IL6ST mutations thus appear to underlie clinical phenocopies through impairment of the IL-6 and IL-11 response pathways.

ObjectivesThe challenges of understanding how interventions influence follow-up medical care are magnified during genomic testing because few patients have received it to date and because the scope of information it provides is complex and often unexpected. We tested a novel strategy for quantifying downstream healthcare utilization after genomic testing to more comprehensively and efficiently identify related services. We also evaluated the effectiveness of different methods for collecting these data.MethodsWe developed a risk-based approach for a trial of newborn genomic sequencing in which we defined primary conditions based on existing diagnoses and family histories of disease and defined secondary conditions based on unexpected findings. We then created patient-specific lists of services associated with managing primary and secondary conditions. Services were quantified based on medical record reviews, surveys, and telephone check-ins with parents.ResultsBy focusing on services that genomic testing would most likely influence in the short-term, we reduced the number of services in our analyses by more than 90% compared with analyses of all observed services. We also identified the same services that were ordered in response to unexpected findings as were identified during expert review and by confirming whether recommendations were completed. Data also showed that quantifying healthcare utilization with surveys and telephone check-ins alone would have missed the majority of attributable services.ConclusionsOur risk-based strategy provides an improved approach for assessing the short-term impact of genomic testing and other interventions on healthcare utilization while conforming as much as possible to existing best-practice recommendations.

BackgroundThe Critical Assessment of Functional Annotation (CAFA) is an ongoing, global, community-driven effort to evaluate and improve the computational annotation of protein function.ResultsHere, we report on the results of the third CAFA challenge, CAFA3, that featured an expanded analysis over the previous CAFA rounds, both in terms of volume of data analyzed and the types of analysis performed. In a novel and major new development, computational predictions and assessment goals drove some of the experimental assays, resulting in new functional annotations for more than 1000 genes. Specifically, we performed experimental whole-genome mutation screening in Candida albicans and Pseudomonas aureginosa genomes, which provided us with genome-wide experimental data for genes associated with biofilm formation and motility. We further performed targeted assays on selected genes in Drosophila melanogaster, which we suspected of being involved in long-term memory.ConclusionWe conclude that while predictions of the molecular function and biological process annotations have slightly improved over time, those of the cellular component have not. Term-centric prediction of experimental annotations remains equally challenging; although the performance of the top methods is significantly better than the expectations set by baseline methods in C. albicans and D. melanogaster, it leaves considerable room and need for improvement. Finally, we report that the CAFA community now involves a broad range of participants with expertise in bioinformatics, biological experimentation, biocuration, and bio-ontologies, working together to improve functional annotation, computational function prediction, and our ability to manage big data in the era of large experimental screens.

RNA therapeutics comprise a diverse group of oligonucleotide-based drugs such as antisense oligonucleotides (ASOs), small interfering RNAs (siRNAs), and short hairpin RNAs (shRNAs) that can be designed to selectively interact with drug targets currently undruggable with small molecule-based drugs or monoclonal antibodies. Furthermore, RNA-based therapeutics have the potential to modulate entire disease pathways, and thereby represent a new modality with unprecedented potential for generating disease-modifying drugs for a wide variety of human diseases, including central nervous system (CNS) disorders. Here, we describe different strategies for delivering RNA drugs to the CNS and review recent advances in clinical development of ASO drugs and siRNA-based therapeutics for the treatment of neurological diseases and neuromuscular disorders. Abbreviations 2’-MOE: 2’-O-(2-methoxyethyl); 2’-O-Me: 2’-O-methyl; 2’-F: 2’-fluoro; AD: Alzheimer's disease; ALS: Amyotrophic lateral sclerosis; ALSFRS-R: Revised Amyotrophic Lateral Sclerosis Functional Rating Scale; ARC: Antibody siRNA Conjugate; AS: Angelman Syndrome; ASGRP: Asialoglycoprotein receptor; ASO: Antisense oligonucleotide; AxD: Alexander Disease; BBB: Blood brain barrier; Bp: Basepair; CNM: Centronuclear myopathies; CNS: Central Nervous System; CPP: Cell-penetrating Peptide; CSF: Cerebrospinal fluid; DMD: Duchenne muscular dystrophy; DNA: Deoxyribonucleic acid; FAP: Familial amyloid polyneuropathy; FALS: Familial amyotrophic lateral sclerosis; FDA: The United States Food and Drug Administration; GalNAc: N-acetylgalactosamine; GoF: Gain of function; hATTR: Hereditary transthyretin amyloidosis; HD: Huntington's disease; HRQOL: health-related quality of life; ICV: Intracerebroventricular; IT: Intrathecal; LNA: Locked nucleic acid; LoF: Loss of function; mRNA: Messenger RNA; MS: Multiple Sclerosis; MSA: Multiple System Atrophy; NBE: New Biological Entity; NCE: New Chemical Entity; NHP: Nonhuman primate; nt: Nucleotide; PD: Parkinson's disease; PNP: Polyneuropathy; PNS: Peripheral nervous system; PS: Phosphorothioate; RISC: RNA-Induced Silencing Complex; RNA: Ribonucleic acid; RNAi: RNA interference; s.c.: Subcutaneous; siRNA: Small interfering RNA; SMA: Spinal muscular atrophy; SMN: Survival motor neuron; TTR: Transthyretin

Encoding the slow skeletal muscle isoform of myosin binding protein-C, MYBPC1 is associated with autosomal dominant and recessive forms of arthrogryposis. The authors describe a novel association for MYBPC1 in four patients from three independent families with skeletal muscle weakness, myogenic tremors, and hypotonia with gradual clinical improvement. The patients carried one of two de novo heterozygous variants in MYBPC1, with the p.Leu263Arg variant seen in three individuals and the p.Leu259Pro variant in one individual. Both variants are absent from controls, well conserved across vertebrate species, predicted to be damaging, and located in the M-motif. Protein modeling studies suggested that the p.Leu263Arg variant affects the stability of the M-motif, whereas the p.Leu259Pro variant alters its structure. In vitro biochemical and kinetic studies demonstrated that the p.Leu263Arg variant results in decreased binding of the M-motif to myosin, which likely impairs the formation of actomyosin cross-bridges during muscle contraction. Collectively, our data substantiate that damaging variants in MYBPC1 are associated with a new form of an early-onset myopathy with tremor, which is a defining and consistent characteristic in all affected individuals, with no contractures. Recognition of this expanded myopathic phenotype can enable identification of individuals with MYBPC1 variants without arthrogryposis.

BackgroundFamily screening of a 48-year-old male with recently diagnosed IgG4-related disease (IgG4-RD) revealed unanticipated elevations in plasma IgG4 in his two healthy teenaged sons.MethodsWe performed gene sequencing, immune cell studies, HLA typing, and analyses of circulating cytotoxic CD4+ T lymphocytes and plasmablasts to seek clues to pathogenesis. DNA from a separate cohort of 99 patients with known IgG4-RD was also sequenced for the presence of genetic variants in a specific gene, FGFBP2.ResultsThe three share a previously unreported heterozygous single base deletion in fibroblast growth factor binding protein type 2 (FGFBP2), which causes a frameshift in the coding sequence. The FGFBP2 protein is secreted by cytotoxic T-lymphocytes and binds fibroblast growth factor. The variant sequence in the FGFBP2 protein is predicted to form a disordered random coil rather than a helical-turn-helix structure, unable to adopt a stable conformation. The proband and the two sons had 5-10-fold higher numbers of circulating cytotoxic CD4 + T cells and plasmablasts compared to matched controls. The three members also share a homozygous missense common variant in FGFBP2 found in heterozygous form in ~40% of the population. This common variant was found in 73% of an independent, well characterized IgG4-RD cohort, showing enrichment in idiopathic IgG4-RD.ConclusionsThe presence of a shared deleterious variant and homozygous common variant in FGFBP2 in the proband and sons strongly implicates this cytotoxic T cell product in the pathophysiology of IgG4-RD. The high prevalence of a common FGFBP2 variant in sporadic IgG4-RD supports the likelihood of participation in disease.

Adenosine deaminase associated with RNA1 (ADAR1) deregulation contributes to therapeutic resistance in many malignancies. Here we show that ADAR1-induced hyper-editing in normal human hematopoietic progenitors impairs miR-26a maturation, which represses CDKN1A expression indirectly via EZH2, thereby accelerating cell-cycle transit. However, in blast crisis chronic myeloid leukemia progenitors, loss of EZH2 expression and increased CDKN1A oppose cell-cycle transit. Moreover, A-to-I editing of both the MDM2 regulatory microRNA and its binding site within the 3' UTR region stabilizes MDM2 transcripts, thereby enhancing blast crisis progenitor propagation. These data reveal a dual mechanism governing malignant transformation of progenitors that is predicated on hyper-editing of cell-cycle-regulatory miRNAs and the 3' UTR binding site of tumor suppressor miRNAs.

Increased levels of the urinary albumin-to-creatinine ratio (UACR) are associated with higher risk of kidney disease progression and cardiovascular events, but underlying mechanisms are incompletely understood. Here, we conduct trans-ethnic (n = 564,257) and European-ancestry specific meta-analyses of genome-wide association studies of UACR, including ancestry- and diabetes-specific analyses, and identify 68 UACR-associated loci. Genetic correlation analyses and risk score associations in an independent electronic medical records database (n = 192,868) reveal connections with proteinuria, hyperlipidemia, gout, and hypertension. Fine-mapping and trans-Omics analyses with gene expression in 47 tissues and plasma protein levels implicate genes potentially operating through differential expression in kidney (including TGFB1, MUC1, PRKCI, and OAF), and allow coupling of UACR associations to altered plasma OAF concentrations. Knockdown of OAF and PRKCI orthologs in Drosophila nephrocytes reduces albumin endocytosis. Silencing fly PRKCI further impairs slit diaphragm formation. These results generate a priority list of genes and pathways for translational research to reduce albuminuria.

Small nucleolar RNAs (snoRNAs) are a diverse group of non-coding RNAs that direct chemical modifications at specific residues on other RNA molecules, primarily on ribosomal RNA (rRNA). SnoRNAs are altered in several cancers; however, their role in cell homeostasis as well as in cellular transformation remains poorly explored. Here, we show that specific subsets of snoRNAs are differentially regulated during the earliest cellular response to oncogenic RASG12V expression. We describe a novel function for one H/ACA snoRNA, SNORA24, which guides two pseudouridine modifications within the small ribosomal subunit, in RAS-induced senescence in vivo. We find that in mouse models, loss of Snora24 cooperates with RASG12V to promote the development of liver cancer that closely resembles human steatohepatitic hepatocellular carcinoma (HCC). From a clinical perspective, we further show that human HCCs with low SNORA24 expression display increased lipid content and are associated with poor patient survival. We next asked whether ribosomes lacking SNORA24-guided pseudouridine modifications on 18S rRNA have alterations in their biophysical properties. Single-molecule Fluorescence Resonance Energy Transfer (FRET) analyses revealed that these ribosomes exhibit perturbations in aminoacyl-transfer RNA (aa-tRNA) selection and altered pre-translocation ribosome complex dynamics. Furthermore, we find that HCC cells lacking SNORA24-guided pseudouridine modifications have increased translational miscoding and stop codon readthrough frequencies. These findings highlight a role for specific snoRNAs in safeguarding against oncogenic insult and demonstrate a functional link between H/ACA snoRNAs regulated by RAS and the biophysical properties of ribosomes in cancer.

ObjectivePrimary sclerosing cholangitis (PSC) is a genetically complex, inflammatory bile duct disease of largely unknown aetiology often leading to liver transplantation or death. Little is known about the genetic contribution to the severity and progression of PSC. The aim of this study is to identify genetic variants associated with PSC disease progression and development of complications.DesignWe collected standardised PSC subphenotypes in a large cohort of 3402 patients with PSC. After quality control, we combined 130 422 single nucleotide polymorphisms of all patients-obtained using the Illumina immunochip-with their disease subphenotypes. Using logistic regression and Cox proportional hazards models, we identified genetic variants associated with binary and time-to-event PSC subphenotypes.ResultsWe identified genetic variant rs853974 to be associated with liver transplant-free survival (p=6.07×10-9). Kaplan-Meier survival analysis showed a 50.9% (95% CI 41.5% to 59.5%) transplant-free survival for homozygous AA allele carriers of rs853974 compared with 72.8% (95% CI 69.6% to 75.7%) for GG carriers at 10 years after PSC diagnosis. For the candidate gene in the region, RSPO3, we demonstrated expression in key liver-resident effector cells, such as human and murine cholangiocytes and human hepatic stellate cells.ConclusionWe present a large international PSC cohort, and report genetic loci associated with PSC disease progression. For liver transplant-free survival, we identified a genome-wide significant signal and demonstrated expression of the candidate gene RSPO3 in key liver-resident effector cells. This warrants further assessments of the role of this potential key PSC modifier gene.

Angiopoietin-like 4 (ANGPTL4) is an endogenous inhibitor of lipoprotein lipase that modulates lipid levels, coronary atherosclerosis risk, and nutrient partitioning. We hypothesize that loss of ANGPTL4 function might improve glucose homeostasis and decrease risk of type 2 diabetes (T2D). We investigate protein-altering variants in ANGPTL4 among 58,124 participants in the DiscovEHR human genetics study, with follow-up studies in 82,766 T2D cases and 498,761 controls. Carriers of p.E40K, a variant that abolishes ANGPTL4 ability to inhibit lipoprotein lipase, have lower odds of T2D (odds ratio 0.89, 95% confidence interval 0.85-0.92, p = 6.3 × 10-10), lower fasting glucose, and greater insulin sensitivity. Predicted loss-of-function variants are associated with lower odds of T2D among 32,015 cases and 84,006 controls (odds ratio 0.71, 95% confidence interval 0.49-0.99, p = 0.041). Functional studies in Angptl4-deficient mice confirm improved insulin sensitivity and glucose homeostasis. In conclusion, genetic inactivation of ANGPTL4 is associated with improved glucose homeostasis and reduced risk of T2D.

Electrocardiographic PR interval measures atrio-ventricular depolarization and conduction, and abnormal PR interval is a risk factor for atrial fibrillation and heart block. Our genome-wide association study of over 92,000 European-descent individuals identifies 44 PR interval loci (34 novel). Examination of these loci reveals known and previously not-yet-reported biological processes involved in cardiac atrial electrical activity. Genes in these loci are over-represented in cardiac disease processes including heart block and atrial fibrillation. Variants in over half of the 44 loci were associated with atrial or blood transcript expression levels, or were in high linkage disequilibrium with missense variants. Six additional loci were identified either by meta-analysis of ~105,000 African and European-descent individuals and/or by pleiotropic analyses combining PR interval with heart rate, QRS interval, and atrial fibrillation. These findings implicate developmental pathways, and identify transcription factors, ion-channel genes, and cell-junction/cell-signaling proteins in atrio-ventricular conduction, identifying potential targets for drug development.

Primary sclerosing cholangitis (PSC) is strongly associated with several human leukocyte antigen (HLA) haplotypes. Due to extensive linkage disequilibrium and multiple polymorphic candidate genes in the HLA complex, identifying the alleles responsible for these associations has proven difficult. We aimed to evaluate whether studying populations of admixed or non-European descent could help in defining the causative HLA alleles. When assessing haplotypes carrying HLA-DRB1*13:01 (hypothesized to specifically increase the susceptibility to chronic cholangitis), we observed that every haplotype in the Scandinavian PSC population carried HLA-DQB1*06:03. In contrast, only 65% of HLA-DRB1*13:01 haplotypes in an admixed/non-European PSC population carried this allele, suggesting that further assessments of the PSC-associated haplotype HLA-DRB1*13:01-DQA1*01:03-DQB1*06:03 in admixed or multi-ethnic populations could aid in identifying the causative allele.

BackgroundEnterohepatic Helicobacter-like organisms (HLO) have been recognized as causes of human infection since 1984, primarily as a cause of bacteremia and cellulitis in immunocompromised individuals, but the spectrum of illness due to HLO has expanded based on numerous reports from Japan.MethodsWe report 4 epidemiologically unrelated immunocompetent patients with HLO bacteremia diagnosed within a 2-year period. Three patients had cellulitis and 1 patient had unexplained fever. 16S ribosomal deoxyribonucleic acid (rDNA) sequence analysis of 2 isolates suggested that they were Helicobacter cinaedi, and whole-genome sequencing showed that they differed only slightly from reference strains.ConclusionsWe believe that this syndrome is more common than reported, but it is easily overlooked because the skin lesions resemble streptococcal cellulitis and respond very rapidly to β-lactam antibiotics, and the organism is difficult to isolate from the blood. All HLO in our series were isolated from blood using the ESP system and were not detected in 2 other widely used commercial blood culture systems.

Exfoliation syndrome (XFS) is the most common known risk factor for secondary glaucoma and a major cause of blindness worldwide. Variants in two genes, LOXL1 and CACNA1A, have previously been associated with XFS. To further elucidate the genetic basis of XFS, we collected a global sample of XFS cases to refine the association at LOXL1, which previously showed inconsistent results across populations, and to identify new variants associated with XFS. We identified a rare protective allele at LOXL1 (p.Phe407, odds ratio (OR) = 25, P = 2.9 × 10-14) through deep resequencing of XFS cases and controls from nine countries. A genome-wide association study (GWAS) of XFS cases and controls from 24 countries followed by replication in 18 countries identified seven genome-wide significant loci (P < 5 × 10-8). We identified association signals at 13q12 (POMP), 11q23.3 (TMEM136), 6p21 (AGPAT1), 3p24 (RBMS3) and 5q23 (near SEMA6A). These findings provide biological insights into the pathology of XFS and highlight a potential role for naturally occurring rare LOXL1 variants in disease biology.

The rapid development of genomic sequencing technologies has decreased the cost of genetic analysis to the extent that it seems plausible that genome-scale sequencing could have widespread availability in pediatric care. Genomic sequencing provides a powerful diagnostic modality for patients who manifest symptoms of monogenic disease and an opportunity to detect health conditions before their development. However, many technical, clinical, ethical, and societal challenges should be addressed before such technology is widely deployed in pediatric practice. This article provides an overview of the Newborn Sequencing in Genomic Medicine and Public Health Consortium, which is investigating the application of genome-scale sequencing in newborns for both diagnosis and screening.

BackgroundA major bottleneck in our understanding of the molecular underpinnings of life is the assignment of function to proteins. While molecular experiments provide the most reliable annotation of proteins, their relatively low throughput and restricted purview have led to an increasing role for computational function prediction. However, assessing methods for protein function prediction and tracking progress in the field remain challenging.ResultsWe conducted the second critical assessment of functional annotation (CAFA), a timed challenge to assess computational methods that automatically assign protein function. We evaluated 126 methods from 56 research groups for their ability to predict biological functions using Gene Ontology and gene-disease associations using Human Phenotype Ontology on a set of 3681 proteins from 18 species. CAFA2 featured expanded analysis compared with CAFA1, with regards to data set size, variety, and assessment metrics. To review progress in the field, the analysis compared the best methods from CAFA1 to those of CAFA2.ConclusionsThe top-performing methods in CAFA2 outperformed those from CAFA1. This increased accuracy can be attributed to a combination of the growing number of experimental annotations and improved methods for function prediction. The assessment also revealed that the definition of top-performing algorithms is ontology specific, that different performance metrics can be used to probe the nature of accurate predictions, and the relative diversity of predictions in the biological process and human phenotype ontologies. While there was methodological improvement between CAFA1 and CAFA2, the interpretation of results and usefulness of individual methods remain context-dependent.

The epigenome orchestrates genome accessibility, functionality, and three-dimensional structure. Because epigenetic variation can impact transcription and thus phenotypes, it may contribute to adaptation. Here, we report 1,107 high-quality single-base resolution methylomes and 1,203 transcriptomes from the 1001 Genomes collection of Arabidopsis thaliana. Although the genetic basis of methylation variation is highly complex, geographic origin is a major predictor of genome-wide DNA methylation levels and of altered gene expression caused by epialleles. Comparison to cistrome and epicistrome datasets identifies associations between transcription factor binding sites, methylation, nucleotide variation, and co-expression modules. Physical maps for nine of the most diverse genomes reveal how transposons and other structural variants shape the epigenome, with dramatic effects on immunity genes. The 1001 Epigenomes Project provides a comprehensive resource for understanding how variation in DNA methylation contributes to molecular and non-molecular phenotypes in natural populations of the most studied model plant.

Formative research suggests that a human embryonic stem cell-specific alternative splicing gene regulatory network, which is repressed by Muscleblind-like (MBNL) RNA binding proteins, is involved in cell reprogramming. In this study, RNA sequencing, splice isoform-specific quantitative RT-PCR, lentiviral transduction, and in vivo humanized mouse model studies demonstrated that malignant reprogramming of progenitors into self-renewing blast crisis chronic myeloid leukemia stem cells (BC LSCs) was partially driven by decreased MBNL3. Lentiviral knockdown of MBNL3 resulted in reversion to an embryonic alternative splice isoform program typified by overexpression of CD44 transcript variant 3, containing variant exons 8-10, and BC LSC proliferation. Although isoform-specific lentiviral CD44v3 overexpression enhanced chronic phase chronic myeloid leukemia (CML) progenitor replating capacity, lentiviral shRNA knockdown abrogated these effects. Combined treatment with a humanized pan-CD44 monoclonal antibody and a breakpoint cluster region - ABL proto-oncogene 1, nonreceptor tyrosine kinase (BCR-ABL1) antagonist inhibited LSC maintenance in a niche-dependent manner. In summary, MBNL3 down-regulation-related reversion to an embryonic alternative splicing program, typified by CD44v3 overexpression, represents a previously unidentified mechanism governing malignant progenitor reprogramming in malignant microenvironments and provides a pivotal opportunity for selective BC LSC detection and therapeutic elimination.

The Undiagnosed Diseases Network (UDN) builds on the successes of the Undiagnosed Diseases Program at the National Institutes of Health (NIH UDP). Through support from the NIH Common Fund, a coordinating center, six additional clinical sites, and two sequencing cores comprise the UDN. The objectives of the UDN are to: (1) improve the level of diagnosis and care for patients with undiagnosed diseases through the development of common protocols designed by an enlarged community of investigators across the network; (2) facilitate research into the etiology of undiagnosed diseases, by collecting and sharing standardized, high-quality clinical and laboratory data including genotyping, phenotyping, and environmental exposure data; and (3) create an integrated and collaborative research community across multiple clinical sites, and among laboratory and clinical investigators, to investigate the pathophysiology of these rare diseases and to identify options for patient management. Broad-based data sharing is at the core of achieving these objectives, and the UDN is establishing the policies and governance structure to support broad data sharing.

There are few better examples of the need for data sharing than in the rare disease community, where patients, physicians, and researchers must search for "the needle in a haystack" to uncover rare, novel causes of disease within the genome. Impeding the pace of discovery has been the existence of many small siloed datasets within individual research or clinical laboratory databases and/or disease-specific organizations, hoping for serendipitous occasions when two distant investigators happen to learn they have a rare phenotype in common and can "match" these cases to build evidence for causality. However, serendipity has never proven to be a reliable or scalable approach in science. As such, the Matchmaker Exchange (MME) was launched to provide a robust and systematic approach to rare disease gene discovery through the creation of a federated network connecting databases of genotypes and rare phenotypes using a common application programming interface (API). The core building blocks of the MME have been defined and assembled. Three MME services have now been connected through the API and are available for community use. Additional databases that support internal matching are anticipated to join the MME network as it continues to grow.

There is increasing evidence from genome-wide association studies for a strong inherited genetic basis of susceptibility to acute lymphoblastic leukaemia (ALL) in children, yet the effects of protein-coding variants on ALL risk have not been systematically evaluated. Here we show a missense variant in CDKN2A associated with the development of ALL at genome-wide significance (rs3731249, P=9.4 × 10(-23), odds ratio=2.23). Functional studies indicate that this hypomorphic variant results in reduced tumour suppressor function of p16(INK4A), increases the susceptibility to leukaemic transformation of haematopoietic progenitor cells, and is preferentially retained in ALL tumour cells. Resequencing the CDKN2A-CDKN2B locus in 2,407 childhood ALL cases reveals 19 additional putative functional germline variants. These results provide direct functional evidence for the influence of inherited genetic variation on ALL risk, highlighting the important and complex roles of CDKN2A-CDKN2B tumour suppressors in leukaemogenesis.

Levels of low-density lipoprotein (LDL) cholesterol, high-density lipoprotein (HDL) cholesterol, triglycerides and total cholesterol are heritable, modifiable risk factors for coronary artery disease. To identify new loci and refine known loci influencing these lipids, we examined 188,577 individuals using genome-wide and custom genotyping arrays. We identify and annotate 157 loci associated with lipid levels at P < 5 × 10(-8), including 62 loci not previously associated with lipid levels in humans. Using dense genotyping in individuals of European, East Asian, South Asian and African ancestry, we narrow association signals in 12 loci. We find that loci associated with blood lipid levels are often associated with cardiovascular and metabolic traits, including coronary artery disease, type 2 diabetes, blood pressure, waist-hip ratio and body mass index. Our results demonstrate the value of using genetic data from individuals of diverse ancestry and provide insights into the biological mechanisms regulating blood lipids to guide future genetic, biological and therapeutic research.

UnlabelledMalassezia commensal yeasts are associated with a number of skin disorders, such as atopic eczema/dermatitis and dandruff, and they also can cause systemic infections. Here we describe the 7.67-Mbp genome of Malassezia sympodialis, a species associated with atopic eczema, and contrast its genome repertoire with that of Malassezia globosa, associated with dandruff, as well as those of other closely related fungi. Ninety percent of the predicted M. sympodialis protein coding genes were experimentally verified by mass spectrometry at the protein level. We identified a relatively limited number of genes related to lipid biosynthesis, and both species lack the fatty acid synthase gene, in line with the known requirement of these yeasts to assimilate lipids from the host. Malassezia species do not appear to have many cell wall-localized glycosylphosphatidylinositol (GPI) proteins and lack other cell wall proteins previously identified in other fungi. This is surprising given that in other fungi these proteins have been shown to mediate interactions (e.g., adhesion and biofilm formation) with the host. The genome revealed a complex evolutionary history for an allergen of unknown function, Mala s 7, shown to be encoded by a member of an amplified gene family of secreted proteins. Based on genetic and biochemical studies with the basidiomycete human fungal pathogen Cryptococcus neoformans, we characterized the allergen Mala s 6 as the cytoplasmic cyclophilin A. We further present evidence that M. sympodialis may have the capacity to undergo sexual reproduction and present a model for a pseudobipolar mating system that allows limited recombination between two linked MAT loci.ImportanceMalassezia commensal yeasts are associated with a number of skin disorders. The previously published genome of M. globosa provided some of the first insights into Malassezia biology and its involvement in dandruff. Here, we present the genome of M. sympodialis, frequently isolated from patients with atopic eczema and healthy individuals. We combined comparative genomics with sequencing and functional characterization of specific genes in a population of clinical isolates and in closely related model systems. Our analyses provide insights into the evolution of allergens related to atopic eczema and the evolutionary trajectory of the machinery for sexual reproduction and meiosis. We hypothesize that M. sympodialis may undergo sexual reproduction, which has important implications for the understanding of the life cycle and virulence potential of this medically important yeast. Our findings provide a foundation for the development of genetic and genomic tools to elucidate host-microbe interactions that occur on the skin and to identify potential therapeutic targets.

BackgroundThe chemokine RANTES (regulated on activation, normal T-cell expressed and secreted)/CCL5 is involved in the pathogenesis of cardiovascular disease in mice, whereas less is known in humans. We hypothesised that its relevance for atherosclerosis should be reflected by associations between CCL5 gene variants, RANTES serum concentrations and protein levels in atherosclerotic plaques and risk for coronary events.Methods and findingsWe conducted a case-cohort study within the population-based MONICA/KORA Augsburg studies. Baseline RANTES serum levels were measured in 363 individuals with incident coronary events and 1,908 non-cases (mean follow-up: 10.2±4.8 years). Cox proportional hazard models adjusting for age, sex, body mass index, metabolic factors and lifestyle factors revealed no significant association between RANTES and incident coronary events (HR [95% CI] for increasing RANTES tertiles 1.0, 1.03 [0.75-1.42] and 1.11 [0.81-1.54]). None of six CCL5 single nucleotide polymorphisms and no common haplotype showed significant associations with coronary events. Also in the CARDIoGRAM study (>22,000 cases, >60,000 controls), none of these CCL5 SNPs was significantly associated with coronary artery disease. In the prospective Athero-Express biobank study, RANTES plaque levels were measured in 606 atherosclerotic lesions from patients who underwent carotid endarterectomy. RANTES content in atherosclerotic plaques was positively associated with macrophage infiltration and inversely associated with plaque calcification. However, there was no significant association between RANTES content in plaques and risk for coronary events (mean follow-up 2.8±0.8 years).ConclusionsHigh RANTES plaque levels were associated with an unstable plaque phenotype. However, the absence of associations between (i) RANTES serum levels, (ii) CCL5 genotypes and (iii) RANTES content in carotid plaques and either coronary artery disease or incident coronary events in our cohorts suggests that RANTES may not be a novel coronary risk biomarker. However, the potential relevance of RANTES levels in platelet-poor plasma needs to be investigated in further studies.

Through an international consortium, we have collected 37 tau- and TAR DNA-binding protein 43 (TDP-43)-negative frontotemporal lobar degeneration (FTLD) cases, and present here the first comprehensive analysis of these cases in terms of neuropathology, genetics, demographics and clinical data. 92% (34/37) had fused in sarcoma (FUS) protein pathology, indicating that FTLD-FUS is an important FTLD subtype. This FTLD-FUS collection specifically focussed on aFTLD-U cases, one of three recently defined subtypes of FTLD-FUS. The aFTLD-U subtype of FTLD-FUS is characterised clinically by behavioural variant frontotemporal dementia (bvFTD) and has a particularly young age of onset with a mean of 41 years. Further, this subtype had a high prevalence of psychotic symptoms (36% of cases) and low prevalence of motor symptoms (3% of cases). We did not find FUS mutations in any aFTLD-U case. To date, the only subtype of cases reported to have ubiquitin-positive but tau-, TDP-43- and FUS-negative pathology, termed FTLD-UPS, is the result of charged multivesicular body protein 2B gene (CHMP2B) mutation. We identified three FTLD-UPS cases, which are negative for CHMP2B mutation, suggesting that the full complement of FTLD pathologies is yet to be elucidated.