Gene length as a biological timer to establish temporal transcriptional regulation

Genome-wide identification and expression analysis of epigenetic regulator gene families in the medicinal mushroom Ganoderma lucidum

Epigenetic regulator (ER) genes, crucial for fungal growth and development, remain largely unexplored in Ganoderma lucidum, a medicinal mushroom valued for its bioactive compounds. This study identified 81 ER genes in G. lucidum, distributed across 12 chromosomes and classified into six families: 3 chromatin remodelers, 4 DNA methyltransferases, 7 histone acetyltransferases, 22 histone deacetylases, 23 histone methyltransferases, and 22 histone demethyltransferases. Comparative and phylogenetic analyses with other species revealed conserved orthologs and species-specific clusters. Gene duplication analysis suggested whole-genome duplication expanded ER gene families, primarily histone demethyltransferases under purifying selection. Additionally, gene structure, motif, and domain analyses revealed family-specific intron/exon organization and conserved domains. Transcriptome profiling across four developmental stages (mycelium, primordia, young and mature fruiting body) revealed dynamic stage-specific expression patterns, suggesting their developmental significance. The result of qRT-PCR validated the expression patterns for 18 ER genes, laying foundation for future research exploring epigenetic regulation in fungal development and bioactive compound production.

Temporal characterisation and electrophysiological implications of TBI-induced serine/threonine kinase activity in mouse cortex

Traumatic brain injury (TBI) remains the leading cause of death and disability worldwide with no existing effective treatment. The early phase after TBI induction triggers numerous molecular cascades to regulate adaptive processes and cortical network activity. Kinases play a particularly prominent role in modifying peptide substrates, which include ion channels, receptors, transcription factors and inflammatory mediators. This study aimed to better understand the post-injury serine/threonine kinome; (1) Which kinases conduct phosphorylation-induced alterations of target peptides following unilateral TBI in mouse cortex? (2) How do these kinases effectuate pathological network hyperexcitability, which has detrimental long-term outcomes? We used a serine/threonine kinase assay at 4 h, 24 h and 72 h post-TBI to identify hyper-/hypo-active/phosphorylated kinases and peptides in the ipsilateral and contralateral cortical hemispheres relative to sham-operated controls. We pharmacologically mimicked the changes seen in ERK1/2 and PKC kinase activity, and using microelectrode array recordings we explored their significant electrophysiological implications on spontaneous and evoked cortical activity. We then used these findings to manipulate key kinase activity changes at 24 h post-TBI to rescue the hyperexcitability that is seen in the contralateral cortical network at this timepoint back to sham level. The contribution of specific downstream peptide target channel/receptor subunits was also shown. We conclude that volatile kinase activity has potent implications on cortical network activity after the injury and that these kinases and/or their peptide substrates should be more seriously considered as therapeutic targets for the clinical treatment of TBI.

Comprehensive blueprint of Salmonella genomic plasticity identifies hotspots for pathogenicity genes

Understanding the dynamic evolution of Salmonella is vital for effective bacterial infection management. This study explores the role of the flexible genome, organised in regions of genomic plasticity (RGP), in shaping the pathogenicity of Salmonella lineages. Through comprehensive genomic analysis of 12,244 Salmonella spp. genomes covering 2 species, 6 subspecies, and 46 serovars, we uncover distinct integration patterns of pathogenicity-related gene clusters into RGP, challenging traditional views of gene distribution. These RGP exhibit distinct preferences for specific genomic spots, and the presence or absence of such spots across Salmonella lineages profoundly shapes strain pathogenicity. RGP preferences are guided by conserved flanking genes surrounding integration spots, implicating their involvement in regulatory networks and functional synergies with integrated gene clusters. Additionally, we emphasise the multifaceted contributions of plasmids and prophages to the pathogenicity of diverse Salmonella lineages. Overall, this study provides a comprehensive blueprint of the pathogenicity potential of Salmonella. This unique insight identifies genomic spots in nonpathogenic lineages that hold the potential for harbouring pathogenicity genes, providing a foundation for predicting future adaptations and developing targeted strategies against emerging human pathogenic strains.

