Open Access

Comparative expression profile of microRNAs in Anopheles anthropophagus midgut after blood-feeding and Plasmodium infection

  • Wenquan Liu1,
  • Zhenhua Hao1,
  • Liyang Huang1,
  • Lingzi Chen1,
  • Qimei Wei1,
  • Liya Cai1 and
  • Shaohui Liang1Email author
Parasites & Vectors201710:86

https://doi.org/10.1186/s13071-017-2027-6

Received: 12 October 2016

Accepted: 9 February 2017

Published: 16 February 2017

Abstract

Background

Anopheles anthropophagus is one of the major vectors of malaria in Asia. MicroRNAs (miRNAs) play important roles in cell development and differentiation as well as in the cellular response to stress and infection. In a former study, we have investigated the global miRNA profiles in relation to sex in An. anthropophagus. However, the miRNAs contributing to the blood-feeding and infection with Plasmodium are still unknown.

Methods

High-throughput sequencing was performed to identify miRNA profiles of An. anthropophagus midguts after blood-feeding and Plasmodium infection. The expression patterns of miRNA in different midgut libraries were compared based on transcripts per million reads (TPM), and further confirmed by Northern blots. Target prediction and pathway analysis were carried out to investigate the role of regulated miRNAs in blood-feeding and Plasmodium infection.

Results

We identified 67 known and 21 novel miRNAs in all three libraries (sugar-feeding, blood-feeding and Plasmodium infection) in An. anthropophagus midguts. Comparing with the sugar-feeding, the experssion of nine (6 known and 3 novel) and ten (9 known and 1 novel) miRNAs were significantly upregulated and downregulated respectively after blood-feeding (P < 0.05, fold change ≥ 2 and TPM ≥ 10). Plasmodium infection induced the expression of thirteen (9 known and 4 novel) and eleven (9 known and 2 novel) miRNAs significantly upregulated and downregulated, respectively, compared with blood-feeding. The representative upregulated miR-92a in blood-feeding and downregulated miR-275 in Plasmodium infection were further confirmed by Northern Blot. Putative targets of these regulated miRNAs were further investigated and classified into their pathways.

Conclusions

This study suggests that miRNAs are involved in the blood-feeding and Plasmodium infection in An. anthropophagus midgut. Further studies of the function of these differential expressed miRNAs will facilitate in better understanding of mosquito biology and anti-parasite immunity.

Keywords

MicroRNA (miRNA) Anopheles anthropophagus Midgut Blood-feeding Plasmodium infection

Background

Malaria, caused by Plasmodium parasites, is a major threat to human health worldwide. According to the reports from WHO, there were 214 million new infection cases of malaria, and about 438,000 malaria deaths worldwide in 2015 [1]. Plasmodium spp. are transmitted to humans by the blood-feeding of Anopheles mosquitoes. To accomplish the life-cycle, the Plasmodium must overcome numerous attacks from the innate immunity system in mosquitoes [2, 3]. Mosquitoes also have developed various mechanisms to confront Plasmodium infection [3, 4]. The midgut represents the first barrier for the survival and development of Plasmodium [46]. A key concept that has emerged from recent studies is the molecular mechanisms of mosquito midgut that negatively and positively modulate the invasion of the parasite [5, 6].

MicroRNAs (miRNAs) are 18–24 nucleotides non-coding RNAs that regulate gene expression at the post-transcriptional level [7]. They are now considered as a key mechanism of gene regulation in many cellular processes including development, differentiation, apoptosis and innate immunity [810]. miRNAs are also involved in the physiological functions of mosquitoes, such as sexual differences and blood-feeding, even in the control of viral and parasitic infections [1116]. Until now, the miRNA profiles of midguts from Anopheles gambiae and Anopheles stephensi have been reported, and several miRNAs expression levels were shown to be altered during the response to blood-feeding and Plasmodium infection [13, 1518]. For example, blood-feeding and Plasmodium infection in An. stephensi revealed regulation of 13 and 16 miRNAs respectively [18]. Anopheles anthropophagus is a species of mosquito that sucks human blood and transmits malaria (Plasmodium vivax and Plasmodium falciparum) as well as Brugia malayi filariasis in Southeast Asia [12, 19, 20]. The role of miRNAs in An. anthropophagus during blood-feeding and Plasmodium infection are still unknown.

In this study, we employed small RNA sequencing to identify miRNA expression profiles from three samples of An. anthropophagus midguts: sugar-feeding, blood-feeding and Plasmodium infection. The differential expression of miRNAs were further analyzed by target prediction and pathway analysis to reveal their roles in blood-feeding and Plasmodium infection. Our results provide novel regulated miRNAs information of An. anthropophagus during blood-feeding and parasite infection. Understanding the functions of these regulated miRNAs will help investigate mosquito biology and control mosquito-borne infectious diseases.

Methods

Mosquito rearing and Plasmodium infection

Anopheles anthropophagus (China wild type strain) was reared and maintained in humidified incubators at 26 ± 1 °C on a 12 h light:dark cycle as described previously [12]. Mice (18–20 g male ICR mice) were used for mosquito blood-feeding and Plasmodium infection experiments. Stock solution of Plasmodium berghei ANKA strain (0.2 ml) were thawed and intraperitoneally (i.p.) injected into ICR mice by using a 1.0 ml tuberculin syringe. When the gametocytes were confirmed, naive 4–5 day-old female mosquitoes were fed on P. berghei-infected or uninfected ICR mice. Mosquitoes were kept at 26 ± 1 °C until dissection.

Mosquito dissections and total RNA isolation

Dissection of female adults midgut was performed 48 h post-blood-feeding on Plasmodium-infected or uninfected ICR mice. Adult female mosquitoes with 10% sugar-feeding were collected as control sample at 48 h. Dissections were performed on ice in RNAlater® Stabilization Solution (Ambion, Austin, U.S.A) and kept on ice. The midguts were dissected from the abdomen as described before [11]. Total RNAs were extracted from dissected tissue using mirVanaTMmiRNA Isolation Kit (Ambion). Quality and quantity of RNA was checked by using denaturalization agar gels and Du530 Spectrophotometer (Beckman, Krefeld, Gemany).

