Next Article in Journal
SMART Drumlines Ineffective in Catching White Sharks in the High Energy Capes Region of Western Australia: Acoustic Detections Confirm That Sharks Are Not Always Amenable to Capture
Previous Article in Journal
Cmcrf1, a Putative Zn2Cys6 Fungal Transcription Factor, Is Involved in Conidiation, Carotenoid Production, and Fruiting Body Development in Cordyceps militaris
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Revised Annotation and Characterization of Novel Aedes albopictus miRNAs and Their Potential Functions in Dengue Virus Infection

1
School of Biological Sciences, Universiti Sains Malaysia, Gelugor 11800, Pulau Pinang, Malaysia
2
Infectomics Cluster, Advanced Medical & Dental Institute, Universiti Sains Malaysia, Kepala Batas 13200, Pulau Pinang, Malaysia
3
Biotechnology Research Institute, Universiti Malaysia Sabah, Jalan UMS, Kota Kinabalu 88400, Sabah, Malaysia
4
Malaysia Genome and Vaccine Institute, National Institutes of Biotechnology Malaysia, Jalan Bangi, Kajang 43000, Selangor, Malaysia
*
Authors to whom correspondence should be addressed.
Biology 2022, 11(10), 1536; https://doi.org/10.3390/biology11101536
Submission received: 26 August 2022 / Revised: 2 October 2022 / Accepted: 9 October 2022 / Published: 19 October 2022

Abstract

:

Simple Summary

Aedes albopictus (Ae. albopictus) is an important vector of the dengue virus. Genetics and molecular studies of virus infection in mosquito vectors are important to uncover the basic biology of the virus. It has been reported that miRNAs are important and possess functional roles in virus infection in Ae. albopictus. Here, we report a comprehensive catalog of miRNAs using the latest genome version of Ae. albopictus. We discovered a total of 72 novel mature miRNAs, 44 of which were differentially expressed in C6/36 cells infected with the dengue virus. Target prediction analysis revealed that the differentially expressed miRNAs were involved in lipid metabolism and protein processing in the endoplasmic reticulum. Results from this study provide a valuable resource for researchers to study miRNAs in this mosquito vector, especially in host–virus interactions.

Abstract

The Asian tiger mosquito, Ae. albopictus, is a highly invasive species that transmits several arboviruses including dengue (DENV), Zika (ZIKV), and chikungunya (CHIKV). Although several studies have identified microRNAs (miRNAs) in Ae. albopictus, it is crucial to extend and improve current annotations with both the newly improved genome assembly and the increased number of small RNA-sequencing data. We combined our high-depth sequence data and 26 public datasets to re-annotate Ae. albopictus miRNAs and found a total of 72 novel mature miRNAs. We discovered that the expression of novel miRNAs was lower than known miRNAs. Furthermore, compared to known miRNAs, novel miRNAs are prone to expression in a stage-specific manner. Upon DENV infection, a total of 44 novel miRNAs were differentially expressed, and target prediction analysis revealed that miRNA-target genes were involved in lipid metabolism and protein processing in endoplasmic reticulum. Taken together, the miRNA annotation profile provided here is the most comprehensive to date. We believed that this would facilitate future research in understanding virus–host interactions, particularly in the role of miRNAs.

1. Introduction

