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

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 .

A molecular atlas of adult C. elegans motor neurons reveals ancient diversity delineated by conserved transcription factor codes

Motor neurons (MNs) constitute an ancient cell type targeted by multiple adult-onset diseases. It is therefore important to define the molecular makeup of adult MNs in animal models and extract organizing principles. Here, we generate a comprehensive molecular atlas of adult Caenorhabditis elegans MNs and a searchable database. Single-cell RNA sequencing of 13,200 cells reveals that ventral nerve cord MNs cluster into 29 molecularly distinct subclasses. Extending C. elegans Neuronal Gene Expression Map and Network (CeNGEN) findings, all MN subclasses are delineated by distinct expression codes of either neuropeptide or transcription factor gene families. Strikingly, combinatorial codes of homeodomain transcription factor genes succinctly delineate adult MN diversity in both C. elegans and mice. Further, molecularly defined MN subclasses in C. elegans display distinct patterns of connectivity. Hence, our study couples the connectivity map of the C. elegans motor circuit with a molecular atlas of its constituent MNs and uncovers organizing principles and conserved molecular codes of adult MN diversity.

Construction of single-cell cross-species chromatin accessibility landscapes with combinatorial-hybridization-based ATAC-seq

Despite recent advances in single-cell genomics, the lack of maps for single-cell candidate cis-regulatory elements (cCREs) in non-mammal species has limited our exploration of conserved regulatory programs across vertebrates and invertebrates. Here, we developed a combinatorial-hybridization-based method for single-cell assay for transposase-accessible chromatin using sequencing (scATAC-seq) named CH-ATAC-seq, enabling the construction of single-cell accessible chromatin landscapes for zebrafish, Drosophila, and earthworms (Eisenia andrei). By integrating scATAC censuses of humans, monkeys, and mice, we systematically identified 152 distinct main cell types and around 0.8 million cell-type-specific cCREs. Our analysis provided insights into the conservation of neural, muscle, and immune lineages across species, while epithelial cells exhibited a higher organ-origin heterogeneity. Additionally, a large-scale gene regulatory network (GRN) was constructed in four vertebrates by integrating scRNA-seq censuses. Overall, our study provides a valuable resource for comparative epigenomics, identifying the evolutionary conservation and divergence of gene regulation across different species.

Intermittent fasting shifts the diurnal transcriptome atlas of transcription factors

Intermittent fasting remains a safe and effective strategy to ameliorate various age-related diseases, but its specific mechanisms are not fully understood. Considering that transcription factors (TFs) determine the response to environmental signals, here, we profiled the diurnal expression of 600 samples across four metabolic tissues sampled every 4 over 24 h from mice placed on five different feeding regimens to provide an atlas of TFs in biological space, time, and feeding regimen. Results showed that 1218 TFs exhibited tissue-specific and temporal expression profiles in ad libitum mice, of which 974 displayed significant oscillations at least in one tissue. Intermittent fasting triggered more than 90% (1161 in 1234) of TFs to oscillate somewhere in the body and repartitioned their tissue-specific expression. A single round of fasting generally promoted TF expression, especially in skeletal muscle and adipose tissues, while intermittent fasting mainly suppressed TF expression. Intermittent fasting down-regulated aging pathway and upregulated the pathway responsible for the inhibition of mammalian target of rapamycin (mTOR). Intermittent fasting shifts the diurnal transcriptome atlas of TFs, and mTOR inhibition may orchestrate intermittent fasting-induced health improvements. This atlas offers a reference and resource to understand how TFs and intermittent fasting may contribute to diurnal rhythm oscillation and bring about specific health benefits.

Single-cell division tracing and transcriptomics reveal cell types and differentiation paths in the regenerating lung

Understanding the molecular and cellular processes involved in lung epithelial regeneration may fuel the development of therapeutic approaches for lung diseases. We combine mouse models allowing diphtheria toxin-mediated damage of specific epithelial cell types and parallel GFP-labeling of functionally dividing cells with single-cell transcriptomics to characterize the regeneration of the distal lung. We uncover cell types, including Krt13+ basal and Krt15+ club cells, detect an intermediate cell state between basal and goblet cells, reveal goblet cells as actively dividing progenitor cells, and provide evidence that adventitial fibroblasts act as supporting cells in epithelial regeneration. We also show that diphtheria toxin-expressing cells can persist in the lung, express specific inflammatory factors, and transcriptionally resemble a previously undescribed population in the lungs of COVID-19 patients. Our study provides a comprehensive single-cell atlas of the distal lung that characterizes early transcriptional and cellular responses to concise epithelial injury, encompassing proliferation, differentiation, and cell-to-cell interactions.

Maternal KLF17 controls zygotic genome activation by acting as a messenger for RNA Pol II recruitment in mouse embryos

Initiation of timely and sufficient zygotic genome activation (ZGA) is crucial for the beginning of life, yet our knowledge of transcription factors (TFs) contributing to ZGA remains limited. Here, we screened the proteome of early mouse embryos after cycloheximide (CHX) treatment and identified maternally derived KLF17 as a potential TF for ZGA genes. Using a conditional knockout (cKO) mouse model, we further investigated the role of maternal KLF17 and found that it promotes embryonic development and full fertility. Mechanistically, KLF17 preferentially binds to promoters and recruits RNA polymerase II (RNA Pol II) in early 2-cell embryos, facilitating the expression of major ZGA genes. Maternal Klf17 knockout resulted in a downregulation of 9% of ZGA genes and aberrant RNA Pol II pre-configuration, which could be partially rescued by introducing exogenous KLF17. Overall, our study provides a strategy for screening essential ZGA factors and identifies KLF17 as a crucial TF in this process.

Myc Dysregulation in Activated Macrophages Initiates Iron-Mediated Lipid Peroxidation that Fuels Type I Interferon and Compromises TB Resistance

A quarter of human population is infected with Mycobacterium tuberculosis, but less than 10% of those infected develop clinical, mostly pulmonary, TB. To dissect mechanisms of susceptibility in immunocompetent individuals, we developed a genetically defined sst1-susceptible mouse model that uniquely reproduces a defining feature of human TB: development of necrotic lung lesions after infection with virulent Mtb. In this study, we explored the connectivity of the sst1-regulated pathways during prolonged macrophage activation with TNF. We determined that the aberrant response of the sst1-susceptible macrophages to TNF was primarily driven by conflicting Myc and antioxidant response pathways that resulted in a coordinated failure to properly sequester intracellular iron and activate ferroptosis inhibitor enzymes. Consequently, iron-mediated lipid peroxidation fueled IFNβ superinduction and sustained the Type I Interferon (IFN-I) pathway hyperactivity that locked the sst1-susceptible macrophages in a state of unresolving stress and compromised their resistance to Mtb. The accumulation of the aberrantly activated, stressed, macrophages within granuloma microenvironment led to the local failure of anti-tuberculosis immunity and tissue necrosis. Our findings suggest a novel link between metabolic dysregulation in macrophages and susceptibility to TB, offering insights into potential therapeutic targets aimed at modulating macrophage function and improving TB control.

Key homeobox transcription factors regulate the development of the firefly's adult light organ and bioluminescence

