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

Genetic and polygenic investigation of heart rate variability to identify biomarkers associated with Anxiety disorders

Given that anxiety disorders (AD) are associated with reduced vagally-mediated heart rate variability (HRV), genetic variants related to HRV may provide insight into anxiety etiology. This study used polygenic risk scores (PRS) to explore the genetic overlap between AD and HRV, and investigated whether HRV-related polymorphisms influence anxiety risk. Resting vagally-mediated HRV was measured using a wearable device in 188 European individuals (AD=101, healthy controls=87). AD PRS was tested for association with resting HRV, and HRV PRS for association with AD. We also investigated 15 significant hits from an HRV genome-wide association study (GWAS) for association with resting HRV and AD and if this association is mediated through resting HRV. The AD PRS and HRV PRS showed nominally significant associations with resting HRV and anxiety disorders, respectively. HRV GWAS variants associated with resting HRV were rs12980262 (NDUFA11), rs2680344 (HCN4), rs4262 and rs180238 (GNG11), and rs10842383 (LINC00477). Mediation analyses revealed that NDUFA11 rs12980262 A-carriers and GNG11 rs180238 and rs4262 C-carriers had higher anxiety risk through lower HRV. This study supports an anxiety-HRV genetic relationship, with HRV-related genetic variants translating to AD. This study encourages exploration of HRV genetics to understand mechanisms and identify novel treatment targets for anxiety.

Quality and applicability of cadaveric donor eyes for molecular biology research: An Indian experience

PURPOSE

Human ocular tissue banking plays an important part in the advancement of translational research for identifying the molecular processes involved in disease etiology and pathogenesis. Timely obtaining a good-quality ocular tissue from a cadaveric donor is exceedingly difficult, especially in remote areas, with a variable transportation time (within 12-24 h), raising concerns about RNA quality and its subsequent applications. Therefore, we assessed the utility of retinal tissues from cadaver donor and enucleated eyes based on the RNA quality and gene expression by real-time polymerase chain reaction (PCR).

SETTINGS AND DESIGN

Prospective study.

METHODS

Retina tissues were separated from the donor/enucleated eyes received in the eye bank within 24 h of death (n = 15) and within an hour from OR (n = 3), respectively, and stored immediately at -80 degree. RNA was isolated using trizol, and the quantity and quality were assessed using Qubit and agarose gel electrophoresis, respectively. QPCR was performed for measuring the expression of different retinal-specific genes. The cellular viability of the retina was assessed by establishing explant primary cell cultures.

STATISTICAL ANALYSIS

The data were calculated as an average of normalised Ct values ± standard error of the mean.

RESULTS

RNA obtained from cadaveric tissues despite being partially degraded showed a uniform strong gene expression of several retinal-specific genes such as PAX6, RHO, TUBB3, CRX , and ALDH1L1 . The primary cultures established from cadaveric tissues showed viable cells.

CONCLUSION

The cadaver donor tissues collected within 24 hours of death can be effectively utilized for gene expression profiling.

Identifying plasma proteomic signatures from health to heart failure, across the ejection fraction spectrum

Circulating proteins may provide insights into the varying biological mechanisms involved in heart failure (HF) with preserved ejection fraction (HFpEF) and reduced ejection fraction (HFrEF). We aimed to identify specific proteomic patterns for HF, by comparing proteomic profiles across the ejection fraction spectrum. We investigated 4210 circulating proteins in 739 patients with normal (Stage A/Healthy) or elevated (Stage B) filling pressures, HFpEF, or ischemic HFrEF (iHFrEF). We found 2122 differentially expressed proteins between iHFrEF-Stage A/Healthy, 1462 between iHFrEF-HFpEF and 52 between HFpEF-Stage A/Healthy. Of these 52 proteins, 50 were also found in iHFrEF vs. Stage A/Healthy, leaving SLITRK6 and NELL2 expressed in lower levels only in HFpEF. Moreover, 108 proteins, linked to regulation of cell fate commitment, differed only between iHFrEF-HFpEF. Proteomics across the HF spectrum reveals overlap in differentially expressed proteins compared to stage A/Healthy. Multiple proteins are unique for distinguishing iHFrEF from HFpEF, supporting the capacity of proteomics to discern between these conditions.

Aging promotes metabolic dysfunction-associated steatotic liver disease by inducing ferroptotic stress

Susceptibility to the biological consequences of aging varies among organs and individuals. We analyzed hepatocyte transcriptomes of healthy young and aged male mice to generate an aging hepatocyte gene signature, used it to deconvolute transcriptomic data from humans and mice with metabolic dysfunction-associated liver disease, validated findings with functional studies in mice and applied the signature to transcriptomic data from other organs to determine whether aging-sensitive degenerative mechanisms are conserved. We discovered that the signature enriches in diseased livers in parallel with degeneration. It is also enriched in failing human hearts, diseased kidneys and pancreatic islets from individuals with diabetes. The signature includes genes that control ferroptosis. Aged mice develop more hepatocyte ferroptosis and liver degeneration than young mice when fed diets that induce metabolic stress. Inhibiting ferroptosis shifts the liver transcriptome of old mice toward that of young mice and reverses aging-exacerbated liver damage, identifying ferroptosis as a tractable, conserved mechanism for aging-related tissue degeneration.

A Novel Macrophage Subpopulation Conveys Increased Genetic Risk of Coronary Artery Disease

BACKGROUND

Coronary artery disease (CAD), the leading cause of death worldwide, is influenced by both environmental and genetic factors. Although over 250 genetic risk loci have been identified through genome-wide association studies, the specific causal variants and their regulatory mechanisms are still largely unknown, particularly in disease-relevant cell types such as macrophages.

METHODS

We utilized single-cell RNA-seq and single-cell multiomics approaches in primary human monocyte-derived macrophages to explore the transcriptional regulatory network involved in a critical pathogenic event of coronary atherosclerosis-the formation of lipid-laden foam cells. The relative genetic contribution to CAD was assessed by partitioning disease heritability across different macrophage subpopulations. Meta-analysis of single-cell RNA-seq data sets from 38 human atherosclerotic samples was conducted to provide high-resolution cross-referencing to macrophage subpopulations in vivo.

RESULTS

We identified 18 782 cis-regulatory elements by jointly profiling the gene expression and chromatin accessibility of >5000 macrophages. Integration with CAD genome-wide association study data prioritized 121 CAD-related genetic variants and 56 candidate causal genes. We showed that CAD heritability was not uniformly distributed and was particularly enriched in the gene programs of a novel CD52-hi lipid-handling macrophage subpopulation. These CD52-hi macrophages displayed significantly less lipoprotein accumulation and were also found in human atherosclerotic plaques. We investigated the cis-regulatory effect of a risk variant rs10488763 on FDX1, implicating the recruitment of AP-1 and C/EBP-β in the causal mechanisms at this locus.

CONCLUSIONS

Our results provide genetic evidence of the divergent roles of macrophage subsets in atherogenesis and highlight lipid-handling macrophages as a key subpopulation through which genetic variants operate to influence disease. These findings provide an unbiased framework for functional fine-mapping of genome-wide association study results using single-cell multiomics and offer new insights into the genotype-environment interactions underlying atherosclerotic disease.

RGC-32 mediates proinflammatory and profibrotic pathways in immune-mediated kidney disease

Systemic lupus erythematosus is an autoimmune disease that results in immune-mediated damage to kidneys and other organs. We investigated the role of response gene to complement-32 (RGC-32), a proinflammatory and profibrotic mediator induced by TGFβ and C5b-9, in nephrotoxic nephritis (NTN), an experimental model that mimics human lupus nephritis. Proteinuria, loss of renal function and kidney histopathology were attenuated in RGC-32 KO NTN mice. RGC-32 KO NTN mice displayed downregulation of the CCL20/CCR6 and CXCL9/CXCR3 ligand/receptor pairs resulting in decreased renal recruitment of IL-17+ and IFNγ+ cells and subsequent decrease in the influx of innate immune cells. RGC-32 deficiency attenuated renal fibrosis as demonstrated by decreased deposition of collagen I, III and fibronectin. Thus, RGC-32 is a unique mediator shared by the Th17 and Th1 dependent proinflammatory and profibrotic pathways and a potential novel therapeutic target in the treatment of immune complex mediated glomerulonephritis such as lupus nephritis.

Job-related exhaustion risk variant in UST is associated with dementia and DNA methylation

Previous genome-wide association and replication study for job-related exhaustion indicated a risk variant, rs13219957 in the UST gene. Epidemiological studies suggest connection of stress-related conditions and dementia risk. Therefore, we first studied association of rs13219957 and register-based incident dementia using survival models in the Finnish National FINRISK study surveys (N = 26,693). The AA genotype of rs13219957 was significantly associated with 40% increased risk of all-cause dementia. Then we analysed the UST locus association with brain pathology in the Vantaa 85+ cohort and found association with tau pathology (Braak stage) but not with amyloid pathology. Finally, in the functional analyses, rs13219957 showed a highly significant association with two DNA methylation sites of UST, and UST expression. Thus, the results suggest a common risk variant for a stress-related condition and dementia. Mechanisms to mediate the connection may include differential DNA methylation and transcriptional regulation of UST.

