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

BACKGROUND: Lycium is an economically and ecologically important genus of shrubs, consisting of approximately 70 species distributed worldwide, 15 of which are located in China. Despite the economic and ecological importance of Lycium, its phylogeny, interspecific relationships, and evolutionary history remain relatively unknown. In this study, we constructed a phylogeny and estimated divergence time based on the chloroplast genomes (CPGs) of 15 species, including subspecies, of the genus Lycium from China. RESULTS: We sequenced and annotated 15 CPGs in this study. Comparative analysis of these genomes from these Lycium species revealed a typical quadripartite structure, with a total sequence length ranging from 154,890 to 155,677 base pairs (bp). The CPGs was highly conserved and moderately differentiated. Through annotation, we identified a total of 128-132 genes. Analysis of the boundaries of inverted repeat (IR) regions showed consistent positioning: the junctions of the IRb/LSC region were located in rps19 in all Lycium species, IRb/SSC between the ycf1 and ndhF genes, and SSC/IRa within the ycf1 gene. Sequence variation in the SSC region exceeded that in the IR region. We did not detect major expansions or contractions in the IR region or rearrangements or insertions in the CPGs of the 15 Lycium species. Comparative analyses revealed five hotspot regions in the CPG: trnR(UCU), atpF-atpH, ycf3-trnS(GGA), trnS(GGA), and trnL-UAG, which could potentially serve as molecular markers. In addition, phylogenetic tree construction based on the CPG indicated that the 15 Lycium species formed a monophyletic group and were divided into two typical subbranches and three minor branches. Molecular dating suggested that Lycium diverged from its sister genus approximately 17.7 million years ago (Mya) and species diversification within the Lycium species of China primarily occurred during the recent Pliocene epoch. CONCLUSION: The divergence time estimation presented in this study will facilitate future research on Lycium, aid in species differentiation, and facilitate diverse investigations into this economically and ecologically important genus.

Complete chloroplast genomes of 13 species of sect. Tuberculata Chang (Camellia L.): genomic features, comparative analysis, and phylogenetic relationships.

BACKGROUND: Sect. Tuberculata belongs to Camellia, and its members are characterized by a wrinkled pericarp and united filaments. All the plants in this group, which are endemic to China, are highly valuable for exploring the evolution of Camellia and have great potential for use as an oil source. However, due to the complex and diverse phenotypes of these species and the difficulty of investigating them in the field, their complex evolutionary history and interspecific definitions have remained largely unelucidated. RESULTS: Therefore, we newly sequenced and annotated 12 chloroplast (cp) genomes and retrieved the published cp genome of Camellia anlungensis Chang in sect. Tuberculata. In this study, comparative analysis of the cp genomes of the thirteen sect. Tuberculata species revealed a typical quadripartite structure characterized by a total sequence length ranging from 156,587 bp to 157,068 bp. The cp.genome arrangement is highly conserved and moderately differentiated. A total of 130 to 136 genes specific to the three types were identified by annotation, including protein-coding genes (coding sequences (CDSs)) (87-91), tRNA genes (35-37), and rRNA genes (8). The total observed frequency ranged from 23,045 (C. lipingensis) to 26,557 (C. anlungensis). IR region boundaries were analyzed to show that the ycf1 gene of C. anlungensis is located in the IRb region, while the remaining species are present only in the IRa region. Sequence variation in the SSC region is greater than that in the IR region, and most protein-coding genes have high codon preferences. Comparative analyses revealed six hotspot regions (tRNA-Thr(GGT)-psbD, psbE-petL, ycf15-tRNA-Leu(CAA), ndhF-rpl32, ndhD, and trnL(CAA)-ycf15) in the cp genomes that could serve as potential molecular markers. In addition, the results of phylogenetic tree construction based on the cp genomes showed that the thirteen sect. Tuberculata species formed a monophyletic group and were divided into two evolutionarily independent clades, confirming the independence of the section. CONCLUSIONS: In summary, we obtained the cp genomes of thirteen sect. Tuberculata plants and performed the first comparative analysis of this group. These results will help us better characterize the plants in this section, deepen our understanding of their genetic characteristics and phylogenetic relationships, and lay the theoretical foundation for their accurate classification, elucidation of their evolutionary changes, and rational development and utilization of this section in the future.

Genome-Wide Identification and Expression Analysis of the COL Gene Family in Hemerocallis citrina Baroni.

