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

Rare variant contribution to the heritability of coronary artery disease

Whole genome sequences (WGS) enable discovery of rare variants which may contribute to missing heritability of coronary artery disease (CAD). To measure their contribution, we apply the GREML-LDMS-I approach to WGS of 4949 cases and 17,494 controls of European ancestry from the NHLBI TOPMed program. We estimate CAD heritability at 34.3% assuming a prevalence of 8.2%. Ultra-rare (minor allele frequency ≤ 0.1%) variants with low linkage disequilibrium (LD) score contribute ~50% of the heritability. We also investigate CAD heritability enrichment using a diverse set of functional annotations: i) constraint; ii) predicted protein-altering impact; iii) cis-regulatory elements from a cell-specific chromatin atlas of the human coronary; and iv) annotation principal components representing a wide range of functional processes. We observe marked enrichment of CAD heritability for most functional annotations. These results reveal the predominant role of ultra-rare variants in low LD on the heritability of CAD. Moreover, they highlight several functional processes including cell type-specific regulatory mechanisms as key drivers of CAD genetic risk.

GWAS of multiple neuropathology endophenotypes identifies new risk loci and provides insights into the genetic risk of dementia

Genome-wide association studies (GWAS) have identified >80 Alzheimer's disease and related dementias (ADRD)-associated genetic loci. However, the clinical outcomes used in most previous studies belie the complex nature of underlying neuropathologies. Here we performed GWAS on 11 ADRD-related neuropathology endophenotypes with participants drawn from the following three sources: the National Alzheimer's Coordinating Center, the Religious Orders Study and Rush Memory and Aging Project, and the Adult Changes in Thought study (n = 7,804 total autopsied participants). We identified eight independent significantly associated loci, of which four were new (COL4A1, PIK3R5, LZTS1 and APOC2). Separately testing known ADRD loci, 19 loci were significantly associated with at least one neuropathology after false-discovery rate adjustment. Genetic colocalization analyses identified pleiotropic effects and quantitative trait loci. Methylation in the cerebral cortex at two sites near APOC2 was associated with cerebral amyloid angiopathy. Studies that include neuropathology endophenotypes are an important step in understanding the mechanisms underlying genetic ADRD risk.

FAVOR-GPT: a generative natural language interface to whole genome variant functional annotations

Motivation

Functional Annotation of genomic Variants Online Resources (FAVOR) offers multi-faceted, whole genome variant functional annotations, which is essential for Whole Genome and Exome Sequencing (WGS/WES) analysis and the functional prioritization of disease-associated variants. A versatile chatbot designed to facilitate informative interpretation and interactive, user-centric summary of the whole genome variant functional annotation data in the FAVOR database is needed.

Results

We have developed FAVOR-GPT, a generative natural language interface powered by integrating large language models (LLMs) and FAVOR. It is developed based on the Retrieval Augmented Generation (RAG) approach, and complements the original FAVOR portal, enhancing usability for users, especially those without specialized expertise. FAVOR-GPT simplifies raw annotations by providing interpretable explanations and result summaries in response to the user's prompt. It shows high accuracy when cross-referencing with the FAVOR database, underscoring the robustness of the retrieval framework.

Availability and implementation

Researchers can access FAVOR-GPT at FAVOR's main website (https://favor.genohub.org).

Cell-type-specific mapping of enhancers and target genes from single-cell multimodal data

Mapping enhancers and target genes in disease-related cell types has provided critical insights into the functional mechanisms of genetic variants identified by genomewide association studies (GWAS). However, most existing analyses rely on bulk data or cultured cell lines, which may fail to identify cell-type-specific enhancers and target genes. Recently, single-cell multimodal data measuring both gene expression and chromatin accessibility within the same cells have enabled the inference of enhancer-gene pairs in a cell-type-specific and context-specific manner. However, this task is challenged by the data's high sparsity, sequencing depth variation, and the computational burden of analyzing a large number of enhancer-gene pairs. To address these challenges, we propose scMultiMap, a statistical method that infers enhancer-gene association from sparse multimodal counts using a joint latent-variable model. It adjusts for technical confounding, permits fast moment-based estimation and provides analytically derived p -values. In systematic analyses of blood and brain data, scMultiMap shows appropriate type I error control, high statistical power with greater reproducibility across independent datasets and stronger consistency with orthogonal data modalities. Meanwhile, its computational cost is less than 1% of existing methods. When applied to single-cell multimodal data from postmortem brain samples from Alzheimer's disease (AD) patients and controls, scMultiMap gave the highest heritability enrichment in microglia and revealed new insights into the regulatory mechanisms of AD GWAS variants in microglia.

A methodology for gene level omics-WAS integration identifies genes influencing traits associated with cardiovascular risks: the Long Life Family Study

