Elucidating the origins of phycocyanobilin biosynthesis and phycobiliproteins.

Terrestrial ecosystems and human societies depend on oxygenic photosynthesis, which began to reshape our atmosphere approximately 2.5 billion years ago. The earliest known organisms carrying out oxygenic photosynthesis are the cyanobacteria, which use large complexes of phycobiliproteins as light-harvesting antennae. Phycobiliproteins rely on phycocyanobilin (PCB), a linear tetrapyrrole (bilin) chromophore, as the light-harvesting pigment that transfers absorbed light energy from phycobilisomes to the chlorophyll-based photosynthetic apparatus. Cyanobacteria synthesize PCB from heme in two steps: A heme oxygenase converts heme into biliverdin IXα (BV), and the ferredoxin-dependent bilin reductase (FDBR) PcyA then converts BV into PCB. In the current work, we examine the origins of this pathway. We demonstrate that PcyA evolved from pre-PcyA proteins found in nonphotosynthetic bacteria and that pre-PcyA enzymes are active FDBRs that do not yield PCB. Pre-PcyA genes are associated with two gene clusters. Both clusters encode bilin-binding globin proteins, phycobiliprotein paralogs that we designate as BBAGs (bilin biosynthesis-associated globins). Some cyanobacteria also contain one such gene cluster, including a BBAG, two V4R proteins, and an iron-sulfur protein. Phylogenetic analysis shows that this cluster is descended from those associated with pre-PcyA proteins and that light-harvesting phycobiliproteins are also descended from BBAGs found in other bacteria. We propose that PcyA and phycobiliproteins originated in heterotrophic, nonphotosynthetic bacteria and were subsequently acquired by cyanobacteria.

phyloBARCODER: A web tool for phylogenetic classification of eukaryote metabarcodes using custom reference databases.

