NOVOPlasty: de novo assembly of organelle genomes from whole genome data

The complete chloroplast genome of Taraxacum albidum (Asteraceae), a Japanese endemic dandelion

Taraxacum albidum, a perennial herb of the Asteraceae family, exhibits both tetraploid and pentaploid in Japan. This study sequenced and characterized the complete chloroplast genome of T. albidum, revealing a 151,451 bp sequence with a typical quadripartite structure, comprising one large single-copy (LSC) region of 84,052 bp, one small single-copy (SSC) region of 18,541 bp, and two inverted repeat (IR) regions, IRa and IRb, each 24,429 bp in length. The chloroplast genome, excluding duplicates, contained 113 unique genes, including 79 protein-coding genes, 30 transfer RNA genes, and four ribosomal RNA genes. The GC content of this genome was 37.7%. Phylogenetic analysis revealed that T. albidum is most closely related to T. mongolicum, with the chloroplast genome sequences being nearly identical, differing by only one nucleotide. These findings suggest that the maternal lineage of T. albidum likely originates from T. mongolicum or its closely related species.

The complete chloroplast genome of Rhododendronambiguum and comparative genomics of related species

Rhododendron Linnaeus, 1753, the largest genus of woody plants in the Northern Hemisphere, includes some of the most significant species in horticulture. Rhododendronambiguum Hemsl, 1911, a member of subsection Triflora Sleumer 1947, exemplifies typical alpine Rhododendron species. The analysis of the complete chloroplast genome of R.ambiguum offers new insights into the evolution of Rhododendron species and enhances the resolution of phylogenetic relationships. This genome is composed of 207,478 base pairs, including a pair of inverted repeats (IRs) of 47,249 bp each, separated by a large single-copy (LSC) region of 110,367 bp and a small single-copy (SSC) region of 2,613 bp. It contains 110 genes: 77 protein-coding genes, 29 tRNAs, four unique rRNAs (4.5S, 5S, 16S, and 23S), with 16 genes duplicated in the IRs. Comparative analyses reveal substantial diversity in the Rhododendron chloroplast genome structures, identifying a fourth variant pattern. Specifically, four highly divergent regions (trnI-rpoB, ndhE-psaC, rpl32-ndhF, rrn16S-trnI) were noted in the intergenic spacers. Additionally, 76 simple sequence repeats were identified. Positive selection signals were detected in four genes (cemA, rps4, rpl16, and rpl14), evidenced by high Ka/Ks ratios. Phylogenetic reconstruction based on two datasets (shared protein-coding genes and complete chloroplast genomes) suggests that R.ambiguum is closely related to R.concinnum Hemsley, 1889. However, the phylogenetic positions of subsection Triflora Pojarkova, 1952 species remain unresolved, indicating that the use of complete chloroplast genomes for phylogenetic research in Rhododendron requires careful consideration. Overall, our findings provide valuable genetic information that will enhance understanding of the evolution, molecular biology, and genetic improvement of Rhododendron spieces.

Chromosome-level genome assembly of marmalade hoverfly Episyrphus balteatus (Diptera: Syrphidae)

Episyrphus balteatus can provide dual ecosystem services including pest control and pollination, which the larvae are excellent predators of aphid pest whereas adults are efficient pollinator. In this study, we assembled a high-quality genome of E. balteatus from northern China geographical population at the chromosome level by using Illumina, PacBio long reads, and Hi-C technologies. The 467.42 Mb genome was obtained from 723 contigs, with a contig N50 of 9.16 Mb and Scaffold N50 of 118.85 Mb, and 90.25% (431.75 Mb) of the assembly was anchored to 4 pseudo-autosomes and one pseudo-heterosome. In total, 14,848 protein-coding genes were annotated, and 95.14% of genes were fully represented in NR, GO, KEGG databases. Besides, we also obtained the mitochondrial genome of E. balteatus of 16, 837 bp in length with 37 typical mitochondrial genes. Overall, this high-quality genome is valuable for evolutionary and genetic studies of E. balteatus and other Syrphidae hoverfly species.

The plastome and phylogenetic status of Cotoneaster rosiflorus (Rosaceae)

Endemic to Taiwan Province, China, Cotoneaster rosiflorus Kun-Cheng Chang & Fu-Yuan Lu 2011 (Rosaceae) holds significant ecological and ornamental importance. Despite its value, research on its molecular data and phylogenetic position has remained limited. In this study, we addressed this gap by sequencing the genome-skimming data, assembling its plastome, and investigating its phylogenetic position. The plastome, spanning 159,449 bp in length, consisted of a large single-copy (87,433 bp), a small single-copy (19,262 bp), and two inverted repeat regions (26,377 bp). We annotated a total of 128 functional genes, including 84 protein-coding genes, 36 transfer genes, and eight ribosomal RNA genes. The phylogenetic results indicated that C. rosiflorus is closely related to C. dammerii, suggesting that C. rosiflorus might have captured its chloroplast from C. dammerii through hybridization and introgression events. This study offered valuable insights for forthcoming phylogenetic and population genetic investigations of Cotoneaster.

Genome sequencing of Coryphaenoides yaquinae reveals convergent and lineage-specific molecular evolution in deep-sea adaptation

Abyssal (3501-6500 m) and hadal (>6500 m) fauna evolve under harsh abiotic stresses, characterized by high hydrostatic pressure, darkness and food shortage, providing unique opportunities to investigate mechanisms underlying environmental adaptation. Genomes of several hadal species have recently been reported. However, the genetic adaptation of deep sea species across a broad spectrum of ocean depths has yet to be thoroughly investigated, due to the challenges imposed by collecting the deep sea species. To elucidate the correlation between genetic innovation and vertical distribution, we generated a chromosome-level genome assembly of the macrourids Coryphaenoides yaquinae, which is widely distributed in the abyssal/hadal zone ranging from 3655 to 7259 m in depth. Genomic comparisons among shallow, abyssal and hadal-living species identified idiosyncratic and convergent genetic alterations underlying the extraordinary adaptations of deep-sea species including light perception, circadian regulation, hydrostatic pressure and hunger tolerance. The deep-sea fishes (Coryphaenoides Sp. and Pseudoliparis swirei) venturing into various ocean depths independently have undergone convergent amino acid substitutions in multiple proteins such as rhodopsin 1, pancreatic and duodenal homeobox 1 and melanocortin 4 receptor which are known or verified in zebrafish to be related with vision adaptation and energy expenditure. Convergent evolution events were also identified in heat shock protein 90 beta family member 1 and valosin-containing protein genes known to be related to hydrostatic pressure adaptation specifically in fishes found around the hadal range. The uncovering of the molecular convergence among the deep-sea species shed new light on the common genetic innovations required for deep-sea adaptation by the fishes.

Daphnia sp. (Branchiopoda: Cladocera) Mitochondrial Genome Gene Rearrangement and Phylogenetic Position Within Branchiopoda

