Mitochondrial DNA markers reveal high genetic diversity and strong genetic differentiation in populations of Dendrolimus kikuchii Matsumura (Lepidoptera: Lasiocampidae)

Cryptic Diversity and Demographic Expansion of Plasmodium knowlesi Malaria Vectors in Malaysia

Although Malaysia is considered free of human malaria, there has been a growing number of Plasmodium knowlesi cases. This alarming trend highlighted the need for our understanding of this parasite and its associated vectors, especially considering the role of genetic diversity in the adaptation and evolution among vectors in endemic areas, which is currently a significant knowledge gap in their fundamental biology. Thus, this study aimed to investigate the genetic diversity of Anopheles balabacensis, Anopheles cracens, Anopheles introlatus, and Anopheles latens-the vectors for P. knowlesi malaria in Malaysia. Based on cytochrome c oxidase 1 (CO1) and internal transcribed spacer 2 (ITS2) markers, the genealogic networks of An. latens showed a separation of the haplotypes between Peninsular Malaysia and Malaysia Borneo, forming two distinct clusters. Additionally, the genetic distances between these clusters were high (2.3-5.2% for CO1) and (2.3-4.7% for ITS2), indicating the likely presence of two distinct species or cryptic species within An. latens. In contrast, no distinct clusters were observed in An. cracens, An. balabacensis, or An. introlatus, implying a lack of pronounced genetic differentiation among their populations. It is worth noting that there were varying levels of polymorphism observed across the different subpopulations, highlighting some levels of genetic variation within these mosquito species. Nevertheless, further analyses revealed that all four species have undergone demographic expansion, suggesting population growth and potential range expansion for these vectors in this region.

Population Genetics Reveals That the Western Tianshan Mountains Populations of Agrilus mali (Coleoptera: Buprestidae) May Have Not been Recently Introduced

Agrilus mali Matsumura is a wood-boring beetle that aggressively attacks species of the genus Malus, that has recently caused serious damage to the wild apple tree M. sieversii (Lebed.) in the western Tianshan Mountains in Xinjiang. It was first detected there in the early 1990s and spread rapidly, being thus considered a regional invasive pest. To explore the possible outbreak mechanism of the local population and characterize the genetic differentiation of A. mali across different regions of China, we used three mitochondrial genes (COI, COII, and CytB) to investigate the genetic diversity and genetic structure of 17 A. mali populations containing 205 individuals collected from five Chinese provinces. Among them, nine populations were from the western Tianshan Mountains. Ultimately, of the 136 pairwise Fst comparisons, 99 showed high genetic differentiation among overall populations, and Tianshan populations exhibited significant differentiation with most of the non-Tianshan populations. Furthermore, A. mali populations represented relatively abundant haplotypes (54 haplotypes). Nine populations from the Tianshan Mountains showed 32 haplotypes (26 of which were unique), displaying relatively high genetic diversity. Additionally, the Mantel test revealed population genetic differentiation among either overall populations or the Tianshan Mountains populations, likely caused by geographical isolation. Phylogenic relationships showed that all populations clustered into three clades, and Tianshan Mountains populations, including CY, occupied one of the three clades. These results suggest that A. mali in the western Tianshan region has possibly been present in the area for a long period, and may not have been introduced recently. Highly frequent gene flows within Tianshan populations are possibly caused by human activities and may enhance the adaptability of A. mali along the western Tianshan Mountains, leading to periodic outbreaks. These findings enhance our understanding of jewel beetle population genetics and provide valuable information for pest management.

Chromosome-Level Genome Assembly Reveals Significant Gene Expansion in the Toll and IMD Signaling Pathways of Dendrolimus kikuchii

