Unique Transcripts Research Articles

Unveiling Key Genes and Unique Transcription Factors Involved in Secondary Cell Wall Formation in Pinus taeda.

Pinus taeda is a key timber species, and extensive research has been conducted on its wood formation. However, a comprehensive investigation into the biosynthetic pathways of lignin, cellulose, and hemicellulose in P. taeda is lacking, resulting in an incomplete understanding of secondary cell wall (SCW) formation in this species. In this study, we systematically analyzed transcriptomic data from previously published sources and constructed detailed pathways for lignin, cellulose, and hemicellulose biosynthesis. We identified 188 lignin-related genes and 78 genes associated with cellulose and hemicellulose biosynthesis. An RT-qPCR highlighted 15 key lignin biosynthesis genes and 13 crucial genes for cellulose and hemicellulose biosynthesis. A STEM analysis showed that most essential enzyme-coding genes clustered into Profile 14, suggesting their significant role in SCW formation. Additionally, we identified seven NAC and six MYB transcription factors (TFs) from atypical evolutionary clades, with distinct expression patterns from those of the previously characterized NAC and MYB genes, indicating potentially unique functions in SCW formation. This research provides the first comprehensive overview of lignin, cellulose, and hemicellulose biosynthetic genes in P. taeda and underscores the importance of non-canonical NAC and MYB TFs, laying a genetic foundation for future studies on SCW regulatory mechanisms.

Activation, regulation, and effects of the mitochondrial unfolded protein response in endothelial cells

Abstract Introduction Endoplasmic reticulum stress and the resulting unfolded protein response (UPRER) play significant roles in endothelial dysfunction and atherosclerosis. However, the understanding of the mitochondrial unfolded protein response (UPRMito) in endothelial dysfunction is limited. This study aims to characterize the activation, regulation, and effects of the UPRMito to provide a comprehensive understanding of endothelial UPRMito. Methods and Results RNA sequencing and quantitative proteomic analysis were employed to investigate the impact of stressors on mitochondrial (MitoBlock-6, CDDO) and endoplasmic (Tunicamycin) protein homeostasis on the transcriptome and proteome of Human Coronary Artery Endothelial cells (HCAEC, Fig.1). The analysis identified 12,745 unique transcripts and 7,891 unique proteins. Mitocarta annotation highlighted the regulatory effect induced by these stressors, particularly on mitochondrial proteins and transcripts (Fig. 2). Geneset enrichment analyses revealed upregulation of pathways regulated by activating transcription factors 4/5 (ATF4/5) and the integrated stress response marker CHOP. Analysis of protein expression changes unveiled that 21 proteins were upregulated by both CDDO and MitoBlock-6, including essential mitochondrial chaperones (HSP10 and HSP60). In contrast, the UPRER stressor Tunicamycin did not induce canonical UPRMito activation. Additionally, 142 proteins were commonly downregulated, affecting crucial biological processes such as mitochondrial translation, mitochondrial gene expression, and ribosome biogenesis. Importantly, NAD+ADP-ribosyltransferase activity was highly downregulated, involving PARP proteins with roles in DNA damage repair. On a functional level, mitochondrial stressors led to dose-dependent alterations in HCAEC viability, apoptosis, and reactive oxygen species formation (ROS) measured by resazurin-based assays, caspase 3/7 activity, and MitoSox staining, respectively. Low doses and short incubations with mitochondrial stressors promoted endothelial cell health through cytoprotective UPRMito pathways, while prolonged incubation led to unresolvable mitochondrial stress (exemplarily shown in Fig. 3). SiRNA-mediated gene silencing of the UPRMito regulators ATF4 and ATF5 promoted mitochondrial ROS formation in cells treated with mitochondrial stressors. These effects were particularly more pronounced after ATF5-Knockdown, suggesting that ATF5 is the main regulator of the UPRMito in HCAEC (Fig. 4). Moreover, high-dose UPRMito stressors led to an impairment in mitochondrial membrane integrity with cytosolic release of pro-inflammatory mtDNA. Conclusion UPRMito and UPRER lead to a different stress response. Depending on the underlying conditions, both stress responses can exert both cytoprotective and cytotoxic effects. Further studies are required to analyze the in vivo relevance and possible therapeutic exploitation of the UPRMito in endothelial dysfunction and atherosclerosis.Figures

