Analysis of Amblyomma americanum microRNAs in response to Ehrlichia chaffeensis infection and their potential role in vectorial capacity

Background: MicroRNAs (miRNAs) represent a subset of small noncoding RNAs and carry tremendous potential for regulating gene expression at the post-transcriptional level. They play pivotal roles in distinct cellular mechanisms including inhibition of bacterial, parasitic, and viral infections via immune response pathways. Intriguingly, pathogens have developed strategies to manipulate the host’s miRNA profile, fostering environments conducive to successful infection. Therefore, changes in an arthropod host’s miRNA profile in response to pathogen invasion could be critical in understanding host-pathogen dynamics. Additionally, this area of study could provide insights into discovering new targets for disease control and prevention. The main objective of the present study is to investigate the functional role of differentially expressed miRNAs upon Ehrlichia chaffeensis, a tick-borne pathogen, infection in tick vector, Amblyomma americanum.

Methods: Small RNA libraries from uninfected and E. chaffeensis-infected Am. americanum midgut and salivary gland tissues were prepared using the Illumina Truseq kit. Small RNA sequencing data was analyzed using miRDeep2 and sRNAtoolbox to identify novel and known miRNAs. The differentially expressed miRNAs were validated using a quantitative PCR assay. Furthermore, a miRNA inhibitor approach was used to determine the functional role of selected miRNA candidates.

Results: The sequencing of small RNA libraries generated >147 million raw reads in all four libraries and identified a total of >250 miRNAs across the four libraries. We identified 23 and 14 differentially expressed miRNAs in salivary glands, and midgut tissues infected with E. chaffeensis, respectively. Three differentially expressed miRNAs (miR-87, miR-750, and miR-275) were further characterized to determine their roles in pathogen infection. Inhibition of target miRNAs significantly decreased the E. chaffeensis load in tick tissues, which warrants more in-depth mechanistic studies.

Conclusions: The current study identified known and novel miRNAs and suggests that interfering with these miRNAs may impact the vectorial capacity of ticks to harbor Ehrlichia. This study identified several new miRNAs for future analysis of their functions in tick biology and tick-pathogen interaction studies.

Introduction

MicroRNAs (miRNAs) are non-coding RNAs with a size ranging from 18-25 nucleotides and play a significant role in post-transcriptional gene regulation (Bartel and Chen, 2004; Bartel, 2009). Latest studies have revealed the significance of miRNAs in arthropod immunity and host-pathogen interactions (Miesen et al., 2016; Momen-Heravi and Bala, 2018). In animals, miRNAs regulate post-transcriptional gene expression by binding to the 3’-untranslated region (3’-UTR), but there are also instances where the miRNA binds to the 5’-untranslated regions (5’-UTR), or coding regions. Perfect complementarity of 2-8 nucleotides at the 5’ end of the miRNA (seed region) is necessary for miRNA regulation, and the remaining sequence of miRNA might carry mismatches or bulges (Bartel, 2009; Rigoutsos, 2009; Schnall-Levin et al., 2010). The miRNA is transcribed as a primary miRNA transcript and processed by Drosha and Pasha into a pre-miRNA. The pre-miRNA is exported to the cytoplasm and processed by Dicer into a mature miRNA, which is then loaded into the microRNA-induced silencing complex (miRISC) and targets the complementary mRNA for degradation (Flynt et al., 2010; Asgari, 2018). Small non-coding RNAs (SncRNAs), including miRNAs, have shown tremendous potential in gene regulation at the post-transcriptional level in animals, plants, and arthropods, including ticks (Carrington and Ambros, 2003; Bartel, 2004; Griffiths-Jones et al., 2008; Lai, 2015). Although more than 800 tick species are present worldwide, ticks are underrepresented in available miRNA resources. Databases such as miRbase have 49 Ixodes scapularis miRNAs and 24 Rhipicephalus microplus miRNAs, while MirGeneDB 2.1 contains 64 Ixodes scapularis miRNAs (Fromm et al., 2022).

The lone-star tick (Amblyomma Americanum) is an aggressive human-biting tick species, a known vector of numerous disease-causing agents, including Ehrlichia chaffeensis, E. ewingii, heartland virus, Bourbon virus, Francisella tularensis, Borrelia lonestari (Sanchez-Vicente and Tokarz, 2023). Lone star tick bites are also known to cause a food allergy, Alpha-Gal Syndrome (AGS), or red meat allergy (Commins et al., 2011; Crispell et al., 2019; Sharma and Karim, 2021; Sharma et al., 2024). Am. americanum ticks are prevalent in the southern United States and have expanded their geographic range to the northeastern United States, and Canada (Stafford et al., 2018; Nelder et al., 2019). E. chaffeensis, a tick-borne Gram-negative obligatory intracellular bacterium, causes a severe flu-like febrile disease called human monocytic ehrlichiosis (HME), a prevalent life-threatening disease (Adams et al., 2017). Ehrlichiosis is an underreported tick-borne disease, and the pathogen infection of E. chaffeensis within the tick vector is a black box, and dynamics of vectorial capacity are largely unknown.

Given the contribution of miRNAs in numerous cellular processes, including development, immunity, and pathogen response in arthropods, the functional characterization of tick miRNAs in tick biology, and host-pathogen interactions remains to be investigated (Alvarez-Garcia and Miska, 2005; Saldaña et al., 2017). Several omics studies have characterized the time-dependent, tissues-dependent blood-meal and pathogen-induced differential gene expression in variety of tick species (Bartel, 2004; Anderson et al., 2008; Karim et al., 2011; Karim and Ribeiro, 2015; Popara et al., 2015; Villar et al., 2015; Antunes et al., 2019; Guizzo et al., 2022; Adegoke et al., 2023, Adegoke et al., 2024). However, studies addressing the role of pathogens in differentially modulating tick’s small RNAs are limited (Artigas-Jerónimo et al., 2019; Hermance et al., 2019; Ramasamy et al., 2020; Kumar et al., 2022). A handful of studies have investigated how miRNAs regulate the tick-pathogen interaction. These studies have shown that tick miRNAs promote the transmission of Anaplasma phagocytophilum and Powassan virus, thereby facilitating infection establishment in Ixodes scapularis (Artigas-Jerónimo et al., 2019; Hermance et al., 2019; Ramasamy et al., 2020). An elegant study demonstrated the tick miRNA-mediated regulation of vertebrate host gene expression at the tick-host interface (Hackenberg et al., 2017). However, regarding miRNA-mediated gene expression, the lone star tick (Am. americanum) is an underrepresented tick species, and there is an urgent need to investigate how E. chaffeensis utilizes tick miRNAs to promote its survival and persistence within the tick vector. This knowledge gap is of great concern, especially considering the increasing threat of the lone star tick to public health significance.

Understanding the molecular interactions between tick vectors and E. chaffeensis, and the characterization of differentially regulated tick miRNAs are needed to develop new approaches to combat tick-borne infections. Identification of new and novel miRNAs by using new sequencing platforms opened up a new avenue of research (Artigas-Jerónimo et al., 2019; Kumar et al., 2022; Luo et al., 2023). In this work, a small RNA sequencing approach was utilized to identify miRNAs induced by E. chaffeensis in A. americanum tissues, and differentially expressed miRNAs were functionally characterized using a miRNA inhibitor approach.

Materials and methods

Ethics statement

All animal experiments were performed in strict accordance with the recommendations in the NIH Guide for the Care and Use of Laboratory Animals. The Institutional Animal Care and Use Committee of the University of Southern Mississippi approved the protocol for blood feeding of field-collected ticks (protocol # 15101501.3).

Ticks and tissue dissections

Adult Am. americanum ticks, both infected and uninfected were prepared as described (Karim et al., 2012). Briefly, lab-grown E. chaffeensis (Arkansas strain) was microinjected in engorged nymphs and 1x DMEM media was injected as a control cohort for uninfected ticks (Karim et al., 2012; Budachetri et al., 2020: Budachetri et al., 2022; Adegoke et al., 2024, Adegoke et al., 2024). The injected engorged nymphs were kept in sterile vials with a piece of filter paper. The engorged nymphs molted into either unfed male or female adult ticks within 2 months. The ticks were maintained under standard conditions as outlined by Patrick and Hair (1975). A qPCR assay was used to determine the infection level of E. chaffeensis in freshly molted ticks (Dunphy et al., 2014). Uninfected and E. chaffeensis adult ticks were infested on a rabbit for blood-feeding (see the experimental design Figure 1). Partially blood-fed uninfected and pathogen-infected female ticks were removed from the rabbit, and tick tissues were dissected as described earlier (Karim et al., 2012), and dissected tissues were directly stored in Trizol (Life Technologies, Carlsbad, CA, USA). Samples were kept at −80°C until use.

Figure 1

Figure 1 Schematic workflow of microRNA study in lone star tick tissues.

RNA extraction and small RNA sequencing

RNA was extracted using Trizol RNA extraction methods from ten pooled tick tissue samples (uninfected, E. chaffeensis-infected midguts, and salivary glands). The integrity of the extracted RNA was determined using the standard spectrophotometric method as described earlier (Kumar et al., 2022). Small RNA library synthesis and small RNA sequencing were outsourced to the University of Mississippi Medical Center (UMMC) Molecular and Genomics Core laboratory. Briefly, four small RNA libraries (uninfected midguts, E. chaffeensis infected midguts, uninfected salivary glands, and E. chaffeensis infected salivary glands) were prepared using the Illumina Truseq Kit according to the manufacturer’s guidelines. Short adapter oligonucleotides were ligated to each end of the small RNAs. Following this, a cDNA copy was made with reverse transcriptase, and polymerase chain reaction (PCR) was utilized to incorporate sample-specific barcodes and Illumina sequencing adapters. The final concentration of all Next Generation Sequencing (NGS) libraries was determined with a Qubit fluorometric assay and a DNA 1000 high-sensitivity chip on an Agilent 2100 Bioanalyzer was used to assess the DNA fragment size for each library followed a purification step by polyacrylamide gel-electrophoresis. The sample libraries were pooled and sequenced on an Illumina Next Seq 500 (single end 36 bases) using TruSeq SBS kit v3 (Illumina) according to the manufacturer’s protocols.

Data analysis

There are currently two computational strategies available for miRNA profiling: 1) the genome-based strategy, which maps small RNA-seq reads to a reference genome and evaluates sequences that generate the characteristic hairpin structure of miRNA precursors (Bortolomeazzi et al., 2019), and 2) the machine-learning-based strategy, where biogenesis features of sequences are extracted based on the available miRNA sequences in microRNA databases such as miRBase (Kozomara and Griffiths-Jones, 2014) and from the analysis of miRNA duplex structures (Vitsios et al., 2017). During the analysis of this data, the Am. americanum genome sequence was not available. Therefore, microRNA data was analyzed by the smallRNAtoolbox webserver (Aparicio-Puerta et al., 2022). Recently, Chou et al. (2023) sequenced the genome of the Am. americanum assembled the long-read sequenced genome and provided a file of a partially annotated genome. The miRDeep2 software package version 2.0.0.8 (Friedländer et al., 2012) was used to predict the novel and known miRNAs in all tick tissue samples using the partially annotated genome of Am. americanum (Supplementary Tables S1-S5). sRNAtoolbox algorithm has not been optimized for novel microRNA discovery or may have limitations in accurately identifying and characterizing novel microRNAs. Since the genome of the tick Am. americanum has not been fully annotated for noncoding RNAs such as rRNAs, tRNAs, snRNAs, snoRNAs, etc., the percentage abundance of these non-coding RNAs in our data has not been determined.

