DnaSP v5: a software for comprehensive analysis of DNA polymorphism data

Phylogenetic and population genetic analyses of Thrips tabaci Lindeman (Thysanoptera: Thripidae) on Allium host in India

Background

Onion thrips (Thrips tabaci) is a complex of cryptic species with subtle morphological differences and distinct genetic backgrounds; thus, species identification using traditional methods remains challenging. The existence of different haplotypes and genotypes within a species can significantly influence various aspects of its biology, including host preference, reproductive capacity, resistance to pesticides, and vector competence for plant viruses. Understanding the genetic diversity and population structure of cryptic species within T. tabaci will not only aid in the development of more effective control strategies tailored to specific genetic variants but also in monitoring population dynamics, tracking invasive species, and implementing quarantine measures to prevent the spread of economically damaging thrips biotypes.

Methods

This study aims to explore intraspecies genetic diversity and molecular evolutionary relationships of the mitochondrial cytochrome oxidase gene subunit I (mtCOI) in T. tabaci populations from India. To capture diversity within the Indian T. tabaci populations, amplicon sequencing was performed for the thrips mtCOI gene from eight diverse localities in India. A total of 48 sequences retrieved for the mtCOI gene from the NCBI Nucleotide database were analysed.

Results

Multiple insertions and deletions were detected at various genomic positions across the populations from different localities, with the highest variation observed in the 300-400 genome position range. Molecular diversity analyses identified 30 haplotypes within the population, with certain subpopulations exhibiting higher gene flow. Analysis of single nucleotide polymorphism patterns within the mtCOI gene across diverse Indian locales revealed significant intrapopulation genetic heterogeneity and its potential repercussions on gene functionality. Elevated F statistics (Fst) values in the northern-western subpopulations suggested high genetic variability, particularly evident in haplotype networks originating mainly from the northern region, notably Delhi. While most populations displayed stable and ancient evolutionary histories, thrips populations from northern, western, and north-eastern regions indicated rapid growth.

The first record of Chironomusnuditarsis Keyl, 1961 from Sevan Lake (Armenia) confirmed by morphology, karyotype and COI gene sequence

Chironomusnuditarsis Keyl, 1961 is recorded from Sevan Lake for the first time. This species is widespread in Europe, the Caucasus, and Siberia. For species identification, we used a comprehensive approach that included morphological, cytogenetic and molecular genetic analyses. Morphological analysis showed a high similarity with the description. Nine chromosome banding sequences ndtA1, ndtA2, ndtB2, ndtC1, ndtD1, ndtE1, ndtF1, ndtG1, and ndtG2 were found. The banding sequences ndtA1, ndtA2, ndtG1, and ndtG2 are species-specific for C.nuditarsis and allow us to accurately distinguish it from the sibling species Ch.curabilis Belyanina, Sigareva et Loginova, 1990. Molecular-genetic analysis of the COI gene sequences has shown low genetic distances of 0.38-0.95% in the sibling species Ch.nuditarsis and Ch.curabilis complex and the insufficiency of using a single COI as a molecular marker for their separation.

Complete plastome genomes of three medicinal heliotropiaceae species: comparative analyses and phylogenetic relationships

BACKGROUND

Heliotropiaceae is a family of the order Boraginales and has over 450 species. The members of the family Heliotropiaceae have been widely reported to be used in traditional medicine Over time, the classification of Heliotropiaceae has remained uncertain and has moved from family to subfamily, or conversely.

RESULTS

In the present study, we sequenced, analyzed, and compared the complete plastomes of Euploca strigosa, Heliotropium arbainense, and Heliotropium longiflorum with the genomes of related taxa. The lengths of the plastomes of E. strigosa, H. arbainense, and H. longiflorum were 155,174 bp, 154,709 bp, and 154,496 bp, respectively. Each plastome consisted of 114 genes: 80 protein-coding genes, 4 ribosomal RNA genes, and 30 transfer RNA genes. The long repeats analysis indicated that reverse, palindromic, complement and forward repeats were all found in the three plastomes. The simple repeats analysis showed that the plastomes of E. strigosa, H. arbainense, and H. longiflorum contained 158, 165, and 151 microsatellites, respectively. The phylogenetic analysis confirmed two major clades in the Boraginales: clade I comprised Boraginaceae, while clade II included Heliotropiaceae, Ehretiaceae, Lennoaceae, and Cordiaceae. Inside the family Heliotropiaceae, E. strigosa is nested within the Heliotropium genus.

CONCLUSIONS

This study expands our knowledge of the evolutionary relationships within Heliotropiaceae and offers useful genetic resources.

Natural variation of domestication-related genes contributed to latitudinal expansion and adaptation in soybean

Soybean is a major source of protein and edible oil worldwide. Originating from the Huang-Huai-Hai region, which has a temperate climate, soybean has adapted to a wide latitudinal gradient across China. However, the genetic mechanisms responsible for the widespread latitudinal adaptation in soybean, as well as the genetic basis, adaptive differentiation, and evolutionary implications of theses natural alleles, are currently lacking in comprehensive understanding. In this study, we examined the genetic variations of fourteen major gene loci controlling flowering and maturity in 103 wild species, 1048 landraces, and 1747 cultivated species. We found that E1, E3, FT2a, J, Tof11, Tof16, and Tof18 were favoured during soybean improvement and selection, which explained 75.5% of the flowering time phenotypic variation. These genetic variation was significantly associated with differences in latitude via the LFMM algorithm. Haplotype network and geographic distribution analysis suggested that gene combinations were associated with flowering time diversity contributed to the expansion of soybean, with more HapA clustering together when soybean moved to latitudes beyond 35°N. The geographical evolution model was developed to accurately predict the suitable planting zone for soybean varieties. Collectively, by integrating knowledge from genomics and haplotype classification, it was revealed that distinct gene combinations improve the adaptation of cultivated soybeans to different latitudes. This study provides insight into the genetic basis underlying the environmental adaptation of soybean accessions, which could contribute to a better understanding of the domestication history of soybean and facilitate soybean climate-smart molecular breeding for various environments.

Mitochondrial Genome Characteristics Reveal Evolution of Acanthopsetta nadeshnyi (Jordan and Starks, 1904) and Phylogenetic Relationships

In the present study, the mitochondrial genomic characteristics of Acanthopsetta nadeshnyi have been reported and have depicted the phylogenetic relationship among Pleuronectidae. Combined with a comparative analysis of 13 PCGs, the TN93 model was used to review the neutral evolution and habitat evolution catalysis of the mitogenome to verify the distancing and purification selectivity of the mitogenome in Pleuronectidae. At the same time, a species differentiation and classification model based on mitogenome analysis data was established. This study is expected to provide a new perspective on the phylogenetic relationship and taxonomic status of A. nadeshnyi and lay a foundation for further exploration of environmental and biological evolutionary mechanisms.

Complete Chloroplast Genome of Megacarpaea megalocarpa and Comparative Analysis with Related Species from Brassicaceae

Megacarpaea megalocarpa, a perennial herbaceous species belonging to the Brassicaceae family, has potential medicinal value. We isolated and characterized the chloroplast (cp) genome of M. megalocarpa and compared it with closely related species. The chloroplast genome displayed a typical quadripartite structure, spanning 154,877 bp, with an overall guanine-cytosine (GC) content of 36.20%. Additionally, this genome contained 129 genes, 105 simple sequence repeats (SSRs), and 48 long repeat sequences. Significantly, the ycf1 gene exhibited a high degree of polymorphism at the small single copy (SSC) region and the inverted repeat a (IRa) boundary. Despite this polymorphism, relative synonymous codon usage (RSCU) values were found to be similar across species, and no large segment rearrangements or inversions were detected. The large single copy (LSC) and SSC regions showed higher sequence variations and nucleotide polymorphisms compared to the IR region. Thirteen distinct hotspot regions were identified as potential molecular markers. Our selection pressure analysis revealed that the protein-coding gene rpl20 is subjected to different selection pressures in various species. Phylogenetic analysis positioned M. megalocarpa within the expanded lineage II of the Brassicaceae family. The estimated divergence time suggests that M. megalocarpa diverged approximately 4.97 million years ago. In summary, this study provides crucial baseline information for the molecular identification, phylogenetic relationships, conservation efforts, and utilization of wild resources in Megacarpaea.

Assembly and comparative analysis of the complete mitochondrial and chloroplast genome of Cyperus stoloniferus (Cyperaceae), a coastal plant possessing saline-alkali tolerance

BACKGROUND

Cyperus stoloniferus is an important species in coastal ecosystems and possesses economic and ecological value. To elucidate the structural characteristics, variation, and evolution of the organelle genome of C. stoloniferus, we sequenced, assembled, and compared its mitochondrial and chloroplast genomes.

RESULTS

