Ultrafast and scalable variant annotation and prioritization with big functional genomics data

Enhanced interactions among gut mycobiomes with the deterioration of glycemic control

BACKGROUND

The gut mycobiome is closely linked to health and disease; however, its role in the progression of type 2 diabetes mellitus (T2DM) remains obscure. Here, a multi-omics approach was employed to explore the role of intestinal fungi in the deterioration of glycemic control.

METHODS

350 participants without hypoglycemic therapies were invited for a standard oral glucose tolerance test to determine their status of glycemic control. The gut mycobiome was identified through internal transcribed spacer sequencing, host genetics were determined by genotyping array, and plasma metabolites were measured with untargeted liquid chromatography mass spectrometry.

FINDINGS

The richness of fungi was higher, whereas its dissimilarity was markedly lower, in participants with T2DM. Moreover, the diversity and composition of fungi were closely associated with insulin sensitivity and pancreatic β-cell functions. With the exacerbation of glycemic control, the co-occurrence network among fungus taxa became increasingly complex, and the complexity of the interaction network was inversely associated with insulin sensitivity. Mendelian randomization analysis further demonstrated that the Archaeorhizomycetes class, Fusarium genus, and Neoascochyta genus were causally linked to impaired glucose metabolism. Furthermore, integrative analysis with metabolomics showed that increased 4-hydroxy-2-oxoglutaric acid, ketoleucine, lysophosphatidylcholine (20:3/0:0), and N-lactoyl-phenylalanine, but decreased lysophosphatidylcholine (O-18:2), functioned as key molecules linking the adverse effect of Fusarium genus on insulin sensitivity.

CONCLUSIONS

Our study uncovers a strong association between disturbance in gut fungi and the progression of T2DM and highlights the potential of targeting the gut mycobiome for the management of T2DM.

FUNDINGS

This study was supported by MOST and NSFC of China.

Genome-Wide Methylation Analysis Reveals a KCNK3-Prominent Causal Cascade on Hypertension

BACKGROUND

Despite advances in understanding hypertension's genetic structure, how noncoding genetic variants influence it remains unclear. Studying their interaction with DNA methylation is crucial to deciphering this complex disease's genetic mechanisms.

METHODS

We investigated the genetic and epigenetic interplay in hypertension using whole-genome bisulfite sequencing. Methylation profiling in 918 males revealed allele-specific methylation and methylation quantitative trait loci. We engineered rs1275988T/C mutant mice using CRISPR (clustered regularly interspaced short palindromic repeats)/Cas9 (CRISPR-associated protein 9), bred them for homozygosity, and subjected them to a high-salt diet. Telemetry captured their cardiovascular metrics. Protein-DNA interactions were elucidated using DNA pull-downs, mass spectrometry, and Western blots. A wire myograph assessed vascular function, and analysis of the Kcnk3 gene methylation highlighted the mutation's role in hypertension.

RESULTS

We discovered that DNA methylation-associated genetic effects, especially in non-cytosine-phosphate-guanine (non-CpG) island and noncoding distal regulatory regions, significantly contribute to hypertension predisposition. We identified distinct methylation quantitative trait locus patterns in the hypertensive population and observed that the onset of hypertension is influenced by the transmission of genetic effects through the demethylation process. By evidence-driven prioritization and in vivo experiments, we unearthed rs1275988 in a cell type-specific enhancer as a notable hypertension causal variant, intensifying hypertension through the modulation of local DNA methylation and consequential alterations in Kcnk3 gene expression and vascular remodeling. When exposed to a high-salt diet, mice with the rs1275988C/C genotype exhibited exacerbated hypertension and significant vascular remodeling, underscored by increased aortic wall thickness. The C allele of rs1275988 was associated with elevated DNA methylation levels, driving down the expression of the Kcnk3 gene by attenuating Nr2f2 (nuclear receptor subfamily 2 group F member 2) binding at the enhancer locus.

CONCLUSIONS

