Population Genetic Approaches Research Articles

Regulatory regions in natural transposable element insertions drive interindividual differences in response to immune challenges in Drosophila

BackgroundVariation in gene expression underlies interindividual variability in relevant traits including immune response. However, the genetic variation responsible for these gene expression changes remains largely unknown. Among the non-coding variants that could be relevant, transposable element insertions are promising candidates as they have been shown to be a rich and diverse source of cis-regulatory elements.ResultsIn this work, we use a population genetics approach to identify transposable element insertions likely to increase the tolerance of Drosophila melanogaster to bacterial infection by affecting the expression of immune-related genes. We identify 12 insertions associated with allele-specific expression changes in immune-related genes. We experimentally validate three of these insertions including one likely to be acting as a silencer, one as an enhancer, and one with a dual role as enhancer and promoter. The direction in the change of gene expression associated with the presence of several of these insertions is consistent with an increased survival to infection. Indeed, for one of the insertions, we show that this is the case by analyzing both natural populations and CRISPR/Cas9 mutants in which the insertion is deleted from its native genomic context.ConclusionsWe show that transposable elements contribute to gene expression variation in response to infection in D. melanogaster and that this variation is likely to affect their survival capacity. Because the role of transposable elements as regulatory elements is not restricted to Drosophila, transposable elements are likely to play a role in immune response in other organisms as well.

Identification of incompletely penetrant variants and interallelic interactions in autosomal recessive disorders by a population-genetic approach.

We aimed to identify incompletely penetrant (IP) variants and interallelic interactions in autosomal recessive disorders by a population-genetic approach. Genotype and clinical data were collected from 9038 patients of European origin with ASL, ATP7B, CAPN3, CFTR, CTNS, DHCR7, GAA, GALNS, GALT, IDUA, MUT, NPHS1, NPHS2, PAH, PKHD1, PMM2, or SLC26A4-related disorders. We calculated the relative allele frequency of each pathogenic variant (n = 1936) to the loss-of-function (LOF) variants of the corresponding gene in the patient ( ) and the general population ( ) and estimated the penetrance of each variant by calculating their ratio: (V/LOF ratio). We classified all variants as null or hypomorphic based on the associated clinical phenotype. We found 25 variants, 29% of the frequent 85 variants, to be underrepresented in the patient population (V/LOF ratio <30% with p < 7.22 × 10-5 ), including 22 novel ones in the ASL, CAPN3, CFTR, GAA, GALNS, PAH, andPKHD1 genes. In contrast to the completely penetrant variants (CP), the majority of the IP variants were hypomorphic (IP: 16/18, 88%; CP: 177/933, 19.0%; p = 5.12 × 10-10 ). Among them, only the NPHS2 R229Q variant was subject to interallelic interactions. The proposed algorithm identifies frequent IP variants and estimates their penetrance and interallelic interactions in large patient cohorts.

Genetic diversity of Huaya India (Melicoccus oliviformis Kunth), a neglected Neotropical fruit crop

Currently, some species of Sapindaceae are important fruit crops worldwide. The Huaya India (Meliccocus oliviformis, Sapindaceae) is a neglected Neotropical fruit tree consumed locally in the Maya Lowlands of Mexico, where it exists in both wild and domesticated forms. Our objective was to evaluate the genetic diversity of the Huaya India in its possible domestication area and thus generate knowledge that serves as the basis for a commercial management. A total of 450 individuals collected from 15 natural vegetation sites and 15 Maya villages, were characterized using nine microsatellite loci and population genetics approaches were applied. STRUCTURE, Neighbor-Joining and PCoA analyses suggested the existence of three main groups: 1) one composed by 14 natural vegetation sites, 2) one integrated by 10 Maya villages plus one natural vegetation site, 3) one composed by five Maya villages. At the species level, genetic differentiation was high (FST = 0.562) and gene flow was low (Nm = 0.395); between genetic groups, differentiation was low and gene flow was high. Genetic diversity was low at the level species (HE = 0.19) and higher in the group composed for only natural vegetation sites. When we considered only two groups (natural vegetation sites vs Maya villages) to explore a possible bottleneck as a consequence of human management, the natural vegetation sites showed higher, and significant, genetic diversity (HE = 0.231) than the Maya villages (HE = 0.152). This study can serve as a basis to develop management strategies for Huaya India in the Maya Lowlands of Mexico, but without compromising its conservation.