In high-altitude (4500 m) freshwater lakes, Daphnia is the apex species and the dominant zooplankton. It frequently dwells in the same lake as the Gammarid. Branchiopoda, a class of Arthropoda, Crustacea, is a relatively primitive group in the subphylum Crustacea, which originated in the Cambrian period of the Paleozoic. The complete mitogenome sequence of Daphnia sp. (Branchiopoda: Cladocera) was sequenced and annotated in this study and deposited in GenBank. The sequence structure of this species was studied by comparing the original sequences with BLAST. In addition, we have also researched the mechanisms of their mitochondrial gene rearrangement by establishing a model. We have used the Bayesian inference [BI] and maximum likelihood [ML] methods to proceed with phylogenetic analysis inference, which generates identical phylogenetic topology that reveals the phylogenetic state of Daphnia. The complete mitogenome of Daphnia sp. shows that it was 15,254 bp in length and included two control regions (CRs) and 37 genes (13 protein-coding genes, 22 tRNAs and two ribosomal RNAs [16S and 12S]). In addition to tRNA-Ser (GCT), other tRNAs have a typical cloverleaf secondary structure. Meanwhile, the mitogenome of Daphnia sp. was clearly rearranged when compared to the mitogenome of typical Daphnia. In a word, we report a newly sequenced mitogenome of Daphnia sp. with a unique rearrangement phenomenon. These results will be helpful for further phylogenetic research and provide a foundation for future studies on the characteristics of the mitochondrial gene arrangement process in Daphnia.

Complete Mitochondrial Genome of Clupeonella abrau (Maliatsky, 1930) (Clupeiformes), an Endemic Freshwater Fish from the Lake Abrau (Russia)

The Abrau sprat (tyulka or sardelka) Clupeonella abrau (Maliatsky, 1930) is an endemic fish of the Lake Abrau (Krasnodar Krai, Russia). The full C. abrau mitochondrial genome (16 650 bp) showed a gene arrangement conserved in Clupeidae and 98.8% similarity with the mitochondrial genome of the related species Black and Caspian Sea sprat C. cultriventris from the Black Sea. The COX1 gene sequence was additionally studied in a museum specimen collected in the Lake Abrau in 1938. Variability in modern Abrau sprat COX1 gene locus was estimated at approximately 0.15%, the difference between C. abrau and C. cultriventris was 1.2%, and the difference between the museum and modern C. abrau specimens from the Lake Abrau was 0.92%. The study confirmed that the Abrau sprat is present in the fish community and is capable of reproducing in the lake. Various scenarios were proposed to explain colonization of the Lake Abrau by C. abrau.

Characterizing the complete mitogenome of Odontothrips phaseoli (Thysanoptera: Thripidae) and its mitochondrial phylogeny

Described originally from Heilongjiang, China, Odontothrips phaseoli is a potential pest of threatening bean plant in northern China. The complete mitochondrial genome of O. phaseoli was sequenced and assembled, with a total length of 15,540 bp. Within this genome, 37 genes have been identified: 13 PCGs, 22 tRNAs, two rRNAs, and two putative control regions. Most PCGs terminate with TAA, while four genes (atp8, nad1, nad2 and nad4) use an incomplete 'T' and nad6 employs TAG as the stop codon. Compared to the mitogenome of the ancestral insect, O. phaseoli displays significant gene rearrangement. However, it retains three conserved gene blocks in common with its related species, Megalurothrips usitatus, both of which belong to the Megalurothrips genus-group. The phylogenetic tree, constructed based on the entire mitogenome dataset of all thrips species available in NCBI, shows that the two species cluster closely together. This alignment might underscore the close link between gene arrangements and the phylogeny relationships.

Complete mitochondrial genome of the edible mushroom Singerocybe alboinfundibuliformis (Clitocybaceae, Agaricales)

Singerocybe alboinfundibuliformis (Seok et al.) Yang, Qin & Takah 2014 is an edible mushroom distributed in several East or Southeast Asian countries. Herein, we report the mitochondrial genome of S. alboinfundibuliformis based on Illumina sequencing data. The overall length of the mitochondrial genome is 64,279 bp, with a GC content of 29.0%. It contains 14 typical protein-coding genes, 27 tRNA genes, two rRNA genes, and 13 intergenic ORFs. Most of these genes (39 out of 56) are transcribed at the forward strand, and few (17 out of 56) are transcribed at the reverse strand. Among these genes, only the rnl gene is invaded by an intron, and all other genes are intron-free. Phylogenetic analysis based on mitochondrial amino acid sequences supports the phylogenetic position of S. alboinfundibuliformis in Clitocybaceae, being close to Lepista sordida (Schumach.) Singer 1951. This study serves as a springboard for future investigation on fungal evolution in Clitocybaceae.

Complete mitochondrial genome of the bird's nest fungus Nidula shingbaensis (Nidulariaceae, Agaricales)

Bird's nest fungi involve six different genera, but only one of these genera (i.e. Cyathus) have available mitochondrial genomes (mitogenomes) to date. In this study, we report the first mitogenome in the genus Nidula with Nidula shingbaensis K. Das & R.L. Zhao 2013 as a representative. The mitogenome is a circular molecule of 65,793 bp with a GC content of 26.2%. There are a total of 43 genes, including 14 typical protein-coding genes, 26 tRNA genes, two rRNA genes, and one free-standing intergenic open reading frame (ORF). Three introns (two in cox1 and one in cob) are present in the mitogenome, with each containing an ORF encoding for a LAGLIDADG endonuclease. Phylogenetic analysis based on mitochondrial amino acid sequences confirms the phylogenetic placement of N. shingbaensis in Nidulariaceae in Agaricales. This study serves as a springboard for future investigation on fungal evolution in Nidulariaceae.

Decoding the Chloroplast Genome of Tetrastigma (Vitaceae): Variations and Phylogenetic Selection Insights

Tetrastigma (Vitaceae) is known for its ornamental, medicinal, and ecological significance. However, the structural and variational characteristics of the Tetrastigma chloroplast genome and their impact on phylogenetic relationships remain underexplored. This study utilized bioinformatics methods to assemble and annotate the chloroplast genomes of 10 Tetrastigma species and compare them with five previously sequenced species. This study analyzed gene composition, simple sequence repeats, and codon usage patterns, revealing a high A/T content, uniquely identified pentanucleotide repeats in five species and several preferred codons. In addition, comparative analyses were conducted of the chloroplast genomes of 15 Tetrastigma species, examining their structural differences and identifying polymorphic hotspots (rps16, rps16-trnQ, trnS, trnD, psbC-trnS-psbZ, accD-psaI, psbE-petL-petG, etc.) suitable for DNA marker development. Furthermore, phylogenetic and selective pressure analyses were performed based on the chloroplast genomes of these 15 Tetrastigma species, validating and elucidating intra-genus relationships within Tetrastigma. Futhermore, several genes under positive selection, such as atpF and accD, were identified, shedding light on the adaptive evolution of Tetrastigma. Utilizing 40 Vitaceae species, the divergence time of Tetrastigma was estimated, clarifying the evolutionary relationships within Tetrastigma relative to other genera. The analysis revealed diverse divergences of Tetrastigma in the Miocene and Pliocene, with possible ancient divergence events before the Eocene. Furthermore, family-level selective pressure analysis identified key features distinguishing Tetrastigma from other genera, showing a higher degree of purifying selection. This research enriches the chloroplast genome data for Tetrastigma and offers new insights into species identification, phylogenetic analysis, and adaptive evolution, enhancing our understanding of the genetic diversity and evolutionary history of these species.

Camellia neriifolia and Camellia ilicifolia (Theaceae) as separate species: evidence from morphology, anatomy, palynology, molecular systematics

BACKGROUND

The systematic status of sect. Tuberculata and its taxonomy have recently attracted considerable attention. However, the different bases for defining the characteristics of sect. Tuberculata has led to many disagreements among the plants in this group. Camellia neriifolia and Camellia ilicifolia have been the subject of taxonomic controversy and have been treated as different species or varieties of the same species. Therefore, it is important to use multiple methods, i.e., integrative taxonomy, to determine the taxonomic status of C. neriifolia and C. ilicifolia. This is the first study to systematically explore the taxonomic position of these two plants on the basis of Morphology, Anatomy, Palynology and Molecular Systematics.

