Gene expression-based biomarkers for Anopheles gambiae age grading

Predicting the age of field Anopheles mosquitoes using mass spectrometry and deep learning

Mosquito-borne diseases like malaria are rising globally, and improved mosquito vector surveillance is needed. Survival of Anopheles mosquitoes is key for epidemiological monitoring of malaria transmission and evaluation of vector control strategies targeting mosquito longevity, as the risk of pathogen transmission increases with mosquito age. However, the available tools to estimate field mosquito age are often approximate and time-consuming. Here, we show a rapid method that combines matrix-assisted laser desorption/ionization-time-of-flight mass spectrometry with deep learning for mosquito age prediction. Using 2763 mass spectra from the head, legs, and thorax of 251 field-collected Anopheles arabiensis mosquitoes, we developed deep learning models that achieved a best mean absolute error of 1.74 days. We also demonstrate consistent performance at two ecological sites in Senegal, supported by age-related protein changes. Our approach is promising for malaria control and the field of vector biology, benefiting other disease vectors like Aedes mosquitoes.

Evaluation of diffuse reflectance spectroscopy for predicting age, species, and cuticular resistance of Anopheles gambiae s.l under laboratory conditions

Mid-infrared spectroscopy (MIRS) combined with machine learning analysis has shown potential for quick and efficient identification of mosquito species and age groups. However, current technology to collect spectra is destructive to the sample and does not allow targeting specific tissues of the mosquito, limiting the identification of other important biological traits such as insecticide resistance. Here, we assessed the use of a non-destructive approach of MIRS for vector surveillance, micro diffuse reflectance spectroscopy (µDRIFT) using mosquito legs to identify species, age and cuticular insecticide resistance within the Anopheles gambiae s.l. complex. These mosquitoes are the major vectors of malaria in Africa and the focus on surveillance in malaria control programs. Legs required significantly less scanning time and showed more spectral consistence compared to other mosquito tissues. Machine learning models were able to identify An. gambiae and An. coluzzii with an accuracy of 0.73, two ages groups (3 and 10 days old) with 0.77 accuracy and we obtained accuracy of 0.75 when identifying cuticular insecticide resistance. Our results highlight the potential of different mosquito tissues and µDRIFT as tools for biological trait identification on mosquitoes that transmit malaria. These results can guide new ways of identifying mosquito traits which can help the creation of innovative surveillance programs by adapting new technology into mosquito surveillance and control tools.

Identification and Characterization of Small RNA Markers of Age in the Blow Fly Cochliomyia macellaria (Fabricius) (Diptera: Calliphoridae)

Blow fly development is important in decomposition ecology, agriculture, and forensics. Much of the impact of these species is from immature samples, thus knowledge of their development is important to enhance or ameliorate their effects. One application of this information is the estimation of immature insect age to provide temporal information for death investigations. While traditional markers of age such as stage and size are generally accurate, they lack precision in later developmental stages. We used miRNA sequencing to measure miRNA expression, throughout development, of the secondary screwworm, Cochliomyia macellaria (Fabricius) (Diptera: Calliphoridae) and identified 217 miRNAs present across the samples. Ten were identified to be significantly differentially expressed in larval samples and seventeen were found to be significantly differentially expressed in intrapuparial samples. Twenty-eight miRNAs were identified to be differentially expressed between sexes. Expression patterns of two miRNAs, miR-92b and bantam, were qPCR-validated in intrapuparial samples; these and likely food-derived miRNAs appear to be stable markers of age in C. macellaria. Our results support the use of miRNAs for developmental markers of age and suggest further investigations across species and under a range of abiotic and biotic conditions.

Back to the Future: Quantifying Wing Wear as a Method to Measure Mosquito Age