Non-canonical retinoid signaling in neural development, regeneration and synaptic function

Canonical retinoid signaling via nuclear receptors and gene regulation is critical for the initiation of developmental processes such as cellular differentiation, patterning and neurite outgrowth, but also mediates nerve regeneration and synaptic functions in adult nervous systems. In addition to canonical transcriptional regulation, retinoids also exert rapid effects, and there are now multiple lines of evidence supporting non-canonical retinoid actions outside of the nucleus, including in dendrites and axons. Together, canonical and non-canonical retinoid signaling provide the precise temporal and spatial control necessary to achieve the fine cellular coordination required for proper nervous system function. Here, we examine and discuss the evidence supporting non-canonical actions of retinoids in neural development and regeneration as well as synaptic function, including a review of the proposed molecular mechanisms involved.

Dissipative scaling of development and aging in multicellular organisms

Evolution, self-replication and ontogenesis are highly dynamic, irreversible and self-organizing processes dissipating energy. While progress has been made to decipher the role of thermodynamics in cellular fission, it is not yet clear how entropic balances shape organism growth and aging. This paper derives a general dissipation theory for the life history of organisms. It implies a self-regulated energy dissipation facilitating exponential growth within a hierarchical and entropy lowering self-organization. The theory predicts ceilings in energy expenditures imposed by geometric constrains, which promote thermal optimality during development, and a dissipative scaling across organisms consistent with ecological scaling laws combining isometric and allometric terms. The theory also illustrates how growing organisms can tolerate damage through continuous extension and production of new dissipative structures low in entropy. However, when organisms reduce their rate of cell division and reach a steady adult state, they become thermodynamically unstable, increase internal entropy by accumulating damage, and age.

Embryo-endometrial interaction associated with the location of the embryo during the mobility phase in mares

Embryo-maternal crosstalk is essential to establish pregnancy, with the equine embryo moving throughout the uterus on days 9-15 (ovulation = day 0) as part of this interaction. We hypothesized that the presence of a mobile embryo induces local changes in the gene expression of the endometrium. On Day 12, the endometrial transcripts were compared among three groups: uterine horn with an embryo (P+, n = 7), without an embryo (P-, n = 7) in pregnant mares, and both uterine horns of nonbred mares (NB, n = 6). We identified 1,101 differentially expressed genes (DEGs) between P+ vs. NB and 1,229 DEGs between P- vs. NB. The genes upregulated in both P+ and P- relative to NB were involved in growth factor pathway and fatty acid activation, while downregulated genes were associated with oxytocin signaling pathway and estrogen receptor signaling. Comparing the transcriptome of P+ to that of P-, we found 59 DEGs, of which 30 genes had a higher expression in P+. These genes are associated with regulating vascular growth factors and the immune system, all known to be essential in early pregnancy. Overall, this study suggests that the mobile embryo influences the endometrial gene expression locally.

Age or lifestyle-induced accumulation of genotoxicity is associated with a length-dependent decrease in gene expression

DNA damage has long been advocated as a molecular driver of aging. DNA damage occurs in a stochastic manner, and is therefore more likely to accumulate in longer genes. The length-dependent accumulation of transcription-blocking damage, unlike that of somatic mutations, should be reflected in gene expression datasets of aging. We analyzed gene expression as a function of gene length in several single-cell RNA sequencing datasets of mouse and human aging. We found a pervasive age-associated length-dependent underexpression of genes across species, tissues, and cell types. Furthermore, we observed length-dependent underexpression associated with UV-radiation and smoke exposure, and in progeroid diseases, Cockayne syndrome, and trichothiodystrophy. Finally, we studied published gene sets showing global age-related changes. Genes underexpressed with aging were significantly longer than overexpressed genes. These data highlight a previously undetected hallmark of aging and show that accumulation of genotoxicity in long genes could lead to reduced RNA polymerase II processivity.

