Search

Background: 'Hongyang' kiwifruit (Actinidia chinensis cv 'Hongyang') is a high-quality variety of A. chinensis with the advantages of high yield, early ripening, and high stress tolerance. Studies have confirmed that the Shaker K+ genes family is involved in plant uptake and distribution of potassium (K+). Results: Twenty-eight Shaker genes were identified and analyzed from the 'Hongyang' kiwifruit (A. chinensis cv 'Hongyang') genome. Subcellular localization results showed that more than one-third of the AcShaker genes were on the cell membrane. Phylogenetic analysis indicated that the AcShaker genes were divided into six subfamilies (I-VI). Conservative model, gene structure, and structural domain analyses showed that AcShaker genes of the same subfamily have similar sequence features and structure. The promoter cis-elements of the AcShaker genes were classified into hormone-associated cis-elements and environmentally stress-associated cis-elements. The results of chromosomal localization and duplicated gene analysis demonstrated that AcShaker genes were distributed on 18 chromosomes, and segmental duplication was the prime mode of gene duplication in the AcShaker family. GO enrichment analysis manifested that the ion-conducting pathway of the AcShaker family plays a crucial role in regulating plant growth and development and adversity stress. Compared with the transcriptome data of the control group, all AcShaker genes were expressed under low-K+stress, and the expression differences of three genes (AcShaker15, AcShaker17, and AcShaker22) were highly significant. The qRT-PCR results showed a high correlation with the transcriptome data, which indicated that these three differentially expressed genes could regulate low-K+ stress and reduce K+ damage in kiwifruit plants, thus improving the resistance to low-K+ stress. Comparison between the salt stress and control transcriptomic data revealed that the expression of AcShaker11 and AcShaker18 genes was significantly different and lower under salt stress, indicating that both genes could be involved in salt stress resistance in kiwifruit. Conclusion: The results showed that 28 AcShaker genes were identified and all expressed under low K+ stress, among which AcShaker22 was differentially and significantly upregulated. The AcShaker22 gene can be used as a candidate gene to cultivate new varieties of kiwifruit resistant to low K+ and provide a reference for exploring more properties and functions of the AcShaker genes. [ABSTRACT FROM AUTHOR]

Background: Potassium (K) is an essential nutrient for plant growth and development. Maize (Zea mays) is a widely planted crops in the world and requires a huge amount of K fertilizer. Arbuscular mycorrhizal fungi (AMF) are closely related to the K uptake of maize. Genetic improvement of maize K utilization efficiency will require elucidating the molecular mechanisms of maize K uptake through the mycorrhizal pathway. Here, we employed transcriptome and gene family analysis to elucidate the mechanism influencing the K uptake and utilization efficiency of mycorrhizal maize. Methods and results: The transcriptomes of maize were studied with and without AMF inoculation and under different K conditions. AM symbiosis increased the K concentration and dry weight of maize plants. RNA sequencing revealed that genes associated with the activity of the apoplast and nutrient reservoir were significantly enriched in mycorrhizal roots under low-K conditions but not under high-K conditions. Weighted gene correlation network analysis revealed that three modules were strongly correlated with K content. Twenty-one hub genes enriched in pathways associated with glycerophospholipid metabolism, glycerolipid metabolism, starch and sucrose metabolism, and anthocyanin biosynthesis were further identified. In general, these hub genes were upregulated in AMF-colonized roots under low-K conditions. Additionally, the members of 14 gene families associated with K obtain were identified (ARF: 38, ILK: 4, RBOH: 12, RUPO: 20, MAPKK: 89, CBL: 14, CIPK: 44, CPK: 40, PIN: 10, MYB: 174, NPF: 79, KT: 19, HAK/HKT/KUP: 38, and CPA: 8) from maize. The transcript levels of these genes showed that 92 genes (ARF:6, CBL:5, CIPK:13, CPK:2, HAK/HKT/KUP:7, PIN:2, MYB:26, NPF:16, RBOH:1, MAPKK:12 and RUPO:2) were upregulated with AM symbiosis under low-K conditions. Conclusions: This study indicated that AMF increase the resistance of maize to low-K stress by regulating K uptake at the gene transcription level. Our findings provide a genome-level resource for the functional assignment of genes regulated by K treatment and AM symbiosis in K uptake-related gene families in maize. This may contribute to elucidate the molecular mechanisms of maize response to low K stress with AMF inoculation, and provided a theoretical basis for AMF application in the crop field. [ABSTRACT FROM AUTHOR]

Background: Invertases (INVs) are key enzymes in sugar metabolism, cleaving sucrose into glucose and fructose and playing an important role in plant development and the stress response, however, the INV gene family in passion fruit has not been systematically reported. Results: In this study, a total of 16 PeINV genes were identified from the passion fruit genome and named according to their subcellular location and chromosome position. These include six cell wall invertase (CWINV) genes, two vacuolar invertase (VINV) genes, and eight neutral/alkaline invertase (N/AINV) genes. The gene structures, phylogenetic tree, and cis-acting elements of PeINV gene family were predicted using bioinformatics methods. Results showed that the upstream promoter region of the PeINV genes contained various response elements; particularly, PeVINV2, PeN/AINV3, PeN/AINV5, PeN/AINV6, PeN/AINV7, and PeN/AINV8 had more response elements. Additionally, the expression profiles of PeINV genes under different abiotic stresses (drought, salt, cold temperature, and high temperature) indicated that PeCWINV5, PeCWINV6, PeVINV1, PeVINV2, PeN/AINV2, PeN/AINV3, PeN/AINV6, and PeN/AINV7 responded significantly to these abiotic stresses, which was consistent with cis-acting element prediction results. Sucrose, glucose, and fructose are main soluble components in passion fruit pulp. The contents of total soluble sugar, hexoses, and sweetness index increased significantly at early stages during fruit ripening. Transcriptome data showed that with an increase in fruit development and maturity, the expression levels of PeCWINV2, PeCWINV5, and PeN/AINV3 exhibited an up-regulated trend, especially for PeCWINV5 which showed highest abundance, this correlated with the accumulation of soluble sugar and sweetness index. Transient overexpression results demonstrated that the contents of fructose, glucose and sucrose increased in the pulp of PeCWINV5 overexpressing fruit. It is speculated that this cell wall invertase gene, PeCWINV5, may play an important role in sucrose unloading and hexose accumulation. Conclusion: In this study, we systematically identified INV genes in passion fruit for the first time and further investigated their physicochemical properties, evolution, and expression patterns. Furthermore, we screened out a key candidate gene involved in hexose accumulation. This study lays a foundation for further study on INV genes and will be beneficial on the genetic improvement of passion fruit breeding. [ABSTRACT FROM AUTHOR]

Background: Alfalfa (Medicago sativa L.) is an essential leguminous forage with high nutrition and strong adaptability. The TIFY family is a plant-specific transcription factor identified in many plants. However, few reports have been reported on the phylogenetic analysis and gene expression profiling of TIFY family genes in alfalfa. Result: A total of 84 TIFY genes belonging to 4 categories were identified in alfalfa, including 58 MsJAZs, 18 MsZMLs, 4 MsTIFYs and 4 MsPPDs, respectively. qRT-PCR data from 8 genes in different tissues revealed that most MsTIFY genes were highly expressed in roots. The expression of MsTIFY14 was up-regulated after different times in both thrips-resistant and susceptible alfalfa after thrips feeding, and the expression of the remaining MsTIFYs had a strong correlation with the time of thrips feeding. Different abiotic stresses, including drought, salt, and cold, could induce or inhibit the expression of MsTIFY genes to varying degrees. In addition, the eight genes were all significantly up-regulated by JA and/or SA. Interestingly, MsTIFY77 was induced considerably by all the biotic, abiotic, or plant hormones (JA or SA) except ABA. Conclusion: Our study identified members of the TIFY gene family in alfalfa and analyzed their structures and possible functions. It laid the foundation for further research on the molecular functions of TIFYs in alfalfa. [ABSTRACT FROM AUTHOR]

Autism Spectrum Disorder (ASD) presents a significant challenge due to its complex genetic basis and associated comorbidities. Among the genes implicated in ASD, SHANK3 has been identified as a critical player, affecting synaptic structure and function. This review examines the role of SHANK3 in ASD, highlighting the genetic diversity and the systemic nature of the disorder. Utilizing animal models, studies have uncovered autism-like behaviours and synaptic dysfunctions linked to SHANK3 deficiency, suggesting potential therapeutic targets. Furthermore, the review delves into the specific gene families associated with ASD, emphasizing the dynamic regulation between translation and transcription processes and the impact of mutations on synaptic translation and proteins. Molecular changes in SHANK3-deficient animal models reveal alterations in protein composition, localization, and transcription, particularly affecting the striatum and involving essential proteins and signalling pathways. Therapeutic strategies, including pharmaceutical compounds and genetic restoration, show promise in addressing the neuropsychiatric symptoms and physiological abnormalities observed in SHANK3-deficient mice. This research not only advances our understanding of ASD's neurobiological basis but also underscores the potential of targeted interventions to mitigate symptoms and improve the quality of life for individuals affected by ASD and related disorders. [ABSTRACT FROM AUTHOR]

