Next Article in Journal
Direct Comparison of Chol-siRNA Polyplexes and Chol-DsiRNA Polyplexes Targeting STAT3 in a Syngeneic Murine Model of TNBC
Previous Article in Journal
CLASH Analyst: A Web Server to Identify In Vivo RNA–RNA Interactions from CLASH Data
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
ncRNA
Volume 8
Issue 1
10.3390/ncrna8010007
ncrna-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

LncRNAs Are Differentially Expressed between Wildtype and Cell Line Strains of African Trypanosomes

by
Hyung Chul Kim
1 and
Emmitt R. Jolly
1,2,3,*
1
Department of Biology, Case Western Reserve University, Cleveland, OH 44106, USA
2
Center for Global Health and Disease, Case Western Reserve University, Cleveland, OH 44106, USA
3
Center for RNA Science and Therapeutics, Case Western Reserve University, Cleveland, OH 44106, USA
*
Author to whom correspondence should be addressed.
Non-Coding RNA 2022, 8(1), 7; https://doi.org/10.3390/ncrna8010007
Submission received: 23 November 2021 / Revised: 2 January 2022 / Accepted: 8 January 2022 / Published: 12 January 2022
Browse Figures
Review Reports Versions Notes

Abstract

:
Trypanosoma brucei is a parasitic protist that causes African sleeping sickness. The establishment of T. brucei cell lines has provided a significant advantage for the majority of T. brucei research. However, these cell lines were isolated and maintained in culture for decades, occasionally accumulating changes in gene expression. Since trypanosome strains have been maintained in culture for decades, it is possible that difference may have accumulated in fast-evolving non-coding RNAs between trypanosomes from the wild and those maintained extensively in cultures. To address this, we compared the lncRNA expression profile of trypanosomes maintained as cultured cell lines (CL) to those extracted from human patients, wildtype (WT). We identified lncRNAs from CL and WT from available transcriptomic data and demonstrate that CL and WT have unique sets of lncRNAs expressed. We further demonstrate that the unique and shared lncRNAs are differentially expressed between CL and WT parasites, and that these lncRNAs are more evenly up-regulated and down-regulated than protein-coding genes. We validated the expression of these lncRNAs using qPCR. Taken together, this study demonstrates that lncRNAs are differentially expressed between cell lines and wildtype T. brucei and provides evidence for potential evolution of lncRNAs, specifically in T. brucei maintained in culture.

1. Introduction

Trypanosoma brucei is the causative agent for African trypanosomiasis, a neglected tropical disease that leads to sleep disturbance and is almost always fatal if untreated [1]. The disease is transmitted by the bite of the tsetse fly, during which the protozoan T. brucei invades the mammalian host and enters the bloodstream. Mammalian host specificity depends on the subspecies of the parasite: T. brucei rhodesiense (Tbr) and T. brucei gambiense (Tbg) infect human hosts, while T. brucei brucei (Tbb) primarily infects cattle. Tbb is unable to survive in humans due to its susceptibility to trypanolytic factors (TLF) present in human sera [2]. Tbr overcomes this susceptibility by expressing serum resistance-associated (SRA) gene [3,4] while Tbg uses multiple independent components to nullify TLF [5,6,7,8,9]. Despite these differences, whole genome sequencing demonstrates that the human-infecting Tbr and cattle-infecting Tbb are highly similar [10].
In addition to causing debilitating disease, T. brucei has a distinctly unique evolutionary lineage and is an interesting research organism [11]. The establishment of a cell line for in vitro study has been crucial for trypanosome research [12,13]. However, the most commonly used cell lines for T. brucei studies are sourced from isolates collected decades ago [14]. In fact, the first trypanosome cell line to have its genome sequenced, Tbb TREU927, was isolated over 50 years before it was sequenced [14,15]. Furthermore, while cell lines allow for an in-depth study at the in vitro level, they may behave very differently from in vivo conditions. Cell lines in humans have been shown to exhibit a significant difference from in vivo samples [16,17], and in similar fashion, different T. brucei cell lines vary in their expression of RNA editing pathways [18]. Additionally, transcriptomic analysis between Tbr from human patients and cell-cultured Tbb shows that cultured trypanosomes have an expression profile that suggests faster proliferation when compared to human blood samples [19]. While the previous transcriptional analysis established differences in protein-coding genes, it did not also account for non-coding RNA.
Long non-coding RNA (lncRNA) is a class of non-coding RNA that is over 200 nucleotides in length with little to no protein-coding potential [20]. LncRNAs have diverse functions in gene regulation and interact with DNA, RNA, proteins, and chromatin [21,22,23,24,25,26,27]. In T. brucei, one lncRNA has been shown to have a functional role in parasite differentiation during the bloodstream stage [28]. LncRNAs have been reported to have faster evolutionary rates than protein-coding genes [29]. While lncRNA sequences and expression patterns evolve more rapidly than protein-coding genes, their tissue-specificity is usually conserved [30]. Although lncRNA sequences are not as conserved as protein-coding regions, they are more conserved than introns and random intergenic regions [31,32]. Thus, lncRNAs offer a unique opportunity to examine how decades of culturing has impacted T. brucei. In our study, we pooled T. brucei RNA-seq datasets from both cell line (CL) and wildtype (WT) sources and used two lncRNA identification tools to identify a high-confidence set of lncRNAs. We then show that the CL and WT strains each has a unique set of lncRNAs and also share common lncRNAs that have different expression patterns compared to protein-coding genes. We further verify the expression of newly identified lncRNA expression using quantitative RT-PCR.

2. Results

2.1. Novel lncRNAs Were Identified in Both CL and WT T. brucei

