Codon usage bias and genetic diversity in chloroplast genomes of Elaeagnus species (Myrtiflorae: Elaeagnaceae)

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 chromosome-level Elaeagnus mollis genome and transcriptomes provide insights into genome evolution, glycerolipid and vitamin E biosynthesis in seeds

Elaeagnus mollis, which has seeds with high lipid and vitamin E contents, is a valuable woody oil plant with potential for utilization. Currently, the biosynthesis and regulation mechanism of glycerolipids and vitamin E are still unknown in E. mollis. Here, we present the chromosome-level reference genome of E. mollis (scaffold N50: ~40.66Mbp, genome size: ~591.48Mbp) by integrating short-read, long-read, and Hi-C sequencing platforms. A total of 36,796 protein-coding sequences, mainly located on 14 proto-chromosomes, were predicted. Additionally, two whole genome duplication (WGD) events were suggested to have occurred ~54.07 and ~35.06 million years ago (MYA), with Elaeagnaceae plants probably experiencing both WGD events. Furthermore, the long terminal retrotransposons in E. mollis were active ~0.23MYA, and one of them was inferred to insert into coding sequence of the negative regulatory lipid synthesis gene, EMF2. Through transcriptomic and metabonomic analysis, key genes contributing to the high lipid and vitamin E levels of E. mollis seeds were identified, while miRNA regulation was also considered. This comprehensive work on the E. mollis genome not only provides a solid theoretical foundation and experimental basis for the efficient utilization of seed lipids and vitamin E, but also contributes to the exploration of new genetic resources.

Comparative Analysis of the Codon Usage Pattern in the Chloroplast Genomes of Gnetales Species

Codon usage bias refers to the preferential use of synonymous codons, a widespread phenomenon found in bacteria, plants, and animals. Codon bias varies among species, families, and groups within kingdoms and between genes within an organism. Codon usage bias (CUB) analysis sheds light on the evolutionary dynamics of various species and optimizes targeted gene expression in heterologous host plants. As a significant order of gymnosperms, species within Gnetales possess extremely high ecological and pharmaceutical values. However, comprehensive analyses of CUB within the chloroplast genomes of Gnetales species remain unexplored. A systematic analysis was conducted to elucidate the codon usage patterns in 13 diverse Gnetales species based on the chloroplast genomes. Our results revealed that chloroplast coding sequences (cp CDSs) in 13 Gnetales species display a marked preference for AT bases and A/T-ending codons. A total of 20 predominantly high-frequency codons and between 2 and 7 optimal codons were identified across these species. The findings from the ENC-plot, PR2-plot, and neutrality analyses suggested that both mutation pressure and natural selection exert influence on the codon bias in these 13 Gnetales species, with natural selection emerging as the predominant influence. Correspondence analysis (COA) demonstrated variation in the codon usage patterns among the Gnetales species and indicated mutation pressure is another factor that could impact CUB. Additionally, our research identified a positive correlation between the measure of idiosyncratic codon usage level of conservatism (MILC) and synonymous codon usage order (SCUO) values, indicative of CUB's potential influence on gene expression. The comparative analysis concerning codon usage frequencies among the 13 Gnetales species and 4 model organisms revealed that Saccharomyces cerevisiae and Nicotiana tabacum were the optimal exogenous expression hosts. Furthermore, the cluster and phylogenetic analyses illustrated distinct patterns of differentiation, implying that codons, even with weak or neutral preferences, could affect the evolutionary trajectories of these species. Our results reveal the characteristics of codon usage patterns and contribute to an enhanced comprehension of evolutionary mechanisms in Gnetales species.

Assembly and comparative analysis of the complete mitochondrial genome of Brassica rapa var. Purpuraria

BACKGROUND

Purple flowering stalk (Brassica rapa var. purpuraria) is a widely cultivated plant with high nutritional and medicinal value and exhibiting strong adaptability during growing. Mitochondrial (mt) play important role in plant cells for energy production, developing with an independent genetic system. Therefore, it is meaningful to assemble and annotate the functions for the mt genome of plants independently. Though there have been several reports referring the mt genome of in Brassica species, the genome of mt in B. rapa var. purpuraria and its functional gene variations when compared to its closely related species has not yet been addressed.