Background: The KT/HAK/KUP is the largest K+ transporter family in plants, playing crucial roles in K+ absorption, transport, and defense against environmental stress. Sweet watermelon is an economically significant horticultural crop belonging to the genus Citrullus, with a high demand for K+ during its growth process. However, a comprehensive analysis of the KT/HAK/KUP gene family in watermelon has not been reported. Results: 14 KT/HAK/KUP genes were identified in the genomes of each of seven Citrullus species. These KT/HAK/KUPs in watermelon were unevenly distributed across seven chromosomes. Segmental duplication is the primary driving force behind the expansion of the KT/HAK/KUP family, subjected to purifying selection during domestication (Ka/Ks < 1), and all KT/HAK/KUPs exhibit conserved motifs and could be phylogenetically classified into four groups. The promoters of KT/HAK/KUPs contain numerous cis-regulatory elements related to plant growth and development, phytohormone response, and stress response. Under K+ deficiency, the growth of watermelon seedlings was significantly inhibited, with cultivated watermelon experiencing greater impacts (canopy width, redox enzyme activity) compared to the wild type. All KT/HAK/KUPs in C. lanatus and C. amarus exhibit specific expression responses to K+-deficiency and drought stress by qRT-PCR. Notably, ClG42_07g0120700/CaPI482276_07g014010 were predominantly expressed in roots and were further induced by K+-deficiency and drought stress. Additionally, the K+ transport capacity of ClG42_07g0120700 under low K+ stress was confirmed by yeast functional complementation assay. Conclusions: KT/HAK/KUP genes in watermelon were systematically identified and analyzed at the pangenome level and provide a foundation for understanding the classification and functions of the KT/HAK/KUPs in watermelon plants. [ABSTRACT FROM AUTHOR]

Background: Phytochrome-interacting factors (PIFs) plays an important role in plants as hubs for intracellular signaling regulation. The PIF gene family has been identified and characterized in many plants, but alfalfa (Medicago sativa L.), an important perennial high-quality legume forage, has not been reported on the PIF gene family. Results: In this study, we presented the identification and characterization of five MsPIF genes in alfalfa (Medicago sativa L.). Phylogenetic analysis indicated that PIFs from alfalfa and other four plant species could be divided into three groups supported by similar motif analysis. The collinearity analysis of the MsPIF gene family showed the presence of two gene pairs, and the collinearity analysis with AtPIFs showed three gene pairs, indicating that the evolutionary process of this family is relatively conservative. Analysis of cis-acting elements in promoter regions of MsPIF genes indicated that various elements were related to light, abiotic stress, and plant hormone responsiveness. Gene expression analyses demonstrated that MsPIFs were primarily expressed in the leaves and were induced by various abiotic stresses. Conclusion: This study conducted genome-wide identification, evolution, synteny analysis, and expression analysis of the PIFs in alfalfa. Our study lays a foundation for the study of the biological functions of the PIF gene family and provides a useful reference for improving abiotic stress resistance in alfalfa. [ABSTRACT FROM AUTHOR]

ABSTRACT Pentatricopeptide repeat (PPR) gene family constitutes one of the largest gene families in plants, which mainly participate in RNA editing and RNA splicing of organellar RNAs, thereby affecting the organellar development. Recently, some evidence elucidated the important roles of PPR proteins in the albino process of plant leaves. However, the functions of PPR genes in the woody mangrove species have not been investigated. In this study, using a typical true mangrove Kandelia obovata, we systematically identified 298 PPR genes and characterized their general features and physicochemical properties, including evolutionary relationships, the subcellular localization, PPR motif type, the number of introns and PPR motifs, and isoelectric point, and so forth. Furthermore, we combined genome‐wide association studies (GWAS) and transcriptome analysis to identify the genetic architecture and potential PPR genes associated with propagule leaves colour variations of K. obovata. As a result, we prioritized 16 PPR genes related to the albino phenotype using different strategies, including differentially expressed genes analysis and genetic diversity analysis. Further analysis discovered two genes of interest, namely Maker00002998 (PLS‐type) and Maker00003187 (P‐type), which were differentially expressed genes and causal genes detected by GWAS analysis. Moreover, we successfully predicted downstream target chloroplast genes (rps14, rpoC1 and rpoC2) bound by Maker00002998 PPR proteins. The experimental verification of RNA editing sites of rps14, rpoC1, and rpoC2 in our previous study and the verification of interaction between Maker00002998 and rps14 transcript using in vitro RNA pull‐down assays revealed that Maker00002998 PPR protein might be involved in the post‐transcriptional process of chloroplast genes. Our result provides new insights into the roles of PPR genes in the albinism mechanism of K. obovata propagule leaves. [ABSTRACT FROM AUTHOR]

Among isopods (order Isopoda), species belonging to the suborder Oniscidea are the most successful colonisers of terrestrial habitats and the largest family in this order is Armadillidae. In this study, the complete mitochondrial genome of an armadillid isopod, Cubaris murina, was assembled and characterised in detail for the first time. The gene order and secondary structure of transfer RNA (tRNA) genes were compared among isopods. The 14,212-bp-long mitogenome consists of 13 protein-coding genes, 20 tRNAs, two ribosomal RNAs, and four non-coding regions (NCRs), each longer than 100 bp. Two tRNA genes (trnG and trnI) were missing and several other tRNA genes exhibited secondary structures that departed from the typical cloverleaf pattern. Deletions and nucleotide insertions in the T- and D-arms and mismatches in the acceptor stem were commonly observed in these tRNAs. Notably, the mitogenome of C. murina exhibited a new gene order containing four tRNA rearrangements, and differences in the location of the NCRs were detected when compared to those reported for other isopods, including oniscid species. Interestingly, the second NCR had a similar number of secondary structures and repetitive elements to those reported in the control region of other crustaceans. The novel gene order, tRNA secondary structures, and NCR organisation reported herein for the mitogenome of C. murina confirms the hyper-variability of mitochondrial genomes in isopods. This study is the first to characterise in detail the mitogenome of the invasive isopod C. murina, belonging to the species-rich family Armadillidae within the Isopoda, and advances our understanding of mitochondrial genome architecture in isopods. [ABSTRACT FROM AUTHOR]

Summary: Spartina alterniflora is an exo‐recretohalophyte Poaceae species that is able to grow well in seashore, but the genomic basis underlying its adaptation to salt tolerance remains unknown. Here, we report a high‐quality, chromosome‐level genome assembly of S. alterniflora constructed through PacBio HiFi sequencing, combined with high‐throughput chromosome conformation capture (Hi‐C) technology and Illumina‐based transcriptomic analyses. The final 1.58 Gb genome assembly has a contig N50 size of 46.74 Mb. Phylogenetic analysis suggests that S. alterniflora diverged from Zoysia japonica approximately 21.72 million years ago (MYA). Moreover, whole‐genome duplication (WGD) events in S. alterniflora appear to have expanded gene families and transcription factors relevant to salt tolerance and adaptation to saline environments. Comparative genomics analyses identified numerous species‐specific genes, significantly expanded genes and positively selected genes that are enriched for 'ion transport' and 'response to salt stress'. RNA‐seq analysis identified several ion transporter genes including the high‐affinity K+ transporters (HKTs), SaHKT1;2, SaHKT1;3 and SaHKT1;8, and high copy number of Salt Overly Sensitive (SOS) up‐regulated under high salt conditions, and the overexpression of SaHKT2;4 in Arabidopsis thaliana conferred salt tolerance to the plant, suggesting specialized roles for S. alterniflora to adapt to saline environments. Integrated metabolomics and transcriptomics analyses revealed that salt stress activate glutathione metabolism, with differential expressions of several genes such as γ‐ECS, GSH‐S, GPX, GST and PCS in the glutathione metabolism. This study suggests several adaptive mechanisms that could contribute our understanding of evolutional basis of the halophyte. [ABSTRACT FROM AUTHOR]

A recent study published in Integrative Zoology presents findings on the chromosome-level genome assembly of the tree sparrow. The researchers identified a burst of new genes with different structures but a similar testis-specific expression pattern in the tree sparrow. These new genes were found to be the result of segmental duplications (SDs), which are large and highly identical chromosomal segments. The study also suggests that the testis-specific expression pattern of these new genes supports the "out of testis" hypothesis, which proposes that new genes emerge and acquire new functions in the testis before being expressed in other tissues. This research provides valuable insights into the role of SDs in genome evolution and the adaptation of tree sparrows to different environments. [Extracted from the article]

Blister beetles (Coleoptera: Meloidae) produce a natural defensive toxin cantharidin (CTD), which has been used for various cancer treatments and other diseases. Currently, the lack of chromosome‐level reference genomes in Meloidae limits further understanding of the mechanism of CTD biosynthesis and environmental adaptation. In this study, the chromosome‐level genome assembly of Mylabris phalerata was generated based on PacBio and Hi‐C sequencing. This reference genome was about 136.68 Mb in size with contig N50 of 9.17 Mb and composed of 12 chromosomes. In comparison to six other Coleoptera insects, M. phalerata exhibited multiple expanded gene families enriched in juvenile hormone (JH) biosynthetic process pathway, farnesol dehydrogenase activity, and cytochrome P450, which may be related to CTD biosynthesis. Consistently, the transcriptomic analysis suggested the "terpenoid backbone biosynthesis" pathway and "the juvenile hormone" as putative core pathways of CTD biosynthesis and presented eight up‐regulated differential expression genes in male adults as candidate genes. It is possible that the restricted feeding niche and lifestyle of M. phalerata were the cause of the gene family's contraction of odorant binding proteins. The ABC transporters (ABCs) related to exporting bound toxins out of the cell and the resistance to the self‐secreted toxins (e.g. CTD) were also contracted, possibly due to other self‐protection strategies in M. phalerata. A foundation of understanding CTD biosynthesis and environmental adaptation of blister beetles will be established by our reference genome and discoveries. [ABSTRACT FROM AUTHOR]