Integrated temporal transcriptional and epigenetic single-cell analysis reveals the intrarenal immune characteristics in an early-stage model of IgA nephropathy during its acute injury.

Kidney inflammation plays a crucial role in the pathogenesis of IgA nephropathy (IgAN), yet the specific phenotypes of immune cells involved in disease progression remain incompletely understood. Utilizing joint profiling through longitudinal single-cell RNA-sequencing (scRNAseq) and single-cell assay for transposase-accessible chromatin sequencing (scATACseq) can provide a comprehensive framework for elucidating the development of cell subset diversity and how chromatin accessibility regulates transcription. We aimed to characterize the dynamic immune cellular landscape at a high resolution in an early IgAN mouse model with acute kidney injury (AKI). A murine model was utilized to mimic 3 immunological states -"immune stability (IS), immune activation (IA) and immune remission (IR)" in early human IgAN-associated glomerulopathy during AKI, achieved through lipopolysaccharide (LPS) injection. Urinary albumin to creatinine ratio (UACR) was measured to further validate the exacerbation and resolution of kidney inflammation during this course. Paired scRNAseq and scATACseq analysis was performed on CD45+ immune cells isolated from kidney tissues obtained from CTRL (healthy vehicle), IS, IA and IR (4 or 5 mice each). The analyses revealed 7 major cell types and 24 clusters based on 72304 single-cell transcriptomes, allowing for the identification and characterization of various immune cell types within each cluster. Our data offer an impartial depiction of the immunological characteristics, as the proportions of immune cell types fluctuated throughout different stages of the disease. Specifically, these analyses also revealed novel subpopulations, such as a macrophage subset (Nlrp1b Mac) with distinct epigenetic features and a unique transcription factor motif profile, potentially exerting immunoregulatory effects, as well as an early subset of Tex distinguished by their effector and cytolytic potential (CX3CR1-transTeff). Furthermore, in order to investigate the potential interaction between immune cells and renal resident cells, we conducted single-cell RNA sequencing on kidney cells obtained from a separate cohort of IS and IA mice without isolating immune cells. These findings underscored the diverse roles played by macrophages and CD8+ T cells in maintaining homeostasis of endothelial cells (ECs) under stress. This study presents a comprehensive analysis of the dynamic changes in immune cell profiles in a model of IgAN, identifying key cell types and their roles and interactions. These findings significantly contribute to the understanding of the pathogenesis of IgAN and may provide potential targets for therapeutic intervention.

Modulation of fungal phosphate homeostasis by the plant hormone strigolactone

Inter-kingdom communication through small molecules is essential to the coexistence of organisms in an ecosystem. In soil communities, the plant root is a nexus of interactions for a remarkable number of fungi and is a source of small-molecule plant hormones that shape fungal compositions. Although hormone signaling pathways are established in plants, how fungi perceive and respond to molecules is unclear because many plant-associated fungi are recalcitrant to experimentation. Here, we develop an approach using the model fungus, Saccharomyces cerevisiae, to elucidate mechanisms of fungal response to plant hormones. Two plant hormones, strigolactone and methyl jasmonate, produce unique transcript profiles in yeast, affecting phosphate and sugar metabolism, respectively. Genetic analysis in combination with structural studies suggests that SLs require the high-affinity transporter Pho84 to modulate phosphate homeostasis. The ability to study small-molecule plant hormones in a tractable genetic system should have utility in understanding fungal-plant interactions.

Sex specific gene expression is present prior to metamorphosis in the sea urchin