A high-quality genome is of significant value when seeking to control forest pests such as Dendrolimus kikuchii, a destructive member of the order Lepidoptera that is widespread in China. Herein, a high quality, chromosome-level reference genome for D. kikuchii based on Nanopore, Pacbio HiFi sequencing and the Hi-C capture system is presented. Overall, a final genome assembly of 705.51 Mb with contig and scaffold N50 values of 20.89 and 24.73 Mb, respectively, was obtained. Of these contigs, 95.89% had unique locations on 29 chromosomes. In silico analysis revealed that the genome contained 15,323 protein-coding genes and 63.44% repetitive sequences. Phylogenetic analyses indicated that D. kikuchii may diverged from the common ancestor of Thaumetopoea. Pityocampa, Thaumetopoea ni, Heliothis virescens, Hyphantria armigera, Spodoptera frugiperda, and Spodoptera litura approximately 122.05 million years ago. Many gene families were expanded in the D. kikuchii genome, particularly those of the Toll and IMD signaling pathway, which included 10 genes in peptidoglycan recognition protein, 19 genes in MODSP, and 11 genes in Toll. The findings from this study will help to elucidate the mechanisms involved in protection of D. kikuchii against foreign substances and pathogens, and may highlight a potential channel to control this pest.

Population genetic structure of the malaria vector Anopheles minimus in Thailand based on mitochondrial DNA markers

BACKGROUND

The malaria vector Anopheles minimus has been influenced by external stresses affecting the survival rate and vectorial capacity of the population. Since An. minimus habitats have continuously undergone ecological changes, this study aimed to determine the population genetic structure and the potential gene flow among the An. minimus populations in Thailand.

METHODS

Anopheles minimus was collected from five malaria transmission areas in Thailand using Centers for Disease Control and Prevention (CDC) light traps. Seventy-nine females from those populations were used as representative samples. The partial mitochondrial cytochrome c oxidase subunit I (COI), cytochrome c oxidase subunit II (COII) and cytochrome b (Cytb) gene sequences were amplified and analyzed to identify species and determine the current population genetic structure. For the past population, we determined the population genetic structure from the 60 deposited COII sequences in GenBank of An. minimus collected from Thailand 20 years ago.

RESULTS

The current populations of An. minimus were genetically divided into two lineages, A and B. Lineage A has high haplotype diversity under gene flow similar to the population in the past. Neutrality tests suggested population expansion of An. minimus, with the detection of abundant rare mutations in all populations, which tend to arise from negative selection.

CONCLUSIONS

This study revealed that the population genetic structure of An. minimus lineage A was similar between the past and present populations, indicating high adaptability of the species. There was substantial gene flow between the eastern and western An. minimus populations without detection of significant gene flow barriers.

From the river to the ocean: mitochondrial DNA analyses provide evidence of spectacled caimans (Caiman crocodilus Linnaeus 1758) mainland–insular dispersal

There are few studies that have investigated the evolutionary history of large vertebrates on islands off the Caribbean coast of South America. Here we use the spectacled caiman (Caiman crocodilus) to investigate among- and within-population patterns of genetic diversity to understand connectivity between island and mainland populations. The spectacled caiman is naturally distributed across Central and South America including the islands of Trinidad and Tobago, which are considered to have the only natural insular populations of the species. Because of this apparent isolation, we sought to determine whether caimans on Trinidad and Tobago comprise a unique lineage and have reduced genetic diversity compared to mainland caimans. We test these hypotheses by using mitochondrial DNA variation to assess the phylogenetic and phylogeographical relationships of the C. crocodilus populations inhabiting these islands within the evolutionary context of the entire spectacled caiman complex. Phylogenetic analyses placed the Trinidad and Tobago samples together with samples from Colombia, Venezuela and Brazil into one well-supported clade, which corresponds to the defined Orinoco/upper Negro lineage. Interestingly, the majority of sequences from Trinidad and Tobago are similar or identical to haplotypes reported from Venezuela and Colombia, supporting the idea of a dispersal process from the Orinoco River to these islands. We discuss the implications of our findings for systematics and the conservation of the species and how these dispersal movements could shape the current phylogeographical structure depicted for C. crocodilus.

High genetic diversity and strong genetic structure of Strongyllodes variegatus populations in oilseed rape production areas of China