Ancient Adaptative Evolution of ACE2 in East Asians.

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has been posing an unprecedented challenge to global public health. SARS-CoV-2 and several other coronaviruses utilize angiotensin-converting enzyme 2 (ACE2) as their entry receptors. The ACE2 gene has been found to experience episodic positive selection across mammals. However, much remains unknown about how the ACE2 gene evolved in human populations. Here, we use population genetics approaches to investigate the evolution of the ACE2 gene in 26 human populations sampled globally. We find the ACE2 gene exhibits an extremely low nucleotide diversity in the East Asian populations. Strong signals of selective sweep are detected in the East Asian populations, but not in the other human populations. The selective sweep in ACE2 is estimated to begin in East Asian populations ∼23,600 years ago. Our study provides novel insights into the evolution of the ACE2 gene within human populations.

Genetic Structure of the Goniopora lobata and G. djiboutiensis Species Complex Is Better Explained by Oceanography Than by Morphological Characteristics

Species delimitation of closely related corals is often challenging due to high intraspecies morphological variation and phenotypic plasticity with a lack of characteristic features and scarcity of relevant molecular markers. Goniopora spp. are one such coralline group, and the species status of Goniopora lobata and Goniopora djiboutiensis, an Indian and Pacific Ocean hermatypic coral species complex, has been questioned on the basis of previous molecular and morphological analyses. To further examine the species boundaries between G. lobata and G. djiboutiensis in Japan, specimens collected from areas spanning from Ryukyu Island to temperate Japanese regions were morphologically identified based on traditional morphological descriptions. Then, the genetic structure of the G. lobata and G. djiboutiensis species complex was examined using six newly developed microsatellite markers. The majority of the collected specimens had intermediate morphologies, and a STRUCTURE analysis using the LOCPRIOR model based on typical G. lobata and G. djiboutiensis morphology indicated that there were no genetic differences between these morphologies. On the other hand, STRUCTURE analysis based on oceanographic regions revealed that there was a genetic break between the temperate and subtropical regions. This weak genetic break corresponded with the Kuroshio-associated barrier, which prevents larval transport between subtropical and temperate regions. This study confirms that the current morphological identification criteria for G. lobata and G. djiboutiensis do not match the existing genetic boundaries and thus the two should be regarded as a species complex. This study also highlighted the robustness of using a multi-locus population genetic approach, including a geographic context, to confirm the species boundaries of closely related species.

Genetic and Morphological Approach for Western Corn Rootworm Resistance Management

The western corn rootworm (WCR), is one of the most serious pests of maize in the United States. In this study, we aimed to find a reliable pattern of difference related to resistance type using population genetic and geometric morphometric approaches. To perform a detailed population genetic analysis of the whole genome, we used single nucleotide polymorphisms (SNPs) markers. For the morphometric analyses, hindwings of the resistant and non-resistant WCR populations from the US were used. Genetic results showed that there were some differences among the resistant US populations. The low value of pairwise FST = 0.0181 estimated suggests a lack of genetic differentiation and structuring among the putative populations genotyped. However, STRUCTURE analysis revealed three genetic clusters. Heterozygosity estimates (HO and HE) over all loci and populations were very similar. There was no exact pattern, and resistance could be found throughout the whole genome. The geometric morphometric results confirmed the genetic results, with the different genetic populations showing similar wing shape. Our results also confirmed that the hindwings of WCR carry valuable genetic information. This study highlights the ability of geometric morphometrics to capture genetic patterns and provides a reliable and low-cost alternative for preliminary estimation of population structure. The combined use of SNPs and geometric morphometrics to detect resistant variants is a novel approach where morphological traits can provide additional information about underlying population genetics, and morphology can retain useful information about genetic structure. Additionally, it offers new insights into an important and ongoing area of pest management on how to prevent or delay pest evolution towards resistant populations, minimizing the negative impacts of resistance.

Identifying Genomic Regions Targeted During Eggplant Domestication Using Transcriptome Data.

