OMA orthology in 2024: improved prokaryote coverage, ancestral and extant GO enrichment, a revamped synteny viewer and more in the OMA Ecosystem.

In this update paper, we present the latest developments in the OMA browser knowledgebase, which aims to provide high-quality orthology inferences and facilitate the study of gene families, genomes and their evolution. First, we discuss the addition of new species in the database, particularly an expanded representation of prokaryotic species. The OMA browser now offers Ancestral Genome pages and an Ancestral Gene Order viewer, allowing users to explore the evolutionary history and gene content of ancestral genomes. We also introduce a revamped Local Synteny Viewer to compare genomic neighborhoods across both extant and ancestral genomes. Hierarchical Orthologous Groups (HOGs) are now annotated with Gene Ontology annotations, and users can easily perform extant or ancestral GO enrichments. Finally, we recap new tools in the OMA Ecosystem, including OMAmer for proteome mapping, OMArk for proteome quality assessment, OMAMO for model organism selection and Read2Tree for phylogenetic species tree construction from reads. These new features provide exciting opportunities for orthology analysis and comparative genomics. OMA is accessible at https://omabrowser.org.

The fitness consequences of genetic divergence between polymorphic gene arrangements.

Inversions restrict recombination when heterozygous with standard arrangements, but often have few noticeable phenotypic effects. Nevertheless, there are several examples of inversions that can be maintained polymorphic by strong selection under laboratory conditions. A long-standing model for the source of such selection is divergence between arrangements with respect to recessive or partially recessive deleterious mutations, resulting in a selective advantage to heterokaryotypic individuals over homokaryotypes. This paper uses a combination of analytical and numerical methods to investigate this model, for the simple case of an autosomal inversion with multiple independent nucleotide sites subject to mildly deleterious mutations. A complete lack of recombination in heterokaryotypes is assumed, as well as constancy of the frequency of the inversion over space and time. It is shown that a significantly higher mutational load will develop for the less frequent arrangement. A selective advantage to heterokaryotypes is only expected when the two alternative arrangements are nearly equal in frequency, so that their mutational loads are very similar in size. The effects of some Drosophila pseudoobscura polymorphic inversions on fitness traits seem to be too large to be explained by this process, although it may contribute to some of the observed effects. Several population genomic statistics can provide evidence for signatures of a reduced efficacy of selection associated with the rarer of two arrangements, but there is currently little published data that are relevant to the theoretical predictions.

Phylomitogenomics elucidates the evolution of symbiosis in Thoracotremata (Decapoda: Cryptochiridae, Pinnotheridae, Varunidae).

Background: Thoracotremata belong to the large group of "true" crabs (infraorder Brachyura), and they exhibit a wide range of physiological and morphological adaptations to living in terrestrial, freshwater and marine habitats. Moreover, the clade comprises various symbiotic taxa (Aphanodactylidae, Cryptochiridae, Pinnotheridae, some Varunidae) that are specialised in living with invertebrate hosts, but the evolutionary history of these symbiotic crabs is still partially unresolved. Methods: Here we assembled and characterised the complete mitochondrial genomes (hereafter mitogenomes) of three gall crab species (Cryptochiridae): Kroppcarcinus siderastreicola, Opecarcinus hypostegus and Troglocarcinus corallicola. A phylogenetic tree of the Thoracotremata was reconstructed using 13 protein-coding genes and two ribosomal RNA genes retrieved from three new gall crab mitogenomes and a further 72 available thoracotreme mitogenomes. Furthermore, we applied a comparative analysis to characterise mitochondrial gene order arrangement, and performed a selection analysis to test for selective pressure of the protein-coding genes in symbiotic Cryptochiridae, Pinnotheridae, and Varunidae (Asthenognathus inaequipes and Tritodynamia horvathi). Results: The results of the phylogenetic reconstruction confirm the monophyly of Cryptochiridae, which clustered separately from the Pinnotheridae. The latter clustered at the base of the tree with robust branch values. The symbiotic varunids A. inaequipes and T. horvathi clustered together in a clade with free-living Varunidae species, highlighting that symbiosis in the Thoracotremata evolved independently on multiple occasions. Different gene orders were detected in symbionts and free-living species when compared with the ancestral brachyuran gene order. Lastly, the selective pressure analysis detected two positively selected sites in the nad6 gene of Cryptochiridae, but the evidence for positive selection in Pinnotheridae and A. inaequipes and T. horvathi was weak. Adaptive evolution of mitochondrial protein-coding genes is perhaps related to the presumably higher energetic demands of a symbiotic lifestyle.