Background

Strongyllodes variegatus (Fairmaire) is a major insect pest of oilseed rape in China. Despite its economic importance, the contribution of its population genetics in the development of any suitable protection control strategy for the management of oilseed rape crops is poorly studied. It is a much urgent need to prevent its spread to the rest of the world.

Results

Using the sequences of mitochondrial DNA cytochrome c oxidase subunit I (COI) and cytochrome b (Cytb) as genetic markers, we analyzed the population genetic diversity and structure of 437 individuals collected from 15 S. variegatus populations located in different oilseed rape production areas in China. In addition, we estimated the demographic history using neutrality test and mismatch distribution analysis. The high level of genetic diversity was detected among the COI and Cytb sequences of S. variegatus. The population structure analyses strongly suggested three distinct genetic and geographical regions in China with limited gene flow. The Mantel test showed that the genetic distance was greatly influenced by the geographical distance. The demographic analyses showed that S. variegatus had experienced population fluctuation during the Pleistocene Epoch, which was likely to be related to the climatic changes.

Conclusion

Overall, these results demonstrate that the strong genetic structure of S. variegatus populations in China, which is attributed by the isolation through the geographical distance among populations, their weak flight capacity and subsequent adaptation to the regional ecological conditions.

Mitochondrial DNA analysis reveals spatial genetic structure and high genetic diversity of Massicus raddei (Blessig) (Coleoptera: Cerambycidae) in China

The oak longhorned beetle (OLB), Massicus raddei (Blessig, 1872) (Coleoptera: Cerambycidae), is widely distributed in Asia (China, the Korean Peninsula, Japan, Vietnam and the Russian Far-East), but pest outbreaks have occurred only in Liaoning Province and Jilin Province of China. In order to explore possible mechanisms of local population outbreaks and characterize the genetic diversity and genetic structure of M. raddei across its range in China, three mitochondrial genes (COI, Cytb, and COII) were sequenced and analyzed for seven M. raddei populations collected from six provinces in China. From these different populations, we found a high haplotype and nucleotide diversity. Haplotype networks and phylogenetic analyses both demonstrate apparent genetic diversification between SC (southern China) and NC (northern China) population groups. A set of 21 pairwise comparisons for Fst (pairwise fixation indices) and Nm (genetic flow index) showed significant genetic differentiation and limited gene flow except for two pairs, Shandong (SD) and Liaoning (LN), and Anhui (AH) and Henan (HN). This pattern suggested that the periodic outbreak of the LN population could not be attributed to the absence of genetic flow with other spatial populations and that regional environmental factors might be responsible. AMOVA (Analysis of molecular variance) showed that the greater molecular genetic variation was among populations. Based on Tajima's D statistic, Fu's Fs, and the mismatch distribution test, we determined that the seven populations sampled were stable and had not experienced any recent population expansion. The fact that all the sampled populations showed only unique haplotypes and lacked shared or ancestral haplotypes, as well as the nonstar-like distribution of haplotype network for concatenated genes, collectively provided powerful evidence of the stable and isolated nature of most populations. The high genetic differentiation and spatial genetic structuring among populations are both likely related to the beetle's moderate flight capacity, regional variation in host tree species and microclimate, as well as the geographic distance between sampling sites.

Comparative mitochondrial genome analysis of Dendrolimus houi (Lepidoptera: Lasiocampidae) and phylogenetic relationship among Lasiocampidae species