Identifying genes and traits that have diverged during domestication provides key information of importance for maintaining and even increasing yield and nutrients in existing crops. A “bottom-up” population genetics approach was used to identify signatures of selection across the eggplant genome, to better understand the process of domestication. RNA-seq data were obtained for 4 wild eggplants (Solanum insanum L.) and 16 domesticated eggplants (S. melongena L.) and mapped to the eggplant genome. Single-nucleotide polymorphism (SNPs) exhibiting signatures of selection in domesticates were identified as those exhibiting high FST between the 2 populations (evidence of significant divergence) and low π for the domesticated population (indicative of a selective sweep). Some of these regions appear to overlap with previously identified quantitative trait loci for domestication traits. Genes in regions of linkage disequilibrium surrounding these SNPs were searched against the Arabidopsis thaliana and tomato genomes to find orthologs. Subsequent gene ontology (GO) enrichment analysis identified over-representation of GO terms related to photosynthesis and response to the environment. This work reveals genomic changes involved in eggplant domestication and improvement, and how this compares to observed changes in the tomato genome, revealing shared chromosomal regions involved in the domestication of both species.

Dispersal and mating patterns determine the fate of naturally dispersed populations: evidence from Bombina orientalis

BackgroundIn contrast to the explosive increase of a population following biological invasion, natural dispersal, i.e., when a population disperses from its original range into a new range, is a passive process that is affected by resources, the environment, and other factors. Natural dispersal is also negatively impacted by genetic drift and the founder effect. Although the fates of naturally dispersed populations are unknown, they can adapt evolutionarily over time to the new environment. Can naturally dispersed populations evolve beneficial adaptive strategies to offset these negative effects to maintain their population in a stable state?ResultsThe current study addressed this question by focusing on the toad Bombina orientalis, the population of which underwent natural dispersal following the Last Glacial Maximum in Northeast Asia. Population genetic approaches were used to determine the genetic structure, dispersal pattern, and mating system of the population of B. orientalis in northeast China (Northern population). The results showed that this northern population of B. orientalis is a typical naturally dispersed population, in which the stable genetic structure and high level of genetic diversity of the population have been maintained through the long-distance biased dispersal behavior of males and the pattern of promiscuity within the population.ConclusionsOur findings suggest that naturally dispersed populations can evolve effective adaptive strategies to maintain a stable population. Different species may have different strategies. The relevance of these maintenance mechanisms for naturally dispersed populations provide a new perspective for further understanding the processes of speciation and evolution.

Effector Avr4 in Phytophthora infestans Escapes Host Immunity Mainly Through Early Termination.

Effector genes play critical roles in the antagonistic interactions between plants and pathogens. However, knowledge of mutation mechanisms and evolutionary processes in effector genes and the contribution of climatic factors to the evolution of effector genes are fragmented but important in sustainable management of plant diseases and securing food supply under changing climates. Here, we used a population genetic approach to explore the evolution of the Avr4 gene in Phytophthora infestans, the causal agent of potato blight. We found that the Avr4 gene exhibited a high genetic diversity generated by point mutation and sequence deletion. Frameshifts caused by a single base-pair deletion at the 194th nucleotide position generate two stop codons, truncating almost the entire C-terminal, which is important for effector function and R4 recognition in all sequences. The effector is under natural selection for adaptation supported by comparative analyses of population differentiation (FST) and isolation-by-distance between Avr4 sequences and simple sequence repeat marker loci. Furthermore, we found that local air temperature was positively associated with pairwise FST in the Avr4 sequences. These results suggest that the evolution of the effector gene is influenced by local air temperature, and the C-terminal truncation is one of the main mutation mechanisms in the P. infestans effector gene to circumvent the immune response of potato plants. The implication of these results to agricultural and natural sustainability in future climate conditions is discussed.

Integration of population and functional genomics to understand mechanisms of artemisinin resistance in Plasmodium falciparum

Resistance to antimalarial drugs, and in particular to the artemisinin derivatives and their partner drugs, threatens recent progress toward regional malaria elimination and eventual global malaria eradication. Population-level studies utilizing whole-genome sequencing approaches have facilitated the identification of regions of the parasite genome associated with both clinical and in vitro drug-resistance phenotypes. However, the biological relevance of genes identified in these analyses and the establishment of a causal relationship between genotype and phenotype requires functional characterization. Here we examined data from population genomic and transcriptomic studies in the context of data generated from recent functional studies, using a new population genetic approach designed to identify potential favored mutations within the region of a selective sweep (iSAFE). We identified several genes functioning in pathways now known to be associated with artemisinin resistance that were supported in early population genomic studies, as well as potential new drug targets/pathways for further validation and consideration for treatment of artemisinin-resistant Plasmodium falciparum. In addition, we establish the utility of iSAFE in identifying positively-selected mutations in population genomic studies, potentially accelerating the time to functional validation of candidate genes.

