A Novel Approach to High-Quality Postmortem Tissue Procurement: The GTEx Project

Autophagopathies: from autophagy gene polymorphisms to precision medicine for human diseases

At a time when complex diseases affect globally 280 million people and claim 14 million lives every year, there is an urgent need to rapidly increase our knowledge into their underlying etiologies. Though critical in identifying the people at risk, the causal environmental factors (microbiome and/or pollutants) and the affected pathophysiological mechanisms are not well understood. Herein, we consider the variations of autophagy-related (ATG) genes at the heart of mechanisms of increased susceptibility to environmental stress. A comprehensive autophagy genomic resource is presented with 263 single nucleotide polymorphisms (SNPs) for 69 autophagy-related genes associated with 117 autoimmune, inflammatory, infectious, cardiovascular, neurological, respiratory, and endocrine diseases. We thus propose the term 'autophagopathies' to group together a class of complex human diseases the etiology of which lies in a genetic defect of the autophagy machinery, whether directly related or not to an abnormal flux in autophagy, LC3-associated phagocytosis, or any associated trafficking. The future of precision medicine for common diseases will lie in our ability to exploit these ATG SNP x environment relationships to develop new polygenetic risk scores, new management guidelines, and optimal therapies for afflicted patients.Abbreviations: ATG, autophagy-related; ALS-FTD, amyotrophic lateral sclerosis-frontotemporal dementia; ccRCC, clear cell renal cell carcinoma; CD, Crohn disease; COPD, chronic obstructive pulmonary disease; eQTL, expression quantitative trait loci; HCC, hepatocellular carcinoma; HNSCC, head and neck squamous cell carcinoma; GTEx, genotype-tissue expression; GWAS, genome-wide association studies; LAP, LC3-associated phagocytosis; LC3-II, phosphatidylethanolamine conjugated form of LC3; LD, linkage disequilibrium; LUAD, lung adenocarcinoma; MAF, minor allele frequency; MAP1LC3/LC3: microtubule associated protein 1 light chain 3; NSCLC, non-small cell lung cancer; OS, overall survival; PtdIns3K CIII, class III phosphatidylinositol 3 kinase; PtdIns3P, phosphatidylinositol-3-phosphate; SLE, systemic lupus erythematosus; SNPs, single-nucleotide polymorphisms; mQTL, methylation quantitative trait loci; ULK, unc-51 like autophagy activating kinase; UTRs, untranslated regions; WHO, World Health Organization.

The role of macrophage scavenger receptor 1 (MSR1) in inflammatory disorders and cancer

Macrophage scavenger receptor 1 (MSR1), also named CD204, holds key inflammatory roles in multiple pathophysiologic processes. Present primarily on the surface of various types of macrophage, this receptor variably affects processes such as atherosclerosis, innate and adaptive immunity, lung and liver disease, and more recently, cancer. As highlighted throughout this review, the role of MSR1 is often dichotomous, being either host protective or detrimental to the pathogenesis of disease. We will discuss the role of MSR1 in health and disease with a focus on the molecular mechanisms influencing MSR1 expression, how altered expression affects disease process and macrophage function, the limited cell signalling pathways discovered thus far, the emerging role of MSR1 in tumour associated macrophages as well as the therapeutic potential of targeting MSR1.

Shared Genetic Regulatory Networks Contribute to Neuropathic and Inflammatory Pain: Multi-Omics Systems Analysis

The mechanisms of chronic pain are complex, and genetic factors play an essential role in the development of chronic pain. Neuropathic pain (NP) and inflammatory pain (IP) are two primary components of chronic pain. Previous studies have uncovered some common biological processes in NP and IP. However, the shared genetic mechanisms remained poorly studied. We utilized multi-omics systematic analyses to investigate the shared genetic mechanisms of NP and IP. First, by integrating several genome-wide association studies (GWASs) with multi-omics data, we revealed the significant overlap of the gene co-expression modules in NP and IP. Further, we uncovered the shared biological pathways, including the previously reported mitochondrial electron transport and ATP metabolism, and stressed the role of genetic factors in chronic pain with neurodegenerative diseases. Second, we identified 24 conservative key drivers (KDs) contributing to NP and IP, containing two well-established pain genes, IL1B and OPRM1, and some novel potential pain genes, such as C5AR1 and SERPINE1. The subnetwork of those KDs highlighted the processes involving the immune system. Finally, gene expression analysis of the KDs in mouse models underlined two of the KDs, SLC6A15 and KCNQ5, with unidirectional regulatory functions in NP and IP. Our study provides strong evidence to support the current understanding of the shared genetic regulatory networks underlying NP and IP and potentially benefit the future common therapeutic avenues for chronic pain.

Composition of the infiltrating immune cells in the brain of healthy individuals: effect of aging

Immune cells infiltrating the central nervous system (CNS) are involved in the defense against invading microbes as well as in the pathogenesis of neuroinflammatory diseases. In these conditions, the presence of several types of immune and inflammatory cells have been demonstrated. However, some studies have also reported low amounts of immune cells that have been detected in the CNS of healthy individuals, but the cell types present have not been systematically analyzed. To do this, we now used brain samples from The Genotype- Tissue Expression (GTEx) project to analyze the relative abundance of 22 infiltrating leukocyte types using a digital cytometry tool (CIBERSORTx). To characterize cell proportions in different parts of the CNS, samples from 13 different anatomic brain regions were used. The data obtained demonstrated that several leukocyte types were present in the CNS. Six leukocyte types (CD4 memory resting T cells, M0 macrophages, plasma cells, CD8 T cells, CD4 memory activated T cells, and monocytes) were present with a proportion higher than 0.05, i.e. 5%. These six cell types were present in most brain regions with only insignificant variation. A consistent association with age was seen with monocytes, CD8 T cells, and follicular helper T cells. Taken together, these data show that several infiltrating immune cell types are present in the non-diseased CNS tissue and that the proportions of infiltrating cells are affected by age in a manner that is consistent with literature on immunosenecence and inflammaging.

Epidemiological evidence for associations between variants in CHRNA genes and risk of lung cancer and chronic obstructive pulmonary disease

Background

Genetic studies have previously reported that single-nucleotide polymorphisms (SNPs) in CHRNA genes (such as CHRNA3, CHRNA4, CHRNA5, or CHRNA3-CHRNA5-CHRNB4 clusters) are linked to the risk of neoplastic and non-neoplastic diseases. However, these conclusions were controversial and no systematic research synopsis has been available. We aimed to synthesize current knowledge of variants in the CHRNA genes on the risk of diseases.

Methods

We systematically searched for publications using PubMed, Medline, and Web of Science on or before 25 August 2021. A total of 1,818 publications were identified, of which 29 were deemed eligible for inclusion that could be used to perform meta-analysis based on at least three data sources to assess whether the morbidity associated with neoplastic and non-neoplastic diseases can be attributed to SNPs in CHRNA genes. To further evaluate the authenticity of cumulative evidence proving significant associations, the present study covered the Venice criteria and false-positive report probability tests. Through the Encyclopedia of DNA Elements (ENCODE) project, we created functional annotations for strong associations.

Results

Meta-analyses were done for nine genetic variants with two diseases {chronic obstructive pulmonary disease (COPD) and lung cancer (LC)}that had at least three data sources. Interestingly, eight polymorphisms were significantly related to changes in the susceptibility COPD and LC (p < 0.05). Of these, strong evidence was assigned to six variants (28 significant associations): CHRNA3 rs1051730, CHRNA3 rs6495309, and CHRNA5 rs16969968 with COPD risk, and CHRNA3 rs1051730, CHRNA3 rs578776, CHRNA3 rs6495309, CHRNA3 rs938682, CHRNA5 rs16969968, and CHRNA5 rs588765 with LC risk; moderate evidence was assigned to five SNPs (12 total associations) with LC or COPD risk. Data from ENCODE and other public databases showed that SNPs with strong evidence may be located in presumptive functional regions.

Conclusions

Our study summarized comprehensive evidence showing that common mutations in CHRNA genes are strongly related to LC and COPD risk. The study also elucidated the vital function of CHRNA genes in genetic predispositions to human diseases.

LRRC15 inhibits SARS-CoV-2 cellular entry in trans

Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2) infection is mediated by the entry receptor angiotensin-converting enzyme 2 (ACE2). Although attachment factors and coreceptors facilitating entry are extensively studied, cellular entry factors inhibiting viral entry are largely unknown. Using a surfaceome CRISPR activation screen, we identified human LRRC15 as an inhibitory attachment factor for SARS-CoV-2 entry. LRRC15 directly binds to the receptor-binding domain (RBD) of spike protein with a moderate affinity and inhibits spike-mediated entry. Analysis of human lung single-cell RNA sequencing dataset reveals that expression of LRRC15 is primarily detected in fibroblasts and particularly enriched in pathological fibroblasts in COVID-19 patients. ACE2 and LRRC15 are not coexpressed in the same cell types in the lung. Strikingly, expression of LRRC15 in ACE2-negative cells blocks spike-mediated viral entry in ACE2+ cell in trans, suggesting a protective role of LRRC15 in a physiological context. Therefore, LRRC15 represents an inhibitory attachment factor for SARS-CoV-2 that regulates viral entry in trans.