Hemerocallis citrina Baroni (H. citrina) is an important specialty vegetable that is not only edible and medicinal but also has ornamental value. However, much remains unknown about the regulatory mechanisms associated with the growth, development, and flowering rhythm of this plant. CO, as a core regulatory factor in the photoperiod pathway, coordinates light and circadian clock inputs to transmit flowering signals. We identified 18 COL genes (HcCOL1-HcCOL18) in the H. citrina cultivar 'Mengzihua' and studied their chromosomal distribution, phylogenetic relationships, gene and protein structures, collinearity, and expression levels in the floral organs at four developmental stages. The results indicate that these genes can be classified into three groups based on phylogenetic analysis. The major expansion of the HcCOL gene family occurred via segmental duplication, and the Ka/Ks ratio indicated that the COL genes of Arabidopsis thaliana, Oryza sativa, Phalaenopsis equestris, and H. citrina were under purifying selection. Many cis-elements, including light response elements, abiotic stress elements, and plant hormone-inducible elements, were distributed in the promoter sequences of the HcCOL genes. Expression analysis of HcCOL genes at four floral developmental stages revealed that most of the HcCOL genes were expressed in floral organs and might be involved in the growth, development, and senescence of the floral organs of H. citrina. This study lays a foundation for the further elucidation of the function of the HcCOL gene in H. citrina and provides a theoretical basis for the molecular design breeding of H. citrina.

Comprehensive comparative analysis and development of molecular markers for Lasianthus species based on complete chloroplast genome sequences.

BACKGROUND: Lasianthus species are widely used in traditional Chinese folk medicine with high medicinal value. However, source materials and herbarium specimens are often misidentified due to morphological characteristics and commonly used DNA barcode fragments are not sufficient for accurately identifying Lasianthus species. To improve the molecular methods for distinguishing among Lasianthus species, we report the complete chloroplast (CP) genomes of Lasianthus attenuatus, Lasianthus henryi, Lasianthus hookeri, Lasianthus sikkimensis, obtained via high-throughput Illumina sequencing. RESULTS: These showed CP genomes size of 160164-160246 bp and a typical quadripartite structure, including a large single-copy region (86675-86848 bp), a small single-copy region (17177-17326 bp), and a pair of inverted repeats (28089-28135 bp). As a whole, the gene order, GC content and IR/SC boundary structure were remarkably similar among of the four Lasianthus CP genomes, the partial gene length and IR, LSC and SSC regions length are still different. The average GC content of the CP genomes was 36.71-36.75%, and a total of 129 genes were detected, including 83 different protein-coding genes, 8 different rRNA genes and 38 different tRNA genes. Furthermore, we compared our 4 complete CP genomes data with publicly available CP genome data from six other Lasianthus species, and we initially screened eleven highly variable region fragments were initially screened. We then evaluated the identification efficiency of eleven highly variable region fragments and 5 regular barcode fragments. Ultimately, we found that the optimal combination fragment' ITS2 + psaI-ycf4' could authenticated the Lasianthus species well. Additionally, the results of genome comparison of Rubiaceae species showed that the coding region is more conservative than the non-coding region, and the ycf1 gene shows the most significant variation. Finally, 49 species of CP genome sequences belonging to 16 genera of the Rubiaceae family were used to construct phylogenetic trees. CONCLUSIONS: Our research is the first to analyze the chloroplast genomes of four species of Lasianthus in detail and we ultimately determined that the combination fragment' ITS2 + psaI-ycf4' is the optimal barcode combination for identifying the genus of Lasianthus. Meanwhile, we gathered the available CP genome sequences from the Rubiaceae and used them to construct the most comprehensive phylogenetic tree for the Rubiaceae family. These investigations provide an important reference point for further studies in the species identification, genetic diversity, and phylogenetic analyses of Rubiaceae species.

De novo assembly of the complete mitochondrial genome of pepino (Solanum muricatum) using PacBio HiFi sequencing: insights into structure, phylogenetic implications, and RNA editing.