Our research reveals new insights into the complex interplay between genetic variations and DNA methylation in hypertension. We underscore hypomethylation's potential in hypertension onset and identify rs1275988 as a causal variant in vascular remodeling. This work advances our understanding of hypertension's molecular mechanisms and encourages personalized health care strategies.

TOPMed imputed genomics enhances genomic atlas of the human proteome in brain, cerebrospinal fluid, and plasma

Comprehensive expression quantitative trait loci studies have been instrumental for understanding tissue-specific gene regulation and pinpointing functional genes for disease-associated loci in a tissue-specific manner. Compared to gene expressions, proteins more directly affect various biological processes, often dysregulated in disease, and are important drug targets. We previously performed and identified tissue-specific protein quantitative trait loci in brain, cerebrospinal fluid, and plasma. We now enhance this work by analyzing more proteins (1,300 versus 1,079) and an almost twofold increase in high quality imputed genetic variants (8.4 million versus 4.4 million) by using TOPMed reference panel. We identified 38 genomic regions associated with 43 proteins in brain, 150 regions associated with 247 proteins in cerebrospinal fluid, and 95 regions associated with 145 proteins in plasma. Compared to our previous study, this study newly identified 12 loci in brain, 30 loci in cerebrospinal fluid, and 22 loci in plasma. Our improved genomic atlas uncovers the genetic control of protein regulation across multiple tissues. These resources are accessible through the Online Neurodegenerative Trait Integrative Multi-Omics Explorer for use by the scientific community.

Epilepsy and psychiatric comorbidities: A bidirectional mendelian randomization study

STUDY OBJECTIVES

Psychiatric comorbidities are relatively common among patients with epilepsy; however, the underlying mechanisms of this association remain largely unknown. The objective of this Mendelian randomization (MR) study was to analyze the genetic correlations and causality underlying these reciprocal associations.

METHODS

Single-nucleotide polymorphisms associated with epilepsy (29,677 controls and 15,212 cases) and seven psychiatric comorbidities (485,436 controls and 269,495 cases) were identified from genome-wide association studies. Causal significance was estimated using inverse variance weighting. Sensitivity analyses included the weighted median, MR-Egger, and MR-PRESSO. The psychiatric comorbidities analyzed in this study included attention deficit hyperactivity disorder (ADHD), autism spectrum disorder, major depressive disorder, bipolar disorder, schizophrenia, obsessive-compulsive disorder (OCD), and anorexia nervosa.

RESULTS

Both forward and reverse genetic associations were observed for the selected psychiatric disorders. Notably, ADHD was significantly associated with an increased risk of generalized epilepsy (odds ratio [OR], 1.09; 95 % confidence interval [CI], 1.01-1.18; p = 0.013). However, MR-PRESSO detected the existence of pleiotropy (p = 0.001). Additionally, focal epilepsy was significantly associated with a higher risk of OCD (OR, 1.44; 95 % CI, 1.08-1.92; p = 0.013), and all sensitivity tests yielded favorably nonsignificant results. There was no significant genetic association between epilepsy and other examined psychiatric disorders. However, due to the detection of pleiotropy by MR-Egger and considerations related to the threshold for genetic instruments, a cautious approach is warranted in interpreting some of the results.

CONCLUSIONS

This study revealed significant genetic causality between focal epilepsy and OCD, as well as between ADHD and generalized epilepsy. However, no casual significance was observed with other psychiatric comorbidities examined. Considering the inherent limitations of MR studies, further research is warranted to definitively clarify these genetic causal associations.

A genome-wide association study based on the China Kadoorie Biobank identifies genetic associations between snoring and cardiometabolic traits

