Assessment of population genetic structure in the arbovirus vector midge, Culicoides brevitarsis (Diptera: Ceratopogonidae), using multi-locus DNA microsatellites

Genetic analyses of Culicoides peregrinus (Diptera: Ceratopogonidae) reveal population differentiation in southern India a hotspot of bluetongue disease outbreaks and lack of gene flow across its range

Among India's 7 potential vectors of bluetongue virus (BTV), Culicoides peregrinus has assumed importance due to the recent isolation of BTV-23, with several severe bluetongue disease (BTD) outbreaks in Southern India and high abundance throughout India. These factors warrant the elucidation of its genetic structure across various physiographic regions. Adult Culicoides were captured from 14 collection locations across different physiographic areas of India. The study used 2 genetic markers, ie mitochondrial cytochrome oxidase-1 (COX-1) and nuclear ribosomal protein-S2B intron region (EPIC). COX-1 sequences of C. peregrinus from the other areas of its distribution were also considered. A high intraspecific genetic distance, ie 2.46% (within India) and 3.5% (across its distribution range) indicates the possibility of cryptic species. This vector has a relatively uniform genetic structure with a slight difference in the southern Indian population. However, the populations from different countries are genetically distinct and exhibit high genetic distances among themselves. There is evidence of demographic expansion with India as the probable ancestral region with 3 distinct lineages: India-Bangladesh, Thailand-China-Malaysia, and Australia. STRUCTURE analysis using EPIC marker (ribosomal protein-S2B intron region) depicted the distribution of 2 alleles in India but no population substructure. Pairwise FST indicated the south Indian population to have low but significant genetic differentiation with other Indian populations. More fine scale studies to ascertain the genetic differentiation of the Indian population along with direct experimental validation of vector competence is imperative to confirm C. peregrinus as the causative vector of BTD outbreaks.

Molecular characterisation of common Culicoides biting midges (Diptera: Ceratopogonidae) in Ireland

BACKGROUND

Biting midges of the genus Culicoides (Diptera: Ceratopogonidae) act as vectors for several arboviruses, including bluetongue virus (BTV) and Schmallenberg virus (SBV), which affect livestock health and productivity. In Ireland, limited genetic data are available regarding the diversity of Culicoides species. This study represents the first attempt to characterise Culicoides in this region using molecular techniques.

METHODS

Adult Culicoides samples were captured using Onderstepoort Veterinary Institute (OVI) traps across six locations in Ireland. Subsequent molecular analyses involved polymerase chain reaction (PCR) and sequencing of the cytochrome oxidase subunit 1 (CO1) and the internal transcriber spacer (ITS) barcoding regions to obtain species identities. In addition, using both markers, we inferred the population genetic structure and potential colonisation pathways of Culicoides obsoletus sensu stricto (s. str.), the major vector species in Ireland.

RESULTS

DNA barcoding facilitated identification of 177 specimens. Eight common Culicoides species were identified through DNA barcoding of CO1 and ITS gene regions. The presence of putative vectors of bluetongue virus (BTV) and Schmallenberg virus (SBV) were also confirmed, including species in the subgenus Avaritia (C. obsoletus s. str., C. scoticus, C. chiopterus, and C. dewulfi) and subgenus Culicoides s. str. (C. pulicaris and C. punctatus). Phylogenetic analysis confirmed the relationship between these vector species and facilitated the placement of Culicoides spp. that could not be identified to species level through DNA barcoding. Haplotype network analysis of C. obsoletus showed that some haplotypes of these species are shared between Continental Europe, the UK, and Ireland, suggesting a possible incursion pathway for this vector.

CONCLUSIONS

DNA barcoding employing a combination of two barcodes, CO1 and ITS, proved effective in identifying Culicoides, especially species within the obsoletus complex, which are difficult to morphologically distinguish. Our findings also suggest that investigation of the population genetic structure of Culicoides spp. could be used to model the potential introduction routes of midge-borne pathogens into the country.

Whole-transcriptome analyses of ovine lung microvascular endothelial cells infected with bluetongue virus

