Search

Myocyte Enhancer Factor 2C (MEF2C) is a transcription factor that plays a crucial role in neurogenesis and synapse development. Genetic studies have identified MEF2C as a gene that influences cognition and risk for neuropsychiatric disorders, including autism spectrum disorder (ASD) and schizophrenia (SCZ). Here, we investigated the involvement of MEF2C in these phenotypes using human-derived neural stem cells (NSCs) and glutamatergic induced neurons (iNs), which represented early and late neurodevelopmental stages. For these cellular models, MEF2C function had previously been disrupted, either by direct or indirect mutation, and gene expression assayed using RNA-seq. We integrated these RNA-seq data with MEF2C ChIP-seq data to identify dysregulated direct target genes of MEF2C in the NSCs and iNs models. Several MEF2C direct target gene-sets were enriched for SNP-based heritability for intelligence, educational attainment and SCZ, as well as being enriched for genes containing rare de novo mutations reported in ASD and/or developmental disorders. These gene-sets are enriched in both excitatory and inhibitory neurons in the prenatal and adult brain and are involved in a wide range of biological processes including neuron generation, differentiation and development, as well as mitochondrial function and energy production. We observed a trans expression quantitative trait locus (eQTL) effect of a single SNP at MEF2C (rs6893807, which is associated with IQ) on the expression of a target gene, BNIP3L. BNIP3L is a prioritized risk gene from the largest genome-wide association study of SCZ and has a function in mitophagy in mitochondria. Overall, our analysis reveals that either direct or indirect disruption of MEF2C dysregulates sets of genes that contain multiple alleles associated with SCZ risk and cognitive function and implicates neuron development and mitochondrial function in the etiology of these phenotypes. Author summary: Schizophrenia is a complex disorder caused by many genes. Current drugs for schizophrenia are only partially effective and do not treat cognitive deficits, which are key factors for explaining disability, leading to unemployment, homelessness and social isolation. Large-scale genetic studies of schizophrenia and cognitive function have been effective at identifying individual SNPs and genes that contribute to these phenotypes but have struggled to immediately uncover the bigger picture of the underlying biology of the disorder. Here we take an individual gene associated with schizophrenia and cognitive function called MEF2C, which on its own is a very important regulator of brain development. We use functional genomics data from studies where MEF2C has been mutated to identify sets of other genes that are influenced by MEF2C in developing and mature neurons. We show that several of these gene-sets are enriched for common variants associated with schizophrenia and cognitive function, and for rare variants that increase risk of various neurodevelopmental disorders. These gene-sets are involved in neuron development and mitochondrial function, providing evidence that these biological processes may be important in the context of the molecular mechanisms that underpin schizophrenia and cognitive function. [ABSTRACT FROM AUTHOR]

Pseudouridine (Ψ) is an ubiquitous RNA modification, present in the tRNAs and rRNAs of species across all domains of life. Conserved pseudouridine synthases modify the mRNAs of diverse eukaryotes, but the modification has yet to be identified in bacterial mRNAs. Here, we report the discovery of pseudouridines in mRNA from E. coli. By testing the mRNA modification capacity of all 11 known pseudouridine synthases, we identify RluA as the predominant mRNA-modifying enzyme. RluA, a known tRNA and 23S rRNA pseudouridine synthase, modifies at least 31 of the 44 high-confidence sites we identified in E. coli mRNAs. Using RNA structure probing data to inform secondary structures, we show that the target sites of RluA occur in a common sequence and structural motif comprised of a ΨURAA sequence located in the loop of a short hairpin. This recognition element is shared with previously identified target sites of RluA in tRNAs and rRNA. Overall, our work identifies pseudouridine in key mRNAs and suggests the capacity of Ψ to regulate the transcripts that contain it. Author summary: While RNAs are initially transcribed using only the nucleotides A, G, C and U, these nucleotides can be enzymatically modified into many different post-transcriptional modifications. tRNAs and rRNAs across all domains of life are modified with pseudouridine, an isoform of uridine that is inserted by highly conserved enzymes. In many eukaryotes, it has been demonstrated that some of these enzymes can also insert pseudouridines in mRNA, where they are poised to impact gene expression through their effects on the transcript. Here we show that protein-coding transcripts in E. coli are also modified with pseudouridine, indicating that mRNA pseudouridylation is also a highly conserved activity. RluA is the main mRNA pseudouridine synthase, introducing the modification into the transcripts of dozens of protein coding genes with high specificity. Its target sites are defined by a combined sequence and secondary structure motif. Two additional enzymes, RluC and RluD, introduce a few additional sites. All three of these enzymes belong to the same sub-family of pseudouridine synthases, and homologs of these also modify mRNAs in humans. Thus, mRNA modification by these enzymes might be a conserved activity with the capacity to impact gene regulation. [ABSTRACT FROM AUTHOR]

Our understanding of the molecular pathways that regulate oogenesis and define cellular identity in the Arthropod female reproductive system and the extent of their conservation is currently very limited. This is due to the focus on model systems, including Drosophila and Daphnia, which do not reflect the observed diversity of morphologies, reproductive modes, and sex chromosome systems. We use single-nucleus RNA and ATAC sequencing to produce a comprehensive single nucleus atlas of the adult Artemia franciscana female reproductive system. We map our data to the Fly Cell Atlas single-nucleus dataset of the Drosophila melanogaster ovary, shedding light on the conserved regulatory programs between the two distantly related Arthropod species. We identify the major cell types known to be present in the Artemia ovary, including germ cells, follicle cells, and ovarian muscle cells. Additionally, we use the germ cells to explore gene regulation and expression of the Z chromosome during meiosis, highlighting its unique regulatory dynamics and allowing us to explore the presence of meiotic sex chromosome silencing in this group. Author summary: Oogenesis is a highly complex process involving multiple cell-types and an extremely orchestrated program that unfolds in the female reproductive system. Despite the extreme diversity of Arthropod reproductive modes and sex determination systems, our current understanding of oogenesis is limited to a few model species. This makes it difficult to study and formulate hypotheses about the evolutionary history, constraints, and importance of the individual elements of this process. To fill this gap, we used single-nucleus expression and chromatin-accessibility data to produce a single-nucleus atlas of the Artemia franciscana female reproductive system. By comparing our dataset to the published Drosophila single-nucleus data (over 400 million years of divergence), we were able to highlight the extreme conservation of several of the molecular pathways of oogenesis and meiosis. We found evidence of global transcriptional quiescence and chromatin condensation in late germ cells, highlighting the conserved role of this repressive stage in arthropod oogenesis. Additionally, we explored the expression patterns of the ZW sex chromosomes during oogenesis. Our data shows that the Z-chromosome is consistently downregulated in germline cells. While this is partly driven by a lack of dosage compensation in the germline, a subset of cells show stronger repression of the Z chromosome. [ABSTRACT FROM AUTHOR]

The pervasiveness of gene expression variation and its contribution to phenotypic variation and evolution is well known. This gene expression variation is context dependent, with differences in regulatory architecture often associated with intrinsic and environmental factors, and is modulated by regulatory elements that can act in cis (linked) or in trans (unlinked) relative to the genes they affect. So far, little is known about how this genetic variation affects the evolution of regulatory architecture among closely related tissues during population divergence. To address this question, we analyzed gene expression in the midgut, hindgut, and Malpighian tubule as well as microbiome composition in the two gut tissues in four Drosophila melanogaster strains and their F1 hybrids from two divergent populations: one from the derived, European range and one from the ancestral, African range. In both the transcriptome and microbiome data, we detected extensive tissue- and genetic background-specific effects, including effects of genetic background on overall tissue specificity. Tissue-specific effects were typically stronger than genetic background-specific effects, although the two gut tissues were not more similar to each other than to the Malpighian tubules. An examination of allele specific expression revealed that, while both cis and trans effects were more tissue-specific in genes expressed differentially between populations than genes with conserved expression, trans effects were more tissue-specific than cis effects. Despite there being highly variable regulatory architecture, this observation was robust across tissues and genetic backgrounds, suggesting that the expression of trans variation can be spatially fine-tuned as well as or better than cis variation during population divergence and yielding new insights into cis and trans regulatory evolution. Author summary: Genetic variants regulating gene expression can act in cis (linked) or in trans (unlinked) relative to the genes they affect and are thought to be important during adaptation because they can spatially and temporally fine-tune gene expression. In this study, we used the fruit fly Drosophila melanogaster to compare gene expression between inbred parental strains and their offspring in order to characterize the basis of gene expression regulation and inheritance. We examined gene expression in three tissues (midgut, hindgut, and Malpighian tubule) and four genetic backgrounds stemming from Europe and the ancestral range in Africa. Additionally, we characterized the bacterial community composition in the two gut tissues. We detected extensive tissue- and genetic background-specific effects on gene expression and bacterial community composition, although tissue-specific effects were typically stronger than genetic background effects. Genes with cis and trans regulatory effects were more tissue-specific than genes with conserved expression, while those with trans effects were more tissue-specific than those with cis effects. These results suggest that the expression of trans variation can be spatially fine-tuned as well as (or better than) cis variation as populations diverge from one another. Our study yields novel insight into the genetic basis of gene regulatory evolution. [ABSTRACT FROM AUTHOR]

Sex determination in the nematode C. elegans is controlled by the master regulator XOL-1 during embryogenesis. Expression of xol-1 is dependent on the ratio of X chromosomes and autosomes, which differs between XX hermaphrodites and XO males. In males, xol-1 is highly expressed and in hermaphrodites, xol-1 is expressed at very low levels. XOL-1 activity is known to be critical for the proper development of C. elegans males, but its low expression was considered to be of minimal importance in the development of hermaphrodite embryos. Our study reveals that XOL-1 plays an important role as a regulator of developmental timing during hermaphrodite embryogenesis. Using a combination of imaging and bioinformatics techniques, we found that hermaphrodite embryos have an accelerated rate of cell division, as well as a more developmentally advanced transcriptional program when xol-1 is lost. Further analyses reveal that XOL-1 is responsible for regulating the timing of initiation of dosage compensation on the X chromosomes, and the appropriate expression of sex-biased transcriptional programs in hermaphrodites. We found that xol-1 mutant embryos overexpress the H3K9 methyltransferase MET-2 and have an altered H3K9me landscape. Some of these effects of the loss of xol-1 gene were reversed by the loss of met-2. These findings demonstrate that XOL-1 plays an important role as a developmental regulator in embryos of both sexes, and that MET-2 acts as a downstream effector of XOL-1 activity in hermaphrodites. Author summary: Various organisms have differing ways of determining, at a molecular level, what the sex of a developing embryo is supposed to be. The two sexes in the nematode C. elegans, hermaphrodite and male, have different numbers of X chromosomes. Hermaphrodites have two X chromosomes and males only have one. This mismatch raises an additional problem as hermaphrodites will have twice the amount of genes expressed from the X compared to males. This is solved by a process called dosage compensation, which equalizes gene expression between the sexes. A molecular sensor called XOL-1 detects the number of X chromosomes in an embryo and kick-starts the process of proper sexual development and/or dosage compensation. XOL-1 was known to be very important for activating male development but was believed to not have any role in the development of the hermaphrodite sex. We show that XOL-1 has some crucial roles in hermaphrodite embryos in controlling the rate of development of the embryo, and regulating the timing of dosage compensation. We also show that MET-2, a protein that deposits repressive methyl marks on DNA-bound histone proteins, is involved in this process. [ABSTRACT FROM AUTHOR]

Coordinated activation and inhibition of F-actin supports the movements of morphogenesis. Understanding the proteins that regulate F-actin is important, since these proteins are mis-regulated in diseases like cancer. Our studies of C. elegans embryonic epidermal morphogenesis identified the GTPase CED-10/Rac1 as an essential activator of F-actin. However, we need to identify the GEF, or Guanine-nucleotide Exchange Factor, that activates CED-10/Rac1 during embryonic cell migrations. The two-component GEF, CED-5/CED-12, is known to activate CED-10/Rac1 to promote cell movements that result in the engulfment of dying cells during embryogenesis, and a later cell migration of the larval Distal Tip Cell. It is believed that CED-5/CED-12 powers cellular movements of corpse engulfment and DTC migration by promoting F-actin formation. Therefore, we tested if CED-5/CED-12 was involved in embryonic migrations, and got a contradictory result. CED-5/CED-12 definitely support embryonic migrations, since their loss led to embryos that died due to failed epidermal cell migrations. However, CED-5/CED-12 inhibited F-actin in the migrating epidermis, the opposite of what was expected for a CED-10 GEF. To address how CED-12/CED-5 could have two opposing effects on F-actin, during corpse engulfment and cell migration, we investigated if CED-12 harbors GAP (GTPase Activating Protein) functions. A candidate GAP region in CED-12 faces away from the CED-5 GEF catalytic region. Mutating a candidate catalytic Arginine in the CED-12 GAP region (R537A) altered the epidermal cell migration function, and not the corpse engulfment function. We interfered with GEF function by interfering with CED-5's ability to bind Rac1/CED-10. Mutating Serine-Arginine in CED-5/DOCK predicted to bind and stabilize Rac1 for catalysis, resulted in loss of both ventral enclosure and corpse engulfment. Genetic and expression studies strongly support that the GAP function likely acts on different GTPases. Thus, we propose CED-5/CED-12 support the cycling of multiple GTPases, by using distinct domains, to both promote and inhibit F-actin nucleation. Author summary: GTPases in their active state promote actin nucleation that drives cellular events, from cell migrations, to cell shape changes, to cell-cell interactions. To function correctly, GTPases need to cycle from the active, GTP-bound state, to the inactive, GDP-bound state. This cycle is supported by Guanine-nucleotide Exchange Factors, or GEFs, that support activation as GDP is switched for GTP, and GTPase-Activating Proteins, or GAPs that support hydrolysis back to the GDP bound state. The Rac1/CED-10 GTPase has a well-studied GEF, CED-5/CED-12, that promotes Rac1 activation during cell engulfment of dying cells. Here we tested if CED-5/CED-12 also functioned as the activator, or GEF for Rac1 during embryonic epidermal cell migrations. Surprisingly, CED-5/CED-12 behaved completely opposite to what was expected during epidermal cell migration. Therefore, we investigated if CED-5/CED-12 could harbor a GAP function. Comparing models of human and C. elegans protein structures suggested a putative GAP region, which we mutated to show that CED-12 likely functions as a GAP. Genetic and gene expression tests identify other GTPases, CDC-42 and RHO-1, as likely targets of this newly uncovered CED-12 GAP function. These findings suggest CED-5/CED-12 can support the cycling of multiple GTPases, and either promote or inhibit F-actin nucleation, using its GEF or GAP functions. [ABSTRACT FROM AUTHOR]