Aedes albopictus (Ae. albopictus) is an important vector of arboviruses including dengue (DENV), Zika (ZIKV), and chikungunya (CHIKV) [1]. This mosquito is a robust and highly invasive species in both temperate and tropical regions of the world [2]. Due to its high adaptability, Ae. albopictus has spread rapidly and is considered a public threat throughout the world; hence, studying the fundamental biology of this mosquito is important for controlling its aggressive spread.
Previous studies in Ae. albopictus have shown the importance of non-coding RNAs (ncRNAs) in their development and virus infection [3,4,5,6,7]. Regulatory ncRNAs in metazoans include small RNAs, which consist of micro-RNAs (miRNAs), short-interfering RNAs (siRNAs), and PIWI-interacting RNAs (piRNAs). Particularly, miRNAs are 18–22 nucleotides (nts) in length and function in the regulation of gene expression by targeting messenger RNAs (mRNAs), resulting in the post-transcriptional silencing of the target genes [8,9]. Furthermore, miRNAs have been shown to play important roles in the development, growth, and infection of animals [10,11,12]. For instance, through a next-generation sequencing approach, multiple studies have reported that the expression profiles of Ae. albopictus miRNAs were altered at different developmental stages, during virus infection, and upon diapause conditions [3,6,13,14,15]. These studies suggest the involvement of miRNAs in the regulation of Ae. albopictus development and virus infection.
High-quality genomes coupled with high-depth RNA-sequencing data are important for accurate and comprehensive annotation of coding and non-coding genes. The alignment of small RNA reads to the reference genome is an essential step to annotate miRNAs in a species. Several studies used the genomes of closely related species as a proxy reference to facilitate the miRNA expression analysis [6,13,16]. For example, previous studies by Su et al. 2017 and Su et al. 2019 used Ae. aegypti miRNAs in miRBase (http://www.mirbase.org/, accessed on 25 Augest 2022) to evaluate the expression of Ae. albopictus miRNAs in the midgut when exposed to (DENV) infection [6,13]. Both studies reported novel Ae. albopictus miRNAs without mapping small RNA reads against the Ae. albopictus genome. Although this approach does not greatly affect the differential expression analysis of highly abundant miRNAs in non-model organisms [17], the level of error in annotation is unknown. By using the genomes of closely related species to discover miRNAs of another organism, we limit ourselves to only miRNAs that are conserved by both species. miRNA target prediction software, such as miRDeep2 [18], will excise potential precursors from the genome of closely related species based on the alignment of small RNA reads. Therefore, the predicted precursor sequences, which are derived from the reference genome of closely related species, may or may not be the true precursors of the newly identified miRNAs.
In this study, we aim to investigate miRNA expression upon DENV1 infection in Ae. Albopictus C6/36 cells. Before identifying the list of differentially expressed miRNAs upon infection, we first re-annotated miRNAs in Ae. albopictus using the updated genome assembly. The most recent and updated genome assembly of Ae. Albopictus was released in 2019 (AalbF2) using PacBio technology [19]. We decided to use this genome version for miRNA annotation. To gain a comprehensive annotation catalog, we combined our C6/36 cell-derived sequencing libraries with 26 publicly available datasets. These datasets encompassed multiple stages of Ae. albopictus development including embryo, larvae, pupae, adult males, sugar-fed/blood-fed females, diapause, and non-diapause pharate larvae. Finally, we predicted target sites of differentially expressed miRNAs to gain insights into their regulatory roles in Ae. albopictus during virus infection.

2. Materials and Methods

2.1. Cell Culture and Virus

Ae. albopictus C6/36 cells were obtained from American Type Culture Collection (ATCC: CRL-1660) and were cultured in Leibovitz’s L-15 medium (Gibco, Waltham, MA, USA, 41300039), supplemented with 10% Fetal Bovine Serum (FBS, Gibco, 10270) and 10% Tryptose Phosphate Broth (TPB) solution (Sigma, St. Louis, MO, USA, T9157). The C6/36 cells were incubated at 25 °C without CO2. BHK-21 cells (ATCC: CCL-10) were cultured at 37 °C in Dulbecco’s modified Eagle Medium (DMEM, Gibco, 11995065), supplemented with 10% FBS (Gibco, 10270) and 5% CO2. Dengue virus serotype 1 (Hawaiian strain) was propagated in C6/36 cells and titered using BHK-21 cells. Determination of DENV1 titer was carried out using 50% tissue culture infectious dose-cytopathic effect (TCID50-CPE) as previously described [5,20,21]. DENV1 (Hawaiian strain) used in this study was a gift from David Perera, University Malaysia Sarawak.

2.2. Virus Infection, RNA Extraction and Sequencing

C6/36 cells were infected with DENV1 at multiplicity of infection (MOI) of 0.25. Three days after infection, RNA extraction was carried out using miRNeasy Mini Kit 50 (Qiagen, Hilden, Germany, 217004) according to the manufacturer’s protocol. Total RNA was then subjected to next-generation sequencing. The RNA-sequencing libraries were prepared using standard Illumina protocols and sequenced using the HiSeq-SE50 platform, generating paired-end reads of 50 bp in size.

2.3. Verification of DENV1 Infection in C6/36 Cells

To verify the presence of DENV1 infection in C6/36 cells, the total RNA of both uninfected and DENV1-infected samples was subjected to cDNA synthesis using the Tetro cDNA synthesis kit (Bioline, London, UK, BIO-65042). DENV1-specific primer was used in cDNA synthesis. The polymerase chain reaction (PCR) assay was performed using DENV1-specific primers [22].

2.4. miRNA Identification

Publicly available small RNA-seq datasets were downloaded from NCBI Sequence Reads Archive (SRA). Prior to downstream analyses, small RNA-seq adapters were clipped using Trimmomatic version 0.38 [23] and low-quality reads were removed. Reads of 18–32 bp size were retained for downstream analysis.
To identify putatively novel miRNAs, clean reads from all small RNA-seq used in this study were pooled and aligned to the Ae. albopictus genome (AalbF2, NCBI) using miRDeep2 [18]. We used miRNAs predicted by Palatini et al. 2020 as known miRNAs [19]. For known miRNA precursors, we used all precursor miRNAs within the arthropod family that were retrieved from miRBase2.1. Predicted precursors that had significant Randfold p-value (p-value < 0.05) were retained. Minimum miRDeep score was set to >4. Then, we retained predicted miRNAs that shared the same seed sequence as arthropod miRNAs. For miRNAs that did not share the same seed sequence, only predicted precursors that had at least 1000 reads of mature miRNAs or those having at least 10 reads of predicted star sequence were retained.

2.5. Differential Expression of miRNA

The quantifier module in miRDeep2 package was used to quantify miRNAs in uninfected and DENV1-infected libraries [18]. In brief, total small RNA reads of both uninfected and DENV1-infected libraries were aligned against precursor miRNAs, and the resulting number of mapped reads indicates the miRNA abundance within each sample. Raw counts generated from the quantifier module were subjected to further analysis in the R/Bioconductor environment using edgeR package [24].

2.6. miRNA Target Prediction and Functional Analysis

Target prediction was performed using miRANDA [25] and RNAhybrid [26] using 3′UTR extracted from the Ae. Albopictus annotation file. miRANDA predictions were performed using the following parameters: score cut-off 140, energy cut-off −20, gap open penalty −9, gap extension penalty −4, seed regions at positions 2–8. RNAhybrid predictions were performed using the following parameters: energy cut-off −20, p-value < 0.1, binding required in miRNA positions 2–7. Functional analysis was conducted using g:Profiler [27]. We used g:SCS threshold for multiple testing corrections.

3. Results and Discussion

3.1. RNA-Sequencing Libraries and Read Mapping

In this study, we generated high-depth small RNA-seq libraries which were derived from DENV1-infected and uninfected Ae. albopictus C6/36 cell lines. All reads generated in this study were deposited in Short Read Archive (SRA) with the accession number SRP193815. A total of 270 million cleaned reads of 18–32 bp in length were acquired from our pooled set of six C6/36 small RNA-seq libraries (Table 1). These 18–32 bp raw reads were aligned to the Ae. albopictus genome (AalbF2, NCBI) and the resulting alignment files were subjected to downstream miRNA prediction analysis. In general, more than 90% of clean reads aligned to the reference genome (Table 1).
To visualize the heterogeneity of mapped reads in each library, the total number of small RNAs and their non-identical reads were plotted according to their size (Figure 1). The length distribution of total small RNAs revealed a major peak of 22 nts in size, but this peak collapsed when the total reads were converted into small RNA tags (Figure 1). The collapsibility of the peaks indicates the existence of miRNA in the libraries, since miRNAs are usually made up of relatively few identical sequences [28]. Unlike the 22 nt peak, peaks corresponding to longer small RNAs (24–32 nt) did not heavily collapse when converted into tags, suggesting the presence of a complex pool of small RNA populations [29].

3.2. Identification of Novel miRNAs in Ae. albopictus

To identify putatively novel miRNAs for Ae. albopictus, we combined reads from six libraries generated in this project with 26 public datasets (Supplemental Table S1), resulting in a total of approximately 350 million reads. The reads were then mapped to the genome (AalbF2, NCBI) and the resulting alignment file was fed into miRDeep2 for novel miRNA discovery. In our miRDeep2 analysis, we used miRNAs predicted from Palatini et al. 2020 as known Ae. albopictus miRNAs. We also defined other known miRNAs as miRNAs that have sequence conservations with miRNAs from other species, especially from the arthropod family curated from miRBase 2.1 (http://www.mirbase.org/, accessed on 25 August 2022). Most conserved miRNAs have been identified by previous studies [6,13,14,15,16,19,30].
A total of 72 novel mature miRNAs, which derived from 88 precursors, were identified. A complete list of both mature and precursor novel miRNAs identified in this study can be found in Supplemental Info S1. Palatini et al. 2020 reported a total of 121 mature miRNAs, which originated from 230 precursors from distinct loci in the genome of Ae. albopictus (AalBF2) [19]. We successfully confirmed the presence of miRNAs previously predicted by Palatini et al., 2020 in the AalbF2 genome version (Supplemental Info S1) [19]. Therefore, the total number of mature miRNAs in Ae. Albopictus is 193.
The majority of mature Ae. albopictus miRNAs were found to be 22 bp in length, while in the case of precursor miRNAs, most of them ranged between 58 and 65 bp (Figure 2). Ae. albopictus miRNAs varied in terms of genomic loci—49.9% and 48.2% of them were located in the intergenic and intronic regions, while those remaining resided within the coding sequence (CDS).
We found that 29 of the novel miRNAs shared the same seed sequence with arthropod miRNAs curated from miRBase 2.1 (http://www.mirbase.org/, accessed on 25 August 2022) (Supplementary Info S1). Despite sharing the same seed sequence, these miRNAs were still considered novel because the predicted precursor sequences were not similar to known precursors. miRNAs interact with their target mRNA via base pairing of its “seed” sequence, which is usually in the position of 2–8 from the 5′ end of the mature miRNAs [31]. Hence, these 29 novel miRNAs may target the same genes as the known miRNAs with which they shared the same seed region.
The number of novel miRNAs detected in this study was higher compared to other Ae. albopictus studies [15,30]. For example, Gu et al. 2013 reported the discovery of 15 novel miRNAs, whereas Liu et al., 2015 and Batz et al., 2017 identified 5 and 10 novel miRNAs in Ae. albopictus, respectively. This finding is presumably due to the high sequencing depth and high percentage of mapped reads in every replicate of C6/36 small RNA libraries generated in this study. It was reported that libraries with a lower number of reads make the discovery of miRNAs, especially the less expressed, more challenging and somewhat limited [32].
We also discovered that multiple locations in the Ae. albopictus genome can give rise to the same mature and precursor miRNAs. A total of 118 distinct locations were found to be the genomic loci of 72 novel miRNAs, suggesting a gene duplication event. This was also true for Ae. albopictus protein-coding genes. The duplication of several members of gene families involved in insecticide-resistance, diapause, immunity, olfaction, and sex determination was observed in the Ae. albopictus genome, contributing to its large genome size [33].
From the standpoint of miRNA biogenesis, the generation of functional novel mature miRNAs requires the transcription of certain genomic loci that are capable of producing RNA whose secondary structures are recognizable by Drosha and Pasha (DGCR8 in vertebrates) in the nucleus, and Dicer and Loquacious in the cytoplasm [9,34,35,36,37,38]. Therefore, the birth of novel miRNA requires a relatively simple prerequisite, which is loci in the genome that evolve to produce RNAs that fold in such a way that miRNA protein machineries can identify and process them, resulting in fully functional mature miRNAs [39,40]. Due to the inherent property of RNA to form imperfect folding (the biogenesis of miRNAs requires an imperfect folding of RNA hairpins), the emergence of novel miRNAs, from the evolutionary viewpoint, is prone to being more feasible than the birth of new coding genes [40].

3.3. Expression of Novel miRNAs in Different Stages of Ae. albopictus Development

We then assessed the dynamics of novel miRNA expression throughout the Ae. albopictus development to have a better understanding of its distribution. We used the public dataset from Gu et al., 2013, which consists of six different stages, namely 0–24 h embryo, larvae, pupae, adult males, blood-fed females and sugar-fed females [30]. In general, the expression level of novel miRNAs was lower than known miRNAs across all six stages (Figure 3a). Next, we determined the specificity of novel miRNA expression across different developmental stages. We calculated the specificity score of each miRNA using the Jensen–Shannon (JS) score, which ranges from zero to one [41]. This specificity metric quantifies the similarity in expression value across developmental stages. miRNAs with a JS score of one indicate the extreme case in which they are only expressed in one specific stage, whereas those having a score close to zero are ubiquitously expressed in all stages with a relatively similar value of expression [41].
Here, we discovered that novel miRNAs were more specific in expression across different developmental stages compared to known miRNAs (Figure 3b). For instance, 23% of known miRNAs had a JS score of 1, while for novel miRNAs, 46% of them had a JS score of 1. Known miRNAs, on the other hand, have a more ubiquitous expression profile. This observation indicates that, unlike known miRNAs, novel miRNAs possess a narrower window of expression than known miRNAs. We speculate that novel miRNAs in Ae. albopictus possibly function as stage-specific tuners and regulators of genes involved in the control of developmental transition. Despite their lower level of expression, the specificity of novel miRNAs expression hinted to the possibility that they may execute distinct functions at specific developmental time-points.

3.4. Ae. albopictus Novel miRNAs Were Differentially Expressed upon DENV Infection

Previous studies have reported that Ae. albopictus miRNAs experienced changes in expression level upon DENV infection [3,6,13,14]. We asked if our newly annotated miRNAs were differentially expressed upon virus infection. We chose DENV serotype 1 (DENV1) in this study. DENV serotype 2 (DENV2) has been used in many genome-wide expression studies either in mammalian or mosquito hosts [42,43,44,45,46,47,48,49,50]. DENV2 has been the center of attention in studies involving DENV, despite the fact that other serotypes are similarly as important as serotype 2. Different serotypes of DENV display different clinical manifestations and levels of severity. A case study in Singapore reported that patients with a DENV1 infection may exhibit more severe illnesses compared to those infected with DENV2 [51]. Furthermore, different symptoms were observed between DENV1- and DENV2-infected patients. This study suggests that different serotypes of DENV can or may exhibit different clinical manifestations and levels of severity [51]. Therefore, studies involving different serotypes of DENV are important to provide insights into molecular mechanisms or virus pathogenesis.
In this study, we discovered that our newly annotated Ae. albopictus miRNAs in C6/36 cells were differentially expressed after 3 days post-infection with DENV1. A total of 143 miRNAs were found to be differentially expressed (p-value < 0.05), whereby 44 of them were novel miRNAs. A full list of differentially expressed miRNAs and their corresponding p-values can be found in Supplemental Table S2. To visualize the relationship between log fold change (FC) and p-value, a volcano plot was plotted (Figure 4). From the plot, we observed that value of log2 FC ranged between −8 and 4.
Avila-Bonilla et al., 2017 performed a small RNA-sequencing of C6/36 cells that were persistently infected with DENV2 to analyze changes in miRNA expression profile (Avila-Bonilla et al., 2017). Several differentially expressed miRNAs in C6/36 cells persistently infected with DENV2, such as miR-190, miR-970, miR-263a, miR-927, and miR-92b, were also found to be differentially expressed in our DENV1-infected C6/36 cells [14]. Therefore, we hypothesized that these miRNAs may be involved in host–DENV interaction in Ae. albopictus, regardless of the virus serotype. Meanwhile, other miRNAs may participate in the serotype-specific antiviral immunity in Ae. albopictus.
It was discovered that the overexpression of miR-927 in C6/36 cells resulted in a higher viral titer and viral genome copy number of DENV2. This result suggests that miR-927 may possess pro-viral activity in Ae. Albopictus cells [52]. Furthermore, the same group also showed that filamin (FLN), which has been associated with the regulation of the Toll pathway, is the direct target of miR-927 by a dual-luciferase gene reporter assay [52,53,54,55]. miR-927 was also found to be differentially expressed in this study, despite the serotype of DENV used being different. This further supports the idea that certain Ae. albopictus miRNAs are involved in host–DENV interaction, regardless of the virus serotype.
Next, we performed a miRNA target prediction to gain further insight into the potential downstream effects of regulation by differentially expressed novel miRNAs in C6/36 cells following DENV1 infection. Undeniably, the identification of miRNA target genes remains a complicated task. In silico analysis and computational methods are important in predicting potential miRNA target genes, despite the outcomes of such approaches varying from one algorithm to another [17]. To alleviate the certainty in the miRNA target prediction, we merged the outcomes from two miRNA target prediction tools, namely miRanda [25] and RNAhybrid [26]. miRNA binding sites derived from both software packages were deemed to be highly confident.
Both miRANDA and RNAhybrid predicted a total of 10,867 genes as targets for differentially expressed miRNAs. Gene Ontology (GO) analyses using g:Profiler [27] showed that these target genes were involved in many types of biological processes, and molecular functions inside cells such as DNA and protein-binding, transport and localization (Figure 5). Meanwhile, KEGG pathway analysis revealed two important pathways: glycerophospholipid metabolism and protein processing in the endoplasmic reticulum (ER). A full list of GO and KEGG pathways from g:Profiler analysis can be found in Supplemental Table S4.
Previous studies have reported that among the cellular factors required for the life cycle of the enveloped virus, including DENV, are lipids [56,57,58]. DENV utilized a specific lipid, bis(monoacylglycerol)phosphate, which functions as a co-factor for viral fusion in late endosomes [58]. DENV entry into host cells is established via endocytosis and uses lipids as factors that trigger the fusion of viral particles into the endosome, hence delivering the viral genome into the cytosol [59,60,61]. One of the critical steps in the DENV life cycle is the translation of the viral genome into a single polypeptide chain by host ribosomes, and this process occurs in the rough endoplasmic reticulum [59,60]. Therefore, post-transcriptional regulation of proteins involved in lipid metabolism and protein processing in ER may be crucial in either facilitating or inhibiting virus replication.
The results presented in this study are exclusively based on immortalized C6/36 cells. One major limitation of cell line research is that it does not reflect the true nature of virus infection and pathogenesis in the whole organism. We expect that miRNA responses in Ae. albopictus adult mosquitoes are more complex than what we discovered here. However, data generated in this study are valuable for forthcoming explorations of miRNA functions in virus infection, especially in Aedes mosquitoes. Our work uncovers a part of the still-growing area of research in host–DENV interaction and offers insights into new research directions to pursue.

4. Conclusions

In summary, we provided a comprehensive genome-wide reannotation of Ae. albopictus miRNAs and their responses in DENV infection in C6/36 cells. Although our understanding of the functional roles of host miRNAs in DENV infection is somewhat limited, results generated in this study provide valuable resources for future research. Further research is required to fully understand the role of specific miRNAs in DENV infection and replication inside mosquito cells.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/biology11101536/s1, Supplemental Info S1: List of known and novel miRNAs in Ae. albopictus genome-wide; Supplemental Table S1: List of public datasets used in this study; Supplemental Table S2: List of differentially expressed miRNAs upon DENV1 infection in C6/36 cells; Supplemental Table S3: List of primers used in this study; Supplemental Table S4: List of GO terms of target genes of the differentially expressed miRNAs.

Author Contributions

Conceptualization: A.A. and G.A.; methodology: A.A; data curation; A.A. and G.A.; writing—original draft preparation: A.A.; writing—review and editing: M.A.H. and G.A.; supervision: G.A. and M.A.Y.; project administration: M.A.H. and G.A.; funding acquisition: G.A. and M.A.Y. All authors have read and agreed to the published version of the manuscript.

Funding

This study was funded by the Universiti Sains Malaysia Research Universiti Grant (1001/PBIOLOGI/811320) and MOSTI Science Fund (305/PBIOLOGI/613238).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

All reads generated in this study are available in Short Read Archive (SRA) with the accession number of SRP193815.

Acknowledgments

We thank David Perera for providing the DENV1 virus (Hawaiian strain).

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Gloria-Soria, A.; Payne, A.F.; Bialosuknia, S.M.; Stout, J.; Mathias, N.; Eastwood, G.; Ciota, A.T.; Kramer, L.D.; Armstrong, P.M. Vector Competence of Aedes albopictus Populations from the Northeastern United States for Chikungunya, Dengue, and Zika Viruses. Am. J. Trop. Med. Hyg. 2021, 104, 1123–1130. [Google Scholar] [CrossRef] [PubMed]
  2. Ryan, S.J.; Carlson, C.J.; Mordecai, E.A.; Johnson, L.R. Global expansion and redistribution of Aedes-borne virus transmission risk with climate change. PLoS Negl. Trop. Dis. 2019, 13, e0007213. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  3. Liu, Y.; Zhou, Y.; Wu, J.; Zheng, P.; Li, Y.; Zheng, X.; Puthiyakunnon, S.; Tu, Z.; Chen, X.-G. The expression profile of Aedes albopictus miRNAs is altered by dengue virus serotype-2 infection. Cell Biosci. 2015, 5, 16. [Google Scholar] [CrossRef] [Green Version]
  4. Liu, W.; Cheng, P.; Zhang, K.; Gong, M.; Zhang, Z.; Zhang, R. Systematic identification and characterization of long noncoding RNAs (lncRNAs) during Aedes albopictus development. PLoS Negl. Trop. Dis. 2022, 16, e0010245. [Google Scholar] [CrossRef]
  5. Azlan, A.; Obeidat, S.M.; Das, K.T.; Yunus, M.A.; Azzam, G. Genome-wide identification of Aedes albopictus long noncoding RNAs and their association with dengue and Zika virus infection. PLoS Negl. Trop. Dis. 2021, 15, e0008351. [Google Scholar] [CrossRef] [PubMed]
  6. Su, J.; Wang, G.; Li, C.; Xing, D.; Yan, T.; Zhu, X.; Liu, Q.; Wu, Q.; Guo, X.; Zhao, T. Screening for differentially expressed miRNAs in Aedes albopictus (Diptera: Culicidae) exposed to DENV-2 and their effect on replication of DENV-2 in C6/36 cells. Parasit. Vectors 2019, 12, 44. [Google Scholar] [CrossRef]
  7. Xu, Y.; Dong, Y.; Xu, Y.; Lai, Z.; Jin, B.; Hao, Y.; Gao, Y.; Sun, Y.; Chen, X.; Gu, J. Differentiation of Long Non-Coding RNA and mRNA Expression Profiles in Male and Female Aedes albopictus. Front. Genet. 2019, 10, 975. [Google Scholar] [CrossRef]
  8. O’Brien, J.; Hayder, H.; Zayed, Y.; Peng, C. Overview of MicroRNA Biogenesis, Mechanisms of Actions, and Circulation. Front. Endocrinol. 2018, 9, 402. [Google Scholar] [CrossRef] [Green Version]
  9. Azlan, A.; Dzaki, N.; Azzam, G. Argonaute: The executor of small RNA function. J. Genet. Genom. 2016, 43, 481–494. [Google Scholar] [CrossRef]
  10. Kunze-Schumacher, H.; Krueger, A. The Role of MicroRNAs in Development and Function of Regulatory T Cells—Lessons for a Better Understanding of MicroRNA Biology. Front. Immunol. 2020, 11, 2185. [Google Scholar] [CrossRef]
  11. Tribolet, L.; Kerr, E.; Cowled, C.; Bean, A.G.D.; Stewart, C.R.; Dearnley, M.; Farr, R.J. MicroRNA Biomarkers for Infectious Diseases: From Basic Research to Biosensing. Front. Microbiol. 2020, 11, 1197. [Google Scholar] [CrossRef] [PubMed]
  12. Wong, R.R.; Abd-Aziz, N.; Affendi, S.; Poh, C.L. Role of microRNAs in antiviral responses to dengue infection. J. Biomed. Sci. 2020, 27, 4. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  13. Su, J.; Li, C.; Zhang, Y.; Yan, T.; Zhu, X.; Zhao, M.; Xing, D.; Dong, Y.; Guo, X.; Zhao, T. Identification of microRNAs expressed in the midgut of Aedes albopictus during dengue infection. Parasit. Vectors 2017, 10, 63. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  14. Avila-Bonilla, R.G.; Yocupicio-Monroy, M.; Marchat, L.A.; De Nova-Ocampo, M.A.; Del Ángel, R.M.; Salas-Benito, J.S. Analysis of the miRNA profile in C6/36 cells persistently infected with dengue virus type 2. Virus Res. 2017, 232, 139–151. [Google Scholar] [CrossRef] [PubMed]
  15. Batz, Z.A.; Goff, A.C.; Armbruster, P.A. MicroRNAs are differentially abundant during Aedes albopictus diapause maintenance but not diapause induction. Insect Mol. Biol. 2017, 26, 721–733. [Google Scholar] [CrossRef]
  16. Etebari, K.; Hussain, M.; Asgari, S. Identification of microRNAs from Plutella xylostella larvae associated with parasitization by Diadegma semiclausum. Insect Biochem. Mol. Biol. 2013, 43, 309–318. [Google Scholar] [CrossRef]
  17. Etebari, K.; Asgari, S. Accuracy of MicroRNA Discovery Pipelines in Non-Model Organisms Using Closely Related Species Genomes. PLoS ONE 2014, 9, e84747. [Google Scholar] [CrossRef] [Green Version]
  18. Friedländer, M.R.; Mackowiak, S.D.; Li, N.; Chen, W.; Rajewsky, N. miRDeep2 accurately identifies known and hundreds of novel microRNA genes in seven animal clades. Nucleic Acids Res. 2012, 40, 37–52. [Google Scholar] [CrossRef] [Green Version]
  19. Palatini, U.; Masri, R.A.; Cosme, L.V.; Koren, S.; Thibaud-Nissen, F.; Biedler, J.K.; Krsticevic, F.; Johnston, J.S.; Halbach, R.; Crawford, J.E.; et al. Improved reference genome of the arboviral vector Aedes albopictus. Genome Biol. 2020, 21, 215. [Google Scholar] [CrossRef]
  20. Atieh, T.; Baronti, C.; de Lamballerie, X.; Nougairède, A. Simple reverse genetics systems for Asian and African Zika viruses. Sci. Rep. 2016, 6, 39384. [Google Scholar] [CrossRef]
  21. Li, J.; Hu, D.; Ding, X.; Chen, Y.; Pan, Y.; Qiu, L.; Che, X. Enzyme-Linked Immunosorbent Assay-Format Tissue Culture Infectious Dose-50 Test for Titrating Dengue Virus. PLoS ONE 2011, 6, e22553. [Google Scholar] [CrossRef]
  22. Johnson, B.W.; Russell, B.J.; Lanciotti, R.S. Serotype-Specific Detection of Dengue Viruses in a Fourplex Real-Time Reverse Transcriptase PCR Assay. J. Clin. Microbiol. 2005, 43, 4977–4983. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  23. Bolger, A.M.; Lohse, M.; Usadel, B. Trimmomatic: A flexible trimmer for Illumina sequence data. Bioinformatics 2014, 30, 2114–2120. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  24. Robinson, M.D.; McCarthy, D.J.; Smyth, G.K. edgeR: A Bioconductor package for differential expression analysis of digital gene expression data. Bioinformatics 2010, 26, 139–140. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  25. Enright, A.J.; John, B.; Gaul, U.; Tuschl, T.; Sander, C.; Marks, D.S. MicroRNA targets in Drosophila. Genome Biol. 2003, 5, R1. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  26. Krüger, J.; Rehmsmeier, M. RNAhybrid: microRNA target prediction easy, fast and flexible. Nucleic Acids Res. 2006, 34, W451–W454. [Google Scholar] [CrossRef] [PubMed]
  27. Raudvere, U.; Kolberg, L.; Kuzmin, I.; Arak, T.; Adler, P.; Peterson, H.; Vilo, J. g:Profiler: A web server for functional enrichment analysis and conversions of gene lists (2019 update). Nucleic Acids Res. 2019, 47, W191–W198. [Google Scholar] [CrossRef] [Green Version]
  28. Armisen, J.; Gilchrist, M.J.; Wilczynska, A.; Standart, N.; Miska, E.A. Abundant and dynamically expressed miRNAs, piRNAs, and other small RNAs in the vertebrate Xenopus tropicalis. Genome Res. 2009, 19, 1766–1775. [Google Scholar] [CrossRef] [Green Version]
  29. Harding, J.L.; Horswell, S.; Heliot, C.; Armisen, J.; Zimmerman, L.B.; Luscombe, N.M.; Miska, E.A.; Hill, C.S. Small RNA profiling of Xenopus embryos reveals novel miRNAs and a new class of small RNAs derived from intronic transposable elements. Genome Res. 2014, 24, 96–106. [Google Scholar] [CrossRef] [Green Version]
  30. Gu, J.; Hu, W.; Wu, J.; Zheng, P.; Chen, M.; James, A.A.; Chen, X.; Tu, Z. miRNA genes of an invasive vector mosquito, Aedes albopictus. PLoS ONE 2013, 8, e67638. [Google Scholar] [CrossRef]
  31. Lai, E.C. Micro RNAs are complementary to 3′ UTR sequence motifs that mediate negative post-transcriptional regulation. Nat. Genet. 2002, 30, 363–364. [Google Scholar] [CrossRef] [PubMed]
  32. Metpally, R.; Nasser, S.; Malenica, I.; Courtright, A.; Carlson, E.; Ghaffari, L.; Villa, S.; Tembe, W.; Van Keuren-Jensen, K. Comparison of Analysis Tools for miRNA High Throughput Sequencing Using Nerve Crush as a Model. Front. Genet. 2013, 4, 20. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  33. Chen, X.-G.; Jiang, X.; Gu, J.; Xu, M.; Wu, Y.; Deng, Y.; Zhang, C.; Bonizzoni, M.; Dermauw, W.; Vontas, J.; et al. Genome sequence of the Asian Tiger mosquito, Aedes albopictus, reveals insights into its biology, genetics, and evolution. Proc. Natl. Acad. Sci. USA 2015, 112, E5907–E5915. [Google Scholar] [CrossRef] [Green Version]
  34. Kim, Y.-K.; Kim, B.; Kim, V.N. Re-evaluation of the roles of DROSHA, Exportin 5, and DICER in microRNA biogenesis. Proc. Natl. Acad. Sci. USA 2016, 113, E1881–E1889. [Google Scholar] [CrossRef] [Green Version]
  35. Naganuma, M.; Tadakuma, H.; Tomari, Y. Single-molecule analysis of processive double-stranded RNA cleavage by Drosophila Dicer-2. Nat. Commun. 2021, 12, 4268. [Google Scholar] [CrossRef] [PubMed]
  36. Denli, A.M.; Tops, B.B.J.; Plasterk, R.H.A.; Ketting, R.F.; Hannon, G.J. Processing of primary microRNAs by the Microprocessor complex. Nature 2004, 432, 231–235. [Google Scholar] [CrossRef]
  37. Gregory, R.I.; Yan, K.-P.; Amuthan, G.; Chendrimada, T.; Doratotaj, B.; Cooch, N.; Shiekhattar, R. The Microprocessor complex mediates the genesis of microRNAs. Nature 2004, 432, 235–240. [Google Scholar] [CrossRef]
  38. Lee, Y.; Ahn, C.; Han, J.; Choi, H.; Kim, J.; Yim, J.; Lee, J.; Provost, P.; Rådmark, O.; Kim, S.; et al. The nuclear RNase III Drosha initiates microRNA processing. Nature 2003, 425, 415–419. [Google Scholar] [CrossRef]
  39. Berezikov, E. Evolution of microRNA diversity and regulation in animals. Nat. Rev. Genet. 2011, 12, 846–860. [Google Scholar] [CrossRef]
  40. Chen, K.; Rajewsky, N. The evolution of gene regulation by transcription factors and microRNAs. Nat. Rev. Genet. 2007, 8, 93–103. [Google Scholar] [CrossRef]
  41. Cabili, M.N.; Trapnell, C.; Goff, L.; Koziol, M.; Tazon-Vega, B.; Regev, A.; Rinn, J.L. Integrative annotation of human large intergenic noncoding RNAs reveals global properties and specific subclasses. Genes Dev. 2011, 25, 1915–1927. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  42. Lim, X.-N.; Shan, C.; Marzinek, J.K.; Dong, H.; Ng, T.S.; Ooi, J.S.G.; Fibriansah, G.; Wang, J.; Verma, C.S.; Bond, P.J.; et al. Molecular basis of dengue virus serotype 2 morphological switch from 29 °C to 37 °C. PLoS Pathog. 2019, 15, e1007996. [Google Scholar] [CrossRef] [Green Version]
  43. Hanley, J.P.; Tu, H.A.; Dragon, J.A.; Dickson, D.M.; del Rio-Guerra, R.; Tighe, S.W.; Eckstrom, K.M.; Selig, N.; Scarpino, S.V.; Whitehead, S.S.; et al. Immunotranscriptomic profiling the acute and clearance phases of a human challenge dengue virus serotype 2 infection model. Nat. Commun. 2021, 12, 3054. [Google Scholar] [CrossRef]
  44. Angleró-Rodríguez, Y.I.; MacLeod, H.J.; Kang, S.; Carlson, J.S.; Jupatanakul, N.; Dimopoulos, G. Aedes aegypti Molecular Responses to Zika Virus: Modulation of Infection by the Toll and Jak/Stat Immune Pathways and Virus Host Factors. Front. Microbiol. 2017, 8, 2050. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  45. Bonizzoni, M.; Dunn, W.A.; Campbell, C.L.; Olson, K.E.; Marinotti, O.; James, A.A. Complex Modulation of the Aedes aegypti Transcriptome in Response to Dengue Virus Infection. PLoS ONE 2012, 7, e50512. [Google Scholar] [CrossRef] [PubMed]
  46. Miesen, P.; Ivens, A.; Buck, A.H.; van Rij, R.P. Small RNA Profiling in Dengue Virus 2-Infected Aedes Mosquito Cells Reveals Viral piRNAs and Novel Host miRNAs. PLoS Negl. Trop. Dis. 2016, 10, e0004452. [Google Scholar] [CrossRef]
  47. Sim, S.; Jupatanakul, N.; Ramirez, J.L.; Kang, S.; Romero-Vivas, C.M.; Mohammed, H.; Dimopoulos, G. Transcriptomic Profiling of Diverse Aedes aegypti Strains Reveals Increased Basal-level Immune Activation in Dengue Virus-refractory Populations and Identifies Novel Virus-vector Molecular Interactions. PLoS Negl. Trop. Dis. 2013, 7, e2295. [Google Scholar] [CrossRef] [Green Version]
  48. Sim, S.; Dimopoulos, G. Dengue Virus Inhibits Immune Responses in Aedes aegypti Cells. PLoS ONE 2010, 5, e10678. [Google Scholar] [CrossRef] [Green Version]
  49. Tsujimoto, H.; Hanley, K.A.; Sundararajan, A.; Devitt, N.P.; Schilkey, F.D.; Hansen, I.A. Dengue virus serotype 2 infection alters midgut and carcass gene expression in the Asian tiger mosquito, Aedes albopictus. PLoS ONE 2017, 12, e0171345. [Google Scholar] [CrossRef] [Green Version]
  50. Zanini, F.; Pu, S.-Y.; Bekerman, E.; Einav, S.; Quake, S.R. Single-cell transcriptional dynamics of flavivirus infection. eLife 2018, 7, e32942. [Google Scholar] [CrossRef]
  51. Yung, C.-F.; Lee, K.-S.; Thein, T.-L.; Tan, L.-K.; Gan, V.C.; Wong, J.G.X.; Lye, D.C.; Ng, L.-C.; Leo, Y.-S. Dengue serotype-specific differences in clinical manifestation, laboratory parameters and risk of severe disease in adults, singapore. Am. J. Trop. Med. Hyg. 2015, 92, 999–1005. [Google Scholar] [CrossRef] [PubMed]
  52. Avila-Bonilla, R.G.; Yocupicio-Monroy, M.; Marchat, L.A.; Pérez-Ishiwara, D.G.; Cerecedo-Mercado, D.A.; del Ángel, R.M.; Salas-Benito, J.S. miR-927 has pro-viral effects during acute and persistent infection with dengue virus type 2 in C6/36 mosquito cells. J. Gen. Virol. 2020, 101, 825–839. [Google Scholar] [CrossRef] [PubMed]
  53. Sharma, C.P.; Ezzell, R.M.; Arnaout, M.A. Direct interaction of filamin (ABP-280) with the beta 2-integrin subunit CD18. J. Immunol. 1995, 154, 3461–3470. [Google Scholar]
  54. Edwards, D.N.; Towb, P.; Wasserman, S.A. An activity-dependent network of interactions links the Rel protein Dorsal with its cytoplasmic regulators. Development 1997, 124, 3855–3864. [Google Scholar] [CrossRef]
  55. Loo, D.T.; Kanner, S.B.; Aruffo, A. Filamin binds to the cytoplasmic domain of the beta1-integrin. Identification of amino acids responsible for this interaction. J. Biol. Chem. 1998, 273, 23304–23312. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  56. Koh, C.; Islam, M.N.; Ye, Y.H.; Chotiwan, N.; Graham, B.; Belisle, J.T.; Kouremenos, K.A.; Dayalan, S.; Tull, D.L.; Klatt, S.; et al. Dengue virus dominates lipid metabolism modulations in Wolbachia-coinfected Aedes aegypti. Commun. Biol. 2020, 3, 518. [Google Scholar] [CrossRef]
  57. Melo, C.F.O.R.; Delafiori, J.; Dabaja, M.Z.; de Oliveira, D.N.; Guerreiro, T.M.; Colombo, T.E.; Nogueira, M.L.; Proenca-Modena, J.L.; Catharino, R.R. The role of lipids in the inception, maintenance and complications of dengue virus infection. Sci. Rep. 2018, 8, 11826. [Google Scholar] [CrossRef]
  58. Zaitseva, E.; Yang, S.-T.; Melikov, K.; Pourmal, S.; Chernomordik, L.V. Dengue Virus Ensures Its Fusion in Late Endosomes Using Compartment-Specific Lipids. PLoS Pathog. 2010, 6, e1001131. [Google Scholar] [CrossRef] [Green Version]
  59. Mukhopadhyay, S.; Kuhn, R.J.; Rossmann, M.G. A structural perspective of the flavivirus life cycle. Nat. Rev. Microbiol. 2005, 3, 13–22. [Google Scholar] [CrossRef]
  60. Rodenhuis-Zybert, I.A.; Wilschut, J.; Smit, J.M. Dengue virus life cycle: Viral and host factors modulating infectivity. Cell. Mol. Life Sci. 2010, 67, 2773–2786. [Google Scholar] [CrossRef]
  61. Latanova, A.; Starodubova, E.; Karpov, V. Flaviviridae Nonstructural Proteins: The Role in Molecular Mechanisms of Triggering Inflammation. Viruses 2022, 14, 1808. [Google Scholar] [CrossRef] [PubMed]
Figure 1. Length distribution of small RNAs in DENV1-infected and uninfected C6/36 cells. Size distribution of small RNA reads in uninfected samples (a) and DENV1-infected samples (b).
Figure 1. Length distribution of small RNAs in DENV1-infected and uninfected C6/36 cells. Size distribution of small RNA reads in uninfected samples (a) and DENV1-infected samples (b).
Biology 11 01536 g001
Figure 2. Size distribution of mature miRNA (a) and precursor miRNA (b).
Figure 2. Size distribution of mature miRNA (a) and precursor miRNA (b).
Biology 11 01536 g002
Figure 3. (a) Expression distribution of lncRNAs and protein-coding genes in Ae. albopictus. (b) Distribution of specificity score of lncRNAs and protein-coding genes in Ae. albopictus.
Figure 3. (a) Expression distribution of lncRNAs and protein-coding genes in Ae. albopictus. (b) Distribution of specificity score of lncRNAs and protein-coding genes in Ae. albopictus.
Biology 11 01536 g003
Figure 4. Volcano plot of differentially expressed miRNAs upon DENV1 infection. Each dot represents miRNAs (novel and known) with their corresponding log2 fold change and −log10 p-value. Red dots represent significant p-value (<0.05), while black dots represent insignificant p-value (>0.05).
Figure 4. Volcano plot of differentially expressed miRNAs upon DENV1 infection. Each dot represents miRNAs (novel and known) with their corresponding log2 fold change and −log10 p-value. Red dots represent significant p-value (<0.05), while black dots represent insignificant p-value (>0.05).
Biology 11 01536 g004
Figure 5. GO bar graphs of target genes of the differentially expressed miRNAs. Only the top 10 GO terms were plotted in the bar graphs.
Figure 5. GO bar graphs of target genes of the differentially expressed miRNAs. Only the top 10 GO terms were plotted in the bar graphs.
Biology 11 01536 g005
Table 1. Mapping statistics of small RNA reads generated in this study.
Table 1. Mapping statistics of small RNA reads generated in this study.
SampleRaw ReadsClean Reads18–32 bp ReadsAlignment Rate (%)
DENV replicate 178933646768132215189388595.73
DENV replicate 277299313743931645085253496.02
DENV replicate 355658966541553834216381995.29
Uninfected replicate 154757609530822344147661696.62
Uninfected replicate 252985076515168474187728795.78
Uninfected replicate 358901861570089734199786895.31
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Azlan, A.; Yunus, M.A.; Abdul Halim, M.; Azzam, G. Revised Annotation and Characterization of Novel Aedes albopictus miRNAs and Their Potential Functions in Dengue Virus Infection. Biology 2022, 11, 1536. https://doi.org/10.3390/biology11101536

AMA Style

Azlan A, Yunus MA, Abdul Halim M, Azzam G. Revised Annotation and Characterization of Novel Aedes albopictus miRNAs and Their Potential Functions in Dengue Virus Infection. Biology. 2022; 11(10):1536. https://doi.org/10.3390/biology11101536

Chicago/Turabian Style

Azlan, Azali, Muhammad Amir Yunus, Mardani Abdul Halim, and Ghows Azzam. 2022. "Revised Annotation and Characterization of Novel Aedes albopictus miRNAs and Their Potential Functions in Dengue Virus Infection" Biology 11, no. 10: 1536. https://doi.org/10.3390/biology11101536

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