Bluetongue virus (BTV) infection induces profound and intricate changes in the transcriptional profile of the host to facilitate its survival and replication. However, there have been no whole-transcriptome studies on ovine lung microvascular endothelial cells (OLMECs) infected with BTV. In this study, we comprehensively analysed the whole-transcriptome sequences of BTV-1 serotype-infected and mock-infected OLMECs and subsequently performed bioinformatics differential analysis. Our analysis revealed 1215 differentially expressed mRNA transcripts, 82 differentially expressed long noncoding RNAs (lncRNAs) transcripts, 63 differentially expressed microRNAs (miRNAs) transcripts, and 42 differentially expressed circular RNAs (circRNAs) transcripts. Annotation from Gene Ontology, enrichment from the Kyoto Encyclopedia of Genes and Genomes, and construction of endogenous competing RNA network analysis revealed that the differentially expressed RNAs primarily participated in viral sensing and signal transduction pathways, antiviral and immune responses, inflammation, and extracellular matrix (ECM)-related pathways. Furthermore, protein‒protein interaction network analysis revealed that BTV may regulate the conformation of ECM receptor proteins and change their biological activity through a series of complex mechanisms. Finally, on the basis of real-time fluorescence quantitative polymerase chain reaction results, the expression trends of the differentially expressed RNA were consistent with the whole-transcriptome sequencing data, such as downregulation of the expression of COL4A1, ITGA8, ITGB5, and TNC and upregulation of the expression of CXCL10, RNASEL, IRF3, IRF7, and IFIHI. This study provides a novel perspective for further investigations of the mechanism of the ECM in the BTV-host interactome and the pathogenesis of lung microvascular endothelial cells.

Immigration and seasonal bottlenecks: high inbreeding despite high genetic diversity in an oscillating population of Culicoides sonorensis (Diptera: Ceratopogonidae)

Most population genetic studies concern spatial genetic differentiation, but far fewer aim at analyzing the temporal genetic changes that occur within populations. Vector species, including mosquitoes and biting midges, are often characterized by oscillating adult population densities, which may affect their dispersal, selection, and genetic diversity over time. Here, we used a population of Culicoides sonorensis from a single site in California to investigate short-term (intra-annual) and long-term (inter-annual) temporal variation in genetic diversity over a 3 yr period. This biting midge species is the primary vector of several viruses affecting both wildlife and livestock, thus a better understanding of the population dynamics of this species can help inform epidemiological studies. We found no significant genetic differentiation between months or years, and no correlation between adult populations and the inbreeding coefficient (FIS). However, we show that repeated periods of low adult abundance during cooler winter months resulted in recurring bottleneck events. Interestingly, we also found a high number of private and rare alleles, which suggests both a large, stable population, as well as a constant influx of migrants from nearby populations. Overall, we showed that the high number of migrants maintains a high level of genetic diversity by introducing new alleles, while this increased diversity is counterbalanced by recurrent bottleneck events potentially purging unfit alleles each year. These results highlight the temporal influences on population structure and genetic diversity in C. sonorensis and provide insight into factors effecting genetic variation that may occur in other vector species with fluctuating populations.

Population Genetic Structure and Population History of the Biting Midge Culicoides mahasarakhamense (Diptera: Ceratopogonidae)

Biting midges of the genus Culicoides Latreille are significant pests and vectors of disease agents transmitted to humans and other animals. Understanding the genetic structure and diversity of these insects is crucial for effective control programs. This study examined the genetic diversity, genetic structure, and demographic history of Culicoides mahasarakhamense, a possible vector of avian haemosporidian parasites and Leishmania martiniquensis, in Thailand. The star-like shape of the median joining haplotype network, a unimodal mismatch distribution, and significant negative values for Tajima's D and Fu's F_S tests indicated that populations had undergone recent expansion. Population expansion time was estimated to be 2000-22,000 years ago. Population expansion may have been triggered by climatic amelioration from cold/dry to warm/humid conditions at the end of the last glaciations, resulting in the increased availability of host blood sources. Population pairwise F_ST revealed that most (87%) comparisons were not genetically different, most likely due to a shared recent history. The exception to the generally low level of genetic structuring is a population from the northern region that is genetically highly different from others. Population isolation in the past and the limitation of ongoing gene flows due to large geographic distance separation are possible explanations for genetic differentiation.

Species delimitation and mitonuclear discordance within a species complex of biting midges

The inability to distinguish between species can be a serious problem in groups responsible for pathogen transmission. Culicoides biting midges transmit many pathogenic agents infecting wildlife and livestock. In North America, the C. variipennis species complex contains three currently recognized species, only one of which is a known vector, but limited species-specific characters have hindered vector surveillance. Here, genomic data were used to investigate population structure and genetic differentiation within this species complex. Single nucleotide polymorphism data were generated for 206 individuals originating from 17 locations throughout the United States and Canada. Clustering analyses suggest the occurrence of two additional cryptic species within this complex. All five species were significantly differentiated in both sympatry and allopatry. Evidence of hybridization was detected in three different species pairings indicating incomplete reproductive isolation. Additionally, COI sequences were used to identify the hybrid parentage of these individuals, which illuminated discordance between the divergence of the mitochondrial and nuclear datasets.

