AnimalTFDB 4.0: a comprehensive animal transcription factor database updated with variation and expression annotations

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.

Bioinformatics tools and resources for cancer and application

ABSTRACT

Tumor bioinformatics plays an important role in cancer research and precision medicine. The primary focus of traditional cancer research has been molecular and clinical studies of a number of fundamental pathways and genes. In recent years, driven by breakthroughs in high-throughput technologies, large-scale cancer omics data have accumulated rapidly. How to effectively utilize and share these data is particularly important. To address this crucial task, many computational tools and databases have been developed over the past few years. To help researchers quickly learn and understand the functions of these tools, in this review, we summarize publicly available bioinformatics tools and resources for pan-cancer multi-omics analysis, regulatory analysis of tumorigenesis, tumor treatment and prognosis, immune infiltration analysis, immune repertoire analysis, cancer driver gene and driver mutation analysis, and cancer single-cell analysis, which may further help researchers find more suitable tools for their research.

Different parts of the mussel Gigantidas haimaensis holobiont responded differently to deep-sea sampling stress

Acute environmental changes cause stress during conventional deep-sea biological sampling without in situ fixation and affect gene expressions of samples collected. However, the degree of influence and underlying mechanisms are hardly investigated. Here, we conducted comparative transcriptomic analyses between in situ and onboard fixed gills and between in situ and onboard fixed mantles of deep-sea mussel Gigantidas haimaensis to assess the effects of incidental sampling stress. Results showed that transcription, translation, and energy metabolism were upregulated in onboard fixed gills and mantles, thereby mobilizing rapid gene expression to tackle the stress. Autophagy and phagocytosis that related to symbiotic interactions between the host and endosymbiont were downregulated in the onboard fixed gills. These findings demonstrated that symbiotic gill and nonsymbiotic mantle responded differently to sampling stress, and symbiosis in the gill was perturbed. Further comparative metatranscriptomic analysis between in situ and onboard fixed gills revealed that stress response genes, peptidoglycan biosynthesis, and methane fixation were upregulated in the onboard fixed endosymbiotic Gammaproteobacteria inside the gills, implying that energy metabolism of the endosymbiont was increased to cope with sampling stress. Furthermore, comparative analysis between the mussel G. haimaensis and the limpet Bathyacmaea lactea transcriptomes resultedidentified six transcription factor orthologs upregulated in both onboard fixed mussel mantles and limpets, including sharply increased early growth response protein 1 and Kruppel-like factor 5. They potentially play key roles in initiating the response of sampled deep-sea macrobenthos to sampling stress. Our results clearly show that in situ fixed biological samples are vital for studying deep-sea environmental adaptation.

Evolutionary origin of the chordate nervous system revealed by amphioxus developmental trajectories

The common ancestor of all vertebrates had a highly sophisticated nervous system, but questions remain about the evolution of vertebrate neural cell types. The amphioxus, a chordate that diverged before the origin of vertebrates, can inform vertebrate evolution. Here we develop and analyse a single-cell RNA-sequencing dataset from seven amphioxus embryo stages to understand chordate cell type evolution and to study vertebrate neural cell type origins. We identified many new amphioxus cell types, including homologues to the vertebrate hypothalamus and neurohypophysis, rooting the evolutionary origin of these structures. On the basis of ancestor-descendant reconstruction of cell trajectories of the amphioxus and other species, we inferred expression dynamics of transcription factor genes throughout embryogenesis and identified three ancient developmental routes forming chordate neurons. We characterized cell specification at the mechanistic level and generated mutant lines to examine the function of five key transcription factors involved in neural specification. Our results show three developmental origins for the vertebrate nervous system: an anterior FoxQ2-dependent mechanism that is deeply conserved in invertebrates, a less-conserved route leading to more posterior neurons in the vertebrate spinal cord and a mechanism for specifying neuromesoderm progenitors that is restricted to chordates. The evolution of neuromesoderm progenitors may have led to a dramatic shift in posterior neural and mesodermal cell fate decisions and the body elongation process in a stem chordate.

In-silico analysis of cattle blood transcriptome to identify lncRNAs and their role during bovine tuberculosis

Long noncoding RNAs (lncRNAs) are RNA molecules with a length greater than 200 nucleotides that do not code for functional proteins. Although, genes play a vital role in immune response against a disease, it is less known that lncRNAs also contribute through gene regulation. Bovine tuberculosis is a significant zoonotic disease caused by Mycobacterium bovis (M. bovis) in cattle. Here, we report the in-silico analysis of the publicly available transcriptomic data of calves infected with M. bovis. A total of 51,812 lncRNAs were extracted across all the samples. A total of 216 genes and 260 lncRNAs were found to be differentially expressed across all the 4 conditions-infected vs uninfected at 8- and 20-week post-infection (WPI), 8 vs 20-WPI of both infected and uninfected. Gene Ontology and Functional annotation showed that 8 DEGs were annotated with immune system GOs and 2 DEGs with REACTOME immune system pathways. Co-expression analysis of DElncRNAs with DEGs revealed the involvement of lncRNAs with the genes annotated with Immune related GOs and pathways. Overall, our study sheds light on the dynamic transcriptomic changes in response to M. bovis infection, particularly highlighting the involvement of lncRNAs with immune-related genes. The identified immune pathways and gene-lncRNA interactions offer valuable insights for further research in understanding host-pathogen interactions and potential avenues for genetic improvement strategies in cattle.

Genome-wide analysis of the cytochrome P450 gene family in Pacific oyster Crassostrea gigas and their expression profiles during gonad development

The cytochrome P450 (CYP) gene superfamily plays a significant role in various physiological processes, producing different compounds such as hormones, fatty acids, and biomolecules. However, little information is known their roles during gonad development in Pacific oyster (Crassostrea gigas). In this study, total of 116 CgCYP (Crassostrea gigas cytochrome P450) genes were identified and their expression pattern was analyzed for the first time. The relative molecular weights of these CgCYP genes ranged from 63.52 to 113.41 kDa, and the length of encoded amino acids ranged from 103 to 993. And total 26 cis-acting elements of these CgCYP genes were identified. GO and KEGG enrichment analysis showed some CgCYP genes are essential for the metabolism of male and female sex hormones. Additionally, expression anslysis showed 69 CgCYP genes were over-expressed in early gonad development and triploid infertile individuals. More importantly, expression levels of CgCYP1, CgCYP15, CgCYP34, CgCYP46, CgCYP69, CgCYP87, CgCYP88, and CgCYP103, were found to be significantly higher in female gonad, suggesting their important roles in female gonad development. The results of this study will provide a better understanding of the CgCYP genes in the gonad development of Pacific oyster.

Polo-like kinase 1 inhibition modulates urinary tract smooth muscle contraction and bladder cell transcriptional programs

The serine/threonine kinase polo-like kinase 1 (PLK1) is a master regulator of cell proliferation and contraction, but its physiological role in the lower urinary tract is unknown. We utilized transcriptomic programs of human bladder smooth muscle cells (hBSMCs), 3D bladder spheroid viability assays, and human ureterovesical junction contractility measurements to elucidate the impacts of PLK1 inhibition. This work reveals PLK1 reduction with the selective inhibitor TAK-960 (500 nM) suppresses high K+-evoked contractions of human urinary smooth muscle ex vivo while decreasing urothelial cell viability. Transcriptomic analysis of hBSMCs treated with TAK-960 shows modulation of cell cycle and contraction pathways, specifically through altered expression of Cys2/His2-type zinc finger transcription factors. In bladder spheroids, PLK1 inhibition also suppresses smooth muscle contraction protein filamin. Taken together, these findings establish PLK1 is a critical governor of urinary smooth muscle contraction and urothelial proliferation with implications for lower urinary tract disorders. Targeting PLK1 pharmacologically may therefore offer therapeutic potential to ameliorate hypercontractility and aberrant growth. Further elucidation of PLK1 signaling networks promises new insights into pathogenesis and much needed treatment advances for debilitating urinary symptoms.

Expression genome-wide association study identifies key regulatory variants enriched with metabolic and immune functions in four porcine tissues

BACKGROUND

Integration of high throughput DNA genotyping and RNA-sequencing data enables the discovery of genomic regions that regulate gene expression, known as expression quantitative trait loci (eQTL). In pigs, efforts to date have been mainly focused on purebred lines for traits with commercial relevance as such growth and meat quality. However, little is known on genetic variants and mechanisms associated with the robustness of an animal, thus its overall health status. Here, the liver, lung, spleen, and muscle transcriptomes of 100 three-way crossbred female finishers were studied, with the aim of identifying novel eQTL regulatory regions and transcription factors (TFs) associated with regulation of porcine metabolism and health-related traits.

RESULTS