Deciphering the evolutionary forces controlling insecticide resistance in malaria vectors remains a prerequisite to designing molecular tools to detect and assess resistance impact on control tools. Here, we demonstrate that a 4.3kb transposon-containing structural variation is associated with pyrethroid resistance in central/eastern African populations of the malaria vector Anopheles funestus. In this study, we analysed Pooled template sequencing data and direct sequencing to identify an insertion of 4.3kb containing a putative retro-transposon in the intergenic region of two P450s CYP6P5-CYP6P9b in mosquitoes of the malaria vector Anopheles funestus from Uganda. We then designed a PCR assay to track its spread temporally and regionally and decipher its role in insecticide resistance. The insertion originates in or near Uganda in East Africa, where it is fixed and has spread to high frequencies in the Central African nation of Cameroon but is still at low frequency in West Africa and absent in Southern Africa. A marked and rapid selection was observed with the 4.3kb-SV frequency increasing from 3% in 2014 to 98% in 2021 in Cameroon. A strong association was established between this SV and pyrethroid resistance in field populations and is reducing pyrethroid-only nets' efficacy. Genetic crosses and qRT-PCR revealed that this SV enhances the expression of CYP6P9a/b but not CYP6P5. Within this structural variant (SV), we identified putative binding sites for transcription factors associated with the regulation of detoxification genes. An inverse correlation was observed between the 4.3kb SV and malaria parasite infection, indicating that mosquitoes lacking the 4.3kb SV were more frequently infected compared to those possessing it. Our findings highlight the underexplored role and rapid spread of SVs in the evolution of insecticide resistance and provide additional tools for molecular surveillance of insecticide resistance. Author summary: Anopheles funestus is an important malaria transmitting mosquito whose control relies on pyrethroid insecticides. The ability to control malaria transmitting mosquitoes is threatened by the emergence of insecticide resistance. Developing effective strategies to prevent, slow, or overcome resistance require an understanding of the underlying genetic mechanisms. In this study, we identified a transposon containing insertion in Anopheles funestus mosquitoes in Central/East Africa but absent in other regions. Tracking its distribution and frequency over the years revealed a strong selection with frequencies reaching fixation in less than 3 years. Phenotype/genotype association studies revealed a strong correlation between pyrethroid resistance and genetic variant genotype but negative correlation with plasmodium parasite infection in mosquitoes. These findings shed light on the often-overlooked role of structural variations in the evolution of insecticide resistance, offering valuable insights and tools for monitoring and addressing this critical issue in malaria control. [ABSTRACT FROM AUTHOR]

Bulky DNA adducts such as those induced by ultraviolet light are removed from the genomes of multicellular organisms by nucleotide excision repair, which occurs through two distinct mechanisms, global repair, requiring the DNA damage recognition-factor XPC (xeroderma pigmentosum complementation group C), and transcription-coupled repair (TCR), which does not. TCR is initiated when elongating RNA polymerase II encounters DNA damage, and thus analysis of genome-wide excision repair in XPC-mutants only repairing by TCR provides a unique opportunity to map transcription events missed by methods dependent on capturing RNA transcription products and thus limited by their stability and/or modifications (5'-capping or 3'-polyadenylation). Here, we have performed eXcision Repair-sequencing (XR-seq) in the model organism Caenorhabditis elegans to generate genome-wide repair maps in a wild-type strain with normal excision repair, a strain lacking TCR (csb-1), and a strain that only repairs by TCR (xpc-1). Analysis of the intersections between the xpc-1 XR-seq repair maps with RNA-mapping datasets (RNA-seq, long- and short-capped RNA-seq) reveal previously unrecognized sites of transcription and further enhance our understanding of the genome of thtableis important model organism. Author summary: Most RNA sequencing technologies are based on capturing either processed RNA or primary transcripts to define the RNA expression profile of cell lines or organisms. These methods suffer from RNA's inherent instability in revealing a complete view of RNA Polymerase II (RNAPII) transcription. Here, we used a C. elegans xpc-1 mutant which only performs transcription-coupled repair (TCR) to generate a complete RNAPII transcription profile by using the eXcision Repair-sequencing (XR-seq) technology to analyze repair of UV-induced DNA damage. XR-seq captures all the excision products from RNAPII-transcribed regions in the form of 21–28 nucleotide-long DNA oligomers and thus in a sense converts the unstable RNA signal to a more stable DNA signal enabling the generation of the most complete RNAPII transcript profile to date including heretofore unknown transcripts. [ABSTRACT FROM AUTHOR]

The piRNA pathway is a conserved germline-specific small RNA pathway that ensures genomic integrity and continued fertility. In C. elegans and other nematodes, Type-I piRNAs are expressed from >10,000 independently transcribed genes clustered within two discrete domains of 1.5 and 3.5 MB on Chromosome IV. Clustering of piRNA genes contributes to their germline-specific expression, but the underlying mechanisms are unclear. We analyze isolated germ nuclei to demonstrate that the piRNA genomic domains are located in a heterochromatin-like environment. USTC (Upstream Sequence Transcription Complex) promotes strong association of nucleosomes throughout piRNA clusters, yet organizes the local nucleosome environment to direct the exposure of individual piRNA genes. Localization of USTC to the piRNA domains depends upon the ATPase chromatin remodeler ISW-1, which maintains high nucleosome density across piRNA clusters and ongoing production of piRNA precursors. Overall, this work provides insight into how chromatin states coordinate transcriptional regulation over large genomic domains, with implications for global genome organization. Author summary: Here we address how gene expression is coordinated over long stretches of DNA that is embedded in a dense chromatin environment. We focus on a large DNA domain that contains thousands of tiny noncoding piRNA genes that are only expressed in germ cells. Our work demonstrates that a specific factor called USTC works with chromatin remodeling factors and with the core transcription complex to activate piRNA gene expression. The ISW-1 chromatin remodeler creates a dense environment that is favorable for USTC to be recruited to piRNA genes. USTC then clears the chromatin away from a space around each piRNA gene that allows for the factor TBP-1 to bind and presumably recruit the core transcriptional machinery to turn on the short piRNA. This work illuminates a mechanism by which both local and long-distance regulatory information can be integrated to coordinate gene expression. [ABSTRACT FROM AUTHOR]

Heterochromatin is critical for maintaining genome stability, especially in flowering plants, where it relies on a feedback loop involving the H3K9 methyltransferase, KRYPTONITE (KYP), and the DNA methyltransferase CHROMOMETHYLASE3 (CMT3). The H3K9 demethylase INCREASED IN BONSAI METHYLATION 1 (IBM1) counteracts the detrimental consequences of KYP-CMT3 activity in transcribed genes. IBM1 expression in Arabidopsis is uniquely regulated by methylation of the 7th intron, allowing it to monitor global H3K9me2 levels. We show the methylated intron is prevalent across flowering plants and its underlying sequence exhibits dynamic evolution. We also find extensive genetic and expression variations in KYP, CMT3, and IBM1 across flowering plants. We identify Arabidopsis accessions resembling weak ibm1 mutants and Brassicaceae species with reduced IBM1 expression or deletions. Evolution towards reduced IBM1 activity in some flowering plants could explain the frequent natural occurrence of diminished or lost CMT3 activity and loss of gene body DNA methylation, as cmt3 mutants in A. thaliana mitigate the deleterious effects of IBM1. Author summary: In flowering plants, IBM1 histone demethylase plays a crucial role in regulating chromatin structure by removing H3K9me2, a modification associated with heterochromatin. This process involves a distinctive mechanism where the methylation of a repetitive sequence within an intron allows IBM1 to monitor and respond to global H3K9me2 levels. We discovered that this intron methylation sensor is widely observed across flowering plants, although the exact sequences involved show considerable variation. This suggests that while the mechanism for controlling IBM1 activity is conserved, it has adapted differently in various plant species. Additionally, our findings include the identification of Arabidopsis thaliana accessions that mimic weak ibm1 mutants and several Brassicaceae species with diminished IBM1 expression or deletions, which correlate with reduced CMT3 activity and gene body DNA methylation. [ABSTRACT FROM AUTHOR]

Lager yeasts are limited to a few strains worldwide, imposing restrictions on flavour and aroma diversity and hindering our understanding of the complex evolutionary mechanisms during yeast domestication. The recent finding of diverse S. eubayanus lineages from Patagonia offers potential for generating new lager yeasts with different flavour profiles. Here, we leverage the natural genetic diversity of S. eubayanus and expand the lager yeast repertoire by including three distinct Patagonian S. eubayanus lineages. We used experimental evolution and selection on desirable traits to enhance the fermentation profiles of novel S. cerevisiae x S. eubayanus hybrids. Our analyses reveal an intricate interplay of pre-existing diversity, selection on species-specific mitochondria, de-novo mutations, and gene copy variations in sugar metabolism genes, resulting in high ethanol production and unique aroma profiles. Hybrids with S. eubayanus mitochondria exhibited greater evolutionary potential and superior fitness post-evolution, analogous to commercial lager hybrids. Using genome-wide screens of the parental subgenomes, we identified genetic changes in IRA2, IMA1, and MALX genes that influence maltose metabolism, and increase glycolytic flux and sugar consumption in the evolved hybrids. Functional validation and transcriptome analyses confirmed increased maltose-related gene expression, influencing greater maltotriose consumption in evolved hybrids. This study demonstrates the potential for generating industrially viable lager yeast hybrids from wild Patagonian strains. Our hybridization, evolution, and mitochondrial selection approach produced hybrids with high fermentation capacity and expands lager beer brewing options. Author summary: Lager beer dominates the global market, accounting for over 90% of commercial beer varieties. The main player in lager fermentation is the yeast Saccharomyces pastorianus, an interspecific hybrid between S. cerevisiae and S. eubayanus. Despite its popularity, the range of flavours and aromas found in lager beers is restricted by the low genetic diversity of available lager strains. Here, we explored if lager yeast profiles can be diversified by leveraging natural isolates of S. eubayanus from Chilean Patagonia. We generated de novo hybrids between S. cerevisiae and three distinct S. eubayanus Patagonian lineages. Through experimental evolution and selection on fermentation traits, we improved the fermentation profiles of the hybrids. We found that mutations in IRA2, IMA1, and MALX genes enhanced their maltose and maltotriose metabolism, resulting in higher ethanol production and unique aroma profiles. Our results also confirm that S. eubayanus mitochondria confer a greater evolutionary potential than S. cerevisiae mitochondria. The current study encourages the use of wild yeast strains to develop new brewing applications to expand the repertoire of de novo lager yeasts. [ABSTRACT FROM AUTHOR]

A quarter of humanity is estimated to have been exposed to Mycobacterium tuberculosis (Mtb) with a 5–10% risk of developing tuberculosis (TB) disease. Variability in responses to Mtb infection could be due to host or pathogen heterogeneity. Here, we focused on host genetic variation in a Peruvian population and its associations with gene regulation in monocyte-derived macrophages and dendritic cells (DCs). We recruited former household contacts of TB patients who previously progressed to TB (cases, n = 63) or did not progress to TB (controls, n = 63). Transcriptomic profiling of monocyte-derived DCs and macrophages measured the impact of genetic variants on gene expression by identifying expression quantitative trait loci (eQTL). We identified 330 and 257 eQTL genes in DCs and macrophages (False Discovery Rate (FDR) < 0.05), respectively. Four genes in DCs showed interaction between eQTL variants and TB progression status. The top eQTL interaction for a protein-coding gene was with FAH, the gene encoding fumarylacetoacetate hydrolase, which mediates the last step in mammalian tyrosine catabolism. FAH expression was associated with genetic regulatory variation in cases but not controls. Using public transcriptomic and epigenomic data of Mtb-infected monocyte-derived dendritic cells, we found that Mtb infection results in FAH downregulation and DNA methylation changes in the locus. Overall, this study demonstrates effects of genetic variation on gene expression levels that are dependent on history of infectious disease and highlights a candidate pathogenic mechanism through pathogen-response genes. Furthermore, our results point to tyrosine metabolism and related candidate TB progression pathways for further investigation. Author summary: Tuberculosis (TB), caused by Mycobacterium tuberculosis (Mtb), is a leading cause of mortality globally. Mtb-exposed individuals have heterogeneous outcomes following infection with Mtb, where some progress to TB, while many remain asymptomatic. We hypothesized that genetic variation could partly explain risk of TB, by regulating the expression of genes involved in the control of Mtb infection. In this study in Peru, we recruited former TB patients and Mtb-exposed household contacts of TB patients who remained asymptomatic to define how genetic variation regulates expression of genes in the innate immune system. We identified 4 genetic variants that have different effects on gene expression based on the TB history of the participant. The lead variant regulated expression of a key enzyme (Fumarylacetoacetate Hydrolase; FAH) in tyrosine catabolism. Knocking out the gene in myeloid cells increased susceptibility to Mtb infection. The results implicate FAH as a candidate host factor involved in TB progression. [ABSTRACT FROM AUTHOR]

Although introns are typically tens to thousands of nucleotides, there are notable exceptions. In flies as well as humans, a small number of genes contain introns that are more than 1000 times larger than typical introns, exceeding hundreds of kilobases (kb) to megabases (Mb). It remains unknown why gigantic introns exist and how cells overcome the challenges associated with their transcription and RNA processing. The Drosophila Y chromosome contains some of the largest genes identified to date: multiple genes exceed 4Mb, with introns accounting for over 99% of the gene span. Here we demonstrate that co-transcriptional splicing of these gigantic Y-linked genes is important to ensure successful transcription: perturbation of splicing led to the attenuation of transcription, leading to a failure to produce mature mRNA. Cytologically, defective splicing of the Y-linked gigantic genes resulted in disorganization of transcripts within the nucleus suggestive of entanglement of transcripts, likely resulting from unspliced long RNAs. We propose that co-transcriptional splicing maintains the length of nascent transcripts of gigantic genes under a critical threshold, preventing their entanglement and ensuring proper gene expression. Our study reveals a novel biological significance of co-transcriptional splicing. Author summary: A small number of genes in our genome are characterized by excessively large introns, sometimes exceeding megabases. It remains poorly understood why such large introns exist or how they are processed during gene expression. By examining genes on Drosophila Y chromosome that contain gigantic introns, it is shown that the transcripts of these genes are co-transcriptionally spliced, and co-transcriptional splicing of these genes is critical for proper gene expression. Perturbation of splicing led to attenuation of transcription specifically in genes with gigantic introns, but not genes with average-sized introns. Whereas transcripts of genes with gigantic introns normally form separate nuclear domains, splicing perturbation led to intermingling of these domains, leading to a proposal that co-transcriptional splicing may function to keep transcript length short to avoid entanglement of transcripts, which may stall the progression of RNA polymerase II. Taken together, these results reveal a novel function of co-transcriptional splicing that specifically impacts the expression of genes with gigantic introns. [ABSTRACT FROM AUTHOR]