Genome-wide association study of the common retinal disorder epiretinal membrane: Significant risk loci in each of three American populations

Epiretinal membrane (ERM) is a common retinal condition characterized by the presence of fibrocellular tissue on the retinal surface, often with visual distortion and loss of visual acuity. We studied European American (EUR), African American (AFR), and Latino (admixed American, AMR) ERM participants in the Million Veteran Program (MVP) for genome-wide association analysis-a total of 38,232 case individuals and 557,988 control individuals. We completed a genome-wide association study (GWAS) in each population separately, and then results were meta-analyzed. Genome-wide significant (GWS) associations were observed in all three populations studied: 31 risk loci in EUR subjects, 3 in AFR, and 2 in AMR, with 48 in trans-ancestry meta-analysis. Many results replicated in the FinnGen sample. Several GWS variants associate to alterations in gene expression in the macula. ERM showed significant genetic correlation to multiple traits. Pathway enrichment analyses implicated collagen and collagen-adjacent mechanisms, among others. This well-powered ERM GWAS identified novel genetic associations that point to biological mechanisms for ERM.

Relationship between sex biases in gene expression and sex biases in autism and Alzheimer's disease

BACKGROUND

Sex differences in the brain may play an important role in sex-differential prevalence of neuropsychiatric conditions.

METHODS

In order to understand the transcriptional basis of sex differences, we analyzed multiple, large-scale, human postmortem brain RNA-Seq datasets using both within-region and pan-regional frameworks.

RESULTS

We find evidence of sex-biased transcription in many autosomal genes, some of which provide evidence for pathways and cell population differences between chromosomally male and female individuals. These analyses also highlight regional differences in the extent of sex-differential gene expression. We observe an increase in specific neuronal transcripts in male brains and an increase in immune and glial function-related transcripts in female brains. Integration with single-nucleus data suggests this corresponds to sex differences in cellular states rather than cell abundance. Integration with case-control gene expression studies suggests a female molecular predisposition towards Alzheimer's disease, a female-biased disease. Autism, a male-biased diagnosis, does not exhibit a male predisposition pattern in our analysis.

CONCLUSION

Overall, these analyses highlight mechanisms by which sex differences may interact with sex-biased conditions in the brain. Furthermore, we provide region-specific analyses of sex differences in brain gene expression to enable additional studies at the interface of gene expression and diagnostic differences.

Meta-Analysis of Genome-Wide Association Studies Reveals Genetic Mechanisms of Supraventricular Arrhythmias

BACKGROUND

Substantial data support a heritable basis for supraventricular tachycardias, but the genetic determinants and molecular mechanisms of these arrhythmias are poorly understood. We sought to identify genetic loci associated with atrioventricular nodal reentrant tachycardia (AVNRT) and atrioventricular accessory pathways or atrioventricular reciprocating tachycardia (AVAPs/AVRT).

METHODS

We performed multiancestry meta-analyses of genome-wide association studies to identify genetic loci for AVNRT (4 studies) and AVAP/AVRT (7 studies). We assessed evidence supporting the potential causal effects of candidate genes by analyzing relations between associated variants and cardiac gene expression, performing transcriptome-wide analyses, and examining prior genome-wide association studies.

RESULTS

Analyses comprised 2384 AVNRT cases and 106 489 referents, and 2811 AVAP/AVRT cases and 1,483 093 referents. We identified 2 significant loci for AVNRT, which implicate NKX2-5 and TTN as disease susceptibility genes. A transcriptome-wide association analysis supported an association between reduced predicted cardiac expression of NKX2-5 and AVNRT. We identified 3 significant loci for AVAP/AVRT, which implicate SCN5A, SCN10A, and TTN/CCDC141. Variant associations at several loci have been previously reported for cardiac phenotypes, including atrial fibrillation, stroke, Brugada syndrome, and electrocardiographic intervals.

CONCLUSIONS

Our findings highlight gene regions associated with ion channel function (AVAP/AVRT), as well as cardiac development and the sarcomere (AVAP/AVRT and AVNRT) as important potential effectors of supraventricular tachycardia susceptibility.

An eQTL-based approach reveals candidate regulators of LINE-1 RNA levels in lymphoblastoid cells

Long interspersed element 1 (LINE-1; L1) are a family of transposons that occupy ~17% of the human genome. Though a small number of L1 copies remain capable of autonomous transposition, the overwhelming majority of copies are degenerate and immobile. Nevertheless, both mobile and immobile L1s can exert pleiotropic effects (promoting genome instability, inflammation, or cellular senescence) on their hosts, and L1's contributions to aging and aging diseases is an area of active research. However, because of the cell type-specific nature of transposon control, the catalogue of L1 regulators remains incomplete. Here, we employ an eQTL approach leveraging transcriptomic and genomic data from the GEUVADIS and 1000Genomes projects to computationally identify new candidate regulators of L1 RNA levels in lymphoblastoid cell lines. To cement the role of candidate genes in L1 regulation, we experimentally modulate the levels of top candidates in vitro, including IL16, STARD5, HSD17B12, and RNF5, and assess changes in TE family expression by Gene Set Enrichment Analysis (GSEA). Remarkably, we observe subtle but widespread upregulation of TE family expression following IL16 and STARD5 overexpression. Moreover, a short-term 24-hour exposure to recombinant human IL16 was sufficient to transiently induce subtle, but widespread, upregulation of L1 subfamilies. Finally, we find that many L1 expression-associated genetic variants are co-associated with aging traits across genome-wide association study databases. Our results expand the catalogue of genes implicated in L1 RNA control and further suggest that L1-derived RNA contributes to aging processes. Given the ever-increasing availability of paired genomic and transcriptomic data, we anticipate this new approach to be a starting point for more comprehensive computational scans for regulators of transposon RNA levels.

A functional variant of ALDH1A2 is associated with hand osteoarthritis in the Chinese population

Genome-wide association study identified common variants within the ALDH1A2 gene as the susceptible loci of hand osteoarthritis (HOA) in UK and Iceland populations. Located in chromosome 15, ALDH1A2 encodes aldehyde dehydrogenase family 1 member A2, which is an enzyme that catalyses the synthesis of retinoic acid from retinaldehyde. Our purposes were to replicate the association of functional variant in ALDH1A2 with the development of HOA in the Chinese population. Variant rs12915901 of ALDH1A2 was genotyped in 872 HOA patients and 1223 healthy controls. Subchondral bone samples were collected from 40 patients who had undergone a trapeziectomy, and the tissue expression of ALDH1A2 was analysed. The chi-square analysis was used to compare the frequency of genotype and risk allele between the HOA cases and controls. The Student t test was used to compare the mRNA expression of ALDH1A2 between patients with genotype AA/AG and those with genotype GG. The frequency of genotype AA was significantly higher in HOA patients than in the controls (7.6% vs. 5.1%, p = .01). The frequency of allele A was significantly higher in the patients than in the controls (28.9% vs. 24.6%, p = .005). The mRNA expression of ALDH1A2 was 1.31-folds higher in patients with genotype GG than in the patients with genotype AA/AG (0.000617 ± 0.000231 vs. 0.000471 ± 0.000198, p = .04). Variant rs12915901 of ALDH1A2 contributed to the susceptibility of HOA in the Chinese population. Allele A of rs12915901 can add to the risk of HOA possibly via down-regulation of ALDH1A2 expression.

Aplp1 interacts with Lag3 to facilitate transmission of pathologic α-synuclein

Pathologic α-synuclein (α-syn) spreads from cell-to-cell, in part, through binding to the lymphocyte-activation gene 3 (Lag3). Here we report that amyloid β precursor-like protein 1 (Aplp1) interacts with Lag3 that facilitates the binding, internalization, transmission, and toxicity of pathologic α-syn. Deletion of both Aplp1 and Lag3 eliminates the loss of dopaminergic neurons and the accompanying behavioral deficits induced by α-syn preformed fibrils (PFF). Anti-Lag3 prevents the internalization of α-syn PFF by disrupting the interaction of Aplp1 and Lag3, and blocks the neurodegeneration induced by α-syn PFF in vivo. The identification of Aplp1 and the interplay with Lag3 for α-syn PFF induced pathology deepens our insight about molecular mechanisms of cell-to-cell transmission of pathologic α-syn and provides additional targets for therapeutic strategies aimed at preventing neurodegeneration in Parkinson's disease and related α-synucleinopathies.

Genomic dissection and mutation-specific target discovery for breast cancer PIK3CA hotspot mutations