Bioinformatics

Bowtie1 was used for indexing the fasta file of the Amblyomma americanum genome. miRDeep2 was used for profiling novel and known miRNAs, and sRNAtoolbox webserver (specific modules – sRNAbench, sRNAde) was used to determine differentially expressed miRNAs. As mentioned above, our small RNA data were analyzed by using a smallRNAtoolbox webserver (Aparicio-Puerta et al., 2022). Briefly, after an initial sequencing quality control step in FastQC (http://www.bioinformatics.babraham.ac.uk/projects/fastqc), preprocessing, mapping, and annotation were mainly conducted in sRNAbench (Aparicio-Puerta et al., 2022) with customized scripts as necessary. Simply, the obtained sequence reads were 36 nucleotides (nt) in length, but many small RNAs were between 27 and 33 nt. When forcing the detection of at least 10 nt of adapter as a typically used minimum length, only RNA molecules of up to 26 nt can be resolved (read length plus minimum adapter length). Therefore, to detect all small RNAs <36 nucleotides, we implemented iterative adapter detection and trimming. First, the adapter was detected in the whole read, and, if not found, it was then searched using iteratively shorter minimum adapter lengths at the 3′ end. After adapter trimming, the reads collapsed into unique reads followed by read count assignment, i.e., counting the number of times that each unique read was sequenced.

Differential expression and normalization

Our experimental design resulted in several possible comparisons: (i) uninfected versus E.chaffeensis infected salivary glands (SGs), (ii) uninfected versus E. chaffeensis infected midgut (MGs), (iii) uninfected MGs versus uninfected SGs, and (iv) E. chaffeensis-infected MGs versus E. chaffeensis-infected SGs. An edgeR tool in the sRNAtool box was used to determine differential miRNAs expression between uninfected and E. chaffeensis-infected tick tissues (Robinson et al., 2010). Briefly, using the differential expression module of sRNAtoolbox (i.e. sRNAde), we generated an expression matrix with the raw read counts for input into edgeR to obtain differential miRNA expression. The edgeR normalizes the data using the trimmed mean of M-values (TMM) method. We also generated an expression matrix with reads per million (RPM)-normalized expression values using the “single assignment” procedure in sRNAbench. As a result, each read mapping multiple times was only assigned once to the miRNA with the highest expression and only affected reads mapping to several different reference sequences, i.e., normally miRNA sequences from the same family. The RPM values were obtained by dividing the read count of a given miRNA by the total number of reads mapped to the miRNA library.

Validation of differentially expressed miRNAs

All miRNAs that were differentially expressed in small RNA sequencing data were validated by qRT-PCR in Am. americanum tick tissues. Mir-X miRNA qRT-PCR TB Green kit (Takara BIO, San Jose, CA, USA) was used for cDNA synthesis and miRNA expression analysis. This kit includes the Mir-X miRNA First-Strand Synthesis Kit, which transforms RNA into complementary DNA (cDNA) and enables the quantification of particular miRNA sequences through real-time PCR. Briefly, utilizing a single-tube method, RNA molecules undergo polyadenylation and reverse transcription by the action of poly(A) polymerase and SMART^® MMLV Reverse Transcriptase, both components of the mRQ Enzyme Mix are supplied with the kit. Subsequently, real-time qPCR is conducted using the TB Green Advantage^® qPCR Premix and mRQ 3’ Primer, in combination with miRNA-specific 5’ primers, to quantify specific miRNA expression. Primers used are listed in Supplementary Table S6. Conditions used for qRT-PCR were initial denaturation of 95°C for 10 mins, then 40 cycles of 95°C for 5 secs, and 60°C for 20 secs.

MicroRNA inhibition assay

We selected three microRNAs, including Aam-miR-87, Aam-miR-750, and Aam-miR-275 for further characterization. A miRNA inhibitor approach works by sterically blocking specific miRNA functions using an oligonucleotide that complements the mature miRNA target (Lennox et al., 2013). These inhibitors were designed and synthesized by Integrated DNA Technologies (IDT, Coralville, IA), which also synthesized non-target negative controls alongside the specific inhibitors. Three groups of E. chaffeensis-infected ticks were injected with miRNA inhibitors for each selected miRNA (Aam-miR-87, Aam-miR-750, and Aam-miR-275), while a fourth group was injected with non-target negative controls. Each group contained 25 female ticks. Using a published study (Ramasamy et al., 2020), we selected a 1.05 nanomoles dose for this assay. Briefly, 1.05 nanomoles of each microRNA inhibitor and negative control were injected into each female tick within their respective groups. The injected ticks were allowed a 48-hour recovery period along with non-treated males (25 females/15 males) in a laboratory incubator, maintained at a temperature range of 23 ± 2˚C with a humidity level of 95%, and under a light cycle of 14 hours of light and 10 hours of darkness. After the recovery period of 48 hrs, all groups of ticks were infested on a sheep in separate stockinet cells to blood feed. The partially blood-fed ticks were removed from the sheep, and tick tissues (MGs, and SGs) were dissected within 4 hrs for downstream analysis. RNA was extracted from each tissue using the Trizol method as described previously (Kumar et al., 2022). In a qPCR assay, the E. chaffeensis infection in individual tick tissues was assessed using disulfide bond formation (dsb) gene primers (Dunphy et al., 2014). The level of miRNA inhibition was determined by using QRT-PCR primers designed by the miRprimer2 algorithm (Busk, 2014), in conjunction with the Mir-X™ miRNA qRT-PCR TB Green^® Kit from TAKARA (San Jose, CA, USA). After confirming miRNA inhibition in tick tissues, the relative E. chaffeensis load compared to the non-target negative control was quantified in the respective tick tissues.

Results and discussion

Profile characteristics of small RNA libraries

The small RNA sequencing yielded a total of >147 million raw reads, comprised of >32 million from uninfected midgut, >32million from the E. chaffeensis-infected midgut, >40million from uninfected salivary glands, and >41 million from the E. chaffeensis-infected salivary glands. Following the adapter trimming process and subsequent removal of short reads (those ≤20 nucleotides (nt) in length), the small RNA reads left for downstream analysis were >41 million from E. chaffeensis-infected samples and >36 million from uninfected samples. The distribution of read length presents an indication of the types of small RNAs present in both E. chaffeensis-infected and uninfected tick tissues (midgut and salivary gland). Sequencing stats including raw reads, adapter cleaned reads, reads in analysis, quality filter reads have been provided in Supplementary Table S7, while processing stats of reads (in percentage) and raw miRNA summary (number of detected miRNAs and its precursors) have been provided correspondingly as Supplementary Table S8 and Supplementary Table S9. Both types of samples, E. chaffeensis-infected and uninfected, exhibited two main peaks at 22 nt (representing miRNAs/siRNAs) and 29 nt (representing piRNAs) (Figure 2). Notably, PIWI-interacting RNAs (piRNAs) constitute a class of small RNAs, which usually vary in size from 26 to 31 nucleotides (Iwasaki et al., 2015; Santos et al., 2023). The piRNAs associate with PIWI proteins, which belong to the Argonaute family of proteins and are active in the testes of mammals. These RNAs play a crucial role in germ cell and stem cell development in invertebrates. One of their primary functions is the silencing of transposable elements (TEs) to protect genomic integrity (Aravin et al., 2007; Brennecke et al., 2007, Brennecke et al., 2008). piRNAs are passed down from the female germline to progeny, ensuring the stability of the genome across generations. When males harboring a specific family of transposon elements (TEs) mate with females devoid of these elements, the females lack the necessary complementary piRNAs to defend their genome resulting in an overabundance of transposons. Such proliferation could potentially lead to sterility, a situation termed hybrid dysgenesis (Erwin et al., 2015). This concept of hybrid dysgenesis holds potential for tick control.

Figure 2

Figure 2 (A) Small RNA sequence length distribution in uninfected (clean) and E. chaffeensis-infected tick tissues (SG, MG). MicroRNAs are twenty-two (22) nucleotides in length. (B) Pie-chart distribution of small RNA reads. There are mainly two characteristic peaks of tick small RNAs (miRNAs/siRNAs-22 nt, piRNAs-29 nt). Aa- Am. americanum, PF-Partially fed (2 days), SG- Salivary glands, MG- Midgut, nt-nucleotide. Ech – E. chaffeensis infected.

In silico mapping of E. chaffeensis infected small RNA sequences to E. chaffeensis Arkansas strain genome

Upon performing in silico mapping of small RNA reads infected with E. chaffeensis to the genome of the E. chaffeensis Arkansas strain (GCF_000013145.1_ASM1314v1_genomic.fna), we detected 1,413 E. chaffeensis sequences in the midgut out of a total of 2,492,851 reads. In addition, we found 3,185 E. chaffeensis sequences in the salivary glands from a total of 19,103,069 reads. From these results, it appears that our samples were infected with E. chaffeensis, and the infection level was low, a hallmark of Ehrlichia infection within the tick vector (Kennedy et al., 2021).

Differentially expressed microRNAs in E. chaffeensis infected tick salivary glands

Our small RNA sequencing analysis identified 360 microRNAs including known and predicted ones in uninfected and E. chaffeensis-infected tick tissues (Supplementary Data-Supplementary Table S1). The list of differentially expressed miRNAs in E. chaffeensis-infected salivary glands (SGs) includes 18 upregulated miRNAs, Aam-miR-5322, Aam-miR-1, Aam-miR-750, Aam-miR-993, Aam-miR-5307, Aam-miR-87, Aam-miR-5315, Aam-miR-133, Aam-let-7, Aam-miR-3931, Aam-miR-263a, Aam-miR-8, Aam-miR-5305, Aam-bantam, Aam-miR-279, Aam-miR-5306, Aam-miR-276, Aam-miR-315) and 5 downregulated miRNAs Aam-miR-12, Aam-miR-2a, Aam-miR-10, Aam-miR-7, Aam-miR-285 (Figure 3).