Background and Aims Carex breviculmis is a perennial herb with good resistance and is widely used for forage production and turf management. It is important in ecology, environmental protection and biodiversity conservation, but faces several challenges due to human activities. However, the absence of genome sequences has limited basic research and the improvement of wild plants. Methods We annotated the genome of C. breviculmis and conducted a systematic analysis to explore its resistance to harsh environments. We also conducted a comparative analysis of Achnatherum splendens , which is similarly tolerant to harsh environments. Key Results The assembled the genome comprises 469.01 Mb, revealing 37 372 genes with a BUSCO completeness score of 99.0 %. The genome has 52.03 % repetitive sequences, primarily influenced by recent LTR insertions that have contributed to its expansion. Phylogenetic analysis suggested that C. breviculmis diverged from C. littledalei ~6.61 million years ago. Investigation of repetitive sequences and expanded gene families highlighted a rapid expansion of tandem duplicate genes, particularly in areas related to sugar metabolism, synthesis of various amino acids, and phenylpropanoid biosynthesis. Additionally, our analysis identified crucial genes involved in secondary metabolic pathways, such as glycolysis, phenylpropanoid biosynthesis and amino acid metabolism, which have undergone positive selection. We reconstructed the sucrose metabolic pathway and identified significant gene expansions, including 16 invertase, 9 sucrose phosphate synthase and 12 sucrose synthase genes associated with sucrose metabolism, which showed varying levels of expansion. Conclusions The expansion of these genes, coupled with subsequent positive selection, contributed to the ability of C. breviculmis to adapt to environmental stressors. This study lays the foundation for future research on the evolution of Carex plants, their environmental adaptations, and potential genetic breeding. [ABSTRACT FROM AUTHOR]

Background: Nuclear factor Y (NF-Y) plays a vital role in numerous biological processes as well as responses to biotic and abiotic stresses. However, its function in ginger (Zingiber officinale Roscoe), a significant medicinal and dietary vegetable, remains largely unexplored. Although the NF-Y family has been thoroughly identified in many plant species, and the function of individual NF-Y TFs has been characterized, there is a paucity of knowledge concerning this family in ginger. Methods: We identified the largest number of NF-Y genes in the ginger genome using two BLASTP methods as part of our ginger genome research project. The conserved motifs of NF-Y proteins were analyzed through this process. To examine gene duplication events, we employed the Multiple Collinearity Scan toolkit (MCScanX). Syntenic relationships of NF-Y genes were mapped using the Dual Synteny Plotter software. Multiple sequence alignments were performed with MUSCLE under default parameters, and the resulting alignments were used to generate a maximum likelihood (ML) phylogenetic tree with the MEGA X program. RNA-seq analysis was conducted on collected samples, and statistical analyses were performed using Sigma Plot v14.0 (SYSTAT Software, USA). Results: In this study, the ginger genome was utilized to identify 36 NF-Y genes (10 ZoNF-YAs, 16 ZoNF-YBs, and 10 ZoNF-YCs), which were renamed based on their chromosomal distribution. Ten distinct motifs were identified within the ZoNF-Y genes, with certain unique motifs being vital for gene function. By analyzing their chromosomal location, gene structure, conserved protein motifs, and gene duplication events, we gained a deeper understanding of the evolutionary characteristics of these ZoNF-Y genes. Detailed analysis of ZoNF-Y gene expression patterns across various tissues, performed through RNA-seq and qRT-PCR, revealed their significant role in regulating ginger rhizome and flower growth and development. Additionally, we identified the ZoNF-Y family genes that responded to abiotic stresses. Conclusion: This study represents the first identification of the ZoNF-Y family in ginger. Our findings contribute to research on evolutionary characteristics and provide a better understanding of the molecular basis for development and abiotic stress response. Furthermore, it lays the foundation for further functional characterization of ZoNF-Y genes with an aim of ginger crop improvement. [ABSTRACT FROM AUTHOR]

Wax biosynthesis is closely controlled by many regulators under different environmental conditions. We have previously shown that the module miR156–SQUAMOSA PROMOTER BINDING PROTEIN-LIKE9 (SPL9)–DEWAX is involved in the diurnal regulation of wax production; however, it was not determined whether other SPLs are also involved in wax synthesis. Here, we report that SPL13 also regulates drought-induced wax production, by directly and indirectly affecting the expression of the two wax biosynthesis genes ECERIFERUM1 (CER1) and CER4 , respectively. In addition, we show that SPL13 together with SPL9 redundantly regulates wax accumulation under both normal and drought stress conditions, and that simultaneous mutation of both genes additively increases cuticle permeability and decreases drought tolerance. However, in contrast to SPL9, SPL13 does not seem to participate in the DEWAX-mediated diurnal regulation of wax production. [ABSTRACT FROM AUTHOR]

Furthering our knowledge of the skin microbiome is essential to understand health and disease in canines. To date, studies into the canine skin microbiome have focused on 16S rRNA high throughput sequencing however, these lack the granularity of species and strain level taxonomic characterisation and their associated functions. The aim of this study was to provide a comprehensive assessment of the skin microbiome by analysing the skin microbiome of 72 healthy adult colony dogs, across four distinct skin sites and four breeds, using metagenomic sequencing. Our analysis revealed that breed and skin site are drivers of variation, and a core group of taxa and genes are present within the skin microbiome of healthy dogs, comprising 230 taxa and 1219 gene families. We identified 15 species within the core microbiome that are represented by more than one strain. The biosynthesis of secondary metabolites pathway was enriched in the core microbiome suggesting the skin microbiome may play a role in colonisation resistance and protection from invading pathogens. Additionally, we uncovered the novelty of the canine skin microbiome and show that further investigation is required to increase the suitability of current databases for metagenomic sequencing of canine skin samples. [ABSTRACT FROM AUTHOR]

The SET domain genes (SDGs) are significant contributors to various aspects of plant growth and development, mainly includes flowering, pollen development, root growth, regulation of the biological clock and branching patterns. To clarify the biological functions of the chrysanthemum SDG family, the SDG family members of four chrysanthemum cultivars and three related wild species were identified; their physical and chemical properties, protein domains and conserved motifs were predicted and analyzed. The results showed that 59, 67, 67, 102, 106, 114, and 123 SDGs were identified from Chrysanthemum nankingense, Chrysanthemum lavandulifolium, Chrysanthemum seticuspe, Chrysanthemum × morifolium cv. 'Hechengxinghuo', 'Zhongshanzigui', 'Quanxiangshuichang' and 'Jinbeidahong', respectively. The SDGs were divided into 5–7 subfamilies by cluster analysis; different conserved motifs were observed in particular families. The SDGs of C. lavandulifolium and C. seticuspe were distributed unevenly on 9 chromosomes. SDG promoters of different species include growth and development, photo-response, stress response and hormone responsive elements, among them, the cis-acting elements related to MeJA response had the largest proportion. The expression of chrysanthemum SDG genes was observed for most variable selected genes which has close association with important Arabidopsis thaliana genes related to flowering regulation. The qPCR results showed that the expression trend of SDG genes varied in different tissues at different growth stages with high expression in the flowering period. The ClSDG29 showed higher expression in the flower and bud tissues, which indicate that ClSDG29 might be associated with flowering regulation in chrysanthemum. In summary, the results of this study can provide a basis for subsequent research on chrysanthemum flowering time regulation. [ABSTRACT FROM AUTHOR]

Background: The transposons of the hAT superfamily are the most widespread transposons ever known. SLEEPER genes encode domesticated transposases from the hAT superfamily, which may have lost their transposable functions during long-term evolution and transformed into host proteins that regulate plant growth and development. Results: This study identified 162 members of the SLEEPER gene family from Brassica napus. These members are widely distributed on 19 chromosomes, mainly in the Cn subgenome, and have promoters with various cis-acting elements related to hormone regulation, abiotic stress, and growth and development regulation. Most of the genes in this family contain similar conserved domains and motifs, and the closer the genes are distributed on evolutionary branches, the more similar their structures are. Transcriptome sequencing performed on tissues at different growth stages from B. napus line 3529 indicated that these genes had different expression patterns, and nearly half of the genes were not detectably expressed in all samples. Conclusions: This study investigated the gene structure, expression patterns, evolutionary features, and gene localization of the SLEEPER family members to confirm the significance of these genes in the growth of B. napus, providing a reference for the study of transposon domestication and outstanding genetic resources for the genetic improvement of B. napus. [ABSTRACT FROM AUTHOR]

Background: Salt Overly Sensitive 1 (SOS1), a plasma membrane Na+/H+ exchanger, is essential for plant salt tolerance. Salt damage is a significant abiotic stress that impacts plant species globally. All living organisms require copper (Cu), a necessary micronutrient and a protein cofactor for many biological and physiological processes. High Cu concentrations, however, may result in pollution that inhibits the growth and development of plants. The function and production of mangrove ecosystems are significantly impacted by rising salinity and copper contamination. Results: A genome-wide analysis and bioinformatics techniques were used in this study to identify 20 SOS1 genes in the genome of Kandelia obovata. Most of the SOS1 genes were found on the plasma membrane and dispersed over 11 of the 18 chromosomes. Based on phylogenetic analysis, KoSOS1s can be categorized into four groups, similar to Solanum tuberosum. Kandelia obovata's SOS1 gene family expanded due to tandem and segmental duplication. These SOS1 homologs shared similar protein structures, according to the results of the conserved motif analysis. The coding regions of 20 KoSOS1 genes consist of amino acids ranging from 466 to 1221, while the exons include amino acids ranging from 3 to 23. In addition, we found that the 2.0 kb upstream promoter region of the KoSOS1s gene contains several cis-elements associated with phytohormones and stress responses. According to the expression experiments, seven randomly chosen genes experienced up- and down-regulation of their expression levels in response to copper (CuCl2) and salt stressors. Conclusions: For the first time, this work systematically identified SOS1 genes in Kandelia obovata. Our investigations also encompassed physicochemical properties, evolution, and expression patterns, thereby furnishing a theoretical framework for subsequent research endeavours aimed at functionally characterizing the Kandelia obovata SOS1 genes throughout the life cycle of plants. [ABSTRACT FROM AUTHOR]