We assembled the mitochondrial and chloroplast genomes of C. stoloniferus. The total length of the mitochondrial genome (mtDNA) was 927,413 bp, with a GC content of 40.59%. It consists of two circular DNAs, including 37 protein-coding genes (PCGs), 22 tRNAs, and five rRNAs. The length of the chloroplast genome (cpDNA) was 186,204 bp, containing 93 PCGs, 40 tRNAs, and 8 rRNAs. The mtDNA and cpDNA contained 81 and 129 tandem repeats, respectively, and 346 and 1,170 dispersed repeats, respectively, both of which have 270 simple sequence repeats. The third high-frequency codon (RSCU > 1) in the organellar genome tended to end at A or U, whereas the low-frequency codon (RSCU < 1) tended to end at G or C. The RNA editing sites of the PCGs were relatively few, with only 9 and 23 sites in the mtDNA and cpDNA, respectively. A total of 28 mitochondrial plastid DNAs (MTPTs) in the mtDNA were derived from cpDNA, including three complete trnT-GGU, trnH-GUG, and trnS-GCU. Phylogeny and collinearity indicated that the relationship between C. stoloniferus and C. rotundus are closest. The mitochondrial rns gene exhibited the greatest nucleotide variability, whereas the chloroplast gene with the greatest nucleotide variability was infA. Most PCGs in the organellar genome are negatively selected and highly evolutionarily conserved. Only six mitochondrial genes and two chloroplast genes exhibited Ka/Ks > 1; in particular, atp9, atp6, and rps7 may have undergone potential positive selection.

CONCLUSION

We assembled and validated the mtDNA of C. stoloniferus, which contains a 15,034 bp reverse complementary sequence. The organelle genome sequence of C. stoloniferus provides valuable genomic resources for species identification, evolution, and comparative genomic research in Cyperaceae.

Species Composition of a Small Mammal Community and Prevalence of Echinococcus spp. in the Alpine Pastoral Area of the Eastern Tibetan Plateau

We aimed to investigate the species composition of a small mammal community and the prevalence of Echinococcus spp. in a typical endemic area of the Tibetan Plateau. One pika and five rodent species were identified based on the morphological characteristics of 1278 small mammal specimens collected during 2014-2019. Detection of Echinococcus DNA in tissue samples from small mammal specimens revealed that Ochotona curzoniae (pika, total prevalence: 6.02%, 26/432), Neodon fuscus (5.91%, 38/643), N. leucurus (2.50%, 3/120), and Alexandromys limnophilus (21.74%, 10/46) were infected by both E. multilocularis and E. shiquicus; Cricetulus longicaudatus (16.67%, 1/6) was infected by E. shiquicus; and no infection was detected in N. irene (0/15). Neodon fuscus and O. curzoniae were the two most abundant small mammal species. There was no significant difference in the prevalence of pika and the overall rodent species assemblage (6.26%, 53/846); however, the larger rodent populations suggested that more attention should be paid to their role in the transmission of echinococcosis in the wildlife reservoir, which has long been underestimated. Moreover, although DNA barcoding provides a more efficient method than traditional morphological methods for identifying large numbers of small mammal samples, commonly used barcodes failed to distinguish the three Neodon species in this study. The close genetic relationships between these species suggest the need to develop more powerful molecular taxonomic tools.

Phylogenetic inferences based on distinct molecular markers reveals a novel Babesia (Babesia pantanalensis nov. sp.) and a Hepatozoon americanum-related genotype in crab-eating foxes (Cerdocyon thous)

Piroplasmids and Hepatozoon spp. Are apicomplexan protozoa that may cause disease in several canid species. The present study aimed to expand the knowledge on the diversity of piroplasmids and Hepatozoon in crab-eating foxes (Cerdocyon thous; n = 12) sampled in the Pantanal of Mato Grosso do Sul State, central-western Brazil. PCR assays based on the 18S rRNA were used as screening. Three (25%) and 11 (91.7%) were positive for piroplasmids and Hepatozoon spp., respectively. Co-infection was found in three C. thous. Phylogenetic analyses based on the near-complete 18S rRNA, cox-1 and hsp70 genes evidenced the occurrence of a novel of Babesia spp. (namely Babesia pantanalensis nov. sp.) closely related to Rangelia vitalii and Babesia sp. 'Coco'. This finding was supported by the genetic divergence analysis which showed (i) high divergence, ranging from 4.17 to 5.62% for 18 S rRNA, 6.16% for hps70 and 4.91-9.25% for cox-1 and (ii) the genotype network (which displayed sequences separated from the previously described Piroplasmida species by median vectors and several mutational events). Also, phylogenetic analysis based on the 18S rRNA gene of Hepatozoon spp. positioned the sequences obtained herein in a clade phylogenetically related to Hepatozoon sp. 'Curupira 2', Hepatozoon sp. detected in domestic and wild canids from Uruguay and Hepatozoon americanum. The present study described Babesia pantanalensis nov sp. and Hepatozoon closely related to H. americanum in crab-eating foxes from Brazil. Moreover, the coinfection by piroplasmids and Hepatozoon sp. for the first time in crab-eating foxes strongly suggesting that this wild canid species potentially acts as a bio-accumulate of hemoprotozoan in wild environment.

SNP Array for Small-Shrimp (Genus Acetes) Origin Determination Using Machine Learning

Accurate origin determination of seafood is crucial for consumer trust and safety. This study was performed to develop a machine learning-based single-nucleotide polymorphism (SNP) analysis technique to determine the origin of Acetes species in salted small-shrimp products. Mitochondrial DNA (COI and 16S rRNA) analysis revealed genetic variations among species and origins. Eight candidate SNPs were identified, six of which were developed into markers for genotyping analysis. Using the developed markers, an SNP array was created and SNP data from salted small-shrimp samples were obtained. Machine learning analysis using a supervised learning algorithm achieved 100% accuracy in classifying the origin of Acetes based on SNP data. This method offers a reliable method for regulatory bodies to combat food fraud and ensure product integrity. The approach can be further improved by expanding the data set to encompass a wider range of species and origins. This study highlights the potential of SNP analysis and machine learning for ensuring seafood authenticity and promoting sustainable practices.

An integrative approach highlights the discrepancy in the genetic, phenotypic, and presumptive taxonomic structure of Phoxinus (Actinopterygii, Leuciscidae, Phoxininae) in Bulgaria

The present drainage network of Bulgaria is the result of a complex Neogene and Quaternary evolution. Karst, which has developed on 23% of the territory, further complicates the hydrological pattern. Fresh waters of Bulgaria drain into the Black Sea and the Aegean Sea basins and can be roughly divided into the Danube (Middle and Lower Danube), non-Danube Black Sea, East Aegean, and West Aegean hydrological regions. Phoxinus, a small leuciscid fish, has a mosaic distribution in all four of these regions, inhabiting small mountainous and semi-mountainous streams. Based on morphology, it was identified as three species, Phoxinus phoxinus in the Danube, Phoxinus strandjae in the non-Danube, and Phoxinus strymonicus in West Aegean region. Later, molecular data revealed Phoxinus csikii and Phoxinus lumaireul in the Middle Danube and P. csikii in the Lower Danube. Phoxinus has been the focus of many studies, showing a high molecular and morphological diversity, which is not entirely consistent with previous morphology-only-based taxonomic concepts. In this study, molecular (a mitochondrial marker and a nuclear marker) and morphological data from both historical and recently sampled collections were analysed to assess the applicability of the integrative approach in Phoxinus. The results showed a significant influence of the complex paleo- and recent hydrology on the currently observed genetic structure of the considered populations and species. Furthermore, the study also demonstrated a strong influence of phenotypic plasticity on the morphological analysis of Phoxinus and the lack of a clear differentiation between P. csikii and P. strandjae. A barcoded specimen was designated as neotype to fix the species named P. strandjae in the current taxonomic concept. Finally, a significant discordance between genetically delimited clades and phenotypic groups did not allow a proper delineation of the species distributed in Bulgaria, demonstrating that more molecular markers are needed for further taxonomic study of the Phoxinus complex.

Plastome structure, phylogeny and evolution of plastid genes in Reevesia (Helicteroideae, Malvaceae)

Reevesia is an eastern Asian-eastern North American disjunction genus in the family Malvaceae s.l. and comprises approximately 25 species. The relationships within the genus are not well understood. Here, 15 plastomes representing 12 Reevesia species were compared, with the aim of better understanding the species circumscription and phylogenetic relationships within the genus and among genera in the family Malvaceae s.l. The 11 newly sequenced plastomes range between 161,532 and 161, 945 bp in length. The genomes contain 114 unique genes, 18 of which are duplicated in the inverted repeats (IRs). Gene content of these plastomes is nearly identical. All the protein-coding genes are under purifying selection in the Reevesia plastomes compared. The top ten hypervariable regions, SSRs, and the long repeats identified are potential molecular markers for future population genetic and phylogenetic studies. Phylogenetic analysis based on the whole plastomes confirmed the monophyly of Reevesia and a close relationship with Durio (traditional Bombacaceae) in subfamily Helicteroideae, but not with the morphologically similar genera Pterospermum and Sterculia (both of traditional Sterculiaceae). Phylogenetic relationships within Reevesia suggested that two species, R. pubescens and R. thyrsoidea, as newly defined, are not monophyletic. Six taxa, R. membranacea, R. xuefengensis, R. botingensis, R. lofouensis, R. longipetiolata and R. pycnantha, are suggested to be recognized.