SNP rs9387478 at ROS1-DCBLD1 Locus is Significantly Associated with Lung Cancer Risk and Poor Survival in Indian Population

OBJECTIVE

Receptor tyrosine kinases (RTK) are relevant therapeutic targets in the treatment of lung cancer. Germline susceptibility variants that influence these RTKs may provide new insights into their regulation.  rs9387478 is located in the genomic interval between two RTK-genes ROS1/DCBLD1, of which ROS1 alterations are implicated in lung carcinogenesis and treatment response while the latter remains poorly understood.

MATERIALS AND METHODS

Venous blood was drawn from 100 control and 231 case subjects. Genotype was scored by restriction fragment length polymorphism (RFLP), PCR amplification followed by HindIII digestion. Logistic regression was applied to compare the association between variables. Survival curve was plotted to draw a correlation between the genotype and overall survival. Also, eQTL and chromatin state changes were analyzed and correlated with the survival of patients using available datasets.

RESULTS

In our population smoking correlated significantly with lung cancer [OR= 2.607] with the presence of the minor allele 'A' enhancing the nicotine dependence [CA (OR=3.23)]. Individuals with homozygous risk allele 'A' had a higher chance of developing lung cancer [OR=2.65] than individuals with CA/CC implying a recessive model of association. Patients with CC/CA genotype had better overall survival than patients with AA genotype [161 days/142 days vs 54 days, p=0.005]. The homozygous risk allele was significantly associated with increased DCBLD1 and ROS1 expression in lung cancer, with enriched active histone marks due to the polymorphism. Interestingly, increased DCBLD1 expression was associated with poor outcomes in lung cancer.

CONCLUSION

Overall, our study provides strong evidence that rs9387478 is significantly associated with both nicotine dependence and lung cancer in our North Indian cohort. The association of the SNP with prognostic genes, DCBLD1 and ROS1 make rs9387478 a promising prognostic marker in the North Indian population. The results obtained are significant, however, the study needs to be performed in a larger sample size.

A convergent malignant phenotype in B-cell acute lymphoblastic leukemia involving the splicing factor SRRM1

A significant proportion of infant B-cell acute lymphoblastic leukemia (B-ALL) patients remains with a dismal prognosis due to yet undetermined mechanisms. We performed a comprehensive multicohort analysis of gene expression, gene fusions, and RNA splicing alterations to uncover molecular signatures potentially linked to the observed poor outcome. We identified 87 fusions with significant allele frequency across patients and shared functional impacts, suggesting common mechanisms across fusions. We further identified a gene expression signature that predicts high risk independently of the gene fusion background and includes the upregulation of the splicing factor SRRM1. Experiments in B-ALL cell lines provided further evidence for the role of SRRM1 on cell survival, proliferation, and invasion. Supplementary analysis revealed that SRRM1 potentially modulates splicing events associated with poor outcomes through protein-protein interactions with other splicing factors. Our findings reveal a potential convergent mechanism of aberrant RNA processing that sustains a malignant phenotype independently of the underlying gene fusion and that could potentially complement current clinical strategies in infant B-ALL.

Genetic and molecular features of seizure-freedom following surgical resections for focal epilepsy: A pilot study

Objective

Seizure outcomes after brain surgery for drug-resistant epilepsy (DRE) are very heterogeneous and difficult to predict with models utilizing the current clinical, imaging, and electrophysiological variables. In this pilot study, we investigated whether genetic and molecular biomarkers (e.g., genomic, transcriptomic) can provide additional insight into differential response to surgery.

Methods

Post-operative seizure-outcomes were collected at last follow-up (>6 months) for 201 adult patients with DRE who underwent surgery between 2004 and 2020. Resected tissue was sent for miRNA sequencing (n = 132) and mRNA sequencing (n = 135). Following the selection of 10 genes (SCN1A, NBEA, PTEN, GABRA1, LGL1, DEPDC5, IL1A, ABCB1, C3, CALHM1), we investigated SNPs in those 10 genes from previously acquired exome sequencing data (n = 106). Logistic regression was performed to test for associations between individual features (mRNAs, miRNAs, and SNPs) and post-operative seizure-outcome with an exploratory FDR P < 0.25 as the threshold for significance. Post-operative time-to-seizure analyses were performed for each SNP using a Cox proportional hazards model.

Results

The majority of patients (83%) had temporal lobe epilepsy. Mean age at surgery was 38.3 years, and 56% were female. Three SNPs (rs10276036, rs11975994, rs1128503) in multi-drug resistance gene, ABCB1, were associated with post-operative seizure outcomes. Patients with alternate alleles in ABCB1 were more likely to be seizure-free at last follow-up (52–56% reduction in seizure recurrence; FDR P = 0.24). All three SNPs were in linkage disequilibrium and highly correlated with each other. Median post-operative time-to-seizure was 63 months for patients with 2 alternate alleles, 24–33 months with 1 alternate allele, and 10–11 months with 0 alternate alleles. These SNPs improved outcome prediction beyond MRI and sex alone. No independent miRNAs or mRNAs were significantly associated with seizure-outcome (P > 0.05). However, pathway analysis identified “cancer drug resistance by drug efflux” (mir-154 and mir-379) as enriched (P = 0.02), supporting the role of drug response genes in post-operative seizure recurrence.

Significance

ABCB1 may have a role in epileptogenesis and surgery outcomes independent of its drug efflux activity necessitating further investigation. SNPs in ABCB1 may serve as independent predictors of post-operative outcome.

ColocQuiaL: a QTL-GWAS colocalization pipeline

SUMMARY

Identifying genomic features responsible for genome-wide association study (GWAS) signals has proven to be a difficult challenge; many researchers have turned to colocalization analysis of GWAS signals with expression quantitative trait loci (eQTL) and splicing quantitative trait loci (sQTL) to connect GWAS signals to candidate causal genes. The ColocQuiaL pipeline provides a framework to perform these colocalization analyses at scale across the genome and returns summary files and locus visualization plots to allow for detailed review of the results. As an example, we used ColocQuiaL to perform colocalization between a recent type 2 diabetes GWAS and Genotype-Tissue Expression (GTEx) v8 single-tissue eQTL and sQTL data.

AVAILABILITY AND IMPLEMENTATION

ColocQuiaL is primarily written in R and is freely available on GitHub: https://github.com/bvoightlab/ColocQuiaL.

Association between SNPs in microRNAs and microRNAs-Machinery Genes with Susceptibility of Leprosy in the Amazon Population

Leprosy is a chronic neurodermatological disease caused by the bacillus Mycobacterium leprae. Recent studies show that SNPs in genes related to miRNAs have been associated with several diseases in different populations. This study aimed to evaluate the association of twenty-five SNPs in genes encoding miRNAs related to biological processes and immune response with susceptibility to leprosy and its polar forms paucibacillary and multibacillary in the Brazilian Amazon. A total of 114 leprosy patients and 71 household contacts were included in this study. Genotyping was performed using TaqMan Open Array Genotyping. Ancestry-informative markers were used to estimate individual proportions of case and control groups. The SNP rs2505901 (pre-miR938) was associated with protection against the development of paucibacillary leprosy, while the SNPs rs639174 (DROSHA), rs636832 (AGO1), and rs4143815 (miR570) were associated with protection against the development of multibacillary leprosy. In contrast, the SNPs rs10739971 (pri-let-7a1), rs12904 (miR200C), and rs2168518 (miR4513) are associated with the development of the paucibacillary leprosy. The rs10739971 (pri-let-7a1) polymorphism was associated with the development of leprosy, while rs2910164 (miR146A) and rs10035440 (DROSHA) was significantly associated with an increased risk of developing multibacillary leprosy.

Identification and validation of a ferroptosis-related gene signature for predicting survival in skin cutaneous melanoma

PURPOSE

Ferroptosis plays a crucial role in the initiation and progression of melanoma. This study developed a robust signature with ferroptosis-related genes (FRGs) and assessed the ability of this signature to predict OS in patients with skin cutaneous melanoma (SKCM).

METHODS

RNA-sequencing data and clinical information of melanoma patients were extracted from TCGA, GEO, and GTEx. Univariate, multivariate, and LASSO regression analyses were conducted to identify the gene signature. A 10 FRG signature was an independent and strong predictor of survival. The predictive performance was assessed using ROC curve. The functions of this gene signature were assessed by GO and KEGG analysis. The statuses of low-risk and high-risk groups according to the gene signature were compared by GSEA. In addition, we investigated the possible relationship of FRGs with immunotherapy efficacy.

RESULTS

A prognostic signature with 10 FRGs (CYBB, IFNG, FBXW7, ARNTL, PROM2, GPX2, JDP2, SLC7A5, TUBE1, and HAMP) was identified by Cox regression analysis. This signature had a higher prediction efficiency than clinicopathological features (AUC = 0.70). The enrichment analyses of DEGs indicated that ferroptosis-related immune pathways were largely enriched. Furthermore, GSEA showed that ferroptosis was associated with immunosuppression in the high-risk group. Finally, immune checkpoints such as PDCD-1 (PD-1), CTLA4, CD274 (PD-L1), and LAG3 were also differential expression in two risk groups.

CONCLUSIONS

The 10 FRGs signature were a strong predictor of OS in SKCM and could be used to predict therapeutic targets for melanoma.

Pan-Cancer Analysis Reveals the Relation between TRMT112 and Tumor Microenvironment