A vast majority of the available T. brucei cell lines were isolated decades ago [14], and have been extensively used for T. brucei research, including transcriptomic analysis. While RNA-seq has been performed several times in CL T. brucei [19,33,34,35,36], WT T. brucei has rarely been investigated using transcriptomic analysis. Despite recent interest in the roles of lncRNAs, a thorough comparison of the differences in these RNAs has not been performed. To compare lncRNAs from CL and WT T. brucei, we used thirty-three RNA-seq datasets from three published transcriptomic studies for lncRNA identification [19,33,34]. An RNA-seq dataset was considered CL if the study used an established T. brucei cell line in culture or in mice. On the other hand, an RNA-seq dataset was considered WT if the source parasites were recently collected from human blood or cerebrospinal fluid (CSF). For the CL studies, any datasets from the procyclic form of T. brucei was discarded, as no WT dataset is available for this life stage. For the WT datasets, any dataset that had lower than 1% alignment rate to the T. brucei TREU927 genome was discarded. This criterion was chosen as the WT samples collected directly from patients were contaminated with host tissue, and the published work showed that 1–10% of the reads were trypanosome-specific [19]. Overall, nine runs of CL and 24 WT datasets were used for this process (Supplementary Data S1). The principal-component analysis showed that the CL datasets and WT datasets form two major clusters with subclusters in each group due to heterogeneous nature of the data (Supplementary Figure S1).
Next, we designed a computational pipeline to identify potential lncRNAs from these datasets at high-confidence level by using two different lncRNA identification tools. After quality checking and trimming the adapters, the RNA-seq datasets were aligned to the T. brucei TREU927 genome and de novo assembled into transcripts. These assemblies were merged, and rRNAs and transcripts matching spliced-leader RNA (SL-RNA) were removed from the assembly. The resulting assembly then underwent the lncRNA identification process using two different tools (Figure 1). One of the tools, FEELnc, chosen is a high-performing program optimized for non-model organisms [37], and the second, CPC2, is a species-neutral lncRNA identification tool that has been used in non-model organisms [38,39,40,41]. For FEELnc, to increase the sensitivity and specificity, we used the shuffle mode of the coding potential calculation with a relaxed ORF definition, multi k-mer frequencies, and increased the mRNA and lncRNA specificity thresholds to 0.99 [37]. Additionally, the mono-exonic lncRNAs were included in this process, as T. brucei largely does not undergo cis-splicing, and previous work in various organisms have identified mono-exonic lncRNAs [42,43,44,45]. The results from FEELnc and CPC2 were then compared to select only the overlapping lncRNAs that were identified by both tools, resulting in a reduction of potentially identified lncRNAs. These high-confidence potential lncRNAs were then quantified using Salmon in mapping-based mode for further differential expression analysis [46].
LncRNAs identified through this process included novel transcripts and transcripts that have already been annotated in the T. brucei TREU927 genome as mRNAs, pseudogenic transcripts, and ncRNAs. Despite being classified as mRNAs, most of these transcripts are annotated as hypothetical proteins or putative genes that have not been experimentally verified. Since the lncRNA identification tools were used with a higher sensitivity and specificity threshold, we decided to retain these transcripts except for genes annotated with more defined protein domains in their gene names. However, we retained any transcripts that have been annotated as putative variant surface glycoprotein (VSG) or expression site-associated gene (ESAG) as most of these transcripts are annotated as pseudogenes, degenerates, or frameshifted. With these criteria, we removed 43 transcripts previously classified as mRNAs from the lncRNA list. In the past, pseudogenes have been regarded as “junk genes” that were considered defunct pieces of protein-coding genes. Nevertheless, some pseudogenes have now been shown to produce lncRNAs that can have further regulatory functions [47,48]. With the possibility of pseudogene-derived lncRNAs, we retained the transcripts that have been previously classified as pseudogenic transcripts. The transcripts that have been previously classified as ncRNA were also retained in this analysis, as their lengths were sufficient to classify them as lncRNAs. Using this approach, we identified 978 potential lncRNAs composing of 535 novel transcripts, 152 transcripts annotated as mRNAs, 246 annotated pseudogenic transcripts, and 45 annotated ncRNAs. With these potential lncRNAs identified, we then investigated whether they were differentially expressed between CL and WT T. brucei.

2.1.1. CL and WT T. brucei lncRNA Expression Patterns Are Distinct

In addition to CL T. brucei being isolated decades ago, most CL T. brucei were collected from ungulates [14]. We speculated that some of the newly identified potential lncRNAs might be expressed exclusively in CL or in WT. To test this, the transcripts were quantified and normalized for gene expression comparison, and only the transcripts with stringent, average normalized count over 10 were included for this analysis. For comparison, we also performed the same analysis for protein-coding genes. Among the 978 potential lncRNAs, 112 (11.4%) transcripts were expressed only in the CL, 82 (8.4%) transcripts were expressed only in the WT, and 457 (46.7%) transcripts were expressed in both CL and WT (Figure 2A). The remaining transcripts were expressed, but at low level in both strains. Compared to the lncRNAs, protein-coding genes showed more conservation. Among 9729 protein-coding transcripts, 185 (1.9%) transcripts were expressed only in the CL, 394 (4.0%) transcripts were expressed only in the WT, and 8426 (86.6%) transcripts were expressed in both CL and WT (Figure 2B). However, if we used a normalized count over 1 for inclusion, we identified 61 (6.2%) lncRNAs expressed only in the CL, 111 (11.3%) lncRNAs expressed only in the WT, and 646 (66.1%) lncRNAs expressed in both CL and WT (Supplementary Figure S2A). For protein-coding transcripts with normalized count over 1, 65 (0.7%) transcripts were expressed only in the CL, 356 (3.7%) transcripts were expressed only in the WT, and 8886 (91.3%) transcripts were expressed in both CL and WT (Supplementary Figure S2B). This analysis suggests that there has been greater divergence of CL and WT strain lncRNAs as compared to the protein-coding genes. The top 15 expressed CL, WT, and shared transcripts are highlighted in Supplementary Data S2.
While these unique sets of lncRNAs can have potentially important roles in either strain, the differential expression of the common sets of lncRNAs can also have a significant role in gene regulation. We used differential expression analysis to compare the log2 fold change (log2FC) of the WT to the CL. Among the 978 identified lncRNAs, 139 (14.2%) transcripts were up-regulated in WT while 133 (13.6%) transcripts were down-regulated in WT when compared to the CL (p < 0.01) (Figure 3A). We then compared this result to the differential expression of protein-coding transcripts. Among the 9729 protein-coding transcripts in T. brucei TREU927 genome, 1155 (11.9%) were up-regulated and 2855 (20.3%) were down-regulated in the WT when compared to the CL (p < 0.01) (Figure 3B). The broadly even up-regulated and down-regulated pattern of lncRNA compared to the largely down-regulated pattern of protein-coding transcripts in WT suggests that the lncRNAs have a more evenly up-regulated and down-regulated expression pattern in both strains when compared to protein-coding genes.

2.1.2. RT-qPCR Assay Confirms the In Silico lncRNA Expression in T. brucei