Small RNA sequencing

The small RNA samples were collected and subsequently sequenced by illumina Hiseq2000 as described before [12]. Briefly, the small RNAs were ligated with RNA adapter followed by reverse transcription using RT primers. Following PCR amplification of the adaptor enriched fragments, the PCR-amplified cDNAs were size-selected using electroelution to obtain the small RNA population with length 119–134 bp. These small RNA libraries were then sequenced using illumina Hiseq2000 (BGI, China). There were 3 biological replicates for each library sequencing.

Computational analysis of small RNA sequencing data

Raw reads generated by high-throughput sequencing was processed as previously described by us [11, 12] with slight modifications. First, low quality reads and reads smaller than 15 nucleotides (nt) were removed from the three small RNA read datasets of midguts, respectively. Clean reads derived from deep sequencing were trimmed and filtered with BOWTIE software to fetch sequences having an appropriate length (15–32 nt). Mature and pre-miRNA sequences of available mosquito species, i.e. Anopheles gambiae and Aedes aegypti, were used as a reference miRNAs database from miRBase v.21. Identification of the novel miRNA was performed by using RNA fold and miRDeep2 as described before [12]. The small RNA datas have been submitted to the NCBI/GEO database with the accession number GSE93545. Transcripts per million reads (TPM) for each miRNA in all three libraries was preformed to generate a comparative analysis of different midgut samples [17, 18]. A heatmap was generated by MeV software.

Northern blot

Northern hybridizations were conducted using digoxigenin-labeled miRCURY LNA probes (Exiqon, Vedbak, Denmark) as described before [12]. Briefly, the midguts of female mosquitoes were collected at 48 h after blood-feeding. The total RNA sample was extracted using mirVana™ miRNA Isolation Kit (Ambion). Total RNA (20 μg) was loaded in 15% denaturing polyacrylamide gels. The RNA gels were transferred to a nylon membrane (Ambion), crosslinked using a UV crosslinker, and prehybridized, then hybridized overnight in the ULTRAhyb-Oligo Hybridization Buffer (Ambion) with the appropriate DIG-labeled probe at 42 °C. After washed three times, the membranes were then incubated for 5 min in development buffer. Substrate (1:100 diluted in development buffer) was applied on to the membranes and incubated in dark for 10 min. Chemiluminescence signal was then measured to detect miRNA on the membrane. Antisense 5′ digoxigenin-labeled miRCURY LNA probes sequences as follow, aan-miR-92a: 5′-TCA GCC GCT GTC ACA CGC ACA G-3′; aan-miR-275: 5′-GAC CAA TCG CCG TCC CCG CCG-3′.

miRNA target prediction and pathway analysis

mRNA targets of regulated miRNAs during blood-feeding and Plasmodium infection were predicted according to the protocol reported before [18, 19]. Briefly, 3′UTR sequence of An. gambiae and Ae. aegypti were downloaded from VectorBase and subjected to target prediction using RNA hybrid tool. Target predictions were carried out based on the following three parameters: (i) the perfect complementarity of the miRNA with the 3′UTR sequence of mRNA targets; (ii) the energy of the miRNA:mRNA target duplex ≤ -20 Kcal/mol; (iii) P-value < 0.05. The selected miRNA:mRNA interaction networks were generated and visualized by Cytoscape.

Statistical analysis

Statistical tests for identifying significant differentially expressed miRNAs were performed using t-test. The P-value cut-off was carried out on the data with the significance selected as 0.05.

Results

Small RNA sequencing analysis

The small RNA libraries were enriched, and separately yielded 17.51 million, 15.87 million and 6.5 million raw reads from the midgut samples of sugar-feeding (SF), blood-feeding (BF) and Plasmodium infection (PI). After filtering for linker sequences, and removing ambiguous reads, high quality clean reads with sizes ranged from 15 to 32 nt were collected, of which 14.45 million (82.56%) reads for SF, 11.75 million (74.02%) reads for BF, and about 3.2 million (50.38%) reads for PI (Table 1). The main populations of small RNA in the length distribution is 20–23 nt (Fig. 1), of which 10.1 million (68.89%) reads for SF, 5.76 million for BF (49.01%) and 1.45 million for PI (45.51%). The 20–23 nt peaks of small RNA are consisted with the expected size of microRNAs. After aligned to the known miRNA and pre-miRNA in miRBase (version 21.0), the miRNA reads from SF, BF and PI midguts which can match to the miRBase are 4.03 million, 1.88 million and 0.21million separately (Table 1).
Table 1

Composition of small RNAs in midguts of sugar-feeding (SF), blood-feeding (BF) and Plasmodium infected blood-feeding (PI) at 48 h

Sample name

Raw reads

Clean reads (%)

20–23 nt reads (%)

Reads matched to known miRNAs

SF

17,510,882

14,457,010 (82.56)

10,104,472 (68.89)

4,032,195

BF

15,875,249

11,750,726 (74.02)

5,759,187 (49.01)

1,878,567

PI

6,511,125

3,180,161 (50.38)

1,447,177 (45.51)

209,907

Total

39,897,256

29,387,897

17,310,836

6,120,669

Fig. 1

Length distribution of small RNA reads in the An. anthropophagus midguts. Female mosquitoes midguts at 48 h sugar feeding (SF), blood-feeding (BF) and Plasmodium infected-blood-feeding (PI)

Identification of miRNAs in midguts of An. anthropophagus

