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

OGU: A Toolbox for Better Utilising Organelle Genomic Data

Organelle genomes serve as crucial datasets for investigating the genetics and evolution of plants and animals, genome diversity, and species identification. To enhance the collection, analysis, and visualisation of such data, we have developed a novel open-source software tool named Organelle Genome Utilities (OGU). The software encompasses three modules designed to streamline the handling of organelle genome data. The data collection module is dedicated to retrieving, validating and organising sequence information. The evaluation module assesses sequence variance using a range of methods, including novel metrics termed stem and terminal phylogenetic diversity. The primer module designs universal primers for downstream applications. Finally, a visualisation pipeline has been developed to present comprehensive insights into organelle genomes across different lineages rather than focusing solely on individual species. The performance, compatibility and stability of OGU have been rigorously evaluated through benchmarking with four datasets, including one million mixed GenBank records, plastid genomic data from the Lamiaceae family, mitochondrial data from rodents, and 308 plastid genomes sourced from various angiosperm families. Based on software capabilities, we identified 30 plastid intergenic spacers. These spacers exhibit a moderate evolutionary rate and offer practical utility comparable to coding regions, highlighting the potential applications of intergenic spacers in organelle genomes. We anticipate that OGU will substantially enhance the efficient utilisation of organelle genomic data and broaden the prospects for related research endeavours.

Comparative analysis of chloroplast genome of Lonicera japonica cv. Damaohua

Lonicera japonica is a well-known medicinal plant, and the Damaohua cultivar is one of the oldest known honeysuckle cultivars in China. The 155,151 bp chloroplast genome of this cultivar was obtained through Illumina sequencing. The genome includes a pair of inverted repeats (IRa and IRb; 23,789 bp each), a large single-copy region (88,924 bp), and a small single-copy (SSC) region (18,649 bp). In total, 127 unique genes were identified: 80 protein-coding, 39 tRNA, and 8 rRNA genes. Only ycf3 contained two introns. Eighty-nine large repetitive sequences and 54 simple sequence repeats were detected. Fifty potential RNA editing sites were predicted. Adaptive evolution analysis revealed that infA, matK, petB, petD, rbcL, rpl16, rpl2, rps3, ycf1, and ycf2 were positively selected, possibly reflecting the specific environmental adaptations of this cultivar. Sequence alignment and analysis revealed several candidate fragments for Lonicera species identification, such as the intergenic regions rpoB-petN, rbcL-accD, and psaA-ycf3. The IR region boundary and phylogenetic analysis revealed that the L. japonica cv. Damaohua chloroplast genome was most closely related to the L. japonica genome, but there were five distinct differences between the two. There are four sites with high variability between L. japonica and L. japonica cv. Damaohua with nucleotide variability (Pi) greater than 0.002, including rps2-rpoC2, atpB-rbcL, ycf1, and ycf1-trnN GUU. The differences between L. japonica and L. japonica cv. Damaohua were further confirmed by the single nucleotide polymorphism sites between these two species. Therefore, this study revealed that the chloroplast genome can serve as a universal super barcode for plant identification, which can identify differences and help distinguish Lonicera japonica from related species. An understanding of Lonicera japonica cv. Damaohua chloroplast genomics and a comparative analysis of Lonicera species will provide a scientific basis for breeding, species identification, systematic evolution analysis, and chloroplast genetic engineering research on medicinal honeysuckle plants.

Deciphering the complex organelle genomes of two Rhododendron species and insights into adaptive evolution patterns in high-altitude

BACKGROUND

The genomes within organelles are crucial for physiological functions such as respiration and photosynthesis and may also contribute to environmental adaptation. However, the limited availability of genetic resources, particularly mitochondrial genomes, poses significant challenges for in-depth investigations.

RESULTS

Here, we explored various assembly methodologies and successfully reconstructed the complex organelle genomes of two Rhododendron species: Rhododendron nivale subsp. boreale and Rhododendron vialii. The mitogenomes of these species exhibit various conformations, as evidenced by long-reads mapping. Notably, only the mitogenome of R. vialii can be depicted as a singular circular molecule. The plastomes of both species conform to the typical quadripartite structure but exhibit elongated inverted repeat (IR) regions. Compared to the high similarity between plastomes, the mitogenomes display more obvious differences in structure, repeat sequences, and codon usage. Based on the analysis of 58 organelle genomes from angiosperms inhabiting various altitudes, we inferred the genetic adaptations associated with high-altitude environments. Phylogenetic analysis revealed partial inconsistencies between plastome- and mitogenome-derived phylogenies. Additionally, evolutionary lineage was determined to exert a greater influence on codon usage than altitude. Importantly, genes such as atp4, atp9, mttB, and clpP exhibited signs of positive selection in several high-altitude species, suggesting a potential link to alpine adaptation.

CONCLUSIONS

We tested the effectiveness of different organelle assembly methods for dealing with complex genomes, while also providing and validating high-quality organelle genomes of two Rhododendron species. Additionally, we hypothesized potential strategies for high-altitude adaptation of organelles. These findings offer a reference for the assembly of complex organelle genomes, while also providing new insights and valuable resources for understanding their adaptive evolution patterns.

The complete mitochondrial genome of Periphyllus koelreuteriae (Takahashi, 1919) (Hemiptera: Aphididae: Chaitophorinae)

We sequenced the complete mitochondrial genome of the aphid Periphyllus koelreuteriae. It is 16,828 bp long and includes 13 protein-coding genes (PCGs), 22 transfer RNA genes (tRNAs), two ribosome RNA genes (rRNAs), a control region, and a repeat region. The control region contains a 273-bp repeat motif repeated 2.3 times. The species-specific repeat region between trnE and trnF comprises two 340-bp repeat units. The maximum-likelihood tree based on all 37 mitochondrial genes indicated a close relationship between P. koelreuteriae and Periphyllus diacerivorus. This study provides a valuable mitogenomic resource for future research on Chaitophorinae phylogeny and P. koelreuteriae diversification.

The complete genome sequence of the chloroplast of Bidens aurea

Bidens aurea (Asteraceae), a native of tropical America is now widespread in Asia and the Americas. We explored the B. aurea chloroplast genome and conducted a phylogenetic analysis. The chloroplast genome was circular, consisting of a large single copy (LSC) of 83,909 base pairs (bp), a small single copy (SSC) of 18,407 bp, and two inverted repeat regions (IR) of 24,729 bp each. Phylogenetic analysis showed that the 19 Bidens taxa were divided into five major clades, and B. aurea was most closely related to two species. Our findings offer a high-quality B. aurea chloroplast genome, aiding DNA barcode development and evolutionary history reconstruction.

New mitochondrial genomes of three whip spider species from the Amazon (Arachnida, Amblypygi) with phylogenetic relationships and comparative analysis

The complete mitochondrial genomes of the whip spiders Charinus carajas, C. ferreus, and Heterophrynus longicornis were sequenced, annotated, and compared with other mitogenomes of whip spiders and arachnids. The three new mitogenomes have the 37 genes usually observed in Metazoa: 13 protein-coding genes (PCGs), 22 transfer RNAs (tRNAs), and two ribosomal RNAs (rRNAs), plus a non-coding control region (CR). Most PCGs presented an ATN start codon, except cox1 in both Charinus species, initiating with TTA. Most PCGs terminated with stop codons TAA or TAG, except nad5 of C. carajas and cox3 of H. longicornis, which presented an incomplete stop codon (T). The Ka/Ks ratios were less than one for all the PCGs, indicating these genes are under purifying selection. All the tRNAs, except for serine 1 (trnS1), had the typical cloverleaf-shaped secondary structure. All the phylogenetic analyses resolved Charinus carajas and C. ferreus as monophyletic groups. Nonetheless, we did not recover the monophyly of Heterophrynus longicornis. The phylogenies under partitioned models did not recover suprageneric taxonomic groups as clades, but the Bayesian inference under the CAT infinite mixture model recovered the family Phrynidae and the superfamily Phrynoidea as monophyletic groups.