Dysregulated epigenetic modifications play a critical role in cancer development where TRMT112 is a member of the transfer RNA (tRNA) methyltransferase family. Till now, no studies have revealed the linkage between TRMT112 expression and diverse types of tumors. Based on TCGA data, we first probed into the relation between TRMT112 and prognosis and the potential role of TRMT112 in tumor microenvironment across 33 types of tumor. TRMT112 presented with increased expression in most cancers, which was significantly prognostic. Furthermore, TRMT112 was associated with tumor-associated fibroblasts in a variety of cancers. Additionally, a positive relationship was identified between TRMT112 expression and multiple tumor-related immune infiltrations, such as dendritic cells, CD8+ T cells, macrophages, CD4+ T cells, neutrophils, and B cells in lung adenocarcinoma and breast invasive carcinoma. In summary, our results suggest that TRMT112 might be a potential prognostic predictor of cancers and involved in regulating multiple cancer-related immune responses to some extent.

Defining novel causal SNPs and linked phenotypes at melanoma-associated loci

A number of genomic regions have been associated with melanoma risk through genome-wide association studies; however, the causal variants underlying the majority of these associations remain unknown. Here, we sequenced either the full locus or the functional regions including exons of 19 melanoma-associated loci in 1959 British melanoma cases and 737 controls. Variant filtering followed by Fisher's exact test analyses identified 66 variants associated with melanoma risk. Sequential conditional logistic regression identified the distinct haplotypes on which variants reside, and massively parallel reporter assays provided biological insights into how these variants influence gene function. We performed further analyses to link variants to melanoma risk phenotypes and assessed their association with melanoma-specific survival. Our analyses replicate previously known associations in the melanocortin 1 receptor (MC1R) and tyrosinase (TYR) loci, while identifying novel potentially causal variants at the MTAP/CDKN2A and CASP8 loci. These results improve our understanding of the architecture of melanoma risk and outcome.

Selecting Genetic Variants and Interactions Associated with Amyotrophic Lateral Sclerosis: A Group LASSO Approach

Amyotrophic lateral sclerosis (ALS) is a multi-system neurodegenerative disease that affects both upper and lower motor neurons, resulting from a combination of genetic, environmental, and lifestyle factors. Usually, the association between single-nucleotide polymorphisms (SNPs) and this disease is tested individually, which leads to the testing of multiple hypotheses. In addition, this classical approach does not support the detection of interaction-dependent SNPs. We applied a two-step procedure to select SNPs and pairwise interactions associated with ALS. SNP data from 276 ALS patients and 268 controls were analyzed by a two-step group LASSO in 2000 iterations. In the first step, we fitted a group LASSO model to a bootstrap sample and a random subset of predictors (25%) from the original data set aiming to screen for important SNPs and, in the second step, we fitted a hierarchical group LASSO model to evaluate pairwise interactions. An in silico analysis was performed on a set of variables, which were prioritized according to their bootstrap selection frequency. We identified seven SNPs (rs16984239, rs10459680, rs1436918, rs1037666, rs4552942, rs10773543, and rs2241493) and two pairwise interactions (rs16984239:rs2118657 and rs16984239:rs3172469) potentially involved in nervous system conservation and function. These results may contribute to the understanding of ALS pathogenesis, its diagnosis, and therapeutic strategy improvement.

Association of the SLC47A1 Gene Variant With Responses to Metformin Monotherapy in Drug-naive Patients With Type 2 Diabetes

CONTEXT

Although metformin is the first-line treatment for type 2 diabetes, the blood sugar-lowering effect of metformin varies among populations. SLC47A1 plays an important role in metformin pharmacokinetics and pharmacodynamics.

OBJECTIVE

We performed a systematic review and meta-analysis to investigate the association between SLC47A1 rs2289669 (G > A) and the metformin response in drug-naive patients with type 2 diabetes.

METHODS

Studies published until January 27, 2022, were retrieved from Cochrane CENTRAL, Embase, PubMed, and Web of Science. Two reviewers independently screened titles, abstracts, and full-text articles. Studies conducted in newly diagnosed or drug-naive patients with type 2 diabetes who received metformin monotherapy were included. A total of 6 studies involving 953 patients were included in this meta-analysis. We extracted the study characteristics and changes in glycated hemoglobin (HbA1c) levels before and after treatment according to the SLC47A1 rs2289669 genotype. Changes in HbA1c levels were analyzed using mean differences (MDs) and 95% CIs. SLC47A1 rs2289669 was associated with changes in HbA1c levels (A carrier vs GG; MD = -0.55; 95% CI, -0.91 to - 0.20; I² = 63%). The sensitivity analysis yielded similar results to the main analysis (MD range, -0.64 to -0.37). When comparing all 3 genotypes, there were significant differences in HbA1c level changes between AA vs GG and GA vs GG, but not in GA vs AA.

CONCLUSION

This meta-analysis showed that SLC47A1 rs2289669 is associated with the glycemic response to metformin in drug-naive patients with type 2 diabetes.

Network Embedding Across Multiple Tissues and Data Modalities Elucidates the Context of Host Factors Important for COVID-19 Infection

COVID-19 is a heterogeneous disease caused by SARS-CoV-2. Aside from infections of the lungs, the disease can spread throughout the body and damage many other tissues, leading to multiorgan failure in severe cases. The highly variable symptom severity is influenced by genetic predispositions and preexisting diseases which have not been investigated in a large-scale multimodal manner. We present a holistic analysis framework, setting previously reported COVID-19 genes in context with prepandemic data, such as gene expression patterns across multiple tissues, polygenetic predispositions, and patient diseases, which are putative comorbidities of COVID-19. First, we generate a multimodal network using the prior-based network inference method KiMONo. We then embed the network to generate a meaningful lower-dimensional representation of the data. The input data are obtained via the Genotype-Tissue Expression project (GTEx), containing expression data from a range of tissues with genomic and phenotypic information of over 900 patients and 50 tissues. The generated network consists of nodes, that is, genes and polygenic risk scores (PRS) for several diseases/phenotypes, as well as for COVID-19 severity and hospitalization, and links between them if they are statistically associated in a regularized linear model by feature selection. Applying network embedding on the generated multimodal network allows us to perform efficient network analysis by identifying nodes close by in a lower-dimensional space that correspond to entities which are statistically linked. By determining the similarity between COVID-19 genes and other nodes through embedding, we identify disease associations to tissues, like the brain and gut. We also find strong associations between COVID-19 genes and various diseases such as ischemic heart disease, cerebrovascular disease, and hypertension. Moreover, we find evidence linking PTPN6 to a range of comorbidities along with the genetic predisposition of COVID-19, suggesting that this kinase is a central player in severe cases of COVID-19. In conclusion, our holistic network inference coupled with network embedding of multimodal data enables the contextualization of COVID-19-associated genes with respect to tissues, disease states, and genetic risk factors. Such contextualization can be exploited to further elucidate the biological importance of known and novel genes for severity of the disease in patients.

NSD1 mutations deregulate transcription and DNA methylation of bivalent developmental genes in Sotos syndrome

Sotos syndrome (SS), the most common overgrowth with intellectual disability (OGID) disorder, is caused by inactivating germline mutations of NSD1, which encodes a histone H3 lysine 36 methyltransferase. To understand how NSD1 inactivation deregulates transcription and DNA methylation (DNAm), and to explore how these abnormalities affect human development, we profiled transcription and DNAm in SS patients and healthy control individuals. We identified a transcriptional signature that distinguishes individuals with SS from controls and was also deregulated in NSD1-mutated cancers. Most abnormally expressed genes displayed reduced expression in SS; these downregulated genes consisted mostly of bivalent genes and were enriched for regulators of development and neural synapse function. DNA hypomethylation was strongly enriched within promoters of transcriptionally deregulated genes: overexpressed genes displayed hypomethylation at their transcription start sites while underexpressed genes featured hypomethylation at polycomb binding sites within their promoter CpG island shores. SS patients featured accelerated molecular aging at the levels of both transcription and DNAm. Overall, these findings indicate that NSD1-deposited H3K36 methylation regulates transcription by directing promoter DNA methylation, partially by repressing polycomb repressive complex 2 (PRC2) activity. These findings could explain the phenotypic similarity of SS to OGID disorders that are caused by mutations in PRC2 complex-encoding genes.

LRLoop: a method to predict feedback loops in cell-cell communication

MOTIVATION

Intercellular communication (i.e. cell-cell communication) plays an essential role in multicellular organisms coordinating various biological processes. Previous studies discovered that feedback loops between two cell types are a widespread and vital signaling motif regulating development, regeneration and cancer progression. While many computational methods have been developed to predict cell-cell communication based on gene expression datasets, these methods often predict one-directional ligand-receptor interactions from sender to receiver cells and are not suitable to identify feedback loops.

RESULTS

Here, we describe ligand-receptor loop (LRLoop), a new method for analyzing cell-cell communication based on bi-directional ligand-receptor interactions, where two pairs of ligand-receptor interactions are identified that are responsive to each other and thereby form a closed feedback loop. We first assessed LRLoop using bulk datasets and found our method significantly reduces the false positive rate seen with existing methods. Furthermore, we developed a new strategy to assess the performance of these methods in single-cell datasets. We used the between-tissue interactions as an indicator of potential false-positive prediction and found that LRLoop produced a lower fraction of between-tissue interactions than traditional methods. Finally, we applied LRLoop to the single-cell datasets obtained from retinal development. We discovered many new bi-directional ligand-receptor interactions among individual cell types that potentially control proliferation, neurogenesis and/or cell fate specification.