Vector biologists have long sought the ability to accurately quantify the age of wild mosquito populations, a metric used to measure vector control efficiency. This has proven difficult due to the difficulties of working in the field and the biological complexities of wild mosquitoes. Ideal age grading techniques must overcome both challenges while also providing epidemiologically relevant age measurements. Given these requirements, the Detinova parity technique, which estimates age from the mosquito ovary and tracheole skein morphology, has been most often used for mosquito age grading despite significant limitations, including being based solely on the physiology of ovarian development. Here, we have developed a modernized version of the original mosquito aging method that evaluated wing wear, expanding it to estimate mosquito chronological age from wing scale loss. We conducted laboratory experiments using adult Anopheles gambiae held in insectary cages or mesocosms, the latter of which also featured ivermectin bloodmeal treatments to change the population age structure. Mosquitoes were age graded by parity assessments and both human- and computational-based wing evaluations. Although the Detinova technique was not able to detect differences in age population structure between treated and control mesocosms, significant differences were apparent using the wing scale technique. Analysis of wing images using averaged left- and right-wing pixel intensity scores predicted mosquito age at high accuracy (overall test accuracy: 83.4%, average training accuracy: 89.7%). This suggests that this technique could be an accurate and practical tool for mosquito age grading though further evaluation in wild mosquito populations is required.

Assessing Near-Infrared Spectroscopy (NIRS) for Evaluation of Aedes aegypti Population Age Structure

Given that older Aedes aegypti (L.) mosquitoes typically pose the greatest risk of pathogen transmission, the capacity to age grade wild Ae. aegypti mosquito populations would be a valuable tool in monitoring the potential risk of arboviral transmission. Here, we compared the effectiveness of near-infrared spectroscopy (NIRS) to age grade field-collected Ae. aegypti with two alternative techniques—parity analysis and transcript abundance of the age-associated gene SCP1. Using lab-reared mosquitoes of known ages from three distinct populations maintained as adults under laboratory or semi-field conditions, we developed and validated four NIRS models for predicting the age of field-collected Ae. aegypti. To assess the accuracy of these models, female Ae. aegypti mosquitoes were collected from Maricopa County, AZ, during the 2017 and 2018 monsoon season, and a subset were age graded using the three different age-grading techniques. For both years, each of the four NIRS models consistently graded parous mosquitoes as significantly older than nulliparous mosquitoes. Furthermore, a significant positive linear association occurred between SCP1 and NIRS age predictions in seven of the eight year/model combinations, although considerable variation in the predicted age of individual mosquitoes was observed. Our results suggest that although the NIRS models were not adequate in determining the age of individual field-collected mosquitoes, they have the potential to quickly and cost effectively track changes in the age structure of Ae. aegypti populations across locations and over time.

Use of transcriptional age grading technique to determine the chronological age of Sri Lankan Aedes aegypti and Aedes albopictus females

BACKGROUND

Aedes aegypti and Ae. albopictus are important vectors of human diseases such as dengue, chikungunya, and zika. In Sri Lanka, they have been responsible for transmitting dengue virus. One of the most important parameters influencing the likelihood of arbovirus transmission is the age structure of the mosquito population. However, mosquito age is difficult to measure with accuracy. This study aims to construct multivariate calibration models using the transcriptional abundance of three age-responsive genes: Ae15848 (calcium-binding protein), Ae8505 (structural component of cuticle), and Ae4274 (fizzy cell cycle/cell division cycle 20).

METHODS

The transcriptional age-grading technique was applied to determine the chronological age of Ae. aegypti and Ae. albopictus female mosquito populations from Sri Lanka using the age-responsive genes Ae15848, Ae8505, and Ae4274. Furthermore, Ae. aegypti samples obtained from colonies reared at two temperatures (23 and 27 °C) were used to investigate the influence of temperature on this age-grading technique. Expression levels of these three genes were quantified using reverse transcription qualitative PCR (qRT-PCR), and results were normalized against the housekeeping gene ribosomal gene S17 (RpS17).

RESULTS

