Mitochondrial DNA signature for range-wide populations of Bicyclus anynana suggests a rapid expansion from recent refugia

Genetic diversity and population genetics of Schistosoma haematobium isolated from children in Lusaka and Siavonga districts, Zambia

Urogenital schistosomiasis remains a pervasive health challenge in rural Zambian communities. This study explores the molecular epidemiology and genetic diversity of Schistosoma haematobium using mitochondrial genes (cox1 and nadh1). Urine samples from 421 children in Siavonga and Lusaka districts, Zambia, were collected between December 2020 and February 2022. Microscopy and DNA extraction facilitated the identification of S. haematobium, followed by amplification, sequencing, and phylogenetic analysis of cox1 and nadh1 genes. Phylogenetic analysis revealed clustering with samples from mainland African countries, emphasizing shared haplotypes. Both mitochondrial genes exhibited substantial diversity, with 5 haplotypes from 37 cox1 sequences and 12 haplotypes from 23 nadh1 sequences. High haplotype diversity (0.621-0.808) and low nucleotide diversity (0.00181-0.03288) were observed. Siavonga and Lusaka districts shared the majority of S. haematobium haplotypes. Molecular variance and genetic differentiation analysis indicated variations within populations rather than between populations (cox1: -0.025, nadh1: 0.01646). These findings suggest a limited differentiation between S. haematobium populations in Siavonga and Lusaka, potentially indicating gene flow. Tajima's test revealed negative values, indicating a departure from neutrality, introduction of rare alleles, and recent population expansion. This study contributes essential insights into S. haematobium population genetics, crucial for effective urogenital schistosomiasis control in Zambia.

Patterns of Genetic Diversity and Gene Flow Associated With an Aridity Gradient in Populations of Common Mole-rats, Cryptomys hottentotus hottentotus

Genetic adaptation is the change of a population toward a phenotype that best fits the present ecological conditions of the environment it inhabits. As environmental conditions change, allele frequencies shift, resulting in different populations of the same species possessing genetic variation and divergent phenotypes. Cooperatively breeding common mole-rats (Cryptomys hottentotus hottentotus) inhabit environments along an aridity gradient in South Africa, which provides an opportunity for local genetic adaptations to occur. Using one mitochondrial gene (cytochrome b) and 3,540 SNP loci across the whole genome, we determined the phylogenetic relationship, population structure and genetic diversity of five populations of C. h. hottentotus located along an aridity gradient. Mitochondrial data identified population-specific clades that were less distinct in the two mesic populations, potentially indicating historical or recent gene flow, or the retention of ancestral haplotypes. Arid and semi-arid populations formed a distinct cluster from the non-arid populations. Genetic diversity and gene flow were higher in arid-dwelling individuals, suggesting greater connectivity and interactions between colonies in arid regions in comparison to mesic ones. Using an Aridity Index, we determined that isolation by environment, rather than isolation by geographical distance, best explains the genetic distance between the populations. Further analyses using target loci may determine if there are differing underlying genetic adaptations among populations of C. h. hottentotus. These analyses could help unravel population differences in response to environmental factors within a subspecies of bathyergid mole-rat and determine the adaptive capacity of this small nonmigratory subterranean rodent species in response to aridification in the face of climate change.

Zygosity-based sex determination in a butterfly drives hypervariability of Masculinizer

Nature has devised many ways of producing males and females. Here, we report on a previously undescribed mechanism for Lepidoptera that functions without a female-specific gene. The number of alleles or allele heterozygosity in a single Z-linked gene (BaMasc) is the primary sex-determining switch in Bicyclus anynana butterflies. Embryos carrying a single BaMasc allele develop into WZ (or Z0) females, those carrying two distinct alleles develop into ZZ males, while (ZZ) homozygotes initiate female development, have mismatched dosage compensation, and die as embryos. Consequently, selection against homozygotes has favored the evolution of spectacular allelic diversity: 205 different coding sequences of BaMasc were detected in a sample of 246 females. The structural similarity of a hypervariable region (HVR) in BaMasc to the HVR in Apis mellifera csd suggests molecular convergence between deeply diverged insect lineages. Our discovery of this primary switch highlights the fascinating diversity of sex-determining mechanisms and underlying evolutionary drivers.

Genotyping and epidemiological distribution of diarrhea-causing isolates of Giardia duodenalis in southeastern part of West Bengal, India

The prevalence and genetic diversity of the protozoan pathogen Giardia duodenalis have been extensively studied worldwide. There is currently a lack of data regarding the genetic variability of the organism in eastern India. Understanding the circulating genotypes and associated risk factors is crucial for effective planning and implementing control measures. Therefore, the objective of the study was to conduct an epidemiological study to determine the prevalence and identify the various genotypes present. This survey adds to our knowledge on the occurrence and distribution of Giardia genotypes in the studied region. The overall prevalence was found to be 6.8%. This parasitic infection was significantly associated with two age groups, i.e., >0-5 years and >5-12 years. Using a multilocus genotyping method, we genotyped 52 human Giardia isolates that were obtained from diarrheal patients. Two distinct assemblages were found in the population-30.8% belonged to assemblage A; 63.5% belonged to assemblage B, prevalent in the population; and 5.7% belonged to a combined assemblage A+B. Sub-assemblage AII was found in 17.3% of the cases, followed by sub-assemblage AI (13.5%). High levels of genetic diversity were found within the population of assemblage B undergoing balancing selection. Overall, the high prevalence of the parasite observed, particularly among children, raises a major concern and necessitates implementation of robust control measures. Furthermore, we report the presence of numerous unique genotypes, circulating in this limited geographical boundary, which can be useful dataset for future studies.

Preliminary population studies of the grassland swallowtail butterfly Euryades corethrus (Lepidoptera, Papilionidae)

Euryades corethrus is a Troidini butterfly (Papilionidae, Papilioninae), endemic to grasslands in southern Brazil, Uruguay, Argentina and Paraguay. Formerly abundant, nowadays it is in the Red list of endangered species for those areas. During its larval stage, it feeds on Aristolochia spp, commonly found in southern grasslands. These native grassland areas are diminishing, being converted to crops and pastures, causing habitat loss for Aristolochia and E. corethrus. This study aimed to assess the genetic diversity, population structure and demographic history of E. corethrus. We sampled eight populations from Rio Grande do Sul, Brazil and based on Cytochrome Oxidase subunit I (COI) molecular marker, our results suggest a low genetic variability between populations, presence of gene flow and, consequently, lack of population structure. A single maternally inherited-genetic marker is insufficient for population-level decisions, but barcoding is a useful tool during early stages of population investigation, bringing out genomic diversity patterns within the target species. Those populations likely faced a bottleneck followed by a rapid expansion during the last glaciation and subsequent stabilization in effective population size. Habitat loss is a threat, which might cause isolation, loss of genetic variability and, ultimately, extinction of E. corethrus if no habitat conservation policy is adopted.