BACKGROUND: Solanum muricatum is an emerging horticultural fruit crop with rich nutritional and antioxidant properties. Although the chromosome-scale genome of this species has been sequenced, its mitochondrial genome sequence has not been reported to date. RESULTS: PacBio HiFi sequencing was used to assemble the circular mitogenome of S. muricatum, which was 433,466 bp in length. In total, 38 protein-coding, 19 tRNA, and 3 rRNA genes were annotated. The reticulate mitochondrial conformations with multiple junctions were verified by polymerase chain reaction, and codon usage, sequence repeats, and gene migration from chloroplast to mitochondrial genome were determined. A collinearity analysis of eight Solanum mitogenomes revealed high structural variability. Overall, 585 RNA editing sites in protein coding genes were identified based on RNA-seq data. Among them, mttB was the most frequently edited (52 times), followed by ccmB (46 times). A phylogenetic analysis based on the S. muricatum mitogenome and those of 39 other taxa (including 25 Solanaceae species) revealed the evolutionary and taxonomic status of S. muricatum. CONCLUSIONS: We provide the first report of the assembled and annotated S. muricatum mitogenome. This information will help to lay the groundwork for future research on the evolutionary biology of Solanaceae species. Furthermore, the results will assist the development of molecular breeding strategies for S. muricatum based on the most beneficial agronomic traits of this species.

In silico analysis of the wheat BBX gene family and identification of candidate genes for seed dormancy and germination.

BACKGROUND: B-box (BBX) proteins are a type of zinc finger proteins containing one or two B-box domains. They play important roles in development and diverse stress responses of plants, yet their roles in wheat remain unclear. RESULTS: In this study, 96 BBX genes were identified in the wheat genome and classified into five subfamilies. Subcellular localization prediction results showed that 68 TaBBXs were localized in the nucleus. Protein interaction prediction analysis indicated that interaction was one way that these proteins exerted their functions. Promoter analysis indicated that TaBBXs may play important roles in light signal, hormone, and stress responses. qRT-PCR analysis revealed that 14 TaBBXs were highly expressed in seeds compared with other tissues. These were probably involved in seed dormancy and germination, and their expression patterns were investigated during dormancy acquisition and release in the seeds of wheat varieties Jing 411 and Hongmangchun 21, showing significant differences in seed dormancy and germination phenotypes. Subcellular localization analysis confirmed that the three candidates TaBBX2-2 A, TaBBX4-2 A, and TaBBX11-2D were nuclear proteins. Transcriptional self-activation experiments further demonstrated that TaBBX4-2A was transcriptionally active, but TaBBX2-2A and TaBBX11-2D were not. Protein interaction analysis revealed that TaBBX2-2A, TaBBX4-2A, and TaBBX11-2D had no interaction with each other, while TaBBX2-2A and TaBBX11-2D interacted with each other, indicating that TaBBX4-2A may regulate seed dormancy and germination by transcriptional regulation, and TaBBX2-2A and TaBBX11-2D may regulate seed dormancy and germination by forming a homologous complex. CONCLUSIONS: In this study, the wheat BBX gene family was identified and characterized at the genomic level by bioinformatics analysis. These observations provide a theoretical basis for future studies on the functions of BBXs in wheat and other species.

The composition and function of bacterial communities in Bombyx mori (Lepidoptera: Bombycidae) changed dramatically with infected fungi: A new potential to culture Cordyceps cicadae.

Cordyceps cicadae (Hypocreales: Cordycipitaceae) is a renowned entomopathogenic fungus used as herbal medicine in China. However, wild C. cicadae resources have been threatened by heavy harvesting. We hypothesised that Bombyx mori L. (Lepidoptera: Bombycidae) could be a new alternative to cultivate C. cicadae due to the low cost of rearing. Bacterial communities are crucial for the formation of Cordyceps and for promoting the production of metabolites. To better understand the bacterial community structure associated with Cordyceps, three Claviciptaceae fungi were used to explore the pathogenicity of the silkworms. Here, fifth-instar silkworms were infected with C. cicadae, Cordyceps cateniannulata (Hypocreales: Cordycipitaceae) and Beauveria bassiana (Hypocreales: Cordycipitaceae). Subsequently, we applied high-throughput sequencing to explore the composition of bacterial communities in silkworms. Our results showed that all three fungi were highly pathogenic to silkworms, which suggests that silkworms have the potential to cultivate Cordyceps. After fungal infection, the diversity of bacterial communities in silkworms decreased significantly, and the abundance of Staphylococcus increased in mummified larvae, which may play a role in the death process when the host suffers infection by entomopathogenic fungi. Furthermore, there were high similarities in the bacterial community composition and function in the C. cicadae and C. cateniannulata infected samples, and the phylogenetic analysis suggested that these similarities may be related to the fungal phylogenetic relationship. Our findings reveal that infection with different entomopathogenic fungi affects the composition and function of bacterial communities in silkworms and that the bacterial species associated with Cordyceps are primarily host dependent, while fungal infection affects bacterial abundance.