The expression of Ae15848 and Ae8505 decreased with the age of mosquitoes and showed the most significant and consistent change while expression of Ae4274 increased with age. The multivariate calibration models showed > 80% correlation between expression of these age-responsive genes and the age of female mosquitoes at both temperatures. At 27 °C the accuracy of age predictions using the models was 2.19 (± 1.66) days and 2.58 (± 2.06) days for Ae. aegypti and Ae. albopictus females, respectively. The accuracy of the model for Ae. aegypti at 23 °C was 3.42 (± 2.74) days.

CONCLUSIONS

An adult rearing temperature difference of 4 °C (23-27 °C) did not significantly affect the age predictions. The calibration models created during this study could be successfully used to estimate the age of wild Ae. aegypti and Ae. albopictus mosquitoes from Sri Lanka.

Mosquito Age Grading and Vector-Control Programmes

An ability to characterize the age of mosquito populations could provide cost-effective and compelling entomological evidence for the potential epidemiological impacts of vector control. The average age of a mosquito population is the most important determinant of vectorial capacity and the likelihood of disease transmission. Yet, despite decades of research, defining the age of a wild-caught mosquito remains a challenging, impractical, and unreliable process. Emerging chemometric and existing transcriptional approaches may overcome many of the limitations of current morphological techniques, but their utility in terms of field-based monitoring programmes remains largely untested. Herein, we review the potential advantages and disadvantages of new and existing age-grading tools in an operational context.

Analysis of natural female post-mating responses of Anopheles gambiae and Anopheles coluzzii unravels similarities and differences in their reproductive ecology

Anopheles gambiae and An. coluzzii, the two most important malaria vectors in sub-Saharan Africa, are recently radiated sibling species that are reproductively isolated even in areas of sympatry. In females from these species, sexual transfer of male accessory gland products, including the steroid hormone 20-hydroxyecdysone (20E), induces vast behavioral, physiological, and transcriptional changes that profoundly shape their post-mating ecology, and that may have contributed to the insurgence of post-mating, prezygotic reproductive barriers. As these barriers can be detected by studying transcriptional changes induced by mating, we set out to analyze the post-mating response of An. gambiae and An. coluzzii females captured in natural mating swarms in Burkina Faso. While the molecular pathways shaping short- and long-term mating-induced changes are largely conserved in females from the two species, we unravel significant inter-specific differences that suggest divergent regulation of key reproductive processes such as egg development, processing of seminal secretion, and mating behavior, that may have played a role in reproductive isolation. Interestingly, a number of these changes occur in genes previously shown to be regulated by the sexual transfer of 20E and may be due to divergent utilization of this steroid hormone in the two species.

Quantitative real-time PCR analysis of Anopheles dirus TEP1 and NOS during Plasmodium berghei infection, using three reference genes

Quantitative reverse transcription PCR (qRT-PCR) has been an integral part of characterizing the immunity of Anopheles mosquitoes towards Plasmodium invasion. Two anti-Plasmodium factors of Anopheles, thioester-containing protein 1 (TEP1) and nitric oxide synthase (NOS), play a role in the refractoriness of Anopheles towards Plasmodium infection and are generally expressed during infection. However, these are less studied in Anopheles dirus, a dominant malaria vector in Southeast Asia. Furthermore, most studies used a single reference gene for normalization during gene expression analysis without proper validation. This may lead to erroneous quantification of expression levels. Therefore, the present study characterized and investigated the expression profiles of TEP1 and NOS of Anopheles dirus during P. berghei infection. Prior to that, the elongation factor 1-alpha (EF1), actin 1 (Act) and ribosomal protein S7 (S7) genes were validated for their suitability as a set of reference genes. TEP1 and NOS expressions in An. dirus were found to be significantly induced after P. berghei infection.

Metabolomics of aging assessed in individual parasitoid wasps