Adult fireflies exhibit unique flashing courtship signals, emitted by specialized light organs, which develop mostly independently from larval light organs during the pupal stage. The mechanisms of adult light organ development have not been thoroughly studied until now. Here we show that key homeobox transcription factors AlABD-B and AlUNC-4 regulate the development of adult light organs and bioluminescence in the firefly Aquatica leii. Interference with the expression of AlAbd-B and AlUnc-4 genes results in undeveloped or non-luminescent adult light organs. AlABD-B regulates AlUnc-4, and they interact with each other. AlABD-B and AlUNC-4 activate the expression of the luciferase gene AlLuc1 and some peroxins. Four peroxins are involved in the import of AlLUC1 into peroxisomes. Our study provides key insights into the development of adult light organs and flash signal control in fireflies.

Denervation alters the secretome of myofibers and thereby affects muscle stem cell lineage progression and functionality

Skeletal muscle function crucially depends on innervation while repair of skeletal muscle relies on resident muscle stem cells (MuSCs). However, it is poorly understood how innervation affects MuSC properties and thereby regeneration of skeletal muscle. Here, we report that loss of innervation causes precocious activation of MuSCs concomitant with the expression of markers of myogenic differentiation. This aberrant activation of MuSCs after loss of innervation is accompanied by profound alterations on the mRNA and protein level. Combination of muscle injury with loss of innervation results in impaired regeneration of skeletal muscle including shifts in myogenic populations concomitant with delayed maturation of regenerating myofibers. We further demonstrate that loss of innervation leads to alterations in myofibers and their secretome, which then affect MuSC behavior. In particular, we identify an increased secretion of Osteopontin and transforming growth factor beta 1 (Tgfb1) by myofibers isolated from mice which had undergone sciatic nerve transection. The altered secretome results in the upregulation of early activating transcription factors, such as Junb, and their target genes in MuSCs. However, the combination of different secreted factors from myofibers after loss of innervation is required to cause the alterations observed in MuSCs after loss of innervation. These data demonstrate that loss of innervation first affects myofibers causing alterations in their secretome which then affect MuSCs underscoring the importance of proper innervation for MuSC functionality and regeneration of skeletal muscle.

Associating growth factor secretions and transcriptomes of single cells in nanovials using SEC-seq

Cells secrete numerous bioactive molecules that are essential for the function of healthy organisms. However, scalable methods are needed to link individual cell secretions to their transcriptional state over time. Here, by developing and using secretion-encoded single-cell sequencing (SEC-seq), which exploits hydrogel particles with subnanolitre cavities (nanovials) to capture individual cells and their secretions, we simultaneously measured the secretion of vascular endothelial growth factor A (VEGF-A) and the transcriptome for thousands of individual mesenchymal stromal cells. Our data indicate that VEGF-A secretion is heterogeneous across the cell population and is poorly correlated with the VEGFA transcript level. The highest VEGF-A secretion occurs in a subpopulation of mesenchymal stromal cells characterized by a unique gene expression signature comprising a surface marker, interleukin-13 receptor subunit alpha 2 (IL13RA2), which allowed the enrichment of this subpopulation. SEC-seq enables the identification of gene signatures linked to specific secretory states, facilitating mechanistic studies, the isolation of secretory subpopulations and the development of means to modulate cellular secretion.

The pleiotropic enhancer enh9 promotes cell proliferation and migration in non-small cell lung cancer via ERMP1 and PD-L1

Enhancers, cis-acting DNA elements for transcriptional regulation, are important regulators of cell identity and disease. However, of the hundreds of thousands of enhancers annotated in the human genome, only a few have been studied for their regulatory mechanisms and functions in cancer progression and therapeutic resistance. Here, we report the pleiotropy of one enhancer (named enh9) in both cell proliferation and migration in non-small cell lung cancer (NSCLC) cells. By integrating multi-genomic data, ERMP1 and PD-L1 were screened out as potential targets of enh9. CUT&Tag sequencing demonstrated that enh9 was involved in the genomic interactions between the transcription factor RELA and the promoters of ERMP1 and PD-L1. In addition, ERMP1 and PD-L1 were validated to be involved in cell proliferation and migration, respectively. Our study fully elucidated the function and transcriptional regulation mechanisms of enh9 in NSCLC. The exploration on enhancers is promising to provide new insights for cancer diagnosis and therapy.

Histone demethylase KDM1A promotes hepatic steatosis and inflammation by increasing chromatin accessibility in NAFLD

Nonalcoholic fatty liver disease (NAFLD) is the most prevalent chronic liver disease without specific Food and Drug Administration-approved drugs. Recent advances suggest that chromatin remodeling and epigenetic alteration contribute to the development of NAFLD. The functions of the corresponding molecular modulator in NAFLD, however, are still elusive. KDM1A, commonly known as lysine-specific histone demethylase 1, has been reported to increase glucose uptake in hepatocellular carcinoma. In addition, a recent study suggests that inhibition of KDM1A reduces lipid accumulation in primary brown adipocytes. We here investigated the role of KDM1A, one of the most important histone demethylases, in NAFLD. In this study, we observed a significant upregulation of KDM1A in NAFLD mice, monkeys, and humans compared to the control group. Based on these results, we further found that the KDM1A can exacerbate lipid accumulation and inflammation in hepatocytes and mice. Mechanistically, KDM1A exerted its effects by elevating chromatin accessibility, subsequently promoting the development of NAFLD. Furthermore, the mutation of KDM1A blunted its capability to promote the development of NAFLD. In summary, our study discovered that KDM1A exacerbates hepatic steatosis and inflammation in NAFLD via increasing chromatin accessibility, further indicating the importance of harnessing chromatin remodeling and epigenetic alteration in combating NAFLD. KDM1A might be considered as a potential therapeutic target in this regard.

An intronic enhancer of Cebpa regulates adipocyte differentiation and adipose tissue development via long-range loop formation

Cebpa is a master transcription factor gene for adipogenesis. However, the mechanisms of enhancer-promoter chromatin interactions controlling Cebpa transcriptional regulation during adipogenic differentiation remain largely unknown. To reveal how the three-dimensional structure of Cebpa changes during adipogenesis, we generated high-resolution chromatin interactions of Cebpa in 3T3-L1 preadipocytes and 3T3-L1 adipocytes using circularized chromosome conformation capture sequencing (4C-seq). We revealed dramatic changes in chromatin interactions and chromatin status at interaction sites during adipogenic differentiation. Based on this, we identified five active enhancers of Cebpa in 3T3-L1 adipocytes through epigenomic data and luciferase reporter assays. Next, epigenetic repression of Cebpa-L1-AD-En2 or -En3 by the dCas9-KRAB system significantly down-regulated Cebpa expression and inhibited adipocyte differentiation. Furthermore, experimental depletion of cohesin decreased the interaction intensity between Cebpa-L1-AD-En2 and the Cebpa promoter and down-regulated Cebpa expression, indicating that long-range chromatin loop formation was mediated by cohesin. Two transcription factors, RXRA and PPARG, synergistically regulate the activity of Cebpa-L1-AD-En2. To test whether Cebpa-L1-AD-En2 plays a role in adipose tissue development, we injected dCas9-KRAB-En2 lentivirus into the inguinal white adipose tissue (iWAT) of mice to suppress the activity of Cebpa-L1-AD-En2. Repression of Cebpa-L1-AD-En2 significantly decreased Cebpa expression and adipocyte size, altered iWAT transcriptome, and affected iWAT development. We identified functional enhancers regulating Cebpa expression and clarified the crucial roles of Cebpa-L1-AD-En2 and Cebpa promoter interaction in adipocyte differentiation and adipose tissue development.