Introduction: Glioblastoma (GBM) invasion studies have focused on coding genes, while few studies evaluate long non-coding RNAs (lncRNAs), transcripts without protein-coding potential, for role in GBM invasion. We leveraged CRISPR-interference (CRISPRi) to evaluate invasive function of GBM-associated lncRNAs in an unbiased functional screen, characterizing and exploring the mechanism of identified candidates. Methods: We implemented a CRISPRi lncRNA loss-of-function screen evaluating association of lncRNA knockdown (KD) with invasion capacity in Matrigel. Top screen candidates were validated using CRISPRi and oligonucleotide(ASO)-mediated knockdown in three tumor lines. Clinical relevance of candidates was assessed via The Cancer Genome Atlas(TCGA) and Genotype-Tissue Expression(GTEx) survival analysis. Mediators of lncRNA effect were identified via differential expression analysis following lncRNA KD and assessed for tumor invasion using knockdown and rescue experiments. Results: Forty-eight lncRNAs were significantly associated with 33–83% decrease in invasion (p<0.01) upon knockdown. The top candidate, LINC03045, identified from effect size and p-value, demonstrated 82.7% decrease in tumor cell invasion upon knockdown, while LINC03045 expression was significantly associated with patient survival and tumor grade(p<0.0001). RNAseq analysis of LINC03045 knockdown revealed that WASF3, previously implicated in tumor invasion studies, was highly correlated with lncRNA expression, while WASF3 KD was associated with significant decrease in invasion. Finally, WASF3 overexpression demonstrated rescue of invasive function lost with LINC03045 KD. Conclusion: CRISPRi screening identified LINC03045, a previously unannotated lncRNA, as critical to GBM invasion. Gene expression is significantly associated with tumor grade and survival. RNA-seq and mechanistic studies suggest that this novel lncRNA may regulate invasion via WASF3. Author summary: Glioblastoma is the most common malignant brain tumor to originate from brain tissue in adults. Most research has focused on the role of genes that code for protein, but genes that do not code for protein can still produce long non-coding RNAs (lncRNAs) that have been demonstrated to have important functions in cancer processes. We utilized a method called CRISPRi to screen for lncRNAs that are important in glioblastoma invasion, which is the capacity of cancer cells to invade into surrounding tissue. With this CRISPRi invasion screen, we identified LINC03045 as a gene that affects glioblastoma invasion. We validated this finding using cultured cancer stem cells and found that knocking down this gene decreased invasion. We also analyzed the public databases for clinical relevance, and found that LINC03045 was associated with patient survival and tumor grade. We also identified the gene WASF3, previously implicated in multiple tumor invasion studies (an element of the JAK-STAT pathway) as a downstream element of LINC03045. Thus, here, we have identified LINC03045 as a critical non-coding gene in glioblastoma invasion, and further characterizing this gene could potentially make it a therapeutic target for glioblastoma therapy in the future. [ABSTRACT FROM AUTHOR]

Long interspersed element 1 (LINE-1; L1) are a family of transposons that occupy ~17% of the human genome. Though a small number of L1 copies remain capable of autonomous transposition, the overwhelming majority of copies are degenerate and immobile. Nevertheless, both mobile and immobile L1s can exert pleiotropic effects (promoting genome instability, inflammation, or cellular senescence) on their hosts, and L1's contributions to aging and aging diseases is an area of active research. However, because of the cell type-specific nature of transposon control, the catalogue of L1 regulators remains incomplete. Here, we employ an eQTL approach leveraging transcriptomic and genomic data from the GEUVADIS and 1000Genomes projects to computationally identify new candidate regulators of L1 RNA levels in lymphoblastoid cell lines. To cement the role of candidate genes in L1 regulation, we experimentally modulate the levels of top candidates in vitro, including IL16, STARD5, HSD17B12, and RNF5, and assess changes in TE family expression by Gene Set Enrichment Analysis (GSEA). Remarkably, we observe subtle but widespread upregulation of TE family expression following IL16 and STARD5 overexpression. Moreover, a short-term 24-hour exposure to recombinant human IL16 was sufficient to transiently induce subtle, but widespread, upregulation of L1 subfamilies. Finally, we find that many L1 expression-associated genetic variants are co-associated with aging traits across genome-wide association study databases. Our results expand the catalogue of genes implicated in L1 RNA control and further suggest that L1-derived RNA contributes to aging processes. Given the ever-increasing availability of paired genomic and transcriptomic data, we anticipate this new approach to be a starting point for more comprehensive computational scans for regulators of transposon RNA levels. Author summary: Transposable elements, or jumping genes, are fragments of DNA that have or once had the ability to mobilize to a new location within our genome. In humans, the most abundant transposable element is LINE-1 (L1), accounting for ~17% of our total DNA. Though L1 is generally repressed in healthy human cells, depression of transposable elements (including L1) has been observed in aging and in aging-associated diseases. Additionally, there is increasing evidence that L1 transcriptional levels may promote features of aging, highlighting the importance of understanding the mechanisms that regulate L1 RNA levels. Here, we computationally identify new candidate regulators of L1 RNA levels, provide experimental evidence that candidate regulators influence L1 RNA levels, and demonstrate that genetic variants associated with differences in L1 RNA levels are co-associated with aging phenotypes. Our approach expands the toolkit that can be used to characterize transposable element regulation and highlights specific genes for further study. Importantly, our results reiterate the notion that L1 levels are linked with aging phenotypes and represent a potential therapeutic target for age-related decline. [ABSTRACT FROM AUTHOR]

Whether single-cell RNA-sequencing (scRNA-seq) captures the same biological information as single-nucleus RNA-sequencing (snRNA-seq) remains uncertain and likely to be context-dependent. Herein, a head-to-head comparison was performed in matched normal-adenocarcinoma human lung samples to assess biological insights derived from scRNA-seq versus snRNA-seq and better understand the cellular transition that occurs from normal to tumoral tissue. Here, the transcriptome of 160,621 cells/nuclei was obtained. In non-tumor lung, cell type proportions varied widely between scRNA-seq and snRNA-seq with a predominance of immune cells in the former (81.5%) and epithelial cells (69.9%) in the later. Similar results were observed in adenocarcinomas, in addition to an overall increase in cell type heterogeneity and a greater prevalence of copy number variants in cells of epithelial origin, which suggests malignant assignment. The cell type transition that occurs from normal lung tissue to adenocarcinoma was not always concordant whether cells or nuclei were examined. As expected, large differential expression of the whole-cell and nuclear transcriptome was observed, but cell-type specific changes of paired normal and tumor lung samples revealed a set of common genes in the cells and nuclei involved in cancer-related pathways. In addition, we showed that the ligand-receptor interactome landscape of lung adenocarcinoma was largely different whether cells or nuclei were evaluated. Immune cell depletion in fresh specimens partly mitigated the difference in cell type composition observed between cells and nuclei. However, the extra manipulations affected cell viability and amplified the transcriptional signatures associated with stress responses. In conclusion, research applications focussing on mapping the immune landscape of lung adenocarcinoma benefit from scRNA-seq in fresh samples, whereas snRNA-seq of frozen samples provide a low-cost alternative to profile more epithelial and cancer cells, and yield cell type proportions that more closely match tissue content. Author summary: Single-cell transcriptomic datasets provide unprecedented opportunities to disentangle the complex tissue microenvironment and cellular origin of cancer. Data are scarce regarding the pros and cons of single-cell RNA sequencing (scRNA-seq) of freshly explanted human tissues over single-nuclei sequencing (snRNA-seq) from the same archived frozen tissues. Lung adenocarcinoma represents a medically valuable case study to compare the biological signal recovered through cells and nuclei sequencing. Here, we sequenced the transcriptome of 160,621 cells/nuclei in paired normal-adenocarcinoma lung samples. Cell type proportions varied widely between scRNA-seq and snRNA-seq with a predominance of immune cells in the former and epithelial cells in the later. Adenocarcinomas were characterized by an increase in cell type heterogeneity and a greater prevalence of malignant epithelial cells in both scRNA-seq and snRNA-seq. The cellular and gene expression transition that occur from normal lung to adenocarcinoma showed common and discordant biological insights whether cells or nuclei were examined. Research applications focussing on mapping the immune landscape of lung cancer benefit from scRNA-seq in fresh samples, whereas snRNA-seq of the same frozen samples provide a low-cost and more flexible alternative to profile more epithelial and cancer cells, and yield cell type proportions that more closely match tissue content. [ABSTRACT FROM AUTHOR]

APOBEC-induced mutations occur in 50% of sequenced human tumors, with APOBEC3A (A3A) being a major contributor to mutagenesis in breast cancer cells. The mechanisms that cause A3A activation and mutagenesis in breast cancers are still unknown. Here, we describe factors that influence basal A3A mRNA transcript levels in breast cancer cells. We found that basal A3A mRNA correlates with A3A protein levels and predicts the amount of APOBEC signature mutations in a panel of breast cancer cell lines, indicating that increased basal transcription may be one mechanism leading to breast cancer mutagenesis. We also show that alteration of ERBB2 expression can drive A3A mRNA levels, suggesting the enrichment of the APOBEC mutation signature in Her2-enriched breast cancer could in part result from elevated A3A transcription. Hierarchical clustering of transcripts in primary breast cancers determined that A3A mRNA was co-expressed with other genes functioning in viral restriction and interferon responses. However, reduction of STAT signaling via inhibitors or shRNA in breast cancer cell lines had only minor impact on A3A abundance. Analysis of single cell RNA-seq from primary tumors indicated that A3A mRNA was highest in infiltrating immune cells within the tumor, indicating that correlations of A3A with STAT signaling in primary tumors may be result from higher immune infiltrates and are not reflective of STAT signaling controlling A3A expression in breast cancer cells. Analysis of ATAC-seq data in multiple breast cancer cell lines identified two transcription factor sites in the APOBEC3A promoter region that could promote A3A transcription. We determined that Rel-A, and Bach1, which have binding sites in these peaks, elevated basal A3A expression. Our findings highlight a complex and variable set of transcriptional activators for A3A in breast cancer cells. Author summary: The protein APOBEC3A (A3A) is a major cause of mutations in breast cancers by changing cytidines in single-stranded DNA to uridine. The amount of A3A transcript is elevated specifically in cancer cells containing A3A-induced mutations, suggesting that up-regulation of A3A transcription may be an underlying cause of cancer mutations. However, little is known about the factors that influence A3A transcript level in breast cancer. Here, we identify three different proteins that impact A3A transcription: Her2, Rel-A and Bach1. Her2 activation is a common oncogenic event in some breast cancer subtypes and likely alters transcription of a variety of genes, including APOBEC3A, through enhanced signaling pathways or cell stress signaling. Contrastingly, Rel-A and Bach1 are transcription factors that we provide evidence that likely impact A3A transcription by binding the A3A promoter. Additionally, we found that correlations between A3A expression and canonical A3A transcriptional regulators, STAT1 and STAT2, in primary tumors stems from immune cells infiltrating into the tumor instead of direct transcriptional regulation in breast cancer cells. Our results demonstrate that A3A transcriptional regulation is highly complex with multiple factors and inputs potentially leading to an increase in mutations. [ABSTRACT FROM AUTHOR]

Chromatin organization plays a crucial role in gene regulation by controlling the accessibility of DNA to transcription machinery. While significant progress has been made in understanding the regulatory role of clock proteins in circadian rhythms, how chromatin organization affects circadian rhythms remains poorly understood. Here, we employed ATAC-seq (Assay for Transposase-Accessible Chromatin with Sequencing) on FAC-sorted Drosophila clock neurons to assess genome-wide chromatin accessibility at dawn and dusk over the circadian cycle. We observed significant oscillations in chromatin accessibility at promoter and enhancer regions of hundreds of genes, with enhanced accessibility either at dusk or dawn, which correlated with their peak transcriptional activity. Notably, genes with enhanced accessibility at dusk were enriched with E-box motifs, while those more accessible at dawn were enriched with VRI/PDP1-box motifs, indicating that they are regulated by the core circadian feedback loops, PER/CLK and VRI/PDP1, respectively. Further, we observed a complete loss of chromatin accessibility rhythms in per01 null mutants, with chromatin consistently accessible at both dawn and dusk, underscoring the critical role of Period protein in driving chromatin compaction during the repression phase at dawn. Together, this study demonstrates the significant role of chromatin organization in circadian regulation, revealing how the interplay between clock proteins and chromatin structure orchestrates the precise timing of biological processes throughout the day. This work further implies that variations in chromatin accessibility might play a central role in the generation of diverse circadian gene expression patterns in clock neurons. Author summary: The intricate organization of chromatin is crucial for orchestrating precise gene expression patterns throughout various biological processes such as development and differentiation. In this study, we adapted ATAC-seq (Assay for Transposase-Accessible Chromatin) to assay chromatin accessibility genome-wide from FAC (Fluorescence-activated cell)-sorted Drosophila clock neurons at dawn and dusk timepoints over the circadian cycle. Our findings reveal that regulatory elements of hundreds of clock-regulated genes undergo diurnal oscillations in chromatin accessibility, with more accessible chromatin during the activation phase leading to higher transcription, and compacted chromatin during the repression phase leading to gene silencing. Motif analysis further elucidated distinct regulatory motifs associated with chromatin accessibility changes. Specifically, genes exhibiting increased accessibility around dusk were enriched with E-boxes, whereas those with enhanced accessibility around dawn showed enrichment for VP-boxes. Notably, our studies unveiled that these diurnal fluctuations in chromatin accessibility within clock neurons are dependent on the presence of a functional clock, as per01 null mutants displayed a loss of these rhythmic patterns. In summary, our studies shed light on the intricate relationship between circadian transcriptional regulation and chromatin accessibility dynamics. These insights underscore the role of chromatin accessibility dynamics in priming the genome for rhythmic gene expression, thereby contributing to the orchestration of circadian rhythms. [ABSTRACT FROM AUTHOR]

The position of the nucleus before it divides during mitosis is variable in different budding yeasts. Studies in the pathogenic intron-rich fungus Cryptococcus neoformans reveal that the nucleus moves entirely into the daughter bud before its division. Here, we report functions of a zinc finger motif containing spliceosome protein C. neoformans Slu7 (CnSlu7) in cell cycle progression. The budding yeast and fission yeast homologs of Slu7 have predominant roles for intron 3' splice site definition during pre-mRNA splicing. Using a conditional knockdown strategy, we show CnSlu7 is an essential factor for viability and is required for efficient cell cycle progression with major role during mitosis. Aberrant nuclear migration, including improper positioning of the nucleus as well as the spindle, were frequently observed in cells depleted of CnSlu7. However, cell cycle delays observed due to Slu7 depletion did not activate the Mad2-dependent spindle assembly checkpoint (SAC). Mining of the global transcriptome changes in the Slu7 knockdown strain identified downregulation of transcripts encoding several cell cycle regulators and cytoskeletal factors for nuclear migration, and the splicing of specific introns of these genes was CnSlu7 dependent. To test the importance of splicing activity of CnSlu7 on nuclear migration, we complemented Slu7 knockdown cells with an intron less PAC1 minigene and demonstrated that the nuclear migration defects were significantly rescued. These findings show that CnSlu7 regulates the functions of diverse cell cycle regulators and cytoskeletal components, ensuring timely cell cycle transitions and nuclear division during mitosis. Author summary: Nuclear position in eukaryotic cells is spatio-temporally regulated during the cell cycle. Unlike in Saccharomyces cerevisiae, where nuclear migration to mother-daughter bud neck precedes mitotic segregation, in Cryptococcus neoformans, the nucleus moves entirely into the daughter bud before division. Transcription dynamics during the C. neoformans cell cycle show periodic expression changes and since all nascent pre-mRNAs must be processed to mRNA, splicing is indispensable for gene expression. Yet, how evolutionarily conserved splicing factors modulate the complex intron-rich C. neoformans transcriptome and their impacts on atypical mitotic nuclear position is unexplored. Here, we show CnSlu7 a zinc finger motif spliceosome factor is essential for viability and for cell cycle progression with a major mitotic role. Live imaging of the nucleus, the spindle and other factors revealed significant abnormalities in nuclear migration from mother to daughter cells, as well as nuclear and spindle position defects upon CnSlu7 depletion. Depletion of CnSlu7 altered the transcript status of various cell cycle players, including those critical for nuclear migration, as their introns were CnSlu7 dependent for splicing. Together, these findings highlight the roles of C. neoformans Slu7 in fine-tuning gene expression levels of transcripts that ensures timely cell cycle progression and spatial control of nuclear position. [ABSTRACT FROM AUTHOR]