Population Genomics of Adaptive Radiations: Exceptionally High Levels of Genetic Diversity and Recombination in an Endemic Spider From the Canary Islands

The spider genus Dysdera has undergone a remarkable diversification in the oceanic archipelago of the Canary Islands, with ~60 endemic species having originated during the 20 million years since the origin of the archipelago. This evolutionary radiation has been accompanied by substantial dietary shifts, often characterised by phenotypic modifications encompassing morphological, metabolic and behavioural changes. Hence, these endemic spiders represent an excellent model for understanding the evolutionary drivers and to pinpoint the genomic determinants underlying adaptive radiations. Recently, we achieved the first chromosome-level genome assembly of one of the endemic species, D. silvatica, providing a high-quality reference sequence for evolutionary genomics studies. Here, we conducted a low coverage-based resequencing study of a natural population of D. silvatica from La Gomera island. Taking advantage of the new high-quality genome, we characterised genome-wide levels of nucleotide polymorphism, divergence and linkage disequilibrium, and inferred the demographic history of this population. We also performed comprehensive genome-wide scans for recent positive selection. Our findings uncovered exceptionally high levels of nucleotide diversity and recombination in this geographically restricted endemic species, indicative of large historical effective population sizes. We also identified several candidate genomic regions that are potentially under positive selection, highlighting relevant biological processes, such as vision and nitrogen extraction as potential adaptation targets. These processes may ultimately drive species diversification in this genus. This pioneering study of spiders that are endemic to an oceanic archipelago lays the groundwork for broader population genomics analyses aimed at understanding the genetic mechanisms driving adaptive radiation in island ecosystems.

Phylogenetics and biogeography of the olive family (Oleaceae)

BACKGROUND AND AIMS

Progress in the systematic studies of the olive family (Oleaceae) during the last two decades provides the opportunity to update its backbone phylogeny and to investigate its historical biogeography. We also aimed to understand the factors underlying the disjunct distribution pattern between East Asia and both West Asia and Europe that is found more commonly in this family than in any other woody plant family.

METHODS

Using a sampling of 298 species out of ~750, the largest in a phylogenetic study of Oleaceae thus far, with a set of 36 plastid and nuclear markers, we reconstructed and dated a new phylogenetic tree based on maximum likelihood and Bayesian methods and checked for any reticulation events. We also assessed the relative support of four competing hypotheses [Qinghai-Tibet Plateau uplift (QTP-only hypothesis); climatic fluctuations (climate-only hypothesis); combined effects of QTP uplift and climate (QTP-climate hypothesis); and no effects (null hypothesis)] in explaining these disjunct distributions.

KEY RESULTS

We recovered all tribes and subtribes within Oleaceae as monophyletic, but uncertainty in the position of tribe Forsythieae remains. Based on this dataset, no reticulation event was detected. Our biogeographical analyses support the QTP-climate hypothesis as the likely main explanation for the East-West Eurasian disjunctions in Oleaceae. Our results also show an earlier origin of Oleaceae at ~86 Mya and the role of Tropical Asia as a main source of species dispersals.

CONCLUSION

Our new family-wide and extensive phylogenetic tree highlights both the stable relationships within Oleaceae, including the polyphyly of the genus Chionanthus, and the need for further systematic studies within the largest and most undersampled genera of the family (Chionanthus and Jasminum). Increased sampling will also help to fine-tune biogeographical analyses across spatial scales and geological times.

The Complete Mitochondrial Genome of Nephropsis grandis: Insights into the Phylogeny of Nephropidae Mitochondrial Genome

The systematic phylogeny of Pleocyemata species, particularly within the family Nephropidae, remains incomplete. In order to enhance the taxonomy and systematics of Nephropidae within the evolutionary context of Pleocyemata, we embarked upon a comprehensive study aiming to elucidate the phylogenetic position of Nephropsis grandis. Consequently, we determined the complete mitochondrial DNA sequence for N. grandis. The circular genome spans a length of 15,344 bp and exhibits a gene composition analogous to that observed in other metazoans, encompassing a comprehensive set of 37 genes. Additionally, the genome features an AT-rich region. The rRNAs exhibited the highest AT content among the 37 genes (70.41%), followed by tRNAs (67.42%) and protein-coding genes (PCGs) (62.76%). The absence of a dihydrouracil arm in trnS1 prevented the formation of the canonical cloverleaf secondary structure. Selective pressure analysis indicated that the PCGs underwent purifying selection. The Ka/Ks ratios for cox1, cox2, cox3, and cob were considerably lower compared to other PCGs, implying strong purifying selection acting upon these particular genes. The mitochondrial gene order in N. grandis was consistent with the reported order in ancestral Pleocyemata. Phylogenetic revealed that N. grandis forms a cluster with the genus Metanephrops, and this cluster further groups with Homarus and the genus Nephrops within the Nephropidae family. These findings provide robust support for N. grandis as an ancestral member of the Nephropidae family. This study highlights the significance of employing complete mitochondrial genomes in phylogenetic analysis and deepens our understanding of the evolution of the Nephropidae family.

The complete chloroplast genome of Erodium cicutarium (Linnaeus) l' Héritier ex Aiton 1789 (Geraniaceae): genome characterization and phylogenetic consideration

Erodium cicutarium is an annual herbaceous plant valued for its applications in traditional medicine. However, the chloroplast genome of E. cicutarium has yet to be reported. In this study, we assembled chloroplast genomes of Erodium cicutarium using Illumina sequencing reads. The chloroplast genome was 114,652 bp long, harbored 111 complete genes, and its overall GC content was 39.1%. In Maximum Likelihood (ML) and Bayesian inference (BI) trees, the 13 Erodium species divided into three main clades, with E. cicutarium and E. carvifolium forming a monophyletic group, suggesting a close relationship between the two species. The E. cicutarium cp genome presented in this study lays a good foundation for the Erodium.

The complete mitochondrial genome of the shining leaf chafer Mimela junii (Duftschmidt, 1805) (Coleoptera: Scarabaeidae)

The complete mitochondrial genome of the shining leaf chafer Mimela junii was sequenced and is herein described. The mitogenome consists of a circular molecule of 16,805 bp, with an overall AT content of 75.7%. It encodes for 13 protein-coding genes (PCGs), 22 transfer RNA genes (tRNAs), two ribosomal RNA genes (rRNAs) and contains a non-coding Control Region (CR) characterized by the presence of tandem repeats. The gene order corresponds to the ancestral Pancrustacea model and mitogenome characteristics are congruous with those of hexapods. In the phylogenetic analysis, M. junii is nested within a paraphyletic Anomala with high support, and is herein associated with Anomala corpulenta with medium/low support.

New mitochondrial gene order arrangements and evolutionary implications in the class Octocorallia

The complete mitochondrial genomes of octocorals typically range from 18.5 kb to 20.5 kb in length and include 14 protein-coding genes (PCGs), two ribosomal RNA genes and one tRNA. To date, seven different gene orders (A-G) have been described, yet comprehensive investigations of the actual number of arrangements, as well as comparative analyses and evolutionary reconstructions of mitochondrial genome evolution within the whole class Octocorallia, have been often overlooked. Here, we considered the complete mitochondrial genomes available for octocorals and explored their structure and gene order variability. Our results updated the actual number of mitochondrial gene order arrangements so far known for octocorals from 7 to 14 and allowed us to explore and preliminarily discuss the role of some of the structural and functional factors in the mitogenomes. We performed comparative mitogenomic analyses on the existing and novel octocoral gene orders, considering different mitogenomic structural features such as genome size, GC percentage, AT and GC skewness. The mitochondrial gene order history mapped on a recently published nuclear loci phylogeny showed that the most common rearrangement events in octocorals are inversions, inverted transpositions and transpositions. Furthermore, gene order rearrangement events were restricted only to some regions of the tree. Overall, different rearrangement events arose independently and from the ancestral and most common gene order, instead of being derived from other rearranged orders. Finally, our data demonstrate how the study of mitochondrial gene orders can be used to explore the evolution of octocorals and in some cases can be used to assess the phylogenetic placement of certain taxa.

