Optimizing RNA-Seq Mapping with STAR

Systematic Characterization of p53-Regulated Long Noncoding RNAs across Human Cancers Reveals Remarkable Heterogeneity among Different Tumor Types

The p53 tumor suppressor protein, a sequence-specific DNA binding transcription factor, regulates the expression of a large number of genes, in response to various forms of cellular stress. Although the protein coding target genes of p53 have been well studied, less is known about its role in regulating long noncoding genes and their functional relevance to cancer. Here we report the genome-wide identification of a large set (>1,000) of long noncoding RNAs (lncRNA), which are putative p53 targets in a colon cancer cell line and in human patient datasets from five different common types of cancer. These lncRNAs have not been annotated by other studies of normal unstressed systems. In the colon cancer cell line, a high proportion of these lncRNAs are uniquely induced by different chemotherapeutic agents that activate p53, whereas others are induced by more than one agent tested. Further, subsets of these lncRNAs independently predict overall and disease-free survival of patients across the five different common cancer types. Interestingly, both genetic alterations and patient survival associated with different lncRNAs are unique to each cancer tested, indicating extraordinary tissue-specific variability in the p53 noncoding response. The newly identified noncoding p53 target genes have allowed us to construct a classifier for tumor diagnosis and prognosis.

IMPLICATIONS

Our results not only identify myriad p53-regulated long noncoding (lncRNA), they also reveal marked drug-induced, as well as tissue- and tumor-specific heterogeneity in these putative p53 targets and our findings have enabled the construction of robust classifiers for diagnosis and prognosis.

AR coactivators, CBP/p300, are critical mediators of DNA repair in prostate cancer

Castration resistant prostate cancer (CRPC) remains an incurable disease stage with ineffective treatments options. Here, the androgen receptor (AR) coactivators CBP/p300, which are histone acetyltransferases, were identified as critical mediators of DNA damage repair (DDR) to potentially enhance therapeutic targeting of CRPC. Key findings demonstrate that CBP/p300 expression increases with disease progression and selects for poor prognosis in metastatic disease. CBP/p300 bromodomain inhibition enhances response to standard of care therapeutics. Functional studies, CBP/p300 cistrome mapping, and transcriptome in CRPC revealed that CBP/p300 regulates DDR. Further mechanistic investigation showed that CBP/p300 attenuation via therapeutic targeting and genomic knockdown decreases homologous recombination (HR) factors in vitro, in vivo, and in human prostate cancer (PCa) tumors ex vivo. Similarly, CBP/p300 expression in human prostate tissue correlates with HR factors. Lastly, targeting CBP/p300 impacts HR-mediate repair and patient outcome. Collectively, these studies identify CBP/p300 as drivers of PCa tumorigenesis and lay the groundwork to optimize therapeutic strategies for advanced PCa via CBP/p300 inhibition, potentially in combination with AR-directed and DDR therapies.

The Long Non-Coding RNA Gene AC027288.3 Plays a Role in Human Endometrial Stromal Fibroblast Decidualization

During the secretory phase of the menstrual cycle, endometrial fibroblast cells begin to change into large epithelial-like cells called decidual cells in a process called decidualization. This differentiation continues more broadly in the endometrium and forms the decidual tissue during early pregnancy. The cells undergoing decidualization as well as the resulting decidual cells, support successful implantation and placentation during early pregnancy. This study was carried out to identify new potentially important long non-coding RNA (lncRNA) genes that may play a role in human endometrial stromal fibroblast cells (hESF) undergoing decidualization in vitro, and several were found. The expression of nine was further characterized. One of these, AC027288.3, showed a dramatic increase in the expression of hESF cells undergoing decidualization. When AC027288.3 expression was targeted, the ability of the cells to undergo decidualization as determined by the expression of decidualization marker protein-coding genes was significantly altered. The most affected markers of decidualization whose expression was significantly reduced were FOXO1, FZD4, and INHBA. Therefore, AC027288.3 may be a major upstream regulator of the WNT-FOXO1 pathway and activin-SMAD3 pathways previously shown as critical for hESF decidualization. Finally, we explored possible regulators of AC027288.3 expression during human ESF decidualization. Expression was regulated by cAMP and progesterone. Our results suggest that AC027288.3 plays a role in hESF decidualization and identifies several other lncRNA genes that may also play a role.

The repertoire and structure of adhesion GPCR transcript variants assembled from publicly available deep-sequenced human samples

Alternative splicing and multiple transcription start and termination sites can produce a diverse repertoire of mRNA transcript variants from a given gene. While the full picture of the human transcriptome is still incomplete, publicly available RNA datasets have enabled the assembly of transcripts. Using publicly available deep sequencing data from 927 human samples across 48 tissues, we quantified known and new transcript variants, provide an interactive, browser-based application Splice-O-Mat and demonstrate its relevance using adhesion G protein-coupled receptors (aGPCRs) as an example. On average, 24 different transcript variants were detected for each of the 33 human aGPCR genes, and several dominant transcript variants were not yet annotated. Variable transcription starts and complex exon-intron structures encode a flexible protein domain architecture of the N- and C termini and the seven-transmembrane helix domain (7TMD). Notably, we discovered the first GPCR (ADGRG7/GPR128) with eight transmembrane helices. Both the N- and C terminus of this aGPCR were intracellularly oriented, anchoring the N terminus in the plasma membrane. Moreover, the assessment of tissue-specific transcript variants, also for other gene classes, in our application may change the evaluation of disease-causing mutations, as their position in different transcript variants may explain tissue-specific phenotypes.

Oncogenic ETS fusions promote DNA damage and proinflammatory responses via pericentromeric RNAs in extracellular vesicles

Aberrant expression of the E26 transformation-specific (ETS) transcription factors characterizes numerous human malignancies. Many of these proteins, including EWS:FLI1 and EWS:ERG fusions in Ewing sarcoma (EwS) and TMPRSS2:ERG in prostate cancer (PCa), drive oncogenic programs via binding to GGAA repeats. We report here that both EWS:FLI1 and ERG bind and transcriptionally activate GGAA-rich pericentromeric heterochromatin. The respective pathogen-like HSAT2 and HSAT3 RNAs, together with LINE, SINE, ERV, and other repeat transcripts, are expressed in EwS and PCa tumors, secreted in extracellular vesicles (EVs), and are highly elevated in plasma of patients with EwS with metastatic disease. High human satellite 2 and 3 (HSAT2,3) levels in EWS:FLI1- or ERG-expressing cells and tumors were associated with induction of G2/M checkpoint, mitotic spindle, and DNA damage programs. These programs were also activated in EwS EV-treated fibroblasts, coincident with accumulation of HSAT2,3 RNAs, proinflammatory responses, mitotic defects, and senescence. Mechanistically, HSAT2,3-enriched cancer EVs induced cGAS-TBK1 innate immune signaling and formation of cytosolic granules positive for double-strand RNAs, RNA-DNA, and cGAS. Hence, aberrantly expressed ETS proteins derepress pericentromeric heterochromatin, yielding pathogenic RNAs that transmit genotoxic stress and inflammation to local and distant sites. Monitoring HSAT2,3 plasma levels and preventing their dissemination may thus improve therapeutic strategies and blood-based diagnostics.

Machine learning on alignment features for parent-of-origin classification of simulated hybrid RNA-seq

BACKGROUND

Parent-of-origin allele-specific gene expression (ASE) can be detected in interspecies hybrids by virtue of RNA sequence variants between the parental haplotypes. ASE is detectable by differential expression analysis (DEA) applied to the counts of RNA-seq read pairs aligned to parental references, but aligners do not always choose the correct parental reference.

RESULTS

We used public data for species that are known to hybridize. We measured our ability to assign RNA-seq read pairs to their proper transcriptome or genome references. We tested software packages that assign each read pair to a reference position and found that they often favored the incorrect species reference. To address this problem, we introduce a post process that extracts alignment features and trains a random forest classifier to choose the better alignment. On each simulated hybrid dataset tested, our machine-learning post-processor achieved higher accuracy than the aligner by itself at choosing the correct parent-of-origin per RNA-seq read pair.

CONCLUSIONS

For the parent-of-origin classification of RNA-seq, machine learning can improve the accuracy of alignment-based methods. This approach could be useful for enhancing ASE detection in interspecies hybrids, though RNA-seq from real hybrids may present challenges not captured by our simulations. We believe this is the first application of machine learning to this problem domain.

Rewiring of RNA methylation by the oncometabolite fumarate in renal cell carcinoma