BACKGROUND

Recent advancements in high-throughput genomics and targeted therapies have provided tremendous potential to identify and therapeutically target distinct mutations associated with cancers. However, to date the majority of targeted therapies are used to treat all functional mutations within the same gene, regardless of affected codon or phenotype.

RESULTS

In this study, we developed a functional genomic analysis workflow with a unique isogenic cell line panel bearing two distinct hotspot PIK3CA mutations, E545K and H1047R, to accurately identify targetable differences between mutations within the same gene. We performed RNA-seq and ATAC-seq and identified distinct transcriptomic and epigenomic differences associated with each PIK3CA hotspot mutation. We used this data to curate a select CRISPR knock out screen to identify mutation-specific gene pathway vulnerabilities. These data revealed AREG as a E545K-preferential target that was further validated through in vitro analysis and publicly available patient databases.

CONCLUSIONS

Using our multi-modal genomics framework, we discover distinct differences in genomic regulation between PIK3CA hotspot mutations, suggesting the PIK3CA mutations have different regulatory effects on the function and downstream signaling of the PI3K complex. Our results demonstrate the potential to rapidly uncover mutation specific molecular targets, specifically AREG and a proximal gene regulatory region, that may provide clinically relevant therapeutic targets. The methods outlined provide investigators with an integrative strategy to identify mutation-specific targets for the treatment of other oncogenic mutations in an isogenic system.

Enrichment of a subset of Neanderthal polymorphisms in autistic probands and siblings

Homo sapiens and Neanderthals underwent hybridization during the Middle/Upper Paleolithic age, culminating in retention of small amounts of Neanderthal-derived DNA in the modern human genome. In the current study, we address the potential roles Neanderthal single nucleotide polymorphisms (SNP) may be playing in autism susceptibility in samples of black non-Hispanic, white Hispanic, and white non-Hispanic people using data from the Simons Foundation Powering Autism Research (SPARK), Genotype-Tissue Expression (GTEx), and 1000 Genomes (1000G) databases. We have discovered that rare variants are significantly enriched in autistic probands compared to race-matched controls. In addition, we have identified 25 rare and common SNPs that are significantly enriched in autism on different ethnic backgrounds, some of which show significant clinical associations. We have also identified other SNPs that share more specific genotype-phenotype correlations but which are not necessarily enriched in autism and yet may nevertheless play roles in comorbid phenotype expression (e.g., intellectual disability, epilepsy, and language regression). These results strongly suggest Neanderthal-derived DNA is playing a significant role in autism susceptibility across major populations in the United States.

Genes Differentially Expressed Across Major Arteries Are Enriched in Endothelial Dysfunction-Related Gene Sets: Implications for Relative Inter-artery Atherosclerosis Risk

Atherosclerosis differs across major arteries. Although the biological basis is not fully understood, limited evidence of genetic differences has been documented. This study, therefore, was aimed to identify differentially expressed genes between clinically relevant major arteries and investigate their enrichment in endothelial dysfunction-related gene sets. A bioinformatic analysis of publicly available gene-level read counts for coronary, aortic, and tibial arteries was performed. Differential gene expression was conducted with DeSeq2 at a false discovery rate of 0.05. Differentially expressed genes were then subjected to over-representation analysis and active-subnetwork-oriented enrichment analysis, both at a false discovery rate of 0.005. Enriched terms common to both analyses were categorized for each contrast into immunity/inflammation-, membrane biology-, lipid metabolism-, and coagulation-related terms, and the top differentially expressed genes validated against Swiss Institute of Bioinformatics' Bgee database. There was mostly upregulation of differentially expressed genes for the coronary/tibial and aorta/tibial contrasts, but milder changes for the coronary/aorta contrast. Transcriptomic differences between coronary or aortic versus tibial samples largely involved immunity/inflammation-, membrane biology-, lipid metabolism-, and coagulation-related genes, suggesting potential to modulate endothelial dysfunction and atherosclerosis. These results imply atheroprone coronary and aortic environments compared with tibial artery tissue, which may explain observed relative inter-artery atherosclerosis risk.

Brain Expression Levels of Commonly Measured Blood Biomarkers of Neurological Damage Differ with Respect to Sex, Race, and Age

It is increasingly evident that blood biomarkers have potential to improve the diagnosis and management of both acute and chronic neurological conditions. The most well-studied candidates, and arguably those with the broadest utility, are proteins that are highly enriched in neural tissues and released into circulation upon cellular damage. It is currently unknown how the brain expression levels of these proteins is influenced by demographic factors such as sex, race, and age. Given that source tissue abundance is likely a key determinant of the levels observed in the blood during neurological pathology, understanding such influences is important in terms of identifying potential clinical scenarios that could produce diagnostic bias. In this study, we leveraged existing mRNA sequencing data originating from 2,642 normal brain specimens harvested from 382 human donors to examine potential demographic variability in the expression levels of genes which code for 28 candidate blood biomarkers of neurological damage. Existing mass spectrometry data originating from 26 additional normal brain specimens harvested from 26 separate human donors was subsequently used to tentatively assess whether observed transcriptional variance was likely to produce corresponding variance in terms of protein abundance. Genes associated with several well-studied or emerging candidate biomarkers including neurofilament light chain (NfL), ubiquitin carboxyl-terminal hydrolase isozyme L1 (UCH-L1), neuron-specific enolase (NSE), and synaptosomal-associated protein 25 (SNAP-25) exhibited significant differences in expression with respect to sex, race, and age. In many instances, these differences in brain expression align well with and provide a mechanistic explanation for previously reported differences in blood levels.

A comparison of RNA-Seq data preprocessing pipelines for transcriptomic predictions across independent studies

BACKGROUND

RNA sequencing combined with machine learning techniques has provided a modern approach to the molecular classification of cancer. Class predictors, reflecting the disease class, can be constructed for known tissue types using the gene expression measurements extracted from cancer patients. One challenge of current cancer predictors is that they often have suboptimal performance estimates when integrating molecular datasets generated from different labs. Often, the quality of the data is variable, procured differently, and contains unwanted noise hampering the ability of a predictive model to extract useful information. Data preprocessing methods can be applied in attempts to reduce these systematic variations and harmonize the datasets before they are used to build a machine learning model for resolving tissue of origins.

RESULTS

We aimed to investigate the impact of data preprocessing steps-focusing on normalization, batch effect correction, and data scaling-through trial and comparison. Our goal was to improve the cross-study predictions of tissue of origin for common cancers on large-scale RNA-Seq datasets derived from thousands of patients and over a dozen tumor types. The results showed that the choice of data preprocessing operations affected the performance of the associated classifier models constructed for tissue of origin predictions in cancer.

CONCLUSION

By using TCGA as a training set and applying data preprocessing methods, we demonstrated that batch effect correction improved performance measured by weighted F1-score in resolving tissue of origin against an independent GTEx test dataset. On the other hand, the use of data preprocessing operations worsened classification performance when the independent test dataset was aggregated from separate studies in ICGC and GEO. Therefore, based on our findings with these publicly available large-scale RNA-Seq datasets, the application of data preprocessing techniques to a machine learning pipeline is not always appropriate.

Psychobiological regulation of plasma and saliva GDF15 dynamics in health and mitochondrial diseases

GDF15 (growth differentiation factor 15) is a marker of cellular energetic stress linked to physical-mental illness, aging, and mortality. However, questions remain about its dynamic properties and measurability in human biofluids other than blood. Here, we examine the natural dynamics and psychobiological regulation of plasma and saliva GDF15 in four human studies representing 4,749 samples from 188 individuals. We show that GDF15 protein is detectable in saliva (8% of plasma concentration), likely produced by salivary glands secretory duct cells. Using a brief laboratory socio-evaluative stressor paradigm, we find that psychosocial stress increases plasma (+3.5-5.9%) and saliva GDF15 (+43%) with distinct kinetics, within minutes. Moreover, saliva GDF15 exhibits a robust awakening response, declining by ~40-89% within 30-45 minutes from its peak level at the time of waking up. Clinically, individuals with genetic mitochondrial OxPhos diseases show elevated baseline plasma and saliva GDF15, and post-stress GDF15 levels in both biofluids correlate with multi-system disease severity, exercise intolerance, and the subjective experience of fatigue. Taken together, our data establish that saliva GDF15 is dynamic, sensitive to psychological states, a clinically relevant endocrine marker of mitochondrial diseases. These findings also point to a shared psychobiological pathway integrating metabolic and mental stress.

Association of IL-18 gene polymorphisms with clinical aspects of hyperlipidemia in middle-aged and early people in the community

Hyperlipidemia is notorious for causing coronary artery disease (CAD). IL-18 is a proinflammtory cytokine that contributes to the pathogenesis of CAD. Previous reports have revealed that genetic polymorphism of IL-18 is associated with its expression level as well as the susceptibility to CAD. In the present study, we aim to investigate the relationship between IL-18 single nucleotide polymorphisms (SNPs) and hyperlipidemia in the Han Chinese population in Taiwan. A total of 580 participants older than 30 were recruited from the community. We collected the demographics, self-reported disease histories, and lifestyles. We also assessed the levels of lipid profiles including total cholesterol (CHOL), triglyceride, low-density lipoprotein cholesterol (LDL-C) and high-density lipoprotein cholesterol. Two SNPs, rs3882891C/A (intron 5) and rs1946518A/C (promoter -607) of IL-18 were elucidated by polymerase chain reaction-restriction fragment length polymorphism (PCR-RFLP) methods. Our results revealed that rs3882891 AA was associated with lower risk of hypercholesterolemia, higher CHOL and LDL-C in subjects (p=0.003, p=0.000 and p=0.005 separately), and rs1946518 CC was associated with hypercholesterolemia, higher CHOL and LDL-C as well (p=0.021, p=0.003 and p=0.001 separately) Furthermore, both SNPs were associated with IL-18 expression level, which was examined by Genotype-Tissue Expression (GTEx) Portal (p=0.042 and 0.016 separately). Finally, the haplotype of IL-18 was subsequently arranged in the order of rs3882891 and rs1946518. The result revealed that the AC haplotype of 2 IL-18 SNPs was also associated with lower risk of hypercholesterolemia, lower levels of CHOL and LDL-C (p=0.01, p=0.001 and 0.003). The current study is the first to report the association between IL-18 SNPs and hyperlipidemia in the Chinese Han population.

Haploinsufficiency of phosphodiesterase 10A activates PI3K/AKT signaling independent of PTEN to induce an aggressive glioma phenotype

Glioblastoma is universally fatal and characterized by frequent chromosomal copy number alterations harboring oncogenes and tumor suppressors. In this study, we analyzed exome-wide human glioblastoma copy number data and found that cytoband 6q27 is an independent poor prognostic marker in multiple data sets. We then combined CRISPR-Cas9 data, human spatial transcriptomic data, and human and mouse RNA sequencing data to nominate PDE10A as a potential haploinsufficient tumor suppressor in the 6q27 region. Mouse glioblastoma modeling using the RCAS/tv-a system confirmed that Pde10a suppression induced an aggressive glioma phenotype in vivo and resistance to temozolomide and radiation therapy in vitro. Cell culture analysis showed that decreased Pde10a expression led to increased PI3K/AKT signaling in a Pten-independent manner, a response blocked by selective PI3K inhibitors. Single-nucleus RNA sequencing from our mouse gliomas in vivo, in combination with cell culture validation, further showed that Pde10a suppression was associated with a proneural-to-mesenchymal transition that exhibited increased cell adhesion and decreased cell migration. Our results indicate that glioblastoma patients harboring PDE10A loss have worse outcomes and potentially increased sensitivity to PI3K inhibition.

Shared genetic architecture between autoimmune disorders and B-cell acute lymphoblastic leukemia: insights from large-scale genome-wide cross-trait analysis

BACKGROUND

To study the shared genetic structure between autoimmune diseases and B-cell acute lymphoblastic leukemia (B-ALL) and identify the shared risk loci and genes and genetic mechanisms involved.

METHODS

Based on large-scale genome-wide association study (GWAS) summary-level data sets, we observed genetic overlaps between autoimmune diseases and B-ALL, and cross-trait pleiotropic analysis was performed to detect shared pleiotropic loci and genes. A series of functional annotation and tissue-specific analysis were performed to determine the influence of pleiotropic genes. The heritability enrichment analysis was used to detect crucial immune cells and tissues. Finally, bidirectional Mendelian randomization (MR) methods were utilized to investigate the casual associations.

RESULTS

Our research highlighted shared genetic mechanisms between seven autoimmune disorders and B-ALL. A total of 73 pleiotropic loci were identified at the genome-wide significance level (P < 5 × 10-8), 16 of which had strong evidence of colocalization. We demonstrated that several loci have been previously reported (e.g., 17q21) and discovered some novel loci (e.g., 10p12, 5p13). Further gene-level identified 194 unique pleiotropic genes, for example IKZF1, GATA3, IKZF3, GSDMB, and ORMDL3. Pathway analysis determined the key role of cellular response to cytokine stimulus, B cell activation, and JAK-STAT signaling pathways. SNP-level and gene-level tissue enrichment suggested that crucial role pleiotropic mechanisms involved in the spleen, whole blood, and EBV-transformed lymphocytes. Also, hyprcoloc and stratified LD score regression analyses revealed that B cells at different developmental stages may be involved in mechanisms shared between two different diseases. Finally, two-sample MR analysis determined causal effects of asthma and rheumatoid arthritis on B-ALL.

CONCLUSIONS

Our research proved shared genetic architecture between autoimmune disorders and B-ALL and shed light on the potential mechanism that might involve in.

The association of cigarette smoking with DNA methylation and gene expression in human tissue samples

Cigarette smoking adversely affects many aspects of human health, and epigenetic responses to smoking may reflect mechanisms that mediate or defend against these effects. Prior studies of smoking and DNA methylation (DNAm), typically measured in leukocytes, have identified numerous smoking-associated regions (e.g., AHRR). To identify smoking-associated DNAm features in typically inaccessible tissues, we generated array-based DNAm data for 916 tissue samples from the GTEx (Genotype-Tissue Expression) project representing 9 tissue types (lung, colon, ovary, prostate, blood, breast, testis, kidney, and muscle). We identified 6,350 smoking-associated CpGs in lung tissue (n = 212) and 2,735 in colon tissue (n = 210), most not reported previously. For all 7 other tissue types (sample sizes 38-153), no clear associations were observed (false discovery rate 0.05), but some tissues showed enrichment for smoking-associated CpGs reported previously. For 1,646 loci (in lung) and 22 (in colon), smoking was associated with both DNAm and local gene expression. For loci detected in both lung and colon (e.g., AHRR, CYP1B1, CYP1A1), top CpGs often differed between tissues, but similar clusters of hyper- or hypomethylated CpGs were observed, with hypomethylation at regulatory elements corresponding to increased expression. For lung tissue, 17 hallmark gene sets were enriched for smoking-associated CpGs, including xenobiotic- and cancer-related gene sets. At least four smoking-associated regions in lung were impacted by lung methylation quantitative trait loci (QTLs) that co-localize with genome-wide association study (GWAS) signals for lung function (FEV1/FVC), suggesting epigenetic alterations can mediate the effects of smoking on lung health. Our multi-tissue approach has identified smoking-associated regions in disease-relevant tissues, including effects that are shared across tissue types.

Creation of signatures and identification of molecular subtypes based on disulfidptosis-related genes for glioblastoma patients' prognosis and immunological activity

BACKGROUND

In recent times, disulfidptosis, an intricate form of cellular demise, has garnered attention due to its impact on prognosis, tumor progression and treatment response. Nevertheless, the exact significance of disulfidptosis-related genes (DisRGs) in glioblastoma (GBM) remains enigmatic.

METHODS

The GEO and TCGA databases provided transcriptional and clinically relevant data on tumor samples, while the GTEx database provided data on healthy tissues. Disulfidptosis-related genes (DisRGs) were procured from previous scholarly investigations. The expression profile of DisRGs was initially scrutinized among patients diagnosed with GBM, subsequent to which their prognostic value was explored. Through consensus clustering, we constructed DisRGs-related clusters and gene subtypes. Our results established that the DisRG-related clusters had differentially expressed genes, resulting in a DisulfidptosisScore model, which had a positive prognostic value.

RESULTS

The differential expression profile of 24 DisRGs between GBM samples and healthy samples was acquired. Through consensus cluster analysis, two distinct disulfidptosis subtypes, namely DisRGcluster A and DisRGcluster B, were identified. Then, the DisulfidptosisScore model including 4 characteristic genes was constructed.Notably, patients with GBM assigned with lower score demonstrated a considerably longer overall survival (OS) compared to those with higher score.

CONCLUSION

We have effectively devised a prognostic model associated with disulfidptosis, presenting autonomous prognostic predictions for patients with GBM. These findings serve as a valuable addition to the current comprehension of disulfidptosis and offer fresh theoretical substantiation for the development of enhanced treatment strategies.

A loss-of-function variant in ZCWPW1 causes human male infertility with sperm head defect and high DNA fragmentation

BACKGROUND

Male infertility is a global health issue. The more causative genes related to human male infertility should be further explored. The essential role of Zcwpw1 in male mouse fertility has been established and the role of ZCWPW1 in human reproduction needs further investigation to verify.

METHODS