Conservation Genetics of the Critically Endangered Southern River Terrapin (Batagur affinis) in Malaysia: Genetic Diversity and Novel Subspecies Distribution Ranges

A population genetics study was carried out on the Southern River terrapin (Batagur affinis) from four places in Peninsular Malaysia: Pasir Gajah, Kemaman (KE), Terengganu; Bukit Pinang (BP), Kedah; Bota Kanan (BK), Perak; and Bukit Paloh, Kuala Berang (KB), Terengganu. The goal of this study is to identify genetic differences in two subspecies of B. affinis in Malaysia. No previous reports were available on the genetic diversity, phylogenetic relationships and matrilineal hereditary structure of these terrapin populations in Malaysia. The sequencing identified 46 single nucleotide polymorphisms that defined six mitochondrial haplotypes in the Southern River terrapins. Tajima's D test and Fu's Fs neutrality tests were performed to evaluate the signatures of recent historical demographic events. Based on the tests, the B. affinis edwardmolli was newly subspecies identified in the west coast-northern region of Kedah state. In addition, the B. affinis edwardmolli in Bukit Paloh, Kuala Berang (KB), Terengganu (Population 4), was shown to have a single maternal lineage compared to other populations. Low genetic diversity, but significant genetic differences, were detected among the studied Southern River terrapin populations.

Genetic Signature of Pinctada fucata Inferred from Population Genomics: Source Tracking of the Invasion in Mischief Reef of Nansha Islands

Among the anthropogenic stresses that marine ecosystems face, biological invasions are one of the major threats. Recently, as a result of increasingly intense anthropogenic disturbance, numerous marine species have been introduced to their non-native ranges. However, many introduced species have uncertain original sources. This prevents the design and establishment of methods for controlling or preventing these introduced species. In the present study, genomic sequencing and population genetic analysis were performed to detect the geographic origin of the introduced Pinctada fucata population in the Mischief Reef of the South China Sea. The results of population genetic structure analysis showed a close relationship between the Mischief Reef introduced population and the Lingshui population, indicating that Lingshui may be the potential geographical origin. Furthermore, lower heterozygosity and nucleotide diversity were observed in the introduced population in Mischief Reef, indicating lower genetic diversity than in other native populations. We also identified some selected genomic regions and genes of the introduced population, including genes related to temperature and salinity tolerance. These genes may play important roles in the adaptation of the introduced population. Our study will improve our understanding of the invasion history of the P. fucata population. Furthermore, the results of the present study will also facilitate further control and prevention of invasion in Mischief Reef, South China Sea.

Honey Bee Habitat Sharing Enhances Gene Flow of the Parasite Nosema ceranae

Host-parasite co-evolution is a process of reciprocal, adaptive genetic change. In natural conditions, parasites can shift to other host species, given both host and parasite genotypes allow this. Even though host-parasite co-evolution has been extensively studied both theoretically and empirically, few studies have focused on parasite gene flow between native and novel hosts. Nosema ceranae is a native parasite of the Asian honey bee Apis cerana, which infects epithelial cells of mid-guts. This parasite successfully switched to the European honey bee Apis mellifera, where high virulence has been reported. In this study, we used the parasite N. ceranae and both honey bee species as model organisms to study the impacts of two-host habitat sharing on parasite diversity and virulence. SNVs (Single Nucleotide Variants) were identified from parasites isolated from native and novel hosts from sympatric populations, as well as novel hosts from a parapatric population. Parasites isolated from native hosts showed the highest levels of polymorphism. By comparing the parasites isolated from novel hosts between sympatric and parapatric populations, habitat sharing with the native host significantly enhanced parasite diversity, suggesting there is continuing gene flow of parasites between the two host species in sympatric populations.

Alternative splicing in seasonal plasticity and the potential for adaptation to environmental change

Seasonal plasticity is accomplished via tightly regulated developmental cascades that translate environmental cues into trait changes. Little is known about how alternative splicing and other posttranscriptional molecular mechanisms contribute to plasticity or how these mechanisms impact how plasticity evolves. Here, we use transcriptomic and genomic data from the butterfly Bicyclus anynana, a model system for seasonal plasticity, to compare the extent of differential expression and splicing and test how these axes of transcriptional plasticity differ in their potential for evolutionary change. Between seasonal morphs, we find that differential splicing affects a smaller but functionally unique set of genes compared to differential expression. Further, we find strong support for the novel hypothesis that spliced genes are more susceptible than differentially expressed genes to erosion of genetic variation due to selection on seasonal plasticity. Our results suggest that splicing plasticity is especially likely to experience genetic constraints that could affect the potential of wild populations to respond to rapidly changing environments.

Morphological and molecular characterization of Castniidae (Lepidoptera) associated to sugarcane in Colombia

The giant sugarcane borer, Telchin licus, has been reported as an economically important sugarcane pest in Colombia; however, its taxonomic status has been scarcely investigated and previous reports offer an ambiguous characterization of both the immature and adult stages. The objective of this work is to identify Telchin species affecting sugarcane and alternative hosts in different departments of the country by integrating molecular analysis and conventional morphology. To date, T. licus has been found in the departments of Caquetá, Casanare, and Meta, while T. atymnius has been found in Antioquia, Caldas, Nariño, and Valle del Cauca. Sugarcane, Musaceae, and Heliconiaceae have been found to be hosts to both species. Additionally, the species T. cacica has also been registered in the department of Nariño, affecting heliconias and plantains. Genetic variation within the species allowed differentiation at the molecular level of subspecies of T. licus and T. atymnius, confirming that the subspecies present in Colombia are T. licus magdalena, T. atymnius humboldti, and T. atymnius atymnius. The haplotype diversity of populations is closely related to their geographical distribution, indicating low gene flow between populations and possible speciation inside the country. Analysis of genetic variance showed significant differences among and within T. atymnius populations, which may suggest a high genetic structure along the regions where it is found and the possible presence of additional subspecies to those previously reported. To understand the geographical and environmental conditions that determine the pest's distribution in Colombia, this information needs to be complemented with ecological considerations of possible geographical isolation and association of alternative hosts.