Rediscovering the unusual, solitary bryozoan Monobryozoon ambulans Remane, 1936: first molecular and new morphological data clarify its phylogenetic position.

BACKGROUND: One of the most peculiar groups of the mostly colonial phylum Bryozoa is the taxon Monobryozoon, whose name already implies non-colonial members of the phylum. Its peculiarity and highly unusual lifestyle as a meiobenthic clade living on sand grains has fascinated many biologists. In particular its systematic relationship to other bryozoans remains a mystery. Despite numerous searches for M. ambulans in its type locality Helgoland, a locality with a long-lasting marine station and tradition of numerous courses and workshops, it has never been reencountered until today. Here we report the first observations of this almost mythical species, Monobryozoon ambulans. RESULTS: For the first time since 1938, we present new modern, morphological analyses of this species as well as the first ever molecular data. Our detailed morphological analysis confirms most previous descriptions, but also ascertains the presence of special ambulatory polymorphic zooids. We consider these as bud anlagen that ultimately consecutively separate from the animal rendering it pseudo-colonial. The remaining morphological data show strong ties to alcyonidioidean ctenostome bryozoans. Our morphological data is in accordance with the phylogenomic analysis, which clusters it with species of Alcyonidium as a sister group to multiporate ctenostomes. Divergence time estimation and ancestral state reconstruction recover the solitary state of M. ambulans as a derived character that probably evolved in the Late Cretaceous. In this study, we also provide the entire mitogenome of M. ambulans, which-despite the momentary lack of comparable data-provides important data of a unique and rare species for comparative aspects in the future. CONCLUSIONS: We were able to provide first sequence data and modern morphological data for the unique bryozoan, M. ambulans, which are both supporting an alcyonidioidean relationship within ctenostome bryozoans.

Macroecological correlates of Darwinian shortfalls across terrestrial vertebrates.

Most described species have not been explicitly included in phylogenetic trees-a problem named the Darwinian shortfall-owing to a lack of molecular and/or morphological data, thus hampering the explicit incorporation of evolution into large-scale biodiversity analyses. We investigate potential drivers of the Darwinian shortfall in tetrapods, a group in which at least one-third of described species still lack phylogenetic data, thus necessitating the imputation of their evolutionary relationships in fully sampled phylogenies. We show that the number of preserved specimens in scientific collections is the main driver of phylogenetic knowledge accumulation, highlighting the major role of biological collections in unveiling novel biodiversity data and the importance of continued sampling efforts to reduce knowledge gaps. Additionally, large-bodied and wide-ranged species, as well as terrestrial and aquatic amphibians and reptiles, are phylogenetically better known. Future efforts should prioritize phylogenetic research on organisms that are narrow-ranged, small-bodied and underrepresented in scientific collections, such as fossorial species. Addressing the Darwinian shortfall will be imperative for advancing our understanding of evolutionary drivers shaping biodiversity patterns and implementing comprehensive conservation strategies.

Molecular identification and phylogenetic relationship of fishes belonging to the Family Danionidae from Brahmaputra Basin, Assam, Northeast India.

BACKGROUD: The Northeast India, being part of two global biodiversity hotspot namely the Indo-Burma and Eastern Himalayan Hotspots supports a wide variety of rich aquatic biodiversity including fishes. The family Danionidae is a widely diverse group inhabiting the upper colder stretches of river although few are abundant in the lower stretches. The persisting similarity in the morphological appearance and body colouration within the members of this family seeks an integrated method to identify the species correctly. METHODS AND RESULTS: In the present study, the mt-DNA barcode was generated for correct identification and confirmation of the species. A total of nine mitochondrial cytochrome c oxidase subunit I gene sequences were generated for each species under the study. The pairwise distance values ranged from 0.09 to 9.11% within species and 9.06-32.71% between species. A neighbour-joining tree was constructed based on the Kimura 2 parameter model. Two major groups were observed where Danioninae formed a sister group to the Chedrinae and Rasborinae. CONCLUSION: The present study is a preliminary work to document and identify the species under the family Danionidae from Brahmaputra basin, Assam, using molecular tools and establish the phylogenetic relationship.

Comparative complete chloroplast genome of Geum japonicum: evolution and phylogenetic analysis.