RESULTS

Extensive specimen reviews and field surveys showed that many differences exist in C. neriifolia and C. ilicifolia, such as the number of trunk (heavily debarked vs. slightly peeling), leaf type (smooth thin leathery, shiny vs. smooth leathery, obscure or slightly shiny), leaf margin (entire vs. serrate), flower type (subsessile vs. sessile), number of styles (3-4 vs. 3), and sepal (ovate vs. round). Moreover, C. neriifolia has a more distinctive faint yellow flower color, and trunk molting was more severe in C. neriifolia than that in C. ilicifolia. In addition, micromorphological analysis of the leaf epidermis showed that the two species differed in the anticlinal wall, stomatal apparatus, and stomatal cluster, and pollen morphology analyses based on pollen size, germination furrow, and polar and equatorial axes showed that they are both distinct from each other. The results of the phylogenetic tree constructed based on the whole chloroplast genome, protein-coding genes, and ITS2 showed that both C. ilicifolia and C. neriifolia were clustered in different branches and gained high support.

CONCLUSIONS

The results combine morphology, anatomy, palynology, and molecular systematics to treat both C. neriifolia and C. ilicifolia as separate species in the sect. Tuberculata, and the species names continue to be used as they were previously. In conclusion, clarifying the taxonomic status of C. neriifolia and C. ilicifolia deepens our understanding of the systematic classification of sect. Tuberculata.

Coptisaustrogaoligongensis (Ranunculaceae), a new species from West Yunnan, China

Based on morphological and plastid data, we have described and confirmed that Coptisaustrogaoligongensis distributed in Tongbiguan Provincial Nature Reserve, Yingjiang County, Yunnan Province, is a new species of Coptis. It is distinctly different from C.teetasubsp.teeta and C.teetasubsp.lohitensis with differences mainly reflected in the following features: former leaf segment lobes contiguous to each other, and lateral segments equal to central one; plants without developed stolons; inflorescences with only 1-3 flowers; petals have short claws. Phylogenetic analysis indicated that C.austrogaoligongensis is a sister to C.teetasubsp.teeta and C.teetasubsp.lohitensis.

First complete mitochondrial genome of Pselliophora (Diptera, Tipulidae): genome description and phylogenetic implications

Pselliophora is widely distributed in Eurasia and China. To explore the characteristics of the mitogenome of Pselliophora and reveal phylogenetic relationships, the mitogenome of Pselliophora bifascipennis Brunetti, 1911 was sequenced and annotated. This is the first complete mitochondrial genome in this genus. Its mitogenome is 15821 bp in length, containing 13 protein-coding genes, 22 tRNA and 2 rRNA genes. Nucleotide compositions of its whole mitogenome are 39.09% for A, 38.49% for T, 13.42% for C, and 9.01% for G. Consistent with previous observations of Tipulidae species, the mitogenome of Pselliophora bifascipennis is highly conserved in gene size, organization and codon usage, and secondary structures of tRNAs. Most tRNAs have the typical clover-leaf structure. The control region is 1006 bp long with an A + T content of 92.7%. Phylogenetic tree analysis using the sequences of the mitochondrial genomes of Pselliophora bifascipennis and other Tipulidae species showed that Pselliophora bifascipennis is closely related to Tanyptera hebeiensis. These two species are grouped on the same branch, which is in accordance with the traditional morphological classification. The results of this study lay a foundation for screening molecular markers of mitochondrion for molecular identification and genetic structure research in Tipulidae species.

Complete genome sequence of strain Lactiplantibacillus plantarum beLP1 a potential probiotic strain with antimicrobial properties isolated from kimchi

The complete genome of the potential probiotic Lactiplantibacillus plantarum strain beLP1, isolated from kimchi in South Korea, was sequenced using Illumina and PacBio technologies. The genome comprises one circular chromosome and one plasmid without antimicrobial resistance genes.

The complete female mitogenome of Potomida semirugata (Lamarck, 1819)

Freshwater mussels perform important ecological functions in ecosystems, such as water filtration and energy cycling. Unlike marine bivalves, freshwater mussels have unique characteristics including internal fertilization and parental care. Some freshwater mussels are facing a high risk of extinction due to several factors such as climate change and habitat loss. Potomida semirugata (Lamarck, 1819) is one of the freshwater mussel species with a high risk of extinction and listed as Endangered in the Red List of the International Union for Conservation of Nature. Here, we present the first F-type mitogenome sequence of P. semirugata. The genome was sequenced on an Illumina high-throughput platform from a P. semirugata specimen collected from the Tersakan River (Turkey). The 16,093 bp mitochondrial genome sequence contains 13 protein-coding genes, 22 transfer RNAs, and two ribosomal RNAs. Phylogenetic analysis placed P. semirugata in the Lamprotulini clade with Potomida littoralis, as expected. Potomida semirugata is a poorly studied species and the genomic resource provided here will contribute to a better understanding of its biological characterization.

Chloroplast genome of plantago major, a medicinal plant in China

Plantago major (Plantaginaceae) is a medicinal plant in Chinese folk culture, known for its famous medicinal components such as plantagomain. In this study, we conducted genome sequencing of P. major using Illumina sequencing technology. The assembled complete chloroplast genome had a length of 165,044 bp, comprising a large single-copy regions (82,963 bp), a small single-copy regions (4,633 bp), and a pair of inverted repeat regions (38,724 bp). A total of 140 genes were detected, including 94 CDS, 38 tRNA, and 8 rRNA. Phylogenetic analysis revealed a close genetic relationship between P. major and P. rigida. These findings provide valuable data for a comprehensive understanding of the biological characteristics of P. major.

Characterization, comparison, and phylogenetic analyses of chloroplast genomes of Euphorbia species

The genus Euphorbia (Euphorbiaceae) has near-cosmopolitan distribution and serves as a significant resource for both ornamental and medicinal purposes. Despite its economic importance, Euphorbia's taxonomy has long been challenged by the intricate nature of morphological traits exhibiting high levels of convergence. While molecular markers are essential for phylogenetic studies, their availability for Euphorbia has been limited. To address this gap, we conducted comparative analyses focusing on the chloroplast (CP) genomes of nine Euphorbia species, incorporating three newly sequenced and annotated accessions. In addition, phylogenetic informativeness and nucleotide diversity were computed to identify candidate markers for phylogenetic analyses among closely related taxa in the genus. Our investigation revealed relatively conserved sizes and structures of CP genomes across the studied species, with notable interspecific variations observed primarily in non-coding regions and IR/SC borders. By leveraging phylogenetic informativeness and nucleotide diversity, we identified rpoB gene as the optimal candidate for species delimitation and shallow-level phylogenetic inference within the genus. Through this comprehensive analysis of CP genomes across multiple taxa, our study sheds light on the evolutionary dynamics and taxonomic intricacies of Euphorbia, offering valuable insights into its CP genome evolution and taxonomy.

Comparative Analysis of Complete Mitochondrial Genomes of Five Chromodorididae Species (Nudibranchia:Doridina)