RESULTS

The mt genome of B. rapa var. purpuraria was assembled through the Illumina and Nanopore sequencing platforms, which revealed a length of 219,775 bp with a typical circular structure. The base composition of the whole B. rapa var. purpuraria mt genome revealed A (27.45%), T (27.31%), C (22.91%), and G (22.32%). 59 functional genes, composing of 33 protein-coding genes (PCGs), 23 tRNA genes, and 3 rRNA genes, were annotated. The sequence repeats, codon usage, RNA editing, nucleotide diversity and gene transfer between the cp genome and mt genome were examined in the B. rapa var. purpuraria mt genome. Phylogenetic analysis show that B. rapa var. Purpuraria was closely related to B. rapa subsp. Oleifera and B. juncea. Ka/Ks analysis reflected that most of the PCGs in the B. rapa var. Purpuraria were negatively selected, illustrating that those mt genes were conserved during evolution.

CONCLUSIONS

The results of our findings provide valuable information on the B.rapa var. Purpuraria genome, which might facilitate molecular breeding, genetic variation and evolutionary researches for Brassica species in the future.

The complete chloroplast genome of Mussaenda pubescens and phylogenetic analysis

The chloroplast (cp) genome sequence of Mussaenda pubescens, a promising resource that is used as a traditional medicine and drink, is important for understanding the phylogenetic relationships among the Mussaenda family and genetic improvement and reservation. This research represented the first comprehensive description of the morphological characteristics of M. pubescens, as well as an analysis of the complete cp genome and phylogenetic relationship. The results indicated a close relationship between M. pubescens and M. hirsutula based on the morphological characteristics of the flower and leaves. The cp was sequenced using the Illumina NovaSeq 6000 platform. The results indicated the cp genome of M. pubescens spanned a total length of 155,122 bp, including a pair of inverted repeats (IRA and IRB) with a length of 25,871 bp for each region, as well as a large single-copy (LSC) region and a small single-copy (SSC) region with lengths of 85,370 bp and 18,010 bp, respectively. The results of phylogenetic analyses demonstrated that species within the same genus displayed a tendency to group closely together. It was suggested that Antirhea, Cinchona, Mitragyna, Neolamarckia, and Uncaria might have experienced an early divergence. Furthermore, M. hirsutula showed a close genetic connection to M. pubescens, with the two species having partially overlapping distributions in China. This study presents crucial findings regarding the identification, evolution, and phylogenetic research on Mussaenda plants, specifically targeting M. pubescens.

Decoding the complete organelle genomic architecture of Stewartia gemmata: an early-diverging species in Theaceae

BACKGROUND

Theaceae, comprising 300 + species, holds significance in biodiversity, economics, and culture, notably including the globally consumed tea plant. Stewartia gemmata, a species of the earliest diverging tribe Stewartieae, is critical to offer insights into Theaceae's origin and evolutionary history.

RESULT

We sequenced the complete organelle genomes of Stewartia gemmata using short/long reads sequencing technologies. The chloroplast genome (158,406 bp) exhibited a quadripartite structure including the large single-copy region (LSC), a small single-copy region (SSC), and a pair of inverted repeat regions (IRs); 114 genes encoded 80 proteins, 30 tRNAs, and four rRNAs. The mitochondrial genome (681,203 bp) exhibited alternative conformations alongside a monocyclic structure: 61 genes encoding 38 proteins, 20 tRNAs, three rRNAs, and RNA editing-impacting genes, including ATP6, RPL16, COX2, NAD4L, NAD5, NAD7, and RPS1. Comparative analyses revealed frequent recombination events and apparent rRNA gene gains and losses in the mitochondrial genome of Theaceae. In organelle genomes, the protein-coding genes exhibited a strong A/U bias at codon endings; ENC-GC3 analysis implies selection-driven codon bias. Transposable elements might facilitate interorganelle sequence transfer. Phylogenetic analysis confirmed Stewartieae's early divergence within Theaceae, shedding light on organelle genome characteristics and evolution in Theaceae.

CONCLUSIONS