High dispersal capacity of Culicoides obsoletus (Diptera: Ceratopogonidae), vector of bluetongue and Schmallenberg viruses, revealed by landscape genetic analyses

BACKGROUND

In the last two decades, recurrent epizootics of bluetongue virus and Schmallenberg virus have been reported in the western Palearctic region. These viruses affect domestic cattle, sheep, goats and wild ruminants and are transmitted by native hematophagous midges of the genus Culicoides (Diptera: Ceratopogonidae). Culicoides dispersal is known to be stratified, i.e. due to a combination of dispersal processes occurring actively at short distances and passively or semi-actively at long distances, allowing individuals to jump hundreds of kilometers.

METHODS

Here, we aim to identify the environmental factors that promote or limit gene flow of Culicoides obsoletus, an abundant and widespread vector species in Europe, using an innovative framework integrating spatial, population genetics and statistical approaches. A total of 348 individuals were sampled in 46 sites in France and were genotyped using 13 newly designed microsatellite markers.

RESULTS

We found low genetic differentiation and a weak population structure for C. obsoletus across the country. Using three complementary inter-individual genetic distances, we did not detect any significant isolation by distance, but did detect significant anisotropic isolation by distance on a north-south axis. We employed a multiple regression on distance matrices approach to investigate the correlation between genetic and environmental distances. Among all the environmental factors that were tested, only cattle density seems to have an impact on C. obsoletus gene flow.

CONCLUSIONS

The high dispersal capacity of C. obsoletus over land found in the present study calls for a re-evaluation of the impact of Culicoides on virus dispersal, and highlights the urgent need to better integrate molecular, spatial and statistical information to guide vector-borne disease control.

Disjunction between canola distribution and the genetic structure of its recently described pest, the canola flower midge (Contarinia brassicola)

Population genomics is a useful tool to support integrated pest management as it can elucidate population dynamics, demography, and histories of invasion. Here, we use a restriction site-associated DNA sequencing approach combined with whole-genome amplification (WGA) to assess genomic population structure of a newly described pest of canola, the diminutive canola flower midge, Contarinia brassicola. Clustering analyses recovered little geographic structure across the main canola production region but differentiated several geographically disparate populations at edges of the agricultural zone. Given a lack of alternative hypotheses for this pattern, we suggest these data support alternative hosts for this species and thus our canola-centric view of this midge as a pest has limited our understanding of its biology. These results speak to the need for increased surveying efforts across multiple habitats and other potential hosts within Brassicaceae to improve both our ecological and evolutionary knowledge of this species and contribute to effective management strategies. We additionally found that use of WGA prior to library preparation was an effective method for increasing DNA quantity of these small insects prior to restriction site-associated DNA sequencing and had no discernible impact on genotyping consistency for population genetic analysis; WGA is therefore likely to be tractable for other similar studies that seek to randomly sample markers across the genome in small organisms.

Genetic Diversity of Culicoides stellifer (Diptera: Ceratopogonidae) in the Southeastern United States Compared With Sequences From Ontario, Canada

Much of the bluetongue (BT) and epizootic hemorrhagic disease (EHD) research in North America focuses on white-tail deer and Culicoides sonorensis (Wirth & Jones) (Diptera: Ceratopogonidae), though several other biting midge species have been suggested as vectors. Culicoides stellifer (Coquillett) has been associated with hosts susceptible to hemorrhagic disease (HD), and more recently, specimens from Florida have tested positive for EHD and BT viral RNA. If C. stellifer is acting as a vector, this could have an impact on the distribution of HD in North America. To determine if gene flow is occurring across the range of C. stellifer within the southeast United States, a mitochondrial haplotype analysis was performed using the COI gene. Our haplotype network showed no population structure in C. stellifer from Florida, Texas, and South Carolina, as the overall genetic divergence between these sites was equal to the genetic divergence within each. We also compared these haplotypes to published sequences of C. stellifer collected in Ontario, Canada. Surprisingly, the genetic diversity of the flies from Ontario was two times greater than what was observed between the southeast U.S. collection sites. This considerable divergence could be evidence of a cryptic species. A better understanding of the connectivity between C. stellifer populations across all of North America will give insight into the distribution of HD. Our results show that gene flow is occurring between sites in the southeastern United States and potentially throughout the eastern distribution of the species.