Transcription and genome integrity

Transcription can cause genome instability by promoting R-loop formation but also act as a mutation-suppressing machinery by sensing of DNA lesions leading to the activation of DNA damage signaling and transcription-coupled repair. Recovery of RNA synthesis following the resolution of repair of transcription-blocking lesions is critical to avoid apoptosis and several new factors involved in this process have recently been identified. Some DNA repair proteins are recruited to initiating RNA polymerases and this may expediate the recruitment of other factors that participate in the repair of transcription-blocking DNA lesions. Recent studies have shown that transcription of protein-coding genes does not always give rise to spliced transcripts, opening the possibility that cells may use the transcription machinery in a splicing-uncoupled manner for other purposes including surveillance of the transcribed genome.

Codon Usage Bias Correlates With Gene Length in Neurodegeneration Associated Genes

Codon usage analysis is a crucial part of molecular characterization and is used to determine the factors affecting the evolution of a gene. The length of a gene is an important parameter that affects the characteristics of the gene, such as codon usage, compositional parameters, and sometimes, its functions. In the present study, we investigated the association of various parameters related to codon usage with the length of genes. Gene expression is affected by nucleotide disproportion. In sixty genes related to neurodegenerative disorders, the G nucleotide was the most abundant and the T nucleotide was the least. The nucleotide T exhibited a significant association with the length of the gene at both the overall compositional level and the first and second codon positions. Codon usage bias (CUB) of these genes was affected by pyrimidine and keto skews. Gene length was found to be significantly correlated with codon bias in neurodegeneration associated genes. In gene segments with lengths below 1,200 bp and above 2,400 bp, CUB was positively associated with length. Relative synonymous CUB, which is another measure of CUB, showed that codons TTA, GTT, GTC, TCA, GGT, and GGA exhibited a positive association with length, whereas codons GTA, AGC, CGT, CGA, and GGG showed a negative association. GC-ending codons were preferred over AT-ending codons. Overall analysis indicated that the association between CUB and length varies depending on the segment size; however, CUB of 1,200–2,000 bp gene segments appeared not affected by gene length. In synopsis, analysis suggests that length of the genes correlates with various imperative molecular signatures including A/T nucleotide disproportion and codon choices. In the present study we additionally evaluated various molecular features and their correlation with different indices of codon usage, like the Codon Adaptation Index (CAI) and Relative Dynonymous Codon Usage (RSCU) of codons. We also considered the impact of gene fragment size on different molecular features in genes related to neurodegeneration. This analysis will aid our understanding of and in potentially modulating gene expression in cases of defective gene functioning in clinical settings.

Does the Pachytene Checkpoint, a Feature of Meiosis, Filter Out Mistakes in Double-Strand DNA Break Repair and as a side-Effect Strongly Promote Adaptive Speciation?

This essay aims to explain two biological puzzles: why eukaryotic transcription units are composed of short segments of coding DNA interspersed with long stretches of non-coding (intron) DNA, and the near ubiquity of sexual reproduction. As is well known, alternative splicing of its coding sequences enables one transcription unit to produce multiple variants of each encoded protein. Additionally, padding transcription units with non-coding DNA (often many thousands of base pairs long) provides a readily evolvable way to set how soon in a cell cycle the various mRNAs will begin being expressed and the total amount of mRNA that each transcription unit can make during a cell cycle. This regulation complements control via the transcriptional promoter and facilitates the creation of complex eukaryotic cell types, tissues, and organisms. However, it also makes eukaryotes exceedingly vulnerable to double-strand DNA breaks, which end-joining break repair pathways can repair incorrectly. Transcription units cover such a large fraction of the genome that any mis-repair producing a reorganized chromosome has a high probability of destroying a gene. During meiosis, the synaptonemal complex aligns homologous chromosome pairs and the pachytene checkpoint detects, selectively arrests, and in many organisms actively destroys gamete-producing cells with chromosomes that cannot adequately synapse; this creates a filter favoring transmission to the next generation of chromosomes that retain the parental organization, while selectively culling those with interrupted transcription units. This same meiotic checkpoint, reacting to accidental chromosomal reorganizations inflicted by error-prone break repair, can, as a side effect, provide a mechanism for the formation of new species in sympatry. It has been a long-standing puzzle how something as seemingly maladaptive as hybrid sterility between such new species can arise. I suggest that this paradox is resolved by understanding the adaptive importance of the pachytene checkpoint, as outlined above.