A total of 67 known and 21 novel miRNAs were identified in the SF, BF and PI midgut libraries (Tables 2 and 3). Both known and novel miRNAs were present in the miRNA databases of An. anthropophagus adult as described before [12]. The distribution and relative abundance of these miRNAs were analyzed according to the set of “abundant” (TPM ≥ 1000) or “rare” (TPM < 10). Most of known and novel miRNAs were found in the SF library (n = 87) with 29 abundant and 20 rare miRNAs. In case of BF library, 26 miRNAs were found to be abundant and 7 miRNAs (miR-193, miR-2944b, miR-307, miR-309, miR-79-5p, miR-965-1 and miR-N6) were not detected. Meanwhile, 27 abundant and 20 rare miRNAs were found in the PI library, of which 7 miRNAs (miR-193, miR-210, miR-252, miR-263b-5p, miR-2944b, miR-307 and miR-87) have no read counts. A novel miRNA miR-N6 is unique in the PI midguts.
Table 2

List of known microRNAs expressed in An. anthropophagus midguts

    

Transcripts per million reads(TPM)

Fold change

Reference miRNA database

No.

Name

Length

Sequence

SF

BF

PI

BF vs SF

PI vs BF

aga

aae

1

aan-bantam

23

UGAGAUCACUUUGAAAGCUGAUU

45,048.5

24,928

74,181

0.553

2.975

Y

Y

2

aan-let-7

21

UGAGGUAGUUGGUUGUAUAGU

1,079.9

646.7

1,199.2

0.598

1.854

Y

Y

3

aan-miR-1*#

22

UGGAAUGUAAAGAAGUAUGGAG

78

20.6

75.5

0.264

3.665

Y

Y

4

aan-miR-10

22

ACCCUGUAGAUCCGAAUUUGUU

60.7

51.9

136

0.855

2.62

Y

Y

5

aan-miR-100#

22

AACCCGUAGAUCCGAACUUGUG

6,270.2

4632

17,444

0.738

3.766

Y

Y

6

aan-miR-1000

21

AUAUUGUCCUGUCACAGCAGU

56.6

129.1

115

2.280

0.89

Y

Y

7

aan-miR-11

22

CAUCACAGUCUGAGUUCUUGCU

3,794.7

1775.7

2,952.7

0.467

1.662

Y

Y

8

aan-miR-1174

21

UCAGAUCUACUUAAUACCCAU

4,040

1,674

3,367.6

0.414

2.011

Y

Y

9

aan-miR-1175-3p*

24

UGAGAUUCUACUUCUCCGACUUAA

7,421.2

2,027

4,887.2

0.273

2.411

Y

Y

10

aan-miR-12*#

23

UGAGUAUUACAUCAGGUACUGGU

5,206.8

1,768.2

13,763.1

0.339

7.783

Y

Y

11

aan-miR-125

22

UCCCUGAGACCCUAACUUGUGA

308.6

257.7

791.3

0.835

3.07

Y

Y

12

aan-miR-13*#

23

UAUCACAGCCAUUUUGACGAGUU

790.5

241.7

1,988.8

0.305

8.228

 

Y

13

aan-miR-133#

22

UUGGUCCCCUUCAACCAGCUGU

96.3

87.3

17.4

0.906

0.199

Y

Y

14

aan-miR-137

22

UAUUGCUUGAGAAUACACGUAG

20.1

14.8

34.9

0.736

2.358

 

Y

15

aan-miR-14#

22

UCAGUCUUUUUCUCUCUCCUAU

8,439

3,009.3

1,009.7

0.356

0.335

Y

Y

16

aan-miR-184

22

UGGACGGAGAACUGAUAAGGGC

6,949.9

4,432

9,130

0.637

2.06

Y

Y

17

aan-miR-1890

22

UGAAAUCUUUGAUUAGGUCUGG

112.6

41.8

80.2

0.371

1.918

Y

Y

18

aan-miR-1891

22

UGAGGAGUUAAUUUGCGUGUUU

3.8

3.2

3.5

0.842

1.093

Y

Y

19

aan-miR-190

24

AGAUAUGUUUGAUAUUCUUGGUUG

47

37.6

34.9

0.8

0.928

Y

Y

20

aan-miR-193

22

UACUGGCCUACUAAGUCCCAAC

0.2

0

0

0

Y

Y

21

aan-miR-210#

19

CUUGUGCGUGUGACAACGG

5.1

11.1

0

2.176

Y

Y

22

aan-miR-252*#

21

UAAGUACUAGUGCCGCAGGAG

59

192

0

3.254

 

Y

23

aan-miR-263a#

23

AAUGGCACUGGAAGAAUUCACGG

36.3

50.8

262.7

1.399

5.17

Y

Y

24

aan-miR-263b-5p

23

CUUGGCACUGGGAGAAUUCACAG

0.2

1.1

0

5.5

 

Y

25

aan-miR-275#

20

UCAGGUACCUGAAGUAGCGC

3,269.4

9,040.2

2,035.9

2.765

0.225

Y

Y

26

aan-miR-276*#

20

AGCGAGGUAUAGAGUUCCUA

23.9

7.4

27.9

0.309

3.770

Y

Y

27

aan-miR-277-5p

22

UAAAUGCACUAUCUGGUACGAC

15.3

7.7

14.9

0.503

1.935

Y

Y

28

aan-miR-278

23

ACGGACGAUAGUCUUCAGCGGCC

55.7

22.8

31.4

0.409

1.377

Y

Y

29

aan-miR-279*

22

UGACUAGAUCCACACUCAUUAA

2,787.5

926.2

1,746.5

0.332

1.885

Y

Y

30

aan-miR-283*

24

CAAUAUCAGCUGGUAAUUCUGGGC

1,231.3

407.5

693.7

0.330

1.702

Y

Y

31

aan-miR-285

22

UAGCACCAUUCGAAAUCAGUAC

10.4

23

18.6

2.211

0.808

 

Y

32

aan-miR-2944a

24

GAAGGAACUUCUGCUGUGAUCUGA

9.9

7.1