Geum japonicum (Rosaceae) has been widely used in China as a traditional herbal medicine due to its high economic and medicinal value. However, the appearance of Geum species is relatively similar, making identification difficult by conventional phenotypic methods, and the studies of genomics and species evolution are lacking. To better distinguish the medicinal varieties and fill this gap, we carried out relevant research on the chloroplast genome of G. japonicum. Results show a typical quadripartite structure of the chloroplast genome of G. japonicum with a length of 156,042 bp. There are totally 131 unique genes in the genome, including 87 protein-coding genes, 36 tRNA genes, and 8 rRNA genes, and there were also 87 SSRs identified and mostly mononucleotide Adenine-Thymine. We next compared the plastid genomes among four Geum species and obtained 14 hypervariable regions, including ndhF, psbE, trnG-UCC, ccsA, trnQ-UUG, rps16, psbK, trnL-UAA, ycf1, ndhD, atpA, petN, rps14, and trnK-UUU. Phylogenetic analysis revealed that G. japonicum is most closely related to Geum aleppicum, and possibly has some evolutionary relatedness with an ancient relic plant Taihangia rupestris. This research enriched the genome resources and provided fundamental insights for evolutionary studies and the phylogeny of Geum.

The complete chloroplast genome sequence and phylogenetic relationship analysis of Eomecon chionantha, one species unique to China.

Eomecon chionantha Hance, an endemic species in China, has a long medical history in Chinese ethnic minority medicine and is known for its anti-inflammatory and analgesic effects. However, studies of E. chionantha are lacking. In this study, we investigated the characteristics of the E. chionantha chloroplast genome and determined the taxonomic position of E. chionantha in Papaveraceae via phylogenetic analysis. In addition, we determined molecular markers to identify E. chionantha at the molecular level by comparing the chloroplast genomes of E. chionantha and its closely related species. The complete chloroplast genomic information indicated that E. chionantha chloroplast DNA (178,808 bp) contains 99 protein-coding genes, 8 rRNAs, and 37 tRNAs. Meanwhile, we were able to identify a total of 54 simple sequence repeats through our analysis. Our findings from the phylogenetic analysis suggest that E. chionantha shares a close relationship with four distinct species, namely Macleaya microcarpa, Coreanomecon hylomeconoides, Hylomecon japonica, and Chelidonium majus. Additionally, using the Kimura two-parameter model, we successfully identified five hypervariable regions (ycf4-cemA, ycf3-trnS-GGA, trnC-GCA-petN, rpl32-trnL-UAG, and psbI-trnS-UGA). To the best of our knowledge, this is the first report of the complete chloroplast genome of E. chionantha, providing a scientific reference for further understanding of E. chionantha from the perspective of the chloroplast genome and establishing a solid foundation for the future identification, taxonomic determination and evolutionary analysis of this species.

Comparative Analysis of the Complete Mitochondrial Genomes of Three Sisoridae (Osteichthyes, Siluriformes) and the Phylogenetic Relationships of Sisoridae.

The family Sisoridae is one of the largest and most diverse Asiatic catfish families, with most species occurring in the water systems of the Qinhai-Tibetan Plateau and East Himalayas. At present, the phylogenetic relationship of the Sisoridae is relatively chaotic. In this study, the mitochondrial genomes (mitogenomes) of three species Creteuchiloglanis kamengensis, Glaridoglanis andersonii, and Exostoma sp. were systematically investigated, the phylogenetic relationships of the family were reconstructed and to determine the phylogenetic position of Exostoma sp. within Sisoridae. The lengths of the mitogenomes' sequences of C. kamengensis, G. andersonii, and Exostoma sp. were 16,589 bp, 16,531 bp, and 16,529 bp, respectively. They all contained one identical control region (D-loop), two ribosomal RNAs (rRNAs), 13 protein-coding genes (PCGs) and 22 transfer RNA (tRNA) genes. We applied two approaches, Bayesian Inference (BI) and Maximum Likelihood (ML), to construct phylogenetic trees. Our findings revealed that the topological structure of both ML and BI trees exhibited significant congruence. Specifically, the phylogenetic tree strongly supports the monophyly of Sisorinae and Glyptosternoids and provides new molecular biological data to support the reconstruction of phylogenetic relationships with Sisoridae. This study is of great scientific value for phylogenetic and genetic variation studies of the Sisoridae.

