Molecular control of endurance training adaptation in male mouse skeletal muscle.

Skeletal muscle has an enormous plastic potential to adapt to various external and internal perturbations. Although morphological changes in endurance-trained muscles are well described, the molecular underpinnings of training adaptation are poorly understood. We therefore aimed to elucidate the molecular signature of muscles of trained male mice and unravel the training status-dependent responses to an acute bout of exercise. Our results reveal that, even though at baseline an unexpectedly low number of genes define the trained muscle, training status substantially affects the transcriptional response to an acute challenge, both quantitatively and qualitatively, in part associated with epigenetic modifications. Finally, transiently activated factors such as the peroxisome proliferator-activated receptor-γ coactivator 1α are indispensable for normal training adaptation. Together, these results provide a molecular framework of the temporal and training status-dependent exercise response that underpins muscle plasticity in training.

Alterations in the hepatocyte epigenetic landscape in steatosis.

Fatty liver disease or the accumulation of fat in the liver, has been reported to affect the global population. This comes with an increased risk for the development of fibrosis, cirrhosis, and hepatocellular carcinoma. Yet, little is known about the effects of a diet containing high fat and alcohol towards epigenetic aging, with respect to changes in transcriptional and epigenomic profiles. In this study, we took up a multi-omics approach and integrated gene expression, methylation signals, and chromatin signals to study the epigenomic effects of a high-fat and alcohol-containing diet on mouse hepatocytes. We identified four relevant gene network clusters that were associated with relevant pathways that promote steatosis. Using a machine learning approach, we predict specific transcription factors that might be responsible to modulate the functionally relevant clusters. Finally, we discover four additional CpG loci and validate aging-related differential CpG methylation. Differential CpG methylation linked to aging showed minimal overlap with altered methylation in steatosis.

Increased chromatin accessibility facilitates intron retention in specific cell differentiation states.

Dynamic intron retention (IR) in vertebrate cells is of widespread biological importance. Aberrant IR is associated with numerous human diseases including several cancers. Despite consistent reports demonstrating that intrinsic sequence features can help introns evade splicing, conflicting findings about cell type- or condition-specific IR regulation by trans-regulatory and epigenetic mechanisms demand an unbiased and systematic analysis of IR in a controlled experimental setting. We integrated matched mRNA sequencing (mRNA-Seq), whole-genome bisulfite sequencing (WGBS), nucleosome occupancy methylome sequencing (NOMe-Seq) and chromatin immunoprecipitation sequencing (ChIP-Seq) data from primary human myeloid and lymphoid cells. Using these multi-omics data and machine learning, we trained two complementary models to determine the role of epigenetic factors in the regulation of IR in cells of the innate immune system. We show that increased chromatin accessibility, as revealed by nucleosome-free regions, contributes substantially to the retention of introns in a cell-specific manner. We also confirm that intrinsic characteristics of introns are key for them to evade splicing. This study suggests an important role for chromatin architecture in IR regulation. With an increasing appreciation that pathogenic alterations are linked to RNA processing, our findings may provide useful insights for the development of novel therapeutic approaches that target aberrant splicing.

A comprehensive approach for genome-wide efficiency profiling of DNA modifying enzymes.

A precise understanding of DNA methylation dynamics is of great importance for a variety of biological processes including cellular reprogramming and differentiation. To date, complex integration of multiple and distinct genome-wide datasets is required to realize this task. We present GwEEP (genome-wide epigenetic efficiency profiling) a versatile approach to infer dynamic efficiencies of DNA modifying enzymes. GwEEP relies on genome-wide hairpin datasets, which are translated by a hidden Markov model into quantitative enzyme efficiencies with reported confidence around the estimates. GwEEP predicts de novo and maintenance methylation efficiencies of Dnmts and furthermore the hydroxylation efficiency of Tets. Its design also allows capturing further oxidation processes given available data. We show that GwEEP predicts accurately the epigenetic changes of ESCs following a Serum-to-2i shift and applied to Tet TKO cells confirms the hypothesized mutual interference between Dnmts and Tets.

Resident memory CD4+ T lymphocytes mobilize from bone marrow to contribute to a systemic secondary immune reaction.