Analysis of the insecticide resistance mechanism in Anopheles culicifacies sensu lato from a malaria-endemic state in India

BACKGROUND

Understanding the dynamics and mechanisms of insecticide resistance in malaria vectors is crucial for vector control activities. The present study investigates the level of insecticide resistance in Anopheles culicifacies and explores the role of two main mechanisms in conferring resistance target site insensitivity and metabolic resistance.

METHODS

A. culicifacies mosquitoes were collected and the voltage-gated sodium channel (VGSC) gene was amplified and sequenced to analyse the knockdown resistance (kdr) mutations. Further, a non-experimental homology model was generated to investigate the effect of kdr mutations on the conformation of protein. Metabolic resistance was determined using bioassay-based resistant and susceptible mosquitoes and the expression levels of the genes CYP6Z1 and GSTe2 were compared between the two groups.

RESULTS

Sequence analysis of the VGSC gene revealed the presence of L1014F (n=48 [17%]), L1014S and V1010L (n=5 [1.7%]) mutations in the study area. In gene expression studies, a significant upregulation of CYP6Z1 in deltamethrin-resistant (fold change 243.62; p=0.02) mosquitoes and that of GSTe2 in dichlorodiphenyltrichloroethane (fold change 403.45; p=0.01) and alpha-cypemethrin resistant (fold change 217.51; p=0.0005) mosquitoes was observed.

CONCLUSIONS

The study revealed that expression of the genes (CYP6Z1 and GSTe2) conferring metabolic resistance play a key role in insecticide resistance in A. culicifacies populations in central India. However, mutations L101F, L10104S and V10101L also have a role to some extent in spreading resistance. GeneBank accession numbers: MW559058, MW559059 and MW559060 Cover Image: Workflow of Chimera-Modeller interface. In the top window of Chimera's multi-align viewer the sequence alignment of VGSC proteins of human (pdb id_6AGF), cockroach (pdb id 5XOM) and A. culicifacies (ACT176122.1) is shown. The dialog box in the middle is of the comparative modelling tool of Modeller. The A. culicifacies sequence is designated as the target while human and cockroach sequences are templates. Upon selection of the template sequences in the dialog box, the structures of the respective proteins are displayed in the Chimera window. As the run is completed, the results are displayed in the form of a list of models with their scores in a table.

Meta-Analysis of Mitochondrial DNA Control Region Diversity to Shed Light on Phylogenetic Relationship and Demographic History of African Sheep (Ovis aries) Breeds

To improve sheep breeding and conservation of genetic resources, the mitochondrial DNA control region (mtDNA CR) of 399 sequences of African indigenous sheep breeds from previously published research articles were meta-analyzed to elucidate their phylogenetic relationship, diversity, and demographic history. A total of 272 haplotypes were found, of which 207 were unique and a high level of mtDNA CR variability was observed. Generally, the number of polymorphic sites, nucleotide and haplotype diversity were high (284, 0.254 ± 0.012 and 0.993 ± 0.002, respectively). The median-joining (MJ) network of haplotypes produced three major haplogroups (A, B and C), with haplogroup B being dominant. A mixture of populations suggests a common matrilineal origin and lack of and/or a weak phylogeographic structure. Mismatch analysis showed recent expansion of North African breeds, whereas East African and continental populations exhibited selection pressures for adaptation. A slight historical genetic difference was also observed between the fat tail and thin tail sheep breeds. However, further investigations are required using more samples and long sequence segments to achieve deeper levels of conclusions on the African sheep phylogenetic relationship. The present meta-analysis results contribute to the general understanding of African native sheep populations for improved management of sheep genetic resources.

DNA barcode analysis of the endangered green turtle (Chelonia mydas) in Mexico1

Technological and analytical advances to study evolutionary biology, ecology, and conservation of green turtles (Chelonia mydas) are realized through molecular approaches including DNA barcoding. We characterized the usefulness of COI DNA barcodes in green turtles in Mexico to better understand genetic divergence and other genetic parameters of this species. We analyzed 63 sequences, including 25 from green turtle field specimens collected from the Gulf of Mexico and from the Mexican Pacific and 38 already present in the Barcode of Life Data Systems (BOLD). A total of 13 haplotypes were identified with four novel haplotypes from the Pacific Ocean and three novel haplotypes from the Atlantic Ocean. Intraspecific distance values among COI gene sequences by two different models were 0.01, demonstrating that there is not a subdivision for green turtle species. Otherwise, the interspecific distance interval ranged from 0.07 to 0.13, supporting a clear subdivision among all sea turtle species. Haplotype and total nucleotide diversity values of the COI gene reflect a medium genetic diversity average. Green turtles of the Mexican Pacific showed common haplotypes to some Australian and Chinese turtles, but different from the haplotypes of the Mexican Atlantic. COI analysis revealed new haplotypes and confirmed that DNA barcodes were useful for evaluation of the population diversity of green turtles in Mexico.

Sleeping with the enemy: unravelling the symbiotic relationships between the scale worm Neopolynoe chondrocladiae (Annelida: Polynoidae) and its carnivorous sponge hosts

The North Atlantic deep-water polynoid worm Neopolynoe chondrocladiae is involved in an exceptional symbiotic relationship with two hosts: the carnivorous sponges Chondrocladia robertballardi and Chondrocladia virgata. While this is an obligate symbiotic relationship, its real nature is unclear. We used a multidisciplinary approach to narrow down the type of symbiotic relationship between symbiont and hosts. Molecular connectivity analyses using COI and 16S suggest that N. chondrocladiae has high potential for dispersal, connecting sites hundreds of kilometres apart, likely aided by oceanographic currents. Microbial analyses on different anatomical parts of five Chondrocladia species suggest that the presence of the worm in C. robertballardi does not affect the microbiome of the sponge. MicroCT analysis on N. chondrocladiae show that it has dorsally oriented parapodia, which might prevent the worm from getting trapped in the sponge. A faecal pellet recovered from the worm suggests that the polynoid feeds on the crustacean prey captured by the sponge, something corroborated by our stable isotope analysis. Light and confocal microscopy images suggest that N. chondrocladiae elytra produce bioluminescence. We propose that the worm might use bioluminescence as a lure for prey (increasing the food available for both the sponge and the polynoid) and thus fuelling a mutualistic relationship.