Zucchini yellow mosaic virus Genomic Sequences from Papua New Guinea: Lack of Genetic Connectivity with Northern Australian or East Timorese Genomes, and New Recombination Findings

Zucchini yellow mosaic virus (ZYMV) isolates were obtained in Papua New Guinea (PNG) from cucumber (Cucumis sativus) or pumpkin (Cucurbita spp.) plants showing mosaic symptoms growing at Kongop in the Mount Hagen District, Western Highlands Province, or Zage in the Goroka District, Eastern Highlands Province. The samples were blotted onto FTA cards, which were sent to Australia, where they were subjected to high-throughput sequencing. When the coding regions of the nine new ZYMV genomic sequences found were compared with those of 64 other ZYMV sequences from elsewhere, they grouped together, forming new minor phylogroup VII within ZYMV's major phylogroup A. Genetic connectivity was lacking between ZYMV genomic sequences from PNG and its neighboring countries, Australia and East Timor; the closest match between a PNG and any other genomic sequence was a 92.8% nucleotide identity with a sequence in major phylogroup A's minor phylogroup VI from Japan. When the RDP5.2 recombination analysis program was used to compare 66 ZYMV sequences, evidence was obtained of 30 firm recombination events involving 41 sequences, and all isolates from PNG were recombinants. There were 21 sequences without recombination events in major phylogroup A, whereas there were only 4 such sequences within major phylogroup B. ZYMV's P1, Cl, N1a-Pro, P3, CP, and NIb regions contained the highest evidence of recombination breakpoints. Following removal of recombinant sequences, seven minor phylogroups were absent (I, III, IV, V, VI, VII, and VIII), leaving only minor phylogroups II and IX. By contrast, when a phylogenetic tree was constructed using recombinant sequences with their recombinationally derived tracts removed before analysis, five previous minor phylogroups remained unchanged within major phylogroup A (II, III, IV, V, and VII) while four formed two new merged phylogroups (I/VI and VIII/IX). Absence of genetic connectivity between PNG, Australian, and East Timorese ZYMV sequences, and the 92.8% nucleotide identity between a PNG sequence and the closest sequence from elsewhere, suggest that a single introduction may have occurred followed by subsequent evolution to adapt to the PNG environment. The need for enhanced biosecurity measures to protect against potentially damaging virus movements crossing the seas separating neighboring countries in this region of the world is discussed.

Genomic analysis of bluetongue virus episystems in Australia and Indonesia

The distribution of bluetongue viruses (BTV) in Australia is represented by two distinct and interconnected epidemiological systems (episystems)—one distributed primarily in the north and one in the east. The northern episystem is characterised by substantially greater antigenic diversity than the eastern episystem; yet the forces that act to limit the diversity present in the east remain unclear. Previous work has indicated that the northern episystem is linked to that of island South East Asia and Melanesia, and that BTV present in Indonesia, Papua New Guinea and East Timor, may act as source populations for new serotypes and genotypes of BTV to enter Australia’s north. In this study, the genomes of 49 bluetongue viruses from the eastern episystem and 13 from Indonesia were sequenced and analysed along with 27 previously published genome sequences from the northern Australian episystem. The results of this analysis confirm that the Australian BTV population has its origins in the South East Asian/Melanesian episystem, and that incursions into northern Australia occur with some regularity. In addition, the presence of limited genetic diversity in the eastern episystem relative to that found in the north supports the presence of substantial, but not complete, barriers to gene flow between the northern and eastern Australian episystems. Genetic bottlenecks between each successive episystem are evident, and appear to be responsible for the reduction in BTV genetic diversity observed in the north to south–east direction.

Genotyping of whole genome amplified reduced representation libraries reveals a cryptic population of Culicoides brevitarsis in the Northern Territory, Australia

BACKGROUND

The advent of genotyping by Next Generation Sequencing has enabled rapid discovery of thousands of single nucleotide polymorphism (SNP) markers and high throughput genotyping of large populations at an affordable cost. Genotyping by sequencing (GBS), a reduced representation library sequencing method, allows highly multiplexed sequencing of genomic subsets. This method has limitations for small organisms with low amounts of genomic DNA, such as the bluetongue virus (BTV) vectors, Culicoides midges.

RESULTS