A profound collaboration between the germline and somatic cells of an organism is the creation of a functional gonad. Here we establish a foundation for studying molecular gonadogenesis in the sea urchin by use of RNA-seq, quantitative mRNA measurements, and in-situ hybridizations throughout the life cycle of the variegated sea urchin, Lytechinus variegatus (Lv). We found through three distinct analyses that the ovary and testis of this echinoderm expresses unique transcripts involved in gametogenesis, and also discovered uncharacterized gene products unique to each gonad. We further developed a pipeline integrating timepoint RNA-seq data throughout development to identify hallmark gene expression in gonads. We found that meiotic and candidate genes involved in sex determination are first expressed surprisingly early during larval growth, and well before metamorphosis. We further discovered that individual larvae express varying amounts of male- or female- hallmarks before metamorphosis, including germline, oocyte, sperm, and meiotic related genes. These distinct male- or female- gonad gene profiles may indicate the onset of, and commitment to, development of a bipotential gonad primordium, and may include metabolic differences, supported by the observation that transcripts involved in glycolysis are highly enriched in the ovary compared to the testis. Together these data support a hypothesis that sex determination is initiated prior to metamorphosis in the sea urchin and that the many uncharacterized genes unique to each gonad type characterized herein may reveal unique pathways and mechanisms in echinoderm reproduction.

Fine mapping of RNA isoform diversity using an innovative targeted long-read RNA sequencing protocol with novel dedicated bioinformatics pipeline.

Solving the structure of mRNA transcripts is a major challenge for both research and molecular diagnostic purposes. Current approaches based on short-read RNA sequencing and RT-PCR techniques cannot fully explore the complexity of transcript structure. The emergence of third-generation long-read sequencing addresses this problem by solving this sequence directly. However, genes with low expression levels are difficult to study with the whole transcriptome sequencing approach. To fix this technical limitation, we propose a novel method to capture transcripts of a gene panel using a targeted enrichment approach suitable for Pacific Biosciences and Oxford Nanopore Technologies platforms. We designed a set of probes to capture transcripts of a panel of genes involved in hereditary breast and ovarian cancer syndrome. We present SOSTAR (iSofOrmS annoTAtoR), a versatile pipeline to assemble, quantify and annotate isoforms from long read sequencing using a new tool specially designed for this application. The significant enrichment of transcripts by our capture protocol, together with the SOSTAR annotation, allowed the identification of 1,231 unique transcripts within the gene panel from the eight patients sequenced. The structure of these transcripts was annotated with a resolution of one base relative to a reference transcript. All major alternative splicing events of the BRCA1 and BRCA2 genes described in the literature were found. Complex splicing events such as pseudoexons were correctly annotated. SOSTAR enabled the identification of abnormal transcripts in the positive controls. In addition, a case of unexplained inheritance in a family with a history of breast and ovarian cancer was solved by identifying an SVA retrotransposon in intron 13 of the BRCA1 gene. We have validated a new protocol for the enrichment of transcripts of interest using probes adapted to the ONT and PacBio platforms. This protocol allows a complete description of the alternative structures of transcripts, the estimation of their expression and the identification of aberrant transcripts in a single experiment. This proof-of-concept opens new possibilities for RNA structure exploration in both research and molecular diagnostics.

Abstract B090: Evaluating the molecular determinants of PDAC racial health disparities