AVAILABILITY AND IMPLEMENTATION

An R package is available at https://github.com/Pinlyu3/LRLoop. The source code can be found at figshare (https://doi.org/10.6084/m9.figshare.20126138.v1). The datasets can be found at figshare (https://doi.org/10.6084/m9.figshare.20126021.v1).

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

A locus at 19q13.31 significantly reduces the ApoE ε4 risk for Alzheimer's Disease in African Ancestry

African descent populations have a lower Alzheimer disease risk from ApoE ε4 compared to other populations. Ancestry analysis showed that the difference in risk between African and European populations lies in the ancestral genomic background surrounding the ApoE locus (local ancestry). Identifying the mechanism(s) of this protection could lead to greater insight into the etiology of Alzheimer disease and more personalized therapeutic intervention. Our objective is to follow up the local ancestry finding and identify the genetic variants that drive this risk difference and result in a lower risk for developing Alzheimer disease in African ancestry populations. We performed association analyses using a logistic regression model with the ApoE ε4 allele as an interaction term and adjusted for genome-wide ancestry, age, and sex. Discovery analysis included imputed SNP data of 1,850 Alzheimer disease and 4,331 cognitively intact African American individuals. We performed replication analyses on 63 whole genome sequenced Alzheimer disease and 648 cognitively intact Ibadan individuals. Additionally, we reproduced results using whole-genome sequencing of 273 Alzheimer disease and 275 cognitively intact admixed Puerto Rican individuals. A further comparison was done with SNP imputation from an additional 8,463 Alzheimer disease and 11,365 cognitively intact non-Hispanic White individuals. We identified a significant interaction between the ApoE ε4 allele and the SNP rs10423769_A allele, (β = -0.54,SE = 0.12,p-value = 7.50x10-6) in the discovery data set, and replicated this finding in Ibadan (β = -1.32,SE = 0.52,p-value = 1.15x10-2) and Puerto Rican (β = -1.27,SE = 0.64,p-value = 4.91x10-2) individuals. The non-Hispanic Whites analyses showed an interaction trending in the "protective" direction but failing to pass a 0.05 significance threshold (β = -1.51,SE = 0.84,p-value = 7.26x10-2). The presence of the rs10423769_A allele reduces the odds ratio for Alzheimer disease risk from 7.2 for ApoE ε4/ε4 carriers lacking the A allele to 2.1 for ApoE ε4/ε4 carriers with at least one A allele. This locus is located approximately 2 mB upstream of the ApoE locus, in a large cluster of pregnancy specific beta-1 glycoproteins on chromosome 19 and lies within a long noncoding RNA, ENSG00000282943. This study identified a new African-ancestry specific locus that reduces the risk effect of ApoE ε4 for developing Alzheimer disease. The mechanism of the interaction with ApoEε4 is not known but suggests a novel mechanism for reducing the risk for ε4 carriers opening the possibility for potential ancestry-specific therapeutic intervention.

Cumulative Evidence for Relationships Between Multiple Variants in the TERT and CLPTM1L Region and Risk of Cancer and Non-Cancer Disease

Background

Genetic studies previously reported that variants in TERT-CLPTM1L genes were related to susceptibility of cancer and non-cancer diseases. However, conclusions were not always concordant.

Methods

We performed meta-analyses to assess correlations between 23 variants within TERT-CLPTM1L region and susceptibility to 12 cancers and 1 non-cancer disease based on data in 109 papers (involving 139,510 cases and 208,530 controls). Two approaches (false-positive report probability test and Venice criteria) were adopted for assessing the cumulative evidence of significant associations. Current study evaluated the potential role of these variants based on data in Encyclopedia of DNA Elements (ENCODE) Project.

Results

Thirteen variants were statistically associated with susceptibility to 11 cancers and 1 non-cancer disease (p < 0.05). Besides, 12 variants with eight cancers and one non-cancer disease were rated as strong evidence (rs2736098, rs401681, and rs402710 in bladder cancer; rs2736100, rs2853691, and rs401681 in esophageal cancer; rs10069690 in gastric cancer; rs2736100 and rs2853676 in glioma; rs2242652, rs2736098, rs2736100, rs2853677, rs31489, rs401681, rs402710, rs465498, and rs4975616 in lung cancer; rs2736100 in idiopathic pulmonary fibrosis and myeloproliferative neoplasms; and rs401681 in pancreatic and skin cancer). According to data from ENCODE and other public databases, 12 variants with strong evidence might fall within putative functional regions.

Conclusions

This paper demonstrated that common variants of TERT-CLPTM1L genes were related to susceptibility to bladder, esophageal, gastric, lung, pancreatic, and skin cancer, as well as to glioma, myeloproliferative neoplasms, and idiopathic pulmonary fibrosis, and, besides, the crucial function of the TERT-CLPTM1L region in the genetic predisposition to human diseases is elucidated.

The topology of genome-scale metabolic reconstructions unravels independent modules and high network flexibility

The topology of metabolic networks is recognisably modular with modules weakly connected apart from sharing a pool of currency metabolites. Here, we defined modules as sets of reversible reactions isolated from the rest of metabolism by irreversible reactions except for the exchange of currency metabolites. Our approach identifies topologically independent modules under specific conditions associated with different metabolic functions. As case studies, the E.coli iJO1366 and Human Recon 2.2 genome-scale metabolic models were split in 103 and 321 modules respectively, displaying significant correlation patterns in expression data. Finally, we addressed a fundamental question about the metabolic flexibility conferred by reversible reactions: “Of all Directed Topologies (DTs) defined by fixing directions to all reversible reactions, how many are capable of carrying flux through all reactions?”. Enumeration of the DTs for iJO1366 model was performed using an efficient depth-first search algorithm, rejecting infeasible DTs based on mass-imbalanced and loopy flux patterns. We found the direction of 79% of reversible reactions must be defined before all directions in the network can be fixed, granting a high degree of flexibility.

Genome-scale metabolic reconstructions represent all biochemical reactions that an organism can accomplish. These reconstructions are complex and often difficult to study in great detail. A way to overcome this limitation is to focus on specific pathways or subsystems. We present a novel method to identify metabolic modules based on the network topology. The method relies on reaction directions and ignores currency metabolites, which artificially connect distant metabolic reactions. In this way, topologically independent modules are built, where inputs and outputs are controlled by irreversible reactions. The method is automatic and unbiased, and, the result is a set of condition specific modules with defined metabolic functions. As a proof-of-concept we generated biologically relevant modules for the E.coli and Human genome-scale metabolic reconstructions supported by transcriptomic data. Finally, we applied the novel approach to study the network flexibility conferred by reversible reactions. In the case of the E. coli model, we found that the direction of 79% of structurally reversible reactions (those not directionally constrained by surrounding irreversible reactions) must be fixed to determine all the reaction directions in the network. Therefore, reversible reactions operate practically independent of each other.

PolarMorphism enables discovery of shared genetic variants across multiple traits from GWAS summary statistics

MOTIVATION

Pleiotropic SNPs are associated with multiple traits. Such SNPs can help pinpoint biological processes with an effect on multiple traits or point to a shared etiology between traits. We present PolarMorphism, a new method for the identification of pleiotropic SNPs from genome-wide association studies (GWAS) summary statistics. PolarMorphism can be readily applied to more than two traits or whole trait domains. PolarMorphism makes use of the fact that trait-specific SNP effect sizes can be seen as Cartesian coordinates and can thus be converted to polar coordinates r (distance from the origin) and theta (angle with the Cartesian x-axis, in the case of two traits). r describes the overall effect of a SNP, while theta describes the extent to which a SNP is shared. r and theta are used to determine the significance of SNP sharedness, resulting in a P-value per SNP that can be used for further analysis.

RESULTS

We apply PolarMorphism to a large collection of publicly available GWAS summary statistics enabling the construction of a pleiotropy network that shows the extent to which traits share SNPs. We show how PolarMorphism can be used to gain insight into relationships between traits and trait domains and contrast it with genetic correlation. Furthermore, pathway analysis of the newly discovered pleiotropic SNPs demonstrates that analysis of more than two traits simultaneously yields more biologically relevant results than the combined results of pairwise analysis of the same traits. Finally, we show that PolarMorphism is more efficient and more powerful than previously published methods.

AVAILABILITY AND IMPLEMENTATION

code: https://github.com/UMCUGenetics/PolarMorphism, results: 10.5281/zenodo.5844193.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Knockout of Sorbin And SH3 Domain Containing 2 (Sorbs2) in Cardiomyocytes Leads to Dilated Cardiomyopathy in Mice

Background

Sorbin and SH3 domain containing 2 (Sorbs2) protein is a cytoskeletal adaptor with an emerging role in cardiac biology and disease; yet, its potential relevance to adult‐onset cardiomyopathies remains underexplored. Sorbs2 global knockout mice display lethal arrhythmogenic cardiomyopathy; however, the causative mechanisms remain unclear. Herein, we examine Sorbs2 dysregulation in heart failure, characterize novel Sorbs2 cardiomyocyte‐specific knockout mice (Sorbs2‐cKO), and explore associations between Sorbs2 genetic variations and human cardiovascular disease.

