GENCODE reference annotation for the human and mouse genomes

Hepatic IRE1α-XBP1 signaling promotes GDF15-mediated anorexia and body weight loss in chemotherapy

Platinum-based chemotherapy drugs can lead to the development of anorexia, a detrimental effect on the overall health of cancer patients. However, managing chemotherapy-induced anorexia and subsequent weight loss remains challenging due to limited effective therapeutic strategies. Growth differentiation factor 15 (GDF15) has recently gained significant attention in the context of chemotherapy-induced anorexia. Here, we report that hepatic GDF15 plays a crucial role in regulating body weight in response to chemo drugs cisplatin and doxorubicin. Cisplatin and doxorubicin treatments induce hepatic Gdf15 expression and elevate circulating GDF15 levels, leading to hunger suppression and subsequent weight loss. Mechanistically, selective activation by chemotherapy of hepatic IRE1α-XBP1 pathway of the unfolded protein response (UPR) upregulates Gdf15 expression. Genetic and pharmacological inactivation of IRE1α is sufficient to ameliorate chemotherapy-induced anorexia and body weight loss. These results identify hepatic IRE1α as a molecular driver of GDF15-mediated anorexia and suggest that blocking IRE1α RNase activity offers a therapeutic strategy to alleviate the adverse anorexia effects in chemotherapy.

Human nonunion tissues display differential gene expression in comparison to physiological fracture callus

The healing of bone fractures can become aberrant and lead to nonunions which in turn have a negative impact on patient health. Understanding why a bone fails to normally heal will enable us to make a positive impact in a patient's life. While we have a wealth of molecular data on rodent models of fracture repair, it is not the same with humans. As such, there is still a lack of information regarding the molecular differences between normal physiological repair and nonunions. This study was designed to address this gap in our molecular knowledge of the human repair process by comparing differentially expressed genes (DEGs) between physiological fracture callus and two different nonunion types, hypertrophic (HNU) and oligotrophic (ONU). RNA sequencing data revealed over ∼18,000 genes in each sample. Using the physiological callus as the control and the nonunion samples as the experimental groups, bioinformatic analyses identified 67 and 81 statistically significant DEGs for HNU and ONU, respectively. Out of the 67 DEGs for the HNU, 34 and 33 were up and down-regulated, respectively. Similarly, out of the 81 DEGs for the ONU, 48 and 33 were up and down-regulated, respectively. Additionally, we also identified common genes between the two nonunion samples; 8 (10.8 %) upregulated and 12 (22.2 %) downregulated. We further identified many biological processes, with several statistically significant ones. Some of these were related to muscle and were common between the two nonunion samples. This study represents the first comprehensive attempt to understand the global molecular events occurring in human nonunion biology. With further research, we can perhaps decipher new molecular pathways involved in aberrant healing of human bone fractures that can be therapeutically targeted.

ONECUT2 is a druggable driver of luminal to basal breast cancer plasticity

PURPOSE

Tumor heterogeneity complicates patient treatment and can be due to transitioning of cancer cells across phenotypic cell states. This process is associated with the acquisition of independence from an oncogenic driver, such as the estrogen receptor (ER) in breast cancer (BC), resulting in tumor progression, therapeutic failure and metastatic spread. The transcription factor ONECUT2 (OC2) has been shown to be a master regulator protein of metastatic castration-resistant prostate cancer (mCRPC) tumors that promotes lineage plasticity to a drug-resistant neuroendocrine (NEPC) phenotype. Here, we investigate the role of OC2 in the dynamic conversion between different molecular subtypes in BC.

METHODS

We analyze OC2 expression and clinical significance in BC using public databases and immunohistochemical staining. In vitro, we perform RNA-Seq, RT-qPCR and western-blot after OC2 enforced expression. We also assess cellular effects of OC2 silencing and inhibition with a drug-like small molecule in vitro and in vivo.

RESULTS

OC2 is highly expressed in a substantial subset of hormone receptor negative human BC tumors and tamoxifen-resistant models, and is associated with poor clinical outcome, lymph node metastasis and heightened clinical stage. OC2 inhibits ER expression and activity, suppresses a gene expression program associated with luminal differentiation and activates a basal-like state at the gene expression level. We also show that OC2 is required for cell growth and survival in metastatic BC models and that it can be targeted with a small molecule inhibitor providing a novel therapeutic strategy for patients with OC2 active tumors.

CONCLUSIONS

The transcription factor OC2 is a driver of BC heterogeneity and a potential drug target in distinct cell states within the breast tumors.

Construction of a novel disulfidptosis-related lncRNAs signature for prognosis prediction and anti-tumor immunity in laryngeal squamous cell carcinoma

Disulfidptosis, an innovative type of controlled cellular death linked to metabolic dysfunction, has garnered attention. However, there is limited knowledge regarding the involvement of disulfidptosisrelated lnRNAs (DRlncRNAs) in laryngeal squamous cell carcinoma (LSCC). The objective of our team in this study seeks to establish a DRlncRNAs signature, investigate their prognostic value in LSCC, and explore their associations with immune cell subpopulations, biological signaling pathways, and exploring implications for drug sensitivity. We accessed LSCC patients' RNA-seq data and pertinent clinical data for subsequent further analysis from The Cancer Genome Atlas (TCGA) portal. A literature search was conducted focusing on disulfidptosis-related genes. Pearson correlation coefficients were calculated to identify DRlncRNAs. Differential expression analysis of lncRNAs was performed. Utilizing univariate Cox regression analysis, we identified disulfidptosis-associated prognostic lncRNAs. The LASSO-Cox regression analysis was employed to refine this set of lncRNAs and construct a disulfidptosis-related lncRNAs signature. Various statistical techniques were employed to appraise model predictive performance. Subsequently, risk groups were stratified based on the risk score derived from the DRlncRNAs signature. The superiority of the risk score in prognostication over traditional clinicopathological features in LSCC patients was demonstrated. Evident distinctions emerged between risk groups, particularly in immune cell subpopulations like activated mast cells, eosinophils, and activated NK cells. Finally, the low-risk group demonstrated reduced IC50 values for specific chemotherapeutics like cisplatin and gemcitabine. The in vitro experiments indicated differential behavior of our DRlncRNAs. The DRlncRNAs signature can serve as a robust biomarker with the ability to predict both prognosis and therapeutic responses among patients with LSCC.

TDP1 mutation causing SCAN1 neurodegenerative syndrome hampers the repair of transcriptional DNA double-strand breaks

TDP1 removes transcription-blocking topoisomerase I cleavage complexes (TOP1ccs), and its inactivating H493R mutation causes the neurodegenerative syndrome SCAN1. However, the molecular mechanism underlying the SCAN1 phenotype is unclear. Here, we generate human SCAN1 cell models using CRISPR-Cas9 and show that they accumulate TOP1ccs along with changes in gene expression and genomic distribution of R-loops. SCAN1 cells also accumulate transcriptional DNA double-strand breaks (DSBs) specifically in the G1 cell population due to increased DSB formation and lack of repair, both resulting from abortive removal of transcription-blocking TOP1ccs. Deficient TDP1 activity causes increased DSB production, and the presence of mutated TDP1 protein hampers DSB repair by a TDP2-dependent backup pathway. This study provides powerful models to study TDP1 functions under physiological and pathological conditions and unravels that a gain of function of the mutated TDP1 protein, which prevents DSB repair, rather than a loss of TDP1 activity itself, could contribute to SCAN1 pathogenesis.

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.

Cross-ancestry atlas of gene, isoform, and splicing regulation in the developing human brain