Metabolic reprogramming is a hallmark of cancer that facilitates changes in many adaptive biological processes. Mutations in the tricarboxylic acid cycle enzyme fumarate hydratase (FH) lead to fumarate accumulation and cause hereditary leiomyomatosis and renal cell cancer (HLRCC). HLRCC is a rare, inherited disease characterized by the development of non-cancerous smooth muscle tumors of the uterus and skin, and an increased risk of an aggressive form of kidney cancer. Fumarate has been shown to inhibit 2-oxoglutarate-dependent dioxygenases (2OGDDs) involved in the hydroxylation of HIF1α, as well as in DNA and histone demethylation. However, the link between fumarate accumulation and changes in RNA post-transcriptional modifications has not been defined. Here, we determine the consequences of fumarate accumulation on the activity of different members of the 2OGDD family targeting RNA modifications. By evaluating multiple RNA modifications in patient-derived HLRCC cell lines, we show that mutation of FH selectively affects the levels of N6-methyladenosine (m6A), while the levels of 5-formylcytosine (f5C) in mitochondrial tRNA are unaffected. This supports the hypothesis of a differential impact of fumarate accumulation on distinct RNA demethylases. The observation that metabolites modulate specific subsets of RNA-modifying enzymes offers new insights into the intersection between metabolism and the epitranscriptome.

Evidence for widespread cytoplasmic structuring into mesoscale condensates

Compartmentalization is an essential feature of eukaryotic life and is achieved both via membrane-bound organelles, such as mitochondria, and membrane-less biomolecular condensates, such as the nucleolus. Known biomolecular condensates typically exhibit liquid-like properties and are visualized by microscopy on the scale of ~1 µm (refs. 1,2). They have been studied mostly by microscopy, examining select individual proteins. So far, several dozen biomolecular condensates are known, serving a multitude of functions, for example, in the regulation of transcription3, RNA processing4 or signalling5,6, and their malfunction can cause diseases7,8. However, it remains unclear to what extent biomolecular condensates are utilized in cellular organization and at what length scale they typically form. Here we examine native cytoplasm from Xenopus egg extract on a global scale with quantitative proteomics, filtration, size exclusion and dilution experiments. These assays reveal that at least 18% of the proteome is organized into mesoscale biomolecular condensates at the scale of ~100 nm and appear to be stabilized by RNA or gelation. We confirmed mesoscale sizes via imaging below the diffraction limit by investigating protein permeation into porous substrates with defined pore sizes. Our results show that eukaryotic cytoplasm organizes extensively via biomolecular condensates, but at surprisingly short length scales.

Use of a Candida albicans SC5314 PacBio HiFi reads dataset to close gaps in the reference genome assembly, reveal a subtelomeric gene family, and produce accurate phased allelic sequences

Candida albicans SC5314 is the most-often used strain for molecular manipulation of the species. The SC5314 reference genome sequence is the result of considerable effort from many scientists and has advanced research into fungal biology and pathogenesis. Although the resource is highly developed and presented in a phased diploid format, the sequence includes gaps and does not extend to the telomeres on its eight chromosome pairs. Accurate SC5314 genome assembly is complicated by the presence of extensive repeated sequences and considerable allelic length variation at some loci. Advances in genome sequencing technology provide the tools to obtain highly accurate long-read data that span even the most-difficult-to-assemble genome regions. Here, we describe derivation of a PacBio HiFi data set and creation of a collapsed haploid telomere-to-telomere assembly of the SC5314 genome (ASM3268872v1) that revealed previously unknown features of the strain. ASM3268872v1 subtelomeric distances were up to 19 kb larger than in the reference genome and revealed a family of highly conserved DNA helicase-encoding genes at 10 of the 16 chromosome ends. We also describe alignments of individual HiFi reads to deduce accurate diploid sequences for the most notoriously difficult-to-assemble C. albicans genes: the agglutinin-like sequence (ALS) gene family. We provide a tutorial that demonstrates how the HiFi reads can be visualized to explore any region of interest. Availability of the HiFi reads data set and the ASM3268872v1 comparative guide assembly will streamline research efforts because accurate diploid sequences can be derived using simple in silico methods rather than time-consuming laboratory-bench approaches.

Cis-regulatory modes of Ultrabithorax inactivation in butterfly forewings

Hox gene clusters encode transcription factors that drive regional specialization during animal development: for example the Hox factor Ubx is expressed in the insect metathoracic (T3) wing appendages and differentiates them from T2 mesothoracic identities. Hox transcriptional regulation requires silencing activities that prevent spurious activation and regulatory crosstalks in the wrong tissues, but this has seldom been studied in insects other than Drosophila, which shows a derived Hox dislocation into two genomic clusters that disjoined Antennapedia (Antp) and Ultrabithorax (Ubx). Here, we investigated how Ubx is restricted to the hindwing in butterflies, amidst a contiguous Hox cluster. By analysing Hi-C and ATAC-seq data in the butterfly Junonia coenia, we show that a Topologically Associated Domain (TAD) maintains a hindwing-enriched profile of chromatin opening around Ubx. This TAD is bordered by a Boundary Element (BE) that separates it from a region of joined wing activity around the Antp locus. CRISPR mutational perturbation of this BE releases ectopic Ubx expression in forewings, inducing homeotic clones with hindwing identities. Further mutational interrogation of two non-coding RNA encoding regions and one putative cis-regulatory module within the Ubx TAD cause rare homeotic transformations in both directions, indicating the presence of both activating and repressing chromatin features. We also describe a series of spontaneous forewing homeotic phenotypes obtained in Heliconius butterflies, and discuss their possible mutational basis. By leveraging the extensive wing specialization found in butterflies, our initial exploration of Ubx regulation demonstrates the existence of silencing and insulating sequences that prevent its spurious expression in forewings.

Increased iron uptake by splenic hematopoietic stem cells promotes TET2-dependent erythroid regeneration

Hematopoietic stem cells (HSCs) are capable of regenerating the blood system, but the instructive cues that direct HSCs to regenerate particular lineages lost to the injury remain elusive. Here, we show that iron is increasingly taken up by HSCs during anemia and induces erythroid gene expression and regeneration in a Tet2-dependent manner. Lineage tracing of HSCs reveals that HSCs respond to hemolytic anemia by increasing erythroid output. The number of HSCs in the spleen, but not bone marrow, increases upon anemia and these HSCs exhibit enhanced proliferation, erythroid differentiation, iron uptake, and TET2 protein expression. Increased iron in HSCs promotes DNA demethylation and expression of erythroid genes. Suppressing iron uptake or TET2 expression impairs erythroid genes expression and erythroid differentiation of HSCs; iron supplementation, however, augments these processes. These results establish that the physiological level of iron taken up by HSCs has an instructive role in promoting erythroid-biased differentiation of HSCs.

Compartments in medulloblastoma with extensive nodularity are connected through differentiation along the granular precursor lineage

Medulloblastomas with extensive nodularity are cerebellar tumors characterized by two distinct compartments and variable disease progression. The mechanisms governing the balance between proliferation and differentiation in MBEN remain poorly understood. Here, we employ a multi-modal single cell transcriptome analysis to dissect this process. In the internodular compartment, we identify proliferating cerebellar granular neuronal precursor-like malignant cells, along with stromal, vascular, and immune cells. In contrast, the nodular compartment comprises postmitotic, neuronally differentiated malignant cells. Both compartments are connected through an intermediate cell stage resembling actively migrating CGNPs. Notably, we also discover astrocytic-like malignant cells, found in proximity to migrating and differentiated cells at the transition zone between the two compartments. Our study sheds light on the spatial tissue organization and its link to the developmental trajectory, resulting in a more benign tumor phenotype. This integrative approach holds promise to explore intercompartmental interactions in other cancers with varying histology.

Efficacy and safety of a SOD1-targeting artificial miRNA delivered by AAV9 in mice are impacted by miRNA scaffold selection

Toxic gain-of-function mutations in superoxide dismutase 1 (SOD1) contribute to approximately 2%-3% of all amyotrophic lateral sclerosis (ALS) cases. Artificial microRNAs (amiRs) delivered by adeno-associated virus (AAV) have been proposed as a potential treatment option to silence SOD1 expression and mitigate disease progression. Primary microRNA (pri-miRNA) scaffolds are used in amiRs to shuttle a hairpin RNA into the endogenous miRNA pathway, but it is unclear whether different primary miRNA (pri-miRNA) scaffolds impact the potency and safety profile of the expressed amiR in vivo. In our process to develop an AAV amiR targeting SOD1, we performed a preclinical characterization of two pri-miRNA scaffolds, miR155 and miR30a, sharing the same guide strand sequence. We report that, while the miR155-based vector, compared with the miR30a-based vector, leads to a higher level of the amiR and more robust suppression of SOD1 in vitro and in vivo, it also presents significantly greater risks for CNS-related toxicities in vivo. Despite miR30a-based vector showing relatively lower potency, it can significantly delay the development of ALS-like phenotypes in SOD1-G93A mice and increase survival in a dose-dependent manner. These data highlight the importance of scaffold selection in the pursuit of highly efficacious and safe amiRs for RNA interference gene therapy.