The complete chloroplast genome sequence of Arytera littoralis (Sapindoideae, Sapindaceae)

Arytera littoralis Blume 1847 is an evergreen small tree belonging to the Sapindaceae family. It is distributed in South China to SE Asia and the Solomon Islands. However, the chloroplast genome of A. littoralis has yet to be reported. In this study, the complete chloroplast genome of Arytera littoralis was determined. The total genome size was 161,091 bp in length, consisting of two inverted repeats (IRs) (28,432 bp) separated by the large single-copy (LSC) (85,737 bp) and small single-copy (SSC) (18,490 bp) regions. The genome contained 133 genes, including 87 protein-coding genes, 37 tRNA genes, and eight rRNA genes. The GC content of the complete chloroplast genome was 37.8%. A maximum-likelihood (ML) phylogenetic tree of A. littoralis and 13 related species from the family Sapindaceae indicated that A. littoralis was close to a clade composed of Sapindus, Nephelium, Litchi, and Dimocarpus. This study will offer essential genetic resources of A. littoralis and provide insights into the phylogeny and evolution of Sapindaceae.

Comparative analysis of the whole mitochondrial genomes of four species in sect. Chrysantha (Camellia L.), endemic taxa in China

BACKGROUND

The sect. Chrysantha Chang of plants with yellow flowers of Camellia species as the "Queen of the Tea Family", most of these species are narrowly distributed endemics of China and are currently listed Grde-II in National Key Protected Wild Plant of China. They are commercially important plants with horticultural medicinal and scientific research value. However, the study of the sect. Chrysantha species genetics are still in its infancy, to date, the mitochondrial genome in sect. Chrysantha has been still unexplored.

RESULTS

In this study, we provide a comprehensive assembly and annotation of the mitochondrial genomes for four species within the sect. Chrysantha. The results showed that the mitochondrial genomes were composed of closed-loop DNA molecules with sizes ranging from 850,836 bp (C. nitidissima) to 1,098,121 bp (C. tianeensis) with GC content of 45.71-45.78% and contained 48-58 genes, including 28-37 protein-coding genes, 17-20 tRNA genes and 2 rRNA genes. We also examined codon usage, sequence repeats, RNA editing and selective pressure in the four species. Then, we performed a comprehensive comparison of their basic structures, GC contents, codon preferences, repetitive sequences, RNA editing sites, Ka/Ks ratios, haplotypes, and RNA editing sites. The results showed that these plants differ little in gene type and number. C. nitidissima has the greatest variety of genes, while C. tianeensis has the greatest loss of genes. The Ka/Ks values of the atp6 gene in all four plants were greater than 1, indicating positive selection. And the codons ending in A and T were highly used. In addition, the RNA editing sites differed greatly in number, type, location, and efficiency. Twelve, six, five, and twelve horizontal gene transfer (HGT) fragments were found in C. tianeensis, Camellia huana, Camellia liberofilamenta, and C. nitidissima, respectively. The phylogenetic tree clusters the four species of sect. Chrysantha plants into one group, and C. huana and C. liberofilamenta have closer affinities.

CONCLUSIONS

In this study, the mitochondrial genomes of four sect. Chrysantha plants were assembled and annotated, and these results contribute to the development of new genetic markers, DNA barcode databases, genetic improvement and breeding, and provide important references for scientific research, population genetics, and kinship identification of sect. Chrysantha plants.

The complete chloroplast genome of Fuchsia standishii J. Harrison, 1840 (Onagraceae) from Yunnan, China

Fuchsia standishii J. Harrison, 1840, a perennial shrub, is renowned for its vividly colored and uniquely shaped blooms, which have an extended flowering season. Commonly cultivated as an ornamental potted plant, it is utilized in traditional Chinese medicine. In this study, we successfully sequenced and assembled the complete chloroplast genome of F. standishii using high-throughput Illumina sequencing technology. The assembled chloroplast genome displays a typical quadripartite structure, with a total length of 156,391 bp. It consists of a pair of inverted repeat regions (IRs), each measuring 25,069 bp, separated by a large single-copy region (LSC) of 87,754 bp and a small single-copy region (SSC) of 18,499 bp. The overall GC content of the genome is 37.60%. The genome includes a total of 129 genes, comprising 84 protein-coding genes, 37 tRNA genes, and eight rRNA genes. Phylogenetic analysis of 17 complete chloroplast genomes revealed that F. standishii forms a monophyletic group with the entire Circaea. This study provides a molecular foundation for future phylogenetic research on Fuchsia.

Complete chloroplast genome sequence of Artemisia littoricola (Asteraceae) from Dokdo Island Korea: genome structure, phylogenetic analysis, and biogeography study

The Asteraceae family, particularly the Artemisia genus, presents taxonomic challenges due to limited morphological characteristics and frequent natural hybridization. Molecular tools, such as chloroplast genome analysis, offer solutions for accurate species identification. In this study, we sequenced and annotated the chloroplast genome of Artemisia littoricola sourced from Dokdo Island, employing comparative analyses across six diverse Artemisia species. Our findings reveal conserved genome structures with variations in repeat sequences and junction boundaries. Notably, the chloroplast genome of A. littoricola spans 150,985 bp, consistent with other Artemisia species, and comprises 131 genes, including 86 protein-coding, 37 tRNA, and 8 rRNA genes. Among these genes, 16 possess a single intron, while clpP and ycf3 exhibit two introns each. Furthermore, 18 genes display duplicated copies within the IR regions. Moreover, the genome possesses 42 Simple Sequence Repeats (SSRs), predominantly abundant in A/T content and located within intergenic spacer regions. The analysis of codon usage revealed that the codons for leucine were the most frequent, with a preference for ending with A/U. While the chloroplast genome exhibited conservation overall, non-coding regions showed lower conservation compared to coding regions, with the Inverted Repeat (IR) region displaying higher conservation than single-copy regions. Phylogenetic analyses position A. littoricola within subgenus Dracunculus, indicating a close relationship with A. scoparia and A. desertorum. Additionally, biogeographic reconstructions suggest ancestral origins in East Asia, emphasizing Mongolia, China (North East and North Central and South Central China), and Korea. This study underscores the importance of chloroplast genomics in understanding Artemisia diversity and evolution, offering valuable insights into taxonomy, evolutionary patterns, and biogeographic history. These findings not only enhance our understanding of Artemisia's intricate biology but also contribute to conservation efforts and facilitate the development of molecular markers for further research and applications in medicine and agriculture.

An acidophilic fungus promotes prey digestion in a carnivorous plant

Leaves of the carnivorous sundew plants (Drosera spp.) secrete mucilage that hosts microorganisms, but whether this microbiota contributes to prey digestion is unclear. We identified the acidophilic fungus Acrodontium crateriforme as the dominant species in the mucilage microbial communities, thriving in multiple sundew species across the global range. The fungus grows and sporulates on sundew glands as its preferred acidic environment, and its presence in traps increased the prey digestion process. A. crateriforme has a reduced genome similar to other symbiotic fungi. During A. crateriforme-Drosera spatulata coexistence and digestion of prey insects, transcriptomes revealed significant gene co-option in both partners. Holobiont expression patterns during prey digestion further revealed synergistic effects in several gene families including fungal aspartic and sedolisin peptidases, facilitating prey digestion in leaves, as well as nutrient assimilation and jasmonate signalling pathway expression. This study establishes that botanical carnivory is defined by adaptations involving microbial partners and interspecies interactions.