KDM6A Regulates Cell Plasticity and Pancreatic Cancer Progression by Noncanonical Activin Pathway

BACKGROUND & AIMS

Inactivating mutations of KDM6A, a histone demethylase, were frequently found in pancreatic ductal adenocarcinoma (PDAC). We investigated the role of KDM6A (lysine demethylase 6A) in PDAC development.

METHODS

We performed a pancreatic tissue microarray analysis of KDM6A protein levels. We used human PDAC cell lines for KDM6A knockout and knockdown experiments. We performed bromouridine sequencing analysis to elucidate the effects of KDM6A loss on global transcription. We performed studies with Ptf1aCre; LSL-KrasG12D; Trp53R172H/+; Kdm6afl/fl or fl/Y, Ptf1aCre; Kdm6afl/fl or fl/Y, and orthotopic xenograft mice to investigate the impacts of Kdm6a deficiency on pancreatic tumorigenesis and pancreatitis.

RESULTS

Loss of KDM6A was associated with metastasis in PDAC patients. Bromouridine sequencing analysis showed up-regulation of the epithelial-mesenchymal transition pathway in PDAC cells deficient in KDM6A. Loss of KDM6A promoted mesenchymal morphology, migration, and invasion in PDAC cells in vitro. Mechanistically, activin A and subsequent p38 activation likely mediated the role of KDM6A loss. Inhibiting either activin A or p38 reversed the effect. Pancreas-specific Kdm6a-knockout mice pancreata showed accelerated PDAC progression, developed a more aggressive undifferentiated type of PDAC, and increased metastases in the background of Kras and p53 mutations. Kdm6a-deficient pancreata in a pancreatitis model had a delayed recovery with increased PDAC precursor lesions compared with wild-type pancreata.

CONCLUSIONS

Loss of KDM6A accelerates PDAC progression and metastasis, most likely by a noncanonical p38-dependent activin A pathway. KDM6A also promotes pancreatic tissue recovery from pancreatitis. Activin A might be used as a therapeutic target for KDM6A-deficient PDACs.

The hunt for a yet unknown: Common molecular signature in some genetically monomorphic enterobacteria

Mark Achtman introduced the term "genetically monomorphic bacteria" (GM bacteria) for some human and plant pathogens. They displayed a great uniformity in terms of their "genetic" properties. This "uniformity" poses a challenge to microbiologists. To address these problems, we used CodonW and IslandViewer 3 as analytical tools and took Escherichia coli, Salmonella, and Shigella strains as a model organisms. We hypothesized that GM bacterium contains a common molecular signature among them. We have found a significant correlation regarding the number of protein-coding genes, predicted highly expressed genes, and the highest length of gene in this regard. On the other hand, the correspondence analysis of pathogenicity-related genes identified by IslandViewer 3 displayed a somewhat unique pattern in GM bacteria. The probable pathogenic genes are clustered into two separate groups, which is a hallmark of some pattern. Similar genes of non-monomorphic pathogenic strain clustered almost similarly, but the clusters are joined together, they are not completely separated. These features, in our considered view, may be considered as codon usages signatures of these bacteria, and E. coli in particular.