Population diversity and virulence characteristics of Cryptococcus neoformans/C. gattii species complexes isolated during the pre-HIV-pandemic era

Cryptococcosis has become a major global health problem since the advent of the HIV pandemic in 1980s. Although its molecular epidemiology is well-defined, using isolates recovered since then, no pre-HIV-pandemic era epidemiological data exist. We conducted a molecular epidemiological study using 228 isolates of the C. neoformans/C. gattii species complexes isolated before 1975. Genotypes were determined by URA5 restriction fragment length polymorphism analysis and multi-locus sequence typing. Population genetics were defined by nucleotide diversity measurements, neutrality tests, and recombination analysis. Growth at 37°C, melanin synthesis, capsule production, and urease activity as virulence factors were quantified. The pre-HIV-pandemic isolates consisted of 186 (81.5%) clinical, 35 (15.4%) environmental, and 7 (3.1%) veterinary isolates. Of those, 204 (89.5%) belonged to C. neoformans VNI (64.0%), VNII (14.9%) and VNIV (10.5%) while 24 (10.5%) belonged to C. gattii VGIII (7.5%), VGI (2.6%) and VGII (0.5%). Among the 47 sequence types (STs) identified, one of VNII and 8 of VNIV were novel. ST5/VNI (23.0%) in C. neoformans and ST75/VGIII (25.0%) in C. gattii were the most common STs in both species complexes. Among C. neoformans, VNIV had the highest genetic diversity (Hd = 0.926) and the minimum recombination events (Rm = 10), and clinical isolates had less genetic diversity (Hd = 0.866) than environmental (Hd = 0.889) and veterinary isolates (Hd = 0.900). Among C. gattii, VGI had a higher nucleotide diversity (π = 0.01436) than in VGIII (π = 0.00328). The high-virulence genotypes (ST5/VNI and VGIIIa/serotype B) did not produce higher virulence factors levels than other genotypes. Overall, high genetic variability and recombination rates were found for the pre-HIV-pandemic era among strains of the C. neoformans/C. gattii species complexes. Whole genome analysis and in vivo virulence studies would clarify the evolution of the genetic diversity and/or virulence of isolates of the C. neoformans/C. gattii species complexes during the pre- and post-HIV-pandemic eras.

Since the beginning of the HIV pandemic in 1980, infections due to isolates of the Cryptococcus neoformans/C. gattii species complexes have caused many deaths worldwide, especially in the HIV-infected population. Annually, approximately one-third, of all AIDS-related deaths,—representing more than 1,000,000 cases,—are caused by cryptococcosis. Since 1980, extensive molecular epidemiological surveys have been conducted, and the VNI molecular type has been found to be responsible for more than 90% of cryptococcosis in HIV patients. Whether the high VNI prevalence is associated with the HIV pandemic remains controversial as information on the isolates of the pre-HIV pandemic era is lacking. Therefore, this study of the molecular epidemiology and in vitro characteristics of the strains from the pre-HIV-pandemic era was undertaken. We found that only 64% of cryptococcosis was caused by VNI, and 9 sequence types existed only in the pre-HIV pandemic era. Unlike what was already known about the strains collected during the HIV pandemic era, ST5 and VGIIIa,—supposedly high virulence genotypes,—did not express higher virulence factors than other genotypes. These results implied that the HIV pandemic altered both the molecular epidemiology and virulence of Cryptococcus neoformans/C. gattii species complexes have been altered during HIV pandemic. However, detailed mechanism of these alteration remains to be deciphered further.

New insights into the coexistence of Contracaecum rudolphii A and Contracaecum rudolphii B (Nematoda: Anisakidae) in Phalacrocorax carbo sinensis from Sardinia: genetic variability and phylogenetic analysis

Contracaecum sp. nematodes are important parasites of fish eating birds that can cause animal health problems. In the present study, specimens of Contracaecum rudolphii sensu lato, from the great cormorant Phalacrocorax carbo sinensis from Sardinia, were characterized based on morphological and molecular data. The morphological analysis allowed to identify all the fourth stage larvae (n = 1918) as Contracaecum sp., and adults, male (n = 5845) and female (n = 8312), as C. rudolphii sensu lato. Population genetics and phylogenetic relationships were inferred based on mitochondrial and nuclear markers. Multiple sequence alignment of the ribosomal internal transcribed spacer showed the coexistence of C. rudolphii A (n = 157), C. rudolphii B (n = 22) and a rare heterozygote of these species. Moreover, mitochondrial markers, namely NADH dehydrogenase subunits I (nad1), cytochrome c oxidase subunit (cox1 and cox2) and small subunit of rRNA (rrnS), showed that the studied C. rudolphii A populations had undergone bottleneck, or founder effect event, subsequent to a rapid population growth and expansion. The observed heterozygote is with a mitochondrial pattern of C. rudolphii B. Although, both Contracaecum species showed high genetic diversity, no genetic structure between localities was detected. Phylogenetic reconstructions supported the paraphyly of the avian Contracaecum species including C. ogmorhini (parasite of otariids).

Host use diversification during range shifts shapes global variation in Lepidopteran dietary breadth

Niche breadths tend to be greater at higher latitudes. This pattern is frequently assumed to emerge from the cumulative effects of multiple, independent local adaptation events along latitudinal environmental gradients, although evidence that generalization is more beneficial at higher-latitude locations remains equivocal. Here I propose an alternative hypothesis: that latitudinal variation in niche breadths emerges as a non-adaptive consequence of range shift dynamics. Based on analysis of a global dataset comprising more than 6,934 globally distributed dietary records from 4,410 Lepidopteran species, this hypothesis receives robust support. Population-level dietary niche breadths are better explained by the relative position of the population within its geographic range and the species' poleward range extent than by the latitude of diet observation. Broader diets are observed closer to poleward range limits and in species that have attained higher latitudes. Moreover, latitudinal variation in diet breadth is more prominent within and among species undergoing rapid, contemporary range shifts than for species with more stable ranges. Together these results suggest that latitudinal patterns in niche breadth represent a transient and emergent property of recent geographic range dynamics and need not require underlying gradients in selective agents or fitness trade-offs. The results have wide-ranging implications for global ecology and for anticipating changes in host use during ongoing distributional shifts of pests and disease vectors.