Neuropsychiatric genome-wide association studies (GWASs), including those for autism spectrum disorder and schizophrenia, show strong enrichment for regulatory elements in the developing brain. However, prioritizing risk genes and mechanisms is challenging without a unified regulatory atlas. Across 672 diverse developing human brains, we identified 15,752 genes harboring gene, isoform, and/or splicing quantitative trait loci, mapping 3739 to cellular contexts. Gene expression heritability drops during development, likely reflecting both increasing cellular heterogeneity and the intrinsic properties of neuronal maturation. Isoform-level regulation, particularly in the second trimester, mediated the largest proportion of GWAS heritability. Through colocalization, we prioritized mechanisms for about 60% of GWAS loci across five disorders, exceeding adult brain findings. Finally, we contextualized results within gene and isoform coexpression networks, revealing the comprehensive landscape of transcriptome regulation in development and disease.

A high-resolution map of functional miR-181 response elements in the thymus reveals the role of coding sequence targeting and an alternative seed match

MicroRNAs (miRNAs) are critical post-transcriptional regulators in many biological processes. They act by guiding RNA-induced silencing complexes to miRNA response elements (MREs) in target mRNAs, inducing translational inhibition and/or mRNA degradation. Functional MREs are expected to predominantly occur in the 3' untranslated region and involve perfect base-pairing of the miRNA seed. Here, we generate a high-resolution map of miR-181a/b-1 (miR-181) MREs to define the targeting rules of miR-181 in developing murine T cells. By combining a multi-omics approach with computational high-resolution analyses, we uncover novel miR-181 targets and demonstrate that miR-181 acts predominantly through RNA destabilization. Importantly, we discover an alternative seed match and identify a distinct set of targets with repeat elements in the coding sequence which are targeted by miR-181 and mediate translational inhibition. In conclusion, deep profiling of MREs in primary cells is critical to expand physiologically relevant targetomes and establish context-dependent miRNA targeting rules.

Analysis of 3760 hematologic malignancies reveals rare transcriptomic aberrations of driver genes

BACKGROUND

Rare oncogenic driver events, particularly affecting the expression or splicing of driver genes, are suspected to substantially contribute to the large heterogeneity of hematologic malignancies. However, their identification remains challenging.

METHODS

To address this issue, we generated the largest dataset to date of matched whole genome sequencing and total RNA sequencing of hematologic malignancies from 3760 patients spanning 24 disease entities. Taking advantage of our dataset size, we focused on discovering rare regulatory aberrations. Therefore, we called expression and splicing outliers using an extension of the workflow DROP (Detection of RNA Outliers Pipeline) and AbSplice, a variant effect predictor that identifies genetic variants causing aberrant splicing. We next trained a machine learning model integrating these results to prioritize new candidate disease-specific driver genes.

RESULTS

We found a median of seven expression outlier genes, two splicing outlier genes, and two rare splice-affecting variants per sample. Each category showed significant enrichment for already well-characterized driver genes, with odds ratios exceeding three among genes called in more than five samples. On held-out data, our integrative modeling significantly outperformed modeling based solely on genomic data and revealed promising novel candidate driver genes. Remarkably, we found a truncated form of the low density lipoprotein receptor LRP1B transcript to be aberrantly overexpressed in about half of hairy cell leukemia variant (HCL-V) samples and, to a lesser extent, in closely related B-cell neoplasms. This observation, which was confirmed in an independent cohort, suggests LRP1B as a novel marker for a HCL-V subclass and a yet unreported functional role of LRP1B within these rare entities.

CONCLUSIONS

Altogether, our census of expression and splicing outliers for 24 hematologic malignancy entities and the companion computational workflow constitute unique resources to deepen our understanding of rare oncogenic events in hematologic cancers.

Functional canonical RNAi in mice expressing a truncated Dicer isoform and long dsRNA

Canonical RNA interference (RNAi) is sequence-specific mRNA degradation guided by small interfering RNAs (siRNAs) made by RNase III Dicer from long double-stranded RNA (dsRNA). RNAi roles include gene regulation, antiviral immunity or defense against transposable elements. In mammals, RNAi is constrained by Dicer's adaptation to produce another small RNA class-microRNAs. However, a truncated Dicer isoform (ΔHEL1) supporting RNAi exists in mouse oocytes. A homozygous mutation to express only the truncated ΔHEL1 variant causes dysregulation of microRNAs and perinatal lethality in mice. Here, we report the phenotype and canonical RNAi activity in DicerΔHEL1/wt mice, which are viable, show minimal miRNome changes, but their endogenous siRNA levels are an order of magnitude higher. We show that siRNA production in vivo is limited by available dsRNA, but not by Protein kinase R, a dsRNA sensor of innate immunity. dsRNA expression from a transgene yields sufficient siRNA levels to induce efficient RNAi in heart and muscle. DicerΔHEL1/wt mice with enhanced canonical RNAi offer a platform for examining potential and limits of mammalian RNAi in vivo.

A platform for rapid patient-derived cutaneous neurofibroma organoid establishment and screening

Localized cutaneous neurofibromas (cNFs) are benign tumors that arise in the dermis of patients affected by neurofibromatosis type 1 syndrome. cNFs are benign lesions: they do not undergo malignant transformation or metastasize. Nevertheless, they can cover a significant proportion of the body, with some individuals developing hundreds to thousands of lesions. cNFs can cause pain, itching, and disfigurement resulting in substantial socio-emotional repercussions. Currently, surgery and laser desiccation are the sole treatment options but may result in scarring and potential regrowth from incomplete removal. To identify effective systemic therapies, we introduce an approach to establish and screen cNF organoids. We optimized conditions to support the ex vivo growth of genomically diverse cNFs. Patient-derived cNF organoids closely recapitulate cellular and molecular features of parental tumors as measured by immunohistopathology, methylation, RNA sequencing, and flow cytometry. Our cNF organoid platform enables rapid screening of hundreds of compounds in a patient- and tumor-specific manner.

Endothelial ROBO4 suppresses PTGS2/COX-2 expression and inflammatory diseases

Accumulating evidence suggests that endothelial cells can be useful therapeutic targets. One of the potential targets is an endothelial cell-specific protein, Roundabout4 (ROBO4). ROBO4 has been shown to ameliorate multiple diseases in mice, including infectious diseases and sepsis. However, its mechanisms are not fully understood. In this study, using RNA-seq analysis, we found that ROBO4 downregulates prostaglandin-endoperoxide synthase 2 (PTGS2), which encodes cyclooxygenase-2. Mechanistic analysis reveals that ROBO4 interacts with IQ motif-containing GTPase-activating protein 1 (IQGAP1) and TNF receptor-associated factor 7 (TRAF7), a ubiquitin E3 ligase. In this complex, ROBO4 enhances IQGAP1 ubiquitination through TRAF7, inhibits prolonged RAC1 activation, and decreases PTGS2 expression in inflammatory endothelial cells. In addition, Robo4-deficiency in mice exacerbates PTGS2-associated inflammatory diseases, including arthritis, edema, and pain. Thus, we reveal the molecular mechanism by which ROBO4 suppresses the inflammatory response and vascular hyperpermeability, highlighting its potential as a promising therapeutic target for inflammatory diseases.

Combining DNA and protein alignments to improve genome annotation with LiftOn

As the number and variety of assembled genomes continues to grow, the number of annotated genomes is falling behind, particularly for eukaryotes. DNA-based mapping tools help to address this challenge, but they are only able to transfer annotation between closely-related species. Here we introduce LiftOn, a homology-based software tool that integrates DNA and protein alignments to enhance the accuracy of genome-scale annotation and to allow mapping between relatively distant species. LiftOn's protein-centric algorithm considers both types of alignments, chooses optimal open reading frames, resolves overlapping gene loci, and finds additional gene copies where they exist. LiftOn can reliably transfer annotation between genomes representing members of the same species, as we demonstrate on human, mouse, honey bee, rice, and Arabidopsis thaliana . It can further map annotation effectively across species pairs as far apart as mouse and rat or Drosophila melanogaster and D. erecta .

Identifying drug targets for schizophrenia through gene prioritization

Background

Schizophrenia genome-wide association studies (GWASes) have identified >250 significant loci and prioritized >100 disease-related genes. However, gene prioritization efforts have mostly been restricted to locus-based methods that ignore information from the rest of the genome.