Resident memory T lymphocytes (TRM ) of epithelial tissues and the Bm protect their host tissue. To what extent these cells are mobilized and contribute to systemic immune reactions is less clear. Here, we show that in secondary immune reactions to the measles-mumps-rubella (MMR) vaccine, CD4+ TRM are mobilized into the blood within 16 to 48 h after immunization in humans. This mobilization of TRM is cognate: TRM recognizing other antigens are not mobilized, unless they cross-react with the vaccine. We also demonstrate through methylome analyses that TRM are mobilized from the Bm. These mobilized cells make significant contribution to the systemic immune reaction, as evidenced by their T-cell receptor Vß clonotypes represented among the newly generated circulating memory T-cells, 14 days after vaccination. Thus, TRM of the Bm confer not only local, but also systemic immune memory.

Epigenomic and transcriptional profiling identifies impaired glyoxylate detoxification in NAFLD as a risk factor for hyperoxaluria.

Epigenetic modifications (e.g. DNA methylation) in NAFLD and their contribution to disease progression and extrahepatic complications are poorly explored. Here, we use an integrated epigenome and transcriptome analysis of mouse NAFLD hepatocytes and identify alterations in glyoxylate metabolism, a pathway relevant in kidney damage via oxalate release-a harmful waste product and kidney stone-promoting factor. Downregulation and hypermethylation of alanine-glyoxylate aminotransferase (Agxt), which detoxifies glyoxylate, preventing excessive oxalate accumulation, is accompanied by increased oxalate formation after metabolism of the precursor hydroxyproline. Viral-mediated Agxt transfer or inhibiting hydroxyproline catabolism rescues excessive oxalate release. In human steatotic hepatocytes, AGXT is also downregulated and hypermethylated, and in NAFLD adolescents, steatosis severity correlates with urinary oxalate excretion. Thus, this work identifies a reduced capacity of the steatotic liver to detoxify glyoxylate, triggering elevated oxalate, and provides a mechanistic explanation for the increased risk of kidney stones and chronic kidney disease in NAFLD patients.

Differential degradation of RNA species by autophagy-related pathways in Arabidopsis.

The plant vacuole recycles proteins and RNA delivered to it by autophagy. In this study, by isolating intact vacuoles from Arabidopsis plants, followed by subsequent RNA purification, and deep sequencing, we provide a comprehensive characterization of Arabidopsis vacuolar RNAome. In the vacuolar RNAome, we detected ribosomal RNAs, transfer RNAs, including those of chloroplast origin, and in addition small RNA types. As autophagy is a main mechanism for the transport of RNA to the vacuole, atg5-1 mutants deficient in autophagy were included in our analysis. We observed severely reduced amounts of most chloroplast-derived RNA species in these mutants. Comparisons with cellular RNA composition provided an indication of possible up-regulation of alternative RNA breakdown pathways. By contrast, vacuolar RNA processing and composition in plants lacking vacuolar ribonuclease 2, involved in cellular RNA homeostasis, only showed minor alterations, possibly because of the presence of further so far unknown vacuolar RNase species. Among the small RNA types, we detected mature miRNAs in all vacuolar preparations but at much lower frequency in atg5-1, raising the possibility of a biological role for vacuolar miRNAs.

Lasp1 regulates adherens junction dynamics and fibroblast transformation in destructive arthritis.

The LIM and SH3 domain protein 1 (Lasp1) was originally cloned from metastatic breast cancer and characterised as an adaptor molecule associated with tumourigenesis and cancer cell invasion. However, the regulation of Lasp1 and its function in the aggressive transformation of cells is unclear. Here we use integrative epigenomic profiling of invasive fibroblast-like synoviocytes (FLS) from patients with rheumatoid arthritis (RA) and from mouse models of the disease, to identify Lasp1 as an epigenomically co-modified region in chronic inflammatory arthritis and a functionally important binding partner of the Cadherin-11/ß-Catenin complex in zipper-like cell-to-cell contacts. In vitro, loss or blocking of Lasp1 alters pathological tissue formation, migratory behaviour and platelet-derived growth factor response of arthritic FLS. In arthritic human TNF transgenic mice, deletion of Lasp1 reduces arthritic joint destruction. Therefore, we show a function of Lasp1 in cellular junction formation and inflammatory tissue remodelling and identify Lasp1 as a potential target for treating inflammatory joint disorders associated with aggressive cellular transformation.

Dysregulation of COVID-19 related gene expression in the COPD lung.