An infertile man with oligoasthenoteratozoospermia phenotype and his parents were recruited from West China Second University Hospital, Sichuan University. A total of 200 healthy Han Chinese volunteers without any evidence of infertility were recruited as normal controls, while an additional 150 infertile individuals were included to assess the prevalence of ZCWPW1 variants in a sporadic male sterile population. The causative gene variant was identified by Whole-exome sequencing and Sanger sequencing. The phenotype of the oligoasthenoteratozoospermia was determined by Papanicolaou staining, immunofluorescence staining and electron microscope. In-vitro experiments, western blot and in-silicon analysis were applied to assess the pathogenicity of the identified variant. Additionally, we examined the influence of the variant on the DNA fragmentation and DNA repair capability by Sperm Chromatin Dispersion and Neutral Comet Assay.

RESULTS

The proband exhibits a phenotype of oligoasthenoteratozoospermia, his spermatozoa show head defects by semen examination, Papanicolaou staining and electron microscope assays. Whole-exome sequencing and Sanger sequencing found the proband carries a homozygous ZCWPW1 variant (c.1064C > T, p. P355L). Immunofluorescence analysis shows a significant decrease in ZCWPW1 expression in the proband's sperm. By exogenous expression with ZCWPW1 mutant plasmid in vitro, the obvious declined expression of ZCWPW1 with the mutation is validated in HEK293T. After being treated by hydroxyurea, MUT-ZCWPW1 transfected cells and empty vector transfected cells have a higher level of γ-H2AX, increased tail DNA and reduced H3K9ac level than WT-ZCWPW1 transfected cells. Furthermore, the Sperm Chromatin Dispersion assay revealed the proband's spermatozoa have high DNA fragmentation.

CONCLUSIONS

It is the first report that a novel homozygous missense mutation in ZCWPW1 caused human male infertility with sperm head defects and high DNA fragmentation. This finding enriches the gene variant spectrum and etiology of oligoasthenoteratozoospermia.

PRAME induces genomic instability in uveal melanoma

PRAME is a CUL2 ubiquitin ligase subunit that is normally expressed in the testis but becomes aberrantly overexpressed in many cancer types in association with aneuploidy and metastasis. Here, we show that PRAME is expressed predominantly in spermatogonia around the time of meiotic crossing-over in coordination with genes mediating DNA double strand break repair. Expression of PRAME in somatic cells upregulates pathways involved in meiosis, chromosome segregation and DNA repair, and it leads to increased DNA double strand breaks, telomere dysfunction and aneuploidy in neoplastic and non-neoplastic cells. This effect is mediated at least in part by ubiquitination of SMC1A and altered cohesin function. PRAME expression renders cells susceptible to inhibition of PARP1/2, suggesting increased dependence on alternative base excision repair pathways. These findings reveal a distinct oncogenic function of PRAME that can be targeted therapeutically in cancer.

Transcription readthrough is prevalent in healthy human tissues and associated with inherent genomic features

Transcription termination is a crucial step in the production of conforming mRNAs and functional proteins. Under cellular stress conditions, the transcription machinery fails to identify the termination site and continues transcribing beyond gene boundaries, a phenomenon designated as transcription readthrough. However, the prevalence and impact of this phenomenon in healthy human tissues remain unexplored. Here, we assessed transcription readthrough in almost 3000 transcriptome profiles representing 23 human tissues and found that 34% of the expressed protein-coding genes produced readthrough transcripts. The production of readthrough transcripts was restricted in genomic regions with high transcriptional activity and was associated with inefficient splicing and increased chromatin accessibility in terminal regions. In addition, we showed that these transcripts contained several binding sites for the same miRNA, unravelling a potential role as miRNA sponges. Overall, this work provides evidence that transcription readthrough is pervasive and non-stochastic, not only in abnormal conditions but also in healthy tissues. This suggests a potential role for such transcripts in modulating normal cellular functions.

A Pregnancy and Childhood Epigenetics Consortium (PACE) meta-analysis highlights potential relationships between birth order and neonatal blood DNA methylation

Higher birth order is associated with altered risk of many disease states. Changes in placentation and exposures to in utero growth factors with successive pregnancies may impact later life disease risk via persistent DNA methylation alterations. We investigated birth order with Illumina DNA methylation array data in each of 16 birth cohorts (8164 newborns) with European, African, and Latino ancestries from the Pregnancy and Childhood Epigenetics Consortium. Meta-analyzed data demonstrated systematic DNA methylation variation in 341 CpGs (FDR adjusted P < 0.05) and 1107 regions. Forty CpGs were located within known quantitative trait loci for gene expression traits in blood, and trait enrichment analysis suggested a strong association with immune-related, transcriptional control, and blood pressure regulation phenotypes. Decreasing fertility rates worldwide with the concomitant increased proportion of first-born children highlights a potential reflection of birth order-related epigenomic states on changing disease incidence trends.

Cerebrospinal Fluid C1-Esterase Inhibitor and Tie-1 Levels Affect Cognitive Performance: Evidence from Proteome-Wide Mendelian Randomization

OBJECTIVE

The association of cerebrospinal fluid (CSF) protein levels with cognitive function in the general population remains largely unexplored. We performed Mendelian randomization (MR) analyses to query which CSF proteins may have potential causal effects on cognitive performance.

METHODS AND ANALYSIS

Genetic associations with CSF proteins were obtained from a genome-wide association study conducted in up to 835 European-ancestry individuals and for cognitive performance from a meta-analysis of GWAS including 257,841 European-ancestry individuals. We performed Mendelian randomization (MR) analyses to test the effect of randomly allocated variation in 154 genetically predicted CSF protein levels on cognitive performance. Findings were validated by performing colocalization analyses and considering cognition-related phenotypes.

RESULTS

Genetically predicted C1-esterase inhibitor levels in the CSF were associated with a better cognitive performance (SD units of cognitive performance per 1 log-relative fluorescence unit (RFU): 0.23, 95% confidence interval: 0.12 to 0.35, p = 7.91 × 10-5), while tyrosine-protein kinase receptor Tie-1 (sTie-1) levels were associated with a worse cognitive performance (-0.43, -0.62 to -0.23, p = 2.08 × 10-5). These findings were supported by colocalization analyses and by concordant effects on distinct cognition-related and brain-volume measures.

CONCLUSIONS

Human genetics supports a role for the C1-esterase inhibitor and sTie-1 in cognitive performance.

Genomic dissection and mutation-specific target discovery for breast cancer PIK3CA hotspot mutations

Background

Recent advancements in high-throughput genomics and targeted therapies have provided tremendous potential to identify and therapeutically target distinct mutations associated with cancers. However, to date the majority of targeted therapies are used to treat all functional mutations within the same gene, regardless of affected codon or phenotype.

Results

In this study, we developed a functional genomic analysis workflow with a unique isogenic cell line panel bearing two distinct hotspot PIK3CA mutations, E545K and H1047R, to accurately identify targetable differences between mutations within the same gene. We performed RNA-seq and ATAC-seq and identified distinct transcriptomic and epigenomic differences associated with each PIK3CA hotspot mutation. We used this data to curate a select CRISPR knock out screen to identify mutation-specific gene pathway vulnerabilities. These data revealed AREG as a E545K-preferential target that was further validated through in vitro analysis and publicly available patient databases.

Conclusions

Using our multi-modal genomics framework, we discover distinct differences in genomic regulation between PIK3CA hotspot mutations, suggesting the PIK3CA mutations have different regulatory effects on the function and downstream signaling of the PI3K complex. Our results demonstrate the potential to rapidly uncover mutation specific molecular targets, specifically AREG and a proximal gene regulatory region, that may provide clinically relevant therapeutic targets. The methods outlined provide investigators with an integrative strategy to identify mutation-specific targets for the treatment of other oncogenic mutations in an isogenic system.

Comprehensive characterization of coding and non-coding single nucleotide polymorphisms of the Myoneurin (MYNN) gene using molecular dynamics simulation and docking approaches

Genome-wide association studies (GWAS) identified a coding single nucleotide polymorphism, MYNN rs10936599, at chromosome 3q. MYNN gene encodes myoneurin protein, which has been associated with several cancer pathogenesis and disease development processes. However, there needed to be a more detailed characterization of this polymorphism's (and other coding and non-coding polymorphisms) structural, functional, and molecular impact. The current study addressed this gap and analyzed different properties of rs10936599 and non-coding SNPs of MYNN via a thorough computational method. The variant, rs10936599, was predicted functionally deleterious by nine functionality prediction approaches, like SIFT, PolyPhen-2, and REVEL, etc. Following that, structural modifications were estimated through the HOPE server and Mutation3D. Moreover, the mutation was found in a conserved and active residue, according to ConSurf and CPORT. Further, the secondary structures were predicted, followed by tertiary structures, and there was a significant deviation between the native and variant models. Similarly, molecular simulation also showed considerable differences in the dynamic pattern of the wildtype and mutant structures. Molecular docking revealed that the variant binds with better docking scores with ligand NOTCH2. In addition to that, non-coding SNPs located at the MYNN locus were retrieved from the ENSEMBL database. These were found to disrupt the transcription factor binding regulatory regions; nonetheless, only two affect miRNA target sites. Again, eight non-coding variants were detected in the testes with normalized expression, whereas HaploReg v4.1 unveiled annotations for non-coding variants. In summary, in silico comprehensive characterization of coding and non-coding single nucleotide polymorphisms of MYNN gene will assist researchers to work on MYNN gene and establish their association with certain types of cancers.