Complete chloroplast genomes of the hemiparasitic genus Cymbaria: Insights into comparative analysis, development of molecular markers, and phylogenetic relationships

The hemiparasitic tribe Cymbarieae (Orobanchaceae) plays a crucial role in elucidating the initial stage of the transition from autotrophism to heterotrophism. However, the complete chloroplast genome of the type genus Cymbaria has yet to be reported. In addition, the traditional Mongolian medicine Cymbaria daurica is frequently subjected to adulteration or substitution because of the minor morphological differences with Cymbaria mongolica. In this study, the complete chloroplast genomes of the two Cymbaria species were assembled and annotated, and those of other published 52 Orobanchaceae species were retrieved for comparative analyses. We found that the Cymbaria chloroplast genomes are characterized by pseudogenization or loss of stress-relevant genes (ndh) and a unique rbcL-matK inversion. Unlike the high variability observed in holoparasites, Cymbaria and other hemiparasites exhibit high similarity to autotrophs in genome size, guanine-cytosine (GC) content, and intact genes. Notably, four pairs of specific DNA barcodes were developed and validated to distinguish the medicinal herb from its adulterants. Phylogenetic analyses revealed that the genus Cymbaria and the Schwalbea-Siphonostegia clade are grouped into the tribe Cymbarieae, which forms a sister clade to the remaining Orobanchaceae parasitic lineages. Moreover, the diversification of monophyletic Cymbaria occurred during the late Miocene (6.72 Mya) in the Mongol-Chinese steppe region. Our findings provide valuable genetic resources for studying the phylogeny of Orobanchaceae and plant parasitism, and genetic tools to validate the authenticity of the traditional Mongolian medicine "Xinba.".

Delimiting species, revealing cryptic diversity in Molytinae (Coleoptera: Curculionidae) weevil through DNA barcoding

The subfamily Molytinae (Coleoptera: Curculionidae), being the second largest group within the family Curculionidae, exhibits a diverse range of hosts and poses a serious threat to agricultural and forestry industries. We used 1,290 cytochrome c oxidase subunit I (COI) barcodes to assess the efficiency of COI barcodes in species differentiation and uncover cryptic species diversity within weevils of Molytinae. The average Kimura 2-parameter distances within species, genus, and subfamily were 2.90%, 11.0%, and 22.26%, respectively, indicating significant genetic differentiation at both levels. Moreover, there exists a considerable degree of overlap between intraspecific (0%-27.50%) and interspecific genetic distances (GDs; 0%-39.30%). The application of Automatic barcode gap discovery, Assemble Species by Automatic Partitioning, Barcode Index Number, Poisson Tree Processes (PTP), Bayesian Poisson Tree Processes (bPTP), and jMOTU resulted in the identification of 279, 275, 494, 322, 320, and 279 molecular operational taxonomic units, respectively. The integration of 6 methods successfully delimited species of Molytinae in 86.6% of all examined morphospecies, surpassing a threshold value of 3% GD (73.0%). A total of 28 morphospecies exhibiting significant intraspecific divergences were assigned to multiple MOTUs, respectively, suggesting the presence of cryptic diversity or population divergence. The identification of cryptic species within certain morphological species in this study necessitates further investigation through comprehensive taxonomic practices in the future.

Natural variation in MORE GRAINS 1 regulates grain number and grain weight in rice

Grain yield is determined mainly by grain number and grain weight. In this study, we identified and characterized MORE GRAINS1 (MOG1), a gene associated with grain number and grain weight in rice (Oryza sativa L.), through map-based cloning. Overexpression of MOG1 increased grain yield by 18.6%-22.3% under field conditions. We determined that MOG1, a bHLH transcription factor, interacts with OsbHLH107 and directly activates the expression of LONELY GUY (LOG), which encodes a cytokinin-activating enzyme and the cell expansion gene EXPANSIN-LIKE1 (EXPLA1), positively regulating grain number per panicle and grain weight. Natural variations in the promoter and coding regions of MOG1 between Hap-LNW and Hap-HNW alleles resulted in changes in MOG1 expression level and transcriptional activation, leading to functional differences. Haplotype analysis revealed that Hap-HNW, which results in a greater number and heavier grains, has undergone strong selection but has been poorly utilized in modern lowland rice breeding. In summary, the MOG1-OsbHLH107 complex activates LOG and EXPLA1 expression to promote cell expansion and division of young panicles through the cytokinin pathway, thereby increasing grain number and grain weight. These findings suggest that Hap-HNW could be used in strategies to breed high-yielding temperate japonica lowland rice.

New Insight into Nucleotide Changes on Irradiated Bactrocera dorsalis (Hendel), A Pest of Horticultural Importance

Bactrocera dorsalis (Hendel) is a major quarantine pest species infesting most of the tropical fruits. Its infestation had significantly reduced and disrupted the export market trade, thus, very crucial to be controlled during the preharvest and postharvest. One of the most sustainable control methods is by using the radiation technique to reduce the pest population, thus curbing the spread of this pest to new geographical areas. The objective of this study was to measure the nucleotide changes in B. dorsalis (larval, pupal and adult stages) which had been irradiated with 50 to 400 Gray, using Gamma Cell Biobeam GM8000 irradiator with Cesium-137 source at the Malaysian Nuclear Agency, Selangor, Malaysia. Data from the treated samples (with and without morphological changes) were analysed using cytochrome oxidase subunit I (COI). The alignment of 59 sequences resulted in 0.92% variables with only four characters that were parsimony informative, and six sites (30, 60, 234, 282, 483 and 589) which had nucleotide changes, but had not been translated to another protein. Low polymorphism was presented on the sample groups, with only four haplotypes, but with high diversity value (Hd) = 0.5885. The phylogeny trees formed soft polytomy in both trees [neighbour joining (NJ) and maximum parsimony (MP)] presenting a mixture of individuals but did not show any significant difference between treatments. This finding concluded that low mutation had occurred on the treated B. dorsalis and this information is very valuable in getting new insight on the survival of B. dorsalis in the horticulture industry.

Phylogeography of Morella nana: The Wumeng Mountains as a natural geographical isolation boundary on the Yunnan-Guizhou Plateau

The Yunnan-Guizhou Plateau (YGP) is characterized by the distinctive isolated habitat of the limestone Karst Islands and features the Wumeng Mountains, which divide the YGP into the two Plateaus of Yunnan and Guizhou. This study aimed to assess the effects of geographic isolation and past climate fluctuation on the distribution of flora in the YGP. To achieve this, we carried out the phylogeographical pattern and genetic structure based on chloroplast and nuclear ribosomal DNA sequence in relation to past (Last Glacial Maximum) and present distributions based on ecological niche modeling for Morella nana, an important wild plant resource and endemic to the YGP once considered a vulnerable species. The results suggest that the genetic and chlorotype network structures of M. nana are divided into at least two groups: cpDNA chlorotype H2 (or dominant nrDNA haplotypes h1 and h2), distributed primarily to the east of the Wumeng Mountains, and cpDNA chlorotypes H1 and H3-H10 (or dominant nrDNA haplotype h2 and h3), distributed to the west of the Wumeng Mountains. A deep genetic split was noted within the two groups to reach 25 steps, especially for the cpDNA fragment variation. This east-west divergence reveals the existence of a natural geographical isolation boundary in the form of the Wumeng Mountains, and supports the existence of at least two glacial refuges during the Quaternary glacial period, along with two genetic diversity center, and at least two large geographic protection units for the important species of M. nana. This study indicates that the phylogeographical pattern of M. nana can be attributed to geographic/environmental isolation caused by the Wumeng Mountains and climate fluctuation during the last glacial maximum, and proposes an effective strategy to protecting this important plant resource.

Origin and evolution of a new tetraploid mangrove species in an intertidal zone