Single-cell profiling of MC1R-inhibited melanocytes

The human red hair color (RHC) trait is caused by increased pheomelanin (red-yellow) and reduced eumelanin (black-brown) pigment in skin and hair due to diminished melanocortin 1 receptor (MC1R) function. In addition, individuals harboring the RHC trait are predisposed to melanoma development. While MC1R variants have been established as causative of RHC and are a well-defined risk factor for melanoma, it remains unclear mechanistically why decreased MC1R signaling alters pigmentation and increases melanoma susceptibility. Here, we use single-cell RNA sequencing (scRNA-seq) of melanocytes isolated from RHC mouse models to define a MC1R-inhibited Gene Signature (MiGS) comprising a large set of previously unidentified genes which may be implicated in melanogenesis and oncogenic transformation. We show that one of the candidate MiGS genes, TBX3, a well-known anti-senescence transcription factor implicated in melanoma progression, binds both E-box and T-box elements to regulate genes associated with melanogenesis and senescence bypass. Our results provide key insights into further mechanisms by which melanocytes with reduced MC1R signaling may regulate pigmentation and offer new candidates of study toward understanding how individuals with the RHC phenotype are predisposed to melanoma.

Cell specification and functional interactions in the pig blastocyst inferred from single-cell transcriptomics and uterine fluids proteomics

The embryonic development of the pig comprises a long in utero pre- and peri-implantation development, which dramatically differs from mice and humans. During this peri-implantation period, a complex series of paracrine signals establishes an intimate dialogue between the embryo and the uterus. To better understand the biology of the pig blastocyst during this period, we generated a large dataset of single-cell RNAseq from early and hatched blastocysts, spheroid and ovoid conceptus and proteomic datasets from corresponding uterine fluids. Our results confirm the molecular specificity and functionality of the three main cell populations. We also discovered two previously unknown subpopulations of the trophectoderm, one characterised by the expression of LRP2, which could represent progenitor cells, and the other, expressing pro-apoptotic markers, which could correspond to the Rauber's layer. Our work provides new insights into the biology of these populations, their reciprocal functional interactions, and the molecular dialogue with the maternal uterine environment.

Optimizing sequence design strategies for perturbation MPRAs: a computational evaluation framework

The advent of perturbation-based massively parallel reporter assays (MPRAs) technique has facilitated the delineation of the roles of non-coding regulatory elements in orchestrating gene expression. However, computational efforts remain scant to evaluate and establish guidelines for sequence design strategies for perturbation MPRAs. In this study, we propose a framework for evaluating and comparing various perturbation strategies for MPRA experiments. Within this framework, we benchmark three different perturbation approaches from the perspectives of alteration in motif-based profiles, consistency of MPRA outputs, and robustness of models that predict the activities of putative regulatory motifs. While our analyses show very similar results across multiple benchmarking metrics, the predictive modeling for the approach involving random nucleotide shuffling shows significant robustness compared with the other two approaches. Thus, we recommend designing sequences by randomly shuffling the nucleotides of the perturbed site in perturbation-MPRA, followed by a coherence check to prevent the introduction of other variations of the target motifs. In summary, our evaluation framework and the benchmarking findings create a resource of computational pipelines and highlight the potential of perturbation-MPRA in predicting non-coding regulatory activities.

Comparative transcriptomics coupled to developmental grading via transgenic zebrafish reporter strains identifies conserved features in neutrophil maturation

Neutrophils are evolutionarily conserved innate immune cells playing pivotal roles in host defense. Zebrafish models have contributed substantially to our understanding of neutrophil functions but similarities to human neutrophil maturation have not been systematically characterized, which limits their applicability to studying human disease. Here we show, by generating and analysing transgenic zebrafish strains representing distinct neutrophil differentiation stages, a high-resolution transcriptional profile of neutrophil maturation. We link gene expression at each stage to characteristic transcription factors, including C/ebp-β, which is important for late neutrophil maturation. Cross-species comparison of zebrafish, mouse, and human samples confirms high molecular similarity of immature stages and discriminates zebrafish-specific from pan-species gene signatures. Applying the pan-species neutrophil maturation signature to RNA-sequencing data from human neuroblastoma patients reveals association between metastatic tumor cell infiltration in the bone marrow and an overall increase in mature neutrophils. Our detailed neutrophil maturation atlas thus provides a valuable resource for studying neutrophil function at different stages across species in health and disease.

Expression network analysis of bovine skin infested with Rhipicephalus australis identifies pro-inflammatory genes contributing to tick susceptibility

The skin is the primary feeding site of ticks that infest livestock animals such as cattle. The highly specialised functions of skin at the molecular level may be a factor contributing to variation in susceptibility to tick infestation; but these remain to be well defined. The aim of this study was to investigate the bovine skin transcriptomic profiles of tick-naïve and tick-infested cattle and to uncover the gene expression networks that influence contrasting phenotypes of host resistance to ticks. RNA-Seq data was obtained from skin of Brangus cattle with high (n = 5) and low (n = 6) host resistance at 0 and 12 weeks following artificial tick challenge with Rhipicephalus australis larvae. No differentially expressed genes were detected pre-infestation between high and low resistance groups, but at 12-weeks there were 229 differentially expressed genes (DEGs; FDR < 0.05), of which 212 were the target of at least 1866 transcription factors (TFs) expressed in skin. Regulatory impact factor (RIF) analysis identified 158 significant TFs (P < 0.05) of which GRHL3, and DTX1 were also DEGs in the experiment. Gene term enrichment showed the significant TFs and DEGs were enriched in processes related to immune response and biological pathways related to host response to infectious diseases. Interferon Type 1-stimulated genes, including MX2, ISG15, MX1, OAS2 were upregulated in low host resistance steers after repeated tick challenge, suggesting dysregulated wound healing and chronic inflammatory skin processes contributing to host susceptibility to ticks. The present study provides an assessment of the bovine skin transcriptome before and after repeated tick challenge and shows that the up-regulation of pro-inflammatory genes is a prominent feature in the skin of tick-susceptible animals. In addition, the identification of transcription factors with high regulatory impact provides insights into the potentially meaningful gene-gene interactions involved in the variation of phenotypes of bovine host resistance to ticks.

Identification of gene biomarkers for brain diseases via multi-network topological semantics extraction and graph convolutional network

BACKGROUND

Brain diseases pose a significant threat to human health, and various network-based methods have been proposed for identifying gene biomarkers associated with these diseases. However, the brain is a complex system, and extracting topological semantics from different brain networks is necessary yet challenging to identify pathogenic genes for brain diseases.

RESULTS

In this study, we present a multi-network representation learning framework called M-GBBD for the identification of gene biomarker in brain diseases. Specifically, we collected multi-omics data to construct eleven networks from different perspectives. M-GBBD extracts the spatial distributions of features from these networks and iteratively optimizes them using Kullback-Leibler divergence to fuse the networks into a common semantic space that represents the gene network for the brain. Subsequently, a graph consisting of both gene and large-scale disease proximity networks learns representations through graph convolution techniques and predicts whether a gene is associated which brain diseases while providing associated scores. Experimental results demonstrate that M-GBBD outperforms several baseline methods. Furthermore, our analysis supported by bioinformatics revealed CAMP as a significantly associated gene with Alzheimer's disease identified by M-GBBD.

CONCLUSION

Collectively, M-GBBD provides valuable insights into identifying gene biomarkers for brain diseases and serves as a promising framework for brain networks representation learning.