Parental selection based on molecular information under a population genetics approach

The correct choice of parents that will compose optimal segregating populations is the key to success for breeding programs. It was postulated the hypothesis that this choice of these parents could be made based on information of molecular markers analyzed in the context of population structure. Ten parental populations were simulated and 45 hybrid combinations were obtained from the dialel crosses. Each population consisted of 200 individuals with 50 independent loci. The populations were evaluated for the Hardy-Weinberg Equilibrium (HWE), Coefficient of Inbreeding (F), Heterozygosity (H), and the Polymorphic Information Content (PIC). Genetic diversity between pairs of parental populations was evaluated using five dissimilarity measures. Values of Mantel correlation were obtained for the pairs of the dissimilarity matrices, and the PIC, H, and F values ​​were obtained in the hybrid combinations. All parental populations were under HWE, and the combination that emerged from this condition was the hybrid 3x5, with only 26% of the loci manifesting HWE. This same hybrid was among those with lower F estimates and higher values ​​of H, which indicated the existence of greater divergence between their parentals. There was agreement on the indication of the more and less divergent hybrid combinations for the dissimilarity measures. This fact is important because the variability, associated with the good average potential, are important criteria for the formation of an initial population in breeding programs of any kind, involving sexual processes.

Effect of social structure and introduction history on genetic diversity and differentiation

Invasive species are a global threat to biodiversity, and understanding their history and biology is a major goal of invasion biology. Population-genetic approaches allow insights into these features, as population structure is shaped by factors such as invasion history (number, origin and age of introductions) and life-history traits (e.g., mating system, dispersal capability). We compared the relative importance of these factors by investigating two closely related ants, Tetramorium immigrans and Tetramorium tsushimae, that differ in their social structure and invasion history in North America. We used mitochondrial DNA sequences and microsatellite alleles to estimate the source and number of introduction events of the two species, and compared genetic structure among native and introduced populations. Genetic diversity of both species was strongly reduced in introduced populations, which also differed genetically from native populations. Genetic differentiation between ranges and the reduction in microsatellite diversity were more severe in the more recently introduced and supercolonial T.tsushimae. However, the loss of mitochondrial haplotype diversity was more pronounced in T.immigrans, which has single-queen colonies and was introduced earlier. Tetramorium immigrans was introduced at least twice from Western Europe to North America and once independently to South America. Its monogyny might have limited genetic diversity per introduction, but new mutations and successive introductions over a long time may have added to the gene pool in the introduced range. Polygyny in T.tsushimae probably facilitated the simultaneous introduction of several queens from a Japanese population to St. Louis, USA. In addition to identifying introduction pathways, our results reveal how social structure can influence the population-genetic consequences of founder events.

On Deleterious Mutations in Perennials: Inbreeding Depression, Mutation Load, and Life-History Evolution.

AbstractIn angiosperms, perennials typically present much higher levels of inbreeding depression than annuals. One hypothesis to explain this pattern stems from the observation that inbreeding depression is expressed across multiple life stages in angiosperms. It posits that increased inbreeding depression in more long-lived species could be explained by differences in the way mutations affect fitness, through the life stages at which they are expressed. In this study, we investigate this hypothesis. We combine a physiological growth model and multilocus population genetics approaches to describe a full genotype-to-phenotype-to-fitness map. We study the behavior of mutations affecting growth or survival and explore their consequences in terms of inbreeding depression and mutation load. Although our results agree with empirical data only within a narrow range of conditions, we argue that they may point us toward the type of traits capable of generating high inbreeding depression in long-lived species-that is, traits under sufficiently strong selection, on which selection decreases sharply as life expectancy increases. Then we study the role deleterious mutations maintained at mutation-selection balance may play in the joint evolution of growth and survival strategies.

Deep divergence between island populations in lichenized fungi