Mitochondrial genome characteristics of six Phylloscopus species and their phylogenetic implication.

The mitochondrial genomes of six Phylloscopus species-small insectivores belonging to the Phylloscopidae family-were obtained using the Illumina sequencing platform. The mitogenomes were closed circular molecules 16,922-17,007 bp in size, containing 13 protein-coding genes, 22 tRNA genes, two rRNA genes, and two control regions (CR1 and remnant CR2). The gene orders were conserved in 35 sampled Phylloscopus mitogenomes in the GenBank database, with a gene rearrangement of cytb-trnT-CR1-trnP-nad6-trnE-remnant CR2-trnF-rrnS. The average base compositions of the six Phylloscopus mitogenomes were 29.43% A, 32.75% C, 14.68% G, and 23.10% T, with the A+T content slightly higher than that of G+C. ATG and TAA were the most frequent initiating and terminating codons, respectively. Several conserved boxes were identified in CR1, including C-string in domain I; F, E, D, and C boxes, as well as bird similarity and B boxes, in domain II; and CSB1 in domain III. Tandem repeats were observed in remnant CR2 of the Phylloscopus fuscatus and Phylloscopus proregulus mitogenomes. A phylogenetic analysis with maximum likelihood (ML) and Bayesian inference (BI) methods, based on 13 protein-coding genes and two rRNA genes, indicated that the Phylloscopus species was divided into two larger clades, with a splitting time approximately 11.06 million years ago (mya). The taxa of Phylloscopus coronatus/Phylloscopus burkii and Phylloscopus inornatus/P. proregulus were located at the basal position of the different clades. The phylogenetic result of the cox1 gene showed that Seicercus was nested within Phylloscopus. The complete set of mitogenomes of the Phylloscopus species provides potentially useful resources for the further exploration of the taxonomic status and phylogenetic history of Phylloscopidae.

Relaxed purifying selection pressure drives accelerated and dynamic gene rearrangements in thrips (Insecta: Thysanoptera) mitochondrial genomes.

Insect mitochondrial genomes (mitogenome) generally present a typical gene order, which is considered as the ancestral arrangement. All sequenced mitogenomes in the Thysanoptera display high levels of gene rearrangement. Due to limited number of thrips mitogenomes sequenced, how gene rearrangement may be shaped by evolution remain unclear. Here, we analyzed 33 thrips mitogenomes, including 14 newly sequenced. These mitogenomes were diverse in organization, nucleotides substitution and gene arrangements. We found 28 highly rearranged gene orders with the breakpoints of gene rearrangements from 25 to 33. Reconstruction of the ancestors mitochondrial gene arrangements states indicated that Tubulifera have more complex pathways than Terebrantia in the gene order evolution. Molecular calibration estimated that divergence of two suborders occurred in the middle Triassic while the radiation of thrips was associated with the arose and flourish of angiosperm. Our evolutionary hypothesis testing suggests that relaxation of selection pressure enabled the early phase of Thysanoptera evolution, followed by a stronger selective pressure fixed diversification. Our analyses found gene inversion increases the nonsynonymous substitution rates and provide an evolutionary hypothesis driving the diverse gene orders.

Comparative Mitogenomic Analyses of Darkling Beetles (Coleoptera: Tenebrionidae) Provide Evolutionary Insights into tRNA-like Sequences.

Tenebrionidae is widely recognized owing to its species diversity and economic importance. Here, we determined the mitochondrial genomes (mitogenomes) of three Tenebrionidae species (Melanesthes exilidentata, Anatolica potanini, and Myladina unguiculina) and performed a comparative mitogenomic analysis to characterize the evolutionary characteristics of the family. The tenebrionid mitogenomes were highly conserved with respect to genome size, gene arrangement, base composition, and codon usage. All protein-coding genes evolved under purifying selection. The largest non-coding region (i.e., control region) showed several unusual features, including several conserved repetitive fragments (e.g., A+T-rich regions, G+C-rich regions, Poly-T tracts, TATA repeat units, and longer repetitive fragments) and tRNA-like structures. These tRNA-like structures can bind to the appropriate anticodon to form a cloverleaf structure, although base-pairing is not complete. We summarized the quantity, types, and conservation of tRNA-like sequences and performed functional and evolutionary analyses of tRNA-like sequences with various anticodons. Phylogenetic analyses based on three mitogenomic datasets and two tree inference methods largely supported the monophyly of each of the three subfamilies (Stenochiinae, Pimeliinae, and Lagriinae), whereas both Tenebrioninae and Diaperinae were consistently recovered as polyphyletic. We obtained a tenebrionid mitogenomic phylogeny: (Lagriinae, (Pimeliinae, ((Tenebrioninae + Diaperinae), Stenochiinae))). Our results provide insights into the evolution and function of tRNA-like sequences in tenebrionid mitogenomes and contribute to our general understanding of the evolution of Tenebrionidae.