Polyploidy is a major factor in the evolution of plants, yet we know little about the origin and evolution of polyploidy in intertidal species. This study aimed to identify the evolutionary transitions in three true-mangrove species of the genus Acanthus distributed in the Indo-West Pacific region. For this purpose, we took an integrative approach that combined data on morphology, cytology, climatic niche, phylogeny, and biogeography of 493 samples from 42 geographic sites. Our results show that the Acanthus ilicifolius lineage distributed east of the Thai-Malay Peninsula possesses a tetraploid karyotype, which is morphologically distinct from that of the lineage on the west side. The haplotype networks and phylogenetic trees for the chloroplast genome and eight nuclear genes reveal that the tetraploid species has two sub-genomes, one each from A. ilicifolius and A . ebracteatus, the paternal and maternal parents, respectively. Population structure analysis also supports the hybrid speciation history of the new tetraploid species. The two sub-genomes of the tetraploid species diverged from their diploid progenitors during the Pleistocene. Environmental niche models revealed that the tetraploid species not only occupied the near-entire niche space of the diploids, but also expanded into novel environments. Our findings suggest that A. ilicifolius species distributed on the east side of the Thai-Malay Peninsula should be regarded as a new species, A. tetraploideus, which originated from hybridization between A. ilicifolius and A. ebracteatus, followed by chromosome doubling. This is the first report of a true-mangrove allopolyploid species that can reproduce sexually and clonally reproduction, which explains the long-term adaptive potential of the species.

Delimiting species, revealing cryptic diversity, and population divergence in Qinghai-Tibet Plateau weevils through DNA barcoding

The Leptomias group represents one of the most diverse taxonomic group of weevils in the Qinghai-Tibet Plateau and its adjacent areas. Despite the potential of hidden diversity, relatively few comprehensive studies have been conducted on species diversity in this taxonomic group. In this study, we performed DNA barcoding analysis for species of the Leptomias group using a comprehensive DNA barcode dataset that included 476 sequences representing 54 morphospecies. Within the dataset, our laboratory contributed 474 sequences, and 390 sequences were newly generated for this study. The average Kimura 2-parameter distances among morphospecies and genera were 0.76% and 19.15%, respectively. In 94.4% of the species, the minimum interspecific distances exceeded the maximum intraspecific distances, indicating the presence of barcode gaps in most species of Leptomias group. The application of Automatic Barcode Gap Discovery, Assemble Species by Automatic Partitioning, Barcode Index Number, Bayesian Poisson tree processes, jMOTU, and Neighbor-joining tree methods revealed 45, 45, 63, 54, and 55 distinct clusters representing single species, respectively. Additionally, a total of four morphospecies, Leptomias kangmarensis, L. midlineatus, L. siahus, and L. sp.9RL, were found to be assigned to multiple subclade each, indicating the geographical divergences and the presence of cryptic diversity. Our findings of this study demonstrate that Qinghai-Tibet Plateau exhibits a higher species diversity of the Leptomias group, and it is imperative to investigate cryptic species within certain morphospecies using integrative taxonomic approaches in future studies. Moreover, the construction of a DNA barcode reference library presented herein establishes a robust foundational dataset to support forthcoming research on weevil taxonomy, phylogenetics, ecology, and evolution.

Natural polymorphisms in ZMET2 encoding a DNA methyltransferase modulate the number of husk layers in maize

DNA methylation affects agronomic traits and the environmental adaptability of crops, but the natural polymorphisms in DNA methylation-related genes and their contributions to phenotypic variation in maize (Zea mays) remain elusive. Here, we show that a polymorphic 10-bp insertion/deletion variant in the 3'UTR of Zea methyltransferase2 (ZMET2) alters its transcript level and accounts for variation in the number of maize husk layers. ZMET2 encodes a chromomethylase and is required for maintaining genome-wide DNA methylation in the CHG sequence context. Disruption of ZMET2 increased the number of husk layers and resulted in thousands of differentially methylated regions, a proportion of which were also distinguishable in natural ZMET2 alleles. Population genetic analyses indicated that ZMET2 was a target of selection and might play a role in the spread of maize from tropical to temperate regions. Our results provide important insights into the natural variation of ZMET2 that confers both global and locus-specific effects on DNA methylation, which contribute to phenotypic diversity in maize.

Integrating phylogenetic, phylogeographic, and morphometric analyses to reveal cryptic lineages within the genus Asaccus (Reptilia: Squamata: Phyllodactylidae) in Iran

The Middle Eastern endemic genus Asaccus comprises Southwest Asian leaf-toed geckos. To date, this genus includes 19 species of leaf-toed geckos (seven in Arabia and 12 in the Zagros Mountains). Despite a recent study on the taxonomy and phylogeny of Asaccus species in Iran, controversies still remain surrounding the phylogeny and phylogeography of the genus. Here, we used an integrative approach to determine the phylogeny and phylogeography of Asaccus species using two mitochondrial genes (12 S and Cyt b), and one nuclear gene (c-mos). Our results uncovered 22 distinct lineages, demonstrating a significant cryptic diversity that challenges the current morphological classifications of these species. Phylogenetic analyses reinforce the monophyly of the Asaccus group, positioning A. montanus as a basal lineage, which supports a deep evolutionary divergence dating back to the Late Oligocene, approximately 27.94 million years ago. This genetic diversity also highlights the impact of historical climatic and geographical changes on species diversification. The findings advocate for an integrative approach combining both molecular and morphological data to resolve species identities accurately, thereby enhancing conservation strategies to protect these genetically distinct lineages.

Genetic characterization and phylogenetic analysis of the Nigella sativa (black seed) plastome

In this study, the complete plastome sequence of Nigella sativa (black seed), was analyzed for the first time. The plastome spans approximately 154,120 bp, comprising four sections: the Large Single-Copy (LSC) (85,538 bp), the Small Single-Copy (SSC) (17,984 bp), and two Inverted Repeat (IR) regions (25,299 bp). A comparative study of N. sativa's plastome with ten other species from various genera in the Ranunculaceae family reveals substantial structural variations. The contraction of the inverted repeat region in N. sativa influences the boundaries of single-copy regions, resulting in a shorter plastome size than other species. When comparing the plastome of N. sativa with those of its related species, significant divergence is observed, particularly except for N. damascena. Among these, the plastome of A. glaucifolium displays the highest average pairwise sequence divergence (0.2851) with N. sativa, followed by A. raddeana (0.2290) and A. coerulea (0.1222). Furthermore, the study identified 12 distinct hotspot regions characterized by elevated Pi values (> 0.1). These regions include trnH-GUG-psbA, matK-trnQ-UUG, psbK-trnR-UCU, atpF-atpI, rpoB-psbD, ycf3-ndhJ, ndhC-cemA, petA-psaJ, trnN-GUU-ndhF, trnV-GAC-rps12, ycf2-trnI-CAU, and ndhA-ycf1. Approximately, 24 tandem and 48 palindromic and forward repeats were detected in N. sativa plastome. The analysis revealed 32 microsatellites with the majority being mononucleotide repeats. In the N. sativa plastome, phenylalanine had the highest number of codons (1982 codons), while alanine was the least common amino acid with 260 codons. A phylogenetic tree, constructed using protein-coding genes, revealed a distinct monophyletic clade comprising N. sativa and N. damascene, closely aligned with the Cimicifugeae tribe and exhibiting robust support. This plastome provides valuable genetic information for precise species identification, phylogenetic resolution, and evolutionary studies of N. sativa.

Population genetic structure and demographic history of short mackerel, Rastrelliger brachysoma, in the Gulf of Thailand

The short mackerel Rastrelliger brachysoma (Bleeker 1851) is an important fish in the Gulf of Thailand (GoT). The biology of this species has been intensively studied, but its genetic diversity is little known. The genetic diversity, population genetic structure, and demographic history of this species in the GoT were studied using complete mt control region sequences. The CR sequences of 455 mackerel samples collected from 23 localities at four fishing grounds revealed 333 haplotypes with haplotype diversity (h) per population, ranging between 0.8933 and 1.000, with an average of 0.9781. In turn, the nucleotide diversity (µ) ranged between 0.0119 ± 0.0060 and 0.0333 ± 0.0174, with an average of 0.0220 ± 0.00059.A haplotype network analysis showed that all sequences segregated into two subgroups named, clade I and clade II. Two clades were separated by 26 mutational steps. Each clade formed star-like clusters with many haplotypes derived from a common haplotype. Moreover, an analysis of molecular variance (AMOVA) revealed no significant differences among the studied localities, suggesting the presence of a single population in the GoT. Pairwise differences between samples from different fishing regions also indicated no population structure. Both Tajima's D and Fu's FS statistics were highly significant for the two clades but nonsignificant for the entire population according to a mismatch distribution analysis. These results confirmed that both clades experienced demographic expansion. The estimated expansion times for clade I and clade II were 1,542.307 years (1.5423 ka BP) and 7,602.541 (7.6025 ka BP) years, respectively.

A Mysterious Asian Firefly Genus, Oculogryphus Jeng, Engel & Yang (Coleoptera, Lampyridae): The First Complete Mitochondrial Genome and Its Phylogenetic Implications