Methods

To more accurately characterize genes involved in schizophrenia etiology, we applied a combination of highly-predictive tools to a published GWAS of 67,390 schizophrenia cases and 94,015 controls. We combined both locus-based methods (fine-mapped coding variants, distance to GWAS signals) and genome-wide methods (PoPS, MAGMA, ultra-rare coding variant burden tests). To validate our findings, we compared them with previous prioritization efforts, known neurodevelopmental genes, and results from the PsyOPS tool.

Results

We prioritized 62 schizophrenia genes, 41 of which were also highlighted by our validation methods. In addition to DRD2 , the principal target of antipsychotics, we prioritized 9 genes that are targeted by approved or investigational drugs. These included drugs targeting glutamatergic receptors ( GRIN2A and GRM3 ), calcium channels ( CACNA1C and CACNB2 ), and GABA B receptor ( GABBR2 ). These also included genes in loci that are shared with an addiction GWAS ( e.g. PDE4B and VRK2 ).

Conclusions

We curated a high-quality list of 62 genes that likely play a role in the development of schizophrenia. Developing or repurposing drugs that target these genes may lead to a new generation of schizophrenia therapies. Rodent models of addiction more closely resemble the human disorder than rodent models of schizophrenia. As such, genes prioritized for both disorders could be explored in rodent addiction models, potentially facilitating drug development.

A Lung Cancer Mouse Model Database

Lung cancer, the leading cause of cancer mortality, exhibits diverse histological subtypes and genetic complexities. Numerous preclinical mouse models have been developed to study lung cancer, but data from these models are disparate, siloed, and difficult to compare in a centralized fashion. Here we established the Lung Cancer Mouse Model Database (LCMMDB), an extensive repository of 1,354 samples from 77 transcriptomic datasets covering 974 samples from genetically engineered mouse models (GEMMs), 368 samples from carcinogen-induced models, and 12 samples from a spontaneous model. Meticulous curation and collaboration with data depositors have produced a robust and comprehensive database, enhancing the fidelity of the genetic landscape it depicts. The LCMMDB aligns 859 tumors from GEMMs with human lung cancer mutations, enabling comparative analysis and revealing a pressing need to broaden the diversity of genetic aberrations modeled in GEMMs. Accompanying this resource, we developed a web application that offers researchers intuitive tools for in-depth gene expression analysis. With standardized reprocessing of gene expression data, the LCMMDB serves as a powerful platform for cross-study comparison and lays the groundwork for future research, aiming to bridge the gap between mouse models and human lung cancer for improved translational relevance.

Regulation of T helper cell differentiation by the interplay between histone modification and chromatin interaction

The development and function of the immune system are controlled by temporospatial gene expression programs, which are regulated by cis-regulatory elements, chromatin structure, and trans-acting factors. In this study, we cataloged the dynamic histone modifications and chromatin interactions at regulatory regions during T helper (Th) cell differentiation. Our data revealed that the H3K4me1 landscape established by MLL4 in naive CD4+ T cells is critical for restructuring the regulatory interaction network and orchestrating gene expression during the early phase of Th differentiation. GATA3 plays a crucial role in further configuring H3K4me1 modification and the chromatin interaction network during Th2 differentiation. Furthermore, we demonstrated that HSS3-anchored chromatin loops function to restrict the activity of the Th2 locus control region (LCR), thus coordinating the expression of Th2 cytokines. Our results provide insights into the mechanisms of how the interplay between histone modifications, chromatin looping, and trans-acting factors contributes to the differentiation of Th cells.

Hyper-physiologic mechanical cues, as an osteoarthritis disease-relevant environmental perturbation, cause a critical shift in set points of methylation at transcriptionally active CpG sites in neo-cartilage organoids

BACKGROUND

Osteoarthritis (OA) is a complex, age-related multifactorial degenerative disease of diarthrodial joints marked by impaired mobility, joint stiffness, pain, and a significant decrease in quality of life. Among other risk factors, such as genetics and age, hyper-physiological mechanical cues are known to play a critical role in the onset and progression of the disease (Guilak in Best Pract Res Clin Rheumatol 25:815-823, 2011). It has been shown that post-mitotic cells, such as articular chondrocytes, heavily rely on methylation at CpG sites to adapt to environmental cues and maintain phenotypic plasticity. However, these long-lasting adaptations may eventually have a negative impact on cellular performance. We hypothesize that hyper-physiologic mechanical loading leads to the accumulation of altered epigenetic markers in articular chondrocytes, resulting in a loss of the tightly regulated balance of gene expression that leads to a dysregulated state characteristic of the OA disease state.

RESULTS

We showed that hyper-physiological loading evokes consistent changes in CpGs associated with expression changes (ML-tCpGs) in ITGA5, CAV1, and CD44, among other genes, which together act in pathways such as anatomical structure morphogenesis (GO:0009653) and response to wound healing (GO:0042060). Moreover, by comparing the ML-tCpGs and their associated pathways to tCpGs in OA pathophysiology (OA-tCpGs), we observed a modest but particular interconnected overlap with notable genes such as CD44 and ITGA5. These genes could indeed represent lasting detrimental changes to the phenotypic state of chondrocytes due to mechanical perturbations that occurred earlier in life. The latter is further suggested by the association between methylation levels of ML-tCpGs mapped to CD44 and OA severity.

CONCLUSION

Our findings confirm that hyper-physiological mechanical cues evoke changes to the methylome-wide landscape of chondrocytes, concomitant with detrimental changes in positional gene expression levels (ML-tCpGs). Since CAV1, ITGA5, and CD44 are subject to such changes and are central and overlapping with OA-tCpGs of primary chondrocytes, we propose that accumulation of hyper-physiological mechanical cues can evoke long-lasting, detrimental changes in set points of gene expression that influence the phenotypic healthy state of chondrocytes. Future studies are necessary to confirm this hypothesis.

Disregarding multimappers leads to biases in the functional assessment of NGS data

BACKGROUND

Standard ChIP-seq and RNA-seq processing pipelines typically disregard sequencing reads whose origin is ambiguous ("multimappers"). This usual practice has potentially important consequences for the functional interpretation of the data: genomic elements belonging to clusters composed of highly similar members are left unexplored.

RESULTS

In particular, disregarding multimappers leads to the underrepresentation in epigenetic studies of recently active transposable elements, such as AluYa5, L1HS and SVAs. Furthermore, this common strategy also has implications for transcriptomic analysis: members of repetitive gene families, such the ones including major histocompatibility complex (MHC) class I and II genes, are under-quantified.

CONCLUSION

Revealing inherent biases that permeate routine tasks such as functional enrichment analysis, our results underscore the urgency of broadly adopting multimapper-aware bioinformatic pipelines -currently restricted to specific contexts or communities- to ensure the reliability of genomic and transcriptomic studies.

Atlantoaxial dislocation in the setting of NMLFS

BACKGROUND

Nablus mask-like facial syndrome (NMFLS) is an extremely rare genetic syndrome characterized by facial dysmorphia as well as developmental delay. In the present report we describe a potential association between non-traumatic atlanto-occipital dislocation and NMFLS in an 11-year old female lacking typical facial features of NMFLS.

CASE DESCRIPTION

An 11-year-old female with autism presented with symptoms of persistent headache and vomiting as well as neck stiffness. Further investigation and CT imaging revealed congenital malformation of the skull base and craniocervical junction with complete posterior subluxation of the left occipital condyle. MRI findings later corroborated the findings on CT.

CONCLUSIONS

The patient was successfully treated with occipitocervical fusion. The findings in this case suggest the possibility that atlanto-occipital instability and generalized occipitocervical may be associated with NMFLS.

Variant-to-function analysis of the childhood obesity chr12q13 locus implicates rs7132908 as a causal variant within the 3' UTR of FAIM2