The complete mitochondrial genome of cattle tick clade C reveals the genetic relationship within Rhipicephalus microplus complex.

Cattle ticks (Rhipicephalus microplus) are important economic ectoparasites causing direct and indirect damage to cattle and leading to severe economic losses in cattle husbandry. It is common knowledge that R. microplus is a species complex including five clades; however, the relationships within the R. microplus complex remain unresolved. In the present study, we assembled the complete mitochondrial genome of clade C by next-generation sequencing and proved its correctness based on long PCR amplification. It was 15,004 bp in length and consisted of 13 protein genes, 22 transfer genes, and two ribosomal genes located in the two strains. There were two copies of the repeat region (pseudo-nad1 and tRNA-Glu). Data revealed that cox1, cox2, and cox3 genes were conserved within R. microplus with small genetic differences. Ka/Ks ratios suggested that 12 protein genes (excluding nad6) may be neutral selection. The genetic and phylogenetic analyses indicated that clade C was greatly close to clade B. Findings in the current study provided more data for the identification and differentiation of the R. microplus complex and made up for the lack of information about R. microplus clade C.

Evolution and chemical diversity of the volatile compounds in Salvia species.

INTRODUCTION: The plant essential oils are composed of volatile compounds and have significant value in preventing and treating neurological diseases, anxiety, depression, among others. The genus Salvia has been shown to be an important medicinal resource, especially the aerial parts of genus Salvia, which are rich in volatile compounds; however, the chemical diversity and distribution patterns of volatile compounds in Salvia species are still unknown. OBJECTIVE: The work is performed to analyse the chemical diversity and distribution patterns of volatile compounds in genus Salvia. METHODS: The genomic single nucleotide polymorphisms (SNPs) combined with gas chromatography-mass spectrometry (GC-MS) were used to explore the evolution and chemical diversity of Salvia volatile compounds. Initially, the genetic relationship of genus Salvia was revealed by phylogenetic tree that was constructed based on SNPs. And then, GC-MS was applied to explore the chemical diversity of volatile compounds. RESULTS: The results indicated that the volatile compounds were mainly monoterpenoids, sesquiterpenoids, and aliphatic compounds. The genomic SNPs divided species involved in this work into four branches. The volatile compounds in the first and second branches were mainly sesquiterpenoids and monoterpenoids, respectively. Species in the third branch contained more aliphatic compounds and sesquiterpenoids. And those in the fourth branch were also rich in monoterpenoids but had relatively simple chemical compositions. CONCLUSION: This study offered new insights into the phylogenetic relationships besides chemistry diversity and distribution pattern of volatile compounds of genus Salvia, providing theoretical guidance for the investigations and development of secondary metabolites.

Comparative mitochondrial genome analysis of three leafhopper species of the genus Abrus Dai & Zhang (Hemiptera: Cicadellidae: Deltocephalinae) from China with phylogenetic implication.

BACKGROUND: The phylogenetic position and classification of Athysanini are poorly defined, as it includes a large group of polyphyletic genera that have historically been assigned to it mainly because they still exhibit the most typical deltocephaline genitalic and external body characters but lack the distinctive characteristics that other tribes possess. The bamboo-feeding leafhopper genus Abrus belong to the tribe Athysanini of subfamily Deltocephalinae, which currently comprises 19 valid described species, and are limited to the Oriental and Palaearctic regions in China. Although the taxonomy of Abrus are well updated, the references on comparative mitogenomic analyses of Abrus species are only known for a single species. In this study, we sequenced and analyzed the complete mitochondrial genomes (mitogenomes) of Abrus daozhenensis Chen, Yang & Li, 2012 (16,391bp) and A. yunshanensis Chen, Yang & Li, 2012 (15,768bp) (Athysanini), and compared with published mitogenome sequence of A. expansivus Xing & Li, 2014 (15,904bp). RESULTS: These Abrus species shared highly conserved mitogenomes with similar gene order to that of the putative ancestral insect with 37 typical genes and a non-coding A + T-rich region. The nucleotide composition of these genomes is highly biased toward A + T nucleotides (76.2%, 76.3%, and 74.7%), AT-skews (0.091 to 0.095, and 0.095), negative GC-skews (- 0.138, - 0.161, and - 0.138), and codon usage. All 22 tRNA genes had typical cloverleaf secondary structures, except for trnS1 (AGN) which lacks the dihydrouridine arm, and distinctively trnG in the mitogenome of A. expansivus lacks the TψC arm. Phylogenetic analyses based on 13 PCGs, 2 rRNA genes, and 22 tRNA genes consistently recovered the monophyletic Opsiini, Penthimiini, Selenocephalini, Scaphoideini, and Athysanini (except Watanabella graminea, previously sequenced species as Chlorotettix nigromaculatus) based on limited available mitogenome sequence data of 37 species. CONCLUSION: At present, Abrus belongs to the tribe Athysanini based on both morphological and molecular datasets, which is strongly supported in present phylogenetic analyses in both BI and ML methods using the six concatenated datasets: amino acid sequences and nucleotides from different combinations of protein-coding genes (PCGs), ribosomal RNA (rRNAs), and transfer RNA (tRNAs). Phylogenetic trees reconstructed herein based on the BI and ML analyses consistently recovered monophylitic Athysanini, except Watanabella graminea (Athysanini) in Opsiini with high support values.