To validate the lncRNA identification pipeline, we performed qPCR analysis on several identified lncRNAs. However, after the data were assembled using Stringtie, we were quite surprised. The data predicted that a significant number of lncRNA were potentially spliced. We questioned the validity of the output, and consequently performed a de novo assembly of Tb protein coding genes where only two spliced genes are already known. Using the same analysis, our data were consistent and supported the currently established two spliced protein coding genes. We selected five lncRNAs with unique sequences for RT-qPCR analysis. Of the five lncRNAs, three were mono-exonic lncRNAs and we identified two as potentially multi-exonic lncRNAs. Tb927.11.10190 (telomerase reverse transcriptase) was used as the reference gene [49]. Expression was detected for all five transcripts, with mono-exonic lncRNAs expression being higher than the reference gene and the multi-exonic lncRNAs expression at similar or lower level than the reference gene (Figure 4). In addition, the melt curves showed one peak for each qPCR sample (Supplementary Figure S3). All five expressed lncRNAs showed very little to no genomic DNA contamination, confirming that the predicted lncRNAs are expressed in vivo.

3. Discussion

In this study, we pooled 33 RNA-seq datasets and used two lncRNA identification tools to identify potential lncRNAs in CL and WT T. brucei. We identified 978 lncRNAs comprised of 536 novel transcripts and 445 transcripts that have been annotated as mRNA, pseudogenic transcript, or ncRNA in the genome. We then provided data that both CL and WT each have a unique set of lncRNAs that are expressed exclusively in each strain, and that lncRNAs are differentially expressed between CL and WT. Further, we compared lncRNA differential expression to protein-coding genes to show that lncRNAs are proportionally more evenly up-regulated and down-regulated while protein-coding genes are largely down-regulated in WT when compared to CL. Finally, we verified the expression of these lncRNAs using RT-qPCR. Our findings support that cell lines have a distinct class of differentially genes when compared to non-cultured samples [16,17]. The differential expression analysis also bolsters the idea that lncRNAs are an ideal standard to assess the effect of in vitro culturing compared to in vivo samples.
A strange finding in the identified lncRNAs is that our data suggest that many of these transcripts are spliced. Trypanosomatids are known to have very few verified cis-spliced transcripts [50,51,52]. Trypanosome transcripts primarily undergo trans-splicing of SL-RNA to the 5′ cap of the polycistronic precursor to mature [50,53]. The known minimum eukaryotic intron length is 30 bp [54]. Our analysis found that 403 out of the 978 potential lncRNAs are multi-exonic by this definition, of which 400 are novel transcripts and three are pseudogenic transcripts. We tested two of the potentially multi-exonic lncRNAs by RT-qPCR that targeted the exon junctions (Figure 4). These data seemed to verify, with melt curves that showed a single peak, that trypanosome may indeed have additional cis-spliced genes among their lncRNA populations (Supplementary Figure S3). Despite the dearth of verified cis-spliced transcripts, the T. brucei TREU927 reference genome actually annotates many multi-exonic transcripts. Out of 9881 transcripts labeled as protein-coding genes, four transcripts have multiple exons. Furthermore, of the 1383 transcripts labeled as pseudogenes, 26 transcripts are annotated with multiple exons. The presence of a predicted higher number of multi-exonic transcripts in the pseudogenic transcripts, which can potentially lead to pseudogene-derived lncRNAs [47,48], suggests that multi-exonic lncRNAs may exist in T. brucei. Furthermore, the RT-qPCR for multi-exonic lncRNAs confirmed their expression despite the fact that the primers were designed to target exon junction that are thousands of base pairs apart. Therefore, the data suggest the existence of multi-exonic lncRNAs in T. brucei. Further investigation is in underway.
Previous transcriptomic analysis in cultured and bloodstream/CSF suggested that cultured trypanosomes potentially have faster division times than human-sourced trypanosomes due to higher levels of RNA-binding proteins, cytoskeletal proteins, mitochondrial proteins, translation initiation factors, and proteasome subunits [19]. Our analysis of protein-coding genes concurs that WT strains have a lower expression of the genes listed above, along with a large proportion of protein-coding genes. However, the differential expression pattern in protein-coding genes differs from that of lncRNAs. While protein-coding genes were largely down-regulated in WT, lncRNAs were proportionally up-regulated and down-regulated (Figure 3). Moreover, we found that both CL and WT strains each express a unique set of lncRNAs: 112 (11.4%) transcripts expressed only in the CL, 82 (8.4%) transcripts expressed only in the WT. Such discrepancy in the expression level and the unique transcripts support the idea that lncRNAs have faster evolutionary rate than protein-coding genes [29].
The existence of unique lncRNAs does not necessarily require that they have unique functions in CL or WT. LncRNAs are more conserved in terms of tissue-specificity [30], and both CL and WT are bloodstream forms of T. brucei. Furthermore, lncRNA conservation has multiple levels: sequence, structure, function, and syntenic context of neighboring genes [55]. Therefore, some of these unique lncRNAs in CL and WT may be conserved at the structural level, or they could be functional counterparts with different sequences and structures. For example, the human lncRNA JPX and the mouse lncRNA Jpx have different sequences and secondary structures, but both bind the same protein partner during X-chromosome inactivation [56]. Likewise, some of the potentially unique lncRNAs identified may have evolved new sequences or completely degenerated through numerous cell culture passages, leaving behind structural remains emblematic of functional cousins. Only one lncRNA has been functionally analyzed in T. brucei. The lncRNA, grumpy, controls differentiation from proliferating slender form to non-proliferating stumpy form in bloodstream T. brucei [28], signifying that lncRNAs are crucial in parasite biology. An important set of “mRNAs” that were retained as lncRNAs were annotated in the genome as degenerate and frame-shifted forms of variant surface glycoproteins (VSGs), which are essential for evasion of the host immune system [57]. In terms of lncRNA conservation, these transcripts could function in regulating VSG-switching to evade the host immune system, and have potential as therapeutic targets.

4. Materials and Methods

4.1. Genomic and Transcriptomic Data Analysis