The ch12q13 locus is among the most significant childhood obesity loci identified in genome-wide association studies. This locus resides in a non-coding region within FAIM2; thus, the underlying causal variant(s) presumably influence disease susceptibility via cis-regulation. We implicated rs7132908 as a putative causal variant by leveraging our in-house 3D genomic data and public domain datasets. Using a luciferase reporter assay, we observed allele-specific cis-regulatory activity of the immediate region harboring rs7132908. We generated isogenic human embryonic stem cell lines homozygous for either rs7132908 allele to assess changes in gene expression and chromatin accessibility throughout a differentiation to hypothalamic neurons, a key cell type known to regulate feeding behavior. The rs7132908 obesity risk allele influenced expression of FAIM2 and other genes and decreased the proportion of neurons produced by differentiation. We have functionally validated rs7132908 as a causal obesity variant that temporally regulates nearby effector genes and influences neurodevelopment and survival.

Early adversity causes sex-specific deficits in perforant pathway connectivity and contextual memory in adolescent mice

BACKGROUND

Early life adversity impairs hippocampal development and function across diverse species. While initial evidence indicated potential variations between males and females, further research is required to validate these observations and better understand the underlying mechanisms contributing to these sex differences. Furthermore, most of the preclinical work in rodents was performed in adult males, with only few studies examining sex differences during adolescence when such differences appear more pronounced. To address these concerns, we investigated the impact of limited bedding (LB), a mouse model of early adversity, on hippocampal development in prepubescent and adolescent male and female mice.

METHODS

RNA sequencing, confocal microscopy, and electron microscopy were used to evaluate the impact of LB and sex on hippocampal development in prepubescent postnatal day 17 (P17) mice. Additional studies were conducted on adolescent mice aged P29-36, which included contextual fear conditioning, retrograde tracing, and ex vivo diffusion magnetic resonance imaging (dMRI).

RESULTS

More severe deficits in axonal innervation and myelination were found in the perforant pathway of prepubescent and adolescent LB males compared to LB female littermates. These sex differences were due to a failure of reelin-positive neurons located in the lateral entorhinal cortex (LEC) to innervate the dorsal hippocampus via the perforant pathway in males, but not LB females, and were strongly correlated with deficits in contextual fear conditioning.

CONCLUSIONS

LB impairs the capacity of reelin-positive cells located in the LEC to project and innervate the dorsal hippocampus in LB males but not female LB littermates. Given the critical role that these projections play in supporting normal hippocampal function, a failure to establish proper connectivity between the LEC and the dorsal hippocampus provides a compelling and novel mechanism to explain the more severe deficits in myelination and contextual freezing found in adolescent LB males.

Disentangling genetic effects on transcriptional and post-transcriptional gene regulation through integrating exon and intron expression QTLs

Expression quantitative trait loci (eQTL) studies typically consider exon expression of genes and discard intronic RNA sequencing reads despite their information on RNA metabolism. Here, we quantify genetic effects on exon and intron levels of genes and their ratio in lymphoblastoid cell lines, revealing thousands of cis-QTLs of each type. While genetic effects are often shared between cis-QTL types, 7814 (47%) are not detected as top cis-QTLs at exon levels. We show that exon levels preferentially capture genetic effects on transcriptional regulation, while exon-intron-ratios better detect those on co- and post-transcriptional processes. Considering all cis-QTL types substantially increases (by 71%) the number of colocalizing variants identified by genome-wide association studies (GWAS). It further allows dissecting the potential gene regulatory processes underlying GWAS associations, suggesting comparable contributions by transcriptional (50%) and co- and post-transcriptional regulation (46%) to complex traits. Overall, integrating intronic RNA sequencing reads in eQTL studies expands our understanding of genetic effects on gene regulatory processes.

Comprehensive characterization of stemness-related lncRNAs in triple-negative breast cancer identified a novel prognostic signature related to treatment outcomes, immune landscape analysis and therapeutic guidance: a silico analysis with in vivo experiments

BACKGROUND

Cancer stem cells (CSCs) and long non-coding RNAs (lncRNAs) are known to play a crucial role in the growth, migration, recurrence, and drug resistance of tumor cells, particularly in triple-negative breast cancer (TNBC). This study aims to investigate stemness-related lncRNAs (SRlncRNAs) as potential prognostic indicators for TNBC patients.

METHODS

Utilizing RNA sequencing data and corresponding clinical information from the TCGA database, and employing Weighted Gene Co-expression Network Analysis (WGCNA) on TNBC mRNAsi sourced from an online database, stemness-related genes (SRGs) and SRlncRNAs were identified. A prognostic model was developed using univariate Cox and LASSO-Cox analysis based on SRlncRNAs. The performance of the model was evaluated using Kaplan-Meier analysis, ROC curves, and ROC-AUC. Additionally, the study delved into the underlying signaling pathways and immune status associated with the divergent prognoses of TNBC patients.

RESULTS

The research identified a signature of six SRlncRNAs (AC245100.6, LINC02511, AC092431.1, FRGCA, EMSLR, and MIR193BHG) for TNBC. Risk scores derived from this signature were found to correlate with the abundance of plasma cells. Furthermore, the nominated chemotherapy drugs for TNBC exhibited considerable variability between different risk score groups. RT-qPCR validation confirmed abnormal expression patterns of these SRlncRNAs in TNBC stem cells, affirming the potential of the SRlncRNAs signature as a prognostic biomarker.

CONCLUSION

The identified signature not only demonstrates predictive power in terms of patient outcomes but also provides insights into the underlying biology, signaling pathways, and immune status associated with TNBC prognosis. The findings suggest the possibility of guiding personalized treatments, including immune checkpoint gene therapy and chemotherapy strategies, based on the risk scores derived from the SRlncRNA signature. Overall, this research contributes valuable knowledge towards advancing precision medicine in the context of TNBC.

Characterizing nascent transcription patterns of PROMPTs, eRNAs, and readthrough transcripts in the ENCODE4 deeply profiled cell lines

Arising as co-products of canonical gene expression, transcription-associated lincRNAs, such as promoter upstream transcripts (PROMPTs), enhancer RNAs (eRNAs), and readthrough (RT) transcripts, are often regarded as byproducts of transcription, although they may be important for the expression of nearby genes. We identified regions of nascent expression of these lincRNA in 16 human cell lines using Bru-seq techniques, and found distinctly regulated patterns of PROMPT, eRNA, and RT transcription using the diverse biochemical approaches in the ENCODE4 deeply profiled cell lines collection. Transcription of these lincRNAs was influenced by sequence-specific features and the local or 3D chromatin landscape. However, these sequence and chromatin features do not describe the full spectrum of lincRNA expression variability we identify, highlighting the complexity of their regulation. This may suggest that transcription-associated lincRNAs are not merely byproducts, but rather that the transcript itself, or the act of its transcription, is important for genomic function.

Heat shock protein 72 supports extracellular matrix production in metastatic mammary tumors

This study identified tumorigenic processes most dependent on murine heat shock protein 72 (HSP72) in the mouse mammary tumor virus-PyMT mammary tumor model, which give rise to spontaneous mammary tumors that exhibit HSP72-dependent metastasis to the lung. RNA-seq expression profiling of Hspa1a/Hspa1b (Hsp72) WT and Hsp72-/- primary mammary tumors discovered significantly lower expression of genes encoding components of the extracellular matrix (ECM) in Hsp72 knockout mammary tumors compared to WT controls. In vitro studies found that genetic or chemical inhibition of HSP72 activity in cultured collagen-expressing human or murine cells also reduces mRNA and protein levels of COL1A1 and several other ECM-encoding genes. In search of a possible mechanistic basis for this relationship, we found HSP72 to support the activation of the tumor growth factor-β-suppressor of mothers against decapentaplegic-3 signaling pathway and evidence of suppressor of mothers against decapentaplegic-3 and HSP72 coprecipitation, suggesting potential complex formation. Human COL1A1 mRNA expression was found to have prognostic value for HER2+ breast tumors over other breast cancer subtypes, suggesting a possible human disease context where targeting HSP72 may have a therapeutic rationale. Analysis of human HER2+ breast tumor gene expression data using a gene set comprising ECM-related gene and protein folding-related gene as an input to the statistical learning algorithm, Galgo, found a subset of these genes that can collectively stratify patients by relapse-free survival, further suggesting a potential interplay between the ECM and protein-folding genes may contribute to tumor progression.

