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

Mineralocorticoid Receptor and Aldosterone: Interaction Between NR3C2 Genetic Variants, Sex, and Age in a Mixed Cohort

CONTEXT

Hypertension, a prevalent cardiovascular risk, often involves dysregulated aldosterone and its interaction with the mineralocorticoid receptor (MR). Experimental designs in animal models and human cohorts have demonstrated a sex and age dependency of aldosterone secretion that expands our pathophysiologic understanding.

OBJECTIVE

This study explores the genetic variation of NR3C2, which encodes MR, in relation to aldosterone, considering age, sex, and race.

METHODS

Incorporating 720 Caucasians and 145 Africans from the HyperPATH cohort, we investigated the impact of rs4835490, a single nucleotide risk allele variant, on aldosterone levels and vasculature.

RESULTS

Notably, a significant association between rs4835490 and plasma aldosterone under liberal salt conditions emerged in individuals of European ancestry (P = .0002). Homozygous carriers of the risk A allele exhibited elevated plasma aldosterone levels (AA = 8.1 ± .9 vs GG = 4.9 ± .5 ng/dL). Additionally, aldosterone activation through posture (P = .025) and urinary excretion (P = .0122) showed notable associations. Moreover, genetic interactions with race, sex, and age were observed. Caucasian females under 50 years displayed higher plasma aldosterone, urine aldosterone, and posture aldosterone with the AA genotype compared to females over 50 years, suggesting a potential connection with menopausal or estrogen influences. Interestingly, such age-dependent interactions were absent in the African cohort.

CONCLUSION

Our study highlights the significance of the NR3C2 genetic variation and its interplay with age, sex, and race in aldosterone activation. The findings point toward an estrogen-modulating effect on MR activation, particularly in women, underlining the role of aldosterone dysregulation in hypertension development. This insight advances our comprehension of hypertension's complexities and opens avenues for personalized interventions. Clinical Trial Registration Number: NCT03029806 (registered January 24, 2017).

Machine learning model identifies genetic predictors of cisplatin-induced ototoxicity in CERS6 and TLR4

BACKGROUND

Cisplatin-induced ototoxicity remains a significant concern in pediatric cancer treatment due to its permanent impact on quality of life. Previously, genetic association analyses have been performed to detect genetic variants associated with this adverse reaction.

METHODS

In this study, a combination of interpretable neural networks and Generative Adversarial Networks (GANs) was employed to identify genetic markers associated with cisplatin-induced ototoxicity. The applied method, BRI-Net, incorporates biological domain knowledge to define the network structure and employs adversarial training to learn an unbiased representation of the data, which is robust to known confounders. Leveraging genomic data from a cohort of 362 cisplatin-treated pediatric cancer patients recruited by the CPNDS (Canadian Pharmacogenomics Network for Drug Safety), this model revealed two statistically significant single nucleotide polymorphisms to be associated with cisplatin-induced ototoxicity.

RESULTS

Two markers within the CERS6 (rs13022792, p-value: 3 × 10-4) and TLR4 (rs10759932, p-value: 7 × 10-4) genes were associated with this cisplatin-induced adverse reaction. CERS6, a ceramide synthase, contributes to elevated ceramide levels, a known initiator of apoptotic signals in mouse models of inner ear hair cells. TLR4, a pattern-recognition protein, initiates inflammation in response to cisplatin, and reduced TLR4 expression has been shown in murine hair cells to confer protection from ototoxicity.

CONCLUSION

Overall, these findings provide a foundation for understanding the genetic landscape of cisplatin-induced ototoxicity, with implications for improving patient care and treatment outcomes.

CAUSALdb2: an updated database for causal variants of complex traits

Unraveling the causal variants from genome wide association studies (GWASs) is pivotal for understanding genetic underpinnings of complex traits and diseases. Despite continuous efforts, tools to refine and prioritize GWAS signals need enhancement to address the direct causal implications of genetic variations. To overcome challenges related to statistical fine-mapping in identifying causal variants, CAUSALdb has been updated with novel features and comprehensive datasets, morphing into CAUSALdb2. This expanded repository integrates 15 057 updated GWAS summary statistics across 10 839 unique traits and implements both LD-based and LD-free fine-mapping approaches, including innovative applications of approximate Bayes Factor and SuSiE. Additionally, by incorporating larger LD reference panels such as TOPMED and UK Biobank, and integrating functional annotations via PolyFun, CAUSALdb2 enhances the accuracy and context of fine-mapping results. The database now supports interrogation of additional causal signals and offers sophisticated visualizations to aid researchers in deciphering complex genetic architectures. By facilitating a deeper and more precise characterisation of causal variants, CAUSALdb2 serves as a crucial tool for advancing the genetic analysis of complex diseases. Available freely, CAUSALdb2 continues to set benchmarks in the post-GWAS era, fostering the development of targeted diagnostics and therapeutics derived from responsible genetic research. Explore these advancements at http://mulinlab.org/causaldb.

Best practices for germline variant and DNA methylation analysis of second- and third-generation sequencing data

This comprehensive review provides insights and suggested strategies for the analysis of germline variants using second- and third-generation sequencing technologies (SGS and TGS). It addresses the critical stages of data processing, starting from alignment and preprocessing to quality control, variant calling, and the removal of artifacts. The document emphasized the importance of meticulous data handling, highlighting advanced methodologies for annotating variants and identifying structural variations and methylated DNA sites. Special attention is given to the inspection of problematic variants, a step that is crucial for ensuring the accuracy of the analysis, particularly in clinical settings where genetic diagnostics can inform patient care. Additionally, the document covers the use of various bioinformatics tools and software that enhance the precision and reliability of these analyses. It outlines best practices for the annotation of variants, including considerations for problematic genetic alterations such as those in the human leukocyte antigen region, runs of homozygosity, and mitochondrial DNA alterations. The document also explores the complexities associated with identifying structural variants and copy number variations, underscoring the challenges posed by these large-scale genomic alterations. The objective is to offer a comprehensive framework for researchers and clinicians, ensuring that genetic analyses conducted with SGS and TGS are both accurate and reproducible. By following these best practices, the document aims to increase the diagnostic accuracy for hereditary diseases, facilitating early diagnosis, prevention, and personalized treatment strategies. This review serves as a valuable resource for both novices and experts in the field, providing insights into the latest advancements and methodologies in genetic analysis. It also aims to encourage the adoption of these practices in diverse research and clinical contexts, promoting consistency and reliability across studies.

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.

A Genome-Wide Association Study Suggests New Susceptibility Loci for Primary Antiphospholipid Syndrome

OBJECTIVE

Primary antiphospholipid syndrome (PAPS) is a rare autoimmune disease characterized by the presence of antiphospholipid antibodies and the occurrence of thrombotic events and pregnancy complications. Our study aimed to identify novel genetic susceptibility loci associated with PAPS.

METHODS

We performed a genome-wide association study comprising 5,485 individuals (482 affected individuals) of European ancestry. Significant and suggestive independent variants from a meta-analysis of approximately 7 million variants were evaluated for functional and biological process enrichment. The genetic risk variability for PAPS in different populations was also assessed. Hierarchical clustering, Mahalanobis distance, and Dirichlet Process Mixtures with uncertainty clustering methods were used to assess genetic similarities between PAPS and other immune-mediated diseases.

RESULTS

We revealed genetic associations with PAPS in a regulatory locus within the HLA class II region near HLA-DRA and in STAT1-STAT4 with a genome-wide level of significance; 34 additional suggestive genetic susceptibility loci for PAPS were also identified. The disease risk allele near HLA-DRA is associated with overexpression of HLA-DRB6, HLA-DRB9, HLA-DQA2, and HLA-DQB2 in immune cells, vascular tissue, and nervous tissue. This association is independent of the association between PAPS and HLA-DRB1*1302. Functional analyses highlighted immune-related pathways in PAPS-associated loci. The comparison with other immune-mediated diseases revealed a close genetic relatedness to neuromyelitis optica, systemic sclerosis, and Sjögren syndrome, suggesting co-localized causal variations close to STAT1-STAT4, TNPO3, and BLK.

CONCLUSION

This study represents a comprehensive large-scale genetic analysis for PAPS and provides new insights into the genetic basis and pathophysiology of this rare disease.

xQTLatlas: a comprehensive resource for human cellular-resolution multi-omics genetic regulatory landscape