The firefly genus Oculogryphus Jeng, Engel & Yang, 2007 is a rare-species group endemic to Asia. Since its establishment, its position has been controversial but never rigorously tested. To address this perplexing issue, we are the first to present the complete mitochondrial sequence of Oculogryphus, using the material of O. chenghoiyanae Yiu & Jeng, 2018 determined through a comprehensive morphological identification. Our analyses demonstrate that its mitogenome exhibits similar characteristics to that of Stenocladius, including a rearranged gene order between trnC and trnW, and a long intergenic spacer (702 bp) between the two rearranged genes, within which six remnants (29 bp) of trnW were identified. Further, we incorporated this sequence into phylogenetic analyses of Lampyridae based on different molecular markers and datasets using ML and BI analyses. The results consistently place Oculogryphus within the same clade as Stenocladius in all topologies, and the gene rearrangement is a synapomorphy for this clade. It suggests that Oculogryphus should be classified together with Stenocladius in the subfamily Ototretinae at the moment. This study provides molecular evidence confirming the close relationship between Oculogryphus and Stenocladius and discovers a new phylogenetic marker helpful in clarifying the monophyly of Ototretinae, which also sheds a new light on firefly evolution.

DNA barcoding reveals cryptic species in the sea slater Ligia italica (Crustacea, Isopoda) from Tunisia

Barcoding studies have provided significant insights into phylogenetic relationships among species belonging to the genus Ligia (Crustacea, Isopoda). Herein the diversity of the Italian sea slater Ligia italica from Tunisia is studied for the first time. Samples were collected from 18 localities in Tunisia, and the analysis included previously published sequences from Italy and Greece available in GenBank. Bayesian and Maximum Likelihood phylogenetic analyses were carried out using a fragment of the mitochondrial COI gene. Putative cryptic species were explored using the 'barcode gap' approach in the software ASAP. A genetic landscape shape analysis was carried out using the program Alleles in Space. The analyses revealed highly divergent and well-supported clades of L. italica dispersed across Tunisia (Clades A1 and A2), Greece (Clade B) and Italy (Clades C1 and C2). High genetic dissimilarity among clades suggested that L. italica constitute a cryptic species complex. Divergence among different L. italica lineages (Clades A, B and C) occurred around 7-4.5 Ma. The detected high genetic distances among clades did not result from atypical mitochondrial DNAs or intracellular infection by Wolbachia bacteria. The complex history of the Mediterranean Sea appears to have played a significant role in shaping the phylogeographic pattern of Ligia italica. Additional morphological and molecular studies are needed to confirm the existence of cryptic species in Ligia italica in Mediterranean.

Mitochondrial DNA and Y chromosome reveal the genetic structure of the native Polish Konik horse population

Polish Konik remains one of the most important horse breeds in Poland. The primitive, native horses with a stocky body and mouse-like coat color are protected by a conservation program, while their Polish population consists of about 3,480 individuals, representing 16 dam and six sire lines. To define the population's genetic structure, mitochondrial DNA and Y chromosome sequence variables were identified. The mtDNA whole hypervariable region analysis was carried out using the Sanger sequencing method on 233 Polish Koniks belonging to all dam lines, while the Y chromosome analysis was performed with the competitive allele-specific PCR genotyping method on 36 horses belonging to all sire lines. The analysis of the mtDNA hypervariable region detected 47 SNPs, which assigned all tested horses to 43 haplotypes. Most dam lines presented more than one haplotype; however, five dam lines were represented by only one haplotype. The haplotypes were classified into six (A, B, E, J, G, R) recognized mtDNA haplogroups, with most horses belonging to haplogroup A, common among Asian horse populations. Y chromosome analysis allocated Polish Koniks in the Crown group, condensing all modern horse breeds, and divided them into three haplotypes clustering with coldblood breeds (28 horses), warmblood breeds (two horses), and Duelmener Pony (six horses). The clustering of all Wicek sire line stallions with Duelmener horses may suggest a historical relationship between the breeds. Additionally, both mtDNA and Y chromosome sequence variability results indicate crossbreeding before the studbooks closure or irregularities in the pedigrees occurred before the DNA testing introduction.

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

OBJECTIVE

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

METHODS

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

RESULTS

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

The chloroplast genome sequences of Ipomoea alba and I. obscura (Convolvulaceae): genome comparison and phylogenetic analysis

Ipomoea species have diverse uses as ornamentals, food, and medicine. However, their genomic information is limited; I. alba and I. obscura were sequenced and assembled. Their chloroplast genomes were 161,353 bp and 159,691 bp, respectively. Both genomes exhibited a quadripartite structure, consisting of a pair of inverted repeat (IR) regions, which are separated by the large single-copy (LSC) and small single-copy (SSC) regions. The overall GC content was 37.5% for both genomes. A total of 104 and 93 simple sequence repeats, 50 large repeats, and 30 and 22 short tandem repeats were identified in the two chloroplast genomes, respectively. G and T were more preferred than C and A at the third base position based on the Parity Rule 2 plot analysis, and the neutrality plot revealed correlation coefficients of 0.126 and 0.105, indicating the influence of natural selection in shaping the codon usage bias in most protein-coding genes (CDS). Genome comparative analyses using 31 selected Ipomoea taxa from Thailand showed that their chloroplast genomes are rather conserved, but the presence of expansion or contraction of the IR region was identified in some of these Ipomoea taxa. A total of five highly divergent regions were identified, including the CDS genes accD, ndhA, and ndhF, as well as the intergenic spacer regions psbI-atpA and rpl32-ccsA. Phylogenetic analysis based on both the complete chloroplast genome sequence and CDS datasets of 31 Ipomoea taxa showed that I. alba is resolved as a group member for series (ser.) Quamoclit, which contains seven other taxa, including I. hederacea, I. imperati, I. indica, I. nil, I. purpurea, I. quamoclit, and I. × sloteri, while I. obscura is grouped with I. tiliifolia, both of which are under ser. Obscura, and is closely related to I. biflora of ser. Pes-tigridis. Divergence time estimation using the complete chloroplast genome sequence dataset indicated that the mean age of the divergence for Ipomoeeae, Argyreiinae, and Astripomoeinae, was approximately 29.99 Mya, 19.81 Mya, and 13.40 Mya, respectively. The node indicating the divergence of I. alba from the other members of Ipomoea was around 10.06 Mya, and the split between I. obscura and I. tiliifolia is thought to have happened around 17.13 Mya. The split between the I. obscura accessions from Thailand and Taiwan is thought to have taken place around 0.86 Mya.

The chloroplast genomes of two medicinal species (Veronica anagallis-aquatica L. and Veronica undulata Wall.) and its comparative analysis with related Veronica species

Veronica anagallis-aquatica L. and Veronica undulata Wall. are widely used ethnomedicinal plants in China. The two species have different clinical efficacies, while their extremely similar morphology and unclear interspecific relationship make it difficult to accurately identify them, leading to increased instances of mixed usage. This article reports on the complete chloroplast genomes sequence of these two species and their related Veronica species to conduct a comparative genomics analysis and phylogenetic construction. The results showed that the chloroplast (cp) genomes of Veronica exhibited typical circular quadripartite structures, with total lengths of 149,386 to 152,319 base pairs (bp), and GC content of 37.9 to 38.1%, and the number of genes was between 129-134. The total number of simple sequence repeats (SSRs) in V. anagallis-aquatica and V. undulata is 37 and 36, while V. arvensis had the highest total number of 56, predominantly characterized by A/T single bases. The vast majority of long repeat sequence types are forward repeats and palindromic repeats. Selective Ka/Ks values showed that three genes were under positive selection. Sequence differences often occur in the non-coding regions of the large single-copy region (LSC) and small single-copy region (SSC), with the lowest sequence variation in the inverted repeat regions (IR). Seven highly variable regions (trnT-GGU-psbD, rps8-rpl16, trnQ-UUG, trnN-GUU-ndhF, petL, ycf3, and ycf1) were detected, which may be potential molecular markers for identifying V. anagallis-aquatica and V. undulata. The phylogenetic tree indicates that there is a close genetic relationship between the genera Veronica and Neopicrorhiza, and V. anagallis-aquatica and V. undulata are sister groups. The molecular clock analysis results indicate that the divergence time of Veronica may occur at ∼ 9.09 Ma, and the divergence time of these two species occurs at ∼ 0.48 Ma. It is speculated that climate change may be the cause of Veronica species diversity and promote the radiation of the genus. The chloroplast genome data of nine Veronica specie provides important insights into the characteristics and evolution of the chloroplast genome of this genus, as well as the phylogenetic relationships of the genus Veronica.

Genetic Variability in Leishmaniasis-Causing Leishmania infantum in Humans and Dogs from North-East Spain