Cellular senescence promotes progenitor cell expansion during axolotl limb regeneration

Axolotl limb regeneration is accompanied by the transient induction of cellular senescence within the blastema, the structure that nucleates regeneration. The precise role of this blastemal senescent cell (bSC) population, however, remains unknown. Here, through a combination of gain- and loss-of-function assays, we elucidate the functions and molecular features of cellular senescence in vivo. We demonstrate that cellular senescence plays a positive role during axolotl regeneration by creating a pro-proliferative niche that supports progenitor cell expansion and blastema outgrowth. Senescent cells impact their microenvironment via Wnt pathway modulation. Further, we identify a link between Wnt signaling and senescence induction and propose that bSC-derived Wnt signals facilitate the proliferation of neighboring cells in part by preventing their induction into senescence. This work defines the roles of cellular senescence in the regeneration of complex structures.

Transcriptome-based identification of tumor-reactive and bystander CD8+ T cell receptor clonotypes in human pancreatic cancer

Pancreatic ductal adenocarcinoma (PDAC) is generally refractory to immune checkpoint blockade, although patients with genetically unstable tumors can show modest therapeutic benefit. We previously demonstrated the presence of tumor-reactive CD8+ T cells in PDAC samples. Here, we charted the tumor-infiltrating T cell repertoire in PDAC by combining single-cell transcriptomics with functional testing of T cell receptors (TCRs) for reactivity against autologous tumor cells. On the basis of a comprehensive dataset including 93 tumor-reactive and 65 bystander TCR clonotypes, we delineated a gene signature that effectively distinguishes between these T cell subsets in PDAC, as well as in other tumor indications. This revealed a high frequency of tumor-reactive TCR clonotypes in three genetically unstable samples. In contrast, the T cell repertoire in six genetically stable PDAC tumors was largely dominated by bystander T cells. Nevertheless, multiple tumor-reactive TCRs were successfully identified in each of these samples, thereby providing a perspective for personalized immunotherapy in this treatment-resistant indication.

Influence of social lifestyles on host-microbe symbioses in the bees

Microbiomes are increasingly recognised as critical for the health of an organism. In eusocial insect societies, frequent social interactions allow for high-fidelity transmission of microbes across generations, leading to closer host-microbe coevolution. The microbial communities of bees with other social lifestyles are less studied, and few comparisons have been made between taxa that vary in social structure. To address this gap, we leveraged a cloud-computing resource and publicly available transcriptomic data to conduct a survey of microbial diversity in bee samples from a variety of social lifestyles and taxa. We consistently recover the core microbes of well-studied corbiculate bees, supporting this method's ability to accurately characterise microbial communities. We find that the bacterial communities of bees are influenced by host location, phylogeny and social lifestyle, although no clear effect was found for fungal or viral microbial communities. Bee genera with more complex societies tend to harbour more diverse microbes, with Wolbachia detected more commonly in solitary tribes. We present a description of the microbiota of Euglossine bees and find that they do not share the "corbiculate core" microbiome. Notably, we find that bacteria with known anti-pathogenic properties are present across social bee genera, suggesting that symbioses that enhance host immunity are important with higher sociality. Our approach provides an inexpensive means of exploring microbiomes of a given taxa and identifying avenues for further research. These findings contribute to our understanding of the relationships between bees and their associated microbial communities, highlighting the importance of considering microbiome dynamics in investigations of bee health.

Testing for Allele-specific Expression from Human Brain Samples

Many single nucleotide polymorphisms (SNPs) identified by genome-wide association studies exert their effects on disease risk as expression quantitative trait loci (eQTL) via allele-specific expression (ASE). While databases for probing eQTLs in tissues from normal individuals exist, one may wish to ascertain eQTLs or ASE in specific tissues or disease-states not characterized in these databases. Here, we present a protocol to assess ASE of two possible target genes (GPNMB and KLHL7) of a known genome-wide association study (GWAS) Parkinson's disease (PD) risk locus in postmortem human brain tissue from PD and neurologically normal individuals. This was done using a sequence of RNA isolation, cDNA library generation, enrichment for transcripts of interest using customizable cDNA capture probes, paired-end RNA sequencing, and subsequent analysis. This method provides increased sensitivity relative to traditional bulk RNAseq-based and a blueprint that can be extended to the study of other genes, tissues, and disease states. Key features • Analysis of GPNMB allele-specific expression (ASE) in brain lysates from cognitively normal controls (NC) and Parkinson's disease (PD) individuals. • Builds on the ASE protocol of Mayba et al. (2014) and extends application from cells to human tissue. • Increased sensitivity by enrichment for desired transcript via RNA CaptureSeq (Mercer et al., 2014). • Optimized for human brain lysates from cingulate gyrus, caudate nucleus, and cerebellum.

Functional Profiling of Soft Tissue Sarcoma Using Mechanistic Models

Soft tissue sarcoma is an umbrella term for a group of rare cancers that are difficult to treat. In addition to surgery, neoadjuvant chemotherapy has shown the potential to downstage tumors and prevent micrometastases. However, finding effective therapeutic targets remains a research challenge. Here, a previously developed computational approach called mechanistic models of signaling pathways has been employed to unravel the impact of observed changes at the gene expression level on the ultimate functional behavior of cells. In the context of such a mechanistic model, RNA-Seq counts sourced from the Recount3 resource, from The Cancer Genome Atlas (TCGA) Sarcoma project, and non-diseased sarcomagenic tissues from the Genotype-Tissue Expression (GTEx) project were utilized to investigate signal transduction activity through signaling pathways. This approach provides a precise view of the relationship between sarcoma patient survival and the signaling landscape in tumors and their environment. Despite the distinct regulatory alterations observed in each sarcoma subtype, this study identified 13 signaling circuits, or elementary sub-pathways triggering specific cell functions, present across all subtypes, belonging to eight signaling pathways, which served as predictors for patient survival. Additionally, nine signaling circuits from five signaling pathways that highlighted the modifications tumor samples underwent in comparison to normal tissues were found. These results describe the protective role of the immune system, suggesting an anti-tumorigenic effect in the tumor microenvironment, in the process of tumor cell detachment and migration, or the dysregulation of ion homeostasis. Also, the analysis of signaling circuit intermediary proteins suggests multiple strategies for therapy.

MiR-662 is associated with metastatic relapse in early-stage breast cancer and promotes metastasis by stimulating cancer cell stemness

BACKGROUND

Breast cancer (BC) metastasis, which often occurs in bone, contributes substantially to mortality. MicroRNAs play a fundamental role in BC metastasis, although microRNA-regulated mechanisms driving metastasis progression remain poorly understood.

METHODS

MiRome analysis in serum from BC patients was performed by TaqMan™ low-density array. MiR-662 was overexpressed following MIMIC-transfection or lentivirus transduction. Animal models were used to investigate the role of miR-662 in BC (bone) metastasis. The effect of miR-662-overexpressing BC cell conditioned medium on osteoclastogenesis was investigated. ALDEFLUOR assays were performed to study BC stemness. RNA-sequencing transcriptomic analysis of miR-662-overexpressing BC cells was performed to evaluate gene expression changes.

RESULTS

High levels of hsa-miR-662 (miR-662) in serum from BC patients, at baseline (time of surgery), were associated with future recurrence in bone. At an early-stage of the metastatic disease, miR-662 could mask the presence of BC metastases in bone by inhibiting the differentiation of bone-resorbing osteoclasts. Nonetheless, metastatic miR-662-overexpressing BC cells then progressed as overt osteolytic metastases thanks to increased stem cell-like traits.

CONCLUSIONS

MiR-662 is involved in BC metastasis progression, suggesting it may be used as a prognostic marker to identify BC patients at high risk of metastasis.

TDP-43 pathology in Drosophila induces glial-cell type specific toxicity that can be ameliorated by knock-down of SF2/SRSF1

Accumulation of cytoplasmic inclusions of TAR-DNA binding protein 43 (TDP-43) is seen in both neurons and glia in a range of neurodegenerative disorders, including amyotrophic lateral sclerosis (ALS), frontotemporal dementia (FTD) and Alzheimer's disease (AD). Disease progression involves non-cell autonomous interactions among multiple cell types, including neurons, microglia and astrocytes. We investigated the effects in Drosophila of inducible, glial cell type-specific TDP-43 overexpression, a model that causes TDP-43 protein pathology including loss of nuclear TDP-43 and accumulation of cytoplasmic inclusions. We report that TDP-43 pathology in Drosophila is sufficient to cause progressive loss of each of the 5 glial sub-types. But the effects on organismal survival were most pronounced when TDP-43 pathology was induced in the perineural glia (PNG) or astrocytes. In the case of PNG, this effect is not attributable to loss of the glial population, because ablation of these glia by expression of pro-apoptotic reaper expression has relatively little impact on survival. To uncover underlying mechanisms, we used cell-type-specific nuclear RNA sequencing to characterize the transcriptional changes induced by pathological TDP-43 expression. We identified numerous glial cell-type specific transcriptional changes. Notably, SF2/SRSF1 levels were found to be decreased in both PNG and in astrocytes. We found that further knockdown of SF2/SRSF1 in either PNG or astrocytes lessens the detrimental effects of TDP-43 pathology on lifespan, but extends survival of the glial cells. Thus TDP-43 pathology in astrocytes or PNG causes systemic effects that shorten lifespan and SF2/SRSF1 knockdown rescues the loss of these glia, and also reduces their systemic toxicity to the organism.