Genetic diversity and phylogenetic relationships of tsetse flies of the palpalis group in Congo Brazzaville based on mitochondrial cox1 gene sequences

BACKGROUND

Despite the morphological characterization established in the 1950s and 1960s, the identity of extant taxa that make up Glossina fuscipes (s.l.) in the Congo remains questionable. Previous claims of overlap between G. fuscipes (believed to be G. f. quanzensis) and G. palpalis palpalis around Brazzaville city further complicate the taxonomic status and population dynamics of the two taxa. This study aimed to determine the phylogenetic relationships between G. fuscipes (s.l.) and G. p. palpalis and to assess genetic variation among G. fuscipes (s.l.) populations in Congo Brazzaville.

METHODS

We collected 263 G. fuscipes (s.l.) from northern and central regions, and 65 G. p. palpalis from southern part of the country. The mitochondrial cytochrome c oxidase subunit 1 (cox1) gene was amplified using taxa-specific primer pairs. Sequence data were analyzed in DnaSP and Arlequin to assess the genetic diversity, differentiation and demographic history of G. fuscipes (s.l.) populations.

RESULTS

The general BLAST analysis yielded a similarity of 99% for G. fuscipes (s.l.) and G. p. palpalis. BLASTn analysis for G. fuscipes (s.l.) showed > 98% identity with GenBank sequences for G. fuscipes (s.l.), with BEMB population showing 100% similarity with G. f. fuscipes. Glossina fuscipes (s.l.) populations showed high haplotype diversity (H = 46, Hd = 0.884), moderate nucleotide diversity ( = 0.012) and moderate (FST = 0.072) to high (FST = 0.152) genetic differentiation. Most of the genetic variation (89.73%) was maintained within populations. The mismatch analysis and neutrality tests indicated recent tsetse population expansions.

CONCLUSIONS

Phylogenetic analysis revealed minor differences between G. fuscipes (s.l.) and G. p. palpalis. Genetic diversity of G. fuscipes (s.l.) was high in the populations sampled except one. Genetic differentiation ranged from moderate to high among subpopulations. There was a restricted gene flow between G. fuscipes (s.l.) populations in the north and central part of the country. Genetic signatures based on cox1 showed recent expansion and recovery of G. fuscipes (s.l.) populations from previous bottlenecks. To fully understand the species distribution limits, we recommend further studies involving a wider sampling scheme including the swampy Mossaka focus for G. fuscipes (s.l.) and the entire range of G. p. palpalis in South Congo.

Phylogeography Approach of Diloboderus abderus (Coleoptera: Melolonthidae) in the Southern Cone of America

Diloboderus abderus (Sturm, 1826) (Coleoptera: Melolonthidae) is a serious soil pest of corn, wheat, oat, and natural and cultivated pastures in Argentina, Paraguay, Uruguay, and southern Brazil. Despite its economic importance, the genetic diversity and population structure of D. abderus remain unknown. We sequenced a fragment of the mitochondrial gene cytochrome oxidase I region (COI), of six populations of D. abderus from the Southern Cone of America. The mtDNA marker revealed a high haplotype diversity, high pairwise F_ST values, and significant genetic variations among populations. No correlation was found between genetic and geographical distances, yet the most common haplotype (Dab01) was present in four out of the six populations. Analysis of molecular variance showed that most of the variation was within populations of D. abderus. Tajima's D and Fu's F_S tests indicated no evidence that D. abderus populations are under recent expansion. Our results indicate that genetic-based traits will likely remain localized or spread slowly, and management strategies need to be undertaken on a small scale.

Population genetic structure and assessment of allochronic divergence in the Macoun’s Arctic (Oeneis macounii) butterfly

Patterns in the genetic variation of species can be used to infer their specific demographic and evolutionary history and provide insight into the general mechanisms underlying population divergence and speciation. The Macoun’s Arctic (Oeneis macounii (W.H. Edwards, 1885); MA) butterfly occurs across Canada and parts of the northern United States in association with jack pine (Pinus banksiana Lamb.) and lodgepole pine (Pinus contorta Douglas ex Loudon). MA’s current distribution is highly fragmented, and the extent of reproductive isolation among allopatric populations is unknown. Furthermore, although MA is biennial, adults emerge every year in some populations. These populations presumably consist of two alternate-year cohorts, providing the opportunity for sympatric divergence via allochronic isolation. Using mitochondrial DNA (mtDNA) and amplified fragment length polymorphism (AFLP) markers, we analyzed MA’s genetic structure to determine the current and historical role of allopatric and allochronic isolation in MA population divergence. Both markers revealed high diversity and a low, but significant, degree of spatial structure and pattern of isolation by distance. Phylogeographic structure was generally absent, with low divergence among mtDNA haplotypes. MA likely exhibits low dispersal and gene flow among most allopatric populations; however, there was no evidence of differentiation resulting from allochronic isolation for sympatric cohorts.

Phylogeographic pattern of the plane leaf miner, Phyllonorycter platani (STAUDINGER, 1870) (Lepidoptera: Gracillariidae) in Europe

Background

The plane leaf miner, Phyllonorycter platani is a widely distributed insect species on plane trees and has a well-documented colonisation history in Europe over the last century. However, phylogeographic data of the species are lacking.

Results

We analysed 284 individuals from 38 populations across Europe, Asia, and North America. A 1242 bp fragment of the mitochondrial COI gene and an 893 bp fragment of the 28S rDNA has been Sanger sequenced. Twenty-four haplotypes were detected on the COI gene, and two alleles were identified on the 28S rDNA. We revealed two distinct clades for both markers reflecting the geographic origins, Asia and Europe. The genetic distance between the two main clades is 2.08% on the COI gene and 0.10% on the nuclear DNA.

An overlapping zone of the two clades was found across Eastern Europe and the Anatolian Peninsula. We detected heterozygote individuals of the 28S rDNA gene in Moldavia, Ukraine and in the southern part of Turkey. These suggest that the two clades can hybridise. Furthermore, the presence of European type homozygote individuals has been confirmed in the southern part of Turkey as well.

Conclusions