Abstract Pancreatic cancer is the 3rd leading cause of cancer-related deaths, with a 5-year survival rate of just 13%. Pancreatic ductal adenocarcinoma (PDAC), accounting for over 90% of pancreatic malignancies, has the highest incidence rate, mortality rate, and shortest survival times in non- Hispanic Black (defined here as African and/or African American ancestry) patients compared to other races. To identify molecular features that may contribute to racial health disparities in PDAC, we used quantitative mass spectrometry to determine the proteomic signatures unique to tumor racial origin in a cohort of 30 Caucasian and 12 African American patients. This analysis of the bulk tumor proteome revealed 183 proteins were differentially expressed between these groups. The protein with the highest expression in tumors from Black patients relative to Whites was Ring Finger Protein 2 (RNF2) that acts as an epigenetic transcriptional repressor in developmental processes pertaining to cellular differentiation. To determine the molecular pathways in PDAC regulated by RNF2, chromatin immunoprecipitation and sequencing (ChIP-seq) was performed on MIA-PaCa2 and a primary PDAC cell line derived from a Black patient to identify binding sites for RNF2 in PDAC cells. From the 1,142 replicated peaks identified, RNF2 peaks were localized over the transcription start site of 644 unique transcripts, including GATA6, which is an important determinant in PDAC subtype delineation with higher GATA6 expression associated with the less aggressive and more chemosensitive Classical subtype. Modulating RNF2 levels/activity by exogenous expression, short hairpin RNA (shRNA) knockdown, and small molecule inhibition confirms that GATA6 is a target of RNF2 repression in PDAC. Multi-omic (transcriptomic, proteomic, and phosphoproteomic) characterization was performed on MIA-PaCa2 cells with CRISPR-Cas9 knockout (KO) of RNF2. Differential gene expression (DE) analysis revealed changes in platinum drug resistance, EMT, and innate immune response, which are characteristic of the aggressive Inflammatory PDAC subtype we previously reported using proteomics. Total and phospho-proteomics corroborated these results identifying changes in DNA repair, inflammation response, TNFα mediated signaling, and apoptotic signaling. These results reflect the enrichment of Inflammatory subtype markers in tumors originating in Black patients. In agreement with these results, RNF2 KO shows a significant reduction in migration, as well as an increase in sensitivity to oxaliplatin. Our work identifies differences in the underlying proteome of PDAC tumors associated with tumor racial origin and describes how elevated expression of RNF2 in PDACs from Black patients contribute to inflammation and subtype delineation. This study delineates the associations between tumor racial origin, subtype specification, and response to therapy that could mitigate PDAC health disparities and be broadly applied across all PDAC patients. Citation Format: Anthony E Johansen-Sallee, Alexa M Barber, Henry C.-H. Law, Vignesh Vudatha, Jose Trevino, Nicholas T Woods. Evaluating the molecular determinants of PDAC racial health disparities [abstract]. In: Proceedings of the AACR Special Conference in Cancer Research: Advances in Pancreatic Cancer Research; 2024 Sep 15-18; Boston, MA. Philadelphia (PA): AACR; Cancer Res 2024;84(17 Suppl_2):Abstract nr B090.

Establishment of single-cell transcriptional states during seed germination

Germination involves highly dynamic transcriptional programs as the cells of seeds reactivate and express the functions necessary for establishment in the environment. Individual cell types have distinct roles within the embryo, so must therefore have cell type-specific gene expression and gene regulatory networks. We can better understand how the functions of different cell types are established and contribute to the embryo by determining how cell type-specific transcription begins and changes through germination. Here we describe a temporal analysis of the germinating Arabidopsis thaliana embryo at single-cell resolution. We define the highly dynamic cell type-specific patterns of gene expression and how these relate to changing cellular function as germination progresses. Underlying these are unique gene regulatory networks and transcription factor activity. We unexpectedly discover that most embryo cells transition through the same initial transcriptional state early in germination, even though cell identity has already been established during embryogenesis. Cells later transition to cell type-specific gene expression patterns. Furthermore, our analyses support previous findings that the earliest events leading to the induction of seed germination take place in the vasculature. Overall, our study constitutes a general framework with which to characterize Arabidopsis cell transcriptional states through seed germination, allowing investigation of different genotypes and other plant species whose seed strategies may differ.

The extensive m5C epitranscriptome of Thermococcus kodakarensis is generated by a suite of RNA methyltransferases that support thermophily