Dendrolimus houi is one of the most common caterpillars infesting Gymnosperm trees, and widely distributed in several countries in Southeast Asia, and exists soley or coexists with several congeners and some Lasiocampidae species in various forest habitats. However, natural hybrids occasionally occur among some closely related species in the same habitat, and host preference, extreme climate stress, and geographic isolation probably lead to their uncertain taxonomic consensus. The mitochondrial DNA (mtDNA) of D. houi was extracted and sequenced by using high-throughput technology, and the mitogenome composition and characteristics were compared and analyzed of these species, then the phylogenetic relationship was constructed using the maximum likelihood method (ML) and the Bayesian method (BI) based on their 13 protein-coding genes (PCGs) dataset, which were combined and made available to download which were combined and made available to download among global Lasiocampidae species data. Mitogenome of D. houi was 15,373 bp in length, with 37 genes, including 13 PCGs, 22 tRNA genes (tRNAs) and 2 rRNA genes (rRNAs). The positions and sequences of genes were consistent with those of most known Lasiocampidae species. The nucleotide composition was highly A+T biased, accounting for ~80% of the whole mitogenome. All start codons of PCGs belonged to typical start codons ATN except for COI which used CGA, and most stop codons ended with standard TAA or TAG, while COI, COII, ND4 ended with incomplete T. Only tRNASer (AGN) lacked DHU arm, while the remainder formed a typical "clover-shaped" secondary structure. For Lasiocampidae species, their complete mitochondrial genomes ranged from 15,281 to 15,570 bp in length, and all first genes started from trnM in the same direction. And base composition was biased toward A and T. Finally, both two methods (ML and BI) separately revealed that the same phylogenetic relationship of D. spp. as ((((D. punctatus + D. tabulaeformis) + D. spectabilis) + D. superans) + (D. kikuchii of Hunan population + D. houi) as in previous research, but results were different in that D. kikuchii from a Yunnan population was included, indicating that different geographical populations of insects have differentiated. And the phylogenetic relationship among Lasiocampidae species was ((((Dendrolimus) + Kunugia) + Euthrix) + Trabala). This provides a better theoretical basis for Lasiocampidae evolution and classification for future research directions.

Lineage Divergence of Dendrolimus punctatus in Southern China Based on Mitochondrial Genome

In southern China, the masson pine caterpillar, Dendrolimus punctatus, has caused serious damage to the Pinus massoniana (Lamb.) pine forests. Here, the whole mitochondrial DNA (mtDNA) was employed to analyze the population evolution of D. punctatus and to understand the process underlying its current phylogenetic pattern. D. punctatus populations within its distribution range in China were categorized into five subgroups: central and eastern China (CEC), southwestern China (SWC), Yibin in Sichuan (SC), Baise in Guangxi (GX), and Luoding in Guangdong (GD), with a high level of haplotype diversity and nucleotide diversity among them. The genetic distances between subgroups are relatively large; however, the genetic distances between populations within the CEC subgroup were relatively small, suggesting that many populations were closely related in this subgroup. The mantel test showed that geographic distance had an important impact on the genetic distance of different geographic populations (r = 0.3633, P < 0.001). The neutrality tests, Bayesian skyline plot, and haplotype network showed that D. punctatus experienced a population expansion around 100,000 years ago. The divergence times of GX/SC, SWC, GD, and CEC were 0.347, 0.236, 0.200, and 0.110 million years ago, respectively. The SWC, CEC, and GD subgroups might have evolved from GX or SC subgroups. The population genetic structure of D. punctatus was closely related to its host tree species, geographic distance among populations, the weak flight capacity, and many eco-environment conditions.

The complete mitochondrial genome of Dendrolimus kikuchii (Lepidoptera: Lasiocampidae)

As a serious forest pest on coniferous trees, Dendrolimus kikuchii has caused widespread concern in China. Here, its complete mitochondrial genome (mitogenome) has been sequenced with 15,382 bp in length. The mitogenome has a base composition of A (40.87%), T (37.83%), C (13.43%), and G (7.87%), and consists of 13 protein-coding genes (PCGs), 22 transfer RNA (tRNA) genes, two ribosomal RNA (rRNA) genes, and an A + T-rich region. The phylogenetic relationships among the lasiocampid species were (Trabala vishnou+ ((Apatelopteryx phenax+ Euthrix laeta) + (Dendrolimus kikuchii+ (D. spectabilis+ (D. tabulaeformis + D. punctatus))))), which were supported by a posterior probability of 1.00 and a high bootstrap value of 100%.