Methods and Results

Bioinformatic analyses show myocardial Sorbs2 mRNA is consistently upregulated in humans with adult‐onset cardiomyopathies and in heart failure models. We generated Sorbs2‐cKO mice and report that they develop progressive systolic dysfunction and enlarged cardiac chambers, and they die with congestive heart failure at about 1 year old. After 3 months, Sorbs2‐cKO mice begin to show atrial enlargement and P‐wave anomalies, without dysregulation of action potential–associated ion channel and gap junction protein expressions. After 6 months, Sorbs2‐cKO mice exhibit impaired contractility in dobutamine‐treated hearts and skinned myofibers, without dysregulation of contractile protein expressions. From our comprehensive survey of potential mechanisms, we found that within 4 months, Sorbs2‐cKO hearts have defective microtubule polymerization and compensatory upregulation of structural cytoskeletal and adapter proteins, suggesting that this early intracellular structural remodeling is responsible for contractile dysfunction. Finally, we identified genetic variants that associate with decreased Sorbs2 expression and human cardiac phenotypes, including conduction abnormalities, atrial enlargement, and dilated cardiomyopathy, consistent with Sorbs2‐cKO mice phenotypes.

Conclusions

Our studies show that Sorbs2 is essential for maintaining structural integrity in cardiomyocytes, likely through strengthening the interactions between microtubules and other cytoskeletal proteins at cross‐link sites.

Comprehensive Profiling Reveals Prognostic and Immunogenic Characteristics of Necroptosis in Soft Tissue Sarcomas

Soft tissue sarcomas (STSs) are heterogeneous malignancies derived from mesenchymal cells. Due to its rarity, heterogeneity, and limited overall response to chemotherapy, STSs represent a therapeutic challenge. Necroptosis is a novel therapeutic strategy for enhancing immunotherapy of cancer. Nevertheless, no research has explored the relationship between necroptosis-related genes (NRGs) and STSs. In this study, differentially expressed NRGs were identified using The Cancer Genome Atlas (TCGA) and The Cancer Genotype-Tissue Expression (GTEx) project. The expression levels of 34 NRGs were significantly different. Several key NRGs were validated using RT-qPCR and our own sequencing data. Patients with STSs were divided into two clusters using consensus cluster analysis, and significant differences were observed in their survival (p=0.002). We found the differentially expressed genes (DEGs) between the two clusters and carried out subsequent analysis. The necroptosis-related gene signatures with 10 key DEGs were identified with a risk score constructed. The prognosis of TCGA-SARC cohort with low necroptosis-related risk score was better (p<0.001). Meanwhile, the low-risk group had a significantly increased immune infiltration. Using the data of GSE17118 and another immunotherapy cohort as external validations, we observed significant survival differences between the two risk groups (p=0.019). The necroptosis-related risk score proved to be an independent prognostic factor, and a nomogram was further established and integrated with other clinical features. Notably, the necroptosis-related gene signature could also act as the prognostic indicator in other malignancies based on pan-cancer analysis. In summary, the study outlines NRGs in STSs and their potential role in prognosis and will be one of the important directions for future research.

Associations of GNAS and RGS Gene Polymorphisms with the Risk of Ritodrine-Induced Adverse Events in Korean Women with Preterm Labor: A Cohort Study

Ritodrine, a β2-adrenergic receptor agonist, is among most commonly prescribed tocolytic agents. This study aimed to evaluate the associations of single nucleotide polymorphisms in GNAS, RGS2, and RGS5 with the risk of ritodrine-induced adverse events (AEs) and develop a risk scoring system to identify high-risk patients. This is the prospective cohort study conducted at the Ewha Woman’s University Mokdong Hospital between January 2010 and October 2016. Pregnant women were included if they were treated with ritodrine for preterm labor with regular uterine contractions (at least 3 every 10 min) and cervical dilation. A total of 6, 3, and 5 single nucleotide polymorphisms (SNPs) of GNAS, RGS2, and RGS5 genes were genotyped and compared in patients with and without ritodrine-induced AEs. A total of 163 patients were included in this study. After adjusting confounders, GNAS rs3730168 (per-allele odds ratio (OR): 2.1; 95% confidence interval (95% CI): 1.0–4.3) and RGS2 rs1152746 (per-allele OR: 2.6, 95% CI: 1.1–6.5) were significantly associated with ritodrine-induced AEs. According to the constructed risk scoring models, patients with 0, 1, 2, 3, 4, and 5 points showed 0%, 13%, 19%, 31%, 46%, and 100% risks of AEs. This study suggested that GNAS and RGS2 polymorphisms could affect the risk of AEs in patients treated with ritodrine.

Clotting factor genes are associated with preeclampsia in high-altitude pregnant women in the Peruvian Andes

Preeclampsia is a multi-organ complication of pregnancy characterized by sudden hypertension and proteinuria that is among the leading causes of preterm delivery and maternal morbidity and mortality worldwide. The heterogeneity of preeclampsia poses a challenge for understanding its etiology and molecular basis. Intriguingly, risk for the condition increases in high-altitude regions such as the Peruvian Andes. To investigate the genetic basis of preeclampsia in a population living at high altitude, we characterized genome-wide variation in a cohort of preeclamptic and healthy Andean families (n = 883) from Puno, Peru, a city located above 3,800 meters of altitude. Our study collected genomic DNA and medical records from case-control trios and duos in local hospital settings. We generated genotype data for 439,314 SNPs, determined global ancestry patterns, and mapped associations between genetic variants and preeclampsia phenotypes. A transmission disequilibrium test (TDT) revealed variants near genes of biological importance for placental and blood vessel function. The top candidate region was found on chromosome 13 of the fetal genome and contains clotting factor genes PROZ, F7, and F10. These findings provide supporting evidence that common genetic variants within coagulation genes play an important role in preeclampsia. A selection scan revealed a potential adaptive signal around the ADAM12 locus on chromosome 10, implicated in pregnancy disorders. Our discovery of an association in a functional pathway relevant to pregnancy physiology in an understudied population of Native American origin demonstrates the increased power of family-based study design and underscores the importance of conducting genetic research in diverse populations.

Chromosome-specific retention of cancer-associated DNA hypermethylation following pharmacological inhibition of DNMT1

The DNA methylation status of the X-chromosome in cancer cells is often overlooked because of computational difficulties. Most of the CpG islands on the X-chromosome are mono-allelically methylated in normal female cells and only present as a single copy in male cells. We treated two colorectal cancer cell lines from a male (HCT116) and a female (RKO) with increasing doses of a DNA methyltransferase 1 (DNMT1)-specific inhibitor (GSK3685032/GSK5032) over several months to remove as much non-essential CpG methylation as possible. Profiling of the remaining DNA methylome revealed an unexpected, enriched retention of DNA methylation on the X-chromosome. Strikingly, the identified retained X-chromosome DNA methylation patterns accurately predicted de novo DNA hypermethylation in colon cancer patient methylomes in the TCGA COAD/READ cohort. These results suggest that a re-examination of tumors for X-linked DNA methylation changes may enable greater understanding of the importance of epigenetic silencing of cancer related genes.

A local human Vδ1 T cell population is associated with survival in nonsmall-cell lung cancer

Murine tissues harbor signature γδ T cell compartments with profound yet differential impacts on carcinogenesis. Conversely, human tissue-resident γδ cells are less well defined. In the present study, we show that human lung tissues harbor a resident Vδ1 γδ T cell population. Moreover, we demonstrate that Vδ1 T cells with resident memory and effector memory phenotypes were enriched in lung tumors compared with nontumor lung tissues. Intratumoral Vδ1 T cells possessed stem-like features and were skewed toward cytolysis and helper T cell type 1 function, akin to intratumoral natural killer and CD8+ T cells considered beneficial to the patient. Indeed, ongoing remission post-surgery was significantly associated with the numbers of CD45RA-CD27- effector memory Vδ1 T cells in tumors and, most strikingly, with the numbers of CD103+ tissue-resident Vδ1 T cells in nonmalignant lung tissues. Our findings offer basic insights into human body surface immunology that collectively support integrating Vδ1 T cell biology into immunotherapeutic strategies for nonsmall cell lung cancer.

RNA analysis of the GALNS transcript reveals novel pathogenic mechanisms associated with Morquio syndrome A

Morquio syndrome A (Mucopolysaccharidosis IVA, MPS IVA) is an autosomal recessive lysosomal storage disorder caused by deficiency of N-acetyl-galactosamine-6-sulfatase (GALNS) which catabolizes the glycosaminoglycans (GAG), keratan sulfate and chondroitin-6-sulfate. Homozygous or compound heterozygous pathogenic variants in the GALNS result in the deficiency of the enzyme and consequent GAG accumulations. DNA sequence and copy number analysis of the GALNS coding region fails to identify biallelic causative pathogenic variants in up to 15% of patients with Morquio syndrome A. RNA transcript analysis was performed to identify pathogenic alterations in two unrelated families with Morquio syndrome A in whom a single heterozygous or no pathogenic alteration was detected by standard analysis of the GALNS gene. RNA sequencing and quantitative expression analysis identified the overabundance of an aberrant GALNS transcript isoform and a reduction of the clinically relevant isoform (NM_000512.4) in the Morquio syndrome A patients from both families. The aberrant isoform (ENST00000568613.1) was produced by alternative splicing and contained intronic sequence that was likely a cryptic exon predicted to result in a reading frame shift and generation of a premature termination codon. These findings indicated that the aberrant splicing is likely the novel molecular defect in our patients. RNA transcript analysis could be useful to identify pathogenic alterations and increase the yield of molecular diagnosis in patients with Morquio syndrome A whose genetic variants are not found by standard sequencing or gene dosage analysis.