We have shown that both post-glacial recolonization and recent expansion events influenced the present genetic structure of P. platani. The genetic patterns revealed at least two refugia during the last ice age: one in the Balkan Peninsula and the other in the Caucasus region. Recent expansion was detected in some European and Central Asian populations. The two main clades (Europe/Asia) show definite genetic differences; however, several hybrid individuals were found in the overlapping zone as well (stretching over Eastern Europe and the Anatolian Peninsula). Discrepancies in mitochondrial and nuclear data indicate introgressions in the southern part of the Anatolian Peninsula.

Electronic supplementary material

The online version of this article (10.1186/s12862-018-1240-z) contains supplementary material, which is available to authorized users.

Phylogeographic analysis reveals high genetic structure with uniform phenotypes in the paper wasp Protonectarina sylveirae (Hymenoptera: Vespidae)

Swarm-founding wasps are endemic and common representatives of neotropical fauna and compose an interesting social tribe of vespids, presenting both complex social characteristics and uncommon traits for a eusocial group, such as the absence of castes with distinct morphology. The paper wasp Protonectarina sylveirae (Saussure) presents a broad distribution from Brazil, Argentina and Paraguay, occurring widespread in the Atlantic rainforest and arboreal Caatinga, being absent in the Amazon region. Given the peculiar distribution among swarm-founding wasps, an integrative approach to reconstruct the evolutionary history of P. sylveirae in a spatial-temporal framework was performed to investigate: the presence of genetic structure and its relationship with the geography, the evolution of distinct morphologic lineages and the possible historical event(s) in Neotropical region, which could explain the observed phylogeographic pattern. Individuals of P. sylveirae were obtained from populations of 16 areas throughout its distribution for DNA extraction and amplification of mitochondrial genes 12S, 16S and COI. Analysis of genetic diversity, construction of haplotype net, analysis of population structure and dating analysis of divergence time were performed. A morphometric analysis was also performed using 8 measures of the body of the adult (workers) to test if there are morphological distinction among populations. Thirty-five haplotypes were identified, most of them exclusively of a group and a high population structure was found. The possibility of genetic divergence because of isolation by distance was rejected. Morphological analysis pointed to a great uniformity in phenotypes, with only a small degree of differentiation between populations of south and the remaining. Divergence time analysis showed a Middle/Late Miocene origin, a period where an extensive marine ingression occurred in South America. Divergence of haplogroups began from the Plio/Pleistocene boundary and the last glacial maximum most likely modeled the current distribution of species, even though it was not the cause of genetic breaks.

Integrative taxonomic methods reveal an incorrect synonymisation of the South African Pseudonereis podocirra (Schmarda) as the widespread Pseudonereis variegata (Grube) from Chile

The nereidid Pseudonereis variegata (Grube, 1866) described from Chile includes 14 synonymised species from 10 type localities with a discontinuous distribution, but no taxonomic or molecular studies have investigated the status of this species outside Chile. Two synonymised species, Mastigonereis podocirra Schmarda, 1861 and Nereis (Nereilepas) stimpsonis Grube, 1866, were described from South Africa and investigated here using morphological examination. MtCOI species delimitation analyses and morphology were used to determine the status of P. variegata in South Africa. Morphological examination revealed that museum and freshly collected specimens from South Africa that conform to the general description of P. variegata are similar to M. podocirra and N. stimpsonis with respect to the consistent absence of homogomph spinigers in the inferior neuropodial fascicle, expanded notopodial ligules and the subterminal attachment of dorsal cirri in posterior parapodia. The synonymy of M. podocirra and N. stimpsonis as P. variegata are rejected and P. podocirra, comb. nov. is reinstated. Morphologically, Pseudonereis podocirra differed from specimens from Chile with regard to the numbers of paragnaths, the absence of homogomph spinigers and changes in parapodial morphology along the body. Independence of these species was further supported by genetic distances, automatic barcode gap discovery and multi-rate Poisson tree process species delimitation analyses of 77 mtCOI sequences. Haplotype network revealed no genetic structuring within the South African populations. http://zoobank.org/urn:lsid:zoobank.org:pub:F0B1A5AF-9CE9-4A43-ACCF-17117E1C2F21

Diversity of tortricid moths in apple orchards: evidence for a cryptic species of Grapholita (Lepidoptera: Tortricidae) from China

Understanding herbivore diversity both at the species and genetic levels is a key to effective pest management. We examined moth samples from multiple locations from a major apple growing region in China. For specimen collection, we used a pheromone trap designed to attract Grapholita molesta (Busck) (Lepidoptera: Tortricidae). Surprisingly, we found a second species captured at high proportions. Its external morphology (e.g., male genitalia and forewing coloration) was the same as for Grapholita funebrana Treitschke (Lepidoptera: Tortricidae) specimens from Europe. However, the barcode sequence of the mitochondrial gene cytochrome oxidase I (COI) diverged markedly between specimens from China and Europe, and the genetic distance value between the specimens from the two regions as estimated using the Juke-Cantor (JC) model amounted to 0.067. These morphological and molecular findings together point to a cryptic species in G. funebrana from China. Further molecular analyses based on COI and COII genes revealed its extremely high genetic diversity, indicating that the origin of this species includes the sampling region. Moreover, molecular data suggest that this species passed through a recent population expansion. This is the first report on a cryptic species in G. funebrana, as well as the first report on its genetic diversity.

Mitochondrial DNA sequence variation in Iranian native dogs

The dog mtDNA diversity picture from wide geographical sampling but from a small number of individuals per region or breed, displayed little geographical correlation and high degree of haplotype sharing between very distant breeds. For a clear picture, we extensively surveyed Iranian native dogs (n = 305) in comparison with published European (n = 443) and Southwest Asian (n = 195) dogs. Twelve haplotypes related to haplogroups A, B and C were shared by Iranian, European, Southwest Asian and East Asian dogs. In Iran, haplotype and nucleotide diversities were highest in east, southeast and northwest populations while western population had the least. Sarabi and Saluki dog populations can be assigned into haplogroups A, B, C and D; Qahderijani and Kurdi to haplogroups A, B and C, Torkaman to haplogroups A, B and D while Sangsari and Fendo into haplogroups A and B, respectively. Evaluation of population differentiation using pairwise F_ST generally revealed no clear population structure in most Iranian dog populations. The genetic signal of a recent demographic expansion was detected in East and Southeast populations. Further, in accordance with previous studies on dog-wolf hybridization for haplogroup d2 origin, the highest number of d2 haplotypes in Iranian dog as compared to other areas of Mediterranean basin suggests Iran as the probable center of its origin. Historical evidence showed that Silk Road linked Iran to countries in South East Asia and other parts of the world, which might have probably influenced effective gene flow within Iran and these regions. The medium nucleotide diversity observed in Iranian dog calls for utilization of appropriate management techniques in increasing effective population size.