Background: The DELLA proteins, a class of GA signaling repressors, belong to the GRAS family of plant-specific nuclear proteins. Members of DELLA gene family encode transcriptional regulators with diverse functions in plant development and abiotic stress responses. To date, DELLAs have been identified in various plant species, such as Arabidopsis thaliana, Malus domestica, Populus trichocarpa, and other land plants. Most information of DELLA family genes was obtained from A. thaliana, whereas little is known about the DELLA gene family in blueberry. Results: In this study, we identified three DELLA genes in blueberry (Vaccinium darrowii, VdDELLA) and provided a complete overview of VdDELLA gene family, describing chromosome localization, protein properties, conserved domain, motif organization, and phylogenetic analysis. Three VdDELLA members, containing two highly conserved DELLA domain and GRAS domain, were distributed across three chromosomes. Additionally, cis-acting elements analysis indicated that VdDELLA genes might play a critical role in blueberry developmental processes, hormone, and stress responses. Expression analysis using quantitative real-time PCR (qRT-PCR) revealed that all of three VdDELLA genes were differentially expressed across various tissues. VdDELLA2 was the most highly expressed VdDELLA in all denoted tissues, with a highest expression in mature fruits. In addition, all of the three VdDELLA genes actively responded to diverse abiotic stresses. Based on qRT-PCR analysis, VdDELLA2 might act as a key regulator in V. darrowii in response to salt stress, whereas VdDELLA1 and VdDELLA2 might play an essential role in cold stress response. Under drought stress, all of three VdDELLA genes were involved in mediating drought response. Furthermore, their transiently co-localization with nuclear markers in A. thaliana protoplasts demonstrated their transcriptional regulator roles. Conclusions: In this study, three VdDELLA genes were identified in V. darrowii genome. Three VdDELLA genes were closely related to the C. moschata DELLA genes, S. lycopersicum DELLA genes, and M. domestica DELLA genes, respectively, indicating their similar biological functions. Expression analysis indicated that VdDELLA genes were highly efficient in blueberry fruit development. Expression patterns under different stress conditions revealed the differentially expressed VdDELLA genes responding to salt, drought, and cold stress. Overall, these results enrich our understanding of evolutionary relationship and potential functions of VdDELLA genes, which provide valuable information for further studies on genetic improvement of the plant yield and plant resistance. [ABSTRACT FROM AUTHOR]

Background: CINNAMYL ALCOHOL DEHYDROGENASE (CAD) catalyzes the NADPH-dependent reduction of cinnamaldehydes into cinnamyl alcohols and is a key enzyme found at the final step of the monolignol pathway. Cinnamyl alcohols and their conjugates are subsequently polymerized in the secondary cell wall to form lignin. CAD genes are typically encoded by multi-gene families and thus traditionally organized into general classifications of functional relevance. Results: In silico analysis of the hexaploid Triticum aestivum genome revealed 47 high confidence TaCAD copies, of which three were determined to be the most significant isoforms (class I) considered bone fide CADs. Class I CADs were expressed throughout development both in RNAseq data sets as well as via qRT-PCR analysis. Of the 37 class II TaCADs identified, two groups were observed to be significantly co-expressed with class I TaCADs in developing tissue and under chitin elicitation in RNAseq data sets. These co-expressed class II TaCADs were also found to be phylogenetically unrelated to a separate clade of class II TaCADs previously reported to be an influential resistance factor to pathogenic fungal infection. Lastly, two groups were phylogenetically identified as class III TaCADs, which possess distinct conserved gene structures. However, the lack of data supporting their catalytic activity for cinnamaldehydes and their bereft transcriptional presence in lignifying tissues challenges their designation and function as CADs. Conclusions: Taken together, our comprehensive transcriptomic analyses suggest that TaCAD genes contribute to overlapping but nonredundant functions during T. aestivum growth and development across a wide variety of agroecosystems and provide tolerance to various stressors. [ABSTRACT FROM AUTHOR]

Background: Hermansky-Pudlak Syndrome (HPS), a rare autosomal recessive disorder, is characterized by oculocutaneous albinism, bleeding diathesis, and sometimes severe lung problems and inflammatory bowel disease. Symptoms include skin and hair pigmentation variations, along with visual impairments. Variants in eleven genes encoding protein complexes essential for membrane trafficking and intracellular endosomal transport pathways underlie various recognized HPS subtypes. This study focuses on HPS-9, a subtype of Hermansky-Pudlak Syndrome caused by a variant in the BLOC1S6 gene, which is a subunit of the BLOC1 complex. In this study, a novel Copy Number Variation (CNV) in the aforementioned gene in an Iranian family is reported. The study aims to better understand the etiology of HPS-9 symptoms by identifying and confirming the variant and determining whether the gene is expressed despite the deletion. There have only been five reports of this syndrome in the literature thus far. Our novel CNV represents a significant contribution to understanding the genetic basis of HPS-9. Results: This study investigates a male patient presenting with albinism. Whole Exome Sequencing (WES) identified a homozygous deletion of approximately 350 bp using CNV analysis. The deletion affects the intronic region of the BLOC1S6 gene, causing uncertainties in defining the exact boundaries due to WES limitations. Primer walking and GAP-PCR techniques were used to define the deletion boundaries. Subsequent assessments of this variant across other family members helped identify homozygous affected members and heterozygous carriers. The absence of BLOC1S6 expression in the affected individual was confirmed through Real-time PCR experiments. These findings underscore the importance of understanding the implications for the patient's healthcare and potential therapeutic approaches. Conclusion: This study introduces a case of Hermansky-Pudlak Syndrome Type 9 (HPS-9) caused by a homozygous deletion in the BLOC1S6 gene. We identified an approximately 7-kb deletion encompassing exon 1 and the intronic region of the gene. The absence of BLOC1S6 expression, confirmed via Real-time PCR, highlights the importance of studying the pathogenicity of the deletion and its impact on the patient's health. Our findings contribute to the sparse knowledge on HPS-9 and underscore the need for further exploration into the genetic causes of this rare disorder. [ABSTRACT FROM AUTHOR]

Background: The mitogen-activated protein kinase (MAPK) cascade is crucial cell signal transduction mechanism that plays an important role in plant growth and development, metabolism, and stress responses. The MAPK cascade includes three protein kinases, MAPK, MAPKK, and MAPKKK. The three protein kinases mediate signaling to downstream response molecules by sequential phosphorylation. The MAPK gene family has been identified and analyzed in many plants, however it has not been investigated in alfalfa. Results: In this study, Medicago sativa MAPK genes (referred to as MsMAPKs) were identified in the tetraploid alfalfa genome. Eighty MsMAPKs were divided into four groups, with eight in group A, 21 in group B, 21 in group C and 30 in group D. Analysis of the basic structures of the MsMAPKs revealed presence of a conserved TXY motif. Groups A, B and C contained a TEY motif, while group D contained a TDY motif. RNA-seq analysis revealed tissue-specificity of two MsMAPKs and tissue-wide expression of 35 MsMAPKs. Further analysis identified MsMAPK members responsive to drought, salt, and cold stress conditions. Two MsMAPKs (MsMAPK70 and MsMAPK75) responds to salt and cold stresses; two MsMAPKs (MsMAPK60 and MsMAPK73) responds to cold and drought stresses; four MsMAPKs (MsMAPK1, MsMAPK33, MsMAPK64 and MsMAPK71) responds to salt and drought stresses; and two MsMAPKs (MsMAPK5 and MsMAPK7) responded to all three stresses. Conclusion: This study comprehensively identified and analysed the alfalfa MAPK gene family. Candidate genes related to abiotic stresses were screened by analysing the RNA-seq data. The results provide key information for further analysis of alfalfa MAPK gene functions and improvement of stress tolerance. [ABSTRACT FROM AUTHOR]

Formin proteins, characterized by the FH2 domain, are critical in regulating actin-driven cellular processes and cytoskeletal dynamics during abiotic stress. However, no genome-wide analysis of the formin gene family has yet to be conducted in the economically significant plant potato (Solanum tuberosum L.). In this study, 26 formin genes were identified and characterized in the potato genome (named as StFH), each containing the typical FH2 domain and distributed across the ten chromosomes. The StFH was categorized into seven subgroups (A-G) and the gene structure and motif analysis demonstrated higher structural similarities within the subgroups. Besides, the StFH exhibited ancestry and functional similarities with Arabidopsis. The Ka/Ks ratio indicated that StFH gene pairs were evolving through purifying selection, with five gene pairs exhibiting segmental duplications and two pairs exhibiting tandem duplications. Subcellular localization analysis suggested that most of the StFH genes were located in the chloroplast and plasma membrane. Moreover, 54 cis-acting regulatory elements (CAREs) were identified in the promoter regions, some of which were associated with stress responses. According to gene ontology analysis, the majority of the StFH genes were involved in biological processes, with 63 out of 74 GO terms affecting actin polymerization. Six major transcription factor families, including bZIP, C2H2, ERF, GATA, LBD, NAC, and HSF, were identified that were involved in the regulation of StFH genes in various abiotic stresses, including drought. Further, the 60 unique microRNAs targeted 24 StFH by regulating gene expression in response to drought stress were identified. The expression of StFH genes in 14 different tissues, particularly in drought-responsive tissues such as root, stem, shoot apex, and leaf, underscores their significance in managing drought stress. RNA-seq analysis of the drought-resistant Qingshu No. 9 variety revealed the potential role of up-regulated genes, including StFH2, StFH10, StFH19, and StFH25, in alleviating drought stress. Overall, these findings provide crucial insights into the response to drought stress in potatoes and can be utilized in breeding programs to develop potato cultivars with enhanced drought-tolerant traits. [ABSTRACT FROM AUTHOR]