14.9

0.717

2.098

 

Y

33

aan-miR-2944b

23

GAAGGAACUCCCGGUGUGAUAUA

1.2

0

0

0

 

Y

34

aan-miR-2a#

23

UAUCACAGCCAGCUUUGAAGAGC

3,687.6

1,745.2

415.5

0.473

0.238

 

Y

35

aan-miR-2b*

24

UAUCACAGCCAGCUUUGAUGAGCU

791.7

269.1

618.2

0.339

2.297

 

Y

36

aan-miR-2c

23

UAUCACAGCCAGCUUUGAUGAGC

657.4

228

435.7

0.346

1.910

 

Y

37

aan-miR-305-5p

24

AUUGUACUUCAUCAGGUGCUCUGG

2,879.6

1,101.7

3,594

0.382

3.262

Y

Y

38

aan-miR-306

22

UCAGGUACUGGAUGACUCUCAG

14,095.1

11,671.1

11,518.1

0.828

0.986

Y

 

39

aan-miR-307

20

UCACAACCUCCUUGAGUGAG

1.3

0

0

0

Y

Y

40

aan-miR-308-5p#

22

AAUCACAGGAGUAUACUGUGAG

85.4

39.7

338.2

0.464

8.518

Y

Y

41

aan-miR-309

22

UCACUGGGCAAAGUUUGUCGCA

0.2

0

3.5

0

Y

 

42

aan-miR-31#

23

UGGCAAGAUGUUGGCAUAGCUGA

27.3

10.7

207

0.392

19.35

 

Y

43

aan-miR-315-5p*

22

UUUUGAUUGUUGCUCAGAAAGC

2.9

12.2

7

4.206

0.573

Y

Y

44

aan-miR-317

24

UGAACACAUCUGGUGGUAUCUCAG

6,342

6,462.2

2,795.8

1.018

0.432

Y

 

45

aan-miR-34

20

UGGCAGUGUGGUUAGCUGGU

584.9

347.7

990.1

0.594

2.847

Y

Y

46

aan-miR-375-1

22

UUUGUUCGUUUGGCUCGAGUUA

10

6.9

20.9

0.69

3.028

Y

Y

47

aan-miR-71-3p

22

UCUCACUACCUUGUCUUUCAUG

279.6

174.3

100.6

0.623

0.577

 

Y

48

aan-miR-79-5p

22

UAAAGCUAGAUUACCAAAGCAU

1.4

0

3.7

0

Y

Y

49

aan-miR-8-5p

23

UAAUACUGUCAGGUAAAGAUGUC

1,990.3

3,340.8

2,245.8

1.678

0.672

Y

Y

50

aan-miR-87

21

GGUGAGCAAAUAUUCAGGUGU

3.3

2.6

0

0.787

0

Y

 

51

aan-miR-927

22

UUUAGAAUUCCUACGCUUUACC

142

156.1

156.9

1.099

1.005

Y

 

52

aan-miR-92a*

22

UAUUGCACUUGUCCCGGCCUAU

2,103.2

17,462.4

8,613.8

8.3

0.493

Y

 

53

aan-miR-92b*#

22

AAUUGCACUUGUCCCGGCCUGC

503.4

2,215.5

366

4.40

0.165

Y

Y

54

aan-miR-932-3p*

23

UCAAUUCCGUAGUGCAUUGCAGU

1.7

10.4

7.4

6.117

0.711

 

Y

55

aan-miR-957#

22

UGAAACCGUCCAAAACUGAGGC

98.6

260.4

31.4

2.640

0.120

Y

Y

56

aan-miR-965-1

22

UAAGCGUAUAGCUUUUCCCAUU

2

0

7.4

0

Y

Y

57

aan-miR-970

21

UCAUAAGACACACGCGGCUAU

4,435.5

2,880.7

1,837.2

0.649

0.637

Y

Y

58

aan-miR-980#

24

CGGCCGUUCAUUGGGUCAUCUAGC

25.5

26.9

4.9

1.054

0.182

 

Y

59

aan-miR-981

22

UUCGUUGUCGACGAAACCUGCA

3.5

7.7

7.4

2.2

0.961

Y

Y

60

aan-miR-988-5p

22

CCCCUUGUUGCAAACCUCACGC

4

2.2

7.4

0.55

3.363

Y

Y

61

aan-miR-989*

21

UGUGAUGUGACGUAGUGGUAC

78

360

314

4.615

0.872

Y

Y

62

aan-miR-996*

21

UGACUAGAUUACAUGCUCGUC

6,082.2

1,855

5,933.3

0.305

3.198

Y

Y

63

aan-miR-998

19

UAGCACCAUGAGAUUCAGC

3.5

0.5

3.7

0.142

7.4

 

Y

64

aan-miR-9a

23

UCUUUGGUUAUCUAGCUGUAUGA

1,308.7

650.5

2,178.8

0.497

3.349

Y

Y

65

aan-miR-9b

24

UCUUUGGUGAUUUUAGCUGUAUGC

2,386.6

1,190.8

3,806.8

0.499

3.196

Y

Y

66

aan-miR-9c

22

UCUUUGGUAUUCUAGCUGUAGA

8,976.6

4,496.3

12,856.5

0.500

2.859

Y

Y

67

aan-miR-iab-4#

22

ACGUAUACUGAAUGUAUCCUGA

14.6

3.3

26.1

0.226

7.909

Y

Y

Abbreviations: SF sugar-feeding, BF blood-feeding, PI Plasmodium infection, TPM transcripts per million reads; Fold change, BF/SF or PI/BF, aan Anopheles anthropophagus, aga Anopheles gambiae, aae Aedes aegypti

*Represents the statistical significance in the differential expression of individual miRNA in BF vs SF (P-value < 0.05)

#Represents the statistical significance in the differential expression of individual miRNA in PI vs BF (P-value < 0.05)

Table 3