Despite the high prevalence of snoring in Asia, little is known about the genetic etiology of snoring and its causal relationships with cardiometabolic traits. Based on 100,626 Chinese individuals, a genome-wide association study on snoring was conducted. Four novel loci were identified for snoring traits mapped on SLC25A21, the intergenic region of WDR11 and FGFR, NAA25, ALDH2, and VTI1A, respectively. The novel loci highlighted the roles of structural abnormality of the upper airway and craniofacial region and dysfunction of metabolic and transport systems in the development of snoring. In the two-sample bi-directional Mendelian randomization analysis, higher body mass index, weight, and elevated blood pressure were causal for snoring, and a reverse causal effect was observed between snoring and diastolic blood pressure. Altogether, our results revealed the possible etiology of snoring in China and indicated that managing cardiometabolic health was essential to snoring prevention, and hypertension should be considered among snorers.

Development of genome-wide association studies on childhood obesity and its indicators: A scoping review and enrichment analysis

BACKGROUND

The progress of genome-wide association studies (GWAS) in childhood obesity and its indicators is challenging and there are differences in genetic studies in children and adults.

OBJECTIVE

To illustrate the history of the development of GWAS in childhood obesity and its indicators and summarize the GWAS loci.

METHODS

PubMed, Web of Science, Embase and GWAS Catalog databases were systematically searched from 1 January 2005 to 19 October 2022 for literature related to GWAS of childhood BMI, body fatness and obesity. The nearest genes were used as positional genes to perform gene set analyses including the enrichment of pathways, tissues and diseases.

RESULTS

Twenty articles published between 2007 and 2021 were included in this scoping review, which identified 116 SNPs reaching genome-wide significance with childhood BMI (n = 50), body fatness (n = 31) and obesity (n = 35). The study populations were European in 16 studies, non-European in three studies (1 East Asian; 1 American; 1 Mexican) and trans-ancestry in one study. Several enriched pathways, tissues and diseases were identified through enrichment analysis of genes associated with childhood obesity and its indicators.

CONCLUSIONS

The innovations in tools and methods enable GWAS to better explore the genetic characteristics of obesity in children and adolescents. However, the number of GWAS in American, Asian and African populations is limited compared to the European population.

FAVOR: functional annotation of variants online resource and annotator for variation across the human genome

Large biobank-scale whole genome sequencing (WGS) studies are rapidly identifying a multitude of coding and non-coding variants. They provide an unprecedented resource for illuminating the genetic basis of human diseases. Variant functional annotations play a critical role in WGS analysis, result interpretation, and prioritization of disease- or trait-associated causal variants. Existing functional annotation databases have limited scope to perform online queries and functionally annotate the genotype data of large biobank-scale WGS studies. We develop the Functional Annotation of Variants Online Resources (FAVOR) to meet these pressing needs. FAVOR provides a comprehensive multi-faceted variant functional annotation online portal that summarizes and visualizes findings of all possible nine billion single nucleotide variants (SNVs) across the genome. It allows for rapid variant-, gene- and region-level queries of variant functional annotations. FAVOR integrates variant functional information from multiple sources to describe the functional characteristics of variants and facilitates prioritizing plausible causal variants influencing human phenotypes. Furthermore, we provide a scalable annotation tool, FAVORannotator, to functionally annotate large-scale WGS studies and efficiently store the genotype and their variant functional annotation data in a single file using the annotated Genomic Data Structure (aGDS) format, making downstream analysis more convenient. FAVOR and FAVORannotator are available at https://favor.genohub.org.

Evaluation of a whole-exome sequencing pipeline and benchmarking of causal germline variant prioritizers

Most causal variants of Mendelian diseases are exonic. Whole-exome sequencing (WES) has become the diagnostic gold standard, but causative variant prioritization constitutes a bottleneck. Here we assessed an in-house sample-to-sequence pipeline and benchmarked free prioritization tools for germline causal variants from WES data. WES of 61 unselected patients with a known genetic disease cause was obtained. Variant prioritizations were performed by diverse tools and recorded to obtain a diagnostic yield when the causal variant was present in the first, fifth, and 10th top rankings. A fraction of causal variants was not captured by WES (8.2%) or did not pass the quality control criteria (13.1%). Most of the applications inspected were unavailable or had technical limitations, leaving nine tools for complete evaluation. Exomiser performed best in the top first rankings, while LIRICAL led in the top fifth rankings. Based on the more conservative top 10th rankings, Xrare had the highest diagnostic yield, followed by a three-way tie among Exomiser, LIRICAL, and PhenIX, then followed by AMELIE, TAPES, Phen-Gen,  AIVar, and VarNote-PAT. Xrare, Exomiser, LIRICAL, and PhenIX are the most efficient options for variant prioritization in real patient WES data.

