AnimalTFDB 3.0: a comprehensive resource for annotation and prediction of animal transcription factors

TFTG: A comprehensive database for human transcription factors and their targets

Transcription factors (TFs) are major contributors to gene transcription, especially in controlling cell-specific gene expression and disease occurrence and development. Uncovering the relationship between TFs and their target genes is critical to understanding the mechanism of action of TFs. With the development of high-throughput sequencing techniques, a large amount of TF-related data has accumulated, which can be used to identify their target genes. In this study, we developed TFTG (Transcription Factor and Target Genes) database (http://tf.liclab.net/TFTG), which aimed to provide a large number of available human TF-target gene resources by multiple strategies, besides performing a comprehensive functional and epigenetic annotations and regulatory analyses of TFs. We identified extensive available TF-target genes by collecting and processing TF-associated ChIP-seq datasets, perturbation RNA-seq datasets and motifs. We also obtained experimentally confirmed relationships between TF and target genes from available resources. Overall, the target genes of TFs were obtained through integrating the relevant data of various TFs as well as fourteen identification strategies. Meanwhile, TFTG was embedded with user-friendly search, analysis, browsing, downloading and visualization functions. TFTG is designed to be a convenient resource for exploring human TF-target gene regulations, which will be useful for most users in the TF and gene expression regulation research.

An association weight matrix identified biological pathways associated with bull fertility traits in a multi-breed population

Using seven indicator traits, we investigated the genetic basis of bull fertility and predicted gene interactions from SNP associations. We used percent normal sperm as the key phenotype for the association weight matrix-partial correlation information theory (AWM-PCIT) approach. Beyond a simple list of candidate genes, AWM-PCIT predicts significant gene interactions and associations for the selected traits. These interactions formed a network of 537 genes: 38 genes were transcription cofactors, and 41 genes were transcription factors. The network displayed two distinct clusters, one with 294 genes and another with 243 genes. The network is enriched in fertility-associated pathways: steroid biosynthesis, p53 signalling, and the pentose phosphate pathway. Enrichment analysis also highlighted gene ontology terms associated with 'regulation of neurotransmitter secretion' and 'chromatin formation'. Our network recapitulates some genes previously implicated in another network built with lower-density genotypes. Sequence-level data also highlights additional candidate genes relevant to bull fertility, such as FOXO4, FOXP3, GATA1, CYP27B1, and EBP. A trio of regulatory genes-KDM5C, LRRK2, and PME-was deemed core to the network because of their overarching connections. This trio probably influences bull fertility through their interaction with genes, both known and unknown as to their role in male fertility. Future studies may target the trio and their target genes to enrich our understanding of male fertility further.

An activity-specificity trade-off encoded in human transcription factors

Transcription factors (TFs) control specificity and activity of gene transcription, but whether a relationship between these two features exists is unclear. Here we provide evidence for an evolutionary trade-off between the activity and specificity in human TFs encoded as submaximal dispersion of aromatic residues in their intrinsically disordered protein regions. We identified approximately 500 human TFs that encode short periodic blocks of aromatic residues in their intrinsically disordered regions, resembling imperfect prion-like sequences. Mutation of periodic aromatic residues reduced transcriptional activity, whereas increasing the aromatic dispersion of multiple human TFs enhanced transcriptional activity and reprogramming efficiency, promoted liquid-liquid phase separation in vitro and more promiscuous DNA binding in cells. Together with recent work on enhancer elements, these results suggest an important evolutionary role of suboptimal features in transcriptional control. We propose that rational engineering of amino acid features that alter phase separation may be a strategy to optimize TF-dependent processes, including cellular reprogramming.

Transcription factor roles in the local adaptation to temperature in the Andean Spiny Toad Rhinella spinulosa

Environmental temperature strongly influences the adaptation dynamics of amphibians, whose limited regulation capabilities render them susceptible to thermal oscillations. A central element of the adaptive strategies is the transcription factors (TFs), which act as master regulators that orchestrate stress responses, enabling species to navigate the fluctuations of their environment skillfully. Our study delves into the intricate relationship between TF expression and thermal adaptation mechanisms in the Rhinella spinulosa populations. We sought to elucidate the dynamic modulations of TF expression in prometamorphic and metamorphic tadpoles that inhabit two thermally contrasting environments (Catarpe and El Tatio Geyser, Chile) and which were exposed to two thermal treatments (25 °C vs. 20 °C). Our findings unravel an intriguing dichotomy in response strategies between these populations. First, results evidence the expression of 1374 transcription factors. Regarding the temperature shift, the Catarpe tadpoles show a multifaceted approach by up-regulating crucial TFs, including fosB, atf7, and the androgen receptor. These dynamic regulatory responses likely underpin the population's ability to navigate thermal fluctuations effectively. In stark contrast, the El Tatio tadpoles exhibit a more targeted response, primarily up-regulating foxc1. This differential expression suggests a distinct focus on specific TFs to mitigate the effects of temperature variations. Our study contributes to understanding the molecular mechanisms governing thermal adaptation responses and highlights the resilience and adaptability of amphibians in the face of ever-changing environmental conditions.

Exploring the oncogenic potential of circSOD2 in clear cell renal cell carcinoma: a novel positive feedback loop

BACKGROUND

Existing studies have found that circular RNAs (circRNAs) act as sponges for micro RNAs (miRNAs) to control downstream genes. However, the specific functionalities and mechanisms of circRNAs in human clear cell renal cell carcinoma (ccRCC) have yet to be thoroughly investigated.

METHODS

Patient cohorts from online databases were used to screen candidate circRNAs, while another cohort from our hospital was obtained for validation. CircSOD2 was identified as a potential oncogenic target, and its relevant characteristics were investigated during ccRCC progression through various assays. A positive feedback loop containing downstream miRNA and its target gene were identified using bioinformatics and validated by luciferase reporter assays, RNA pull-down, and high-throughput sequencing.

RESULTS

CircSOD2 expression was elevated in tumor samples and significantly correlated with overall survival (OS) and the tumor stage of ccRCC patients, which appeared in the enhanced proliferation, invasion, and migration of tumor cells. Through competitive binding to circSOD2, miR-532-3p can promote the expression of PAX5 and the progression of ccRCC, and such regulation can be salvaged by miR-532-3p inhibitor.

CONCLUSION

A novel positive feedback loop, PAX5/circSOD2/miR-532-3p/PAX5 was identified in the study, indicating that the loop may play an important role in the diagnosis and prognostic prediction in ccRCC patients.

Exploring the potential biological significance of KDELR family genes in lung adenocarcinoma

The Lys-Asp-Glu-Leu receptor (KDELR) family genes play critical roles in a variety of biological processes in different tumors. Our study aimed to provide a comprehensive analysis of the potential roles of KDELRs in lung adenocarcinoma (LUAD). Utilizing data from The Cancer Genome Atlas (TCGA) and Gene Expression Omnibus (GEO) database, as well as clinical samples, we conducted a series of analyses and validations using R software tools and various online resources. The results showed that KDELR family genes and proteins were highly expressed and associated with a poor prognosis of LUAD. Promoter hypomethylation and the competing endogenous RNA (ceRNA) network of PCAT6/hsa-miR-326/KDELR1 might be potential causes of aberrant KDELR1 overexpression in LUAD. Three key Transcription factors (TFs) (SPI1, EP300, and MAZ) and a TFs-miRNAs-KDELRs network (involving 11 TFs) might be involved in modulating KDELRs expression abnormalities. Gene Set Enrichment Analysis (GSEA) indicated enrichment of genes highly expressing KDELR1, KDELR2, and KDELR3 in MTORC1_SIGNALING, P53_PATHWAY, and ANGIOGENESIS. Negative correlations between KDELRs expression and CD8 + T cell infiltration, as well as CTLA-4 expression. Our multiple analyses suggested that the KDELRs are important signaling molecules in LUAD. These results provided novel insights for developing prognostic markers and novel therapies of LUAD.

MOI is a comprehensive database collecting processed multi-omics data associated with viral infection

Viral infections pose significant public health challenges, exemplified by the global impact of COVID-19 caused by SARS-CoV-2. Understanding the intricate molecular mechanisms governing virus-host interactions is pivotal for effective intervention strategies. Despite the burgeoning multi-omics data on viral infections, a centralized database elucidating host responses to viruses remains lacking. In response, we have developed a comprehensive database named 'MOI' (available at http://www.fynn-guo.cn/ ), specifically designed to aggregate processed Multi-Omics data related to viral Infections. This meticulously curated database serves as a valuable resource for conducting detailed investigations into virus-host interactions. Leveraging high-throughput sequencing data and metadata from PubMed and Gene Expression Omnibus (GEO), MOI comprises over 3200 viral-infected samples, encompassing human and murine infections. Standardized processing pipelines ensure data integrity, including bulk RNA sequencing (RNA-seq), single-cell RNA-seq (scRNA-seq), Chromatin Immunoprecipitation sequencing (ChIP-seq), and Assay for Transposase-Accessible Chromatin using sequencing (ATAC-seq). MOI offers user-friendly interfaces presenting comprehensive cell marker tables, gene expression data, and epigenetic landscape charts. Analytical tools for DNA sequence conversion, FPKM calculation, differential gene expression, and Gene Ontology (GO)/ Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment enhance data interpretation. Additionally, MOI provides 16 visualization plots for intuitive data exploration. In summary, MOI serves as a valuable repository for researchers investigating virus-host interactions. By centralizing and facilitating access to multi-omics data, MOI aims to advance our understanding of viral pathogenesis and expedite the development of therapeutic interventions.

Ras-dependent activation of BMAL2 regulates hypoxic metabolism in pancreatic cancer

To identify drivers of malignancy in human pancreatic ductal adenocarcinoma (PDAC), we performed regulatory network analysis on a large collection of expression profiles from laser capture microdissected samples of PDAC and benign precursors. We discovered that BMAL2 plays a role in the initiation, progression, post resection survival, and KRAS activity in PDAC. Functional analysis of BMAL2 target genes led us to hypothesize that it plays a role in regulating the response to hypoxia, a critical but poorly understood feature of PDAC physiology. Knockout of BMAL2 in multiple human PDAC cell lines revealed effects on viability and invasion, particularly under hypoxic conditions. Loss of BMAL2 also affected glycolysis and other metabolic processes. We found that BMAL2 directly regulates hypoxia-responsive target genes. We also found that BMAL2 is necessary for the stabilization of HIF1A upon exposure to hypoxia, but destabilizes HIF2A under hypoxia. These data demonstrate that BMAL2 is a master transcriptional regulator of hypoxia responses in PDAC and may serve as a long-sought molecular switch that distinguishes HIF1A- and HIF2A-dependent modes of hypoxic metabolism.

Quantitative proteomics reveals the dynamic proteome landscape of zebrafish embryos during the maternal-to-zygotic transition

Maternal-to-zygotic transition (MZT) is central to early embryogenesis. However, its underlying molecular mechanisms are still not well described. Here, we revealed the expression dynamics of 5,000 proteins across four stages of zebrafish embryos during MZT, representing one of the most systematic surveys of proteome landscape of the zebrafish embryos during MZT. Nearly 700 proteins were differentially expressed and were divided into six clusters according to their expression patterns. The proteome expression profiles accurately reflect the main events that happen during the MZT, i.e., zygotic genome activation (ZGA), clearance of maternal mRNAs, and initiation of cellular differentiation and organogenesis. MZT is modulated by many proteins at multiple levels in a collaborative fashion, i.e., transcription factors, histones, histone-modifying enzymes, RNA helicases, and P-body proteins. Significant discrepancies were discovered between zebrafish proteome and transcriptome profiles during the MZT. The proteome dynamics database will be a valuable resource for bettering our understanding of MZT.

scRank infers drug-responsive cell types from untreated scRNA-seq data using a target-perturbed gene regulatory network

Cells respond divergently to drugs due to the heterogeneity among cell populations. Thus, it is crucial to identify drug-responsive cell populations in order to accurately elucidate the mechanism of drug action, which is still a great challenge. Here, we address this problem with scRank, which employs a target-perturbed gene regulatory network to rank drug-responsive cell populations via in silico drug perturbations using untreated single-cell transcriptomic data. We benchmark scRank on simulated and real datasets, which shows the superior performance of scRank over existing methods. When applied to medulloblastoma and major depressive disorder datasets, scRank identifies drug-responsive cell types that are consistent with the literature. Moreover, scRank accurately uncovers the macrophage subpopulation responsive to tanshinone IIA and its potential targets in myocardial infarction, with experimental validation. In conclusion, scRank enables the inference of drug-responsive cell types using untreated single-cell data, thus providing insights into the cellular-level impacts of therapeutic interventions.

Cholesterol efflux from C1QB-expressing macrophages is associated with resistance to chimeric antigen receptor T cell therapy in primary refractory diffuse large B cell lymphoma

Chimeric antigen receptor T (CAR-T) cell therapy has demonstrated promising efficacy in early trials for relapsed/refractory diffuse large B cell lymphoma (DLBCL). However, its efficacy in treating primary refractory DLBCL has not been comprehensively investigated, and the underlying resistance mechanisms remain unclear. Here, we report the outcomes of a phase I, open-label, single-arm clinical trial of relmacabtagene autoleucel (relma-cel), a CD19-targeted CAR-T cell product, with safety and efficacy as primary endpoints. Among the 12 enrolled patients, 8 experienced grade 4 hematologic toxicity of treatment-emergent adverse event. No grade ≥3 cytokine release syndrome or neurotoxicity occurred. Single-cell RNA sequencing revealed an increase proportion of C1QB-expressing macrophages in patients with progressive disease before CAR-T cell therapy. Cholesterol efflux from M2 macrophages was found to inhibit CAR-T cells cytotoxicity by inducing an immunosuppressive state in CD8+ T cells, leading to their exhaustion. Possible interactions between macrophages and CD8+ T cells, mediating lipid metabolism (AFR1-FAS), immune checkpoint activation, and T cell exhaustion (LGALS9-HAVCR2, CD86-CTLA4, and NECTIN2-TIGIT) were enhanced during disease progression. These findings suggest that cholesterol efflux from macrophages may trigger CD8+ T cell exhaustion, providing a rationale for metabolic reprogramming to counteract CAR-T treatment failure. Chinadrugtrials.org.cn identifier: CTR20200376.

Deciphering the molecular landscape: evolutionary progression from gynecomastia to aggressive male breast cancer

BACKGROUND

Gynecomastia denotes the benign proliferation of glandular breast tissue and stands as a recognized risk factor for male breast cancer. Nonetheless, the underlying carcinogenic mechanisms orchestrating the progression from gynecomastia to cancer remain poorly understood.

METHODS

This study employed single-cell RNA sequencing (scRNA-seq) to meticulously dissect the cellular landscape of gynecomastia and unravel potential associations with male breast cancer at a single-cell resolution. Pseudotime and evolutionary analyses were executed to delineate the distinct features characterizing gynecomastia and male breast cancer. The TCGA database, along with cell-cell communication analysis and immunohistochemistry staining, was harnessed to validate differential gene expression, specifically focusing on CD13.

RESULT

From the copy number variation profiles and evolutionary tree, we inferred shared mutation characteristics (18p+ and 18q+) underpinning both conditions. The developmental trajectory unveiled an intriguing overlap between gynecomastia and malignant epithelial cells. Moreover, the differential gene CD13 emerged as a common denominator in both gynecomastia and male breast cancer when compared with normal mammary tissue. Cell-cell interaction analysis and communication dynamics within the tumor microenvironment spotlighted distinctions between CD13+ and CD13- subsets, with the former exhibiting elevated expression of FGFR1-FGF7.

CONCLUSIONS

Our investigation provides novel insights into the evolutionary progression from gynecomastia to male breast cancer, shedding light on the pivotal role of CD13 in driving this transition. The identification of CD13 as a potential therapeutic target suggests the feasibility of CD13-targeted interventions, specifically tailored for male breast cancer treatment.

Identifying Reproducible Transcription Regulator Coexpression Patterns with Single Cell Transcriptomics

The proliferation of single cell transcriptomics has potentiated our ability to unveil patterns that reflect dynamic cellular processes, rather than cell type compositional effects that emerge from bulk tissue samples. In this study, we leverage a broad collection of single cell RNA-seq data to identify the gene partners whose expression is most coordinated with each human and mouse transcription regulator (TR). We assembled 120 human and 103 mouse scRNA-seq datasets from the literature (>28 million cells), constructing a single cell coexpression network for each. We aimed to understand the consistency of TR coexpression profiles across a broad sampling of biological contexts, rather than examine the preservation of context-specific signals. Our workflow therefore explicitly prioritizes the patterns that are most reproducible across cell types. Towards this goal, we characterize the similarity of each TR's coexpression within and across species. We create single cell coexpression rankings for each TR, demonstrating that this aggregated information recovers literature curated targets on par with ChIP-seq data. We then combine the coexpression and ChIP-seq information to identify candidate regulatory interactions supported across methods and species. Finally, we highlight interactions for the important neural TR ASCL1 to demonstrate how our compiled information can be adopted for community use.

An oncogenic enhancer promotes melanoma progression via regulating ETV4 expression

BACKGROUND

Enhancers are important gene regulatory elements that promote the expression of critical genes in development and disease. Aberrant enhancer can modulate cancer risk and activate oncogenes that lead to the occurrence of various cancers. However, the underlying mechanism of most enhancers in cancer remains unclear. Here, we aim to explore the function and mechanism of a crucial enhancer in melanoma.

METHODS

Multi-omics data were applied to identify an enhancer (enh17) involved in melanoma progression. To evaluate the function of enh17, CRISPR/Cas9 technology were applied to knockout enh17 in melanoma cell line A375. RNA-seq, ChIP-seq and Hi-C data analysis integrated with luciferase reporter assay were performed to identify the potential target gene of enh17. Functional experiments were conducted to further validate the function of the target gene ETV4. Multi-omics data integrated with CUT&Tag sequencing were performed to validate the binding profile of the inferred transcription factor STAT3.

RESULTS

An enhancer, named enh17 here, was found to be aberrantly activated and involved in melanoma progression. CRISPR/Cas9-mediated deletion of enh17 inhibited cell proliferation, migration, and tumor growth of melanoma both in vitro and in vivo. Mechanistically, we identified ETV4 as a target gene regulated by enh17, and functional experiments further support ETV4 as a target gene that is involved in cancer-associated phenotypes. In addition, STAT3 acts as a transcription factor binding with enh17 to regulate the transcription of ETV4.

CONCLUSIONS

Our findings revealed that enh17 plays an oncogenic role and promotes tumor progression in melanoma, and its transcriptional regulatory mechanisms were fully elucidated, which may open a promising window for melanoma prevention and treatment.

Single-cell transcriptomic atlas reveals increased regeneration in diseased human inner ear balance organs

Mammalian inner ear hair cell loss leads to permanent hearing and balance dysfunction. In contrast to the cochlea, vestibular hair cells of the murine utricle have some regenerative capacity. Whether human utricular hair cells regenerate in vivo remains unknown. Here we procured live, mature utricles from organ donors and vestibular schwannoma patients, and present a validated single-cell transcriptomic atlas at unprecedented resolution. We describe markers of 13 sensory and non-sensory cell types, with partial overlap and correlation between transcriptomes of human and mouse hair cells and supporting cells. We further uncover transcriptomes unique to hair cell precursors, which are unexpectedly 14-fold more abundant in vestibular schwannoma utricles, demonstrating the existence of ongoing regeneration in humans. Lastly, supporting cell-to-hair cell trajectory analysis revealed 5 distinct patterns of dynamic gene expression and associated pathways, including Wnt and IGF-1 signaling. Our dataset constitutes a foundational resource, accessible via a web-based interface, serving to advance knowledge of the normal and diseased human inner ear.

O-GlcNAcylation controls pro-fibrotic transcriptional regulatory signaling in myofibroblasts

Tissue injury causes activation of mesenchymal lineage cells into wound-repairing myofibroblasts (MFs), whose uncontrolled activity ultimately leads to fibrosis. Although this process is triggered by deep metabolic and transcriptional reprogramming, functional links between these two key events are not yet understood. Here, we report that the metabolic sensor post-translational modification O-linked β-D-N-acetylglucosaminylation (O-GlcNAcylation) is increased and required for myofibroblastic activation. Inhibition of protein O-GlcNAcylation impairs archetypal myofibloblast cellular activities including extracellular matrix gene expression and collagen secretion/deposition as defined in vitro and using ex vivo and in vivo murine liver injury models. Mechanistically, a multi-omics approach combining proteomic, epigenomic, and transcriptomic data mining revealed that O-GlcNAcylation controls the MF transcriptional program by targeting the transcription factors Basonuclin 2 (BNC2) and TEA domain transcription factor 4 (TEAD4) together with the Yes-associated protein 1 (YAP1) co-activator. Indeed, inhibition of protein O-GlcNAcylation impedes their stability leading to decreased functionality of the BNC2/TEAD4/YAP1 complex towards promoting activation of the MF transcriptional regulatory landscape. We found that this involves O-GlcNAcylation of BNC2 at Thr455 and Ser490 and of TEAD4 at Ser69 and Ser99. Altogether, this study unravels protein O-GlcNAcylation as a key determinant of myofibroblastic activation and identifies its inhibition as an avenue to intervene with fibrogenic processes.

Quantitative estimates of the regulatory influence of long non-coding RNAs on global gene expression variation using TCGA breast cancer transcriptomic data

Long non-coding RNAs (lncRNAs) have received attention in recent years for their regulatory roles in diverse biological contexts including cancer, yet large gaps remain in our understanding of their mechanisms and global maps of their targets. In this work, we investigated a basic unanswered question of lncRNA systems biology: to what extent can gene expression variation across individuals be attributed to lncRNA-driven regulation? To answer this, we analyzed RNA-seq data from a cohort of breast cancer patients, explaining each gene's expression variation using a small set of automatically selected lncRNA regulators. A key aspect of this analysis is that it accounts for confounding effects of transcription factors (TFs) as common regulators of a lncRNA-mRNA pair, to enrich the explained gene expression for lncRNA-mediated regulation. We found that for 16% of analyzed genes, lncRNAs can explain more than 20% of expression variation. We observed 25-50% of the putative regulator lncRNAs to be in 'cis' to, i.e., overlapping or located proximally to the target gene. This led us to quantify the global regulatory impact of such cis-located lncRNAs, which was found to be substantially greater than that of trans-located lncRNAs. Additionally, by including statistical interaction terms involving lncRNA-protein pairs as predictors in our regression models, we identified cases where a lncRNA's regulatory effect depends on the presence of a TF or RNA-binding protein. Finally, we created a high-confidence lncRNA-gene regulatory network whose edges are supported by co-expression as well as a plausible mechanism such as cis-action, protein scaffolding or competing endogenous RNAs. Our work is a first attempt to quantify the extent of gene expression control exerted globally by lncRNAs, especially those located proximally to their regulatory targets, in a specific biological (breast cancer) context. It also marks a first step towards systematic reconstruction of lncRNA regulatory networks, going beyond the current paradigm of co-expression networks, and motivates future analyses assessing the generalizability of our findings to additional biological contexts.

Association between CACNA1D polymorphisms and hypospadias in a southern Chinese population

BACKGROUND

Hypospadias is a congenital genitourinary malformation, with the etiology remaining complex and poorly understood. Despite several genes have been identified to be associated with the risk of hypospadias, current understanding of the susceptibility loci for hypospadias yet remained largely improved. The CACNA1D gene encodes calcium voltage-gated channel subunit alpha 1d and may be involved in androgen signaling. However, the genetic susceptibility of CACNA1D associated with hypospadias has yet been addressed.

OBJECTIVE

To evaluate the association between CACNA1D polymorphisms and the susceptibility to hypospadias.

METHODS

In this study, we accessed the association between two potential regulatory SNPs (rs3774491 and rs898415) within CACNA1D and hypospadias in a cohort of southern Chinese population which comprised of 740 cases and 948 healthy individuals. Both SNP and haplotypic associations were evaluated. Bioinformatic analysis of the regulatory abilities of the CACNA1D SNPs were carried out by utilizing public ChIP-seq and DNase-seq data. The expression of Cacna1d in mouse external genitalia and testis was evaluated by qPCR.

RESULTS

We found that the allele C in rs3774491 and allele G in rs898415 were significantly associated with an increased risk of hypospadias, especially for proximal hypospadias. Further model-based genotypic analyses showed that these association were prominent in additive model and recessive models. Bioinformatic analyses indicated that both SNPs were colocalized with DNase and multiple histone marker across multiple tissues, suggesting the regulatory potentials for these variants. Cacna1d is detectable in both testis and external genitalia of mouse, but the expression level was more prominent in testis than that in external genitalia, suggesting tissue-specific differences in its expression.

CONCLUSION

Our findings provide evidence for CACNA1D as a novel predisposing gene for hypospadias, shedding new light on the genetic basis of malformation of urinary tract. Further investigations are warranted to elucidate the functional implication of CACNA1D underlying the development of hypospadias.

LEVEL OF EVIDENCE

N/A.

Systematical identification of key genes and regulatory genetic variants associated with prognosis of esophageal squamous cell carcinoma

Esophageal squamous cell carcinoma (ESCC) stands as a highly lethal malignancy characterized by pronounced recurrence and metastasis, resulting in a bleak 5-year survival rate. Despite extensive investigations, encompassing genome-wide association studies, the identification of robust prognostic markers has remained elusive. In this study, leveraging four independent data sets comprising 404 ESCC patients, we conducted a systematic analysis to unveil pivotal genes influencing overall survival. our meta-analysis identified 278 genes significantly associated with ESCC prognosis. Further exploration of the prognostic landscape involved an examination of expression quantitative trait loci for these genes, leading to the identification of six tag single nucleotide polymorphisms predictive of overall survival in a cohort of 904 ESCC patients. Notably, functional annotation spotlighted rs11227223, residing in the enhancer region of nuclear paraspeckle assembly transcript 1 (NEAT1), as a crucial variant likely exerting a substantive biological role. Through a series of biochemistry experiments, we conclusively demonstrated that the rs11227223-T allele, indicative of a poorer prognosis, augmented NEAT1 expression. Our results underscore the substantive role of NEAT1 and its regulatory variant in prognostic predictions for ESCC. This comprehensive analysis not only advances our comprehension of ESCC prognosis but also unveils a potential avenue for targeted interventions, offering promise for enhanced clinical outcomes.

Complete characterization of the yak testicular development using accurate full-length transcriptome sequencing

Alternative splicing is a prevalent phenomenon in testicular tissues. Due to the low assembly accuracy of short-read RNA sequencing technology in analyzing post-transcriptional regulatory events, full-length (FL) transcript sequencing is highly demanded to accurately determine FL splicing variants. In this study, we performed FL transcriptome sequencing of testicular tissues from 0.5, 1.5, 2.5, and 4-year-old yaks and 4-year-old cattle-yaks using Oxford Nanopore Technologies. The obtained sequencing data were predicted to have 47,185 open reading frames (ORFs), including 26,630 complete ORFs, detected 7645 fusion transcripts, 15,355 alternative splicing events, 25,798 simple sequence repeats, 7628 transcription factors, and 35,503 long non-coding RNAs. A total of 40,038 novel transcripts were obtained from the sequencing data, and the proportion was almost close to the number of known transcripts identified. Structural analysis and functional annotation of these novel transcripts resulted in the successful annotation of 9568 transcripts, with the highest and lowest annotation numbers in the Nr and KOG databases, respectively. Weighted gene co-expression network analysis revealed the key regulatory pathways and hub genes at various stages of yak testicular development. Our findings enhance our comprehension of transcriptome complexity, contribute to genome annotation refinement, and provide foundational data for further investigations into male sterility in cattle-yaks.

Linc00673-V3 positively regulates autophagy by promoting Smad3-mediated LC3B transcription in NSCLC

Since its first discovery, long noncoding RNA Linc00673 has been linked to carcinogenesis and metastasis of various human cancers. Linc00673 had five transcriptional isoforms and their biological functions remained to be explored. Here we have reported that Linc00673-V3, one of the isoforms of Linc00673, promoted non-small cell lung cancer chemoresistance, and increased Linc00673-V3 expression level was associated with enhanced autophagy. Mechanistically, we discerned the existence of a stem-loop configuration engendered by the 1-100-nt and 2200-2275-nt fragments within Linc00673-V3. This structure inherently interacted with Smad3, thereby impeding its ubiquitination and subsequent degradation orchestrated by E3 ligase STUB1. The accumulation of Smad3 contributed to autophagy via up-regulation of LC3B transcription and ultimately conferred chemoresistance in NSCLC. Our results revealed a novel transcriptional regulation network between Linc00673-V3, Smad3, and LC3B, which provided an important insight into the interplay between autophagy regulation and non-canonical function of Smad3. Furthermore, the results from in vivo experiments suggested Linc00673-V3 targeted antisense oligonucleotide as a promising therapeutic strategy to overcome chemotherapy resistance in NSCLC.

A Comparative Full-Length Transcriptome Analysis Using Oxford Nanopore Technologies (ONT) in Four Tissues of Bovine Origin

The transcriptome complexity and splicing patterns in male and female cattle are ambiguous, presenting a substantial obstacle to genomic selection programs that seek to improve productivity, disease resistance, and reproduction in cattle. A comparative transcriptomic analysis using Oxford Nanopore Technologies (ONT) was conducted in bovine testes (TESTs), ovaries (OVAs), muscles (MUSCs), and livers (LIVs). An average of 5,144,769 full-length reads were obtained from each sample. The TESTs were found to have the greatest number of alternative polyadenylation (APA) events involved in processes such as sperm flagellum development and fertilization in male reproduction. In total, 438 differentially expressed transcripts (DETs) were identified in the LIVs in a comparison of females vs. males, and 214 DETs were identified in the MUSCs between females and males. Additionally, 14,735, 36,347, and 33,885 DETs were detected in MUSC vs. LIV, MUSC vs. TEST, and OVA vs. TEST comparisons, respectively, revealing the complexity of the TEST. Gene Set Enrichment Analysis (GSEA) showed that these DETs were mainly involved in the "spermatogenesis", "flagellated sperm motility", "spermatid development", "reproduction", "reproductive process", and "microtubule-based movement" KEGG pathways. Additional studies are necessary to further characterize the transcriptome in different cell types, developmental stages, and physiological conditions in bovines and ascertain the functions of the novel transcripts.

Unveiling the role of HP1α-HDAC1-STAT1 axis as a therapeutic target for HP1α-positive intrahepatic cholangiocarcinoma

BACKGROUND

Intrahepatic cholangiocarcinoma (ICCA) is a heterogeneous group of malignant tumors characterized by high recurrence rate and poor prognosis. Heterochromatin Protein 1α (HP1α) is one of the most important nonhistone chromosomal proteins involved in transcriptional silencing via heterochromatin formation and structural maintenance. The effect of HP1α on the progression of ICCA remained unclear.

METHODS

The effect on the proliferation of ICCA was detected by experiments in two cell lines and two ICCA mouse models. The interaction between HP1α and Histone Deacetylase 1 (HDAC1) was determined using Electrospray Ionization Mass Spectrometry (ESI-MS) and the binding mechanism was studied using immunoprecipitation assays (co-IP). The target gene was screened out by RNA sequencing (RNA-seq). The occupation of DNA binding proteins and histone modifications were predicted by bioinformatic methods and evaluated by Cleavage Under Targets and Tagmentation (CUT & Tag) and Chromatin immunoprecipitation (ChIP).

RESULTS

HP1α was upregulated in intrahepatic cholangiocarcinoma (ICCA) tissues and regulated the proliferation of ICCA cells by inhibiting the interferon pathway in a Signal Transducer and Activator of Transcription 1 (STAT1)-dependent manner. Mechanistically, STAT1 is transcriptionally regulated by the HP1α-HDAC1 complex directly and epigenetically via promoter binding and changes in different histone modifications, as validated by high-throughput sequencing. Broad-spectrum HDAC inhibitor (HDACi) activates the interferon pathway and inhibits the proliferation of ICCA cells by downregulating HP1α and targeting the heterodimer. Broad-spectrum HDACi plus interferon preparation regimen was found to improve the antiproliferative effects and delay ICCA development in vivo and in vitro, which took advantage of basal activation as well as direct activation of the interferon pathway. HP1α participates in mediating the cellular resistance to both agents.

CONCLUSIONS

HP1α-HDAC1 complex influences interferon pathway activation by directly and epigenetically regulating STAT1 in transcriptional level. The broad-spectrum HDACi plus interferon preparation regimen inhibits ICCA development, providing feasible strategies for ICCA treatment. Targeting the HP1α-HDAC1-STAT1 axis is a possible strategy for treating ICCA, especially HP1α-positive cases.

CUX1 regulates human hematopoietic stem cell chromatin accessibility via the BAF complex

CUX1 is a homeodomain-containing transcription factor that is essential for the development and differentiation of multiple tissues. CUX1 is recurrently mutated or deleted in cancer, particularly in myeloid malignancies. However, the mechanism by which CUX1 regulates gene expression and differentiation remains poorly understood, creating a barrier to understanding the tumor-suppressive functions of CUX1. Here, we demonstrate that CUX1 directs the BAF chromatin remodeling complex to DNA to increase chromatin accessibility in hematopoietic cells. CUX1 preferentially regulates lineage-specific enhancers, and CUX1 target genes are predictive of cell fate in vivo. These data indicate that CUX1 regulates hematopoietic lineage commitment and homeostasis via pioneer factor activity, and CUX1 deficiency disrupts these processes in stem and progenitor cells, facilitating transformation.

Placental Transcriptome Analysis in Connection with Low Litter Birth Weight Phenotype (LBWP) Sows

It is possible to identify sub-populations of sows in every pig herd that consistently give birth to low birth weight (BW) piglets, irrespective of the litter size. A previous study from our group demonstrated that placental development is a main factor affecting the litter birth weight phenotype (LBWP) in sows, thereby impacting the BW of entire litters, but the biological and molecular pathways behind this phenomenon are largely unknown. The aim of this study was to investigate the differential gene expression in placental tissues at day 30 of gestation between low LBWP (LLBWP) vs. high LBWP (HLBWP) sows from a purebred Large White maternal line. Using mRNA sequencing, we found 45 differentially expressed genes (DEGs) in placental tissues of LLBWP and HLBWP sows. Furthermore, (GO) enrichment of upregulated DEGs predicted that there were two biological processes significantly related to cornification and regulation of cell population proliferation. To better understand the molecular interaction between cell proliferation and cornification, we conducted transcriptional factor binding site (TFBS) prediction analysis. The results indicated that a highly significant TFBS was located at the 5' upstream of all four upregulated genes (CDSN, DSG3, KLK14, KRT17), recognized by transcription factors EGR4 and FOSL1. Our findings provide novel insight into how transcriptional regulation of two different biological processes interact in placental tissues of LLBWP sows.

Integrated transcriptomic and metabolomic analysis provides insights into cold tolerance in lettuce (Lactuca sativa L.)

The popular leafy vegetable lettuce (Lactuca sativa L.) is susceptible to cold stress during the growing season, which slows growth rate, causes leaf yellowing and necrosis, and reduced yield and quality. In this study, transcriptomic and metabolomic analyses of two cold-resistant lettuce cultivars (GWAS-W42 and F11) and two cold-sensitive lettuce cultivars (S13K079 and S15K058) were performed to identify the mechanisms involved in the cold response of lettuce. Overall, transcriptome analysis identified 605 differentially expressed genes (DEGs), including significant enrichment of genes involved in the flavonoid and flavonol (CHS, CHI, F3H, FLS, CYP75B1, HCT, etc.) biosynthetic pathways related to oxidation-reduction and catalytic activity. Untargeted metabolomic analysis identified fifteen flavonoid metabolites and 28 other metabolites potentially involved in the response to cold stress; genistein, quercitrin, quercetin derivatives, kaempferol derivatives, luteolin derivatives, apigenin and their derivatives accumulate at higher levels in cold-resistant cultivars. Moreover, MYBs, bHLHs, WRKYs and Dofs also play positive role in the low temperature response, which affected the expression of structural genes contributing to the variation of metabolites between the resistant and sensitive. These results provide valuable evidence that the metabolites and genes involved in the flavonoid biosynthetic pathway play important roles in the response of lettuce to cold stress.

Screening of drug targets for tuberculosis on the basis of transcription factor regulatory network and mRNA sequencing technology

Background

Tuberculosis is a worldwide epidemic disease, posing a serious threat to human health. To find effective drug action targets for Mycobacterium tuberculosis, differentially expressed genes in tuberculosis patients and healthy people were screened by mRNA sequencing in this study. A total of 556 differentially expressed genes in tuberculosis patients and healthy people were screened out by mRNA sequencing technology. 26 transcription factors and 66 corresponding target genes were screened out in the AnimalTFDB 3.0 database, and a transcription factor regulatory network was constructed.

Results

Three key transcription factors (TP53, KLF5 and GATA2) and one key gene (AKT1) were screened as new potential drug targets and diagnostic targets for tuberculosis by MCODE cluster analysis, and the key genes and key transcription factors were verified by RT-PCR. Finally, we constructed the and a key factor and KEGG signaling pathway regulatory network to clarify the possible molecular pathogenesis of tuberculosis.

Conclusion

This study suggested M. tuberculosis may activate the AKT1 gene expression by regulating transcription factors TP53, KLF5, and GATA2, thus activating the B cell receptor signaling pathway to induce the infection and invasion of M. tuberculosis. AKT1, TP53, KLF5, and GATA2 can be used as new potential drug targets for tuberculosis.

Early-stage lung cancer is driven by a transitional cell state dependent on a KRAS-ITGA3-SRC axis

Glycine-12 mutations in the GTPase KRAS (KRASG12) are an initiating event for development of lung adenocarcinoma (LUAD). KRASG12 mutations promote cell-intrinsic rewiring of alveolar type-II progenitor (AT2) cells, but to what extent such changes interplay with lung homeostasis and cell fate pathways is unclear. Here, we generated single-cell RNA-seq (scRNA-seq) profiles from AT2-mesenchyme organoid co-cultures, mice, and stage-IA LUAD patients, identifying conserved regulators of AT2 transcriptional dynamics and defining the impact of KRASG12D mutation with temporal resolution. In AT2WT organoids, we found a transient injury/plasticity state preceding AT2 self-renewal and AT1 differentiation. Early-stage AT2KRAS cells exhibited perturbed gene expression dynamics, most notably retention of the injury/plasticity state. The injury state in AT2KRAS cells of patients, mice, and organoids was distinguishable from AT2WT states via altered receptor expression, including co-expression of ITGA3 and SRC. The combination of clinically relevant KRASG12D and SRC inhibitors impaired AT2KRAS organoid growth. Together, our data show that an injury/plasticity state essential for lung repair is co-opted during AT2 self-renewal and LUAD initiation, suggesting that early-stage LUAD may be susceptible to interventions that target specifically the oncogenic nature of this cell state.

LncRNA BlncAD1 Modulates Bovine Adipogenesis by Binding to MYH10, PI3K/Akt Signaling Pathway, and miR-27a-5p/CDK6 Axis

Research on adipogenesis will help to improve the meat quality of livestock. Long noncoding RNAs (lncRNAs) are involved in mammalian adipogenesis as epigenetic modulators. In this study, we analyzed lncRNA expression during bovine adipogenesis and detected 195 differentially expressed lncRNAs, including lncRNA BlncAD1, which was significantly upregulated in mature bovine adipocytes. Gain- and loss-of-function experiments confirmed that BlncAD1 promoted the proliferation, apoptosis, and differentiation of bovine preadipocytes. RNA pull-down revealed that the nonmuscle myosin 10 (MYH10) is a potential binding protein of BlncAD1. Then, we elucidated that loss of BlncAD1 caused increased ubiquitination of MYH10, which confirmed that BlncAD1 regulates adipogenesis by enhancing the stability of the MYH10 protein. Western blotting was used to demonstrate that BlncAD1 activated the PI3K/Akt signaling pathway. Bioinformatic analysis and dual-luciferase reporter assays indicated that BlncAD1 competitively absorbed miR-27a-5p. The overexpression and interference of miR-27a-5p in bovine preadipocytes displayed that miR-27a-5p inhibited proliferation, apoptosis, and differentiation. Further results suggested that miR-27a-5p targeted the CDK6 gene and that BlncAD1 controlled the proliferation of bovine preadipocytes by modulating the miR-27a-5p/CDK6 axis. This study revealed the complex mechanisms of BlncAD1 underlying bovine adipogenesis for the first time, which would provide useful information for genetics and breeding improvement of Chinese beef cattle.

Causal gene regulatory analysis with RNA velocity reveals an interplay between slow and fast transcription factors

Single-cell expression dynamics, from differentiation trajectories or RNA velocity, have the potential to reveal causal links between transcription factors (TFs) and their target genes in gene regulatory networks (GRNs). However, existing methods either overlook these expression dynamics or necessitate that cells be ordered along a linear pseudotemporal axis, which is incompatible with branching trajectories. We introduce Velorama, an approach to causal GRN inference that represents single-cell differentiation dynamics as a directed acyclic graph of cells, constructed from pseudotime or RNA velocity measurements. Additionally, Velorama enables the estimation of the speed at which TFs influence target genes. Applying Velorama, we uncover evidence that the speed of a TF's interactions is tied to its regulatory function. For human corticogenesis, we find that slow TFs are linked to gliomas, while fast TFs are associated with neuropsychiatric diseases. We expect Velorama to become a critical part of the RNA velocity toolkit for investigating the causal drivers of differentiation and disease.

GRouNdGAN: GRN-guided simulation of single-cell RNA-seq data using causal generative adversarial networks

We introduce GRouNdGAN, a gene regulatory network (GRN)-guided reference-based causal implicit generative model for simulating single-cell RNA-seq data, in silico perturbation experiments, and benchmarking GRN inference methods. Through the imposition of a user-defined GRN in its architecture, GRouNdGAN simulates steady-state and transient-state single-cell datasets where genes are causally expressed under the control of their regulating transcription factors (TFs). Training on six experimental reference datasets, we show that our model captures non-linear TF-gene dependencies and preserves gene identities, cell trajectories, pseudo-time ordering, and technical and biological noise, with no user manipulation and only implicit parameterization. GRouNdGAN can synthesize cells under new conditions to perform in silico TF knockout experiments. Benchmarking various GRN inference algorithms reveals that GRouNdGAN effectively bridges the existing gap between simulated and biological data benchmarks of GRN inference algorithms, providing gold standard ground truth GRNs and realistic cells corresponding to the biological system of interest.

Glioblastoma stem cells deliver ABCB4 transcribed by ATF3 via exosomes conferring glioblastoma resistance to temozolomide

Glioblastoma stem cells (GSCs) play a key role in glioblastoma (GBM) resistance to temozolomide (TMZ) chemotherapy. With the increase in research on the tumour microenvironment, exosomes secreted by GSCs have become a new focus in GBM research. However, the molecular mechanism by which GSCs affect drug resistance in GBM cells via exosomes remains unclear. Using bioinformatics analysis, we identified the specific expression of ABCB4 in GSCs. Subsequently, we established GSC cell lines and used ultracentrifugation to extract secreted exosomes. We conducted in vitro and in vivo investigations to validate the promoting effect of ABCB4 and ABCB4-containing exosomes on TMZ resistance. Finally, to identify the transcription factors regulating the transcription of ABCB4, we performed luciferase assays and chromatin immunoprecipitation-quantitative PCR. Our results indicated that ABCB4 is highly expressed in GSCs. Moreover, high expression of ABCB4 promoted the resistance of GSCs to TMZ. Our study found that GSCs can also transmit their highly expressed ABCB4 to differentiated glioma cells (DGCs) through exosomes, leading to high expression of ABCB4 in these cells and promoting their resistance to TMZ. Mechanistic studies have shown that the overexpression of ABCB4 in GSCs is mediated by the transcription factor ATF3. In conclusion, our results indicate that GSCs can confer resistance to TMZ in GBM by transmitting ABCB4, which is transcribed by ATF3, through exosomes. This mechanism may lead to drug resistance and recurrence of GBM. These findings contribute to a deeper understanding of the mechanisms underlying drug resistance in GBM and provide novel insights into its treatment.

NPAS2, transcriptionally activated by ARRB1, promotes the malignant behaviours of lung adenocarcinoma cells and regulates the reprogramming of glucose metabolism

Lung adenocarcinoma (LUAD) is a serious threat to public health and is accompanied by increased morbidity and mortality worldwide. Neuronal PAS domain protein2 (NPAS2) has been confirmed as an oncogene in LUAD; however, little is known about its molecular mechanism. Here, the expression level of NPAS2 was detected in LUAD cell lines and 16HBE cells. Gain- and loss-of-function experiments were performed. Cell Counting Kit-8, colony formation, flow cytometry, wound-healing and Transwell assays were conducted to assess cell proliferation, apoptosis, migration and invasion, respectively. Reprogramming of glucose metabolism was evaluated via oxygen consumption rate (OCR), complexes activities, lactic production and glucose consumption. The expression of critical proteins was examined by western blot. We demonstrated aberrant upregulation of NPAS2 and β-arrestin-1 (ARRB1) in LUAD cell lines. ARRB1 was found to be a critical transcription factor of NPAS2 with binding sites within the promoter region of NPAS2, thereby causing its transcriptional activation. Functional experiments revealed that NPAS2 depletion significantly inhibited the malignant behaviours of A549 cells by suppressing cell proliferation, migration, invasion and epithelial-mesenchymal transition and promoting cell apoptosis. Meanwhile, NPAS2 depletion increased OCR and activities of complexes (I, II, III and V), and reduced lactic acid production and glucose uptake in A549 cells, indicating that NPAS2 depletion inhibited aerobic glycolysis, accompanied by reduced expression of glycolytic enzymes. However, the changes caused by NPAS2 knockdown were partly restored by ARRB1 overexpression. In conclusion, our study suggests that ARRB1 could transcriptionally activate NPAS2, facilitating malignant activities and glycolysis, and ultimately promoting the progression of LUAD, proving a novel therapeutic strategy for the treatment of LUAD.

CENPA promotes glutamine metabolism and tumor progression by up-regulating SLC38A1 in endometrial cancer

Glutamine addiction is a significant hallmark of metabolic reprogramming in tumors and is crucial to the progression of cancer. Nevertheless, the regulatory mechanisms of glutamine metabolism in endometrial cancer (EC) remains elusive. In this research, we found that elevated expression of CENPA and solute carrier family 38 member 1 (SLC38A1) were firmly associated with worse clinical stage and unfavorable outcomes in EC patients. In addition, ectopic overexpression or silencing of CENPA could either enhance or diminish glutamine metabolism and tumor progression in EC. Mechanistically, CENPA directly regulated the transcriptional activity of the target gene, SLC38A1, leading to enhanced glutamine uptake and metabolism, thereby promoting EC progression. Notably, a prognostic model utilizing the expression levels of CENPA and SLC38A1 genes independently emerged as a prognostic factor for EC. More importantly, CENPA and SLC38A1 were significantly elevated and positively correlated, as well as indicative of poor prognosis in multiple cancers. In brief, our study confirmed that CENPA is a critical transcription factor involved in glutamine metabolism and tumor progression through modulating SLC38A1. This revelation suggests that targeting CENPA could be an appealing therapeutic approach to address pan-cancer glutamine addiction.

Activation of the FOXM1/ASF1B/PRDX3 axis confers hyperproliferative and antioxidative stress reactivity to gastric cancer

Nucleosome assembly during DNA replication is dependent on histone chaperones. Recent studies suggest that dysregulated histone chaperones contribute to cancer progression, including gastric cancer (GC). Further studies are required to explore the prognostic and therapeutic implications of histone chaperones and their mechanisms of action in GC progression. Here we identified histone chaperone ASF1B as a potential biomarker for GC proliferation and prognosis. ASF1B was significantly upregulated in GC, which was associated with poor prognosis. In vitro and in vivo experiments demonstrated that the inhibition of ASF1B suppressed the malignant characteristics of GC, while overexpression of ASF1B had the opposite effect. Mechanistically, transcription factor FOXM1 directly bound to the ASF1B-promoter region, thereby regulating its transcription. Treatment with thiostrepton, a FOXM1 inhibitor, not only suppressed ASF1B expression, but also inhibited GC progression. Furthermore, ASF1B regulated the mitochondrial protein peroxiredoxin 3 (PRDX3) transcription in a FOXM1-dependent manner. The crucial role of ASF1B-regulated PRDX3 in GC cell proliferation and oxidative stress balance was also elucidated. In summary, our study suggests that the FOXM1-ASF1B-PRDX3 axis is a potential therapeutic target for treating GC.

Genetic analysis of PALB2 gene WD40 domain in canine mammary tumour patients

BACKGROUND

DNA repair mechanisms are essential for tumorigenesis and disruption of HR mechanism is an important predisposing factor of human breast cancers (BC). PALB2 is an important part of the HR. There are similarities between canine mammary tumours (CMT) and BCs. As its human counterpart, PALB2 mutations could be a predisposing factor of CMT.

OBJECTIVES

In this study, we aimed to investigate the impacts of PALB2 variants on tumorigenesis and canine mammary tumor (CMT) malignancy.

METHODS

We performed Sanger sequencing to detect germline mutations in the WD40 domain of the canine PALB2 gene in CMT patients. We conducted in silico analysis to investigate the variants, and compared the germline PALB2 mutations in humans that cause breast cancer (BC) with the variants detected in dogs with CMT.

RESULTS

We identified an intronic (c.3096+8C>G) variant, two exonic (p.A1050V and p.R1354R) variants, and a 3' UTR variant (c.4071T>C). Of these, p.R1354R and c.4071T>C novel variants were identified for the first time in this study. We found that the p.A1050V mutation had a significant effect. However, we could not determine sufficient similarity due to the differences in nucleotide/amino acid sequences between two species. Nonetheless, possible variants of human sequences in the exact location as their dog counterparts are associated with several cancer types, implying that the variants could be crucial for tumorigenesis in dogs. Our results did not show any effect of the variants on tumor malignancy.

CONCLUSIONS

The current project is the first study investigating the relationship between the PALB2 gene WD40 domain and CMTs. Our findings will contribute to a better understanding of the pathogenic mechanism of the PALB2 gene in CMTs. In humans, variant positions in canines have been linked to cancer-related phenotypes such as familial BC, endometrial tumor, and hereditary cancer predisposition syndrome. The results of bioinformatics analyses should be investigated through functional tests or case-control studies.

Sexually dimorphic control of affective state processing and empathic behaviors

Recognizing the affective states of social counterparts and responding appropriately fosters successful social interactions. However, little is known about how the affective states are expressed and perceived and how they influence social decisions. Here, we show that male and female mice emit distinct olfactory cues after experiencing distress. These cues activate distinct neural circuits in the piriform cortex (PiC) and evoke sexually dimorphic empathic behaviors in observers. Specifically, the PiC → PrL pathway is activated in female observers, inducing a social preference for the distressed counterpart. Conversely, the PiC → MeA pathway is activated in male observers, evoking excessive self-grooming behaviors. These pathways originate from non-overlapping PiC neuron populations with distinct gene expression signatures regulated by transcription factors and sex hormones. Our study unveils how internal states of social counterparts are processed through sexually dimorphic mechanisms at the molecular, cellular, and circuit levels and offers insights into the neural mechanisms underpinning sex differences in higher brain functions.

In search of the ratio of miRNA expression as robust biomarkers for constructing stable diagnostic models among multi-center data

MicroRNAs (miRNAs) are promising biomarkers for the early detection of disease, and many miRNA-based diagnostic models have been constructed to distinguish patients and healthy individuals. To thoroughly utilize the miRNA-profiling data across different sequencing platforms or multiple centers, the models accounting the batch effects were demanded for the generalization of medical application. We conducted transcription factor (TF)-mediated miRNA-miRNA interaction network analysis and adopted the within-sample expression ratios of miRNA pairs as predictive markers. The ratio of the expression values between each miRNA pair turned out to be stable across multiple data sources. A genetic algorithm-based classifier was constructed to quantify risk scores of the probability of disease and discriminate disease states from normal states in discovery, with a validation dataset for COVID-19, renal cell carcinoma, and lung adenocarcinoma. The predictive models based on the expression ratio of interacting miRNA pairs demonstrated good performances in the discovery and validation datasets, and the classifier may be used accurately for the early detection of disease.

Identification of different myofiber types in pigs muscles and construction of regulatory networks

BACKGROUND

Skeletal muscle is composed of muscle fibers with different physiological characteristics, which plays an important role in regulating skeletal muscle metabolism, movement and body homeostasis. The type of skeletal muscle fiber directly affects meat quality. However, the transcriptome and gene interactions between different types of muscle fibers are not well understood.

RESULTS

In this paper, we selected 180-days-old Large White pigs and found that longissimus dorsi (LD) muscle was dominated by fast-fermenting myofibrils and soleus (SOL) muscle was dominated by slow-oxidizing myofibrils by frozen sections and related mRNA and protein assays. Here, we selected LD muscle and SOL muscle for transcriptomic sequencing, and identified 312 differentially expressed mRNA (DEmRs), 30 differentially expressed miRNA (DEmiRs), 183 differentially expressed lncRNA (DElRs), and 3417 differentially expressed circRNA (DEcRs). The ceRNA network included ssc-miR-378, ssc-miR-378b-3p, ssc-miR-24-3p, XR_308817, XR_308823, SMIM8, MAVS and FOS as multiple core nodes that play important roles in muscle development. Moreover, we found that different members of the miR-10 family expressed differently in oxidized and glycolytic muscle fibers, among which miR-10a-5p was highly expressed in glycolytic muscle fibers (LD) and could target MYBPH gene mRNA. Therefore, we speculate that miR-10a-5p may be involved in the transformation of muscle fiber types by targeting the MYHBP gene. In addition, PPI analysis of differentially expressed mRNA genes showed that ACTC1, ACTG2 and ACTN2 gene had the highest node degree, suggesting that this gene may play a key role in the regulatory network of muscle fiber type determination.

CONCLUSIONS

We can conclude that these genes play a key role in regulating muscle fiber type transformation. Our study provides transcriptomic profiles and ceRNA interaction networks for different muscle fiber types in pigs, providing reference for the transformation of pig muscle fiber types and the improvement of meat quality.

Single-nucleus transcriptomics and chromatin accessibility analysis of musk gland development in Chinese forest musk deer (Moschus berezovskii)

Musk secreted by male forest musk deer (Moschus berezovskii) musk glands is an invaluable component of medicine and perfume. Musk secretion depends on musk gland maturation; however, the mechanism of its development remains elusive. Herein, using single cell multiome ATAC + gene expression coupled with several bioinformatic analyses, a dynamic transcriptional cell atlas of musk gland development was revealed, and key genes and transcription factors affecting its development were determined. Twelve cell types, including two different types of acinar cells (Clusters 0 and 10) were identified. Single-nucleus RNA and single-nucleus ATAC sequencing analyses revealed that seven core target genes associated with musk secretion (Hsd17b2, Acacb, Lss, Vapa, Aldh16a1, Aldh7a1, and Sqle) were regulated by 12 core transcription factors (FOXO1, CUX2, RORA, RUNX1, KLF6, MGA, NFIC, FOXO3, ETV5, NR3C1, HSF4, and MITF) during the development of Cluster 0 acinar cells. Kyoto Encyclopedia of Genes and Genomes enrichment showed significant changes in the pathways associated with musk secretion during acinar cell development. Gene set variation analysis also revealed that certain pathways associated with musk secretion were enriched in 6-year-old acinar cells. A gene co-expression network was constructed during acinar cell development to provide a precise understanding of the connections between transcription factors, genes, and pathways. Finally, intercellular communication analysis showed that intercellular communication is involved in musk gland development. This study provides crucial insights into the changes and key factors underlying musk gland development, which serve as valuable resources for studying musk secretion mechanisms and promoting the protection of this endangered species.

SRRM2 splicing factor modulates cell fate in early development

Embryo development is an orchestrated process that relies on tight regulation of gene expression to guide cell differentiation and fate decisions. The Srrm2 splicing factor has recently been implicated in developmental disorders and diseases, but its role in early mammalian development remains unexplored. Here, we show that Srrm2 dosage is critical for maintaining embryonic stem cell pluripotency and cell identity. Srrm2 heterozygosity promotes loss of stemness, characterised by the coexistence of cells expressing naive and formative pluripotency markers, together with extensive changes in gene expression, including genes regulated by serum-response transcription factor (SRF) and differentiation-related genes. Depletion of Srrm2 by RNA interference in embryonic stem cells shows that the earliest effects of Srrm2 heterozygosity are specific alternative splicing events on a small number of genes, followed by expression changes in metabolism and differentiation-related genes. Our findings unveil molecular and cellular roles of Srrm2 in stemness and lineage commitment, shedding light on the roles of splicing regulators in early embryogenesis, developmental diseases and tumorigenesis.

Construction of gastric cancer prognostic signature based on the E26 transcription factor and the identification of novel oncogene ELK3

The aim of the present study was to investigate the function of 29 E26 (ETS) transcription factor families in gastric cancer (GC) and determine their association with prognosis. Our analysis of the expression of the ETS family revealed that 28 genes were dysregulated in GC, and that their expression was associated with multiple clinicopathological features (P<0.05). Based on the expression signature of the ETS family, consensus clustering was performed to generate two gastric cancer subtypes. These subtypes exhibited differences in overall survival (OS, P = 0.161), disease-free survival (DFS, P<0.05) and GC grade (P<0.01). Functional enrichment analysis of the target genes associated with the ETS family indicated that these genes primarily contribute to functions that facilitate tumor progression. A systematic statistical analysis was used to construct a prognostic model related to OS and DFS in association with the ETS family. This model demonstrated that the maximum area under the curve (AUC) values for predicting OS and DFS were 0.729 and 0.670, respectively, establishing ETS as an independent prognostic factor for GC Furthermore, a nomogram was created from the prognostic signature, and its predictive accuracy was confirmed by a calibration curve. Finally, the expression and prognostic significance of the six genes comprising the model were also examined. Among these, ELK3 was found to be significantly overexpressed in GC clinical samples. Subsequent in vitro and in vivo studies verified that ELK3 regulates GC proliferation and metastasis, highlighting its potential as a therapeutic target for gastric cancer.

Annelid adult cell type diversity and their pluripotent cellular origins

Many annelids can regenerate missing body parts or reproduce asexually, generating all cell types in adult stages. However, the putative adult stem cell populations involved in these processes, and the diversity of cell types generated by them, are still unknown. To address this, we recover 75,218 single cell transcriptomes of the highly regenerative and asexually-reproducing annelid Pristina leidyi. Our results uncover a rich cell type diversity including annelid specific types as well as novel types. Moreover, we characterise transcription factors and gene networks that are expressed specifically in these populations. Finally, we uncover a broadly abundant cluster of putative stem cells with a pluripotent signature. This population expresses well-known stem cell markers such as vasa, piwi and nanos homologues, but also shows heterogeneous expression of differentiated cell markers and their transcription factors. We find conserved expression of pluripotency regulators, including multiple chromatin remodelling and epigenetic factors, in piwi+ cells. Finally, lineage reconstruction analyses reveal computational differentiation trajectories from piwi+ cells to diverse adult types. Our data reveal the cell type diversity of adult annelids by single cell transcriptomics and suggest that a piwi+ cell population with a pluripotent stem cell signature is associated with adult cell type differentiation.

Polyploidy Promotes Hypertranscription, Apoptosis Resistance, and Ciliogenesis in Cancer Cells and Mesenchymal Stem Cells of Various Origins: Comparative Transcriptome In Silico Study

Mesenchymal stem cells (MSC) attract an increasing amount of attention due to their unique therapeutic properties. Yet, MSC can undergo undesirable genetic and epigenetic changes during their propagation in vitro. In this study, we investigated whether polyploidy can compromise MSC oncological safety and therapeutic properties. For this purpose, we compared the impact of polyploidy on the transcriptome of cancer cells and MSC of various origins (bone marrow, placenta, and heart). First, we identified genes that are consistently ploidy-induced or ploidy-repressed through all comparisons. Then, we selected the master regulators using the protein interaction enrichment analysis (PIEA). The obtained ploidy-related gene signatures were verified using the data gained from polyploid and diploid populations of early cardiomyocytes (CARD) originating from iPSC. The multistep bioinformatic analysis applied to the cancer cells, MSC, and CARD indicated that polyploidy plays a pivotal role in driving the cell into hypertranscription. It was evident from the upregulation of gene modules implicated in housekeeping functions, stemness, unicellularity, DNA repair, and chromatin opening by means of histone acetylation operating via DNA damage associated with the NUA4/TIP60 complex. These features were complemented by the activation of the pathways implicated in centrosome maintenance and ciliogenesis and by the impairment of the pathways related to apoptosis, the circadian clock, and immunity. Overall, our findings suggest that, although polyploidy does not induce oncologic transformation of MSC, it might compromise their therapeutic properties because of global epigenetic changes and alterations in fundamental biological processes. The obtained results can contribute to the development and implementation of approaches enhancing the therapeutic properties of MSC by removing polyploid cells from the cell population.

Intergenic sequences harboring potential enhancer elements contribute to Axenfeld-Rieger syndrome by regulating PITX2

Recent studies have uncovered that noncoding sequence variants may relate to Axenfeld-Rieger syndrome (ARS), a rare developmental anomaly with genetic heterogeneity. However, how these genomic regions are functionally and structurally associated with ARS is still unclear. In this study, we performed genome-wide linkage analysis and whole-genome sequencing in a Chinese family with ARS and identified a heterozygous deletion of about 570 kb (termed LOH-1) in the intergenic sequence between paired-like homeodomain transcription factor 2 (PITX2) and family with sequence similarity 241 member A. Knockout of LOH-1 homologous sequences caused ARS phenotypes in mice. RNA-Seq and real-time quantitative PCR revealed a significant reduction in Pitx2 gene expression in LOH-1-/- mice, while forkhead box C1 expression remained unchanged. ChIP-Seq and bioinformatics analysis identified a potential enhancer region (LOH-E1) within LOH-1. Deletion of LOH-E1 led to a substantial downregulation of the PITX2 gene. Mechanistically, we found a sequence (hg38 chr4:111,399,594-111,399,691) that is on LOH-E1 could regulate PITX2 by binding to RAD21, a critical component of the cohesin complex. Knockdown of RAD21 resulted in reduced PITX2 expression. Collectively, our findings indicate that a potential enhancer sequence that is within LOH-1 may regulate PITX2 expression remotely through cohesin-mediated loop domains, leading to ARS when absent.

Fructose improves titanium dioxide nanoparticles induced alterations in developmental competence of mouse oocytes

AIMS

We investigated the effects of intraperitoneal injections of titanium dioxide nanoparticles (TiO2 NPs, 100 mg/kg) for 5 consecutive days on the developmental competence of murine oocytes. Furthermore, study the effects of TiO2 NPs on antioxidant and oxidative stress biomarkers, as well as their effects on expression of apoptotic and hypoxia inducing factor-1α (HIF1A) protein translation. Moreover, the possible ameliorating effects of intraperitoneal injections of fructose (2.75 mM/ml) was examined.

MATERIALS AND METHODS

Thirty sexually mature (8-12 weeks old; ~ 25 g body weight) female mice were used for the current study. The female mice were assigned randomly to three treatment groups: Group1 (G1) mice were injected intraperitoneal (ip) with deionized water for 5 consecutive days; Group 2 (G2) mice were injected ip with TiO2 NPs (100 mg/kg BW) for 5 consecutive days; Group 3 (G3) mice were injected ip with TiO2 NPs (100 mg/kg BW + fructose (2.75 mM) for 5 consecutive days.

RESULTS

Nano-titanium significantly decreased expression of GSH, GPx, and NO, expression of MDA and TAC increased. The rates of MI, MII, GVBD and degenerated oocytes were significantly less for nano-titanium treated mice, but the rate of activated oocytes was significantly greater than those in control oocytes. TiO2 NPs significantly increased expression of apoptotic genes (BAX, Caspase 3 and P53) and HIF1A. Intraperitoneal injection of fructose (2.75 mM/kg) significantly alleviated the detrimental effects of TiO2 NPs. Transmission electron microscopy indicated that fructose mitigated adverse effects of TiO2 NPs to alter the cell surface of murine oocytes.

CONCLUSION

Results of this study suggest that the i/p infusion of fructose for consecutive 5 days enhances development of murine oocytes and decreases toxic effects of TiO2 NPs through positive effects on oxidative and antioxidant biomarkers in cumulus-oocyte complexes and effects to inhibit TiO2-induced increases in expression of apoptotic and hypoxia inducing factors.

Evolved changes in phenotype across skeletal muscles in deer mice native to high altitude

The cold and hypoxic conditions at high altitude necessitate high metabolic O2 demands to support thermogenesis while hypoxia reduces O2 availability. Skeletal muscles play key roles in thermogenesis, but our appreciation of muscle plasticity and adaptation at high altitude has been hindered by past emphasis on only a small number of muscles. We examined this issue in deer mice (Peromyscus maniculatus). Mice derived from both high-altitude and low-altitude populations were born and raised in captivity and then acclimated as adults to normoxia or hypobaric hypoxia (12 kPa O2 for 6-8 wk). Maximal activities of citrate synthase (CS), cytochrome c oxidase (COX), β-hydroxyacyl-CoA dehydrogenase (HOAD), hexokinase (HK), pyruvate kinase (PK), and lactate dehydrogenase (LDH) were measured in 20 muscles involved in shivering, locomotion, body posture, ventilation, and mastication. Principal components analysis revealed an overall difference in muscle phenotype between populations but no effect of hypoxia acclimation. High-altitude mice had greater activities of mitochondrial enzymes and/or lower activities of PK or LDH across many (but not all) respiratory, limb, core and mastication muscles compared with low-altitude mice. In contrast, chronic hypoxia had very few effects across muscles. Further examination of CS in the gastrocnemius showed that population differences in enzyme activity stemmed from differences in protein abundance and mRNA expression but not from population differences in CS amino acid sequence. Overall, our results suggest that evolved increases in oxidative capacity across many skeletal muscles, at least partially driven by differences in transcriptional regulation, may contribute to high-altitude adaptation in deer mice.NEW & NOTEWORTHY Most previous studies of muscle plasticity and adaptation in high-altitude environments have focused on a very limited number of skeletal muscles. Comparing high-altitude versus low-altitude populations of deer mice, we show that a large number of muscles involved in shivering, locomotion, body posture, ventilation, and mastication exhibit greater mitochondrial enzyme activities in the high-altitude population. Therefore, evolved increases in mitochondrial oxidative capacity across skeletal muscles contribute to high-altitude adaptation.

A role for Retinoblastoma 1 in hindbrain morphogenesis by regulating GBX family

The hindbrain, which develops from the anterior end of the neural tube expansion, can differentiate into the metencephalon and myelencephalon, with varying sizes and functions. The midbrain-hindbrain boundary (MHB) and hindbrain myelencephalon/ventral midline (HMVM) are known to be the source of the progenitors for the anterior hindbrain and myelencephalon, respectively. However, the molecular networks regulating hindbrain morphogenesis in these structures remain unclear. In this study, we show that retinoblastoma 1 (rb1) is highly expressed at the MHB and HMVM in zebrafish. Knocking out rb1 in mice and zebrafish results in an enlarged hindbrain due to hindbrain neuronal hyperproliferation. Further study reveals that Rb1 controls the hindbrain morphogenesis by suppressing the expression of Gbx1/Gbx2, essential transcription factors for hindbrain development, through its binding to E2f3/Hdac1, respectively. Interestingly, we find that Gbx1 and Gbx2 are expressed in different types of hindbrain neurons, suggesting distinct roles in hindbrain morphogenesis. In summary, our study clarifies the specific role of RB1 in hindbrain neural cell proliferation and morphogenesis by regulating the E2f3-Gbx1 axis and the Hdac1-Gbx2 axis. These findings provide a research paradigm for exploring the differential proliferation of neurons in various brain regions.

Zinc finger protein 263 upregulates interleukin 33 and suppresses autophagy to accelerate the malignant progression of non-small cell lung cancer

PURPOSE

Non-small cell lung cancer (NSCLC) is a complex disease that remains a major public health concern worldwide. One promising avenue for NSCLC treatment is the targeting of transcription factors that regulate key pathways involved in cancer progression. In this study, we investigated the role of the transcription factor ZNF263 in NSCLC and its impact on the regulation of IL33, apoptosis, and autophagy.

METHODS

Levels of ZNF263 in tissues and cell lines were identified, after which the effects of its knockdown on cellular malignant behaviors, apoptosis and autophagy were assessed. Based on bioinformatics analysis, ZNF263 was found to bind to IL33 promoter, their mutual relationship was confirmed, as well as the role of IL33 in the regulation of ZNF263. The involvement of ZNF263 in the growth of xenograft tumors was assessed using tumor-bearing nude mouse models.

RESULTS

Experimental results revealed that ZNF263 was upregulated in NSCLC tissue samples and cell lines. Its expression level is positively correlated with cellular malignant behaviors. We further demonstrated that ZNF263 upregulated IL33 expression, which, in turn, promoted the proliferation and migration, inhibited apoptosis and autophagy in NSCLC cells. Furthermore, ZNF263 knockdown reduced the growth of xenograft tumors in nude mice.

CONCLUSION

This finding suggests that the inhibition of ZNF263 or IL33 may represent a novel therapeutic strategy for NSCLC. Importantly, our results highlight the crucial role of transcription factors in NSCLC and their potential as therapeutic targets.

DELVE: feature selection for preserving biological trajectories in single-cell data

Single-cell technologies can measure the expression of thousands of molecular features in individual cells undergoing dynamic biological processes. While examining cells along a computationally-ordered pseudotime trajectory can reveal how changes in gene or protein expression impact cell fate, identifying such dynamic features is challenging due to the inherent noise in single-cell data. Here, we present DELVE, an unsupervised feature selection method for identifying a representative subset of molecular features which robustly recapitulate cellular trajectories. In contrast to previous work, DELVE uses a bottom-up approach to mitigate the effects of confounding sources of variation, and instead models cell states from dynamic gene or protein modules based on core regulatory complexes. Using simulations, single-cell RNA sequencing, and iterative immunofluorescence imaging data in the context of cell cycle and cellular differentiation, we demonstrate how DELVE selects features that better define cell-types and cell-type transitions. DELVE is available as an open-source python package: https://github.com/jranek/delve .