Using genome and transcriptome data from African-ancestry female participants to identify putative breast cancer susceptibility genes

African-ancestry (AA) participants are underrepresented in genetics research. Here, we conducted a transcriptome-wide association study (TWAS) in AA female participants to identify putative breast cancer susceptibility genes. We built genetic models to predict levels of gene expression, exon junction, and 3' UTR alternative polyadenylation using genomic and transcriptomic data generated in normal breast tissues from 150 AA participants and then used these models to perform association analyses using genomic data from 18,034 cases and 22,104 controls. At Bonferroni-corrected P < 0.05, we identified six genes associated with breast cancer risk, including four genes not previously reported (CTD-3080P12.3, EN1, LINC01956 and NUP210L). Most of these genes showed a stronger association with risk of estrogen-receptor (ER) negative or triple-negative than ER-positive breast cancer. We also replicated the associations with 29 genes reported in previous TWAS at P < 0.05 (one-sided), providing further support for an association of these genes with breast cancer risk. Our study sheds new light on the genetic basis of breast cancer and highlights the value of conducting research in AA populations.

Integrated transcriptomics and histopathology approach identifies a subset of rejected donor livers with potential suitability for transplantation

BACKGROUND

Liver transplantation is an effective treatment for liver failure. There is a large unmet demand, even as not all donated livers are transplanted. The clinical selection criteria for donor livers based on histopathological evaluation and liver function tests are variable. We integrated transcriptomics and histopathology to characterize donor liver biopsies obtained at the time of organ recovery. We performed RNA sequencing as well as manual and artificial intelligence-based histopathology (10 accepted and 21 rejected for transplantation).

RESULTS

We identified two transcriptomically distinct rejected subsets (termed rejected-1 and rejected-2), where rejected-2 exhibited a near-complete transcriptomic overlap with the accepted livers, suggesting acceptability from a molecular standpoint. Liver metabolic functional genes were similarly upregulated, and extracellular matrix genes were similarly downregulated in the accepted and rejected-2 groups compared to rejected-1. The transcriptomic pattern of the rejected-2 subset was enriched for a gene expression signature of graft success post-transplantation. Serum AST, ALT, and total bilirubin levels showed similar overlapping patterns. Additional histopathological filtering identified cases with borderline scores and extensive molecular overlap with accepted donor livers.

CONCLUSIONS

Our integrated approach identified a subset of rejected donor livers that are likely suitable for transplantation, demonstrating the potential to expand the pool of transplantable livers.

regionalpcs: improved discovery of DNA methylation associations with complex traits

We have developed the regional principal components (rPCs) method, a novel approach for summarizing gene-level methylation. rPCs address the challenge of deciphering complex epigenetic mechanisms in diseases like Alzheimer's disease (AD). In contrast to traditional averaging, rPCs leverage principal components analysis to capture complex methylation patterns across gene regions. Our method demonstrated a 54% improvement in sensitivity over averaging in simulations, offering a robust framework for identifying subtle epigenetic variations. Applying rPCs to the AD brain methylation data in ROSMAP, combined with cell type deconvolution, we uncovered 838 differentially methylated genes associated with neuritic plaque burden-significantly outperforming conventional methods. Integrating methylation quantitative trait loci (meQTL) with genome-wide association studies (GWAS) identified 17 genes with potential causal roles in AD, including MS4A4A and PICALM. Our approach is available in the Bioconductor package regionalpcs, opening avenues for research and facilitating a deeper understanding of the epigenetic landscape in complex diseases.

The ABC's of Alzheimer risk gene ABCA7

Alzheimer's disease (AD) is a growing problem worldwide. Since ABCA7's identification as a risk gene, it has been extensively researched for its role in the disease. We review its recently characterized structure and what the mechanistic insights teach us about its function. We furthermore provide an overview of identified ABCA7 mutations, their presence in different ancestries and protein domains and how they might cause AD. For ABCA7 PTC variants and a VNTR expansion, haploinsufficiency is proposed as the most likely mode-of-action, although splice events could further influence disease risk. Overall, the need to better understand expression of canonical ABCA7 and its isoforms in disease is indicated. Finally, ABCA7's potential functions in lipid metabolism, phagocytosis, amyloid deposition, and the interplay between these three, is described. To conclude, in this review, we provide a comprehensive overview and discussion about the current knowledge on ABCA7 in AD, and what research questions remain. HIGHLIGHTS: Alzheimer's risk-increasing variants in ABCA7 can be found in up to 7% of AD patients. We review the recently characterized protein structure of ABCA7. We present latest insights in genetics, expression patterns, and functions of ABCA7.

JAK-STAT signaling maintains homeostasis in T cells and macrophages

Immune cells need to sustain a state of constant alertness over a lifetime. Yet, little is known about the regulatory processes that control the fluent and fragile balance that is called homeostasis. Here we demonstrate that JAK-STAT signaling, beyond its role in immune responses, is a major regulator of immune cell homeostasis. We investigated JAK-STAT-mediated transcription and chromatin accessibility across 12 mouse models, including knockouts of all STAT transcription factors and of the TYK2 kinase. Baseline JAK-STAT signaling was detected in CD8+ T cells and macrophages of unperturbed mice-but abrogated in the knockouts and in unstimulated immune cells deprived of their normal tissue context. We observed diverse gene-regulatory programs, including effects of STAT2 and IRF9 that were independent of STAT1. In summary, our large-scale dataset and integrative analysis of JAK-STAT mutant and wild-type mice uncovered a crucial role of JAK-STAT signaling in unstimulated immune cells, where it contributes to a poised epigenetic and transcriptional state and helps prepare these cells for rapid response to immune stimuli.

Oxidative phosphorylation is a pivotal therapeutic target of fibrodysplasia ossificans progressiva

Heterotopic ossification (HO) is a non-physiological bone formation where soft tissue progenitor cells differentiate into chondrogenic cells. In fibrodysplasia ossificans progressiva (FOP), a rare genetic disease characterized by progressive and systemic HO, the Activin A/mutated ACVR1/mTORC1 cascade induces HO in progenitors in muscle tissues. The relevant biological processes aberrantly regulated by activated mTORC1 remain unclear, however. RNA-sequencing analyses revealed the enrichment of genes involved in oxidative phosphorylation (OXPHOS) during Activin A-induced chondrogenesis of mesenchymal stem cells derived from FOP patient-specific induced pluripotent stem cells. Functional analyses showed a metabolic transition from glycolysis to OXPHOS during chondrogenesis, along with increased mitochondrial biogenesis. mTORC1 inhibition by rapamycin suppressed OXPHOS, whereas OXPHOS inhibitor IACS-010759 inhibited cartilage matrix formation in vitro, indicating that OXPHOS is principally involved in mTORC1-induced chondrogenesis. Furthermore, IACS-010759 inhibited the muscle injury-induced enrichment of fibro/adipogenic progenitor genes and HO in transgenic mice carrying the mutated human ACVR1. These data indicated that OXPHOS is a critical downstream mediator of mTORC1 signaling in chondrogenesis and therefore is a potential FOP therapeutic target.

Chromosome evolution screens recapitulate tissue-specific tumor aneuploidy patterns