The Long Life Family Study (LLFS) enrolled 4953 participants in 539 pedigrees displaying exceptional longevity. To identify genetic mechanisms that affect cardiovascular risks in the LLFS population, we developed a multi-omics integration pipeline and applied it to 11 traits associated with cardiovascular risks. Using our pipeline, we aggregated gene-level statistics from rare-variant analysis, GWAS, and gene expression-trait association by Correlated Meta-Analysis (CMA). Across all traits, CMA identified 64 significant genes after Bonferroni correction (p ≤ 2.8 × 10-7), 29 of which replicated in the Framingham Heart Study (FHS) cohort. Notably, 20 of the 29 replicated genes do not have a previously known trait-associated variant in the GWAS Catalog within 50 kb. Thirteen modules in Protein-Protein Interaction (PPI) networks are significantly enriched in genes with low meta-analysis p-values for at least one trait, three of which are replicated in the FHS cohort. The functional annotation of genes in these modules showed a significant over-representation of trait-related biological processes including sterol transport, protein-lipid complex remodeling, and immune response regulation. Among major findings, our results suggest a role of triglyceride-associated and mast-cell functional genes FCER1A, MS4A2, GATA2, HDC, and HRH4 in atherosclerosis risks. Our findings also suggest that lower expression of ATG2A, a gene we found to be associated with BMI, may be both a cause and consequence of obesity. Finally, our results suggest that ENPP3 may play an intermediary role in triglyceride-induced inflammation. Our pipeline is freely available and implemented in the Nextflow workflow language, making it easily runnable on any compute platform ( https://nf-co.re/omicsgenetraitassociation ).

Identification of 16 novel Alzheimer's disease susceptibility loci using multi-ancestry meta-analyses of clinical Alzheimer's disease and AD-by-proxy cases from four whole genome sequencing datasets

Alzheimer's disease (AD) is the most prevalent form of dementia. While many AD-associated genetic determinants have been previously identified, few studies have analyzed individuals of non-European ancestry. Here, we describe a multi-ancestry genome-wide association study of clinically-diagnosed AD and AD-by-proxy using whole genome sequencing data from NIAGADS, NIMH, UKB, and All of Us (AoU) consisting of 49,149 cases (12,074 clinically-diagnosed and 37,075 AD-by-proxy) and 383,225 controls. Nearly half of NIAGADS and AoU participants are of non-European ancestry. For clinically-diagnosed AD, we identified 14 new loci - five common (FBN2,/SCL27A6, AC090115.1, DYM, KCNG1/AL121785.1, TIAM1) and nine rare (VWA5B1, RNU6-755P/LMX1A, MOB1A, MORC1-AS1, LINC00989, PDE4D, RNU2-49P/CDO1, NEO1, and SLC35G3/AC022916.1). Meta-analysis of UKB and AoU AD-by-proxy cases yielded two new rare loci (RPL23/LASP1 and CEBPA/ AC008738.6) which were also nominally significant in NIAGADS. In summary, we provide evidence for 16 novel AD loci and advocate for more studies using WGS-based GWAS of diverse cohorts.

VCAT: an integrated variant function annotation tools

The development of sequencing technology has promoted discovery of variants in the human genome. Identifying functions of these variants is important for us to link genotype to phenotype, and to diagnose diseases. However, it usually requires researchers to visit multiple databases. Here, we presented a one-stop webserver for variant function annotation tools (VCAT, https://biomed.nscc-gz.cn/zhaolab/VCAT/ ) that is the first one connecting variant to functions via the epigenome, protein, drug and RNA. VCAT is also the first one to make all annotations visualized in interactive charts or molecular structures. VCAT allows users to upload data in VCF format, and download results via a URL. Moreover, VCAT has annotated a huge number (1,262,041,068) of variants collected from dbSNP, 1000 Genomes projects, gnomAD, ICGC, TCGA, and HPRC Pangenome project. For these variants, users are able to searcher their functions, related diseases and drugs from VCAT. In summary, VCAT provides a one-stop webserver to explore the potential functions of human genomic variants including their relationship with diseases and drugs.

Whole genome sequencing based analysis of inflammation biomarkers in the Trans-Omics for Precision Medicine (TOPMed) consortium

Inflammation biomarkers can provide valuable insight into the role of inflammatory processes in many diseases and conditions. Sequencing based analyses of such biomarkers can also serve as an exemplar of the genetic architecture of quantitative traits. To evaluate the biological insight, which can be provided by a multi-ancestry, whole-genome based association study, we performed a comprehensive analysis of 21 inflammation biomarkers from up to 38 465 individuals with whole-genome sequencing from the Trans-Omics for Precision Medicine (TOPMed) program (with varying sample size by trait, where the minimum sample size was n = 737 for MMP-1). We identified 22 distinct single-variant associations across 6 traits-E-selectin, intercellular adhesion molecule 1, interleukin-6, lipoprotein-associated phospholipase A2 activity and mass, and P-selectin-that remained significant after conditioning on previously identified associations for these inflammatory biomarkers. We further expanded upon known biomarker associations by pairing the single-variant analysis with a rare variant set-based analysis that further identified 19 significant rare variant set-based associations with 5 traits. These signals were distinct from both significant single variant association signals within TOPMed and genetic signals observed in prior studies, demonstrating the complementary value of performing both single and rare variant analyses when analyzing quantitative traits. We also confirm several previously reported signals from semi-quantitative proteomics platforms. Many of these signals demonstrate the extensive allelic heterogeneity and ancestry-differentiated variant-trait associations common for inflammation biomarkers, a characteristic we hypothesize will be increasingly observed with well-powered, large-scale analyses of complex traits.

A Large-Scale Genome-Wide Study of Gene-Sleep Duration Interactions for Blood Pressure in 811,405 Individuals from Diverse Populations

Although both short and long sleep duration are associated with elevated hypertension risk, our understanding of their interplay with biological pathways governing blood pressure remains limited. To address this, we carried out genome-wide cross-population gene-by-short-sleep and long-sleep duration interaction analyses for three blood pressure traits (systolic, diastolic, and pulse pressure) in 811,405 individuals from diverse population groups. We discover 22 novel gene-sleep duration interaction loci for blood pressure, mapped to 23 genes. Investigating these genes' functional implications shed light on neurological, thyroidal, bone metabolism, and hematopoietic pathways that necessitate future investigation for blood pressure management that caters to sleep health lifestyle. Non-overlap between short sleep (12) and long sleep (10) interactions underscores the plausible nature of distinct influences of both sleep duration extremes in cardiovascular health. Several of our loci are specific towards a particular population background or sex, emphasizing the importance of addressing heterogeneity entangled in gene-environment interactions, when considering precision medicine design approaches for blood pressure management.

Genome-wide association study identifies high-impact susceptibility loci for HCC in North America

BACKGROUND AND AIMS

Despite the substantial impact of environmental factors, individuals with a family history of liver cancer have an increased risk for HCC. However, genetic factors have not been studied systematically by genome-wide approaches in large numbers of individuals from European descent populations (EDP).

APPROACH AND RESULTS

We conducted a 2-stage genome-wide association study (GWAS) on HCC not affected by HBV infections. A total of 1872 HCC cases and 2907 controls were included in the discovery stage, and 1200 HCC cases and 1832 controls in the validation. We analyzed the discovery and validation samples separately and then conducted a meta-analysis. All analyses were conducted in the presence and absence of HCV. The liability-scale heritability was 24.4% for overall HCC. Five regions with significant ORs (95% CI) were identified for nonviral HCC: 3p22.1, MOBP , rs9842969, (0.51, [0.40-0.65]); 5p15.33, TERT , rs2242652, (0.70, (0.62-0.79]); 19q13.11, TM6SF2 , rs58542926, (1.49, [1.29-1.72]); 19p13.11 MAU2 , rs58489806, (1.53, (1.33-1.75]); and 22q13.31, PNPLA3 , rs738409, (1.66, [1.51-1.83]). One region was identified for HCV-induced HCC: 6p21.31, human leukocyte antigen DQ beta 1, rs9275224, (0.79, [0.74-0.84]). A combination of homozygous variants of PNPLA3 and TERT showing a 6.5-fold higher risk for nonviral-related HCC compared to individuals lacking these genotypes. This observation suggests that gene-gene interactions may identify individuals at elevated risk for developing HCC.

CONCLUSIONS

Our GWAS highlights novel genetic susceptibility of nonviral HCC among European descent populations from North America with substantial heritability. Selected genetic influences were observed for HCV-positive HCC. Our findings indicate the importance of genetic susceptibility to HCC development.

FastBiCmrMLM: a fast and powerful compressed variance component mixed logistic model for big genomic case-control genome-wide association study

Large sample datasets have been regarded as the primary basis for innovative discoveries and the solution to missing heritability in genome-wide association studies. However, their computational complexity cannot consider all comprehensive effects and all polygenic backgrounds, which reduces the effectiveness of large datasets. To address these challenges, we included all effects and polygenic backgrounds in a mixed logistic model for binary traits and compressed four variance components into two. The compressed model combined three computational algorithms to develop an innovative method, called FastBiCmrMLM, for large data analysis. These algorithms were tailored to sample size, computational speed, and reduced memory requirements. To mine additional genes, linkage disequilibrium markers were replaced by bin-based haplotypes, which are analyzed by FastBiCmrMLM, named FastBiCmrMLM-Hap. Simulation studies highlighted the superiority of FastBiCmrMLM over GMMAT, SAIGE and fastGWA-GLMM in identifying dominant, small α (allele substitution effect), and rare variants. In the UK Biobank-scale dataset, we demonstrated that FastBiCmrMLM could detect variants as small as 0.03% and with α ≈ 0. In re-analyses of seven diseases in the WTCCC datasets, 29 candidate genes, with both functional and TWAS evidence, around 36 variants identified only by the new methods, strongly validated the new methods. These methods offer a new way to decipher the genetic architecture of binary traits and address the challenges outlined above.

Variant effect predictors: a systematic review and practical guide

Large-scale association analyses using whole-genome sequence data have become feasible, but understanding the functional impacts of these associations remains challenging. Although many tools are available to predict the functional impacts of genetic variants, it is unclear which tool should be used in practice. This work provides a practical guide to assist in selecting appropriate tools for variant annotation. We conducted a MEDLINE search up to November 10, 2023, and included tools that are applicable to a broad range of phenotypes, can be used locally, and have been recently updated. Tools were categorized based on the types of variants they accept and the functional impacts they predict. Sequence Ontology terms were used for standardization. We identified 118 databases and software packages, encompassing 36 variant types and 161 functional impacts. Combining only three tools, namely SnpEff, FAVOR, and SparkINFERNO, allows predicting 99 (61%) distinct functional impacts. Thirty-seven tools predict 89 functional impacts that are not supported by any other tool, while 75 tools predict pathogenicity and can be used within the ACMG/AMP guidelines in a clinical context. We launched a website allowing researchers to select tools based on desired variants and impacts. In summary, more than 100 tools are already available to predict approximately 160 functional impacts. About 60% of the functional impacts can be predicted by the combination of three tools. Unexpectedly, recent tools do not predict more impacts than older ones. Future research should allow predicting the functionality of so far unsupported variant types, such as gene fusions.URL: https://cardio-care.shinyapps.io/VEP_Finder/ .Registration: OSF Registries on November 10, 2023, https://osf.io/s2gct .

The potential impact of polymorphisms in METTL3 gene on knee osteoarthritis susceptibility

Objective

This study was aimed to explore the correlation between METTL3 polymorphisms and susceptibility to knee osteoarthritis (KOA).

Methods

The relationship of five single nucleotide polymorphisms (SNPs) in the METTL3 gene with the susceptibility of KOA was analyzed through multinomial logistic regression analysis in this a case-control study. Genotyping was performed on 228 KOA patients and 252 unaffected individuals from South China based on the TaqMan method. The MDR software (version 3.0.2) was utilized for the analysis of SNP interactions.

Results

Out of the five SNPs examined, the T > G change in the METTL3 gene at the rs1061026 locus increased the risk of KOA, while rs1139130 A > G and rs1263802 C > T variants were found to be linked with a reduced risk of developing KOA with statistical significance. The rs1061027 A > C and rs1263801 C > G variants did not show significant association (p＞0.05). The rs1061026 TG/GG genotype showed a significant correlation with an increased risk of KOA in the following subgroups: the males, individuals with a BMI ranging from 24 to 28, smokers, those who were not engaged in physical exercise (PE), patients who had experienced KOA symptoms for eight years or longer, and those without a family history of the disease or reported swelling. On the other hand, the rs1139130 AG/GG genotype demonstrated a protective effect against KOA among the females, individuals with a BMI greater than or equal to 24, a unilateral KOA, or a KOA duration of 8 years or less, non-smokers, non-alcohol drinkers, those who were not engaged in PE, and those who had no injury or family history, or no experience of knee swelling. Additionally, it was observed that the rs1263802 CT/TT genotypes showed a protective effect among patients without a history of injury. Furthermore, individuals with the haplotypes GAT, GGC, TAT, and TGC were found to have a significantly lower susceptibility to KOA compared to the reference haplotype TAC.

Conclusions

The METTL3 gene variant rs1061026 could increase the risk of KOA, whereas the variants of rs1139130 as well as rs1263802 might exert a protective effect against KOA. These variants could potentially function as susceptibility markers for KOA among the population from South China.

GPN-MSA: an alignment-based DNA language model for genome-wide variant effect prediction

Whereas protein language models have demonstrated remarkable efficacy in predicting the effects of missense variants, DNA counterparts have not yet achieved a similar competitive edge for genome-wide variant effect predictions, especially in complex genomes such as that of humans. To address this challenge, we here introduce GPN-MSA, a novel framework for DNA language models that leverages whole-genome sequence alignments across multiple species and takes only a few hours to train. Across several benchmarks on clinical databases (ClinVar, COSMIC, OMIM), experimental functional assays (DMS, DepMap), and population genomic data (gnomAD), our model for the human genome achieves outstanding performance on deleteriousness prediction for both coding and non-coding variants.

Variants located in intron 6 of SMN1 lead to misdiagnosis in genetic detection and screening for SMA

Accurate genetic diagnosis is necessary for guiding the treatment of spinal muscular atrophy (SMA). An updated consensus for the diagnosis and management of SMA was published in 2018. However, clinicians should remain alert to some pitfalls of genetic testing that can occur when following a routine diagnosis. In this study, we report the diagnosis of three unrelated individuals who were initially misdiagnosed as carrying a homozygous deletion of SMN1 exon 7. MLPA (P060 and P021) and qPCR were used to detect the copy number of SMN. SMN1 variants were identified by SMN1 clone and next-generation sequencing (NGS). Transcription of SMN1 variants was detected using qRT-PCR and ex vivo splicing analysis. Among the three individuals, one was identified as a patient with SMA carrying a heterozygous deletion and a pathogenic variant (c.835-17_835-14delCTTT) of SMN1, one was a healthy carrier only carrying a heterozygous deletion of SMN1 exon 7, and the third was a patient with nemaline myopathy 2 carrying a heterozygous deletion of SMN1 exon 7. The misdiagnosis of these individuals was attributed to the presence of the c.835-17_835-14delCTTT or c.835-17C > G variants in SMN1 intron 6, which affect the amplification of SMN1 exon 7 during MLPA-P060 and qPCR testing. However, MLPA-P021 and NGS analyses were unaffected by these variants. These results support that additional detection methods should be employed in cases where the SMN1 copy number is ambiguous to minimize the misdiagnosis of SMA.

Mitochondrial and Nuclear DNA Variants in Amyotrophic Lateral Sclerosis: Enrichment in the Mitochondrial Control Region and Sirtuin Pathway Genes in Spinal Cord Tissue

Amyotrophic Lateral Sclerosis (ALS) is a progressive disease with prevalent mitochondrial dysfunctions affecting both upper and lower motor neurons in the motor cortex, brainstem, and spinal cord. Despite mitochondria having their own genome (mtDNA), in humans, most mitochondrial genes are encoded by the nuclear genome (nDNA). Our study aimed to simultaneously screen for nDNA and mtDNA genomes to assess for specific variant enrichment in ALS compared to control tissues. Here, we analysed whole exome (WES) and whole genome (WGS) sequencing data from spinal cord tissues, respectively, of 6 and 12 human donors. A total of 31,257 and 301,241 variants in nuclear-encoded mitochondrial genes were identified from WES and WGS, respectively, while mtDNA reads accounted for 73 and 332 variants. Despite technical differences, both datasets consistently revealed a specific enrichment of variants in the mitochondrial Control Region (CR) and in several of these genes directly associated with mitochondrial dynamics or with Sirtuin pathway genes within ALS tissues. Overall, our data support the hypothesis of a variant burden in specific genes, highlighting potential actionable targets for therapeutic interventions in ALS.

SUMMIT-FA: a new resource for improved transcriptome imputation using functional annotations

Transcriptome-wide association studies (TWAS) integrate gene expression prediction models and genome-wide association studies (GWAS) to identify gene-trait associations. The power of TWAS is determined by the sample size of GWAS and the accuracy of the expression prediction model. Here, we present a new method, the Summary-level Unified Method for Modeling Integrated Transcriptome using Functional Annotations (SUMMIT-FA), which improves gene expression prediction accuracy by leveraging functional annotation resources and a large expression quantitative trait loci (eQTL) summary-level dataset. We build gene expression prediction models in whole blood using SUMMIT-FA with the comprehensive functional database MACIE and eQTL summary-level data from the eQTLGen consortium. We apply these models to GWAS for 24 complex traits and show that SUMMIT-FA identifies significantly more gene-trait associations and improves predictive power for identifying "silver standard" genes compared to several benchmark methods. We further conduct a simulation study to demonstrate the effectiveness of SUMMIT-FA.

A Large-Scale Genome-Wide Study of Gene-Sleep Duration Interactions for Blood Pressure in 811,405 Individuals from Diverse Populations

Although both short and long sleep duration are associated with elevated hypertension risk, our understanding of their interplay with biological pathways governing blood pressure remains limited. To address this, we carried out genome-wide cross-population gene-by-short-sleep and long-sleep duration interaction analyses for three blood pressure traits (systolic, diastolic, and pulse pressure) in 811,405 individuals from diverse population groups. We discover 22 novel gene-sleep duration interaction loci for blood pressure, mapped to genes involved in neurological, thyroidal, bone metabolism, and hematopoietic pathways. Non-overlap between short sleep (12) and long sleep (10) interactions underscores the plausibility of distinct influences of both sleep duration extremes in cardiovascular health. With several of our loci reflecting specificity towards population background or sex, our discovery sheds light on the importance of embracing granularity when addressing heterogeneity entangled in gene-environment interactions, and in therapeutic design approaches for blood pressure management.

Computational Tools to Assist in Analyzing Effects of the SERPINA1 Gene Variation on Alpha-1 Antitrypsin (AAT)

In the rapidly advancing field of bioinformatics, the development and application of computational tools to predict the effects of single nucleotide variants (SNVs) are shedding light on the molecular mechanisms underlying disorders. Also, they hold promise for guiding therapeutic interventions and personalized medicine strategies in the future. A comprehensive understanding of the impact of SNVs in the SERPINA1 gene on alpha-1 antitrypsin (AAT) protein structure and function requires integrating bioinformatic approaches. Here, we provide a guide for clinicians to navigate through the field of computational analyses which can be applied to describe a novel genetic variant. Predicting the clinical significance of SERPINA1 variation allows clinicians to tailor treatment options for individuals with alpha-1 antitrypsin deficiency (AATD) and related conditions, ultimately improving the patient's outcome and quality of life. This paper explores the various bioinformatic methodologies and cutting-edge approaches dedicated to the assessment of molecular variants of genes and their product proteins using SERPINA1 and AAT as an example.

Multi-omics Integration Identifies Genes Influencing Traits Associated with Cardiovascular Risks: The Long Life Family Study

The Long Life Family Study (LLFS) enrolled 4,953 participants in 539 pedigrees displaying exceptional longevity. To identify genetic mechanisms that affect cardiovascular risks in the LLFS population, we developed a multi-omics integration pipeline and applied it to 11 traits associated with cardiovascular risks. Using our pipeline, we aggregated gene-level statistics from rare-variant analysis, GWAS, and gene expression-trait association by Correlated Meta-Analysis (CMA). Across all traits, CMA identified 64 significant genes after Bonferroni correction (p ≤ 2.8×10-7), 29 of which replicated in the Framingham Heart Study (FHS) cohort. Notably, 20 of the 29 replicated genes do not have a previously known trait-associated variant in the GWAS Catalog within 50 kb. Thirteen modules in Protein-Protein Interaction (PPI) networks are significantly enriched in genes with low meta-analysis p-values for at least one trait, three of which are replicated in the FHS cohort. The functional annotation of genes in these modules showed a significant over-representation of trait-related biological processes including sterol transport, protein-lipid complex remodeling, and immune response regulation. Among major findings, our results suggest a role of triglyceride-associated and mast-cell functional genes FCER1A, MS4A2, GATA2, HDC, and HRH4 in atherosclerosis risks. Our findings also suggest that lower expression of ATG2A, a gene we found to be associated with BMI, may be both a cause and consequence of obesity. Finally, our results suggest that ENPP3 may play an intermediary role in triglyceride-induced inflammation. Our pipeline is freely available and implemented in the Nextflow workflow language, making it easily runnable on any compute platform (https://nf-co.re/omicsgenetraitassociation).

Genetic influences on circulating retinol and its relationship to human health

Retinol is a fat-soluble vitamin that plays an essential role in many biological processes throughout the human lifespan. Here, we perform the largest genome-wide association study (GWAS) of retinol to date in up to 22,274 participants. We identify eight common variant loci associated with retinol, as well as a rare-variant signal. An integrative gene prioritisation pipeline supports novel retinol-associated genes outside of the main retinol transport complex (RBP4:TTR) related to lipid biology, energy homoeostasis, and endocrine signalling. Genetic proxies of circulating retinol were then used to estimate causal relationships with almost 20,000 clinical phenotypes via a phenome-wide Mendelian randomisation study (MR-pheWAS). The MR-pheWAS suggests that retinol may exert causal effects on inflammation, adiposity, ocular measures, the microbiome, and MRI-derived brain phenotypes, amongst several others. Conversely, circulating retinol may be causally influenced by factors including lipids and serum creatinine. Finally, we demonstrate how a retinol polygenic score could identify individuals more likely to fall outside of the normative range of circulating retinol for a given age. In summary, this study provides a comprehensive evaluation of the genetics of circulating retinol, as well as revealing traits which should be prioritised for further investigation with respect to retinol related therapies or nutritional intervention.

Carrier frequencies, trends, and geographical distribution of hearing loss variants in China: The pooled analysis of 2,161,984 newborns

The aim of this study is to comprehensively investigate the prevalence and distribution patterns of three common genetic variants associated with hearing loss (HL) in Chinese neonatal population. Methods: Prior to June 30, 2023, an extensive search and screening process was conducted across multiple literature databases. R software was utilized for conducting meta-analyses, cartography, and correlation analyses. Results: Firstly, our study identified a total of 99 studies meeting the inclusion criteria. Notably, provinces such as Qinghai, Tibet, Jilin, and Heilongjiang lack large-scale genetic screening data for neonatal deafness. Secondly, in Chinese newborns, the carrier frequencies of GJB2 variants (c.235delC, c.299_300delAT) were 1.63 % (95 %CI 1.52 %-1.76 %) and 0.33 % (95 %CI 0.30 %-0.37 %); While SLC26A4 variants (c.919-2A > G, c.2168A > G) exhibited carrier rates of 0.95 % (95 %CI 0.86 %-1.04 %) and 0.17 % (95 %CI 0.15 %-0.19 %); Additionally, Mt 12S rRNA m.1555 A > G variant was found at a rate of 0.24 % (95 % CI 0.22 %-0.26 %). Thirdly, the mutation rate of GJB2 c.235delC was higher in the east of the Heihe-Tengchong line, whereas the mutation rate of Mt 12S rRNA m.1555 A > G variant exhibited the opposite pattern. Forthly, no significant correlation exhibited the opposite pattern of GJB2 variants, but there was a notable correlation among SLC26A4 variants. Lastly, strong regional distribution correlations were evident between mutation sites from different genes, particularly between SLC26A4 (c.919-2A > G and c.2168A > G) and GJB c.299_300delAT. Conclusions: The most prevalent deafness genes among Chinese neonates were GJB2 c.235delC variant, followed by SLC26A4 c.919-2A > G variant. These gene mutation rates exhibit significant regional distribution characteristics. Consequently, it is imperative to enhance genetic screening efforts to reduce the incidence of deafness in high-risk areas.

NCAD v1.0: a database for non-coding variant annotation and interpretation

The application of whole genome sequencing is expanding in clinical diagnostics across various genetic disorders, and the significance of non-coding variants in penetrant diseases is increasingly being demonstrated. Therefore, it is urgent to improve the diagnostic yield by exploring the pathogenic mechanisms of variants in non-coding regions. However, the interpretation of non-coding variants remains a significant challenge, due to the complex functional regulatory mechanisms of non-coding regions and the current limitations of available databases and tools. Hence, we develop the non-coding variant annotation database (NCAD, http://www.ncawdb.net/), encompassing comprehensive insights into 665,679,194 variants, regulatory elements, and element interaction details. Integrating data from 96 sources, spanning both GRCh37 and GRCh38 versions, NCAD v1.0 provides vital information to support the genetic diagnosis of non-coding variants, including allele frequencies of 12 diverse populations, with a particular focus on the population frequency information for 230,235,698 variants in 20,964 Chinese individuals. Moreover, it offers prediction scores for variant functionality, five categories of regulatory elements, and four types of non-coding RNAs. With its rich data and comprehensive coverage, NCAD serves as a valuable platform, empowering researchers and clinicians with profound insights into non-coding regulatory mechanisms while facilitating the interpretation of non-coding variants.

The causal effect of HbA1c on white matter brain aging by two-sample Mendelian randomization analysis

Background

Poor glycemic control with elevated levels of hemoglobin A1c (HbA1c) is associated with increased risk of cognitive impairment, with potentially varying effects between sexes. However, the causal impact of poor glycemic control on white matter brain aging in men and women is uncertain.

Methods

We used two nonoverlapping data sets from UK Biobank cohort: gene-outcome group (with neuroimaging data, (N = 15,193; males/females: 7,101/8,092)) and gene-exposure group (without neuroimaging data, (N = 279,011; males/females: 122,638/156,373)). HbA1c was considered the exposure and adjusted "brain age gap" (BAG) was calculated on fractional anisotropy (FA) obtained from brain imaging as the outcome, thereby representing the difference between predicted and chronological age. The causal effects of HbA1c on adjusted BAG were studied using the generalized inverse variance weighted (gen-IVW) and other sensitivity analysis methods, including Mendelian randomization (MR)-weighted median, MR-pleiotropy residual sum and outlier, MR-using mixture models, and leave-one-out analysis.

Results

We found that for every 6.75 mmol/mol increase in HbA1c, there was an increase of 0.49 (95% CI = 0.24, 0.74; p-value = 1.30 × 10-4) years in adjusted BAG. Subgroup analyses by sex and age revealed significant causal effects of HbA1c on adjusted BAG, specifically among men aged 60-73 (p-value = 2.37 × 10-8).

Conclusion

Poor glycemic control has a significant causal effect on brain aging, and is most pronounced among older men aged 60-73 years, which provides insights between glycemic control and the susceptibility to age-related neurodegenerative diseases.

VarCards2: an integrated genetic and clinical database for ACMG-AMP variant-interpretation guidelines in the human whole genome

VarCards, an online database, combines comprehensive variant- and gene-level annotation data to streamline genetic counselling for coding variants. Recognising the increasing clinical relevance of non-coding variations, there has been an accelerated development of bioinformatics tools dedicated to interpreting non-coding variations, including single-nucleotide variants and copy number variations. Regrettably, most tools remain as either locally installed databases or command-line tools dispersed across diverse online platforms. Such a landscape poses inconveniences and challenges for genetic counsellors seeking to utilise these resources without advanced bioinformatics expertise. Consequently, we developed VarCards2, which incorporates nearly nine billion artificially generated single-nucleotide variants (including those from mitochondrial DNA) and compiles vital annotation information for genetic counselling based on ACMG-AMP variant-interpretation guidelines. These annotations include (I) functional effects; (II) minor allele frequencies; (III) comprehensive function and pathogenicity predictions covering all potential variants, such as non-synonymous substitutions, non-canonical splicing variants, and non-coding variations and (IV) gene-level information. Furthermore, VarCards2 incorporates 368 820 266 documented short insertions and deletions and 2 773 555 documented copy number variations, complemented by their corresponding annotation and prediction tools. In conclusion, VarCards2, by integrating over 150 variant- and gene-level annotation sources, significantly enhances the efficiency of genetic counselling and can be freely accessed at http://www.genemed.tech/varcards2/.

An in-depth association analysis of genetic variants within nicotine-related loci: Meeting in middle of GWAS and genetic fine-mapping

In the last two decades of Genome-wide association studies (GWAS), nicotine-dependence-related genetic loci (e.g., nicotinic acetylcholine receptor - nAChR subunit genes) are among the most replicable genetic findings. Although GWAS results have reported tens of thousands of SNPs within these loci, further analysis (e.g., fine-mapping) is required to identify the causal variants. However, it is computationally challenging for existing fine-mapping methods to reliably identify causal variants from thousands of candidate SNPs based on the posterior inclusion probability. To address this challenge, we propose a new method to select SNPs by jointly modeling the SNP-wise inference results and the underlying structured network patterns of the linkage disequilibrium (LD) matrix. We use adaptive dense subgraph extraction method to recognize the latent network patterns of the LD matrix and then apply group LASSO to select causal variant candidates. We applied this new method to the UK biobank data to identify the causal variant candidates for nicotine addiction. Eighty-one nicotine addiction-related SNPs (i.e.,-log(p) > 50) of nAChR were selected, which are highly correlated (average r2>0.8) although they are physically distant (e.g., >200 kilobase away) and from various genes. These findings revealed that distant SNPs from different genes can show higher LD r2 than their neighboring SNPs, and jointly contribute to a complex trait like nicotine addiction.

A statistical framework for powerful multi-trait rare variant analysis in large-scale whole-genome sequencing studies

Large-scale whole-genome sequencing (WGS) studies have improved our understanding of the contributions of coding and noncoding rare variants to complex human traits. Leveraging association effect sizes across multiple traits in WGS rare variant association analysis can improve statistical power over single-trait analysis, and also detect pleiotropic genes and regions. Existing multi-trait methods have limited ability to perform rare variant analysis of large-scale WGS data. We propose MultiSTAAR, a statistical framework and computationally-scalable analytical pipeline for functionally-informed multi-trait rare variant analysis in large-scale WGS studies. MultiSTAAR accounts for relatedness, population structure and correlation among phenotypes by jointly analyzing multiple traits, and further empowers rare variant association analysis by incorporating multiple functional annotations. We applied MultiSTAAR to jointly analyze three lipid traits (low-density lipoprotein cholesterol, high-density lipoprotein cholesterol and triglycerides) in 61,861 multi-ethnic samples from the Trans-Omics for Precision Medicine (TOPMed) Program. We discovered new associations with lipid traits missed by single-trait analysis, including rare variants within an enhancer of NIPSNAP3A and an intergenic region on chromosome 1.

hipFG: high-throughput harmonization and integration pipeline for functional genomics data

SUMMARY

Preparing functional genomic (FG) data with diverse assay types and file formats for integration into analysis workflows that interpret genome-wide association and other studies is a significant and time-consuming challenge. Here we introduce hipFG (Harmonization and Integration Pipeline for Functional Genomics), an automatically customized pipeline for efficient and scalable normalization of heterogenous FG data collections into standardized, indexed, rapidly searchable analysis-ready datasets while accounting for FG datatypes (e.g. chromatin interactions, genomic intervals, quantitative trait loci).

AVAILABILITY AND IMPLEMENTATION

hipFG is freely available at https://bitbucket.org/wanglab-upenn/hipFG. A Docker container is available at https://hub.docker.com/r/wanglab/hipfg.

Rare variants in long non-coding RNAs are associated with blood lipid levels in the TOPMed whole-genome sequencing study

Long non-coding RNAs (lncRNAs) are known to perform important regulatory functions in lipid metabolism. Large-scale whole-genome sequencing (WGS) studies and new statistical methods for variant set tests now provide an opportunity to assess more associations between rare variants in lncRNA genes and complex traits across the genome. In this study, we used high-coverage WGS from 66,329 participants of diverse ancestries with measurement of blood lipids and lipoproteins (LDL-C, HDL-C, TC, and TG) in the National Heart, Lung, and Blood Institute (NHLBI) Trans-Omics for Precision Medicine (TOPMed) program to investigate the role of lncRNAs in lipid variability. We aggregated rare variants for 165,375 lncRNA genes based on their genomic locations and conducted rare-variant aggregate association tests using the STAAR (variant-set test for association using annotation information) framework. We performed STAAR conditional analysis adjusting for common variants in known lipid GWAS loci and rare-coding variants in nearby protein-coding genes. Our analyses revealed 83 rare lncRNA variant sets significantly associated with blood lipid levels, all of which were located in known lipid GWAS loci (in a ±500-kb window of a Global Lipids Genetics Consortium index variant). Notably, 61 out of 83 signals (73%) were conditionally independent of common regulatory variation and rare protein-coding variation at the same loci. We replicated 34 out of 61 (56%) conditionally independent associations using the independent UK Biobank WGS data. Our results expand the genetic architecture of blood lipids to rare variants in lncRNAs.

Translating pharmacogenomic sequencing data into drug response predictions-How to interpret variants of unknown significance

The rapid development of sequencing technologies during the past 20 years has provided a variety of methods and tools to interrogate human genomic variations at the population level. Pharmacogenes are well known to be highly polymorphic and a plethora of pharmacogenomic variants has been identified in population sequencing data. However, so far only a small number of these variants have been functionally characterized regarding their impact on drug efficacy and toxicity and the significance of the vast majority remains unknown. It is therefore of high importance to develop tools and frameworks to accurately infer the effects of pharmacogenomic variants and, eventually, aggregate the effect of individual variations into personalized drug response predictions. To address this challenge, we here first describe the technological advances, including sequencing methods and accompanying bioinformatic processing pipelines that have enabled reliable variant identification. Subsequently, we highlight advances in computational algorithms for pharmacogenomic variant interpretation and discuss the added value of emerging strategies, such as machine learning and the integrative use of omics techniques that have the potential to further contribute to the refinement of personalized pharmacological response predictions. Lastly, we provide an overview of experimental and clinical approaches to validate in silico predictions. We conclude that the iterative feedback between computational predictions and experimental validations is likely to rapidly improve the accuracy of pharmacogenomic prediction models, which might soon allow for an incorporation of the entire pharmacogenetic profile into personalized response predictions.

Whole Genome Sequencing Based Analysis of Inflammation Biomarkers in the Trans-Omics for Precision Medicine (TOPMed) Consortium

Inflammation biomarkers can provide valuable insight into the role of inflammatory processes in many diseases and conditions. Sequencing based analyses of such biomarkers can also serve as an exemplar of the genetic architecture of quantitative traits. To evaluate the biological insight, which can be provided by a multi-ancestry, whole-genome based association study, we performed a comprehensive analysis of 21 inflammation biomarkers from up to 38,465 individuals with whole-genome sequencing from the Trans-Omics for Precision Medicine (TOPMed) program. We identified 22 distinct single-variant associations across 6 traits - E-selectin, intercellular adhesion molecule 1, interleukin-6, lipoprotein-associated phospholipase A2 activity and mass, and P-selectin - that remained significant after conditioning on previously identified associations for these inflammatory biomarkers. We further expanded upon known biomarker associations by pairing the single-variant analysis with a rare variant set-based analysis that further identified 19 significant rare variant set-based associations with 5 traits. These signals were distinct from both significant single variant association signals within TOPMed and genetic signals observed in prior studies, demonstrating the complementary value of performing both single and rare variant analyses when analyzing quantitative traits. We also confirm several previously reported signals from semi-quantitative proteomics platforms. Many of these signals demonstrate the extensive allelic heterogeneity and ancestry-differentiated variant-trait associations common for inflammation biomarkers, a characteristic we hypothesize will be increasingly observed with well-powered, large-scale analyses of complex traits.

Rare Variants in Primary Immunodeficiency Genes and Their Functional Partners in Severe COVID-19

The development of severe COVID-19, which is a complex multisystem disease, is thought to be associated with many genes whose action is modulated by numerous environmental and genetic factors. In this study, we focused on the ideas of the omnigenic model of heritability of complex traits, which assumes that a small number of core genes and a large pool of peripheral genes expressed in disease-relevant tissues contribute to the genetics of complex traits through interconnected networks. We hypothesized that primary immunodeficiency disease (PID) genes may be considered as core genes in severe COVID-19, and their functional partners (FPs) from protein-protein interaction networks may be considered as peripheral near-core genes. We used whole-exome sequencing data from patients aged ≤ 45 years with severe (n = 9) and non-severe COVID-19 (n = 11), and assessed the cumulative contribution of rare high-impact variants to disease severity. In patients with severe COVID-19, an excess of rare high-impact variants was observed at the whole-exome level, but maximal association signals were detected for PID + FP gene subsets among the genes intolerant to LoF variants, haploinsufficient and essential. Our exploratory study may serve as a model for new directions in the research of host genetics in severe COVID-19.

Pathogenic variants of human GABRA1 gene associated with epilepsy: A computational approach

Critical for brain development, neurodevelopmental and network disorders, the GABRA1 gene encodes for the α1 subunit, an abundantly and developmentally expressed subunit of heteropentameric gamma-aminobutyric acid A receptors (GABAARs) mediating primary inhibition in the brain. Mutations of the GABAAR subunit genes including GABRA1 gene are associated with epilepsy, a group of syndromes, characterized by unprovoked seizures and diagnosed by integrative approach, that involves genetic testing. Despite the diagnostic use of genetic testing, a large fraction of the GABAAR subunit gene variants including the variants of GABRA1 gene is not known in terms of their molecular consequence, a challenge for precision and personalized medicine. Addressing this, one hundred thirty-seven GABRA1 gene variants of unknown clinical significance have been extracted from the ClinVar database and computationally analyzed for pathogenicity. Eight variants (L49H, P59L, W97R, D99G, G152S, V270G, T294R, P305L) are predicted as pathogenic and mapped to the α1 subunit's extracellular domain (ECD), transmembrane domains (TMDs) and extracellular linker. This is followed by the integration with relevant data for cellular pathology and severity of the epilepsy syndromes retrieved from the literature. Our results suggest that the pathogenic variants in the ECD of GABRA1 (L49H, P59L, W97R, D99G, G152S) will probably manifest decreased surface expression and reduced current with mild epilepsy phenotypes while V270G, T294R in the TMDs and P305L in the linker between the second and the third TMDs will likely cause reduced cell current with severe epilepsy phenotypes. The results presented in this study provides insights for clinical genetics and wet lab experimentation.

The sequence kernel association test for multicategorical outcomes

Disease heterogeneity is ubiquitous in biomedical and clinical studies. In genetic studies, researchers are increasingly interested in understanding the distinct genetic underpinning of subtypes of diseases. However, existing set-based analysis methods for genome-wide association studies are either inadequate or inefficient to handle such multicategorical outcomes. In this paper, we proposed a novel set-based association analysis method, sequence kernel association test (SKAT)-MC, the sequence kernel association test for multicategorical outcomes (nominal or ordinal), which jointly evaluates the relationship between a set of variants (common and rare) and disease subtypes. Through comprehensive simulation studies, we showed that SKAT-MC effectively preserves the nominal type I error rate while substantially increases the statistical power compared to existing methods under various scenarios. We applied SKAT-MC to the Polish breast cancer study (PBCS), and identified gene FGFR2 was significantly associated with estrogen receptor (ER)+ and ER- breast cancer subtypes. We also investigated educational attainment using UK Biobank data (N = 127 , 127 $N=127,127$ ) with SKAT-MC, and identified 21 significant genes in the genome. Consequently, SKAT-MC is a powerful and efficient analysis tool for genetic association studies with multicategorical outcomes. A freely distributed R package SKAT-MC can be accessed at https://github.com/Zhiwen-Owen-Jiang/SKATMC.

Deciphering the causal relationship between blood pressure and regional white matter integrity: A two-sample Mendelian randomization study

Elevated arterial blood pressure (BP) is a common risk factor for cerebrovascular and cardiovascular diseases, but no causal relationship has been established between BP and cerebral white matter (WM) integrity. In this study, we performed a two-sample Mendelian randomization (MR) analysis with individual-level data by defining two nonoverlapping sets of European ancestry individuals (genetics-exposure set: N = 203,111; mean age = 56.71 years, genetics-outcome set: N = 16,156; mean age = 54.61 years) from UK Biobank to evaluate the causal effects of BP on regional WM integrity, measured by fractional anisotropy of diffusion tensor imaging. Two BP traits: systolic and diastolic blood pressure were used as exposures. Genetic variant was carefully selected as instrumental variable (IV) under the MR analysis assumptions. We existing large-scale genome-wide association study summary data for validation. The main method used was a generalized version of inverse-variance weight method while other MR methods were also applied for consistent findings. Two additional MR analyses were performed to exclude the possibility of reverse causality. We found significantly negative causal effects (FDR-adjusted p < .05; every 10 mmHg increase in BP leads to a decrease in FA value by .4% ~ 2%) of BP traits on a union set of 17 WM tracts, including brain regions related to cognitive function and memory. Our study extended the previous findings of association to causation for regional WM integrity, providing insights into the pathological processes of elevated BP that might chronically alter the brain microstructure in different regions.

An allelic-series rare-variant association test for candidate-gene discovery

Allelic series are of candidate therapeutic interest because of the existence of a dose-response relationship between the functionality of a gene and the degree or severity of a phenotype. We define an allelic series as a collection of variants in which increasingly deleterious mutations lead to increasingly large phenotypic effects, and we have developed a gene-based rare-variant association test specifically targeted to identifying genes containing allelic series. Building on the well-known burden test and sequence kernel association test (SKAT), we specify a variety of association models covering different genetic architectures and integrate these into a Coding-Variant Allelic-Series Test (COAST). Through extensive simulations, we confirm that COAST maintains the type I error and improves the power when the pattern of coding-variant effect sizes increases monotonically with mutational severity. We applied COAST to identify allelic-series genes for four circulating-lipid traits and five cell-count traits among 145,735 subjects with available whole-exome sequencing data from the UK Biobank. Compared with optimal SKAT (SKAT-O), COAST identified 29% more Bonferroni-significant associations with circulating-lipid traits, on average, and 82% more with cell-count traits. All of the gene-trait associations identified by COAST have corroborating evidence either from rare-variant associations in the full cohort (Genebass, n = 400,000) or from common-variant associations in the GWAS Catalog. In addition to detecting many gene-trait associations present in Genebass by using only a fraction (36.9%) of the sample, COAST detects associations, such as that between ANGPTL4 and triglycerides, that are absent from Genebass but that have clear common-variant support.

The Impact of Genomic Variation on Function (IGVF) Consortium

Our genomes influence nearly every aspect of human biology from molecular and cellular functions to phenotypes in health and disease. Human genetics studies have now associated hundreds of thousands of differences in our DNA sequence ("genomic variation") with disease risk and other phenotypes, many of which could reveal novel mechanisms of human biology and uncover the basis of genetic predispositions to diseases, thereby guiding the development of new diagnostics and therapeutics. Yet, understanding how genomic variation alters genome function to influence phenotype has proven challenging. To unlock these insights, we need a systematic and comprehensive catalog of genome function and the molecular and cellular effects of genomic variants. Toward this goal, the Impact of Genomic Variation on Function (IGVF) Consortium will combine approaches in single-cell mapping, genomic perturbations, and predictive modeling to investigate the relationships among genomic variation, genome function, and phenotypes. Through systematic comparisons and benchmarking of experimental and computational methods, we aim to create maps across hundreds of cell types and states describing how coding variants alter protein activity, how noncoding variants change the regulation of gene expression, and how both coding and noncoding variants may connect through gene regulatory and protein interaction networks. These experimental data, computational predictions, and accompanying standards and pipelines will be integrated into an open resource that will catalyze community efforts to explore genome function and the impact of genetic variation on human biology and disease across populations.

Identification and validation of risk score model based on gene set activity as a diagnostic biomarker for endometriosis

Objective

The enigmatic nature of Endometriosis (EMS) pathogenesis necessitates investigating alterations in signaling pathway activity to enhance our comprehension of the disease's characteristics.

Methods

Three published gene expression profiles (GSE11691, GSE25628, and GSE7305 datasets) were downloaded, and the "combat" algorithm was employed for batch correction, gene expression difference analysis, and pathway enrichment difference analysis. The protein-protein interaction (PPI) network was constructed to identify core genes, and the relative enrichment degree of gene sets was evaluated. The Lasso regression model identified candidate gene sets with diagnostic value, and a risk scoring diagnostic model was constructed for further validation on the GSE86534 and GSE5108 datasets. CIBERSORT was used to assess the composition of immune cells in EMS, and the correlation between EMS diagnostic value gene sets and immune cells was evaluated.

Results

A total of 568 differentially expressed genes were identified between eutopic and ectopic endometrium, with 10 core genes in the PPI network associated with cell cycle regulation. Inflammation-related pathways, including cytokine-receptor signaling and chemokine signaling pathways, were significantly more active in ectopic endometrium compared to eutopic endometrium. Diagnostic gene sets for EMS, such as homologous recombination, base excision repair, DNA replication, P53 signaling pathway, adherens junction, and SNARE interactions in vesicular transport, were identified. The risk score's area under the curve (AUC) was 0.854, as indicated by the receiver operating characteristic (ROC) curve, and the risk score's diagnostic value was validated by the validation cohort. Immune cell infiltration analysis revealed correlations between the risk score and Macrophages M2, Plasma cells, resting NK cells, activated NK cells, and regulatory T cells.

Conclusion

The risk scoring diagnostic model, based on pathway activity, demonstrates high diagnostic value and offers novel insights and strategies for the clinical diagnosis and treatment of Endometriosis.

Rare variants in long non-coding RNAs are associated with blood lipid levels in the TOPMed Whole Genome Sequencing Study

Long non-coding RNAs (lncRNAs) are known to perform important regulatory functions. Large-scale whole genome sequencing (WGS) studies and new statistical methods for variant set tests now provide an opportunity to assess the associations between rare variants in lncRNA genes and complex traits across the genome. In this study, we used high-coverage WGS from 66,329 participants of diverse ancestries with blood lipid levels (LDL-C, HDL-C, TC, and TG) in the National Heart, Lung, and Blood Institute (NHLBI) Trans-Omics for Precision Medicine (TOPMed) program to investigate the role of lncRNAs in lipid variability. We aggregated rare variants for 165,375 lncRNA genes based on their genomic locations and conducted rare variant aggregate association tests using the STAAR (variant-Set Test for Association using Annotation infoRmation) framework. We performed STAAR conditional analysis adjusting for common variants in known lipid GWAS loci and rare coding variants in nearby protein coding genes. Our analyses revealed 83 rare lncRNA variant sets significantly associated with blood lipid levels, all of which were located in known lipid GWAS loci (in a ±500 kb window of a Global Lipids Genetics Consortium index variant). Notably, 61 out of 83 signals (73%) were conditionally independent of common regulatory variations and rare protein coding variations at the same loci. We replicated 34 out of 61 (56%) conditionally independent associations using the independent UK Biobank WGS data. Our results expand the genetic architecture of blood lipids to rare variants in lncRNA, implicating new therapeutic opportunities.

COVID-19 severity: does the genetic landscape of rare variants matter?

Rare variants affecting host defense against pathogens may be involved in COVID-19 severity, but most rare variants are not expected to have a major impact on the course of COVID-19. We hypothesized that the accumulation of weak effects of many rare functional variants throughout the exome may contribute to the overall risk in patients with severe disease. This assumption is consistent with the omnigenic model of the relationship between genetic and phenotypic variation in complex traits, according to which association signals tend to spread across most of the genome through gene regulatory networks from genes outside the major pathways to disease-related genes. We performed whole-exome sequencing and compared the burden of rare variants in 57 patients with severe and 29 patients with mild/moderate COVID-19. At the whole-exome level, we observed an excess of rare, predominantly high-impact (HI) variants in the group with severe COVID-19. Restriction to genes intolerant to HI or damaging missense variants increased enrichment for these classes of variants. Among various sets of genes, an increased signal of rare HI variants was demonstrated predominantly for primary immunodeficiency genes and the entire set of genes associated with immune diseases, as well as for genes associated with respiratory diseases. We advocate taking the ideas of the omnigenic model into account in COVID-19 studies.

Multitrait genome-wide analyses identify new susceptibility loci and candidate drugs to primary sclerosing cholangitis

Primary sclerosing cholangitis (PSC) is a rare autoimmune bile duct disease that is strongly associated with immune-mediated disorders. In this study, we implemented multitrait joint analyses to genome-wide association summary statistics of PSC and numerous clinical and epidemiological traits to estimate the genetic contribution of each trait and genetic correlations between traits and to identify new lead PSC risk-associated loci. We identified seven new loci that have not been previously reported and one new independent lead variant in the previously reported locus. Functional annotation and fine-mapping nominated several potential susceptibility genes such as MANBA and IRF5. Network-based in silico drug efficacy screening provided candidate agents for further study of pharmacological effect in PSC.

Quantifying portable genetic effects and improving cross-ancestry genetic prediction with GWAS summary statistics

Polygenic risk scores (PRS) calculated from genome-wide association studies (GWAS) of Europeans are known to have substantially reduced predictive accuracy in non-European populations, limiting their clinical utility and raising concerns about health disparities across ancestral populations. Here, we introduce a statistical framework named X-Wing to improve predictive performance in ancestrally diverse populations. X-Wing quantifies local genetic correlations for complex traits between populations, employs an annotation-dependent estimation procedure to amplify correlated genetic effects between populations, and combines multiple population-specific PRS into a unified score with GWAS summary statistics alone as input. Through extensive benchmarking, we demonstrate that X-Wing pinpoints portable genetic effects and substantially improves PRS performance in non-European populations, showing 14.1%-119.1% relative gain in predictive R2 compared to state-of-the-art methods based on GWAS summary statistics. Overall, X-Wing addresses critical limitations in existing approaches and may have broad applications in cross-population polygenic risk prediction.

Powerful, scalable and resource-efficient meta-analysis of rare variant associations in large whole genome sequencing studies

Meta-analysis of whole genome sequencing/whole exome sequencing (WGS/WES) studies provides an attractive solution to the problem of collecting large sample sizes for discovering rare variants associated with complex phenotypes. Existing rare variant meta-analysis approaches are not scalable to biobank-scale WGS data. Here we present MetaSTAAR, a powerful and resource-efficient rare variant meta-analysis framework for large-scale WGS/WES studies. MetaSTAAR accounts for relatedness and population structure, can analyze both quantitative and dichotomous traits and boosts the power of rare variant tests by incorporating multiple variant functional annotations. Through meta-analysis of four lipid traits in 30,138 ancestrally diverse samples from 14 studies of the Trans Omics for Precision Medicine (TOPMed) Program, we show that MetaSTAAR performs rare variant meta-analysis at scale and produces results comparable to using pooled data. Additionally, we identified several conditionally significant rare variant associations with lipid traits. We further demonstrate that MetaSTAAR is scalable to biobank-scale cohorts through meta-analysis of TOPMed WGS data and UK Biobank WES data of ~200,000 samples.