Sequencing and analysis of the complete mitochondrial genomes of Toona sinensis and Toona ciliata reveal evolutionary features of Toona.

BACKGROUND: Toona is a critical genus in the Meliaceae, and the plants of this group are an asset for both restorative and restorative purposes, the most flexible of which are Toona sinensis and Toona ciliata. To concentrate on the advancement of mitochondrial(Mt) genome variety in T.sinensis and T.ciliata, the Mt genomes of the two species were sequenced in high throughput independently, after de novo assembly and annotation to construct a Mt genome map for comparison in genome structure. Find their repetitive sequences and analyze them in comparison with the chloroplast genome, along with Maximum-likelihood(ML) phylogenetic analysis with 16 other relatives. RESULTS: (1) T. sinensis and T.ciliata are both circular structures with lengths of 683482 bp and 68300 bp, respectively. They share a high degree of similarity in encoding genes and have AT preferences. All of them have the largest Phe concentration and are the most frequently used codons. (2) Both of their Mt genome are highly preserved in terms of structural and functional genes, while the main variability is reflected in the length of tRNA, the number of genes, and the value of RSCU. (3) T. siniensis and T. ciliata were detected to have 94 and 87 SSRs, respectively, of which mononucleotides accounted for the absolute proportion. Besides, the vast majority of their SSRs were found to be poly-A or poly-T. (4)10 and 11 migrating fragments were identified in the comparison with the chloroplast genome, respectively. (5) In the ML evolutionary tree, T.sinensis and T.ciliata clustered individually into a small branch with 100% support, reflecting two species of Toona are very similarly related to each other. CONCLUSIONS: This research provides a basis for the exploitation of T.sinensis and T.ciliata in terms of medicinal, edible, and timber resources to avoid confusion; at the same time, it can explore the evolutionary relationship between the Toona and related species, which does not only have an important practical value, but also provides a theoretical basis for future hybrid breeding of forest trees, molecular markers, and evolutionary aspects of plants, which has great scientific significance.

Comparative analyses and phylogenetic relationships of thirteen Pholidota species (Orchidaceae) inferred from complete chloroplast genomes.

BACKGROUND: The orchid genus Pholidota Lindl. ex Hook. is economically important as some species has long been used in traditional medicine. However, the systematic status of the genus and intergeneric relationships inferred from previous molecular studies are unclear due to insufficient sampling and lack of informative sites. So far, only limited genomic information has been available. The taxonomy of Pholidota remains unresolved and somewhat controversial. In this study, the complete chloroplast (cp.) genomes of thirteen Pholidota species were sequenced and analyzed to gain insight into the phylogeny of Pholidota and mutation patterns in their cp. genomes. RESULTS: All examined thirteen Pholidota cp. genomes exhibited typical quadripartite circular structures, with the size ranging from 158,786 to 159,781 bp. The annotation contained a total of 135 genes in each cp. genome, i.e., 89 protein-coding genes, 38 tRNA genes, and eight rRNA genes. The codon usage analysis indicated the preference of A/U-ending codons. Repeat sequence analysis identified 444 tandem repeats, 322 palindromic repeats and 189 dispersed repeats. A total of 525 SSRs, 13,834 SNPs and 8,630 InDels were detected. Six mutational hotspots were identified as potential molecular markers. These molecular markers and highly variable regions are expected to facilitate future genetic and genomic studies. Our phylogenetic analyses confirmed the polyphyletic status of the genus Pholidota, with species grouped into four main clades: Pholidota s.s. was resolved as the sister to a clade containing species of Coelogyne; the other two clades clustered together with species of Bulleyia and Panisea, respectively; species P. ventricosa was placed at the basal position, deviated from all other species. CONCLUSION: This is the first study to comprehensively examine the genetic variations and systematically analyze the phylogeny and evolution of Pholidota based on plastid genomic data. These findings contribute to a better understanding of plastid genome evolution of Pholidota and provide new insights into the phylogeny of Pholidota and its closely related genera within the subtribe Coelogyninae. Our research has laid the foundation for future studies on the evolutionary mechanisms and classification of this economically and medicinally important genus.