Whole chromosome and arm-level copy number alterations occur at high frequencies in tumors, but their selective advantages, if any, are poorly understood. Here, utilizing unbiased whole chromosome genetic screens combined with in vitro evolution to generate arm- and subarm-level events, we iteratively selected the fittest karyotypes from aneuploidized human renal and mammary epithelial cells. Proliferation-based karyotype selection in these epithelial lines modeled tissue-specific tumor aneuploidy patterns in patient cohorts in the absence of driver mutations. Hi-C-based translocation mapping revealed that arm-level events usually emerged in multiples of two via centromeric translocations and occurred more frequently in tetraploids than diploids, contributing to the increased diversity in evolving tetraploid populations. Isogenic clonal lineages enabled elucidation of pro-tumorigenic mechanisms associated with common copy number alterations, revealing Notch signaling potentiation as a driver of 1q gain in breast cancer. We propose that intrinsic, tissue-specific proliferative effects underlie tumor copy number patterns in cancer.

Transcription-coupled DNA-protein crosslink repair by CSB and CRL4CSA-mediated degradation

DNA-protein crosslinks (DPCs) arise from enzymatic intermediates, metabolism or chemicals like chemotherapeutics. DPCs are highly cytotoxic as they impede DNA-based processes such as replication, which is counteracted through proteolysis-mediated DPC removal by spartan (SPRTN) or the proteasome. However, whether DPCs affect transcription and how transcription-blocking DPCs are repaired remains largely unknown. Here we show that DPCs severely impede RNA polymerase II-mediated transcription and are preferentially repaired in active genes by transcription-coupled DPC (TC-DPC) repair. TC-DPC repair is initiated by recruiting the transcription-coupled nucleotide excision repair (TC-NER) factors CSB and CSA to DPC-stalled RNA polymerase II. CSA and CSB are indispensable for TC-DPC repair; however, the downstream TC-NER factors UVSSA and XPA are not, a result indicative of a non-canonical TC-NER mechanism. TC-DPC repair functions independently of SPRTN but is mediated by the ubiquitin ligase CRL4CSA and the proteasome. Thus, DPCs in genes are preferentially repaired in a transcription-coupled manner to facilitate unperturbed transcription.

Multiplex, single-cell CRISPRa screening for cell type specific regulatory elements

CRISPR-based gene activation (CRISPRa) is a promising therapeutic approach for gene therapy, upregulating gene expression by targeting promoters or enhancers in a tissue/cell-type specific manner. Here, we describe an experimental framework that combines highly multiplexed perturbations with single-cell RNA sequencing (sc-RNA-seq) to identify cell-type-specific, CRISPRa-responsive cis- regulatory elements and the gene(s) they regulate. Random combinations of many gRNAs are introduced to each of many cells, which are then profiled and partitioned into test and control groups to test for effect(s) of CRISPRa perturbations of both enhancers and promoters on the expression of neighboring genes. Applying this method to a library of 493 gRNAs targeting candidate cis- regulatory elements in both K562 cells and iPSC-derived excitatory neurons, we identify gRNAs capable of specifically upregulating intended target genes and no other neighboring genes within 1 Mb, including gRNAs yielding upregulation of six autism spectrum disorder (ASD) and neurodevelopmental disorder (NDD) risk genes in neurons. A consistent pattern is that the responsiveness of individual enhancers to CRISPRa is restricted by cell type, implying a dependency on either chromatin landscape and/or additional trans- acting factors for successful gene activation. The approach outlined here may facilitate large-scale screens for gRNAs that activate therapeutically relevant genes in a cell type-specific manner.

Worldwide study of the taste of bitter medicines and their modifiers

The bitter taste of medicines hinders patient compliance, but not everyone experiences these difficulties because people worldwide differ in their bitterness perception. To better understand how people from diverse ancestries perceive medicines and taste modifiers, 338 adults, European and recent US and Canada immigrants from Asia, South Asia, and Africa, rated the bitterness intensity of taste solutions on a 100-point generalized visual analog scale and provided a saliva sample for genotyping. The taste solutions were five medicines, tenofovir alafenamide (TAF), moxifloxacin, praziquantel, amodiaquine, and propylthiouracil (PROP), and four other solutions, TAF mixed with sucralose (sweet, reduces bitterness) or 6-methylflavone (tasteless, reduces bitterness), sucralose alone, and sodium chloride alone. Bitterness ratings differed by ancestry for two of the five drugs (amodiaquine and PROP) and for TAF mixed with sucralose. Genetic analysis showed that people with variants in one bitter receptor variant gene (TAS2R38) reported PROP was more bitter than did those with a different variant (p= 7.6e-19) and that people with either an RIMS2 or a THSD4 genotype found sucralose more bitter than did others (p=2.6e-8, p=7.9e-11, resp.). Our findings may help guide the formulation of bad-tasting medicines to meet the needs of those most sensitive to them.

Systematic mapping of TF-mediated cell fate changes by a pooled induction coupled with scRNA-seq and multi-omics approaches

Transcriptional regulation controls cellular functions through interactions between transcription factors (TFs) and their chromosomal targets. However, understanding the fate conversion potential of multiple TFs in an inducible manner remains limited. Here, we introduce iTF-seq as a method for identifying individual TFs that can alter cell fate toward specific lineages at a single-cell level. iTF-seq enables time course monitoring of transcriptome changes, and with biotinylated individual TFs, it provides a multi-omics approach to understanding the mechanisms behind TF-mediated cell fate changes. Our iTF-seq study in mouse embryonic stem cells identified multiple TFs that trigger rapid transcriptome changes indicative of differentiation within a day of induction. Moreover, cells expressing these potent TFs often show a slower cell cycle and increased cell death. Further analysis using bioChIP-seq revealed that GCM1 and OTX2 act as pioneer factors and activators by increasing gene accessibility and activating the expression of lineage specification genes during cell fate conversion. iTF-seq has utility in both mapping cell fate conversion and understanding cell fate conversion mechanisms.

Protocol to process crosslinking and immunoprecipitation data into annotated binding sites

Here, we present a protocol for using Skipper, a pipeline designed to process crosslinking and immunoprecipitation (CLIP) data into annotated binding sites. We describe steps for partitioning annotated transcript regions and fitting data to a beta-binomial model to call windows of enriched binding. From raw CLIP data, we detail how users can map reproducible RNA-binding sites to call enriched windows and perform downstream analysis. This protocol supports optional customizations for different use cases. For complete details on the use and execution of this protocol, please refer to Boyle et al.1.

Common risk alleles for schizophrenia within the major histocompatibility complex predict white matter microstructure

Recent research has highlighted the role of complement genes in shaping the microstructure of the brain during early development, and in contributing to common allele risk for Schizophrenia. We hypothesised that common risk variants for schizophrenia within complement genes will associate with structural changes in white matter microstructure within tracts innervating the frontal lobe. Results showed that risk alleles within the complement gene set, but also intergenic alleles, significantly predict axonal density in white matter tracts connecting frontal cortex with parietal, temporal and occipital cortices. Specifically, risk alleles within the Major Histocompatibility Complex region in chromosome 6 appeared to drive these associations. No significant associations were found for the orientation dispersion index. These results suggest that changes in axonal packing - but not in axonal coherence - determined by common risk alleles within the MHC genomic region - including variants related to the Complement system - appear as a potential neurobiological mechanism for schizophrenia.

Massively parallel screen uncovers many rare 3' UTR variants regulating mRNA abundance of cancer driver genes

Understanding the function of rare non-coding variants represents a significant challenge. Using MapUTR, a screening method, we studied the function of rare 3' UTR variants affecting mRNA abundance post-transcriptionally. Among 17,301 rare gnomAD variants, an average of 24.5% were functional, with 70% in cancer-related genes, many in critical cancer pathways. This observation motivated an interrogation of 11,929 somatic mutations, uncovering 3928 (33%) functional mutations in 155 cancer driver genes. Functional MapUTR variants were enriched in microRNA- or protein-binding sites and may underlie outlier gene expression in tumors. Further, we introduce untranslated tumor mutational burden (uTMB), a metric reflecting the amount of somatic functional MapUTR variants of a tumor and show its potential in predicting patient survival. Through prime editing, we characterized three variants in cancer-relevant genes (MFN2, FOSL2, and IRAK1), demonstrating their cancer-driving potential. Our study elucidates the function of tens of thousands of non-coding variants, nominates non-coding cancer driver mutations, and demonstrates their potential contributions to cancer.