The effects of inversion polymorphisms on patterns of neutral genetic diversity.

The strong reduction in the frequency of recombination in heterozygotes for an inversion and a standard gene arrangement causes the arrangements to become partially isolated genetically, resulting in sequence divergence between them and changes in the levels of neutral variability at nucleotide sites within each arrangement class. Previous theoretical studies on the effects of inversions on neutral variability have assumed either that the population is panmictic or that it is divided into 2 populations subject to divergent selection. Here, the theory is extended to a model of an arbitrary number of demes connected by migration, using a finite island model with the inversion present at the same frequency in all demes. Recursion relations for mean pairwise coalescent times are used to obtain simple approximate expressions for diversity and divergence statistics for an inversion polymorphism at equilibrium under recombination and drift, and for the approach to equilibrium following the sweep of an inversion to a stable intermediate frequency. The effects of an inversion polymorphism on patterns of linkage disequilibrium are also examined. The reduction in effective recombination rate caused by population subdivision can have significant effects on these statistics. The theoretical results are discussed in relation to population genomic data on inversion polymorphisms, with an emphasis on Drosophila melanogaster. Methods are proposed for testing whether or not inversions are close to recombination-drift equilibrium, and for estimating the rate of recombinational exchange in heterozygotes for inversions; difficulties involved in estimating the ages of inversions are also discussed.

Mitochondrial Genome Evolution in Annelida-A Systematic Study on Conservative and Variable Gene Orders and the Factors Influencing its Evolution.

The mitochondrial genomes of Bilateria are relatively conserved in their protein-coding, rRNA, and tRNA gene complement, but the order of these genes can range from very conserved to very variable depending on the taxon. The supposedly conserved gene order of Annelida has been used to support the placement of some taxa within Annelida. Recently, authors have cast doubts on the conserved nature of the annelid gene order. Various factors may influence gene order variability including, among others, increased substitution rates, base composition differences, structure of noncoding regions, parasitism, living in extreme habitats, short generation times, and biomineralization. However, these analyses were neither done systematically nor based on well-established reference trees. Several focused on only a few of these factors and biological factors were usually explored ad-hoc without rigorous testing or correlation analyses. Herein, we investigated the variability and evolution of the annelid gene order and the factors that potentially influenced its evolution, using a comprehensive and systematic approach. The analyses were based on 170 genomes, including 33 previously unrepresented species. Our analyses included 706 different molecular properties, 20 life-history and ecological traits, and a reference tree corresponding to recent improvements concerning the annelid tree. The results showed that the gene order with and without tRNAs is generally conserved. However, individual taxa exhibit higher degrees of variability. None of the analyzed life-history and ecological traits explained the observed variability across mitochondrial gene orders. In contrast, the combination and interaction of the best-predicting factors for substitution rate and base composition explained up to 30% of the observed variability. Accordingly, correlation analyses of different molecular properties of the mitochondrial genomes showed an intricate network of direct and indirect correlations between the different molecular factors. Hence, gene order evolution seems to be driven by molecular evolutionary aspects rather than by life history or ecology. On the other hand, variability of the gene order does not predict if a taxon is difficult to place in molecular phylogenetic reconstructions using sequence data or not. We also discuss the molecular properties of annelid mitochondrial genomes considering canonical views on gene evolution and potential reasons why the canonical views do not always fit to the observed patterns without making some adjustments. [Annelida; compositional biases; ecology; gene order; life history; macroevolution; mitochondrial genomes; substitution rates.].

From comparative gene content and gene order to ancestral contigs, chromosomes and karyotypes.