The role of salinity on genome-wide DNA methylation dynamics in European sea bass gills

Epigenetic modifications, like DNA methylation, generate phenotypic diversity in fish and ultimately lead to adaptive evolutionary processes. Euryhaline marine species that migrate between salinity-contrasted habitats have received little attention regarding the role of salinity on whole-genome DNA methylation. Investigation of salinity-induced DNA methylation in fish will help to better understand the potential role of this process in salinity acclimation. Using whole-genome bisulfite sequencing, we compared DNA methylation patterns in European sea bass (Dicentrarchus labrax) juveniles in seawater and after freshwater transfer. We targeted the gill as a crucial organ involved in plastic responses to environmental changes. To investigate the function of DNA methylation in gills, we performed RNAseq and assessed DNA methylome-transcriptome correlations. We showed a negative correlation between gene expression levels and DNA methylation levels in promoters, first introns and first exons. A significant effect of salinity on DNA methylation dynamics with an overall DNA hypomethylation in freshwater-transferred fish compared to seawater controls was demonstrated. This suggests a role of DNA methylation changes in salinity acclimation. Genes involved in key functions as metabolism, ion transport and transepithelial permeability (junctional complexes) were differentially methylated and expressed between salinity conditions. Expression of genes involved in mitochondrial metabolism (tricarboxylic acid cycle) was increased, whereas the expression of DNA methyltransferases 3a was repressed. This study reveals novel links between DNA methylation, mainly in promoters and first exons/introns, and gene expression patterns following salinity change.

Early clonal extinction in glioblastoma progression revealed by genetic barcoding

Glioblastoma progression in its early stages remains poorly understood. Here, we transfer PDGFB and genetic barcodes in mouse brain to initiate gliomagenesis and enable direct tracing of glioblastoma evolution from its earliest possible stage. Unexpectedly, we observe a high incidence of clonal extinction events and progressive divergence in clonal sizes, even after the acquisition of malignant phenotype. Computational modeling suggests these dynamics result from clonal-based cell-cell competition. Through bulk and single-cell transcriptome analyses, coupled with lineage tracing, we reveal that Myc transcriptional targets have the strongest correlation with clonal size imbalances. Moreover, we show that the downregulation of Myc expression is sufficient to drive competitive dynamics in intracranially transplanted gliomas. Our findings provide insights into glioblastoma evolution that are inaccessible using conventional retrospective approaches, highlighting the potential of combining clonal tracing and transcriptomic analyses in this field.

The genetic basis of variation in immune defense against Lysinibacillus fusiformis infection in Drosophila melanogaster

The genetic causes of phenotypic variation often differ depending on the population examined, particularly if the populations were founded by relatively small numbers of genotypes. Similarly, the genetic causes of phenotypic variation among similar traits (resistance to different xenobiotic compounds or pathogens) may also be completely different or only partially overlapping. Differences in genetic causes for variation in the same trait among populations suggests context dependence for how selection acts on those traits. Similarities in the genetic causes of variation for different traits, on the other hand, suggests pleiotropy which would also influence how natural selection shapes variation in a trait. We characterized immune defense against a natural Drosophila pathogen, the Gram-positive bacterium Lysinibacillus fusiformis, in three different populations and found almost no overlap in the genetic architecture of variation in survival post infection. However, when comparing our results to a similar experiment with the fungal pathogen, B. bassiana, we found a convincing shared QTL peak for both pathogens. This peak contains the Bomanin cluster of Drosophila immune effectors. Loss of function mutants and RNAi knockdown experiments confirms a role of some of these genes in immune defense against both pathogens. This suggests that natural selection may act on the entire cluster of Bomanin genes (and the linked region under the QTL) or specific peptides for specific pathogens.

TDP-43 pathology in Drosophila induces glial-cell type specific toxicity that can be ameliorated by knock-down of SF2/SRSF1

Accumulation of cytoplasmic inclusions of TAR-DNA binding protein 43 (TDP-43) is seen in both neurons and glia in a range of neurodegenerative disorders, including amyotrophic lateral sclerosis (ALS), frontotemporal dementia (FTD) and Alzheimer's disease (AD). Disease progression involves non-cell autonomous interactions among multiple cell types, including neurons, microglia and astrocytes. We investigated the effects in Drosophila of inducible, glial cell type-specific TDP-43 overexpression, a model that causes TDP-43 protein pathology including loss of nuclear TDP-43 and accumulation of cytoplasmic inclusions. We report that TDP-43 pathology in Drosophila is sufficient to cause progressive loss of each of the 5 glial sub-types. But the effects on organismal survival were most pronounced when TDP-43 pathology was induced in the perineural glia (PNG) or astrocytes. In the case of PNG, this effect is not attributable to loss of the glial population, because ablation of these glia by expression of pro-apoptotic reaper expression has relatively little impact on survival. To uncover underlying mechanisms, we used cell-type-specific nuclear RNA sequencing to characterize the transcriptional changes induced by pathological TDP-43 expression. We identified numerous glial cell-type specific transcriptional changes. Notably, SF2/SRSF1 levels were found to be decreased in both PNG and in astrocytes. We found that further knockdown of SF2/SRSF1 in either PNG or astrocytes lessens the detrimental effects of TDP-43 pathology on lifespan, but extends survival of the glial cells. Thus TDP-43 pathology in astrocytes or PNG causes systemic effects that shorten lifespan and SF2/SRSF1 knockdown rescues the loss of these glia, and also reduces their systemic toxicity to the organism.

A modERN resource: identification of Drosophila transcription factor candidate target genes using RNAi