From the reference human genome to human pangenome: Premise, promise and challenge

The Reference Human Genome remains the single most important resource for mapping genetic variations and assessing their impact. However, it is monophasic, incomplete and not representative of the variation that exists in the population. Given the extent of ethno-geographic diversity and the consequent diversity in clinical manifestations of these variations, population specific references were developed overtime. The dramatically plummeting cost of sequencing whole genomes and the advent of third generation long range sequencers allowing accurate, error free, telomere-to-telomere assemblies of human genomes present us with a unique and unprecedented opportunity to develop a more composite standard reference consisting of a collection of multiple genomes that capture the maximal variation existing in the population, with the deepest annotation possible, enabling a realistic, reliable and actionable estimation of clinical significance of specific variations. The Human Pangenome Project thus is a logical next step promising a more accurate and global representation of genomic variations. The pangenome effort must be reciprocally complemented with precise variant discovery tools and exhaustive annotation to ensure unambiguous clinical assessment of the variant in ethno-geographical context. Here we discuss a broad roadmap, the challenges and way forward in developing a universal pangenome reference including data visualization techniques and integration of prior knowledge base in the new graph based architecture and tools to submit, compare, query, annotate and retrieve relevant information from the pangenomes. The biggest challenge, however, will be the ethical, legal and social implications and the training of human resource to the new reference paradigm.

QTLbase2: an enhanced catalog of human quantitative trait loci on extensive molecular phenotypes