List of novel microRNAs expressed in An. anthropophagus midguts

    

Transcripts per million reads(TPM)

Fold change

No.

Name

Length

Sequence

SF

BF

PI

BF vs SF

PI vs BF

1

aan-miR-N1#

18

CGCUGCAGUACUGGCGCC

37.2

17.5

306.8

0.470

17.53

2

aan-miR-N2

21

AUCCGGUGAUAGGCUGACCCG

50.6

26.5

45.3

0.523

1.709

3

aan-miR-N3*

20

UUAGAAUGUGGAAUCUGUUU

27.9

7.4

24.4

0.265

3.297

4

aan-miR-N4

23

UUGGUGUUAUAUCUUACAGUGAG

3,160.5

2,462.1

3,824.2

0.779

1.553

5

aan-miR-N5

23

UUGGUGUUAUAUCUUACAGUGAG

9.8

4.2

10.5

0.428

2.5

6

aan-miR-N6#

20

UAUCACAGCCAGCUUUGAAG

0

0

149.9

7

aan-miR-N7

22

UGCAUUCAGUGGGGCGGUCGUG

4.2

4.8

3.5

1.142

0.729

8

aan-miR-N8

22

UGUUAACUGUAAGACUGUGUCG

5.1

1.1

7

0.215

6.363

9

aan-miR-N9

21

UAGCACCAUGAGAUUCAGCUC

794,875

861,632.1

704,445.5

1.083

0.817

10

aan-miR-N10*

22

UCAAUUCCGUAGUGCAUUGCAGU

80.5

273.1

491.5

3.392

1.799

11

aan-miR-N11

22

UUGGUGUUAUAUCUUACAGUGAG

1,479.9

592.8

348.6

0.400

0.588

12

aan-miR-N12

22

GUAGGCCGGCGGAAACUACUUGC

351.2

132.3

142.9

0.376

1.080

13

aan-miR-N13

23

UUGGCCGGUACGGGCUGACCGGGC

49.5

29.6

55.8

0.597

1.885

14

aan-miR-N14#

22

UGAACCGGCGUAGCGUGAAAGCA

8,461.5

3,769.9

18,898.1

0.445

5.012

15

aan-miR-N15

21

CUAAGUACUAGUGCCGCAGGAG

3,368.1

2,159.9

6,177.4

0.641

2.860

16

aan-miR-N16*#

19

UUAGAAUGUGGAAUCUGUUU

755.7

4,010.7

324.2

5.307

0.080

17

aan-miR-N17

22

UGAAAUCUUUGAUUAGGUCUGG

59.3

60.3

125.5

1.016

2.081

18

aan-miR-N18#

18

UAUCAGCGGUAGUUACCUG

22.7

13.8

3.5

0.607

0.253

19

aan-miR-N19

20

GUGCAUUGUAGUUGCAUUGCA

3,771.7

1,716.9

3,099.1

0.455

1.805

20

aan-miR-N20#

18

GUUGCUGUCCGCUGAAGCA

59.8

25.4

209.2

0.424

8.236

21

aan-miR-N21*

23

UGGCAAGAUGUUGGCAUAGCAGCU

14.2

97.6

134.1

6.873

1.373

Abbreviations: N the abbreviation of Novel, SF sugar feeding, BF blood feeding, PI Plasmodium infection, TPM transcripts per million reads; Fold change, BF/SF or PI/ BF, aan Anopheles anthropophagus

*Represents the statistical significance in the differential expression of individual miRNA in BF vs SF (P-value < 0.05)

#Represents the statistical significance in the differential expression of individual miRNA in PI vs BF (P-value < 0.05)

Differential expression of miRNAs after blood-feeding and Plasmodium infection

To investigate the miRNAs regulated in blood-feeding and Plasmodium infection, the TPM value of individual miRNA between BF and SF, PI and BF were compared in the form of the fold change (Fig. 2). The set of significantly regulated miRNAs were selected on the basis of their P-value < 0.05, fold change ≥ 2 and TPM ≥10 (Additional file 1: Table S1 and Additional file 2: Table S2). With respect to BF library, six known miRNAs (miR-252, miR-315-5p, miR-92a, miR-92b, miR-932-3p and miR-989) and three novel miRNAs (miR-N10, miR-N16 and miR-N21) showed significantly upregulated expression compared with SF. Especially, the expression levels for miR-92a, miR-932-3p, miR-N16 and miR-N21 were enhanced more than 5-fold after blood-feeding. Meanwhile, the expression of 10 miRNAs (miR-1, miR-1175-3p, miR-12, miR-13, miR-276, miR-279, miR-283, miR-2b, miR-996 and miR-N3) was found to be significantly downregulated after blood-feeding (Tables 2 and 3).
Fig. 2

Heatmap of differential expressed miRNAs in the An. anthropophagus midgut. Expression profile of known miRNAs (a) and novel miRNAs (b) in sugar-feeding midgut (SF), blood-feeding midgut (BF) and Plasmodium-infected midguts (PI). Colour scale-bar from light green to dark red indicates relative increase in miRNA expression

In the case of the PI library, 9 known miRNAs (miR-1, miR100, miR-12, miR-13, miR-263a, miR-276, miR-308-5p, miR-31, miR-iab-4) and 4 novel miRNAs (miR-N1, miR-N6, miR-N14 and miR-N20) exhibited significantly upregulated expression. Notably, the TPM value of aan-miR-N16 was almost 150-fold enhanced in PI midgut compared with BF midgut. Meanwhile, the expression of miR-31 was enhanced more than 10-fold after Plasmodium infection. On the other hand, the expression levels of 9 known (miR-133, miR-14, miR-210, miR-252, miR-2a, miR-275, miR-92b, miR-957 and miR-980) and 2 novel miRNAs (miR-N16 and miR-N18) were significantly downregulated in the Plasmodium infection midgut. In particular, the TPM value for miR-252 was 191 in BF midgut but zero in the PI midgut (Tables 2 and 3).