Understanding how genetic variants influence molecular phenotypes in different cellular contexts is crucial for elucidating the molecular and cellular mechanisms behind complex traits, which in turn has spurred significant advances in research into molecular quantitative trait locus (xQTL) at the cellular level. With the rapid proliferation of data, there is a critical need for a comprehensive and accessible platform to integrate this information. To meet this need, we developed xQTLatlas (http://www.hitxqtl.org.cn/), a database that provides a multi-omics genetic regulatory landscape at cellular resolution. xQTLatlas compiles xQTL summary statistics from 151 cell types and 339 cell states across 55 human tissues. It organizes these data into 20 xQTL types, based on four distinct discovery strategies, and spans 13 molecular phenotypes. Each entry in xQTLatlas is meticulously annotated with comprehensive metadata, including the origin of the tissue, cell type, cell state and the QTL discovery strategies utilized. Additionally, xQTLatlas features multiscale data exploration tools and a suite of interactive visualizations, facilitating in-depth analysis of cell-level xQTL. xQTLatlas provides a valuable resource for deepening our understanding of the impact of functional variants on molecular phenotypes in different cellular environments, thereby facilitating extensive research efforts.

Identification and verification of disulfidptosis-related genes in sepsis-induced acute lung injury

Background

Sepsis-induced acute lung injury (ALI) is a common and serious complication of sepsis that eventually progresses to life-threatening hypoxemia. Disulfidptosis is a newly discovered type of cell death associated with the pathogenesis of different diseases. This study investigated the potential association between sepsis-induced acute lung injury and disulfidptosis by bioinformatics analysis.

Methods

In order to identify differentially expressed genes (DEGs) linked to sepsis, we screened appropriate data sets from the GEO database and carried out differential analysis. The key genes shared by DEGs and 39 disulfidptosis-related genes were identified: ACSL4 and MYL6 mRNA levels of key genes were detected in different datasets. We then used a series of bioinformatics analysis techniques, such as immune cell infiltration analysis, protein-protein interaction (PPI) network, genetic regulatory network, and receiver operating characteristic (ROC), to investigate the possible relationship between key genes and sepsis. Then, experimental verification was obtained for changes in key genes in sepsis-induced acute lung injury. Finally, to investigate the relationship between genetic variants of MYL6 or ACSL4 and sepsis, Mendelian randomization (MR) analysis was applied.

Results

Two key genes were found in this investigation: myosin light chain 6 (MYL6) and Acyl-CoA synthetase long-chain family member 4 (ACSL4). We verified increased mRNA levels of key genes in training datasets. Immune cell infiltration analysis showed that key genes were associated with multiple immune cell levels. Building the PPI network between MYL6 and ACSL4 allowed us to determine that their related genes had distinct biological functions. The co-expression genes of key genes were involved in different genetic regulatory networks. In addition, both the training and validation datasets confirmed the diagnostic capabilities of key genes by using ROC curves. Additionally, both in vivo and in vitro experiments confirmed that the mRNA levels of ACSL4 and MYL6 in sepsis-induced acute lung injury were consistent with the results of bioinformatics analysis. Finally, MR analysis revealed a causal relationship between MYL6 and sepsis.

Conclusion

We have discovered and confirmed that the key genes ACSL4 and MYL6, which are linked to disulfidptosis in sepsis-induced acute lung injury, may be useful in the diagnosis and management of septic acute lung injury.

A functional variant rs10409772 in FUT6 promoter regulates colorectal cancer progression through PKA/CREB signaling

Fucosyltransferase 6 (FUT6) is overexpressed in colorectal cancer tissue according to TCGA samples and immunohistochemistry results of a tissue microarray. FUT6 effects cell migration, tumor formation and proliferation of colorectal cancer cells in different essays. FUT6 promotes cancer cell proliferation in vitro and colorectal tumorigenesis in vivo by upregulating PKA/CREB pathway activation. Moreover, FUT6 expression is regulated by rs10409772 shown in the luciferase essays, a single nucleotide polymorphism in the promoter of FUT6. Our study suggests that elevated expression of FUT6 promotes PKA/CREB signaling, which in turn augments colorectal carcinogenesis, indicating a potential therapeutic target for colorectal cancer patients with increased FUT6 expression.

Estimation of genetic variation in vitiligo associated genes: Population genomics perspective

BACKGROUND

Vitiligo is an auto-immune progressive depigmentation disorder of the skin due to loss of melanocytes. Genetic risk is one of the important factors for development of vitiligo. Preponderance of vitiligo in certain ethnicities is known which can be analysed by understanding the distribution of allele frequencies across normal populations. Earlier GWAS identified 108 risk alleles for vitiligo in Europeans and East Asians. In this study, 64 of these risk alleles were used for analysing their enrichment and depletion across populations (1000 Genomes Project and IndiGen) with reference to 1000 Genomes dataset. Genetic risk scores were calculated and Fisher's exact test was performed to understand statistical significance of their variation in each population with respect to 1000 Genomes dataset as reference. In addition to SNPs reported in GWAS, significant variation in allele frequencies of 1079 vitiligo-related genes were also analysed. Two-tailed Chi-square test and Bonferroni's multiple adjustment values along with fixation index (≥ 0.5) and minimum allele frequency (≥ 0.05) were calculated and used to prioritise the variants based on pairwise comparison across populations.

RESULTS

Risk alleles rs1043101 and rs10768122 belong to 3 prime UTR of glutamate receptor gene SLC1A2 are found to be highly enriched in the South Asian population when compared with the 'global normal' population. Intron variant rs4766578 (ATXN2) was found to be deleted in SAS, EAS and AFR and enriched in EUR and AMR1. This risk allele is found to be under positive selection in SAS, AMR1 and EUR. From the ancillary vitiligo gene list, nonsynonymous variant rs16891982 was found to be enriched in the European and the Admixed American populations and depleted in all others. rs2279238 and rs11039155 belonging to the LXR-α gene involved in regulation of metalloproteinase 2 and 9 (melanocyte precursors) were found to be associated with vitiligo in the North Indian population (in earlier study).

CONCLUSION

The differential enrichment/depletion profile of the risk alleles provides insight into the underlying inter-population variations. This would provide clues towards prioritisation of SNPs associated with vitiligo thereby elucidating its preponderance in different ethnic groups.

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.

Epigenetic associations of GPNMB rs199347 variant with alcohol consumption in Parkinson's disease

Introduction

Alcohol consumption can induce a neuroinflammatory response and contribute to the progression of neurodegeneration. However, its association with Parkinson's disease (PD), the second most common neurodegenerative disorder, remains undetermined. Recent studies suggest that the glycoprotein non-metastatic melanoma protein B (GPNMB) is a potential biomarker for PD. We evaluated the association of rs199347, a variant of the GPNMB gene, with alcohol consumption and methylation upstream of GPNMB.

Methods

We retrieved genetic and DNA methylation data obtained from participants enrolled in the Taiwan Biobank (TWB) between 2008 and 2016. After excluding individuals with incomplete or missing information about potential PD risk factors, we included 1,357 participants in our final analyses. We used multiple linear regression to assess the association of GPNMB rs199347 and chronic alcohol consumption (and other potential risk factors) with GPNMB cg17274742 methylation.

Results

There was no difference between the distribution of GPNMB rs199347 genotypes between chronic alcohol consumers and the other study participants. A significant interaction was observed between the GPNMB rs199347 variant and alcohol consumption (p = 0.0102) concerning cg17274742 methylation. Compared to non-chronic alcohol consumers with the AA genotype, alcohol drinkers with the rs199347 GG genotype had significantly lower levels (hypomethylation) of cg17274742 (p = 0.0187).

Conclusion

Alcohol consumption among individuals with the rs199347 GG genotype was associated with lower levels of cg17274742 methylation, which could increase expression of the GPNMB gene, an important neuroinflammatory-related risk gene for PD.

REEV: review, evaluate and explain variants

In the era of high throughput sequencing, special software is required for the clinical evaluation of genetic variants. We developed REEV (Review, Evaluate and Explain Variants), a user-friendly platform for clinicians and researchers in the field of rare disease genetics. Supporting data was aggregated from public data sources. We compared REEV with seven other tools for clinical variant evaluation. REEV (semi-)automatically fills individual ACMG criteria facilitating variant interpretation. REEV can store disease and phenotype data related to a case to use these for phenotype similarity measures. Users can create public permanent links for individual variants that can be saved as browser bookmarks and shared. REEV may help in the fast diagnostic assessment of genetic variants in a clinical as well as in a research context. REEV (https://reev.bihealth.org/) is free and open to all users and there is no login requirement.

Association of NID2 SNPs with Glioma Risk and Prognosis in the Chinese Population

Glioma is the most common primary intracranial tumor with high mortality and poor prognosis. The purpose of this study was to investigate how single-nucleotide polymorphisms (SNPs) of the NID2 gene affect glioma risk and prognosis. Four candidate SNPs of NID2 in 529 glioma patients and 478 healthy controls were successfully genotyped by Agena MassARRAY mass spectrometer. Logistic regression was utilized to assess the associations between NID2 SNPs and glioma risk under different genetic models. Furthermore, the relationship between risk-related SNPs in NID2 and the prognosis of glioma patients was explored through Kaplan-Meier (KM) survival curve and Cox proportional hazard regression analysis. The results showed that rs11846847 (OR 1.24, p = 0.017) and rs1874569 (OR 1.22, p = 0.026) were significantly associated with an increased risk of glioma, and rs11846847 also had a risk-increasing effect on glioma in participants ≤ 40 years old. The interaction model of rs11846847 and rs1874569 could be more suitable for forecasting glioma risk. We also discovered a significant association between rs1874569 and poor prognosis in glioma patients (HR 1.32, p = 0.039) and especially CC genotype was relevant to shorter overall survival (OS) and progression-free survival (PFS) in patients with high-grade glioma. Additionally, the study demonstrated that gross total resection or chemotherapy improve glioma prognosis in the Chinese Han population. This study is the first to provide evidence for the association of NID2 SNPs with glioma risk and prognosis, suggesting that NID2 variants might be potential factors for glioma.

Identification of an novel genetic variant associated with osteoporosis: insights from the Taiwan Biobank Study

Purpose

The purpose of this study was to identify new independent significant SNPs associated with osteoporosis using data from the Taiwan Biobank (TWBB).

Material and Methods

The dataset was divided into discovery (60%) and replication (40%) subsets. Following data quality control, genome-wide association study (GWAS) analysis was performed, adjusting for sex, age, and the top 5 principal components, employing the Scalable and Accurate Implementation of the Generalized mixed model approach. This was followed by a meta-analysis of TWBB1 and TWBB2. The Functional Mapping and Annotation (FUMA) platform was used to identify osteoporosis-associated loci. Manhattan and quantile-quantile plots were generated using the FUMA platform to visualize the results. Independent significant SNPs were selected based on genome-wide significance (P < 5 × 10-8) and independence from each other (r2 < 0.6) within a 1 Mb window. Positional, eQTL(expression quantitative trait locus), and Chromatin interaction mapping were used to map SNPs to genes.

Results

A total of 29 084 individuals (3154 osteoporosis cases and 25 930 controls) were used for GWAS analysis (TWBB1 data), and 18 918 individuals (1917 cases and 17 001 controls) were utilized for replication studies (TWBB2 data). We identified a new independent significant SNP for osteoporosis in TWBB1, with the lead SNP rs76140829 (minor allele frequency = 0.055, P-value = 1.15 × 10-08). Replication of the association was performed in TWBB2, yielding a P-value of 6.56 × 10-3. The meta-analysis of TWBB1 and TWBB2 data demonstrated a highly significant association for SNP rs76140829 (P-value = 7.52 × 10-10). In the positional mapping of rs76140829, 6 genes (HABP2, RP11-481H12.1, RNU7-165P, RP11-139 K1.2, RP11-57H14.3, and RP11-214 N15.5) were identified through chromatin interaction mapping in mesenchymal stem cells.

Conclusions

Our GWAS analysis using the Taiwan Biobank dataset unveils rs76140829 in the VTI1A gene as a key risk variant associated with osteoporosis. This finding expands our understanding of the genetic basis of osteoporosis and highlights the potential regulatory role of this SNP in mesenchymal stem cells.

Novel genetic association of the Furin gene polymorphism rs1981458 with COVID-19 severity among Indian populations

SARS CoV-2, the causative agent for the ongoing COVID-19 pandemic, it enters the host cell by activating the ACE2 receptor with the help of two proteasesi.e., Furin and TMPRSS2. Therefore, variations in these genes may account for differential susceptibility and severity between populations. Previous studies have shown that the role of ACE2 and TMPRSS2 gene variants in understanding COVID-19 susceptibility among Indian populations. Nevertheless, a knowledge gap exists concerning the COVID-19 susceptibility of Furin gene variants among diverse South Asian ethnic groups. Investigating the role of Furin gene variants and their global phylogeographic structure is essential to comprehensively understanding COVID-19 susceptibility in these populations. We have used 450 samples from diverse Indian states and performed linear regression to analyse the Furin gene variant's with COVID-19 Case Fatality Rate (CFR) that could be epidemiologically associated with disease severity outcomes. Associated genetic variants were further evaluated for their expression and regulatory potential through various Insilco analyses. Additionally, we examined the Furin gene using next-generation sequencing (NGS) data from 393 diverse global samples, with a particular emphasis on South Asia, to investigate its Phylogeographic structure among diverse world populations. We found a significant positive association for the SNP rs1981458 with COVID-19 CFR (p < 0.05) among diverse Indian populations at different timelines of the first and second waves. Further, QTL and other regulatory analyses showed various significant associations for positive regulatory roles of rs1981458 and Furin gene, mainly in Immune cells and virus infection process, highlighting their role in host immunity and viral assembly and processing. The Furin protein-protein interaction suggested that COVID-19 may contribute to Pulmonary arterial hypertension via a typical inflammation mechanism. The phylogeographic architecture of the Furin gene demonstrated a closer genetic affinity of South Asia with West Eurasian populations. Therefore, it is worth proposing that for the Furin gene, the COVID-19 susceptibility of South Asians will be more similar to the West Eurasian population. Our previous studies on the ACE2 and TMPRSS2 genes showed genetic affinity of South Asian with East Eurasians and West Eurasians, respectively. Therefore, with the collective information from these three important genes (ACE2, TMPRSS2 and Furin) we modelled COVID-19 susceptibilityof South Asia in between these two major ancestries with an inclination towards West Eurasia. In conclusion, this study, for the first time, concluded the role of rs1981458 in COVID-19 severity among the Indian population and outlined its regulatory potential.This study also highlights that the genetic structure for COVID-19 susceptibilityof South Asia is distinct, however, inclined to the West Eurasian population. We believe this insight may be utilised as a genetic biomarker to identify vulnerable populations, which might be directly relevant for developing policies and allocating resources more effectively during an epidemic.

Causal relationship between gut microbiota and chronic renal failure: a two-sample Mendelian randomization study

Background

Observational studies and some experimental investigations have indicated that gut microbiota are closely associated with the incidence and progression of chronic renal failure. However, the causal relationship between gut microbiota and chronic renal failure remains unclear. The present study employs a two-sample Mendelian randomization approach to infer the causal relationship between gut microbiota and chronic renal failure at the genetic level. This research aims to determine whether there is a causal effect of gut microbiota on the risk of chronic renal failure, aiming to provide new evidence to support targeted gut therapy for the treatment of chronic renal failure.

Methods

Employing genome-wide association study (GWAS) data from the public MiBioGen and IEU OpenGWAS platform, a two-sample Mendelian randomization analysis was conducted. The causal relationship between gut microbiota and chronic renal failure was inferred using five different methods: Inverse Variance Weighted, MR-Egger, Weighted Median, Simple Mode, and Weighted Mode. The study incorporated sensitivity analyses that encompassed evaluations for pleiotropy and heterogeneity. Subsequently, the results of the Mendelian randomization analysis underwent a stringent correction for multiple testing, employing the False Discovery Rate method to enhance the validity of our findings.

Results

According to the results from the Inverse Variance Weighted method, seven bacterial genera show a significant association with the outcome variable chronic renal failure. Of these, Ruminococcus (gauvreauii group) (OR = 0.82, 95% CI = 0.71-0.94, p = 0.004) may act as a protective factor against chronic renal failure, while the genera Escherichia-Shigella (OR = 1.22, 95% CI = 1.08-1.38, p = 0.001), Lactococcus (OR = 1.1, 95% CI = 1.02-1.19, p = 0.013), Odoribacter (OR = 1.23, 95% CI = 1.03-1.49, p = 0.026), Enterorhabdus (OR = 1.14, 95% CI = 1.00-1.29, p = 0.047), Eubacterium (eligens group) (OR = 1.18, 95% CI = 1.02-1.37, p = 0.024), and Howardella (OR = 1.18, 95% CI = 1.09-1.28, p < 0.001) may be risk factors for chronic renal failure. However, after correction for multiple comparisons using False Discovery Rate, only the associations with Escherichia-Shigella and Howardella remain significant, indicating that the other genera have suggestive associations. Sensitivity analyses did not reveal any pleiotropy or heterogeneity.

Conclusion

Our two-sample Mendelian randomization study suggests that the genera Escherichia-Shigella and Howardella are risk factors for chronic renal failure, and they may serve as potential targets for future therapeutic interventions. However, the exact mechanisms of action are not yet clear, necessitating further research to elucidate their precise roles fully.

A perspective on genetic and polygenic risk scores-advances and limitations and overview of associated tools

Polygenetic Risk Scores are used to evaluate an individual's vulnerability to developing specific diseases or conditions based on their genetic composition, by taking into account numerous genetic variations. This article provides an overview of the concept of Polygenic Risk Scores (PRS). We elucidate the historical advancements of PRS, their advantages and shortcomings in comparison with other predictive methods, and discuss their conceptual limitations in light of the complexity of biological systems. Furthermore, we provide a survey of published tools for computing PRS and associated resources. The various tools and software packages are categorized based on their technical utility for users or prospective developers. Understanding the array of available tools and their limitations is crucial for accurately assessing and predicting disease risks, facilitating early interventions, and guiding personalized healthcare decisions. Additionally, we also identify potential new avenues for future bioinformatic analyzes and advancements related to PRS.

Association of Vitamin D Deficiency and Vitamin D Receptor Gene Polymorphisms with Type 1 Diabetes Risk: A South Indian Familial Study

Objective

Vitamin D is a potent immune modulator and is associated with autoimmune diseases, including type 1 diabetes (T1D). The vitamin D levels and its receptor gene polymorphisms together in T1D are not yet investigated in the South Indian population. The present study focused on exploring the significance of vitamin D levels and vitamin D receptor (VDR) gene polymorphisms with the risk of developing T1D in the South Indian population.

Methods

Patients with T1D and unaffected first-degree relatives (FDRs) were included in this study. Genotyping of VDR polymorphisms at four different loci (FokI- F/f, BsmI- B/b, TaqI- T/t, and ApaI- A/a) was assessed through the amplification refractive mutation system-polymerase chain reaction method. Serum vitamin D levels were measured in 98 T1D patients and 75 age- and sex-matched siblings.

Results

A total of 120 patients with T1D and 214 FDRs were included. Vitamin D deficiency (VDD) was observed in a higher proportion of T1D patients than in controls (52% vs. 32%; p<0.03). The frequency of the FokI-FF genotype was significantly higher [odds ratio (OR)=1.66; p<0.03] in T1D patients conferring a susceptible association with the disease. Nevertheless, the increased frequency of heterozygous Ff genotype (OR=0.57; p<0.02) among controls may confer a protective association with T1D. Furthermore, the transmission disequilibrium test revealed over-transmission of ApaI-A (T: U=15/5; p<0.006) and BsmI-B alleles (T: U=17/5; p<0.01) and under-transmission of BsmI-b/ApaI-a/TaqI-T haplotype (T: U=5.4/14.4; p=0.04) from parents to T1D patients.

Conclusion

The present study concludes that VDD is the major contributing risk factor to T1D development in the South Indian population. Furthermore, the FokI-FF genotype, BsmI-B, and ApaI-A alleles were positively associated with T1D. In contrast, the FokI-Ff genotype and BsmI-b/ApaI-a/TaqI-T haplotype were negatively associated with T1D.

Genetic susceptibility and causal pathway analysis of eye disorders coexisting in multiple sclerosis

Introduction

The comorbidity of optic neuritis with multiple sclerosis has been well recognized. However, the causal association between multiple sclerosis and optic neuritis, as well as other eye disorders, remains incompletely understood. To address these gaps, we investigated the genetically relationship between multiple sclerosis and eye disorders, and explored potential drugs.

Methods

In order to elucidate the genetic susceptibility and causal links between multiple sclerosis and eye disorders, we performed two-sample Mendelian randomization analyses to examine the causality between multiple sclerosis and eye disorders. Additionally, causal single-nucleotide polymorphisms were annotated and searched for expression quantitative trait loci data. Pathway enrichment analysis was performed to identify the possible mechanisms responsible for the eye disorders coexisting with multiple sclerosis. Potential therapeutic chemicals were also explored using the Cytoscape.

Results

Mendelian randomization analysis revealed that multiple sclerosis increased the incidence of optic neuritis while reducing the likelihood of concurrent of cataract and macular degeneration. Gene Ontology enrichment analysis implicated that lymphocyte proliferation, activation and antigen processing as potential contributors to the pathogenesis of eye disorders coexisting with multiple sclerosis. Furthermore, pharmaceutical agents traditionally employed for allograft rejection exhibited promising therapeutic potential for the eye disorders coexisting with multiple sclerosis.

Discussion

Multiple sclerosis genetically contributes to the development of optic neuritis while mitigating the concurrent occurrence of cataract and macular degeneration. Further research is needed to validate these findings and explore additional mechanisms underlying the comorbidity of multiple sclerosis and eye disorders.

Exploring noncoding variants in genetic diseases: from detection to functional insights

Previous studies on genetic diseases predominantly focused on protein-coding variations, overlooking the vast noncoding regions in the human genome. The development of high-throughput sequencing technologies and functional genomics tools has enabled the systematic identification of functional noncoding variants. These variants can impact gene expression, regulation, and chromatin conformation, thereby contributing to disease pathogenesis. Understanding the mechanisms that underlie the impact of noncoding variants on genetic diseases is indispensable for the development of precisely targeted therapies and the implementation of personalized medicine strategies. The intricacies of noncoding regions introduce a multitude of challenges and research opportunities. In this review, we introduce a spectrum of noncoding variants involved in genetic diseases, along with research strategies and advanced technologies for their precise identification and in-depth understanding of the complexity of the noncoding genome. We will delve into the research challenges and propose potential solutions for unraveling the genetic basis of rare and complex 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/.

A genome-wide association study suggests new susceptibility loci for primary antiphospholipid syndrome

Objectives

Primary antiphospholipid syndrome (PAPS) is a rare autoimmune disease characterized by the presence of antiphospholipid antibodies and the occurrence of thrombotic events and pregnancy complications. Our study aimed to identify novel genetic susceptibility loci associated with PAPS.

Methods

We performed a genome-wide association study comprising 5,485 individuals (482 affected individuals) of European ancestry. Significant and suggestive independent variants from a meta-analysis of approximately 7 million variants were evaluated for functional and biological process enrichment. The genetic risk variability for PAPS in different populations was also assessed. Hierarchical clustering, Mahalanobis distance, and Dirichlet Process Mixtures with uncertainty clustering methods were used to assess genetic similarities between PAPS and other immune-mediated diseases.

Results

We revealed genetic associations with PAPS in a regulatory locus within the HLA class II region near HLA-DRA and in STAT4 with a genome-wide level of significance. 34 additional suggestive genetic susceptibility loci for PAPS were also identified. The disease risk allele in the HLA class II locus is associated with overexpression of HLA-DRB6 , HLA-DRB9 , HLA-DPB2 , HLA-DQA2 and HLA-DQB2 , and is independent of the association between PAPS and HLA-DRB1*1302 . Functional analyses highlighted immune and nervous system related pathways in PAPS-associated loci. The comparison with other immune-mediated diseases revealed a close genetic relatedness to neuromyelitis optica, systemic sclerosis, and Sjögren's syndrome, suggesting colocalized causal variations close to STAT4 , TNPO3 , and BLK .

Conclusions

This study represents a comprehensive large-scale genetic analysis for PAPS and provides new insights into the genetic basis and pathophysiology of this rare disease.

Genetic polymorphisms in CYP2B6 may be associated with lung cancer risk in the Chinese Han population

BACKGROUND

Our study aimed to elucidate the association between single nucleotide polymorphisms (SNPs) in CYP2B6 gene and susceptibility to lung cancer (LC).

METHODS

Five SNPs in CYP2B6 were genotyped in Chinese Han population (507 cases and 505 controls) utilizing Agena MassARRAY. The relationship between these SNPs and LC susceptibility was assessed using odds ratios, 95% confidence intervals, and χ2 tests. Additionally, multifactor dimensionality reduction was employed to analyze SNP-SNP interactions. Bioinformatics methods were applied to investigate the function of these SNPs.

RESULTS

We found that rs2099361 was associated with an increased susceptibility to LC in the codominant model (OR = 1.31, p = 0.045). Stratification analysis revealed the allele G at rs4803418 and the allele T at rs4803420 of CYP2B6 (BMI >24 kg/m2) were significantly linked to decreased susceptibility of LC. Conversely, the allele C at rs12979270 (BMI >24 kg/m2) showed increased susceptibility to LC. Moreover, a robust redundant relationship between rs12979270 and rs4803420 was identified in the study. According to the VannoPortal database, we found that rs4803420, rs12979270 and rs2099361 may modulate the binding affinity of LMNB1, SP1 and HDAC2, respectively.

CONCLUSIONS

Our results suggest that SNPs in the CYP2B6 gene play crucial roles in LC susceptibility.

Association Analyses of Predicted Loss-of-Function Variants Prioritized 15 Genes as Blood Pressure Regulators

BACKGROUND

Hypertension, clinically defined by elevated blood pressure (BP), is an important cause of mortality and morbidity worldwide. Many risk factors for hypertension are known, including a positive family history, which suggests that genetics contribute to interindividual BP variation. Genome-wide association studies (GWAS) have identified > 1000 loci associated with BP, yet the identity of the genes responsible for these associations remains largely unknown.

METHODS

To pinpoint genes that causally affect variation of BP in humans, we analyzed predicted loss-of-function (pLoF) variants in the UK Biobank whole-exome sequencing dataset (n = 454,709 participants, 6% non-European ancestry). We analyzed genetic associations between systolic or diastolic BP (SBP/DBP) and single pLoF variants (additive and recessive genetic models) as well as with the burden of very rare pLoF variants (minor allele frequency [MAF] < 0.01%).

RESULTS

Single pLoF variants in 10 genes were associated with BP (ANKDD1B, ENPEP, PNCK, BTN3A2, C1orf145 [OBSCN-AS1], CASP9, DBH, KIAA1161 [MYORG], OR4X1, and TMC3). We also found a burden of rare pLoF variants in 5 additional genes associated with BP (TTN, NOS3, FES, SMAD6, COL21A1). Except for PNCK, which is located on the X-chromosome, these genes map near variants previously associated with BP by GWAS, validating the study of pLoF variants to prioritize causal genes at GWAS loci.

CONCLUSIONS

Our study highlights 15 genes that likely modulate BP in humans, including 5 genes that harbour pLoF variants associated with lower BP.

Ocular and neural genes jointly regulate the visuospatial working memory in ADHD children

OBJECTIVE

Working memory (WM) deficits have frequently been linked to attention deficit hyperactivity disorder (ADHD). Despite previous studies suggested its high heritability, its genetic basis, especially in ADHD, remains unclear. The current study aimed to comprehensively explore the genetic basis of visual-spatial working memory (VSWM) in ADHD using wide-ranging genetic analyses.

METHODS

The current study recruited a cohort consisted of 802 ADHD individuals, all met DSM-IV ADHD diagnostic criteria. VSWM was assessed by Rey-Osterrieth complex figure test (RCFT), which is a widely used psychological test include four memory indexes: detail delayed (DD), structure delayed (SD), structure immediate (SI), detail immediate (DI). Genetic analyses were conducted at the single nucleotide polymorphism (SNP), gene, pathway, polygenic and protein network levels. Polygenic Risk Scores (PRS) were based on summary statistics of various psychiatric disorders, including ADHD, autism spectrum disorder (ASD), major depressive disorder (MDD), schizophrenia (SCZ), obsessive compulsive disorders (OCD), and substance use disorder (SUD).

RESULTS

Analyses at the single-marker level did not yield significant results (5E-08). However, the potential signals with P values less than E-05 and their mapped genes suggested the regulation of VSWM involved both ocular and neural system related genes, moreover, ADHD-related genes were also involved. The gene-based analysis found RAB11FIP1, whose encoded protein modulates several neurodevelopment processes and visual system, as significantly associated with DD scores (P = 1.96E-06, Padj = 0.036). Candidate pathway enrichment analyses (N = 53) found that forebrain neuron fate commitment significantly enriched in DD (P = 4.78E-04, Padj = 0.025), and dopamine transport enriched in SD (P = 5.90E-04, Padj = 0.031). We also observed a significant negative relationship between DD scores and ADHD PRS scores (P = 0.0025, Empirical P = 0.048).

CONCLUSIONS

Our results emphasized the joint contribution of ocular and neural genes in regulating VSWM. The study reveals a shared genetic basis between ADHD and VSWM, with GWAS indicating the involvement of ADHD-related genes in VSWM. Additionally, the PRS analysis identifies a significant relationship between ADHD-PRS and DD scores. Overall, our findings shed light on the genetic basis of VSWM deficits in ADHD, and may have important implications for future research and clinical practice.

A candidate locus in the renalase gene and susceptibility to blood pressure responses to the dietary salt

BACKGROUND

High dietary salt confers a risk of elevating blood pressure (BP) and the development of hypertension. BP to salt intake may be determined in part by individual genetic predisposition. Identifying these genetic underpinnings will enhance our understanding of the biological mechanisms of BP regulation. This study aims to assess the genetic association with salt sensitivity of BP (SSBP) within two well-phenotyped multinational cohorts.

METHODS

A total of 720 white participants from the HyperPATH consortium program were selected and genotyped using a multiethnic genotyping array. Individuals consumed two study diets containing high (>200 mEq/day) and low (<10 mEq/day) sodium content, after which SSBP, aldosterone, and plasma renin activity (PRA) were assessed in a controlled inpatient research setting.

RESULTS

A top signal (rs10887801; beta = 4.57, P  = 5.03E - 07) at the renalase gene ( RNLS ) region was significantly associated with SSBP. We also identified seven single nucleotide variants with linkage disequilibrium to the top signal at this region that comprised a significant haplotype (TCTTAGTT, P  = 0.00081). Homozygous carriers of the T-risk allele of the key single nucleotide variant had higher SSBP ( P  ≤ 0.00001) and lower PRA ( P  = 0.0076) compared with the nonrisk allele.

CONCLUSION

We identified significant associations between genetic variants of the RNLS gene and BP responses to dietary salt intervention and PRA that suggest susceptibility to volume-driven hypertension. These findings may contribute to a better understanding of the genetic mechanisms underlying BP regulation, support the role of RNLS in the pathogenesis of SSBP, and identify individuals who may be at risk from excess dietary salt intake.

Causal relationships between gut microbiota and programmed cell death protein 1/programmed cell death-ligand 1: A bidirectional Mendelian randomization study

Background

Multiple clinical studies have indicated that the gut microbiota influences the effects of immune checkpoint blockade (ICB) therapy comprising PD-1/PD-L1 inhibitors, but the causal relationship is unclear. Because of numerous confounders, many microbes related to PD-1/PD-L1 have not been identified. This study aimed to determine the causal relationship between the microbiota and PD-1/PD-L1 and identify possible biomarkers for ICB therapy.

Method

We used bidirectional two-sample Mendelian randomization with two different thresholds to explore the potential causal relationship between the microbiota and PD-1/PD-L1 and species-level microbiota GWAS to verify the result.

Result

In the primary forward analysis, genus_Holdemanella showed a negative correlation with PD-1 [βIVW = -0.25; 95% CI (-0.43 to -0.07); PFDR = 0.028] and genus_Prevotella9 showed a positive correlation with PD-1 [βIVW = 0.2; 95% CI (0.1 to 0.4); PFDR = 0.027]; order_Rhodospirillales [βIVW = 0.2; 95% CI (0.1 to 0.4); PFDR = 0.044], family_Rhodospirillaceae [βIVW = 0.2; 95% CI (0 to 0.4); PFDR = 0.032], genus_Ruminococcaceae_UCG005 [βIVW = 0.29; 95% CI (0.08 to 0.5); PFDR = 0.028], genus_Ruminococcus_gnavus_group [βIVW = 0.22; 95% CI (0.05 to 0.4); PFDR = 0.029], and genus_Coprococcus_2 [βIVW = 0.4; 95% CI (0.1 to 0.6); PFDR = 0.018] were positively correlated with PD-L1; and phylum_Firmicutes [βIVW = -0.3; 95% CI (-0.4 to -0.1); PFDR = 0.031], family_ClostridialesvadinBB60group [βIVW = -0.31; 95% CI (-0.5 to -0.11), PFDR = 0.008], family_Ruminococcaceae [βIVW = -0.33; 95% CI (-0.58 to -0.07); PFDR = 0.049], and genus_Ruminococcaceae_UCG014 [βIVW = -0.35; 95% CI (-0.57 to -0.13); PFDR = 0.006] were negatively correlated with PD-L1. The one significant species in further analysis was species_Parabacteroides_unclassified [βIVW = 0.2; 95% CI (0-0.4); PFDR = 0.029]. Heterogeneity (P &gt; 0.05) and pleiotropy (P &gt; 0.05) analyses confirmed the robustness of the MR results.

Immunogenetic variations predict immune-related adverse events for PD-1/PD-L1 inhibitors

BACKGROUND

PD-1/PD-L1 inhibitors have brought remarkable benefits but can cause profound immune-related adverse events (irAEs). The host immunogenetic background is likely to play a role in irAE susceptibility. In this study, we aimed to identify potential immunogenetic biomarkers to predict irAEs.

METHODS

Patients with solid tumours receiving PD-1/PD-L1 blockade were recruited and followed up. Genes considered pivotal contributors to tumour-immunity and autoimmune diseases were screened out via protein-protein interaction network and Cytoscape. Consequently, thirty-nine variants in eighteen genes were genotyped using the multiplex genotyping assay. Association analysis between genetic variants and irAEs as well as irAEs-free survival was performed.

RESULTS

Four immunogenetic variants as predictive biomarkers of irAEs were identified. The C allele of Mitogen-Activated Protein Kinase 1 (MAPK1) rs3810610 (odds ratio [OR] = 1.495, 95% confidence interval [CI] = 1.093-2.044, P = 0.012) was a risk predictor while the A allele of PTPRC rs6428474 (OR = 0.717, 95% CI = 0.521-0.987, P = 0.041) was a protective factor for all-grade irAEs. The A allele of ADAD1 rs17388568 (OR = 2.599, 95% CI = 1.355-4.983, P = 0.003) increased the risk while the G allele of IL6 rs1800796 (OR = 0.425, 95% CI = 0.205-0.881, P = 0.018) protected patients from high-grade irAEs. Significant immunogenetic variants reached a similar tendency in PD-1 blockade or lung cancer subgroups. In multivariate Cox regression analysis, the MAPK1 rs3810610 was an independent factor regarding all-grade irAEs-free survival (CC versus CT or TT: hazard ratio [HR] = 0.71, 95% CI = 0.52-0.99, P = 0.042). ADAD1 rs17388568 (AA versus AG or GG: HR = 0.11, 95% CI = 0.025-0.49, P = 0.004) and IL6 rs1800796 (GG or GC versus CC: HR = 3.10, 95% CI = 1.315-7.29, P = 0.01) were independent variables for high-grade irAEs-free survival.

CONCLUSION

We first identified several immunogenetic polymorphisms associated with irAEs and irAEs-free survival in PD-1/PD-L1 blockade-treated tumour patients, and they may serve as potential predictive biomarkers.

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.

Pancan-MNVQTLdb: systematic identification of multi-nucleotide variant quantitative trait loci in 33 cancer types

Multi-nucleotide variants (MNVs) are defined as clusters of two or more nearby variants existing on the same haplotype in an individual. Recent studies have identified millions of MNVs in human populations, but their functions remain largely unknown. Numerous studies have demonstrated that single-nucleotide variants could serve as quantitative trait loci (QTLs) by affecting molecular phenotypes. Therefore, we propose that MNVs can also affect molecular phenotypes by influencing regulatory elements. Using the genotype data from The Cancer Genome Atlas (TCGA), we first identified 223 759 unique MNVs in 33 cancer types. Then, to decipher the functions of these MNVs, we investigated the associations between MNVs and six molecular phenotypes, including coding gene expression, miRNA expression, lncRNA expression, alternative splicing, DNA methylation and alternative polyadenylation. As a result, we identified 1 397 821 cis-MNVQTLs and 402 381 trans-MNVQTLs. We further performed survival analysis and identified 46 173 MNVQTLs associated with patient overall survival. We also linked the MNVQTLs to genome-wide association studies (GWAS) data and identified 119 762 MNVQTLs that overlap with existing GWAS loci. Finally, we developed Pancan-MNVQTLdb (http://gong_lab.hzau.edu.cn/mnvQTLdb/) for data retrieval and download. Pancan-MNVQTLdb will help decipher the functions of MNVs in different cancer types and be an important resource for genetic and cancer research.

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.

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.

In silico analysis of upstream variants in Brazilian patients with Familial Hypercholesterolemia

Familial hypercholesterolemia (FH) is a prevalent autosomal genetic disease associated with increased risk of early cardiovascular events and death due to chronic exposure to very high levels of low-density lipoprotein cholesterol (LDL-c). Pathogenic variants in the coding regions of LDLR, APOB and PCSK9 account for most FH cases, and variants in non-coding regions maybe involved in FH as well. Variants in the upstream region of LDLR, APOB and PCSK9 were screened by targeted next-generation sequencing and their effects were explored using in silico tools. Twenty-five patients without pathogenic variants in FH-related genes were selected. 3 kb upstream regions of LDLR, APOB and PCSK9 were sequenced using the AmpliSeq (Illumina) and Miseq Reagent Nano Kit v2 (Illumina). Sequencing data were analyzed using variant discovery and functional annotation tools. Potentially regulatory variants were selected by integrating data from public databases, published data and context-dependent regulatory prediction score. Thirty-four single nucleotide variants (SNVs) in upstream regions were identified (6 in LDLR, 15 in APOB, and 13 in PCSK9). Five SNVs were prioritized as potentially regulatory variants (rs934197, rs9282606, rs36218923, rs538300761, g.55038486A>G). APOB rs934197 was previously associated with increased rate of transcription, which in silico analysis suggests that could be due to reducing binding affinity of a transcriptional repressor. Our findings highlight the importance of variant screening outside of coding regions of all relevant genes. Further functional studies are necessary to confirm that prioritized variants could impact gene regulation and contribute to the FH phenotype.

Performance evaluation of differential splicing analysis methods and splicing analytics platform construction

A proportion of previously defined benign variants or variants of uncertain significance in humans, which are challenging to identify, may induce an abnormal splicing process. An increasing number of methods have been developed to predict splicing variants, but their performance has not been completely evaluated using independent benchmarks. Here, we manually sourced ∼50 000 positive/negative splicing variants from > 8000 studies and selected the independent splicing variants to evaluate the performance of prediction methods. These methods showed different performances in recognizing splicing variants in donor and acceptor regions, reminiscent of different weight coefficient applications to predict novel splicing variants. Of these methods, 66.67% exhibited higher specificities than sensitivities, suggesting that more moderate cut-off values are necessary to distinguish splicing variants. Moreover, the high correlation and consistent prediction ratio validated the feasibility of integration of the splicing prediction method in identifying splicing variants. We developed a splicing analytics platform called SPCards, which curates splicing variants from publications and predicts splicing scores of variants in genomes. SPCards also offers variant-level and gene-level annotation information, including allele frequency, non-synonymous prediction and comprehensive functional information. SPCards is suitable for high-throughput genetic identification of splicing variants, particularly those located in non-canonical splicing regions.

Association between Polymorphisms in CFH, ARMS2, CFI, and C3 Genes and Response to Anti-VEGF Treatment in Neovascular Age-Related Macular Degeneration

Neovascular age-related macular degeneration (nAMD) is the leading cause of vision loss in the elderly. The gold standard of nAMD treatment is intravitreal injections of vascular endothelial growth factor (VEGF) inhibitors. Genetic factors may influence the response to anti-VEGF therapy and result in a high degree of response variability. The aim of the study was to evaluate the association of the polymorphisms in genes related to the complement system (rs2285714-CFI, rs10490924-ARMS2, rs2230199-C3, rs800292-CFH, and rs6677604-CFH) with nAMD its clinical features and optical coherent tomography (OCT) biomarkers of treatment response to anti-VEGF therapy. Genotyping by allele-specific PCR was performed in 193 AMD patients and 147 age-matched controls. A prospective study of the dynamics of changes in OCT biomarkers during aflibercept treatment included 110 treatment-naive patients. Allele T rs10490924 was associated with the increased risk of nAMD. For both rs800292 and rs6677604, carriage of the A allele was protective and decreased the nAMD risk. Associations of rs2230199 with central retinal thickness (CRT) and intraretinal cysts were revealed. The height of pigment epithelium detachment and the height of neuroretinal detachment were significantly higher in carriers of the minor allele of rs2285714, both at baseline and during treatment. The reduction of CRT was associated with higher CRT at baseline and the presence of the T allele of rs2285714. By the end of one-year follow-up the patients homozygous for the minor allele rs2285714 had significantly higher odds of the presence of anastomoses and loops and active neovascular membrane. Furthermore, minor allele carriers had decreased levels of complement factor I level in aqueous humor but not in the plasma, which may be due to the influence of rs2285714 on tissue-specific splicing. Our results suggest that the severity of AMD macular lesions is associated with rs2285714 and rs2230199 polymorphisms, which could be explained by their high regulatory potential. Patients with the minor allele of rs2285714 respond worse to antiangiogenic therapy.

ANANASTRA: annotation and enrichment analysis of allele-specific transcription factor binding at SNPs

We present ANANASTRA, https://ananastra.autosome.org, a web server for the identification and annotation of regulatory single-nucleotide polymorphisms (SNPs) with allele-specific binding events. ANANASTRA accepts a list of dbSNP IDs or a VCF file and reports allele-specific binding (ASB) sites of particular transcription factors or in specific cell types, highlighting those with ASBs significantly enriched at SNPs in the query list. ANANASTRA is built on top of a systematic analysis of allelic imbalance in ChIP-Seq experiments and performs the ASB enrichment test against background sets of SNPs found in the same source experiments as ASB sites but not displaying significant allelic imbalance. We illustrate ANANASTRA usage with selected case studies and expect that ANANASTRA will help to conduct the follow-up of GWAS in terms of establishing functional hypotheses and designing experimental verification.

Integrative genetic and immune cell analysis of plasma proteins in healthy donors identifies novel associations involving primary immune deficiency genes

BACKGROUND

Blood plasma proteins play an important role in immune defense against pathogens, including cytokine signaling, the complement system, and the acute-phase response. Recent large-scale studies have reported genetic (i.e., protein quantitative trait loci, pQTLs) and non-genetic factors, such as age and sex, as major determinants to inter-individual variability in immune response variation. However, the contribution of blood-cell composition to plasma protein heterogeneity has not been fully characterized and may act as a mediating factor in association studies.

METHODS

Here, we evaluated plasma protein levels from 400 unrelated healthy individuals of western European ancestry, who were stratified by sex and two decades of life (20-29 and 60-69 years), from the Milieu Intérieur cohort. We quantified 229 proteins by Luminex in a clinically certified laboratory and their levels of variation were analyzed together with 5.2 million single-nucleotide polymorphisms. With respect to non-genetic variables, we included 254 lifestyle and biochemical factors, as well as counts of seven circulating immune cell populations measured by hemogram and standardized flow cytometry.

RESULTS

Collectively, we found 152 significant associations involving 49 proteins and 20 non-genetic variables. Consistent with previous studies, age and sex showed a global, pervasive impact on plasma protein heterogeneity, while body mass index and other health status variables were among the non-genetic factors with the highest number of associations. After controlling for these covariates, we identified 100 and 12 pQTLs acting in cis and trans, respectively, collectively associated with 87 plasma proteins and including 19 novel genetic associations. Genetic factors explained the largest fraction of the variability of plasma protein levels, as compared to non-genetic factors. In addition, blood-cell fractions, including leukocytes, lymphocytes, monocytes, neutrophils, eosinophils, basophils, and platelets, had a larger contribution to inter-individual variability than age and sex and appeared as confounders of specific genetic associations. Finally, we identified new genetic associations with plasma protein levels of five monogenic Mendelian disease genes including two primary immunodeficiency genes (Ficolin-3 and FAS).

CONCLUSIONS

Our study identified novel genetic and non-genetic factors associated to plasma protein levels which may inform health status and disease management.

The 2022 Nucleic Acids Research database issue and the online molecular biology database collection

The 2022 Nucleic Acids Research Database Issue contains 185 papers, including 87 papers reporting on new databases and 85 updates from resources previously published in the Issue. Thirteen additional manuscripts provide updates on databases most recently published elsewhere. Seven new databases focus specifically on COVID-19 and SARS-CoV-2, including SCoV2-MD, the first of the Issue's Breakthrough Articles. Major nucleic acid databases reporting updates include MODOMICS, JASPAR and miRTarBase. The AlphaFold Protein Structure Database, described in the second Breakthrough Article, is the stand-out in the protein section, where the Human Proteoform Atlas and GproteinDb are other notable new arrivals. Updates from DisProt, FuzDB and ELM comprehensively cover disordered proteins. Under the metabolism and signalling section Reactome, ConsensusPathDB, HMDB and CAZy are major returning resources. In microbial and viral genomes taxonomy and systematics are well covered by LPSN, TYGS and GTDB. Genomics resources include Ensembl, Ensembl Genomes and UCSC Genome Browser. Major returning pharmacology resource names include the IUPHAR/BPS guide and the Therapeutic Target Database. New plant databases include PlantGSAD for gene lists and qPTMplants for post-translational modifications. The entire Database Issue is freely available online on the Nucleic Acids Research website (https://academic.oup.com/nar). Our latest update to the NAR online Molecular Biology Database Collection brings the total number of entries to 1645. Following last year's major cleanup, we have updated 317 entries, listing 89 new resources and trimming 80 discontinued URLs. The current release is available at http://www.oxfordjournals.org/nar/database/c/.