To reconstruct the ancestral genome of a set of phylogenetically related descendant species, we use the RACCROCHE pipeline for organizing a large number of generalized gene adjacencies into contigs and then into chromosomes. Separate reconstructions are carried out for each ancestral node of the phylogenetic tree for focal taxa. The ancestral reconstructions are monoploids; they each contain at most one member of each gene family constructed from descendants, ordered along the chromosomes. We design and implement a new computational technique for solving the problem of estimating the ancestral monoploid number of chromosomes x. This involves a "g-mer" analysis to resolve a bias due long contigs, and gap statistics to estimate x. We find that the monoploid number of all the rosid and asterid orders is [Formula: see text]. We show that this is not an artifact of our method by deriving [Formula: see text] for the metazoan ancestor.

Evolution of mitogenomic gene order in Orthoptera.

Mitochondrial gene order has contributed to the elucidation of evolutionary relationships in several animal groups. It generally has found its application as a phylogenetic marker for deep nodes. Yet, in Orthoptera limited research has been performed on the gene order, although the group represents one of the oldest insect orders. We performed a comprehensive study on mitochondrial genome rearrangements (MTRs) within Orthoptera in the context of mitogenomic sequence-based phylogeny. We used 280 published mitogenome sequences from 256 species, including three outgroup species, to reconstruct a molecular phylogeny. Using a heuristic approach, we assigned MTR scenarios to the edges of the phylogenetic tree and reconstructed ancestral gene orders to identify possible synapomorphies in Orthoptera. We found all types of MTRs in our dataset: inversions, transpositions, inverse transpositions, and tandem-duplication/random loss events (TDRL). Most of the suggested MTRs were in single and unrelated species. Out of five MTRs which were unique in subgroups of Orthoptera, we suggest four of them to be synapomorphies; those were in the infraorder Acrididea, in the tribe Holochlorini, in the subfamily Pseudophyllinae, and in the two families Phalangopsidae and Gryllidae or their common ancestor (leading to the relationship ((Phalangopsidae + Gryllidae) + Trigonidiidae)). However, similar MTRs have been found in distant insect lineages. Our findings suggest convergent evolution of specific mitochondrial gene orders in several species, deviant from the evolution of the mitogenome DNA sequence. As most MTRs were detected at terminal nodes, a phylogenetic inference of deeper nodes based on MTRs is not supported. Hence, the marker does not seem to aid resolving the phylogeny of Orthoptera, but adds further evidence for the complex evolution of the whole group, especially at the genetic and genomic levels. The results indicate a high demand for more research on patterns and underlying mechanisms of MTR events in Orthoptera.

The Complete Mitochondrial Genome and Gene Arrangement of the Enigmatic Scaphopod Pictodentalium vernedei.

The enigmatic scaphopods, or tusk shells, are a small and rare group of molluscs whose phylogenomic position among the Conchifera is undetermined, and the taxonomy within this class also needs revision. Such work is hindered by there only being a very few mitochondrial genomes in this group that are currently available. Here, we present the assembly and annotation of the complete mitochondrial genome from Dentaliida Pictodentalium vernedei, whose mitochondrial genome is 14,519 bp in size, containing 13 protein-coding genes, 22 tRNA genes and two rRNA genes. The nucleotide composition was skewed toward A-T, with a 71.91% proportion of AT content. Due to the mitogenome-based phylogenetic analysis, we defined P. vernedei as a sister to Graptacme eborea in Dentaliida. Although a few re-arrangements occurred, the mitochondrial gene order showed deep conservation within Dentaliida. Yet, such a gene order in Dentaliida largely diverges from Gadilida and other molluscan classes, suggesting that scaphopods have the highest degree of mitogenome arrangement compared to other molluscs.

The metazoan landscape of mitochondrial DNA gene order and content is shaped by selection and affects mitochondrial transcription.