Regulation of transcription is a fundamental process that allows bacteria to respond to external stimuli with appropriate timing and magnitude of response. In the soil bacterium Bacillus subtilis, transcriptional regulation is at the core of developmental processes needed for cell survival. Gene expression in cells transitioning from exponential phase to stationary phase is under the control of a group of transcription factors called transition state regulators (TSRs). TSRs influence numerous developmental processes including the decision between biofilm formation and motility, genetic competence, and sporulation, but the extent to which TSRs influence bacterial physiology remains to be fully elucidated. Here, we demonstrate two TSRs, ScoC and AbrB, along with the MarR-family transcription factor PchR negatively regulate production of the iron chelator pulcherrimin in B. subtilis. Genetic analysis of the relationship between the three transcription factors indicate that all are necessary to limit pulcherrimin production during exponential phase and influence the rate and total amount of pulcherrimin produced. Similarly, expression of the pulcherrimin biosynthesis gene yvmC was found to be under control of ScoC, AbrB, and PchR and correlated with the amount of pulcherrimin produced by each background. Lastly, our in vitro data indicate a weak direct role for ScoC in controlling pulcherrimin production along with AbrB and PchR. The layered regulation by two distinct regulatory systems underscores the important, role for pulcherrimin in B. subtilis physiology. Author summary: Regulation of gene expression is important for survival in ever changing environments. In the soil bacterium Bacillus subtilis, key developmental processes are controlled by overlapping networks of transcription factors, some of which are termed transition state regulators (TSRs). Despite decades of research, the scope of how TSRs influence B. subtilis physiology is still unclear. We found that three transcription factors, two of which are TSRs, converge to inhibit production of the iron-chelator pulcherrimin. Only when all three are missing is pulcherrimin production elevated. Finally, we demonstrate that expression of pulcherrimin biosynthesis genes occurs via direct and indirect regulation by the trio of transcription factors. Due to its iron chelating ability, pulcherrimin has been characterized as a modulator of niche development with antioxidant properties. Thus, our findings that TSRs control pulcherrimin, concurrently with other developmental phenotypes, provides new insight into how TSRs impact B. subtilis and its interaction with the environment. [ABSTRACT FROM AUTHOR]

CgHog1, terminal kinase of the high-osmolarity glycerol signalling pathway, orchestrates cellular response to multiple external stimuli including surplus-environmental iron in the human fungal pathogen Candida glabrata (Cg). However, CgHog1 substrates remain unidentified. Here, we show that CgHog1 adversely affects Cg adherence to host stomach and kidney epithelial cells in vitro, but promotes Cg survival in the iron-rich gastrointestinal tract niche. Further, CgHog1 interactome and in vitro phosphorylation analysis revealed CgSub2 (putative RNA helicase) to be a CgHog1 substrate, with CgSub2 also governing iron homeostasis and host adhesion. CgSub2 positively regulated EPA1 (encodes a major adhesin) expression and host adherence via its interactor CgHtz1 (histone H2A variant). Notably, both CgHog1 and surplus environmental iron had a negative impact on CgSub2-CgHtz1 interaction, with CgHTZ1 or CgSUB2 deletion reversing the elevated adherence of Cghog1Δ to epithelial cells. Finally, the surplus-extracellular iron led to CgHog1 activation, increased CgSub2 phosphorylation, elevated CgSub2-CgHta (canonical histone H2A) interaction, and EPA1 transcriptional activation, thereby underscoring the iron-responsive, CgHog1-induced exchange of histone partners of CgSub2. Altogether, our work mechanistically defines how CgHog1 couples Epa1 adhesin expression with iron abundance, and point towards specific chromatin composition modification programs that probably aid fungal pathogens align their adherence to iron-rich (gut) and iron-poor (blood) host niches. Author summary: The human gut can serve as a source of disseminated infections caused by its resident microbes. Candida glabrata (Cg) is one such opportunistic fungal pathogen that is isolated from both iron-rich (human gut) and iron-poor (blood) host niches. Cg is presumed to translocate across the epithelial cell barrier of the gastrointestinal tract into the bloodstream. However, how Cg colonizes and traverses through the host gut, and adapts to the iron-rich gut milieu, is unknown. We show that Cg ingeniously attunes expression of the major cell surface protein Epa1, which mediates adherence to the host epithelial cells, to the environmental iron abundance, by exchanging the interacting partner (histone variant CgHtz1) of the CgSub2 RNA helicase. We also show that the conserved mitogen-activated protein kinase CgHog1, that may control this partner exchange by phosphorylating CgSub2, is essential for Cg survival in the mouse gut and bloodstream infections. Altogether, our work yields molecular insights into how nutrient fluctuations in the mammalian host lead to the chromatin-based transcriptional rewiring of a key virulence gene. [ABSTRACT FROM AUTHOR]

The Integrator is a multi-subunit protein complex that catalyzes the maturation of snRNA transcripts via 3' cleavage, a step required for snRNA incorporation with snRNP for spliceosome biogenesis. Here we developed a GFP based in vivo snRNA misprocessing reporter as a readout of Integrator function and performed a genome-wide RNAi screen for Integrator regulators. We found that loss of the Argonaute encoding csr-1 gene resulted in widespread 3' misprocessing of snRNA transcripts that is accompanied by a significant increase in alternative splicing. Loss of the csr-1 gene down-regulates the germline expression of Integrator subunits 4 and 6 and is accompanied by a reduced protein translation efficiency of multiple Integrator catalytic and non-catalytic subunits. Through isoform and motif mutant analysis, we determined that CSR-1's effect on snRNA processing is dependent on its catalytic slicer activity but does not involve the CSR-1a isoform. Moreover, mRNA-sequencing revealed high similarity in the transcriptome profile between csr-1 and Integrator subunit knockdown via RNAi. Together, our findings reveal CSR-1 as a new regulator of the Integrator complex and implicate a novel role of this Argonaute protein in snRNA 3' processing. Author summary: Small nuclear RNA molecules are an important structural component of the spliceosome, which is the cellular machinery that performs RNA splicing. Proper RNA splicing is important for ensuring proteins are accurately produced within the cell, and this process is broadly important for contributing to basic physiological processes such as development. Inside the cell, small nuclear RNA is processed by a protein complex called the Integrator, and mutations to the Integrator affect RNA splicing and can lead to neurodevelopmental defects. Here, we used the roundworm Caenorhabditis elegans to describe how an Argonaute encoding gene called csr-1 is required for maintaining Integrator protein expression and that loss of csr-1 gene expression contributes to the dysregulation of small nuclear RNA processing. These results provide new insights into our understanding of fundamental factors that regulate small nuclear RNA processing in cells, which are directly important in RNA splicing control. [ABSTRACT FROM AUTHOR]

Fuchs endothelial corneal dystrophy (FECD) is an age-related cause of vision loss, and the most common repeat expansion-mediated disease in humans characterised to date. Up to 80% of European FECD cases have been attributed to expansion of a non-coding CTG repeat element (termed CTG18.1) located within the ubiquitously expressed transcription factor encoding gene, TCF4. The non-coding nature of the repeat and the transcriptomic complexity of TCF4 have made it extremely challenging to experimentally decipher the molecular mechanisms underlying this disease. Here we comprehensively describe CTG18.1 expansion-driven molecular components of disease within primary patient-derived corneal endothelial cells (CECs), generated from a large cohort of individuals with CTG18.1-expanded (Exp+) and CTG 18.1-independent (Exp-) FECD. We employ long-read, short-read, and spatial transcriptomic techniques to interrogate expansion-specific transcriptomic biomarkers. Interrogation of long-read sequencing and alternative splicing analysis of short-read transcriptomic data together reveals the global extent of altered splicing occurring within Exp+ FECD, and unique transcripts associated with CTG18.1-expansions. Similarly, differential gene expression analysis highlights the total transcriptomic consequences of Exp+ FECD within CECs. Furthermore, differential exon usage, pathway enrichment and spatial transcriptomics reveal TCF4 isoform ratio skewing solely in Exp+ FECD with potential downstream functional consequences. Lastly, exome data from 134 Exp- FECD cases identified rare (minor allele frequency <0.005) and potentially deleterious (CADD>15) TCF4 variants in 7/134 FECD Exp- cases, suggesting that TCF4 variants independent of CTG18.1 may increase FECD risk. In summary, our study supports the hypothesis that at least two distinct pathogenic mechanisms, RNA toxicity and TCF4 isoform-specific dysregulation, both underpin the pathophysiology of FECD. We anticipate these data will inform and guide the development of translational interventions for this common triplet-repeat mediated disease. Author summary: Fuchs endothelial corneal dystrophy (FECD) leads to vision loss and is the most common repeat expansion-mediated disease characterised to date; most individuals with FECD harbour a non-coding CTG repeat expansion within the gene TCF4. FECD represents an important paradigm for other rare and devastating neurological repeat-mediated diseases, given its relatively mild and tissue-specific nature. Intriguingly, despite TCF4 being ubiquitously expressed, individuals with FECD only experience corneal disease, and the biological reason for this tissue-specificity remains elusive. Here, we use tissue from 31 individuals with FECD to perform complementary long-read, short-read and spatial transcriptomic analyses to enhance our understanding of mechanisms underpinning this disease. These data highlight that at least two mechanisms, RNA toxicity and TCF4 isoform dysregulation, underlie the disease state in affected corneal cells. Furthermore, TCF4 isoform skewing, with evidence of downregulation, suggests this mechanism in part may explain the unique vulnerability of the cornea. In addition, 7/134 FECD expansion negative cases were identified to harbour rare and potentially deleterious TCF4 variants, further supporting the hypothesis that dysregulation of TCF4 may be key to FECD pathophysiology. Biological insights presented here will guide the development of personalised FECD therapies and may inform the development of repeat-expansion mediated therapies more broadly. [ABSTRACT FROM AUTHOR]

Transcriptional cis-regulatory modules, e.g., enhancers, control the time and location of metazoan gene expression. While changes in enhancers can provide a powerful force for evolution, there is also significant deep conservation of enhancers for developmentally important genes, with function and sequence characteristics maintained over hundreds of millions of years of divergence. Not well understood, however, is how the overall regulatory composition of a locus evolves, with important outstanding questions such as how many enhancers are conserved vs. novel, and to what extent are the locations of conserved enhancers within a locus maintained? We begin here to address these questions with a comparison of the respective single-minded (sim) loci in the two dipteran species Drosophila melanogaster (fruit fly) and Aedes aegypti (mosquito). sim encodes a highly conserved transcription factor that mediates development of the arthropod embryonic ventral midline. We identify two enhancers in the A. aegypti sim locus and demonstrate that they function equivalently in both transgenic flies and transgenic mosquitoes. One A. aegypti enhancer is highly similar to known Drosophila counterparts in its activity, location, and autoregulatory capability. The other differs from any known Drosophila sim enhancers with a novel location, failure to autoregulate, and regulation of expression in a unique subset of midline cells. Our results suggest that the conserved pattern of sim expression in the two species is the result of both conserved and novel regulatory sequences. Further examination of this locus will help to illuminate how the overall regulatory landscape of a conserved developmental gene evolves. Author summary: The expression patterns and roles of genes, especially those involved in core developmental processes, are often conserved over vast evolutionary distances. Paradoxically, the DNA sequences surrounding these genes, which contain the cis-regulatory sequences (enhancers) that regulate gene expression, tend to be highly diverged. The manner and extent to which enhancers are functionally conserved, and how the overall organization of regulatory sequences within a locus is preserved or restructured, is not well understood. In this paper, we investigate these questions by identifying enhancers controlling expression of a master nervous system regulatory gene named sim in the mosquito Aedes aegypti, and comparing their functions and locations to those in the well-characterized sim locus of the fruit fly Drosophila melanogaster. Our results suggest that the two species generate identical patterns of sim expression through a mix of conserved and novel regulatory sequences. Continued exploration of the sim locus in these two species will help to build a comprehensive picture of how a regulatory locus for a master developmental regulator has evolved. [ABSTRACT FROM AUTHOR]

Tumor-associated macrophages (TAM) subtypes have been shown to impact cancer prognosis and resistance to immunotherapy. However, there is still a lack of systematic investigation into their molecular characteristics and clinical relevance in different cancer types. Single-cell RNA sequencing data from three different tumor types were used to cluster and type macrophages. Functional analysis and communication of TAM subpopulations were performed by Gene Ontology-Biological Process and CellChat respectively. Differential expression of characteristic genes in subpopulations was calculated using zscore as well as edgeR and Wilcoxon rank sum tests, and subsequently gene enrichment analysis of characteristic genes and anti-PD-1 resistance was performed by the REACTOME database. We revealed the heterogeneity of TAM, and identified eleven subtypes and their impact on prognosis. These subtypes expressed different molecular functions respectively, such as being involved in T cell activation, apoptosis and differentiation, or regulating viral bioprocesses or responses to viruses. The SPP1 pathway was identified as a critical mediator of communication between TAM subpopulations, as well as between TAM and epithelial cells. Macrophages with high expression of SPP1 resulted in poorer survival. By in vitro study, we showed SPP1 mediated the interactions between TAM clusters and between TAM and tumor cells. SPP1 promoted the tumor-promoting ability of TAM, and increased PDL1 expression and stemness of tumor cells. Inhibition of SPP1 attenuated N-cadherin and β-catenin expression and the activation of AKT and STAT3 pathway in tumor cells. Additionally, we found that several subpopulations could decrease the sensitivity of anti-PD-1 therapy in melanoma. SPP1 signal was a critical pathway of communication between macrophage subtypes. Some specific macrophage subtypes were associated with immunotherapy resistance and prognosis in some cancer types. Author summary: Macrophages are natural immune cell-specialized, phagocytic cells, and many studies have been conducted to analyze the functional role of macrophages, but integrating macrophages from multiple cancers to analyze the molecular functions and interactions of macrophages is poorly known.We clustered and typed macrophages using single-cell RNA sequencing data from three different tumor types. Gene Ontology-Biological Processes and Cell Chat were used for functional analysis and communication of TAM subpopulations, respectively. We found heterogeneity of TAM and identified 11 subtypes and their impact on prognosis.The SPP1 pathway is recognized as a key mediator of communication between TAM subpopulations as well as between TAM and epithelial cells, and macrophages with high expression of SPP1 have a lower survival rate. Through in vitro studies, we also found that SPP1 mediated interactions between TAM clusters and between TAM and tumor cells. SPP1 promoted the tumor-promoting ability of TAM and increased PDL1 expression and stemness of tumor cells. In addition, we performed gene enrichment analysis using the REACTOME database and identified subpopulations that decrease the sensitivity of melanoma patients to anti-PD-1 therapy. Some specific macrophage subtypes are associated with immunotherapy resistance and prognosis in certain cancer types. [ABSTRACT FROM AUTHOR]