Background: SET domain-containing histone lysine methyltransferases (HKMTs) and JmjC domain-containing histone demethylases (JHDMs) are essential for maintaining dynamic changes in histone methylation across parasite development and infection. However, information on the HKMTs and JHDMs in human pathogenic piroplasms, such as Babesia duncani and Babesia microti, and in veterinary important pathogens, including Babesia bigemina, Babesia bovis, Theileria annulata and Theileria parva, is limited. Results: A total of 38 putative KMTs and eight JHDMs were identified using a comparative genomics approach. Phylogenetic analysis revealed that the putative KMTs can be divided into eight subgroups, while the JHDMs belong to the JARID subfamily, except for BdJmjC1 (BdWA1_000016) and TpJmjC1 (Tp Muguga_02g00471) which cluster with JmjC domain only subfamily members. The motifs of SET and JmjC domains are highly conserved among piroplasm species. Interspecies collinearity analysis provided insight into the evolutionary duplication events of some SET domain and JmjC domain gene families. Moreover, relative gene expression analysis by RT‒qPCR demonstrated that the putative KMT and JHDM gene families were differentially expressed in different intraerythrocytic developmental stages of B. duncani, suggesting their role in Apicomplexa parasite development. Conclusions: Our study provides a theoretical foundation and guidance for understanding the basic characteristics of several important piroplasm KMT and JHDM families and their biological roles in parasite differentiation. [ABSTRACT FROM AUTHOR]

Phenylalanine ammonia lyase (PAL) is a widely studied enzyme in plant biology due to its role in connecting primary metabolism to secondary phenylpropanoid metabolism, significantly influencing plant growth, development, and stress response. Although PAL genes have been extensively studied in various plant species but their exploration in cucumber has been limited. This study successfully identified 11 CsPAL genes in Cucumis sativus (cucumber). These CsPAL genes were categorized based on their conserved sequences revealing patterns through MEME analysis and multiple sequence alignment. Interestingly, cis-elements related to stress were found in the promoter regions of CsPAL genes, indicating their involvement in responding to abiotic stress. Furthermore, these gene's promoters contained components associated with light, development and hormone responsiveness. This suggests that they may have roles in hormone developmental processes. MicroRNAs were identified as a key regulators for the CsPAL genes, playing a crucial role in modulating their expression. This discovery underscores the complex regulatory network involved in the plant's response to various stress conditions. The influence of these microRNAs further highlights the complicated mechanisms that plants use to manage stress. Gene expression patterns were analyzed using RNA-seq data. The significant upregulation of CsPAL9 during HT3h (heat stress for 3 h) and the heightened upregulation of both CsPAL9 and CsPAL7 under HT6h (heat stress for 6 h) in the transcriptome study suggest a potential role for these genes in cucumber's tolerance to heat stress. This comprehensive investigation aims to enhance our understanding of the PAL gene family's versatility, offering valuable insights for advancements in cucumber genetics. [ABSTRACT FROM AUTHOR]

Potato is one of the four staple food crops in the world. It has a wide range of cultivation, high yield, and high nutritional value. Enhancing the photosynthesis of potato is particularly important as it leads to an increase in the potato yield. The light-harvesting pigment-binding protein complex is very important for plant photosynthesis. We identified 12 Stlhcb gene family members from the potato variety "Atlantic" using transcriptome sequencing and bioinformatics. The proteins encoded by the Stlhcb gene family have between 3358 and 4852 atomic number, a relative molecular weight between 24060.16 and 34624.54 Da, and an isoelectric point between 4.99 and 8.65. The RT-qPCR results showed that the 12 Stlhcb genes were expressed in a tissue-specific and time-dependent fashion under low light. The relative expression of the Stlhcb genes in the leaves was significantly higher than that in the stems and roots, and the relative expression of these genes first increased and then decreased with the prolongation of light exposure time. The Stcp24 gene with the highest expression was cloned, and an expression vector was constructed. A subcellular localization analysis was performed in tobacco and an overexpression experiment was performed in potato using an Agrobacterium-mediated method. The subcellular localization analysis showed that the protein encoded by Stcp24 was located in chloroplasts as expected. Overexpression of Stcp24 in transgenic potato increased the yield of potatoes and the content of chlorophyll a and b; increased the net photosynthetic rate, transpiration rate, stomatal conductance, electron transport efficiency, and semi-saturated light intensity; and promoted photosynthesis and plant growth. This study provides a reference for the study of the function of the potato light-harvesting pigment-binding protein gene family. It lays a foundation for further study of the mechanism of the photosynthesis of potato, improvement of the light energy utilization of potato, and molecular breeding of potato. [ABSTRACT FROM AUTHOR]

Background and Objectives Materials and Methods Results Conclusion Missense variants in exon 7 of the ABO gene can lead to the formation of cisAB alleles. These alleles encode glycosyltransferases (GTs) capable of synthesizing both A and B antigens. In this study, we report the discovery of a novel cisAB allele and characterize it at molecular, protein and serological levels.Blood and DNA samples from the proband and seven relatives were examined using standard and modified ABO phenotyping, polymerase chain reaction‐restriction fragment length polymorphism and ABO gene sequencing. We assessed the impact of the p.Leu324Ser variant on the protein structure of the mutant GT using bioinformatics tools.Molecular tests revealed a c.971T>C (p.Leu324Ser) variant in the ABO gene in five of the eight individuals. This variant results in a GT that produces more A antigens and fewer B antigens. Bioinformatics analysis suggests that the amino acid substitution (p.Leu324Ser) could potentially affect enzymatic activity and specificity of the GT.We identified a novel cisAB allele resulting from a c.971T>C variant in the ABO gene. This variant led to the expression of an ABweak phenotype. [ABSTRACT FROM AUTHOR]

Background: Hylurgus ligniperda (Coleoptera: Curculionidae) is a worldwide forest quarantine pest. It is widely distributed, has many host tree species, and possesses strong adaptability. To explore its environmental adaptability and the related molecular mechanisms, we conducted chromosome-level genome sequencing and analyzed the transcriptome under different environmental factors, identifying key expressed genes. Results: We employed PacBio, Illumina, and Hi-C sequencing techniques to assemble a 520 Mb chromosomal-level genome of H. ligniperda, obtaining an N50 of 39.97 Mb across 138 scaffolds. A total of 10,765 protein-coding genes were annotated after repeat masking. Fourteen chromosomes were identified, among which Hyli14 was determined to be the sex chromosome. Survival statistics were tested over various growth periods under high temperature and low humidity conditions. The maximum survival period of adults reached 292 days at 25 °C, 65% relative humidity. In comparison, the maximum survival period was 14 days under 35 °C, 65% relative humidity, and 106 days under 25°C, 40% relative humidity. This indicated that environmental stress conditions significantly reduced adults' survival period. We further conducted transcriptome analysis to screen for potentially influential differentially expressed genes, such as CYP450 and Histone. Subsequently, we performed gene family analysis to gain insights into their functions and interactions, such as CYP450 and Histone. CYP450 genes affected the detoxification metabolism of enzymes in the Cytochrome P450 pathway to adapt to different environments. Histone genes are involved in insect hormone biosynthesis and longevity-regulating pathways in H. ligniperda to adapt to environmental stress. Conclusions: The genome at the chromosome level of H. ligniperda was assembled for the first time. The mortality of H. ligniperda increased significantly at 35 ℃, 65% RH, and 25 ℃, 40% RH. CYP450 and Histone genes played an important role in response to environmental stress. This genome offers a substantial genetic resource for investigating the molecular mechanisms behind beetle invasion and spread. [ABSTRACT FROM AUTHOR]

SUMMARY 12‐Oxo‐phytodienoic acid reductases (OPRs) perform vital functions in plants. However, few studies have been reported in sugarcane (Saccharum spp.), and it is of great significance to systematically investigates it in sugarcane. Here, 61 ShOPRs, 32 SsOPRs, and 36 SoOPRs were identified from R570 (Saccharum spp. hybrid cultivar R570), AP85‐441 (Saccharum spontaneum), and LA‐purple (Saccharum officinarum), respectively. These OPRs were phylogenetically classified into four groups, with close genes similar structures. During evolution, OPR gene family was mainly expanded via whole‐genome duplications/segmental events and predominantly underwent purifying selection, while sugarcane OPR genes may function differently in response to various stresses. Further, ScOPR2, a tissue‐specific OPR, which was localized in cytoplasm and cell membrane and actively response to salicylic acid (SA), methyl jasmonate, and smut pathogen (Sporisorium scitamineum) stresses, was cloned from sugarcane. In addition, both its transient overexpression and stable overexpression enhanced the resistance of transgenic plants to pathogen infection, most probably through activating pathogen‐associated molecular pattern/pattern‐recognition receptor‐triggered immunity, producing reactive oxygen species, and initiating mitogen‐activated protein kinase cascade. Subsequently, the transmission of SA and hypersensitive reaction were triggered, which stimulated the transcription of defense‐related genes. These findings provide insights into the function of ScOPR2 gene for disease resistance. [ABSTRACT FROM AUTHOR]