Deciphering the fine-scale molecular mechanisms that shape the genetic effects at disease-associated loci from genome-wide association studies (GWAS) remains challenging. The key avenue is to identify the essential molecular phenotypes that mediate the causal variant and disease under particular biological conditions. Therefore, integrating GWAS signals with context-specific quantitative trait loci (QTLs) (such as different tissue/cell types, disease states, and perturbations) from extensive molecular phenotypes would present important strategies for full understanding of disease genetics. Via persistent curation and systematic data processing of large-scale human molecular trait QTLs (xQTLs), we updated our previous QTLbase database (now QTLbase2, http://mulinlab.org/qtlbase) to comprehensively analyze and visualize context-specific QTLs across 22 molecular phenotypes and over 95 tissue/cell types. Overall, the resource features the following major updates and novel functions: (i) 960 more genome-wide QTL summary statistics from 146 independent studies; (ii) new data for 10 previously uncompiled QTL types; (iii) variant query scope expanded to fit 195 QTL datasets based on whole-genome sequencing; (iv) supports filtering and comparison of QTLs for different biological conditions, such as stimulation types and disease states; (v) a new linkage disequilibrium viewer to facilitate variant prioritization across tissue/cell types and QTL types.

Echtvar: compressed variant representation for rapid annotation and filtering of SNPs and indels

Germline and somatic variants within an individual or cohort are interpreted with information from large cohorts. Annotation with this information becomes a computational bottleneck as population sets grow to terabytes of data. Here, we introduce echtvar, which efficiently encodes population variants and annotation fields into a compressed archive that can be used for rapid variant annotation and filtering. Most variants, represented by chromosome, position and alleles are encoded into 32-bits-half the size of previous encoding schemes and at least 4 times smaller than a naive encoding. The annotations, stored separately within the same archive, are also encoded and compressed. We show that echtvar is faster and uses less space than existing tools and that it can effectively reduce the number of candidate variants. We give examples on germ-line and somatic variants to document how echtvar can facilitate exploratory data analysis on genetic variants. Echtvar is available at https://github.com/brentp/echtvar under an MIT license.

Scalable approaches for functional analyses of whole-genome sequencing non-coding variants

Non-coding genetic variants outside of protein-coding genome regions play an important role in genetic and epigenetic regulation. It has become increasingly important to understand their roles, as non-coding variants often make up the majority of top findings of genome-wide association studies (GWAS). In addition, the growing popularity of disease-specific whole-genome sequencing (WGS) efforts expands the library of and offers unique opportunities for investigating both common and rare non-coding variants, which are typically not detected in more limited GWAS approaches. However, the sheer size and breadth of WGS data introduce additional challenges to predicting functional impacts in terms of data analysis and interpretation. This review focuses on the recent approaches developed for efficient, at-scale annotation and prioritization of non-coding variants uncovered in WGS analyses. In particular, we review the latest scalable annotation tools, databases and functional genomic resources for interpreting the variant findings from WGS based on both experimental data and in silico predictive annotations. We also review machine learning-based predictive models for variant scoring and prioritization. We conclude with a discussion of future research directions which will enhance the data and tools necessary for the effective functional analyses of variants identified by WGS to improve our understanding of disease etiology.

Interrogating cell type-specific cooperation of transcriptional regulators in 3D chromatin

Context-specific activities of transcription regulators (TRs) in the nucleus modulate spatiotemporal gene expression precisely. Using the largest ChIP-seq data and chromatin loops in the human K562 cell line, we initially interrogated TR cooperation in 3D chromatin via a graphical model and revealed many known and novel TRs manipulating context-specific pathways. To explore TR cooperation across broad tissue/cell types, we systematically leveraged large-scale open chromatin profiles, computational footprinting, and high-resolution chromatin interactions to investigate tissue/cell type-specific TR cooperation. We first delineated a landscape of TR cooperation across 40 human tissue/cell types. Network modularity analyses uncovered the commonality and specificity of TR cooperation in different conditions. We also demonstrated that TR cooperation information can better interpret the disease-causal variants identified by genome-wide association studies and recapitulate cell states during neural development. Our study characterizes shared and unique patterns of TR cooperation associated with the cell type specificity of gene regulation in 3D chromatin.

•

Computational inference of transcriptional regulator (TR) cooperation in 3D chromatin

•

A landscape of 3D TR cooperation across 40 human tissue/cell types

•

TR cooperation can better interpret the disease-causal variants identified by GWAS

•

Cooperation of certain TRs shapes context-specific gene regulation in cell development

VannoPortal: multiscale functional annotation of human genetic variants for interrogating molecular mechanism of traits and diseases

Interpreting the molecular mechanism of genomic variations and their causal relationship with diseases/traits are important and challenging problems in the human genetic study. To provide comprehensive and context-specific variant annotations for biologists and clinicians, here, by systematically integrating over 4TB genomic/epigenomic profiles and frequently-used annotation databases from various biological domains, we develop a variant annotation database, called VannoPortal. In general, the database has following major features: (i) systematically integrates 40 genome-wide variant annotations and prediction scores regarding allele frequency, linkage disequilibrium, evolutionary signature, disease/trait association, tissue/cell type-specific epigenome, base-wise functional prediction, allelic imbalance and pathogenicity; (ii) equips with our recent novel index system and parallel random-sweep searching algorithms for efficient management of backend databases and information extraction; (iii) greatly expands context-dependent variant annotation to incorporate large-scale epigenomic maps and regulatory profiles (such as EpiMap) across over 33 tissue/cell types; (iv) compiles many genome-scale base-wise prediction scores for regulatory/pathogenic variant classification beyond protein-coding region; (v) enables fast retrieval and direct comparison of functional evidence among linked variants using highly interactive web panel in addition to plain table; (vi) introduces many visualization functions for more efficient identification and interpretation of functional variants in single web page. VannoPortal is freely available at http://mulinlab.org/vportal.