The genetic structure of hypoderatid mites (Actinotrichida: Astigmata) parasitizing great cormorant (Phalacrocorax carbo) during host post-breeding dispersal in Milicz, SW Poland

The astigmatid family Hypoderatidae includes over 80 mite species exhibiting peculiar life cycle. Deutonymphs are parasitic instars inhabiting subcutaneous or visceral tissues of birds and rodents, while all other instars are nidicolous forms. In this study we investigated genetic diversity of deutonymphs of two hypoderatid species, Neottialges evansi and Phalacrodectes gaudi, collected from 16 individuals of the great cormorant Phalacrocorax carbo in SW Poland during host post-breeding dispersal. The initial alternative hypotheses were: (1) populations of deutonymphs in both mite species found in tissues of particular bird individuals form genetically distinct populations, (2) mites are panmictic among hosts. The topologies of NJ phylogenetic trees and median-joining haplotype networks reconstructed for COI haplotypes revealed lack of hostdependent genetic structure in populations of N. evansi and P. gaudi. Furthermore, high haplotype diversity (Hd) and low nucleotide diversity (Pi) prove high genetic differentiation of both mite species. We concluded, that unlimited dispersal of mites among cormorant specimens could be explained by host specific breeding behavior: free mating between cormorants originating from different European populations and low contribution of reunited pairs in breeding colonies in subsequent breeding seasons, reuse of the same nest material by different members of the colony, and host behavior called prospecting.

Genetic Polymorphism of Aedes albopictus Population Inferred From ND5 Gene Variabilities In Subang Jaya, Malaysia

This study was performed to establish the genetic variability of Aedes albopictus within Subang Jaya, Selangor, Malaysia, by using the nicotinamide adenine dinucleotide dehydrogenase 5 subunit (ND5) mitochondrial DNA (mtDNA) marker. A total of 90 samples were collected from 9 localities within an area of the Subang Jaya Municipality. Genetic variability was determined through the amplification and sequencing of a fragment of the ND5 gene. Eight distinct mtDNA haplotypes were identified. The evolutionary relationship of the local haplotypes alongside 28 reference strains was used to construct a phylogram, the analysis of which revealed low genetic differentiation in terms of both nucleotide and haplotype diversity. Bayesian method was used to infer the phylogenetic tree, revealing a unique relationship between local isolates. The study corroborates the reliability of ND5 to identify distinct lineages for polymorphism-based studies and supplements the existing body of knowledge regarding its genetic diversity. This in turn could potentially aid existing vector control strategies to help mitigate the risk and spread of the dengue virus.

Differentiation in putative male sex pheromone components across and within populations of the African butterfly Bicyclus anynana as a potential driver of reproductive isolation

Sexual traits are often the most divergent characters among closely related species, suggesting an important role of sexual traits in speciation. However, to prove this, we need to show that sexual trait differences accumulate before or during the speciation process, rather than being a consequence of it. Here, we contrast patterns of divergence among putative male sex pheromone (pMSP) composition and the genetic structure inferred from variation in the mitochondrial cytochrome oxidase 1 and nuclear CAD loci in the African butterfly Bicyclus anynana (Butler, 1879) to determine whether the evolution of "pheromonal dialects" occurs before or after the differentiation process. We observed differences in abundance of some shared pMSP components as well as differences in the composition of the pMSP among B. anynana populations. In addition, B. anynana individuals from Kenya displayed differences in the pMSP composition within a single population that appeared not associated with genetic differences. These differences in pMSP composition both between and within B. anynana populations were as large as those found between different Bicyclus species. Our results suggest that "pheromonal dialects" evolved within and among populations of B. anynana and may therefore act as precursors of an ongoing speciation process.

Mitochondrial population structure and post-glacial dispersal of longnose sucker Catostomus catostomus in Labrador, Canada: evidence for multiple refugial origins and limited ongoing gene flow