Mitochondrial genome is an important molecular marker for exploring the phylogenetic relationships of species and revealing molecular evolution. In the present study, 5 mitogenomes of Chromodorididae (Chromodoris lochi, Chromodoris colemani, Chromodoris elisabethina, Chromodoris annae and Hypselodoris whitei) were systemically investigated. The lengths of the mitogenomes sequences were 14248 bp, 14257 bp, 14252 bp, 14254 bp and 14856 bp, respectively. Most protein-coding genes (PCGs) were initiated with the common ATG codon and terminated with the TAA and TAG. We calculated Ka/Ks values for all 13 PCGs of Chromodorididae species, all ratios were less than 1, indicating selection by purification. Phylogenetic relationships were constructed by Bayesian inference (BI) and maximum likelihood (ML) methods based on all complete genomes of 50 species, primarily from the family Chromodorididae (Doridina) and 2 outgroups. This phylogenetic tree provided further additional references for the classification of the suborder Doridina. Gene rearrangement suggested a more conserved pattern of gene sequences in the superfamily Chromodoridoidea. These results and newly sequenced will contribute to a better understanding of Chromodorididae and provide reference for further phylogenetic studies.

Experimental evaluation of genetic variability based on DNA metabarcoding from the aquatic environment: Insights from the Leray COI fragment

Intraspecific genetic variation is important for the assessment of organisms' resistance to changing environments and anthropogenic pressures. Aquatic DNA metabarcoding provides a non-invasive method in biodiversity research, including investigations at the within-species level. Through the analysis of eDNA samples collected from the Peter the Great Gulf of the Japan Sea, in this study, we aimed to evaluate the identification of Amplicon Sequence Variants (ASVs) in marine eDNA among abundant species of the Zostera sp. community: Hexagrammos octogrammus, Pholidapus dybowskii (Teleostei: Perciformes), and Pandalus latirostris (Arthropoda: Decapoda). These species were collected from two distant locations to produce mock communities and gather aquatic eDNA both on the community and individual level. Our approach highlights the efficacy of eDNA metabarcoding in capturing haplotypic diversity and the potential for this methodology to track genetic diversity accurately, contributing to conservation efforts and ecosystem management. Additionally, our results elucidate the impact of nuclear mitochondrial DNA segments (NUMTs) on the reliability of metabarcoding data, indicating the necessity for cautious interpretation of such data in ecological studies. Moreover, we analyzed 83 publicly available COI sequence datasets from common groups of multicellular organisms (Mollusca, Echinodermata, Crustacea, Polychaeta, and Actinopterygii). The results reflect the decrease in population diversity that arises from using the metabarcode compared to the COI barcode.

Genomic insights into hybrid zone formation: The role of climate, landscape, and demography in the emergence of a novel hybrid lineage

Population demographic changes, alongside landscape, geographic and climate heterogeneity, can influence the timing, stability and extent of introgression where species hybridise. Thus, quantifying interactions across diverged lineages, and the relative contributions of interspecific genetic exchange and selection to divergence at the genome-wide level is needed to better understand the drivers of hybrid zone formation and maintenance. We used seven latitudinally arrayed transects to quantify the contributions of climate, geography and landscape features to broad patterns of genetic structure across the hybrid zone of Populus trichocarpa and P. balsamifera and evaluated the demographic context of hybridisation over time. We found genetic structure differed among the seven transects. While ancestry was structured by climate, landscape features influenced gene flow dynamics. Demographic models indicated a secondary contact event may have influenced contemporary hybrid zone formation with the origin of a putative hybrid lineage that inhabits regions with higher aridity than either of the ancestral groups. Phylogenetic relationships based on chloroplast genomes support the origin of this hybrid lineage inferred from demographic models based on the nuclear data. Our results point towards the importance of climate and landscape patterns in structuring the contact zones between P. trichocarpa and P. balsamifera and emphasise the value whole genome sequencing can have to advancing our understanding of how neutral processes influence divergence across space and time.

Characterization of the complete mitochondrial genome of Mucor indicus Lendn. 1930 (Mucorales: Mucoraceae), isolated from the wine fermentation system

Mucor indicus Lendn. 1930 has been widely used in food fermentation; however, its mitochondrial genome characteristics are not well understood. In this study, the complete mitochondrial genome of M. indicus was obtained, which was 61,400 bp in length with a GC content of 33%. The M. indicus mitochondrial genome was found to contain 14 core protein-coding genes, four free-standing ORFs, 18 intronic ORFs, 26 tRNAs, and two rRNA genes. Phylogenetic trees were generated for 25 early-differentiated fungi using the Bayesian inference (BI) method, which demonstrated that M. indicus is closely related to Mucor piriformis. This study provides useful information for the classification and evolution of Mucor species or other early-differentiated fungi.

The absence of canonical respiratory complex I subunits in male-type mitogenomes of three Donax species

Bivalves are an extraordinary class of animals in which species with a doubly uniparental inheritance (DUI) of mitochondrial DNA have been described. DUI is characterized as a mitochondrial homoplasmy of females and heteroplasmy of male individuals where F-type mitogenomes are passed to the progeny with mother egg cells and divergent M-type mitogenomes are inherited with fathers sperm cells. However, in most cases only male individuals retain divergent mitogenome inherited with spermatozoa. Additionally, in many of bivalves, unique mitochondrial features, like additional genes, gene duplication, gene extensions, mitochondrial introns, and recombination, were observed. In this study, we sequenced and assembled male-type mitogenomes of three Donax species. Comparative analysis of mitochondrial sequences revealed a lack of all seven NADH dehydrogenase subunits as well as the presence of three long additional open reading frames lacking identifiable homology to any of the existing genes.

The complete mitochondrial genome of a Biwa goby, Gymnogobius isaza (Tanaka, 1916)

We determined the complete mitochondrial DNA sequence of a Biwa goby, Gymnogobius isaza (Tanaka, 1916) using next-generation sequencing methods. The composition of its mitogenome is the same as that observed in most other vertebrates, comprising of 13 protein-coding genes, 22 tRNA genes, two rRNA genes, and two control regions. Our molecular phylogenetic analysis confirmed the close phylogenetic relationship between G. isaza and G. petschiliensis. This mitogenome information will be useful for distribution surveys using environmental DNA and the development of conservation strategies for this species.

Twenty-nine newly sequenced genomes and a comprehensive genome dataset for the insect endosymbiont Buchnera

Most phloem-feeding insects face nutritional deficiency and rely on their intracellular symbionts to provide nutrients, and most of endosymbiont genomes have undergone reduction. However, the study of genome reduction processes of endosymbionts has been constrained by the limited availability of genome data from different insect lineages. The obligate relationship between aphids and Buchnera aphidicola (hereafter Buchnera) makes them a classic model for studying insect-endosymbiont interaction. Here, we report 29 newly sequenced Buchnera genomes from 11 aphid subfamilies, and a comprehensive dataset based on 90 Buchnera genomes from 14 aphid subfamilies. The dataset shows a significant genomic difference of Buchnera among different aphid lineages. The dataset exhibits a more balanced distribution of Buchnera (from 14 aphid subfamilies) genome sizes, ranging from 400 kb to 600 kb, which can illustrate the genome reduction process of Buchnera. The new genome data provide valuable insights into the microevolutionary processes leading to genomic reduction of insect endosymbionts.

Characterization of the complete chloroplast genome of the rare medicinal plant: Mandragora caulescens (Solanaceae)

In this study, we assembled high-quality chloroplast genomes of Mandragora caulescens through a reference-guided approach using high-throughput Illumina sequencing reads. The resulting chloroplast genome assembly displayed a typical quadripartite structural organization, comprising a large single-copy (LSC) region of 85,233 bp, two inverted repeat (IR) regions of 25,685 bp each, and a small single-copy (SSC) region of 18,207 bp. The chloroplast genome harbored 141 complete genes, and its overall GC content was 38.0%. In maximum-likelihood (ML) and Bayesian inference (BI) trees, the 19 Solanaceae species formed a monophyletic group, dividing into two main clades. M. caulescens and Nicandra physalodes formed a monophyletic group, suggesting a close relationship between the two species. The M. caulescens cp genome presented in this study lays a good foundation for further genetic and genomic studies of the Solanaceae.