Serotonin G Protein-Coupled Receptor-Based Biosensing Modalities in Yeast

Serotonin is a key neurotransmitter involved in numerous physiological processes and serves as an important precursor for manufacturing bioactive indoleamines and alkaloids used in the treatment of human pathologies. In humans, serotonin sensing and signaling can occur by 12 G protein-coupled receptors (GPCRs) coupled to Gα proteins. In yeast, human serotonin GPCRs coupled to Gα proteins have previously been shown to function as whole-cell biosensors of serotonin. However, systematic characterization of serotonin biosensing modalities between variant serotonin GPCRs and application thereof for high-resolution serotonin quantification is still awaiting. To systematically assess GPCR signaling in response to serotonin, we characterized reporter gene expression at two different pHs of a 144-sized library encoding all 12 human serotonin GPCRs in combination with 12 different Gα proteins engineered in yeast. From this screen, we observed changes in the biosensor sensitivities of >4 orders of magnitude. Furthermore, adopting optimal biosensing designs and pH conditions enabled high-resolution high-performance liquid chromatography-validated sensing of serotonin produced in yeast. Lastly, we used the yeast platform to characterize 19 serotonin GPCR polymorphisms found in human populations. While major differences in signaling were observed among the individual polymorphisms when studied in yeast, a cross-comparison of selected variants in mammalian cells showed both similar and disparate results. Taken together, our study highlights serotonin biosensing modalities of relevance to both biotechnological and potential human health applications.

The Less Well-Known Little Brothers: The SLC9B/NHA Sodium Proton Exchanger Subfamily—Structure, Function, Regulation and Potential Drug-Target Approaches

The SLC9 gene family encodes Na+/H+ exchangers (NHEs), a group of membrane transport proteins critically involved in the regulation of cytoplasmic and organellar pH, cell volume, as well as systemic acid-base and volume homeostasis. NHEs of the SLC9A subfamily (NHE 1–9) are well-known for their roles in human physiology and disease. Much less is known about the two members of the SLC9B subfamily, NHA1 and NHA2, which share higher similarity to prokaryotic NHEs than the SLC9A paralogs. NHA2 (also known as SLC9B2) is ubiquitously expressed and has recently been shown to participate in renal blood pressure and electrolyte regulation, insulin secretion and systemic glucose homeostasis. In addition, NHA2 has been proposed to contribute to the pathogenesis of polycystic kidney disease, the most common inherited kidney disease in humans. NHA1 (also known as SLC9B1) is mainly expressed in testis and is important for sperm motility and thus male fertility, but has not been associated with human disease thus far. In this review, we present a summary of the structure, function and regulation of expression of the SLC9B subfamily members, focusing primarily on the better-studied SLC9B paralog, NHA2. Furthermore, we will review the potential of the SLC9B subfamily as drug targets.

Single-nucleus cross-tissue molecular reference maps toward understanding disease gene function

Understanding gene function and regulation in homeostasis and disease requires knowledge of the cellular and tissue contexts in which genes are expressed. Here, we applied four single-nucleus RNA sequencing methods to eight diverse, archived, frozen tissue types from 16 donors and 25 samples, generating a cross-tissue atlas of 209,126 nuclei profiles, which we integrated across tissues, donors, and laboratory methods with a conditional variational autoencoder. Using the resulting cross-tissue atlas, we highlight shared and tissue-specific features of tissue-resident cell populations; identify cell types that might contribute to neuromuscular, metabolic, and immune components of monogenic diseases and the biological processes involved in their pathology; and determine cell types and gene modules that might underlie disease mechanisms for complex traits analyzed by genome-wide association studies.

Regulation of A-to-I RNA editing and stop codon recoding to control selenoprotein expression during skeletal myogenesis

Selenoprotein N (SELENON), a selenocysteine (Sec)-containing protein with high reductive activity, maintains redox homeostasis, thereby contributing to skeletal muscle differentiation and function. Loss-of-function mutations in SELENON cause severe neuromuscular disorders. In the early-to-middle stage of myoblast differentiation, SELENON maintains redox homeostasis and modulates endoplasmic reticulum (ER) Ca2+ concentration, resulting in a gradual reduction from the middle-to-late stages due to unknown mechanisms. The present study describes post-transcriptional mechanisms that regulate SELENON expression during myoblast differentiation. Part of an Alu element in the second intron of SELENON pre-mRNA is frequently exonized during splicing, resulting in an aberrant mRNA that is degraded by nonsense-mediated mRNA decay (NMD). In the middle stage of myoblast differentiation, ADAR1-mediated A-to-I RNA editing occurs in the U1 snRNA binding site at 5' splice site, preventing Alu exonization and producing mature mRNA. In the middle-to-late stage of myoblast differentiation, the level of Sec-charged tRNASec decreases due to downregulation of essential recoding factors for Sec insertion, thereby generating a premature termination codon in SELENON mRNA, which is targeted by NMD.

Elucidating Human Milk Oligosaccharide biosynthetic genes through network-based multi-omics integration

Human Milk Oligosaccharides (HMOs) are abundant carbohydrates fundamental to infant health and development. Although these oligosaccharides were discovered more than half a century ago, their biosynthesis in the mammary gland remains largely uncharacterized. Here, we use a systems biology framework that integrates glycan and RNA expression data to construct an HMO biosynthetic network and predict glycosyltransferases involved. To accomplish this, we construct models describing the most likely pathways for the synthesis of the oligosaccharides accounting for >95% of the HMO content in human milk. Through our models, we propose candidate genes for elongation, branching, fucosylation, and sialylation of HMOs. Our model aggregation approach recovers 2 of 2 previously known gene-enzyme relations and 2 of 3 empirically confirmed gene-enzyme relations. The top genes we propose for the remaining 5 linkage reactions are consistent with previously published literature. These results provide the molecular basis of HMO biosynthesis necessary to guide progress in HMO research and application with the goal of understanding and improving infant health and development.

Human milk oligosaccharides are fundamental to infant health. Here the authors deploy a multi-omics systems biology approach to elucidate their biosynthetic network, including the associated enzymes and likely structures of ambiguous oligosaccharides.

Ceramides and Acute Kidney Injury

Altered lipid metabolism is a characteristic feature and potential driving factor of acute kidney injury (AKI). Of the lipids that accumulate in injured renal tissues, ceramides are potent regulators of metabolism and cell fate. Up-regulation of ceramide synthesis is a common feature shared across several AKI etiologies in vitro and in vivo. Furthermore, ceramide accumulation is an early event in the natural history of AKI that precedes cell death and organ dysfunction. Emerging evidence suggests that inhibition of ceramide accumulation may improve renal outcomes in several models of AKI. This review examines the landscape of ceramide metabolism and regulation in the healthy and injured kidney. Furthermore, we discuss the body of literature regarding ceramides as therapeutic targets for AKI and consider potential mechanisms by which ceramides drive kidney pathogenesis.

Comprehensive Molecular Analyses of a TNF Family-Based Gene Signature as a Potentially Novel Prognostic Biomarker for Cervical Cancer

Background

Increasing evidence suggests that tumour necrosis factor (TNF) family genes play important roles in cervical cancer (CC). However, whether TNF family genes can be used as prognostic biomarkers of CC and the molecular mechanisms of TNF family genes remain unclear.

Methods

A total of 306 CC and 13 normal samples were obtained from The Cancer Genome Atlas (TCGA) and Genotype-Tissue Expression (GTEx) databases. We identified differentially expressed TNF family genes between CC and normal samples and subjected them to univariate Cox regression analysis for selecting prognostic TNF family genes. Least absolute shrinkage and selection operator (LASSO) regression and multivariate Cox regression analyses were performed to screen genes to establish a TNF family gene signature. Gene set enrichment analysis (GSEA) was performed to investigate the biological functions of the TNF family gene signature. Finally, methylation and copy number variation data of CC were used to analyse the potential molecular mechanisms of TNF family genes.

Results

A total of 26 differentially expressed TNF family genes were identified between the CC and normal samples. Next, a TNF family gene signature, including CD27, EDA, TNF, TNFRSF12A, TNFRSF13C, and TNFRSF9 was constructed based on univariate Cox, LASSO, and multivariate Cox regression analyses. The TNF family gene signature was related to age, pathological stages M and N, and could predict patient survival independently of clinical factors. Moreover, KEGG enrichment analysis suggested that the TNF family gene signature was mainly involved in the TGF-β signaling pathway, and the TNF family gene signature could affect the immunotherapy response. Finally, we confirmed that the mRNA expressions of CD27, TNF, TNFRSF12A, TNFRSF13C, and TNFRSF9 were upregulated in CC, while that of EDA was downregulated. The mRNA expressions of CD27, EDA, TNF, TNFRSF12A, TNFRSF13C, and TNFRSF9 might be influenced by gene methylation and copy number variation.