BACKGROUND: Chronic obstructive pulmonary disease (COPD) patients are at increased risk of poor outcome from Coronavirus disease (COVID-19). Early data suggest elevated Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2) receptor angiotensin converting enzyme 2 (ACE2) expression, but relationships to disease phenotype and downstream regulators of inflammation in the Renin-Angiotensin system (RAS) are unknown. We aimed to determine the relationship between RAS gene expression relevant to SARS-CoV-2 infection in the lung with disease characteristics in COPD, and the regulation of newly identified SARS-CoV-2 receptors and spike-cleaving proteases, important for SARS-CoV-2 infection. METHODS: We quantified gene expression using RNA sequencing of epithelial brushings and bronchial biopsies from 31 COPD and 37 control subjects. RESULTS: ACE2 gene expression (log2-fold change (FC)) was increased in COPD compared to ex-smoking (HV-ES) controls in epithelial brushings (0.25, p = 0.042) and bronchial biopsies (0.23, p = 0.050), and correlated with worse lung function (r = - 0.28, p = 0.0090). ACE2 was further increased in frequent exacerbators compared to infrequent exacerbators (0.51, p = 0.00045) and associated with use of ACE inhibitors (ACEi) (0.50, p = 0.0034), having cardiovascular disease (0.23, p = 0.048) or hypertension (0.34, p = 0.0089), and inhaled corticosteroid use in COPD subjects in bronchial biopsies (0.33, p = 0.049). Angiotensin II receptor type (AGTR)1 and 2 expression was decreased in COPD bronchial biopsies compared to HV-ES controls with log2FC of -0.26 (p = 0.033) and - 0.40, (p = 0.0010), respectively. However, the AGTR1:2 ratio was increased in COPD subjects compared with HV-ES controls, log2FC of 0.57 (p = 0.0051). Basigin, a newly identified potential SARS-CoV-2 receptor was also upregulated in both brushes, log2FC of 0.17 (p = 0.0040), and bronchial biopsies, (log2FC of 0.18 (p = 0.017), in COPD vs HV-ES. Transmembrane protease, serine (TMPRSS)2 was not differentially regulated between control and COPD. However, various other spike-cleaving proteases were, including TMPRSS4 and Cathepsin B, in both epithelial brushes (log2FC of 0.25 (p = 0.0012) and log2FC of 0.56 (p = 5.49E-06), respectively) and bronchial biopsies (log2FC of 0.49 (p = 0.00021) and log2FC of 0.246 (p = 0.028), respectively). CONCLUSION: This study identifies key differences in expression of genes related to susceptibility and aetiology of COVID-19 within the COPD lung. Further studies to understand the impact on clinical course of disease are now required.

Epigenetic control of region-specific transcriptional programs in mouse cerebellar and cortical astrocytes.

Astrocytes from the cerebral cortex (CTX) and cerebellum (CB) share basic molecular programs, but also form distinct spatial and functional subtypes. The regulatory epigenetic layers controlling such regional diversity have not been comprehensively investigated so far. Here, we present an integrated epigenome analysis of methylomes, open chromatin, and transcriptomes of astroglia populations isolated from the cortex or cerebellum of young adult mice. Besides a basic overall similarity in their epigenomic programs, cortical astrocytes and cerebellar astrocytes exhibit substantial differences in their overall open chromatin structure and in gene-specific DNA methylation. Regional epigenetic differences are linked to differences in transcriptional programs encompassing genes of region-specific transcription factor networks centered around Lhx2/Foxg1 in CTX astrocytes and the Zic/Irx families in CB astrocytes. The distinct epigenetic signatures around these transcription factor networks point to a complex interconnected and combinatorial regulation of region-specific transcriptomes. These findings suggest that key transcription factors, previously linked to temporal, regional, and spatial control of neurogenesis, also form combinatorial networks important for astrocytes. Our study provides a valuable resource for the molecular basis of regional astrocyte identity and physiology.

A gene therapy for inherited blindness using dCas9-VPR-mediated transcriptional activation.

Catalytically inactive dCas9 fused to transcriptional activators (dCas9-VPR) enables activation of silent genes. Many disease genes have counterparts, which serve similar functions but are expressed in distinct cell types. One attractive option to compensate for the missing function of a defective gene could be to transcriptionally activate its functionally equivalent counterpart via dCas9-VPR. Key challenges of this approach include the delivery of dCas9-VPR, activation efficiency, long-term expression of the target gene, and adverse effects in vivo. Using dual adeno-associated viral vectors expressing split dCas9-VPR, we show efficient transcriptional activation and long-term expression of cone photoreceptor-specific M-opsin (Opn1mw) in a rhodopsin-deficient mouse model for retinitis pigmentosa. One year after treatment, this approach yields improved retinal function and attenuated retinal degeneration with no apparent adverse effects. Our study demonstrates that dCas9-VPR-mediated transcriptional activation of functionally equivalent genes has great potential for the treatment of genetic disorders.

Neuronal deficiency of p38α-MAPK ameliorates symptoms and pathology of APP or Tau-transgenic Alzheimer's mouse models.

Alzheimer's disease (AD) is the leading cause of dementia with very limited therapeutic options. Amyloid ß (Aß) and phosphorylated Tau (p-Tau) are key pathogenic molecules in AD. P38α-MAPK is specifically activated in AD lesion sites. However, its effects on AD pathogenesis, especially on p-Tau-associated brain pathology, and the underlying molecular mechanisms remain unclear. We mated human APP-transgenic mice and human P301S Tau-transgenic mice with mapk14-floxed and neuron-specific Cre-knock-in mice. We observed that deletion of p38α-MAPK specifically in neurons improves the cognitive function of both 9-month-old APP and Tau-transgenic AD mice, which is associated with decreased Aß and p-Tau load in the brain. We further used next-generation sequencing to analyze the gene transcription in brains of p38α-MAPK deficient and wild-type APP-transgenic mice, which indicated that deletion of p38α-MAPK regulates the transcription of calcium homeostasis-related genes, especially downregulates the expression of grin2a, a gene encoding NMDAR subunit NR2A. Cell culture experiments further verified that deletion of p38α-MAPK inhibits NMDA-triggered calcium influx and neuronal apoptosis. Our systemic studies of AD pathogenic mechanisms using both APP- and Tau-transgenic mice suggested that deletion of neuronal p38α-MAPK attenuates AD-associated brain pathology and protects neurons in AD pathogenesis. This study supports p38α-MAPK as a novel target for AD therapy.

Feeding exogenous dsRNA interferes with endogenous sRNA accumulation in Paramecium.

Supply of exogenous dsRNA (exo-dsRNA), either by injection or by feeding, is a fast and powerful alternative to classical knockout studies. The biotechnical potential of feeding techniques is evident from the numerous studies focusing on oral administration of dsRNA to control pests and viral infection in crops/animal farming. We aimed to dissect the direct and indirect effects of exo-dsRNA feeding on the endogenous short interfering RNA (endo-siRNA) populations of the free-living ciliate Paramecium. We introduced dsRNA fragments against Dicer1 (DCR1), involved in RNA interference (RNAi) against exo- and few endo-siRNAs, and an RNAi unrelated gene, ND169. Any feeding, even the control dsRNA, diminishes genome wide the accumulation of endo-siRNAs and mRNAs. This cannot be explained by direct off-target effects and suggests mechanistic overlaps of the exo- and endo-RNAi mechanisms. Nevertheless, we observe a stronger down-regulation of mRNAs in DCR1 feeding compared with ND169 knockdown. This is likely due to the direct involvement of DCR1 in endo-siRNA accumulation. We further observed a cis-regulatory effect on mRNAs that overlap with phased endo-siRNAs. This interference of exo-dsRNA with endo-siRNAs warrants further investigations into secondary effects in target species/consumers, risk assessment of dsRNA feeding applications, and environmental pollution with dsRNA.

Targeted De-Methylation of the FOXP3-TSDR Is Sufficient to Induce Physiological FOXP3 Expression but Not a Functional Treg Phenotype.

CD4+ regulatory T cells (Tregs) are key mediators of immunological tolerance and promising effector cells for immuno-suppressive adoptive cellular therapy to fight autoimmunity and chronic inflammation. Their functional stability is critical for their clinical utility and has been correlated to the demethylated state of the TSDR/CNS2 enhancer element in the Treg lineage transcription factor FOXP3. However, proof for a causal contribution of the TSDR de-methylation to FOXP3 stability and Treg induction is so far lacking. We here established a powerful transient-transfection CRISPR-Cas9-based epigenetic editing method for the selective de-methylation of the TSDR within the endogenous chromatin environment of a living cell. The induced de-methylated state was stable over weeks in clonal T cell proliferation cultures even after expression of the editing complex had ceased. Epigenetic editing of the TSDR resulted in FOXP3 expression, even in its physiological isoform distribution, proving a causal role for the de-methylated TSDR in FOXP3 regulation. However, successful FOXP3 induction was not associated with a switch towards a functional Treg phenotype, in contrast to what has been reported from FOXP3 overexpression approaches. Thus, TSDR de-methylation is required, but not sufficient for a stable Treg phenotype induction. Therefore, targeted demethylation of the TSDR may be a critical addition to published in vitro Treg induction protocols which so far lack FOXP3 stability.

Identification and Characterization of Post-activated B Cells in Systemic Autoimmune Diseases.

Autoimmune diseases (AID) such as systemic lupus erythematosus (SLE), primary Sjögren's syndrome (pSS), and rheumatoid arthritis (RA) are chronic inflammatory diseases in which abnormalities of B cell function play a central role. Although it is widely accepted that autoimmune B cells are hyperactive in vivo, a full understanding of their functional status in AID has not been delineated. Here, we present a detailed analysis of the functional capabilities of AID B cells and dissect the mechanisms underlying altered B cell function. Upon BCR activation, decreased spleen tyrosine kinase (Syk) and Bruton's tyrosine kinase (Btk) phosphorylation was noted in AID memory B cells combined with constitutive co-localization of CD22 and protein tyrosine phosphatase (PTP) non-receptor type 6 (SHP-1) along with hyporesponsiveness to TLR9 signaling, a Syk-dependent response. Similar BCR hyporesponsiveness was also noted specifically in SLE CD27- B cells together with increased PTP activities and increased transcripts for PTPN2, PTPN11, PTPN22, PTPRC, and PTPRO in SLE B cells. Additional studies revealed that repetitive BCR stimulation of normal B cells can induce BCR hyporesponsiveness and that tissue-resident memory B cells from AID patients also exhibited decreased responsiveness immediately ex vivo, suggesting that the hyporesponsive status can be acquired by repeated exposure to autoantigen(s) in vivo. Functional studies to overcome B cell hyporesponsiveness revealed that CD40 co-stimulation increased BCR signaling, induced proliferation, and downregulated PTP expression (PTPN2, PTPN22, and receptor-type PTPs). The data support the conclusion that hyporesponsiveness of AID and especially SLE B cells results from chronic in vivo stimulation through the BCR without T cell help mediated by CD40-CD154 interaction and is manifested by decreased phosphorylation of BCR-related proximal signaling molecules and increased PTPs. The hyporesponsiveness of AID B cells is similar to a form of functional anergy.

Unique and assay specific features of NOMe-, ATAC- and DNase I-seq data.

Chromatin accessibility maps are important for the functional interpretation of the genome. Here, we systematically analysed assay specific differences between DNase I-seq, ATAC-seq and NOMe-seq in a side by side experimental and bioinformatic setup. We observe that most prominent nucleosome depleted regions (NDRs, e.g. in promoters) are roboustly called by all three or at least two assays. However, we also find a high proportion of assay specific NDRs that are often 'called' by only one of the assays. We show evidence that these assay specific NDRs are indeed genuine open chromatin sites and contribute important information for accurate gene expression prediction. While technically ATAC-seq and DNase I-seq provide a superb high NDR calling rate for relatively low sequencing costs in comparison to NOMe-seq, NOMe-seq singles out for its genome-wide coverage allowing to not only detect NDRs but also endogenous DNA methylation and as we show here genome wide segmentation into heterochromatic B domains and local phasing of nucleosomes outside of NDRs. In summary, our comparisons strongly suggest to consider assay specific differences for the experimental design and for generalized and comparative functional interpretations.

Genome-Wide Analysis of the Nucleosome Landscape in Individuals with Coffin-Siris Syndrome.

The switch/sucrose non-fermenting (SWI/SNF) complex is an ATP-dependent chromatin remodeller that regulates the spacing of nucleosomes and thereby controls gene expression. Heterozygous mutations in genes encoding subunits of the SWI/SNF complex have been reported in individuals with Coffin-Siris syndrome (CSS), with the majority of the mutations in ARID1B. CSS is a rare congenital disorder characterized by facial dysmorphisms, digital anomalies, and variable intellectual disability. We hypothesized that mutations in genes encoding subunits of the ubiquitously expressed SWI/SNF complex may lead to alterations of the nucleosome profiles in different cell types. We performed the first study on CSS-patient samples and investigated the nucleosome landscapes of cell-free DNA (cfDNA) isolated from blood plasma by whole-genome sequencing. In addition, we studied the nucleosome landscapes of CD14+ monocytes from CSS-affected individuals by nucleosome occupancy and methylome-sequencing (NOMe-seq) as well as their expression profiles. In cfDNA of CSS-affected individuals with heterozygous ARID1B mutations, we did not observe major changes in the nucleosome profile around transcription start sites. In CD14+ monocytes, we found few genomic regions with different nucleosome occupancy when compared to controls. RNA-seq analysis of CD14+ monocytes of these individuals detected only few differentially expressed genes, which were not in proximity to any of the identified differential nucleosome-depleted regions. In conclusion, we show that heterozygous mutations in the human SWI/SNF subunit ARID1B do not have a major impact on the nucleosome landscape or gene expression in blood cells. This might be due to functional redundancy, cell-type specificity, or alternative functions of ARID1B.

Transcriptomics of a KDELR1 knockout cell line reveals modulated cell adhesion properties.

KDEL receptors (KDELRs) represent transmembrane proteins of the secretory pathway which regulate the retention of soluble ER-residents as well as retrograde and anterograde vesicle trafficking. In addition, KDELRs are involved in the regulation of cellular stress response and ECM degradation. For a deeper insight into KDELR1 specific functions, we characterised a KDELR1-KO cell line (HAP1) through whole transcriptome analysis by comparing KDELR1-KO cells with its respective HAP1 wild-type. Our data indicate more than 300 significantly and differentially expressed genes whose gene products are mainly involved in developmental processes such as cell adhesion and ECM composition, pointing out to severe cellular disorders due to a loss of KDELR1. Impaired adhesion capacity of KDELR1-KO cells was further demonstrated through in vitro adhesion assays, while collagen- and/or laminin-coating nearly doubled the adhesion property of KDELR1-KO cells compared to wild-type, confirming a transcriptional adaptation to improve or restore the cellular adhesion capability. Perturbations within the secretory pathway were verified by an increased secretion of ER-resident PDI and decreased cell viability under ER stress conditions, suggesting KDELR1-KO cells to be severely impaired in maintaining cellular homeostasis.

Transcriptome Kinetics of Saccharomyces cerevisiae in Response to Viral Killer Toxin K1.

The K1 A/B toxin secreted by virus-infected Saccharomyces cerevisiae strains kills sensitive cells via disturbance of cytoplasmic membrane functions. Despite decades of research, the mechanisms underlying K1 toxicity and immunity have not been elucidated yet. In a novel approach, this study aimed to characterize transcriptome changes in K1-treated sensitive yeast cells in a time-dependent manner. Global transcriptional profiling revealed substantial cellular adaptations in target cells resulting in 1,189 differentially expressed genes in total. Killer toxin K1 induced oxidative, cell wall and hyperosmotic stress responses as well as rapid down-regulation of transcription and translation. Essential pathways regulating energy metabolism were also significantly affected by the toxin. Remarkably, a futile cycle of the osmolytes trehalose and glycogen was identified probably representing a critical feature of K1 intoxication. In silico analysis suggested several transcription factors involved in toxin-triggered signal transduction. The identified transcriptome changes provide valuable hints to illuminate the still unknown molecular events leading to K1 toxicity and immunity implicating an evolutionarily conserved response at least initially counteracting ionophoric toxin action.

Exogenous RNAi mechanisms contribute to transcriptome adaptation by phased siRNA clusters in Paramecium.

Extensive research has characterized distinct exogenous RNAi pathways interfering in gene expression during vegetative growth of the unicellular model ciliate Paramecium. However, role of RNAi in endogenous transcriptome regulation, and environmental adaptation is unknown. Here, we describe the first genome-wide profiling of endogenous sRNAs in context of different transcriptomic states (serotypes). We developed a pipeline to identify, and characterize 2602 siRNA producing clusters (SRCs). Our data show no evidence that SRCs produce miRNAs, and in contrast to other species, no preference for strand specificity of siRNAs. Interestingly, most SRCs overlap coding genes and a separate group show siRNA phasing along the entire open reading frame, suggesting that the mRNA transcript serves as a source for siRNAs. Integrative analysis of siRNA abundance and gene expression levels revealed surprisingly that mRNA and siRNA show negative as well as positive associations. Two RNA-dependent RNA Polymerase mutants, RDR1 and RDR2, show a drastic loss of siRNAs especially in phased SRCs accompanied with increased mRNA levels. Importantly, most SRCs depend on both RDRs, reminiscent to primary siRNAs in the RNAi against exogenous RNA, indicating mechanistic overlaps between exogenous and endogenous RNAi contributing to flexible transcriptome adaptation.