The complete mitochondrial genome of Membranipora villosa Hincks, 1880 (Bryozoa: Gymnolaemata: Cheilostomatida): phylogenetic relationship of two kelp-encrusting bryozoans within the suborder Membraniporina

The two commonest kelp-encrusting bryozoans, Membranipora villosa and M. membranacea, are difficult to distinguish morphologically. Molecular studies of M. villosa should thus be helpful for the identification of both species because the mitogenome of M. membranacea was already sequenced. The complete mitogenome of M. villosa collected from Sinjido was determined in this study through Illumina NovaSeq sequencing. Maximum-likelihood (ML) analysis was based on concatenated 13 protein-coding genes dataset from nine bryozoan species. The mitogenome length was 15,407 bp, and its gene arrangement was similar to those of the mitogenome of other membraniporids, having 13 PCGs, two ribosomal RNAs, and 22 tRNAs. It had an overall A + T content of 63.7% (29.7% A, 16.7% C, 19.6% G, and 34.0% T). M. villosa and M. membranacea showed sequence differences of 20% for the total length of mitogenome and 16.1.% for 13 PCGs. Molecular data definitely consider them to be separate species. Phylogenetic analyses based on the amino acids of 13 PCGs indicated that M. villosa has the closest relationship with another kelp-encrusting bryozoan, M. membranacea of membraniporids. The phylogenetic position of genera and families within the suborder Membraniporina coincides with the Bayesian phylogenetic analysis of the mixed concatenated alignment consisting of three partitions.

Complete mitogenome and intra-family comparative mitogenomics showed distinct position of Pama Croaker Otolithoides pama

The Pama Croaker, Otolithoides pama, is an economically important fish species in Bangladesh. Intra-family similarities in morphology and typical barcode sequences of cox1 create ambiguities in its identification. Therefore, morphology and the complete mitochondrial genome of O. pama, and comparative mitogenomics within the family Sciaenidae have been studied. Extracted genomic DNA was subjected to Illumina-based short read sequencing for De-Novo mitogenome assembly. The complete mitogenome of O. pama (Accession: OQ784575.1) was 16,513 bp, with strong AC biasness and strand asymmetry. Relative synonymous codon usage (RSCU) among 13 protein-coding genes (PCGs) of O. pama was also analyzed. The studied mitogenomes including O. pama exhibited consistent sizes and gene orders, except for the genus Johnius which possessed notably longer mitogenomes with unique gene rearrangements. Different genetic distance metrics across 30 species of Sciaenidae family demonstrated 12S rRNA and the control region (CR) as the most conserved and variable regions, respectively, while most of the PCGs undergone a purifying selection. Different phylogenetic trees were congruent with one another, where O. pama was distinctly placed. This study would contribute to distinguishing closely related fish species of Sciaenidae family and can be instrumental in conserving the genetic diversity of O. pama.

Characterization and phylogenetic analysis of the complete mitochondrial genome of Saccharomycopsis fibuligera (lindner) Klocker 1907 (saccharomycetales: saccharomycopsidaceae)

Saccharomycopsis fibuligera (Lindner) Klocker 1907 is frequently employed in the fermentation of metabolites such as citric acid, ethanol, mannitol, and pyruvate. Its heat tolerance and alcohol-producing capabilities during fermentation make it a desirable option for bread and wine production. To date, the mitochondrial genome of S. fibuligera has not been sequenced. In the present study, we obtained the full mitochondrial genome of S. fibuligera, which is 57,302 bp long and has a GC content of 24.40%. This genome contained 14 core protein-coding genes, 3 independent ORFs, 21 intronic ORFs, 25 tRNAs, and 2 rRNA genes. By utilizing the Bayesian inference phylogenetic method, we constructed phylogenetic trees for 24 Saccharomycotina fungi, which indicated that S. fibuligera is closely related to S. capsularis.

The complete chloroplast genome of Dicliptera tinctoria (Nees) Kostel. and comparative analysis of chloroplast genomes in Acanthaceae

Dicliptera tinctoria is a member of Acanthaceae, which has a wide distribution and contains potentially medicinal species, and exhibited pharmaceutical potentials. This study sequenced and characterized the complete chloroplast genome of Dicliptera tinctoria. The newly sequenced cpDNA of D. tinctoria was 150,733 bp in length and had a typical quadripartite structure consisting of a large single copy (LSC, 82,895 bp), a small single copy (SSC, 17,249 bp), and two inverted repeat (IRs, 25,295 bp each) regions. This genome also contained 80 protein-coding genes, 30 transfer RNAs, and four ribosomal RNAs, which is identical to other chloroplast genomes in Acanthaceae family. Nucleotides diversity analysis among chloroplast genomes of Acanthaceae species revealed eight hypervariable regions, including trnK_UUU-matK, trnC_GCA-petN, accD, rps12-clpP, rps3-rps19, ycf1-ndhF, ccsA-ndhD, and ycf1. Phylogenetic analysis revealed the paraphyly of Dicliptera species and monophyly in four Acanthaceae subfamilies. These results provide an overview of genomic variations in Acanthaceae chloroplast genome, which is helpful for further genomic studies.

The complete mitochondrial genome of Zaomma eriococci (hymenoptera: encyrtidae)

The complete mitochondrial genome of the Zaomma eriococci (Ferrière, 1955) (Hymenoptera: Encyrtidae) was obtained through next-generation sequencing, making the first reported complete mitochondrial genome of the genus Zaomma. The mitochondrial genome is 15,648 bp in length and includes 37 classical eukaryotic mitochondrial genes along with an A + T rich region. All 13 protein-coding genes (PCGs) initiate with typical ATN codons. Of these, 10 PCG genes terminate with TAA, while three terminate with TAG. Additionally, there are 22 tRNA genes, ranging in size from 62 to 70 bp. The maximum likelihood phylogenetic tree was constructed based on 13 PCGs, indicates that Z. eriococci is closely related to Tassonia gloriae. This mitochondrial genome will serve as a valuable molecular resource for species identification, genetic analysis, and comparative genomic studies of Z. eriococci, contributing to the growing collection of mitochondrial genomes within the family Encyrtidae.

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.

Phylogenetic analysis and divergence time estimation of Lycium species in China based on the chloroplast genomes

BACKGROUND

Lycium is an economically and ecologically important genus of shrubs, consisting of approximately 70 species distributed worldwide, 15 of which are located in China. Despite the economic and ecological importance of Lycium, its phylogeny, interspecific relationships, and evolutionary history remain relatively unknown. In this study, we constructed a phylogeny and estimated divergence time based on the chloroplast genomes (CPGs) of 15 species, including subspecies, of the genus Lycium from China.

RESULTS

We sequenced and annotated 15 CPGs in this study. Comparative analysis of these genomes from these Lycium species revealed a typical quadripartite structure, with a total sequence length ranging from 154,890 to 155,677 base pairs (bp). The CPGs was highly conserved and moderately differentiated. Through annotation, we identified a total of 128-132 genes. Analysis of the boundaries of inverted repeat (IR) regions showed consistent positioning: the junctions of the IRb/LSC region were located in rps19 in all Lycium species, IRb/SSC between the ycf1 and ndhF genes, and SSC/IRa within the ycf1 gene. Sequence variation in the SSC region exceeded that in the IR region. We did not detect major expansions or contractions in the IR region or rearrangements or insertions in the CPGs of the 15 Lycium species. Comparative analyses revealed five hotspot regions in the CPG: trnR(UCU), atpF-atpH, ycf3-trnS(GGA), trnS(GGA), and trnL-UAG, which could potentially serve as molecular markers. In addition, phylogenetic tree construction based on the CPG indicated that the 15 Lycium species formed a monophyletic group and were divided into two typical subbranches and three minor branches. Molecular dating suggested that Lycium diverged from its sister genus approximately 17.7 million years ago (Mya) and species diversification within the Lycium species of China primarily occurred during the recent Pliocene epoch.