Leishmania infantum is the primary cause of visceral and cutaneous leishmaniasis in the European Mediterranean region. Subspecies-level characterization of L. infantum aids epidemiological studies by offering insights into the evolution and geographical distribution of the parasite and reservoir identity. In this study, conducted in north-east Spain, 26 DNA samples of L. infantum were analyzed, comprising 21 from 10 humans and 5 from 5 dogs. Minicircle kinetoplast DNA (kDNA) polymerase chain reaction assays using primers MC1 and MC2, followed by sequencing, were employed to assess intraspecific genetic variability. Single-nucleotide polymorphism (SNP) analysis detected seven genotypes (G1, G2, G12*-G15*, and G17*), with five being reported for the first time (*). The most prevalent was the newly described G13 (54%), while the other currently identified genotypes were predominantly found in single samples. The in silico restriction fragment length polymorphism (RFLP) method revealed five genotypes (B, F, N, P, and W), one of them previously unreported (W). Genotype B was the most prevalent (85%), comprising three SNP genotypes (G1, G2, and G13), whereas the other RFLP genotypes were associated with single SNP genotypes. These kDNA genotyping methods revealed significant intraspecific genetic diversity in L. infantum, demonstrating their suitability for fingerprinting and strain monitoring.

Complete Organelle Genome of the Desiccation-Tolerant (DT) Moss Tortula atrovirens and Comparative Analysis of the Pottiaceae Family

Tortula atrovirens (Sm.) Lindb. is an important component of biological soil crusts and possesses an extraordinary tolerance against desiccation in dryland habitats. However, knowledge of the organelle genome of this desiccation-tolerant (DT) moss is still lacking. Here, we assembled the first reported Tortula organelle genome and conducted a comprehensive analysis within the Pottiaceae family. T. atrovirens exhibited the second largest chloroplast genome (129,646 bp) within the Pottiaceae, whereas its mitogenome (105,877 bp) and those of other mosses were smaller in size compared to other land plants. The chloroplast and mitochondrial genomes of T. atrovirens were characterized by the expansion of IR boundaries and the absence of homologous recombination-mediated by large repeats. A total of 57 RNA editing sites were detected through mapping RNA-seq data. Moreover, the gene content and order were highly conserved among the Pottiaceae organelle genomes. Phylogenetic analysis showed that bryophytes are paraphyletic, with their three lineages (hornworts, mosses, and liverworts) and vascular plants forming successive sister clades. Timmiella anomala is clearly separated from the monophyletic Pottiaceae, and T. atrovirens is closely related to Syntrichia filaris within the Pottioideae. In addition, we detected four hypervariable regions for candidate-molecular markers. Our findings provide valuable insights into the organelle genomes of T. atrovirens and the evolutionary relationships within the Pottiaceae family, facilitating future discovery of DT genetic resources from bryophytes.

Structure, gene composition, divergence time and phylogeny analysis of the woody desert species Neltuma alba, Neltuma chilensis and Strombocarpa strombulifera

Neltuma alba (Algarrobo blanco), Neltuma chilensis (Algarrobo Chileno) and Strombocarpa strombulifera (Fortuna) are some of the few drought resistant trees and shrubs found in small highly fragmented populations, throughout the Atacama Desert. We reconstructed their plastid genomes using de novo assembly of paired-end reads from total genomic DNA. We found that the complete plastid genomes of N. alba and N. chilensis are larger in size compared to species of the Strombocarpa genus. The Strombocarpa species presented slightly more GC content than the Neltuma species. Therefore, we assume that Strombocarpa species have been exposed to stronger natural selection than Neltuma species. We observed high variation values in the number of cpSSRs (chloroplast simple sequence repeats) and repeated elements among Neltuma and Strombocarpa species. The p-distance results showed a low evolutionary divergence within the genus Neltuma, whereas a high evolutionary divergence was observed between Strombocarpa species. The molecular divergence time found in Neltuma and Strombocarpa show that these genera diverged in the late Oligocene. With this study we provide valuable information about tree species that provide important ecosystem services in hostile environments which can be used to determine these species in the geographically isolated communities, and keep the highly fragmented populations genetically healthy.

Diversity of major histocompatibility complex of II B gene and mate choice in a monogamous and long-lived seabird, the Little Auk (Alle alle)

The major histocompatibility complex (MHC) plays a key role in the adaptive immune system of vertebrates, and is known to influence mate choice in many species. In birds, the MHC has been extensively examined but mainly in galliforms and passerines while other taxa that represent specific ecological and evolutionary life-histories, like seabirds, are underexamined. Here, we characterized diversity of MHC Class II B exon 2 in a colonial pelagic seabird, the Little Auk (or Dovekie Alle alle). We further examined whether MHC variation could be maintained through balancing selection and disassortative mating. We found high polymorphism at the genotyped MHC fragment, characterizing 99 distinct alleles across 140 individuals from three populations. The alleles frequencies exhibited a similar skewed distribution in both sexes, with the four most commonly occurring alleles representing approximately 35% of allelic variation. The results of a Bayesian site-by-site selection analysis suggest evidence of balancing selection and no direct evidence for MHC-dependent disassortative mating preferences in the Little Auk. The latter result might be attributed to the high overall polymorphism of the examined fragment, which itself may be maintained by the large population size of the species.

The complete mitochondrial genomes of five critical phytopathogenic Bipolaris species: features, evolution, and phylogeny

In the present study, three mitogenomes from the Bipolaris genus (Bipolaris maydis, B. zeicola, and B. oryzae) were assembled and compared with the other two reported Bipolaris mitogenomes (B. oryzae and B. sorokiniana). The five mitogenomes were all circular DNA molecules, with lengths ranging from 106,403 bp to 135,790 bp. The mitogenomes of the five Bipolaris species mainly comprised the same set of 13 core protein-coding genes (PCGs), two rRNAs, and a certain number of tRNAs and unidentified open reading frames (ORFs). The PCG length, AT skew and GC skew showed large variability among the 13 PCGs in the five mitogenomes. Across the 13 core PCGs tested, nad6 had the least genetic distance among the 16 Pleosporales species we investigated, indicating that this gene was highly conserved. In addition, the Ka/Ks values for all 12 core PCGs (excluding rps3) were < 1, suggesting that these genes were subject to purifying selection. Comparative mitogenomic analyses indicate that introns were the main factor contributing to the size variation of Bipolaris mitogenomes. The introns of the cox1 gene experienced frequent gain/loss events in Pleosporales species. The gene arrangement and collinearity in the mitogenomes of the five Bipolaris species were almost highly conserved within the genus. Phylogenetic analysis based on combined mitochondrial gene datasets showed that the five Bipolaris species formed well-supported topologies. This study is the first report on the mitogenomes of B. maydis and B. zeicola, as well as the first comparison of mitogenomes among Bipolaris species. The findings of this study will further advance investigations into the population genetics, evolution, and genomics of Bipolaris species.

Plastid phylogenomics contributes to the taxonomic revision of taxa within the genus Sanicula L. and acceptance of two new members of the genus

Introduction

The genus Sanicula L. is a taxonomically complicated taxa within Apiaceae, as its high variability in morphology. Although taxonomists have performed several taxonomic revisions for this genus, the interspecific relationships and species boundaries have not been satisfactorily resolved, especially for those endemic to China. This study mainly focused on S. giraldii var. ovicalycina, S. tienmuensis var. pauciflora, and S. orthacantha var. stolonifera and also described two new members of the genus.

Methods

We newly sequenced sixteen plastomes from nine Sanicula species. Combined with eleven plastomes previously reported by us and one plastome downloaded, we performed a comprehensively plastid phylogenomics analysis of 21 Sanicula taxa.

Results and Discussion

The comparative results showed that 21 Sanicula plastomes in their structure and features were highly conserved and further justified that two new species were indeed members of Sanicula. Nevertheless, eleven mutation hotspot regions were still identified. Phylogenetic analyses based on plastome data and the ITS sequences strongly supported that these three varieties were clearly distant from three type varieties. The results implied that these three varieties should be considered as three independent species, which were further justified by their multiple morphological characters. Therefore, revising these three varieties into three independent species was reasonable and convincing. Moreover, we also identified and described two new Sanicula species (S. hanyuanensis and S. langaoensis) from Sichuan and Shanxi, China, respectively. Based on their distinct morphological characteristics and molecular phylogenetic analysis, two new species were included in Sanicula. In summary, our study impelled the revisions of Sanicula members and improved the taxonomic system of the genus.

Chloroplast Genomes of Vitis flexuosa and Vitis amurensis: Molecular Structure, Phylogenetic, and Comparative Analyses for Wild Plant Conservation

The chloroplast genome plays a crucial role in elucidating genetic diversity and phylogenetic relationships. Vitis vinifera L. (grapevine) is an economically important species, prompting exploration of wild genetic resources to enhance stress resilience. We meticulously assembled the chloroplast genomes of two Korean Vitis L. species, V. flexuosa Thunb. and V. amurensis Rupr., contributing valuable data to the Korea Crop Wild Relatives inventory. Through exhaustive specimen collection spanning diverse ecological niches across South Korea, we ensured comprehensive representation of genetic diversity. Our analysis, which included rigorous codon usage bias assessment and repeat analysis, provides valuable insights into amino acid preferences and facilitates the identification of potential molecular markers. The assembled chloroplast genomes were subjected to meticulous annotation, revealing divergence hotspots enriched with nucleotide diversity, thereby presenting promising candidates for DNA barcodes. Additionally, phylogenetic analysis reaffirmed intra-genus relationships and identified related crops, shedding light on evolutionary patterns within the genus. Comparative examination with chloroplast genomes of other crops uncovered conserved sequences and variable regions, offering critical insights into genetic evolution and adaptation. Our study advances the understanding of chloroplast genomes, genetic diversity, and phylogenetic relationships within Vitis species, thereby laying a foundation for enhancing grapevine genetic diversity and resilience to environmental challenges.

