Identifying Canadian mosquito species through DNA barcodes

Morphometric and molecular characterization of Channa marulius from Riverine system of Punjab, Pakistan

OBJECTIVE

Channidae family, are major freshwater fish species amongst the local aquatic fauna of Pakistan, while, there is limited availability of local data on their molecular identification and phylogenetic analysis.

METHODS

Channa species were collected from four different geographical sites in the tertiary of Punjab province on the Indus and Chenab rivers of Pakistan. Morphometric records and molecular techniques were used to determine the intraspecific variations among populations of Channa marulius. Mitochondrial DNA was extracted from the flesh of C. marulius, while, COI gene was used for molecular identification and variation levels were estimated by using Principal Component Analysis.

RESULTS

Data recorded on the basis of morphometric parameters clearly divided the C. marulius of different locations into two distinct categories, which accounted for a cumulative variability of 97.6%. Non-significance (P < 0.05) among the C. marulius showed that it contains a unique control haplotype localized within the sub-population. The intra-species distance ranged from 0.000 to 0.001 for four different populations, in contrast, the sequences retrieved from the NCBI database exhibited a range span of 0.000-0.003, while, sequence diversity ranged from 0.000 to 0.006 for this intra-specific comparison. The cladogram was also constructed for C. marulius of different geographical locations for observation of phylogenetic relationship. The conclusion drawn from the phylogenetic analysis of C. marulius populations used in this study, contributes significantly to the understanding of genetic variations within populations of this species. The findings provide valuable insight to devise conservation strategies in fisheries management programs in Pakistan.

Population Genetics of Culex tritaeniorhynchus (Diptera: Culicidae) in Türkiye

PURPOSE

Mosquitoes are important vectors of pathogens that can affect humans and animals. Culex tritaeniorhynchus is an important vector of arboviruses such as Japanese encephalitis virus, West Nile virus among various human and animal communities. These diseases are of major public health concern and can have huge economic and health burdens in prevalent countries. Although populations of this important mosquito species have been detected in the Mediterranean and Aegean regions of Türkiye; little is known about its population structure. Our study is to examine the population genetics and genetic composition of Cx. tritaeniorhynchus mosquitoes collected from several localities using cytochrome oxidase subunit I (COI) and the NADH dehydrogenase subunit 5 genes (ND5). This is the first extensive study of Cx. tritaeniorhynchus in the mainland Türkiye with sampling spanning many of provinces.

METHODS

In this study, DNA extraction, amplification of mitochondrial COI and ND5 genes and population genetic analyses were performed on ten geographic populations of Culex tritaeniorhynchus in the Aegean and Mediterranean region of Türkiye.

RESULTS

Between 2019 and 2020, 96 samples were collected from 10 geographic populations in the Aegean and Mediterranean regions; they were molecularly analyzed and 139 sequences (50 sequence for COI and 89 sequence for ND5) were used to determine the population structure and genetic diversity. For ND5 gene region, the samples produced 24 haplotypes derived from 15 variable sites and for COI gene region, 43 haplotypes were derived from 17 variable sites. The haplotype for both gene regions was higher than nucleotide diversity. Haplotype phylogeny revealed two groups present in all populations. AMOVA test results show that the geographical populations were the same for all gene regions. Results suggest that Cx. tritaeniorhynchus is a native population in Türkiye, the species is progressing towards speciation and there is no genetic differentiation between provinces and regions.

CONCLUSION

This study provides useful information on the molecular identifcation and genetic diversity of Cx. tritaeniorhynchus; these results are important to improve mosquito control programs.

The monsoon-associated equine South African pointy mosquito 'Aedes caballus'; the first comprehensive record from southeastern Iran with a description of ecological, morphological, and molecular aspects

The equine South African pointy vector mosquito, Aedes caballus, poses a significant threat to human health due to its capacity for transmitting arboviruses. Despite favorable climate for its existence in southeast Iran, previous records of this species in the area have indicated very low abundance. This comprehensive field and laboratory study aimed to assess its current adult population status in this region, utilizing a combination of ecological, morphological and molecular techniques. Four distinct types of traps were strategically placed in three fixed and two variable mosquito sampling sites in the southern strip of Sistan and Baluchistan Province. Subsequently, DNA was extracted from trapped mosquitoes and subjected to PCR amplification using the molecular markers COI, ITS2, and ANT. In total, 1734 adult Ae. caballus specimens were collected from rural areas, with the majority being captured by CO2-baited bednet traps. A notable increase in the abundance of this species was observed following rainfall in February. The genetic analysis revealed multiple haplotypes based on COI and ITS2 sequences, with COI gene divergence at 0.89%, and ITS2 sequence divergence at 1.6%. This suggests that previous challenges in morphological identification may have led to misidentifications, with many adults previously classified as Ae. vexans potentially being Ae. caballus. The findings of this study hold significant implications for public health authorities, providing valuable insights for integrated and targeted vector control and disease management efforts.

Genetic variability of Aedes aegypti (Diptera: Culicidae) in El Salvador and Honduras: presence of a widespread haplotype and implications for mosquito control

BACKGROUND

This study examined population genetics of Aedes aegypti in El Salvador and Honduras, two adjacent countries in Central America. Aedes aegypti is associated with yellow fever, dengue, chikungunya, and Zika. Each year, thousands of cases of dengue are typically reported in El Salvador and Honduras.

METHODS

In El Salvador, collections were obtained from five Departments. In Honduras, samples were obtained from six municipalities in four Departments. Mitochondrial DNA cytochrome oxidase I (COI) was sequenced, and consensus sequences were combined with available sequences from El Salvador to determine haplotype number, haplotype diversity, nucleotide diversity, and Tajima's D. A haplotype network was produced to examine the relationship between genotypes.

RESULTS

In El Salvador, there were 17 haplotypes, while in Honduras there were 4 haplotypes. In both El Salvador and Honduras, Haplotype 1 is most abundant and widespread. In El Salvador, haplotype H2 was also widespread in 10 of 11 sampled municipalities, but it was not present in Honduras. The capital of El Salvador (San Salvador) and the eastern region of ES had the highest haplotype diversity of regions sampled.

CONCLUSIONS

Haplotype 1 and H2 each belong to different phylogenetic lineages of Ae. aegypti. The most geographically widespread haplotype (H1) may have been present the longest and could be a remnant from previous eradication programs. These data may contribute to future control programs for Ae. aegypti in the two countries.

A DNA barcode reference library of Croatian mosquitoes (Diptera: Culicidae): implications for identification and delimitation of species, with notes on the distribution of potential vector species

BACKGROUND

Mosquitoes pose a risk to human health worldwide, and correct species identification and detection of cryptic species are the most important keys for surveillance and control of mosquito vectors. In addition to traditional identification based on morphology, DNA barcoding has recently been widely used as a complementary tool for reliable identification of mosquito species. The main objective of this study was to create a reference DNA barcode library for the Croatian mosquito fauna, which should contribute to more accurate and faster identification of species, including cryptic species, and recognition of relevant vector species.

METHODS

Sampling was carried out in three biogeographical regions of Croatia over six years (2017-2022). The mosquitoes were morphologically identified; molecular identification was based on the standard barcoding region of the mitochondrial COI gene and the nuclear ITS2 region, the latter to identify species within the Anopheles maculipennis complex. The BIN-RESL algorithm assigned the COI sequences to the corresponding BINs (Barcode Index Number clusters) in BOLD, i.e. to putative MOTUs (Molecular Operational Taxonomic Units). The bPTP and ASAP species delimitation methods were applied to the genus datasets in order to verify/confirm the assignment of specimens to specific MOTUs.

RESULTS

A total of 405 mosquito specimens belonging to six genera and 30 morphospecies were collected and processed. Species delimitation methods assigned the samples to 31 (BIN-RESL), 30 (bPTP) and 28 (ASAP) MOTUs, with most delimited MOTUs matching the morphological identification. Some species of the genera Culex, Aedes and Anopheles were assigned to the same MOTUs, especially species that are difficult to distinguish morphologically and/or represent species complexes. In total, COI barcode sequences for 34 mosquito species and ITS2 sequences for three species of the genus Anopheles were added to the mosquito sequence database for Croatia, including one individual from the Intrudens Group, which represents a new record for the Croatian mosquito fauna.

CONCLUSION

We present the results of the first comprehensive study combining morphological and molecular identification of most mosquito species present in Croatia, including several invasive and vector species. With the exception of some closely related species, this study confirmed that DNA barcoding based on COI provides a reliable basis for the identification of mosquito species in Croatia.

Large-Scale Cytochrome C Oxidase Subunit I Gene Data Analysis for the Development of a Multiplex Polymerase Chain Reaction Test Capable of Identifying Biting Midge Vector Species and Haplotypes (Diptera: Ceratopogonidae) of the Culicoides Subgenus Avaritia Fox, 1955