Non-invasive prediction of preeclampsia using the maternal plasma cell-free DNA profile and clinical risk factors

Background

Preeclampsia (PE) is a pregnancy complication defined by new onset hypertension and proteinuria or other maternal organ damage after 20 weeks of gestation. Although non-invasive prenatal testing (NIPT) has been widely used to detect fetal chromosomal abnormalities during pregnancy, its performance in combination with maternal risk factors to screen for PE has not been extensively validated. Our aim was to develop and validate classifiers that predict early- or late-onset PE using the maternal plasma cell-free DNA (cfDNA) profile and clinical risk factors.

Methods

We retrospectively collected and analyzed NIPT data of 2,727 pregnant women aged 24-45 years from four hospitals in China, which had previously been used to screen for fetal aneuploidy at 12 + 0 ~ 22 + 6 weeks of gestation. According to the diagnostic criteria for PE and the time of diagnosis (34 weeks of gestation), a total of 143 early-, 580 late-onset PE samples and 2,004 healthy controls were included. The wilcoxon rank sum test was used to identify the cfDNA profile for PE prediction. The Fisher's exact test and Mann-Whitney U-test were used to compare categorical and continuous variables of clinical risk factors between PE samples and healthy controls, respectively. Machine learning methods were performed to develop and validate PE classifiers based on the cfDNA profile and clinical risk factors.

Results

By using NIPT data to analyze cfDNA coverages in promoter regions, we found the cfDNA profile, which was differential cfDNA coverages in gene promoter regions between PE and healthy controls, could be used to predict early- and late-onset PE. Maternal age, body mass index, parity, past medical histories and method of conception were significantly differential between PE and healthy pregnant women. With a false positive rate of 10%, the classifiers based on the combination of the cfDNA profile and clinical risk factors predicted early- and late-onset PE in four datasets with an average accuracy of 89 and 80% and an average sensitivity of 63 and 48%, respectively.

Conclusion

Incorporating cfDNA profiles in classifiers might reduce performance variations in PE models based only on clinical risk factors, potentially expanding the application of NIPT in PE screening in the future.

Improved multi-ancestry fine-mapping identifies cis-regulatory variants underlying molecular traits and disease risk

Multi-ancestry statistical fine-mapping of cis-molecular quantitative trait loci (cis-molQTL) aims to improve the precision of distinguishing causal cis-molQTLs from tagging variants. However, existing approaches fail to reflect shared genetic architectures. To solve this limitation, we present the Sum of Shared Single Effects (SuShiE) model, which leverages LD heterogeneity to improve fine-mapping precision, infer cross-ancestry effect size correlations, and estimate ancestry-specific expression prediction weights. We apply SuShiE to mRNA expression measured in PBMCs (n=956) and LCLs (n=814) together with plasma protein levels (n=854) from individuals of diverse ancestries in the TOPMed MESA and GENOA studies. We find SuShiE fine-maps cis-molQTLs for 16% more genes compared with baselines while prioritizing fewer variants with greater functional enrichment. SuShiE infers highly consistent cis-molQTL architectures across ancestries on average; however, we also find evidence of heterogeneity at genes with predicted loss-of-function intolerance, suggesting that environmental interactions may partially explain differences in cis-molQTL effect sizes across ancestries. Lastly, we leverage estimated cis-molQTL effect-sizes to perform individual-level TWAS and PWAS on six white blood cell-related traits in AOU Biobank individuals (n=86k), and identify 44 more genes compared with baselines, further highlighting its benefits in identifying genes relevant for complex disease risk. Overall, SuShiE provides new insights into the cis-genetic architecture of molecular traits.

Updating mRNA variants of the human RSK4 gene and their expression in different stressed situations

We determined RNA spectrum of the human RSK4 (hRSK4) gene (also called RPS6KA6) and identified 29 novel mRNA variants derived from alternative splicing, which, plus the NCBI-documented ones and the five we reported previously, totaled 50 hRSK4 RNAs that, by our bioinformatics analyses, encode 35 hRSK4 protein isoforms of 35-762 amino acids. Many of the mRNAs are bicistronic or tricistronic for hRSK4. The NCBI-normalized NM_014496.5 and the protein it encodes are designated herein as the Wt-1 mRNA and protein, respectively, whereas the NM_001330512.1 and the long protein it encodes are designated as the Wt-2 mRNA and protein, respectively. Many of the mRNA variants responded differently to different situations of stress, including serum starvation, a febrile temperature, treatment with ethanol or ethanol-extracted clove buds (an herbal medicine), whereas the same stressed situation often caused quite different alterations among different mRNA variants in different cell lines. Mosifloxacin, an antibiotics and also a functional inhibitor of hRSK4, could inhibit the expression of certain hRSK4 mRNA variants. The hRSK4 gene likely uses alternative splicing as a handy tool to adapt to different stressed situations, and the mRNA and protein multiplicities may partly explain the incongruous literature on its expression and comports.

ERBB3 Overexpression is Enriched in Diverse Patient Populations with Castration-sensitive Prostate Cancer and is Associated with a Unique AR Activity Signature

PURPOSE

Despite successful clinical management of castration-sensitive prostate cancer (CSPC), the 5-year survival rate for men with castration-resistant prostate cancer is only 32%. Combination treatment strategies to prevent disease recurrence are increasing, albeit in biomarker-unselected patients. Identifying a biomarker in CSPC to stratify patients who will progress on standard-of-care therapy could guide therapeutic strategies.

EXPERIMENTAL DESIGN

Targeted deep sequencing was performed for the University of Illinois (UI) cohort (n = 30), and immunostaining was performed on a patient tissue microarray (n = 149). Bioinformatic analyses identified pathways associated with biomarker overexpression (OE) in the UI cohort, consolidated RNA sequencing samples accessed from Database of Genotypes and Phenotypes (n = 664), and GSE209954 (n = 68). Neutralizing antibody patritumab and ectopic HER3 OE were utilized for functional mechanistic experiments.

RESULTS

We identified ERBB3 OE in diverse patient populations with CSPC, where it was associated with advanced disease at diagnosis. Bioinformatic analyses showed a positive correlation between ERBB3 expression and the androgen response pathway despite low dihydrotestosterone and stable expression of androgen receptor (AR) transcript in Black/African American men. At the protein level, HER3 expression was negatively correlated with intraprostatic androgen in Black/African American men. Mechanistically, HER3 promoted enzalutamide resistance in prostate cancer cell line models and HER3-targeted therapy resensitized therapy-resistant prostate cancer cell lines to enzalutamide.

CONCLUSIONS

In diverse patient populations with CSPC, ERBB3 OE was associated with high AR signaling despite low intraprostatic androgen. Mechanistic studies demonstrated a direct link between HER3 and enzalutamide resistance. ERBB3 OE as a biomarker could thus stratify patients for intensification of therapy in castration-sensitive disease, including targeting HER3 directly to improve sensitivity to AR-targeted therapies.

Transcriptomic analysis reveals Cilostazol's role in ameliorating cardiovascular disease: Inhibition of monocyte-to-macrophage differentiation and reduction of endothelial cell reactive oxygen species production

Background

Cardiovascular diseases (CVDs) are the leading global cause of death, with atherosclerosis as the primary cause. Chronic inflammation, endothelial dysfunction, and the role of molecules like nitric oxide and reactive oxygen species are crucial in this context. Our previous research indicated that cilostazol and ginkgo biloba extract could enhance the ability of endothelial cells to dissolve blood clots, but the effects of cilostazol on monocytes remain unexplored.

Method

This study utilized peripheral blood mononuclear cells from 10 healthy donors, treated ex vivo with cilostazol. RNA-sequencing, over-representation analysis, xCell stromal cell analysis, and Gene Set Enrichment Analysis were employed to investigate the gene expression changes and biological pathways affected by cilostazol treatment.

Results