Repeated migration, interbreeding and bottlenecking shaped the phylogeography of the selfing grass Brachypodium stacei

Brachypodium stacei is the most ancestral lineage in the genus Brachypodium, a model system for grass functional genomics. B. stacei shows striking and sometimes contradictory biological and evolutionary features, including a high selfing rate yet extensive admixture, an ancient Miocene origin yet with recent evolutionary radiation, and adaptation to different dry climate conditions in its narrow distribution range. Therefore, it constitutes an ideal system to study these life history traits. We studied the phylogeography of 17 native circum-Mediterranean B. stacei populations (39 individuals) using genome-wide RADseq SNP data and complete plastome sequences. Nuclear SNP data revealed the existence of six distinct genetic clusters, low levels of intra-population genetic diversity and high selfing rates, albeit with signatures of admixture. Coalescence-based dating analysis detected a recent split between crown lineages in the Late Quaternary. Plastome sequences showed incongruent evolutionary relationships with those recovered by the nuclear data, suggesting interbreeding and chloroplast capture events between genetically distant populations. Demographic and population dispersal coalescent models identified an ancestral origin of B. stacei in the western-central Mediterranean islands, followed by an early colonization of the Canary Islands and two independent colonization events of the eastern Mediterranean region through long-distance dispersal and bottleneck events as the most likely evolutionary history. Climate niche data identified three arid niches of B. stacei in the southern Mediterranean region. Our findings indicate that the phylogeography of B. stacei populations was shaped by recent radiations, frequent extinctions, long-distance dispersal events, occasional interbreeding, and adaptation to local climates.

Characterization of the complete mitochondrial genome of the medical fungus Ganoderma resinaceum Boud., 1889 (Polyporales: Ganodermataceae)

The medical mushroom Ganoderma resinaceum Boud., 1889, is of great interest in pharmacy due to its diverse functional active ingredients. However, the mitochondrial genome of G. resinaceum remains unexplored. Here, we present the complete mitochondrial genome of G. resinaceum, which spans 67,458 bp and has a GC content of 25.65%. This genome encompasses 15 core protein-coding genes, 8 independent ORFs, 15 intronic ORFs, 27 tRNAs, and 2 rRNA genes. Through phylogenetic analysis using Bayesian inference (BI), we elucidated the evolutionary relationships among 34 Basidiomycota fungi, revealing distinct clades and indicating a close relationship between G. resinaceum and G. subamboinense.

The complete chloroplast genomes and phylogenetic analysis of Exbucklandia longipetala and Exbucklandia populnea (Hamamelidaceae)

Exbucklandia longipetala and Exbucklandia populnea are two evergreen trees of the genus Exbucklandia in the family Hamamelidaceae. In this study, the complete chloroplast genomes of E. longipetala and E. populnea were sequenced, assembled, and annotated. The total lengths of the chloroplast genomes were 160,723 bp and 160,744 bp, respectively, and both had a GC content of 38.1%. The complete chloroplast genomes of these two species had typical quadripartite structures: LSC region (88,972 bp and 88,989 bp), SSC region (18,907 bp and 18,911 bp) and a pair of inverted repeats both of 26,422 bp. Both species contained 114 unique genes, including 80 protein-coding genes, 30 tRNA genes, and 4 rRNA genes. Phylogenetic analysis indicated that E. longipetala and E. populnea are sister species to each other. Our results provide useful genetic resources for further studies on the origin and evolution of Hamamelidaceae.

Multiple origins of freshwater invasion and parental care reflecting ancient vicariances in the bivalve family Cyrenidae (Mollusca)

Habitat transitions in living organisms are key innovations often coupled with species diversification after their successful adaptation to new environment. The Cyrenidae is among the most well-known heterodont bivalve groups that have successfully invaded freshwater systems from brackish water environments and display diverse lineage-specific developmental modes. Phylogenetic and molecular clock-based divergence time analyses using 12 complete mitochondrial genome sequences suggest that Cyrenidae species independently colonized freshwater habitats during three distinct spatial and geological periods: one from the American continents approximately in the Early Jurassic and the two others from Australasian/East Asian continents in the Early/Middle Cretaceous and the Paleogene-Neogene boundary, respectively. This study provides significant insight into the temporal and spatial patterns of multiple freshwater invasions, aligning with ancient vicariance events inferred from different geological timelines of plate tectonics. Additionally, mitogenome phylogeny confirms the earlier hypothesis of the repeated parallel evolution of parental care system within this bivalve group.

The mitochondrial genome sequence of Manchurian Hare (Lepus mandshuricus)

The Manchurian hare (Lepus mandshuricus) is widely distributed in eastern Russia and northeastern China, but due to limited research, its taxonomic status remains somewhat ambiguous. The mitochondrial genome of the Manchurian hare was 16,705 bp in length, which was consisted of 13 protein-coding genes (PCGs), 22 transfer RNA (tRNA) genes, 2 ribosomal RNA (rRNA) genes and one control region. The overall nucleotide composition is 31.7% A, 29.4% T, 13.3% G, and 25.6% C, indicating a high AT content. Phylogenetic analysis reveals a closer relationship of the Manchurian hare with the Korean hare (Lepus coreanus) and the Iberian hare (Lepus granatensis), while its relationship with the Hainan rabbit (Lepus hainanus), European hare (Lepus europaeus) and the snowshoe hare (Lepus americanus) is more distant. The mitochondrial genome of the Manchurian hare is of vital importance for the phylogenetic analysis of lagomorphs and provides valuable data for deeper evolutionary inquiries.

Species limits and hybridization in Andean leaf-eared mice (Phyllotis)

Leaf-eared mice (genus Phyllotis) are among the most widespread and abundant small mammals in the Andean Altiplano, but species boundaries and distributional limits are often poorly delineated due to sparse survey data from remote mountains and high-elevation deserts. Here we report a combined analysis of mitochondrial DNA variation and whole-genome sequence (WGS) variation in Phyllotis mice to delimit species boundaries, to assess the timescale of diversification of the group, and to examine evidence for interspecific hybridization. Estimates of divergence dates suggest that most diversification of Phyllotis occurred during the past 3 million years. Consistent with the Pleistocene Aridification hypothesis, our results suggest that diversification of Phyllotis largely coincided with climatically induced environmental changes in the mid- to late Pleistocene. Contrary to the Montane Uplift hypothesis, most diversification in the group occurred well after the major phase of uplift of the Central Andean Plateau. Species delimitation analyses revealed surprising patterns of cryptic diversity within several nominal forms, suggesting the presence of much undescribed alpha diversity in the genus. Results of genomic analyses revealed evidence of ongoing hybridization between the sister species Phyllotis limatus and P. vaccarum and suggest that the contemporary zone of range overlap between the two species represents an active hybrid zone.

Chloroplast genomes of Simarouba Aubl., molecular evolution and comparative analyses within Sapindales

Simarouba, a neotropical genus in the family Simaroubaceae, currently lacks comprehensive genomic data in existing databases. This study aims to fill this gap by providing genomic resources for three Simarouba species, S. amara, S. versicolor, and S. glauca. It also aims to perform comparative molecular evolutionary analyses in relation to other species within the order Sapindales. The analysis of these three Simarouba species revealed the presence of the typical quadripartite structure expected in plastomes. However, some pseudogenization events were identified in the psbC, infA, rpl22, and ycf1 genes. In particular, the CDS of the psbC gene in S. amara was reduced from 1422 bp to 584 bp due to a premature stop codon. Nucleotide diversity data pointed to gene and intergenic regions as promising candidates for species and family discrimination within the group, specifically matK, ycf1, ndhF, rpl32, petA-psbJ, and trnS-trnG. Selection signal analyses showed strong evidence for positive selection on the rpl23 gene. Phylogenetic analyses indicated that S. versicolor and S. glauca have a closer phylogenetic relationship than S. amara. We provide chloroplast genomes of three Simaruba species and use them to elucidate plastome evolution, highlight the presence of pseudogenization, and identify potential DNA barcode regions.