Gene Size Matters: An Analysis of Gene Length in the Human Genome

While it is expected for gene length to be associated with factors such as intron number and evolutionary conservation, we are yet to understand the connections between gene length and function in the human genome. In this study, we show that, as expected, there is a strong positive correlation between gene length, transcript length, and protein size as well as a correlation with the number of genetic variants and introns. Among tissue-specific genes, we find that the longest transcripts tend to be expressed in the blood vessels, nerves, thyroid, cervix uteri, and the brain, while the smallest transcripts tend to be expressed in the pancreas, skin, stomach, vagina, and testis. We report, as shown previously, that natural selection suppresses changes for genes with longer transcripts and promotes changes for genes with smaller transcripts. We also observe that genes with longer transcripts tend to have a higher number of co-expressed genes and protein-protein interactions, as well as more associated publications. In the functional analysis, we show that bigger transcripts are often associated with neuronal development, while smaller transcripts tend to play roles in skin development and in the immune system. Furthermore, pathways related to cancer, neurons, and heart diseases tend to have genes with longer transcripts, with smaller transcripts being present in pathways related to immune responses and neurodegenerative diseases. Based on our results, we hypothesize that longer genes tend to be associated with functions that are important in the early development stages, while smaller genes tend to play a role in functions that are important throughout the whole life, like the immune system, which requires fast responses.

A report on DNA sequence determinants in gene expression

The biased usage of nucleotides in coding sequence and its correlation with gene expression has been observed in several studies. A complex set of interactions between genes and other components of the expression system determine the amount of proteins produced from coding sequences. It is known that the elongation rate of polypeptide chain is affected by both codon usage bias and specific amino acid compositional constraints. Therefore, it is of interest to review local DNA-sequence elements and other positional as well as combinatorial constraints that play significant role in gene expression.

Comparative and selection sweep analysis of CNV was associated to litter size in Dazu black goats

This study aims to identify the relative Copy number variation (CNV) associated with the litter size of Dazu black goats based on the unpublished CNV analytical results of our previously published sequencing data, in which the litter-size groups were classified into extreme low- and high-yield groups. Firstly, to compare the existence of valuable CNV in Dazu black goats with different fertility levels with mixed pools. We obtained 4992 and 4888 CNVs from the HY and LY, which overlapping 1461 genes, and classified on the original CNV type. Three genes [LOC108633278, PPP1R12A, and YIPF4] were observed in the intersection between the HY deletion and the LY duplication groups. Secondly, on individuals level, we identified a novel candidate CNV (Chr1_50215501, F_ST = 0.148, V_ST = 0.347) from 214 autosomal credible CNVs to be significant with litter size in the Dazu black goat, which located in the CBLB gene. This finding indicates the CBLB gene may affect the litter size of the Dazu black goats through structural variations, and Chr1_50215501 can be an effective genetic marker for marker-assisted selection breeding, and this study was also helps understand the molecular mechanism related to the goat litter size.

Characterization of novel primary miRNA transcription units in human cells using Bru-seq nascent RNA sequencing

MicroRNAs (miRNAs) are key contributors to gene regulatory networks. Because miRNAs are processed from RNA polymerase II transcripts, insight into miRNA regulation requires a comprehensive understanding of the regulation of primary miRNA transcripts. We used Bru-seq nascent RNA sequencing and hidden Markov model segmentation to map primary miRNA transcription units (TUs) across 32 human cell lines, allowing us to describe TUs encompassing 1443 miRNAs from miRBase and 438 from MirGeneDB. We identified TUs for 61 miRNAs with an unknown CAGE TSS signal for MirGeneDB miRNAs. Many primary transcripts containing miRNA sequences failed to generate mature miRNAs, suggesting that miRNA biosynthesis is under both transcriptional and post-transcriptional control. In addition to constitutive and cell-type specific TU expression regulated by differential promoter usage, miRNA synthesis can be regulated by transcription past polyadenylation sites (transcriptional read through) and promoter divergent transcription (PROMPTs). We identified 197 miRNA TUs with novel promoters, 97 with transcriptional read-throughs and 3 miRNA TUs that resemble PROMPTs in at least one cell line. The miRNA TU annotation data resource described here reveals a greater complexity in miRNA regulation than previously known and provides a framework for identifying cell-type specific differences in miRNA transcription in cancer and cell transition states.