An expression genome-wide association study with 535,896 genotypes and the expression of 12,680 genes in liver, 13,310 genes in lung, 12,650 genes in spleen, and 12,595 genes in muscle resulted in 4,293, 10,630, 4,533, and 6,871 eQTL regions for each of these tissues, respectively. Although only a small fraction of the eQTLs were annotated as cis-eQTLs, these presented a higher number of polymorphisms per region and significantly stronger associations with their target gene compared to trans-eQTLs. Between 20 and 115 eQTL hotspots were identified across the four tissues. Interestingly, these were all enriched for immune-related biological processes. In spleen, two TFs were identified: ERF and ZNF45, with key roles in regulation of gene expression.

CONCLUSIONS

This study provides a comprehensive analysis with more than 26,000 eQTL regions identified that are now publicly available. The genomic regions and their variants were mostly associated with tissue-specific regulatory roles. However, some shared regions provide new insights into the complex regulation of genes and their interactions that are involved with important traits related to metabolism and immunity.

Tumour vasculature at single-cell resolution

Tumours can obtain nutrients and oxygen required to progress and metastasize through the blood supply1. Inducing angiogenesis involves the sprouting of established vessel beds and their maturation into an organized network2,3. Here we generate a comprehensive atlas of tumour vasculature at single-cell resolution, encompassing approximately 200,000 cells from 372 donors representing 31 cancer types. Trajectory inference suggested that tumour angiogenesis was initiated from venous endothelial cells and extended towards arterial endothelial cells. As neovascularization elongates (through angiogenic stages SI, SII and SIII), APLN+ tip cells at the SI stage (APLN+ TipSI) advanced to TipSIII cells with increased Notch signalling. Meanwhile, stalk cells, following tip cells, transitioned from high chemokine expression to elevated TEK (also known as Tie2) expression. Moreover, APLN+ TipSI cells not only were associated with disease progression and poor prognosis but also hold promise for predicting response to anti-VEGF therapy. Lymphatic endothelial cells demonstrated two distinct differentiation lineages: one responsible for lymphangiogenesis and the other involved in antigen presentation. In pericytes, endoplasmic reticulum stress was associated with the proangiogenic BASP1+ matrix-producing pericytes. Furthermore, intercellular communication analysis showed that neovascular endothelial cells could shape an immunosuppressive microenvironment conducive to angiogenesis. This study depicts the complexity of tumour vasculature and has potential clinical significance for anti-angiogenic therapy.

Deciphering lineage specification during early embryogenesis in mouse gastruloids using multilayered proteomics

Gastrulation is a critical stage in embryonic development during which the germ layers are established. Advances in sequencing technologies led to the identification of gene regulatory programs that control the emergence of the germ layers and their derivatives. However, proteome-based studies of early mammalian development are scarce. To overcome this, we utilized gastruloids and a multilayered mass spectrometry-based proteomics approach to investigate the global dynamics of (phospho) protein expression during gastruloid differentiation. Our findings revealed many proteins with temporal expression and unique expression profiles for each germ layer, which we also validated using single-cell proteomics technology. Additionally, we profiled enhancer interaction landscapes using P300 proximity labeling, which revealed numerous gastruloid-specific transcription factors and chromatin remodelers. Subsequent degron-based perturbations combined with single-cell RNA sequencing (scRNA-seq) identified a critical role for ZEB2 in mouse and human somitogenesis. Overall, this study provides a rich resource for developmental and synthetic biology communities endeavoring to understand mammalian embryogenesis.

Summary data-based Mendelian randomization and single-cell RNA sequencing analyses identify immune associations with low-level LGALS9 in sepsis

Sepsis is one of the major challenges in intensive care units, characterized by the complexity of the host immune status. To gain a deeper understanding of the pathogenesis of sepsis, it is crucial to study the phenotypic changes in immune cells and their underlying molecular mechanisms. We conducted Summary data-based Mendelian randomization analysis by integrating genome-wide association studies data for sepsis with expression quantitative trait locus data, revealing a significant decrease in the expression levels of 17 biomarkers in sepsis patients. Furthermore, based on single-cell RNA sequencing data, we elucidated potential molecular mechanisms at single-cell resolution and identified that LGALS9 inhibition in sepsis patients leads to the activation and differentiation of monocyte and T-cell subtypes. These findings are expected to assist researchers in gaining a more in-depth understanding of the immune dysregulation in sepsis.

Advances in computational and experimental approaches for deciphering transcriptional regulatory networks: Understanding the roles of cis-regulatory elements is essential, and recent research utilizing MPRAs, STARR-seq, CRISPR-Cas9, and machine learning has yielded valuable insights

Understanding the influence of cis-regulatory elements on gene regulation poses numerous challenges given complexities stemming from variations in transcription factor (TF) binding, chromatin accessibility, structural constraints, and cell-type differences. This review discusses the role of gene regulatory networks in enhancing understanding of transcriptional regulation and covers construction methods ranging from expression-based approaches to supervised machine learning. Additionally, key experimental methods, including MPRAs and CRISPR-Cas9-based screening, which have significantly contributed to understanding TF binding preferences and cis-regulatory element functions, are explored. Lastly, the potential of machine learning and artificial intelligence to unravel cis-regulatory logic is analyzed. These computational advances have far-reaching implications for precision medicine, therapeutic target discovery, and the study of genetic variations in health and disease.

KLF15-activated MARCH2 boosts cell proliferation and epithelial-mesenchymal transition and presents diagnostic significance for hepatocellular carcinoma

PURPOSE

Membrane associated ubiquitin ligase MARCH2 majorly involves in inflammation response and protein trafficking. However, its comprehensive role in hepatocellular carcinoma (HCC) is largely unknown.

METHODS

Firstly, multiple bioinformatic analyses were applied to determine MARCH2 mRNA level, its expression comparison in diverse molecular and immune subtypes, and diagnostic value in HCC. Subsequently, RNA-seq, real-time quantitative PCR, immunohistochemistry and cell proliferation assay are used to explore the epithelial-mesenchymal transition (EMT) and proliferation by gene-silencing or overexpressing in cultured HCC cells or in vivo xenograft. Moreover, dual luciferase reporter assay and immunoblotting are delved into verify the transcription factor that activating MARCH2 promoter.

RESULTS

Multiple bioinformatic analyses demonstrate that MARCH2 is upregulated in multiple cancer types and exhibits startling diagnostic value as well as distinct molecular and immune subtypes in HCC. RNA-seq analysis reveals MARCH2 may promote EMT, cell proliferation and migration in HepG2 cells. Furthermore, overexpression of MARCH2 triggers EMT and significantly enhances HCC cell migration, proliferation and colony formation in a ligase activity-dependent manner. Additionally, above observations are validated in the HepG2 mice xenografts. For up-stream mechanism, transcription factor KLF15 is highly expressed in HCC and activates MARCH2 expression.

CONCLUSION

KLF15 activated MARCH2 triggers EMT and serves as a fascinating biomarker for precise diagnosis of HCC. Consequently, MARCH2 emerges as a promising candidate for target therapy in cancer management.

Transcriptional response to an alternative diet on liver, muscle, and rumen of beef cattle

Feed cost represents a major economic determinant within cattle production, amounting to an estimated 75% of the total variable costs. Consequently, comprehensive approaches such as optimizing feed utilization through alternative feed sources, alongside the selection of feed-efficient animals, are of great significance. Here, we investigate the effect of two diets, traditional corn-grain fed and alternative by-product based, on 14 phenotypes related to feed, methane emission and production efficiency and on multi-tissue transcriptomics data from liver, muscle, and rumen wall, derived from 52 Nellore bulls, 26 on each diet. To this end, diets were contrasted at the level of phenotype, gene expression, and gene-phenotype network connectivity. As regards the phenotypic level, at a P value < 0.05, significant differences were found in favour of the alternative diet for average daily weight gain at finishing, dry matter intake at finishing, methane emission, carcass yield and subcutaneous fat thickness at the rib-eye muscle area. In terms of the transcriptional level of the 14,776 genes expressed across the examined tissues, we found 487, 484, and 499 genes differentially expressed due to diet in liver, muscle, and rumen, respectively (P value < 0.01). To explore differentially connected phenotypes across both diet-based networks, we focused on the phenotypes with the largest change in average number of connections within diets and tissues, namely methane emission and carcass yield, highlighting, in particular, gene expression changes involving SREBF2, and revealing the largest differential connectivity in rumen and muscle, respectively. Similarly, from examination of differentially connected genes across diets, the top-ranked most differentially connected regulators within each tissue were MEOX1, PTTG1, and BASP1 in liver, muscle, and rumen, respectively. Changes in gene co-expression patterns suggest activation or suppression of specific biological processes and pathways in response to dietary interventions, consequently impacting the phenotype. The identification of genes that respond differently to diets and their associated phenotypic effects serves as a crucial stepping stone for further investigations, aiming to build upon our discoveries. Ultimately, such advancements hold the promise of improving animal welfare, productivity, and sustainability in livestock farming.