Gene expression can be influenced by genetic variants that are closely linked to the expressed gene (cis eQTLs) and variants in other parts of the genome (trans eQTLs). We created a multiparental mapping population by sampling genotypes from a single natural population of Mimulus guttatus and scored gene expression in the leaves of 1,588 plants. We find that nearly every measured gene exhibits cis regulatory variation (91% have FDR < 0.05). cis eQTLs are usually allelic series with three or more functionally distinct alleles. The cis locus explains about two thirds of the standing genetic variance (on average) but varies among genes and tends to be greatest when there is high indel variation in the upstream regulatory region and high nucleotide diversity in the coding sequence. Despite mapping over 10,000 trans eQTL / affected gene pairs, most of the genetic variance generated by trans acting loci remains unexplained. This implies a large reservoir of trans acting genes with subtle or diffuse effects. Mapped trans eQTLs show lower allelic diversity but much higher genetic dominance than cis eQTLs. Several analyses also indicate that trans eQTLs make a substantial contribution to the genetic correlations in expression among different genes. They may thus be essential determinants of "gene expression modules," which has important implications for the evolution of gene expression and how it is studied by geneticists. Author summary: Mimulus guttatus (yellow monkeyflower) is a model for the study of quantitative trait evolution in natural populations. Research has focused mainly on whole organism traits like flower size or herbivore resistance, but the level of expression of a gene is also a quantitative trait. In this study, we dissect leaf transcriptome variation using a breeding design that estimates the contribution of individual loci to expression variation (eQTLs). We find rough agreement to the "oligogenic model" of inheritance where a major locus (the cis regulatory region) generates much of the genetic variation in the population. Associations studies usually characterize genetic effects as binary (e.g. the two alternatives at a single nucleotide polymorphism or "SNP"), but this description is insufficient for Mimulus. Most loci exhibit multiple, and in some cases, a continuum of alleles. We find that trans eQTLs have different features than cis eQTLs, both in terms of the diversity and genetic dominance of alleles. These genetic features of eQTLs are critical determinants of the "G matrix," the genetic variances and covariances among all genes which determine how gene expression will evolve under selection in response to changing environmental conditions. Our finding of large effect sizes and high allelic diversity suggests that the G matrix may be surprisingly malleable, even on ecological timescales. [ABSTRACT FROM AUTHOR]

RNA sequencing (RNA-Seq) is widely used to capture transcriptome dynamics across tissues, biological entities, and conditions. Currently, few or no methods can handle multiple biological variables (e.g., tissues/ phenotypes) and their interactions simultaneously, while also achieving dimension reduction (DR). We propose INSIDER, a general and flexible statistical framework based on matrix factorization, which is freely available at https://github.com/kai0511/insider. INSIDER decomposes variation from different biological variables and their interactions into a shared low-rank latent space. Particularly, it introduces the elastic net penalty to induce sparsity while considering the grouping effects of genes. It can achieve DR of high-dimensional data (of > = 3 dimensions), as opposed to conventional methods (e.g., PCA/NMF) which generally only handle 2D data (e.g., sample × expression). Besides, it enables computing 'adjusted' expression profiles for specific biological variables while controlling variation from other variables. INSIDER is computationally efficient and accommodates missing data. INSIDER also performed similarly or outperformed a close competing method, SDA, as shown in simulations and can handle complex missing data in RNA-Seq data. Moreover, unlike SDA, it can be used when the data cannot be structured into a tensor. Lastly, we demonstrate its usefulness via real data analysis, including clustering donors for disease subtyping, revealing neuro-development trajectory using the BrainSpan data, and uncovering biological processes contributing to variables of interest (e.g., disease status and tissue) and their interactions. Author summary: The RNA-Seq studies aim to understand the contributing role of gene activities to phenotypes, especially human disorders, reveal important biological processes (BPs) that contribute to phenotypes, and understand the involvement of tissues (e.g., brain structures/regions) in disease development. However, our ability to understand these key issues in high-dimensional settings is hampered by the limitations of currently available methods for analyzing RNA-seq data. To address these issues, we propose a novel statistical method to jointly model variation from multiple biological variables (e.g., donor, tissue, phenotype), to introduce interactions (tissue × phenotype), and to reduce dimensionality. The favorable properties of our method are discussed below. Its advantages over one closest existing computational tool and its wide applicability are shown via simulation and applications. A computationally efficient implementation of our method is provided. [ABSTRACT FROM AUTHOR]

RpoS is an alternative sigma factor needed for the induction of the general stress response in many gammaproteobacteria. Tight regulation of RpoS levels and activity is required for bacterial growth and survival under stress. In Escherichia coli, various stresses lead to higher levels of RpoS due to increased translation and decreased degradation. During non-stress conditions, RpoS is unstable, because the adaptor protein RssB delivers RpoS to the ClpXP protease. RpoS degradation is prevented during stress by the sequestration of RssB by anti-adaptors, each of which is induced in response to specific stresses. Here, we examined how the stabilization of RpoS is reversed during recovery of the cell from stress. We found that RpoS degradation quickly resumes after recovery from phosphate starvation, carbon starvation, and when transitioning from stationary phase back to exponential phase. This process is in part mediated by the anti-adaptor IraP, known to promote RpoS stabilization during phosphate starvation via the sequestration of adaptor RssB. The rapid recovery from phosphate starvation is dependent upon a feedback loop in which RpoS transcription of rssB, encoding the adaptor protein, plays a critical role. Crl, an activator of RpoS that specifically binds to and stabilizes the complex between the RNA polymerase and RpoS, is also required for the feedback loop to function efficiently, highlighting a critical role for Crl in restoring RpoS basal levels. Author summary: In their native environments, bacteria are exposed to constant changes in nutrient availability, as well as other biotic and abiotic stressors. To adjust to these changes, bacteria must rewire gene expression to adapt to or avoid stress-induced damage. A key player in the global response to general stresses is the alternative sigma factor RpoS, a promoter specificity -determining subunit of RNA polymerase. RpoS levels increase with stress, due to increased translation and stabilization of the otherwise unstable RpoS protein. Here, we examine how the cell restores homeostasis after the stress has passed. We show that a negative feedback loop in which RpoS regulates the transcription of an adaptor for proteolysis poises the cell to rapidly resume RpoS degradation upon the exit from stress. [ABSTRACT FROM AUTHOR]

A number of studies have demonstrated that epigenetic factors regulate plant developmental timing in response to environmental changes. However, we still have an incomplete view of how epigenetic factors can regulate developmental events such as organogenesis, and the transition from cell division to cell expansion, in plants. The small number of cell types and the relatively simple developmental progression required to form the Arabidopsis petal makes it a good model to investigate the molecular mechanisms driving plant organogenesis. In this study, we investigated how the RABBIT EARS (RBE) transcriptional repressor maintains the downregulation of its downstream direct target, TCP5, long after RBE expression dissipates. We showed that RBE recruits the Groucho/Tup1-like corepressor TOPLESS (TPL) to repress TCP5 transcription in petal primordia. This process involves multiple layers of changes such as remodeling of chromatin accessibility, alteration of RNA polymerase activity, and histone modifications, resulting in an epigenetic memory that is maintained through multiple cell divisions. This memory functions to maintain cell divisions during the early phase of petal development, and its attenuation in a cell division-dependent fashion later in development enables the transition from cell division to cell expansion. Overall, this study unveils a novel mechanism by which the memory of an epigenetic state, and its cell-cycle regulated decay, acts as a timer to precisely control organogenesis. Author summary: Epigenetic changes, such as histone modifications or DNA methylation, can modify gene expression without affecting the genetic code. Such alterations can persist for many cell generations, and often are induced by various environmental inputs in both plants and animals. By contrast, little is known of how intrinsic epigenetic changes can act to regulate the temporal progression of development. Here we showed that the expression of RABBIT EARS (RBE) represses the transcription of its target, TCP5, through inducing a variety of epigenetic changes at the TCP5 locus. These changes dissipated over a period of days, allowing for eventual activation of TCP5 expression and a concomitant shift from cell division to cell expansion during Arabidopsis petal development. The gradual loss of repressive epigenetic marks at the TCP5 locus represents an internal timing mechanism to temporally control the development of the Arabidopsis petal. [ABSTRACT FROM AUTHOR]

Bacteria have evolved complex transcriptional regulatory networks, as well as many diverse regulatory strategies at the RNA level, to enable more efficient use of metabolic resources and a rapid response to changing conditions. However, most RNA-based regulatory mechanisms are not well conserved across different bacterial species despite controlling genes important for virulence or essential biosynthetic processes. Here, we characterize the activity of, and assess the fitness benefit conferred by, twelve cis-acting regulatory RNAs (including several riboswitches and a T-box), in the opportunistic pathogen Streptococcus pneumoniae TIGR4. By evaluating native locus mutants of each regulator that result in constitutively active or repressed expression, we establish that growth defects in planktonic culture are associated with constitutive repression of gene expression, while constitutive activation of gene expression is rarely deleterious. In contrast, in mouse nasal carriage and pneumonia models, strains with either constitutively active and repressed gene expression are significantly less fit than matched control strains. Furthermore, two RNA-regulated pathways, FMN synthesis/transport and pyrimidine synthesis/transport display exceptional sensitivity to mis-regulation or constitutive gene repression in both planktonic culture and in vivo environments. Thus, despite lack of obvious phenotypes associated with constitutive gene expression in vitro, the fitness benefit conferred on bacteria via fine-tuned metabolic regulation through cis-acting regulatory RNAs is substantial in vivo, and therefore easily sufficient to drive the evolution and maintenance of diverse RNA regulatory mechanisms. Author summary: Regulation of gene expression is a nearly universal trait across all organisms that enables efficient utilization of resources, rapid reaction to a changing environment, or development of specialized subpopulations. However, the mechanisms that enable gene control are extraordinarily diverse. In bacteria, there are numerous regulatory mechanisms that act at the RNA level. Such mechanisms are frequently responsible for fine-tuning expression levels of critical metabolic processes or virulence factors, but are often not well conserved across divergent bacterial species. In this work, we seek to evaluate the selective pressures driving the evolution and maintenance of RNA cis-regulatory mechanisms. Using the opportunistic pathogen Streptococcus pneumoniae as a model, we measured the fitness of strains carrying native locus mutations to RNA regulators that result in constitutively active and repressed gene expression. We find that only constitutively repressed mutants show significant growth defects under stringent in vitro culture conditions, but both constitutively active and repressed mutants are significantly less fit within mouse models of infection and carriage. We also find that specific pathways such as pyrimidine and FMN synthesis/transport are very sensitive to mis-regulation and constitutive repression respectively. This study demonstrates that while individual variants may show mild fitness defects, when considered in aggregate, these defects correspond to significantly depressed fitness values that are more than sufficient to drive evolutionary trajectory for organisms with large population sizes such as S. pneumoniae. Thus, RNA-regulators responsible for fine-tuning metabolic genes are likely critical for S. pneumoniae success in a host. [ABSTRACT FROM AUTHOR]

TOP2 inhibitors (TOP2i) are effective drugs for breast cancer treatment. However, they can cause cardiotoxicity in some women. The most widely used TOP2i include anthracyclines (AC) Doxorubicin (DOX), Daunorubicin (DNR), Epirubicin (EPI), and the anthraquinone Mitoxantrone (MTX). It is unclear whether women would experience the same adverse effects from all drugs in this class, or if specific drugs would be preferable for certain individuals based on their cardiotoxicity risk profile. To investigate this, we studied the effects of treatment of DOX, DNR, EPI, MTX, and an unrelated monoclonal antibody Trastuzumab (TRZ) on iPSC-derived cardiomyocytes (iPSC-CMs) from six healthy females. All TOP2i induce cell death at concentrations observed in cancer patient serum, while TRZ does not. A sub-lethal dose of all TOP2i induces limited cellular stress but affects calcium handling, a function critical for cardiomyocyte contraction. TOP2i induce thousands of gene expression changes over time, giving rise to four distinct gene expression response signatures, denoted as TOP2i early-acute, early-sustained, and late response genes, and non-response genes. There is no drug- or AC-specific signature. TOP2i early response genes are enriched in chromatin regulators, which mediate AC sensitivity across breast cancer patients. However, there is increased transcriptional variability between individuals following AC treatments. To investigate potential genetic effects on response variability, we first identified a reported set of expression quantitative trait loci (eQTLs) uncovered following DOX treatment in iPSC-CMs. Indeed, DOX response eQTLs are enriched in genes that respond to all TOP2i. Next, we identified 38 genes in loci associated with AC toxicity by GWAS or TWAS. Two thirds of the genes that respond to at least one TOP2i, respond to all ACs with the same direction of effect. Our data demonstrate that TOP2i induce thousands of shared gene expression changes in cardiomyocytes, including genes near SNPs associated with inter-individual variation in response to DOX treatment and AC-induced cardiotoxicity. Author summary: Anthracycline drugs such as Doxorubicin are effective treatments for breast cancer; however, they can cause cardiotoxicity in some women. It is unclear whether women would experience the same toxicity for all drugs in this class, or whether specific drugs would be better tolerated in specific individuals. We used an in vitro system of induced pluripotent stem cell-derived cardiomyocytes from six healthy females to test the effects of five breast cancer drugs on cell heath and global gene expression. We identified a strong shared cellular and gene expression response to drugs from the same class. However, there is more variation in gene expression levels between individuals following treatment with each anthracycline compared to untreated cells. We found that many genes in regions previously associated with Doxorubicin-induced cardiotoxicity in cancer patients, respond to at least two drugs in the class. This suggests that drugs in the same class induce similar effects on an individual's heart. This work contributes to our understanding of how drug response, in the context of off-target effects, varies across individuals. [ABSTRACT FROM AUTHOR]