Complete mitochondrial genome sequence of the white root rot pathogen Dematophora necatrix (Xylariaceae: Xylariales)

The mitochondrial genome of Dematophora necatrix is 121,350 base pairs in length with a G + C content of 30.19%. Phylogenetic analysis showed that D. necatrix grouped with other members of the Xylariaceae, with which its mitogenome also shares a broadly similar architecture and gene content. The D. necatrix mitogenome contains 14 protein-coding and 26 tRNA-encoding genes, as well as one copy each of the rnl, rns, rps3 and nat1 genes. However, as much as 80% of this genome is intronic or non-coding. This is likely due to expansions and rearrangements caused by the large number of group I introns and the homing endonucleases and reverse-transcriptases they encode. Our study thus provides a valuable foundation from which to explore the mitochondrion's role in the biology of D. necatrix, and also serves as a resource for investigating the pathogen's population biology and general ecology.

Characterization and phylogenetic analysis of the Talaromyces liani (kamyschko) Yilmaz, Frisvad & Samson, 2014 (Eurotiales: trichocomaceae) mitochondrial genome

The filamentous fungus Talaromyces liani (Kamyschko) Yilmaz, Frisvad & Samson, 2014, has attracted considerable interest in biotechnology due to its diverse industrial applications and physiological characteristics. However, the mitochondrial genome of T. liani remains uncharacterized. Here, we present the complete mitochondrial genome of T. liani, comprising 38,000 bp with a GC content of 24.61%. This genome includes 15 core protein-coding genes, 4 independent ORFs, 6 intronic ORFs, 26 tRNAs, and 2 rRNA genes. Phylogenetic analysis using Bayesian inference (BI) revealed the evolutionary relationships among 15 fungi from Eurotiales, strongly supporting distinct clades and indicating that T. liani most closely related to T. pinophilus.

Comparative plastome assembly of the yellow ironweed (Verbesina alternifolia) using Nanopore and Illumina reads

Chloroplast genomes (plastomes) represent a very important source of valuable information for phylogenetic and biogeographic reconstructions. The use of short reads (as those produced from Illumina sequencing), along with de novo read assembly, has been considered the "gold standard" for plastome reconstruction. However, short reads often cannot reconstruct long repetitive regions in chloroplast genomes. Long Nanopore (ONT) reads can help bridging long repetitive regions but are by far more error-prone than those produced by Illumina sequencing. Verbesina is the largest genus of tribe Heliantheae (Asteraceae) and includes species of economic importance as ornamental or as invasive weeds. However, no complete chloroplast genomes have been published yet for the genus. We utilized Illumina and Nanopore sequencing data and different assembly strategies to reconstruct the plastome of Verbesina alternifolia and evaluated the usefulness of the Nanopore assemblies. The two plastome sequence assemblages, one obtained with the Nanopore sequencing and the other inferred with Illumina reads, were identical, except for missing bases in homonucleotide regions. The best-assembled plastome of V. alternifolia was 152,050 bp in length and contained 80, 29, and four unique protein-coding genes, tRNAs, and rRNAs, respectively. When used as reference for mapping Illumina reads, all plastomes performed similarly. In a phylogenetic analysis including 28 other plastomes from closely related taxa (from the Heliantheae alliance), the two Verbesina chloroplast genomes grouped together and were nested among the other members of the tribe Heliantheae s.str. Our study highlights the usefulness of the Nanopore technology for assembling rapidly and cost-effectively chloroplast genomes, especially in taxonomic groups with paucity of publicly available plastomes.

Characterisation of the complete chloroplast genome of Solanum tuberosum cv. White Lady

Potato (Solanum tuberosum) is considered worldwide as one of the most important non-cereal food crops. As a result of its adaptability and worldwide production area, potato displays a vast phenotypical variability as well as genomic diversity. Chloroplast genomes have long been a core issue in plant molecular evolution and phylogenetic studies, and have an important role in revealing photosynthetic mechanisms, metabolic regulations and the adaptive evolution of plants. We sequenced the complete chloroplast genome of the Hungarian cultivar White Lady, which is 155 549 base pairs (bp) in length and is characterised by the typical quadripartite structure composed of a large- and small single-copy region (85 991 bp and 18 374 bp, respectively) interspersed by two identical inverted repeats (25 592 bp). The genome consists of 127 genes of which 82 are protein-coding, eight are ribosomal RNAs and 37 are transfer RNAs. The overall gene content and distribution of the genes on the White Lady chloroplast was the same as found in other potato chloroplasts. The alignment of S. tuberosum chloroplast genome sequences resulted in a highly resolved tree, with 10 out of the 13 nodes recovered having bootstrap values over 90%. By comparing the White Lady chloroplast genome with available S. tuberosum sequences we found that gene content and synteny are highly conserved. The new chloroplast sequence can support further studies of genetic diversity, resource conservation, evolution and applied agricultural research. The new sequence can support further potato genetic diversity and evolutionary studies, resource conservation, and also applied agricultural research.

Characterization and trans-generation dynamics of mitogene pool in the silver carp (Hypophthalmichthys molitrix)

Multicopied mitogenome are prone to mutation during replication often resulting in heteroplasmy. The derived variants in a cell, organ, or an individual animal constitute a mitogene pool. The individual mitogene pool is initiated by a small fraction of the egg mitogene pool. However, the characteristics and relationship between them has not yet been investigated. This study quantitatively analyzed the heteroplasmy landscape, genetic loads, and selection strength of the mitogene pool of egg and hatchling in the silver carp (Hypophthalmichthys molitrix) using high-throughput resequencing. The results showed heteroplasmic sites distribute across the whole mitogenome in both eggs and hatchlings. The dominant substitution was Transversion in eggs and Transition in hatching accounting for 95.23%±2.07% and 85.38%±6.94% of total HP sites, respectively. The total genetic loads were 0.293±0.044 in eggs and 0.228±0.022 in hatchlings (P=0.048). The dN/dS ratio was 58.03±38.98 for eggs and 9.44±3.93 for hatchlings (P=0.037). These results suggest that the mitogenomes were under strong positive selection in eggs with tolerance to variants with deleterious effects, while the selection was positive but much weaker in hatchlings showing marked quality control. Based on these findings, we proposed a trans-generation dynamics model to explain differential development mode of the two mitogene pool between oocyte maturation and ontogenesis of offspring. This study sheds light on significance of mitogene pool for persistence of populations and subsequent integration in ecological studies and conservation practices.

The plastomes of Lepismium cruciforme (Vell.) Miq and Schlumbergera truncata (Haw.) Moran reveal tribe-specific rearrangements and the first loss of the trnT-GGU gene in Cactaceae

BACKGROUND

Recent studies have revealed atypical features in the plastomes of the family Cactaceae, the largest lineage of succulent species adapted to arid and semi-arid regions. Most plastomes sequenced to date are from short-globose and cylindrical cacti, while little is known about plastomes of epiphytic cacti. Published cactus plastomes reveal reduction and complete loss of IRs, loss of genes, pseudogenization, and even degeneration of tRNA structures. Aiming to contribute with new insights into the plastid evolution of Cactaceae, particularly within the tribe Rhipsalideae, we de novo assembled and analyzed the plastomes of Lepismium cruciforme and Schlumbergera truncata, two South American epiphytic cacti.