Conclusion

Our study is the first to demonstrate that CD27, EDA, TNF, TNFRSF12A, TNFRSF13C, and TNFRSF9 might be used as prognostic biomarkers of CC and are associated with the immunotherapy response of CC.

Alternative splicing regulation of membrane trafficking genes during myogenesis

Alternative splicing transitions occur during organ development, and, in numerous diseases, splicing programs revert to fetal isoform expression. We previously found that extensive splicing changes occur during postnatal mouse heart development in genes encoding proteins involved in vesicle-mediated trafficking. However, the regulatory mechanisms of this splicing-trafficking network are unknown. Here, we found that membrane trafficking genes are alternatively spliced in a tissue-specific manner, with striated muscles exhibiting the highest levels of alternative exon inclusion. Treatment of differentiated muscle cells with chromatin-modifying drugs altered exon inclusion in muscle cells. Examination of several RNA-binding proteins revealed that the poly-pyrimidine tract binding protein 1 (PTBP1) and quaking regulate splicing of trafficking genes during myogenesis, and that removal of PTBP1 motifs prevented PTBP1 from binding its RNA target. These findings enhance our understanding of developmental splicing regulation of membrane trafficking proteins which might have implications for muscle disease pathogenesis.

The biological behavior of tRNA-derived fragment tRF-Leu-AAG in pancreatic cancer cells

Pancreatic cancer (PC) is a life-threatening cancer with increasing incidence in developed countries. Reports indicate that tRNA-derived fragments (tRFs) are possible therapeutic targets and biomarkers for cancer treatment. Nonetheless, the effect of tRF-Leu-AAG on PC is unclear. This study aims to explore the role of tRF-Leu-AAG and upstream frameshift mutant 1 (UPF1) in the development of PC and its potential underlying mechanisms. High-throughput second-generation sequencing techniques were used to detect the expression of tRFs in cancerous and adjacent normal tissues from PC patients. The role of tRF-Leu-AAG proliferation in PC cells was investigated via the Cell Counting Kit-8 (CCK8) assay. The effect of tRF-Leu-AAG on the invasion and migration ability of PC cells was also determined by the transwell assay. Thereafter, the downstream target genes of tRF-Leu-AAG were comprehensively predicted using bioinformatics analysis databases. We also used the Dual-Luciferase Reporter assay to assess the nexus between tRF-Leu-AAG and UPF1. Eventually, Western Blot was used to validate the expression of UPF1 in PC cells. A total of 33 tRF expressions significantly varied from PC patients. RT-qPCR confirmed that the expression of tRF-Leu-AAG was observably up-regulated in PC cells as compared to the control cells. Importantly, knockdown of tRF-Leu-AAG observably inhibited cell proliferation, migration, and invasion. Furthermore, according to the predicted frameshift database results, the UPF1 acted as downstream target genes for tRF-Leu-AAG and significantly down-regulated UPF1 expression.

Enhancing Discovery of Genetic Variants for Posttraumatic Stress Disorder Through Integration of Quantitative Phenotypes and Trauma Exposure Information

BACKGROUND

Posttraumatic stress disorder (PTSD) is heritable and a potential consequence of exposure to traumatic stress. Evidence suggests that a quantitative approach to PTSD phenotype measurement and incorporation of lifetime trauma exposure (LTE) information could enhance the discovery power of PTSD genome-wide association studies (GWASs).

METHODS

A GWAS on PTSD symptoms was performed in 51 cohorts followed by a fixed-effects meta-analysis (N = 182,199 European ancestry participants). A GWAS of LTE burden was performed in the UK Biobank cohort (N = 132,988). Genetic correlations were evaluated with linkage disequilibrium score regression. Multivariate analysis was performed using Multi-Trait Analysis of GWAS. Functional mapping and annotation of leading loci was performed with FUMA. Replication was evaluated using the Million Veteran Program GWAS of PTSD total symptoms.

RESULTS

GWASs of PTSD symptoms and LTE burden identified 5 and 6 independent genome-wide significant loci, respectively. There was a 72% genetic correlation between PTSD and LTE. PTSD and LTE showed largely similar patterns of genetic correlation with other traits, albeit with some distinctions. Adjusting PTSD for LTE reduced PTSD heritability by 31%. Multivariate analysis of PTSD and LTE increased the effective sample size of the PTSD GWAS by 20% and identified 4 additional loci. Four of these 9 PTSD loci were independently replicated in the Million Veteran Program.

CONCLUSIONS

Through using a quantitative trait measure of PTSD, we identified novel risk loci not previously identified using prior case-control analyses. PTSD and LTE have a high genetic overlap that can be leveraged to increase discovery power through multivariate methods.

Associations of carotid intima media thickness with gene expression in whole blood and genetically predicted gene expression across 48 tissues

Carotid intima media thickness (cIMT) is a biomarker of subclinical atherosclerosis and a predictor of future cardiovascular events. Identifying associations between gene expression levels and cIMT may provide insight to atherosclerosis etiology. Here, we use two approaches to identify associations between mRNA levels and cIMT: differential gene expression analysis in whole blood and S-PrediXcan. We used microarrays to measure genome-wide whole blood mRNA levels of 5647 European individuals from four studies. We examined the association of mRNA levels with cIMT adjusted for various potential confounders. Significant associations were tested for replication in three studies totaling 3943 participants. Next, we applied S-PrediXcan to summary statistics from a cIMT genome-wide association study (GWAS) of 71 128 individuals to estimate the association between genetically determined mRNA levels and cIMT and replicated these analyses using S-PrediXcan on an independent GWAS on cIMT that included 22 179 individuals from the UK Biobank. mRNA levels of TNFAIP3, CEBPD and METRNL were inversely associated with cIMT, but these associations were not significant in the replication analysis. S-PrediXcan identified associations between cIMT and genetically determined mRNA levels for 36 genes, of which six were significant in the replication analysis, including TLN2, which had not been previously reported for cIMT. There was weak correlation between our results using differential gene expression analysis and S-PrediXcan. Differential expression analysis and S-PrediXcan represent complementary approaches for the discovery of associations between phenotypes and gene expression. Using these approaches, we prioritize TNFAIP3, CEBPD, METRNL and TLN2 as new candidate genes whose differential expression might modulate cIMT.

Shared Genetic Basis and Causal Relationship Between Television Watching, Breakfast Skipping and Type 2 Diabetes: Evidence From a Comprehensive Genetic Analysis

Background

Epidemiological investigations have established unhealthy lifestyles, such as excessive leisurely sedentary behavior (especially TV/television watching) and breakfast skipping, increase the risk of type 2 diabetes (T2D), but the causal relationship is unclear. We aimed to understand how single nucleotide variants contribute to the co-occurrence of unhealthy lifestyles and T2D, thereby providing meaningful insights into disease mechanisms.

Methods

Combining summary statistics from genome-wide association studies (GWAS) on TV watching (N = 422218), breakfast skipping (N = 193860) and T2D (N = 159208) in European pedigrees, we conducted comprehensive pairwise genetic analysis, including high-definition likelihood (HDL-method), cross-phenotype association studies (CPASSOC), GWAS-eQTL colocalization analysis and transcriptome-wide association studies (TWAS), to understand the genetic overlap between them. We also performed bidirectional two-sample Mendelian randomization (MR) analysis for causal inference using genetic instrumental variables, and two-step MR mediation analysis was used to assess any effects explained by body mass index, lipid traits and glycemic traits.

Results

HDL-method showed that T2D shared a strong genetic correlation with TV watching (rg = 0.26; P = 1.63×10-29) and skipping breakfast (rg = 0.15; P =2.02×10-6). CPASSOC identifies eight independent SNPs shared between T2D and TV watching, including one novel shared locus. TWAS and CPASSOC showed that shared genes were enriched in lung, esophageal, adipose, and thyroid tissues and highlighted potential shared regulatory pathways for lipoprotein metabolism, pancreatic β-cell function, cellular senescence and multi-mediator factors. MR showed TV watching had a causal effect on T2D (βIVW = 0.629, PIVW = 1.80×10-10), but no significant results were observed between breakfast skipping and T2D. Mediation analysis provided evidence that body mass index, fasting glucose, hemoglobin A1c and high-density lipoprotein are potential factors that mediate the causal relationship between TV and T2D.

Conclusions

Our findings provide strong evidence of shared genetics and causation between TV watching and T2D and facilitate our identification of common genetic architectures shared between them.

The Power of Single-Cell RNA Sequencing in eQTL Discovery

Genome-wide association studies have successfully mapped thousands of loci associated with complex traits. During the last decade, functional genomics approaches combining genotype information with bulk RNA-sequencing data have identified genes regulated by GWAS loci through expression quantitative trait locus (eQTL) analysis. Single-cell RNA-Sequencing (scRNA-Seq) technologies have created new exciting opportunities for spatiotemporal assessment of changes in gene expression at the single-cell level in complex and inherited conditions. A growing number of studies have demonstrated the power of scRNA-Seq in eQTL mapping across different cell types, developmental stages and stimuli that could be obscured when using bulk RNA-Seq methods. In this review, we outline the methodological principles, advantages, limitations and the future experimental and analytical considerations of single-cell eQTL studies. We look forward to the explosion of single-cell eQTL studies applied to large-scale population genetics to take us one step closer to understanding the molecular mechanisms of disease.