Figure 3

Figure 3 In silico differential expression of predicted microRNAs in E. chaffeensis infected partially fed salivary gland relative to partially fed clean salivary glands. EdgeR was used for differential expression analysis. miRNAs with a Log2 fold-change expression > |0.8|, p-value < 0.05 were considered significantly differentially expressed.

Differentially expressed microRNAs in E. chaffeensis infected tick midgut

Likewise, in the E. chaffeensis infected tick midgut (MG), seven miRNAs were upregulated in Ehrlichia infection, Aam-miR-5309, Aam-miR-79, Aam-miR-1, Aam-miR-275, Aam-miR-315, Aam-miR-278, Aam-miR-263a and 7 downregulated miRNAs, Aam-miR-153, Aam-miR-12, Aam-miR-7, Aam-miR-124, Aam-miR-5314, Aam-miR-5310, Aam-miR-285 (Figure 4). Heat map representation of DE miRNAs (Supplementary Figure S1, Supplementary Figure S2) has also been provided as Supplementary Information). All these miRNAs are listed in miRbase (22.1) and have also been identified in Ixodes scapularis or Rhipicephalus microplus tick species. It is necessary to investigate the roles of these differentially expressed miRNAs in E. chaffeensis-infected tick tissues. It is crucial to gain a deeper understanding through these differentially expressed miRNAs of how E. chaffeensis manipulates the expression of tick microRNAs to ensure its survival, persistence, and transmission. This knowledge holds immense potential in the development of innovative tools to effectively block the transmission of tick-borne pathogens.

Figure 4

Figure 4 In silico differential expression of predicted microRNAs in E. chaffeensis infected partially fed midgut relative to the partially fed clean midgut. EdgeR was used for differential expression analysis. miRNAs with a Log2 fold-change expression > |0.8|, p-value < 0.05 were considered significantly differentially expressed.

In this study, we found miR-1 upregulated in E. chaffeensis infected tick tissues, i.e. SG, and MG in comparison to uninfected ones. It is noteworthy that miR-1 often exhibits increased expression during pathogen infections. miR-1 belongs to a conserved family, which includes miR-7 and miR-34, and is conserved across various organisms such as fruit flies, shrimps, and humans (Takane et al., 2010). It participates in analogous pathways in these organisms, including development and apoptosis, and its upregulation has also been identified during stressful conditions (Huang et al., 2012). In mosquitoes, miR-1 is similarly upregulated during Plasmodium infection (Huang and Zhang, 2012), and facilitates West Nile Virus infection (Hussain et al., 2012). For Bombyx mori, the Nucleo-polyhedrosis virus (NPV) releases miR-1 inside the host to regulate its target RAN (exportin 5, co-factor), a player in miRNA biogenesis responsible for transporting pre-miRNA from the nucleus to the cytoplasm (Singh et al., 2012). Moreover, during Listeria infection in macrophages, miR-1 stimulates IFN-γ-dependent activation of the innate immune response (Xu et al., 2019). Based on these findings, we speculate that miR-1 could be upregulated in E. chaffeensis-infected tick tissues to trigger an immune response against E. chaffeensis.

The upregulation of miR-87 in Ehrlichia-infected SGs points out its putative role in ticks’ innate immune response. Earlier work on other arthropods, including Manduca sexta and Aedes albopictus hint at its potential role in interfering with innate immunity, particularly IMD and toll receptor signaling pathways (Zhang et al., 2014; Liu et al., 2015; Avila-Bonilla et al., 2017). Its in silico predicted targets in Aedes albopictus include TOLL pathway signaling Ser/Thr. Kinase, Toll-like receptor TOLL 1A, Class A Scavenger receptor with Se-Protease domain, and Galectin (Liu et al., 2015), as well as putative TLR 5b (Avila-Bonilla et al., 2017). In Manduca sexta, its predicted target is FADD, an adaptor protein involved in DISC formation (Zhang et al., 2014). Based on the published studies, we hypothesize that Ehrlichia infection differentially regulates the miR-87 to inhibit the Toll pathway for its survival in the tick vector.

The upregulation of miR-bantam in Ehrlichia-infected SGs also suggests its putative role in the pathogen infection of E. chaffeensis. Interestingly, miR-bantam, a conserved microRNA, exhibits high expression levels in insects. For example, it is highly expressed in Drosophila and participates in several key cellular processes such as cell proliferation, apoptosis, development, and the circadian rhythm. In Drosophila, miR-bantam serves two major roles: inhibiting apoptosis by down-regulating the apoptotic gene hid (Brennecke et al., 2003) and promoting cell proliferation by targeting genes like mad (Robins et al., 2005). The inhibition of apoptosis might act as a survival mechanism for E.chaffeensis within tick salivary gland cells. Our hypothesis warrants in-depth future studies to determine the role of miR-bantam in the vectorial capacity of the tick vector.

Another miRNA candidate, mir-79 was also significantly upregulated Ehrlichia infected MG. Earlier studies have described its role in various pathways, including immunity, cell differentiation, neurogenesis, and apoptosis. It has also been implicated in cancer and disease caused by viral infection (Yuva-Aydemir et al., 2011; Pederson et al., 2013; Seddiki et al., 2013; Ouyang et al., 2015; Dong et al., 2017; Artigas-Jerónimo et al., 2019). It is known that mir-79 disrupts the JNK pathway by targeting its component genes pvr (CG8222) and puc (CG7850) (Fullaondo and Lee, 2012). The JNK pathway is an immune response pathway against Gram-negative bacterial pathogens (Bond and Foley, 2009). In the current study, mir-79 is upregulated in the E. chaffeensis-infected midgut. Similarly, in ticks infected with Anaplasma phagocytophilum (a Gram-negative bacterial pathogen), mir-79 was found upregulated, thereby facilitating infection by targeting the Roundabout protein 2 pathway (Robo2) (Artigas-Jerónimo et al., 2019). The upregulation of mir-79 could potentially be a mechanism used by E. chaffeensis to evade the tick’s immune system.

Our in-silico data reveals the downregulation of miR-5310 in the midgut of E. chaffeensis-infected ticks. MiR-5310, a tick-specific miRNA (Barrero et al., 2011), may play a role in tick feeding, as it was shown to be downregulated in Rhipicephalus microplus tick larvae exposed to host odor without being allowed to feed (Barrero et al., 2011). Another study showed its downregulation in Anaplasma phagocytophilum-infected nymphs compared to unfed, uninfected nymphs (Ramasamy et al., 2020). Previous research has also shown the modulation of signaling events upon A. phagocytophilum infection of ticks (Neelakanta et al., 2010; Sultana et al., 2010; Khanal et al., 2017; Taank et al., 2017; Turck et al., 2019; Ramasamy et al., 2020). Thus, in the case of the E. chaffeensis-infected tick midgut in the present study, miR-5310’s speculated role might be to modulate signaling events to protect E. chaffeensis. Its potential role could also be in tick feeding, as shown by its downregulation in Rhipicephalus microplus tick larvae exposed to host odor without being allowed to feed.

In our data, miR-133 is upregulated in the salivary glands of E. chaffeensis-infected ticks. According to a recent study, the infection of ticks with the pathogen Anaplasma phagocytophilum results in the downregulation of tick microRNA-133 (miR-133), leading to the induction of the Ixodes scapularis organic anion transporting polypeptide (isoatp4056) gene expression, which is critical for the pathogen’s survival within the vector and its transmission to the vertebrate host (Ramasamy et al., 2020). Therefore, the upregulation of miR-133 in E. chaffeensis-infected tick salivary glands in our study might suggest the downregulation of organic anion transporting polypeptide (isoatp4056) gene expression, potentially inhibiting E. chaffeensis survival and transmission. This is a hypothetical explanation for the observed miR-133 upregulation, and further investigation is necessary for confirmation.