Common genetic variants in the repressive GATA-family transcription factor (TF) TRPS1 locus are associated with breast cancer risk, and luminal breast cancer cell lines are particularly sensitive to TRPS1 knockout. We introduced an inducible degron tag into the native TRPS1 locus within a luminal breast cancer cell line to identify the direct targets of TRPS1 and determine how TRPS1 mechanistically regulates gene expression. We acutely deplete over 80 percent of TRPS1 from chromatin within 30 minutes of inducing degradation. We find that TRPS1 regulates transcription of hundreds of genes, including those related to estrogen signaling. TRPS1 directly regulates chromatin structure, which causes estrogen receptor alpha (ER) to redistribute in the genome. ER redistribution leads to both repression and activation of dozens of ER target genes. Downstream from these primary effects, TRPS1 depletion represses cell cycle-related gene sets and reduces cell doubling rate. Finally, we show that high TRPS1 activity, calculated using a gene expression signature defined by primary TRPS1-regulated genes, is associated with worse breast cancer patient prognosis. Taken together, these data suggest a model in which TRPS1 modulates the genomic distribution of ER, both activating and repressing transcription of genes related to cancer cell fitness. Author summary: Breast cancer is the most common cancer among women. The majority of cases are luminal, which tend to be estrogen receptor alpha (ER)-positive. ER is well-studied among transcription factors (TFs) because it is ligand-activated. This allows for the rapid induction of ER activity and the identification of primary estrogen-responsive genes. Most TFs have not been so extensively characterized, because their activity is not so rapidly perturbable. TRPS1 is an atypical GATA family TF that is associated with corepressor complexes and transcriptional repression. Here, we use an inducible degron tag system to rapidly deplete endogenously-tagged TRPS1 in luminal breast cancer cells within 30 minutes. We find that TRPS1 directly decreases local chromatin accessibility. This decreases ER binding intensity at TRPS1-proximal ER binding sites. As an indirect effect, ER binding intensity distal to TRPS1 increases in intensity. These effects on ER binding are associated with changes in ER target gene transcription, repressing or activating these genes in concordance with the effect on ER binding intensity. Our work is consistent with a model in which TFs either exclusively activate or exclusively repress transcription of their direct target genes, with indirect primary response genes changing due to the redistribution of limiting activating TFs or coactivators. [ABSTRACT FROM AUTHOR]

Elucidating gene function is a major goal in biology, especially among non-model organisms. However, doing so is complicated by the fact that molecular conservation does not always mirror functional conservation, and that complex relationships among genes are responsible for encoding pathways and higher-order biological processes. Co-expression, a promising approach for predicting gene function, relies on the general principal that genes with similar expression patterns across multiple conditions will likely be involved in the same biological process. For Cryptococcus neoformans, a prevalent human fungal pathogen greatly diverged from model yeasts, approximately 60% of the predicted genes in the genome lack functional annotations. Here, we leveraged a large amount of publicly available transcriptomic data to generate a C. neoformans Co-Expression Network (CryptoCEN), successfully recapitulating known protein networks, predicting gene function, and enabling insights into the principles influencing co-expression. With 100% predictive accuracy, we used CryptoCEN to identify 13 new DNA damage response genes, underscoring the utility of guilt-by-association for determining gene function. Overall, co-expression is a powerful tool for uncovering gene function, and decreases the experimental tests needed to identify functions for currently under-annotated genes. Author summary: A central problem in genetics is the connection between genotype and phenotype. Computational approaches to predict gene function can be especially useful for non-model organisms where extensive functional testing has not yet been performed. Co-expression to predict gene function is based on the principle that genes that share similar expression patterns across multiple environmental conditions or perturbations are likely to be involved in the same biological process. Here, we collected transcriptomic data from the Cryptococcus neoformans field and built a robust co-expression network for predicting gene function, especially biological process information. Not only are we able to use this network for retrospective analysis of known gene clusters, but we are also able to make prospective predictions about gene function, including the well-studied processes of capsule and ergosterol biosynthesis. We also discovered a new role for 13 genes in the response to DNA damaging agents, showing that co-expression can reveal new players in conserved biological processes. [ABSTRACT FROM AUTHOR]

Changes in gene regulation have long been appreciated as a driving force of adaptive evolution, however the relative contributions of cis- and trans-acting changes to gene regulation over short evolutionary timescales remain unclear. Instances of recent, parallel phenotypic evolution provide an opportunity to assess whether parallel patterns are seen at the level of gene expression, and to assess the relative contribution of cis- and trans- changes to gene regulation in the early stages of divergence. Here, we studied gene expression in liver and brown adipose tissue in two wild-derived strains of house mice that independently adapted to cold, northern environments, and we compared them to a strain of house mice from a warm, tropical environment. To investigate gene regulatory evolution, we studied expression in parents and allele-specific expression in F1 hybrids of crosses between warm-adapted and cold-adapted strains. First, we found that the different cold-adapted mice showed both unique and shared changes in expression, but that the proportion of shared changes (i.e. parallelism) was greater than expected by chance. Second, we discovered that expression evolution occurred largely at tissue-specific and cis-regulated genes, and that these genes were over-represented in parallel cases of evolution. Finally, we integrated the expression data with scans for selection in natural populations and found substantial parallelism in the two northern populations for genes under selection. Furthermore, selection outliers were associated with cis-regulated genes more than expected by chance; cis-regulated genes under selection influenced phenotypes such as body size, immune functioning, and activity level. These results demonstrate that parallel patterns of gene expression in mice that have independently adapted to cold environments are driven largely by tissue-specific and cis-regulatory changes, providing insight into the mechanisms of adaptive gene regulatory evolution at the earliest stages of divergence. Author summary: The parallel movement of organisms into novel environments provides an opportunity to understand the molecular basis of adaptation and the repeatability of this process. Mutations affecting the expression of genes are known to underlie much of adaptive evolution. Such mutations can arise in cis- (near the gene of interest) or in trans- (at a distant locus), but the relative contribution of these different kinds of changes to adaptation is poorly understood, especially during very recent divergence. Here, we compared evolved gene expression differences between a warm-adapted house mouse strain and two different cold-adapted strains that have independently evolved similar phenotypic traits, such as increased body size and decreased extremity length during the last few hundred years. Using crosses between warm-adapted and cold-adapted mice, we found that mutations acting in a context specific manner (cis-regulatory and tissue-specific changes) predominate expression divergence and are more likely to be involved in parallel evolution. We used population level selection scans in wild animals to identify regions of the genome under selection and combined these findings with the gene expression data to identify candidate genes underlying adaptation to novel environments. Together, our work describes the gene regulatory dynamics of rapid environmental adaptation, and the repeatability of these patterns over multiple instances of adaptation. [ABSTRACT FROM AUTHOR]

Genome evolution is partly driven by the mobility of transposable elements (TEs) which often leads to deleterious effects, but their activity can also facilitate genetic novelty and catalyze local adaptation. We explored how the intraspecific diversity of TE polymorphisms might contribute to the broad geographic success and adaptive capacity of the emerging oil crop Thlaspi arvense (field pennycress). We classified the TE inventory based on a high-quality genome assembly, estimated the age of retrotransposon TE families and comprehensively assessed their mobilization potential. A survey of 280 accessions from 12 regions across the Northern hemisphere allowed us to quantify over 90,000 TE insertion polymorphisms (TIPs). Their distribution mirrored the genetic differentiation as measured by single nucleotide polymorphisms (SNPs). The number and types of mobile TE families vary substantially across populations, but there are also shared patterns common to all accessions. Ty3/Athila elements are the main drivers of TE diversity in T. arvense populations, while a single Ty1/Alesia lineage might be particularly important for transcriptome divergence. The number of retrotransposon TIPs is associated with variation at genes related to epigenetic regulation, including an apparent knockout mutation in BROMODOMAIN AND ATPase DOMAIN-CONTAINING PROTEIN 1 (BRAT1), while DNA transposons are associated with variation at the HSP19 heat shock protein gene. We propose that the high rate of mobilization activity can be harnessed for targeted gene expression diversification, which may ultimately present a toolbox for the potential use of transposition in breeding and domestication of T. arvense. Author summary: Transposable elements (TEs) are often considered genomic parasites, but they can also generate phenotypic novelty that helps organisms to adapt to new environments. To understand how TEs might contribute to phenotypic diversity and adaptive potential in the emerging oilseed crop Thlaspi arvense (field pennycress), we examined the dynamics of TE variation in a geographically diverse sample of this species. By surveying almost 300 wild accessions from North America and Eurasia we discovered over 90,000 polymorphic TE insertions. We identified not only genetic factors that vary between populations and that are associated with TE mobilization, but also TE families that are most likely to generate genetic diversity of interest to breeders. [ABSTRACT FROM AUTHOR]

Genome-wide association studies (GWASs) have achieved remarkable success in associating thousands of genetic variants with complex traits. However, the presence of linkage disequilibrium (LD) makes it challenging to identify the causal variants. To address this critical gap from association to causation, many fine-mapping methods have been proposed to assign well-calibrated probabilities of causality to candidate variants, taking into account the underlying LD pattern. In this manuscript, we introduce a statistical framework that incorporates expression quantitative trait locus (eQTL) information to fine-mapping, built on the sum of single-effects (SuSiE) regression model. Our new method, SuSiE2, connects two SuSiE models, one for eQTL analysis and one for genetic fine-mapping. This is achieved by first computing the posterior inclusion probabilities (PIPs) from an eQTL-based SuSiE model with the expression level of the candidate gene as the phenotype. These calculated PIPs are then utilized as prior inclusion probabilities for risk variants in another SuSiE model for the trait of interest. By prioritizing functional variants within the candidate region using eQTL information, SuSiE2 improves SuSiE by increasing the detection rate of causal SNPs and reducing the average size of credible sets. We compared the performance of SuSiE2 with other multi-trait fine-mapping methods with respect to power, coverage, and precision through simulations and applications to the GWAS results of Alzheimer's disease (AD) and body mass index (BMI). Our results demonstrate the better performance of SuSiE2, both when the in-sample linkage disequilibrium (LD) matrix and an external reference panel is used in inference. Author summary: Genome-wide association studies (GWASs) have proven powerful in detecting genetic variants associated with complex traits. However, there are challenges in distinguishing the causal variants from other variants strongly correlated with them. To better identify causal SNPs, many fine-mapping methods have been proposed to assign well-calibrated probabilities of causality to candidate variants. We introduce a statistical framework that incorporates expression quantitative trait locus (eQTL) information to fine-mapping, which can improve the accuracy and efficiency of association studies by prioritizing functional variants within the risk genes before evaluating the causation. Our new fine-mapping framework, SuSiE2, connects two sum of single-effects (SuSiE) models, one for eQTL analysis and one for genetic fine-mapping. The posterior inclusion probabilities from an eQTL-based SuSiE model are utilized as prior inclusion probabilities for risk variants in another SuSiE model for the trait of interest. Through simulations and real data analyses focused on body mass index and Alzheimer's disease, we demonstrate that SuSiE2 improves fine-mapping results by increasing statistical power, having appropriate coverage, and reducing the average size of credible sets. [ABSTRACT FROM AUTHOR]

Eukaryotic chromatin is organized into either silenced heterochromatin or relaxed euchromatin regions, which controls the accessibility of transcriptional machinery and thus regulates gene expression. In fission yeast, Schizosaccharomyces pombe, Set1 is the sole H3K4 methyltransferase and is mainly enriched at the promoters of actively transcribed genes. In contrast, Clr4 methyltransferase initiates H3K9 methylation, which has long been regarded as a hallmark of heterochromatic silencing. Lsd1 and Lsd2 are two highly conserved H3K4 and H3K9 demethylases. As these histone-modifying enzymes perform critical roles in maintaining histone methylation patterns and, consequently, gene expression profiles, cross-regulations among these enzymes are part of the complex regulatory networks. Thus, elucidating the mechanisms that govern their signaling and mutual regulations remains crucial. Here, we demonstrated that C-terminal truncation mutants, lsd1-ΔHMG and lsd2-ΔC, do not compromise the integrity of the Lsd1/2 complex but impair their chromatin-binding capacity at the promoter region of target genomic loci. We identified protein-protein interactions between Lsd1/2 and Raf2 or Swd2, which are the subunits of the Clr4 complex (CLRC) and Set1-associated complex (COMPASS), respectively. We showed that Clr4 and Set1 modulate the protein levels of Lsd1 and Lsd2 in opposite ways through the ubiquitin-proteasome-dependent pathway. During heat stress, the protein levels of Lsd1 and Lsd2 are upregulated in a Set1-dependent manner. The increase in protein levels is crucial for differential gene expression under stress conditions. Together, our results support a cross-regulatory model by which Set1 and Clr4 methyltransferases control the protein levels of Lsd1/2 demethylases to shape the dynamic chromatin landscape. Author summary: Histone-modifying enzymes make covalent modifications to histones. These modifications act like chemical tags that can either tighten or loosen the DNA structure, affecting whether genes are turned on or off. In fission yeast, Set1 methylates histone H3 lysine 4, which marks loosely packed DNA and is associated with gene activation; while Clr4 methylates histone H3 lysine 9, which represents a hallmark of the tightly packed DNA and is associated with gene silencing. Here, we show a regulatory relationship between these two enzymes and with two lysine-specific demethylases (Lsd1/2), which can remove the methyl tags added by Clr4 and Set1. Clr4 and Set1 have opposite effects on Lsd1 and Lsd2 protein levels. Clr4 reduces the levels of Lsd1/2, while Set1 promotes their stability. This control over the levels of Lsd1/2 is achieved through the ubiquitin-proteasome-dependent pathway. By studying these interactions, we have uncovered a novel regulatory mechanism that helps fission yeast maintain a balanced level of these active/repressive histone-modifying enzymes. Understanding these intricate regulatory networks offers important insights into gene control and enhances our comprehension of the complex interactions within cells. [ABSTRACT FROM AUTHOR]