METHODS AND RESULTS

Our data reveal many gene losses in both plastomes and the first loss of functionality of the trnT-GGU gene in Cactaceae. The trnT-GGU is a pseudogene in L. cruciforme plastome and appears to be degenerating in the tribe Rhipsalideae. Although the plastome structure is conserved among the species of the tribe Rhipsalideae, with tribe-specific rearrangements, we mapped around 200 simple sequence repeats and identified nine nucleotide polymorphism hotspots, useful to improve the phylogenetic resolutions of the Rhipsalideae. Furthermore, our analysis indicated high gene divergence and rapid evolution of RNA editing sites in plastid protein-coding genes in Cactaceae.

CONCLUSIONS

Our findings show that some characteristics of the Rhipsalideae tribe are conserved, such as plastome structure with IRs containing only the ycf2 and two tRNA genes, structural degeneration of the trnT-GGU gene and ndh complex, and lastly, pseudogenization of rpl33 and rpl23 genes, both plastid translation-related genes.

Mitochondrial genome study of Camellia oleifera revealed the tandem conserved gene cluster of nad5-nads in evolution

Camellia oleifera is a kind of high-quality oil supply species. Its seeds contain rich unsaturated fatty acids and antioxidant active ingredients, which is a kind of high-quality edible oil. In this study, we used bioinformatics methods to decipher a hexaploid Camellia oil tree's mitochondrial (mt) genome based on second-generation sequencing data. A 709,596 bp circular map of C. oleifera mt genome was found for the first time. And 74 genes were annotated in the whole genome. Mt genomes of C. oleifera and three Theaceae species had regions with high similarity, including gene composition and gene sequence. At the same time, five conserved gene pairs were found in 20 species. In all of the mt genomes, most of nad genes existed in tandem pairs. In addition, the species classification result, which, according to the gene differences in tandem with nad5 genes, was consistent with the phylogenetic tree. These initial results provide a valuable basis for the further researches of Camellia oleifera and a reference for the systematic evolution of plant mt genomes.

Tracing Homopolymers in Oikopleura dioica's Mitogenome

Oikopleura dioica is a planktonic tunicate (Appendicularia class) found extensively across the marine waters of the globe. The genome of a single male individual collected from Okinawa, Japan was sequenced using the single-molecule PacBio Hi-Fi method and assembled with NOVOLoci. The mitogenome is 39,268 bp long, featuring a large control region of around 22,000 bp. We annotated the proteins atp6, cob, cox1, cox2, cox3, nad1, nad4, and nad5, and found one more open reading frame that did not match any known gene. This study marks the first complete mitogenome assembly for an appendicularian, and reveals that A and T homopolymers cumulatively account for nearly half of its length. This reference sequence will be an asset for environmental DNA and phylogenetic studies.

Characterization of the complete mitochondrial genome of a newly discovered torrent catfish, Liobagrus geumgangensis, and their phylogenetic relationships

BACKGROUND

A new Liobagrus fish was reported from the Korean Peninsula, but research on this taxon is lacking. Moreover, existing research on the mitogenome of the genus Liobagrus in Korea is very limited, and no studies have been conducted on structural characteristics of transfer RNA (tRNA) or gene order comparisons between taxa; instead, research has been restricted to basic phylogeny.

OBJECTIVE

The complete mitochondrial genome of Liobagrus geumgangensis was analyzed for the first time. We then aimed to reconstruct the phylogenetic relationships of the genus Liobagrus and estimate the divergence time of speciation events.

METHODS

We used a dissected fin clip from an adult of Liobagrus geumgangensis. Genomic DNA was extracted and analyzed with whole genome sequencing (WGS) and assembled by the NOVOPlasty method. The mitogenome sequence was annotated, and a genome map, tRNA structure, and phylogenetic tree were constructed using maximum likelihood analysis. In addition, divergence time was estimated.

RESULTS

The mitochondrial genome was 16,522 bp in length and comprised 37 genes. The overall base composition was 30.5% A, 25.5% T, 28.4% C, and 15.7% G. Most tRNAs exhibited the typical clover leaf shape, except trnS1. Phylogenetic analysis revealed that Liobagrus geumgangensis clustered within a clade with four other Liobagrus species exclusive to the southern region of the Korean Peninsula. Its divergence was estimated to have occurred during the late Miocene.

CONCLUSION

Characteristics of Liobagrus geumgangensis mitogenome were consistent with those of other torrent catfish species. Time scale estimation revealed distinct groupings, with some distributed across mainland Asia and others in the southern region of the Korean Peninsula. Notably, the Korean Peninsula group was identified as its own lineage, comprising entirely endemic species.

The complete plastome of Saxifraga giraldiana Engler and its phylogenetic analysis

Saxifraga giraldiana Engler is a common subalpine and alpine plant belonging to Saxifragaceae. However, the genetic diversity of this species has remained to be explored. In this study, we have assembled and characterized the complete chloroplast genome of S. giraldiana, filling this knowledge gap and uncovering its genetic composition. The chloroplast genome is 147,267 bp long and contains 131 genes, including 85 protein-coding genes, 38 tRNA genes, and eight rRNA genes. Furthermore, we have performed a phylogenetic analysis of 19 representative species within Saxifraga. As a result, we have found that S. giraldiana, together with S. implicans and S. stellariifolia, forms a monophyletic group. These findings have implications for the conservation and utilization of S. giraldiana.

The complete chloroplast genome sequence of red raspberry (Rubus idaeus L.) and phylogenetic analysis

Red raspberries, Rubus idaeus L. 1753 are famous fruits which possess high value bioactive compounds. In this study, we report the complete chloroplast genome of R. idaeus, it displayed a typical quadripartite structure with 155687 bp in length. The genome encodes 127 genes including 79 protein coding genes, 8 rRNA genes and 40 tRNA genes, the overall GC content is 37.2%. Phylogenetic analysis revealed a close relationship between R. idaeus and R. sachalinensis in Section Malaehobatus.

Complete chloroplast genome of endangered species Dipterocarpus retusus Blume (Dipterocarpaceae) and its phylogenetic implications

Dipterocarpus retusus Blume is an endangered species on the IUCN Red List. In this study, we reported the complete chloroplast (cp) genome of D. retusus (GenBank accession number: OP271853). The cp genome was 154,303 bp long, with a large single-copy (LSC) region of 85,586 bp and a small single-copy (SSC) region of 20,273 bp separated by a pair of inverted repeats (IRs) of 24,222 bp. It encodes 128 genes, including 84 protein-coding genes, 36 tRNA genes, and eight ribosomal RNA genes. We also reconstructed the cp genome phylogeny of Dipterocarpus, which indicated D. retusus was closely related with the sympatric species D. gracilis. This study may contribute valuable information to the phylogenetic relationships within the genus Dipterocarpus.

First genome assembly of the order Strepsiptera using PacBio HiFi reads reveals a miniature genome

Twisted-wing insects (Strepsiptera) are an enigmatic order of parasites with unusual life histories and striking sexual dimorphism. Males emerge from hosts as free-living winged adults, while females from most species remain as endoparasites that retain larval traits. Due to scarce genomic data and phylogenetic controversies, Strepsiptera was only recently placed as the closest living relative to beetles (Coleoptera). Here, we report the first PacBio HiFi genome assembly of the strepsipteran Xenos peckii (Xenidae). This de novo assembly size is 72.1 Mb, with a BUSCO score of 87.4%, N50 of 7.3 Mb, 23.4% GC content, and 38.41% repeat content. We identified 8 contigs that contain >75% of the assembly and reflect the haploid chromosome number reported from karyotypic data, and 3 contigs that exhibit sex chromosome coverage patterns. Additionally, the mitochondrial genome is 16,111 bp long and has 37 genes. This long-read assembly for Strepsiptera reveals a miniature genome and provides a unique tool to understand complex genome evolution associated with a parasitic lifestyle and extreme sexual dimorphism.