Two hundred and eighty-seven longnose sucker Catostomus catostomus were collected from 14 lakes in Labrador, 52 from three lakes in Ontario, 43 from two lakes in British Columbia and 32 from a lake in Yukon; a total of 414 in all. The resulting 34 haplotypes (20 in Labrador) contained moderate haplotypic diversity (h = 0·657) and relatively low nucleotide diversity (π = 3·730 × 10(-3) . Mean ϕST (0·453, P < 0·05) over all populations revealed distinct genetic structuring among C. catostomus populations across Canada, based on province, which was validated by the analysis and spatial analysis of molecular variance (c. 80% variation between provinces). These results probably reflect the historical imprint of recolonization from different refugia and possibly indicate limited ongoing gene flow within provinces. A haplotype network revealed one major and two minor clades within Labrador that were assigned to the Atlantic, Beringian and Mississippian refugia, respectively, with tests of neutrality and mismatch distribution indicative of a recent population expansion in Labrador, dated between c. 3500 and 8300 years ago.

Genetic diversity and genetic structure of the striped field mouse Apodemus agrarius coreae (Muridae, Rodentia) in Korea

The aim of this study was to investigate the genetic diversity and genetic structure of the striped field mouse Apodemus agrarius coreae in Korea. The Korean A. a. coreae is characterized by high levels of haplotype diversity (Hd=0.967) and low levels of nucleotide diversity (π=0.00683). Haplogroup 1 is well separated from the haplotypes of the neighboring regions of the Korean Peninsula, while the other haplogroups are closely related to those from the Russian Far East. Thus, further investigations are required to confirm the validity of the subspecies status of A. a. coreae by implementing additional morphological characters as well as genetic data from the populations present in the Korean Peninsula and its neighboring countries. Haplogroup 1 includes most Korean haplotypes and forms a star-like haplotype network structure, which reveals relatively low levels of sequence divergence and high frequency of unique mutations (only few mutations are shared in most of the haplotype nodes). The results indicate that the haplotypes of Haplogroup 1 might have experienced population expansion since their migration into Korea, which was further corroborated with negative results of neutrality tests for Korean population of A. a. coreae.

Differences in speciation progress in feather mites (Analgoidea) inhabiting the same host: the case of Zachvatkinia and Alloptes living on arctic and long-tailed skuas

Recent molecular phylogenetic analyses have revealed that some apparently oligoxenous feather mite species are in fact monoxenous cryptic species with little morphological differentiation. In this study we analyzed two species, Zachvatkinia isolata (Avenzoariidae) and Alloptes (Sternalloptes) stercorarii (Alloptidae) which prefer different parts of the plumage of two sister species of birds: arctic skua (Stercorarius parasiticus) and long-tailed skua (S. longicaudus) breeding on tundra in the High Arctic archipelago of Svalbard. Given that there are no reports about hybridization events between the host species, we expected that both skuas would have a species-specific acarofauna. The genetic distances among DNA-barcode sequences (COI and 28S rDNA), phylogenetic tree topologies, and haplotype networks of the COI sequences of mites suggested extensive gene flow in Z. isolata between and within populations inhabiting both skua species, whereas the Alloptes populations were host specific and sufficiently genetically separated as to warrant species-level status. The discrepancy in the genetic structure of Alloptes and Zachvatkinia populations suggests frequent but transient contacts between the two skua species in which the probability of mite exchange is much higher for Zachvatkinia, which is present in high numbers and inhabits exposed parts of primary flight feathers, than for the less abundant Alloptes that lives primarily in more protected and inaccessible parts of the plumage. We discuss the possible nature of these contacts between host species and the area(s) where they might take place. The star-like structures in the haplotype network as well as high haplotype diversity and low nucleotide diversity observed in Z. isolata are concordant with the known dispersal strategy of feather mites: vertical colonization of new host individuals followed by rapid growth of founder populations.

Genetic fragmentation in India's third longest river system, the Narmada

India’s third longest river, the Narmada, is studied here for the potential effects on native fish populations of river fragmentation due to various barriers including dams and a waterfall. The species we studied include a cyprinid fish, Catla catla, and a mastacembelid, Mastacembelus armatus, both of which are found in the Narmada. Our goal was to use DNA sequence information from the D-loop region of the mitochondrial DNA to explore how this fragmentation could impact the genetic structure of these fish populations. Our results clearly show that these barriers can contribute to the fragmentation of the genetic structure of these fish communities, Furthermore, these barriers enhance the effects of natural isolation by distance and the asymmetry of dispersal flows. This may be a slow process, but it can create significant isolation and result in genetic disparity. In particular, populations furthest upstream having low migration rates could be even more subject to genetic impoverishment. This study serves as a first report of its kind for a river system on the Indian subcontinent. The results of this study also emphasize the need for appropriate attention towards the creation of fish passages across the dams and weirs that could help in maintaining biodiversity.

Assembly and annotation of full mitochondrial genomes for the corn rootworm species, Diabrotica virgifera virgifera and Diabrotica barberi (Insecta: Coleoptera: Chrysomelidae), using Next Generation Sequence data

Complete mitochondrial genomes for two corn rootworm species, Diabrotica virgifera virgifera (16,747 bp) and Diabrotica barberi (16,632; Insecta: Coleoptera: Chrysomelidae), were assembled from Illumina HiSeq2000 read data. Annotation indicated that the order and orientation of 13 protein coding genes (PCGs), and 22 tRNA and 2 rRNA sequences were in typical of insect mitochondrial genomes. Non-standard nad4 and cox3 stop codons were composed of single T nucleotides and likely completed by adenylation, and atypical TTT start codons was predicted for both D. v. virgifera and D. barberinad1 genes. The D. v. virgifera and D. barberi haplotypes showed 819 variable nucleotide positions within PCG regions (7.36% divergence), which suggest that speciation may have occurred ~3.68 million years ago assuming a linear rate of short-term substitution. Phylogenetic analyses of Coleopteran MtD genome show clustering based on family level, and may have the capacity to resolve the evolutionary history within this Order of insects.

Conservation genetics and the implication for recovery of the endangered Mitchell's satyr Butterfly, Neonympha mitchellii mitchellii

The modern delineation of taxonomic groups is often aided by analyses of molecular data, which can also help inform conservation biology. Two subspecies of the butterfly Neonympha mitchellii are classified as federally endangered in the United States: Neonympha mitchellii mitchellii, the Mitchell's satyr, and Neonympha mitchellii francisi, the Saint Francis's satyr. The recent discovery of additional disjunct populations of N. mitchellii in the southeastern US could have important implications for both legal protection and management decisions. We elucidated the relationships among 48 individuals representing 5 N. mitchellii populations using 6 molecular markers (5 nuclear and 1 mitochondrial) under a variety of analytical frameworks. Phylogenetic analysis resulted in moderately supported clades that corresponded with the geographic region where samples originated. Clustering analyses identified 3 groups, wherein the 2 named subspecies formed separate clusters. Coalescent analyses indicated evolutionary divergence between N. m. mitchellii and all other populations but weakly supported divergence among N. m. francisi and the recently discovered populations. Hence, the 2 currently accepted subspecies were clearly different from one another, but the recently discovered populations could not be completely distinguished from N. m. francisi or each other. We propose that N. m. mitchellii and N. m. francisi continue to be managed as separate endangered species.

Little gene flow between domestic silkmoth Bombyx mori and its wild relative Bombyx mandarina in Japan, and possible artificial selection on the CAD gene of B. mori

We analyzed PCR-amplified carbamoyl-phosphate synthetase 2, aspartate transcarbamylase, and dihydroorotase (CAD) gene fragments from 146 Bombyx mori native strains and found extremely low levels of DNA polymorphism. Two haplotypes were identified, one of which was predominant. CAD haplotype analysis of 42 samples of Japanese B. mandarina revealed four haplotypes. No common haplotype was shared between the two species and at least five base substitutions were detected. This result was suggestive of low levels of gene flow between the two species. The nucleotide diversity (π) scores of the two samples differed markedly: lower π values were estimated for B. mori native strains than Japanese B. mandarina. We further analyzed 12 Chinese B. mandarina derived from seven areas of China, including Taiwan. The results clearly indicated that the π score was ~80-fold greater in Chinese B. mandarina than in B. mori. The extremely low level of DNA polymorphism in B. mori compared to its wild relatives suggested that the CAD gene itself or its tightly linked regions are possible targets for silkworm domestication.