The emergence of culicoid-transmitted bluetongue and Schmallenberg viruses in several European countries demonstrated the ability of indigenous biting midge species to transmit pathogens. Entomologic research programs identified members of the Obsoletus Group (Culicoides subgenus Avaritia) as keyplayers in disease epidemiology in Europe. However, morphological identification of potential vectors is challenging due to the recent discovery of new genetic variants (haplotypes) of C. obsoletus sensu stricto (s.s.), forming distinct clades. In this study, 4422 GenBank entries of the mitochondrial cytochrome c oxidase subunit I (COI) gene of subgenus Avaritia members of the genus Culicoides were analyzed to develop a conventional multiplex PCR, capable of detecting all vector species and clades of the Western Palearctic in this subgenus. Numerous GenBank entries incorrectly assigned to a species were identified, analyzed and reassigned. The results suggest that the three C. obsoletus clades represent independent species, whereas C. montanus should rather be regarded as a genetic variant of C. obsoletus s.s. Based on these findings, specific primers were designed and validated with DNA material from field-caught biting midges which achieved very high diagnostic sensitivity (100%) when compared to an established reference PCR (82.6%).

Outline-based geometric morphometrics: Wing cell differences for mosquito vector classification in the Tanaosri mountain range, Thailand

Recent studies have revealed taxonomic signals within the wing cells of certain mosquito species. In our study, wing cell differentiation among mosquito vectors from the Tanaosri mountain range in Thailand was evaluated using the outline-based geometric morphometric (GM) approach. Our focus was on four specific wing cells for GM analysis: the wing contour (external cell), the second submarginal cell (internal cell 1), the first posterior cell (internal cell 2), and the third posterior cell (internal cell 3). Before proceeding with the GM approach, the identity of seven mosquito genera and 21 species was confirmed using molecular techniques. Our validated classification tests demonstrated that the performance of mosquito species classification varies according to genus. Notably, three Aedes species exhibited the highest accuracy for both internal cell 2 and internal cell 3, each registering a score of 93.20 %. In the case of two Mansonia species, the wing contour displayed a remarkable accuracy of 98.57 %. Consequently, we suggest the use of the outline-based GM approach, particularly focusing on the wing contour, for differentiating Mansonia annulifera and Mansonia uniformis. In contrast, the highest accuracy for classifying Culex species was found in internal cell 1, at 75.51 %, highlighting the challenges due to similarities in wing cells within this genus. These findings provide a guideline for future applications of the outline-based GM approach, focusing on wing cells, as an alternative method to classify mosquito vector species.

Mitochondrial and microbial diversity of the invasive mosquito vector species Culex tritaeniorhynchus across its extensive inter-continental geographic range

Background

Culex (Cx.) tritaeniorhynchus is an invasive mosquito species with an extensive and expanding inter-continental distribution, currently reported across Asia, Africa, the Middle East, Europe and now Australia. It is an important vector of medical and veterinary pathogens which cause significant morbidity and mortality in human and animal populations. Across regions endemic for Japanese encephalitis virus (JEV), Cx. tritaeniorhynchus is considered the major vector and has also been shown to contribute to the transmission of several other zoonotic arboviruses including Rift Valley fever virus (RVFV) and West Nile virus (WNV).

Methods

In this study, we used laboratory vector competence experiments to determine if Cx. tritaeniorhynchus from a Southern European population were competent JEV vectors. We also obtained samples from multiple geographically dispersed Cx. tritaeniorhynchus populations from countries within Europe, Africa, Eurasia and Asia to perform phylogenetic analysis to measure the level of mitochondrial divergence using the cytochrome oxidase subunit 1 ( CO1) gene. We also undertook bacterial 16S rRNA gene amplicon sequencing to determine microbial diversity and used multi-locus sequence typing (MLST) to determine any evidence for the presence of strains of the naturally occurring endosymbiotic bacterium Wolbachia.

Results

Cx. tritaeniorhynchus from a Greek population were shown be be competent vectors of JEV with high levels of virus present in saliva. We found a signficant level of mitochondrial genetic diversity using the mosquito CO1 gene between geographically dispersed populations. Furthermore, we report diverse microbiomes identified by 16S rRNA gene amplicon sequencing within and between geographical populations. Evidence for the detection of the endosymbiotic bacteria Wolbachia was confirmed using Wolbachia-specific PCR and MLST.

Conclusions

This study enhances our understanding of the diversity of Cx. tritaeniorhynchus and the associated microbiome across its inter-continental range and highlights the need for greater surveillance of this invasive vector species in Europe.

The insupportable validity of mosquito subspecies (Diptera: Culicidae) and their exclusion from culicid classification

Beginning about 80 years ago, the recognition of morphological varieties of mosquitoes was gradually replaced by the recognition of subspecies. As an examination of revisionary and detailed taxonomic studies of mosquitoes clearly shows, subspecies are untenable concepts which have been synonymized with nominotypical forms or recognized as distinct species. Thus, from our perspective, subspecies is not a functional or practical taxonomic rank. Consequently, in this study we critically assessed the taxonomic status of the 120 nominal taxa distinguished as subspecies before now to determine whether they should be recognized as separate species or synonymous names. As a result, 96 subspecies are formally elevated to specific rank, 22 are relegated to synonymy with nominotypical forms, one is considered a nomen dubium, one a species inquirenda and the names of four nominal species regarded as synonyms are revalidated. The subspecies and their new status are listed in a conspectus. The revalidated species include Anopheles argentinus (Brèthes, 1912), from synonymy with An. pseudopunctipennis Theobald, 1901c; An. peruvianus Tamayo, 1907, from synonymy with An. pseudopunctipennis as nomen dubium; Culex major Edwards, 1935, from synonymy with Cx. annulioris consimilis Newstead, 1907; and Trichoprosopon trichorryes (Dyar & Knab, 1907), from synonymy with Tr. compressum Lutz, 1905. Additionally, the type locality of Anopheles sergentii Theobald, 1907 is restricted to El Outaya, Biskra Province, Algeria. A complete list of species to be retained, added to or removed from the Encyclopedia of Life, with a few corrections, is provided.

Molecular and morphometric differentiation of secondary filariasis vector Coquillettidia mosquitoes (Diptera: Culicidae) in Thailand

Coquillettidia mosquitoes are important nuisance-biting pests and a vector of brugian filariasis in Thailand. However, comprehensive information about these mosquitoes remains unavailable such as molecular and morphometric differences among species. The lack of vector knowledge on Coquillettidia species could affect future disease control. This study aims to investigate differences in molecular variations based on mitochondrial cytochrome oxidase subunit I (COI) gene and wing geometric traits of three Coquillettidia species, namely Cq. crassipes, Cq. nigrosignata, and Cq. ochracea in Thailand. The results of molecular analyses revealed the differences among three Coquillettidia species. The genetic difference measure based on the Kimura two-parameter model among three Coquillettidia species showed low intraspecific distances (0%-3.05%) and large interspecific distances (10.10%-12.41%). The values of intra- and inter-genetic differences of three Coquillettidia species did not overlap which showed the existence of a barcoding gap indicating the efficiency of the identification based on the COI gene. As with molecular analysis, the landmark-based geometric morphometrics approach based on wing shape analysis indicated three distinct species groups which were supported by the high total performance score of cross-validated classification (97.16%). These results provide the first evidence of taxonomic signal based on molecular and wing geometric differences to support species identification and biological variations of Coquillettidia mosquitoes in Thailand for understanding these rare vector mosquitoes in depth and leading to effective further mosquito control.

DNA Barcodes of Mansonia (Mansonia) Blanchard, 1901 (Diptera, Culicidae)

Females of the genus Mansonia feed on the blood of humans, livestock, and other vertebrates to develop their eggs. The females' biting behavior may cause severe disturbance to blood hosts, with a negative impact on public health and economics. Certain species have been identified as potential or effective disease vectors. The accurate species identification of field-collected specimens is of paramount importance for the success of monitoring and control strategies. Mansonia (Mansonia) morphological species boundaries are blurred by patterns of intraspecific heteromorphism and interspecific isomorphism. DNA barcodes can help to solve taxonomic controversies, especially if combined with other molecular tools. We used cytochrome c oxidase subunit I (COI) gene 5' end (DNA barcode) sequences to identify 327 field-collected specimens of Mansonia (Mansonia) spp. The sampling encompassed males and females collected from three Brazilian regions and previously assigned to species based on their morphological characteristics. Eleven GenBank and BOLD sequences were added to the DNA barcode analyses. Initial morphospecies assignments were mostly corroborated by the results of five clustering methods based on Kimura two-parameter distance and maximum likelihood phylogeny. Five to eight molecular operational taxonomic units may represent taxonomically unknown species. The first DNA barcode records for Mansonia fonsecai, Mansonia iguassuensis, and Mansonia pseudotitillans are presented.

High Levels of Diversity in Anopheles Subgenus Kerteszia Revealed by Species Delimitation Analyses

The Anopheles subgenus Kerteszia is a poorly understood group of mosquitoes that includes several species of medical importance. Although there are currently twelve recognized species in the subgenus, previous studies have shown that this is likely to be an underestimate of species diversity. Here, we undertake a baseline study of species delimitation using the barcode region of the mtDNA COI gene to explore species diversity among a geographically and taxonomically diverse range of Kerteszia specimens. Beginning with 10 of 12 morphologically identified Kerteszia species spanning eight countries, species delimitation analyses indicated a high degree of cryptic diversity. Overall, our analyses found support for at least 28 species clusters within the subgenus Kerteszia. The most diverse taxon was Anopheles neivai, a known malaria vector, with eight species clusters. Five other species taxa showed strong signatures of species complex structure, among them Anopheles bellator, which is also considered a malaria vector. There was some evidence for species structure within An. homunculus, although the results were equivocal across delimitation analyses. The current study, therefore, suggests that species diversity within the subgenus Kerteszia has been grossly underestimated. Further work will be required to build on this molecular characterization of species diversity and will rely on genomic level approaches and additional morphological data to test these species hypotheses.