Comparison of plastid genomes and ITS of two sister species in Gentiana and a discussion on potential threats for the endangered species from hybridization.

BACKGROUND: Gentiana rigescens Franchet is an endangered medicinal herb from the family Gentianaceae with medicinal values. Gentiana cephalantha Franchet is a sister species to G. rigescens possessing similar morphology and wider distribution. To explore the phylogeny of the two species and reveal potential hybridization, we adopted next-generation sequencing technology to acquire their complete chloroplast genomes from sympatric and allopatric distributions, as along with Sanger sequencing to produce the nrDNA ITS sequences. RESULTS: The plastid genomes were highly similar between G. rigescens and G. cephalantha. The lengths of the genomes ranged from 146,795 to 147,001 bp in G. rigescens and from 146,856 to 147,016 bp in G. cephalantha. All genomes consisted of 116 genes, including 78 protein-coding genes, 30 tRNA genes, four rRNA genes and four pseudogenes. The total length of the ITS sequence was 626 bp, including six informative sites. Heterozygotes occurred intensively in individuals from sympatric distribution. Phylogenetic analysis was performed based on chloroplast genomes, coding sequences (CDS), hypervariable sequences (HVR), and nrDNA ITS. Analysis based on all the datasets showed that G. rigescens and G. cephalantha formed a monophyly. The two species were well separated in phylogenetic trees using ITS, except for potential hybrids, but were mixed based on plastid genomes. This study supports that G. rigescens and G. cephalantha are closely related, but independent species. However, hybridization was confirmed to occur frequently between G. rigescens and G. cephalantha in sympatric distribution owing to the lack of stable reproductive barriers. Asymmetric introgression, along with hybridization and backcrossing, may probably lead to genetic swamping and even extinction of G. rigescens. CONCLUSION: G. rigescens and G. cephalantha are recently diverged species which might not have undergone stable post-zygotic isolation. Though plastid genome shows obvious advantage in exploring phylogenetic relationships of some complicated genera, the intrinsic phylogeny was not revealed because of matrilineal inheritance here; nuclear genomes or regions are hence crucial for uncovering the truth. As an endangered species, G. rigescens faces serious threats from both natural hybridization and human activities; therefore, a balance between conservation and utilization of the species is extremely critical in formulating conservation strategies.

Genome-wide identification and expression analysis of the AUX/IAA gene family in turnip (Brassica rapa ssp. rapa).

BACKGROUND: Auxin/indoleacetic acid (AUX/IAA) genes encoding short-lived proteins participate in AUX signaling transduction and play crucial roles in plant growth and development. Although the AUX/IAA gene family has been identified in many plants, a systematic analysis of AUX/IAA genes in Brassica rapa ssp. rapa has not yet been reported. RESULTS: We performed a comprehensive genome-wide analysis and found 89 AUX/IAA genes in turnip based on the conserved AUX/IAA domain (pfam02309). Phylogenetic analysis of AUX/IAA genes from turnip, Arabidopsis, and cabbage revealed that these genes cluster into six subgroups (A1, A2, A3, A4, B1, and B2). The motif distribution was also conservative among the internal members of the clade. Enhanced yellow fluorescent protein (EYFP) signals of BrrIAA-EYFPs showed that BrrIAA members functioned as nucleoproteins. Moreover, transcriptional analysis revealed that the expression patterns of AUX/IAA genes in turnip were tissue-dependent. Because orthologs have similar biological functions and interaction networks in plant growth and development, BrrIAA66 in turnip possibly played a role in embryo axis formation, vascular development, lateral root formation, and floral organ development by interacting with BrrARF19 and BrrTIR1. CONCLUSION: These results provide a theoretical basis for further investigation of BrrAUX/IAA genes and lay the foundation for functional analysis of BrrIAA66 in turnip.