Breed and timepoint-based analysis of chicken harderian gland transcriptome during Newcastle disease virus challenge

Introduction: Newcastle disease is a highly infectious disease caused by the Newcastle Disease Virus (NDV) and has a devastating financial impact on the global chicken industry. It was previously established that Leghorn and Fayoumi breeds of chicken exhibit variable resistance against NDV infection. The harderian gland is the less studied tissue of the chicken, known to play an essential role in the immune response. Methods: Our previous study, we reported differential gene expression and long noncoding RNAs (lncRNAs) between challenged and non-challenged chickens in the Harderian gland transcriptomic data. Now, we report the analysis of the same data studying the differential expression patterns between Leghorn and Fayoumi and between different timepoints during disease. First, the pipeline FHSpipe was used for identification of lncRNAs, followed by differential expression analysis by edgeR (GLM), functional annotation by OmicsBox, co-expression analysis using WGCNA and finally validation of selected lncRNAs and co-expressing genes using qRT-PCR. Results: Here, we observed that Leghorn showed a higher number of upregulated immune-related genes than Fayoumi in timepoint-based analysis, especially during the initial stages. Surprisingly, Fayoumi, being comparatively resistant, showed little difference between challenged and non-challenged conditions and different time points of the challenge. The breed-based analysis, which compared Leghorn with Fayoumi in both challenged and non-challenged conditions separately, identified several immune-related genes and positive co-expressing cis lncRNAs to be upregulated in Fayoumi when compared to Leghorn in both challenged and non-challenged conditions. Discussion: The current study shows that Leghorn, being comparatively more susceptible to NDV than Fayoumi, showed several immune-related genes and positive co-expressing cis lncRNAs upregulated in challenged Leghorn when compared to non-challenged Leghorn and also in different timepoints during challenge. While, breed-based analysis showed that there were more upregulated immune genes and positive cis-lncRNAs in Fayoumi than Leghorn. This result clearly shows that the differences in the expression of genes annotated with immune-related GO terms and pathways, i.e., immune-related genes and the co-expressing cis-lncRNAs between Leghorn and Fayoumi, and their role in the presence of differences in the resistance of Leghorn and Fayoumi chicken against NDV. Conclusion: These immune-genes and cis-lncRNAs could play a role in Fayoumi being comparatively more resistant to NDV than Leghorn. Our study elucidated the importance of lncRNAs during the host defense against NDV infection, paving the way for future research on the mechanisms governing the genetic improvement of chicken breeds.

Comprehensive transcriptome sequencing of silkworm Midguts: Uncovering extensive isoform diversity and alternative splicing in BmNPV-Sensitive and BmNPV-resistant strains

The silkworm, Bombyx mori, stands out as one of the few economically valuable insects within the realm of model organisms. However, Bombyx mori nucleopolyhedrovirus (BmNPV) poses a significant threat, decreasing the quality and quantity of silkworm cocoons. Over the past few decades, a multitude of researchers has delved into the mechanisms that underlie silkworm resistance to BmNPV, employing diverse methodologies and approaching the problem from various angles. Despite this extensive research, the role of alternative splicing (AS) in the silkworm's response to BmNPV infection has been largely unexplored. This study leveraged both third-generation (Oxford Nanopore Technologies) and second-generation (Illumina) high-throughput sequencing technologies to meticulously identify and analyze AS patterns in the context of BmNPV response, utilizing two distinct silkworm strains-the susceptible strain 306 and the resistant strain NB. Consequently, we identified five crucial genes (Dsclp, LOC692903, LOC101743583, LOC101742498, LOC101743809) that are linked to the response to BmNPV infection through AS and differential expression. Additionally, a thorough comparative analysis was conducted on their diverse transcriptomic expression profiles, including alternative polyadenylation, simple sequence repeats, and transcription factors.

Structure-based learning to predict and model protein-DNA interactions and transcription-factor co-operativity in cis-regulatory elements

Transcription factor (TF) binding is a key component of genomic regulation. There are numerous high-throughput experimental methods to characterize TF-DNA binding specificities. Their application, however, is both laborious and expensive, which makes profiling all TFs challenging. For instance, the binding preferences of ∼25% human TFs remain unknown; they neither have been determined experimentally nor inferred computationally. We introduce a structure-based learning approach to predict the binding preferences of TFs and the automated modelling of TF regulatory complexes. We show the advantage of using our approach over the classical nearest-neighbor prediction in the limits of remote homology. Starting from a TF sequence or structure, we predict binding preferences in the form of motifs that are then used to scan a DNA sequence for occurrences. The best matches are either profiled with a binding score or collected for their subsequent modeling into a higher-order regulatory complex with DNA. Co-operativity is modelled by: (i) the co-localization of TFs and (ii) the structural modeling of protein-protein interactions between TFs and with co-factors. We have applied our approach to automatically model the interferon-β enhanceosome and the pioneering complexes of OCT4, SOX2 (or SOX11) and KLF4 with a nucleosome, which are compared with the experimentally known structures.

The non-canonical poly(A) polymerase FAM46C promotes erythropoiesis

The post-transcriptional regulation of mRNA is a crucial component of gene expression. The disruption of this process has detrimental effects on the normal development and gives rise to various diseases. Searching for novel post-transcriptional regulators and exploring their roles are essential for understanding development and disease. Through a multimodal analysis of red blood cell trait genome-wide association studies (GWAS) and transcriptomes of erythropoiesis, we identify FAM46C, a non-canonical RNA poly(A) polymerase, as a necessary factor for proper red blood cell development. FAM46C is highly expressed in the late stages of the erythroid lineage, and its developmental upregulation is controlled by an erythroid-specific enhancer. We demonstrate that FAM46C stabilizes mRNA and regulates erythroid differentiation in a polymerase activity-dependent manner. Furthermore, we identify transcripts of lysosome and mitochondria components as highly confident in vivo targets of FAM46C, which aligns with the need of maturing red blood cells for substantial clearance of organelles and maintenance of cellular redox homeostasis. In conclusion, our study unveils a unique role of FAM46C in positively regulating lysosome and mitochondria components, thereby promoting erythropoiesis.

Trans-omic profiling uncovers molecular controls of early human cerebral organoid formation

Defining the molecular networks orchestrating human brain formation is crucial for understanding neurodevelopment and neurological disorders. Challenges in acquiring early brain tissue have incentivized the use of three-dimensional human pluripotent stem cell (hPSC)-derived neural organoids to recapitulate neurodevelopment. To elucidate the molecular programs that drive this highly dynamic process, here, we generate a comprehensive trans-omic map of the phosphoproteome, proteome, and transcriptome of the exit of pluripotency and neural differentiation toward human cerebral organoids (hCOs). These data reveal key phospho-signaling events and their convergence on transcriptional factors to regulate hCO formation. Comparative analysis with developing human and mouse embryos demonstrates the fidelity of our hCOs in modeling embryonic brain development. Finally, we demonstrate that biochemical modulation of AKT signaling can control hCO differentiation. Together, our data provide a comprehensive resource to study molecular controls in human embryonic brain development and provide a guide for the future development of hCO differentiation protocols.

GATA binding protein 2 mediated ankyrin repeat domain containing 26 high expression in myeloid-derived cell lines

BACKGROUND

Thrombocytopenia 2, an autosomal dominant inherited disease characterized by moderate thrombocytopenia, predisposition to myeloid malignancies and normal platelet size and function, can be caused by 5'-untranslated region (UTR) point mutations in ankyrin repeat domain containing 26 (ANKRD26). Runt related transcription factor 1 (RUNX1) and friend leukemia integration 1 (FLI1) have been identified as negative regulators of ANKRD26. However, the positive regulators of ANKRD26 are still unknown.

AIM

To prove the positive regulatory effect of GATA binding protein 2 (GATA2) on ANKRD26 transcription.

METHODS

Human induced pluripotent stem cells derived from bone marrow (hiPSC-BM) and urothelium (hiPSC-U) were used to examine the ANKRD26 expression pattern in the early stage of differentiation. Then, transcriptome sequencing of these iPSCs and three public transcription factor (TF) databases (Cistrome DB, animal TFDB and ENCODE) were used to identify potential TF candidates for ANKRD26. Furthermore, overexpression and dual-luciferase reporter experiments were used to verify the regulatory effect of the candidate TFs on ANKRD26. Moreover, using the GENT2 platform, we analyzed the relationship between ANKRD26 expression and overall survival in cancer patients.

RESULTS