BCLAF1 binds SPOP to stabilize PD-L1 and promotes the development and immune escape of hepatocellular carcinoma

Interaction between programmed death-1 (PD-1) ligand 1 (PD-L1) on tumor cells and PD-1 on T cells allows tumor cells to evade T cell-mediated immune surveillance. Strategies targeting PD-1/PD-L1 have shown clinical benefits in a variety of cancers. However, limited response rates in hepatocellular carcinoma (HCC) have prompted us to investigate the molecular regulation of PD-L1. Here, we identify B cell lymphoma-2-associated transcription factor 1 (BCLAF1) as a key PD-L1 regulator in HCC. Specifically, BCLAF1 interacts with SPOP, an E3 ligase that mediates the ubiquitination and degradation of PD-L1, thereby competitively inhibiting SPOP-PD-L1 interaction and subsequent ubiquitination and degradation of PD-L1. Furthermore, we determined an SPOP-binding consensus (SBC) motif mediating the BCLAF1-SPOP interaction on BCLAF1 protein and mutation of BCLAF1-SBC motif disrupts the regulation of the SPOP-PD-L1 axis. In addition, BCLAF1 expression was positively correlated with PD-L1 expression and negatively correlated with biomarkers of T cell activation, including CD3 and CD8, as well as with the level of immune cell infiltration in HCC tissues. Besides, BCLAF1 depletion leads to a significant reduction of PD-L1 expression in vitro, and this reduction of PD-L1 promoted T cell-mediated cytotoxicity. Notably, overexpression of BCLAF1 sensitized tumor cells to checkpoint therapy in an in vitro HCC cells-Jurkat cells co-culture model, whereas BCLAF1-SBC mutant decreased tumor cell sensitivity to checkpoint therapy, suggesting that BCLAF1 and its SBC motif serve as a novel therapeutic target for enhancing anti-tumor immunity in HCC.

Repaglinide restrains HCC development and progression by targeting FOXO3/lumican/p53 axis

PURPOSE

The recent focus on the roles of N-linked glycoproteins in carcinogenesis across various malignancies has prompted our exploration of aberrantly expressed glycoproteins responsible for HCC progression and potential therapeutic strategy.

METHODS

Mass spectrometry was applied to initially identify abnormally expressed glycoproteins in HCC, which was further assessed by immunohistochemistry (IHC) staining. The role of selected glycoprotein on HCC development and underlying mechanism was systematically investigated by colony formation, mouse xenograft, RNA-sequencing and western blot assays, etc. Chromatin immunoprecipitation (ChIP) and luciferase assays were performed to explore potential transcription factors (TFs) of selected glycoprotein. The regulation of repaglinide (RPG) on expression of lumican and downstream effectors was assessed by western blot and IHC, while its impact on malignant phenotypes of HCC was explored through in vitro and in vivo analyses, including a murine NASH-HCC model established using western diet and carbon tetrachloride (CCl4).

RESULTS

Lumican exhibited upregulation in both serum and tumor tissue, with elevated expression associated with an inferior prognosis in HCC patients. Knockdown of lumican resulted in significantly reduced growth of HCC in vitro and in vivo. Mechanically, lumican promoted HCC malignant phenotypes by inhibiting the p53/p21 signaling pathway. Forkhead Box O3 (FOXO3) was identified as the TF of lumican that transcriptionally enhanced its expression. Without silencing FOXO3, RPG blocked the binding of FOXO3 to the promoter region of lumican, thereby inhibiting the activation of lumican/p53/p21 axis. Mice treated with RPG developed fewer and smaller HCCs than those in the control group at 24 weeks after establishment.

CONCLUSION

Our results indicate that RPG prevented the development and progression of HCC via alteration of FOXO3/lumican/p53 axis.

Transcription factor ATMIN facilitates chemoresistance in nasopharyngeal carcinoma

Despite that the docectaxel-cisplatin-5-fluorouracil (TPF) induction chemotherapy has greatly improved patients' survival and became the first-line treatment for advanced nasopharyngeal carcinoma (NPC), not all patients could benefit from this therapy. The mechanism underlying the TPF chemoresistance remains unclear. Here, by analyzing gene-expression microarray data and survival of patients who received TPF chemotherapy, we identify transcription factor ATMIN as a chemoresistance gene in response to TPF chemotherapy in NPC. Mass spectrometry and Co-IP assays reveal that USP10 deubiquitinates and stabilizes ATMIN protein, resulting the high-ATMIN expression in NPC. Knockdown of ATMIN suppresses the cell proliferation and facilitates the docetaxel-sensitivity of NPC cells both in vitro and in vivo, while overexpression of ATMIN exerts the opposite effect. Mechanistically, ChIP-seq combined with RNA-seq analysis suggests that ATMIN is associated with the cell death signaling and identifies ten candidate target genes of ATMIN. We further confirm that ATMIN transcriptionally activates the downstream target gene LCK and stabilizes it to facilitate cell proliferation and docetaxel resistance. Taken together, our findings broaden the insight into the molecular mechanism of chemoresistance in NPC, and the USP10-ATMIN-LCK axis provides potential therapeutic targets for the management of NPC.

METTL14-mediated m6A epitranscriptomic modification contributes to chemotherapy-induced neuropathic pain by stabilizing GluN2A expression via IGF2BP2

Epigenetics is a biological process that modifies and regulates gene expression, affects neuronal function, and contributes to pain. However, the mechanism by which epigenetics facilitates and maintains chronic pain is poorly understood. We aimed to determine whether N6-methyladenosine (m6A) specifically modified by methyltransferase-like 14 (METTL14) alters neuronal activity and governs pain by sensitizing the GluN2A subunit of the N-methyl-d-aspartate receptor (NMDAR) in the dorsal root ganglion (DRG) neurons in a model of chemotherapy-induced neuropathic pain (CINP). Using dot blotting, immunofluorescence, gain/loss-of-function, and behavioral assays, we found that m6A levels were upregulated in L4-L6 DRG neurons in CINP in a DBP/METTL14-dependent manner, which was also confirmed in human DRGs. Blocking METTL14 reduced m6A methylation and attenuated pain hypersensitivity. Mechanistically, METTL14-mediated m6A modification facilitated the synaptic plasticity of DRG neurons by enhancing the GluN2A subunit of NMDAR, and inhibiting METTL14 blocked this effect. In contrast, overexpression of METTL14 upregulated m6A modifications, enhanced presynaptic NMDAR activity in DRG neurons, and facilitated pain sensation. Our findings reveal a previously unrecognized mechanism of METTL14-mediated m6A modification in DRG neurons to maintain neuropathic pain. Targeting these molecules may provide a new strategy for pain treatment.

A single-cell time-lapse of mouse prenatal development from gastrula to birth

The house mouse (Mus musculus) is an exceptional model system, combining genetic tractability with close evolutionary affinity to humans1,2. Mouse gestation lasts only 3 weeks, during which the genome orchestrates the astonishing transformation of a single-cell zygote into a free-living pup composed of more than 500 million cells. Here, to establish a global framework for exploring mammalian development, we applied optimized single-cell combinatorial indexing3 to profile the transcriptional states of 12.4 million nuclei from 83 embryos, precisely staged at 2- to 6-hour intervals spanning late gastrulation (embryonic day 8) to birth (postnatal day 0). From these data, we annotate hundreds of cell types and explore the ontogenesis of the posterior embryo during somitogenesis and of kidney, mesenchyme, retina and early neurons. We leverage the temporal resolution and sampling depth of these whole-embryo snapshots, together with published data4-8 from earlier timepoints, to construct a rooted tree of cell-type relationships that spans the entirety of prenatal development, from zygote to birth. Throughout this tree, we systematically nominate genes encoding transcription factors and other proteins as candidate drivers of the in vivo differentiation of hundreds of cell types. Remarkably, the most marked temporal shifts in cell states are observed within one hour of birth and presumably underlie the massive physiological adaptations that must accompany the successful transition of a mammalian fetus to life outside the womb.