Metabolomics studies of low-biomass organisms, such as small insects, have previously relied on the pooling of biological samples to overcome detection limits, particularly using NMR. We show that the differentiation of metabolite profiles of individual 1 mg parasitoid wasps of different ages is possible when using a modified sample preparation and a combination of untargeted NMR and LC-MS based metabolomics. Changes were observed between newly emerged and older wasps in glycerolipids, amino acids and circulatory sugars. This advance in chemical profiling has important implications for the study of the behaviour and ecology of parasitoids and many other species of small organisms because predictions and observations are typically made at the level of the individual. Thus, the metabolomic state of low-biomass individuals can now be related to their behaviour and ecological performance. We discuss specifically the utility of age-related metabolomic profiling but our new approach can be applied to a wide range of biological research.

Estimating malaria transmission from humans to mosquitoes in a noisy landscape

A basic quantitative understanding of malaria transmission requires measuring the probability a mosquito becomes infected after feeding on a human. Parasite prevalence in mosquitoes is highly age-dependent, and the unknown age-structure of fluctuating mosquito populations impedes estimation. Here, we simulate mosquito infection dynamics, where mosquito recruitment is modelled seasonally with fractional Brownian noise, and we develop methods for estimating mosquito infection rates. We find that noise introduces bias, but the magnitude of the bias depends on the ‘colour' of the noise. Some of these problems can be overcome by increasing the sampling frequency, but estimates of transmission rates (and estimated reductions in transmission) are most accurate and precise if they combine parity, oocyst rates and sporozoite rates. These studies provide a basis for evaluating the adequacy of various entomological sampling procedures for measuring malaria parasite transmission from humans to mosquitoes and for evaluating the direct transmission-blocking effects of a vaccine.

Molecular biomarkers for chronological age in animal ecology

The chronological age of an individual animal predicts many of its biological characteristics, and these in turn influence population-level ecological processes. Animal age information can therefore be valuable in ecological research, but many species have no external features that allow age to be reliably determined. Molecular age biomarkers provide a potential solution to this problem. Research in this area of molecular ecology has so far focused on a limited range of age biomarkers. The most commonly tested molecular age biomarker is change in average telomere length, which predicts age well in a small number of species and tissues, but performs poorly in many other situations. Epigenetic regulation of gene expression has recently been shown to cause age-related modifications to DNA and to cause changes in abundance of several RNA types throughout animal lifespans. Age biomarkers based on these epigenetic changes, and other new DNA-based assays, have already been applied to model organisms, humans and a limited number of wild animals. There is clear potential to apply these marker types more widely in ecological studies. For many species, these new approaches will produce age estimates where this was previously impractical. They will also enable age information to be gathered in cross-sectional studies and expand the range of demographic characteristics that can be quantified with molecular methods. We describe the range of molecular age biomarkers that have been investigated to date and suggest approaches for developing the newer marker types as age assays in nonmodel animal species.

Candidate biomarkers for mosquito age-grading identified by label-free quantitative analysis of protein expression in Aedes albopictus females

We applied a "shotgun" approach based on nanoliquid chromatography-high resolution mass spectrometry associated to label free quantification (LFQ) to identify proteins varying with age, independently from the physiological state, in Aedes albopictus, a mosquito species which in the last decades invaded temperate regions in North America and Europe, creating concerns for associated high nuisance and risk of arbovirus transmission. The combined "shotgun" and LFQ approach was shown to be highly suitable to simultaneously compare several biological samples, as needed in a study aimed to analyze different age-groups and physiological states of adult mosquito females. The results obtained represent the first wide-scale analysis of protein expression in Ae. albopictus females: >1000 and 665 proteins were identified from few micrograms of crude protein extracts of mosquito heads and thoraxes, respectively. Six of these proteins were shown to significantly vary from 2- to 16-day-old females, independently from their physiological state (i.e. virgin, mated, host-seeking, blood-fed, and gravid).

BIOLOGICAL SIGNIFICANCE