CONCLUSION

The divergence time estimation presented in this study will facilitate future research on Lycium, aid in species differentiation, and facilitate diverse investigations into this economically and ecologically important genus.

Impatiensyingjingensis (Balsaminaceae), a new species from Sichuan, China

This study describes Impatiensyingjingensis X.Q. Song, B.N. Song & Biao Yang, sp. nov., a new species collected from the Yingjing area of the Giant Panda National Park. This new species is distributed at an altitude of 1400-2100 m, with a plant height of 30-130 cm. The flowers are purple-red or light purple red, with 3-9 flowers on each inflorescence and the dorsal auricle of the lateral united petals is thread-like and about 2 cm long, differing significantly from other species of Impatiens. Furthermore, molecular data, as well as micro-morphological evidence under SEM (of pollens), also support the establishment of the new species.

Evolutionary and phylogenetic insights from the mitochondrial genomic analysis of Diceraeus melacanthus and D. furcatus (Hemiptera: Pentatomidae)

The mitochondrial genomes of D. melacanthus and D. furcatus were sequenced and used to investigate the phylogenetic relationships with 54 species of Pentatomidae. Their mitogenomes are 17,197 and 15,444 bp-long, respectively, including 13 protein-coding genes (PCGs), 2 ribosomal RNA genes, and 22/21 transfer RNA genes, with conserved gene arrangement. Leu, Lys, and Ser were the most common amino acids in their PCGs. PCGs evolutionary analysis indicated their mitogenomes are under purifying selection, and the most conserved genes are from the cytochrome complex, reinforcing their suitability as markers for molecular taxonomy. We identified 490 mtSSRs in 56 Pentatomidae species, with large variation and a positive correlation between mtSSR number and genome size. Three mtSSRs were identified in each Diceraeus species. Only the mtSSR in the nad6 (D. melacanthus) and nad4 (D. furcatus) appear to have application as molecular markers for species characterization. Phylogenetic analysis confirmed the monophyly of Pentatomidae. However, our analysis challenged the monophyly of Pentatominae and Podopinae. We also detected unexpected relationships among some tribes and genera, highlighting the complexity of the internal taxonomic structure of Pentatomidae. Both Diceraeus species were grouped in the same clade with the remaining Carpocorini analyzed.

The complete mitochondrial genome of Potomida acarnanica (Kobelt, 1879)

Freshwater mussels (Bivalvia, Unionida) play essential roles in the well-functioning of ecosystems, even providing essential services to humans. However, these bivalves face numerous threats (e.g. habitat loss and fragmentation, pollution, introduction of invasive species, and climate change) which have already led to the extinction of many populations. This underscores the need to fully characterize the biology of these species, particularly those, such as Potomida acarnanica, that are still poorly studied. This study presents the first mitogenome of P. acarnanica (Kobelt, 1879), an endemic species of Greece with a distribution limited to only two river basins. The mitochondrial genome of a P. acarnanica specimen, collected at Pamisos River (Peloponnese, Greece), was sequenced by Illumina high-throughput sequencing. This mitogenome (16,101 bp) is characterized by 13 protein-coding genes, 22 transfer RNA and 2 ribosomal RNA genes. The size of this mitogenome is within the range of another Potomida mitogenome already published for the species Potomida littoralis. In the phylogenetic inference, P. acarnanica was recovered as monophyletic with P. littoralis mitogenome in the Lamprotulini tribe, as expected. This genomic resource will assist in genetically characterizing the species, potentially benefiting future evolutionary studies and conservation efforts.

The First Genome-Wide Survey of Shortbelly Eel (Dysomma anguillare Barnard, 1923) to Provide Genomic Characteristics, Microsatellite Markers and Complete Mitogenome Information

Dysomma anguillare is a demersal eel widespread distributing in tropical waters of the Indo-West Pacific and Atlantic. As an important component of the coastal fishery and marine ecosystem, the lack of genomic information for this species severely restricts the progress of relevant researches. In this study, the abecedarian genome-wide characteristics and phylogenetic relationships analyses were carried out based on next-generation sequencing (NGS) platform. The revised genome size was approximately 1 310 Mb, with the largest scaffold length reaching 23 878 bp through K-mer (K = 17) method. The heterozygosity, repetitive rate and average GC content were about 0.94%, 51.93% and 42.23%, respectively. A total of 1 160 104 microsatellite motifs were identified from the de novo assembled genome of D. anguillare, in which dinucleotide repeats accounted for the largest proportion (592 234, 51.05%), the highest occurrence frequency (14.58%) as well as the largest relative abundance (379.27/Mb). The high-polymorphic and moderate-polymorphic loci composed around 73% of the total single sequence repeats (SSRs), showing a latent capacity for subsequent population genetic structure and genetic diversity appraisal researches. Another byproduct of whole-genome sequencing, the double-stranded and circular mitogenome (16 690 bp) was assembled to investigate the evolutionary relationships of D. anguillare. The phylogenic tree constructed with maximum likelihood (ML) method showed that D. anguillare was closely related to Synaphobranchidae species, and the molecular systematic results further supported classical taxonomy status of D. anguillare.

Characterization and Molecular Phylogenetic Analysis of Subfamily Erebinae (Lepidoptera: Noctuoidea: Erebidae) Using Five Complete Mitochondrial Genomes