Pig-eRNAdb: a comprehensive enhancer and eRNA dataset of pigs

Enhancers and the enhancer RNAs (eRNAs) have been strongly implicated in regulations of transcriptions. Based the multi-omics data (ATAC-seq, ChIP-seq and RNA-seq) from public databases, Pig-eRNAdb is a dataset that comprehensively integrates enhancers and eRNAs for pigs using the machine learning strategy, which incorporates 82,399 enhancers and 37,803 eRNAs from 607 samples across 15 tissues of pigs. This user-friendly dataset covers a comprehensive depth of enhancers and eRNAs annotation for pigs. The coordinates of enhancers and the expression patterns of eRNAs are downloadable. Besides, thousands of regulators on eRNAs, the target genes of eRNAs, the tissue-specific eRNAs, and the housekeeping eRNAs are also accessible as well as the sequence similarity of eRNAs with humans. Moreover, the tissue-specific eRNA-trait associations encompass 652 traits are also provided. It will crucially facilitate investigations on enhancers and eRNAs with Pig-eRNAdb as a reference dataset in pigs.

Simulated microgravity altered the gene expression profiles and inhibited the proliferation of Kupffer cells in the early phase by downregulating LMO2 and EZH2

Microgravity is a primary challenge that need to overcome, when human travel to space. Our study provided evidence that Kupffer cells (KCs) are sensitive to simulated microgravity (SMG), and no similar research report has been found in the literature. Using transcriptome sequencing technology, it was showed that 631 genes were upregulated and 801 genes were downregulated in KCs after treatment under SMG for 3 days. The GO analysis indicated that the proliferation of KCs was affected when exposed to SMG for 3 days. CCK-8 assay confirmed that the proliferation of KCs was inhibited in the third day under the environment of SMG. Furthermore, we identified 8 key genes that affect the proliferation of KCs and predicted 2 transcription factors (TFs) that regulate the 8 key genes. Significantly, we found that microgravity could affect the expression of LMO2 and EZH2 to reduce the transcription of Racgap1, Ccna2, Nek2, Aurka, Plk1, Haus4, Cdc20, Bub1b, which resulting in the reduction in KCs proliferation. These finding suggested that the inhibition of KCs proliferation under microgravity may influence the homeostasis of liver, and LMO2 and EZH2 can be the targets in management of KCs' disturbance in the future practice of space medicine.

Deciphering the role of alternative splicing as a potential regulator in fat-tail development of sheep: a comprehensive RNA-seq based study

Although research on alternative splicing (AS) has been widely conducted in mammals, no study has investigated the splicing profiles of genes involved in fat-tail formation in sheep. Here, for the first time, a comprehensive study was designed to investigate the profile of AS events and their involvement in fat-tail development of sheep. In total, 45 RNA-Seq samples related to seven different studies, which have compared the fat-tailed vs thin-tailed sheep breeds, were analyzed. Two independent tools, rMATS and Whippet, along with a set of stringent filters were applied to identify differential AS (DAS) events between the breeds per each study. Only DAS events that were detected by both tools as well as in at least three datasets with the same ΔPSI trend (percent spliced in), were considered as the final high-confidence set of DAS genes. Final results revealed 130 DAS skipped exon events (69 negative and 61 positive ΔPSI) belonged to 124 genes. Functional enrichment analysis highlighted the importance of the genes in the underlying molecular mechanisms of fat metabolism. Moreover, protein-protein interaction network analysis revealed that DAS genes are significantly connected. Of DAS genes, five transcription factors were found that were enriched in the biological process associated with lipid metabolism like "Fat Cell Differentiation". Further investigations of the findings along with a comprehensive literature review provided a reliable list of candidate genes that may potentially contribute to fat-tail formation including HSD11B1, SIRT2, STRN3 and TCF7L2. Based on the results, it can be stated that the AS patterns may have evolved, during the evolution of sheep breeds, as another layer of regulation to contribute to biological complexity by reprogramming the gene regulatory networks. This study provided the theoretical basis of the molecular mechanisms behind the sheep fat-tail development in terms of AS.

Mapping Expression Quantitative Trait Loci Targeting Candidate Genes for Pregnancy in Beef Cows

Despite collective efforts to understand the complex regulation of reproductive traits, no causative genes and/or mutations have been reported yet. By integrating genomics and transcriptomics data, potential regulatory mechanisms may be unveiled, providing opportunities to dissect the genetic factors governing fertility. Herein, we identified regulatory variants from RNA-Seq data associated with gene expression regulation in the uterine luminal epithelial cells of beef cows. We identified 4676 cis and 7682 trans eQTLs (expression quantitative trait loci) affecting the expression of 1120 and 2503 genes, respectively (FDR < 0.05). These variants affected the expression of transcription factor coding genes (71 cis and 193 trans eQTLs) and genes previously reported as differentially expressed between pregnant and nonpregnant cows. Functional over-representation analysis highlighted pathways related to metabolism, immune response, and hormone signaling (estrogen and GnRH) affected by eQTL-regulated genes (p-value ≤ 0.01). Furthermore, eQTLs were enriched in QTL regions for 13 reproduction-related traits from the CattleQTLdb (FDR ≤ 0.05). Our study provides novel insights into the genetic basis of reproductive processes in cattle. The underlying causal mechanisms modulating the expression of uterine genes warrant further investigation.

HHCDB: a database of human heterochromatin regions

Heterochromatin plays essential roles in eukaryotic genomes, such as regulating genes, maintaining genome integrity and silencing repetitive DNA elements. Identifying genome-wide heterochromatin regions is crucial for studying transcriptional regulation. We propose the Human Heterochromatin Chromatin Database (HHCDB) for archiving heterochromatin regions defined by specific or combined histone modifications (H3K27me3, H3K9me2, H3K9me3) according to a unified pipeline. 42 839 743 heterochromatin regions were identified from 578 samples derived from 241 cell-types/cell lines and 92 tissue types. Genomic information is provided in HHCDB, including chromatin location, gene structure, transcripts, distance from transcription start site, neighboring genes, CpG islands, transposable elements, 3D genomic structure and functional annotations. Furthermore, transcriptome data from 73 single cells were analyzed and integrated to explore cell type-specific heterochromatin-related genes. HHCDB affords rich visualization through the UCSC Genome Browser and our self-developed tools. We have also developed a specialized online analysis platform to mine differential heterochromatin regions in cancers. We performed several analyses to explore the function of cancer-specific heterochromatin-related genes, including clinical feature analysis, immune cell infiltration analysis and the construction of drug-target networks. HHCDB is a valuable resource for studying epigenetic regulation, 3D genomics and heterochromatin regulation in development and disease. HHCDB is freely accessible at http://hhcdb.edbc.org/.

Deciphering the Role of Rapidly Evolving Conserved Elements in Primate Brain Development and Exploring Their Potential Involvement in Alzheimer's Disease

