Identification and characterization of the expression profile of microRNAs in Anopheles anthropophagus
© Liu et al.; licensee BioMed Central Ltd. 2014
Received: 26 December 2013
Accepted: 25 March 2014
Published: 1 April 2014
Anopheles anthropophagus, one of the most important mosquito-borne disease vectors in Asia, mainly takes blood meals from humans and transmits both malaria and filariae. MicroRNAs (miRNAs) are small non-coding RNAs, and play a critical role in many cellular processes, including development, differentiation, apoptosis and innate immunity.
We investigated the global miRNA expression profile of male and female adults of A. anthropophagus using illumina Hiseq2000 sequencing combined with Northern blot.
By using the miRNAs of the closely-related species Anopheles gambiae and Aedes aegypti as reference, we obtained 102 miRNAs candidates out of 12.43 million raw sequencing reads for male and 16.51 million reads for female, with 81 of them found as known miRNAs in An. gambiae and/or Ae. aegypti, and the remaining 21 miRNAs were considered as novel. By analyzing the revised read count of miRNAs in male and female, 29 known miRNAs show sexual difference expression: >2-fold in the read count of the same miRNAs in male and female. Especially for miR-989, which is highly expressed in the female mosquitoes, but shows almost no detected expression in male mosquitoes, indicating that miR-989 may be involved in the physiological activity of female mosquito adults. The expression of four miRNAs in different growth stages of mosquito were further identified by Northern blot. Several miRNAs show the stage-specific expression, of which miR-2943 only expressed in the egg stage, suggesting that miR-2943 may be associated with the development of mosquito eggs.
The present study represents the first global characterization of An. anthropophagus miRNAs in sexual differences and stage-specific functions. A better understanding of the functions of these miRNAs will offer new insights in mosquito biology and has implications for the effective control of mosquito-borne infectious diseases.
Mosquitoes are important vectors in the transmission of viruses and parasitic diseases, which represent an important global public health problem, especially in tropical and subtropical areas [1, 2]. Due to the lack of an efficient vaccine and drug against the pathogens including plasmodium and dengue virus, mosquito control has been the primary strategy for preventing mosquito-borne disease [1, 3]. Anopheles anthropophagus (An. anthropophagus), which belongs to the Hyrcanus group of the genus Anopheles, was identified as a major vector for malaria and filariasis in Asia .
MicroRNAs (miRNAs) are small (approximately 18–24 nucleotides in length), non-coding RNAs that are produced by all animals, plants and some viruses [5–7]. They regulate the expression of cellular genes by guiding the RNA-induced silencing complex (RISC) to target the 3′ UTR of mRNAs for cleavage or translational repression [8, 9], and play a critical role in many cellular processes, including development, differentiation, apoptosis and innate immunity [10–14]. Undoubtedly, they are also implicated in the physiological functions of mosquitoes, such as sexual difference and blood feeding, even in the control of viral and parasitic infection [12–15].
Until now, thousands of miRNAs have been reported in animals, plants and viruses (miRBase: http://www.mirbase.org) . Many miRNAs are highly conserved across divergent species while others are specific to a particular evolutionary lineage [13, 15]. The miRNAs of some members of the mosquitoes, including Anopheles gambiae, Aedes aegypti, Aedes albopictus, Culex quinquefasciatus and Anopheles stephensi have been reported [13, 15–21]. However, there was no miRNAs identified from An. anthropophagus despite its important role in transmitting parasitic disease.
In this study, we performed the first systematic analysis of miRNAs in An. Anthropophagus by using high throughput sequencing and bioinformatics approaches. Due to differences in feeding behaviours between the female and male adult mosquitoes, and female adults playing an important role in transmitting pathogens . The miRNA expression in the sexual difference and specific- stage was substantially analyzed to study the potential role of miRNAs in development and physiological activity. Meanwhile, a number of novel An. anthropophagus mosquito-specific miRNAs were discovered. Our analysis also offered insights into the evolution of conserved and lineage-specific miRNAs in mosquitoes, and has implications for the effective control of mosquito-borne infectious diseases.
Source material and ethics statement
Anopheles anthropophagus (China wild type strain originally from Jiangsu Institute of parasitic diseases prevention, Jiangsu, P.R. China) were reared in a humidified insectary at 26 ± 1°C on a 12 hour light: dark cycle. Adult mosquitoes were kept in a 30 × 30 × 40 cm screened cages and provided constant access to water and glucose-soaked sponges.
ICR mice (Animal Experiment Centre of Wenzhou Medical University) were used in this study to offer a blood meal. The procedure was handled in accordance with good animal practices required by the Animal Ethics Procedures and Guidelines of the People’s Republic of China
Total RNA isolation and small RNA library preparation
Total RNA was prepared from 100 adult male or female An. anthropophagus mosquitoes using trizol (Invitrogen) according to the manufacturer’s protocol. The isolation methods of total RNA and small RNA were unbiased in each sample. All samples were ground in liquid nitrogen and the quality of RNA was detected by using denaturalization agar gels and Du-530 Spectrophotometer (Beckman, Gemany).
The RNA smaller than 200 bp were enriched with the mirVana miRNA isolation kit (Ambion, USA). The small RNA samples were sent to Genergy Bio. (Shanghai, China) for small RNA cloning. The population of miRNAs with a length of 15–30 bp was passively eluted from polyacrylamide gels. The RNA was then precipitated with ethanol and dissolved in water. The small RNAs collected had a poly(A)- tail added to their 3'–OH by poly-(A) polymerase. The 5′-phosphate of the small RNAs were ligated to an RNA adapter. First-strand cDNA synthesis was then performed using an oligo(dT)-linker primer and MMLV-RNase H reverse transcriptase (Promega, USA). The resulting cDNAs were PCR amplified to ~ 25 μg/μl.
High-throughput sequencing and computational analysis An. anthropophagus miRNAs
Primers used for PCR amplification were designed for amplicon sequencing according to the instructions of illumina Hiseq2000 (BGI, China). The PCR-amplified cDNAs were size-selected using electroelution to obtain products of 119–134 bp. These cDNAs were then sequenced by illumina Hiseq2000. Adaptors, low quality reads and reads smaller than 18 nucleotides (nt) were firstly removed from the total small RNA read datasets of male and female adults, respectively. No publically available genome is currently accessible for An. anthropophagus, the genome of the related mosquito and the closely related species in genetic distance An. gambiae and Ae. aegypti were used as a reference genome. The clean read datasets were blasted with BOWTIE software according to the following criteria: a 5′ and 3′ linker match of at least 15 nt and an appropriate length (18–28 nt). The pre-miRNAs and mature miRNAs in the miRBase v.20.0 were searched with BLAST software to identify An. anthropophagus miRNAs. Rfam (10.1) database (http://rfam.sanger.ac.uk/)was used to remove non-miRNAs, including rRNA, tRNA, snRNA, snoRNA.
To identify novel mosquito miRNAs, we used a combination of miRDeep2  and randfold  to ask whether non-annotated sequences mapping to the mosquito genomes demonstrated folding properties of pre-miRNAs hairpins. Each novel miRNA follows both expression and biogenesis criteria set forth for identifying new miRNAs, which include (i) a small RNA of appropriate and discrete length (19–24 nt), (ii) arising from one arm of a hairpin precursor, (iii) presence of the star strand, and (iv) evolutionary conservation [5, 25]. The software RNA-fold was used to predict and analyze the structure of new miRNAs [13, 15].
Sample collections from different developmental stages of An. anthropophagus for northern blots, are briefly described below. Embryo collections were made at 0–24 hours after placing a damp collection cup within a cage. In the larval sample collection, we did not separate early and late larval samples and used one mixed larval sample (I-IV instars) instead. Pupa samples were collected from a pool of varied ages. Adults one to five days following eclosion were collected.
All samples were either directly processed for RNA isolation or flash frozen on liquid nitrogen immediately following collection, then stored at −80°C. The isolation method of total RNA is unbiased for each sample as described above. The quantity of total RNA used for each sample is ~ 20 μg. Northern blots were carried out according to the previous reports . Briefly, total RNA were loaded onto 15% denaturing polyacrylamide gels, and run beside 19 and 23 nucleotide long ssDNA markers. The RNA gels were transferred to Bright- Star-Plus nylon membranes (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. Wash conditions were the same as described before . Antisense 5′ digoxigenin-labelled miRCURY LNA probes were prepared by Exiqon (Vedbaek, Denmark). Probe sequences were as follow, aan-miR-184: CCCTTATCAGTTCTCCGTCCA; aan-miR-989: GTACCACTACGTCACATCACA; aan-miR-998: GAGCTGAATCTCATG GTGCTA; aan-miR-2943: TTTGCCTGCAAGTGCCTACTTAA; aan-miR-1000: ACTG CTGTGTCAGGCAATAT.
Results and discussion
High-throughput sequencing of small RNAs
There is also an elevated population of small RNAs 25–30 nt in length, which may represent piwi-interacting RNAs (piRNAs). piRNAs are known to derive from a Dicer-independent manner from single-stranded precursors and suppress repetitive sequences including transposable elements [26, 27]. Recently piRNAs have been shown to be involved in the innate immunity of mosquitoes [26, 27]. Further investigation of the piRNA pathway in An. anthropophagus mosquitoes will improve our understanding of how this important vector species may defend against pathogens.
Most miRNAs are orthologs of known mosquito miRNAs
The miRNA sequencing reads from male and female adults of An. anthropophagus mosquitoes were aligned to the known miRNA strands present in miRBase (version 20.0). For the pre-miRNAs and mature miRNA library, 8,169,859 clean reads from the male adult library corresponded to 7,481 pre-miRNAs (0.55%) and 3,381 mature miRNAs (0.25%); 12,946,645 clean reads from the female adult library corresponded to 8,680 pre-miRNAs (0.4%) and 3,997 mature miRNAs (0.18%) (Figure 1A).
An. anthropophagus, An. gambiae and Ae. Aegypti are important mosquito-borne disease vectors in tropical and sub-tropical areas. An. gambiae and An. anthropophagus transmit plasmodium and cause malaria, and Ae. Aegypti is the primary vector for yellow fever and dengue fever [13, 18, 28]. The genome sequences of An. gambiae and Ae. Aegypti have been published [13, 18, 20]. However, An. anthropophagus genome is not yet sequenced. Since miRNA sequences are highly conserved between species [13, 15], we aligned sequencing reads to known mature miRNAs and the pre-miRNA library of An. gambiae (aga) and Ae. aegypti (aae). 81 known miRNA expressions were detected in An. anthropophagus (Additional file 1: Table S1), of which 62 miRNAs strands mapped with 100% identity to the library of An. gambiae, and 73 miRNAs strands mapped with 100% identity to the library of Ae. aegypti. 53 miRNAs strands were found to be consistent in these three different mosquito species. Normally, some miRNAs can be found both at 5p and 3p of the same pre-miRNAs . So, the locations of An. anthropophagus miRNAs were subsequently analyzed. We identified 14 miRNAs strands, which located both at of the 5p and 3p of the same pre-miRNAs. However, the miRNAs -3p were only found in the library of Ae. Aegypti. Therefore, these 14 miRNAs and the miRNAs only found in the library of Ae. Aegypti were located according to the genome of Ae. Aegypti, (Additional file 1: Table S1). The other miRNAs were located according to the genome of An. gambiae (Additional file 1: Table S1).
Identification of novel mosquito miRNAs
The novel microRNAs (miRNAs) identified from An. anthropophagus adults
Location at genomeb
Shown in Figure 2 are the sequence alignments of 3 novel miRNAs sequences discovered in this study and the hairpins they form. Sequence blast analysis suggests that miR-N6, miR-N17 and miR-N18 are separately found to be orthologous genes with miR-2a, miR-2942 and miR-2944b in Ae. Aegypti (Figure 2). miR-N4, −N5, −N8, −N9, −N10, −N11, −N12 and miR-N15 are found the orthologous genes in Drosophila melanogaster or Bombyx mori. It is suggested that most miRNAs are relatively conserved in arthropods . miR-N19 is found to be the orthologous gene with miR-33, which is a family of microRNA precursors in several animal species, including humans . Unexpectedly, miR-N21 is found to be the orthologous gene with miR-31, which has been characterised as a tumour suppressor miRNA, with its levels varying in breast cancer cells according to the metastatic state of the tumor [31, 32]. No matches were found for miR-N7.
miRNAs show sexual difference expression in adults
Different expression levels of microRNAs (miRNAs) identified from the male and female adults of An. anthropophagus
The read count of male adulta
The read count of female adultb
The revised read count of male adultc
The ratio of female/maled
Confirmation of mosquito miRNAs
Potential impact of trizol in the loss of low GC miRNAs
Trizol is commonly used for RNA isolation [15, 29, 33]. Kim et al. (2012) found that the structured small RNAs with low GC content are recovered inefficiently when a small number of cells (<500,000) are used for RNA isolation with trizol . Since small RNA needs to base pair with other longer RNA for effective precipitation, it is plausible that the higher concentration of starting material contains more large RNA and is helpful for precipitating the small RNA . In our study, 100 male or female mosquitoes were collected for the RNA isolation, and the quantities of total RNA were ~ 111 μg for male and ~ 159 μg for female. In addition, the quantities of total RNA for egg, larvae and pupae were also more than 100 μg in the trizol extraction. These high quantities of total RNA will help to prevent the small RNA for being lost in the isolation with trizol. Meanwhile, the GC content of aan-miR-1000、-184、-989 and aan-miR-2943 in our study were compared with the low GC miRNA (hsa-miR-141 and miR-200c) in Kim’s paper . As shown in Additional file 3: Table S2, the GC contents of aan-miR-1000 and miR-2943 are very contiguous to miR-141 and miR-200c, which is about 43.5%. So, aan-miR-1000 and miR-2943 also belong to low GC miRNAs. The aan-miR-1000 isolated by trizol can be found in both male and female adults in the sequencing experiments, which further demonstrated that the trizol-related impact on low GC miRNA isolation in our study has been limited in sequencing experiments.
In order to certify that the potential of trizol bias issues were also avoided in Northern experiments, the isolation of RNA from new batches of eggs, larvae, pupae and adult mosquitoes (female and male) were repeated with mirVana miRNA isolation kit (Ambion, USA) rather than trizol. The Northern blot results of miR-1000, −184, −2943 and miR-989 with mirVana miRNA isolation kit are similar to the patterns with trizol isolation (Additional file 4: Figure S2). As described before, miR-1000 and miR-2943 belong to the low GC miRNAs family. miR-1000 is present in embryonic, larval, pupal, and adult stages, and the expression level in male adults is higher than those in female adults by comparing with the blot signal. In contrast, miR-2943 is only present in the embryonic stage. The Northern blot results indicated that the similar expression patterns of the low GC miRNAs were detected both in mirVana miRNA kit isolation and trizol method, and the potential trizol bias issues were overcome in the Northern experiments.
In the present study, the miRNA expression profile in An. anthropophagus were investigated and 102 miRNAs including 81 known and 21 new miRNAs were identified from An. anthropophagus. Furthermore, we also investigated the expression difference of miRNAs in male and female adults as well as different developmental stages, indicating several miRNAs exhibited sexual difference and stage-specific differences. Although the loss of low GC miRNAs was overcome in this study, we should pay attention to the trizol bias issue in future experiments. The present study represented the first global characterization of An. anthropophagus miRNAs, which provides novel resources for better understanding of the biology of the mosquitoes, which, in turn, has implications for the effective control of mosquito-borne infectious diseases.
This work was supported by the National Natural Science Foundation of China (81071389).
- Liu H, Gao X, Liang G: Newly recognized mosquito-associated viruses in mainland China, in the last two decades. Virol J. 2011, 68: 1-13.Google Scholar
- Smith LD, Dushoff J, McKenzie EF: The risk of a mosquito-borne infection in a heterogeneous environment. PLoS Biol. 2004, 11: e368.View ArticleGoogle Scholar
- Hancock AP, Sinkins PS, Godfray J, Charles H: Strategies for introducing wolbachia to reduce transmission of mosquito-borne diseases. PLoS Negl Trop Dis. 2011, 4: e1024-View ArticleGoogle Scholar
- Gao Q, Beebe NW, Cooper RD: Molecular identification of the malaria vectors Anopheles anthropophagus and Anopheles sinensis (Diptera: Culicidae) in central China using polymerase chain reaction and appraisal of their position within the Hyrcanus group. J Med Entomol. 2004, 1: 5-11.View ArticleGoogle Scholar
- Bartel DP: MicroRNAs: genomics, biogenesis, mechanism, and function. Cell. 2004, 116: 281-297. 10.1016/S0092-8674(04)00045-5.View ArticlePubMedGoogle Scholar
- Kim VN, Han J, Siomi MC: Biogenesis of small RNAs in animals. Nat Rev Mol Cell Biol. 2009, 10: 126-139. 10.1038/nrm2632.View ArticlePubMedGoogle Scholar
- Zhang R, Wang YQ, Su B: Molecular evolution of a primate-specific microRNA family. Mol Biol Evol. 2008, 25: 1493-1502. 10.1093/molbev/msn094.View ArticlePubMedGoogle Scholar
- Kim YK, Kim VN: Processing of intronic microRNAs. EMBO J. 2007, 26: 775-783. 10.1038/sj.emboj.7601512.PubMed CentralView ArticlePubMedGoogle Scholar
- Bushati N: MicroRNA functions. Ann Rev Cell Dev Biol. 2007, 23: 175-205. 10.1146/annurev.cellbio.23.090506.123406.View ArticleGoogle Scholar
- Miska EA: How microRNAs control cell division, differentiation and death. Curr Opin Genet Dev. 2005, 15: 563-568. 10.1016/j.gde.2005.08.005.View ArticlePubMedGoogle Scholar
- Liu N, Okamura K, Tyler DM, Phillips MD, Chung WJ, Lai EC: The evolution and functional diversification of animal microRNA genes. Cell Res. 2008, 18: 985-996. 10.1038/cr.2008.278.PubMed CentralView ArticlePubMedGoogle Scholar
- Zhao GH, Xu MJ, Zhu XQ: Identification and characterization of microRNAs in Baylisascaris schroederi of the giant panda. Parasit Vectors. 2013, 6: 216-10.1186/1756-3305-6-216.PubMed CentralView ArticlePubMedGoogle Scholar
- 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-10.1186/1471-2164-10-581.PubMed CentralView ArticlePubMedGoogle Scholar
- Xu MJ, Zhou DH, Nisbet AJ, Huang SY, Fan YF, Zhu XQ: Characterization of mouse brain microRNAs after infection with cyst-forming Toxoplasma gondii. Parasit Vectors. 2013, 1: 154.View ArticleGoogle Scholar
- 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-10.1186/1471-2164-11-119.PubMed CentralView ArticlePubMedGoogle Scholar
- Griffiths-Jones S, Saini HK, van Dongen S, Enright AJ: miRBase: tools for microRNA genomics. Nucleic Acids Res. 2008, 36: D154-D158. 10.1093/nar/gkn221.PubMed CentralView ArticlePubMedGoogle Scholar
- Gu J, Hu W, Wu J, Zheng P, Chen M, James AA, Chen X, Tu Z: miRNA genes of an invasive vector mosquito, Aedes albopictus. PLoS One. 2013, 7: e67638.View ArticleGoogle Scholar
- Winter F, Edaye S, Hüttenhofer A, Brune C: Anopheles gambiae miRNAs as actors of defence reaction against plasmodium invasion. Nucleic Acids Res. 2007, 35: 6953-6962. 10.1093/nar/gkm686.PubMed CentralView ArticlePubMedGoogle Scholar
- Puthiyakunnon S, Yao Y, Li Y, Gu J, Peng H, Chen X: Functional characterization of three MicroRNAs of the Asian tiger mosquito, Aedes albopictus. Parasit Vectors. 2013, 1: 230-View ArticleGoogle Scholar
- Mead EA, Tu Z: Cloning, characterization, and expression of microRNAs from the Asian malaria mosquito, Anopheles stephensi. BMC Genomics. 2008, 9: 244-10.1186/1471-2164-9-244.PubMed CentralView ArticlePubMedGoogle Scholar
- Chatterjee R, Chaudhuri K: An approach for the identification of microRNA with an application to Anopheles gambiae. Acta Biochim Pol. 2006, 2: 303-309.Google Scholar
- 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-851. 10.1007/s00436-009-1470-5.View ArticlePubMedGoogle Scholar
- Friedlander MR, Chen W, Adamidi C, Maaskola J, Einspanier R, Knespel S, Rajewsky N: Discovering microRNAs from deep sequencing data using miRDeep. Nat Biotechnol. 2008, 26: 407-415. 10.1038/nbt1394.View ArticlePubMedGoogle Scholar
- Zuker M: Mfold web server for nucleic acid folding and hybridization prediction. Nucleic Acids Res. 2003, 31: 3406-3415. 10.1093/nar/gkg595.PubMed CentralView ArticlePubMedGoogle Scholar
- Ruby JG, Stark A, Johnston WK, Kellis M, Bartel DP, Lai EC: Evolution, biogenesis, expression, and target predictions of a substantially expanded set of Drosophila microRNAs. Genome Res. 2007, 17: 1850-1864. 10.1101/gr.6597907.PubMed CentralView ArticlePubMedGoogle Scholar
- Léger P, Lara E, Jagla B, Sismeiro O, Mansuroglu Z, Coppée JY, Bonnefoy E, Bouloy M: Dicer-2- and piwi-mediated RNA interference in Rift Valley fever virus-infected mosquito cells. J Virol. 2013, 3: 1631-1648.View ArticleGoogle Scholar
- Morazzani EM, Wiley MR, Murreddu MG, Adelman ZN, Myles KM: Production of virus-derived ping-pong-dependent piRNA-like small RNAs in the mosquito soma. PLoS Pathog. 2012, 1: e1002470-View ArticleGoogle Scholar
- Sinka ME, Bangs MJ, Manguin S, Coetzee M, Mbogo CM, Hemingway J, Patil AP, Temperley WH, Gething PW, Kabaria CW, Okara RM, Van Boeckel T, Godfray HC, Harbach RE, Hay SI: 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-10.1186/1756-3305-3-117.PubMed CentralView ArticlePubMedGoogle Scholar
- Xu MJ, Wang CR, Huang SY, Fu JH, Zhou DH, Chang QC, Zheng X, Zhu XQ: Identification and characterization of microRNAs in the pancreatic fluke Eurytrema pancreaticum. Parasit Vectors. 2013, 6: 25-10.1186/1756-3305-6-25.PubMed CentralView ArticlePubMedGoogle Scholar
- Iwakiri Y: A role of miR-33 for cell cycle progression and cell proliferation. Cell Cycle. 2012, 6: 1057-1058.View ArticleGoogle Scholar
- Valastyan S, Weinberg RA: miR-31: a crucial overseer of tumor metastasis and other emerging roles. Cell Cycle. 2010, 11: 2124-2129.View ArticleGoogle Scholar
- Augoff K, Das M, Bialkowska K, McCue B, Plow EF, Sossey-Alaoui K: miR-31 is a broad regulator of β1-integrin expression and function in cancer cells. Mol Cancer Res. 2011, 11: 1500-1508.View ArticleGoogle Scholar
- Ai L, Xu MJ, Chen MX, Zhang YN, Chen SH, Guo J, Cai YC, Zhou XN, Zhu XQ, Chen JX: Characterization of microRNAs in Taenia saginata of zoonotic significance by Solexa deep sequencing and bioinformatics analysis. Parasitol Res. 2012, 6: 2373-2378.View ArticleGoogle Scholar
- Kim YK, Yeo J, Kim B, Ha M, Kim VN: Short structured RNAs with low GC content are selectively lost during extraction from a small number of cells. Mol Cell. 2012, 46: 893-895. 10.1016/j.molcel.2012.05.036.View ArticlePubMedGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. 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.