Let-7, an evolutionarily conserved microRNA in bilateral animals, plays a role in developmental regulation, such as molting and metamorphosis in arthropods, and can disrupt innate immunity by targeting the antimicrobial peptide diptericin (Pasquinelli et al., 2000; Carrington and Ambros, 2003; Garbuzov and Tatar, 2010; Hertel et al., 2012; Ling et al., 2014). Recent studies have also suggested its role in the molting of Hyalomma asiaticum ticks by targeting the ecdysteroid receptor (ECR), a part of the 20E signaling pathway (Sempere et al., 2003; Wu et al., 2019). In this study, let-7’s upregulation in tick salivary glands implies a possible role in targeting the antimicrobial peptide diptericin, potentially allowing E. chaffeensis to evade innate immunity. This may represent a mechanism for E. chaffeensis’s survival and successful transmission in tick salivary glands, but further investigation is required. It should be noted that diptericin inhibits Gram-negative bacteria by disrupting membrane integrity.

In this study, miR-275 is upregulated in the midgut of E. chaffeensis-infected ticks. MiR-275 directly targets and positively regulates the sarco/endoplasmic reticulum Ca2+ adenosine triphosphatase (SERCA), an active player in transporting Ca2+ from the cytosol to the sarco/endoplasmic reticulum (ER) in mosquito guts (Zhao et al., 2017). The transportation of Ca++ from the cytoplasm to the ER is required for the spreading process of Ehrlichia canis (Alves et al., 2014), suggesting that E. chaffeensis may modulate tick machinery via upregulation of miR-275. However, a follow-up study is necessary for confirmation. It’s also worth noting that miR-275 was found to be crucial for blood digestion and egg development in the mosquito Aedes aegypti (Bryant et al., 2010).

In our data set, miR-750 is upregulated in E. chaffeensis-infected tick salivary glands. Past studies have suggested its role in innate immunity, hormone signaling, and stress response (Nunes et al., 2013; Rebijith et al., 2016; Queiroz et al., 2020; Kanoksinwuttipong et al., 2022). A recent study indicated that upregulated miR-750 suppresses its target, the sarcoplasmic calcium-binding protein (Scp), and inhibits apoptosis, thus contributing to pathogen propagation (Kanoksinwuttipong et al., 2022). Given these previous studies, the possible role of miR-750 in inhibiting apoptosis and promoting E. chaffeensis propagation in tick salivary glands can be speculated. This could represent a mechanism by which E. chaffeensis avoids cellular apoptosis and propagates for effective transmission in the tick salivary glands. The roles of all the differentially expressed miRNAs mentioned above are listed in Table 1, along with other necessary details.

Table 1

Table 1 List of differentially expressed microRNAs detected in Ehrlichia chaffeensis infected tick tissues (salivary gland, midgut) and their putative roles.

miR-8 is upregulated in E. chaffeensis-infected tick salivary glands. Its role in innate immunity is unknown so far, but our results have indicated its possible role against E. chaffeensis infection, but further work is required for validation. It is a conserved miRNA, and previous studies have shown its role in development and reproduction. Its expression was found up-regulated in Aedes aegypti during the pupation stage, with its highest expression levels observed in the mid-pupal period. Previous work by Bryant et al. (2010) revealed the upregulation of miR-8 in the fat body of blood-fed female mosquitoes, suggesting a potential regulatory function in the reproduction of Ae. aegypti. In contrast to Ae. aegypti, miR-8 shows abundant expression at different developmental stages of Anopheles stephensi (Feng et al., 2018) and is equally expressed in uninfected and infected Ae. albopictus saliva following CHIKV infection (Maharaj et al., 2015).

miR-263a was also found upregulated in E. chaffeensis infected tick salivary glands indicating its possible role in either innate immunity or transmission, which needs to be revealed by further work. Previous studies have shown its role in development. It was found highly expressed in uninfected and infected Ae. aegypti saliva (Maharaj et al., 2015), and is amongst the most highly expressed miRNAs across developmental stages in many mosquito species (Hu et al., 2015).

miR-12 was found downregulated in E. chaffeensis infected tick salivary glands indicating its possible role in activating immune pathways, and its downregulation might have a probable effect on Ehrlichia survival, and therefore successful transmission. Further work is required to validate this hypothesis. Although its immunological role is unknown in previous studies. Although previous work has shown its role in affecting Wolbachia density in mosquitoes (Osei-Amo et al., 2012). The preferential expression of miR-12 in Anopheles gambiae occurs in the thorax of both males and females, predominantly in midguts and twice as much in their heads (Winter et al., 2007), and its targets are DNA replication licensing factor (MCM6) and monocarboxylate transporter (MCT1) genes, as validated in A.aegypti, by which it affects Wolbachia density in host cells (Osei-Amo et al., 2012). Our in silico data indicates that miR-279 is upregulated in the salivary glands of ticks infected with Ehrlichia chaffeensis, suggesting it may have a role in the pathogen’s survival or its transmission. These hypotheses, still in the realm of speculation, warrant further study to be substantiated. Additionally, a recent study has proposed that miR-279 might influence the resistance to the insecticide deltamethrin. It does this by regulating the expression of its target gene CYP325BB1, which codes for the enzyme cytochrome P450 325bb1, in the mosquito species Culex pipiens pallens (Li et al., 2021). Given that the functions of microRNAs are conserved, this research work might offer substantial insights into the mechanisms driving acaricide resistance, which would be essential in developing new and effective strategies for tick control in the future.

Validation of in silico differentially expressed microRNAs by qRT-PCR

The expression levels of differentially expressed miRNAs were validated using qRT-PCR assays on E. chaffeensis-infected and uninfected tick tissues (Figure 5). The qRT-PCR patterns of the differentially expressed miRNAs were consistent with the next-generation sequencing (NGS) results for the majority of evaluated miRNAs. However, inconsistencies between the NGS and qRT-PCR data patterns were detected. These discrepancies may arise from the different methodologies used to quantify miRNA expression (Saldaña et al., 2017). The comparison of the differentially regulated microRNAs upon blood-feeding or infection of pathogen using Illumina small RNA sequencing and qRT-PCR assay showed inconsistent expression data (Supplementary Figure S3). These expression inconsistenciesbetween the two methodologies can be attributed to several factors. Illumina small RNA sequencing and qRT-PCR are fundamentally different approaches with varying sensitivities, specificities, and biases. Small RNA sequencing is a high-throughput method that can capture a broad range of microRNAs, while qRT-PCR is a targeted approach that amplifies and quantifies specific microRNAs of interest. Variations in library preparation, sequencing chemistry, primer design, and amplification efficiencies can contribute to differences in the quantification of microRNA expression levels between the two techniques. Additionally, differences in normalization strategies can alos affect the quantification of microRNA expression levels and lead to discrepancies in the observed expression patterns between the two techniques. Normalization methods such as total read count normalization, spike-in controls, or reference gene normalization may introduce variability in the expression measurements. Factors such as RNA degradation, PCR amplification biases, contamination, and technical variability can impact the accuracy and reproducibility of microRNA expression measurements.

Figure 5

Figure 5 qPCR validation of differentially expressed miRNAs (in silico) in E. chaffeensis infected and partially fed tick tissues. (A) Midgut (B) Salivary glands. Expression of miRNAs was normalized with clean and partially fed tick tissues (indicated as 1 on the y-axis). At least three biological replicates were used in each of the experiments.

miRNA inhibition in tick tissues reduced E. chaffeensis load

Three differentially expressed miRNAs Aam-miR-87, Aam-miR-750, and Aam-miR-275 were selected for the miRNA inhibition assay in tick tissues based on their putative role in tick immune responses (Table 1). As depicted in Figure 2 and Figure 3, Aam-miR-87 and Aam-miR-750 were found to be upregulated in the salivary glands of Ehrlichia chaffeensis-infected ticks, while Aam-miR-275 was upregulated in the midgut of ticks infected with Ehrlichia chaffeensis. The results of our miRNA inhibitory experiments revealed that suppressing these microRNAs individually reduced the E. chaffeensis load in tick tissues suggesting a role for miRNAs in tick immunity (Figure 6), although further validation is necessary. It is noteworthy to point out the pathogen infection levels in the samples. All samples were infected with E. chaffeensis, but the infection level of E. chaffeensis was not synchronized (Figures 6B, C). It is a technical challenge to have a synchronized number of E. chaffeensis in each sample. Subsequent studies will explore the targets of these microRNAs to gain a deeper understanding of their importance in pathways essential for the survival of E. chaffeensis. Additional functional studies will examine how miRNAs and their specific targets impact pathways that affect tick vector competence. Further investigation is also needed to explore additional differentially expressed miRNAs identified in this study, which could potentially offer valuable insights for preventing E. chaffeensis infection.

Figure 6

Figure 6 microRNA inhibition reduces E. chaffeensis load in tick tissues. (A) microRNA inhibition in tick tissues (MG and SG). miR-87 and miR-750 were inhibited by ~75% and ~80% in tick salivary glands while miR-275 was inhibited by ~63% in tick midgut expression of microRNAs in negative control tick tissues were given a normalized fold expression value of 1, as represented by the dashed line. Tick microRNA inhibition resulted in E. chaffeensis load reduction (dsb gene) in (B) tick salivary glands and (C) midguts. At least three biological replicates were used in each of the experiments. Statistically significant change (P<0.05) is indicated by an asterisk (*).

Limitations and conclusion