The Venus flytrap sea anemone Actinoscyphia liui inhabits the nutrient‐limited deep ocean in the tropical western Pacific. Compared with most other sea anemones, it has undergone a distinct modification of body shape similar to that of the botanic flytrap. However, the molecular mechanism by which such a peculiar sea anemone adapts to a deep‐sea oligotrophic environment is unknown. Here, we report the chromosomal‐level genome of A. liui constructed from PacBio and Hi‐C data. The assembled genome is 522 Mb in size and exhibits a continuous scaffold N50 of 58.4 Mb. Different from most other sea anemones, which typically possess 14–18 chromosomes per haplotype, A. liui has only 11. The reduced number of chromosomes is associated with chromosome fusion, which likely represents an adaptive strategy to economize energy in oligotrophic deep‐sea environments. Comparative analysis with other deep‐sea sea anemones revealed adaptive evolution in genes related to cellular autophagy (TMBIM6, SESN1, SCOCB and RPTOR) and mitochondrial energy metabolism (MDH1B and KAD2), which may aid in A. liui coping with severe food scarcity. Meanwhile, the genome has undergone at least two rounds of expansion in gene families associated with fast synaptic transmission, facilitating rapid responses to water currents and prey. Positive selection was detected on putative phosphorylation sites of muscle contraction‐related proteins, possibly further improving feeding efficiency. Overall, the present study provides insights into the molecular adaptation to deep‐sea oligotrophic environments and sheds light upon potential effects of a novel morphology on the evolution of Cnidaria. [ABSTRACT FROM AUTHOR]

Biological invasions pose a significant threat to ecosystems, disrupting local biodiversity and ecosystem functions. The genomic underpinnings of invasiveness, however, are still largely unknown, making it difficult to predict and manage invasive species effectively. The common reed (Phragmites australis) is a dominant grass species in wetland ecosystems and has become particularly invasive when transferred from Europe to North America. Here, we present a high-quality gap-free, telomere-to-telomere genome assembly of Phragmites australis consisting of 24 pseudochromosomes and a B chromosome. Fully phased subgenomes demonstrated considerable subgenome dominance and revealed the divergence of diploid progenitors approximately 30.9 million years ago. Comparative genomics using chromosome-level scaffolds for three other lineages and a previously published draft genome assembly of an invasive lineage revealed that gene family expansions in the form of tandem duplications may have contributed to the invasiveness of the lineage. This study sheds light on the genome evolution of Arundinoideae grasses and suggests that genetic drivers, such as gene family expansions and tandem duplications, may underly the processes of biological invasion in plants. These findings provide a crucial step toward understanding and managing the genetic basis of invasiveness in plant species. We provide a complete, telomere to telomere gap-free genome assembly for common reed that includes its B chromosome. Comparative genomic analysis revealed that tandem duplications may have contributed to the invasiveness potential of the species. [ABSTRACT FROM AUTHOR]

Chrysosplenium sinicum, a traditional Tibetan medicinal plant, can successfully thrive in low-light environments for long periods of time. To investigate the adaptive evolution of shade plants in low-light environments, we generated a chromosome-scale genome assembly (~320 Mb) for C. sinicum by combining PacBio sequencing and Hi-C technologies. Based on our results, gene families related to photosynthesis and cell respiration greatly expanded and evolved in C. sinicum genome due to intracellular DNA transfer from organelle genome to nuclear genome. Under positive selective pressure, adaptive evolution of light-harvesting complex II (LHCII) component protein CsLhcb1s resulted in the expansion of threonine residues at the phosphorylation site of STN7 kinase, potentially establishing a crucial genomic foundation for enhancing C. sinicum's adaptability in low-light environments. Through transcriptome and metabolome analysis, we identified chrysosplenol and chrysosplenoside as predominant flavonoid metabolites of C. sinicum and predicted their synthesis pathways. In addition, analysis of alternative splicing (AS) revealed that AS events help regulate state transition and flavonoid biosynthesis. The present study provides new insights into the genomes of shade plants exposed to low-light conditions and adaptive evolution of these genomes; in addition, the results improve our current knowledge on the biosynthetic and regulatory processes of chrysosplenol and chrysosplenoside. A chromosome-level genome assembly for Chrysosplenium sinicum combined with transcriptomic and metabolomic data provides insights into the adaptive evolution of shade plants as well as the pathways responsible for biosynthesis of chrysosplenol and chrysosplenoside. [ABSTRACT FROM AUTHOR]

Background: Deoxynivalenol (DON) is a type B trichothecene mycotoxin that is commonly found in cereals and grains worldwide. The presence of this fungal secondary-metabolite raises public-health concerns at both the agriculture and food industry level. Recently, we have shown that DON has a negative impact on gut integrity, a feature also noticed for Campylobacter (C.) jejuni. We further demonstrated that DON increased the load of C. jejuni in the gut and inner organs. In contrast, feeding the less toxic DON metabolite deepoxy-deoxynivalenol (DOM-1) to broilers reduced the Campylobacter load in vivo. Consequently, it can be hypothesized that DON and DOM-1 have a direct effect on the growth profile of C. jejuni. The aim of the present study was to further resolve the nature of this interaction in vitro by co-incubation and RNA-sequencing. Results: The co-incubation of C. jejuni with DON resulted in significantly higher bacterial growth rates from 30 h of incubation onwards. On the contrary, the co-incubation of C. jejuni with DOM-1 reduced the CFU counts, indicating that this DON metabolite might contribute to reduce the burden of C. jejuni in birds, altogether confirming in vivo data. Furthermore, the transcriptomic profile of C. jejuni following incubation with either DON or DOM-1 differed. Co-incubation of C. jejuni with DON significantly increased the expression of multiple genes which are critical for Campylobacter growth, particularly members of the Flagella gene family, frr (ribosome-recycling factor), PBP2 futA-like (Fe3+ periplasmic binding family) and PotA (ATP-binding subunit). Flagella are responsible for motility, biofilm formation and host colonization, which may explain the high Campylobacter load in the gut of DON-fed broiler chickens. On the contrary, DOM-1 downregulated the Flagella gene family and upregulated ribosomal proteins. Conclusion: The results highlight the adaptive mechanisms involved in the transcriptional response of C. jejuni to DON and its metabolite DOM-1, based on the following effects: (a) ribosomal proteins; (b) flagellar proteins; (c) engagement of different metabolic pathways. The results provide insight into the response of an important intestinal microbial pathogen against DON and lead to a better understanding of the luminal or environmental acclimation mechanisms in chickens. [ABSTRACT FROM AUTHOR]

Heat shock protein 90 (HSP90) is important for many organisms, including plants. Based on the whole genome information, the gene number, gene structure, evolutionary relationship, protein structure, and active site of the HSP90 gene family in Rosa chinensis and Rubus idaeus were determined, and the expression of the HSP90 gene under salt, and drought stresses in two rose varieties Wangxifeng and Sweet Avalanche were analyzed. Six and eight HSP90 genes were identified from R. chinensis and Ru. idaeus, respectively. Phylogenetic analysis revealed that the analyzed genes were divided into two Groups and four subgroups (Classes 1a, 1b, 2a, and 2b). Although members within the same classes displayed highly similar gene structures, while the gene structures and conserved domains of Group 1 (Class 1a and 1b) and the Group 2 (Class 2a and 2b) are different. Tandem and segmental duplication genes were found in Ru. idaeus, but not in R. chinensis, perhaps explaining the difference in HSP90 gene quantity in the two analyzed species. Analysis of cis-acting elements revealed abundant abiotic stress, photolight-response, and hormone-response elements in R. chinensis HSP90s. qRT-PCR analysis suggested that RcHSP90-1–1, RcHSP90-5–1 and RcHSP90-6–1 in Sweet Avalanche and Wangxifeng varieties played important regulatory roles under salt and drought stress. The analysis of protein structure and active sites indicate that the potential different roles of RcHSP90-1–1, RcHSP90-5–1, and RcHSP90-6–1 in salt and drought stresses may come from the differences of corresponding protein structures and activation sites. These data will provide information for the breeding of rose varieties with high stress resistance. [ABSTRACT FROM AUTHOR]

Evolution of unisexual flowers involves extreme changes in floral development. Spinach is one of the species to discern the formation and evolution of dioecy. MADS-box gene family is involved in regulation of floral organ identity and development and in many other plant developmental processes. However, there is no systematic analysis of MADS-box family genes in spinach. A comprehensive genome-wide analysis and transcriptome profiling of MADS-box genes were undertaken to understand their involvement in unisexual flower development at different stages in spinach. In total, 54 MADS-box genes found to be unevenly located across 6 chromosomes and can be divided into type I and type II genes. Twenty type I MADS-box genes are subdivided into Mα, Mβ and Mγ subgroups. While thirty-four type II SoMADSs consist of 3 MIKC*, and 31 MIKCC -type genes including sixteen floral homeotic MADS-box genes that are orthologous to the proposed Arabidopsis ABCDE model of floral organ identity determination, were identified in spinach. Gene structure, motif distribution, physiochemical properties, gene duplication and collinearity analyses for these genes are performed in detail. Promoters of both types of SoMADS genes contain mainly MeJA and ABA response elements. Expression profiling indicated that MIKCc genes exhibited more dynamic and intricate expression patterns compared to M-type genes and the majority of type-II genes AP1, SVP, and SOC1 sub-groups showed female flower-biased expression profiles, suggesting their role in carpel development, while PI showed male-biased expression throughout flower developmental stages, suggesting their role in stamen development. These results provide genomic resources and insights into spinach dioecious flower development and expedite spinach improvement. [ABSTRACT FROM AUTHOR]