Transcription factors (TFs) play a key role in development and in cellular responses to the environment by activating or repressing the transcription of target genes in precise spatial and temporal patterns. In order to develop a catalog of target genes of Drosophila melanogaster TFs, the modERN consortium systematically knocked down the expression of TFs using RNAi in whole embryos followed by RNA-seq. We generated data for 45 TFs which have 18 different DNA-binding domains and are expressed in 15 of the 16 organ systems. The range of inactivation of the targeted TFs by RNAi ranged from log2fold change -3.52 to +0.49. The TFs also showed remarkable heterogeneity in the numbers of candidate target genes identified, with some generating thousands of candidates and others only tens. We present detailed analysis from five experiments, including those for three TFs that have been the focus of previous functional studies (ERR, sens, and zfh2) and two previously uncharacterized TFs (sens-2 and CG32006), as well as short vignettes for selected additional experiments to illustrate the utility of this resource. The RNA-seq datasets are available through the ENCODE DCC (http://encodeproject.org) and the Sequence Read Archive (SRA). TF and target gene expression patterns can be found here: https://insitu.fruitfly.org. These studies provide data that facilitate scientific inquiries into the functions of individual TFs in key developmental, metabolic, defensive, and homeostatic regulatory pathways, as well as provide a broader perspective on how individual TFs work together in local networks during embryogenesis.

Splicing accuracy varies across human introns, tissues and age

Alternative splicing impacts most multi-exonic human genes. Inaccuracies during this process may have an important role in ageing and disease. Here, we investigated mis-splicing using RNA-sequencing data from ~14K control samples and 42 human body sites, focusing on split reads partially mapping to known transcripts in annotation. We show that mis-splicing occurs at different rates across introns and tissues and that these splicing inaccuracies are primarily affected by the abundance of core components of the spliceosome assembly and its regulators. Using publicly available data on short-hairpin RNA-knockdowns of numerous spliceosomal components and related regulators, we found support for the importance of RNA-binding proteins in mis-splicing. We also demonstrated that age is positively correlated with mis-splicing, and it affects genes implicated in neurodegenerative diseases. This in-depth characterisation of mis-splicing can have important implications for our understanding of the role of splicing inaccuracies in human disease and the interpretation of long-read RNA-sequencing data.

Pseudo-temporal dynamics of chemoresistant triple negative breast cancer cells reveal EGFR/HER2 inhibition as synthetic lethal during mid-neoadjuvant chemotherapy

In the absence of targetable hormonal axes, chemoresistance for triple-negative breast cancer (TNBC) often compromises patient outcomes. To investigate the underlying tumor dynamics, we performed trajectory analysis on the single-nuclei RNA-seq (snRNA-seq) of chemoresistant tumor clones during neoadjuvant chemotherapy (NAC). It revealed a common tumor trajectory across multiple patients with HER2-like expansions during NAC. Genome-wide CRISPR-Cas9 knock-out on mammary epithelial cells revealed chemosensitivity-promoting knock-outs were up-regulated along the tumor trajectory. Furthermore, we derived a consensus gene signature of TNBC chemoresistance by comparing the trajectory transcriptome with chemoresistant transcriptomes from TNBC cell lines and poor prognosis patient samples to predict FDA-approved drugs, including afatinib (pan-HER inhibitor), targeting the consensus signature. We validated the synergistic efficacy of afatinib and paclitaxel in chemoresistant TNBC cells and confirmed pharmacological suppression of the consensus signature. The study provides a dynamic model of chemoresistant tumor transcriptome, and computational framework for pharmacological intervention.

A multi-omic analysis reveals the esophageal dysbiosis as the predominant trait of eosinophilic esophagitis

BACKGROUND

Eosinophilic esophagitis (EoE) is a chronic immune-mediated rare disease, characterized by esophageal dysfunctions. It is likely to be primarily activated by food antigens and is classified as a chronic disease for most patients. Therefore, a deeper understanding of the pathogenetic mechanisms underlying EoE is needed to implement and improve therapeutic lines of intervention and ameliorate overall patient wellness.

METHODS

RNA-seq data of 18 different studies on EoE, downloaded from NCBI GEO with faster-qdump ( https://github.com/ncbi/sra-tools ), were batch-corrected and analyzed for transcriptomics and metatranscriptomics profiling as well as biological process functional enrichment. The EoE TaMMA web app was designed with plotly and dash. Tabula Sapiens raw data were downloaded from the UCSC Cell Browser. Esophageal single-cell raw data analysis was performed within the Automated Single-cell Analysis Pipeline. Single-cell data-driven bulk RNA-seq data deconvolution was performed with MuSiC and CIBERSORTx. Multi-omics integration was performed with MOFA.

RESULTS

The EoE TaMMA framework pointed out disease-specific molecular signatures, confirming its reliability in reanalyzing transcriptomic data, and providing new EoE-specific molecular markers including CXCL14, distinguishing EoE from gastroesophageal reflux disorder. EoE TaMMA also revealed microbiota dysbiosis as a predominant characteristic of EoE pathogenesis. Finally, the multi-omics analysis highlighted the presence of defined classes of microbial entities in subsets of patients that may participate in inducing the antigen-mediated response typical of EoE pathogenesis.

CONCLUSIONS

Our study showed that the complex EoE molecular network may be unraveled through advanced bioinformatics, integrating different components of the disease process into an omics-based network approach. This may implement EoE management and treatment in the coming years.

Transcriptomic comparison of seeds and silique walls from two rapeseed genotypes with contrasting seed oil content

Silique walls play pivotal roles in contributing photoassimilates and nutrients to fuel seed growth. However, the interaction between seeds and silique walls impacting oil biosynthesis is not clear during silique development. Changes in sugar, fatty acid and gene expression during Brassica napus silique development of L192 with high oil content and A260 with low oil content were investigated to identify key factors affecting difference of their seed oil content. During the silique development, silique walls contained more hexose and less sucrose than seeds, and glucose and fructose contents in seeds and silique walls of L192 were higher than that of A260 at 15 DAF, and sucrose content in the silique walls of L192 were lower than that of A260 at three time points. Genes related to fatty acid biosynthesis were activated over time, and differences on fatty acid content between the two genotypes occurred after 25 DAF. Genes related to photosynthesis expressed more highly in silique walls than in contemporaneous seeds, and were inhibited over time. Gene set enrichment analysis suggested photosynthesis were activated in L192 at 25 and 35 DAF in silique walls and at both 15 and 35 DAF in the seed. Expressions of sugar transporter genes in L192 was higher than that in A260, especially at 35 DAF. Expressions of genes related to fatty acid biosynthesis, such as BCCP2s, bZIP67 and LEC1s were higher in L192 than in A260, especially at 35 DAF. Meanwhile, genes related to oil body proteins were expressed at much lower levels in L192 than in A260. According to the WGCNA results, hub modules, such as ME.turquoise relative to photosynthesis, ME.green relative to embryo development and ME.yellow relative to lipid biosynthesis, were identified and synergistically regulated seed development and oil accumulation. Our results are helpful for understanding the mechanism of oil accumulation of seeds in oilseed rape for seed oil content improvement.

Comparison of transcriptome profiles between medulloblastoma primary and recurrent tumors uncovers novel variance effects in relapses

Nowadays medulloblastoma (MB) tumors can be treated with risk-stratified approaches with up to 80% success rate. However, disease relapses occur in approximately 30% of patients and successful salvage treatment strategies at relapse remain scarce. Acquired copy number changes or TP53 mutations are known to occur frequently in relapses, while methylation profiles usually remain highly similar to those of the matching primary tumors, indicating that in general molecular subgrouping does not change during the course of the disease. In the current study, we have used RNA sequencing data to analyze the transcriptome profiles of 43 primary-relapse MB pairs in order to identify specific molecular features of relapses within various tumor groups. Gene variance analysis between primary and relapse samples demonstrated the impact of age in SHH-MB: the changes in gene expression relapse profiles were more pronounced in the younger patients (< 10 years old), which were also associated with increased DNA aberrations and somatic mutations at relapse probably driving this effect. For Group 3/4 MB transcriptome data analysis uncovered clear sets of genes either active or decreased at relapse that are significantly associated with survival, thus could be potential predictive markers. In addition, deconvolution analysis of bulk transcriptome data identified progression-associated differences in cell type enrichment. The proportion of undifferentiated progenitors increased in SHH-MB relapses with a concomitant decrease of differentiated neuron-like cells, while in Group 3/4 MB relapses cell cycle activity increases and differentiated neuron-like cells proportion decreases as well. Thus, our findings uncovered significant transcriptome changes in the molecular signatures of relapsed MB and could be potentially useful for further clinical purposes.

Gene expression profiling of human macrophages after graphene oxide and graphene nanoplatelets treatment reveals particle-specific regulation of pathways

Graphene and its derivatives are attractive materials envisaged to enable a wealth of novel applications in many fields including energy, electronics, composite materials or health. A comprehensive understanding of the potential adverse effects of graphene-related materials (GRM) in humans is a prerequisite to the safe use of these promising materials. Here, we exploited gene expression profiling to identify transcriptional responses and toxicity pathways induced by graphene oxide (GO) and graphene nanoplatelets (GNP) in human macrophages. Primary human monocyte-derived macrophages (MDM) and a human macrophage cell line, i.e. differentiated THP-1 cells, were exposed to 5 or 20 μg/mL GO and GNP for 6 and 24 h to capture early and more persistent acute responses at realistic or slightly overdose concentrations. GO and GNP induced time-, dose- and macrophage type-specific differential expression of a substantial number of genes with some overlap between the two GRM types (up to 384 genes (9.6%) or 447 genes (20.4%) in THP-1 or MDM, respectively) but also a high number of genes exclusively deregulated from each material type. Furthermore, GRM responses on gene expression were highly different from those induced by inflammogenic material crystalline quartz (maximum of 64 (2.3%) or 318 (11.3%) common genes for MDM treated with 20 μg/mL GO and GNP, respectively). Further bioinformatics analysis revealed that GNP predominantly activated genes controlling inflammatory and apoptotic pathways whereas GO showed only limited inflammatory responses. Interestingly, both GRM affected the expression of genes related to antigen processing and presentation and in addition, GO activated pathways of neutrophil activation, degranulation and immunity in MDM. Overall, this study provides an extensive resource of potential toxicity mechanisms for future safety assessment of GRM in more advanced model systems to verify if the observed changes in gene expression in human macrophages could lead to long-term consequences on human health.

In vivo genome editing using novel AAV-PHP variants rescues motor function deficits and extends survival in a SOD1-ALS mouse model

CRISPR-based gene editing technology represents a promising approach to deliver therapies for inherited disorders, including amyotrophic lateral sclerosis (ALS). Toxic gain-of-function superoxide dismutase 1 (SOD1) mutations are responsible for ~20% of familial ALS cases. Thus, current clinical strategies to treat SOD1-ALS are designed to lower SOD1 levels. Here, we utilized AAV-PHP.B variants to deliver CRISPR-Cas9 guide RNAs designed to disrupt the human SOD1 (huSOD1) transgene in SOD1^G93A mice. A one-time intracerebroventricular injection of AAV.PHP.B-huSOD1-sgRNA into neonatal H11^Cas9 SOD1^G93A mice caused robust and sustained mutant huSOD1 protein reduction in the cortex and spinal cord, and restored motor function. Neonatal treatment also reduced spinal motor neuron loss, denervation at neuromuscular junction (NMJ) and muscle atrophy, diminished axonal damage and preserved compound muscle action potential throughout the lifespan of treated mice. SOD1^G93A treated mice achieved significant disease-free survival, extending lifespan by more than 110 days. Importantly, a one-time intrathecal or intravenous injection of AAV.PHP.eB-huSOD1-sgRNA in adult H11^Cas9 SOD1^G93A mice, immediately before symptom onset, also extended lifespan by at least 170 days. We observed substantial protection against disease progression, demonstrating the utility of our CRISPR editing preclinical approach for target evaluation. Our approach uncovered key parameters (e.g., AAV capsid, Cas9 expression) that resulted in improved efficacy compared to similar approaches and can also serve to accelerate drug target validation.

Altered tRNA processing is linked to a distinct and unusual La protein in Tetrahymena thermophila

Nascent pre-tRNAs are transcribed by RNA polymerase III and immediately bound by La proteins on the UUU-3'OH sequence, using a tandem arrangement of the La motif and an adjacent RNA recognition motif-1 (RRM1), resulting in protection from 3'-exonucleases and promotion of pre-tRNA folding. The Tetrahymena thermophila protein Mlp1 has been previously classified as a genuine La protein, despite the predicted absence of the RRM1. We find that Mlp1 functions as a La protein through binding of pre-tRNAs, and affects pre-tRNA processing in Tetrahymena thermophila and when expressed in fission yeast. However, unlike in other examined eukaryotes, depletion of Mlp1 results in 3'-trailer stabilization. The 3'-trailers in Tetrahymena thermophila are uniquely short relative to other examined eukaryotes, and 5'-leaders have evolved to disfavour pre-tRNA leader/trailer pairing. Our data indicate that this variant Mlp1 architecture is linked to an altered, novel mechanism of tRNA processing in Tetrahymena thermophila.

IntroVerse: a comprehensive database of introns across human tissues

Dysregulation of RNA splicing contributes to both rare and complex diseases. RNA-sequencing data from human tissues has shown that this process can be inaccurate, resulting in the presence of novel introns detected at low frequency across samples and within an individual. To enable the full spectrum of intron use to be explored, we have developed IntroVerse, which offers an extensive catalogue on the splicing of 332,571 annotated introns and a linked set of 4,679,474 novel junctions covering 32,669 different genes. This dataset has been generated through the analysis of 17,510 human control RNA samples from 54 tissues provided by the Genotype-Tissue Expression Consortium. IntroVerse has two unique features: (i) it provides a complete catalogue of novel junctions and (ii) each novel junction has been assigned to a specific annotated intron. This unique, hierarchical structure offers multiple uses, including the identification of novel transcripts from known genes and their tissue-specific usage, and the assessment of background splicing noise for introns thought to be mis-spliced in disease states. IntroVerse provides a user-friendly web interface and is freely available at https://rytenlab.com/browser/app/introverse.

Identification of candidate genes associated with resistance against race 0 of Colletotrichum lentis in Lens ervoides

Resistance to anthracnose caused by the fungal pathogen Colletotrichum lentis was explored through transcriptome sequencing over a period of 24 to 96 h post-inoculation (hpi) of the partially resistant recombinant inbred lines (RIL) LR-66-528 and susceptible LR-66-524 of the crop wild relative Lens ervoides population LR-66. The development of infection vesicles and primary hyphae by C. lentis were significantly higher on susceptible RIL LR-66-524 compared to partially resistant LR-66-528 at 24 and 48 hpi, but exponential trends in fungal growth were observed between 24 to 96 hpi in both RILs. Comparison of inoculated with mock-inoculated samples revealed 3091 disease responsive genes, among which 477 were differentially expressed between the two RILs. These were clustered into six expression clusters with genes that had either high or low expression in one of the RILs. Differentially expressed genes (DEGs) were functionally annotated and included genes coding LRR and NB-ARC domain disease resistance proteins, protein detoxification, LRR receptor-like kinase family proteins, and wall-associated Ser/Thr Kinases. DEGs were compared to genes in previously published anthracnose resistance QTLs mapped in LR-66 and revealed 22 DEGs located in 3 QTLs. Expression of 21 DEGs was validated using RT-qPCR confirming expression trends in RNA-seq.

Venus: An efficient virus infection detection and fusion site discovery method using single-cell and bulk RNA-seq data

Early and accurate detection of viruses in clinical and environmental samples is essential for effective public healthcare, treatment, and therapeutics. While PCR detects potential pathogens with high sensitivity, it is difficult to scale and requires knowledge of the exact sequence of the pathogen. With the advent of next-gen single-cell sequencing, it is now possible to scrutinize viral transcriptomics at the finest possible resolution–cells. This newfound ability to investigate individual cells opens new avenues to understand viral pathophysiology with unprecedented resolution. To leverage this ability, we propose an efficient and accurate computational pipeline, named Venus, for virus detection and integration site discovery in both single-cell and bulk-tissue RNA-seq data. Specifically, Venus addresses two main questions: whether a tissue/cell type is infected by viruses or a virus of interest? And if infected, whether and where has the virus inserted itself into the human genome? Our analysis can be broken into two parts–validation and discovery. Firstly, for validation, we applied Venus on well-studied viral datasets, such as HBV- hepatocellular carcinoma and HIV-infection treated with antiretroviral therapy. Secondly, for discovery, we analyzed datasets such as HIV-infected neurological patients and deeply sequenced T-cells. We detected viral transcripts in the novel target of the brain and high-confidence integration sites in immune cells. In conclusion, here we describe Venus, a publicly available software which we believe will be a valuable virus investigation tool for the scientific community at large.

GPNMB confers risk for Parkinson's disease through interaction with α-synuclein

Many risk loci for Parkinson's disease (PD) have been identified by genome-wide association studies (GWASs), but target genes and mechanisms remain largely unknown. We linked the GWAS-derived chromosome 7 locus (sentinel single-nucleotide polymorphism rs199347) to GPNMB through colocalization analyses of expression quantitative trait locus and PD risk signals, confirmed by allele-specific expression studies in the human brain. In cells, glycoprotein nonmetastatic melanoma protein B (GPNMB) coimmunoprecipitated and colocalized with α-synuclein (aSyn). In induced pluripotent stem cell-derived neurons, loss of GPNMB resulted in loss of ability to internalize aSyn fibrils and develop aSyn pathology. In 731 PD and 59 control biosamples, GPNMB was elevated in PD plasma, associating with disease severity. Thus, GPNMB represents a PD risk gene with potential for biomarker development and therapeutic targeting.

Identification of circular RNA biomarkers for Pien Tze Huang treatment of CCl4‑induced liver fibrosis using RNA‑sequencing

Pien Tze Huang (PZH), a common hepatoprotective Traditional Chinese Medicine that has been found to be an effective treatment for carbon tetrachloride-induced hepatic damage, including liver fibrosis. Circular RNAs (circRNAs) serve a crucial role in regulating gene expression levels via circRNA/micro (mi)RNA/mRNA networks in several human diseases and biological processes. However, whether circRNAs are involved in the underlying mechanism of the therapeutic effects of PZH on liver fibrosis remains unclear. Therefore, the aim of the present study was to investigate these effects using circRNA expression profiles from PZH-treated fibrotic livers in model mice. A case-control study on >59,476 circRNAs from CCl₄-induced (control group, n=6) and PZH-treated (case group, n=6) mice was performed using circRNA sequencing in liver tissues. PZH treatment resulted in the differential expression of 91 circRNAs, including 58 upregulated and 33 downregulated circRNAs. Furthermore, the construction of competing endogenous networks also indicated that differentially expressed circRNAs acted as miRNA sponges. Gene Ontology and Kyoto Encyclopedia of Genes and Genomes pathway enrichment analysis of miRNA targets demonstrated that PZH-affected circRNAs were mainly involved in biological processes such as ‘positive regulation of fibroblast proliferation’, ‘cellular response to interleukin-1’ and ‘regulation of DNA-templated transcription in response to stress’ and in a number of important pathways, such as ‘TNF signaling pathway’, ‘PI3K-Akt signaling pathway’, ‘IL-17 signaling pathway’ and ‘MAPK signaling pathway’. To further validate the bioinformatics data, reverse transcription–quantitative PCR was performed on seven miRNA targets in a human hepatic stellate LX-2 cell model. The results suggested that seven of the miRNAs exhibited regulatory patterns that were consistent with those of the transcriptome sequencing results. Kaplan-Meier survival analysis demonstrated that the expression levels of dihydrodiol dehydrogenase and solute carrier family 7, member 11 gene were significantly associated with patient survival, 269 patients with liver hepatocellular carcinoma from The Cancer Genome Atlas database. To the best of our knowledge, this was the first study to provide evidence that PZH affects circRNA expression levels, which may serve important roles in PZH-treated fibrotic liver through the regulation of functional gene expression. In conclusion, the present study provided new insights into the mechanism underlying the pathogenesis of liver fibrosis and identified potential novel, efficient, therapeutic targets against liver injury.

Adenoviral-based vaccine promotes neoantigen-specific CD8+ T cell stemness and tumor rejection

Upon chronic antigen exposure, CD8⁺ T cells become exhausted, acquiring a dysfunctional state correlated with the inability to control infection or tumor progression. In contrast, stem-like CD8⁺ T progenitors maintain the ability to promote and sustain effective immunity. Adenovirus (Ad)-vectored vaccines encoding tumor neoantigens have been shown to eradicate large tumors when combined with anti-programmed cell death protein 1 (αPD-1) in murine models; however, the mechanisms and translational potential have not yet been elucidated. Here, we show that gorilla Ad vaccine targeting tumor neoepitopes enhances responses to αPD-1 therapy by improving immunogenicity and antitumor efficacy. Single-cell RNA sequencing demonstrated that the combination of Ad vaccine and αPD-1 increased the number of murine polyfunctional neoantigen-specific CD8⁺ T cells over αPD-1 monotherapy, with an accumulation of Tcf1⁺ stem-like progenitors in draining lymph nodes and effector CD8⁺ T cells in tumors. Combined T cell receptor (TCR) sequencing analysis highlighted a broader spectrum of neoantigen-specific CD8⁺ T cells upon vaccination compared to αPD-1 monotherapy. The translational relevance of these data is supported by results obtained in the first 12 patients with metastatic deficient mismatch repair (dMMR) tumors vaccinated with an Ad vaccine encoding shared neoantigens. Expansion and diversification of TCRs were observed in post-treatment biopsies of patients with clinical response, as well as an increase in tumor-infiltrating T cells with an effector memory signature. These findings indicate a promising mechanism to overcome resistance to PD-1 blockade by promoting immunogenicity and broadening the spectrum and magnitude of neoantigen-specific T cells infiltrating tumors.

IL-17A Increases Doxorubicin Efficacy in Triple Negative Breast Cancer

Due to lack of targetable receptors and intertumoral heterogeneity, triple negative breast cancer (TNBC) remains particularly difficult to treat. Doxorubicin (DOX) is typically used as nonselective neoadjuvant chemotherapy, but the diversity of treatment efficacy remains unclear. Comparable to variability in clinical response, an experimental model of TNBC using a 4T1 syngeneic mouse model was found to elicit a differential response to a seven-day treatment regimen of DOX. Single-cell RNA sequencing identified an increase in T cells in tumors that responded to DOX treatment compared to tumors that continued to grow uninhibited. Additionally, compared to resistant tumors, DOX sensitive tumors contained significantly more CD4 T helper cells (339%), γδ T cells (727%), Naïve T cells (278%), and activated CD8 T cells (130%). Furthermore, transcriptional profiles of tumor infiltrated T cells in DOX responsive tumors revealed decreased exhaustion, increased chemokine/cytokine expression, and increased activation and cytotoxic activity. γδ T cell derived IL-17A was identified to be highly abundant in the sensitive tumor microenvironment. IL-17A was also found to directly increase sensitivity of TNBC cells in combination with DOX treatment. In TNBC tumors sensitive to DOX, increased IL-17A levels lead to a direct effect on cancer cell responsiveness and chronic stimulation of tumor infiltrated T cells leading to improved chemotherapeutic efficacy. IL-17A’s role as a chemosensitive cytokine in TNBC may offer new opportunities for treating chemoresistant breast tumors and other cancer types.

Transcriptome Profiling of Transposon-Derived Long Non-coding RNAs Response to Hormone in Strawberry Fruit Development

Strawberry is an economically grown horticulture crop required for fruit consumption. The ripening of its fruit is a complex biological process regulated by various hormones. Abscisic acid (ABA) is a critical phytohormone involved in fruit ripening. However, little is known about the long non-coding RNAs (LncRNAs), especially transposon-derived LncRNA (TE-lncRNA), response to hormones during fruit ripening in octoploid strawberry. In the study, the transcriptome data of developing strawberry fruits treated with ABA and its inhibitor Nordihydroguaiaretic acid (NGDA) were analyzed to identify responsive LncRNAs and coding genes. A total of 14,552 LncRNAs were identified, including 8,617 transposon-derived LncRNAs (TE-LncRNAs), 412 LncRNAs (282 TE-LncRNAs), and 382 ABA-sensitive LncRNAs (231 TE-LncRNAs). Additionally, a weighted co-expression network analysis constructed 27 modules containing coding RNAs and LncRNAs. Seven modules, including "MEdarkorange" and "MElightyellow" were significantly correlated with ABA/NDGA treatments, resulting in 247 hub genes, including 21 transcription factors and 22 LncRNAs (15 TE-LncRNAs). Gene ontology enrichment analysis further revealed that ABA/NDGA-responsive modules, including LncRNAs, were associated with various metabolic pathways involved in strawberry fruit development and ripening, including lipid metabolism, organic acid metabolism, and phenylpropanoid metabolism. The current study identifies many high-confidence LncRNAs in strawberry, with a percentage of them being ABA pathway-specific and 22 hub-responsive LncRNAs, providing new insight into strawberry or other Rosaceae crop fruit ripening.

Dynamic 3D genome reorganization during development and metabolic stress of the porcine liver

Liver development is a complex process that is regulated by a series of signaling pathways. Three-dimensional (3D) chromatin architecture plays an important role in transcriptional regulation; nonetheless, its dynamics and role in the rapid transition of core liver functions during development and obesity-induced metabolic stress remain largely unexplored. To investigate the dynamic chromatin architecture during liver development and under metabolic stress, we generated high-resolution maps of chromatin architecture for porcine livers across six major developmental stages (from embryonic day 38 to the adult stage) and under a high-fat diet-induced obesity. The characteristically loose chromatin architecture supports a highly plastic genome organization during early liver development, which fundamentally contributes to the rapid functional transitions in the liver after birth. We reveal the multi-scale reorganization of chromatin architecture and its influence on transcriptional regulation of critical signaling processes during liver development, and show its close association with transition in hepatic functions (i.e., from hematopoiesis in the fetus to metabolism and immunity after birth). The limited changes in chromatin structure help explain the observed metabolic adaptation to excessive energy intake in pigs. These results provide a global overview of chromatin architecture dynamics associated with the transition of physiological liver functions between prenatal development and postnatal maturation, and a foundational resource that allows for future in-depth functional characterization.

Chimeric transcripts observed in non-canonical FGFR2 fusions with partner genes' breakpoint located in intergenic region in intrahepatic cholangiocarcinoma

Intrahepatic cholangiocarcinoma (ICC) is a fatal bile duct cancer with dismal prognosis and limited therapeutic options. FGFR family fusion have been identified in many diseases, and FGFR2 fusion is a validated oncogenic driver in ICC. At present, a variety of fusion forms have been reported, including gene-gene, gene-intergenic, and intergenic-intergenic fusion. Here, by performing RNA- and DNA-sequencing analysis, FGFR2 fusions were found in 10.1% of ICC, including 4 gene-intergenic fusions. We confirmed that the non-canonical rearrangements can generate chimeric transcripts, and used conventional splicing mechanism to explain the event. Our study provides possible target therapy for these 4 patients and possibility analysis scheme for similar situation.

Recombination landscape dimorphism and sex chromosome evolution in the dioecious plant Rumex hastatulus

There is growing evidence from diverse taxa for sex differences in the genomic landscape of recombination, but the causes and consequences of these differences remain poorly understood. Strong recombination landscape dimorphism between the sexes could have important implications for the dynamics of sex chromosome evolution because low recombination in the heterogametic sex can favour the spread of sexually antagonistic alleles. Here, we present a sex-specific linkage map and revised genome assembly of Rumex hastatulus and provide the first evidence and characterization of sex differences in recombination landscape in a dioecious plant. We present data on significant sex differences in recombination, with regions of very low recombination in males covering over half of the genome. This pattern is evident on both sex chromosomes and autosomes, suggesting that pre-existing differences in recombination may have contributed to sex chromosome formation and divergence. Our analysis of segregation distortion suggests that haploid selection due to pollen competition occurs disproportionately in regions with low male recombination. We hypothesize that sex differences in the recombination landscape have contributed to the formation of a large heteromorphic pair of sex chromosomes in R. hastatulus, but more comparative analyses of recombination will be important to investigate this hypothesis further. This article is part of the theme issue 'Sex determination and sex chromosome evolution in land plants'.

Transcriptome in paraffin samples for the diagnosis and prognosis of adrenocortical carcinoma

DESIGN

Molecular classification is important for the diagnosis and prognosis of adrenocortical tumors (ACT). Transcriptome profiles separate adrenocortical adenomas 'C2' from carcinomas, and identify two groups of carcinomas 'C1A' and 'C1B', of poor and better prognosis respectively. However, many ACT cannot be profiled because of improper or absent freezing procedures, a mandatory requirement so far. The main aim was to determine transcriptome profiles on formalin-fixed paraffin-embedded (FFPE) samples, using the new 3'-end RNA-sequencing technology. A secondary aim was to demonstrate the ability of this technique to explore large FFPE archives, by focusing on the rare oncocytic ACT variants.

METHODS

We included 131 ACT: a training cohort from Cochin hospital and an independent validation cohort from Wuerzburg hospital. The 3' transcriptome was generated from FFPE samples using QuantSeq (Lexogen, Vienna, Austria) and NextSeq500 (Illumina, San Diego, CA, USA).

RESULTS

In the training cohort, unsupervised clustering identified three groups: 'C1A' aggressive carcinomas (n = 28, 29%), 'C1B' more indolent carcinomas (n = 28, 29%), and 'C2' adenomas (n = 39, 41%). The prognostic value of FFPE transcriptome was confirmed in the validation cohort (5-year OS: 26% in 'C1A' (n = 26) and 100% in 'C1B' (n = 10), P = 0.003). FFPE transcriptome was an independent prognostic factor in a multivariable model including tumor stage and Ki-67 (OS HR: 7.5, P = 0.01). Oncocytic ACT (n = 19) did not form any specific cluster. Oncocytic carcinomas (n = 6) and oncocytic ACT of uncertain malignant potential (n = 4) were all in 'C1B'.

CONCLUSIONS

The 3' RNA-sequencing represents a convenient solution for determining ACT molecular class from FFPE samples. This technique should facilitate routine use and large retrospective studies.

Cellular Model of Malignant Transformation of Primary Human Astrocytes Induced by Deadhesion/Readhesion Cycles

Astrocytoma is the most common and aggressive tumor of the central nervous system. Genetic and environmental factors, bacterial infection, and several other factors are known to be involved in gliomagenesis, although the complete underlying molecular mechanism is not fully understood. Tumorigenesis is a multistep process involving initiation, promotion, and progression. We present a human model of malignant astrocyte transformation established by subjecting primary astrocytes from healthy adults to four sequential cycles of forced anchorage impediment (deadhesion). After limiting dilution of the surviving cells obtained after the fourth deadhesion/readhesion cycle, three clones were randomly selected, and exhibited malignant characteristics, including increased proliferation rate and capacity for colony formation, migration, and anchorage-independent growth in soft agar. Functional assay results for these clonal cells, including response to temozolomide, were comparable to U87MG—a human glioblastoma-derived cell lineage—reinforcing malignant cell transformation. RNA-Seq analysis by next-generation sequencing of the transformed clones relative to the primary astrocytes revealed upregulation of genes involved in the PI3K/AKT and Wnt/β-catenin signaling pathways, in addition to upregulation of genes related to epithelial–mesenchymal transition, and downregulation of genes related to aerobic respiration. These findings, at a molecular level, corroborate the change in cell behavior towards mesenchymal-like cell dedifferentiation. This linear progressive model of malignant human astrocyte transformation is unique in that neither genetic manipulation nor treatment with carcinogens are used, representing a promising tool for testing combined therapeutic strategies for glioblastoma patients, and furthering knowledge of astrocytoma transformation and progression.

Comprehensive Transcriptome Analysis Uncovers Hub Long Non-coding RNAs Regulating Potassium Use Efficiency in Nicotiana tabacum

Potassium (K) is the essential element for plant growth. It is one of the critical factors that determine crop yield, quality, and especially leaf development in tobacco. However, the molecular mechanism of potassium use efficiency (KUE), especially non-coding RNA, is still unknown. In this study, tobacco seedlings were employed, and their hydro-cultivation with K treatments of low and sufficient concentrations was engaged. Physiological analysis showed that low potassium treatment could promote malondialdehyde (MDA) accumulation and antioxidant enzyme activities such as peroxidase (POD), ascorbate-peroxidase (APX). After transcriptomic analysis, a total of 10,585 LncRNA transcripts were identified, and 242 of them were significantly differently expressed under potassium starvation. Furthermore, co-expression networks were constructed and generated 78 potential regulation modules in which coding gene and LncRNAs are involved and functional jointly. By further module-trait analysis and module membership (MM) ranking, nine modules, including 616 coding RNAs and 146 LncRNAs, showed a high correlation with K treatments, and 20 hub K-responsive LncRNAs were finally predicted. Following gene ontology (GO) analysis, the results showed potassium starvation inducing the pathway of antioxidative stress which is consistent with the physiology result mentioned above. Simultaneously, a part of detected LncRNAs, such as MSTRG.6626.1, MSTRG.11330.1, and MSTRG.16041.1, were co-relating with a bench of MYB, C3H, and NFYC transcript factors in response to the stress. Overall, this research provided a set of LncRNAs that respond to K concentration from starvation and sufficient supply. Simultaneously, the regulation network and potential co-functioning genes were listed as well. This massive dataset would serve as an outstanding clue for further study in tobacco and other plant species for nutrient physiology and molecular regulation mechanism.

Peering Into Candida albicans Pir Protein Function and Comparative Genomics of the Pir Family

The fungal cell wall, comprised primarily of protein and polymeric carbohydrate, maintains cell structure, provides protection from the environment, and is an important antifungal drug target. Pir proteins (proteins with internal repeats) are linked to cell wall β-1,3-glucan and are best studied in Saccharomyces cerevisiae. Sequential deletion of S. cerevisiae PIR genes produces strains with increasingly notable cell wall damage. However, a true null mutant lacking all five S. cerevisiae PIR genes was never constructed. Because only two PIR genes (PIR1, PIR32) were annotated in the Candida albicans genome, the initial goal of this work was to construct a true Δpir/Δpir null strain in this species. Unexpectedly, the phenotype of the null strain was almost indistinguishable from its parent, leading to the search for other proteins with Pir function. Bioinformatic approaches revealed nine additional C. albicans proteins that share a conserved Pir functional motif (minimally DGQ). Examination of the protein sequences revealed another conserved motif (QFQFD) toward the C-terminal end of each protein. Sequence similarities and presence of the conserved motif(s) were used to identify a set of 75 proteins across 16 fungal species that are proposed here as Pir proteins. The Pir family is greatly expanded in C. albicans and C. dubliniensis compared to other species and the orthologs are known to have specialized function during chlamydospore formation. Predicted Pir structures showed a conserved core of antiparallel beta-sheets and sometimes-extensive loops that contain amino acids with the potential to form linkages to cell wall components. Pir phylogeny demonstrated emergence of specific ortholog groups among the fungal species. Variation in gene expression patterns was noted among the ortholog groups during growth in rich medium. PIR allelic variation was quite limited despite the presence of a repeated sequence in many loci. Results presented here demonstrate that the Pir family is larger than previously recognized and lead to new hypotheses to test to better understand Pir proteins and their role in the fungal cell wall.

Transcriptomic analysis of papillary thyroid cancer focused on immune-subtyping, oncogenic fusion, and recurrence

Objectives

Thyroid cancer is the most common endocrine tumor, with rapidly increasing incidence worldwide. However, its transcriptomic characteristics associated with immunological signatures, driver fusions, and recurrence markers remain unclear. We aimed to investigate the transcriptomic characteristics of advanced papillary thyroid cancer.

Methods

This study included 282 papillary thyroid cancer tumor samples and 155 normal samples from Chungnam National University Hospital and Seoul National University Hospital. Transcriptomic quantification was determined by high-throughput RNA sequencing. We investigated the associations of clinical parameters and molecular signatures using RNA sequencing. We validated predictive biomarkers using the Cancer Genome Atlas database.

Results

Through a comparison of differentially expressed genes, gene sets, and pathways in papillary thyroid cancer compared to normal tumor-adjacent tissue, we found increased immune signaling associated with cytokines or T cells and decreased thyroid hormone synthetic pathways. In addition, patients with recurrence presented increased CD8+ T-cell and Th1-cell signatures. Interestingly, we found differentially overexpressed genes related to immune-escape signaling such as CTLA4, IDO1, LAG3, and PDCD1 in advanced papillary thyroid cancer with a low thyroid differentiation score. Fusion analysis showed that the PI3K and mitogen-activated protein kinase (MAPK) signaling pathways were regulated differently according to the RET fusion partner genes (CCDC6 or NCOA4). Finally, we identified HOXD9 as a novel molecular biomarker that predicts the recurrence of thyroid cancer in addition to known risk factors (tumor size, lymph node metastasis, and extrathyroidal extension).

Conclusion

We identified a high association with immune-escape signaling in the immune-hot group with aggressive clinical characteristics among Korean thyroid cancer patients. Moreover, RET fusion differentially regulated PI3K and MAPK signaling depending on the partner gene of RET, and HOXD9 was found to be a recurrence marker for advanced papillary thyroid cancer.