The Muscleblind-like (MBNL) family is a highly conserved set of RNA-binding proteins (RBPs) that regulate RNA metabolism during the differentiation of various animal tissues. Functional insufficiency of MBNL affects muscle and central nervous system development, and contributes to the myotonic dystrophies (DM), a set of incurable multisystemic disorders. Studies on the regulation of MBNL genes are essential to provide insight into the gene regulatory networks controlled by MBNL proteins and to understand how dysregulation within these networks causes disease. In this study, we demonstrate the evolutionary conservation of an autoregulatory mechanism that governs the function of MBNL proteins by generating two distinct protein isoform types through alternative splicing. Our aim was to further our understanding of the regulatory principles that underlie this conserved feedback loop in a whole-organismal context, and to address the biological significance of the respective isoforms. Using an alternative splicing reporter, our studies show that, during development of the Caenorhabditis elegans central nervous system, the orthologous mbl-1 gene shifts production from long protein isoforms that localize to the nucleus to short isoforms that also localize to the cytoplasm. Using isoform-specific CRISPR/Cas9-generated strains, we showed that expression of short MBL-1 protein isoforms is required for healthy neuromuscular function and neurodevelopment, while expression of long MBL-1 protein isoforms is dispensable, emphasizing a key role for cytoplasmic functionalities of the MBL-1 protein. Furthermore, RNA-seq and lifespan analyses indicated that short MBL-1 isoforms are crucial regulators of miRNA expression and, in consequence, required for normal lifespan. In conclusion, this study provides support for the disruption of cytoplasmic RNA metabolism as a contributor in myotonic dystrophy and paves the way for further exploration of miRNA regulation through MBNL proteins during development and in disease models. Author summary: Alternative splicing is a molecular process by which one gene can produce multiple protein isoforms. These isoforms often have different properties, and normal cell function and development requires expression of the appropriate isoform. The regulation of alternative splicing is carried out by a set of specialized proteins, one of which is the Muscleblind-like protein (MBNL). This protein is essential during the development and differentiation of many animal tissues. We found that this protein regulates the alternative splicing of its own gene in animals as divergent as nematode and human. We showed that this leads to the timely expression of MBNL protein isoforms that localize to the cytoplasm of the cell. In the cytoplasm, these isoforms acquire new functions unrelated to splicing and regulate molecular processes that control the expression of many other genes. These functions are required for normal development of the nervous system and lifespan in nematodes. This finding is particularly important because it helps us understand the regulation and function of MBNL proteins in animal development as well as identify potential disease mechanisms in myotonic dystrophies, neurodegenerative disorders associated with MBNL dysfunction. [ABSTRACT FROM AUTHOR]

Genetic networks are surprisingly robust to perturbations caused by new mutations. This robustness is conferred in part by compensation for loss of a gene's activity by genes with overlapping functions, such as paralogs. Compensation occurs passively when the normal activity of one paralog can compensate for the loss of the other, or actively when a change in one paralog's expression, localization, or activity is required to compensate for loss of the other. The mechanisms of active compensation remain poorly understood in most cases. Here we investigate active compensation for the loss or reduction in expression of the Saccharomyces cerevisiae gene TDH3 by its paralog TDH2. TDH2 is upregulated in a dose-dependent manner in response to reductions in TDH3 by a mechanism requiring the shared transcriptional regulators Gcr1p and Rap1p. TDH1, a second and more distantly related paralog of TDH3, has diverged in its regulation and is upregulated by another mechanism. Other glycolytic genes regulated by Rap1p and Gcr1p show changes in expression similar to TDH2, suggesting that the active compensation by TDH3 paralogs is part of a broader homeostatic response mediated by shared transcriptional regulators. Author summary: Living things have a remarkable ability to modify their molecular and cellular processes to reduce the negative impacts of genetic and environmental perturbations. Here, we examine how cells of baker's yeast respond to mutations changing the expression of a key metabolic gene, which affects cell division rates (i.e., fitness). We find that the expression of other metabolic genes, including a gene with ancestrally similar sequence, function, and regulation (a paralog), changes its expression in ways complementary to the mutated gene, reducing the impact of the mutation on fitness. This type of compensatory expression change has also been reported for other genes, but the molecular mechanisms behind such changes are generally unknown. In this case, we find that upregulation of the paralog is mediated by transcription factors that also directly regulate the mutated gene. This observation suggests that ancestrally shared regulation of paralogous genes can play a key role in homeostasis and might help maintain genes in the genome following gene duplication. [ABSTRACT FROM AUTHOR]

For understanding the evolutionary mechanism of sexually selected exaggerated traits, it is essential to uncover its molecular basis. By using broad-horned flour beetle that has male-specific exaggerated structures (mandibular horn, head horn and gena enlargement), we investigated the transcriptomic and functional characters of sex-biased genes. Comparative transcriptome of male vs. female prepupal heads elucidated 673 sex-biased genes. Counter-intuitively, majority of them were female-biased (584 genes), and GO enrichment analysis showed cell-adhesion molecules were frequently female-biased. This pattern motivated us to hypothesize that female-biased transcripts (i.e. the transcripts diminished in males) may play a role in outgrowth formation. Potentially, female-biased genes may act as suppressors of weapon structure. In order to test the functionality of female-biased genes, we performed RNAi-mediated functional screening for top 20 female-biased genes and 3 genes in the most enriched GO term (cell-cell adhesion, fat1/2/3, fat4 and dachsous). Knockdown of one transcription factor, zinc finger protein 608 (zfp608) resulted in the formation of male-like gena in females, supporting the outgrowth suppression function of this gene. Similarly, knockdown of fat4 induced rudimental, abnormal mandibular horn in female. fat1/2/3RNAi, fat4RNAi and dachsousRNAi males exhibited thick and/or short mandibular horns and legs. These cell adhesion molecules are known to regulate tissue growth direction and known to be involved in the weapon formation in Scarabaeoidea beetles. Functional evidence in phylogenetically distant broad-horned flour beetle suggest that cell adhesion genes are repeatedly deployed in the acquisition of outgrowth. In conclusion, this study clarified the overlooked functions of female-biased genes in weapon development. Author summary: A diverse range of animals have evolved the male sexual weapons that provide an advantage in male-male competition. Understanding the genetic mechanisms behind the sexually dimorphic expression of these weapons is crucial for comprehending their evolution. However, this area of study remains largely unexplored. By comparing the gene expression in the heads of male and female prepupae, we identified 673 sex-biased genes, with 584 genes being highly expressed in females. Through knockdown experiments targeting 23 female-biased genes, we discovered that four genes play a role in determining the size and expression of male-exaggerated traits. Intriguingly, two of these genes, including a transcription factor and a cell adhesion molecule, were found to suppress weapon development in females, implying that down-regulation of these genes in males enhances weapon development. These findings illuminate the often-overlooked function of female-biased genes in male weapon development and propose the idea that the evolution of weapons involves the down-regulation of certain genes in males. [ABSTRACT FROM AUTHOR]

Mycobacterium abscessus is intrinsically resistant to antibiotics effective against other pathogenic mycobacteria largely due to the drug-induced expression of genes that confer resistance. WhiB7 is a major hub controlling the induction of resistance to ribosome-targeting antibiotics. It activates the expression of >100 genes, 7 of which are known determinants of drug resistance; the function of most genes within the regulon is however unknown, but some conceivably encode additional mechanisms of resistance. Furthermore, the hierarchy of gene expression within the regulon, if any, is poorly understood. In the present work we have identified 56 WhiB7 binding sites using chromatin immunoprecipitation sequencing (CHIP-Seq) which accounts for the WhiB7-dependent upregulation of 72 genes, and find that M. abscessus WhiB7 functions exclusively as a transcriptional activator at promoters recognized by σA/σB. We have investigated the role of 18 WhiB7 regulated genes in drug resistance. Our results suggest that while some genes within the regulon (eg. erm41, hflX, eis2 and the ABCFs) play a major role in resistance, others make smaller contributions (eg. MAB_4324c and MAB_1409c) and the observed hypersensitivity ΔMabwhiB7 is a cumulative effect of these individual contributions. Moreover, our CHIP-Seq data implicate additional roles of WhiB7 induced genes beyond antibiotic resistance. Finally, we identify a σH-dependent network in aminoglycoside and tigecycline resistance which is induced upon drug exposure and is further activated by WhiB7 demonstrating the existence of a crosstalk between components of the WhiB7-dependent and -independent circuits. Author summary: The induction of multiple genes that confer resistance to structurally diverse ribosome-targeting antibiotics is funneled through the induction of a single transcriptional activator, WhiB7, by antibiotic- stalled ribosomes. This poses a severe restriction in M. abscessus therapy as treatment with one ribosome-targeting antibiotic confers resistance to all other ribosome-targeting antibiotics. Here we uncover the intricacies of the WhiB7 regulatory circuit, identify three previously unknown determinants of aminoglycoside resistance and unveil a communication between WhiB7 dependent and independent components. This not only expands our understanding of the antibiotic resistance potential of M. abscessus but can also inform the development of much needed therapeutic options. [ABSTRACT FROM AUTHOR]

The Candida albicans genome contains between ten and fifteen distinct TLO genes that all encode a Med2 subunit of Mediator. In order to investigate the biological role of Med2/Tlo in C. albicans we deleted all fourteen TLO genes using CRISPR-Cas9 mutagenesis. ChIP-seq analysis showed that RNAP II localized to 55% fewer genes in the tloΔ mutant strain compared to the parent, while RNA-seq analysis showed that the tloΔ mutant exhibited differential expression of genes required for carbohydrate metabolism, stress responses, white-opaque switching and filamentous growth. Consequently, the tloΔ mutant grows poorly in glucose- and galactose-containing media, is unable to grow as true hyphae, is more sensitive to oxidative stress and is less virulent in the wax worm infection model. Reintegration of genes representative of the α-, β- and γ-TLO clades resulted in the complementation of the mutant phenotypes, but to different degrees. TLOα1 could restore phenotypes and gene expression patterns similar to wild-type and was the strongest activator of glycolytic and Tye7-regulated gene expression. In contrast, the two γ-TLO genes examined (i.e., TLOγ5 and TLOγ11) had a far lower impact on complementing phenotypic and transcriptomic changes. Uniquely, expression of TLOβ2 in the tloΔ mutant stimulated filamentous growth in YEPD medium and this phenotype was enhanced when Tloβ2 expression was increased to levels far in excess of Med3. In contrast, expression of reintegrated TLO genes in a tloΔ/med3Δ double mutant background failed to restore any of the phenotypes tested, suggesting that complementation of these Tlo-regulated processes requires a functional Mediator tail module. Together, these data confirm the importance of Med2/Tlo in a wide range of C. albicans cellular activities and demonstrate functional diversity within the gene family which may contribute to the success of this yeast as a coloniser and pathogen of humans. Author summary: Fungi are an important cause of infectious disease in humans. One of the most potent causes of fungal infections is the single-celled yeast species, Candida albicans, which usually lives in balance as part of the microbial communities that inhabit the gastrointestinal and vaginal tract. However, under certain conditions, for example in immunocompromised hosts, the yeast can overgrow and cause a range of infections ranging from superficial to life-threatening diseases. Our understanding of this organism and how it causes disease is still rudimentary. One of the most unique aspects of the genome of this species is the presence of an expanded family of fourteen TLO genes that encode the Med2 component of Mediator, a multi-peptide complex involved in the transcription of genes. Using CRISPR-Cas9-mediated mutagenesis we have successfully deleted all fourteen C. albicans TLO genes. Molecular and phenotypic analysis of the mutant and strains expressing representative TLO genes, directly confirms that the TLO/MED2 gene family is required for a wide range of important cell functions, including carbohydrate metabolism, cell morphology, stress tolerance and virulence, and that there is functional diversity within the proteins encoded by individual TLO genes. [ABSTRACT FROM AUTHOR]

Haploid males of hymenopteran species produce gametes through an abortive meiosis I followed by meiosis II that can either be symmetric or asymmetric in different species. Thus, one spermatocyte could give rise to two spermatids with either equal or unequal amounts of cytoplasm. It is currently unknown what molecular features accompany these postmeiotic sperm cells especially in species with asymmetric meiosis II such as bees. Here we present testis single-cell RNA sequencing datasets from the honeybee (Apis mellifera) drones of 3 and 14 days after emergence (3d and 14d). We show that, while 3d testes exhibit active, ongoing spermatogenesis, 14d testes only have late-stage spermatids. We identify a postmeiotic bifurcation in the transcriptional roadmap during spermatogenesis, with cells progressing toward the annotated spermatids (SPT) and small spermatids (sSPT), respectively. Despite an overall similarity in their transcriptomic profiles, sSPTs express the fewest genes and the least RNA content among all the sperm cell types. Intriguingly, sSPTs exhibit a relatively high expression level for Hymenoptera-restricted genes and a high mutation load, suggesting that the special meiosis II during spermatogenesis in the honeybee is accompanied by phylogenetically young gene activities. Author summary: Previous histological studies suggest that spermatogenesis is concluded before adult emergence in the honeybee. We complement this histological atlas on testis development and spermatogenesis with a comprehensive dataset from the honeybee adult testes, including single-cell RNA, long-read RNA, bulk mRNA and DNA-seq data. Through the use of gene orthology between the honeybee (Apis mellifera) and Drosophila melanogaster, we show that in newly emerged drones, spermatogenesis is an ongoing process and meiosis II results in two separable clusters of spermatid cells. One of them is identified to contain small spermatids, which are consistent with previous anatomical detection of polar body-like cells in close association with regular spermatids. The characteristic features of small spermatids include an overall low transcription level, high expression of evolutionarily young genes, an enrichment of gene activities regulated by ETS-domain transcription factors, and a high level of mutation load. This work establishes the foundation for future investigations into molecular mechanisms of sperm quality assurance in the honeybee. [ABSTRACT FROM AUTHOR]

A prominent source of mutation in cancer is single-stranded DNA cytosine deamination by cellular APOBEC3 enzymes, which results in signature C-to-T and C-to-G mutations in TCA and TCT motifs. Although multiple enzymes have been implicated, reports conflict and it is unclear which protein(s) are responsible. Here we report the development of a selectable system to quantify genome mutation and demonstrate its utility by comparing the mutagenic activities of three leading candidates—APOBEC3A, APOBEC3B, and APOBEC3H. The human cell line, HAP1, is engineered to express the thymidine kinase (TK) gene of HSV-1, which confers sensitivity to ganciclovir. Expression of APOBEC3A and APOBEC3B, but not catalytic mutant controls or APOBEC3H, triggers increased frequencies of TK mutation and similar TC-biased cytosine mutation profiles in the selectable TK reporter gene. Whole genome sequences from independent clones enabled an analysis of thousands of single base substitution mutations and extraction of local sequence preferences with APOBEC3A preferring YTCW motifs 70% of the time and APOBEC3B 50% of the time (Y = C/T; W = A/T). Signature comparisons with breast tumor whole genome sequences indicate that most malignancies manifest intermediate percentages of APOBEC3 signature mutations in YTCW motifs, mostly between 50 and 70%, suggesting that both enzymes contribute in a combinatorial manner to the overall mutation landscape. Although the vast majority of APOBEC3A- and APOBEC3B-induced single base substitution mutations occur outside of predicted chromosomal DNA hairpin structures, whole genome sequence analyses and supporting biochemical studies also indicate that both enzymes are capable of deaminating the single-stranded loop regions of DNA hairpins at elevated rates. These studies combine to help resolve a long-standing etiologic debate on the source of APOBEC3 signature mutations in cancer and indicate that future diagnostic and therapeutic efforts should focus on both APOBEC3A and APOBEC3B. Author summary: A large source of mutation in cancer is attributable to single-stranded DNA cytosine to uracil (C-to-U) deamination by APOBEC3 family enzymes. Chromosomal DNA U's directly template the insertion of adenine (A), which results in C/G-to-T/A mutations. Chromosomal U's also trigger error-prone DNA repair processes, which less directly result in C/G-to-T/A mutations, C/G-to-G/C mutations, and more complex events such as insertion/deletion mutations. We report a human cell line for cancer mutation studies, where mutations can be identified selectively or by genomic DNA sequencing. Only APOBEC3A and APOBEC3B proved capable of inflicting the aforementioned types of single base substitution mutations. Interestingly, 70% and 50% of the single base substitution mutations attributable to APOBEC3A and APOBEC3B, respectively, occurred in TCA and TCT trinucleotide motifs preceded by a pyrimidine (T or C). This is significant because most APOBEC3-positive breast tumors exhibit an intermediate percentage of signature mutations in these broader tetranucleotide motifs, which supports the idea that both APOBEC3A and APOBEC3B contribute in a combinatorial manner to the overall mutation landscape in these tumors. [ABSTRACT FROM AUTHOR]