Mitochondrial DNA (mtDNA) harbors essential genes in most metazoans, yet the regulatory impact of the multiple evolutionary mtDNA rearrangements has been overlooked. Here, by analyzing mtDNAs from ~8000 metazoans we found high gene content conservation (especially of protein and rRNA genes), and codon preferences for mtDNA-encoded tRNAs across most metazoans. In contrast, mtDNA gene order (MGO) was selectively constrained within but not between phyla, yet certain gene stretches (ATP8-ATP6, ND4-ND4L) were highly conserved across metazoans. Since certain metazoans with different MGOs diverge in mtDNA transcription, we hypothesized that evolutionary mtDNA rearrangements affected mtDNA transcriptional patterns. As a first step to test this hypothesis, we analyzed available RNA-seq data from 53 metazoans. Since polycistron mtDNA transcripts constitute a small fraction of the steady-state RNA, we enriched for polycistronic boundaries by calculating RNA-seq read densities across junctions between gene couples encoded either by the same strand (SSJ) or by different strands (DSJ). We found that organisms whose mtDNA is organized in alternating reverse-strand/forward-strand gene blocks (mostly arthropods), displayed significantly reduced DSJ read counts, in contrast to organisms whose mtDNA genes are preferentially encoded by one strand (all chordates). Our findings suggest that mtDNA rearrangements are selectively constrained and likely impact mtDNA regulation.

Insights into chloroplast genome structure, intraspecific variation, and phylogeny of Cyclamen species (Myrsinoideae).

Species from the flowering plant genus Cyclamen are popular amongst consumers. In particular Cyclamen persicum Mill. has been significantly used commercially, and certain small flowering species such as Cyclamen hederifolium and Cyclamen coum are gradually growing in popularity in the potted flower market. Here, the chloroplast genomes of nine Cyclamen samples including four Cyclamen species and five varieties of C. hederifolium were sequenced for genome structure comparison, White green septal striped leaves related gene screening and DNA molecular markers were developed for phylogenetic analysis. In comparing Cyclamen species' chloroplast genomes, gene content and gene order were found to be highly similar with the length of genomes ranging from 151,626 to 153,058 bp. The chloroplast genome of Cyclamen has 128 genes, including 84 protein-coding genes, 36 transfer RNA genes, and 8 ribosomal RNA genes. Based on intraspecific variation, seven hotspots, including three genes and four intergenic regions, were identified as variable markers for downstream species delimitation and interspecific relationship analyses. Moreover, a phylogenetic tree constructed with complete chloroplast genomes, revealed that Cyclamen are monophyletic with Lysimachia as the closest neighbor. Phylogenetic analyses of the 14 Cyclamen species with the seven variable regions showed five distinct clades within this genus. The highly supported topologies showed these seven regions may be used as candidate DNA barcode sequences to distinguish Cyclamen species. White green septal striped leaves is common in C. hederifolium, however the molecular mechanism of this has not yet been described. Here, we find that the intergenic region rps4-trnT-UGU seems related to white green septal striped leaves.

The complete mitochondrial genome and novel gene arrangement in Nesodiprion zhejiangensis Zhou & Xiao (Hymenoptera: Diprionidae).

The complete mitochondrial genome (mitogenome) of the sawfly, Nesodiprion zhejiangensis Zhou & Xiao, was sequenced, assembled, and deposited in GenBank (Accession Number: OM501121). The 15,660 bp N. zhejiangensis mitogenome encodes for 2 ribosomal RNAs (rrnL and rrnS), 22 transfer RNAs (tRNAs), 13 protein-coding genes (PCGs), and an AT-rich region of 450 bp in length. The nucleotide composition is biased toward adenine and thymine (A + T = 81.8%). Each PCG is initiated by an ATN codon, except for cox2, which starts with a TTG. Of 13 PCGs, 9 have a TAA termination codon, while the remainder terminate with a TAG or a single T. All tRNAs have the classic cloverleaf structure, except for the dihydrouridine (DHU) arm of tRNAval, which forms a simple loop. There are 49 helices belonging to 6 domains in rrnL and 30 helices belonging to 4 domains in rrnS. In comparison to the ancestral architecture, N. zhejiangensis has the most rearranged mitogenome in Symphyta, in which rearrangement events of local inversion and transposition are identified in three gene clusters. Specifically, the main hotspot of gene rearrangement occurred between rrnS and trnY, and rearranged from rrnS-(AT-rich region)-I-Q-M-nd2-W-C-Y to rrnS-Q-W-C-nd2-I-M-(AT-rich region)-Y, involving a local inversion event of a large gene cluster and transposition events of some tRNAs. Transposition of trnA and trnR (rearranged from A-R to R-A) was observed at the nd3-nd5 gene junction while shuffling of trnP and trnT (rearranged from T-P to P-T) occurred at the nd4l-nd6 gene junction. While illegitimate inter-mtDNA recombination might explain the opposite orientations of transcription between rrnS and trnY, transposition events of tRNA in some gene blocks can be accounted for by the tandem duplication/random loss (TDRL) model. Our phylogenetic analysis suggests that N. zhejiangensis is closely related to congeneric species N. biremis and N. japonicus, which together form a sister lineage with the European pine sawfly, Neodiprion sertifer.