Herpesviruses and their genetic diversity in the blood virome of healthy individuals: effect of aging

BACKGROUND

As we age, the functioning of the human immune system declines. The results of this are increases in morbidity and mortality associated with infectious diseases, cancer, cardiovascular disease, and neurodegenerative disease in elderly individuals, as well as a weakened vaccination response. The aging of the immune system is thought to affect and be affected by the human virome, the collection of all viruses present in an individual. Persistent viral infections, such as those caused by certain herpesviruses, can be present in an individual for long periods of time without any overt pathology, yet are associated with disease in states of compromised immune function. To better understand the effects on human health of such persistent viral infections, we must first understand how the human virome changes with age. We have now analyzed the composition of the whole blood virome of 317 individuals, 21-70 years old, using a metatranscriptomic approach. Use of RNA sequencing data allows for the unbiased detection of RNA viruses and active DNA viruses.

RESULTS

The data obtained showed that Epstein-Barr virus (EBV) was the most frequently expressed virus, with other detected viruses being herpes simplex virus 1, human cytomegalovirus, torque teno viruses, and papillomaviruses. Of the 317 studied blood samples, 68 (21%) had EBV expression, whereas the other detected viruses were only detected in at most 6 samples (2%). We therefore focused on EBV in our further analyses. Frequency of EBV detection, relative EBV RNA abundance and the genetic diversity of EBV was not significantly different between age groups (21-59 and 60-70 years old). No significant correlation was seen between EBV RNA abundance and age. Deconvolution analysis revealed a significant difference in proportions of activated dendritic cells, macrophages M1, and activated mast cells between EBV expression positive and negative individuals.

CONCLUSIONS

As it is likely that the EBV RNA quantified in this work is derived from reactivation of the latent EBV virus, these data suggest that age does not affect the rate of reactivation nor the genetic landscape of EBV. These findings offer new insight on the genetic diversity of a persistent EBV infection in the long-term.

Transcriptional, epigenetic and metabolic signatures in cardiometabolic syndrome defined by extreme phenotypes

BACKGROUND

This work is aimed at improving the understanding of cardiometabolic syndrome pathophysiology and its relationship with thrombosis by generating a multi-omic disease signature.

METHODS/RESULTS

We combined classic plasma biochemistry and plasma biomarkers with the transcriptional and epigenetic characterisation of cell types involved in thrombosis, obtained from two extreme phenotype groups (morbidly obese and lipodystrophy) and lean individuals to identify the molecular mechanisms at play, highlighting patterns of abnormal activation in innate immune phagocytic cells. Our analyses showed that extreme phenotype groups could be distinguished from lean individuals, and from each other, across all data layers. The characterisation of the same obese group, 6 months after bariatric surgery, revealed the loss of the abnormal activation of innate immune cells previously observed. However, rather than reverting to the gene expression landscape of lean individuals, this occurred via the establishment of novel gene expression landscapes. NETosis and its control mechanisms emerge amongst the pathways that show an improvement after surgical intervention.

CONCLUSIONS

We showed that the morbidly obese and lipodystrophy groups, despite some differences, shared a common cardiometabolic syndrome signature. We also showed that this could be used to discriminate, amongst the normal population, those individuals with a higher likelihood of presenting with the disease, even when not displaying the classic features.

Genetic variants associated mRNA stability in lung

Background

Expression quantitative trait loci (eQTLs) analyses have been widely used to identify genetic variants associated with gene expression levels to understand what molecular mechanisms underlie genetic traits. The resultant eQTLs might affect the expression of associated genes through transcriptional or post-transcriptional regulation. In this study, we attempt to distinguish these two types of regulation by identifying genetic variants associated with mRNA stability of genes (stQTLs).

Results

Here, we presented a computational framework that takes advantage of recently developed methods to infer the mRNA stability of genes based on RNA-seq data and performed association analysis to identify stQTLs. Using the Genotype-Tissue Expression (GTEx) lung RNA-Seq data, we identified a total of 142,801 stQTLs for 3942 genes and 186,132 eQTLs for 4751 genes from 15,122,700 genetic variants for 13,476 genes on the autosomes, respectively. Interestingly, our results indicated that stQTLs were enriched in the CDS and 3’UTR regions, while eQTLs are enriched in the CDS, 3’UTR, 5’UTR, and upstream regions. We also found that stQTLs are more likely than eQTLs to overlap with RNA binding protein (RBP) and microRNA (miRNA) binding sites. Our analyses demonstrate that simultaneous identification of stQTLs and eQTLs can provide more mechanistic insight on the association between genetic variants and gene expression levels.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08405-y.

Irx5 and transient outward K+ currents contribute to transmural contractile heterogeneities in the mouse ventricle

Previous studies have established that fast transmural gradients of transient outward K⁺ current (I_to,f) correlate with regional differences in action potential (AP) profile and excitation-contraction coupling (ECC) with high I_to,f expression in the epimyocardium (EPI) being associated with short APs and low contractility and vice versa. Herein, we investigated the effects of disrupted I_to,f gradient on contractile properties using mouse models of Irx5 knockout (Irx5-KO) for selective I_to,f elevation in the endomyocardium (ENDO) of the left ventricle (LV) and Kcnd2 ablation (K_V4.2-KO) for selective I_to,freduction in the EPI. Irx5-KO mice exhibited decreased global LV contractility in association with reductions in cell shortening and Ca²⁺ transient amplitudes in isolated ENDO but not EPI cardiomyocytes. Moreover, transcriptional profiling revealed that the primary effect of Irx5 ablation on ECC-related genes was to increase I_to,f gene expression (i.e. Kcnd2 and Kcnip2) in the ENDO, but not the EPI. Indeed, K_V4.2-KO mice showed selective increases in cell shortening and Ca²⁺ transients in isolated EPI cardiomyocytes, leading to enhanced ventricular contractility and mice lacking both Irx5 and Kcnd2 displayed elevated ventricular contractility comparable to K_V4.2-KO mice. Our findings demonstrate that the transmural electromechanical heterogeneities in the healthy ventricles depend on the Irx5-dependent I_to,f gradients. These observations provide a useful framework for assessing the molecular mechanisms underlying the alterations in contractile heterogeneity seen in the diseased heart.

Meta-GWAS Reveals Novel Genetic Variants Associated with Urinary Excretion of Uromodulin

BACKGROUND

Uromodulin, the most abundant protein excreted in normal urine, plays major roles in kidney physiology and disease. The mechanisms regulating the urinary excretion of uromodulin remain essentially unknown.

METHODS

We conducted a meta-analysis of genome-wide association studies for raw (uUMOD) and indexed to creatinine (uUCR) urinary levels of uromodulin in 29,315 individuals of European ancestry from 13 cohorts. We tested the distribution of candidate genes in kidney segments and investigated the effects of keratin-40 (KRT40) on uromodulin processing.

RESULTS

Two genome-wide significant signals were identified for uUMOD: a novel locus (P 1.24E-08) over the KRT40 gene coding for KRT40, a type 1 keratin expressed in the kidney, and the UMOD-PDILT locus (P 2.17E-88), with two independent sets of single nucleotide polymorphisms spread over UMOD and PDILT. Two genome-wide significant signals for uUCR were identified at the UMOD-PDILT locus and at the novel WDR72 locus previously associated with kidney function. The effect sizes for rs8067385, the index single nucleotide polymorphism in the KRT40 locus, were similar for both uUMOD and uUCR. KRT40 colocalized with uromodulin and modulating its expression in thick ascending limb (TAL) cells affected uromodulin processing and excretion.

CONCLUSIONS

Common variants in KRT40, WDR72, UMOD, and PDILT associate with the levels of uromodulin in urine. The expression of KRT40 affects uromodulin processing in TAL cells. These results, although limited by lack of replication, provide insights into the biology of uromodulin, the role of keratins in the kidney, and the influence of the UMOD-PDILT locus on kidney function.

TDP-43 represses cryptic exon inclusion in the FTD-ALS gene UNC13A

A hallmark pathological feature of the neurodegenerative diseases amyotrophic lateral sclerosis (ALS) and frontotemporal dementia (FTD) is the depletion of RNA-binding protein TDP-43 from the nucleus of neurons in the brain and spinal cord1. A major function of TDP-43 is as a repressor of cryptic exon inclusion during RNA splicing2-4. Single nucleotide polymorphisms in UNC13A are among the strongest hits associated with FTD and ALS in human genome-wide association studies5,6, but how those variants increase risk for disease is unknown. Here we show that TDP-43 represses a cryptic exon-splicing event in UNC13A. Loss of TDP-43 from the nucleus in human brain, neuronal cell lines and motor neurons derived from induced pluripotent stem cells resulted in the inclusion of a cryptic exon in UNC13A mRNA and reduced UNC13A protein expression. The top variants associated with FTD or ALS risk in humans are located in the intron harbouring the cryptic exon, and we show that they increase UNC13A cryptic exon splicing in the face of TDP-43 dysfunction. Together, our data provide a direct functional link between one of the strongest genetic risk factors for FTD and ALS (UNC13A genetic variants), and loss of TDP-43 function.