MicroRNAs (miRNAs) regulate gene expression by base-pairing to target sequences in messenger RNAs (mRNAs) and recruiting factors that induce translational repression and mRNA decay. In animals, nucleotides 2–8 at the 5' end of the miRNA, called the seed region, are often necessary and sometimes sufficient for functional target interactions. MiRNAs that contain identical seed sequences are grouped into families where individual members have the potential to share targets and act redundantly. A rare exception seemed to be the miR-238/239ab family in Caenorhabditis elegans, as previous work indicated that loss of miR-238 reduced lifespan while deletion of the miR-239ab locus resulted in enhanced longevity and thermal stress resistance. Here, we re-examined these potentially opposing roles using new strains that individually disrupt each miRNA sister. We confirmed that loss of miR-238 is associated with a shortened lifespan but could detect no longevity or stress phenotypes in animals lacking miR-239a or miR-239b, individually or in combination. Additionally, dozens of genes were mis-regulated in miR-238 mutants but almost no gene expression changes were detected in either miR-239a or miR-239b mutants compared to wild type animals. We present evidence that the lack of redundancy between miR-238 and miR-239ab is independent of their sequence differences; miR-239a or miR-239b could substitute for the longevity role of miR-238 when expressed from the miR-238 locus. Altogether, these studies disqualify miR-239ab as negative regulators of aging and demonstrate that expression, not sequence, dictates the specific role of miR-238 in promoting longevity. Author summary: MicroRNAs (miRNAs) are tiny non-coding RNAs that function in diverse biological pathways. To exert their regulatory influence, miRNAs bind to specific target RNAs through partial base-pairing. A critical aspect of this miRNA-target engagement is the seed sequence, nucleotides 2–8 of the miRNA. MiRNAs that share seed sequences are grouped into families and presumed to have similar functions. Yet, other factors, such as non-seed sequences in the miRNA and its expression level, can also contribute to target regulation and result in distinct roles for miRNAs within a family. To better understand how miRNA family members can have specific functions, we focused on miR-238 and its sisters, miR-239a and miR-239b, because these miRNAs had previously been reported to play opposing longevity roles in the nematode C. elegans. Using new genetic tools, we found that loss of miR-238 alone leads to the misregulation of many genes and a reduced lifespan. However, the lack of miR-239a, miR-239b, or both sisters had almost no effect on gene expression or longevity compared to wild type animals. Strikingly, though, miR-239a or miR-239b could substitute for the aging role of miR-238 when expressed from the miR-238 locus. Thus, expression, not sequence, is the predominant distinguishing feature of mir-238 that bestows upon it a role in aging not shared with the other family members. [ABSTRACT FROM AUTHOR]

Hirschsprung disease (HSCR) is associated with deficiency of the receptor tyrosine kinase RET, resulting in loss of cells of the enteric nervous system (ENS) during fetal gut development. The major contribution to HSCR risk is from common sequence variants in RET enhancers with additional risk from rare coding variants in many genes. Here, we demonstrate that these RET enhancer variants specifically alter the human fetal gut development program through significant decreases in gene expression of RET, members of the RET-EDNRB gene regulatory network (GRN), other HSCR genes, with an altered transcriptome of 2,382 differentially expressed genes across diverse neuronal and mesenchymal functions. A parsimonious hypothesis for these results is that beyond RET's direct effect on its GRN, it also has a major role in enteric neural crest-derived cell (ENCDC) precursor proliferation, its deficiency reducing ENCDCs with relative expansion of non-ENCDC cells. Thus, genes reducing RET proliferative activity can potentially cause HSCR. One such class is the 23 RET-dependent transcription factors enriched in early gut development. We show that their knockdown in human neuroblastoma SK-N-SH cells reduces RET and/or EDNRB gene expression, expanding the RET-EDNRB GRN. The human embryos we studied had major remodeling of the gut transcriptome but were unlikely to have had HSCR: thus, genetic or epigenetic changes in addition to those in RET are required for aganglionosis. Author summary: Common enhancer variants are pervasive in the human genome and has been associated with multiple traits and disorders. Yet our knowledge of how these variants individually or commutatively affect gene expression leading to final phenotype is incomplete. Hirschsprung disease (HSCR) a common developmental disorder of the enteric nervous system is caused due to genetic dysregulation of multiple genes primary among them RET, which contain multiple disease associated enhancer variants. Since HSCR is a developmental disease, we used human fetal gut tissue containing both disease associated and non-associated RET enhancer variants to measure gene expression changes. We demonstrated that HSCR associated enhancer variants cumulatively decrease RET expression which consequently leads to changes in expression of thousands of genes expressed in enteric nervous system as well as surrounding gut tissue like the smooth muscle. Among these affected genes are a repertoire of transcription factors (TFs), whose knockdown leads to changes in RET transcription, highlighting bidirectional feedback between RET and multiple TFs. Our study helps to better understand how common RET enhancer variants routinely discovered in HSCR patients affects the transcriptional program during gut development potentially leading to the observed pathophysiology at birth and has implications for many congenital disorders in general and gastrointestinal motility disorders in particular. [ABSTRACT FROM AUTHOR]

Activating Transcription Factor 4 (ATF4) is an important regulator of gene expression in stress responses and developmental processes in many cell types. Here, we catalogued ATF4 binding sites in the human genome and identified overlaps with trait-associated genetic variants. We probed these genetic variants for allelic regulatory activity using a massively parallel reporter assay (MPRA) in HepG2 hepatoma cells exposed to tunicamycin to induce endoplasmic reticulum stress and ATF4 upregulation. The results revealed that in the majority of cases, the MPRA allelic activity of these SNPs was in agreement with the nucleotide preference seen in the ATF4 binding motif from ChIP-Seq. Luciferase and electrophoretic mobility shift assays in additional cellular models further confirmed ATF4-dependent regulatory effects for the SNPs rs532446 (GADD45A intronic; linked to hematological parameters), rs7011846 (LPL upstream; myocardial infarction), rs2718215 (diastolic blood pressure), rs281758 (psychiatric disorders) and rs6491544 (educational attainment). CRISPR-Cas9 disruption and/or deletion of the regulatory elements harboring rs532446 and rs7011846 led to the downregulation of GADD45A and LPL, respectively. Thus, these SNPs could represent examples of GWAS genetic variants that affect gene expression by altering ATF4-mediated transcriptional activation. Author summary: Genome-wide association studies (GWAS) have gathered a large amount of data connecting genetic variations (mostly single-nucleotide polymorphisms (SNPs)) to phenotypic traits, including diseases, but the molecular mechanisms through which the SNPs lead to phenotypic effects are often obscure. In this paper, we compiled a whole-genome catalogue of chromatin loci interacting with activating transcription factor 4 (ATF4), an important regulator of gene expression in cellular stress responses and other biological processes, and identified trait-associated SNPs overlapping ATF4 binding motifs. We then evaluated the consequences of more than 500 ATF4 binding motif-altering SNPs on allelic transcriptional activity using a massively parallel reporter assay in hepatoma cells in unstressed and stressful conditions. Further characterization of five selected SNPs in additional cellular models confirmed ATF4-dependent regulatory effects for the SNPs rs532446 (GADD45A intronic; linked to blood cell parameters), rs7011846 (LPL upstream; myocardial infarction), rs2718215 (diastolic blood pressure), rs281758 (psychiatric disorders) and rs6491544 (educational attainment). CRISPR/Cas9 technology was used to confirm the effect of rs532446 region on GADD45A and rs7011846 region on LPL gene expression. In conclusion, the data indicate that allelic variants of rs532446, rs7011846 and several other GWAS trait-linked SNPs located in ATF4 binding motifs differ in the ability to activate ATF4-mediated gene expression. [ABSTRACT FROM AUTHOR]

Imprinted genes are subject to germline epigenetic modification resulting in parental-specific allelic silencing. Although genomic imprinting is thought to be important for maternal behaviour, this idea is based on serendipitous findings from a small number of imprinted genes. Here, we undertook an unbiased systems biology approach, taking advantage of the recent delineation of specific neuronal populations responsible for controlling parental care, to test whether imprinted genes significantly converge to regulate parenting behaviour. Using single-cell RNA sequencing datasets, we identified a specific enrichment of imprinted gene expression in a recognised "parenting hub", the galanin-expressing neurons of the preoptic area. We tested the validity of linking enriched expression in these neurons to function by focusing on MAGE family member L2 (Magel2), an imprinted gene not previously linked to parenting behaviour. We confirmed expression of Magel2 in the preoptic area galanin expressing neurons. We then examined the parenting behaviour of Magel2-null(+/p) mice. Magel2-null mothers, fathers and virgin females demonstrated deficits in pup retrieval, nest building and pup-directed motivation, identifying a central role for this gene in parenting. Finally, we show that Magel2-null mothers and fathers have a significant reduction in POA galanin expressing cells, which in turn contributes to a reduced c-Fos response in the POA upon exposure to pups. Our findings identify a novel imprinted gene that impacts parenting behaviour and, moreover, demonstrates the utility of using single-cell RNA sequencing data to predict gene function from expression and in doing so here, have identified a purposeful role for genomic imprinting in mediating parental behaviour. Author summary: Genomic imprinting is a fascinating phenomenon that affects a small sub-group of the approximately 22,000 found in mammals. Unlike most genes that are equally expressed from both inherited parental copies (or alleles), so called imprinted genes are only expressed from one inherited allele, and this is usually fixed so that some imprinted genes are only active from the maternal copy, whereas others are only active from the paternal copy. This silencing of one of the parental copies makes genomic imprinting and evolutionary conundrum and the best way to understand why imprinted genes exist is to investigate the physiologies upon which they impact. Here we investigated imprinted gene expression in the brain circuitry that controls parental behaviours in mammals. We show that as a group the imprinted genes are disproportionately represented in the gene expression profile of the key neurons in this circuitry. We then tested this approach by showing that loss expression of a gene called Magel2 that was one of those imprinted genes identified in this brain circuitry, leads to deficits in parental behaviour in mice. Taken together with previous work, our findings indicate that genomic imprinting plays a particularly important role in the control of parenting behaviour. [ABSTRACT FROM AUTHOR]

RNA-RNA interactions play a crucial role in regulating gene expression and various biological processes, but identifying these interactions on a transcriptomic scale remains a challenge. To address this, we have developed a new biochemical technique called pCp-biotin labelled RNA hybrid and ultraviolet crosslinking and immunoprecipitation (lhCLIP) that enables the transcriptome-wide identification of intra- and intermolecular RNA-RNA interactions mediated by a specific RNA-binding protein (RBP). Using lhCLIP, we have uncovered a diverse landscape of intermolecular RNA interactions recognized by hnRNPK in human cells, involving all major classes of noncoding RNAs (ncRNAs) and mRNA. Notably, hnRNPK selectively binds with snRNA U4, U11, and U12, and shapes the secondary structure of these snRNAs, which may impact RNA splicing. Our study demonstrates the potential of lhCLIP as a user-friendly and widely applicable method for discovering RNA-RNA interactions mediated by a particular protein of interest and provides a valuable tool for further investigating the role of RBPs in gene expression and biological processes. Author summary: RBPs play a crucial role in post-transcriptional regulation of RNAs by modulating RNA-RNA interactions, encompassing processes like pre-mRNA splicing and mRNA degradation. Existing techniques for exploring RNA interactions associated with specific RBPs require some troublesome experimental procedures, including the use of radioisotopes to visualize RBP-RNA complexes and the introduction of exogenous expression of the target protein. In this context, we introduce lhCLIP, a novel method designed for capturing RNA-RNA interactions mediated by a specific protein. Diverging from other approaches reliant on radioisotopes, lhCLIP eliminates the need for their use. Moreover, its heightened sensitivity enables the capture of RNA-RNA interactions mediated by endogenously expressed RBPs, making it an intuitive and competitive option for researchers. We utilized lhCLIP to investigate RNA hybrids that are directly bound by hnRNPK, unveiling a complex array of intra- and intermolecular RNA interactions with the potential to exert influence over RNA splicing. [ABSTRACT FROM AUTHOR]

UPF-1-UPF-2-UPF-3 complex-orchestrated nonsense-mediated mRNA decay (NMD) is a well-characterized eukaryotic cellular surveillance mechanism that not only degrades aberrant transcripts to protect the integrity of the transcriptome but also eliminates normal transcripts to facilitate appropriate cellular responses to physiological and environmental changes. Here, we describe the multifaceted regulatory roles of the Neurospora crassa UPF complex in catalase-3 (cat-3) gene expression, which is essential for scavenging H2O2-induced oxidative stress. First, losing UPF proteins markedly slowed down the decay rate of cat-3 mRNA. Second, UPF proteins indirectly attenuated the transcriptional activity of cat-3 gene by boosting the decay of cpc-1 and ngf-1 mRNAs, which encode a well-studied transcription factor and a histone acetyltransferase, respectively. Further study showed that under oxidative stress condition, UPF proteins were degraded, followed by increased CPC-1 and NGF-1 activity, finally activating cat-3 expression to resist oxidative stress. Together, our data illustrate a sophisticated regulatory network of the cat-3 gene mediated by the UPF complex under physiological and H2O2-induced oxidative stress conditions. Author summary: The UPF complex-mediated NMD pathway is an evolutionarily conserved mRNA surveillance mechanism in eukaryotic cells, which degrades both mutated mRNAs to protect the transcriptome integrity and normal mRNAs to facilitate appropriate cellular responses. In fission yeast, mutant strains lacking UPF proteins are extremely sensitive to H2O2-induced oxidative stress. Here, we discovered a completely opposite H2O2-resistance mechanism in upf deletion mutants of the filamentous fungus N. crassa. Under physiological conditions, the N. crassa UPF complex inhibits the expression of the cat-3 gene to maintain redox homeostasis. In the presence of H2O2, UPF proteins are degraded and thus allow for the expression of the cat-3 gene, eventually ensuring adequate catalase activity to deal with oxidative stress. Our results do not only open up the post-transcriptional regulatory mechanism of catalase gene, but also identify the UPF complex as key node in response to oxidative stress. [ABSTRACT FROM AUTHOR]