RNAs are often modified to invoke new activities. While many modifications are limited in frequency, restricted to non-coding RNAs, or present only in select organisms, 5-methylcytidine (m5C) is abundant across diverse RNAs and fitness-relevant across Domains of life, but the synthesis and impacts of m5C have yet to be fully investigated. Here, we map m5C in the model hyperthermophile, Thermococcus kodakarensis. We demonstrate that m5C is ~25x more abundant in T. kodakarensis than human cells, and the m5C epitranscriptome includes ~10% of unique transcripts. T. kodakarensis rRNAs harbor tenfold more m5C compared to Eukarya or Bacteria. We identify at least five RNA m5C methyltransferases (R5CMTs), and strains deleted for individual R5CMTs lack site-specific m5C modifications that limit hyperthermophilic growth. We show that m5C is likely generated through partial redundancy in target sites among R5CMTs. The complexity of the m5C epitranscriptome in T. kodakarensis argues that m5C supports life in the extremes.

Unravelling transcriptional responses of the willow to Fusarium kuroshium infection

Fusarium kuroshium and the ambrosia beetle Euwallacea kuroshio attack several tree species, including the arroyo willow (Salix lasiolepis). Despite the importance of this pest, there is no molecular information about this plant-fungus interaction. Using an in vitro pathosystem of S. lasiolepis leaves – F. kurosihum and the RNA-seq approach, this study analyzed for the first time the transcriptome of both organisms at 4-, 8-, and 12-days post-inoculation. De novo assembly was done to generate the willow transcriptome, and a total of 361,2151 unique transcripts (UniGenes) were found. 4,583 UniGenes responded to the infection, and the posterior analysis showed that it contains the molecular elements to develop both pattern- and effector-triggered immunity responses (PTI and ETI, respectively). Significant expression changes of genes involved in an oxidative burst, Ca2+ signaling, and mutually antagonistic interaction between salicylic acid- and jasmonic acid/ethylene-dependent signaling suggest that the willow recognized the fungus and triggered downstream responses. Additionally, chromatin remodelers and non-coding RNA play a significant role in plant defense. Concerning F. kuroshium, 178 induced genes that show homology with genes involved in the pathogenicity in other ascomycetes fungi plus transport proteins, ribosomal proteins, and uncharacterized proteins with specific domains were identified. Collectively these data provide a basis for future research on the molecular response of this transcendental tree during biotic stress, specifically during Fusarium dieback disease.

Granulosa cells provide transcriptomic information on ovarian follicle dynamics in southern white rhinoceros

Much remains unknown about the reproductive physiology of southern white rhinoceros (SWR) and the effect of ovarian stimulation prior to ovum pickup (OPU) have not been fully elucidated. Granulosa cells (GC) provide valuable insight into follicle growth and oocyte maturation status. The goals of this study were to evaluate transcriptomic changes in GC from three stages of follicle development and to identify biomarkers possibly associated with follicular growth and maturation as a result of ovarian stimulation. GC collected from SWRs following OPU were assigned stages based upon follicle size. Total RNA was isolated, and cDNA libraries were prepared and sequenced on a NovaSeq 6000. All bioinformatics analyses were performed utilizing the Galaxy web platform. Reads were aligned to CerSimCot1.0, and the manual curation was performed with EquCab3.0. Overall, 39,455 transcripts (21,612 genes) were identified across follicle stages, and manual curation yielded a 61% increase in gene identification from the original annotation. Granulosa cells from preovulatory follicles expressed the highest number of unique transcripts. The following seven biomarkers were determined based upon cluster analysis and patterns of expression: COL1A1, JMY, FBXW11, NRG1, TMPO, MACIR and COL4A1. These data can be used to potentially evaluate the effects of different ovarian stimulation protocols on follicle dynamics, improve OPU results, and support conservation efforts in this species.

Distinct granzyme k expression in immune cells: a single-cell rna-seq meta-analysis.