Furthermore, several miRNA exhibited different expression patterns in BF and PI midguts. For example, four miRNAs including miR-1, miR-12, miR-13 and miR-276 were significantly downregulated in the BF midguts but upregulated in the PI midguts. Conversely, the expression of miR-252 and miR-N16 showed significant upregulation in the BF midguts but were downregulated in the PI midguts (Tables 2 and 3). These data suggest that the same miRNAs may play different roles in the blood-feeding and Plasmodium infection midguts.

Confirmation of mosquito miRNAs

Having identified the microRNAs expression patterns of An. anthropophagus midguts after blood-feeding and Plasmodium infection, Northern hybridization was performed to validate some of these miRNAs. The representative abundant miRNA (miR-275 and miR-92a) were selected in the Northern blot analysis. Total RNA from SF, BF and PI 48 h were probed using locked nucleic acid (LNA) probe. The expression patterns of miR-275 and miR-92a in SF, BF and PI are shown in Fig. 3. The northern signals at ~20 nt indicated that the expression level of miR-92a is upregulated to over 5-fold after blood-feeding, and miR-275 is downregulated to more than 3-fold after Plasmodium infection, which is consistent with the sequencing results.
Fig. 3

Northern blot of differentially expressed miRNAs in the An. anthropophagus midguts. Abbreviations: SF, sugar-feeding midgut; BF, blood-feeding midgut; PI, Plasmodium-infected midgut

Target prediction and network analysis

Target prediction were carried out to understand the putative function of regulated miRNAs by using RNA hybrid (P-value < 0.05). All the mRNA targets are derived from the orthologous genes of An. gambiae in the Vector Base. Maximum numbers of mRNA targets were predicted for miR-92a (n = 744). The targets were further analyzed by KOBAS, and a total of 38 different pathways were regulated by miRNAs after blood-feeding and Plasmodium infection (Fig. 4). In the upregulated miRNA after blood-feeding, miR-252 and miR-92a were identified to target oxidative phosphorylation and peroxisome pathway. miR-92a was also found to target glycolysis, proteasome, ribosome and TGF-beta signaling pathway. Meanwhile, ribosome was commonly targeted by the downregulated miRNAs including miR-13, miR-279 and miR-2b in BF midguts. After Plasmodium infection, endocytosis and fructose mannose metabolism pathway is specifically targeted by miR-31. Spliceosome was common between miR100 and miR-14. In the downregulated miRNA, miR-210 and miR-275 targeted RNA transport and purine metabolism separately (Fig. 4).
Fig. 4

KOBAS analysis of miRNA targets. miRNA targets involved in blood-feeding and Plasmodium infection were predicted by RNA hybrid. Pie chart represents the pathway targeted by the miRNAs which are listed in the pie area

The interaction networks of miRNA:mRNA targets were generated. In BF midguts, AGAP000348-RA and AGAP000120-RA were targeted by downregulated miR-1. A total of 493 mRNA targets were regulated by upregulated miR-989, of which AGAP007839-RA, AGAP0004451-RA and AGAP008345-RA were targeted by two or more upregulated and downregulated miRNAs (Additional file 3: Figure S1). Furthermore, 573 mRNAs were targeted by upregulated miR-31 in the PI midguts. Three mRNA (AGAP000204-RA, AGAP009077-RA and AGAP010306-RA) were targeted by downregulated miR-957 (Additional file 4: Figure S2). Further investigations are important to understand the role of these mRNAs during blood-feeding and Plasmodium infection.

Discussion

Anopheles mosquitoes are the main vector for the transmission of malaria, and the invasion of midguts is one of the most critical steps for the survival and development of Plasmodium [21]. Because of the lack of an adaptive immune system, RNA interference (RNAi) is the most important and primary defense employed by mosquitoes to protect themselves from pathogens [22, 23]. miRNAs are the most important part of RNA interference, and some of them have been proven to be involved in fighting against virus and parasite infections [2124]. In a previous study, we have investigated the miRNA profiles of An. anthropophagus, and discovered that several miRNAs exhibit sexual differences and stage-specific functions [12]. However, their role in Plasmodium infection is poor understood. This study was conducted to identify and elucidate role of miRNAs after blood-feeding and Plasmodium infection in An. anthropophagus midguts.

In the present study, the expression profiles of miRNA at 48 h after blood-feeding and Plasmodium infection were investigated using small RNA sequencing; sugar-feeding midguts were taken as the control. Compared with sugar feeding, there is a significant downregulation for the percent of miRNA library (20–23 nt) after blood-feeding and Plasmodium infection. The similar downregulation expression profiling of miRNAs was also reported in the parasitized blood-feeding midguts of Anopheles stephensi [18] suggested that pathogens including Plasmodium and viruses invade mosquitoes midguts by inhibiting or downregulating the miRNA expression [2224].

Blood-feeding is a critical physiological activity for the mosquito and its ability to transmit disease [25]. Female mosquitoes take blood meals to carry out egg production, and acquire pathogens such as malaria parasites and dengue viruses from an infected host [2628]. Furthermore, mosquito biological characteristics are affected by triggering a cascade of gene regulatory events in the midgut after blood-feeding [28, 29]. In this study, the correlation between blood-feeding and miRNA expression profiles in An. anthropophagus midguts was investigated. We indentified nine significantly upregulated and ten downregulated expression miRNAs in the blood-feeding midguts compared with sugar feeding. For example, miR-92a exhibited a significant enhanced expression level after blood-feeding according to the sequencing and Northern blot results. Li et al. [11] reported that blood-feeding can induce the expression of miR-92a in the midgut of Ae. aegypti females. While Wolbachia infection can downregulate the expression of miR-92a in mosquito cell [30]. In our study, the expression of miR-92a is significantly upregulated in blood-feeding. By target prediction and networks analysis, miR92a was found to target 744 mRNAs and several pathways including the oxidative phosphorylation, proteasome, ribosome and TGF-beta singaling pathway. Our study may shed light on the possible roles of miR92a in An. anthropophagus mosquito physiology related to blood-feeding.