Mosquito-borne diseases, such as malaria, dengue and other arboviroses, are a persistent cause of global mortality and morbidity, affecting hundreds of thousands of people. Billions of people living in tropical areas are at risk of being bitten every day by an infective mosquito female and the spread of tropical species such as Aedes albopictus to temperate areas is creating alarm in the northern hemisphere. Mosquito longevity is a critical factor affecting mosquito-borne pathogen transmission cycles and the mosquito capacity to transmit pathogens. However, large scale analyses of the age structure of mosquito field populations is hampered by the lack of optimal age-grading approaches. Our findings open new perspectives for the development of reliable, simple and cheap protein-based assays to age-grade Ae. albopictus females and, most likely, other mosquito species of higher medical relevance, such as the main dengue vector, Aedes aegypti, and the major Afrotropical malaria vectors. These assays would greatly contribute to epidemiological studies aimed at defining the actual vectorial capacity of a given mosquito species. Moreover, they would be very valuable in assessing the effectiveness of mosquito control interventions based on the relative ratio between young and old individuals before and after the intervention.

Proteomic changes occurring in the malaria mosquitoes Anopheles gambiae and Anopheles stephensi during aging

The age of mosquitoes is a crucial determinant of their ability to transmit pathogens and their resistance to insecticides. We investigated changes to the abundance of proteins found in heads and thoraces of the malaria mosquitoes Anopheles gambiae and Anopheles stephensi as they aged. Protein expression changes were assessed using two-dimensional difference gel electrophoresis and the identity of differentially expressed proteins was determined by using either matrix-assisted laser desorption ionization tandem time-of-flight mass spectrometry or capillary high-pressure liquid chromatography coupled with a linear ion-trap (LTQ)-Orbitrap XL hybrid mass spectrometer. Protein biomarkers were validated by semi quantitative Western blot analysis. Nineteen and nine age dependent protein spots were identified for A. stephensi and A. gambiae, respectively. Among the proteins down-regulated with age were homologs of ADF/Cofilin, cytochome c1, heat shock protein-70 and eukaryotic translation initiation factor 5A (eIF5a). Proteins up-regulated with age included probable methylmalonate-semialdehyde dehydrogenase, voltage-dependent anion-selective channel and fructose bisphosphate aldolase. Semi quantitative Western blot analysis confirmed expression patterns observed by 2-D DIGE for eIF5a and ADF/Cofilin. Further work is recommended to determine whether these biomarkers are robust to infection, blood feeding and insecticide resistance. Robust biomarkers could then be incorporated into rapid diagnostic assays for ecological and epidemiological studies.

BIOLOGICAL SIGNIFICANCE

In this study, we have identified several proteins with characteristic changes in abundance in both A. gambiae and A. stephensi during their aging process. These changes may highlight underlying mechanisms beneath the relationship between mosquito age and factors affecting Plasmodium transmission and mosquito control. The similarity of changes in protein abundance between these species and the primary dengue vector Aedes aegypti, has revealed conserved patterns of aging-specific protein regulation.

Using a near-infrared spectrometer to estimate the age of anopheles mosquitoes exposed to pyrethroids

We report on the accuracy of using near-infrared spectroscopy (NIRS) to predict the age of Anopheles mosquitoes reared from wild larvae and a mixed age-wild adult population collected from pit traps after exposure to pyrethroids. The mosquitoes reared from wild larvae were estimated as <7 or ≥7 d old with an overall accuracy of 79%. The age categories of Anopheles mosquitoes that were not exposed to the insecticide papers were predicted with 78% accuracy whereas the age categories of resistant, susceptible and mosquitoes exposed to control papers were predicted with 82%, 78% and 79% accuracy, respectively. The ages of 85% of the wild-collected mixed-age Anopheles were predicted by NIRS as ≤8 d for both susceptible and resistant groups. The age structure of wild-collected mosquitoes was not significantly different for the pyrethroid-susceptible and pyrethroid-resistant mosquitoes (P = 0.210). Based on these findings, NIRS chronological age estimation technique for Anopheles mosquitoes may be independent of insecticide exposure and the environmental conditions to which the mosquitoes are exposed.