Granzymes are essential serine proteases in cytotoxic T cells and natural killer (NK) cells, with GZMK's expression being less understood. This study aims to uncover GZMK expression profiles across various immune cell types using single-cell RNA sequencing meta-analysis. This study aims to uncover GZMK expression profiles across various immune cell types using single-cell RNA sequencing meta-analysis. We conducted a meta-analysis using cellxgene, an interactive data exploration platform developed by the Chan Zuckerberg Initiative. We focused on mature T cells, NK cells, B cells, and NKT cells. We also checked transcription factor binding sites at the granzyme gene promoter regions using JASPAR. Comparative analysis was also done using mouse single-cell RNA sequencing data. GZMK was the most lowly expressed in NK cells and mature NKT cells in most tissues except for colon and lymph nodes. In mature T cells, GZMK is similarly or more highly expressed than other granzymes. HBCA data revealed weak expression of GZMK in NK cells but strong expression in effector memory CD8-positive, alpha-beta T cells. Combined data shows no significant difference in GZMK expression between cell types. Subtype analysis shows that GZMK expression was higher in CD16-negative, CD56-bright NK cells when compared to CD16-positive, CD56-dim NK cells. We also identified unique transcription factor binding sites for GZMK. While this pattern in mouse data with low Gzmk expression in NK cells and higher T cells was repeated. GZMK expression is distinctively regulated among immune cells and tissues, with unique promoter regions and transcription factor binding sites contributing to this differential expression. These insights into GZMK's role in immune function and regulation offer potential therapeutic targets.

Characterizing the splice map of Turkey Hemorrhagic Enteritis Virus

BackgroundHemorrhagic enteritis, caused by Turkey Hemorrhagic Enteritis Virus (THEV), is a disease affecting turkey poults characterized by immunosuppression and bloody diarrhea. An avirulent THEV strain that retains the immunosuppressive ability is used as a live vaccine. Characterizing the splice map of THEV is an essential step that would allow studies of individual genes mediating its immunosuppressive functions. We used RNA sequencing to characterize the splice map of THEV for the first time, providing key insights into the THEV gene expression and mRNA structures.MethodsAfter infecting a turkey B-cell line with the vaccine strain, samples in triplicates were collected at 4-, 12-, 24-, and 72-hours post-infection. Total RNA was extracted, and poly-A-tailed mRNA sequenced. Reads were mapped to the THEV genome after trimming and transcripts assembled with StringTie. We performed PCR of THEV cDNA, cloned the PCR products, and used Sanger sequencing to validate all identified splice junctions.ResultsResearchers previously annotated the THEV genome as encoding 23 open reading frames (ORFs). We identified 29 spliced transcripts from our RNA sequencing data, all containing novel exons although some exons matched some previously annotated ORFs. The three annotated splice junctions were also corroborated by our data. During validation we identified five additional unique transcripts, a subset of which were further validated by 3’ rapid amplification of cDNA ends (3’ RACE). Thus, we report that the genome of THEV contains 34 transcripts with the coding capacity for all annotated ORFs. However, we found six of the previously annotated ORFs to be truncated ORFs on the basis of the identification of an in-frame upstream start codon or the detection of additional coding exons. We also identified three of the annotated ORFs with longer or shorter isoforms, and seven novel unannotated ORFs that could potentially be translated; although it is beyond the scope of this manuscript to investigate whether they are translated.ConclusionsSimilar to human adenoviruses, all THEV transcripts are spliced and organized into five transcription units under the control of their cognate promoters. The genes are expressed under temporal regulation and THEV also produces multiple distinctly spliced transcripts that code for the same protein. Studies of the newly identified potential proteins should be urgently performed as these proteins may have roles in THEV-induced immunosuppression. Also, knowing the splicing of THEV genes should be invaluable to future research focusing on studying THEV genes, as this will allow accurate cloning of the mRNAs.

Proliferative arrest induces neuron differentiation and innate immune responses in control and Creutzfeldt-Jakob Disease agent infected rat septal neurons.