Genetic diversity and population structure of Physella acuta (Gastropoda: Physidae) in Thailand using mitochondrial gene markers: COI and 16S rDNA

Physella acuta is a freshwater snail native to North America. Understanding the phylogeography and genetic structure of P. acuta will help elucidate its evolution. In this study, we used mitochondrial (COI and 16S rDNA) and nuclear (ITS1) markers to identify the species and examine its genetic diversity, population structure, and demographic history of P. acuta in Thailand. Phylogenetic and network analyses of P. acuta in Thailand pertained to clade A, which exhibits a global distribution. Analysis of the genetic structure of the population revealed that the majority of pairwise comparisons showed no genetic dissimilarity. An isolation-by-distance test indicates no significant correlation between genetic and geographical distances among P. acuta populations, suggesting that gene flow is not restricted by distance. Demographic history and haplotype network analyses suggest a population expansion of P. acuta, as evidenced by the star-like structure detected in the median-joining network. Based on these results, we concluded that P. acuta in Thailand showed gene flow and recent population expansion. Our findings provide fundamental insights into the genetic variation of P. acuta in Thailand.

Association between Variants of the TRPV1 Gene and Body Composition in Sub-Saharan Africans

In humans, the transient receptor potential vanilloid 1 (TRPV1) gene is activated by exogenous (e.g., high temperatures, irritating compounds such as capsaicin) and endogenous (e.g., endocannabinoids, inflammatory factors, fatty acid metabolites, low pH) stimuli. It has been shown to be involved in several processes including nociception, thermosensation, and energy homeostasis. In this study, we investigated the association between TRPV1 gene variants, sensory perception (to capsaicin and PROP), and body composition (BMI and bioimpedance variables) in human populations. By comparing sequences deposited in worldwide databases, we identified two haplotype blocks (herein referred to as H1 and H2) that show strong stabilizing selection signals (MAF approaching 0.50, Tajima's D > +4.5) only in individuals with sub-Saharan African ancestry. We therefore studied the genetic variants of these two regions in 46 volunteers of sub-Saharan descent and 45 Italian volunteers (both sexes). Linear regression analyses showed significant associations between TRPV1 diplotypes and body composition, but not with capsaicin perception. Specifically, in African women carrying the H1-b and H2-b haplotypes, a higher percentage of fat mass and lower extracellular fluid retention was observed, whereas no significant association was found in men. Our results suggest the possible action of sex-driven balancing selection at the non-coding sequences of the TRPV1 gene, with adaptive effects on water balance and lipid deposition.

The Analysis, Description, and Examination of the Maize LAC Gene Family's Reaction to Abiotic and Biotic Stress

Laccase (LAC) is a diverse group of genes found throughout the plant genome essential for plant growth and the response to stress by converting monolignin into intricate lignin formations. However, a comprehensive investigation of maize laccase has not yet been documented. A bioinformatics approach was utilized in this research to conduct a thorough examination of maize (Zea mays L.), resulting in the identification and categorization of 22 laccase genes (ZmLAC) into six subfamilies. The gene structure and motifs of each subgroup were largely consistent. The distribution of the 22 LAC genes was uneven among the maize chromosomes, with the exception of chromosome 9. The differentiation of the genes was based on fragment replication, and the differentiation time was about 33.37 million years ago. ZmLAC proteins are primarily acidic proteins. There are 18 cis-acting elements in the promoter sequences of the maize LAC gene family associated with growth and development, stress, hormones, light response, and stress response. The analysis of tissue-specific expression revealed a high expression of the maize LAC gene family prior to the V9 stage, with minimal expression at post-V9. Upon reviewing the RNA-seq information from the publicly available transcriptome, it was discovered that ZmLAC5, ZmLAC10, and ZmLAC17 exhibited significant expression levels when exposed to various biotic and abiotic stress factors, suggesting their crucial involvement in stress responses and potential value for further research. This study offers an understanding of the functions of the LAC genes in maize's response to biotic and abiotic stress, along with a theoretical basis for comprehending the molecular processes at play.

Hapsolutely: a user-friendly tool integrating haplotype phasing, network construction, and haploweb calculation

Motivation

Haplotype networks are a routine approach to visualize relationships among alleles. Such visual analysis of single-locus data is still of importance, especially in species diagnosis and delimitation, where a limited amount of sequence data usually are available and sufficient, along with other datasets in the framework of integrative taxonomy. In diploid organisms, this often requires separating (phasing) sequences with heterozygotic positions, and typically separate programs are required for phasing, reformatting of input files, and haplotype network construction. We therefore developed Hapsolutely, a user-friendly program with an ergonomic graphical user interface that integrates haplotype phasing from single-locus sequences with five approaches for network/genealogy reconstruction.

Results

Among the novel options implemented, Hapsolutely integrates phasing and graphical reconstruction steps of haplotype networks, supports input of species partition data in the common SPART and SPART-XML formats, and calculates and visualizes haplowebs and fields for recombination, thus allowing graphical comparison of allele distribution and allele sharing among subsets for the purpose of species delimitation. The new tool has been specifically developed with a focus on the workflow in alpha-taxonomy, where exploring fields for recombination across alternative species partitions may help species delimitation.

Availability and implementation

Hapsolutely is written in Python, and integrates code from Phase, SeqPHASE, and PopART in C++ and Haxe. Compiled stand-alone executables for MS Windows and Mac OS along with a detailed manual can be downloaded from https://www.itaxotools.org; the source code is openly available on GitHub (https://github.com/iTaxoTools/Hapsolutely).

Unveiling the genetic structure of pig population in a Himalayan state Uttarakhand through microsatellite and mitochondrial DNA analyses

This study traced the maternal lineage of the domestic swine populations using mitochondrial DNA control region markers and genetic diversity using microsatellite markers in Uttarakhand, an Indian state situated at the foothills of the world's youngest (geo-dynamically sensitive) mountain system, "the Himalayas". Analysis of 68 maternally unrelated individuals revealed 20 haplotypes. The maternal signature of the Pacific, Southeast Asian, European, and ubiquitously distributed Chinese haplotypes was present in Uttarakhand's domestic pig population. The D3 haplotype reported in wild pigs from North India was also identified in 47 domestic samples. A unique gene pool, UKD (Uttarakhand Domestic), as another lineage specific to this region has been proposed. Genotypes were analyzed, using 13 sets of microsatellite markers. The observed (Ho) and expected (He) heterozygosities were 0.83 ± 0.02 and 0.84 ± 0.01, respectively. The average polymorphic information content value of 0.83 ± 0.01 indicated the high informativeness of the marker. The overall mean FIS value for all the microsatellite markers was low (F = 0.04, P < 0.01). Seven loci deviated from Hardy-Weinberg equilibrium (HWE) at a significant level (p < 0.05). Two clusters were identified, indicating overlapping populations. These results suggested that though belonging to different maternal lineages, the traditional management practices in Uttarakhand have allowed for genetic mixing and the sharing of genetic material among pig populations. It could contribute to increased genetic diversity but might also result in the loss of distinct genetic characteristics or breed purity of the local breeds if not carefully managed.

Comparative Mitogenomics and Phylogenetic Implications for Nine Species of the Subfamily Meconematinae (Orthoptera: Tettigoniidae)

Currently, the subfamily Meconematinae encompasses 1029 species, but whole-mitochondrial-genome assemblies have only been made available for 13. In this study, the whole mitochondrial genomes (mitogenomes) of nine additional species in the subfamily Meconematinae were sequenced. The size ranged from 15,627 bp to 17,461 bp, indicating double-stranded circular structures. The length of the control region was the main cause of the difference in mitochondrial genome length among the nine species. All the mitogenomes including 13 protein-coding genes (PCGs), 22 transfer RNA genes (tRNAs), 2 ribosomal RNA genes (rRNAs) and a control region (CR). The majority strand encoded 23 genes, and the minority strand encoded 14 genes. A phylogenetic analysis reaffirmed the monophyletic status of each subfamily, but the monophysitism of Xizicus, Xiphidiopsis and Phlugiolopsis was not supported.

"Point by point" source: The Chinese pine plantations in North China by evidence from mtDNA

The geographical variation and domestication of tree species are an important part of the theory of forest introduction, and the tracing of the germplasm is the theoretical basis for the establishment of high-quality plantations. Chinese pine (Pinus tabuliformis Carr.) is an important native timber tree species widely distributed in northern China, but it is unclear exactly where germplasm of the main Chinese pine plantation populations originated. Here, using two mtDNA markers, we analyzed 796 individuals representing 35 populations (matR marker), and 873 individuals representing 38 populations (nad5-1 marker) of the major natural and artificial populations in northern China, respectively (Shanxi, Hebei and Liaoning provinces). The results confirmed that the core position of natural SX* populations ("*" means natural population) in the Chinese pine populations of northern China, the genetic diversity of HB and LN plantations was higher than that of natural SX* populations, and there was a large difference in genetic background within the groups of SX* and LN, HB showed the opposite. More importantly, we completed the "point by point" tracing of the HB and LN plantings. The results indicated that almost all HB populations originated from SX* (GDS*, ZTS*, GCS*, and THS*), which resulted in homogeneity of the genetic background of HB populations. Most of germplasm of the LN plantations originated from LN* (ZJS* and WF*), and the other part originated from GDS* (SX*), resulting in the large differences in the genetic background within the LN group. Our results provided a reliable theoretical basis for the scientific allocation, management, and utilization of Chinese pine populations in northern China, and for promoting the high-quality establishment of Chinese pine plantations.

Eutrichophilus cordiceps Mjöberg, 1910 (Ischnocera: Trichodectidae) in Spiny Tree Porcupines (Coendou villosus): New locality records and the first molecular evidence of association with Bartonella sp

The chewing louse genus Eutrichophilus Mjöberg has 19 species only associated with porcupines (Rodentia: Erethizontidae). Of these species, E. cercolabes, E. cordiceps, E. emersoni, E. minor, E. moojeni, and E. paraguayensis have been recorded in Brazil. In the present study, we report E. cordiceps for the first time in the São Paulo State (Bauru Municipality) and for the second time in the Santa Catarina State (Lages Municipality), providing scanning electron images and light microscopy for the eggs, as well as the first molecular data (18S rRNA) for the genus. Additionally, Bartonella sp. was detected for the first time in this chewing lice species.

Comparative mitogenome research revealed the phylogenetics and evolution of the superfamily Tenebrionoidea (Coleoptera: Polyphage)

Despite the worldwide distribution and rich diversity of the superfamily Tenebrionoidea, the knowledge of the mitochondrial genomes (mtgenome) characteristics of the superfamily is still very limited, and its phylogenetics and evolution remain unresolved. In the present study, we newly sequenced mtgenomes from 19 species belonging to Tenebrionoidea, and a total of 90 mitochondrial genomes from 16 families of Tenebrionoidea were used for phylogenetic analysis. There exist 37 genes for all 82 species of complete mtgenomes of 16 families investigated, and their characteristics are identical as reported mtgenomes of other Tenebrionoids. The Ka/Ks analysis suggests that all 13 PCGs have undergone a strong purifying selection. The phylogenetic analysis suggests the monophyly of Mordellidae, Meloidae, Oedemeridae, Pyrochroidae, Salpingidae, Scraptiidae, Lagriidae, and Tenebrionidae, and the Mordellidae is close to the Ripiphoridae. The "Tenebrionidae clade" and "Meloidae clade" are monophyletic, and both of them are sister groups. In the "Meloidae clade," Meloidae is close to Anthicidae. In the "Tenebrionidae clade," the family Lagriidae and Tenebrionidae are sister groups. The divergence time analysis suggests that Tenebrionoidea originated in the late Jurassic, Meloidae Mordellidae, Lagriidae, and Tenebrionidae in the Cretaceous, Oedemeridae in Paleogene. The work lays a base for the study of mtgenome, phylogenetics, and evolution of the superfamily Tenebrionoidea.

Post-mortem detection of hemoplasmas (hemotropic Mycoplasma spp.) in South American fur seal (Arctocephalus australis) sampled in Rio Grande do Sul State, southern Brazil

Hemotropic mycoplasmas are bacteria that attaches to erythrocytes surface, which some species presents zoonotic concerns. In the suborder Pinnipedia, genera Otaria and Arctocephalus are prominent in Brazil. This study investigated the occurrence of hemoplasmas in Arctocephalus sp. and Otaria flavescens found dead along the coast of a Southern Brazilian State. DNA from 135 spleen samples were extracted and subjected to conventional PCR protocols, targeting the 16 S rRNA and 23 S rRNA gene. Three (2.22 %) Arctocephalus australis were positive in the 16 S rRNA gene, and no samples amplified in the 23 S rRNA gene. Samples from this study clustered with Zalophus californianus and Arctocephalus tropicalis mycoplasmas on a Bayesian phylogenetic analysis. Genetic diversity analysis suggested distinct genotypes, indicating A. australis as a new host for hemoplasma, and also a potential putative novel hemoplasma genotype. These findings raises future awareness for pinnipeds conservation, and adds Mycoplasma spp. to be taken into consideration when clinically evaluating rescued animals.

Identification of fish specimens of the Tocantins River, Brazil, using DNA barcoding

The fish fauna of the Tocantins River possesses many endemic species; however, it is little studied in molecular terms and is quite threatened by the construction of several hydroelectric dams. Therefore, the objective of this study was to identify the ichthyofauna of the Tocantins River using DNA barcoding. For this, collections were carried out in five points of this river, which resulted in the capture of 725 individuals from which partial sequences of the cytochrome oxidase subunit I (COI) gene were obtained for genetic analysis. A total of 443 haplotypes were recovered with the mean intraspecific K2P genetic distance of 1.82%. Altogether, 138 species were identified based on morphological criteria, which was a quantity that was much lower than that indicated by the four molecular methods (assemble species by automatic partitioning [ASAP], barcode index number [BIN], generalized mixed Yule coalescent (GMYC), and Bayesian Poisson tree processes [bPTP]) through which 152-157 molecular entities were identified. In all, 41 unique BINs were obtained based on the data generated in the BOLDSystems platform. According to the result indicated by ASAP (species delimitation approach considered the most appropriate in the present study), there was an increase of 17 molecular entities (12.32%), when compared to the number of species identified through their morphological criteria, as it can show cryptic diversity, candidates for new species, and misidentifications. There were 21 incongruities indicated between the different identification approaches for species. Therefore, it is suggested that these taxonomic problems be cautiously evaluated by experts to solve such taxonomic issues.

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

PURPOSE

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

METHODS

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

RESULTS

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

CONCLUSION

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

Bacterial burden and molecular characterization of Coxiella burnetii in shedding pregnant and postpartum ewes from Saint Kitts

This study aimed to evaluate the bacterial burden and perform molecular characterization of Coxiella burnetii during shedding in pregnant (vaginal, mucus and feces) and postpartum (vaginal mucus, feces and milk) ewes from Saint Kitts. Positive IS1111 DNA (n=250) for C. burnetii samples from pregnant (n=87) and postpartum (n=74) Barbados Blackbelly ewes in a previous investigation were used for this study. Vaginal mucus (n=118), feces (n=100), and milk (n=32) positive IS1111 C. burnetii-DNA were analysed by real time qPCR (icd gene). For molecular characterization of C. burnetii, selected (n=10) IS1111 qPCR positive samples were sequenced for fragments of the IS1111 element and the 16 S rRNA gene. nBLAST, phylogenetic and haplotype analyses were performed. Vaginal mucus, feces and milk had estimated equal amounts of bacterial DNA (icd copies), and super spreaders were detected within the fecal samples. C. burnetii haplotypes had moderate to high diversity, were ubiquitous worldwide and similar to previously described in ruminants and ticks and humans.

Contrasting population genetic structure of three semi-terrestrial brachyuran crabs on the coast of the Japanese archipelago

Anthropogenic activities have reduced ecotones between the ocean and land, which is likely to threaten the population of brackish-water brachyuran crabs. To assess the current status of these crabs, we examine the population genetic structures of three semi-terrestrial brachyuran crabs widely distributed along the coast of the Japan and to clarify factors determining their genetic structures. We collected 184 Orisarma dehaani, 252 Chiromantes haematocheir, and 151 Helice tridens crabs from 36 localities of the Japanese archipelago. Genome-wide SNP data from these crabs were analyzed using MIG-seq. Bayesian clustering of STRUCTURE and DAPC analysis were used to identify genetically disturbed populations and to visualize genetic differentiation between local populations. Genetic population structure showed clear differentiation between populations on the Pacific coast of the Tohoku region and on other Japanese coasts in O. dehaani, but not in C. haematocheir or H. tridens. The inbreeding coefficient of O. dehaani was significantly higher on the Pacific coast of the Tohoku region compared to other Japanese coasts. C. haematocheir and H. tridens had homogeneous genetic structures along the Japanese coast, but showed genetic differentiation of a local population at their range limits. Thus, O. dehaani showed little gene flow and clear genetic differentiation between populations in the Tohoku Pacific region and those on other Japanese coasts due to ocean currents. Although such a regional differentiation was not found in C. haematocheir and H. tridens, one population of C. haematocheir was genetically isolated at the edge of its distribution range and likely vulnerable to environmental changes.