Myzomyia and Pyretophorus series of Anopheles mosquitoes acting as probable vectors of the goat malaria parasite Plasmodium caprae in Thailand

Unlike malaria parasites in humans, non-human primates, rodents, and birds, ungulate malaria parasites and their vectors have received little attention. As a result, understanding of the hosts, vectors, and biology of ungulate malaria parasites has remained limited. In this study, we aimed to identify the vectors of the goat malaria parasite Plasmodium caprae. A total of 1019 anopheline and 133 non-anopheline mosquitoes were collected from goat farms in Thailand, where P. caprae-infected goats were discovered. Anopheline mosquitoes were identified using molecular biological methods that target the cytochrome c oxidase subunit 1 (cox1), the cytochrome c oxidase subunit 2 (cox2) genes, and the internal transcribed spacer 2 (ITS2) region. Pool and individual mosquitoes were tested for P. caprae using the head-thorax parts that contain the salivary glands, with primers targeting three genetic markers including cytochrome b, cytochrome c oxidase subunit 1, and 18S small subunit ribosomal RNA genes. Additionally, goat blood samples were collected concurrently with mosquito surveys and screened to determine the status of malaria infection. This study revealed nine mosquito species belonging to six groups on goat farms, including Hyrcanus, Barbirostris, Subpictus, Funestus, Tessellatus, and Annularis. The DNA of P. caprae was detected in Anopheles subpictus and Anopheles aconitus. This is the first time An. subpictus and An. aconitus have been implicated as probable vectors of P. caprae.

Modification of the Folmer primers for the cytochrome c oxidase gene facilitates identification of mosquitoes

BACKGROUND

Accurate identification of mosquito species is essential for the development and optimization of strategies to control mosquitoes and mosquito-borne diseases. Problems with the morphological identification of mosquito species have led to the use of molecular identification techniques, in particular the Folmer cytochrome c oxidase subunit I (COI) PCR system (FCOS), originally designed to identify a range of other invertebrates.

METHODS

As there can be difficulties identifying mosquitoes using FCOS, we re-evaluated the FCOS primers and developed a new COI-based SYBR PCR (the Auburn COI system-AUCOS) to improve the molecular identification of mosquitoes. Sequence data in GenBank for 33 species from 10 genera of mosquitoes were used to develop our AUCOS primers. Two molecular assays (AUCOS, FCOS) and morphological identification were carried out on mosquitoes collected from the field in Auburn, Alabama (USA) and on Saint Kitts.

RESULTS

With a convenience sample of individual mosquitoes comprising 19 species from six genera in Saint Kitts (n = 77) and Auburn (n = 48), our AUCOS provided higher-quality sequence data than FCOS. It also proved more sensitive than FCOS, successfully amplifying 67.5% (85/126) as opposed to 16.7% (21/126) of the samples. The species determined by morphology, or genus with damaged samples, matched that as determined by AUCOS for 84.9% (62/73) of the samples. Morphological classification was confirmed by FCOS with 81.0% (17/21) of samples producing utilizable sequences. While both FCOS and AUCOS correctly identified all the Aedes, Anopheles, Deinocerites, and Uranotaenia species in the study, identification of Culex species was less successful with both methods: 50.0% (3/6) by FCOS and 35.7% (5/14) by AUCOS.

CONCLUSIONS

The AUCOS DNA barcoding system for mosquito species described in this study is superior to the existing FCOS for the identification of mosquito species. As AUCOS and FCOS amplify the same variable region of the COI, the large amount of existing data on GenBank can be used to identify mosquito species with sequences produced by either PCR.

Mitochondrial DNA barcoding of mosquito species (Diptera: Culicidae) in Thailand

The correct identification of mosquito species is important for effective mosquito vector control. However, the standard morphological identification of mosquito species based on the available keys is not easy with specimens in the field due to missing or damaged morphological features during mosquito collections, often leading to the misidentification of morphologically indistinguishable. To resolve this problem, we collected mosquito species across Thailand to gather genetic information, and evaluated the DNA barcoding efficacy for mosquito species identification in Thailand. A total of 310 mosquito samples, representing 73 mosquito species, were amplified using mitochondrial cytochrome c oxidase subunit I (COI) primers. The average maximum intraspecific genetic variation of the 73 mosquito species was 1% ranged from 0–5.7%. While, average minimum interspecific genetic variation (the distance to the nearest neighbour) of the 73 mosquito species was 7% ranged from 0.3–12.9%. The identification of success rates based on the “Best Match,” “Best Close Match,” and “All Species Barcodes” methods were 97.7%, 91.6%, and 81%, respectively. Phylogenetic analyses of Anopheles COI sequences demonstrated a clear separation between almost all species (except for those between An. baimaii and An. dirus), with high bootstrap support values (97%–99%). Furthermore, phylogenetic analyses revealed potential sibling species of An. annularis, An. tessellatus, and An. subpictus in Thailand. Our results indicated that DNA barcoding is an effective molecular approach for the accurate identification of mosquitoes in Thailand.

Description of the mitogenome and phylogeny of Aedes spp. (Diptera: Culicidae) from the Amazon region

The genus Aedes (Diptera: Culicidae) includes species of great epidemiological relevance, particularly involved in transmission cycles of leading arboviruses in the Brazilian Amazon region, such as the Zika virus (ZIKV), Dengue virus (DENV), Yellow fever virus (YFV), and Chikungunya virus (CHIKV). We report here the first putatively complete sequencing of the mitochondrial genomes of Brazilian populations of the species Aedes albopictus, Aedes scapularis and Aedes serratus. The sequences obtained showed an average length of 14,947 bp, comprising 37 functional subunits, typical in animal mitochondria (13 PCGs, 22 tRNA, and 2 rRNA). The phylogeny reconstructed by Maximum likelihood method, based on the concatenated sequences of all 13 PCGs produced at least two non-directly related groupings, composed of representatives of the subgenus Ochlerotatus and Stegomyia of the genus Aedes. The data and information produced here may be useful for future taxonomic and evolutionary studies of the genus Aedes, as well as the Culicidae family.

DNA Barcoding of Potential Mosquito Disease Vectors (Diptera, Culicidae) in Jazan Region, Saudi Arabia

The conventional morphological characterization of mosquito species remains heavily used for species identification in Jazan, Saudi Arabia. It requires substantial expertise and time, as well as having difficulty in confirming identity of morphologically similar species. Therefore, to establish a reliable and accurate identification system that can be applied to understanding spatial distribution of local mosquito species from the Jazan region, DNA barcoding was explored as an integrated tool for mosquito species identification. In this study, 44 adult mosquito specimens were analyzed, which contain 16 species belong to three genera of potential mosquito disease vectors (Aedes, Anopheles, and Culex). The specimens were collected from the Jazan region located in southwest Saudi Arabia. These included old and preserved mosquito voucher specimens. In addition, we assessed the genetic distance based on the generated mitochondrial partial COI DNA barcodes to detect cryptic diversity across these taxa. Nine mosquito species belonging to three genera were successfully barcoded and submitted to GenBank, namely: Aedes aegypti, Aedes caspius, Aedes vexans, Aedes vittatus, Anopheles arabiensis, Culex pipiens, Culex quinquefasciatus, Culex sitiens, and Culex tritaeniorhynchus. Of these nine species, Aedes vexans, Aedes vittatus, Culex sitiens, and Culex tritaeniorhynchus were registered in GenBank for the first time from Saudi Arabia. The DNA barcodes generated a 100% match to known barcodes of these mosquito species, that also matched with the morphological identification. Ae. vexans was found to be either a case of cryptic species (subspecies) or a new species from the region. However, more research has to be conducted to prove the latter. This study directly contributes to the development of a molecular reference library of mosquito species from the Jazan region and Saudi Arabia. The library is essential for confirmation of species in support of existing mosquito surveillance and control programmes.

Myzorhynchus series of Anopheles mosquitoes as potential vectors of Plasmodium bubalis in Thailand