Computationally defined and in vitro validated putative genomic safe harbour loci for transgene expression in human cells

Selection of the target site is an inherent question for any project aiming for directed transgene integration. Genomic safe harbour (GSH) loci have been proposed as safe sites in the human genome for transgene integration. Although several sites have been characterised for transgene integration in the literature, most of these do not meet criteria set out for a GSH and the limited set that do have not been characterised extensively. Here, we conducted a computational analysis using publicly available data to identify 25 unique putative GSH loci that reside in active chromosomal compartments. We validated stable transgene expression and minimal disruption of the native transcriptome in three GSH sites in vitro using human embryonic stem cells (hESCs) and their differentiated progeny. Furthermore, for easy targeted transgene expression, we have engineered constitutive landing pad expression constructs into the three validated GSH in hESCs.

Asthma exacerbations and eosinophilia in the UK Biobank: a genome-wide association study

Background

Asthma exacerbations reflect disease severity, affect morbidity and mortality, and may lead to declining lung function. Inflammatory endotypes (e.g. T2-high (eosinophilic)) may play a key role in asthma exacerbations. We aimed to assess whether genetic susceptibility underlies asthma exacerbation risk and additionally tested for an interaction between genetic variants and eosinophilia on exacerbation risk.

Methods

UK Biobank data were used to perform a genome-wide association study of individuals with asthma and at least one exacerbation compared to individuals with asthma and no history of exacerbations. Individuals with asthma were identified using self-reported data, hospitalisation data and general practitioner records. Exacerbations were identified as either asthma-related hospitalisation, general practitioner record of asthma exacerbation or an oral corticosteroid burst prescription. A logistic regression model adjusted for age, sex, smoking status and genetic ancestry via principal components was used to assess the association between genetic variants and asthma exacerbations. We sought replication for suggestive associations (p<5×10-6) in the GERA cohort.

Results

In the UK Biobank, we identified 11 604 cases and 37 890 controls. While no variants reached genome-wide significance (p<5×10-8) in the primary analysis, 116 signals were suggestively significant (p<5×10-6). In GERA, two single nucleotide polymorphisms (rs34643691 and rs149721630) replicated (p<0.05), representing signals near the NTRK3 and ABCA13 genes.

Conclusions

Our study has identified reproducible associations with asthma exacerbations in the UK Biobank and GERA cohorts. Confirmation of these findings in different asthma subphenotypes in diverse ancestries and functional investigation will be required to understand their mechanisms of action and potentially inform therapeutic development.

Gene Co-Expression Analysis of Multiple Brain Tissues Reveals Correlation of FAM222A Expression with Multiple Alzheimer's Disease-Related Genes

Background

Alzheimer's disease (AD) is the most common form of dementia in the elderly marked by central nervous system (CNS) neuronal loss and amyloid plaques. FAM222A, encoding an amyloid plaque core protein, is an AD brain atrophy susceptibility gene that mediates amyloid-β aggregation. However, the expression interplay between FAM222A and other AD-related pathway genes is unclear.

Objective

Our goal was to study FAM222A's whole-genome co-expression profile in multiple tissues and investigate its interplay with other AD-related genes.

Methods

We analyzed gene expression correlations in Genotype-Tissue Expression (GTEx) tissues to identify FAM222A co-expressed genes and performed functional enrichment analysis on identified genes in CNS system.

Results

Genome-wide gene expression profiling identified 673 genes significantly correlated with FAM222A (p < 2.5×10-6) in 48 human tissues, including 298 from 13 CNS tissues. Functional enrichment analysis revealed that FAM222A co-expressed CNS genes were enriched in multiple AD-related pathways. Gene co-expression network analysis for identified genes in each brain region predicted other disease associated genes with similar biological function. Furthermore, co-expression of 25 out of 31 AD-related pathways genes with FAM222A was replicated in brain samples from 107 aged subjects from the Aging, Dementia and TBI Study.

Conclusion

This gene co-expression study identified multiple AD-related genes that are associated with FAM222A, indicating that FAM222A and AD-associated genes can be active simultaneously in similar biological processes, providing evidence that supports the association of FAM222A with AD.

Molecular Regulation of the CNS by Vitamin D

Vitamin D is a lipid-soluble vitamin that can be found in some foods. It is also produced endogenously (in the presence of ultraviolet light), transported through the blood to the targets organs and this is the reason to consider vitamin D as a hormone. It is known that vitamin D has genomic and non-genomic effects. This review is focused mainly on the vitamin D receptors, the importance of vitamin D as a neuromodulator, the role of vitamin D in the pathophysiology of devastating neurological disorders such as Alzheimer's disease, multiple sclerosis, amyotrophic lateral sclerosis, Parkinson's disease and the benefit of vitamin D and its derivates in alleviating these disorders.

ncRNAs: an unexplored cellular defense mechanism in leprosy

Leprosy is an infectious disease primarily caused by the obligate intracellular parasite Mycobacterium leprae. Although it has been considered eradicated in many countries, leprosy continues to be a health issue in developing nations. Besides the social stigma associated with it, individuals affected by leprosy may experience nerve damage leading to physical disabilities if the disease is not properly treated or early diagnosed. Leprosy is recognized as a complex disease wherein socioenvironmental factors, immune response, and host genetics interact to contribute to its development. Recently, a new field of study called epigenetics has emerged, revealing that the immune response and other mechanisms related to infectious diseases can be influenced by noncoding RNAs. This review aims to summarize the significant advancements concerning non-coding RNAs in leprosy, discussing the key perspectives on this novel approach to comprehending the pathophysiology of the disease and identifying molecular markers. In our view, investigations on non-coding RNAs in leprosy hold promise and warrant increased attention from researches in this field.

An eQTL-based Approach Reveals Candidate Regulators of LINE-1 RNA Levels in Lymphoblastoid Cells

Long interspersed element 1 (L1) are a family of autonomous, actively mobile transposons that occupy ~17% of the human genome. A number of pleiotropic effects induced by L1 (promoting genome instability, inflammation, or cellular senescence) have been observed, and L1's contributions to aging and aging diseases is an area of active research. However, because of the cell type-specific nature of transposon control, the catalogue of L1 regulators remains incomplete. Here, we employ an eQTL approach leveraging transcriptomic and genomic data from the GEUVADIS and 1000Genomes projects to computationally identify new candidate regulators of L1 RNA levels in lymphoblastoid cell lines. To cement the role of candidate genes in L1 regulation, we experimentally modulate the levels of top candidates in vitro, including IL16, STARD5, HSDB17B12, and RNF5, and assess changes in TE family expression by Gene Set Enrichment Analysis (GSEA). Remarkably, we observe subtle but widespread upregulation of TE family expression following IL16 and STARD5 overexpression. Moreover, a short-term 24-hour exposure to recombinant human IL16 was sufficient to transiently induce subtle, but widespread, upregulation of L1 subfamilies. Finally, we find that many L1 expression-associated genetic variants are co-associated with aging traits across genome-wide association study databases. Our results expand the catalogue of genes implicated in L1 RNA control and further suggest that L1-derived RNA contributes to aging processes. Given the ever-increasing availability of paired genomic and transcriptomic data, we anticipate this new approach to be a starting point for more comprehensive computational scans for transposon transcriptional regulators.

Genome-Wide Association Study of Blood Mercury in European Pregnant Women and Children

Mercury has high industrial utility and is present in many products, and environmental contamination and occupational exposure are widespread. There are numerous biological systems involved in the absorption, metabolism, and excretion of Hg, and it is possible that some systems may be impacted by genetic variation. If so, genotype may affect tissue concentrations of Hg and subsequent toxic effects. Genome-wide association testing was performed on blood Hg samples from pregnant women of the Avon Longitudinal Study of Parents and Children (n = 2893) and children of the Human Early Life Exposome (n = 1042). Directly-genotyped single-nucleotide polymorphisms (SNPs) were imputed to the Haplotype Reference Consortium r1.1 panel of whole genotypes and modelled againstlog-transformed Hg. Heritability was estimated using linkage disequilibrium score regression. The heritability of Hg was estimated as 24.0% (95% CI: 16.9% to 46.4%) in pregnant women, but could not be determined in children. There were 16 SNPs associated with Hg in pregnant women above a suggestive p-value threshold (p < 1 × 10-5), and 21 for children. However, no SNP passed this threshold in both studies, and none were genome-wide significant (p < 5 × 10-8). SNP-Hg associations were highly discordant between women and children, and this may reflect differences in metabolism, a gene-age interaction, or dose-response effects. Several suggestive variants had plausible links to Hg metabolism, such as rs146099921 in metal transporter SLC39A14, and two variants (rs28618224, rs7154700) in potassium voltage-gated channel genes. The findings would benefit from external validation, as suggestive results may contain both true associations and false positives.