Midgut represents the first barrier for the pathogens to establish infection in mosquitoes [2, 4]. Pathogens such as parasites and endosymbiotic bacteria are known to alter host miRNA profiles for their invasion and development [15, 31]. For example, Wolbachia can induce the expression of miR-12 in An. aegypti mosquito cells to maintan the persistence of infection [31]. In our study, we indentified 13 upregulated miRNAs including miR-12. This finding suggests that miR-12 is also involved in Plasmodium infection, and the role of miR-12 in Plasmodium infection of An. anthropophagus need to be further elucidated. Meanwhile, we also found downregulated expression for 11 miRNA in Plasmodium infection midguts of An. anthropophagus, six of which including miR-14, miR-2a, miR-92b, miR-957, miR-980 and miR-275, were also shown to be downregulated in Plasmodium-infected Anopheles stephensi [18]. miR-275 is indispensable for blood digestion and egg development in the mosquito Ae. aegypti [32]. In this study, we confirmed that miR-275 is also involved in Plasmodium infection of An. anthropophagus by sequencing and further Northern blot analysis.

Conclusions

In conclusion, our study provides the significant experimental data on the expression profile of microRNAs in An. anthropophagus midgut after blood-feeding and Plasmodium infection. Differentially expressed miRNA in SF, BF and PI were identified by small RNA sequencing, and further validated by Northern blot. By comparative analyzing differentially expression levels of the miRNA in sugar-feeding, blood-feeding and Plasmodium infection, we found several significant miRNAs involved in the interaction of mosquito host and parasite Plasmodium. Our study provides novel regulated miRNAs information of An. anthropophagus during blood-feeding and parasite infection. Elucidation of regulated miRNA functions will provide strong foundation for better understanding of the biology of the mosquitoes and mosquito-parasite interactions.

Abbreviations

BF: 

Blood-feeding

FPKM: 

Fragments Per Kilobase of transcript per Million mapped reads

miRNAs: 

microRNAs

PI: 

Plasmodium infection

RNAi: 

RNA interference

SF: 

Sugar feeding

Declarations

Acknowledgements

This work was supported by the National Natural Science Foundation of China (81071389), the Commonweal Technology Application Project of Zhejiang Province (2015C33101) and the Technology Project of Wenzhou City (Y20140481).

Funding

The National Natural Science Foundation of China (81071389), the Commonweal Technology Application Project of Zhejiang Province (2015C33101) and the Technology Project of Wenzhou City (Y20140481).

Availability of data and materials

The datasets supporting the conclusions of this article are available in the Gene Expression Omnibus repository, [GSE93545, https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE93545].

Authors’ contributions

SHL and WQL conceived and designed the study. WQL, ZHH, LYH, LZC, QMW and LYC performed the experiments, analyzed the data and drafted the manuscript. SHL critically revised the manuscript. All authors read and approved the final manuscript.

Competing interests

The authors declare that they have no competing interests.

Consent for publication

Not applicable.

Ethics approval

Animal experiments were performed in accordance with National animal ethics guidelines of the Government of P. R. China after approval by the Animal Ethics Committees of Wenzhou Medical University (Permit number: SYXKZHE20150009).

Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.

Authors’ Affiliations

(1)
Department of Parasitology, Wenzhou Medical University