Macaronesia is characterized by a high degree of endemism and represents a noteworthy system to study the evolutionary history of populations and species. Here, we compare the population-genetic structure in three lichen-forming fungi, the widespread Lobaria pulmonaria and two Macaronesian endemics, L. immixta and L. macaronesica, based on microsatellites. We utilize population genetic approaches to explore population subdivision and evolutionary history of these taxa on the Canary Islands, Madeira, Azores, and the western Iberian Peninsula. A common feature in all species was the deep divergence between populations on the Azores, a pattern expected by the large geographic distance among islands. For both endemic species, there was a major split between archipelagos. In contrast, in the widespread L. pulmonaria, divergent individuals were distributed across multiple archipelagos, suggesting a complex evolutionary history involving repeated migration between islands and mainland.

Estimating divergence times from DNA sequences.

The patterns of genetic variation within and among individuals and populations can be used to make inferences about the evolutionary forces that generated those patterns. Numerous population genetic approaches have been developed in order to infer evolutionary history. Here, we present the “Two-Two (TT)” and the “Two-Two-outgroup (TTo)” methods; two closely related approaches for estimating divergence time based in coalescent theory. They rely on sequence data from two haploid genomes (or a single diploid individual) from each of two populations. Under a simple population-divergence model, we derive the probabilities of the possible sample configurations. These probabilities form a set of equations that can be solved to obtain estimates of the model parameters, including population split times, directly from the sequence data. This transparent and computationally efficient approach to infer population divergence time makes it possible to estimate time scaled in generations (assuming a mutation rate), and not as a compound parameter of genetic drift. Using simulations under a range of demographic scenarios, we show that the method is relatively robust to migration and that the TTo method can alleviate biases that can appear from drastic ancestral population size changes. We illustrate the utility of the approaches with some examples, including estimating split times for pairs of human populations as well as providing further evidence for the complex relationship among Neandertals and Denisovans and their ancestors.

Large genetic diversity and strong positive selection in F-box and GPCR genes among the wild isolates of Caenorhabditis elegans.

The F-box and chemosensory GPCR (csGPCR) gene families are greatly expanded in nematodes, including the model organism Caenorhabditis elegans, compared with insects and vertebrates. However, the intraspecific evolution of these two gene families in nematodes remain unexamined. In this study, we analyzed the genomic sequences of 330 recently sequenced wild isolates of C. elegans using a range of population genetics approaches. We found that F-box and csGPCR genes, especially the Srw family csGPCRs, showed much more diversity than other gene families. Population structure analysis and phylogenetic analysis divided the wild strains into eight non-Hawaiian and three Hawaiian subpopulations. Some Hawaiian strains appeared to be more ancestral than all other strains. F-box and csGPCR genes maintained a great amount of the ancestral variants in the Hawaiian subpopulation and their divergence among the non-Hawaiian subpopulations contributed significantly to population structure. F-box genes are mostly located at the chromosomal arms and high recombination rate correlates with their large polymorphism. Moreover, using both neutrality tests and extended haplotype homozygosity analysis, we identified signatures of strong positive selection in the F-box and csGPCR genes among the wild isolates, especially in the non-Hawaiian population. Accumulation of high-frequency-derived alleles in these genes was found in non-Hawaiian population, leading to divergence from the ancestral genotype. In summary, we found that F-box and csGPCR genes harbor a large pool of natural variants, which may be subjected to positive selection. These variants are mostly mapped to the substrate-recognition domains of F-box proteins and the extracellular and intracellular regions of csGPCRs, possibly resulting in advantages during adaptation by affecting protein degradation and the sensing of environmental cues, respectively.

Molecular systematics and biogeography of an Australian soil‐burrowing cockroach with polymorphic males, Geoscapheus dilatatus (Blattodea: Blaberidae)