Chromodomain helicase DNA-binding protein (CHD) gene family, an ATP (adenosine triphosphate) -dependent chromatin remodeler family, is involved in multiple developmental process and tumor development. However, there have been none pan-cancer analyses of this family. The expression levels, survival profiles, mutation profiles and immune infiltration of the CHD family genes from TCGA and TARGET database were analyzed using online tools or R packages. Interestingly, all types of CHD gene expressions were associated with the prognosis of Neuroblastoma, Acute lymphoblastic leukemia-Phase 3 and Acute Myeloid Leukemia (All P < 0.05). Knock down of CHD7 and CHD9 in K562 (human erythromyeloblastoid leukemia) and HEC-1-B (human endometrial adenocarcinoma) cells significantly inhibit cell proliferation and migration (P < 0.05). Proliferation, colony formation and migration assays were performed in CHD7 and CHD9 knockdown K562 and HBC-1-B cell lines. Mechanisms were also analyzed by PPI and GO ontology for our experiments. Histone modification, especially the methylation of H3K4, might be involved in CHD7 and CHD9 related oncogenesis. Through bioinformatic analysis, we showed CHD genes significantly affected the prognosis of different tumor types, including childhood tumor. Our findings provide new insights into the function and mechanism of CHD gene family, especially in CHD7 and CHD9. [ABSTRACT FROM AUTHOR]

The increasing frequency and severity of heatwaves will intensify stress on plants. Given regional variation in heatwave exposure and expected differences in thermal tolerance between species it is unlikely that all plant species will be affected equally by climate change. However, little is currently known about variation in the responses of plants to heat stress, or how those responses differ among closely related species adapted to different environments. Here we quantify the response of 17 Acacia species (175 RNA‐seq libraries), from across Australia's diverse biomes, to a multi‐day experimental heatwave treatment to identify variation in transcriptomic and physiological responses to heat stress. Genes with known heat response functions showed consistent responses across Acacia species. Up to 10% of all genes and over 100 gene families showed significant clinal variation in the magnitude of their expression plasticity across species. Specifically, gene families linked to the temperature stress response were overrepresented among significant relationships with home range temperature conditions. Gene expression responses seen on the first day of the heatwave were more frequently associated with home range climates, while expression responses by day four were more commonly related to photosystem II acclimation. Comparative transcriptomics on non‐model species has the potential to provide key information on stress response plasticity, especially when linked with our understanding of model species. Our study indicates that the pressing challenge to identifying potentially vulnerable species to climate change could be benefited by the further exploration of clinal variation in transcriptome plasticity. [ABSTRACT FROM AUTHOR]

DNA-binding one zinc finger (DOF) transcription factors are crucial plant-specific regulators involved in growth, development, signal transduction, and abiotic stress response generation. However, the genome-wide identification and characterization of AcDOF genes and their regulatory elements in kiwifruit (Actinidia chinensis) has not been thoroughly investigated. In this study, we screened the kiwifruit genome database and identified 42 AcDOF genes (AcDOF1 to AcDOF42). Phylogenetic analysis facilitated the categorization of these genes into five subfamilies (DOF-a, DOF-b, DOF-c, DOF-d, and DOF-e). We further analyzed the motifs, conserved domains, gene structures, and collinearity of the AcDOFgene family. Gene ontology (GO) enrichment analysis indicated significant enrichment in the "flower development" term and the "response to abiotic stress" category. Promoter prediction analysis revealed numerous cis-regulatory elements related to responses to light, hormones, and low-temperature and drought stress in AcDOF promoters. RNA-seq expression profiles demonstrated the tissue-specific expression of AcDOF genes. Quantitative real-time PCR results showed that six selected genes (AcDOF04, AcDOF09, AcDOF11, AcDOF13, AcDOF21, and AcDOF22) were differentially induced by abscisic acid (ABA), methyl jasmonate (MeJA), and cold, salt, and drought stresses, with AcDOF22 specifically expressed at high levels in drought-tolerant cultivars. Further experiments indicated that transient AcDOF22 overexpression in kiwifruit leaf disks reduced water loss and chlorophyll degradation. Additionally, AcDOF22 was localized to the nucleus and exhibited transcriptional activation, enhancing drought resistance by activating the downstream drought marker gene AcDREB2A. These findings lay the foundation for elucidating the molecular mechanisms of drought resistance in kiwifruit and offer new insights into drought-resistant breeding. [ABSTRACT FROM AUTHOR]

Sulfate transporters (SULTRs) are essential for the transport and absorption of sulfate in plants and serve as critical transport proteins within the sulfur metabolism pathway, significantly influencing plant growth, development, and stress adaptation. A bioinformatics analysis of SULTR genes in soybean was performed, resulting in the identification and classification of twenty-eight putative GmSULTRs into four distinct groups. In this study, the characteristics of the 28 GmSULTR genes, including those involved in collinearity, gene structure, protein motifs, cis-elements, tissue expression patterns, and the response to abiotic stress and plant hormone treatments, were systematically analyzed. This study focused on conducting a preliminary functional analysis of the GmSULTR3;1a gene, wherein a high expression level of GmSULTR3;1a in the roots, stems, and leaves was induced by a sulfur deficiency and GmSULTR3;1a improved the salt tolerance. A further functional characterization revealed that GmSULTR3;1a-overexpressing soybean hairy roots had higher SO42−, GSH, and methionine (Met) contents compared with the wild-type (WT) plant. These results demonstrate that the overexpression of GmSULTR3;1a may promote the sulfur assimilation metabolism and increase the content of sulfur-containing amino acids in plants. [ABSTRACT FROM AUTHOR]

The BTB (Broad-complex, tramtrack, and bric-a-brac) gene family, characterized by a highly conserved BTB domain, is implicated in a spectrum of biological processes, encompassing growth and development, as well as stress responses. Characterization and functional studies of BTB genes in poplar are still limited, especially regarding their response to hormones and biotic/abiotic stresses. In this study, we conducted an HMMER search in conjunction with BLASTp and identified 95 BTB gene models in Populus trichocarpa. Through domain motif and phylogenetic relationship analyses, these proteins were classified into eight families, NPH3, TAZ, Ankyrin, only BTB, BACK, Armadillo, TPR, and MATH. Collinearity analysis of poplar BTB genes with homologs in six other species elucidated evolutionary relationships and functional conservations. RNA-seq analysis of five tissues of poplar identified BTB genes as playing a pivotal role during developmental processes. Comprehensive RT-qPCR analysis of 11 BTB genes across leaves, roots, and xylem tissues revealed their responsive expression patterns under diverse hormonal and biotic/abiotic stress conditions, with varying degrees of regulation observed in the results. This study marks the first in-depth exploration of the BTB gene family in poplar, providing insights into the potential roles of BTB genes in hormonal regulation and response to stress. [ABSTRACT FROM AUTHOR]

CCHC-type zinc finger proteins (CCHC-ZFPs), ubiquitous across plant species, are integral to their growth, development, hormonal regulation, and stress adaptation. Roses (Rosa sp.), as one of the most significant and extensively cultivated ornamentals, account for more than 30% of the global cut-flower market. Despite its significance, the CCHC gene family in roses (Rosa sp.) remains unexplored. This investigation identified and categorized 41 CCHC gene members located on seven chromosomes of rose into 14 subfamilies through motif distribution and phylogenetic analyses involving ten additional plant species, including Ginkgo biloba, Ostreococcus lucimarinus, Arabidopsis thaliana, and others. This study revealed that dispersed duplication likely plays a crucial role in the diversification of the CCHC genes, with the Ka/Ks ratio suggesting a history of strong purifying selection. Promoter analysis highlighted a rich presence of cis-acting regulatory elements linked to both abiotic and biotic stress responses. Differential expression analysis under drought conditions grouped the 41 CCHC gene members into five distinct clusters, with those in group 4 exhibiting pronounced regulation in roots and leaves under severe drought. Furthermore, virus-induced gene silencing (VIGS) of the RcCCHC25 member from group 4 compromised drought resilience in rose foliage. This comprehensive analysis lays the groundwork for further investigations into the functional dynamics of the CCHC gene family in rose physiology and stress responses. [ABSTRACT FROM AUTHOR]

Cytokinins (CKs) are a group of phytohormones that are involved in plant growth, development, and disease resistance. The isopentenyl transferase (IPT) and cytokinin oxidase/dehydrogenase (CKX) families comprise key enzymes controlling CK biosynthesis and degradation. However, an integrated analysis of these two gene families in radish has not yet been explored. In this study, 13 RsIPT and 12 RsCKX genes were identified and characterized, most of which had four copies in Brassica napus and two copies in radish and other diploid Brassica species. Promoter analysis indicated that the genes contained at least one phytohormone or defense and stress responsiveness cis-acting element. RsIPTs and RsCKXs were expanded through segmental duplication. Moreover, strong purifying selection drove the evolution of the two gene families. The expression of the RsIPT and RsCKX genes distinctly showed diversity in different tissues and developmental stages of the root. Expression profiling showed that RsCKX1-1/1-2/1-3 was significantly upregulated in club-resistant materials during primary infection, suggesting their vital function in clubroot resistance. The interaction network of CKX proteins with similar 3D structures also reflected the important role of RsCKX genes in disease resistance. This study provides a foundation for further functional study on the IPT and CKX genes for clubroot resistance improvement in Raphanus. [ABSTRACT FROM AUTHOR]