Sctensor detects many-to-many cell-cell interactions from single cell RNA-sequencing data

BACKGROUND

Complex biological systems are described as a multitude of cell-cell interactions (CCIs). Recent single-cell RNA-sequencing studies focus on CCIs based on ligand-receptor (L-R) gene co-expression but the analytical methods are not appropriate to detect many-to-many CCIs.

RESULTS

In this work, we propose scTensor, a novel method for extracting representative triadic relationships (or hypergraphs), which include ligand-expression, receptor-expression, and related L-R pairs.

CONCLUSIONS

Through extensive studies with simulated and empirical datasets, we have shown that scTensor can detect some hypergraphs that cannot be detected using conventional CCI detection methods, especially when they include many-to-many relationships. scTensor is implemented as a freely available R/Bioconductor package.

MiR-181a targets STING to drive PARP inhibitor resistance in BRCA- mutated triple-negative breast cancer and ovarian cancer

BACKGROUND

Poly (ADP-ribose) polymerase inhibitors (PARPi) are approved for the treatment of BRCA-mutated breast cancer (BC), including triple-negative BC (TNBC) and ovarian cancer (OvCa). A key challenge is to identify the factors associated with PARPi resistance; although, previous studies suggest that platinum-based agents and PARPi share similar resistance mechanisms.

METHODS

Olaparib-resistant (OlaR) cell lines were analyzed using HTG EdgeSeq miRNA Whole Transcriptomic Analysis (WTA). Functional assays were performed in three BRCA-mutated TNBC cell lines. In-silico analysis were performed using multiple databases including The Cancer Genome Atlas, the Genotype-Tissue Expression, The Cancer Cell Line Encyclopedia, Genomics of Drug Sensitivity in Cancer, and Gene Omnibus Expression.

RESULTS

High miR-181a levels were identified in OlaR TNBC cell lines (p = 0.001) as well as in tumor tissues from TNBC patients (p = 0.001). We hypothesized that miR-181a downregulates the stimulator of interferon genes (STING) and the downstream proinflammatory cytokines to mediate PARPi resistance. BRCA1 mutated TNBC cell lines with miR-181a-overexpression were more resistant to olaparib and showed downregulation in STING and the downstream genes controlled by STING. Extracellular vesicles derived from PARPi-resistant TNBC cell lines horizontally transferred miR-181a to parental cells which conferred PARPi-resistance and targeted STING. In clinical settings, STING levels were positively correlated with interferon gamma (IFNG) response scores (p = 0.01). In addition, low IFNG response scores were associated with worse response to neoadjuvant treatment including PARPi for high-risk HER2 negative BC patients (p = 0.001). OlaR TNBC cell lines showed resistance to platinum-based drugs. OvCa cell lines resistant to platinum showed resistance to olaparib. Knockout of miR-181a significantly improved olaparib sensitivity in OvCa cell lines (p = 0.001).

CONCLUSION

miR-181a is a key factor controlling the STING pathway and driving PARPi and platinum-based drug resistance in TNBC and OvCa. The miR-181a-STING axis can be used as a potential marker for predicting PARPi responses in TNBC and OvCa tumors.

Human endogenous retroviruses of the HERV-K (HML-2) family are expressed in the brain of healthy individuals and modify the composition of the brain-infiltrating immune cells

Human endogenous retroviruses (HERVs) are remnants of ancient retroviral infections in the human genome. RNA expression of individual HERVs has frequently been observed in various pathologic conditions, but some activity can also be seen in healthy individuals, e.g. in the blood. To quantitate the basal expression levels of HERVs in the brain, we now used high-throughput sequencing-based metagenomic analysis to characterize the expression profiles of the HERV-K (HML-2) family proviruses in different brain regions of healthy brain tissue. To this end, RNA-seq data from the Genotype-Tissue Expression (GTEx) project was used. The GTEx project is a public resource to study tissue-specific gene expression and regulation, consisting of a large selection of sequenced samples from different tissues. The GTEx data used in this study consisted of 378 samples taken from 13 brain regions from 55 individuals. The data demonstrated that out of 99 intact proviruses in the family 58 were expressed, but the expression profiles were highly divergent and there were no significant differences in the expression profiles between the various anatomic regions of the brain. It is known that the brain contains a variety of infiltrating immune cells, which are probably of great importance both in the normal defense mechanisms as well as in the various pathogenic processes. Digital cytometry (CIBERSORTx) was used to quantify the proportions of the infiltrating immune cells in the same brain samples. Six most abundant (>5 % of the total population) cell types were observed to be CD4 memory resting T cells, M0 macrophages, plasma cells, CD8 T cells, CD4 memory activated T cells, and monocytes. Analysis of the correlations between the individual HERVs and infiltrating cell types indicated that a cluster of 6 HERVs had a notable correlation signature between T cell type infiltrating cell proportions and HERV RNA expression intensity. The correlations between inflammatory type infiltrating cells were negative or weak. Taken together, these data indicate that the expression of HERVs is associated with a T cell type immunity.

A Comprehensive Evaluation of the Effects of RNA-Editing Proteins ADAR and ADARB1 on the Expression of the Drug-Metabolizing Enzymes in HepaRG Cells

Two RNA-editing proteins, the adenosine deaminase acting on RNA, ADAR, and ADARB1, broadly regulate gene expression in editing-dependent and editing-independent manners. Previous studies showed that the expression of the drug-metabolizing cytochrome P450s (P450s) and UDP glucuronosyltransferases (UGTs) changes upon knockdown (KD) of ADAR or ADARB1 in different hepatic cell lines. To systematically survey the effects of these two ADARs on the expression of P450s and UGTs, we used small interfering RNA in HepaRG cells and tested the association between the expression of the P450s and ADARs in a liver sample cohort (n = 246). KD of ADAR increased the expression of the CYP3As and CYP2C9 and reduced the expression of the others, whereas KD of ADARB1 reduced the expression of nearly all genes tested. ADAR KD also suppressed the induction of most P450s, whereas ADARB1 KD had mixed effects depending on the inducer/gene combination. P450 expression was positively associated with both ADARs in liver samples, consistent with the KD results. However, after adjusting for the expression of transcription factors (TFs) known to regulate P450 expression, the associations disappeared, indicating that the effects of ADAR or ADARB1 primarily occur through TFs. Moreover, we found that the expression of normally spliced CYP3A5 transcripts is increased in both KDs, indicating a direct effect of the ADARs on promoting the usage of the cryptic splice site generated by CYP3A5*3. Taken together, our results revealed the nonoverlapping regulatory effects of ADAR and ADARB1 and supported their broad roles in controlling the expression of drug-metabolizing enzymes in the liver. SIGNIFICANCE STATEMENT: Here, this study systematically surveyed the roles of ADAR and ADARB1 in both basal and induced expression of drug-metabolizing enzymes and assessed their coexpression in liver samples. This study's results support that ADAR and ADARB1 regulate the expression of the drug-metabolizing enzymes in the liver, suggesting that factors affecting ADAR expression also have the potential to impact drug metabolism.

A pan-tissue survey of mosaic chromosomal alterations in 948 individuals

Genetic mutations accumulate in an organism's body throughout its lifetime. While somatic single-nucleotide variants have been well characterized in the human body, the patterns and consequences of large chromosomal alterations in normal tissues remain largely unknown. Here, we present a pan-tissue survey of mosaic chromosomal alterations (mCAs) in 948 healthy individuals from the Genotype-Tissue Expression project, augmenting RNA-based allelic imbalance estimation with haplotype phasing. We found that approximately a quarter of the individuals carry a clonally-expanded mCA in at least one tissue, with incidence strongly correlated with age. The prevalence and genome-wide patterns of mCAs vary considerably across tissue types, suggesting tissue-specific mutagenic exposure and selection pressures. The mCA landscapes in normal adrenal and pituitary glands resemble those in tumors arising from these tissues, whereas the same is not true for the esophagus and skin. Together, our findings show a widespread age-dependent emergence of mCAs across normal human tissues with intricate connections to tumorigenesis.

Using autopsies to dissect COVID-19 pathogenesis