The study identified specific gene sets and pathways that were enriched or reduced in response to cilostazol treatment, providing insights into its effects on monocytes and potential therapeutic applications in CVD. The analysis also revealed the potential impact of cilostazol on the stromal cell compartment, further broadening our understanding of its multifaceted role.

Conclusion

The findings offer a nuanced understanding of the advantages and mechanisms of cilostazol in CVD, uncovering novel therapeutic targets and strategies to enhance the clinical application of cilostazol and contributing to the broader implications of this therapy in cardiovascular health.

Evolution of T cells in the cancer-resistant naked mole-rat

Naked mole-rats (NMRs) are best known for their extreme longevity and cancer resistance, suggesting that their immune system might have evolved to facilitate these phenotypes. Natural killer (NK) and T cells have evolved to detect and destroy cells infected with pathogens and to provide an early response to malignancies. While it is known that NMRs lack NK cells, likely lost during evolution, little is known about their T-cell subsets in terms of the evolution of the genes that regulate their function, their clonotypic diversity, and the thymus where they mature. Here we find, using single-cell transcriptomics, that NMRs have a large circulating population of γδT cells, which in mice and humans mostly reside in peripheral tissues and induce anti-cancer cytotoxicity. Using single-cell-T-cell-receptor sequencing, we find that a cytotoxic γδT-cell subset of NMRs harbors a dominant clonotype, and that their conventional CD8 αβT cells exhibit modest clonotypic diversity. Consistently, perinatal NMR thymuses are considerably smaller than those of mice yet follow similar involution progression. Our findings suggest that NMRs have evolved under a relaxed intracellular pathogenic selective pressure that may have allowed cancer resistance and longevity to become stronger targets of selection to which the immune system has responded by utilizing γδT cells.

GRAS1 non-coding RNA protects against DNA damage and cell death by binding and stabilizing NKAP

Non-coding RNA (ncRNA) gene products are involved in diverse biological processes including splicing, epigenetic regulation, gene expression, proliferation, and metabolism. The biological mechanisms by which ncRNAs contribute to cell survival remain poorly understood. We found that the Growth Regulator Antisense 1 (GRAS1) long non-coding RNA (lncRNA) transcript promotes growth in multiple human cell types by protecting against DNA damage. Knockdown of GRAS1 induced DNA damage and cell death, along with significant expression changes in DNA damage response, intrinsic apoptotic signaling, and cellular response to environmental stimulus genes. Extensive DNA damage occurred after GRAS1 knockdown, with numerous double strand breaks occurring in each cell. The number of cells undergoing apoptosis and with fragmented nuclei increased significantly after GRAS1 knockdown. We used RNA antisense purification and mass spectrometry (RAP-MS) to identify the NF-κB activating protein (NKAP) as a direct protein interaction partner of GRAS1 lncRNA. NKAP protein was degraded after GRAS1 knockdown, in a proteasome-dependent manner. Overexpression of GRAS1 or NKAP mitigated the DNA damage effects of GRAS1 knockdown. In summary, GRAS1 and NKAP directly interact to protect against DNA damage and cell death in multiple human cell lines.

Multi-omic QTL mapping in early developmental tissues reveals phenotypic and temporal complexity of regulatory variants underlying GWAS loci

Most GWAS loci are presumed to affect gene regulation, however, only ∼43% colocalize with expression quantitative trait loci (eQTLs). To address this colocalization gap, we identify eQTLs, chromatin accessibility QTLs (caQTLs), and histone acetylation QTLs (haQTLs) using molecular samples from three early developmental (EDev) tissues. Through colocalization, we annotate 586 GWAS loci for 17 traits by QTL complexity, QTL phenotype, and QTL temporal specificity. We show that GWAS loci are highly enriched for colocalization with complex QTL modules that affect multiple elements (genes and/or peaks). We also demonstrate that caQTLs and haQTLs capture regulatory variations not associated with eQTLs and explain ∼49% of the functionally annotated GWAS loci. Additionally, we show that EDev-unique QTLs are strongly depleted for colocalizing with GWAS loci. By conducting one of the largest multi-omic QTL studies to date, we demonstrate that many GWAS loci exhibit phenotypic complexity and therefore, are missed by traditional eQTL analyses.

Widespread variation in molecular interactions and regulatory properties among transcription factor isoforms

Most human Transcription factors (TFs) genes encode multiple protein isoforms differing in DNA binding domains, effector domains, or other protein regions. The global extent to which this results in functional differences between isoforms remains unknown. Here, we systematically compared 693 isoforms of 246 TF genes, assessing DNA binding, protein binding, transcriptional activation, subcellular localization, and condensate formation. Relative to reference isoforms, two-thirds of alternative TF isoforms exhibit differences in one or more molecular activities, which often could not be predicted from sequence. We observed two primary categories of alternative TF isoforms: "rewirers" and "negative regulators", both of which were associated with differentiation and cancer. Our results support a model wherein the relative expression levels of, and interactions involving, TF isoforms add an understudied layer of complexity to gene regulatory networks, demonstrating the importance of isoform-aware characterization of TF functions and providing a rich resource for further studies.

MMP-9-dependent proteolysis of the histone H3 N-terminal tail: a critical epigenetic step in driving oncogenic transcription and colon tumorigenesis

Matrix metalloproteinase 9 (MMP-9) is a member of the MMP family and has been recently identified as a nuclear protease capable of clipping histone H3 N-terminal tails (H3NT). This MMP-9-dependent H3NT proteolysis is critical for establishing an active state of gene transcription during osteoclast differentiation and melanoma development. However, whether H3NT cleavage by MMP-9 plays a similar role in other cellular events has not been explored. Here, we dissect the functional contribution of MMP-9-dependent H3NT clipping to colonic tumorigenesis by using a combination of genome-wide transcriptome data, ChIP/ChIPac-qPCR, CRISPR/dCas9 gene-targeting system, and in vivo xenograft models. We show that MMP-9 is overexpressed in colon cancer cells and catalyzes H3NT proteolysis to drive transcriptional activation of growth stimulatory genes. Our studies using knockdown and inhibition approaches clearly indicate that MMP-9 mediates transcriptional activation and promotes colonic tumorigenesis in a manner dependent on its protease activity toward H3NT. Remarkably, artificial H3NT proteolysis at target gene promoters with dCAS9-MMP-9 is sufficient for establishing their transcriptional competence in colon cancer cells, underscoring the importance of MMP-9-dependent H3NT proteolysis per se in the transactivation process. Our data establish new functions and mechanisms for MMP-9 in driving the oncogenic transcription program in colon cancer through H3NT proteolysis, and demonstrate how this epigenetic pathway can be exploited as a potential therapeutic target for cancer treatment.

FuncPEP v2.0: An Updated Database of Functional Short Peptides Translated from Non-Coding RNAs

Over the past decade, there have been reports of short novel functional peptides (less than 100 aa in length) translated from so-called non-coding RNAs (ncRNAs) that have been characterized using mass spectrometry (MS) and large-scale proteomics studies. Therefore, understanding the bivalent functions of some ncRNAs as transcripts that encode both functional RNAs and short peptides, which we named ncPEPs, will deepen our understanding of biology and disease. In 2020, we published the first database of functional peptides translated from non-coding RNAs-FuncPEP. Herein, we have performed an update including the newly published ncPEPs from the last 3 years along with the categorization of host ncRNAs. FuncPEP v2.0 contains 152 functional ncPEPs, out of which 40 are novel entries. A PubMed search from August 2020 to July 2023 incorporating specific keywords was performed and screened for publications reporting validated functional peptides derived from ncRNAs. We did not observe a significant increase in newly discovered functional ncPEPs, but a steady increase. The novel identified ncPEPs included in the database were characterized by a wide array of molecular and physiological parameters (i.e., types of host ncRNA, species distribution, chromosomal density, distribution of ncRNA length, identification methods, molecular weight, and functional distribution across humans and other species). We consider that, despite the fact that MS can now easily identify ncPEPs, there still are important limitations in proving their functionality.