This study highlights the differential expression of miRNAs in E. chaffeensis-infected tick tissues, which could significantly influence the survival, colonization, and transmission of E. chaffeensis. Additionally, these miRNAs may play a role in the tick’s immune response against the pathogen. Selected microRNAs miR-87, miR-750, and miR-275 have shown promising results against E. chaffeensis survival or colonization in ticks. Further investigation of other differentially expressed miRNAs through miRNA inhibitory experiments is needed to explore these aspects. These tick-specific and E. chaffeensis specific differentially expressed miRNAs could provide potent avenues for treating or inhibiting E. chaffeensis.

Here it is noteworthy to mention that analyzing microRNA sequencing data (from Illumina i.e. short reads) within the context of a long-read sequenced genome may pose challenges in read mapping, alignment, and data integration due to the differences in read lengths and sequencing technologies. Aligning short Illumina reads to a long-read sequenced genome could result in decreased mapping efficiency, especially in regions with structural variations or repetitive elements. This may impact the accuracy of microRNA expression profiling and annotation. While long-read sequences may better capture genomic heterogeneity and structural variations compared to short reads, and provide valuable insights into genome architecture, they may also introduce challenges in accurately identifying and characterizing microRNAs, particularly in regions of high complexity or variation. Long-read sequencing technologies may have significantly higher error rates, which can pose challenges in accurately reconstructing the genome, mapping microRNA sequences, and distinguishing true variations from sequencing errors.

Data availability statement

The raw small RNA sequences were deposited into the NCBI Sequence Read Archive (SRA) repository under the BioProject ID PRJNA992656.

Ethics statement

The animal study was approved by The University of Southern Mississippi IACUC. The study was conducted in accordance with the local legislation and institutional requirements.

Author contributions

DK: Conceptualization, Formal analysis, Investigation, Methodology, Validation, Visualization, Writing – original draft, Writing – review & editing. KB: Investigation, Methodology, Resources, Writing – review & editing. YR: Investigation, Resources, Writing – review & editing. SK: Conceptualization, Formal analysis, Funding acquisition, Investigation, Project administration, Resources, Supervision, Visualization, Writing – original draft, Writing – review & editing.

Funding

The author(s) declare financial support was received for the research, authorship, and/or publication of this article. This research was principally supported by the NIH NIAID Awards #R15AI167013, #R21AI175885, and # R01AI163857. We thank the UMMC Molecular and Genomics facility, supported by the NIHNIGMS (#P20GM103476 & P20GM144041). The funders played no role in the study design, data collection, analysis, publication, decision, or manuscript preparation.

Acknowledgments