Ungulate malaria parasites and their vectors are among the least studied when compared to other medically important species. As a result, a thorough understanding of ungulate malaria parasites, hosts, and mosquito vectors has been lacking, necessitating additional research efforts. This study aimed to identify the vector(s) of Plasmodium bubalis. A total of 187 female mosquitoes (133 Anopheles spp., 24 Culex spp., 24 Aedes spp., and 6 Mansonia spp. collected from a buffalo farm in Thailand where concurrently collected water buffalo samples were examined and we found only Anopheles spp. samples were P. bubalis positive. Molecular identification of anopheline mosquito species was conducted by sequencing of the PCR products targeting cytochrome c oxidase subunit 1 (cox1), cytochrome c oxidase subunit 2 (cox2), and internal transcribed spacer 2 (ITS2) markers. We observed 5 distinct groups of anopheline mosquitoes: Barbirostris, Hyrcanus, Ludlowae, Funestus, and Jamesii groups. The Barbirostris group (Anopheles wejchoochotei or Anopheles campestris) and the Hyrcanus group (Anopheles peditaeniatus) were positive for P. bubalis. Thus, for the first time, our study implicated these anopheline mosquito species as probable vectors of P. bubalis in Thailand.

Molecular phylogeny of Anopheles nivipes based on mtDNA-COII and mosquito diversity in Cambodia-Laos border

BACKGROUND

Few studies have been conducted to investigate the distribution of mosquito vectors and the population structure of secondary vectors in the border region of Cambodia-Laos. The aim of this work was to study the mosquito diversity and molecular phylogeny of secondary vectors, i.e., Anopheles nivipes in this area.

METHODS

1440 adult mosquitoes were trapped in the Cambodia-Laos border. mtDNA-COII were amplified and sequenced from 53 An. nivipes DNA samples. Together with COII sequences deposited in GenBank, a total of 86 COII sequences were used for examining population variations, genetic differentiation, spatial population structure, population expansion, and gene flow patterns.

RESULTS

The adult mosquitoes were classified into 5 genera and 27 species in this border region. The predominant genera were Culex (60.07%, 865/1440) and Anopheles (31.25%, 450/1440), and the major Anopheles species were An. nivipes (73.56%, 331/450) and Anopheles maculatus (14.22%, 64/450). Based on sequences analysis of COII, a high level of genetic differentiation was reported in two Northwest India (Cheema and Bathinda, Punjab) and Cambodia-Laos (Siem Pang, Stung treng) populations (F_ST = 0.97824, 0.97343, P < 0.05), as well as lower gene flow (Nm = 0.01112, 0.01365) in the An. nivipes populations. Phylogenetic analysis and SAMOVA revealed a gene barrier restricting gene flow among three An. nivipes populations. Mantel test suggested a significant correlation between geography and gene distance in all An. nivipes populations (Z = 44,983.1865, r = 0.5575, P = 0.0070). Neutrality test and Mismatch distribution revealed a recent population expansion of An. nivipes in the Cambodia-Laos population.

CONCLUSIONS

Anopheles nivipes was one of the major Anopheles species in the Cambodia-Laos border. Based on sequences analysis of COII, a genetic barrier between Cambodia-Laos and two Indian populations was found, and a recent population expanding or selecting of An. nivipes occurred in the Cambodia-Laos population, suggesting that COII might be an effective marker for describing the molecular phylogeny of An. nivipes. Further investigation and continuous surveillance of An. nivipes are warranted in this region.

Integrated Approaches for the Identification of Mosquitoes (Diptera: Culicidae) from the Volcanoes of Central America Physiographic Subprovince of the State of Chiapas, Mexico

Nowadays, there is a lack of information on the mosquito's fauna and DNA barcoding sequence reference library from many areas in Mexico, including the Volcanoes of Central America physiographic subprovince in the state of Chiapas. Consequently, a survey was undertaken to delineate the mosquito (Diptera: Culicidae) fauna in this region across different seasons using different collecting techniques. All species were identified by morphology and DNA barcoding, and their ecological features were also defined. In total, 62 taxa were morphologically examined, 60 of these were successfully identified based on morphological characteristics, but two were unable to be identified at the species level. The genera Aedes, Anopheles, Culex, and Wyeomyia are the most diverse among mosquito genera collected and include several species of medical and veterinary importance. Ecological characteristics of the immature habitats indicated that they were grouped into four categories namely, (1) large water bodies at ground level, (2) small and shady phytotelmata (e.g., tree holes and bamboo internodes), (3) large phytotelmata (e.g., plant leaves and axis bromeliad), and (4) artificial containers. The cytochrome c oxidase subunit I (COI) DNA barcoding sequences successfully separated the majority of these species, although specific species showed >2% intraspecific genetic divergences.

Improving the DNA Barcode Library of Mosquito Species With New Identifications and Discoveries in North-Central Argentina

The family Culicidae is represented by 244 species in Argentina, many of them with epidemiological importance. DNA barcodes are effective tools for identifying mosquito species, for knowing genetic variability, and for establishing phylogenetic relationships. This work aims to explore mosquito diversity employing different species delimitation approaches and to establish formally a DNA barcode library for the Argentinian mosquito fauna. Barcode fragments of 80 specimens of Argentinian mosquitoes of 28 species of the genera Aedeomyia Theobald (Diptera: Culicidae), Anopheles Meigen (Diptera: Culicidae), Coquillettidia Dyar (Diptera: Culicidae), Culex L. (Diptera: Culicidae), Haemagogus Williston (Diptera: Culicidae), Mansonia Blanchard (Diptera: Culicidae), Nyssorhynchus Blanchard (Diptera: Culicidae), Ochlerotatus Lynch-Arribálzaga (Diptera: Culicidae), Psorophora Robinneau-Desvoidy (Diptera: Culicidae) and Uranotaenia Lynch-Arribálzaga (Diptera: Culicidae) were sequenced. Another 82 sequences were obtained from public databases to establish the phylogenetic relationships using Maximum Likelihood and Bayesian Inference, and the species boundaries based on three approaches (ABGD, GMYC, and mPTP). Sixteen of the 28 species sequenced were recovered as monophyletic, of which 12 were also recognized as molecular operational taxonomic units according to the three methodologies. The disparity between morphology and barcode-based identifications could be explained by synonymy, species complexes occurrence, hybridization, incomplete lineage sorting, or the effect of the geographical scale of sampling. Twenty of the 28 sequenced species are new barcodes for Argentina and 11 are the first for science. This increases from 31 to 52 (12.7 to 21.31%) and from six to 10 (28.57 to 47.62%) the number of species and genera, respectively, with barcode sequences in Argentina. New species records are provided.

Pragmatic Applications and Universality of DNA Barcoding for Substantial Organisms at Species Level: A Review to Explore a Way Forward

DNA barcodes are regarded as hereditary succession codes that serve as a recognition marker to address several queries relating to the identification, classification, community ecology, and evolution of certain functional traits in organisms. The mitochondrial cytochrome c oxidase 1 (CO1) gene as a DNA barcode is highly efficient for discriminating vertebrate and invertebrate animal species. Similarly, different specific markers are used for other organisms, including ribulose bisphosphate carboxylase (rbcL), maturase kinase (matK), transfer RNA-H and photosystem II D1-ApbsArabidopsis thaliana (trnH-psbA), and internal transcribed spacer (ITS) for plant species; 16S ribosomal RNA (16S rRNA), elongation factor Tu gene (Tuf gene), and chaperonin for bacterial strains; and nuclear ITS for fungal strains. Nevertheless, the taxon coverage of reference sequences is far from complete for genus or species-level identification. Applying the next-generation sequencing approach to the parallel acquisition of DNA barcode sequences could greatly expand the potential for library preparation or accurate identification in biodiversity research. Overall, this review articulates on the DNA barcoding technology as applied to different organisms, its universality, applicability, and innovative approach to handling DNA-based species identification.

First Description of the Mitogenome and Phylogeny of Culicinae Species from the Amazon Region

The Culicidae family is distributed worldwide and comprises about 3587 species subdivided into the subfamilies Anophelinae and Culicinae. This is the first description of complete mitochondrial DNA sequences from Aedes fluviatilis, Aedeomyia squamipennis, Coquillettidia nigricans, Psorophora albipes, and Psorophora ferox. The mitogenomes showed an average length of 15,046 pb and 78.02% AT content, comprising 37 functional subunits (13 protein coding genes, 22 tRNAs, and two rRNAs). The most common start codons were ATT/ATG, and TAA was the stop codon for all PCGs. The tRNAs had the typical leaf clover structure, except tRNA^Ser1. Phylogeny was inferred by analyzing the 13 PCGs concatenated nucleotide sequences of 48 mitogenomes. Maximum likelihood and Bayesian inference analysis placed Ps. albipes and Ps. ferox in the Janthinosoma group, like the accepted classification of Psorophora genus. Ae. fluviatilis was placed in the Aedini tribe, but was revealed to be more related to the Haemagogus genus, a result that may have been hampered by the poor sampling of Aedes sequences. Cq. nigricans clustered with Cq. chrysonotum, both related to Mansonia. Ae. squamipennis was placed as the most external lineage of the Culicinae subfamily. The yielded topology supports the concept of monophyly of all groups and ratifies the current taxonomic classification.

Multiple Novel Clades of Anopheline Mosquitoes Caught Outdoors in Northern Zambia

Residual vector populations that do not come in contact with the most frequently utilized indoor-directed interventions present major challenges to global malaria eradication. Many of these residual populations are mosquito species about which little is known. As part of a study to assess the threat of outdoor exposure to malaria mosquitoes within the Southern and Central Africa International Centers of Excellence for Malaria Research, foraging female anophelines were collected outside households in Nchelenge District, northern Zambia. These anophelines proved to be more diverse than had previously been reported in the area. In order to further characterize the anopheline species, sequencing and phylogenetic approaches were utilized. Anopheline mosquitoes were collected from outdoor light traps, morphologically identified, and sent to Johns Hopkins Bloomberg School of Public Health for sequencing. Sanger sequencing from 115 field-derived samples yielded mitochondrial COI sequences, which were aligned with a homologous 488 bp gene segment from known anophelines (n = 140) retrieved from NCBI. Nuclear ITS2 sequences (n = 57) for at least one individual from each unique COI clade were generated and compared against NCBI’s nucleotide BLAST database to provide additional evidence for taxonomical identity and structure. Molecular and morphological data were combined for assignment of species or higher taxonomy. Twelve phylogenetic groups were characterized from the COI and ITS2 sequence data, including the primary vector species Anopheles funestus s.s. and An. gambiae s.s. An unexpectedly large proportion of the field collections were identified as An. coustani and An. sp. 6. Six phylogenetic groups remain unidentified to species-level. Outdoor collections of anopheline mosquitoes in areas frequented by people in Nchelenge, northern Zambia, proved to be extremely diverse. Morphological misidentification and underrepresentation of some anopheline species in sequence databases confound efforts to confirm identity of potential malaria vector species. The large number of unidentified anophelines could compromise the malaria vector surveillance and malaria control efforts not only in northern Zambia but other places where surveillance and control are focused on indoor-foraging and resting anophelines. Therefore, it is critical to continue development of methodologies that allow better identification of these populations and revisiting and cleaning current genomic databases.

DNA Barcoding of Mosquitoes from the Pantanos de Centla Biosphere Reserve, Southeastern Mexico

Accurate identification of mosquito species is essential to support programs that involve the study of distribution and mosquito control. Numerous mosquito species are difficult to identify based only on morphological characteristics, due to the morphological similarities in different life stages and large numbers of some species that are members of morphologically similar species complexes. In the present study, the mosquitoes collected in the Pantanos de Centla Biosphere Reserve, southeastern Mexico, were evaluated using a combination of morphological and molecular approaches (mitochondrial cytochrome c oxidase subunit I [COI] DNA barcode). A total of 1,576 specimens of 10 genera and 35 species, mostly adult stages, were collected. A total of 225 COI DNA barcode sequences were analyzed; most species formed well-supported groups in the neighbor joining, maximum likelihood, and Bayesian inference trees. The intraspecific Kimura 2-parameter (K2P) genetic distance averaged 1.52%. An intraspecific K2P distance of 6.20% was observed in Anopheles crucians s.l., while a deep split was identified in Culex erraticus and Cx. conspirator. This study showed that COI DNA barcodes offer a reliable approach to support mosquito species identification in Mexico.

Genetic diversity and population structure of Phlebotomus argentipes: Vector of Leishmania donovani in Sri Lanka

Phlebotomus argentipes is the vector of Leishmania donovani which causes the disease leishmaniasis, a neglected tropical disease and a growing health problem in Sri Lanka. A proper understanding of the population genetic structure of sand fly vectors is considered important prior to planning and implementation of a successful vector control program. Thus, the present study was conducted to determine the population genetic structure of sand fly vectors in Sri Lanka. Two mitochondrial genes namely Cytochrome c oxidase subunit 1 (Cox 1) and Cytochrome b (Cytb), and the internal transcribed spacer 2 (ITS2) region from the nuclear ribosomal DNA were used for molecular characterization. Analyses included maximum likelihood method, network analysis and DNA polymorphisms. The outcome revealed unique sequences of all genomic regions studied except the cox 1 gene had a relationship with sand flies isolated previously from Sri Lanka, India and Israel and cytb gene of 4 sand flies that aligned with those isolated earlier from Sri Lanka and 3 from Madagascar. Furthermore, cox 1 gene and ITS 2 region analyses based on F_ST values indicated a possible gene flow between the study sites whereas cytb gene analysis favoured the existence of genetically distinct populations of P. argentipes in each of the study sites. Poor population differentiation of P. argentipes, a possible consequence of a gene flow, is indeed of concern due to the risk imposed by promoting the spread of functionally important phenotypes such as insecticide resistance across the country, making future vector control efforts challenging.

A DNA barcode reference library of Portuguese mosquitoes

Mosquitoes are important biological vectors of pathogens and species identification in all life stages is the first step for effective monitoring and control of mosquito-borne diseases. Molecular methods for species identification have been developed over the last years to overcome the limitations of the taxonomic identification based on morphology. DNA barcoding, using a fragment of the mitochondrial cytochrome oxidase I (COI) gene, can be used for species identification but a reliable and comprehensive reference database of verified sequences is required. In this study, we aimed to generate a DNA barcode reference library for the identification of mosquito species from Portuguese mosquito fauna, including most relevant vector species. Mosquitoes captured under the National Vector Surveillance Program (REVIVE) were processed for DNA extraction, COI gene fragment amplification and sequencing. Nighty-eight barcode sequences were obtained, representing 26 species and 6 genera. Sequences were submitted to GenBank and BOLD and were used for validation of phenetic classification. Barcode Index Number (BIN) assignment and Automatic Barcode Gap Discovery (ABGD) were used and clustered COI sequences into twenty-five molecular operational taxonomic units (MOTUs). This is the first comprehensive study that combines morphological and molecular identification of most mosquito species present in Portugal aiming to offer a reliable framework for mosquito species identification.

Mosquito Identification From Bulk Samples Using DNA Metabarcoding: a Protocol to Support Mosquito-Borne Disease Surveillance in Canada

Approximately 80 species of mosquitoes (Diptera: Culicidae) have been documented in Canada. Exotic species such as Aedes albopictus (Skuse) (Diptera: Culicidae) are becoming established. Recently occurring endemic mosquito-borne diseases (MBD) in Canada including West-Nile virus (WNV) and Eastern Equine Encephalitis (EEE) are having significant public health impacts. Here we explore the use of DNA metabarcoding to identify mosquitoes from CDC light-trap collections from two locations in eastern Canada. Two primer pairs (BF2-BR2 and F230) were used to amplify regions of the cytochrome c oxidase subunit I (CO1) gene. High throughput sequencing was conducted using an Illumina MiSeq platform and GenBank-based species identification was applied using a QIIME 1.9 bioinformatics pipeline. From a site in southeastern Ontario, Canada, 26 CDC light trap collections of 72 to >300 individual mosquitoes were used to explore the capacity of DNA metabarcoding to identify and quantify captured mosquitoes. The DNA metabarcoding method identified 33 species overall while 24 species were identified by key. Using replicates from each trap, the dried biomass needed to identify the majority of species was determined to be 76 mg (equivalent to approximately 72 mosquitoes), and at least two replicates from the dried biomass would be needed to reliably detect the majority of species in collections of 144-215 mosquitoes and three replicates would be advised for collections with >215 mosquitoes. This study supports the use of DNA metabarcoding as a mosquito surveillance tool in Canada which can help identify the emergence of new mosquito-borne disease potential threats.

Newly Emerging Pest in China, Rhynchaenusmaculosus (Coleoptera: Curculionidae): Morphology and Molecular Identification with DNA Barcoding

The oak flea weevil, Rhynchaenusmaculosus Yang et Zhang 1991, is a newly emerging pest that severely damages oak (genus Quercus) in China. The first R. maculosus outbreak occurred in 2020 and caused spectacular damage to all oak forests in Jilin province, northeast China. The lack of key morphological characters complicates the identification of this native pest, especially in larva and pupa stages. This is problematic because quick and accurate species identification is crucial for early monitoring and intervention during outbreaks. Here, we provided the first detailed morphological description of R. maculosus at four life stages. Additionally, we used DNA barcodes from larva and pupa specimens collected from three remote locations for molecular identification. The average pairwise divergence of all sequences in this study was 0.51%, well below the 2% to 3% (K-2-parameter) threshold set for one species. All sample sequences matched the R. maculosus morphospecies (KX657706.1 and KX657707.1), with 99.23% to 100% (sequence identity, E value: 0.00) matching success. The tree based on barcodes placed the specimens into the Rhynchaenus group, and the phylogenetic relationship between 62 sequences (30 samples and 32 from GeneBank) had high congruence with the morphospecies taxa. The traditional DNA barcodes were successfully transformed into quick response codes with larger coding capacity for information storage. The results showed that DNA barcoding is reliable for R. maculosus identification. The integration of molecular and morphology-based methods contributes to accurate species identification of this newly emerging oak pest.

Updating the bionomy and geographical distribution of Anopheles (Nyssorhynchus) albitarsis F: A vector of malaria parasites in northern South America

Anopheles albitarsis F is a putative species belonging to the Albitarsis Complex, recognized by rDNA, mtDNA, partial white gene, and microsatellites sequences. It has been reported from the island of Trinidad, Venezuela and Colombia, and incriminated as a vector of malaria parasites in the latter. This study examined mitochondrially encoded cytochrome c oxidase I (MT-CO1) sequences of An. albitarsis F from malaria-endemic areas in Colombia and Venezuela to understand its relations with other members of the Complex, revised and update the geographical distribution and bionomics of An. albitarsis F and explore hypotheses to explain its phylogenetic relationships and geographical expansion. Forty-five MT-CO1 sequences obtained in this study were analyzed to estimate genetic diversity and possible evolutionary relationships. Sequences generated 37 haplotypes clustered in a group where the genetic divergence of Venezuelan populations did not exceed 1.6% with respect to Colombian samples. Anopheles albitarsis F (π = 0.013) represented the most recent cluster located closer to An. albitarsis I (π = 0.009). Barcode gap was detected according to Albitarsis Complex lineages previously reported (threshold 0.014–0.021). Anopheles albitarsis F has a wide distribution in northern South America and might play an important role in the transmission dynamics of malaria due to its high expansion capacity. Future studies are required to establish the southern distribution of An. albitarsis F in Venezuela, and its occurrence in Guyana and Ecuador.

The efficacy of DNA barcoding in the classification, genetic differentiation, and biodiversity assessment of benthic macroinvertebrates

Macroinvertebrates have been recognized as key ecological indicators of aquatic environment and are the most commonly used approaches for water quality assessment. However, species identification of macroinvertebrates (especially of aquatic insects) proves to be very difficult due to the lack of taxonomic expertise in some regions and can become time-consuming. In this study, we evaluated the feasibility of DNA barcoding for the classification of benthic macroinvertebrates and investigated the genetic differentiation in seven orders (Insecta: Ephemeroptera, Plecoptera, Trichoptera, Diptera, Hemiptera, Coleoptera, and Odonata) from four large transboundary rivers of northwest China and further explored its potential application to biodiversity assessment. A total of 1,144 COI sequences, belonging to 176 species, 112 genera, and 53 families were obtained and analyzed. The barcoding gap analysis showed that COI gene fragment yielded significant intra- and interspecific divergences and obvious barcoding gaps. NJ phylogenetic trees showed that all species group into monophyletic species clusters whether from the same population or not, except two species (Polypedilum. laetum and Polypedilum. bullum). The distance-based (ABGD) and tree-based (PTP and MPTP) methods were utilized for grouping specimens into Operational Taxonomic Units (OTUs) and delimiting species. The ABGD, PTP, and MPTP analysis were divided into 177 (p = .0599), 197, and 195 OTUs, respectively. The BIN analysis generated 186 different BINs. Overall, our study showed that DNA barcoding offers an effective framework for macroinvertebrate species identification and sheds new light on the biodiversity assessment of local macroinvertebrates. Also, the construction of DNA barcode reference library of benthic macroinvertebrates in Eurasian transboundary rivers provides a solid backup for bioassessment studies of freshwater habitats using modern high-throughput technologies in the near future.

Wing Morphometrics of Aedes Mosquitoes from North-Eastern France

BACKGROUND

In the context of the increasing circulation of arboviruses, a simple, fast and reliable identification method for mosquitoes is needed. Geometric morphometrics have proven useful for mosquito classification and have been used around the world on known vectors such as Aedes albopictus. Morphometrics applied on French indigenous mosquitoes would prove useful in the case of autochthonous outbreaks of arboviral diseases.

METHODS

We applied geometric morphometric analysis on six indigenous and invasive species of the Aedes genus in order to evaluate its efficiency for mosquito classification.

RESULTS

Six species of Aedes mosquitoes (Ae. albopictus, Ae. cantans, Ae. cinereus, Ae. sticticus, Ae. japonicus and Ae. rusticus) were successfully differentiated with Canonical Variate Analysis of the Procrustes dataset of superimposed coordinates of 18 wing landmarks.

CONCLUSIONS

Geometric morphometrics are effective tools for the rapid, inexpensive and reliable classification of at least six species of the Aedes genus in France.

Molecular Characterization of Mosquito Diversity in the Balearic Islands

Several outbreaks of mosquito-borne diseases have taken place in Europe in recent years. In Spain, both active and passive surveillance have demonstrated that dengue and West Nile viruses are currently circulating, and seven autochthonous dengue cases have been reported in the last 2 yr. The effectiveness of vector control programs largely depends on the accuracy of the taxonomic identification of the species. However, in Spain, identification almost completely relies on the use of morphological keys to characterize the mosquito fauna. This study investigates the congruence between molecular and morphological species boundaries in 13 Spanish mosquito taxa. The Cytochrome c oxidase subunit I (COI) gene region was sequenced from 60 adult specimens collected in Mallorca, plus several representatives from other Spanish regions for comparative purposes. Phylogenetic relationships were established using Bayesian and maximum-likelihood approaches. Using three species delimitation algorithms (ABGD, mPTP, and GMYC), we found strong evidence for cryptic speciation within Anopheles algeriensis Theobald, a widespread mosquito in the Mediterranean basin. We also delimited the Mallorcan rock pool mosquito Aedes mariae (Sergent & Sergent), from mainland European populations. Finally, we found difficulties in the use of wing characters in species keys to distinguish Culiseta annulata (Schrk) from Culiseta subochrea (Edwards). Given that these species are vectors of pathogens of medical relevance and have veterinary importance, their accurate taxonomic identification is essential in European vector surveillance programs.

Phylogenetic Network of Mitochondrial COI Gene Sequences Distinguishes 10 Taxa Within the Neotropical Albitarsis Group (Diptera: Culicidae), Confirming the Separate Species Status of Anopheles albitarsis H (Diptera: Culicidae) and Revealing a Novel Lineage, Anopheles albitarsis J

The Neotropical Albitarsis Group is a complex assemblage of essentially isomorphic species which currently comprises eight recognized species-five formally described (Anopheles albitarsis Lynch-Arribalzaga, An. deaneorum Rosa-Freitas, An. janconnae Wilkerson and Sallum, An. marajoara Galvao and Damasceno, An. oryzalimnetes Wilkerson and Motoki) and three molecularly assigned (An. albitarsis F, G & I)-and one mitochondrial lineage (An. albitarsis H). To further explore species recognition within this important group, 658 base pairs of the mitochondrial DNA cytochrome oxidase subunit I (COI) were analyzed from 988 specimens from South America. We conducted statistical parsimony network analysis, generated estimates of haplotype, nucleotide, genetic differentiation, divergence time, and tested the effect of isolation by distance (IBD). Ten clusters were identified, which confirmed the validity of the eight previously determined species, and confirmed the specific status of the previous mitochondrial lineage An. albitarsis H. High levels of diversity were highlighted in two samples from Pará (= An. albitarsis J), which needs further exploration through additional sampling, but which may indicate another cryptic species. The highest intra-specific nucleotide diversity was observed in An. deaneorum, and the lowest in An. marajoara. Significant correlation between genetic and geographical distance was observed only in An. oryzalimnetes and An. albitarsis F. Divergence time within the Albitarsis Group was estimated at 0.58-2.25 Mya, during the Pleistocene. The COI barcode region was considered an effective marker for species recognition within the Albitarsis Group and a network approach was an analytical method to discriminate among species of this group.

An Integrated Molecular Approach to Untangling Host-Vector-Pathogen Interactions in Mosquitoes (Diptera: Culicidae) From Sylvan Communities in Mexico

There are ~240 species of Culicidae in Mexico, of which some are vectors of arthropod-borne viruses such as Zika virus, dengue virus, chikungunya virus, and West Nile virus. Thus, the identification of mosquito feeding preferences is paramount to understanding of vector–host–pathogen interactions that, in turn, can aid the control of disease outbreaks. Typically, DNA and RNA are extracted separately for animal (insects and blood meal hosts) and viral identification, but this study demonstrates that multiple organisms can be analyzed from a single RNA extract. For the first time, residual DNA present in standard RNA extracts was analyzed by DNA barcoding in concert with Sanger and next-generation sequencing (NGS) to identify both the mosquito species and the source of their meals in blood-fed females caught in seven sylvan communities in Chiapas State, Mexico. While mosquito molecular identification involved standard barcoding methods, the sensitivity of blood meal identification was maximized by employing short primers with NGS. In total, we collected 1,634 specimens belonging to 14 genera, 25 subgenera, and 61 morphospecies of mosquitoes. Of these, four species were new records for Mexico (Aedes guatemala, Ae. insolitus, Limatus asulleptus, Trichoprosopon pallidiventer), and nine were new records for Chiapas State. DNA barcode sequences for >300 bp of the COI gene were obtained from 291 specimens, whereas 130 bp sequences were recovered from another 179 specimens. High intraspecific divergence values (>2%) suggesting cryptic species complexes were observed in nine taxa: Anopheles eiseni (5.39%), An. pseudopunctipennis (2.79%), Ae. podographicus (4.05%), Culex eastor (4.88%), Cx. erraticus (2.28%), Toxorhynchites haemorrhoidalis (4.30%), Tr. pallidiventer (4.95%), Wyeomyia adelpha/Wy. guatemala (7.30%), and Wy. pseudopecten (4.04%). The study increased the number of mosquito species known from 128 species to 138 species for Chiapas State, and 239 for Mexico as a whole. Blood meal analysis showed that Aedes angustivittatus fed on ducks and chicken, whereas Psorophora albipes fed on humans. Culex quinquefasciatus fed on diverse hosts including chicken, human, turkey, and Mexican grackle. No arbovirus RNA was detected by reverse transcriptase–polymerase chain reaction in the surveyed specimens. This study demonstrated, for the first time, that residual DNA present in RNA blood meal extracts can be used to identify host vectors, highlighting the important role of molecular approaches in both vector identification and revealing host–vector–pathogen interactions.

Molecular species delimitation reveals high diversity in the mosquito Anopheles tessellatus Theobald, 1901 (Diptera, Culicidae) across its range

Anopheles tessellatus is a potentially important vector found across South, East and Southeast Asia. While it was formerly considered a formidable vector of human Plasmodium and filarial parasites in the Maldives, and of lesser importance as a vector of human Plasmodium in Sri Lanka and parts of Indonesia, it is currently of little or unknown health importance in many other parts of its range. This study describes the genetic diversity and evolutionary relationships among An. tessellatus populations in nine Asian countries at the COI gene using maximum-likelihood and Bayesian phylogenetic inference tree and cluster-based species delimitation approaches. These analyses reveal exceptional levels of genetic diversity in An. tessellatus populations across its known range, and identify up to six putative species in the newly determined Tessellatus Complex. The existence of such cryptic diversity has potentially important consequences for vector management and disease control. Differences in the ecologies and life histories among these species may have considerable impact on vectorial capacity and may go some way towards explaining why An. tessellatus s.l. has such varying degrees of public health importance across its range.

Foretelling the Phenotype of a Genomic Sequence

Estimating phenotypic features (physical and biochemical traits) in a biological organism from their genomic sequence alone and/or environmental conditions has major applications in anthropological paleontology and criminal forensics, for example. To what extent do genomic sequences generally and causally determine phenotypic features of organisms, environmental conditions aside? We present results of two studies, one in blackfly (Insecta:Diptera:Simuliidae) larvae in two species (Simulium ignescens and S. tunja) with four phenotypic features, including the area and spot pattern of the cephalic apotome (in the form of a latin cross on the dorsal side of the head), the postgenal cleft (area under the head on the ventral side) and general body color in larva specimens; the second in strains of Arabidopsis thaliana. They establish that a substantial component of these phenotypic features (over 75 percent) are at least logically inferable, if not causally determined, by genomic fragments alone, despite the fact that these phenotypic features are not 100 percent determined entirely by genetic traits. These results suggest that it is possible to infer the genetic contribution in the determination of specific phenotypic features of a biological organism, without recourse to the causal chain of metabolomics and proteomic events leading to them from genomic sequences.

Molecular and Morphological Identification of Dermanyssoid Mites (Parasitiformes: Mesostigmata: Dermanyssoidea) Causatives of a Parasitic Outbreak on Captive Snakes

A parasitic outbreak caused by dermanyssoid mites in a herpetarium of the Metropolitan area of the Valley of Mexico is revealed. This outbreak was caused by Hemilaelaps triangulus (Ewing), but a second mite species, Ophionyssus natricis (Gervais), was found in low abundance. The parasitic load is analyzed, and the morphological and molecular diagnostic characters to identify each of the two species involved are given. A barcode analysis is presented, and two more molecular markers are presented and analyzed. Hemilaelaps triangulus is recorded for the first time in Mexico, and this is the first record of massive infestation on captive snakes caused by ixodorhynchid mites, and DNA sequences of ixodorhynchid mites are publicly available for the first time.

Comparative Phylogeography and Integrative Taxonomy of Ochlerotatus caspius (Dipera: Culicidae) and Ochlerotatus dorsalis

Given that accurately identifying pathogen vectors is vital for designing efficient mosquito control programs based on the proper surveillance of the epidemiologically important species, it has been suggested the complementary use of independently evolving genes and morphometric traits as a reliable approach for the characterization and delimitation of related species. Hence, we examined the spatial distribution of COI mtDNA and ITS2 rDNA variation from the historical perspective of Ochlerotatus caspius (Pallas, 1771) and O. dorsalis (Meigen, 1830), while simultaneously testing the utility of the two markers in integrative species delimitation when combined with phenotypic character analyses of larvae and adults. Despite the striking difference in haplotype diversity (high in COI mtDNA, low in ITS2 rDNA), no evident phylogeographic structure was apparent in the Palearctic O. caspius. The Holarctic O. dorsalis species was subdivided into two highly distinctive COI mtDNA phylogroups which corresponded to the Nearctic and Palearctic regions. Strong support for the independence of the two allopatric evolutionary lineages suggested that geographical barrier and climatic changes during Pleistocene caused vicariance of the ancestral range. COI mtDNA reliably distinguished O. caspius and O. dorsalis, while ITS2 rDNA yet again lacked the proper resolution for solving this problem. An integrative approach based on the larval and adult morphological traits have varying taxonomic applications due to their differential diagnostic values. Thus, by the implementation of an integrative taxonomic approach, we successfully detected species borders between the two epidemiologically relevant species and uncovered the presence of cryptic diversity within O. dorsalis.

Identification of mosquitoes (Diptera: Culicidae) from Mexico State, Mexico using morphology and COI DNA barcoding

Mosquitoes are commonly identified to species level using morphological traits, but complementary methods for identification are often necessary when specimens are collected as immature stages, stored inadequately, or when delineation of species complexes is problematic. DNA-barcoding using the mitochondrial cytochrome c oxidase subunit 1 (COI) gene is one such tool used for the morphological identification of species. A comprehensive entomological survey of mosquito species in Mexico State identified by COI DNA barcoding and morphology is documented in this paper. Specimens were collected from all the physiographic provinces in Mexico State between 2017 and 2019. Overall, 2,218 specimens were collected from 157 localities representing both subfamilies Anophelinae and Culicinae. A species checklist that consists of 6 tribes, 10 genera, 20 subgenera, and 51 species, 35 of which are new records for Mexico State, is provided. Three hundred and forty-two COI sequences of 46 species were analysed. Mean intraspecific and interspecific distances ranged between 0% to 3.9% and from 1.2% to 25.3%, respectively. All species groups were supported by high bootstraps values in a Neighbour-Joining analysis, and new COI sequences were generated for eight species: Aedes chionotum Zavortink, Ae. vargasi Schick, Ae. gabriel Schick, Ae. guerrero Berlin, Ae. ramirezi Vargas and Downs, Haemagogus mesodentatus Komp and Kumm, Culex restrictor Dyar and Knab, and Uranotaenia geometrica Theobald. This study provides a detailed inventory of the Culicidae from Mexico State and discusses the utility of DNA barcoding as a complementary tool for accurate mosquito species identification in Mexico.

Modeling host-feeding preference and molecular systematics of mosquitoes in different ecological niches in Canada

Several mosquito-borne viruses (mobovirus) cause infections in Canada. Ecological data on mosquito species and host range in Canada remains elusive. The main aim of the current study is to determine the host range and molecular systematics of mosquito species in Canada. Mosquitoes were collected using BG-Sentinel traps and aspirators at 10 trapping sites in Canada during 2018 and 2019. Mosquitoes collected were identified via morphology and molecular techniques. Mosquito sequences were aligned by MUSCLE algorithm and evolutionary systematics were drawn using MEGA and SDT software. Moreover, the source of blood meals was identified using a DNA barcoding technique. A total of 5,708 female mosquitoes over 34 different taxa were collected. DNA barcodes and evolutionary tree analysis confirmed the identification of mosquito species in Canada. Of the total collected samples, 201 specimens were blood-fed female mosquitoes in 20 different taxa. Four mosquito species represented about half (51.47%) of all collected blood-fed specimens: Aede cinereus (39 specimens, 19.11%), Aedes triseriatus (23, 11.27%), Culex pipiens (22, 10.78%), and Anopheles punctipennis (21, 10.29%). The most common blood meal sources were humans (49 mosquito specimens, 24% of all blood-fed mosquito specimen), pigs (44, 21.5%), American red squirrels (28, 13.7%), white-tailed deers (28, 13.7%), and American crows (16, 7.8%). Here, we present the first analysis of the host-feeding preference of different mosquito species in Canada via molecular techniques. Our results on mosquito distribution and behavior will aid in the development of effective mitigation and control strategies to prevent or reduce human/animal health issues in regards to moboviruses.

DNA barcoding and species delimitation of butterflies (Lepidoptera) from Nigeria

Accurate identification of species is a prerequisite for successful biodiversity management and further genetic studies. Species identification techniques often require both morphological diagnostics and molecular tools, such as DNA barcoding, for correct identification. In particular, the use of the subunit I of the mitochondrial cytochrome c oxidase (COI) gene for DNA barcoding has proven useful in species identification for insects. However, to date, no studies have been carried out on the DNA barcoding of Nigerian butterflies. We evaluated the utility of DNA barcoding applied for the first time to 735 butterfly specimens from southern Nigeria. In total, 699 DNA barcodes, resulting in a record of 116 species belonging to 57 genera, were generated. Our study sample comprised 807 DNA barcodes based on sequences generated from our current study and 108 others retrieved from BOLD. Different molecular analyses, including genetic distance-based evaluation (Neighbor-Joining, Maximum Likelihood and Bayesian trees) and species delimitation tests (TaxonDNA, Automated Barcode Gap Discovery, General Mixed Yule-Coalescent, and Bayesian Poisson Tree Processes) were performed to accurately identify and delineate species. The genetic distance-based analyses resulted in 163 well-separated clusters consisting of 147 described and 16 unidentified species. Our findings indicate that about 90.20% of the butterfly species were explicitly discriminated using DNA barcodes. Also, our field collections reported the first country records of ten butterfly species-Acraea serena, Amauris cf. dannfelti, Aterica galena extensa, Axione tjoane rubescens, Charaxes galleyanus, Papilio lormieri lormeri, Pentila alba, Precis actia, Precis tugela, and Tagiades flesus. Further, DNA barcodes revealed a high mitochondrial intraspecific divergence of more than 3% in Bicyclus vulgaris vulgaris and Colotis evagore. Furthermore, our result revealed an overall high haplotype (gene) diversity (0.9764), suggesting that DNA barcoding can provide information at a population level for Nigerian butterflies. The present study confirms the efficiency of DNA barcoding for identifying butterflies from Nigeria. To gain a better understanding of regional variation in DNA barcodes of this biogeographically complex area, future work should expand the DNA barcode reference library to include all butterfly species from Nigeria as well as surrounding countries. Also, further studies, involving relevant genetic and eco-morphological datasets, are required to understand processes governing mitochondrial intraspecific divergences reported in some species complexes.

A simple PCR-based method for the rapid and accurate identification of spider mites (Tetranychidae) on cassava

The morphological identification of mites entails great challenges. Characteristics such as dorsal setae and aedeagus are widely used, but they show variations between populations, and the technique is time consuming and demands specialized taxonomic expertise that is difficult to access. A successful alternative has been to exploit a region of the mitochondrial cytochrome oxidase I (COI) gene to classify specimens to the species level. We analyzed the COI sequences of four mite species associated with cassava and classified them definitively by detailed morphological examinations. We then developed an identification kit based on the restriction fragment length polymorphism-polymerase chain reaction of subunit I of the COI gene focused on the three restriction enzymes AseI, MboII, and ApoI. This set of enzymes permitted the simple, accurate identification of Mononychellus caribbeanae, M. tanajoa, M. mcgregori, and Tetranychus urticae, rapidly and with few resources. This kit could be a vital tool for the surveillance and monitoring of mite pests in cassava crop protection programs in Africa, Asia, and Latin America.

Mitochondrial genome sequencing and phylogeny of Haemagogus albomaculatus, Haemagogus leucocelaenus, Haemagogus spegazzinii, and Haemagogus tropicalis (Diptera: Culicidae)

The genus Haemagogus (Diptera: Culicidae) comprises species of great epidemiological relevance, involved in transmission cycles of the Yellow fever virus and other arboviruses in South America. So far, only Haemagogus janthinomys has complete mitochondrial sequences available. Given the unavailability of information related to aspects of the evolutionary biology and molecular taxonomy of this genus, we report here, the first sequencing of the mitogenomes of Haemagogus albomaculatus, Haemagogus leucocelaenus, Haemagogus spegazzinii, and Haemagogus tropicalis. The mitogenomes showed an average length of 15,038 bp, average AT content of 79.3%, positive AT-skews, negative GC-skews, and comprised 37 functional subunits (13 PCGs, 22 tRNA, and 02 rRNA). The PCGs showed ATN as start codon, TAA as stop codon, and signs of purifying selection. The tRNAs had the typical leaf clover structure, except tRNA^Ser1. Phylogenetic analyzes of Bayesian inference and Maximum Likelihood, based on concatenated sequences from all 13 PCGs, produced identical topologies and strongly supported the monophyletic relationship between the Haemagogus and Conopostegus subgenera, and corroborated with the known taxonomic classification of the evaluated taxa, based on external morphological aspects. The information produced on the mitogenomes of the Haemagogus species evaluated here may be useful in carrying out future taxonomic and evolutionary studies of the genus.

New records and DNA barcoding of deer flies, Chrysops (Diptera: Tabanidae) in Thailand

Chrysops spp. or deer flies (Diptera: Tabanidae) are hematophagous flies of medical and veterinary importance and some species are important vectors of Trypanosoma evansi, the causative agent of surra in Thailand. However, data regarding deer fly species and their molecular identification are limited. Accurate species identification will indicate the appropriate control measures. In this study, an entomological survey of deer flies from different sites in Thailand between May 2018 and June 2019 were conducted. In addition, mitochondrial cytochrome oxidase subunit I (COI) barcoding region was used for species identification. A total of 82 females were collected and 6 species were identified. Of these, three species are new records for Thailand: C. designatus, C. fuscomarginalis and C. vanderwulpi bringing the species total found in Thailand to nine. The COI sequences revealed an intraspecific divergence of 0.0%-2.65% and an interspecific divergence of 7.03%-13.47%. Phylogenetic analysis showed that all deer fly species were clearly separated into distinct clusters according to morphologically identified species. These results indicated that COI barcodes were capable in discriminating between deer fly species on the basis of the barcoding gap and phylogenetic analysis. Therefore, DNA barcoding is a valuable tool for species identification of deer flies in Thailand.

Longer is Not Always Better: Optimizing Barcode Length for Large-Scale Species Discovery and Identification

New techniques for the species-level sorting of millions of specimens are needed in order to accelerate species discovery, determine how many species live on earth, and develop efficient biomonitoring techniques. These sorting methods should be reliable, scalable, and cost-effective, as well as being largely insensitive to low-quality genomic DNA, given that this is usually all that can be obtained from museum specimens. Mini-barcodes seem to satisfy these criteria, but it is unclear how well they perform for species-level sorting when compared with full-length barcodes. This is here tested based on 20 empirical data sets covering ca. 30,000 specimens (5500 species) and six clade-specific data sets from GenBank covering ca. 98,000 specimens ($>$20,000 species). All specimens in these data sets had full-length barcodes and had been sorted to species-level based on morphology. Mini-barcodes of different lengths and positions were obtained in silico from full-length barcodes using a sliding window approach (three windows: 100 bp, 200 bp, and 300 bp) and by excising nine mini-barcodes with established primers (length: 94-407 bp). We then tested whether barcode length and/or position reduces species-level congruence between morphospecies and molecular operational taxonomic units (mOTUs) that were obtained using three different species delimitation techniques (Poisson Tree Process, Automatic Barcode Gap Discovery, and Objective Clustering). Surprisingly, we find no significant differences in performance for both species- or specimen-level identification between full-length and mini-barcodes as long as they are of moderate length ($>$200 bp). Only very short mini-barcodes (<200 bp) perform poorly, especially when they are located near the 5$^\prime$ end of the Folmer region. The mean congruence between morphospecies and mOTUs was ca. 75% for barcodes $>$200 bp and the congruent mOTUs contain ca. 75% of all specimens. Most conflict is caused by ca. 10% of the specimens that can be identified and should be targeted for re-examination in order to efficiently resolve conflict. Our study suggests that large-scale species discovery, identification, and metabarcoding can utilize mini-barcodes without any demonstrable loss of information compared to full-length barcodes. [DNA barcoding; metabarcoding; mini-barcodes; species discovery.].

Mitochondrial metabolic genes provide phylogeographic relationships of global collections of Aedes aegypti (Diptera: Culicidae)

Phylogeographic relationships among global collections of the mosquito Aedes aegypti were evaluated using the mitochondrial Cytochrome C Oxidase 1 (CO1) and NADH dehydrogenase subunit 4 (ND4) genes including new sequences from Sri Lanka. Phylogeographic analysis estimated that Ae. aegypti arose as a species ~614 thousand years ago (kya) in the late Pleistocene. At 545 kya an “early” East African clade arose that continued to differentiate in East Africa, and eventually gave rise to three lineages one of which is distributed throughout all tropical and subtropical regions, a second that contains Southeast Asian/Sri Lankan mosquitoes and a third that contains mostly New World mosquitoes. West African collections were not represented in this early clade. The late clade continued to differentiate throughout Africa and gave rise to a lineage that spread globally. The most recent branches of the late clade are represented by South-East Asia and India/Pakistan collections. Analysis of migration rates suggests abundant gene flow between India/Pakistan and the rest of the world with the exception of Africa.

Molecular Identification and Geometric Morphometric Analysis of Haematobosca aberrans (Diptera: Muscidae)

The genus Haematobosca Bezzi, 1907 (Diptera: Muscidae) contains haematophagous flies of veterinary importance. A new fly species of this genus was recognised from northern Thailand based on morphological characters and described as Haematobosca aberrans Pont, Duvallet & Changbunjong, 2020. In the present study, the mitochondrial cytochrome c oxidase I (COI) gene was used to confirm the morphological identification of H. aberrans. In addition, landmark-based geometric morphometrics was used to determine sexual dimorphism. The molecular analysis was conducted with 10 COI sequences. The results showed that all sequences were 100% identical. The sequence was not highly similar to reference sequences from GenBank and did not match any identified species from Barcode of Life Data Systems (BOLD). Phylogenetic analysis clearly differentiated this species from other species within the subfamily Stomoxyinae. For geometric morphometric analysis, a total of 16 wing pictures were analysed using the landmark-based approach. The results showed significant differences in wing shape between males and females, with a cross-validated classification score of 100%. The allometric analysis showed that wing shape has no correlation with size. Therefore, the COI gene is effective in species identification of H. aberrans, and geometric morphometrics is also effective in determining sexual dimorphism.