The T. brucei TREU927 genome sequence and annotation were downloaded from TriTrypDB [58]. The most recent version, release 54, of the genome was used for this analysis.
Thirty-three RNA-seq datasets from three studies were used for this study: (1) four runs from cultured bloodstream form T. brucei RNA-seq datasets (PRJEB19907) [33], (2) five runs of CL T. brucei grown in mice and collected from blood and fat tissue (PRJEB11801) [34], and (3) twenty-four runs selected from WT T. brucei collected from human patients (PRJEB18524, PRJEB18523, and PRJEB23278) [19]. For the WT datasets, sixteen runs were directly collected from humans, either from blood or CSF, while eight runs were collected from human patients, injected into rats, and collected back again. The four cultured blood stream forms are T. brucei brucei EATRO 1125, clone AnTat1.1 strain while the five CL parasites were T. brucei brucei EATRO 1125 Antat1.1E 90–13 strain. The WT parasites were T. brucei rhodesiense collected from patients from Kaberamaido district of North-Eastern Uganda. FastQC was used on all datasets to check for the dataset quality [59] (https://www.bioinformatics.babraham.ac.uk/projects/fastqc/, accessed on 8 January 2020), and Trimmomatics was used to trim any adapters present in the datasets [60]. Next, the dataset was aligned to the T. brucei TREU927 genome using STAR [61] (https://github.com/alexdobin/STAR, accessed on 2 October 2020) and de novo assembled into transcripts using Stringtie [62]. Any RNA-seq runs with less than 1% alignment rate to the genome were discarded, resulting in the number of datasets mentioned above. The individual assemblies were then merged using Stringtie with the genome as the reference to generate a single GTF file for downstream analysis.

4.2. Identification of Long Noncoding RNAs

The merged assembly was processed through riboPicker to identify and remove any rRNAs [63] (http://ribopicker.sourceforge.net/, accessed on 5 June 2020). The following databases were used for riboPicker: SILVA rRNA database for large subunit (version 132) and small subunit (version 138) rRNA sequences [64] (https://www.arb-silva.de/, accessed on 5 June 2020) and Rfam database for 5S and 5.8S subunit rRNA sequences (release 14.2) [65] (https://rfam.org/, accessed on 5 June 2020). Since the SILVA small subunit data file was too large to index with riboPicker, it was broken into six smaller sized files for indexing.
After the rRNA filtering, the assembly was processed with two different lncRNA identification tools: FEELnc and CPC2. The FEELnc coding potential calculation module was performed on shuffle mode with relaxed ORF definition and multi k-mer frequencies to increase sensitivity of the module [37] (https://github.com/tderrien/FEELnc, accessed on 16 March 2020). To decrease the false positive rate due to the optimized sensitivity, the specificity thresholds for mRNA and lncRNA were each increased to 0.99. CPC2 was run on the default setting [38] (http://cpc2.gao-lab.org/, accessed on 10 September 2020). The transcript ID was used to find the lncRNA transcripts that were shared by the results from both lncRNA identification tools.

4.3. Statistical Analyses

The transcripts were quantified using the mapping-mode of Salmon [46] with the identified lncRNAs as reference transcripts. The protein-coding gene quantification was done with all annotated mRNA from the T. brucei TREU927 genome with the 152 transcripts identified as lncRNAs removed from the reference. Statistical analyses were performed with R environment (v. 4.0.5 [66]) with DESeq2 [67], readr [68], and tximport [69] libraries loaded. The log2 fold change shrinkage was performed with the apeglm [70] package. The p values used for all analyses were adjusted for multiple testing.

4.4. Animals and Parasites

Bloodstream form T. brucei brucei Lister 427 were cultured in T75 vented cap culture flasks (NEST, Wuxi, China) in complete HMI-9 media (Axenia Biologix, Dixon, CA, USA) containing 10% FBS (VWR, Radnor, PA, USA), 10% Serum Plus medium supplement (MilliporeSigma, Burlington, MA, USA), and 1X penicillin-streptomycin (Thermo Fisher, Waltham, MA, USA). The trypanosomes were maintained in a humidified incubator at 37 °C and 5% CO2, and they were sub-cultured every 48 h. After checking for parasite density for log growth phase, the trypanosomes were collected, pelleted, and flash frozen.

4.5. RNA Extraction and cDNA Generation

RNA was extracted from frozen pellets of T. brucei using the Direct-zol RNA Miniprep Kit (Zymo Research, Irvine, CA, USA) with a minor modification. To decrease the genomic DNA contamination, on-column DNase I treatment was performed twice. RNA concentration and quality were assessed on a Nanodrop 8000 Spectrophotometer (Thermo Scientific, Waltham, MA, USA). The cDNA generation was performed with Superscript II Reverse Transcriptase kit (Thermo Scientific, Waltham, MA, USA) with oligo dT primer following the manufacturer’s protocol.

4.6. RT-qPCR Assays

The qRT-PCR reactions were performed with Power SYBR Green RNA-to-CT 1-Step Kit (Applied Biosystems, Foster City, CA, USA). using two Q-qPCR machines (QuantaBio, Beverly, MA, USA). The no-RT reactions were set up to determine DNA contamination in the RNA samples. For reference gene, Tb927.11.10190 was chosen as it was determined to be an ideal reference gene for T. brucei qPCR studies [49]. The experiments were performed in triplicate and ΔCT method was used to analyze the transcript expression levels relative to the reference gene. The primer sequences are included in Supplementary Data S3.

Supplementary Materials

The following are available online at https://www.mdpi.com/article/10.3390/ncrna8010007/s1, Figure S1: PCA plot of RNA-seq files used for the analysis. Figure S2: Venn diagram of (A) identified T. brucei lncRNA and (B) protein-coding gene expression pattern in CL and WT strains. Only the transcripts with average >1 DESeq2 normalized count were included for this analysis. Figure S3: Melt curves of the RT-qPCR used to confirm lncRNA expression. TbLNC00418, TbLINC00030, and TbLINC00054 are mono-exonic transcripts. TbLNC00355, and TbLNC00041 are multi-exonic transcripts. For multi-exonic transcripts, the exon junctions were targeted for RT-qPCR analysis. Data S1: Table of RNA-seq datasets used for this experiment. All procyclic datasets and dataset with <1% alignment rate to T. brucei TREU927 genome were discarded from the study. Data S2: List of lncRNAs highly expressed in CL, highly expressed in WT, and highly differentially expressed; Data S3: Primer sequences used for qPCR.

Author Contributions

Conceptualization, H.C.K. and E.R.J.; methodology, H.C.K. and E.R.J.; software, H.C.K.; validation, H.C.K.; formal analysis, H.C.K. and E.R.J.; investigation, H.C.K. and E.R.J.; resources, E.R.J.; data curation, H.C.K.; writing—original draft preparation, H.C.K.; writing—review and editing, E.R.J.; visualization, H.C.K.; supervision, E.R.J.; project administration, E.R.J.; funding acquisition, E.R.J. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

All data analyzed in this study can be found in NCBI’s Sequence Read Archive under Accession numbers PRJEB19907, PRJEB11801, PRJEB18523, PRJEB18524, and PRJEB23278, and further included in Supplementary Data S1.

Acknowledgments

Research reported in this publication was supported by CWRU funds provided to ERJ. T. brucei brucei Lister 427 was graciously provided by Zac Mackey at Virginia Tech. (Blacksburg, VA, USA).

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Kennedy, P.G.E. Clinical features, diagnosis, and treatment of human African trypanosomiasis (sleeping sickness). Lancet Neurol. 2013, 12, 186–194. [Google Scholar] [CrossRef]
  2. Pays, E.; Vanhollebeke, B. Human innate immunity against African trypanosomes. Curr. Opin. Immunol. 2009, 21, 493–498. [Google Scholar] [CrossRef] [PubMed]
  3. De Greef, C.; Imberechts, H.; Matthyssens, G.; Van Meirvenne, N.; Hamers, R. A gene expressed only in serum-resistant variants of Trypanosoma brucei rhodesiense. Mol. Biochem. Parasitol. 1989, 36, 169–176. [Google Scholar] [CrossRef]
  4. Van Xong, H.; Vanhamme, L.; Chamekh, M.; Chimfwembe, C.E.; Van Den Abbeele, J.; Pays, A.; Van Meirvenne, N.; Hamers, R.; De Baetselier, P.; Pays, E. A VSG Expression Site–Associated Gene Confers Resistance to Human Serum in Trypanosoma rhodesiense. Cell 1998, 95, 839–846. [Google Scholar] [CrossRef] [Green Version]
  5. Capewell, P.; Clucas, C.; DeJesus, E.; Kieft, R.; Hajduk, S.; Veitch, N.; Steketee, P.C.; Cooper, A.; Weir, W.; MacLeod, A. The TgsGP gene is essential for resistance to human serum in Trypanosoma brucei gambiense. PLoS Pathog. 2013, 9, e1003686. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  6. Capewell, P.; Veitch, N.J.; Turner, C.M.R.; Raper, J.; Berriman, M.; Hajduk, S.L.; MacLeod, A. Differences between Trypanosoma brucei gambiense groups 1 and 2 in their resistance to killing by trypanolytic factor 1. PLoS Negl. Trop. Dis. 2011, 5, e1287. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  7. DeJesus, E.; Kieft, R.; Albright, B.; Stephens, N.A.; Hajduk, S.L. A single amino acid substitution in the group 1 Trypanosoma brucei gambiense haptoglobin-hemoglobin receptor abolishes TLF-1 binding. PLoS Pathog. 2013, 9, e1003317. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  8. Kieft, R.; Capewell, P.; Turner, C.M.R.; Veitch, N.J.; MacLeod, A.; Hajduk, S. Mechanism of Trypanosoma brucei gambiense (group 1) resistance to human trypanosome lytic factor. Proc. Natl. Acad. Sci. USA 2010, 107, 16137–16141. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  9. Uzureau, P.; Uzureau, S.; Lecordier, L.; Fontaine, F.; Tebabi, P.; Homblé, F.; Grélard, A.; Zhendre, V.; Nolan, D.P.; Lins, L.; et al. Mechanism of Trypanosoma brucei gambiense resistance to human serum. Nature 2013, 501, 430–434. [Google Scholar] [CrossRef]
  10. Sistrom, M.; Evans, B.; Benoit, J.; Balmer, O.; Aksoy, S.; Caccone, A. De Novo Genome Assembly Shows Genome Wide Similarity between Trypanosoma brucei brucei and Trypanosoma brucei rhodesiense. PLoS ONE 2016, 11, e0147660. [Google Scholar] [CrossRef] [Green Version]
  11. Baldauf, S.L.; Roger, A.J.; Wenk-Siefert, I.; Doolittle, W.F. A Kingdom-Level Phylogeny of Eukaryotes Based on Combined Protein Data. Science 2000, 290, 972. [Google Scholar] [CrossRef]
  12. Hirumi, H.; Doyle, J.J.; Hirumi, K. African trypanosomes: Cultivation of animal-infective Trypanosoma brucei in vitro. Science 1977, 196, 992. [Google Scholar] [CrossRef]
  13. Hirumi, H.; Hirumi, K. Continuous Cultivation of Trypanosoma brucei Blood Stream Forms in a Medium Containing a Low Concentration of Serum Protein without Feeder Cell Layers. J. Parasitol. 1989, 75, 985–989. [Google Scholar] [CrossRef] [Green Version]
  14. Gibson, W. The origins of the trypanosome genome strains Trypanosoma brucei brucei TREU 927, T. b. gambiense DAL 972, T. vivax Y486 and T. congolense IL3000. Parasites Vectors 2012, 5, 71. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  15. Berriman, M.; Ghedin, E.; Hertz-Fowler, C.; Blandin, G.; Renauld, H.; Bartholomeu Daniella, C.; Lennard Nicola, J.; Caler, E.; Hamlin Nancy, E.; Haas, B.; et al. The Genome of the African Trypanosome Trypanosoma brucei. Science 2005, 309, 416–422. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  16. Januszyk, M.; Rennert, R.C.; Sorkin, M.; Maan, Z.N.; Wong, L.K.; Whittam, A.J.; Whitmore, A.; Duscher, D.; Gurtner, G.C. Evaluating the Effect of Cell Culture on Gene Expression in Primary Tissue Samples Using Microfluidic-Based Single Cell Transcriptional Analysis. Microarrays 2015, 4, 540–550. [Google Scholar] [CrossRef]
  17. Zaitseva, M.; Vollenhoven, B.J.; Rogers, P.A.W. In vitro culture significantly alters gene expression profiles and reduces differences between myometrial and fibroid smooth muscle cells. Mol. Hum. Reprod. 2006, 12, 187–207. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  18. Kirby, L.E.; Koslowsky, D. Cell-line specific RNA editing patterns in Trypanosoma brucei suggest a unique mechanism to generate protein variation in a system intolerant to genetic mutations. Nucleic Acids Res. 2020, 48, 1479–1493. [Google Scholar] [CrossRef]
  19. Mulindwa, J.; Leiss, K.; Ibberson, D.; Kamanyi Marucha, K.; Helbig, C.; Melo do Nascimento, L.; Silvester, E.; Matthews, K.; Matovu, E.; Enyaru, J.; et al. Transcriptomes of Trypanosoma brucei rhodesiense from sleeping sickness patients, rodents and culture: Effects of strain, growth conditions and RNA preparation methods. PLoS Negl. Trop. Dis. 2018, 12, e0006280. [Google Scholar] [CrossRef] [PubMed]
  20. Kapranov, P.; Cheng, J.; Dike, S.; Nix, D.A.; Duttagupta, R.; Willingham, A.T.; Stadler, P.F.; Hertel, J.; Hackermüller, J.; Hofacker, I.L.; et al. RNA Maps Reveal New RNA Classes and a Possible Function for Pervasive Transcription. Science 2007, 316, 1484. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  21. Wang, F.; Chainani, P.; White, T.; Yang, J.; Liu, Y.; Soibam, B. Deep learning identifies genome-wide DNA binding sites of long non-coding RNAs. RNA Biol. 2018, 15, 1468–1476. [Google Scholar] [CrossRef]
  22. Jadaliha, M.; Gholamalamdari, O.; Tang, W.; Zhang, Y.; Petracovici, A.; Hao, Q.; Tariq, A.; Kim, T.G.; Holton, S.E.; Singh, D.K.; et al. A natural antisense lncRNA controls breast cancer progression by promoting tumor suppressor gene mRNA stability. PLoS Genet. 2018, 14, e1007802. [Google Scholar] [CrossRef] [PubMed]
  23. Rinn, J.L.; Kertesz, M.; Wang, J.K.; Squazzo, S.L.; Xu, X.; Brugmann, S.A.; Goodnough, L.H.; Helms, J.A.; Farnham, P.J.; Segal, E.; et al. Functional demarcation of active and silent chromatin domains in human HOX loci by non-coding RNAs. Cell 2007, 129, 1311–1323. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  24. Huarte, M.; Guttman, M.; Feldser, D.; Garber, M.; Koziol, M.J.; Kenzelmann-Broz, D.; Khalil, A.M.; Zuk, O.; Amit, I.; Rabani, M.; et al. A large intergenic non-coding RNA induced by p53 mediates global gene repression in the p53 response. Cell 2010, 142, 409–419. [Google Scholar] [CrossRef] [Green Version]
  25. Khalil, A.M.; Guttman, M.; Huarte, M.; Garber, M.; Raj, A.; Rivea Morales, D.; Thomas, K.; Presser, A.; Bernstein, B.E.; van Oudenaarden, A.; et al. Many human large intergenic non-coding RNAs associate with chromatin-modifying complexes and affect gene expression. Proc. Natl. Acad. Sci. USA 2009, 106, 11667–11672. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  26. Beletskii, A.; Hong, Y.K.; Pehrson, J.; Egholm, M.; Strauss, W.M. PNA interference mapping demonstrates functional domains in the non-coding RNA Xist. Proc. Natl. Acad. Sci. USA 2001, 98, 9215–9220. [Google Scholar] [CrossRef] [Green Version]
  27. Saghafi, T.; Taheri, R.A.; Parkkila, S.; Emameh, R.Z. Phytochemicals as Modulators of Long Non-Coding RNAs and Inhibitors of Cancer-Related Carbonic Anhydrases. Int. J. Mol. Sci. 2019, 20, 2939. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  28. Guegan, F.; Bento, F.; Neves, D.; Sequeira, M.; Notredame, C.; Figueiredo, L.M. A long non-coding RNA controls parasite differentiation in African trypanosomes. bioRxiv 2020. [Google Scholar] [CrossRef]
  29. Bu, D.; Luo, H.; Jiao, F.; Fang, S.; Tan, C.; Liu, Z.; Zhao, Y. Evolutionary annotation of conserved long non-coding RNAs in major mammalian species. Sci. China Life Sci. 2015, 58, 787–798. [Google Scholar] [CrossRef] [Green Version]
  30. Necsulea, A.; Soumillon, M.; Warnefors, M.; Liechti, A.; Daish, T.; Zeller, U.; Baker, J.C.; Grützner, F.; Kaessmann, H. The evolution of lncRNA repertoires and expression patterns in tetrapods. Nature 2014, 505, 635–640. [Google Scholar] [CrossRef]
  31. Ponjavic, J.; Ponting, C.P.; Lunter, G. Functionality or transcriptional noise? Evidence for selection within long non-coding RNAs. Genome Res. 2007, 17, 556–565. [Google Scholar] [CrossRef] [Green Version]
  32. Guttman, M.; Amit, I.; Garber, M.; French, C.; Lin, M.F.; Feldser, D.; Huarte, M.; Zuk, O.; Carey, B.W.; Cassady, J.P.; et al. Chromatin signature reveals over a thousand highly conserved large non-coding RNAs in mammals. Nature 2009, 458, 223–227. [Google Scholar] [CrossRef]
  33. Naguleswaran, A.; Doiron, N.; Roditi, I. RNA-Seq analysis validates the use of culture-derived Trypanosoma brucei and provides new markers for mammalian and insect life-cycle stages. BMC Genom. 2018, 19, 227. [Google Scholar] [CrossRef] [Green Version]
  34. Trindade, S.; Rijo-Ferreira, F.; Carvalho, T.; Pinto-Neves, D.; Guegan, F.; Aresta-Branco, F.; Bento, F.; Young, S.A.; Pinto, A.; Van Den Abbeele, J.; et al. Trypanosoma brucei Parasites Occupy and Functionally Adapt to the Adipose Tissue in Mice. Cell Host Microbe 2016, 19, 837–848. [Google Scholar] [CrossRef] [Green Version]
  35. Jensen, B.C.; Ramasamy, G.; Vasconcelos, E.J.R.; Ingolia, N.T.; Myler, P.J.; Parsons, M. Extensive stage-regulation of translation revealed by ribosome profiling of Trypanosoma brucei. BMC Genom. 2014, 15, 911. [Google Scholar] [CrossRef] [Green Version]
  36. Kolev, N.G.; Franklin, J.B.; Carmi, S.; Shi, H.; Michaeli, S.; Tschudi, C. The transcriptome of the human pathogen Trypanosoma brucei at single-nucleotide resolution. PLoS Pathog. 2010, 6, e1001090. [Google Scholar] [CrossRef] [Green Version]
  37. Wucher, V.; Legeai, F.; Hédan, B.; Rizk, G.; Lagoutte, L.; Leeb, T.; Jagannathan, V.; Cadieu, E.; David, A.; Lohi, H.; et al. FEELnc: A tool for long non-coding RNA annotation and its application to the dog transcriptome. Nucleic Acids Res. 2017, 45, e57. [Google Scholar] [CrossRef] [Green Version]
  38. Kang, Y.-J.; Yang, D.-C.; Kong, L.; Hou, M.; Meng, Y.-Q.; Wei, L.; Gao, G. CPC2: A fast and accurate coding potential calculator based on sequence intrinsic features. Nucleic Acids Res. 2017, 45, W12–W16. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  39. Vanamamalai, V.K.; Garg, P.; Kolluri, G.; Gandham, R.K.; Jali, I.; Sharma, S. Transcriptomic analysis to infer key molecular players involved during host response to NDV challenge in Gallus gallus (Leghorn & Fayoumi). Sci. Rep. 2021, 11, 8486. [Google Scholar] [PubMed]
  40. Xue, J.; Lv, Q.; Khas, E.; Bai, C.; Ma, B.; Li, W.; Cao, Q.; Fan, Z.; Ao, C. Tissue-specific regulatory mechanism of LncRNAs and methylation in sheep adipose and muscle induced by Allium mongolicum Regel extracts. Sci. Rep. 2021, 11, 9186. [Google Scholar] [CrossRef] [PubMed]
  41. Salabi, F.; Jafari, H.; Navidpour, S.; Sadr, A.S. Systematic and computational identification of Androctonus crassicauda long non-coding RNAs. Sci. Rep. 2021, 11, 4720. [Google Scholar] [CrossRef]
  42. Kim, H.C.; Khalil, A.M.; Jolly, E.R. LncRNAs in molluscan and mammalian stages of parasitic schistosomes are developmentally-regulated and coordinately expressed with protein-coding genes. RNA Biol. 2020, 17, 805–815. [Google Scholar] [CrossRef] [Green Version]
  43. Akay, A.; Jordan, D.; Navarro, I.C.; Wrzesinski, T.; Ponting, C.P.; Miska, E.A.; Haerty, W. Identification of functional long non-coding RNAs in C. elegans. BMC Biol. 2019, 17, 14. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  44. Wang, P.; Dai, L.; Ai, J.; Wang, Y.; Ren, F. Identification and functional prediction of cold-related long non-coding RNA (lncRNA) in grapevine. Sci. Rep. 2019, 9, 6638. [Google Scholar] [CrossRef]
  45. Zhang, T.; Zhang, X.; Han, K.; Zhang, G.; Wang, J.; Xie, K.; Xue, Q.; Fan, X. Analysis of long non-coding RNA and mRNA using RNA sequencing during the differentiation of intramuscular preadipocytes in chicken. PLoS ONE 2017, 12, e0172389. [Google Scholar]
  46. Patro, R.; Duggal, G.; Love, M.I.; Irizarry, R.A.; Kingsford, C. Salmon provides fast and bias-aware quantification of transcript expression. Nat. Methods 2017, 14, 417. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  47. Ma, H.; Ma, T.; Chen, M.; Zou, Z.; Zhang, Z. The pseudogene-derived long non-coding RNA SFTA1P suppresses cell proliferation, migration, and invasion in gastric cancer. Biosci. Rep. 2018, 38, BSR20171193. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  48. Milligan, M.J.; Lipovich, L. Pseudogene-derived lncRNAs: Emerging regulators of gene expression. Front. Genet. 2015, 5, 476. [Google Scholar] [CrossRef] [Green Version]
  49. Brenndörfer, M.; Boshart, M. Selection of reference genes for mRNA quantification in Trypanosoma brucei. Mol. Biochem. Parasitol. 2010, 172, 52–55. [Google Scholar] [CrossRef] [PubMed]
  50. Clayton, C. Regulation of gene expression in trypanosomatids: Living with polycistronic transcription. Open Biol. 2019, 9, 190072. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  51. Mair, G.; Shi, H.; Li, H.; Djikeng, A.; Aviles, H.O.; Bishop, J.R.; Falcone, F.H.; Gavrilescu, C.; Montgomery, J.L.; Santori, M.I.; et al. A new twist in trypanosome RNA metabolism: Cis-splicing of pre-mRNA. RNA 2000, 6, 163–169. [Google Scholar] [CrossRef] [Green Version]
  52. Chikne, V.; Gupta, S.K.; Doniger, T.; Cohen-Chalamish, S.; Waldman Ben-Asher, H.; Kolet, L.; Yahia, N.H.; Unger, R.; Ullu, E.; Kolev, N.G.; et al. The Canonical Poly (A) Polymerase PAP1 Polyadenylates Non-Coding RNAs and Is Essential for snoRNA Biogenesis in Trypanosoma brucei. J. Mol. Biol. 2017, 429, 3301–3318. [Google Scholar] [CrossRef]
  53. Siegel, T.N.; Gunasekera, K.; Cross, G.A.M.; Ochsenreiter, T. Gene expression in Trypanosoma brucei: Lessons from high-throughput RNA sequencing. Trends Parasitol. 2011, 27, 434–441. [Google Scholar] [CrossRef] [Green Version]
  54. Piovesan, A.; Caracausi, M.; Ricci, M.; Strippoli, P.; Vitale, L.; Pelleri, M.C. Identification of minimal eukaryotic introns through GeneBase, a user-friendly tool for parsing the NCBI Gene databank. DNA Res. 2015, 22, 495–503. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  55. Diederichs, S. The four dimensions of non-coding RNA conservation. Trends Genet. 2014, 30, 121–123. [Google Scholar] [CrossRef] [PubMed]
  56. Karner, H.; Webb, C.-H.; Carmona, S.; Liu, Y.; Lin, B.; Erhard, M.; Chan, D.; Baldi, P.; Spitale, R.C.; Sun, S. Functional Conservation of LncRNA JPX Despite Sequence and Structural Divergence. J. Mol. Biol. 2020, 432, 283–300. [Google Scholar] [CrossRef] [PubMed]
  57. Stijlemans, B.; Caljon, G.; Van Den Abbeele, J.; Van Ginderachter, J.A.; Magez, S.; De Trez, C. Immune Evasion Strategies of Trypanosoma brucei within the Mammalian Host: Progression to Pathogenicity. Front. Immunol. 2016, 7, 233. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  58. Aslett, M.; Aurrecoechea, C.; Berriman, M.; Brestelli, J.; Brunk, B.P.; Carrington, M.; Depledge, D.P.; Fischer, S.; Gajria, B.; Gao, X.; et al. TriTrypDB: A functional genomic resource for the Trypanosomatidae. Nucleic Acids Res. 2010, 38, D457–D462. [Google Scholar] [CrossRef] [Green Version]
  59. Andrews, S. FastQC: A Quality Control Tool for High Throughput Sequence Data. 2010. Available online: https://www.bioinformatics.babraham.ac.uk/projects/fastqc/ (accessed on 8 January 2020).
  60. Bolger, A.M.; Lohse, M.; Usadel, B. Trimmomatic: A flexible trimmer for Illumina sequence data. Bioinformatics 2014, 30, 2114–2120. [Google Scholar] [CrossRef] [Green Version]
  61. Dobin, A.; Davis, C.A.; Schlesinger, F.; Drenkow, J.; Zaleski, C.; Jha, S.; Batut, P.; Chaisson, M.; Gingeras, T.R. STAR: Ultrafast universal RNA-seq aligner. Bioinformatics 2013, 29, 15–21. [Google Scholar] [CrossRef]
  62. Pertea, M.; Pertea, G.M.; Antonescu, C.M.; Chang, T.-C.; Mendell, J.T.; Salzberg, S.L. StringTie enables improved reconstruction of a transcriptome from RNA-seq reads. Nat. Biotechnol. 2015, 33, 290. [Google Scholar] [CrossRef] [Green Version]
  63. Schmieder, R.; Lim, Y.W.; Edwards, R. Identification and removal of ribosomal RNA sequences from metatranscriptomes. Bioinformatics 2012, 28, 433–435. [Google Scholar] [CrossRef] [Green Version]
  64. Quast, C.; Pruesse, E.; Yilmaz, P.; Gerken, J.; Schweer, T.; Yarza, P.; Peplies, J.; Glöckner, F.O. The SILVA ribosomal RNA gene database project: Improved data processing and web-based tools. Nucleic Acids Res. 2013, 41, D590–D596. [Google Scholar] [CrossRef]
  65. Kalvari, I.; Argasinska, J.; Quinones-Olvera, N.; Nawrocki, E.P.; Rivas, E.; Eddy, S.R.; Bateman, A.; Finn, R.D.; Petrov, A.I. Rfam 13.0: Shifting to a genome-centric resource for non-coding RNA families. Nucleic Acids Res. 2018, 46, D335–D342. [Google Scholar] [CrossRef]
  66. R Core Team. R: A Language and Environment for Statistical Computing; R Foundation for Statistical Computing: Vienna, Austria, 2019. [Google Scholar]
  67. Love, M.I.; Huber, W.; Anders, S. Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biol. 2014, 15, 550. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  68. Wickham, H.; Hester, J. Readr: Read Rectangular Text Data; Version 2.0.2; R Foundation for Statistical Computing: Vienna, Austria, 2020. [Google Scholar]
  69. Soneson, C.; Love, M.; Robinson, M. Differential analyses for RNA-seq: Transcript-level estimates improve gene-level inferences. F1000Research 2016, 4, 1521. [Google Scholar] [CrossRef] [PubMed]
  70. Zhu, A.; Ibrahim, J.G.; Love, M.I. Heavy-tailed prior distributions for sequence count data: Removing the noise and preserving large differences. Bioinformatics 2018, 35, 2084–2092. [Google Scholar] [CrossRef] [PubMed]
Ncrna 08 00007 g001 550
Figure 1. Schematic of the computational pipeline designed to identify lncRNAs from pooled T. brucei RNA-seq datasets using two lncRNA identification tools.
Figure 1. Schematic of the computational pipeline designed to identify lncRNAs from pooled T. brucei RNA-seq datasets using two lncRNA identification tools.
Ncrna 08 00007 g001
Ncrna 08 00007 g002 550
Figure 2. Venn diagram of (A) identified T. brucei lncRNA and (B) protein-coding gene expression pattern in CL and WT strains. Only the transcripts with average >10 DESeq2 normalized count were included for this analysis.
Figure 2. Venn diagram of (A) identified T. brucei lncRNA and (B) protein-coding gene expression pattern in CL and WT strains. Only the transcripts with average >10 DESeq2 normalized count were included for this analysis.
Ncrna 08 00007 g002
Ncrna 08 00007 g003 550
Figure 3. LncRNAs are differentially expressed between CL and WT T. brucei. Red dots indicate transcripts that are significantly differentially expressed (p < 0.01), and the log fold change of WT transcripts are compared to that of CL strain. (A) MA-plot of the identified lncRNAs in WT compared to CL. (B) MA-plot of protein-coding transcripts in WT compared to CL.
Figure 3. LncRNAs are differentially expressed between CL and WT T. brucei. Red dots indicate transcripts that are significantly differentially expressed (p < 0.01), and the log fold change of WT transcripts are compared to that of CL strain. (A) MA-plot of the identified lncRNAs in WT compared to CL. (B) MA-plot of protein-coding transcripts in WT compared to CL.
Ncrna 08 00007 g003
Ncrna 08 00007 g004 550
Figure 4. RT-qPCR confirms the expression of identified lncRNAs. Five lncRNAs with unique sequences were selected for expression confirmation in CL T. brucei, and their relative expression level was calculated compared to the reference gene. TbLNC00418, TbLINC00030, and TbLINC00054 are mono-exonic transcripts. TbLNC00355, and TbLNC00041 are multi-exonic transcripts. Tb927.11.10190 is reference gene.
Figure 4. RT-qPCR confirms the expression of identified lncRNAs. Five lncRNAs with unique sequences were selected for expression confirmation in CL T. brucei, and their relative expression level was calculated compared to the reference gene. TbLNC00418, TbLINC00030, and TbLINC00054 are mono-exonic transcripts. TbLNC00355, and TbLNC00041 are multi-exonic transcripts. Tb927.11.10190 is reference gene.
Ncrna 08 00007 g004
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Kim, H.C.; Jolly, E.R. LncRNAs Are Differentially Expressed between Wildtype and Cell Line Strains of African Trypanosomes. Non-Coding RNA 2022, 8, 7. https://doi.org/10.3390/ncrna8010007

AMA Style

Kim HC, Jolly ER. LncRNAs Are Differentially Expressed between Wildtype and Cell Line Strains of African Trypanosomes. Non-Coding RNA. 2022; 8(1):7. https://doi.org/10.3390/ncrna8010007

Chicago/Turabian Style

Kim, Hyung Chul, and Emmitt R. Jolly. 2022. "LncRNAs Are Differentially Expressed between Wildtype and Cell Line Strains of African Trypanosomes" Non-Coding RNA 8, no. 1: 7. https://doi.org/10.3390/ncrna8010007

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop