Selection Pressure Analysis Research Articles

Molecular evolution of Phytocyanin gene and analysis of expression at different coloring periods in apple (Malus domestica)

BackgroundPC (phytocyanin) is a class of copper-containing electron transfer proteins closely related to plant photosynthesis, abiotic stress responses growth and development in plants, and regulation of the expression of some flavonoids and phenylpropanoids, etc., however, compared with other plants, the PC gene family has not been systematically characterized in apple.ResultsA total of 59 MdPC gene members unevenly distributed across 12 chromosomes were identified at the genome-wide level. The proteins of the MdPC family were classified into four subfamilies based on differences in copper binding sites and glycosylation sites: Apple Early nodulin-like proteins (MdENODLs), Apple Uclacyanin-like proteins (MdUCLs), Apple Stellacyanin-like proteins (MdSCLs), and Apple Plantacyanin-like proteins (MdPLCLs). Some MdPC members with similar gene structures and conserved motifs belong to the same group or subfamily. The internal collinearity analysis revealed 14 collinearity gene pairs among members of the apple MdPC gene. Interspecific collinearity analysis showed that apple had 31 and 35 homologous gene pairs with strawberry and grape, respectively. Selection pressure analysis indicated that the MdPC gene was under purifying selection. Prediction of protein interactions showed that MdPC family members interacted strongly with the Nad3 protein. GO annotation results indicated that the MdPC gene also regulated the biosynthesis of phenylpropanoids. Chip data analysis showed that (MdSCL3, MdSCL7 and MdENODL27) were highly expressed in mature fruits and peels. Many cis-regulatory elements related to light response, phytohormones, abiotic stresses and flavonoid biosynthetic genes regulation were identified 2000 bp upstream of the promoter of the MdPC gene, and qRT-PCR results showed that gene members in Group IV (MdSCL1/3, MdENODL27) were up-regulated at all five stages of apple coloring, but the highest expression was observed at the DAF13 (day after fruit bag removal) stage. The gene members in Group II (MdUCL9, MdPLCL3) showed down-regulated or lower expression in the first four stages of apple coloring but up-regulated and highest expression in the DAF 21 stage.ConclusionHerein, one objective of these findings is to provide valuable information for understanding the structure, molecular evolution, and expression pattern of the MdPC gene, another major objective in this study was designed to lay the groundwork for further research on the molecular mechanism of PC gene regulation of apple fruit coloration.

Sequencing of the Complete Mitochondrial Genome of the Big Brown Mactra Clam, Mactra grandis (Venerida: Mactridae).

Mitochondrial genomes are playing an increasingly important role in molluscan taxonomy, germplasm, and evolution studies. The first complete mitochondrial genome of the commercial big brown mactra clam, Mactra grandis, was characterized using Illumina next-generation sequencing in this study. The 17,289 bp circular genome has a typical gene organization of 13 protein-coding genes (PCGs), 2 rRNAs, and 22 tRNAs, with an obvious (A + T)-bias of 64.54%. All PCGs exhibited a homogeneous bias in nucleotide composition with a (A + T)-bias, a positive GC skew, and a negative AT skew. Results of phylogenetic analysis showed that Mactra grandis was most closely related to Mactra cygnus. The functional gene arrangement of the two species was identical but different from other Mactra species. The congeneric relationships among Mactra species were demonstrated by genetic distance analysis. Additionally, the selective pressure analysis suggested that cox1 was highly efficient for discriminating closely related species in genus Mactra, while nad2 was the most appropriate marker for population genetic analysis.

PCR amplification, sequencing, and in silico analysis of holandric genes TSPYand SRY in Nili Ravi buffalo

The present work has been designed to study the evolutionary perspectives of two important holandric genes, namely, testis-specific protein, Y-encoded (TSPY), and sex determining region, Y-encoded (SRY) in Nili Ravi buffalo. These Y-chromosomal genes are responsible for male sex development. The evolutionary aspects of these genes have not been studied in the Nili Ravi buffaloes from India. A study was carried out on the peripheral blood samples of three adult Nili Ravi bulls belonging to university farm. DNA was extracted using the phenol-chloroform isoamyl alcohol (PCI) technique and reported primer sequences were used to amplify the genes (partially) using a thermocycler. Three biological samples were analyzed using Sanger sequencing and were aligned using Clustal Omega. Separate phylogenetic trees for the genes were constructed using MEGA7 software. The selection pressure analysis and evolutionary divergence were individually studied for each of the genes. The current study reports the partial TSPY and SRY gene sequences of Indian Nili Ravi buffalo. Besides, the molecular evolution analysis indicates that the genes have experienced some purifying selection during the evolutionary process. Molecular characterization of these holandric genes (TSPY and SRY) in Nili Ravi will find applications in sex determination in Nili Ravi and other buffaloes.

Exploring the evolution of CHS gene family in plants.

Chalcone synthase (CHS) is a key enzyme that catalyzes the first committed step of flavonoid biosynthetic pathway. It plays a vital role not only in maintaining plant growth and development, but also in regulating plant response to environmental hazards. However, the systematic phylogenomic analysis of CHS gene family in a wide range of plant species has not been reported yet. To fill this knowledge gap, a large-scale investigation of CHS genes was performed in 178 plant species covering green algae to dicotyledons. A total of 2,011 CHS and 293 CHS-like genes were identified and phylogenetically divided into four groups, respectively. Gene distribution patterns across the plant kingdom revealed the origin of CHS can be traced back to before the rise of algae. The gene length varied largely in different species, while the exon structure was relatively conserved. Selection pressure analysis also indicated the conserved features of CHS genes on evolutionary time scales. Moreover, our synteny analysis pinpointed that, besides genome-wide duplication and tandem duplication, lineage specific transposition events also occurred in the evolutionary trajectory of CHS gene family. This work provides novel insights into the evolution of CHS gene family and may facilitate further research to better understand the regulatory mechanism of traits relating to flavonoid biosynthesis in diverse plants.

Selective Pressure and Evolution of SARS-CoV-2 Lineages BF.7 and BQ.1.1 Circulating in Italy from July to December 2022.

In this work, we studied the selective pressure and evolutionary analysis on the SARS-CoV-2 BF.7 and BQ.1.1 lineages circulating in Italy from July to December 2022. Two different datasets were constructed: the first comprised 694 SARS-CoV-2 BF.7 lineage sequences and the second comprised 734 BQ.1.1 sequences, available in the Italian COVID-19 Genomic (I-Co-Gen) platform and GISAID (last access date 15 December 2022). Alignments were performed with MAFFT v.7 under the Galaxy platform. The HYPHY software was used to study the selective pressure. Four positively selected sites (two in nsp3 and two in the spike) were identified in the BF.7 dataset, and two (one in ORF8 and one in the spike gene) were identified in the BQ.1.1 dataset. Mutation analysis revealed that R408S and N440K are very common in the spike of the BF.7 genomes, as well as L452R among BQ.1.1. N1329D and Q180H in nsp3 were found, respectively, at low and rare frequencies in BF.7, while I121L and I121T were found to be rare in ORF8 for BQ.1.1. The positively selected sites may have been driven by the selection for increased viral fitness, under circumstances of defined selective pressure, as well by host genetic factors.

Phylogenetic relationships, selective pressure and molecular markers development of six species in subfamily Polygonoideae based on complete chloroplast genomes

The subfamily Polygonoideae encompasses a diverse array of medicinal and horticultural plants that hold significant economic value. However, due to the lack of a robust taxonomy based on phylogenetic relationships, the classification within this family is perplexing, and there is also a scarcity of reports on the chloroplast genomes of many plants falling under this classification. In this study, we conducted a comprehensive analysis by sequencing and characterizing the complete chloroplast genomes of six Polygonoideae plants, namely Pteroxygonum denticulatum, Pleuropterus multiflorus, Pleuropterus ciliinervis, Fallopia aubertii, Fallopia dentatoalata, and Fallopia convolvulus. Our findings revealed that these six plants possess chloroplast genomes with a typical quadripartite structure, averaging 162,931 bp in length. Comparative chloroplast analysis, codon usage analysis, and repetitive sequence analysis demonstrated a high level of conservation within the chloroplast genomes of these plants. Furthermore, phylogenetic analysis unveiled a distinct clade occupied by P. denticulatum, while P. ciliinrvis displayed a closer relationship to the three plants belonging to the Fallopia genus. Selective pressure analysis based on maximum likelihood trees showed that a total of 14 protein-coding genes exhibited positive selection, with psbB and ycf1 having the highest number of positive amino acid sites. Additionally, we identified four molecular markers, namely petN-psbM, psal-ycf4, ycf3-trnS-GGA, and trnL-UAG-ccsA, which exhibit high variability and can be utilized for the identification of these six plants.

Genome-Wide Analysis of MYB Gene Family in Chrysanthemum ×morifolium Provides Insights into Flower Color Regulation.

MYBs constitute the second largest transcription factor (TF) superfamily in flowering plants with substantial structural and functional diversity, which have been brought into focus because they affect flower colors by regulating anthocyanin biosynthesis. Up to now, the genomic data of several Chrysanthemum species have been released, which provides us with abundant genomic resources for revealing the evolution of the MYB gene family in Chrysanthemum species. In the present study, comparative analyses of the MYB gene family in six representative species, including C. lavandulifolium, C. seticuspe, C. ×morifolium, Helianthus annuus, Lactuca sativa, and Arabidopsis thaliana, were performed. A total of 1104 MYBs, which were classified into four subfamilies and 35 lineages, were identified in the three Chrysanthemum species (C. lavandulifolium, C. seticuspe, and C. ×morifolium). We found that whole-genome duplication and tandem duplication are the main duplication mechanisms that drove the occurrence of duplicates in CmMYBs (particularly in the R2R3-MYB subfamily) during the evolution of the cultivated chrysanthemums. Sequence structure and selective pressure analyses of the MYB gene family revealed that some of R2R3-MYBs were subjected to positive selection, which are mostly located on the distal telomere segments of the chromosomes and contain motifs 7 and 8. In addition, the gene expression analysis of CmMYBs in different organs and at various capitulum developmental stages of C. ×morifolium indicated that CmMYBS2, CmMYB96, and CmMYB109 might be the negative regulators for anthocyanin biosynthesis. Our results provide the phylogenetic context for research on the genetic and functional evolution of the MYB gene family in Chrysanthemum species and deepen our understanding of the regulatory mechanism of MYB TFs on the flower color of C. ×morifolium.

Identification of Dof transcription factors in Dendrobium huoshanense and expression pattern under abiotic stresses.

Introduction: DNA-binding with one finger (Dof) transcription factors (TFs) are a unique family of TFs found in higher plants that regulate plant responses to light, hormones, and abiotic stresses. The specific involvement of Dof genes in the response to environmental stresses remains unknown in D. huoshanense. Methods: A total of 22 Dof family genes were identified from the D. huoshanense genome. Results: Chromosome location analysis showed that DhDof genes were distributed on 12 chromosomes, with the largest number of Dof genes located on chromosome 8. The phylogenetic tree revealed that DhDofs could be categorized into 11 distinct subgroups. In addition to the common groups, DhDof4, DhDof5, DhDof17, and the AtDof1.4 ortholog were clustered into the B3 subgroup. Group E was a newly identified branch, among which DhDof6, DhDof7, DhDof8, and DhDof9 were in an independent branch. The conserved motifs and gene structure revealed the differences in motif number and composition of DhDofs. The dof domain near the N-terminus was highly conserved and contained a C2-C2-type zinc finger structure linked with four cysteines. Microsynteny and interspecies collinearity revealed gene duplication events and phylogenetic tree among DhDofs. Large-scale gene duplication had not occurred among the DhDofs genes and only in one pair of genes on chromosome 13. Synteny blocks were found more often between D. huoshanense and its relatives and less often between Oryza sativa and Arabidopsis thaliana. Selection pressure analysis indicated that DhDof genes were subject to purifying selection. Expression profiles and correlation analyses revealed that the Dof gene under hormone treatments showed several different expression patterns. DhDof20 and DhDof21 had the highest expression levels and were co-expressed under MeJA induction. The cis-acting element analysis revealed that each DhDof had several regulatory elements involved in plant growth as well as abiotic stresses. qRT-PCR analysis demonstrated that DhDof2 was the main ABA-responsive gene and DhDof7 was the main cold stress-related gene. IAA suppressed the expression of some Dof candidates, and SA inhibited most of the candidate genes. Discussion: Our results may provide new insights for the further investigation of the Dof genes and the screening of the core stress-resistance genes.

Comparative mitochondrial genomics of endemic Mexican vesper yellow bats genus Rhogeessa (Chiroptera: Vespertilionidae) and insights into internal relationships in the family Vespertilionidae

In the species-rich family Vespertilionidae, vesper yellow bats in the genus Rhogeessa include eleven species, three of them endemic to Mexico. These insectivorous bats provide important ecosystem services, including pest control. Even though some aspects of their biology are well- known, only a few genomic resources are available for these species, which limits our understanding of their biology. In this study, we assembled and annotated the mitochondrial genome of four species: R. aenea, R. genowaysi, R. mira, and R. parvula. We generated a phylomitogenomic hypothesis based on translated protein-coding genes for a total of 52 species in the family Vespertilionidae and examined the phylogenetic position of the genus Rhogeessa and species within the family. The AT-rich mitogenomes of R. aenea, R. genowaysi, R. mira, and R. parvula are 16,763, 16,781, 16,807, and 16,794 pb in length, respectively. Each studied mitogenome encodes 13 Protein Coding Genes (PCGs), 22 transfer RNA genes, and 2 rRNA genes, and contains a putative control region (CR). All tRNAs exhibit a ‘cloverleaf’ secondary structure, except tRNA-Serine-1 that lacked the DHU arm in all studied mitogenomes. Selective pressure analyses indicated that all protein-coding genes are exposed to purifying selection. The phylomitogenomic analysis supported the monophyletic status of the family Vespertilionidae, confirmed the placement of Rhogeessa within the tribe Antrozoini, and clarified phylogenetic relationships within and among subfamilies and tribes in this family. Our results indicate that phylomitogenomics are useful to explore the evolutionary history of vesper bats. The assembly and comprehensive analysis of mitochondrial genomes offer the potential to generate molecular references and resources beneficial for genetic analyses aimed at understanding the ecology and evolution of these remarkable bats.

Molecular characterization and evolutionary dynamics of influenza A(H1N1) strains isolated from 2015 to 2017 in North India

Background and Objectives: The influenza A(H1N1) virus is known for large outbreaks, epidemics and pandemics world- wide owing to its genome plasticity which evolves constantly. In the year 2015 and then in 2017, India witnessed an upsurge in cases. Materials and Methods: The study was carried out in this period (2015-2017) with samples from 5 states across north India. The hemagglutinin 1 (HA1) and non-structural 1 (NS1) gene segments of the viral genome were characterised by phylogenetic analysis, selection pressure analysis, prediction of potential glycosylation sites and phylodynamic analysis of the study strains. Results: The study strains belonged to genogroup 6B. A total of 12 mutations were observed, half of which were located on the key receptor binding region of the HA1 protein. Established virulence markers D222G, S183P were observed in 2017 samples. Acquisition of an extra glycosylation site was observed in few strains from 2017 and 2016. Selection pressure analysis found the average dN/dS (v) ratio of 0.2106 and few codon sites in particular showed significant evidence of being under negative selection. Conclusion: The genogroup 6B continues to be the dominant circulating strain in Indian subcontinent region however the presence of pathogenic mutations in the 2017 strains from north India underlines the importance of continued molecular surveillance.

Unveiling kiwifruit TCP genes: evolution, functions, and expression insights

ABSTRACT The TEOSINTE-BRANCHED1/CYCLOIDEA/PROLEFERATING-CELL-FACTORS (TCP) gene family is a plant-specific transcriptional factor family involved in leaf morphogenesis and senescence, lateral branching, hormone crosstalk, and stress responses. To date, a systematic study on the identification and characterization of the TCP gene family in kiwifruit has not been reported. Additionally, the function of kiwifruit TCPs in regulating kiwifruit responses to the ethylene treatment and bacterial canker disease pathogen (Pseudomonas syringae pv. actinidiae, Psa) has not been investigated. Here, we identified 40 and 26 TCP genes in Actinidia chinensis (Ac) and A. eriantha (Ae) genomes, respectively. The synteny analysis of AcTCPs illustrated that whole-genome duplication accounted for the expansion of the TCP family in Ac. Phylogenetic, conserved domain, and selection pressure analysis indicated that TCP family genes in Ac and Ae had undergone different evolutionary patterns after whole-genome duplication (WGD) events, causing differences in TCP gene number and distribution. Our results also suggested that protein structure and cis-element architecture in promoter regions of TCP genes have driven the function divergence of duplicated gene pairs. Three and four AcTCP genes significantly affected kiwifruit responses to the ethylene treatment and Psa invasion, respectively. Our results provided insight into general characters, evolutionary patterns, and functional diversity of kiwifruit TCPs.

Genome-wide analysis and expression profiling of the HD-ZIP gene family in kiwifruit

The homeodomain-leucine zipper (HD-Zip) gene family plays a pivotal role in plant development and stress responses. Nevertheless, a comprehensive characterization of the HD-Zip gene family in kiwifruit has been lacking. In this study, we have systematically identified 70 HD-Zip genes in the Actinidia chinensis (Ac) genome and 55 in the Actinidia eriantha (Ae) genome. These genes have been categorized into four subfamilies (HD-Zip I, II, III, and IV) through rigorous phylogenetic analysis. Analysis of synteny patterns and selection pressures has provided insights into how whole-genome duplication (WGD) or segmental may have contributed to the divergence in gene numbers between these two kiwifruit species, with duplicated gene pairs undergoing purifying selection. Furthermore, our study has unveiled tissue-specific expression patterns among kiwifruit HD-Zip genes, with some genes identified as key regulators of kiwifruit responses to bacterial canker disease and postharvest processes. These findings not only offer valuable insights into the evolutionary and functional characteristics of kiwifruit HD-Zips but also shed light on their potential roles in plant growth and development.

Comparative plastome analysis of the sister genera Ceratocephala and Myosurus (Ranunculaceae) reveals signals of adaptive evolution to arid and aquatic environments

BackgroundExpansion and contraction of inverted repeats can cause considerable variation of plastid genomes (plastomes) in angiosperms. However, little is known about whether structural variations of plastomes are associated with adaptation to or occupancy of new environments. Moreover, adaptive evolution of angiosperm plastid genes remains poorly understood. Here, we sequenced the complete plastomes for four species of xerophytic Ceratocephala and hydrophytic Myosurus, as well as Ficaria verna. By an integration of phylogenomic, comparative genomic, and selection pressure analyses, we investigated evolutionary patterns of plastomes in Ranunculeae and their relationships with adaptation to dry and aquatic habitats.ResultsOwing to the significant contraction of the boundary of IRA/LSC towards the IRA, plastome sizes and IR lengths of Myosurus and Ceratocephala are smaller within Ranunculeae. Compared to other Ranunculeae, the Myosurus plastome lost clpP and rps16, one copy of rpl2 and rpl23, and one intron of rpoC1 and rpl16, and the Ceratocephala plastome added an infA gene and lost one copy of rpl2 and two introns of clpP. A total of 11 plastid genes (14%) showed positive selection, two genes common to Myosurus and Ceratocephala, seven in Ceratocephala only, and two in Myosurus only. Four genes showed strong signals of episodic positive selection. The rps7 gene of Ceratocephala and the rpl32 and ycf4 genes of Myosurus showed an increase in the rate of variation close to 3.3 Ma.ConclusionsThe plastomic structure variations as well as the positive selection of two plastid genes might be related to the colonization of new environments by the common ancestor of Ceratocephala and Myosurus. The seven and two genes under positive selection might be related to the adaptation to dry and aquatic habitats in Ceratocephala and Myosurus, respectively. Moreover, intensified aridity and frequent sea-level fluctuations, as well as global cooling, might have favored an increased rate of change in some genes at about 3.3 Ma, associated with adaptation to dry and aquatic environments, respectively. These findings suggest that changing environments might have influenced structural variations of plastomes and fixed new mutations arising on some plastid genes owing to adaptation to specific habitats.

Investigation of porcine circovirus type 2 and porcine circovirus type 3 infections based on dual TaqMan fluorescent quantitative PCR method and genetic evolutionary analysis of these two viruses.

Porcine circovirus type 2 (PCV2) is the pathogen of Porcine Circovirus Associated Diseases. Porcine circovirus type 3 (PCV3) is a novel porcine circovirus associated with porcine dermatitis and nephropathy syndrome (PDNS) and reproductive failure. PCV2 is clearly pathogenic, while the pathogenicity of PCV3 remains controversial, so it is crucial to monitor the prevalence of PCV2 and PCV3 in healthy and diseased pigs to investigate the effects of PCV3 and PCV2 on the health status of pigs. Here, we developed a PCV2 and PCV3 dual TaqMan quantitative PCR (qPCR) method to test samples from healthy and diseased pigs, to clarify the differences in the positive rates and viral copy numbers of PCV2 and PCV3, and to analyze the genetic evolution and molecular characterization of the viral genomes obtained with sequence alignment and phylogenetic analysis, homology and structural analysis of Cap proteins, and selection pressure analysis. We successfully established a dual TaqMan qPCR method for PCV2 and PCV3 with good repeatability, specificity and sensitivity. In total, 1,385 samples from 15 Chinese provinces were tested with the established qPCR. The total positive rates were 37.47% for PCV3 and 57.95% for PCV2, and the coinfection rate for was 25.49%. The positive rates of PCV3 and PCV2 in 372 healthy pigs were 15.05 and 69.89%, respectively, and the coinfection rate was 12.90%. The positive rates of PCV3 and PCV2 in 246 diseased pigs were 55.69 and 83.33%, respectively, and the coinfection rate was 47.97%. Eighteen PCV3 genomes and 64 PCV2 genomes were identified, including nine each of the PCV3a-1 and PCV3b genotypes, eight of PCV2a, 16 of PCV2b, and 40 of PCV2d. The amino acid identity within the PCV3 Cap proteins was 94.00-100.0%, whereas the PCV2 Cap proteins showed an identity of 81.30-100.0%. PCV3 Cap was most variable at amino acid sites 24, 27, 77, 104 and 150, whereas PCV2 Cap had 10-13 unique sites of variation between genotypes. These results clarify the prevalence and variations of PCV2 and PCV3 in healthy and diseased pigs, which will provide a basis for the prevention and control of the two viral infections.

Genomic Evolution and Selective Pressure Analysis of a Novel Porcine Sapovirus in Shanghai, China.

Porcine sapovirus (PoSaV) is one of the most significant pathogens causing piglet diarrhea, and one with limited genetic characterization. In this study, the prevalence, infection pattern, and genetic evolution of porcine sapovirus were elucidated in detail. The positive rate of PoSaV was 10.1% (20/198), with dual, triple, and quadruple infections of 45%, 40%, and 5%, respectively. To further explore the viral composition in the PoSaV-positive diarrhea feces, metagenomic sequencing was carried out. The results confirmed that RNA viruses accounted for a higher proportion (55.47%), including the two primary viruses of PoSaV (21.78%) and porcine astrovirus (PAstV) (24.54%) in the tested diarrhea feces samples. Afterward, a full-length sequence of the PoSaV isolate was amplified and named SHCM/Mega2023, and also given the identifier of GenBank No. PP388958. Phylogenetic analysis identified the prevalent PoSaV strain SHCM/Mega2023 in the GIII genogroup, involving a recombinant event with MK962338 and KT922089, with the breakpoint at 2969-5132 nucleotides (nt). The time tree revealed that the GIII genogroup exhibits the widest divergence time span, indicating a high likelihood of viral recombination. Moreover, SHCM/Mega2023 had three nucleotide "RPL" insertions at the 151-153 nt site in the VP2 gene, compared to the other GIII strains. Further selective pressure calculations demonstrate that the whole genome of the SHCM/Mega2023 strain was under purifying selection (dN/dS < 1), with seven positively selected sites in the VP1 protein, which might be related to antigenicity. In conclusion, this study presents a novel genomic evolution of PoSaV, offering valuable insights into antigenicity and for vaccine research.

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.

Comparative genomics and phylogenomics of the genus Glycyrrhiza (Fabaceae) based on chloroplast genomes.

Glycyrrhiza (Fabaceae) species are rich in metabolites and widely used in medicine. Research on the chloroplast genome of Glycyrrhiza is important for understanding its phylogenetics, biogeography, genetic diversity, species identification, and medicinal properties. In this study, comparative genomics and phylogenomics of Glycyrrhiza were analyzed based on the chloroplast genome. The chloroplast genomes of six Glycyrrhiza species were obtained using various assembly and annotation tools. The final assembled chloroplast genome sizes for the six Glycyrrhiza species ranged from 126,380 bp to 129,115 bp, with a total of 109-110 genes annotated. Comparative genomics results showed that the chloroplast genomes of Glycyrrhiza showed typically lacking inverted repeat regions, and the genome length, structure, GC content, codon usage, and gene distribution were highly similar. Bioinformatics analysis revealed the presence of 69-96 simple sequence repeats and 61-138 long repeats in the chloroplast genomes. Combining the results of mVISTA and nucleotide diversity, four highly variable regions were screened for species identification and relationship studies. Selection pressure analysis indicated overall purifying selection in the chloroplast genomes of Glycyrrhiza, with a few positively selected genes potentially linked to environmental adaptation. Phylogenetic analyses involving all tribes of Fabaceae with published chloroplast genomes elucidated the evolutionary relationships, and divergence time estimation estimated the chronological order of species differentiations within the Fabaceae family. The results of phylogenetic analysis indicated that species from the six subfamilies formed distinct clusters, consistent with the classification scheme of the six subfamilies. In addition, the inverted repeat-lacking clade in the subfamily Papilionoideae clustered together, and it was the last to differentiate. Co-linear analysis confirmed the conserved nature of Glycyrrhiza chloroplast genomes, and instances of gene rearrangements and inversions were observed in the subfamily Papilionoideae.

Analysis of environmental selection pressure of superoxide dismutase in deep-sea sea cucumber

Analysis of environmental selection pressure of superoxide dismutase in deep-sea sea cucumber

Transcriptomic divergence of the Rheum palmatum complex derived from top-geoherb and non-geoherb areas provides the insights into geoherbalism properties of rhubarb

Geoherb usually represents high-quality medicinal herbs with better clinical therapeutic effects, and elucidating the geoherbalism is essential for the quality improvement of traditional Chinese Medicine. However, few researches were conducted to clarify the geoherbalism based on a large scale of transcriptomics. In the present study, we compared the transcriptomes of Rheum palmatum complex derived from top-geoherb and non-geoherb areas to show the geoherbalism properties of rhubarb. A total of 412.32 Gb clean reads were obtained with unigene numbers of 100,615 after assembly. Based on the obtained transcriptome datasets, key enzyme-encoding genes involved in the anthraquinones biosynthesis were also obtained. We also found that 21 anthraquinone-related unigenes were differentially expressed between two different groups, and some of these DEGs were correlated to the content accumulation of five free anthraquinones, indicating that the gene expression profiles may promote the geoherbalism formation of rhubarb. In addition, the selective pressure analyses indicated that most paired orthologous genes between these two groups were subject to negative selection, and only a low proportion of orthologs under positive selection were detected. Functional annotation analyses indicated that these positive-selected genes related to the functions such as gene expression, substance transport, stress response and metabolism, indicating that discrepant environment also enhanced the formation of geoherbalism. Our study not only provided insights for the genetic mechanism of geoherbalism of rhubarb, but also laid more genetic cues for the future rhubarb germplasms improvement and utilization.

Integrative Analysis of the DUF668 Gene Family in Nicotiana tabacum to Excavate Their Potential Roles in Abiotic Stress Responses

The domain of the unknown function (DUF) gene families assumes pivotal roles in plant metabolic and stress responses. However, our comprehension of the tobacco DUF668 (NtDUF668) gene family and its specific reactions to heavy metal, drought, and salt stresses remain circumscribed. In the current investigation, a comprehensive genome-wide analysis of the NtDUF668 gene family was undertaken utilizing bioinformatics tools. The results unveiled a total of 20 members in the NtDUF668 gene family, denominated NtDUF668-01 to NtDUF668-20. Phylogenetic analyses indicated a closer genetic relationship of DUF668 genes between Nicotiana tabacum and Ipomoea batatas. The examination of gene structure and conservative motifs revealed a bifurcation into two major Clades, aligning with previous studies on DUF668 gene families from various plant species, emphasizing its highly conserved evolutionary mechanism across plants. The exploration of promoter regions of NtDUF668 genes revealed a plethora of cis-acting elements associated with abiotic and biotic stresses, light signaling, and phytohormones. Gene duplication events and selection pressure analysis disclosed the segmental duplication and strong purifying selection pressure during the evolution of NtDUF668 genes. Syntenic analysis indicated a relatively conserved evolutionary mechanism of DUF668 gene families within dicotyledons. Tissue-specific expression analysis suggested that NtDUF668 family members are potentially involved in root development, floral organ formation, and abscission. The expression patterns and qRT–PCR analysis of NtDUF668 genes implied the potentially functional involvements of NtDUF668s in response to multiple abiotic stresses. Furthermore, the stress-triggered member NtDUF668-08 exhibited specific nuclear localization. In conclusion, this genome-wide analysis illuminates the composition, phylogenetic relationships, and potential roles of the NtDUF668 gene family in abiotic stress responses. The identified candidate genes, particularly NtDUF668-08, warrant further research for functional investigation.