Although previous studies have identified human-specific accelerated regions as playing a key role in the recent evolution of the human brain, the characteristics and cellular functions of rapidly evolving conserved elements (RECEs) in ancestral primate lineages remain largely unexplored. Here, based on large-scale primate genome assemblies, we identify 888 RECEs that have been highly conserved in primates that exhibit significantly accelerated substitution rates in the ancestor of the Simiiformes. This primate lineage exhibits remarkable morphological innovations, including an expanded brain mass. Integrative multiomic analyses reveal that RECEs harbor sequences with potential cis-regulatory functions that are activated in the adult human brain. Importantly, genes linked to RECEs exhibit pronounced expression trajectories in the adult brain relative to the fetal stage. Furthermore, we observed an increase in the chromatin accessibility of RECEs in oligodendrocytes from individuals with Alzheimer's disease (AD) compared to that of a control group, indicating that these RECEs may contribute to brain aging and AD. Our findings serve to expand our knowledge of the genetic underpinnings of brain function during primate evolution.

Construction of hub transcription factor-microRNAs-messenger RNA regulatory network in recurrent implantation failure

PURPOSE

Recurrent implantation failure (RIF) affects up to 10% of in vitro fertilization (IVF) patients worldwide. However, the pathogenesis of RIF remains unclear. This study was aimed at identifying hub transcription factors (TFs) of RIF in bioinformatics approaches.

METHODS

The GSE111974 (mRNA), GSE71332 (miRNA), and GSE103465 (mRNA) datasets were downloaded from the Gene Expression Omnibus database from human endometrial tissue using R version 4.2.1 and used to identify differentially expressed TFs (DETFs), differentially expressed miRNAs, and differentially expressed genes for RIF, respectively. DETFs were subjected to functional enrichment analysis and the protein-protein interaction network analysis using the Search Tool for the Retrieval of Interacting Genes (version 11.5) database. Hub TFs were identified using the cytoHubb plug-in, after which a hub TF-miRNA-mRNA network was constructed using Cytoscape v3.8.2.

RESULTS

Fifty-seven DETFs were identified, in which Gene Ontology analysis revealed to be mainly involved in the regulation of transcription. Kyoto Encyclopedia of Genes and Genomes pathway analysis suggested that DETFs were enriched in transcriptional misregulation in cancer, aldosterone synthesis and secretion, AMPK signaling pathway, and cGMP-PKG signaling pathway. EOMES, NKX2-1, and POU5F1 were identified as hub TFs, and a hub TF-miRNA-mRNA regulatory network was constructed using these three hub TFs, four miRNAs, and four genes.

CONCLUSION

Collectively, we identified three promising molecular biomarkers for the diagnosis of RIF, which may further be potential therapeutic targets. This study provides novel insights into the molecular mechanisms underlying RIF. However, further experiments are required to verify these results.

Cross-species single-cell landscape of vertebrate pineal gland

The pineal gland has evolved from a photoreceptive organ in fish to a neuroendocrine organ in mammals. This study integrated multiple daytime single-cell RNA-seq datasets from the pineal glands of zebrafish, rats, and monkeys, providing a detailed examination of the evolutionary transition at single-cell resolution. We identified key factors responsible for the anatomical and functional transformation of the pineal gland. We retrieved and integrated daytime single-cell transcriptomic datasets from the pineal glands of zebrafish, rats, and monkeys, resulting in a total of 22 431 cells after rigorous quality filtering. Comparative analysis was then conducted to elucidate the evolution of pineal cells, their photosensitivity, their role in melatonin production, and the signaling processes within the glands of these species. Our analysis identified distinct cellular compositions of the pineal gland in zebrafish, rats, and monkeys. Zebrafish photoreceptors exhibited comprehensive phototransduction gene expression, while specific genes, including transducin (Gngt1, Gnb3, and Gngt2) and phosducin (Pdc), were consistently present in mammalian pinealocytes. We found transcriptional similarities between the pineal gland and retina, underscoring shared evolutionary and functional pathways. Zebrafish displayed unique light-responsive circadian gene activity compared to rats and monkeys. Key ligand-receptor interactions were identified, especially involving MDK and PTN, influencing melatonin synthesis across species. Furthermore, we observed species-specific GPCR (G protein-coupled receptors) expressions related to melatonin synthesis and their alignment with retinal expressions. Our findings also highlighted specific transcription factors (TFs) and regulatory networks associated with pineal gland evolution and function. Our study provides a detailed analysis of the pineal gland's evolution from fish to mammals. We identified key transcriptional changes and controls that highlight the gland's functional diversity. Notably, we found significant ligand-receptor interactions influencing melatonin synthesis and demonstrated parallels between pineal and retinal expressions. These insights enhance our understanding of the pineal gland's role in phototransduction, melatonin production, and circadian rhythms in vertebrates.

Tumor cell SPTBN1 inhibits M2 polarization of macrophages by suppressing CXCL1 expression

Tumor-associated macrophages (TAMs) are the most abundant immune cells in the tumor microenvironment, and the M2-type TAMs can promote tumor growth, invasion and angiogenesis, and suppress antitumor immune responses. It has been reported that spectrin beta, non-erythrocytic 1 (SPTBN1) may inhibit the infiltration of macrophages in Sptbn1+/-  mouse liver, but whether tumor SPTBN1 affects TAMs polarization remains unclear. This study investigated the effect and mechanism of tumor cell SPTBN1 on polarization and migration of TAMs in hepatoma and breast cancer. By analyzing tumor immune databases, we found a negative correlation between SPTBN1 and abundance of macrophages and myeloid-derived suppressor cells (MDSCs) in the tumor microenvironment. By reverse transcription-quantitative real-time PCR assays and cell migration assays, the migration and M2 polarization of macrophages were enhanced by the culture medium from hepatocellular carcinoma cell line PLC/PRF/5, SNU449, and breast cancer cell line MDA-MB-231 with SPTBN1 suppression, which could be reversed by CXCL1 neutralizing antibody MAB275. Meanwhile, the ability of migration and colony formation of PLC/PRF/5, SNU449, and MDA-MB-231 cells were promoted when coculture with M2 macrophages. We also found that SPTBN1 regulated CXCL1 through p65 by cytoplasmic-nuclear protein isolation experiments and ChIP-qPCR. Our data suggest that tumor cell SPTBN1 inhibits migration and M2-type polarization of TAMs by reducing the expression and secretion of CXCL1 via inhibiting p65 nuclear localization.

LncHLEF promotes hepatic lipid synthesis through miR-2188-3p/GATA6 axis and encoding peptides and enhances intramuscular fat deposition via exosome

Long noncoding RNAs (lncRNAs) have emergingly been implicated in mammalian lipid metabolism. However, their biological functions and regulatory mechanisms underlying adipogenesis remain largely elusive in chicken. Here, we systematically characterized the genome-wide full-length lncRNAs in the livers of pre- and peak-laying hens, and identified a novel intergenic lncRNA, lncHLEF, an RNA macromolecule with a calculated molecular weight of 433 kDa. lncHLEF was primarily distributed in cytoplasm of chicken hepatocyte and significantly up-regulated in livers of peak-laying hens. Functionally, lncHLEF could promote hepatocyte lipid droplet formation, triglycerides and total cholesterol contents. Mechanistically, lncHLEF could not only serve as a competitive endogenous RNA to modulate miR-2188-3p/GATA6 axis, but also encode three small functional polypeptides that directly interact with ACLY protein to enable its stabilization. Importantly, adeno-associated virus-mediated liver-specific lncHLEF overexpression resulted in increased hepatic lipid synthesis and intramuscular fat (IMF) deposition, but did not alter abdominal fat (AbF) deposition. Furthermore, hepatocyte lncHLEF could be delivered into intramuscular and abdominal preadipocytes via hepatocyte-secreted exosome to enhance intramuscular preadipocytes differentiation without altering abdominal preadipocytes differentiation. In conclusion, this study revealed that the lncHLEF could promote hepatic lipid synthesis through two independent regulatory mechanisms, and could enhance IMF deposition via hepatocyte-adipocyte communications mediated by exosome.