Evolution of mitochondrial and nuclear genomes in Pennatulacea.

We examine the phylogeny of sea pens using sequences of whole mitochondrial genomes and the nuclear ribosomal cluster generated through low coverage Illumina sequencing. Taxon sampling includes 30 species in 19 genera representing 13 families. Ancestral state reconstruction shows that most sea pen mitochondrial genomes have the ancestral gene order, and that Pennatulacea with diverse gene orders are found in a single clade. The monophyly of Pennatulidae and Protoptilidae are rejected by both the mitochondrial and nuclear dataset, while the mitochondrial dataset further rejects monophyly of Virgulariidae, and the nuclear dataset rejects monophyly of Kophobelemnidae. We show discordance between nuclear ribosomal gene cluster phylogenies and whole mitochondrial genome phylogenies and highlight key Pennatulacea taxa that could be included in cnidarian genome-wide studies to better resolve the sea pen tree of life. We further illustrate how well frequently sequenced markers capture the overall diversity of the mitochondrial genome and the nuclear ribosomal genes in sea pens.

Taxonomic status and phylogenetic relationship of Anomura (Crustacea: Decapoda) based on mitochondrial sequences and gene order rearrangements.

The complete mitochondrial genome (mitogenome) contributed crucial information, which could improve the better understanding of molecular phylogenetic analysis, evolution, and gene rearrangements. However, only a few Coenobitidae and Albuneidae mitogenomes have been described, and controversies about the Anomuran phylogeny remain. Here, we determined three Coenobitidae species (Coenobita perlatus, Coenobita rugosus, Coenobita clypeatus) and one Albuneidae species (Emerita talpoida). As the representative species of genera or families, and even the first species of the genus, supplementing these four species is essential to study the phylogenetic relationships within Anomura. Most of the four novel mitogenomes have the typical 37 genes, except for E. talpoida, which lacked trnF. Through combining these new data with 31 Anomuran mitogenomes from Genbank, different gene rearrangement patterns and internal phylogenetic relationships of Anomura were investigated. In our phylogeny analysis, all the Anomuran species were clustered into one group, and the polyphyly of Paguroidea was well supported. Moreover, various peculiar mitochondrial gene orders (MGOs) were summarized among Anomura and preliminary determined their rearrangement mechanisms through CREx. Twenty different MGOs were suggested in our research, and we were focused seven MGO patterns (Pattern A-J) in the process of discussing the gene rearrangement mechanism. Currently making full use of the large taxon sampling, our results provide valuable information on the evolutionary status of Anomuran species and are available for systematic rearrangement and phylogenetic analyses.

Mitogenome of Xya pfaendleri (Orthoptera: Caelifera): Its comparative description and phylogenetic position within Tridactylidea.

We report the comparative examination of the complete mitochondrial genome of the pygmy mole cricket Xya pfaendleri (Orthoptera: Caelifera: Tridactylidae). The mitogenome consists of 13 protein-coding regions, 22 tRNAs, two rRNAs, and one control region, following the gene order of the ancestral pancrustacean mitogenome. The length of the mitogenome in Xya pfaendleri is 15352 bp. The start and stop codons of the protein-coding genes exhibit the general pattern observed in orthopterans. The data indicate that the pattern of gene overlapping/intergenic sequences exhibits a significant phylogenetic signal. A phylogenetic tree inferred using 12 mitogenomes (seven belonging to Tridactylidea, three to Acrididea, and two to Ensifera) confirms the sister group relationship of Acrididea and Tridactylidea. The relationship among the families of Tridactylidea is Cylindrachetidae + (Ripipterygidae + Tridactylidae). The mitogenome sequences of Xya and Tridactylus constitute a single clade, sharing a last common ancestor 94 million years ago, and rendering the first genus paraphyletic. The present preliminary data suggest that we still have much to learn about the evolution and diversity of Tridactylidea.

OrchidBase 5.0: updates of the orchid genome knowledgebase.