The emergence of SARS-CoV-2 has resulted in millions of deaths as a result of COVID-19. Suitable models were missing at the beginning of the pandemic, and studies investigating disease pathogenesis relied on patients who had succumbed to COVID-19. Since then, autopsies of patients have substantially contributed to our understanding of the pathogenesis of COVID-19 and associated major organ complications. Here we summarize how autopsies have complemented experimental studies, mainly in animal models, and how they have facilitated critical knowledge of COVID-19 to improve daily clinical practice and develop therapeutic interventions. Employing advanced histopathologic and molecular genetic methods in post-mortem tissues, the COVID-19 pandemic has highlighted the importance of autopsies for virology research and clinical practice in current and emerging infectious diseases.

Identification of novel genetic risk factors of dilated cardiomyopathy: from canine to human

BACKGROUND

Dilated cardiomyopathy (DCM) is a life-threatening heart disease and a common cause of heart failure due to systolic dysfunction and subsequent left or biventricular dilatation. A significant number of cases have a genetic etiology; however, as a complex disease, the exact genetic risk factors are largely unknown, and many patients remain without a molecular diagnosis.

METHODS

We performed GWAS followed by whole-genome, transcriptome, and immunohistochemical analyses in a spontaneously occurring canine model of DCM. Canine gene discovery was followed up in three human DCM cohorts.

RESULTS

Our results revealed two independent additive loci associated with the typical DCM phenotype comprising left ventricular systolic dysfunction and dilatation. We highlight two novel candidate genes, RNF207 and PRKAA2, known for their involvement in cardiac action potentials, energy homeostasis, and morphology. We further illustrate the distinct genetic etiologies underlying the typical DCM phenotype and ventricular premature contractions. Finally, we followed up on the canine discoveries in human DCM patients and discovered candidate variants in our two novel genes.

CONCLUSIONS

Collectively, our study yields insight into the molecular pathophysiology of DCM and provides a large animal model for preclinical studies.

Investigating the potential mechanism of quercetin against cervical cancer

BACKGROUND

Cervical cancer is emerging as a potential target of increased susceptibility to coronavirus disease-2019 (COVID-19), leading to compromised survival rates. Despite this critical link, efficacious anti-cervical cancer/COVID-19 interventions remain limited. Quercetin, known for its efficacy against both cancer and viral infections, holds promise as a therapeutic agent. This study aims to elucidate quercetin's anti-cervical cancer/COVID-19 mechanisms and potential targets.

METHODS

We initiated our investigation with differential gene expression analysis using cervical cancer transcriptome data from The Cancer Genome Atlas (TCGA) and The Genotype-Tissue Expression (GTEx), focusing on intersections with COVID-19-related genes. Network pharmacology was employed to identify the shared targets between cervical cancer/COVID-19 DEGs and quercetin's targets. Subsequently, Cox proportional hazards analyses were employed to establish a risk score based on these genes. Molecular docking techniques were applied to predict quercetin's therapeutic targets and mechanisms for mitigating cervical cancer and COVID-19.

RESULTS

Our findings unveiled 45 potential quercetin targets with anti-cervical cancer/COVID-19 actions. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analyses highlighted significant enrichment in immune pathways and COVID-19-related pathways. A refined risk score model, comprising PLA2G7, TNF, TYK2, F2, and NRP1, effectively stratified cervical cancer patients into distinct risk groups. Importantly, molecular docking analyses illuminated quercetin's remarkable binding affinity to the primary protease of the coronavirus.

CONCLUSIONS

In summation, our study suggests that quercetin holds promise as a potential therapeutic agent for mitigating coronavirus function, specifically through its interaction with the primary protease. This research offers novel insights into exploring COVID-19 susceptibility and enhancing survival in cervical cancer patients.

Characterization of Intercalated Cell Markers KIT and LINC01187 in Chromophobe Renal Cell Carcinoma and Other Renal Neoplasms

Introduction. Chromophobe renal cell carcinoma (chromophobe RCC) is the third major subcategory of renal tumors after clear cell RCC and papillary RCC, accounting for approximately 5% of all RCC subtypes. Other oncocytic neoplasms seen commonly in surgical pathology practice include the eosinophilic variant of chromophobe RCC, renal oncocytoma, and low-grade oncocytic unclassified RCC. Methods. In our recent next-generation sequencing based study, we nominated a lineage-specific novel biomarker LINC01187 (long intergenic non-protein coding RNA 1187) which was found to be enriched in chromophobe RCC. Like KIT (cluster of differentiation 117; CD117), a clinically utilized chromophobe RCC related biomarker, LINC01187 is expressed in intercalated cells of the nephron. In this follow-up study, we performed KIT immunohistochemistry and LINC01187 RNA in situ hybridization (RNA-ISH) on a cohort of chromophobe RCC and other renal neoplasms, characterized the expression patterns, and quantified the expression signals of the two biomarkers in both primary and metastatic settings. Results. LINC01187, in comparison to KIT, exhibits stronger and more uniform expression within tumors while maintaining temporal and spatial consistency. LINC01187 also is devoid of intra-tumoral heterogeneous expression pattern, a phenomenon commonly noted with KIT. Conclusions. LINC01187 expression can augment the currently utilized KIT assay and help facilitate easy microscopic analyses in routine surgical pathology practice.

Integrated immunological analysis of single-cell and bulky tissue transcriptomes reveals the role of interactions between M0 macrophages and naïve CD4+ T cells in the immunosuppressive microenvironment of cervical cancer

In recent decades, the incidence and mortality of cervical cancer have declined in developed countries due to the implementation of screening and vaccination programs. However, cervical cancer remains one of the major culprits of cancer-related deaths in young women. Current studies have found that immune cell-related intercellular communication in the tumor microenvironment has a large impact on the construction of the immunosuppressive microenvironment. In this study, we performed a comprehensive immune analysis on bulk RNA-seq and scRNA-seq data obtained from cervical cancer and revealed that two highly plastic cell populations, M0 macrophages and naïve CD4+ T cells, were significantly correlated with prognosis and clinical phenotypes. Notably, signaling between M0 macrophages and naïve CD4+ T cells as well as intracellular transcription factor activity were significantly altered in the tumor state. Furthermore, we identified overlapping genes between the transcription factor target genes of M0 macrophages or naïve CD4+ T cells and the differentially expressed genes in each type of cell, and these overlapping genes were subsequently subjected to an analysis using the LASSO regression model. Finally, we generated a score index that was significantly associated with the clinical prognosis of cervical cancer. In conclusion, interventions to improve the communication between M0 macrophages and naïve CD4+ T cells may help to improve the immunosuppressive microenvironment of cervical cancer and prevent immune evasion. The relevant molecular mechanisms need to be further validated by experimental and cohort studies.

Relationship between sex biases in gene expression and sex biases in autism and Alzheimer's disease

Sex differences in the brain may play an important role in sex-differential prevalence of neuropsychiatric conditions. In order to understand the transcriptional basis of sex differences, we analyzed multiple, large-scale, human postmortem brain RNA-seq datasets using both within-region and pan-regional frameworks. We find evidence of sex-biased transcription in many autosomal genes, some of which provide evidence for pathways and cell population differences between chromosomally male and female individuals. These analyses also highlight regional differences in the extent of sex-differential gene expression. We observe an increase in specific neuronal transcripts in male brains and an increase in immune and glial function-related transcripts in female brains. Integration with single-cell data suggests this corresponds to sex differences in cellular states rather than cell abundance. Integration with case-control gene expression studies suggests a female molecular predisposition towards Alzheimer's disease, a female-biased disease. Autism, a male-biased diagnosis, does not exhibit a male predisposition pattern in our analysis. Finally, we provide region specific analyses of sex differences in brain gene expression to enable additional studies at the interface of gene expression and diagnostic differences.

Identification of human genetic variants modulating the course of COVID-19 infection with importance in other viral infections

Introduction: COVID-19 has been a major focus of scientific research since early 2020. Due to its societal, economic, and clinical impact worldwide, research efforts aimed, among other questions, to address the effect of host genetics in susceptibility and severity of COVID-19. Methods: We, therefore, performed next-generation sequencing of coding and regulatory regions of 16 human genes, involved in maintenance of the immune system or encoding receptors for viral entry into the host cells, in a subset of 60 COVID-19 patients from the General Hospital Tešanj, Bosnia and Herzegovina, classified into three groups of clinical conditions of different severity ("mild," "moderate," and "severe"). Results: We confirmed that the male sex and older age are risk factors for severe clinical picture and identified 13 variants on seven genes (CD55, IL1B, IL4, IRF7, DDX58, TMPRSS2, and ACE2) with potential functional significance, either as genetic markers of modulated susceptibility to SARS-CoV-2 infection or modifiers of the infection severity. Our results include variants reported for the first time as potentially associated with COVID-19, but further research and larger patient cohorts are required to confirm their effect. Discussion: Such studies, focused on candidate genes and/or variants, have a potential to answer the questions regarding the effect of human genetic makeup on the expected infection outcome. In addition, loci we identified here were previously reported to have clinical significance in other diseases and viral infections, thus confirming a general, broader significance of COVID-19-related research results following the end of the pandemic period.