Rat post-mitotic septal (SEP) neurons, engineered to conditionally proliferate at 33°C, differentiate when arrested at 37.5°C and can be maintained for weeks without cytotoxic effects. Nine independent cDNA libraries were made to follow arrest-induced neural differentiation and innate immune responses in normal (Nl) uninfected and CJ agent infected SEP cells. Proliferating Nl versus latently infected (CJ-) cells showed few RNA-seq differences. However arrest induced major changes. Normal cells displayed a plethora of anti-proliferative transcripts. Additionally, known neuron differentiation transcripts, e.g., Agtr2, Neuregulin-1, GDF6, SFRP4 and Prnp were upregulated. These Nl neurons also displayed many activated IFN innate immune genes, e.g., OAS1, RTP4, ISG20, GTB4, CD80 and cytokines, complement, and clusterin (CLU) that binds to misfolded proteins. In contrast, arrested highly infectious CJ+ cells (10 logs/gm) downregulated many replication controls. Furthermore, arrested CJ+ cells suppressed neuronal differentiation transcripts, including Prnp which is essential for CJ agent infection. CJ+ cells also enhanced IFN stimulated pathways, and analysis of the 342 CJ+ unique transcripts revealed additional innate immune and anti-viral-linked transcripts, e.g., Il17, ISG15, and RSAD2 (viperin). These data show: 1) innate immune transcripts are produced by normal neurons during differentiation; 2) CJ infection can enhance and expand anti-viral responses; 3) latent CJ infection epigenetically imprints many proliferative pathways to thwart complete arrest. CJ+ brain microglia, white blood cells and intestinal myeloid cells with shared transcripts may be stimulated to educe latent CJD infections that can be clinically silent for >30 years.

Alternative TSS use is widespread in Cryptococcus fungi in response to environmental cues and regulated genome-wide by the transcription factor Tur1.

Alternative transcription start site (TSS) usage regulation has been identified as a major means of gene expression regulation in metazoans. However, in fungi, its impact remains elusive as its study has thus far been restricted to model yeasts. Here, we first re-analyzed TSS-seq data to define genuine TSS clusters in 2 species of pathogenic Cryptococcus. We identified 2 types of TSS clusters associated with specific DNA sequence motifs. Our analysis also revealed that alternative TSS usage regulation in response to environmental cues is widespread in Cryptococcus, altering gene expression and protein targeting. Importantly, we performed a forward genetic screen to identify a unique transcription factor (TF) named Tur1, which regulates alternative TSS (altTSS) usage genome-wide when cells switch from exponential phase to stationary phase. ChiP-Seq and DamID-Seq analyses suggest that at some loci, the role of Tur1 might be direct. Tur1 has been previously shown to be essential for virulence in C. neoformans. We demonstrated here that a tur1Δ mutant strain is more sensitive to superoxide stress and phagocytosed more efficiently by macrophages than the wild-type (WT) strain.

Cellular zinc status alters chromatin accessibility and binding of p53 to DNA.

Zn2+ is an essential metal required by approximately 850 human transcription factors. How these proteins acquire their essential Zn2+ cofactor and whether they are sensitive to changes in the labile Zn2+ pool in cells remain open questions. Using ATAC-seq to profile regions of accessible chromatin coupled with transcription factor enrichment analysis, we examined how increases and decreases in the labile zinc pool affect chromatin accessibility and transcription factor enrichment. We found 685 transcription factor motifs were differentially enriched, corresponding to 507 unique transcription factors. The pattern of perturbation and the types of transcription factors were notably different at promoters versus intergenic regions, with zinc-finger transcription factors strongly enriched in intergenic regions in elevated Zn2+ To test whether ATAC-seq and transcription factor enrichment analysis predictions correlate with changes in transcription factor binding, we used ChIP-qPCR to profile six p53 binding sites. We found that for five of the six targets, p53 binding correlates with the local accessibility determined by ATAC-seq. These results demonstrate that changes in labile zinc alter chromatin accessibility and transcription factor binding to DNA.