Comparative analysis of the complete mitogenomes of Camellia sinensis var. sinensis and C. sinensis var. assamica provide insights into evolution and phylogeny relationship

Introduction

Among cultivated tea plants (Camellia sinensis), only four mitogenomes for C. sinensis var. assamica (CSA) have been reported so far but none for C. sinensis var. sinensis (CSS). Here, two mitogenomes of CSS (CSSDHP and CSSRG) have been sequenced and assembled.

Methods

Using a combination of Illumina and Nanopore data for the first time. Comparison between CSS and CSA mitogenomes revealed a huge heterogeneity.

Results

The number of the repetitive sequences was proportional to the mitogenome size and the repetitive sequences dominated the intracellular gene transfer segments (accounting for 88.7%- 92.8% of the total length). Predictive RNA editing analysis revealed that there might be significant editing in NADH dehydrogenase subunit transcripts. Codon preference analysis showed a tendency to favor A/T bases and T was used more frequently at the third base of the codon. ENc plots analysis showed that the natural selection play an important role in shaping the codon usage bias, and Ka/Ks ratios analysis indicated Nad1 and Sdh3 genes may have undergone positive selection. Further, phylogenetic analysis shows that six C. sinensis clustered together, with the CSA and CSS forming two distinct branches, suggesting two different evolutionary pathway.

Discussion

Altogether, this investigation provided an insight into evolution and phylogeny relationship of C. sinensis mitogenome, thereby enhancing comprehension of the evolutionary patterns within C. sinensis species.

Complete plastome sequence of Narcissus pseudonarcissus L., one of the most iconic European plants

Narcissus pseudonarcissus L. is one of the most iconic plants of the European flora. It is a species of great horticultural interest, but also an endangered and protected plant in the wild as a consequence of loss of natural habitats. Complete plastid genome was assembled from next-generation sequencing data obtaining a circular genome of 160,008 bp long assembly. It comprises a pair of inverted repeat regions, a large single-copy region (108,400 bp), and a small single-copy region (16,434 bp). It encodes 131 genes, including 87 protein coding genes, 37 tRNA genes and seven rRNA genes. Phylogeny showed the strict relationship between N. pseudonarcissus and Narcissus poeticus L. The complete plastome will provide a useful genetic resource for future conservation programmes, phylogenetic studies and horticultural applications.

The complete chloroplast genome of Semiaquilegia danxiashanensis (Ranunculaceae), a rare species endemic to Danxia landform in Guangdong, China

Semiaquilegia danxiashanensis is currently known only from the type locality, Danxia Mountain, characterized by its spectacular red sandstone cliffscape. In this study, we assembled the complete chloroplast genome sequence of S. danxiashanensis and inferred its phylogenetic relationships. Total length of the chloroplast genome was 160,548 bp, with an overall GC content of 39%. The chloroplast genome had typical quadripartite structure and contained one LSC region (89,882 bp) and one SSC region (17,386 bp), which were separated by two IRs regions (26,640 bp, respectively). It comprised 133 genes, including 84 protein coding genes, 41 tRNA genes and eight rRNA genes. The maximum likelihood phylogenetic analysis indicated that S. danxiashanensis was sister to S. adoxoides; meanwhile, Semiaquilegia was closely related to both Urophysa and Aquilegia in Ranunculaceae. This study sheds light on the evolutionary history of Semiaquilegia and provides preliminary data for future comparative analysis of chloroplast genomes.

The chloroplast genome structure of Turpinia affinis (Staphyleaceae)

The species Turpinia affinis Merr. et Perry 1941 is widely distributed throughout southwestern China. In folk medicine, this species is often used as a substitute for the Chinese medicine Turpiniae Folium, whose legal origin is T. arguta (Lindl.) Seem. In order to ascertain the relationship between these two species, the chloroplast genome of T. affinis was aequenced and assembled, resulting in a typical quadripartite molecule with a length of 160,769 base pairs and an overall GC content of 37.3%. Additionally, 131 genes were annotated, comprising 86 protein-coding genes, 37 tRNA genes, and 8 rRNA genes. A maximum likelihood analysis demonstrated that the Turpinia species form a monophyletic clade, with T. affinis positioned as the sister taxon to the clade comprising the remaining species within the genus. This outcome enhances the genomic data for the genus Turpinia and will contribute to further investigations into phylogenetics, evolution, and sustainable resource utilization.

Two complete mitochondrial genomes of Boulenophrys (Anura: Megophryidae: Megophryinae): characteristics and phylogenetic implications

The Chinese horned toads, Boulenophrys boettgeri (Boulenger, 1899) and Boulenophrys kuatunensis (Pope, 1929), are two captivating species within the family Megophryidae, which inhabit the mountainous streams in the Eastern of China. In this study, two new complete mitochondrial genomes of B. boettgeri and B. kuatunensis were sequenced, assembled, and annotated using next-generation sequencing. The length of mitochondrial genomes of B. boettgeri and B. kuatunensis was 16,597 and 17,921 bp, respectively, with both containing 13 protein coding genes, 22 tRNA genes, two rRNA genes, and one putative control region. Phylogenetic relationships based on protein-coding mitochondrial genes showed that the two Boulenophrys species formed a cluster with other Boulenophrys species. The two new sequences provide valuable insights into the mitochondrial genomes of these two species, offering important data for understanding the phylogenetic relationships of the genus Boulenophrys.

Genome sequencing and assembly of Indian major carp, Cirrhinus mrigala (Hamilton, 1822)

The freshwater aquaculture species, Cirrhinus mrigala, commonly known as mrigal, holds significance in the carp-dominated aquaculture system, globally. Despite constituting 1.08% of the total freshwater aquaculture production, mrigal is the third most important Indian major carp. However, its genome and associated information is not available. This study aims to address this gap by generation high quality genome assembly using PacBio long reads, Illumina short reads and Hi-C scaffolding. The characterization of assembled highly contiguous genome, 1.057 Gb in size, revealed 39,091 genes with functional annotations. The orthology analysis based on direct orthologs and single copy ortholgs places C. mrigala in a distinct position within the Otophysi clade. Additionally, the study delves into Hox gene clusters, identifying 38 Hox genes distributed in seven clusters. The present genomic information offers potential applications for sustainable aquaculture management, including selection programs for economic traits.

Phylogeny, adaptive evolution, and taxonomy of Acronema (Apiaceae): evidence from plastid phylogenomics and morphological data

Introduction

The genus Acronema, belonging to Apiaceae, includes approximately 25 species distributed in the high-altitude Sino-Himalayan region from E Nepal to SW China. This genus is a taxonomically complex genus with often indistinct species boundaries and problematic generic delimitation with Sinocarum and other close genera, largely due to the varied morphological characteristics.

Methods

To explore the phylogenetic relationships and clarify the limits of the genus Acronema and its related genera, we reconstructed a reliable phylogenetic framework with high support and resolution based on two molecular datasets (plastome data and ITS sequences) and performed morphological analyses.

Results

Both phylogenetic analyses robustly supported that Acronema was a non-monophyletic group that fell into two clades: Acronema Clade and East-Asia Clade. We also newly sequenced and assembled sixteen Acronema complete plastomes and performed comprehensively comparative analyses for this genus. The comparative results showed that the plastome structure, gene number, GC content, codon bias patterns were high similarity, but varied in borders of SC/IR and we identified six different types of SC/IR border. The SC/IR boundaries of Acronema chienii were significantly different from the other Acronema members which was consistent with the type VI pattern in the genus Tongoloa. We also identified twelve potential DNA barcode regions (ccsA, matK, ndhF, ndhG, psaI, psbI, rpl32, rps15, ycf1, ycf3, psaI-ycf4 and psbM-trnD) for species identification in Acronema. The molecular evolution of Acronema was relatively conservative that only one gene (petG) was found to be under positive selection (ω = 1.02489).