Intronic CNVs and gene expression variation in human populations

Introns can be extraordinarily large and they account for the majority of the DNA sequence in human genes. However, little is known about their population patterns of structural variation and their functional implication. By combining the most extensive maps of CNVs in human populations, we have found that intronic losses are the most frequent copy number variants (CNVs) in protein-coding genes in human, with 12,986 intronic deletions, affecting 4,147 genes (including 1,154 essential genes and 1,638 disease-related genes). This intronic length variation results in dozens of genes showing extreme population variability in size, with 40 genes with 10 or more different sizes and up to 150 allelic sizes. Intronic losses are frequent in evolutionarily ancient genes that are highly conserved at the protein sequence level. This result contrasts with losses overlapping exons, which are observed less often than expected by chance and almost exclusively affect primate-specific genes. An integrated analysis of CNVs and RNA-seq data showed that intronic loss can be associated with significant differences in gene expression levels in the population (CNV-eQTLs). These intronic CNV-eQTLs regions are enriched for intronic enhancers and can be associated with expression differences of other genes showing long distance intron-promoter 3D interactions. Our data suggests that intronic structural variation of protein-coding genes makes an important contribution to the variability of gene expression and splicing in human populations.

Defective transcription elongation in a subset of cancers confers immunotherapy resistance

The nature and role of global transcriptional deregulations in cancers are not fully understood. We report that a large proportion of cancers have widespread defects in mRNA transcription elongation (TE). Cancers with TE defects (TEdeff) display spurious transcription and defective mRNA processing of genes characterized by long genomic length, poised promoters and inducible expression. Signaling pathways regulated by such genes, such as pro-inflammatory response pathways, are consistently suppressed in TEdeff tumors. Remarkably, TEdeff correlates with the poor response and outcome in immunotherapy, but not chemo- or targeted therapy, -treated renal cell carcinoma and metastatic melanoma patients. Forced pharmacologic or genetic induction of TEdeff in tumor cells impairs pro-inflammatory response signaling, and imposes resistance to the innate and adaptive anti-tumor immune responses and checkpoint inhibitor therapy in vivo. Therefore, defective TE is a previously unknown mechanism of tumor immune resistance, and should be assessed in cancer patients undergoing immunotherapy.

What does time mean in development?

Biology is dynamic. Timescales range from frenetic sub-second ion fluxes and enzymatic reactions to the glacial millions of years of evolutionary change. Falling somewhere in the middle of this range are the processes we usually study in development: cell division and differentiation, gene expression, cell-cell signalling, and morphogenesis. But what sets the tempo and manages the order of developmental events? Are the order and tempo different between species? How is the sequence of multiple events coordinated? Here, we discuss the importance of time for developing embryos, highlighting the necessity for global as well as cell-autonomous control. New reagents and tools in imaging and genomic engineering, combined with in vitro culture, are beginning to offer fresh perspectives and molecular insight into the origin and mechanisms of developmental time.

The role of H3K79 methylation in transcription and the DNA damage response

Chromatin plays a critical role in organizing and protecting DNA. However, chromatin acts as an impediment for transcription and DNA repair. Histone modifications, such as H3K79 methylation, promote transcription and genomic stability by enhancing transcription elongation and by serving as landing sites for proteins involved in the DNA damage response. This review summarizes the current understanding of the role of H3K79 methylation in transcription, how it affects genome stability and opportunities to develop impactful therapeutic interventions for cancer.