This study employed the GBS method to isolate SNP markers de novo from whole genome amplified Culicoides brevitarsis genomic DNA. The individuals were collected from regions representing two different Australian patterns of BTV strain distribution: the Northern Territory (NT) and the east coast. We isolated 8145 SNPs using GBS. Phylogenetic analysis conducted using the filtered 3263 SNPs revealed the presence of a distinct C. brevitarsis sub-population in the NT and this was confirmed by analysis of mitochondrial DNA. Two loci showed a very strong signal for selection and were unique to the NT population. Bayesian analysis with STRUCTURE indicated a possible two-population cluster.

CONCLUSIONS

The results suggest that genotyping vectors with high density markers in combination with biological and environmental data is useful. However, more extensive sampling over a wider spatial and temporal range is needed. The presence of sub-structure in populations and loci under natural selection indicates the need for further investigation of the role of vectors in shaping the two Australian systems of BTV transmission. The described workflow is transferable to genotyping of small, non-model organisms, including arthropod vectors of pathogens of economic and medical importance.

Population Genetic Structure and Potential Incursion Pathways of the Bluetongue Virus Vector Culicoides brevitarsis (Diptera: Ceratopogonidae) in Australia

Culicoides brevitarsis is a vector of the bluetongue virus (BTV), which infects sheep and cattle. It is an invasive species in Australia with an assumed Asian/South East Asian origin. Using one mitochondrial marker (i.e., part of the cytochrome oxidase subunit I gene) and six nuclear markers, we inferred population genetic structure and possible incursion pathways for Australian C. brevitarsis. Nine mitochondrial haplotypes, with low nucleotide sequence diversity (0.0-0.7%) among these, were identified in a sample of 70 individuals from seven sites. Both sets of markers revealed a homogeneous population structure, albeit with evidence of isolation by distance and two genetically distinct clusters distributed along a north-to-south cline. No evidence of a cryptic species complex was found. The geographical distribution of the mitochondrial haplotypes is consistent with at least two incursion pathways into Australia since the arrival of suitable livestock hosts. By contrast, 15 mitochondrial haplotypes, with up to four times greater nucleotide sequence diversity (0.0-2.9%) among these, were identified in a sample of 16 individuals of the endemic C. marksi (sampled from a site in South Australia and another in New South Wales). A phylogenetic tree inferred using the mitochondrial marker revealed that the Australian and Japanese samples of C. brevitarsis are as evolutionarily different from one another as some of the other Australian species (e.g., C. marksi, C. henryi, C. pallidothorax) are. The phylogenetic tree placed four of the species endemic to Australia (C. pallidothorax, C. bundyensis, C. marksi, C. henryi) in a clade, with a fifth such species (C. bunrooensis) sharing a common ancestor with that clade and a clade comprising two Japanese species (C. verbosus, C. kibunensis).

Delineation of the population genetic structure of Culicoides imicola in East and South Africa

Background

Culicoides imicola Kieffer, 1913 is the main vector of bluetongue virus (BTV) and African horse sickness virus (AHSV) in Sub-Saharan Africa. Understanding the population genetic structure of this midge and the nature of barriers to gene flow will lead to a deeper understanding of bluetongue epidemiology and more effective vector control in this region.

Methods

A panel of 12 DNA microsatellite markers isolated de novo and mitochondrial DNA were utilized in a study of C. imicola populations from Africa and an outlier population from the Balearic Islands. The DNA microsatellite markers and mitochondrial DNA were also used to examine a population of closely related C. bolitinos Meiswinkel midges.

Results

The microsatellite data suggest gene flow between Kenya and south-west Indian Ocean Islands exist while a restricted gene flow between Kenya and South Africa C. imicola populations occurs. Genetic distance correlated with geographic distance by Mantel test. The mitochondrial DNA analysis results imply that the C. imicola populations from Kenya and south-west Indian Ocean Islands (Madagascar and Mauritius) shared haplotypes while C. imicola population from South Africa possessed private haplotypes and the highest nucleotide diversity among the African populations. The Bayesian skyline plot suggested a population growth.

Conclusions

The gene flow demonstrated by this study indicates a potential risk of introduction of new BTV serotypes by wind-borne infected Culicoides into the Islands. Genetic similarity between Mauritius and South Africa may be due to translocation as a result of human-induced activities; this could impact negatively on the livestock industry. The microsatellite markers isolated in this study may be utilised to study C. bolitinos, an important vector of BTV and AHSV in Africa and identify sources of future incursions.

Electronic supplementary material

The online version of this article (doi:10.1186/s13071-015-1277-4) contains supplementary material, which is available to authorized users.