Discussion

The gene petG is one of the genes involved in the transmission of photosynthetic electron chains during photosynthesis, which plays a crucial role in the process of photosynthesis in plants. This is also a manifestation of the adaptive evolution of plants in high-altitude areas to the environment. In conclusion, our study provides novel insights into the plastome adaptive evolution, phylogeny, and taxonomy of genus Acronema.

Plastome structure and phylogenetic relationships of genus Hydrocotyle (apiales): provide insights into the plastome evolution of Hydrocotyle

BACKGROUND

The genus Hydrocotyle Tourn. ex L. is a key group for further study on the evolution of Apiales, comprising around 170 species globally. Previous studies mainly focused on separate sections and provided much information about this genus, but its infrageneric relationships are still confusing. In addition, the genetic basis of its adaptive evolution remains poorly understood. To investigate the phylogeny and evolution of the genus, we selected ten representative species covering two of three diversity distribution centers and exhibiting rich morphology diversity. Comparative plastome analysis was conducted to clarify the structural character of Hydrocotyle plastomes. Positive selection analyses were implemented to assess the evolution of the genus. Phylogenetic inferences with protein-coding sequences (CDS) of Hydrocotyle and 17 related species were also performed.

RESULTS

Plastomes within Hydrocotyle were generally conservative in structure, gene order, and size. A total of 14 regions (rps16-trnK, trnQ-rps16, atpI-atpH, trnC-petN-psbM, ycf3-trnS, accD-psaI-ycf4, petA-psbJ, rps12-rpl20, rpl16 intron, rps3-rpl16 intron, rps9-rpl22, ndhF-rpl32, ndhA intron, and ycf1a) were recognized as hotspot regions within the genus, which suggested to be promising DNA barcodes for global phylogenetic analysis of Hydrocotyle. The ycf15 gene was suggested to be a protein-coding gene for Hydrocotyle species, and it could be used as a DNA barcode to identify Hydrocotyle. In phylogenetic analysis, three monophyletic clades (Clade I, II, III) were identified with evidence of rapid radiation speciation within Clade I. The selective pressure analysis detected that six CDS genes (ycf1b, matK, atpF, accD, rps14, and psbB) of Hydrocotyle species were under positive selection. Within the genus, the last four genes were conservative, suggesting a relation to the unique evolution of the genus in Apiales. Seven genes (atpE, matK, psbH, ycf1a, ycf1b, rpoA, and ycf2) were detected to be under some degree of positive selection in different taxa within the genus Hydrocotyle, indicating their role in the adaptive evolution of species.

CONCLUSIONS

Our study offers new insights into the phylogeny and adaptive evolution of Hydrocotyle. The plastome sequences could significantly enhance phylogenetic resolution and provide genomic resources and potential DNA markers useful for future studies of the genus.

The mitochondrial genome of Huananpotamon koatenense (Rathbun, 1904) (Brachyura, Potamidae) and phylogenetic analysis

Freshwater crabs play essential roles in the well-functioning of the inland aquatic ecosystems. However, due to the lack of sufficient molecular resources, the study of freshwater crabs has been greatly hindered. In this study, the mitochondrial genome of Huananpotamon koatenense, a freshwater crab endemic to China, was sequenced for the first time. This mitogenome sequence is 15,528 bp long, and contains 13 protein-coding genes, 2 rRNA genes and 22 tRNA genes. Phylogenetic analyses based on 25 mitogenomes showed that H. koatenense was clustered with the known congeneric species of H. lichuanense.

The complete mitochondrial genome of the basidiomycetous fungus, Tinctoporellus epimiltinus strain RS1

Tinctoporellus epimiltinus is widely known as a wood-decaying fungus. In the present study, we identified the complete mitochondrial genome of this species using next-generation sequencing technology. Our findings revealed that the genomic structure is a circular molecule with a size of 51,878 bp. Consistent with most Basidiomycota species, it consists of 14 core protein-coding genes, one ribosomal protein gene (rps3), 26 transfer RNA genes, and small and large ribosomal RNA (rns and rnl) genes. Seven additional open reading frames were identified. These included two sequences similar to DNA polymerases, an endonuclease-like sequence, and four hypothetical proteins. The mitochondrial genome exhibited a nucleotide composition of A (36.24%), C (12.04%), G (13.18%), and T (38.55%), resulting in a 25.21% GC content. A phylogenetic tree constructed using the combined mitochondrial gene dataset provided insight into the phylogenetic relationships of this species within the context of Basidiomycota and its members.

New taxonomic insights for Brazilian Syrbatus Reitter (Coleoptera: Staphylinidae: Pselaphinae), including three new species and their mitochondrial genomes

Here we present a taxonomic treatment for the Brazilian species of Syrbatus (Reitter, 1882), including the description of three new species (Syrbatus moustache Asenjo & Valois sp. nov., Syrbatus obsidian Asenjo & Valois sp. nov. and Syrbatus superciliata Asenjo & Valois sp. nov.) from the Quadrilátero Ferrífero (Minas Gerais, Brazil). In addition, we designated lectotypes for the Brazilian species of species-group 2, Syrbatus centralis (Raffray, 1898), Syrbatus hetschkoi (Reitter, 1888), Syrbatus hiatusus (Reitter, 1888), Syrbatus transversalis (Raffray, 1898), and Syrbatus trinodulus (Schaufuss, 1887), besides recognizing the holotype for Syrbatus brevispinus (Reitter, 1882), Syrbatus bubalus (Raffray, 1898), and Syrbatus grouvellei (Raffray, 1898). The mitochondrial genomes (mitogenomes) of the three new species are presented, for which we present the phylogenetic placement among Staphylinidae with previously published data.

Complete mitochondrial genome of Acanthochitona defilippii (Polyplacophora: Chitonida) from South Korea

The chiton (Polyplacophora) occupies a significant position in molluscan evolutionary history as one of the most primitive groups within the phylum Mollusca. Acanthochitona defilippii (Tapparone-Canefri 1874) (Chitonida: Acanthochitonidae) is a commonly found intertidal chiton species in South Korea. In this study, we characterized the complete mitochondrial genome of A. defilippii (14,999 bp long), comprising 13 protein-coding genes (PCGs), 22 transfer RNA genes, two ribosomal RNA genes, and an A + T rich region (166 bp). The base composition is as follows: 31.82% for A, 11.63% for C, 16.69% for G, and 39.86% for T. We reconstructed a maximum likelihood (ML) tree to elucidate phylogenetic relationships among the eight chitonid families using the nucleotide sequences of all PCGs. The ML tree revealed that A. defilippii clustered with Acanthochitona avicula (BP 100) within the family Acanthochitonidae. Acanthochitonidae formed a sister group with Mopaliidae. The results could provide a valuable understanding the phylogenetic relationships of chitonid species.

The complete mitochondrial genome of Heptathela kimurai (Araneae: Heptathelidae)

Heptathela kimurai (Kishida, 1920) is a spider that belongs to the family Heptathelidae which is a basial lineage of spiders. The molecular information of ancestral species belonging to families like Heptathelidae is comparatively limited when compared to spider species from derived families. Here we present the complete mitochondrial genome sequence (mtDNA) of H. kimurai. The sequence was obtained using massively parallel sequencing technology. The circular genome was 14,224 bp in length, and the AT content was 69.53%. The H. kimurai mitochondrial genome contains 13 protein-coding genes (PCGs), 21 tRNA genes, and 2 rRNA genes. The majority of PCGs were found in the heavy strand.