In this study, the complete mitogenomes of Sympis rufibasis, Lacera noctilio, Oxyodes scrobiculata, Mocis undata, and Artena dotata were newly sequenced to bring up-to-date the database using the next-generation sequencing methods. The gene order of all sequenced mitogenomes was identical consisting of 13 protein-coding genes (PCGs), 22 transfer RNA genes (tRNAs), two ribosomal RNA genes (rRNAs), and a non-coding A+T-rich region, which were common to other Lepidopteran insects. All protein-coding genes (PCGs) initiated with a canonical ATN codon and ended with TAN or an incomplete stop codon, single T. The A+T-rich region of S. rufibasis, L. noctilio, O. scrobiculata, M. undata, and A. dotata are 406 bp, 462 bp, 372 bp, 410 bp, and 406 bp long, respectively, containing number of characteristics that are distinctive to Noctuoidea moths. We analyzed concatenated amino acid sequences of protein-coding genes not including rRNAs, using Maximum Likelihood and Bayesian Inference methods. The phylogenetic analyses indicated that the tribe relationships within Erebinae were reconstructed as (Sypnini+((Erebini 1+Poaphilini 1)+((Euclidiini+Catocalini+(Hypopyrini+Erebini 2))+((Hulodini+(Poaphilini 2+Ophiusini))))). Phylogenetic analyses supported and confirmed the monophyly of the subfamilies' relationships as follows: (Hypeninae+Lymantriinae)+((Scoliopterginae+((Calpinae+Erebinae)+((Herminiinae+Aganainae)+Arctiinae)))) within Erebidae.

Refining the phylogeny and taxonomy of the apple tribe Maleae (Rosaceae): insights from phylogenomic analyses of 563 plastomes and a taxonomic synopsis of Photinia and its allies in the Old World

This study addresses the longstanding absence of a comprehensive phylogenetic backbone for the apple tribe Maleae, a deficiency attributed to limited taxon and marker sampling. We conducted an extensive taxon sampling, incorporating 563 plastomes from a diverse range of 370 species encompassing 26 presently recognized genera. Employing a range of phylogenetic inference methods, including RAxML and IQ-TREE2 for Maximum Likelihood (ML) analyses, we established a robust phylogenetic framework for the Maleae tribe. Our phylogenomic investigations provided compelling support for three major clades within Maleae. By integrating nuclear phylogenetic data with morphological and chromosomal evidence, we propose an updated infra-tribal taxonomic system, comprising subtribe Malinae Reveal, subtribe Lindleyinae Reveal, and subtribe Vauqueliniinae B.B.Liu (subtr. nov.). Plastid phylogenetic analysis also confirmed the monophyly of most genera, except for Amelanchier, Malus, Sorbus sensu lato, and Stranvaesia. In addition, we present a comprehensive taxonomic synopsis of Photinia and its morphological allies in the Old World, recognizing 27 species and ten varieties within Photinia, three species and two varieties within Stranvaesia, and two species and three varieties within Weniomeles. Furthermore, we also lectotypified 12 names and made two new combinations, Photiniamicrophylla (J.E.Vidal) B.B.Liu and Weniomelesatropurpurea (P.L.Chiu ex Z.H.Chen & X.F.Jin) B.B.Liu.

The first complete mitochondrional genome of Anopheles gibbinsi using a skimming sequencing approach

Mosquitoes belonging to the genus Anopheles are the only vectors of human malaria. Anopheles gibbinsi has been linked to malaria transmission in Kenya, with recent collections in Zambia reporting the mosquito species exhibiting zoophilic and exophilic behavioral patterns with occasional contact with humans. Given the paucity of genetic data, and challenges to identification and molecular taxonomy of the mosquitoes belonging to the Anopheles genus; we report the first complete mitochondrial genome of An. gibbinsi using a genome skimming approach. An Illumina Novaseq 6000 platform was used for sequencing, the length of the mitochondrial genome was 15401 bp, with 78.5% AT content comprised of 37 genes. Phylogenetic analysis by maximum likelihood using concatenation of the 13 protein coding genes demonstrated that An. marshallii was the closest relative based on existing sequence data. This study demonstrates that the skimming approach is an inexpensive and efficient approach for mosquito species identification and concurrent taxonomic rectification, which may be a useful alternative for generating reference sequence data for evolutionary studies among the Culicidae.

Which contributes more to the relict flora distribution pattern in East Asia, geographical processes or climate change? New evidence from the phylogeography of Rehderodendron kwangtungense

BACKGROUND

Relict species are important for enhancing the understanding of modern biogeographic distribution patterns. Although both geological and climatic changes since the Cenozoic have affected the relict flora in East Asia, the contributions of geographical processes remain unclear. In this study, we employed restriction-site associated DNA sequencing (RAD-seq) and shallow genome sequencing data, in conjunction with ecological niche modeling (ENM), to investigate the spatial genetic patterns and population differentiation history of the relict species Rehderodendron kwangtungense Chun.

RESULTS

A total of 138 individuals from 16 populations were collected, largely covering the natural distribution of R. kwangtungense. The genetic diversity within the R. kwangtungense populations was extremely low (HO = 0.048 ± 0.019; HE = 0.033 ± 0.011). Mantel tests revealed isolation-by-distance pattern (R2 = 0.38, P < 0.001), and AMOVA analysis showed that the genetic variation of R. kwangtungense occurs mainly between populations (86.88%, K = 7). Between 23 and 21 Ma, R. kwangtungense underwent a period of rapid differentiation that coincided with the rise of the Himalayas and the establishment of the East Asian monsoon. According to ENM and population demographic history, the suitable area and effective population size of R. kwangtungense decreased sharply during the glacial period and expanded after the last glacial maximum (LGM).

CONCLUSION

Our study shows that the distribution pattern of southern China mountain relict flora may have developed during the panplain stage between the middle Oligocene and the early Miocene. Then, the flora later fragmented under the force of orogenesis, including intermittent uplift during the Cenozoic Himalayan orogeny and the formation of abundant rainfall associated with the East Asian monsoon. The findings emphasized the predominant role of geographical processes in shaping relict plant distribution patterns.

Sequence comparison of the mitochondrial genomes of Plesionika species (Caridea: Pandalidae), gene rearrangement and phylogenetic relationships of Caridea

Background

Despite the Caridean shrimps' vast species richness and ecological diversity, controversies persist in their molecular classification. Within Caridea, the Pandalidae family exemplifies significant taxonomic diversity. As of June 25, 2023, GenBank hosts only nine complete mitochondrial genomes (mitogenomes) for this family. The Plesionika genus within Pandalidae is recognized as polyphyletic. To improve our understanding of the mitogenome evolution and phylogenetic relationships of Caridea, this study introduces three novel mitogenome sequences from the Plesionika genus: P.  ortmanni, P. izumiae and P. lophotes.

Methods

The complete mitochondrial genomes of three Plesionika species were sequenced utilizing Illumina's next-generation sequencing (NGS) technology. After assembling and annotating the mitogenomes, we conducted structural analyses to examine circular maps, sequence structure characteristics, base composition, amino acid content, and synonymous codon usage frequency. Additionally, phylogenetic analysis was performed by integrating existing mitogenome sequences of true shrimp available in GenBank.

Results

The complete mitogenomes of the three Plesionika species encompass 37 canonical genes, comprising 13 protein-coding genes (PCGs), 22 transfer RNAs (tRNAs), two ribosomal RNAs (rRNAs), and one control region (CR). The lengths of these mitogenomes are as follows: 15,908 bp for P. ortmanni, 16,074 bp for P. izumiae and 15,933 bp for P. lophotes. Our analyses extended to their genomic features and structural functions, detailing base composition, gene arrangement, and codon usage. Additionally, we performed selection pressure analysis on the PCGs of all Pandalidae species available in Genbank, indicating evolutionary purification selection acted on the PCGs across Pandalidae species. Compared with the ancestral Caridea, translocation of two tRNA genes, i.e., trnP or trnT, were found in the two newly sequenced Plesionika species-P. izumiae and P. lophotes. We constructed a phylogenetic tree of Caridea using the sequences of 13 PCGs in mitogenomes. The results revealed that family Pandalidae exhibited robust monophyly, while genus Plesionika appeared to be a polyphyletic group.

Conclusions

Gene rearrangements within the Pandalidae family were observed for the first time. Furthermore, a significant correlation was discovered between phylogenetics of the Caridea clade and arrangement of mitochondrial genes. Our findings offer a detailed exploration of Plesionika mitogenomes, laying a crucial groundwork for subsequent investigations into genetic diversity, phylogenetic evolution, and selective breeding within this genus.

The complete chloroplast genome of Illicium simonsii Maxim. (Illiciaceae), a species with important medicinal properties

Illicium simonsii Maxim (1888) is a medicinal species of the genus Illicium in the Illiciaceae family. It is commonly used to cure gastro-frigid vomiting, cystic hernia, gas pains in the chest, and scabies as folk medicine. To utilize its resources efficiently, the complete chloroplast genome of I. simonsii was sequenced, assembled, and annotated by using high-throughput sequencing data. The complete chloroplast genome was 143,038 bp in length, with a large single-copy region (LSC) of 101,094 bp, a short single-copy region (SSC) of 20,070 bp, and a pair of inverted repeats (IRs) of 21,874 bp. A total of 113 genes were annotated, including 79 protein-coding genes, 30 tRNA genes, and four rRNA genes. The phylogenetic tree exhibited that I. simonsii and Illicium burmanicum form a sister group, and were nested in the monophyletic clade of the Illicium genus.

Mitochondrial genomic characteristics and phylogenetic analysis of a brewing fungus, Rhizopus microsporus Tiegh. 1875 (Mucorales: Rhizopodaceae)

Rhizopus microsporus Tiegh. 1875 is widely used in a variety of industries, such as brewing, wine making, baking, and medicine production, as it has the capability to break down proteins and generate surface-active agents. To date, the mitochondrial genome features of early evolved fungi from the Rhizopus genus have not been extensively studied. Our research obtained a full mitochondrial genome of R. microsporus species, which was 43,837 bp in size and had a GC content of 24.93%. This genome contained 14 core protein-coding genes, 3 independent ORFs, 7 intronic ORFs, 24 tRNAs, and 2 rRNA genes. Through the use of the BI phylogenetic inference method, we were able to create phylogenetic trees for 25 early differentiation fungi which strongly supported the major clades; this indicated that R. microsporus is most closely related to Rhizopus oryzae.

Comparative and phylogenetic analyses of Loranthaceae plastomes provide insights into the evolutionary trajectories of plastome degradation in hemiparasitic plants

BACKGROUND

The lifestyle transition from autotrophy to heterotrophy often leads to extensive degradation of plastomes in parasitic plants, while the evolutionary trajectories of plastome degradation associated with parasitism in hemiparasitic plants remain poorly understood. In this study, phylogeny-oriented comparative analyses were conducted to investigate whether obligate Loranthaceae stem-parasites experienced higher degrees of plastome degradation than closely related facultative root-parasites and to explore the potential evolutionary events that triggered the 'domino effect' in plastome degradation of hemiparasitic plants.

RESULTS

Through phylogeny-oriented comparative analyses, the results indicate that Loranthaceae hemiparasites have undergone varying degrees of plastome degradation as they evolved towards a heterotrophic lifestyle. Compared to closely related facultative root-parasites, all obligate stem-parasites exhibited an elevated degree plastome degradation, characterized by increased downsizing, gene loss, and pseudogenization, thereby providing empirical evidence supporting the theoretical expectation that evolution from facultative parasitism to obligate parasitism may result in a higher degree of plastome degradation in hemiparasites. Along with infra-familial divergence in Loranthaceae, several lineage-specific gene loss/pseudogenization events occurred at deep nodes, whereas further independent gene loss/pseudogenization events were observed in shallow branches.

CONCLUSIONS

The findings suggest that in addition to the increasing levels of nutritional reliance on host plants, cladogenesis can be considered as another pivotal evolutionary event triggering the 'domino effect' in plastome degradation of hemiparasitic plants. These findings provide new insights into the evolutionary trajectory of plastome degradation in hemiparasitic plants.

An atlas of spider development at single-cell resolution provides new insights into arthropod embryogenesis

Spiders are a diverse order of chelicerates that diverged from other arthropods over 500 million years ago. Research on spider embryogenesis, particularly studies using the common house spider Parasteatoda tepidariorum, has made important contributions to understanding the evolution of animal development, including axis formation, segmentation, and patterning. However, we lack knowledge about the cells that build spider embryos, their gene expression profiles and fate. Single-cell transcriptomic analyses have been revolutionary in describing these complex landscapes of cellular genetics in a range of animals. Therefore, we carried out single-cell RNA sequencing of P. tepidariorum embryos at stages 7, 8 and 9, which encompass the establishment and patterning of the body plan, and initial differentiation of many tissues and organs. We identified 20 cell clusters, from 18.5 k cells, which were marked by many developmental toolkit genes, as well as a plethora of genes not previously investigated. We found differences in the cell cycle transcriptional signatures, suggestive of different proliferation dynamics, which related to distinctions between endodermal and some mesodermal clusters, compared with ectodermal clusters. We identified many Hox genes as markers of cell clusters, and Hox gene ohnologs were often present in different clusters. This provided additional evidence of sub- and/or neo-functionalisation of these important developmental genes after the whole genome duplication in an arachnopulmonate ancestor (spiders, scorpions, and related orders). We also examined the spatial expression of marker genes for each cluster to generate a comprehensive cell atlas of these embryonic stages. This revealed new insights into the cellular basis and genetic regulation of head patterning, hematopoiesis, limb development, gut development, and posterior segmentation. This atlas will serve as a platform for future analysis of spider cell specification and fate, and studying the evolution of these processes among animals at cellular resolution.

The complete chloroplast genome of Durio zibethinus L. cultivar Ri6 (Helicteroideae, Malvaceae)

Durian, a member of the Malvaceae family, is famous for its delicious fruits, which have strong scents and are rich in nutrients. In this study, we sequenced and characterized the complete chloroplast genome of Durio zibethinus L. 1774 cultivar Ri6, a popular durian cultivar in Vietnam, using the Illumina Hiseq platform. The results showed a circular chloroplast genome composed of a large single copy of 96,115 bp, a small single copy of 20,819 bp, and two inverted repeat regions of 24,185 bp. This genome consisted of 79 protein-coding genes, 30 transfer RNA genes, and four ribosomal RNA genes. The overall GC content of this genome was 35.7%. Phylogenetic analysis inferred from 78 protein-coding regions revealed monophyly of Durio species and a close relationship between D. zibethinus cultivar Ri6 and cultivar Mongthong. This study provides essential information for further studies examining genetic population, breedings, and species identification among Durio taxa and cultivars.

Chromosome-level genome assembly of the European green woodpecker Picus viridis

The European green woodpecker, Picus viridis, is a widely distributed species found in the Western Palearctic region. Here, we assembled a highly contiguous genome assembly for this species using a combination of short- and long-read sequencing and scaffolded with chromatin conformation capture (Hi-C). The final genome assembly was 1.28 Gb and features a scaffold N50 of 37 Mb and a scaffold L50 of 39.165 Mb. The assembly incorporates 89.4% of the genes identified in birds in OrthoDB. Gene and repetitive content annotation on the assembly detected 15,805 genes and a ∼30.1% occurrence of repetitive elements, respectively. Analysis of synteny demonstrates the fragmented nature of the P. viridis genome when compared to the chicken (Gallus gallus). The assembly and annotations produced in this study will certainly help for further research into the genomics of P. viridis and the comparative evolution of woodpeckers. Five historical and seven contemporary samples have been resequenced and may give insights on the population history of this species.

Characterization of the complete chloroplast genome of Helicteres hirsuta Lour. 1790 (Helicteriodeae: Malvaceae)

Helicteres hirsuta Lour. 1790 is a precious medicinal plant species, especially for treating chronic liver diseases. Genomic data on H. hirsuta are limited. Therefore, this current study aimed to characterize the chloroplast genome of H. hirsuta and reconstruct the phylogenetic relationship among Helicteroideae taxa. Consequently, the complete chloroplast genome of H. hirsuta was 163,404 bp in length and contained 113 unique genes (79 protein-coding genes, 30 tRNA genes, and four rRNA genes). Notably, two introns of clpP gene of H. hirsuta were lost in comparison to that of other Helicteroideae species. The phylogenetic tree based on chloroplast genomes of eleven Helicteroideae species revealed that H. hirsuta was closely related to Reevesia species. In conclusion, our study described the first complete chloroplast genome of H. hirsuta, which is essential for tracing evolutionary history in the Helicteroideae subfamily.