This work utilized the Magnolia High-Performance Computing at USM (http://magnolia.usm.edu).

Conflict of interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

The author(s) declared that they were an editorial board member of Frontiers, at the time of submission. This had no impact on the peer review process and the final decision.

Publisher’s note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Supplementary material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fcimb.2024.1427562/full#supplementary-material

References

Adams, D. A., Thomas, K. R., Jajosky, R. A., Foster, L., Baroi, G., Sharp, P., et al. (2017). Summary of notifiable infectious diseases and conditions - United States 2015. MMWR Morb Mortal Wkly Rep. 64, 1–143. doi: 10.15585/mmwr.mm6453a1

CrossRef Full Text | Google Scholar

Adegoke, A., Hanson, J., Smith, R., Karim, S. (2023). Ehrlichia chaffeensis co-opts phagocytic hemocytes for systemic dissemination in the Lone Star tick, Amblyomma americanum. Preprint. bioRxiv. doi: 10.1101/2023.08.17.553720

CrossRef Full Text | Google Scholar

Adegoke, A., Ribeiro, J. M. C., Smith, R. C., Karim, S. (2024). Tick innate immune responses to hematophagy and Ehrlichia infection at single-cell resolution. Front. Immunol. 14. doi: 10.3389/fimmu.2023.1305976

PubMed Abstract | CrossRef Full Text | Google Scholar

Alvarez-Garcia, I., Miska, E. A. (2005). MicroRNA functions in animal development and human disease. Development 132, 4653–4662. doi: 10.1242/dev.02073

PubMed Abstract | CrossRef Full Text | Google Scholar

Alves, R. N., Levenhagen, M. A., Levenhagen, M. M., Rieck, S. E., Labruna, M. B., Beletti, M. E. (2014). The spreading process of Ehrlichiacanis in macrophages is dependent on actin cytoskeleton, calcium and iron influx and lysosomal evasion. Vet. Microbiol. 168, 442–446. doi: 10.1016/j.vetmic.2013.11.030

PubMed Abstract | CrossRef Full Text | Google Scholar

Anderson, J. M., Sonenshine, D. E., Valenzuela, J. G. (2008). Exploring the mialome of ticks: an annotated catalogue of midgut transcripts from the hard tick, Dermacentor variabilis (Acari: Ixodidae). BMC Genomics 9, 552. doi: 10.1186/1471-2164-9-552

PubMed Abstract | CrossRef Full Text | Google Scholar

Antunes, S., Couto, J., Ferrolho, J., Sanches, G. S., Merino Charrez, J. O., de la Cruz Hernandez, N., et al. (2019). Transcriptome and Proteome Response of Rhipicephalus annulatus Tick Vector to Babesia bigemina Infection. Front. Physiol. 10. doi: 10.3389/fphys.2019.00318

PubMed Abstract | CrossRef Full Text | Google Scholar

Aparicio-Puerta, E., Gómez-Martín, C., Giannoukakos, S., Medina, J. M., Scheepbouwer, C., García-Moreno, A., et al. (2022). sRNAbench and sRNAtoolbox 2022 update: accurate miRNA and sncRNA profiling for model and non-model organisms. Nucleic Acids Res. 50, W710–W717. doi: 10.1093/nar/gkac363

PubMed Abstract | CrossRef Full Text | Google Scholar

Aravin, A. A., Hannon, G. J., Brennecke, J. (2007). The Piwi-piRNA pathway provides an adaptive defense in the transposon arms race. Science 318 (5851), 761–764. doi: 10.1126/science.1146484

PubMed Abstract | CrossRef Full Text | Google Scholar

Artigas-Jerónimo, S., Alberdi, P., Villar Rayo, M., Cabezas-Cruz, A., Prados, P. J. E., Mateos-Hernández, L., et al. (2019). Anaplasma phagocytophilum modifies tick cell microRNA expression and upregulates isc-mir-79 to facilitate infection by targeting the Roundabout protein 2 pathway. Sci. Rep. 9 (1), 9073. doi: 10.1038/s41598-019-45658-2

PubMed Abstract | CrossRef Full Text | Google Scholar

Asgari, S. (2018). “Chapter two - microRNAs as regulators of insect host–pathogen interactions and immunity,” in Advances in insect physiology, vol. 55 . Ed. Smagghe, G. (Cambridge, MA, USA: Academic Press), 19–45.

Google Scholar

Avila-Bonilla, R. G., Yocupicio-Monroy, M., Marchat, L. A., De Nova-Ocampo, M. A., Del Angel, R. M., Salas-Benito, J. S. (2017). Analysis of the miRNA profile in C6/36 cells persistently infected with dengue virus type 2. Virus Res. 232, 139–151. doi: 10.1016/j.virusres.2017.03.005

PubMed Abstract | CrossRef Full Text | Google Scholar

Barrero, R. A., Keeble-Gagnère, G., Zhang, B., Moolhuijzen, P., Ikeo, K., Tateno, Y., et al. (2011). Evolutionary conserved microRNAs are ubiquitously expressed compared to tick-specific miRNAs in the cattle tick Rhipicephalus (Boophilus) microplus. BMC Genomics. 12, 328. doi: 10.1186/1471-2164-12-328

PubMed Abstract | CrossRef Full Text | Google Scholar

Bartel, D. P. (2004). MicroRNAs: genomics, biogenesis, mechanism, and function. Cell 116, 281–297. doi: 10.1016/S0092-8674(04)00045-5

PubMed Abstract | CrossRef Full Text | Google Scholar

Bartel, D. P. (2009). MicroRNAs: target recognition and regulatory functions. Cell. 136, 215–233. doi: 10.1016/j.cell.2009.01.002

PubMed Abstract | CrossRef Full Text | Google Scholar

Bartel, D. P., Chen, C. Z. (2004). Micromanagers of gene expression: The Potentially Widespread Influence of Metazoan MicroRNAs. Nat. Rev. Genet. 5, 396–400. doi: 10.1038/nrg1328

PubMed Abstract | CrossRef Full Text | Google Scholar

Bond, D., Foley, E. (2009). A quantitative RNAi screen for JNK modifiers identifies Pvr as a novel regulator of Drosophila immune signaling. PLoSPathog. 5, e1000655. doi: 10.1371/journal.ppat.1000655

CrossRef Full Text | Google Scholar

Bortolomeazzi, M., Gaffo, E., Bortoluzzi, S. (2019). A survey of software tools for microRNA discovery and characterization using RNA-seq. Brief Bioinf. 20, 918–930. doi: 10.1093/bib/bbx148

CrossRef Full Text | Google Scholar

Brennecke, J., Aravin, A. A., Stark, A., Dus, M., Kellis, M., Sachidanandam, R., et al. (2007). Discrete small RNA-generating loci as master regulators of transposon activity in Drosophila. Cell 128 (6), 1089–1103. doi: 10.1016/j.cell.2007.01.043

PubMed Abstract | CrossRef Full Text | Google Scholar

Brennecke, J., Hipfner, D. R., Stark, A., Russell, R. B., Cohen, S. M. (2003). Bantam encodes a developmentally regulated microRNA that controls cell proliferation and regulates the proapoptotic gene hid in drosophila. Cell 113, 25–36. doi: 10.1016/S0092-8674(03)00231-9

PubMed Abstract | CrossRef Full Text | Google Scholar

Brennecke, J., Malone, C. D., Aravin, A. A., Sachidanandam, R., Stark, A., Hannon, G. J. (2008). An epigenetic role for maternally inherited piRNAs in transposon silencing. Science (New York, N.Y.) 322 (5906), 1387–1392. doi: 10.1126/science.1165171

PubMed Abstract | CrossRef Full Text | Google Scholar

Bryant, B., Macdonald, W., Raikhel, A. S. (2010). microRNA miR-275 is indispensable for blood digestion and egg development in the mosquito Aedes aEgypti. Proc. Natl. Acad. Sci. U S A. 107, 22391–22398. doi: 10.1073/pnas.1016230107

PubMed Abstract | CrossRef Full Text | Google Scholar

Budachetri, K., Lin, M., Chien, R. C., Zhang, W., Brock, G. N., Rikihisa, Y. (2022). Efficacy and immune correlates of OMP-1B and virB2-4 vaccines for protection of dogs from tick transmission of ehrlichia chaffeensis. mBio. 13, e0214022. doi: 10.1128/mbio.02140-22

PubMed Abstract | CrossRef Full Text | Google Scholar

Budachetri, K., Teymournejad, O., Lin, M., Yan, Q., Mestres-Villanueva, M., Brock, G. N., et al. (2020). An entry-triggering protein of ehrlichia is a new vaccine candidate against tick-borne human monocytic ehrlichiosis. mBio 11, 1–13. doi: 10.1128/mBio.00895-20

CrossRef Full Text | Google Scholar

Busk, P. K. (2014). A tool for design of primers for microRNA-specific quantitative RT-qPCR. BMC Bioinf. 15, 29. doi: 10.1186/1471-2105-15-29

CrossRef Full Text | Google Scholar

Carrington, J. C., Ambros, V. (2003). Role of microRNAs in plant and animal development. Science 301, 336–338. doi: 10.1126/science.1085242

PubMed Abstract | CrossRef Full Text | Google Scholar

Chou, S., Poskanzer, K. E., Rollins, M., Thuy-Boun, P. S. (2023). De novo assembly of a long-read Amblyomma americanum tick genome. doi: 10.57844/arcadia-9b6j-q683

CrossRef Full Text | Google Scholar

Commins, S. P., James, H. R., Kelly, L. A., Pochan, S. L., Workman, L. J., Perzanowski, M. S., et al. (2011). The relevance of tick bites to the production of IgE antibodies to the mammalian oligosaccharide galactose-α-1,3-galactose. J. Allergy Clin. Immunol. 127, 1286–93.e6. doi: 10.1016/j.jaci.2011.02.019

PubMed Abstract | CrossRef Full Text | Google Scholar

Crispell, G., Commins, S. P., Archer-Hartman, S. A., Choudhary, S., Dharmarajan, G., Azadi, P., et al. (2019). Discovery of alpha-gal-containing antigens in North American tick species believed to induce red meat allergy. Front. Immunol. 10, 1056. doi: 10.3389/fimmu.2019.01056

PubMed Abstract | CrossRef Full Text | Google Scholar

Dong, C., Sun, X., Guan, Z., Zhang, M., Duan, M. (2017). Modulation of infuenza A virus replication by microRNA-9 through targeting MCPIP1. J. Med. Virol. 89, 41–48. doi: 10.1002/jmv.24604

PubMed Abstract | CrossRef Full Text | Google Scholar

Dunphy, P. S., Luo, T., McBride, J. W. (2014). Ehrlichia chaffeensis exploits host SUMOylation pathways to mediate effector-host interactions and promote intracellular survival. Infection Immun. 82, 4154–4168. doi: 10.1128/IAI.01984-14

CrossRef Full Text | Google Scholar

Erwin, A. A., Galdos, M. A., Wickersheim, M. L., Harrison, C. C., Marr, K. D., Colicchio, J. M., et al. (2015). piRNAs are associated with diverse transgenerational effects on gene and transposon expression in a hybrid dysgenic syndrome of D. virilis. PLoS Genet. 11, e1005332. doi: 10.1371/journal.pgen.1005332

PubMed Abstract | CrossRef Full Text | Google Scholar

Feng, X., Zhou, S., Wang, J., Hu, W. (2018). microRNA profiles and functions in mosquitoes. PLoS Negl. Trop. Dis. 12, e0006463. doi: 10.1371/journal.pntd.0006463

PubMed Abstract | CrossRef Full Text | Google Scholar

Flynt, A. S., Greimann, J. C., Chung, W. J., Lima, C. D., Lai, E. C. (2010). MicroRNA biogenesis via splicing and exosome-mediated trimming in Drosophila. Mol. Cell 38, 900–907. doi: 10.1016/j.molcel.2010.06.014

PubMed Abstract | CrossRef Full Text | Google Scholar

Friedländer, M. R., Mackowiak, S. D., Li, N., Chen, W., Rajewsky, N. (2012). miRDeep2 accurately identifies known and hundreds of novel microRNA genes in seven animal clades. Nucleic Acids Res. 40, 37–52. doi: 10.1093/nar/gkr688

PubMed Abstract | CrossRef Full Text | Google Scholar

Fromm, B., Høye, E., Domanska, D., Zhong, X., Aparicio-Puerta, E., Ovchinnikov, V., et al. (2022). MirGeneDB 2.1: toward a complete sampling of all major animal phyla. Nucleic Acids Res. 50, D204–D210. doi: 10.1093/nar/gkab1101

PubMed Abstract | CrossRef Full Text | Google Scholar

Fullaondo, A., Lee, S. Y. (2012). Identification of putative miRNA involved in Drosophila melanogaster immune response. Dev. Comp. Immunol. 36, 267–273. doi: 10.1016/j.dci.2011.03.034

PubMed Abstract | CrossRef Full Text | Google Scholar

Garbuzov, A., Tatar, M. (2010). Hormonal regulation of Drosophila microRNA let-7 and miR-125 that target innate immunity. Fly (Austin). 4, 306–311. doi: 10.4161/fly.4.4.13008

PubMed Abstract | CrossRef Full Text | Google Scholar

Griffiths-Jones, S., Saini, H. K., van Dongen, S., Enright, A. J. (2008). miRBase: tools for microRNA genomics. Nucleic Acids Res. 36, D154–D158. doi: 10.1093/nar/gkm952

PubMed Abstract | CrossRef Full Text | Google Scholar

Guizzo, M. G., Dolezelikova, K., Neupane, S., Frantova, H., Hrbatova, A., Pafco, B., et al. (2022). Characterization and manipulation of the bacterial community in the midgut of Ixodes ricinus. Parasit Vectors. 15, 248. doi: 10.1186/s13071-022-05362-z

PubMed Abstract | CrossRef Full Text | Google Scholar

Hackenberg, M., Langenberger, D., Schwarz, A., Erhart, J., Kotsyfakis, M. (2017). In silico target network analysis of de novo-discovered, tick saliva-specific microRNAs reveals important combinatorial effects in their interference with vertebrate host physiology. RNA. 23, 1259–1269. doi: 10.1261/rna.061168.117

PubMed Abstract | CrossRef Full Text | Google Scholar

Hermance, M. E., Widen, S. G., Wood, T. G., Thangamani, S. (2019). Ixodes scapularis salivary gland microRNAs are differentially expressed during Powassan virus transmission. Sci. Rep. 9, 1–17. doi: 10.1038/s41598-019-49572-5

PubMed Abstract | CrossRef Full Text | Google Scholar

Hertel, J., Bartschat, S., Wintsche, A., Otto, C., Stadler, P. F. (2012). Evolution of the let-7 microRNA family. RNA Biol. 9, 231–241. doi: 10.4161/rna.18974

PubMed Abstract | CrossRef Full Text | Google Scholar

Hu, W., Criscione, F., Liang, S., Tu, Z. (2015). MicroRNAs of two medically important mosquito species: Aedes aEgypti and Anopheles stephensi. Insect Mol. Biol. 24, 240–252. doi: 10.1111/imb.12152

PubMed Abstract | CrossRef Full Text | Google Scholar

Huang, T., Xu, D., Zhang, X. (2012). Characterization of host microRNAs that respond to DNA virus infection in a crustacean. BMC Genomics 13, 1–10. doi: 10.1186/1471-2164-13-159

PubMed Abstract | CrossRef Full Text | Google Scholar

Huang, T., Zhang, X. (2012). Functional analysis of a crustacean microRNA in host-virus interactions. J. Virol. 86, 12997–13004. doi: 10.1128/JVI.01702-12

PubMed Abstract | CrossRef Full Text | Google Scholar

Hussain, M., Torres, S., Schnettler, E., Funk, A., Grundhoff, A., Pijlman, G. P., et al. (2012). West Nile virus encodes a microRNA-like small RNA in the 30 untranslated region which up-regulates GATA4 mRNA and facilitates virus replication in mosquito cells. Nucleic Acids Res. 40, 2210–2223. doi: 10.1093/nar/gkr848

PubMed Abstract | CrossRef Full Text | Google Scholar

Iwasaki, Y. W., Siomi, M. C., Siomi, H. (2015). PIWI-interacting RNA: its biogenesis and functions. Annu. Rev. Biochem. 84, 405–433. doi: 10.1146/annurev-biochem-060614-034258

PubMed Abstract | CrossRef Full Text | Google Scholar

Kanoksinwuttipong, N., Jaree, P., Somboonwiwat, K. (2022). Shrimp pmo-miR-750 regulates the expression of sarcoplasmic calcium-binding protein facilitating virus infection in Penaeus monodon. Fish Shellfish Immunol. 129, 74–84. doi: 10.1016/j.fsi.2022.08.046

PubMed Abstract | CrossRef Full Text | Google Scholar

Karim, S., Browning, R., Ali, L., Truhett, R. (2012). Laboratory-infected ehrlichiachaffeensis female adult amblyommaamericanum salivary glands reveal differential gene expression. J. Med. Entomology 49, 547–554. doi: 10.1603/ME11214

CrossRef Full Text | Google Scholar

Karim, S., Ribeiro, J. M. (2015). An Insight into the Sialome of the Lone Star Tick, Amblyomma americanum, with a Glimpse on Its Time Dependent Gene Expression. PLoS One 10, e0131292. doi: 10.1371/journal.pone.0131292

PubMed Abstract | CrossRef Full Text | Google Scholar

Karim, S., Singh, P., Ribeiro, J. M. (2011). A deep insight into the sialotranscriptome of the gulf coast tick, Amblyomma maculatum. PLoS One 6, e28525. doi: 10.1371/journal.pone.0028525

PubMed Abstract | CrossRef Full Text | Google Scholar

Kennedy, A. C., BCE1, Marshall, E. (2021). Lone star ticks (Amblyomma americanum):: an emerging threat in delaware. Dela J. Public Health 7, 66–71. doi: 10.32481/djph.2021.01.013

PubMed Abstract | CrossRef Full Text | Google Scholar

Khanal, S., Sultana, H., Catravas, J. D., Carlyon, J. A., Neelakanta, G. (2017). Anaplasmaphagocytophilum infection modulates expression of megakaryocyte cell cycle genes through phosphatidylinositol-3-kinase signaling. PLoS One 12, e0182898. doi: 10.1371/journal.pone.0182898

PubMed Abstract | CrossRef Full Text | Google Scholar

Kozomara, A., Griffiths-Jones, S. (2014). miRBase: annotating high confidence microRNAs using deep sequencing data. Nucleic Acids Res. 42, D68–D73. doi: 10.1093/nar/gkt1181

PubMed Abstract | CrossRef Full Text | Google Scholar

Kumar, D., Downs, L. P., Embers, M., Flynt, A. S., Karim, S. (2022). Identification of microRNAs in the lyme disease vector. Ixodes scapularis. Int. J. Mol. Sci. 23, 5565. doi: 10.3390/ijms23105565

CrossRef Full Text | Google Scholar

Lai, E. C. (2015). Two decades of miRNA biology: lessons and challenges. RNA. 21, 675–677. doi: 10.1261/rna.051193.115

PubMed Abstract | CrossRef Full Text | Google Scholar

Lennox, K. A., Owczarzy, R., Thomas, D. M., Walder, J. A., Behlke, M. A. (2013). Improved performance of anti-miRNA oligonucleotides using a novel non-nucleotide modifier. Mol. Ther. Nucleic Acids 2, e117. doi: 10.1038/mtna.2013.46

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, X., Hu, S., Zhang, H., Yin, H., Wang, H., Zhou, D., et al. (2021). MiR-279-3p regulates deltamethrin resistance through CYP325BB1 in Culex pipiens pallens. Parasites Vectors 14, 528. doi: 10.1186/s13071-021-05033-5

PubMed Abstract | CrossRef Full Text | Google Scholar

Ling, L., Ge, X., Li, Z., Zeng, B., Xu, J., Aslam, A. F. M., et al. (2014). MicroRNA Let-7 regulates molting and metamorphosis in the silkworm, Bombyx mori. Insect Biochem. Mol. Biol. 53, 13–21. doi: 10.1016/j.ibmb.2014.06.011

PubMed Abstract | CrossRef Full Text | Google Scholar

Liu, W., Hao, Z., Huang, L., Chen, L., Wei, Q., Cai, L., et al. (2017). Comparative expression profile of microRNAs in Anopheles anthropophagus midgut after blood-feeding and Plasmodium infection. Parasites Vectors. 10 (1), 86. doi: 10.1186/s13071-017-2027-6

PubMed Abstract | CrossRef Full Text | Google Scholar

Liu, Y., Zhou, Y., Wu, J., Zheng, P., Li, Y., Zheng, X., et al. (2015). The expression profile of Aedes albopictus miRNAs is altered by dengue virus serotype-2 infection. Cell Biosci. 5, 1–11. doi: 10.1186/s13578-015-0009-y

PubMed Abstract | CrossRef Full Text | Google Scholar

Lucas, K. J., Roy, S., Ha, J., Gervaise, A. L., Kokoza, V. A., Raikhel, A. S. (2015). MicroRNA-8 targets the Wingless signaling pathway in the female mosquito fat body to regulate reproductive processes. Proc. Natl. Acad. Sci. U S A. 112, 1440–1445. doi: 10.1073/pnas.1424408112

PubMed Abstract | CrossRef Full Text | Google Scholar

Luo, Y., Peng, L., Shan, W., Sun, M., Luo, L., Liang, W. (2023). Machine learning in the development of targeting microRNAs in human disease. Front. Genet. 13. doi: 10.3389/fgene.2022.1088189

CrossRef Full Text | Google Scholar

Maharaj, P. D., Widen, S. G., Huang, J., Wood, T. G., Thangamani, S. (2015). Discovery of mosquito saliva microRNAs during CHIKV infection. PLoS Negl. Trop. Dis. 9, e0003386. doi: 10.1371/journal.pntd.0003386

PubMed Abstract | CrossRef Full Text | Google Scholar

Miesen, P., Ivens, A., Buck, A. H., van Rij, R. P. (2016). Small RNA Profiling in Dengue Virus 2-Infected Aedes Mosquito Cells Reveals Viral piRNAs and Novel Host miRNAs. PLoSNegl Trop. Dis. 10, e0004452. doi: 10.1371/journal.pntd.0004452

CrossRef Full Text | Google Scholar

Momen-Heravi, F., Bala, S. (2018). miRNA regulation of innate immunity. J. Leukoc. Biol. doi: 10.1002/JLB.3MIR1117-459R

PubMed Abstract | CrossRef Full Text | Google Scholar

Neelakanta, G., Sultana, H., Fish, D., Anderson, J. F., Fikrig, E. (2010). Anaplasmaphagocytophilum induces Ixodes scapularis ticks to express an antifreeze glycoprotein gene that enhances their survival in the cold. J. Clin. Invest. 120, 3179–3190. doi: 10.1172/JCI42868

PubMed Abstract | CrossRef Full Text | Google Scholar

Nelder, M. P., Russell, C. B., Clow, K. M., Johnson, S., Weese, J. S., Cronin, K., et al. (2019). Occurrence and distribution of Ambylomma americanum as determined by passive surveillance in Ontario, Canada, (1999-2016). Ticks Tick Borne Dis. 10 (1), 146–155. doi: 10.1016/j.ttbdis.2018.10.001

PubMed Abstract | CrossRef Full Text | Google Scholar

Nunes, F. M. F., Ihle, K. E., Mutti, N. S., Simões, Z. L. P., Amdam, G. V. (2013). The gene vitellogenin affects microRNA regulation in honey bee (Apis mellifera) fat body and brain. J. Exp. Biol. 216, 3724–3732. doi: 10.1242/jeb.089243

PubMed Abstract | CrossRef Full Text | Google Scholar

Osei-Amo, S., Hussain, M., O'Neill, S. L., Asgari, S. (2012). Wolbachia-induced aae-miR-12 miRNA negatively regulates the expression of MCT1 and MCM6 genes in Wolbachia-infected mosquito cell line. PLoS One 7, e50049. doi: 10.1371/journal.pone.0050049

PubMed Abstract | CrossRef Full Text | Google Scholar

Ouyang, W., Wang, Y. S., Du, X. N., Liu, H. J., Zhang, H. B. (2015). gga-miR-9* inhibits IFN production in antiviral innate immunity by targeting interferon regulatory factor 2 to promote IBDV replication. Vet. Microbiol. 178, 41–49. doi: 10.1016/j.vetmic.2015.04.023

PubMed Abstract | CrossRef Full Text | Google Scholar

Pasquinelli, A. E., Reinhart, B. J., Slack, F., Martindale, M. Q., Kuroda, M. I., Maller, B., et al. (2000). Conservation of the sequence and temporal expression of let-7 heterochronic regulatory RNA. Nature. 408, 86–89. doi: 10.1038/35040556

PubMed Abstract | CrossRef Full Text | Google Scholar

Patrick, C. D., Hair, J. A. (1975). Laboratory rearing procedures and equipment for multi-host ticks (Acarina: Ixodidae). J. Med. Entomol. 12, 389–390. doi: 10.1093/jmedent/12.3.389

PubMed Abstract | CrossRef Full Text | Google Scholar

Pedersen, M. E., Snieckute, G., Kagias, K., Nehammer, C., Multhaupt, H. A., Couchman, J. R., et al. (2013). An epidermal microRNA regulates neuronal migration through control of the cellular glycosylation state. Science 341 (6152), 1404–1408. doi: 10.1126/science.1242528

PubMed Abstract | CrossRef Full Text | Google Scholar

Popara, M., Villar, M., de la Fuente, J. (2015). Proteomics characterization of tick-host-pathogen interactions. Methods Mol. Biol. 1247, 513–527. doi: 10.1007/978-1-4939-2004-4_34

PubMed Abstract | CrossRef Full Text | Google Scholar

Queiroz, F. R., Portilho, L. G., Jeremias, W. J., Babá, É. H., do Amaral, L. R., Silva, L. M., et al. (2020). Deep sequencing of small RNAs reveals the repertoire of miRNAs and piRNAs in Biomphalaria glabrata. Mem Inst Oswaldo Cruz. 115, e190498. doi: 10.1590/0074-02760190498

PubMed Abstract | CrossRef Full Text | Google Scholar

Ramasamy, E., Taank, V., Anderson, J. F., Sultana, H., Neelakanta, G. (2020). Repression of tick microRNA-133 induces organic anion transporting polypeptide expression critical for Anaplasmaphagocytophilum survival in the vector and transmission to the vertebrate host. PLoS Genet. 16, e1008856. doi: 10.1371/journal.pgen.1008856

PubMed Abstract | CrossRef Full Text | Google Scholar

Rebijith, K. B., Asokan, R., Hande, H. R., Krishna Kumar, N. K. (2016). The First Report of miRNAs from a Thysanopteran Insect, Thrips palmiKarny Using High-Throughput Sequencing. PLoS One 11, e0163635. doi: 10.1371/journal.pone.0163635

PubMed Abstract | CrossRef Full Text | Google Scholar

Rigoutsos, I. (2009). New tricks for animal microRNAS: targeting of amino acid coding regions at conserved and nonconserved sites. Cancer Res. 69, 3245–3248. doi: 10.1158/0008-5472.CAN-09-0352

PubMed Abstract | CrossRef Full Text | Google Scholar

Robins, H., Li, Y., Padgett, R. W. (2005). Incorporating structure to predict microRNA targets. Proc. Natl. Acad. Sci. U S A. 102, 4006–4009. doi: 10.1073/pnas.0500775102

PubMed Abstract | CrossRef Full Text | Google Scholar

Robinson, M. D., McCarthy, D. J., Smyth, G. K. (2010). edgeR: a Bioconductor package for differential expression analysis of digital gene expression data. Bioinformatics 26, 139–140. doi: 10.1093/bioinformatics/btp616

PubMed Abstract | CrossRef Full Text | Google Scholar

Saldaña, M., Etebari, K., Hart, C. E., Widen, S. G., Wood, T. G., Thangamani, S., et al. (2017). Zika virus alters the microRNA expression profile and elicits an RNAi response in Aedes aEgypti mosquitoes. PLoSNegl. Trop. Dis. 11, e0005760. doi: 10.1371/journal.pntd.0005760

CrossRef Full Text | Google Scholar

Sanchez-Vicente, S., Tokarz, R. (2023). Tick-borne co-infections: Challenges in molecular and serologic diagnoses. Pathogens. 12, 1371. doi: 10.3390/pathogens12111371

PubMed Abstract | CrossRef Full Text | Google Scholar

Santos, D., Feng, M., Kolliopoulou, A., Taning, C. N. T., Sun, J., Swevers, L. (2023). What are the functional roles of piwi proteins and piRNAs in insects? Insects 14, 187. doi: 10.3390/insects14020187

PubMed Abstract | CrossRef Full Text | Google Scholar

Schnall-Levin, M., Zhao, Y., Perrimon, N., Berger, B. (2010). Conserved microRNA targeting in Drosophila is as widespread in coding regions as in 3'UTRs. Proc. Natl. Acad. Sci. U.S.A. 107, 15751–15756. doi: 10.1073/pnas.1006172107

PubMed Abstract | CrossRef Full Text | Google Scholar

Seddiki, N., Phetsouphanh, C., Swaminathan, S., Xu, Y., Rao, S., Li, J., et al. (2013). Te microRNA-9/B-lymphocyte-induced maturation protein-1/IL-2 axis is diferentially regulated in progressive HIV infection. Eur. J. Immunol. 43 (2), 510–520. doi: 10.1002/eji.201242695

PubMed Abstract | CrossRef Full Text | Google Scholar

Sempere, L. F., Sokol, N. S., Dubrovsky, E. B., Berger, E. M., Ambros, V. (2003). Temporal regulation of microRNA expression in Drosophila melanogaster mediated by hormonal signals and broad-Complex gene activity. Dev. Biol. 259, 9–18. doi: 10.1016/S0012-1606(03)00208-2

PubMed Abstract | CrossRef Full Text | Google Scholar

Sharma, S. R., Choudhary, S. K., Vorobiov, J., Commins, S. P., Karim, S. (2024). Tick bite-induced alpha-gal syndrome and immunologic responses in an alpha-gal deficient murine model. Front Immunol. 14, 1336883. doi: 10.3389/fimmu.2023.1336883

PubMed Abstract | CrossRef Full Text | Google Scholar

Sharma, S. R., Karim, S. (2021). Tick saliva and the alpha-gal syndrome: Finding a needle in a haystack. Front. Cell Infect. Microbiol. 11, 680264. doi: 10.3389/fcimb.2021.680264

PubMed Abstract | CrossRef Full Text | Google Scholar

Singh, C. P., Singh, J., Nagaraju, J. (2012). A baculovirus-encoded microRNA (miRNA) suppresses its host miRNA biogenesis by regulating the exportin-5 cofactor ran. J. Virol. doi: 10.1128/JVI.00064-12

PubMed Abstract | CrossRef Full Text | Google Scholar

Stafford, K. C., 3rd, Molaei, G., Little, E. A. H., Paddock, C. D., Karpathy, S. E., Labonte, A. M. (2018). Distribution and establishment of the lone star tick in connecticut and implications for range expansion and public health. J. Med. Entomol. 55, 1561–1568. doi: 10.1093/jme/tjy115

PubMed Abstract | CrossRef Full Text | Google Scholar

Sultana, H., Neelakanta, G., Kantor, F. S., Malawista, S. E., Fish, D., Montgomery, R. R., et al. (2010). Anaplasmaphagocytophilum induces actin phosphorylation to selectively regulate gene transcription in Ixodes scapularis ticks. J. Exp. Med. 207, 1727–1743. doi: 10.1084/jem.20100276

PubMed Abstract | CrossRef Full Text | Google Scholar

Taank, V., Dutta, S., Dasgupta, A., Steeves, T. K., Fish, D., Anderson, J. F., et al. (2017). Human rickettsial pathogen modulates arthropod organic anion transporting polypeptide and tryptophan pathway for its survival in ticks. Sci. Rep. 7 (1), 13256. doi: 10.1038/s41598-017-13559-x

PubMed Abstract | CrossRef Full Text | Google Scholar

Takane, K., Fujishima, K., Watanabe, Y., Sato, A., Saito, N., Tomita, M., et al. (2010). Computational prediction and experimental validation of evolutionarily conserved microRNA target genes in bilaterian animals. BMC Genomics 11, 101. doi: 10.1186/1471-2164-11-101

PubMed Abstract | CrossRef Full Text | Google Scholar

Turck, J. W., Taank, V., Neelakanta, G., Sultana, H. (2019). Ixodes scapularis Src tyrosine kinase facilitates Anaplasmaphagocytophilum survival in its arthropod vector. Ticks Tick-Borne Diseases. 10, 838–847. doi: 10.1016/j.ttbdis.2019.04.002

CrossRef Full Text | Google Scholar

Villar, M., Ayllon, N., Alberdi, P., Moreno, A., Moreno, M., Tobes, R., et al. (2015). Integrated metabolomics, transcriptomics and proteomics identifies metabolic pathways affected by anaplasmaphagocytophilum infection in tick cells. Mol. Cell. Proteomics: MCP. 14, 3154–3172. doi: 10.1074/mcp.M115.051938

CrossRef Full Text | Google Scholar

Vitsios, D. M., Kentepozidou, E., Quintais, L., Benito-Gutiérrez, E., van Dongen, S., Davis, M. P., et al. (2017). Mirnovo: genome-free prediction of microRNAs from small RNA sequencing data and single-cells using decision forests. Nucleic Acids Res. 45, e177–e177. doi: 10.1093/nar/gkx836

PubMed Abstract | CrossRef Full Text | Google Scholar

Winter, F., Edaye, S., Hüttenhofer, A., Brunel, C. (2007). Anopheles Gambiae miRNAs as actors of defence reaction against Plasmodium invasion. Nucleic Acids Res. 35, 6953–6962. doi: 10.1093/nar/gkm686

PubMed Abstract | CrossRef Full Text | Google Scholar

Wu, F., Luo, J., Chen, Z., Ren, Q., Xiao, R., Liu, W., et al. (2019). MicroRNA let-7 regulates the expression of ecdysteroid receptor (ECR) in Hyalommaasiaticum (Acari: Ixodidae) ticks. Parasit Vectors. 12 (1), 235. doi: 10.1186/s13071-019-3488-6

PubMed Abstract | CrossRef Full Text | Google Scholar

Xu, H., Jiang, Y., Xu, X., Su, X., Liu, Y., Ma, Y., et al. (2019). Inducible degradation of lncRNA Sros1 promotes IFN-gamma-mediated activation of innate immune responses by stabilizing Stat1 mRNA. Nat. Immunol. 20, 1621–1630. doi: 10.1038/s41590-019-0542-7

PubMed Abstract | CrossRef Full Text | Google Scholar

Yuva-Aydemir, Y., Simkin, A., Gascon, E., Gao, F. B. (2011). MicroRNA-9: functional evolution of a conserved small regulatory RNA. RNA Biol. 8, 557–564. doi: 10.4161/rna.8.4.16019

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhang, X., Zheng, Y., Jagadeeswaran, G., Ren, R., Sunkar, R., Jiang, H. (2014). Identification of conserved and novel microRNAs in Manduca sexta and their possible roles in the expression regulation of immunity-related genes. Insect. Biochem. Mol. Biol. 47, 12–22. doi: 10.1016/j.ibmb.2014.01.008

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhao, B., Lucas, K. J., Saha, T. T., Ha, J., Ling, L., Kokoza, V. A., et al. (2017). MicroRNA-275 targets sarco/endoplasmic reticulum Ca²⁺ adenosine triphosphatase (SERCA) to control key functions in the mosquito gut. PLoS Genet. 13, e1006943. doi: 10.1371/journal.pgen.1006943

PubMed Abstract | CrossRef Full Text | Google Scholar