In hiPSC-BMs and hiPSC-Us, we found that the transcription levels of ANKRD26 varied in the absence of RUNX1 and FLI1. We sequenced hiPSC-BM and hiPSC-U and identified 68 candidate TFs for ANKRD26. Together with three public TF databases, we found that GATA2 was the only candidate gene that could positively regulate ANKRD26. Using dual-luciferase reporter experiments, we showed that GATA2 directly binds to the 5'-UTR of ANKRD26 and promotes its transcription. There are two identified binding sites of GATA2 that are located 2 kb upstream of the TSS of ANKRD26. In addition, we discovered that high ANKRD26 expression is always related to a more favorable prognosis in breast and lung cancer patients.

CONCLUSION

We first discovered that the transcription factor GATA2 plays a positive role in ANKRD26 transcription and identified its precise binding sites at the promoter region, and we revealed the importance of ANKRD26 in many tissue-derived cancers.

Hematopoietic stem cell niche generation and maintenance are distinguishable by an epitranscriptomic program

The niche is typically considered as a pre-established structure sustaining stem cells. Therefore, the regulation of its formation remains largely unexplored. Whether distinct molecular mechanisms control the establishment versus maintenance of a stem cell niche is unknown. To address this, we compared perinatal and adult bone marrow mesenchymal stromal cells (MSCs), a key component of the hematopoietic stem cell (HSC) niche. MSCs exhibited enrichment in genes mediating m6A mRNA methylation at the perinatal stage and downregulated the expression of Mettl3, the m6A methyltransferase, shortly after birth. Deletion of Mettl3 from developing MSCs but not osteoblasts led to excessive osteogenic differentiation and a severe HSC niche formation defect, which was significantly rescued by deletion of Klf2, an m6A target. In contrast, deletion of Mettl3 from MSCs postnatally did not affect HSC niche. Stem cell niche generation and maintenance thus depend on divergent molecular mechanisms, which may be exploited for regenerative medicine.

Single nucleotide polymorphisms in the KRT82 promoter region modulate irregular thickening and patchiness in the dorsal skin of New Zealand rabbits

BACKGROUND

While rabbits are used as models in skin irritation tests, the presence of irregular patches and thickening on the dorsal skin can affect precise evaluation. In this study, genes associated with patchiness or non-patchiness on the dorsal skin of New Zealand rabbits were investigated to identify potential regulators of the patchiness phenotype.

RESULTS

The results showed that parameters associated with hair follicles (HFs), such as HF density, skin thickness, and HF depth, were augmented in rabbits with the patchiness phenotype relative to the non-patchiness phenotype. A total of 592 differentially expressed genes (DEGs) were identified between the two groups using RNA-sequencing. These included KRT72, KRT82, KRT85, FUT8, SOX9, and WNT5B. The functions of the DEGs were investigated by GO and KEGG enrichment analyses. A candidate gene, KRT82, was selected for further molecular function verification. There was a significant positive correlation between KRT82 expression and HF-related parameters, and KRT82 overexpression and knockdown experiments with rabbit dermal papilla cells (DPCs) showed that it regulated genes related to skin and HF growth and development. Investigation of single nucleotide polymorphisms (SNPs) in the exons and promoter region of KRT82 identified four SNPs in the promoter region but none in the exons. The G.-631G > T, T.-696T > C, G.-770G > T and A.-873 A > C alleles conformed to the Hardy - Weinberg equilibrium, and three identified haplotypes showed linkage disequilibrium. Luciferase reporter assays showed that the core promoter region of KRT82 was located in the - 600 to - 1200 segment, in which the four SNPs were located.

CONCLUSIONS

The morphological characteristics of the patchiness phenotype were analyzed in New Zealand rabbits and DEGs associated with this phenotype were identified by RNA-sequencing. The biological functions of the gene KRT82 associated with this phenotype were analyzed, and four SNPs were identified in the promoter region of the gene. These findings suggest that KRT82 may be a potential biomarker for the breeding of experimental New Zealand rabbits.

Deciphering infected cell types, hub gene networks and cell-cell communication in infectious bronchitis virus via single-cell RNA sequencing

Infectious bronchitis virus (IBV) is a coronavirus that infects chickens, which exhibits a broad tropism for epithelial cells, infecting the tracheal mucosal epithelium, intestinal mucosal epithelium, and renal tubular epithelial cells. Utilizing single-cell RNA sequencing (scRNA-seq), we systematically examined cells in renal, bursal, and tracheal tissues following IBV infection and identified tissue-specific molecular markers expressed in distinct cell types. We evaluated the expression of viral RNA in diverse cellular populations and subsequently ascertained that distal tubules and collecting ducts within the kidney, bursal mucosal epithelial cells, and follicle-associated epithelial cells exhibit susceptibility to IBV infection through immunofluorescence. Furthermore, our findings revealed an upregulation in the transcription of proinflammatory cytokines IL18 and IL1B in renal macrophages as well as increased expression of apoptosis-related gene STAT in distal tubules and collecting duct cells upon IBV infection leading to renal damage. Cell-to-cell communication unveiled potential interactions between diverse cell types, as well as upregulated signaling pathways and key sender-receiver cell populations after IBV infection. Integrating single-cell data from all tissues, we applied weighted gene co-expression network analysis (WGCNA) to identify gene modules that are specifically expressed in different cell populations. Based on the WGCNA results, we identified seven immune-related gene modules and determined the differential expression pattern of module genes, as well as the hub genes within these modules. Our comprehensive data provides valuable insights into the pathogenesis of IBV as well as avian antiviral immunology.

Defining a Haplotype Encompassing the LCORL-NCAPG Locus Associated with Increased Lean Growth in Beef Cattle

Numerous studies have shown genetic variation at the LCORL-NCAPG locus is strongly associated with growth traits in beef cattle. However, a causative molecular variant has yet to be identified. To define all possible candidate variants, 34 Charolais-sired calves were whole-genome sequenced, including 17 homozygous for a long-range haplotype associated with increased growth (QQ) and 17 homozygous for potential ancestral haplotypes for this region (qq). The Q haplotype was refined to an 814 kb region between chr6:37,199,897-38,014,080 and contained 218 variants not found in qq individuals. These variants include an insertion in an intron of NCAPG, a previously documented mutation in NCAPG (rs109570900), two coding sequence mutations in LCORL (rs109696064 and rs384548488), and 15 variants located within ATAC peaks that were predicted to affect transcription factor binding. Notably, rs384548488 is a frameshift variant likely resulting in loss of function for long isoforms of LCORL. To test the association of the coding sequence variants of LCORL with phenotype, 405 cattle from five populations were genotyped. The two variants were in complete linkage disequilibrium. Statistical analysis of the three populations that contained QQ animals revealed significant (p < 0.05) associations with genotype and birth weight, live weight, carcass weight, hip height, and average daily gain. These findings affirm the link between this locus and growth in beef cattle and describe DNA variants that define the haplotype. However, further studies will be required to define the true causative mutation.

The chromatin accessibility dynamics during cell fate specifications in zebrafish early embryogenesis

Chromatin accessibility plays a critical role in the regulation of cell fate decisions. Although gene expression changes have been extensively profiled at the single-cell level during early embryogenesis, the dynamics of chromatin accessibility at cis-regulatory elements remain poorly studied. Here, we used a plate-based single-cell ATAC-seq method to profile the chromatin accessibility dynamics of over 10 000 nuclei from zebrafish embryos. We investigated several important time points immediately after zygotic genome activation (ZGA), covering key developmental stages up to dome. The results revealed key chromatin signatures in the first cell fate specifications when cells start to differentiate into enveloping layer (EVL) and yolk syncytial layer (YSL) cells. Finally, we uncovered many potential cell-type specific enhancers and transcription factor motifs that are important for the cell fate specifications.

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

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

Loci associated with cave-derived traits concentrate in specific regions of the Mexican cavefish genome

A major goal of modern evolutionary biology is connecting phenotypic evolution with its underlying genetic basis. The Mexican cavefish (Astyanax mexicanus), a characin fish species comprised of a surface ecotype and a cave-derived ecotype, is well suited as a model to study the genetic mechanisms underlying adaptation to extreme environments. Here we map 206 previously published quantitative trait loci (QTL) for cave-derived traits in A. mexicanus to the newest version of the surface fish genome assembly, AstMex3. This analysis revealed that QTL cluster in the genome more than expected by chance, and this clustering is not explained by the distribution of genes in the genome. To investigate whether certain characteristics of the genome facilitate phenotypic evolution, we tested whether genomic characteristics, such as highly mutagenic CpG sites, are reliable predictors of the sites of trait evolution but did not find any significant trends. Finally, we combined the QTL map with previously collected expression and selection data to identify a list of 36 candidate genes that may underlie the repeated evolution of cave phenotypes, including rgrb which is predicted to be involved in phototransduction. We found this gene has disrupted exons in all non-hybrid cave populations but intact reading frames in surface fish. Overall, our results suggest specific "evolutionary hotspots" in the genome may play significant roles in driving adaptation to the cave environment in Astyanax mexicanus and demonstrate how this compiled dataset can facilitate our understanding of the genetic basis of repeated evolution in the Mexican cavefish.

Unraveling Axonal Transcriptional Landscapes: Insights from iPSC-Derived Cortical Neurons and Implications for Motor Neuron Degeneration

Neuronal function and pathology are deeply influenced by the distinct molecular profiles of the axon and soma. Traditional studies have often overlooked these differences due to the technical challenges of compartment specific analysis. In this study, we employ a robust RNA-sequencing (RNA-seq) approach, using microfluidic devices, to generate high-quality axonal transcriptomes from iPSC-derived cortical neurons (CNs). We achieve high specificity of axonal fractions, ensuring sample purity without contamination. Comparative analysis revealed a unique and specific transcriptional landscape in axonal compartments, characterized by diverse transcript types, including protein-coding mRNAs, ribosomal proteins (RPs), mitochondrial-encoded RNAs, and long non-coding RNAs (lncRNAs). Previous works have reported the existence of transcription factors (TFs) in the axon. Here, we detect a subset of previously unreported TFs specific to the axon and indicative of their active participation in transcriptional regulation. To investigate transcripts and pathways essential for central motor neuron (MN) degeneration and maintenance we analyzed KIF1C-knockout (KO) CNs, modeling hereditary spastic paraplegia (HSP), a disorder associated with prominent length-dependent degeneration of central MN axons. We found that several key factors crucial for survival and health were absent in KIF1C-KO axons, highlighting a possible role of these also in other neurodegenerative diseases. Taken together, this study underscores the utility of microfluidic devices in studying compartment-specific transcriptomics in human neuronal models and reveals complex molecular dynamics of axonal biology. The impact of KIF1C on the axonal transcriptome not only deepens our understanding of MN diseases but also presents a promising avenue for exploration of compartment specific disease mechanisms.

Gastrulation-stage gene expression in Nipbl+/- mouse embryos foreshadows the development of syndromic birth defects

In animal models, Nipbl deficiency phenocopies gene expression changes and birth defects seen in Cornelia de Lange syndrome, the most common cause of which is Nipbl haploinsufficiency. Previous studies in Nipbl+/- mice suggested that heart development is abnormal as soon as cardiogenic tissue is formed. To investigate this, we performed single-cell RNA sequencing on wild-type and Nipbl+/- mouse embryos at gastrulation and early cardiac crescent stages. Nipbl+/- embryos had fewer mesoderm cells than wild-type and altered proportions of mesodermal cell subpopulations. These findings were associated with underexpression of genes implicated in driving specific mesodermal lineages. In addition, Nanog was found to be overexpressed in all germ layers, and many gene expression changes observed in Nipbl+/- embryos could be attributed to Nanog overexpression. These findings establish a link between Nipbl deficiency, Nanog overexpression, and gene expression dysregulation/lineage misallocation, which ultimately manifest as birth defects in Nipbl+/- animals and Cornelia de Lange syndrome.

CircAKT3 alleviates postoperative cognitive dysfunction by stabilizing the feedback cycle of miR-106a-5p/HDAC4/MEF2C axis in hippocampi of aged mice

Circular RNAs (circRNAs) have garnered significant attention in the field of neurodegenerative diseases including Alzheimer's diseases due to their covalently closed loop structure. However, the involvement of circRNAs in postoperative cognitive dysfunction (POCD) is still largely unexplored. To identify the genes differentially expressed between non-POCD (NPOCD) and POCD mice, we conducted the whole transcriptome sequencing initially in this study. According to the expression profiles, we observed that circAKT3 was associated with hippocampal neuronal apoptosis in POCD mice. Moreover, we found that circAKT3 overexpression reduced apoptosis of hippocampal neurons and alleviated POCD. Subsequently, through bioinformatics analysis, our data showed that circAKT3 overexpression in vitro and in vivo elevated the abundance of miR-106a-5p significantly, resulting in a decrease of HDAC4 protein and an increase of MEF2C protein. Additionally, this effect of circAKT3 was blocked by miR-106a-5p inhibitor. Interestingly, MEF2C could activate the transcription of miR-106a-5p promoter and form a positive feedback loop. Therefore, our findings revealed more potential modulation ways between circRNA-miRNA and miRNA-mRNA, providing different directions and targets for preclinical studies of POCD.

MKRN1 regulates the expression profiles and transcription factor activity in HeLa cells inhibition suppresses cervical cancer cell progression

Cervical cancer is one of the most common gynecologic malignancies worldwide, necessitating the identification of novel biomarkers and therapeutic targets. This study aimed to investigate the significance of MKRN1 in cervical cancer and explore its potential as a diagnostic marker and therapeutic target. The results indicated that MKRN1 expression was up-regulated in cervical cancer tissues and correlated with advanced tumor stage, higher grade, and poor patient survival. Functional studies demonstrated that targeting MKRN1 effectively inhibited cell proliferation, migration, and invasion, highlighting its critical role in tumor progression and metastasis. Moreover, the knockdown of MKRN1 resulted in altered expression patterns of six transcription factor-encoding genes, revealing its involvement in gene regulation. Co-expression network analysis unveiled complex regulatory mechanisms underlying the effects of MKRN1 knockdown on gene expression. Furthermore, the results suggested that MKRN1 might serve as a diagnostic marker for personalized treatment strategies and a therapeutic target to inhibit tumor growth, metastasis, and overcome drug resistance. The development of MKRN1-targeted interventions might hold promise for advancing personalized medicine approaches in cervical cancer treatment. Further research is warranted to validate these findings, elucidate underlying mechanisms, and translate these insights into improved management and outcomes for cervical cancer patients.

Mechanisms of harmful effects of Microcystis aeruginosa on a brackish water organism Moina mongolica based on physiological and transcriptomic responses

To investigate the detrimental impacts of cyanobacterial bloom, specifically Microcystis aeruginosa, on brackish water ecosystems, the study used Moina mongolica, a cladoceran species, as the test organism. In a chronic toxicology experiment, the survival and reproductive rates of M. mongolica were assessed under M. aeruginosa stress. It was observed that the survival rate of M. mongolica fed with M. aeruginosa significantly decreased with time and their reproduction rate dropped to zero, while the control group remained maintained stable and normal reproduction. To further explore the underlying molecular mechanisms of the effects of M. aeruginosa on M. mongolica, we conducted a transcriptomic analysis on newly hatched M. mongolica cultured under different food conditions for 24 h. The results revealed significant expression differences in 572 genes, with 233 genes significantly up-regulated and 339 genes significantly down-regulated. Functional analysis of these differentially expressed genes identified six categories of physiological functional changes, including nutrition and metabolism, oxidative phosphorylation, neuroimmunology, cuticle and molting, reproduction, and programmed cell death. Based on these findings, we outlined the basic mechanisms of microcystin toxicity. The discovery provides critical insights into the mechanisms of Microcystis toxicity on organisms and explores the response mechanisms of cladocerans under the stress of Microcystis.

The Mediterranean mussel Mytilus galloprovincialis: a novel model for developmental studies in mollusks

A model organism in developmental biology is defined by its experimental amenability and by resources created for the model system by the scientific community. For the most powerful invertebrate models, the combination of both has already yielded a thorough understanding of developmental processes. However, the number of developmental model systems is still limited, and their phylogenetic distribution heavily biased. Members of one of the largest animal lineages, the Spiralia, for example, have long been neglected. In order to remedy this shortcoming, we have produced a detailed developmental transcriptome for the bivalve mollusk Mytilus galloprovincialis, and have expanded the list of experimental protocols available for this species. Our high-quality transcriptome allowed us to identify transcriptomic signatures of developmental progression and to perform a first comparison with another bivalve mollusk: the Pacific oyster Crassostrea gigas. To allow co-labelling studies, we optimized and combined protocols for immunohistochemistry and hybridization chain reaction to create high-resolution co-expression maps of developmental genes. The resources and protocols described here represent an enormous boost for the establishment of Mytilus galloprovincialis as an alternative model system in developmental biology.

Gastrulation-stage gene expression in Nipbl +/- mouse embryos foreshadows the development of syndromic birth defects

In animal models, Nipbl-deficiency phenocopies gene expression changes and birth defects seen in Cornelia de Lange Syndrome (CdLS), the most common cause of which is Nipbl-haploinsufficiency. Previous studies in Nipbl+/- mice suggested that heart development is abnormal as soon as cardiogenic tissue is formed. To investigate this, we performed single-cell RNA-sequencing on wildtype (WT) and Nipbl+/- mouse embryos at gastrulation and early cardiac crescent stages. Nipbl+/- embryos had fewer mesoderm cells than WT and altered proportions of mesodermal cell subpopulations. These findings were associated with underexpression of genes implicated in driving specific mesodermal lineages. In addition, Nanog was found to be overexpressed in all germ layers, and many gene expression changes observed in Nipbl+/- embryos could be attributed to Nanog overexpression. These findings establish a link between Nipbl-deficiency, Nanog overexpression, and gene expression dysregulation/lineage misallocation, which ultimately manifest as birth defects in Nipbl+/- animals and CdLS.

Teaser

Gene expression changes during gastrulation of Nipbl-deficient mice shed light on early origins of structural birth defects.

Combined Transcriptomics and Metabolomics Identify Regulatory Mechanisms of Porcine Vertebral Chondrocyte Development In Vitro

Porcine body length is closely related to meat production, growth, and reproductive performance, thus playing a key role in the profitability of the pork industry. Cartilage development is critical to longitudinal elongation of individual vertebrae. This study isolated primary porcine vertebral chondrocytes (PVCs) to clarify the complex mechanisms of elongation. We used transcriptome and target energy metabolome technologies to confirm crucial genes and metabolites in primary PVCs at different differentiation stages (0, 4, 8, and 12 days). Pairwise comparisons of the four stages identified 4566 differentially expressed genes (DEGs). Time-series gene cluster and functional analyses of these DEGs revealed four clusters related to metabolic processes, cartilage development, vascular development, and cell cycle regulation. We constructed a transcriptional regulatory network determining chondrocyte maturation. The network indicated that significantly enriched transcription factor (TF) families, including zf-C2H2, homeobox, TF_bZIP, and RHD, are important in cell cycle and differentiation processes. Further, dynamic network biomarker (DNB) analysis revealed that day 4 was the tipping point for chondrocyte development, consistent with morphological and metabolic changes. We found 24 DNB DEGs, including the TFs NFATC2 and SP7. Targeted energy metabolome analysis showed that most metabolites were elevated throughout chondrocyte development; notably, 16 differentially regulated metabolites (DRMs) were increased at three time points after cell differentiation. In conclusion, integrated metabolome and transcriptome analyses highlighted the importance of amino acid biosynthesis in chondrocyte development, with coordinated regulation of DEGs and DRMs promoting PVC differentiation via glucose oxidation. These findings reveal the regulatory mechanisms underlying PVC development and provide an important theoretical reference for improving pork production.

Peak Scores Significantly Depend on the Relationships between Contextual Signals in ChIP-Seq Peaks

Chromatin immunoprecipitation followed by massively parallel DNA sequencing (ChIP-seq) is a central genome-wide method for in vivo analyses of DNA-protein interactions in various cellular conditions. Numerous studies have demonstrated the complex contextual organization of ChIP-seq peak sequences and the presence of binding sites for transcription factors in them. We assessed the dependence of the ChIP-seq peak score on the presence of different contextual signals in the peak sequences by analyzing these sequences from several ChIP-seq experiments using our fully enumerative GPU-based de novo motif discovery method, Argo_CUDA. Analysis revealed sets of significant IUPAC motifs corresponding to the binding sites of the target and partner transcription factors. For these ChIP-seq experiments, multiple regression models were constructed, demonstrating a significant dependence of the peak scores on the presence in the peak sequences of not only highly significant target motifs but also less significant motifs corresponding to the binding sites of the partner transcription factors. A significant correlation was shown between the presence of the target motifs FOXA2 and the partner motifs HNF4G, which found experimental confirmation in the scientific literature, demonstrating the important contribution of the partner transcription factors to the binding of the target transcription factor to DNA and, consequently, their important contribution to the peak score.

KnockTF 2.0: a comprehensive gene expression profile database with knockdown/knockout of transcription (co-)factors in multiple species

Transcription factors (TFs), transcription co-factors (TcoFs) and their target genes perform essential functions in diseases and biological processes. KnockTF 2.0 (http://www.licpathway.net/KnockTF/index.html) aims to provide comprehensive gene expression profile datasets before/after T(co)F knockdown/knockout across multiple tissue/cell types of different species. Compared with KnockTF 1.0, KnockTF 2.0 has the following improvements: (i) Newly added T(co)F knockdown/knockout datasets in mice, Arabidopsis thaliana and Zea mays and also an expanded scale of datasets in humans. Currently, KnockTF 2.0 stores 1468 manually curated RNA-seq and microarray datasets associated with 612 TFs and 172 TcoFs disrupted by different knockdown/knockout techniques, which are 2.5 times larger than those of KnockTF 1.0. (ii) Newly added (epi)genetic annotations for T(co)F target genes in humans and mice, such as super-enhancers, common SNPs, methylation sites and chromatin interactions. (iii) Newly embedded and updated search and analysis tools, including T(co)F Enrichment (GSEA), Pathway Downstream Analysis and Search by Target Gene (BLAST). KnockTF 2.0 is a comprehensive update of KnockTF 1.0, which provides more T(co)F knockdown/knockout datasets and (epi)genetic annotations across multiple species than KnockTF 1.0. KnockTF 2.0 facilitates not only the identification of functional T(co)Fs and target genes but also the investigation of their roles in the physiological and pathological processes.

scGRN: a comprehensive single-cell gene regulatory network platform of human and mouse

Gene regulatory networks (GRNs) are interpretable graph models encompassing the regulatory interactions between transcription factors (TFs) and their downstream target genes. Making sense of the topology and dynamics of GRNs is fundamental to interpreting the mechanisms of disease etiology and translating corresponding findings into novel therapies. Recent advances in single-cell multi-omics techniques have prompted the computational inference of GRNs from single-cell transcriptomic and epigenomic data at an unprecedented resolution. Here, we present scGRN (https://bio.liclab.net/scGRN/), a comprehensive single-cell multi-omics gene regulatory network platform of human and mouse. The current version of scGRN catalogs 237 051 cell type-specific GRNs (62 999 692 TF-target gene pairs), covering 160 tissues/cell lines and 1324 single-cell samples. scGRN is the first resource documenting large-scale cell type-specific GRN information of diverse human and mouse conditions inferred from single-cell multi-omics data. We have implemented multiple online tools for effective GRN analysis, including differential TF-target network analysis, TF enrichment analysis, and pathway downstream analysis. We also provided details about TF binding to promoters, super-enhancers and typical enhancers of target genes in GRNs. Taken together, scGRN is an integrative and useful platform for searching, browsing, analyzing, visualizing and downloading GRNs of interest, enabling insight into the differences in regulatory mechanisms across diverse conditions.

TFCheckpoint database update, a cross-referencing system for transcription factors from human, mouse and rat

Prior knowledge about DNA-binding transcription factors (dbTFs), transcription co-regulators (coTFs) and general transcriptional factors (GTFs) is crucial for the study and understanding of the regulation of transcription. This is reflected by the many publications and database resources describing knowledge about TFs. We previously launched the TFCheckpoint database, an integrated resource focused on human, mouse and rat dbTFs, providing users access to a comprehensive overview of these proteins. Here, we describe TFCheckpoint 2.0 (https://www.tfcheckpoint.org/index.php), comprising 13 collections of dbTFs, coTFs and GTFs. TFCheckpoint 2.0 provides an easy and versatile cross-referencing system for users to view and download collections that may otherwise be cumbersome to find, compare and retrieve.

Computational identification of key transcription factors for embryonic and postnatal Sox2+ dental epithelial stem cell

While many reptiles can replace their tooth throughout life, human loss the tooth replacement capability after formation of the permanent teeth. It was thought that the difference in tooth regeneration capability depends on the persistence of a specialized dental epithelial structure, the dental lamina that contains dental epithelial stem cells (DESC). Currently, we know very little about DESC such as what genes are expressed or its chromatin accessibility profile. Multiple markers of DESC have been proposed such as Sox2 and Lgr5 . Few single cell RNA-seq experiments have been performed previously, but no obvious DESC cluster was identified in these scRNA-seq datasets, possible due to that the expression level of DESC markers such as Sox2 and Lgr5 is too low or the percentage of DESC is too low in whole tooth. We utilize a mouse line Sox2-GFP to enrich Sox2+ DESC and use Smart-Seq2 protocol and ATAC-seq protocol to generate transcriptome profile and chromatin accessibility profile of P2 Sox2+ DESC. Additionally, we generate transcriptome profile and chromatin accessibility profile of E11.5 Sox2+ dental lamina cells. With transcriptome profile and chromatin accessibility profile, we systematically identify potential key transcription factors for E11.5 Sox2+ cells and P2 Sox2+ cells. We identified transcription factors including Pitx2, Id3, Pitx1, Tbx1, Trp63, Nkx2-3, Grhl3, Dlx2, Runx1, Nfix, Zfp536 , etc potentially formed the core transcriptional regulatory networks of Sox2+ DESC in both embryonic and postnatal stages.

Novel enhancer mediates the RPL36A-HNRNPH2 readthrough loci and GLA gene expressions associated with fabry disease

Fabry disease (FD) is a rare genetic condition caused by mutations in the GLA gene, located on the X chromosome in the RPL36-HNRNPH2 readthrough genomic region. This gene produces an enzyme called alpha-galactosidase A (α-Gal A). When the enzyme does not function properly due to the mutations, it causes harmful substances called globotriaosylceramide (Gb3) and globotriaosylsphingosine (lyso-Gb3) to build up in the body's lysosomes. This accumulation can damage the kidneys, heart, eyes, and nervous system. Recent studies have shown that the RPL36A-HNRNPH2 readthrough loci, which include RPL36A and HNRNPH2 genes, as well as the regulatory sequence known as the GLA-HNRNPH2 bidirectional promoter, may also play a role in FD. However, the involvement of enhancer RNAs (eRNAs) in FD is still poorly understood despite their known role in various diseases. To investigate this further, we studied an RPL36A enhancer called GH0XJ101390 and showed its genomic setting in the RPL36-HNRNPH2 readthrough region; the eRNA is rich in Homotypic Clusters of TFBSs (HCTs) type and hosts a CpG Island (CGI). To test the functional correlation further with GLA, RPL36A, and HNRNPH2, we used siRNAs to knock down GH0XJ101390 in human kidney embryonic cells 293T. The results showed a significant decrease in RPL36A and GLA expression and a non-significant decrease in HNRNPH2 expression. These findings could have important implications for understanding the regulatory mechanisms of GH0XJ101390 and its potential role in FD. A better understanding of these mechanisms may improve diagnostic and therapeutic methods for FD, which could ultimately benefit patients with this rare condition.

iCAZyGFADB: an insect CAZyme and gene function annotation database

With the continuous upgrading of high-throughput sequencing technology, a large amount of biological genome data has been deciphered and published. The research on functional genes of biological genomes urgently needs a collection of service websites with user-friendly and full annotation functions for a variety of gene function annotation tools. In this study, iCAZyGFADB, which is a database website integrating nine gene function annotation tools, was perfectly developed to meet the needs of biological genome functional annotation. Its nine gene function annotation tools were Carbohydrate-Active Enzymes (CAZyme) annotation, Gene Ontology (GO) annotation, Kyoto Encyclopedia of Genes and Genomes (KEGG) annotation, Cluster of Orthologous Gene (COG) annotation, Evolutionary Genealogy of Genes: Non-supervised Orthologous Groups (eggNOG) annotation, SwissProt annotation, Pfam annotation, KOG annotation and Animal Transcription Factor DataBase (AnimalTFDB) annotation. It has three advantages. First, it is superior to gene function annotation of other biological cloud analysis platforms and runs very fast. Second, all gene annotation functions of the website are free and open to users. Third, it can annotate eight gene functions (GO, KEGG, COG, eggNOG, SwissProt, Pfam, KOG and AnimalTFDB annotation) of a single species at the same time, while other cloud platforms do not have the ability or need to charge to open for users to complete the annotation of eight gene functions at the same time. Moreover, the development and operation of our database will provide great help for gene function annotation research and significantly improve the efficiency of genome function research and reduce the cost of bioinformatics analysis. Genomic functional annotation researchers can access this database through the following website: http://www.icazygfadb.org.cn/. Database URL:  http://www.icazygfadb.org.cn/.

Temporally specific gene expression and chromatin remodeling programs regulate a conserved Pdyn enhancer

Neuronal and behavioral adaptations to novel stimuli are regulated by temporally dynamic waves of transcriptional activity, which shape neuronal function and guide enduring plasticity. Neuronal activation promotes expression of an immediate early gene (IEG) program comprised primarily of activity-dependent transcription factors, which are thought to regulate a second set of late response genes (LRGs). However, while the mechanisms governing IEG activation have been well studied, the molecular interplay between IEGs and LRGs remain poorly characterized. Here, we used transcriptomic and chromatin accessibility profiling to define activity-driven responses in rat striatal neurons. As expected, neuronal depolarization generated robust changes in gene expression, with early changes (1 hr) enriched for inducible transcription factors and later changes (4 hr) enriched for neuropeptides, synaptic proteins, and ion channels. Remarkably, while depolarization did not induce chromatin remodeling after 1 hr, we found broad increases in chromatin accessibility at thousands of sites in the genome at 4 hr after neuronal stimulation. These putative regulatory elements were found almost exclusively at non-coding regions of the genome, and harbored consensus motifs for numerous activity-dependent transcription factors such as AP-1. Furthermore, blocking protein synthesis prevented activity-dependent chromatin remodeling, suggesting that IEG proteins are required for this process. Targeted analysis of LRG loci identified a putative enhancer upstream of Pdyn (prodynorphin), a gene encoding an opioid neuropeptide implicated in motivated behavior and neuropsychiatric disease states. CRISPR-based functional assays demonstrated that this enhancer is both necessary and sufficient for Pdyn transcription. This regulatory element is also conserved at the human PDYN locus, where its activation is sufficient to drive PDYN transcription in human cells. These results suggest that IEGs participate in chromatin remodeling at enhancers and identify a conserved enhancer that may act as a therapeutic target for brain disorders involving dysregulation of Pdyn.

Deletions of singular U1 snRNA gene significantly interfere with transcription and 3'-end mRNA formation

Small nuclear RNAs (snRNAs) are structural and functional cores of the spliceosome. In metazoan genomes, each snRNA has multiple copies/variants, up to hundreds in mammals. However, the expressions and functions of each copy/variant in one organism have not been systematically studied. Focus on U1 snRNA genes, we investigated all five copies in Drosophila melanogaster using two series of constructed strains. Analyses of transgenic flies that each have a U1 promoter-driven gfp revealed that U1:21D is the major and ubiquitously expressed copy, and the other four copies have specificities in developmental stages and tissues. Mutant strains that each have a precisely deleted copy of U1-gene exhibited various extents of defects in fly morphology or mobility, especially deletion of U1:82Eb. Interestingly, splicing was changed at limited levels in the deletion strains, while large amounts of differentially-expressed genes and alternative polyadenylation events were identified, showing preferences in the down-regulation of genes with 1-2 introns and selection of proximal sites for 3'-end polyadenylation. In vitro assays suggested that Drosophila U1 variants pulled down fewer SmD2 proteins compared to the canonical U1. This study demonstrates that all five U1-genes in Drosophila have physiological functions in development and play regulatory roles in transcription and 3'-end formation.

Transcriptomic landscape of posterior regeneration in the annelid Platynereis dumerilii

BACKGROUND

Restorative regeneration, the capacity to reform a lost body part following amputation or injury, is an important and still poorly understood process in animals. Annelids, or segmented worms, show amazing regenerative capabilities, and as such are a crucial group to investigate. Elucidating the molecular mechanisms that underpin regeneration in this major group remains a key goal. Among annelids, the nereididae Platynereis dumerilii (re)emerged recently as a front-line regeneration model. Following amputation of its posterior part, Platynereis worms can regenerate both differentiated tissues of their terminal part as well as a growth zone that contains putative stem cells. While this regeneration process follows specific and reproducible stages that have been well characterized, the transcriptomic landscape of these stages remains to be uncovered.

RESULTS

We generated a high-quality de novo Reference transcriptome for the annelid Platynereis dumerilii. We produced and analyzed three RNA-sequencing datasets, encompassing five stages of posterior regeneration, along with blastema stages and non-amputated tissues as controls. We included two of these regeneration RNA-seq datasets, as well as embryonic and tissue-specific datasets from the literature to produce a Reference transcriptome. We used this Reference transcriptome to perform in depth analyzes of RNA-seq data during the course of regeneration to reveal the important dynamics of the gene expression, process with thousands of genes differentially expressed between stages, as well as unique and specific gene expression at each regeneration stage. The study of these genes highlighted the importance of the nervous system at both early and late stages of regeneration, as well as the enrichment of RNA-binding proteins (RBPs) during almost the entire regeneration process.

CONCLUSIONS

In this study, we provided a high-quality de novo Reference transcriptome for the annelid Platynereis that is useful for investigating various developmental processes, including regeneration. Our extensive stage-specific transcriptional analysis during the course of posterior regeneration sheds light upon major molecular mechanisms and pathways, and will foster many specific studies in the future.

Integrated Multi-Tissue Transcriptome Profiling Characterizes the Genetic Basis and Biomarkers Affecting Reproduction in Sheep (Ovis aries)

The heritability of litter size in sheep is low and controlled by multiple genes, but the research on its related genes is not sufficient. Here, to explore the expression pattern of multi-tissue genes in Chinese native sheep, we selected 10 tissues of the three adult ewes with the highest estimated breeding value in the early study of the prolific Xinggao sheep population. The global gene expression analysis showed that the ovary, uterus, and hypothalamus expressed the most genes. Using the Uniform Manifold Approximation and Projection (UMAP) cluster analysis, these samples were clustered into eight clusters. The functional enrichment analysis showed that the genes expressed in the spleen, uterus, and ovary were significantly enriched in the Ataxia Telangiectasia Mutated Protein (ATM) signaling pathway, and most genes in the liver, spleen, and ovary were enriched in the immune response pathway. Moreover, we focus on the expression genes of the hypothalamic-pituitary-ovarian axis (HPO) and found that 11,016 genes were co-expressed in the three tissues, and different tissues have different functions, but the oxytocin signaling pathway was widely enriched. To further explore the differences in the expression genes (DEGs) of HPO in different sheep breeds, we downloaded the transcriptome data in the public data, and the analysis of DEGs (Xinggao sheep vs. Sunite sheep in Hypothalamus, Xinggao sheep vs. Sunite sheep in Pituitary, and Xinggao sheep vs. Suffolk sheep in Ovary) revealed the neuroactive ligand-receptor interactions. In addition, the gene subsets of the transcription factors (TFs) of DEGs were identified. The results suggest that 51 TF genes and the homeobox TF may play an important role in transcriptional variation across the HPO. Altogether, our study provided the first fundamental resource to investigate the physiological functions and regulation mechanisms in sheep. This important data contributes to improving our understanding of the reproductive biology of sheep and isolating effecting molecular markers that can be used for genetic selection in sheep.

Identifying Novel Drug Targets for Epilepsy Through a Brain Transcriptome-Wide Association Study and Protein-Wide Association Study with Chemical-Gene-Interaction Analysis

Epilepsy is a severe neurological condition affecting 50-65 million individuals worldwide that can lead to brain damage. Nevertheless, the etiology of epilepsy remains poorly understood. Meta-analyses of genome-wide association studies involving 15,212 epilepsy cases and 29,677 controls of the ILAE Consortium cohort were used to conduct transcriptome-wide association studies (TWAS) and protein-wide association studies (PWAS). Furthermore, a protein-protein interaction (PPI) network was generated using the STRING database, and significant epilepsy-susceptible genes were verified using chip data. Chemical-related gene set enrichment analysis (CGSEA) was performed to determine novel drug targets for epilepsy. TWAS analysis identified 21,170 genes, of which 58 were significant (TWASfdr < 0.05) in ten brain regions, and 16 differentially expressed genes were verified based on mRNA expression profiles. The PWAS identified 2249 genes, of which 2 were significant (PWASfdr < 0.05). Through chemical-gene set enrichment analysis, 287 environmental chemicals associated with epilepsy were identified. We identified five significant genes (WIPF1, IQSEC1, JAM2, ICAM3, and ZNF143) that had causal relationships with epilepsy. CGSEA identified 159 chemicals that were significantly correlated with epilepsy (Pcgsea < 0.05), such as pentobarbital, ketone bodies, and polychlorinated biphenyl. In summary, we performed TWAS, PWAS (for genetic factors), and CGSEA (for environmental factors) analyses and identified several epilepsy-associated genes and chemicals. The results of this study will contribute to our understanding of genetic and environmental factors for epilepsy and may predict novel drug targets.

Temporally specific gene expression and chromatin remodeling programs regulate a conserved Pdyn enhancer

Neuronal and behavioral adaptations to novel stimuli are regulated by temporally dynamic waves of transcriptional activity, which shape neuronal function and guide enduring plasticity. Neuronal activation promotes expression of an immediate early gene (IEG) program comprised primarily of activity-dependent transcription factors, which are thought to regulate a second set of late response genes (LRGs). However, while the mechanisms governing IEG activation have been well studied, the molecular interplay between IEGs and LRGs remain poorly characterized. Here, we used transcriptomic and chromatin accessibility profiling to define activity-driven responses in rat striatal neurons. As expected, neuronal depolarization generated robust changes in gene expression, with early changes (1 h) enriched for inducible transcription factors and later changes (4 h) enriched for neuropeptides, synaptic proteins, and ion channels. Remarkably, while depolarization did not induce chromatin remodeling after 1 h, we found broad increases in chromatin accessibility at thousands of sites in the genome at 4 h after neuronal stimulation. These putative regulatory elements were found almost exclusively at non-coding regions of the genome, and harbored consensus motifs for numerous activity-dependent transcription factors such as AP-1. Furthermore, blocking protein synthesis prevented activity-dependent chromatin remodeling, suggesting that IEG proteins are required for this process. Targeted analysis of LRG loci identified a putative enhancer upstream of Pdyn (prodynorphin), a gene encoding an opioid neuropeptide implicated in motivated behavior and neuropsychiatric disease states. CRISPR-based functional assays demonstrated that this enhancer is both necessary and sufficient for Pdyn transcription. This regulatory element is also conserved at the human PDYN locus, where its activation is sufficient to drive PDYN transcription in human cells. These results suggest that IEGs participate in chromatin remodeling at enhancers and identify a conserved enhancer that may act as a therapeutic target for brain disorders involving dysregulation of Pdyn.

Dynamic alterations in metabolomics and transcriptomics associated with intestinal fibrosis in a 2,4,6-trinitrobenzene sulfonic acid-induced murine model

BACKGROUND & AIMS

Intestinal fibrosis is a common and severe complication of inflammatory bowel disease without clear pathogenesis. Abnormal expression of host genes and metabolic perturbations might associate with the onset of intestinal fibrosis. In this study, we aimed to investigate the relationship between the development of intestinal fibrosis and the dynamic alterations in both fecal metabolites and host gene expression.

METHODS

We induced intestinal fibrosis in a murine model using 2,4,6-trinitrobenzene sulfonic acid (TNBS). TNBS-treated or control mice were sacrificed after 4 and 6 weeks of intervention; alterations in colonic genes and fecal metabolites were determined by transcriptomics and metabolomics, respectively. Differential, tendency, enrichment, and correlation analyses were performed to assess the relationship between host genes and fecal metabolites.

RESULTS

RNA-sequencing analysis revealed that 679 differential genes with enduring changes were mainly enriched in immune response-related signaling pathways and metabolism-related biological processes. Among them, 15 lipid metabolism-related genes were closely related to the development of intestinal fibrosis. Moreover, the fecal metabolic profile was significantly altered during intestinal fibrosis development, especially the lipid metabolites. Particularly, dynamic perturbations in lipids were strongly associated with alterations in lipid metabolism-related genes expression. Additionally, six dynamically altered metabolites might serve as biomarkers to identify colitis-related intestinal fibrosis in the murine model.

CONCLUSIONS

Intestinal fibrosis in colitis mice might be related to dynamic changes in gene expression and metabolites. These findings could provide new insights into the pathogenesis of intestinal fibrosis.

SMRT Sequencing Technology Was Used to Construct the Batocera horsfieldi (Hope) Transcriptome and Reveal Its Features

Batocera horsfieldi (Hope) (Coleoptera: Cerambycidae) is an important forest pest in China that mainly infests timber and economic forests. This pest primarily causes plant tissue to necrotize, rot, and eventually die by feeding on the woody parts of tree trunks. To gain a deeper understanding of the genetic mechanism of B. horsfieldi, this study employed single-molecule real-time sequencing (SMRT) and Illumina RNA-seq technologies to conduct full-length transcriptome sequencing of the insect. Total RNA extracted from male and female adults was mixed and subjected to SMRT sequencing, generating a complete transcriptome. Transcriptome analysis, prediction of long non-coding RNA (lncRNA), coding sequences (CDs), analysis of simple sequence repeats (SSR), prediction of transcription factors, and functional annotation of transcripts were performed in this study. The collective 20,356,793 subreads (38.26 G, clean reads) were generated, including 432,091 circular consensus sequences and 395,851 full-length non-chimera reads. The full-length non-chimera reads (FLNC) were clustered and redundancies were removed, resulting in 39,912 consensus reads. SSR and ANGEL software v3.0 were used for predicting SSR and CDs. In addition, four tools were used for annotating 6058 lncRNAs, identifying 636 transcription factors. Furthermore, a total of 84,650 transcripts were functionally annotated in seven different databases. This is the first time that the full-length transcriptome of B. horsfieldi has been obtained using SMRT sequencing. This provides an important foundation for investigating the gene regulation underlying the interaction between B. horsfieldi and its host plants through gene editing in the future and provides a scientific basis for the prevention and control of B. horsfieldi.