Containing the largest number of species, the orchid family provides not only materials for studying plant evolution and environmental adaptation, but economically and culturally important ornamental plants for human society. Previously, we collected genome and transcriptome information of Dendrobium catenatum, Phalaenopsis equestris, and Apostasia shenzhenica which belong to two different subfamilies of Orchidaceae, and developed user-friendly tools to explore the orchid genetic sequences in the OrchidBase 4.0. The OrchidBase 4.0 offers the opportunity for plant science community to compare orchid genomes and transcriptomes and retrieve orchid sequences for further study.In the year 2022, two whole-genome sequences of Orchidoideae species, Platanthera zijinensis and Platanthera guangdongensis, were de novo sequenced, assembled and analyzed. In addition, systemic transcriptomes from these two species were also established. Therefore, we included these datasets to develop the new version of OrchidBase 5.0. In addition, three new functions including synteny, gene order, and miRNA information were also developed for orchid genome comparisons and miRNA characterization.OrchidBase 5.0 extended the genetic information to three orchid subfamilies (including five orchid species) and provided new tools for orchid researchers to analyze orchid genomes and transcriptomes. The online resources can be accessed at https://cosbi.ee.ncku.edu.tw/orchidbase5/.

Characterization of the Evolutionary Pressure on Anisodus tanguticus Maxim. with Complete Chloroplast Genome Sequence.

Anisodus tanguticus Maxim. (Solanaceae), a traditional endangered Tibetan herb, is endemic to the Qinghai-Tibet Plateau. Here, we report the de novo assembled chloroplast (cp) genome sequences of A. tanguticus (155,765 bp). The cp contains a pair of inverted repeated (IRa and IRb) regions of 25,881 bp that are separated by a large single copy (LSC) region (86,516 bp) and a small single copy SSC (17,487 bp) region. A total of 132 functional genes were annotated in the cp genome, including 87 protein-coding genes, 37 tRNA genes, and 8 rRNA genes. Moreover, 199 simple sequence repeats (SSR) and 65 repeat structures were detected. Comparative plastome analyses revealed a conserved gene order and high similarity of protein-coding sequences. The A. tanguticus cp genome exhibits contraction and expansion, which differs from Przewalskia tangutica and other related Solanaceae species. We identified 30 highly polymorphic regions, mostly belonging to intergenic spacer regions (IGS), which may be suitable for the development of robust and cost-effective markers for inferring the phylogeny of the genus Anisodus and family Solanaceae. Analysis of the Ka/Ks ratios of the Hyoscyameae tribe revealed significant positive selection exerted on the cemA, rpoC2, and clpP genes, which suggests that protein metabolism may be an important strategy for A. tanguticus and other species in Hyoscyameae in adapting to the adverse environment on the Qinghai-Tibetan Plateau. Phylogenetic analysis revealed that A. tanguticus clustered closer with Hyoscyamus niger than P. tangutica. Our results provide reliable genetic information for future exploration of the taxonomy and phylogenetic evolution of the Hyoscyameae tribe and related species.

Complete Mitogenome of Oreolalax Omeimontis Reveals Phylogenetic Status and Novel Gene Arrangement of Archaeobatrachia.

Species of the genus Oreolalax displayed crucial morphological characteristics of vertebrates transitioning from aquatic to terrestrial habitats; thus, they can be regarded as a representative vertebrate genus for this landing phenomenon. But the present phylogenetic status of Oreolalax omeimontis has been controversial with morphological and molecular approaches, and specific gene rearrangements were discovered in all six published Oreolalax mitogenomes, which are rarely observed in Archaeobatrachia. Therefore, this study determined the complete mitogenome of O. omeimontis with the aim of identifying its precise phylogenetic position and novel gene arrangement in Archaeobatrachia. Phylogenetic analysis with Bayesian inference and maximum likelihood indicates O. omeimontis is a sister group to O. lichuanensis, which is consistent with previous phylogenetic analysis based on morphological characteristics, but contrasts with other studies using multiple gene fragments. Moreover, although the duplication of trnM occurred in all seven Oreolalax species, the translocation of trnQ and trnM occurred differently in O. omeimontis to the other six, and this unique rearrangement would happen after the speciation of O. omeimontis. In general, this study sheds new light on the phylogenetic relationships and gene rearrangements of Archaeobatrachia.