References

  1. World Malaria Report 2015. http://www.who.int/malaria/publications/world-malaria-report-2015/en/. Accessed on 3 Jan. 2016
  2. Sinden RE, Alavi Y, Raine JD. Mosquito-malaria interactions: a reappraisal of the concepts of susceptibility and refractoriness. Insect Biochem Mol Biol. 2004;34:625–9.View ArticlePubMedGoogle Scholar
  3. Meister S, Koutsos AC, Christophides GK. The Plasmodium parasite - a ‘new’ challenge for insect innate immunity. Int J Parasitol. 2004;34:1473–82.View ArticlePubMedGoogle Scholar
  4. Whitten MM, Shiao SH, Levashina EA. Mosquito midguts and malaria: cell biology, compartmentalization and immunology. Parasite Immunol. 2006;28:121–30.View ArticlePubMedGoogle Scholar
  5. Zhang G, Niu G, Franca CM, Dong Y, Wang X, Butler NS, et al. Anopheles Midgut FREP1 mediates Plasmodium invasion. J Biol Chem. 2015;290:16490–501.View ArticlePubMedPubMed CentralGoogle Scholar
  6. Nikolaeva D, Draper SJ, Biswas S. Toward the development of effective transmission-blocking vaccines for malaria. Expert Rev Vaccines. 2015;14:653–80.View ArticlePubMedGoogle Scholar
  7. Bartel DP. MicroRNAs: genomics, biogenesis, mechanism, and function. Cell. 2004;116:281–97.View ArticlePubMedGoogle Scholar
  8. Miska EA. How microRNAs control cell division, differentiation and death. Curr Opin Genet Dev. 2005;15:563–8.View ArticlePubMedGoogle Scholar
  9. Zhang R, Wang YQ, Su B. Molecular evolution of a primate-specific microRNA family. Mol Biol Evol. 2008;25:1493–502.View ArticlePubMedGoogle Scholar
  10. Bushati N. MicroRNA Functions. Ann Rev Cell Dev Biol. 2007;23:175–205.View ArticleGoogle Scholar
  11. Li S, Mead AE, Liang SH, Tu ZJ. Direct sequencing and expression analysis of a large number of miRNAs in Aedes aegypti and a multi-species survey of novel mosquito miRNAs. BMC Genomics. 2009;10:581.View ArticlePubMedPubMed CentralGoogle Scholar
  12. Liu W, Huang H, Xing C, Li C, Tan F, Liang S. Identification and characterization of the expression profile of microRNAs in Anopheles anthropophagus. Parasit Vectors. 2014;7:159.View ArticlePubMedPubMed CentralGoogle Scholar
  13. Winter F, Edaye S, Hüttenhofer A, Brunel C. Anopheles gambiae miRNAs as actors of defence reaction against Plasmodium invasion. Nucleic Acids Res. 2007;35(20):6953–62.View ArticlePubMedPubMed CentralGoogle Scholar
  14. Skalsky LR, Vanlandingham LD, Scholle F, Higgs S, Cullen RB. Identification of microRNAs expressed in two mosquito vectors, Aedes albopictus and Culex quinquefasciatus. BMC Genomics. 2010;11:119.View ArticlePubMedPubMed CentralGoogle Scholar
  15. Biryukova I, Ye T, Levashina E. Transcriptome-wide analysis of microRNA expression in the malaria mosquito Anopheles gambiae. BMC Genomics. 2014;15:557.View ArticlePubMedPubMed CentralGoogle Scholar
  16. Castellano L, Rizzi E, Krell J, Di Cristina M, Galizi R, Mori A, et al. The germline of the malaria mosquito produces abundant miRNAs, endo-siRNAs, piRNAs and 29-nt small RNAs. BMC Genomics. 2015;16:100.View ArticlePubMedPubMed CentralGoogle Scholar
  17. Jain S, Rana V, Tridibes A, Sunil S, Bhatnagar RK. Dynamic expression of miRNAs across immature and adult stages of the malaria mosquito Anopheles stephensi. Parasit Vectors. 2015;8:179.View ArticlePubMedPubMed CentralGoogle Scholar
  18. Jain S, Rana V, Shrinet J, Sharma A, Tridibes A, Sunil S, Bhatnagar RK. Blood feeding and Plasmodium infection alters the miRNome of Anopheles stephensi. PLoS One. 2014;9:e98402.View ArticlePubMedPubMed CentralGoogle Scholar
  19. Shrinet J, Jain S, Jain J, Bhatnagar RK, Sunil S. Next generation sequencing reveals regulation of distinct Aedes microRNAs during chikungunya virus development. PLoS Negl Trop Dis. 2014;8(1):e2616.View ArticlePubMedPubMed CentralGoogle Scholar
  20. Geng YJ, Gao ST, Huang DN, Zhao YR, Liu JP, Li XH, Zhang RL. Differentially expressed genes between female and male adult Anopheles anthropophagus. Parasitol Res. 2009;105:843–51.Google Scholar
  21. Sinka ME, Bangs MJ, Manguin S, Coetzee M, Mbogo CM, Hemingway J, et al. The dominant Anopheles vectors of human malaria in Africa, Europe and the Middle East: occurrence data, distribution maps and bionomic précis. Parasit Vectors. 2010;3:117.View ArticlePubMedPubMed CentralGoogle Scholar
  22. Carol DB. Mosquito RNAi is the major innate immune pathway controlling arbovirus infection and transmission. Future Microbiol. 2011;6:265–77.View ArticleGoogle Scholar
  23. Cirimotich CM, Dong Y, Garver LS, Sim S, Dimopoulos G. Mosquito immune defenses against Plasmodium infection. Dev Comp Immunol. 2010;34:387–95.View ArticlePubMedGoogle Scholar
  24. Campbell LC, Harrison T, Hess MA, Ebel DG. MicroRNA levels are modulated in Aedes aegypti following exposure to Dengue-2. Insect Mol Biol. 2014;23(1):132–9.View ArticlePubMedGoogle Scholar
  25. Dennison NJ, BenMarzouk-Hidalgo OJ, Dimopoulos G. MicroRNA-regulation of Anopheles gambiae immunity to Plasmodium falciparum infection and midgut microbiota. Dev Comp Immunol. 2015;49:170–8.View ArticlePubMedGoogle Scholar
  26. Cator LJ, Pietri JE, Murdock CC, Ohm JR, Lewis EE, Read AF, et al. Immune response and insulin signalling alter mosquito feeding behaviour to enhance malaria transmission potential. Sci Rep. 2015;5:11947.View ArticlePubMedPubMed CentralGoogle Scholar
  27. Harrington LC, Edman JD, Scott TW. Why do female Aedes aegypti (Diptera: Culicidae) feed preferentially and frequently on human blood? J Med Entomol. 2001;38:411–22.View ArticlePubMedGoogle Scholar
  28. Ferguson HM, Read AF. Mosquito appetite for blood is stimulated by Plasmodium chabaudi infections in themselves and their vertebrate hosts. Malar J. 2004;3:12.View ArticlePubMedPubMed CentralGoogle Scholar
  29. Putnam JL, Scott TW. Blood-feeding behavior of dengue-2 virus-infected Aedes aegypti. Am J Trop Med Hyg. 1995;52:225–7.PubMedGoogle Scholar
  30. Mayoral JG, Etebari K, Hussain M, Khromykh AA, Asgari S. Wolbachia infection modifies the profile, shuttling and structure of microRNAs in a mosquito cell line. PLoS One. 2014;9:e96107.View ArticlePubMedPubMed CentralGoogle Scholar
  31. Osei-Amo S, Hussain M, O’Neill SL, Asgari S. Wolbachia-induced aae-miR-12 miRNA negatively regulates the expression of MCT1 and MCM6 genes in Wolbachia-infected mosquito cell line. PLoS One. 2012;7(11):e50049.View ArticlePubMedPubMed CentralGoogle Scholar
  32. Bryant B, Macdonald W, Raikhel AS. microRNA miR-275 is indispensable for blood digestion and egg development in the mosquito Aedes aegypti.Proc Natl Acad Sci USA. 2010;107(52):22391–8.Google Scholar

Copyright

© The Author(s). 2017

Comments

By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.

Advertisement