The evolution of the MADS-box gene family is essential for the rapid differentiation of floral organs and fruit types in angiosperms. Two key processes drive the evolution of gene families: gene duplication and functional differentiation. Duplicated copies provide the material for variation, while advantageous mutations can confer new functions on gene copies. In this study, we selected the Rosaceae family, which includes a variety of fruit types and flower organs, as well as species that existed before and after whole-genome duplication (WGD). The results indicate that different fruit types are associated with different copies of MADS-box gene family duplications and WGD events. While most gene copies derived from WGD have been lost, MADS-box genes not only retain copies derived from WGD but also undergo further gene duplication. The sequences, protein structures, and expression patterns of these gene copies have undergone significant differentiation. This work provides a clear example of MADS-box genes in the context of gene duplication and functional differentiation, offering new insights into the evolution of fruit types and floral organs. [ABSTRACT FROM AUTHOR]

Leaves are a key forage part for livestock, and the aging of leaves affects forage biomass and quality. Preventing or delaying premature leaf senescence leads to an increase in pasture biomass accumulation and an improvement in alfalfa quality. NAC transcription factors have been reported to affect plant growth and abiotic stress responses. In this study, 48 NAC genes potentially associated with leaf senescence were identified in alfalfa under dark or salt stress conditions. A phylogenetic analysis divided MsNACs into six subgroups based on similar gene structure and conserved motif. These MsNACs were unevenly distributed in 26 alfalfa chromosomes. The results of the collinearity analysis show that all of the MsNACs were involved in gene duplication. Some cis-acting elements related to hormones and stress were screened in the 2-kb promoter regions of MsNACs. Nine of the MsNAC genes were subjected to qRT-PCR to quantify their expression and Agrobacterium-mediated transient expression to verify their functions. The results indicate that Ms.gene031485, Ms.gene032313, Ms.gene08494, and Ms.gene77666 might be key NAC genes involved in alfalfa leaf senescence. Our findings extend the understanding of the regulatory function of MsNACs in leaf senescence. [ABSTRACT FROM AUTHOR]

The R2R3-MYB gene family represents a widely distributed class of plant transcription factors. This gene family plays an important role in many aspects of plant growth and development. However, the characterization of R2R3-MYB genes present in the genome of Coptis teeta has not been reported. Here, we describe the bioinformatic identification and characterization of 88 R2R3-MYB genes in this species, and the identification of members of the R2R3-MYB gene family in species within the order Ranales most closely related to Coptis teeta. The CteR2R3-MYB genes were shown to exhibit a higher degree of conservation compared to those of A. thaliana, as evidenced by phylogeny, conserved motifs, gene structure, and replication event analyses. Cis-acting element analysis confirmed the involvement of CteR2R3-MYB genes in a variety of developmental processes, including growth, cell differentiation, and reproduction mediated by hormone synthesis. In addition, through homology comparisons with the equivalent gene family in A. thaliana, protein regulatory network prediction and transcriptome data analysis of floral organs across three time periods of flower development, 17 candidate genes were shown to exhibit biased expression in two floral phenotypes of C. teeta. This suggests their potential involvement in floral development (anther development) in this species. [ABSTRACT FROM AUTHOR]

The indeterminate domain proteins (IDD proteins) play essential roles in the growth and development of various plant tissues and organs across different developmental stages, but members of this gene family have not yet been characterized in foxtail millet (Setaria italica). To have a comprehensive understanding of the IDD gene family in foxtail millet, we performed a genome-wide characterization and haplotypic variation analysis of the IDD gene family in foxtail millet. In this study, sixteen IDD genes were identified across the reference genome of Yugu1, a foxtail millet cultivar. Phylogenetic analysis revealed that the Setaria italica IDD (SiIDD) proteins were clustered into four groups together with IDD proteins from Arabidopsis thaliana (dicot) and Oryza sativa (monocot). Conserved protein motif and gene structure analyses revealed that the closely clustered SiIDD genes were highly conserved within each subgroup. Furthermore, chromosomal location analysis showed that the SiIDD genes were unevenly distributed on nine chromosomes of foxtail millet and shared collinear relationships with IDD genes of other grass species. Transcriptional analysis revealed that the SiIDD genes differed greatly in their expression patterns, and paralogous genes shared similar expression patterns. In addition, superior haplotypes for two SiIDD genes (SiIDD8 and SiIDD14) were identified to correlate with traits of early heading date, and high thousand seed weight and molecular markers were designed for SiIDD8 and SiIDD14 to distinguish different haplotypes for breeding. Taken together, the results of this study provide useful information for further functional investigation of SiIDD genes, and the superior haplotypes of SiIDD8 and SiIDD14 will be particularly beneficial for improving heading date and yield of foxtail millet in breeding programs. [ABSTRACT FROM AUTHOR]

The plerocercoid larva of Spirometra mansoni can cause a parasitic zoonosis—sparganosis. Malate dehydrogenase (MDH) plays a very important role in the life activities of parasites. However, little is known about the MDH family in S. mansoni. We identified eight new MDH members in S. mansoni in this study. Clustering analysis divided SmMDHs into two groups and revealed patterns similar to the conserved motif organization. RT–qPCR suggested that five MDHs were highly expressed in the mature proglottid and that three MDHs were highly expressed in the gravid proglottid. Phylogenetic analysis revealed that SmMDHs contain both conserved family members and members in the process of further diversification. rSmMDH has an NAD binding domain, a dimer interface and a substrate binding domain. Natural SmMDH was immunolocalized in the tissues and follicles around the uterus in the mature or gravid proglottid and eggshells. The maximum forward and reverse reaction activities of rSmMDH were observed at pH 8.5 and 9.0, respectively. The optimum temperature for enzyme activity was 37 °C in the forward reaction and 40 °C in the reverse reaction. These results lay the foundation for studying the molecular functions and mechanisms of MDHs in S. mansoni and related taxa. [ABSTRACT FROM AUTHOR]

Evolution of unisexual flowers involves extreme changes in floral development. Spinach is one of the species to discern the formation and evolution of dioecy. MADS-box gene family is involved in regulation of floral organ identity and development and in many other plant developmental processes. However, there is no systematic analysis of MADS-box family genes in spinach. A comprehensive genome-wide analysis and transcriptome profiling of MADS-box genes were undertaken to understand their involvement in unisexual flower development at different stages in spinach. In total, 54 MADS-box genes found to be unevenly located across 6 chromosomes and can be divided into type I and type II genes. Twenty type I MADS-box genes are subdivided into Mα, Mβ and Mγ subgroups. While thirty-four type II SoMADSs consist of 3 MIKC*, and 31 MIKCC -type genes including sixteen floral homeotic MADS-box genes that are orthologous to the proposed Arabidopsis ABCDE model of floral organ identity determination, were identified in spinach. Gene structure, motif distribution, physiochemical properties, gene duplication and collinearity analyses for these genes are performed in detail. Promoters of both types of SoMADS genes contain mainly MeJA and ABA response elements. Expression profiling indicated that MIKCc genes exhibited more dynamic and intricate expression patterns compared to M-type genes and the majority of type-II genes AP1, SVP, and SOC1 sub-groups showed female flower-biased expression profiles, suggesting their role in carpel development, while PI showed male-biased expression throughout flower developmental stages, suggesting their role in stamen development. These results provide genomic resources and insights into spinach dioecious flower development and expedite spinach improvement. [ABSTRACT FROM AUTHOR]

Background: Tubulins are major components of the eukaryotic cytoskeletons that are crucial in many cellular processes. Ciliated protists comprise one of the oldest eukaryotic lineages possessing cilia over their cell surface and assembling many diverse microtubular structures. As such, ciliates are excellent model organisms to clarify the origin and evolution of tubulins in the early stages of eukaryote evolution. Nonetheless, the evolutionary history of the tubulin subfamilies within and among ciliate classes is unclear. Results: We analyzed the evolutionary pattern of ciliate tubulin gene family based on genomes/transcriptomes of 60 species covering 10 ciliate classes. Results showed: (1) Six tubulin subfamilies (α_Tub, β_Tub, γ_Tub, δ_Tub, ε_Tub, and ζ_Tub) originated from the last eukaryotic common ancestor (LECA) were observed within ciliates. Among them, α_Tub, β_Tub, and γ_Tub were present in all ciliate species, while δ_Tub, ε_Tub, and ζ_Tub might be independently lost in some species. (2) The evolutionary history of the tubulin subfamilies varied. Evolutionary history of ciliate γ_Tub, δ_Tub, ε_Tub, and ζ_Tub showed a certain degree of consistency with the phylogeny of species after the divergence of ciliate classes, while the evolutionary history of ciliate α_Tub and β_Tub varied among different classes. (3) Ciliate α- and β-tubulin isoforms could be classified into an "ancestral group" present in LECA and a "divergent group" containing only ciliate sequences. Alveolata-specific expansion events probably occurred within the "ancestral group" of α_Tub and β_Tub. The "divergent group" might be important for ciliate morphological differentiation and wide environmental adaptability. (4) Expansion events of the tubulin gene family appeared to be consistent with whole genome duplication (WGD) events in some degree. More Paramecium-specific tubulin expansions were detected than Tetrahymena-specific ones. Compared to other Paramecium species, the Paramecium aurelia complex underwent a more recent WGD which might have experienced more tubulin expansion events. Conclusions: Evolutionary history among different tubulin gene subfamilies seemed to vary within ciliated protists. And the complex evolutionary patterns of tubulins among different ciliate classes might drive functional diversification. Our investigation provided meaningful information for understanding the evolution of tubulin gene family in the early stages of eukaryote evolution. [ABSTRACT FROM AUTHOR]