Integrated analysis of genes and long non-coding RNAs in trachea transcriptome to decipher the host response during Newcastle disease challenge in different breeds of chicken

Newcastle disease is a highly infectious economically devastating disease caused by Newcastle disease Virus in Chicken (Gallus gallus). Leghorn and Fayoumi are two breeds which show differential resistance patterns towards NDV. This study aims to identify the differentially expressed genes and lncRNAs during NDV challenge which could play a potential role in this differential resistance pattern. A total of 552 genes and 1580 lncRNAs were found to be differentially expressing. Of them, 52 genes were annotated with both Immune related pathways and Gene ontologies. We found that most of these genes were upregulated in Leghorn between normal and challenged chicken but several were down regulated between different timepoints after NDV challenge, while Fayoumi showed no such downregulation. We also observed that higher number of positively correlating lncRNAs was found to be downregulated along with these genes. This shows that although Leghorn is showing higher number of differentially expressed genes in challenged than in non-challenged, most of them were downregulated during the disease between different timepoints. With this we hypothesize that the downregulation of immune related genes and co-expressing lncRNAs could play a significant role behind the Leghorn being comparatively susceptible breed than Fayoumi. The computational pipeline is available at https://github.com/Venky2804/FHSpipe.

Dynamic chromatin remodeling in cycling human endometrium at single-cell level

Estrogen-dependent proliferation followed by progesterone-dependent differentiation of the endometrium culminates in a short implantation window. We performed single-cell assay for transposase-accessible chromatin with sequencing on endometrial samples obtained across the menstrual cycle to investigate the regulation of temporal gene networks that control embryo implantation. We identify uniquely accessible chromatin regions in all major cellular constituents of the endometrium, delineate temporal patterns of coordinated chromatin remodeling in epithelial and stromal cells, and gain mechanistic insights into the emergence of a receptive state through integrated analysis of enriched transcription factor (TF) binding sites in dynamic chromatin regions, chromatin immunoprecipitation sequencing analyses, and gene expression data. We demonstrate that the implantation window coincides with pervasive cooption of transposable elements (TEs) into the regulatory chromatin landscape of decidualizing cells and expression of TE-derived transcripts in a spatially defined manner. Our data constitute a comprehensive map of the chromatin changes that control TF activities in a cycling endometrium at cellular resolution.

Single-Cell Analyses Reveal Necroptosis's Potential Role in Neuron Degeneration and Show Enhanced Neuron-Immune Cell Interaction in Parkinson's Disease Progression

Parkinson's disease (PD) is a common neuron degenerative disease among the old, characterized by uncontrollable movements and an impaired posture. Although widely investigated on its pathology and treatment, the disease remains incompletely understood. Single-cell RNA sequencing (scRNA-seq) has been applied to the area of PD, providing valuable data for related research. However, few works have taken deeper insights into the causes of neuron death and cell-cell interaction between the cell types in the brain. Our bioinformatics analyses revealed necroptosis-related genes (NRGs) enrichment in neuron degeneration and selecting the cells by NRGs levels showed two subtypes within the main degenerative cell types in the midbrain. NRG-low subtype was largely replaced by NRG-high subtype in the patients, indicating the striking change of cell state related to necroptosis in PD progression. Moreover, we carried out cell-cell interaction analyses between cell types and found that microglia (MG)'s interaction strength with glutamatergic neuron (GLU), GABAergic neuron (GABA), and dopaminergic neuron (DA) was significantly upregulated in PD. Also, MG show much stronger interaction with NRG-high subtypes and a stronger cell killing function in PD samples. Additionally, we identified CLDN11 as a novel interaction pattern specific to necroptosis neurons and MG. We also found LEF1 and TCF4 as key transcriptional regulators in neuron degeneration. These findings suggest that MG were significantly overactivated in PD patients to clear abnormal neurons, especially the NRG-high cells, explaining the neuron inflammation in PD. Our analyses provide insights into the causes of neuron death and inflammation in PD from single-cell resolution, which could be seriously considered in clinical trials.

Mutant GGGGCC RNA prevents YY1 from binding to Fuzzy promoter which stimulates Wnt/β-catenin pathway in C9ALS/FTD

The GGGGCC hexanucleotide repeat expansion mutation in the chromosome 9 open reading frame 72 (C9orf72) gene is a major genetic cause of amyotrophic lateral sclerosis and frontotemporal dementia (C9ALS/FTD). In this study, we demonstrate that the zinc finger (ZF) transcriptional regulator Yin Yang 1 (YY1) binds to the promoter region of the planar cell polarity gene Fuzzy to regulate its transcription. We show that YY1 interacts with GGGGCC repeat RNA via its ZF and that this interaction compromises the binding of YY1 to the FuzzyYY1 promoter sites, resulting in the downregulation of Fuzzy transcription. The decrease in Fuzzy protein expression in turn activates the canonical Wnt/β-catenin pathway and induces synaptic deficits in C9ALS/FTD neurons. Our findings demonstrate a C9orf72 GGGGCC RNA-initiated perturbation of YY1-Fuzzy transcriptional control that implicates aberrant Wnt/β-catenin signalling in C9ALS/FTD-associated neurodegeneration. This pathogenic cascade provides a potential new target for disease-modifying therapy.

An in vitro model of acute horizontal basal cell activation reveals dynamic gene regulatory networks underlying the acute activation phase

Horizontal basal cells (HBCs) activate only in response to severe olfactory epithelium (OE) injury. This activation is mediated by the loss of the transcription factor TP63. Using the compound phorbol 12-myristate 13-acetate (PMA), we find that we can model the process of acute HBC activation. First, we find that PMA treatment induces a rapid loss in TP63 protein and induces the expression of HOPX and the nuclear translocation of RELA, previously identified to mediate HBC activation. Using bulk RNA sequencing, we find that PMA-treated HBCs pass through various stages of acute activation identifiable by transcriptional regulatory signatures that mimic stages identified in vivo . These temporal stages are associated with varying degrees of engraftment and differentiation potential in transplantation assays. Together, this data shows that our model can model physiologically relevant features of HBC activation and identifies new candidates for mechanistic testing.

Differential developmental blueprints of organ-intrinsic nervous systems