We developed phyloBARCODER (https://github.com/jun-inoue/phyloBARCODER), a new web tool that can identify short DNA sequences to the species level using metabarcoding. phyloBARCODER estimates phylogenetic trees based on uploaded anonymous DNA sequences and reference sequences from databases. Without such phylogenetic contexts, alternative, similarity-based methods independently identify species names and anonymous sequences of the same group by pairwise comparisons between queries and database sequences, with the caveat that they must match exactly or very closely. By putting metabarcoding sequences into a phylogenetic context, phyloBARCODER accurately identifies (1) species or classification of query sequences and (2) anonymous sequences associated with the same species or even with populations of query sequences, with clear and accurate explanations. Version 1 of phyloBARCODER stores a database comprising all eukaryotic mitochondrial gene sequences. Moreover, by uploading their own databases, phyloBARCODER users can conduct species identification specialized for sequences obtained from a local geographic region or those of non-mitochondrial genes, e.g., ITS or rbcL.

HGTphyloDetect: facilitating the identification and phylogenetic analysis of horizontal gene transfer.

BACKGROUND: Horizontal gene transfer (HGT) is an important driver in genome evolution, gain-of-function, and metabolic adaptation to environmental niches. Genome-wide identification of putative HGT events has become increasingly practical, given the rapid growth of genomic data. However, existing HGT analysis toolboxes are not widely used, limited by their inability to perform phylogenetic reconstruction to explore potential donors, and the detection of HGT from both evolutionarily distant and closely related species. RESULTS: In this study, we have developed HGTphyloDetect, which is a versatile computational toolbox that combines high-throughput analysis with phylogenetic inference, to facilitate comprehensive investigation of HGT events. Two case studies with Saccharomyces cerevisiae and Candida versatilis demonstrate the ability of HGTphyloDetect to identify horizontally acquired genes with high accuracy. In addition, HGTphyloDetect enables phylogenetic analysis to illustrate a likely path of gene transmission among the evolutionarily distant or closely related species. CONCLUSIONS: The HGTphyloDetect computational toolbox is designed for ease of use and can accurately find HGT events with a very low false discovery rate in a high-throughput manner. The HGTphyloDetect toolbox and its related user tutorial are freely available at https://github.com/SysBioChalmers/HGTphyloDetect.

Deciphering the genetic basis of resistome and virulome diversity among multidrug-resistant Mycoplasma hominis.

Mycoplasma hominis, a commensal bacterium that commonly inhabits the genital tract, leading to infections in both the genitourinary and extragenital regions. However, the antimicrobial resistance and pathogenic mechanisms of M. hominis isolated from extra-urogenital cystic abscess is largely unknown. This study reports the genomic epidemiological characteristics of a M. hominis isolate recovered from a pelvic abscess sample in China. Genomic DNA was extracted and sequenced using Illumina HiSeq X Ten platform. De novo assembly was performed and in silico analysis was accomplished by multiple bioinformatics tools. For phylogenomic analysis, publicly available M. hominis genomes were retrieved from NCBI GenBank database. Whole genome sequencing data showed that the genome size of M. hominis MH4246 was calculated as 679,746 bp, with 558 protein-coding sequences and a G + C content of 26.9%. M. hominis MH4246 is resistant to fluoroquinolones and macrolides, harboring mutations in the quinolone resistance-determining regions (QRDRs) (GyrA S153L, ParC S91I and ParE V417I) and 23S rRNA gene (G280A, C1500T, T1548C and T2218C). Multiple virulence determinants, such as tuf, hlyA, vaa, oppA, MHO_0730 and alr genes, were identified. Phylogenetic analysis showed that the closest relative of M. hominis MH4246 was the strain MH-1 recovered from China, which differed by 3490 SNPs. Overall, this study contributes to the comprehension of genomic characteristics, antimicrobial resistance patterns, and the mechanisms underlying the pathogenicity of this pathogen.

Assembly and comparative analysis of the complete mitochondrial genome of white towel gourd (Luffa cylindrica).

Mitochondria play an important role in the energy production of plant cells through independent genetic systems. This study has aimed to assemble and annotate the functions of the mitochondrial (mt) genome of Luffa cylindrica. The mt genome of L. cylindrica contained two chromosomes with lengths of 380,879 bp and 67,982 bp, respectively. Seventy-seven genes including 39 protein-coding genes, 34 tRNA genes, 3 rRNA genes, and 1 pseudogene, were identified. About 90.63% of the codons ended with A or U bases, and 98.63% of monomers contained A/T, which contributed to the high A/T content (55.91%) of the complete mt genome. Six genes (ATP8, CCMFC, NAD4, RPL10, RPL5 and RPS4) showed positive selection. Phylogenetic analysis indicates that L. cylindrica is closely related to L. acutangula. The present results provide the mt genome of L. cylindrica, which may facilitate possible genetic variation, evolutionary, and molecular breeding studies of L. cylindrica.

Characterizing the etiology of recurrent tuberculosis using whole genome sequencing-Alaska, USA, 2008-2020.

BACKGROUND: Understanding the etiology of recurrent tuberculosis (rTB) is important for effective TB control. Prior to the advent of whole genome sequencing (WGS), attributing rTB to relapse or reinfection using genetic information was complicated by the limited resolution of conventional genotyping methods. METHODS: We applied a systematic method of evaluating whole genome single nucleotide polymorphism (wgSNP) distances and results of phylogenetic analyses to characterize the etiology of rTB in American Indian and Alaska Native (AIAN) persons in Alaska during 2008-2020. We contextualized our findings through descriptive analyses of surveillance data and results of a literature search for investigations that characterized rTB etiology using WGS. RESULTS: The percentage of TB cases in AIAN persons in Alaska classified as recurrent episodes (11.8%) was three times the national percentage (3.9%). Of 38 recurrent episodes included in genetic analyses, we attributed 25 (65.8%) to reinfection based on wgSNP distances and phylogenetic analyses; this proportion was the highest among 16 published point estimates identified through the literature search. By comparison, we attributed 11 of 38 (28.9%) and 6 of 38 (15.8%) recurrent episodes to reinfection based on wgSNP distances alone and on conventional genotyping methods, respectively. CONCLUSIONS: WGS and attribution criteria involving genetic distances and patterns of relatedness can provide an effective means of elucidating rTB etiology. Our findings indicate that rTB occurs at high proportions among AIAN persons in Alaska and is frequently attributable to reinfection, reinforcing the importance of active surveillance and control measures to limit the spread of TB disease in Alaskan AIAN communities.

GERDH: an interactive multi-omics database for cross-species data mining in horticultural crops.

Horticultural plants contribute immensely to the quality of human's life. The rapid development of omics studies on horticultural plants has resulted in large volumes of valuable growth- and development-related data. Genes that are essential for growth and development are highly conserved in evolution. Cross-species data mining reduces the impact of species heterogeneity and has been extensively used for conserved gene identification. Owing to the lack of a comprehensive database for cross-species data mining using multi-omics data from all horticultural plant species, the current resources in this field are far from satisfactory. Here, we introduce GERDH (https://dphdatabase.com), a database platform for cross-species data mining among horticultural plants, based on 12 961 uniformly processed publicly available omics libraries from more than 150 horticultural plant accessions, including fruits, vegetables and ornamental plants. Important and conserved genes that are essential for a specific biological process can be obtained by cross-species analysis module with interactive web-based data analysis and visualization. Moreover, GERDH is equipped with seven online analysis tools, including gene expression, in-species analysis, epigenetic regulation, gene co-expression, enrichment/pathway and phylogenetic analysis. By interactive cross-species analysis, we identified key genes contributing to postharvest storage. By gene expression analysis, we explored new functions of CmEIN3 in flower development, which was validated by transgenic chrysanthemum analysis. We believe that GERDH will be a useful resource for key gene identification and will allow for omics big data to be more available and accessible to horticultural plant community members.

Fern cell walls and the evolution of arabinogalactan proteins in streptophytes.

Significant changes have occurred in plant cell wall composition during evolution and diversification of tracheophytes. As the sister lineage to seed plants, knowledge on the cell wall of ferns is key to track evolutionary changes across tracheophytes and to understand seed plant-specific evolutionary innovations. Fern cell wall composition is not fully understood, including limited knowledge of glycoproteins such as the fern arabinogalactan proteins (AGPs). Here, we characterize the AGPs from the leptosporangiate fern genera Azolla, Salvinia, and Ceratopteris. The carbohydrate moiety of seed plant AGPs consists of a galactan backbone including mainly 1,3- and 1,3,6-linked pyranosidic galactose, which is conserved across the investigated fern AGPs. Yet, unlike AGPs of angiosperms, those of ferns contained the unusual sugar 3-O-methylrhamnose. Besides terminal furanosidic arabinose, Ara (Araf), the main linkage type of Araf in the ferns was 1,2-linked Araf, whereas in seed plants 1,5-linked Araf is often dominating. Antibodies directed against carbohydrate epitopes of AGPs supported the structural differences between AGPs of ferns and seed plants. Comparison of AGP linkage types across the streptophyte lineage showed that angiosperms have rather conserved monosaccharide linkage types; by contrast bryophytes, ferns, and gymnosperms showed more variability. Phylogenetic analyses of glycosyltransferases involved in AGP biosynthesis and bioinformatic search for AGP protein backbones revealed a versatile genetic toolkit for AGP complexity in ferns. Our data reveal important differences across AGP diversity of which the functional significance is unknown. This diversity sheds light on the evolution of the hallmark feature of tracheophytes: their elaborate cell walls.

Substrate selectivity and unique sequence signatures in SWEET/semiSWEET homologs of four taxonomic groups: Sequence analysis and phylogenetic studies.

The recently discovered SWEET (Sugar Will Eventually be Exported Transporter) proteins are involved in the selective transport of monosaccharides and disaccharides. The prokaryotic counterparts, semiSWEETs, form dimers with each monomer forming a triple-helix transmembrane bundle (THB). The longer eukaryotic SWEETs have seven transmembrane helices with two THBs and a linker helix. Structures of semiSWEETs/SWEETs have been determined experimentally. Experimental studies revealed the role of plant SWEETs in vital physiological processes and identified residues responsible for substrate selectivity. However, SWEETs/semiSWEETs from metazoans and bacteria are not characterized. In this study, we used structure-based sequence alignment and compared more than 2000 SWEET/semiSWEETs from four different taxonomic groups. Conservation of residue/chemical property was examined at all positions. Properties of clades/subclades of phylogenetic trees from each taxonomic group were analyzed. Conservation pattern of known residues in the selectivity-filter was used to predict the substrate preference of plant SWEETs and some clusters of metazoans and bacteria. Some residues at the gating and substrate-binding regions, pore-facing positions and at the helix-helix interface are conserved across all taxonomic groups. Conservation of polar/charged residues at specific pore-facing positions, helix-helix interface and in loops seems to be unique for plant SWEETs. Overall, the number of conserved residues is less in metazoan SWEETs. Plant and metazoan SWEETs exhibit high conservation of four and three proline residues respectively in "proline tetrad." Further experimental studies can validate the predicted substrate selectivity and significance of conserved polar/charged/aromatic residues at structurally and functionally important positions of SWEETs/semiSWEETs in plants, metazoans and bacteria.

Vertical Transmission of Hepatitis C Virus Among Women with a History of Injection Opioid Use.

We evaluated vertical transmission and linkage to care in women with HCV and history of injection drug use employing co-localized testing and treatment. Transmission occurred in 1 of 23 infants, with mother-infant genetic distance of 1.26%. Rates for infant testing, maternal linkage and cure were 77%, 52%, and 100%, respectively.

Solanum aculeatissimum and Solanum torvum chloroplast genome sequences: a comparative analysis with other Solanum chloroplast genomes.

BACKGROUND: Solanum aculeatissimum and Solanum torvum belong to the Solanum species, and they are essential plants known for their high resistance to diseases and adverse conditions. They are frequently used as rootstocks for grafting and are often crossbred with other Solanum species to leverage their resistance traits. However, the phylogenetic relationship between S. aculeatissimum and S. torvum within the Solanum genus remains unclear. Therefore, this paper aims to sequence the complete chloroplast genomes of S. aculeatissimum and S. torvum and analyze them in comparison with 29 other previously published chloroplast genomes of Solanum species. RESULTS: We observed that the chloroplast genomes of S. aculeatissimum and S. torvum possess typical tetrameric structures, consisting of one Large Single Copy (LSC) region, two reverse-symmetric Inverted Repeats (IRs), and one Small Single Copy (SSC) region. The total length of these chloroplast genomes ranged from 154,942 to 156,004 bp, with minimal variation. The highest GC content was found in the IR region, while the lowest was in the SSC region. Regarding gene content, the total number of chloroplast genes and CDS genes remained relatively consistent, ranging from 128 to 134 and 83 to 91, respectively. Nevertheless, there was notable variability in the number of tRNA genes and rRNAs. Relative synonymous codon usage (RSCU) analysis revealed that both S. aculeatissimum and S. torvum preferred codons that utilized A and U bases. Analysis of the IR boundary regions indicated that contraction and expansion primarily occurred at the junction between SSC and IR regions. Nucleotide polymorphism analysis and structural variation analysis demonstrated that chloroplast variation in Solanum species mainly occurred in the LSC and SSC regions. Repeat sequence analysis revealed that A/T was the most frequent base pair in simple repeat sequences (SSR), while Palindromic and Forward repeats were more common in long sequence repeats (LSR), with Reverse and Complement repeats being less frequent. Phylogenetic analysis indicated that S. aculeatissimum and S. torvum belonged to the same meristem and were more closely related to Cultivated Eggplant. CONCLUSION: These findings enhance our comprehension of chloroplast genomes within the Solanum genus, offering valuable insights for plant classification, evolutionary studies, and potential molecular markers for species identification.

Sequence characteristics, genetic diversity and phylogenetic analysis of the Cucurbita ficifolia (Cucurbitaceae) chloroplasts genome.

BACKGROUND: Curcubita ficifolia Bouché (Cucurbitaceae) has high value as a food crop and medicinal plant, and also has horticultural value as rootstock for other melon species. China is home to many different cultivars, but the genetic diversity of these resources and the evolutionary relationships among them, as well as the differences between C. ficifolia and other Cucurbita species, remain unclear. RESULTS: We investigated the chloroplast (cp) genomes of 160 C. ficifolia individuals from 31 populations in Yunnan, a major C. ficifolia production area in China. We found that the cp genome of C. ficifolia is ~151 kb and contains 128 genes, of which 86 are protein coding genes, 34 encode tRNA, and eight encode rRNAs. We also identified 64 SSRs, mainly AT repeats. The cp genome was found to contain a total of 204 SNP and 57 indels, and a total of 21 haplotypes were found in the 160 study individuals. The reverse repeat (IR) region of C. ficifolia contained a few differences compared with this region in the six other Cucurbita species. Sequence difference analysis demonstrated that most of the variable regions were concentrated in the single copy (SC) region. Moreover, the sequences of the coding regions were found to be more similar among species than those of the non-coding regions. The phylogenies reconstructed from the cp genomes of 61 representative species of Cucurbitaceae reflected the currently accepted classification, in which C. ficifolia is sister to the other Cucurbita species, however, different interspecific relationships were found between Cucurbita species. CONCLUSIONS: These results will be valuable in the classification of C. ficifolia genetic resources and will contribute to our understanding of evolutionary relationships within the genus Cucurbita.

Investigating pigeon circovirus infection in a pigeon farm: molecular detection, phylogenetic analysis and complete genome analysis.

BACKGROUND: Pigeon circovirus infections in pigeons (Columba livia domestica) have been reported worldwide. Pigeons should be PiCV-free when utilized as qualified experimental animals. However, pigeons can be freely purchased as experimental animals without any clear guidelines to follow. Herein, we investigated the status quo of PiCV infections on a pigeon farm in Beijing, China, which provides pigeons for experimental use. RESULTS: PiCV infection was verified in at least three types of tissues in all forty pigeons tested. A total of 29 full-length genomes were obtained and deposited in GenBank. The whole genome sequence comparison among the 29 identified PiCV strains revealed nucleotide homologies of 85.8-100%, and these sequences exhibited nucleotide homologies of 82.7-98.9% as compared with those of the reference sequences. The cap gene displayed genetic diversity, with a wide range of amino acid homologies ranging from 64.5% to 100%. Phylogenetic analysis of the 29 full-genome sequences revealed that the PiCV strains in this study could be further divided into four clades: A (17.2%), B (10.4%), C (37.9%) and D (34.5%). Thirteen recombination events were also detected in 18 out of the 29 PiCV genomes obtained in this study. Phylogenetic research using the rep and cap genes verified the recombination events, which occurred between clades A/F, A/B, C/D, and B/D among the 18 PiCV strains studied. CONCLUSIONS: In conclusion, PiCV infection, which is highly genetically varied, is extremely widespread on pigeon farms in Beijing. These findings indicate that if pigeons are to be used as experimental animals, it is necessary to evaluate the impact of PiCV infection on the results.

Phylogeny and expression patterns of ERF genes that are potential reproductive inducers in hybrid larch.

BACKGROUND: Larch is an important component of northern forests and a major cultivated tree species in restoration of forest cover using improved seed material. In recent years, the continuous low seed production has severely affected the production of improved variety seedlings and natural regeneration. However, research on the reproductive growth of gymnosperms is extremely scarce. RESULTS: In this study, based on differential transcriptome analysis of two asexual reproductive phases, namely high-yield and low-yield, we further screened 5 ERF family genes that may affect the reproductive development of larch. We analyzed their genetic relationships and predicted their physicochemical properties. The expression patterns of these genes were analyzed in different tissues, developmental stages, hormone treatments, and environmental conditions in hybrid larch. CONCLUSION: The results showed that all 5 genes were induced by low temperature and ABA, and their expression patterns in different tissues suggested a suppressive role in the development of female cones in larch. Among them, LkoERF3-like1 and LkoERF071 may be involved in the flowering age pathway. This study enriches the scarce research on reproductive development in gymnosperms and provides a theoretical basis and research direction for regulating the reproductive development of larch in seed orchards.

Genome-wide characterization, functional analysis, and expression profiling of the Aux/IAA gene family in spinach.

BACKGROUND: The auxin/indole-3-acetic acid (Aux/IAA) gene family is a crucial element of the auxin signaling pathway, significantly influencing plant growth and development. Hence, we conducted a comprehensive investigation of Aux/IAAs gene family using the Sp75 and Monoe-Viroflay genomes in spinach. RESULTS: A total of 24 definitive Aux/IAA genes were identified, exhibiting diverse attributes in terms of amino acid length, molecular weight, and isoelectric points. This diversity underscores potential specific roles within the family, such as growth regulation and stress response. Structural analysis revealed significant variations in gene length and molecular weight. These variations indicate distinct roles within the Aux/IAA gene family. Chromosomal distribution analysis exhibited a dispersed pattern, with chromosomes 4 and 1 hosting the highest and lowest numbers of Aux/IAA genes, respectively. Phylogenetic analysis grouped the identified genes into distinct clades, revealing potential evolutionary relationships. Notably, the phylogenetic tree highlighted specific gene clusters suggesting shared genetic ancestry and potential functional synergies within spinach. Expression analysis under NAA treatment unveiled gene-specific and time-dependent responses, with certain genes exhibiting distinct temporal expression patterns. Specifically, SpoIAA5 displayed a substantial increase at 2 h post-NAA treatment, while SpoIAA7 and SpoIAA9 demonstrated continuous rises, peaking at the 4-hour time point. CONCLUSIONS: These observations indicate a complex interplay of gene-specific and temporal regulation in response to auxin. Moreover, the comparison with other plant species emphasized both shared characteristics and unique features in Aux/IAA gene numbers, providing insights into the evolutionary dynamics of this gene family. This comprehensive characterization of Aux/IAA genes in spinach not only establishes the foundation for understanding their specific functions in spinach development but also provides a valuable resource for experimental validation and further exploration of their roles in the intricate network of auxin signaling pathways.

Assembly and comparative analysis of the initial complete mitochondrial genome of Primulina hunanensis (Gesneriaceae): a cave-dwelling endangered plant.

BACKGROUND: Primulina hunanensis, a troglobitic plant within the Primulina genus of Gesneriaceae family, exhibits robust resilience to arid conditions and holds great horticultural potential as an ornamental plant. The work of chloroplast genome (cpDNA) has been recently accomplished, however, the mitochondrial genome (mtDNA) that is crucial for plant evolution has not been reported. RESULTS: In this study, we sequenced and assembled the P. hunanensis complete mtDNA, and elucidated its evolutionary and phylogenetic relationships. The assembled mtDNA spans 575,242 bp with 43.54% GC content, encompassing 60 genes, including 37 protein-coding genes (PCGs), 20 tRNA genes, and 3 rRNA genes. Notably, high number of repetitive sequences in the mtDNA and substantial sequence translocation from chloroplasts to mitochondria were observed. To determine the evolutionary and taxonomic positioning of P. hunanensis, a phylogenetic tree was constructed using mitochondrial PCGs from P. hunanensis and 32 other taxa. Furthermore, an exploration of PCGs relative synonymous codon usage, identification of RNA editing events, and an investigation of collinearity with closely related species were conducted. CONCLUSIONS: This study reports the initial assembly and annotation of P. hunanensis mtDNA, contributing to the limited mtDNA repository for Gesneriaceae plants and advancing our understanding of their evolution for improved utilization and conservation.

Comparative analysis of the organelle genomes of Aconitum carmichaelii revealed structural and sequence differences and phylogenetic relationships.

In this study, we conducted an assembly and analysis of the organelle genomes of Aconitum carmichaelii. Our investigation encompassed the examination of organelle genome structures, gene transfer events, and the environmental selection pressures affecting A. carmichaelii. The results revealed distinct evolutionary patterns in the organelle genomes of A. carmichaelii. Especially, the plastome exhibited a more conserved structure but a higher nucleotide substitution rate (NSR), while the mitogenome displayed a more complex structure with a slower NSR. Through homology analysis, we identified several instances of unidirectional protein-coding genes (PCGs) transferring from the plastome to the mitogenome. However, we did not observe any events which genes moved from the mitogenome to the plastome. Additionally, we observed multiple transposable element (TE) fragments in the organelle genomes, with both organelles showing different preferences for the type of nuclear TE insertion. Divergence time estimation suggested that rapid differentiation occurred in Aconitum species approximately 7.96 million years ago (Mya). This divergence might be associated with the reduction in CO2 levels and the significant uplift of the Qinghai-Tibet Plateau (QTP) during the late Miocene. Selection pressure analysis indicated that the dN/dS values of both organelles were less than 1, suggested that organelle PCGs were subject to purification selection. However, we did not detect any positively selected genes (PSGs) in Subg. Aconitum and Subg. Lycoctonum. This observation further supports the idea that stronger negative selection pressure on organelle genes in Aconitum results in a more conserved amino acid sequence. In conclusion, this study contributes to a deeper understanding of organelle evolution in Aconitum species and provides a foundation for future research on the genetic mechanisms underlying the structure and function of the Aconitum plastome and mitogenome.

Elucidating the multichromosomal structure within the Brasenia schreberi mitochondrial genome through assembly and analysis.

Brasenia schreberi, a plant species traditionally utilized in Chinese medicine and cuisine, represents an early evolutionary stage among flowering plants (angiosperms). While the plastid genome of this species has been published, its mitochondrial genome (mitogenome) has not been extensively explored, with a notable absence of thorough comparative analyses of its organellar genomes. In our study, we had assembled the entire mitogenome of B. schreberi utilizing the sequencing data derived from both Illumina platform and Oxford Nanopore. The B. schreberi mitogenome mostly exists as six circular DNA molecules, with the largest being 628,257 base pairs (bp) and the smallest 110,220 bp, amounting to 1.49 megabases (Mb). Then we annotated the mitogenome of B. schreberi. The mitogenome encompasses a total of 71 genes: 40 of these are coding proteins genes (PCGs), 28 are genes for transfer RNA (tRNA), and the remaining 3 are genes for ribosomal RNA (rRNA). In the analysis of codon usage, we noted a unique codon preference specific to each amino acid. The most commonly used codons exhibited an average RSCU of 1.36, indicating a noticeable bias in codon selection. In the repeat sequence analysis, a total of 553 simple sequence repeats (SSRs) were identified, 1,822 dispersed repeats (comprising 1,015 forward and 807 palindromic repeats), and 608 long terminal repeats (LTRs). Additionally, in the analysis of homologous sequences between organelle genomes, we detected 38 homologous sequences derived from the plastid genome, each exceeding 500 bp, within the B. schreberi mitochondrial genome. Notably, ten tRNA genes (trnC-GCA, trnM-CAU, trnI-CAU, trnQ-UUG, trnN-GUU, trnT-GGU, trnW-CCA, trnA-UGC, trnI-GAU, and trnV-GAC) appear to have been completely transferred from the chloroplast to the mitogenome. Utilizing the Deepred-mt to predict the RNA editing sites in the mitogenome, we have identified 675 high-quality RNA editing sites in the 40 mitochondrial PCGs. In the final stage of our study, we performed an analysis of colinearity and inferred the phylogenetic relationship of B. schreberi with other angiosperms, utilizing the mitochondrial PCGs as a basis. The results showed that the non-coding regions of the B. schreberi mitogenome are characterized by an abundance of repetitive sequences and exogenous sequences, and B. schreberi is more closely related with Euryale ferox.

Genome-wide analysis of filamentous temperature-sensitive H protease (ftsH) gene family in soybean.

BACKGROUND: The filamentous temperature-sensitive H protease (ftsH) gene family belongs to the ATP-dependent zinc metalloproteins, and ftsH genes play critical roles in plant chloroplast development and photosynthesis. RESULTS: In this study, we performed genome-wide identification and a systematic analysis of soybean ftsH genes. A total of 18 GmftsH genes were identified. The subcellular localization was predicted to be mainly in cell membranes and chloroplasts, and the gene structures, conserved motifs, evolutionary relationships, and expression patterns were comprehensively analyzed. Phylogenetic analysis of the ftsH gene family from soybean and various other species revealed six distinct clades, all of which showed a close relationship to Arabidopsis thaliana. The members of the GmftsH gene family were distributed on 13 soybean chromosomes, with intron numbers ranging from 3 to 15, 13 pairs of repetitive segment. The covariance between these genes and related genes in different species of Oryza sativa, Zea mays, and Arabidopsis thaliana was further investigated. The transcript expression data revealed that the genes of this family showed different expression patterns in three parts, the root, stem, and leaf, and most of the genes were highly expressed in the leaves of the soybean plants. Fluorescence-based real-time quantitative PCR (qRT-PCR) showed that the expression level of GmftsH genes varied under different stress treatments. Specifically, the genes within this family exhibited various induction levels in response to stress conditions of 4℃, 20% PEG-6000, and 100 mmol/L NaCl. These findings suggest that the GmftsH gene family may play a crucial role in the abiotic stress response in soybeans. It was also found that the GmftsH7 gene was localized on the cell membrane, and its expression was significantly upregulated under 4 ℃ treatment. In summary, by conducting a genome-wide analysis of the GmftsH gene family, a strong theoretical basis is established for future studies on the functionality of GmftsH genes. CONCLUSIONS: This research can potentially serve as a guide for enhancing the stress tolerance characteristics of soybean.

In silico analysis identified bZIP transcription factors genes responsive to abiotic stress in Alfalfa (Medicago sativa L.).

BACKGROUND: Alfalfa (Medicago sativa L.) is the most cultivated forage legume around the world. Under a variety of growing conditions, forage yield in alfalfa is stymied by biotic and abiotic stresses including heat, salt, drought, and disease. Given the sessile nature of plants, they use strategies including, but not limited to, differential gene expression to respond to environmental cues. Transcription factors control the expression of genes that contribute to or enable tolerance and survival during periods of stress. Basic-leucine zipper (bZIP) transcription factors have been demonstrated to play a critical role in regulating plant growth and development as well as mediate the responses to abiotic stress in several species, including Arabidopsis thaliana, Oryza sativa, Lotus japonicus and Medicago truncatula. However, there is little information about bZIP transcription factors in cultivated alfalfa. RESULT: In the present study, 237 bZIP genes were identified in alfalfa from publicly available sequencing data. Multiple sequence alignments showed the presence of intact bZIP motifs in the identified sequences. Based on previous phylogenetic analyses in A. thaliana, alfalfa bZIPs were similarly divided and fell into 10 groups. The physico-chemical properties, motif analysis and phylogenetic study of the alfalfa bZIPs revealed high specificity within groups. The differential expression of alfalfa bZIPs in a suite of tissues indicates that bZIP genes are specifically expressed at different developmental stages in alfalfa. Similarly, expression analysis in response to ABA, cold, drought and salt stresses, indicates that a subset of bZIP genes are also differentially expressed and likely play a role in abiotic stress signaling and/or tolerance. RT-qPCR analysis on selected genes further verified these differential expression patterns. CONCLUSIONS: Taken together, this work provides a framework for the future study of bZIPs in alfalfa and presents candidate bZIPs involved in stress-response signaling.