Abstract
Mycoses induced by C.albicans or A.fumigatus can cause important host damage either by deficient or exaggerated immune response. Regulation of chemokine and cytokine signaling plays a crucial role for an adequate inflammation, which can be modulated by vitamins A and D. Non-coding RNAs (ncRNAs) as transcription factors or cis-acting antisense RNAs are known to be involved in gene regulation. However, the processes during fungal infections and treatment with vitamins in terms of therapeutic impact are unknown. We show that in monocytes both vitamins regulate ncRNAs involved in amino acid metabolism and immune system processes using comprehensive RNA-Seq analyses. Compared to protein-coding genes, fungi and bacteria induced an expression change in relatively few ncRNAs, but with massive fold changes of up to 4000. We defined the landscape of long-ncRNAs (lncRNAs) in response to pathogens and observed variation in the isoforms composition for several lncRNA following infection and vitamin treatment. Most of the involved antisense RNAs are regulated and positively correlated with their sense protein-coding genes. We investigated lncRNAs with stimulus specific immunomodulatory activity as potential marker genes: LINC00595, SBF2-AS1 (A.fumigatus) and RP11-588G21.2, RP11-394l13.1 (C.albicans) might be detectable in the early phase of infection and serve as therapeutic targets in the future.
Similar content being viewed by others
Introduction
Non-coding RNAs (ncRNAs) play essential roles in regulating protein-coding genes and gained an enormous medical relevance in the last decades. Although only less than 3% of the human genome codes for proteins1,2, it is believed that at least 80% of the human genome has a function2. Among these regions, which are mainly transcribed, the ncRNAs have been shown to regulate several basic cellular processes, such as transcription, processing, and translation. They are involved in chromosome structure, DNA replication, gene regulation, genome defense and protein transport. During the last decade, miRNAs received major attention as regulators of various cell processes. However, currently there are 2,468 ncRNA families described (Rfam v.12.0)3, of which 750 are known to exist in human, covering 1,853 genomic regions (http://rfam.sanger.ac.uk). Additionally, the number of long non-coding RNAs (lncRNAs, >200ânt, mostly with exon-intron structure) in humans is estimated to be 15,9414. They can be distinguished by their loci relative to protein annotations: intergenic, sense-intronic, sense-overlapping and antisense.
Vitamins A and D have been shown to play a regulatory role in immunity, prevent cancer development and can be used as therapeutical option in different malignancies5,6,7,8. However, the mechanisms have not been described so far. Although various publications report ncRNAs to be differentially expressed during infection diseases in human9,10,11, there is no knowledge about the role of ncRNAs during fungal infection and only a few studies have been reported about treatment with vitamins A and D. In the following, we summarize ncRNAs which have been brought in relation to vitamin supply or fungal infection.
NcRNAs and atRA
All-trans-retinoic acid (atRA), the active metabolite of vitamin A, significantly regulates protein-coding genes and to an even higher fraction ncRNAs in human nCCIT cells12. In addition to micro RNAs (miRNAs) being regulated by atRA13,14, several lncRNAs have been reported to be influenced by atRA in neuronal and cancer cells. In rats, atRA is involved in hippocampal functions by up-regulating LINC_RBE15 and is able to restore NEAT1 expression in NB4 leukemia cell differentiation16. In the same cells and APL cell lines, a positive correlation of atRA and HOTAIRM1 has been shown. When HOTAIRM1 is regulated by atRA, less cell cycle genes are expressed, and genes involved in DNA replication as well as the expression of HOXD13 and HOXA11 are influenced17,18,19. Retinoic acid activates HOXD-AS1 (inducing the PI3K/Akt pathway) and the apoptosis repressor HOXA-AS2 in peripheral blood neutrophils20,21. Moreover, it has been shown that a reduced level of atRA leads to an induction of leukocyte activation, maturation and defense response19.
AtRA is also known to significantly down-regulate Malat1 expression during the differentiation of EML cells22. Furthermore atRA regulates LINC-RSAS in mammalian brain, which is associated with aging processes23 and lincRNA-EPS, which promotes terminal differentiation and inhibits apoptosis of mature erythrocytes18. In contrast to lncRNAs being regulated by various conditions, the retinoic acid signaling has been shown to be controlled by the lncRNAs LL18/NANCI and LL34, during endoderm development24.
NcRNAs and vitamin D
Several lncRNAs have been reported to be affected by vitamin D in different cell lines. In mouse keratinocytes, a knockout of the the vitamin D receptor (VDR) leads to a down-regulation of tumor suppressor lncRNAs (Kcnq1ot1, lincRNA-p21, Foxn2-as, Gtl2-as, H10-as) and an increase in the expression of oncogenes (H19, HOTTIP, Nespas, mHOTAIR, Malat1 and interacting partner SRA)25,26.
VDR directly down-regulates miR-181a, co-regulates mir-106b, and finally leads to an up-regulation of p27. It targets mir-627 that in turn targets JMJD1A, miR-98 and let-7a-2. MiR-125b inhibits VDR expression and in turn the receptor down-regulates miR-125b27.
NcRNAs during fungal infection
Although the knowledge regarding the regulatory impact of ncRNAs during fungal infection is limited, in Arabidopsis thaliana a few lncRNAs and lincRNAs have been shown to play an important role in antifungal immune response28. Also multiple miRNAs have been shown to be key elements in the respiratory epithelium during Candida infection29. The role of ncRNAs in humans during fungal infection and the influence of vitamins has not been described so far. On protein-coding level only few studies have been reported30,31.
As mentioned above, Vitamins A and D have been shown to play a regulatory role in immunity and preventing cancer where they have been used as therapeutic option5,6,7,8. However, the mechanisms have not been described so far. Here, we raise the hypothesis that ncRNAs might contribute to this effect.
In this paper, we find that vitamins A and D could directly regulate ncRNAs involved in immunity against fungal infections. We explore the effect of different pathogens on the expression of lncRNA in human monocytes. We describe specific effects during fungal infections when compared to bacterial infection with E. coli, which was used as outgroup in this study. Moreover, we explore the impact of atRA and vitamin D administration to obtain a systematic and comprehensive overview of vitamin regulated ncRNAs during infection. Many lncRNAs have been shown to be regulators of inflammatory chemokine and cytokine signaling. Using RNA-Seq we find both vitamins to regulate lncRNAs involved in amino acid metabolism and immune system processes, such as pro-inflammation and apoptosis. We present response specific and treatment specific results on different levels, for example single genes which are potential biomarkers and show structurally conserved regions and alternative splicing events.
Materials and Methods
Preparation of fungi and bacteria
Overnight cultures of C. albicans (SC5314) in YPD medium were washed three times with PBS and resuspended at 108 yeasts/ml in RPMI 1640 GlutaMAX medium (Gibco, UK) supplemented with 10% fetal bovine serum (FBS; Biochrom, Germany).
A. fumigatus (AF293) was grown on AMM plates at 30â°C for 6 d. Conidiospores were harvested by rinsing the plates with water +0.05% Tween-20 (Sigma-Aldrich, Germany) and filtered through 70-μm and 30-μm pre-separation filters (Miltenyi Biotec, UK) to obtain a single-cell suspension. Conidia were then washed twice in PBS and resuspended at 107 conidia/ml in RPMI 1640 GlutaMAX medium supplemented with 10% FBS (Biochrom, Germany). Germlings were obtained by incubation of conidia at 37â°C under continuous shaking for 6â8âh. They were then centrifuged and resuspended at 108 cells/ml in fresh RPMI 1640 GlutaMAX medium supplemented with 10% FBS.
Overnight culture of Escherichia coli (isolate 018:K1:H7) in LB medium was washed three times in PBS and resuspended in RPMI 1640 GlutaMAX medium supplemented with 10% FBS. The concentration of bacteria was adjusted to 109âcfu/ml. All pathogens were heat-killed by incubation at 65â°C for 30âmin and immediately used for stimulation assays.
Monocyte isolation
Human monocytes were isolated from 500âml fresh whole blood (drawn within 1âh before use) of healthy male donors. Blood was layered onto an equal volume of 1-Step Polymorphs (Accurate Chemical & Scientific Corporation, USA) and centrifuged at 650âÃâg for 35âmin. After centrifugation, the peripheral blood mononuclear cells (PBMCs) were collected, and normal osmolarity was restored by adding an equal volume of 0.45% cold NaCl. After erythrocyte lysis using a hypotonic buffer, cells were washed twice in cold PBS and counted using a Neubauer chamber. Cell viability of >95% was assessed by trypan blue staining. Monocytes were isolated from the PBMCs using the monocyte isolation kit II and quadro-MACS (Miltenyi Biotec, UK), following manufacturerâs instructions. Purity of the monocytes (>90%) was assessed as described in Klassert et al.32.
Stimulation assays
Monocytes were resuspended at 5âÃâ106 cells/ml in RPMI 1640 GlutaMAX medium (Gibco, UK) supplemented with 10% FBS (Biochrom, Germany) and 1% Penicillin/Streptomycin (Thermo Fisher Scientific, USA). They were seeded on 6-well plates (VWR International, Germany) and allowed to equilibrate at 37â°C and 5% CO2 for 2âh. Cells were then pre-incubated with 1âμM atRA (Sigma-Aldrich (Germany) or 1α,25 (OH)2D3 (Sigma-Aldrich (Germany) for 30âmin. Then, the heat-killed pathogens were added at a pathogen:host ratio of 1:1 for C. albicans yeast and A. fumigatus germ tubes, and 10:1 in case of E. coli stimulation. After 6âh of incubation at 37â°C and 5% CO2, cell viability >90% was assessed by trypan blue staining, and the monocytes were harvested for RNA isolation. The whole experimental workflow is depicted in Supplementary Fig. S1 by Klassert et al.32.
In total, we had four different immune-stimulatory settings (w/o infection, A. fumigatus infection, C. albicans infection and E. coli infection), in each of which we aimed to address the effect of vitamin A (atRA) or vitamin D supplementation.
RNA sequencing
RNA was isolated from 5âÃâ106 monocytes using the RNeasy Mini Kit (Qiagen, Germany). An additional step was included to remove the residual genomic DNA using DNaseI (Qiagen, Germany). Quantity and quality of total RNA was assessed using a Nanodrop ND-1000 spectrophotometer (Thermo Fisher Scientific, USA) and a Tape Station 2200 (Agilent Technologies, USA), respectively. Strand-specific RNA libraries were constructed using the Dynabeads mRNA DIRECT Micro Purification Kit and the Ion Total RNA-Seq Kit v2.0 (Thermo Fisher Scientific, USA). Library templates were clonally amplified on Ion Sphere particles using the Ion PI Hi-Q Chef Kit and Ion Chef instrument (Thermo Fisher Scientific, USA), loaded onto Ion PI Chips and sequenced on an Ion Proton Sequencer (Thermo Fisher Scientific, USA). For sequencing, in total 48 samples were multiplexed on 12 chips. The raw sequence data in fastq format are stored in the Sequence Read Archive (SRA) at National Center for Biotechnology Information (NCBI) and can be accessed at NCBI homepage (accession number: SRP076532).
Read preprocessing and rRNA filtering
Raw reads in fastq format were quality controlled with FastQC33 (v0.11.3) and trimmed with a window size of 10 using Prinseq34 (v0.20.3) (Qââ¥â20, lengthââ¥â20). The resulting reads were aligned to reference rRNA databases using SortMeRNA35 (v2.0). After quality control, up to 19 242 435 reads with a length of 20â373âbp and a GC content of ~50% were obtained for each sample (see Data availability section for external supplement S1).
Mapping
The quality controlled and rRNA cleaned cleaned reads were aligned to reference genomes using Segemehl36 (v0.2.0). This mapping was performed against the human genome version GRCh38, downloaded from Ensembl (release 80). Indices were build together with the corresponding pathogen genomes, depending on the samples to be mapped. For E. coli, the complete genome of strain K-12 substr. MG1655 (NC_000913.3) was downloaded from NCBI. Genomes of A. fumigatus and C. albicans were obtained from aspergillusgenome.org (21.05.2015) and candidagenome.org (Ca21, 28.05.2015). All mappings were performed with default parameters and the âsplits option of Segemehl to allow for multiple spliced read alignments.
Differential gene expression
We used HTSeq-Count37 (v0.6.0) to quantify strand-specific and unique mapped reads on exon level. As reference, the full human Ensembl annotation (version GRCh38.80), including both protein- and non-coding genes, was used. Overall, 60 604 annotated features, comprising 19 825 protein-coding and 24 150 non-coding genes, were included in the annotation.
The raw read counts and rRNA cleaned read counts from each sample were normalized based on the library size and tested for significantly differentially expressed genes (adjusted p-valueââ¥â0.05) using DESeq238 (v1.10.1).
Gene filtering
The DESeq2 results of different comparisons were additionally filtered by calculating TPM values for each gene in each sample as shown in ref. 39.
where ci is the raw read count of gene i, li is the cumulative exon length of gene i and N is the number of all genes in the given annotation.
For each gene, we calculated four mean TPM values (TPMM), based on the 12 samples corresponding to control or one of the three different infection types. For protein-coding genes, we used a TPMMâ>â5 to filter the final gene sets for further analyses. For ncRNAs, we used a minimum basemean read count of 10.
Identification of different RNA isoforms
In order to identify different isoforms of antisense RNAs and lncRNAs, we investigated the expression patterns of the individual exons of these RNAs. Therefore, we adjusted our DESeq2 pipeline to calculate differentially expressed exons instead of whole genes Results were filtered by an adjusted p-value (â¤0.05) and by basemean expression quantifier (â¥40).
Principal component analyses
Principal component analyses (PCA) were performed subsequent to the corresponding DESeq2 runs on selected gene subsets (e.g. only protein-coding genes or only lncRNAs) and with different variance cut-offs in R.
3-dimensional PCAs were plotted with the scatterplot3d package in R.
Identification of co-regulated genes
To compare the effect of atRA and vitamin D during different infections, log2 fold changes (FC) as computed by DEseq2 were visualized using scatter plots in R.
The scatterplots were overlaid with contour plots for a two-dimensional kernel estimate (kde2d; MASS package) using the default parameters. Outliers were labeled with the respective gene names.
Box plots of certain gene expression patterns were visualized with the help of ref. 40.
Heat maps
Gene-wise hierarchical clustering was performed on variance stabilized read counts to build heat maps of selected differentially expressed gene sets (adjusted p valueââ¤â0.05). The input matrices for each heat map were scaled on rows to visualize changes in expression on gene level.
Secondary structure prediction
The prediction of conserved secondary structures for potential marker genes is based on the multiz7way whole genome alignment, downloaded from UCSC (Octobre 2016) in multiple alignment format (MAF). Next to the human genome it contains the species Pan troglodytes, Macaca mulatta, Canis lupus familiaris, Mus musculus, Rattus norvegicus and Monodelphis domestica. According to the exons of human lncRNAs, alignments were sliced exon-wise and concatenated afterwards. To detect common secondary structures alignments were folded with RNAalifold (v2.1.8) from the Vienna RNA package41 following a window wise strategy (step width 50, window 300). Resulting structures were evaluated by hand.
Ethics statement
The blood of healthy male donors was drawn after written informed consent. This is in accordance with the Declaration of Helsinki, all protocols were approved by the Ethics Committee of the University Hospital Jena (permit number: 3639-12/12).
Data Availability
External supplement: http://www.rna.uni-jena.de/supplements/fungi_infection/.
Results
The role of ncRNAs in humans during fungal infection and the influence of vitamins has not been described so far. Here, we investigate the transcriptional response of human monocytes after 6âh of infection with A. fumigatus, C. albicans and E. coli and, in addition, the influence of atRA and vitamin D during infection using RNA-Seq. We compared the three infections with control samples, as well as the three infections among each other. Furthermore, we determined the effect of vitamins in each type of infection. We particularly focused on antisense RNAs as potential cis-regulatory elements of related sense located protein coding genes and lincRNAs. Antisense RNAs and lincRNAs represent 92.2% of annotated lncRNAs and 56.2% of annotated ncRNAs in the human genome.
Few highly expressed ncRNAs with strong effect in response to infection
We observed that about 5% of the lncRNAs and antisense RNAs were differentially expressed with high significance in at least one of the comparisons (see Table 1). In comparison to the 24% of highly significantly regulated protein-coding mRNAs32, the number of differentially expressed lncRNAs is low. However, these few ncRNAs are much stronger regulated than the protein-coding mRNAs (see Fig. 1; real fold changes up to 4000). Moreover, most of the strongly regulated ncRNAs are present in multiple related organisms with compensatory mutations to keep their potential secondary structures (see Data availability section for external supplement S10). This supports the hypothesis of a few strongly regulated ncRNAs to influence a large set of protein-coding genes.
Interestingly, we see far more up-regulated than down-regulated non-coding RNAs during infection (see Fig. 1, plots in first column and first row). See also Klassert et al.32, where the number of protein-coding genes being up- and down-regulated are comparable. The expression pattern of ncRNAs in human monocytes during C. albicans and E. coli infection are similar (see Fig. 1, plot 2.4 and 4.2, with generally lower fold changes), whereas A. fumigatus infection shows a less intensive impact in human monocytes (see Fig. 1, plots 2.3, 3.4, 3.2 and 4.3). In general, the influence of vitamin supply on protein-coding mRNAs and ncRNAs is not as large as for the pathogens itself (see Fig. 2). At a first sight, the pathogen dependent transcription profiles upon vitamin A supply are comparable to those observed upon vitamin D supply. Nevertheless when we address the specific impact of both vitamins in each stimulatory setting, we could identify different patterns of vitamin dependent transcriptional regulation (see Data availability section for SFig. 4).
Only few ncRNAs significantly expressed in response to atRA supply
Using principal component analysis (PCA) on the top variant 500 protein-coding and non-coding genes in the entire datasets, we were able to detect a strong impact of each pathogen. Nevertheless, the strong atRA specific effect observed in protein-coding mRNAs could not be observed for most of the ncRNA biotypes (see Fig. 3). This result is in line with Fig. 1, where we could observe only a modest regulation of ncRNAs when compared to protein-coding mRNAs during infection (see y-axis of Fig. 3A,B). For each subset of ncRNAs (lncRNAs Fig. 3C, sense RNAs Fig. 3D, and antisense RNAs Fig. 3E), we observe a similar pattern upon vitamin stimulation. We see only weak differences between control, vitamin D-treated and atRA-treated samples. We could only identify a slight atRA-dependent impact on antisense RNAs, which was consistent across all stimulatory settings (see Fig. 3E). This effect could be validated using a strict variance cutoff of â¥2 (see Fig. 3F), which reveals about 30 genes affected by atRA supply. About half of the 30 significantly expressed antisense ncRNAs are known to be related to immune response or infection and inflammation. Among them, we find AC073072.5 and RP11â536K7.5 to be antisense to the protein-based infection key player IL6 and receptor IL2RA, respectively32. Other antisense RNAs remain with unknown function.
With this analysis we are able to support the general assumption that a few regulatory ncRNAs are affected by treatment, but a tremendous number of protein-coding genes might be regulated by these few ncRNAs. One may also argue that only a few ncRNAs and an independent massive number of protein-coding genes are affected by atRA supply.
Antisense RNAs involved in basic cell pathways
The general function of antisense RNAs is unknown. However, they have been proposed to interact with the protein-coding mRNAs on the sensestrand42,43,44. Therefore, we used the potential sense mRNA-targets of all antisense RNAs, to explore their putative impact on known pathways. Expression data of antisense RNAs were artificially transferred to corresponding mRNAs. Pathway enrichment analysis of these protein-coding genes revealed several pathways which clustered together according to their ontology, thereby suggesting an impact of antisense RNAs on functions such as immune response, inflammation and apoptosis (see Fig. 4A). Interestingly, these biological functions were also highlighted by Klassert et al., when using the expression data of all protein-coding genes to perform gene ontology enrichment- and pathway analyses32.
We find antisense RNAs being generally down-regulated during infection with C. albicans and E. coli, but not for A. fumigatus. Protein-coding genes involved in amino acid metabolism mainly have antisense RNAs down-regulated in all tested infections, whereas antisense RNAs involved in other basic cellular processes seem generally up-regulated. For many antisense RNAs, we observe a surprising positive correlation to the expression levels of the sense mRNAs (see Fig. 5). Such a positive correlation of sense and antisense RNA transcription has been former reported45.
One of the significantly enriched regulated pathways from KEGG is the âcytokine-cytokine-receptor-interactionâ pathway46,47. The most significantly involved antisense RNAs are displayed in Fig. 4B. Generally, we see only a marginal effect during A. fumigatus infection, in contrast to the dramatic up-regulation of cytokines and their receptors during C. albicans and E. coli infection. Also on single gene level we confirm a more similar reaction of monocytes to C. albicans and E. coli than between the two fungal pathogens. AC073072.5 (antisense to IL6) is more than 100 times up-regulated during C. albicans infection. Also, IL12A-AS1 is up-regulated during C. albicans infection. However, this effect is reduced by atRA supply.
In general, protein-coding genes for cytokines and their receptors, as described in ref. 32, are highly up-regulated during infection and down-regulated by atRA supply. Interestingly most of their antisense RNAs show a similarly strong behavior as described above for IL6 and AC073072.5 (see Fig. 5), which has never been reported in previous studies. However, a few antisense RNAs show an inverse transcription pattern as their corresponding cytokines, e.g. IL2RBâRP1-151814.6, IL21RâIL21R-AS1, and TNFSF12âRP11-186B7.7.
Retinol metabolism modulated by atRA and vitamin D
Since atRA has been shown to regulate immune response32, it is highly notorious that atRA-metabolism can be found among the top regulated pathways in our dataset (see Data availability section for SFig. 7). Therefore we investigated this pathway in more detail. We observed that the administration of atRA and vitamin D leads to an up-regulation of the enzymes SDR dehydrogenase/reductase member 4 and 9 (DHRS4,DHRS) as well as retinol dehydrogenase (RDH) and alcohol dehydrogenase (ADH) which are all responsible for the conversion of Retinol into Retinal as initial step for the generation of atRA. Interestingly, in absence of atRA supply, we can observe that pathogens inhibit the transcription of these Retinol metabolizing enzymes (see Fig. 6). Also the antisense RNA DHRS4-AS1 is down-regulated by pathogens, whereas upon vitamin supply a slight up-regulation can be observed. Thus, we observe a direct correlation between the antisense RNA and the DHRS gene family members. In contrast, the antisense RNA RP11-696N14 shows no consistent pattern suggesting an independent function. We conclude that atRA is able to regulate its own metabolism at transcriptional level, regulating not only protein-coding genes, but also antisense RNAs with putative functions in this pathway. The regulation of the retinol metabolism constitutes an additional example of counter-active effects between vitamin and pathogen stimulations.
NcRNAs as potential marker genes
In our analysis, we find the ncRNAs known from the literature provided in the external supplement STable 3 (see Data availability section), to be involved in either fungal infections of plants or in human cancer development during atRA supply, to be only marginally affected. We find other previously annotated ncRNAs without an assigned function to be more relevant for fungal infection and upon vitamin supply. For a complete list of potential marker genes see Data availability section for external supplement STable 5). The function of novel ncRNAs is difficult to determine in silico. Secondary structures, with compensatory mutations and evolutionary conserved regions as potential protein interaction sites as provided in the external supplement S10 (see Data availability section). Genomically flanking (protein-coding) genes are described below to hint to a possible function.
Figure 7 depicts ncRNAs as potential marker genes for various scenarios. The following functional assumptions are taken from GeneVisible, if not indicated differently. A. fumigatus-specific ncRNAs LINC00595 and SBF2-AS1 are described to be involved in inflammation and immune response. The first gene is located within a ribosomal protein cluster and near POLR3A, possibly associated with the translation machinery. The latter gene has been described just recently to promote proliferation in non-small cell lung cancer (NSCLC)48.
We identified RP11-588G21.2 and RP11-394I13.1 to be C. albicans-infection specific genes and CTD-3128G10.6 and SRP to be potential marker genes for E. coli infection. RP11-588G21.2 has been reported to possibly being involved in the maturation process of mature immune cells49.
SRP has been found in T-cell exosomes, suggesting that their amount might be increased during infection50.
Among the different infection types, we identified RP11-13A1.1 and LINC00856 to be specific for fungal challenge. For infections in general, we propose the following marker genes: RP11-384O8.1, RP11-662I13.2, and LINC01181. None of them has been reported previously to be of interest.
The ncRNAs AC011899.9 and RP11-483P21.2 appear to be specific for atRA supply (up-regulated by atRA in all stimulatory settings), whereas the infection-specific transcription of RP11-662I13.2 and XXbac-BPG249D20.9 are inhibited by atRA (without an effect in control samples). AC011899.9 and RP11-483P21.2 have previously only been reported to be a marker for mature beta-cells (like RP11-588G21.2)49. The concept of genes being up-regulated during infection and then inhibited by atRA supply has been also observed with protein-coding genes32 and these ncRNAs seem to follow this line.
The potential marker genes for higher vitamin D levels, C20orf197 and RP11-401P9.4, have been previously reported to be involved in tumorogenesis51. As one of only a few genes, we see LINC01270 being up-regulated during infection. Vitamin D supply increases this effect. LINC00900 is significantly up-regulated by atRA and down-regulated by vitamin D.
Isoforms of antisense RNAs and lncRNAs during vitamin D supply
Retinoic acid receptors (RAR) and retinoid X receptors (RXR) have three isotypes α, β, and γ, which are highly conserved between human and mice, suggesting specific functions for each receptor type52,53. LncRNAs can consist of several exons and undergo alternative splicing events leading to different isoforms54,55,56,57. These transcript variations may be expressed under various conditions, fulfilling a variety of different functions, similar to protein isotypes. We identified several alternative splicing events for antisense RNAs and lncRNAs within our samples (see STable 8). As one example, we found the lncRNA AATBC to show different isoforms during E. coli infection and treatment with vitamin D or atRA (see Fig. 8). AATBC has been recently described to be involved in suppression of cell proliferation, as well as apoptosis induction in bladder cancer58. Interestingly, many additional lncRNA isoforms were identified upon vitamin D treatment, especially during E. coli infection (see Fig. 8).
Conclusion
Retinoids and vitamin D are known to play a regulatory role in development and physiology and are critical molecules in medical therapies. The mechanisms for this regulations are widely unknown. In this manuscript, we found lncRNAs to presumably contribute to their regulatory mechanisms. The number of differentially expressed lncRNAs is low compared to protein-coding mRNAs, but these lncRNAs are massively regulated (with fold changes higher than 4000). LncRNAs are regulated by atRA and vitamin D and are potentially involved in immunity, inflammation and apoptosis. The genes we have identified to be most important, and thus to be potential marker genes, have not been described previously and need further testing for validation of their putative function. We propose potential marker genes for infections in general and for specific fungal infections, which might be used as therapeutic targets in the future.
For the first time, the landscape of both the protein-coding mRNAs32 and non-coding RNAs was analyzed during fungal infection using high throughput sequencing of monocytes. Surprisingly, at ncRNA level, the immune responses during C. albicans and E. coli infection were more similar than during A. fumigatus infection. We revealed global systemic effects as well as detailed insights into the regulation of immune related proteins by lncRNAs, in particular antisense RNAs.
Although we did not find a single, but rather a set of key genes regulated by vitamins during infection, we could not identify common enrichments for response element motifs RARE (DR1, DR5) and VDRE in potential enhancer regions as reported before for proteins (see Data availability section for external supplement S10)59. However we demonstrate that vitamins counteract the transcriptional regulation of ncRNAs by A. fumigatus, C. albicans and the bacterium E. coli, thereby influencing the extent of the pro-inflammatory response. Cytokines, chemokines and apoptosis related antisense RNA genes are generally down-regulated by atRA. Retinoic acids are known to activate HOXD-AS1 (inducing the PI3K/Akt pathway) and the apoptosis repressor HOXA-AS2 in peripheral blood neutrophils20,21. We were able to verify these findings for human monocytes upon fungal infection shedding new light on lncRNAs. Furthermore, potential key marker genes in response to atRA include AC011899.9 and RP11-483P21.2 and RP11-588G21.2.
In conclusion, we have found several ncRNAs regulated by both vitamins during infection, which might play an important role regulating protein-coding genes. Moreover, we have identified potential marker genes for A. fumigatus (LINC00595, SBF2-AS1) and C. albicans (RP11-588G21.2, RP11-394l13.1) infections. The significance and fold change of expression of these genes is high enough, to assume their detectability in the early phase of fungal infection.
Additional Information
Accession codes: The raw Ion Proton sequence data in fastq format are stored in the Sequence Read Archive (SRA) at National Center for Biotechnology Information (NCBI) and can be accessed at NCBI homepage (accession number: SRP076532).
How to cite this article: Riege, K. et al. Massive Effect on LncRNAs in Human Monocytes During Fungal and Bacterial Infections and in Response to Vitamins A and D. Sci. Rep. 7, 40598; doi: 10.1038/srep40598 (2017).
Publisher's note: Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
References
ENCODE Project Consortium. Identification and analysis of functional elements in 1% of the human genome by the ENCODE pilot project. Nature 447, 799â816 (2007).
ENCODE Project Consortium. An integrated encyclopedia of DNA elements in the human genome. Nature 489, 57â74 (2012).
Nawrocki, E. P. et al. Rfam 12.0: updates to the RNA families database. Nucleic Acids Res 43, D130âD137 (2015).
Harrow, J. et al. GENCODE: the reference human genome annotation for The ENCODE Project. Genome Res 22, 1760â1774 (2012).
Doldo, E. et al. Vitamin A, cancer treatment and prevention: the new role of cellular retinol binding proteins. BioMed research international 2015, 624627 (2015).
Ordóñez-Morán, P., Dafflon, C., Imajo, M., Nishida, E. & Huelsken, J. HOXA5 counteracts stem cell traits by inhibiting Wnt signaling in colorectal cancer. Cancer cell 28, 815â829 (2015).
Chung, M., Lee, J., Terasawa, T., Lau, J. & Trikalinos, T. A. Vitamin D with or without calcium supplementation for prevention of cancer and fractures: an updated meta-analysis for the U.S. Preventive Services Task Force. Annals of internal medicine 155, 827â838 (2011).
Manson, J. E. et al. The VITamin D and OmegA-3 trial(vital): rationale and design of a large randomized controlled trial of vitamin D and marine omega-3 fatty acid supplements for the primary prevention of cancer and cardiovascular disease. Contemporary clinical trials 33, 159â171 (2012).
Cao, J. The functional role of long non-coding RNAs and epigenetics. Biol Proced Online 16, 11 (2014).
Cao, X. Self-regulation and cross-regulation of pattern-recognition receptor signalling in health and disease. Nat Rev Immunol 16, 35â50 (2015).
Carpenter, S. et al. A long noncoding RNA mediates both activation and repression of immune response genes. Science 341, 789â792 (2013).
Cawley, S. et al. Unbiased mapping of transcription factor binding sites along human chromosomes 21 and 22 points to widespread regulation of noncoding RNAs. Cell 116, 499â509 (2004).
Mallardo, M., Poltronieri, P. & DâUrso, O. F. Non-protein coding RNA biomarkers and differential expression in cancers: a review. J Exp Clin Cancer Res 27, 19 (2008).
Rossi, A. et al. Non-coding RNAs change their expression profile after retinoid induced differentiation of the promyelocytic cell line NB4. BMC Res Notes 3, 24 (2010).
Kour, S. & Rath, P. C. All-trans retinoic acid induces expression of a novel intergenic long noncoding RNA in adult rat primary hippocampal neurons. J Mol Neurosci 58, 266â276 (2016).
Zeng, C. et al. Inhibition of long non-coding RNA NEAT1 impairs myeloid differentiation in acute promyelocytic leukemia cells. BMC Cancer 14, 693 (2014).
Zhang, X. et al. A myelopoiesis-associated regulatory intergenic noncoding RNA transcript within the human HOXA cluster. Blood 113, 2526â2534 (2009).
Fatica, A. Noncoding RNAs in acute myeloid leukemia: From key regulators to clinical players. Scientifica (Cairo) 2012, 925758 (2012).
Zhang, X., Weissman, S. M. & Newburger, P. E. Long intergenic non-coding RNA HOTAIRM1 regulates cell cycle progression during myeloid maturation in NB4 human promyelocytic leukemia cells. RNA Biol 11, 777â787 (2014).
Yarmishyn, A. A. et al. HOXD-AS1 is a novel lncRNA encoded in HOXD cluster and a marker of neuroblastoma progression revealed via integrative analysis of noncoding transcriptome. BMC Genomics 15 Suppl 9, S7 (2014).
Zhao, H., Zhang, X., Frazao, J. B., Condino-Neto, A. & Newburger, P. E. HOX antisense lincRNA HOXA-AS2 is an apoptosis repressor in all trans retinoic acid treated NB4 promyelocytic leukemia cells. J Cell Biochem 114, 2375â2383 (2013).
Ma, X. Y. et al. Malat1 as an evolutionarily conserved lncRNA, plays a positive role in regulating proliferation and maintaining undifferentiated status of early-stage hematopoietic cells. BMC Genomics 16, 676 (2015).
Kour, S. & Rath, P. C. Age-related expression of a repeat-rich intergenic long noncoding RNA in the rat brain. Mol Neurobiol (2016).
Herriges, M. J. et al. Long noncoding RNAs are spatially correlated with transcription factors and regulate lung development. Genes Dev 28, 1363â1379 (2014).
Jiang, Y. J. & Bikle, D. D. LncRNA: a new player in 1 alpha, 25 (OH)(2) vitamin D(3) /VDR protection against skin cancer formation. Exp Dermatol 23, 147â150 (2014).
Jiang, Y. J. & Bikle, D. D. LncRNA profiling reveals new mechanism for VDR protection against skin cancer formation. J Steroid Biochem Mol Biol 144 Pt A, 87â90 (2014).
Long, M. D., Sucheston-Campbell, L. E. & Campbell, M. J. Vitamin D receptor and RXR in the post-genomic era. J Cell Physiol 230, 758â766 (2015).
Zhu, Q. H., Stephen, S., Taylor, J., Helliwell, C. A. & Wang, M. B. Long noncoding RNAs responsive to Fusarium oxysporum infection in Arabidopsis thaliana . New Phytol 201, 574â584 (2014).
Muhammad, S. A. et al. MicroRNA expression profiling of human respiratory epithelium affected by invasive candida infection. PLoS One 10, e0136454 (2015).
Khoo, A.-L. et al. 1,25-dihydroxyvitamin D3 modulates cytokine production induced by candida albicans: Impact of seasonal variation of immune responses. Journal of Infectious Diseases 203, 122â130 (2011).
Klassert, T. E. et al. Modulatory role of vitamin A on the Candida albicans-induced immune response in human monocytes. Medical Microbiology and Immunology 203, 415â424 (2014).
Klassert, T. E. et al. Differential effect of vitamins A and D on the transcriptional landscape of human monocytes during infection. Sci. Rep. 7, 40599; doi: 10.1038/srep40599 (2017).
Andrews, S. et al. FastQC: A quality control tool for high throughput sequence data. Reference Source (2010).
Schmieder, R. & Edwards, R. Quality control and preprocessing of metagenomic datasets. Bioinformatics 27, 863â864 (2011).
Kopylova, E., Noé, L. & Touzet, H. SortMeRNA: fast and accurate filtering of ribosomal RNAs in metatranscriptomic data. Bioinformatics 28, 3211â3217 (2012).
Hoffmann, S. et al. Fast mapping of short sequences with mismatches, insertions and deletions using index structures. PLoS Comput Biol 5, e1000502 (2009).
Anders, S., Pyl, P. T. & Huber, W. HTSeq â a Python framework to work with high-throughput sequencing data. Bioinformatics btu638 (2014).
Love, M. I., Huber, W. & Anders, S. Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome biology 15, 1â21 (2014).
Li, B., Ruotti, V., Stewart, R. M., Thomson, J. A. & Dewey, C. N. RNA-Seq gene expression estimation with read mapping uncertainty. Bioinformatics 26, 493â500 (2010).
Huntley, M. A. et al. ReportingTools: an automated result processing and presentation toolkit for high-throughput genomic analyses. Bioinformatics 29, 3220â3221 (2013).
Lorenz, R. et al. ViennaRNA package 2.0. Algorithms Mol Biol 6, 26 (2011).
Thomason, M. K. & Storz, G. Bacterial antisense RNAs: how many are there, and what are they doing? Annu Rev Genet 44, 167â188 (2010).
Chen, J. et al. Over 20% of human transcripts might form sense-antisense pairs. Nucleic Acids Res 32, 4812â4820 (2004).
Kiyosawa, H. et al. Antisense transcripts with FANTOM2 clone set and their implications for gene regulation. Genome Res 13, 1324â1334 (2003).
Schwartz, J. C. et al. Antisense transcripts are targets for activating small RNAs. Nature structural & molecular biology 15, 842â848 (2008).
Kanehisa, M. & Goto, S. KEGG: kyoto encyclopedia of genes and genomes. Nucleic Acids Res 28, 27â30 (2000).
Kanehisa, M., Sato, Y., Kawashima, M., Furumichi, M. & Tanabe, M. KEGG as a reference resource for gene and protein annotation. Nucleic Acids Res 44, D457âD462 (2015).
Lv, J. et al. High expression of long non-coding RNA SBF2-AS1 promotes proliferation in non-small cell lung cancer. J Exp Clin Cancer Res 35, 75 (2016).
Melton, D. A. & Hrvatin, S. inventors; Markers for mature beta-cells and methods of using the same. United States patent US 20140329704 (2014 Nov 06).
Nolte, E. N. et al. Deep sequencing of rna from immune cell-derived vesicles uncovers the selective incorporation of small non-coding RNA biotypes with potential regulatory functions. Nucleic acids research 40, 9272â9285 (2012).
Zhang, H. et al. Long non-coding RNA expression profiles of hepatitis C virus-related dysplasia and hepatocellular carcinoma. Oncotarget 6, 43770â43778 (2015).
Parrado, A. et al. Retinoic acid receptor alpha1 variants, RARalpha1DeltaB and RARalpha1DeltaBC, define a new class of nuclear receptor isoforms. Nucleic acids research 29, 4901â4908 (2001).
Miano, J. M. & Berk, B. C. Retinoids: versatile biological response modifiers of vascular smooth muscle phenotype. Circ Res 87, 355â362 (2000).
Clark, M. B. et al. Genome-wide analysis of long noncoding RNA stability. Genome research 22, 885â898 (2012).
Gibb, E. A., Brown, C. J. & Lam, W. L. The functional role of long non-coding RNA in human carcinomas. Molecular cancer 10, 38 (2011).
Derrien, T. et al. The GENCODE v7 catalog of human long noncoding RNAs: analysis of their gene structure, evolution, and expression. Genome research 22, 1775â1789 (2012).
Lee, J. T. Epigenetic regulation by long noncoding RNAs. Science (New York, N.Y.) 338, 1435â1439 (2012).
Zhao, Fengjin & Lin, Tianxin & He, Wang & Han, Jinli & Zhu, Dingjun & Hu, Kaishun & Li, Weicong & Zheng, Zaosong & Huang, Jian & Xie, Wenlian . Knockdown of a novel lincRNA AATBC suppresses proliferation and induces apoptosis in bladder cancer. Oncotarget 6, 1064 (2015).
Lalevee, S. et al. Genome-wide in silico identification of new conserved and functional retinoic acid receptor response elements(direct repeats separated by 5âbp). Journal of Biological Chemistry 286, 33322â33334 (2011).
Hulsen, T., de Vlieg, J. & Alkema, W. BioVennâa web application for the comparison and visualization of biological lists using area-proportional Venn diagrams. BMC genomics 9, 488 (2008).
Acknowledgements
This work was supported by the German Research Foundation (DFG). Collaborative Research Center/Transregio 124 - âPathogenic fungi and their human host: Networks of Interactionâ, subproject A5 (T.K., H.S.), subproject B3 (K.R., M.M.), Carl-Zeiss Stiftung (F.H.), DFG project MA-5082/1 (K.R., M.M.) and the German Federal Ministry of Education and Research (BMBF 03Z2JN22 to H.S.). The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript. We want to thank Simone Tänzer and Moira Walters for excellent technical assistance.
Author information
Authors and Affiliations
Contributions
K.R., M.H., T.E.K., H.S. and M.M. conceived and designed the experiments. T.E.K. and J.B. conducted the experiments. M.H., K.R., E.B., M.C., N.M., M.S. and B.V. analysed the data. T.E.K. and H.S. contributed reagents/materials/analysis tools. K.R., M.H., F.H. and M.M. wrote the manuscript. All authors reviewed the manuscript.
Corresponding author
Ethics declarations
Competing interests
The authors declare no competing financial interests.
Rights and permissions
This work is licensed under a Creative Commons Attribution 4.0 International License. The images or other third party material in this article are included in the articleâs Creative Commons license, unless indicated otherwise in the credit line; if the material is not included under the Creative Commons license, users will need to obtain permission from the license holder to reproduce the material. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/
About this article
Cite this article
Riege, K., Hölzer, M., Klassert, T. et al. Massive Effect on LncRNAs in Human Monocytes During Fungal and Bacterial Infections and in Response to Vitamins A and D. Sci Rep 7, 40598 (2017). https://doi.org/10.1038/srep40598
Received:
Accepted:
Published:
DOI: https://doi.org/10.1038/srep40598
This article is cited by
-
Tracing vitamins on the long non-coding lane of the transcriptome: vitamin regulation of LncRNAs
Genes & Nutrition (2024)
-
Exploring the hidden hot world of long non-coding RNAs in thermophilic fungus using a robust computational pipeline
Scientific Reports (2024)
-
Molecular pathways identified from single nucleotide polymorphisms demonstrate mechanistic differences in systemic lupus erythematosus patients of Asian and European ancestry
Scientific Reports (2023)
-
Antifungal Effect of Long Noncoding RNA 9708-1 in the Vulvovaginal Candidiasis Murine Model
Mycopathologia (2021)
-
Genome-wide analysis of differentially expressed profiles of mRNAs, lncRNAs and circRNAs in chickens during Eimeria necatrix infection
Parasites & Vectors (2020)