The organ-intrinsic nervous system is a major interface between visceral organs and the brain, mediating important sensory and regulatory functions in the body-brain axis and serving as critical local processors for organ homeostasis. Molecularly, anatomically, and functionally, organ-intrinsic neurons are highly specialized for their host organs. However, the underlying mechanism that drives this specialization is largely unknown. Here, we describe the differential strategies utilized to achieve organ-specific organization between the enteric nervous system (ENS) 1 and the intrinsic cardiac nervous system (ICNS) 2 , a neuronal network essential for heart performance but poorly characterized. Integrating high-resolution whole-embryo imaging, single-cell genomics, spatial transcriptomics, proteomics, and bioinformatics, we uncover that unlike the ENS which is highly mobile and colonizes the entire gastrointestinal (GI) tract, the ICNS uses a rich set of extracellular matrix (ECM) genes that match with surrounding heart cells and an intermediate dedicated neuronal progenitor state to stabilize itself for a 'beads-on-the-necklace' organization on heart atria. While ICNS- and ENS-precursors are genetically similar, their differentiation paths are influenced by their host-organs, leading to distinct mature neuron types. Co-culturing ENS-precursors with heart cells shifts their identity towards the ICNS and induces the expression of heart-matching ECM genes. Our cross-organ study thus reveals fundamental principles for the maturation and specialization of organ-intrinsic neurons.

In vivo protein turnover rates in varying oxygen tensions nominate MYBBP1A as a mediator of the hyperoxia response

Oxygen deprivation and excess are both toxic. Thus, the body's ability to adapt to varying oxygen tensions is critical for survival. While the hypoxia transcriptional response has been well studied, the post-translational effects of oxygen have been underexplored. In this study, we systematically investigate protein turnover rates in mouse heart, lung, and brain under different inhaled oxygen tensions. We find that the lung proteome is the most responsive to varying oxygen tensions. In particular, several extracellular matrix (ECM) proteins are stabilized in the lung under both hypoxia and hyperoxia. Furthermore, we show that complex 1 of the electron transport chain is destabilized in hyperoxia, in accordance with the exacerbation of associated disease models by hyperoxia and rescue by hypoxia. Moreover, we nominate MYBBP1A as a hyperoxia transcriptional regulator, particularly in the context of rRNA homeostasis. Overall, our study highlights the importance of varying oxygen tensions on protein turnover rates and identifies tissue-specific mediators of oxygen-dependent responses.

HECW1 induces NCOA4-regulated ferroptosis in glioma through the ubiquitination and degradation of ZNF350

Tumor suppression by inducing NCOA4-mediated ferroptosis has been shown to be feasible in a variety of tumors, including gliomas. However, the regulatory mechanism of ferroptosis induced by NCOA4 in glioma has not been studied deeply. HECW1 and ZNF350 are involved in the biological processes of many tumors, but their specific effects and mechanisms on glioma are still unclear. In this study, we found that HECW1 decreased the survival rate of glioma cells and enhanced iron accumulation, lipid peroxidation, whereas ZNF350 showed the opposite effect. Mechanistically, HECW1 directly regulated the ubiquitination and degradation of ZNF350, eliminated the transcriptional inhibition of NCOA4 by ZNF350, and ultimately activated NCOA4-mediated iron accumulation, lipid peroxidation, and ferroptosis. We demonstrate that HECW1 induces ferroptosis and highlight the value of HECW1 and ZNF350 in the prognostic evaluation of patients with glioma. We also elucidate the mechanisms underlying the HECW1/ZNF350/NCOA4 axis and its regulation of ferroptosis. Our findings enrich the understanding of ferroptosis and provide potential treatment options for glioma patients.

Computational Characterization of Undifferentially Expressed Genes with Altered Transcription Regulation in Lung Cancer

A transcriptome profiles the expression levels of genes in cells and has accumulated a huge amount of public data. Most of the existing biomarker-related studies investigated the differential expression of individual transcriptomic features under the assumption of inter-feature independence. Many transcriptomic features without differential expression were ignored from the biomarker lists. This study proposed a computational analysis protocol (mqTrans) to analyze transcriptomes from the view of high-dimensional inter-feature correlations. The mqTrans protocol trained a regression model to predict the expression of an mRNA feature from those of the transcription factors (TFs). The difference between the predicted and real expression of an mRNA feature in a query sample was defined as the mqTrans feature. The new mqTrans view facilitated the detection of thirteen transcriptomic features with differentially expressed mqTrans features, but without differential expression in the original transcriptomic values in three independent datasets of lung cancer. These features were called dark biomarkers because they would have been ignored in a conventional differential analysis. The detailed discussion of one dark biomarker, GBP5, and additional validation experiments suggested that the overlapping long non-coding RNAs might have contributed to this interesting phenomenon. In summary, this study aimed to find undifferentially expressed genes with significantly changed mqTrans values in lung cancer. These genes were usually ignored in most biomarker detection studies of undifferential expression. However, their differentially expressed mqTrans values in three independent datasets suggested their strong associations with lung cancer.

Undifferentially Expressed CXXC5 as a Transcriptionally Regulatory Biomarker of Breast Cancer

This work hypothesizes that some genes undergo radically changed transcription regulations (TRs) in breast cancer (BC), but don't show differential expressions for unknown reasons. The TR of a gene is quantitatively formulated by a regression model between the expression of this gene and multiple transcription factors (TFs). The difference between the predicted and real expression levels of a gene in a query sample is defined as the mqTrans value of this gene, which quantitatively reflects its regulatory changes. This work systematically screens the undifferentially expressed genes with differentially expressed mqTrans values in 1036 samples across five datasets and three ethnic groups. This study calls the 25 genes satisfying the above hypothesis in at least four datasets as dark biomarkers, and the strong dark biomarker gene CXXC5 (CXXC Finger Protein 5) is even supported by all the five independent BC datasets. Although CXXC5 does not show differential expressions in BC, its transcription regulations show quantitative associations with BCs in diversified cohorts. The overlapping long noncoding RNAs (lncRNAs) may have contributed their transcripts to the expression miscalculations of dark biomarkers. The mqTrans analysis serves as a complementary view of the transcriptome-based detections of biomarkers that are ignored by many existing studies.

Homeobox-D 1 and FTO form a transcriptional-epigenetic feedback loop to promote head and neck cancer proliferation

Transcription factors (TFs) and N6-methyladenosine (m6A) modifiers are involved in tumor progression through transcriptional regulation and posttranscriptional regulation of genes, respectively. However, the crosstalk and role of these two types of gene expression regulators in head and neck squamous cell carcinoma (HNSC) remains poorly understood. In this study, we demonstrate that the TF homeobox-D1 (HOXD1) and the m6A demethylase fat mass and obesity-associated protein (FTO) form a positive feedback loop to promote cell proliferation and survival in HNSC. Clinically, HOXD1 expression is dysregulated in multiple cancer types and is associated with worse prognosis in patients with HNSC, stomach adenocarcinoma, uterine corpus endometrial carcinoma, and pheochromocytoma and paraganglioma. Mechanistically, FTO is overexpressed in HNSC tumor samples and positively regulates HOXD1 expression in an m6A-dependent manner. Functionally, deficiency of HOXD1 relieved the resistance of HNSC cells to apoptosis and arrested tumor cells at the G0/G1 phase, thereby inhibiting cell growth, whereas overexpression of HOXD1 caused the opposite effect. Furthermore, HOXD1 activates the transcription of the oncogenic factor FTO by directly targeting its promoter. Downregulation of FTO mimicked the biological effect of HOXD1 knockdown on HNSC. Importantly, overexpression of HOXD1 significantly rescued the proliferation inhibition and apoptosis promotion of HNSC cells induced by deficiency of FTO. Together, our findings reveal HOXD1 as a novel prognostic predictor and a potential target for HNSC, providing mechanistic insights into the role of the HOXD1-FTO circuit in this cancer.