We studied the detailed characterization of organelle genomes, including genome structure, composition, and repeated sequences, along with the identification of lateral gene transfer (LGT) events and complexities. The discovery of a large number of repetitive sequences and simple sequence repeats (SSRs) has led to new insights into molecular phylogenetic markers. Decoding the Stewartia gemmata organellar genome provides valuable genomic resources for further studies in tea plant phylogenomics and evolutionary biology.

Complete chloroplast genome of the Malus baccata var. gracilis provides insights into the evolution and phylogeny of Malus species

Malus baccata (L.) var. gracilis (Rehd.) has high ornamental value and breeding significance, and comparative chloroplast genome analysis was applied to facilitate genetic breeding for desired traits and resistance and provide insight into the phylogeny of this genus. Using data from whole-genome sequencing, a tetrameric chloroplast genome with a length of 159,992 bp and a total GC content of 36.56% was constructed. The M. baccata var. gracilis chloroplast genome consists of a large single-copy sequence (88,100 bp), a short single-copy region (19,186 bp), and two inverted repeat regions, IRa (26,353 bp) and IRb (26,353 bp). This chloroplast genome contains 112 annotated genes, including 79 protein-coding genes (nine multicopy), 29 tRNA genes (eight multicopy), and four rRNA genes (all multicopy). Calculating the relative synonymous codon usage revealed a total of 32 high-frequency codons, and the codons exhibited a biased usage pattern towards A/U as the ending nucleotide. Interspecific sequence comparison and boundary analysis revealed significant sequence variation in the vast single-copy region, as well as generally similar expansion and contraction of the SSC and IR regions for 10 analyzed Malus species. M. baccata var. gracilis and Malus hupehensis were grouped together into one branch based on phylogenetic analysis of chloroplast genome sequences. The chloroplast genome of Malus species provides an important foundation for species identification, genetic diversity analysis, and Malus chloroplast genetic engineering. Additionally, the results can facilitate the use of pendant traits to improve apple tree shape.

Codon usage characterization and phylogenetic analysis of the mitochondrial genome in Hemerocallis citrina

BACKGROUND

Hemerocallis citrina Baroni is a traditional vegetable crop widely cultivated in eastern Asia for its high edible, medicinal, and ornamental value. The phenomenon of codon usage bias (CUB) is prevalent in various genomes and provides excellent clues for gaining insight into organism evolution and phylogeny. Comprehensive analysis of the CUB of mitochondrial (mt) genes can provide rich genetic information for improving the expression efficiency of exogenous genes and optimizing molecular-assisted breeding programmes in H. citrina.

RESULTS

Here, the CUB patterns in the mt genome of H. citrina were systematically analyzed, and the possible factors shaping CUB were further evaluated. Composition analysis of codons revealed that the overall GC (GCall) and GC at the third codon position (GC3) contents of mt genes were lower than 50%, presenting a preference for A/T-rich nucleotides and A/T-ending codons in H. citrina. The high values of the effective number of codons (ENC) are indicative of fairly weak CUB. Significant correlations of ENC with the GC3 and codon counts were observed, suggesting that not only compositional constraints but also gene length contributed greatly to CUB. Combined ENC-plot, neutrality plot, and Parity rule 2 (PR2)-plot analyses augmented the inference that the CUB patterns of the H. citrina mitogenome can be attributed to multiple factors. Natural selection, mutation pressure, and other factors might play a major role in shaping the CUB of mt genes, although natural selection is the decisive factor. Moreover, we identified a total of 29 high-frequency codons and 22 optimal codons, which exhibited a consistent preference for ending in A/T. Subsequent relative synonymous codon usage (RSCU)-based cluster and mt protein coding gene (PCG)-based phylogenetic analyses suggested that H. citrina is close to Asparagus officinalis, Chlorophytum comosum, Allium cepa, and Allium fistulosum in evolutionary terms, reflecting a certain correlation between CUB and evolutionary relationships.

CONCLUSIONS

There is weak CUB in the H. citrina mitogenome that is subject to the combined effects of multiple factors, especially natural selection. H. citrina was found to be closely related to Asparagus officinalis, Chlorophytum comosum, Allium cepa, and Allium fistulosum in terms of their evolutionary relationships as well as the CUB patterns of their mitogenomes. Our findings provide a fundamental reference for further studies on genetic modification and phylogenetic evolution in H. citrina.