High-density resolution of the Kaposi's sarcoma associated herpesvirus transcriptome identifies novel transcript isoforms generated by long-range transcription and alternative splicing.

Kaposi's sarcoma-associated herpesvirus is the etiologic agent of Kaposi's sarcoma and two B-cell malignancies. Recent advancements in sequencing technologies have led to high resolution transcriptomes for several human herpesviruses that densely encode genes on both strands. However, for KSHV progress remained limited due to the overall low percentage of KSHV transcripts, even during lytic replication. To address this challenge, we have developed a target enrichment method to increase the KSHV-specific reads for both short- and long-read sequencing platforms. Furthermore, we combined this approach with the Transcriptome Resolution through Integration of Multi-platform Data (TRIMD) pipeline developed previously to annotate transcript structures. TRIMD first builds a scaffold based on long-read sequencing and validates each transcript feature with supporting evidence from Illumina RNA-Seq and deepCAGE sequencing data. Our stringent innovative approach identified 994 unique KSHV transcripts, thus providing the first high-density KSHV lytic transcriptome. We describe a plethora of novel coding and non-coding KSHV transcript isoforms with alternative untranslated regions, splice junctions and open-reading frames, thus providing deeper insights on gene expression regulation of KSHV. Interestingly, as described for Epstein-Barr virus, we identified transcription start sites that augment long-range transcription and may increase the number of latency-associated genes potentially expressed in KS tumors.

Integrating Iso-seq and RNA-seq data for the reannotation of the greater amberjack genome

The greater amberjack is a very important fishery species with high commercial value, and it is distributed worldwide. Transcriptome-based studies on S. dumerili have been limited by an inadequate reference genome and a lack of well-annotated full-length transcripts. In this study, a total of 12 tissues from juvenile and adult fish both sexes were collected for next-generation RNA sequencing (RNA-seq) and full-length isoform sequencing (Iso-seq). For Iso-seq, a total of 163,218, 149,716, and 189,169 high-quality unique transcript sequences were obtained, with an N50 of 5,441, 5,255, and 5,939, from juvenile, adult male and adult female S. dumerili, respectively. We integrated the Iso-seq and RNA-seq data to construct a comprehensive gene annotation and systematically profiled the dynamics of gene expression across the 12 tissues. Our gene models had greater detail and accuracy than those from NCBI and Ensembl, with more precise polyA locations. These resources serve as a foundation for functional genomic studies and provide valuable insights into the molecular mechanisms underlying the development, reproduction and commercial traits of amberjack.

Cryo-EM structures of Thogoto virus polymerase reveal unique RNA transcription and replication mechanisms among orthomyxoviruses

Influenza viruses and thogotoviruses account for most recognized orthomyxoviruses. Thogotoviruses, exemplified by Thogoto virus (THOV), are capable of infecting humans using ticks as vectors. THOV transcribes mRNA without the extraneous 5′ end sequences derived from cap-snatching in influenza virus mRNA. Here, we report cryo-EM structures to characterize THOV polymerase RNA synthesis initiation and elongation. The structures demonstrate that THOV RNA transcription and replication are able to start with short dinucleotide primers and that the polymerase cap-snatching machinery is likely non-functional. Triggered by RNA synthesis, asymmetric THOV polymerase dimers can form without the involvement of host factors. We confirm that, distinctive from influenza viruses, THOV-polymerase RNA synthesis is weakly dependent of the host factors ANP32A/B/E in human cells. This study demonstrates varied mechanisms in RNA synthesis and host factor utilization among orthomyxoviruses, providing insights into the mechanisms behind thogotoviruses’ broad-infectivity range.

Deciphering novel TCF4-driven mechanisms underlying a common triplet repeat expansion-mediated disease.

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.