AbstractAn iconic group of arid‐adapted insects is the Australian soil‐burrowing cockroaches (Blaberidae: Geoscapheinae), large, wingless insects that evolved burrowing behaviour and associated forms in parallel from wood feeding ancestors in the subfamily Panesthiinae. A particularly problematic taxon within the Geoscapheinae is Geoscapheus dilatatus (Saussure, 1864), which might represent a species complex and whose delimitation has been complicated for decades by the species harbouring polymorphic males. Males can be divided into two main morphs: individuals possessing horn‐like protrusions on the anterior margin of the pronotum (‘tuberculate’) and those without these characters (‘non‐tuberculate’). A less common, third form consists of individuals that possess tubercles but are far larger than other tuberculate males and occur solely to the north of the species' distribution (‘atypical’ tuberculates). Here, we make use of whole mitochondrial genomes and nuclear ribosomal RNA data from individuals across the range of G. dilatatus to conduct the first phylogenetic analysis of this species to date. We recover all tuberculate males (including atypical forms) as monophyletic and the derived form of G. dilatatus, having evolved only once in this species, whereas non‐tuberculate forms are paraphyletic. Fossil‐calibrated molecular clock analysis revealed that the divergence between these two forms occurred during the late Miocene approximately 6.7 Mya, concurrent with an expansion of the continent's drier biomes. Environmental niche modelling suggests that tuberculate male forms are more climatically tolerant than their more restricted non‐tuberculate counterparts, and both forms' predicted fundamental niches are strongly limited by rainfall. Three species delimitation analyses implemented here failed to consistently delimit G. dilatatus beyond a single species. Ultimately, population genetic approaches paired with additional sampling will be necessary to determine these findings more concretely, but at present, we do not consider the results presented here sufficient to delimit G. dilatatus based on morphological differences found in the species' polymorphic males alone.

Comprehensive genotyping of a Brazilian cassava (Manihot esculenta Crantz) germplasm bank: insights into diversification and domestication

Key messageBrazilian cassava diversity was characterized through population genetics and clustering approaches, highlighting contrasted genetic groups and spatial genetic differentiation.Cassava (Manihot esculenta Crantz) is a major staple root crop of the tropics, originating from the Amazonian region. In this study, 3354 cassava landraces and modern breeding lines from the Embrapa Cassava Germplasm Bank (CGB) were characterized. All individuals were subjected to genotyping-by-sequencing (GBS), identifying 27,045 single-nucleotide polymorphisms (SNPs). Identity-by-state and population structure analyses revealed a unique set of 1536 individuals and 10 distinct genetic groups with heterogeneous linkage disequilibrium (LD). On this basis, a density of 1300–4700 SNP markers were selected for large-effect quantitative trait loci (QTL) detection. Identified genetic groups were further characterized for population genetics parameters including minor allele frequency (MAF), observed heterozygosity ({H}_{o}), effective population size estimate widehat{{(N}_{e}}) and polymorphism information content (PIC). Selection footprints and introgressions of M. glaziovii were detected. Spatial population structure analysis revealed five ancestral populations related to distinct Brazilian ecoregions. Estimation of historical relationships among identified populations suggests an early population split from Amazonian to Atlantic forest and Caatinga ecoregions and active gene flows. This study provides a thorough genetic characterization of ex situ germplasm resources from cassava’s center of origin, South America, with results shedding light on Brazilian cassava characteristics and its biogeographical landscape. These findings support and facilitate the use of genetic resources in modern breeding programs including implementation of association mapping and genomic selection strategies.

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

BackgroundIn 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.MethodsHere, 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.ResultsWe 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.ConclusionsThe 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.Graphical

Genetic structure and population demographics of Hypnoidus bicolor (Coleoptera: Elateridae) in the Canadian Prairies.

Following banning of the pesticide lindane in most counties, wireworms (i.e., the soil-living larval stages of click beetles) have become major pests of a variety of economically important field crops. Hypnoidus bicolor is a common pest species in the Canadian Prairies. However, little is known about its life history, which impedes the development of effective integrated pest management (IPM) strategies. Population genetic approaches have the potential to assist in the development of IPM. We sequenced a 622-bp fragment of the COX1 gene from 326 H. bicolor wireworm and click beetles collected from 13 localities on the Canadian Prairies. Two genetically distinct (>4.66% sequence divergence) clades were identified, suggesting that they may be part of a species complex. Clade A predominated and increased in prevalence the further east samples were collected, whereas the opposite was true for clade B. Clade B appears to be comprised of two mitochondrial DNA groups, however, one group was represented by only one haplotype. Both clades were characterized by uneven gene flow among populations with low levels of regional genetic structuring. Clade A appeared to have undergone population and range expansions, which may coincide with the advent of intensive agriculture practices in the prairies. Knowledge of species composition and population structure is important for the development of effective IPM strategies but is often lacking for wireworms. Our study fills these knowledge gaps for a predominant pest species in the prairies, H. bicolor, by providing robust evidence for cryptic forms and characterizing its dispersal patterns and population dynamics. © 2021 Society of Chemical Industry.