MCScanX: a toolkit for detection and evolutionary analysis of gene synteny and collinearity

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

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 identification of YABBY gene family and its expression pattern analysis in Astragalus mongholicus

During plant growth and development, the YABBY gene plays a crucial role in the morphological structure, hormone signaling, stress resistance, crop breeding, and agricultural production of plant lateral organs, leaves, flowers, and fruits. Astragalus mongholicus is a perennial herbaceous plant in the legume family, widely used worldwide due to its high medicinal and edible value. However, there have been no reports of the YABBY gene family in A. mongholicus. This study used bioinformatics methods, combined with databases and analysis websites, to systematically analyze the AmYABBY gene family in the entire genome of A. mongholicus and verified its expression patterns in different tissues of A. mongholicus through transcriptome data and qRT-PCR experiments. A total of seven AmYABBY genes were identified, which can be divided into five subfamilies and distributed on three chromosomes. Two pairs of AmYABBY genes may be involved in fragment duplication on three chromosomes. All AmYABBY proteins have a zinc finger YABBY domain, and members of the same group have similar motif composition and intron - exon structure. In the promoter region of the genes, light-responsive and MeJa-response cis-elements are dominant. AmYABBY is highly expressed in stems and leaves, especially AmYABBY1, AmYABBY2, and AmYABBY3, which play important roles in the growth and development of stems and leaves. The AmYABBY gene family regulates the growth and development of A. mongholicus. In summary, this study provides a theoretical basis for in-depth research on the function of the AmYABBY gene and new insights into the molecular response mechanism of the growth and development of the traditional Chinese medicine A. mongholicus.

Genome-wide identification of the Hsp70 gene family in Penaeus chinensis and their response to environmental stress

The heat shock protein 70 (HSP70) gene family plays a crucial role in the response of organisms to environmental stress. However, it has not been systematically characterized in shrimp. In this study, we identified 25 PcHsp70 genes in the Penaeus chinensis genome. The encoded proteins were categorized into six subgroups based on phylogenetic relationships. Tandem duplication was the main driver of amplification in the PcHsp70 family, and the genes have experienced strong purifying selection during evolution. Transcriptome data analysis revealed that the 25 PcHsp70 members have different expression patterns in shrimp under conditions of low temperature, low salinity, and white spot syndrome virus infection. Among them, PcHsp70.11 was significantly induced under all three stress conditions, suggesting that this gene plays an important role in response to environmental stress in P. chinensis. To the best of our knowledge, this is the first study to systematically analyze the Hsp70 gene family in shrimp. The results provide important information on shrimp Hsp70s, contributing to a better understanding of the role of these genes in environmental stress and providing a basis for further functional studies.

Identification of cotton PIP5K genes and role of GhPIP5K9a in primary root development

Phosphatidylinositol 4 phosphate 5-kinase (PIP5K) is crucial for the phosphatidylinositol (PI) signaling pathway. It plays a significant role in plant growth and development, as well as stress response. However, its effects on cotton are unknown. This study identified PIP5K genes from four cotton species and conducted bioinformatic analyses, with a particular emphasis on the functions of GhPIP5K9a in primary roots. The results showed that cotton PIP5Ks were classified into four subgroups. Analysis of gene structure and motif composition showed obvious conservation within each subgroup. Synteny analysis suggested that the PIP5K gene family experienced significant expansion due to both whole-genome duplication (WGD) and segmental duplication. Transcriptomic data analysis revealed that the majority of GhPIP5K genes had the either low or undetectable levels of expression. Moreover, GhPIP5K9a is highly expressed in the root and was located in plasmalemma. Suppression of GhPIP5K9a transcripts resulted in longer primary roots, longer primary root cells and increased auxin polar transport-related genes expression, and decreased abscisic acid (ABA) content, indicating that GhPIP5K9a negatively regulates cotton primary root growth. This study lays the foundation for further exploration of the role of the PIP5K genes in cotton.

Syntrichia ruralis: emerging model moss genome reveals a conserved and previously unknown regulator of desiccation in flowering plants

Water scarcity, resulting from climate change, poses a significant threat to ecosystems. Syntrichia ruralis, a dryland desiccation-tolerant moss, provides valuable insights into survival of water-limited conditions. We sequenced the genome of S. ruralis, conducted transcriptomic analyses, and performed comparative genomic and transcriptomic analyses with existing genomes and transcriptomes, including with the close relative S. caninervis. We took a genetic approach to characterize the role of an S. ruralis transcription factor, identified in transcriptomic analyses, in Arabidopsis thaliana. The genome was assembled into 12 chromosomes encompassing 21 169 protein-coding genes. Comparative analysis revealed copy number and transcript abundance differences in known desiccation-associated gene families, and highlighted genome-level variation among species that may reflect adaptation to different habitats. A significant number of abscisic acid (ABA)-responsive genes were found to be negatively regulated by a MYB transcription factor (MYB55) that was upstream of the S. ruralis ortholog of ABA-insensitive 3 (ABI3). We determined that this conserved MYB transcription factor, uncharacterized in Arabidopsis, acts as a negative regulator of an ABA-dependent stress response in Arabidopsis. The new genomic resources from this emerging model moss offer novel insights into how plants regulate their responses to water deprivation.

Terpene synthases GhTPS6 and GhTPS47 participate in resistance to Verticillium dahliae in upland cotton

Terpene synthases (TPSs) are enzymes responsible for catalyzing the production of diverse terpenes, the largest class of secondary metabolites in plants. Here, we identified 107 TPS gene loci encompassing 92 full-length TPS genes in upland cotton (Gossypium hirsutum L.). Phylogenetic analysis showed they were divided into six subfamilies. Segmental duplication and tandem duplication events contributed greatly to the expansion of TPS gene family, particularly the TPS-a and TPS-b subfamilies. Expression profile analysis screened out that GhTPSs may mediate the interaction between cotton and Verticillium dahliae. Three-dimensional structures and subcellular localizations of the two selected GhTPSs, GhTPS6 and GhTPS47, which belong to the TPS-a subfamily, demonstrated similarity in protein structures and nucleus and cytoplasm localization. Virus-induced gene silencing (VIGS) of the two GhTPSs yielded plants characterized by increased wilting and chlorosis, more severe vascular browning, and higher disease index than control plants. Additionally, knockdown of GhTPS6 and GhTPS47 led to the down-regulation of cotton terpene synthesis following V. dahliae infection, indicating that these two genes may positively regulate resistance to V. dahliae through the modulation of disease-resistant terpene biosynthesis. Overall, our study represents a comprehensive analysis of the G. hirsutum TPS gene family, revealing their potential roles in defense responses against Verticillium wilt.

Phyllostachys edulis argonaute genes function in the shoot architecture

Argonaute (AGO) proteins are the core components of the RNA-induced silencing complexes (RISC) in the cytoplasm and nucleus, and are necessary for the development of plant shoot meristem, which gives rise to the above-ground plant body. In this study, we identified 23 Phyllostachys edulis AGO genes (PhAGOs) that were distributed unequally on the 14 unmapped scaffolds. Gene collinearity and phylogeny analysis showed that the innovation of PhAGO genes was mainly due to dispersed duplication and whole-genome duplication, which resulted in the enlarged PhAGO family. PhAGO genes were expressed in a temporal-spatial expression pattern, and they encoded proteins differently localized in the cytoplasm and/or nucleus. Overexpression of the PhAGO2 and PhAGO4 genes increased the number of tillers or leaves in Oryza sativa and affected the shoot architecture of Arabidopsis thaliana. These results provided insight into the fact that PhAGO genes play important roles in plant development.

Comparative genomic and transcriptomic analyses provide new insight into symbiotic host specificity

Host specificity plays important roles in expanding the host range of rhizobia, while the genetic information responsible for host specificity remains largely unexplored. In this report, the roots of four symbiotic systems with notable different symbiotic phenotypes and the control were studied at four different post-inoculation time points by RNA sequencning (RNA-seq). The differentially expressed genes (DEGs) were divided into "found only in soybean or Lotus," "only expressed in soybean or Lotus," and "expressed in both hosts" according to the comparative genomic analysis. The distributions of enriched function ontology (GO) terms and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways vary significantly in different symbiotic systems. Host specific genes account for the majority of the DEGs involved in response to stimulus, associated with plant-pathogen interaction pathways, and encoding resistance (R) proteins, the symbiotic nitrogen fixation (SNF) proteins and the target proteins in the SNF-related modules. Our findings provided molecular candidates for better understanding the mechanisms of symbiotic host-specificity.

Genome-wide identification of MATE and ALMT genes and their expression profiling in mungbean (Vigna radiata L.) under aluminium stress

The Multidrug and toxic compound extrusion (MATE) and aluminium activated malate transporter (ALMT) gene families are involved in response to aluminium (Al) stress. In this study, we identified 48 MATE and 14 ALMT gene families in Vigna radiata genome and classified into 5 (MATE) and 3 (ALMT) clades by phylogenetic analysis. All the VrMATE and VrALMT genes were distributed across mungbean chromosomes. Tandem duplication was the main driving force for evolution and expansion of MATE gene family. Collinearity of mungbean with soybean indicated that MATE gene family is closely linked to Glycine max. Eight MATE transporters in clade 2 were found to be associated with previously characterized Al tolerance related MATEs in various plant species. Citrate exuding motif (CEM) was present in seven VrMATEs of clade 2. Promoter analysis revealed abundant plant hormone and stress responsive cis-elements. Results from quantitative real time-polymerase chain reaction (qRT-PCR) revealed that VrMATE19, VrMATE30 and VrALMT13 genes were markedly up-regulated at different time points under Al stress. Overall, this study offers a new direction for further molecular characterization of the MATE and ALMT genes in mungbean for Al tolerance.

Genome-wide analysis of the SWEET gene family in Hemerocallis citrina and functional characterization of HcSWEET4a in response to salt stress

Sugars will be eventually effluxed transporters (SWEETs) have been confirmed to play diverse physiological roles in plant growth, development and stress response. However, the characteristics and functions of the SWEET genes in Hemerocallis citrina remain unclear and poorly elucidated. In this study, the whole genome of Hemerocallis citrina was utilized to conduct bioinformatics analysis and a total of 19 HcSWEET genes were successfully identified. Analysis of the physicochemical properties indicated dominant differences among these HcSWEETs. A phylogenetic analysis revealed that HcSWEET proteins can be divided into 4 clades ranging from Clade I to IV, where proteins within the same clade exhibited shared conserved motifs and gene structures. Five to six exons were contained in the majority of HcSWEET genes, which were unevenly distributed across 11 chromosomes. The gene duplication analysis showed the presence of 4 gene pairs. Comparative syntenic maps revealed that the HcSWEET gene family might present more closed homology in monocotyledons than dicotyledons. Cis-acting element analysis of HcSWEET genes indicated key responsiveness to various hormones, light, and stresses. Additionally, transcriptome sequencing analysis suggested that most HcSWEET genes had a relatively higher expression in roots, and HcSWEET4a was significantly up-regulated under salt stress. Overexpression further verified the possibility that HcSWEET4a was involved in response to salt stress, which provides novel insights and facilitates in-depth studies of the functional analysis of HcSWEETs in resistance to abiotic stress.

Combined genome and transcriptome provides insight into the genetic evolution of an edible halophyte Suaeda salsa adaptation to high salinity

Suaeda salsa L. is a typical halophyte with high value as a vegetable. Here, we report a 447.98 Mb, chromosomal-level genome of S. salsa, assembled into nine pseudomolecules (contig N50 = 1.36 Mb) and annotated with 27,927 annotated protein-coding genes. Most of the assembled S. salsa genome, 58.03%, consists of transposable elements. Some gene families including HKT1, NHX, SOS and CASP related to salt resistance were significantly amplified. We also observed expansion of genes encoding protein that bind the trace elements Zn, Fe, Cu and Mn, and genes related to flavonoid and α-linolenic acid metabolism. Many expanded genes were significantly up-regulated under salinity, which might have contributed to the acquisition of salt tolerance in S. salsa. Transcriptomic data showed that high salinity markedly up-regulated salt-resistance related genes, compared to low salinity. Abundant metabolic pathways of secondary metabolites including flavonoid, unsaturated fatty acids and selenocompound were enriched, which indicates that the species is a nutrient-rich vegetable. Particularly worth mentioning is that there was no significant difference in the numbers of cis-elements in the promoters of salt-related and randomly selected genes in S. salsa when compared with Arabidopsis thaliana, which may affirm that plant salt tolerance is a quantitative rather than a qualitative trait in terms of promoter evolution. Our findings provide deep insight into the adaptation of halophytes to salinity from a genetic evolution perspective.

Evolution of FLOWERING LOCUS T-like genes in angiosperms: a core Lamiales-specific diversification

Plant life history is determined by two transitions, germination and flowering time, in which the phosphatidylethanolamine-binding proteins (PEBPs) FLOWERING LOCUS T (FT) and TERMINAL FLOWER1 (TFL1) play key regulatory roles. Compared with the highly conserved TFL1-like genes, FT-like genes vary significantly in copy numbers in gymnosperms, and monocots within the angiosperms, while sporadic duplications can be observed in eudicots. Here, via a systematic analysis of the PEBPs in angiosperms with a special focus on 12 representative species featuring high-quality genomes in the order Lamiales, we identified a successive lineage-specific but systematic expansion of FT-like genes in the families of core Lamiales. The first expansion event generated FT1-like genes mainly via a core Lamiales-specific whole-genome duplication (cL-WGD), while a likely random duplication produced the FT2-like genes in the lineages containing Scrophulariaceae and the rest of the core Lamiales. Both FT1- and FT2-like genes were further amplified tandemly in some families. These expanded FT-like genes featured highly diverged expression patterns and structural variation, indicating functional diversification. Intriguingly, some core Lamiales contained the relict MOTHER OF FT AND TFL1 like 2 (MFT2) that probably expanded in the common ancestor of angiosperms. Our data showcase the highly dynamic lineage-specific expansion of the FT-like genes, and thus provide important and fresh evolutionary insights into the gene regulatory network underpinning flowering time diversity in Lamiales and, more generally, in angiosperms.

Genome-wide identification and expression analysis of the PP2C gene family in Apocynum venetum and Apocynum hendersonii

BACKGROUND

Protein phosphatase class 2 C (PP2C) is the largest protein phosphatase family in plants. Members of the PP2C gene family are involved in a variety of physiological pathways in plants, including the abscisic acid signalling pathway, the regulation of plant growth and development, etc., and are capable of responding to a wide range of biotic and abiotic stresses, and play an important role in plant growth, development, and response to stress. Apocynum is a perennial persistent herb, divided into Apocynum venetum and Apocynum hendersonii. It mainly grows in saline soil, deserts and other harsh environments, and is widely used in saline soil improvement, ecological restoration, textiles and medicine. A. hendersonii was found to be more tolerant to adverse conditions. The main purpose of this study was to investigate the PP2C gene family and its expression pattern under salt stress and to identify important candidate genes related to salt tolerance.

RESULTS

In this study, 68 AvPP2C genes and 68 AhPP2C genes were identified from the genomes of A. venetum and A. hendersonii, respectively. They were classified into 13 subgroups based on their phylogenetic relationships and were further analyzed for their subcellular locations, gene structures, conserved structural domains, and cis-acting elements. The results of qRT-PCR analyses of seven AvPP2C genes and seven AhPP2C genes proved that they differed significantly in gene expression under salt stress. It has been observed that the PP2C genes in A. venetum and A. hendersonii exhibit different expression patterns. Specifically, AvPP2C2, 6, 24, 27, 41 and AhPP2C2, 6, 24, 27, 42 have shown significant differences in expression under salt stress. This indicates that these genes may play a crucial role in the salt tolerance mechanism of A. venetum and A. hendersonii.

CONCLUSIONS

In this study, we conducted a genome-wide analysis of the AvPP2C and AhPP2C gene families in Apocynum, which provided a reference for further understanding the functional characteristics of these genes.

Understanding the role of P-type ATPases in regulating pollen fertility and development in pigeonpea

The P-type ATPase superfamily genes are the cation and phospholipid pumps that transport ions across the membranes by hydrolyzing ATP. They are involved in a diverse range of functions, including fundamental cellular events that occur during the growth of plants, especially in the reproductive organs. The present work has been undertaken to understand and characterize the P-type ATPases in the pigeonpea genome and their potential role in anther development and pollen fertility. A total of 59 P-type ATPases were predicted in the pigeonpea genome. The phylogenetic analysis classified the ATPases into five subfamilies: eleven P1B, eighteen P2A/B, fourteen P3A, fifteen P4, and one P5. Twenty-three pairs of P-type ATPases were tandemly duplicated, resulting in their expansion in the pigeonpea genome during evolution. The orthologs of the reported anther development-related genes were searched in the pigeonpea genome, and the expression profiling studies of specific genes via qRT-PCR in the pre- and post-meiotic anther stages of AKCMS11A (male sterile), AKCMS11B (maintainer) and AKPR303 (fertility restorer) lines of pigeonpea was done. Compared to the restorer and maintainer lines, the down-regulation of CcP-typeATPase22 in the post-meiotic anthers of the male sterile line might have played a role in pollen sterility. Furthermore, the strong expression of CcP-typeATPase2 in the post-meiotic anthers of restorer line and CcP-typeATPase46, CcP-typeATPase51, and CcP-typeATPase52 in the maintainer lines, respectively, compared to the male sterile line, clearly indicates their potential role in developing male reproductive organs in pigeonpea.

Complete mitochondrial genome of Melia azedarach L., reveals two conformations generated by the repeat sequence mediated recombination

Melia azedarach is a species of enormous value of pharmaceutical industries. Although the chloroplast genome of M. azedarach has been explored, the information of mitochondrial genome (Mt genome) remains surprisingly limited. In this study, we used a hybrid assembly strategy of BGI short-reads and Nanopore long-reads to assemble the Mt genome of M. azedarach. The Mt genome of M. azedarach is characterized by two circular chromosomes with 350,142 bp and 290,387 bp in length, respectively, which encodes 35 protein-coding genes (PCGs), 23 tRNA genes, and 3 rRNA genes. A pair of direct repeats (R1 and R2) were associated with genome recombination, resulting in two conformations based on the Sanger sequencing and Oxford Nanopore sequencing. Comparative analysis identified 19 homologous fragments between Mt and chloroplast genome, with the longest fragment of 12,142 bp. The phylogenetic analysis based on PCGs were consist with the latest classification of the Angiosperm Phylogeny Group. Notably, a total of 356 potential RNA editing sites were predicted based on 35 PCGs, and the editing events lead to the formation of the stop codon in the rps10 gene and the start codons in the nad4L and atp9 genes, which were verified by PCR amplification and Sanger sequencing. Taken together, the exploration of M. azedarach gap-free Mt genome provides a new insight into the evolution research and complex mitogenome architecture.

A reference-grade genome of the xerophyte Ammopiptanthus mongolicus sheds light on its evolution history in legumes and drought-tolerance mechanisms

Plants that grow in extreme environments represent unique sources of stress-resistance genes and mechanisms. Ammopiptanthus mongolicus (Leguminosae) is a xerophytic evergreen broadleaf shrub native to semi-arid and desert regions; however, its drought-tolerance mechanisms remain poorly understood. Here, we report the assembly of a reference-grade genome for A. mongolicus, describe its evolutionary history within the legume family, and examine its drought-tolerance mechanisms. The assembled genome is 843.07 Mb in length, with 98.7% of the sequences successfully anchored to the nine chromosomes of A. mongolicus. The genome is predicted to contain 47 611 protein-coding genes, and 70.71% of the genome is composed of repetitive sequences; these are dominated by transposable elements, particularly long-terminal-repeat retrotransposons. Evolutionary analyses revealed two whole-genome duplication (WGD) events at 130 and 58 million years ago (mya) that are shared by the genus Ammopiptanthus and other legumes, but no species-specific WGDs were found within this genus. Ancestral genome reconstruction revealed that the A. mongolicus genome has undergone fewer rearrangements than other genomes in the legume family, confirming its status as a "relict plant". Transcriptomic analyses demonstrated that genes involved in cuticular wax biosynthesis and transport are highly expressed, both under normal conditions and in response to polyethylene glycol-induced dehydration. Significant induction of genes related to ethylene biosynthesis and signaling was also observed in leaves under dehydration stress, suggesting that enhanced ethylene response and formation of thick waxy cuticles are two major mechanisms of drought tolerance in A. mongolicus. Ectopic expression of AmERF2, an ethylene response factor unique to A. mongolicus, can markedly increase the drought tolerance of transgenic Arabidopsis thaliana plants, demonstrating the potential for application of A. mongolicus genes in crop improvement.

Genome-wide and functional analysis of late embryogenesis abundant (LEA) genes during dormancy and sprouting periods of kernel consumption apricots (P. armeniaca L. × P. sibirica L.)

Late embryogenesis abundant (LEA) proteins play a crucial role in protecting cells from stress, making them potential contributors to abiotic stress tolerance. This study focuses on apricot (P. armeniaca L. × P. sibirica L.), where a comprehensive genome-wide analysis identified 54 LEA genes, categorized into eight subgroups based on phylogenetic relationships. Synteny analysis revealed 14 collinear blocks containing LEA genes between P. armeniaca × P. sibirica and Arabidopsis thaliana, with an additional 9 collinear blocks identified between P. armeniaca × P. sibirica and poplar. Examination of gene structure and conserved motifs indicated that these subgroups exhibit consistent exon-intron patterns and shared motifs. The expansion and duplication of LEA genes in P. armeniaca × P. sibirica were driven by whole-genome duplication (WGD), segmental duplication, and tandem duplication events. Expression analysis, utilizing RNA-seq data and quantitative real-time RT-PCR (qRT-PCR), indicated induction of PasLEA2-20, PasLEA3-2, PasLEA6-1, Pasdehydrin-3, and Pasdehydrin-5 in flower buds during dormancy and sprouting phases. Coexpression network analysis linked LEA genes with 15 cold-resistance genes. Remarkably, during the four developmental stages of flower buds in P. armeniaca × P. sibirica - physiological dormancy, ecological dormancy, sprouting period, and germination stage - the expression patterns of all PasLEAs coexpressed with cold stress-related genes remained consistent. Protein-protein interaction networks, established using Arabidopsis orthologs, emphasized connections between PasLEA proteins and cold resistance pathways. Overexpression of certain LEA genes in yeast and Arabidopsis conferred advantages under cold stress, including increased pod length, reduced bolting time and flowering time, improved survival and seed setting rates, elevated proline accumulation, and enhanced antioxidative enzymatic activities. Furthermore, these overexpressed plants exhibited upregulation of genes related to flower development and cold resistance. The Y1H assay confirmed that PasGBF4 and PasDOF3.5 act as upstream regulatory factors by binding to the promoter region of PasLEA3-2. PasDOF2.4, PasDnaJ2, and PasAP2 were also found to bind to the promoter of Pasdehydrin-3, regulating the expression levels of downstream genes. This comprehensive study explores the evolutionary relationships among PasLEA genes, protein interactions, and functional analyses during various stages of dormancy and sprouting in P. armeniaca × P. sibirica. It offers potential targets for enhancing cold resistance and manipulating flower bud dormancy in this apricot hybrid.

Genome-wide identification and analysis of monocot-specific chimeric jacalins (MCJ) genes in Maize (Zea mays L.)

BACKGROUND

The monocot chimeric jacalins (MCJ) proteins, which contain a jacalin-related lectin (JRL) domain and a dirigent domain (DIR), are specific to Poaceae. MCJ gene family is reported to play an important role in growth, development and stress response. However, their roles in maize have not been thoroughly investigated.

RESULTS

In this study, eight MCJ genes in the maize genome (designated as ZmMCJs) were identified, which displayed unequal distribution across four chromosomes. Phylogenetic relationships between the ZmMCJs were evident through the identification of highly conserved motifs and gene structures. Analysis of transcriptome data revealed distinct expression patterns among the ZmMCJ genes, leading to their classification into four different modules, which were subsequently validated using RT-qPCR. Protein structures of the same module are found to be relatively similar. Subcellular localization experiments indicated that the ZmMCJs are mainly located on the cell membrane. Additionally, hemagglutination and inhibition experiments show that only part of the ZmMCJs protein has lectin activity, which is mediated by the JRL structure, and belongs to the mannose-binding type. The cis-acting elements in the promoter region of ZmMCJ genes predicted their involvement response to phytohormones, such as abscisic acid and jasmonic acid. This suggests that ZmMCJ genes may play a significant role in both biotic and abiotic stress responses.

CONCLUSIONS

Overall, this study adds new insights into our understanding of the gene-protein architecture, evolutionary characteristics, expression profiles, and potential functions of MCJ genes in maize.

Functional and evolutionary comparative analysis of the DIR gene family in Nicotiana tabacum L. and Solanum tuberosum L

BACKGROUND

The dirigent (DIR) genes encode proteins that act as crucial regulators of plant lignin biosynthesis. In Solanaceae species, members of the DIR gene family are intricately related to plant growth and development, playing a key role in responding to various biotic and abiotic stresses. It will be of great application significance to analyze the DIR gene family and expression profile under various pathogen stresses in Solanaceae species.

RESULTS

A total of 57 tobacco NtDIRs and 33 potato StDIRs were identified based on their respective genome sequences. Phylogenetic analysis of DIR genes in tobacco, potato, eggplant and Arabidopsis thaliana revealed three distinct subgroups (DIR-a, DIR-b/d and DIR-e). Gene structure and conserved motif analysis showed that a high degree of conservation in both exon/intron organization and protein motifs among tobacco and potato DIR genes, especially within members of the same subfamily. Total 8 pairs of tandem duplication genes (3 pairs in tobacco, 5 pairs in potato) and 13 pairs of segmental duplication genes (6 pairs in tobacco, 7 pairs in potato) were identified based on the analysis of gene duplication events. Cis-regulatory elements of the DIR promoters participated in hormone response, stress responses, circadian control, endosperm expression, and meristem expression. Transcriptomic data analysis under biotic stress revealed diverse response patterns among DIR gene family members to pathogens, indicating their functional divergence. After 96 h post-inoculation with Ralstonia solanacearum L. (Ras), tobacco seedlings exhibited typical symptoms of tobacco bacterial wilt. The qRT-PCR analysis of 11 selected NtDIR genes displayed differential expression pattern in response to the bacterial pathogen Ras infection. Using line 392278 of potato as material, typical symptoms of potato late blight manifested on the seedling leaves under Phytophthora infestans infection. The qRT-PCR analysis of 5 selected StDIR genes showed up-regulation in response to pathogen infection. Notably, three clustered genes (NtDIR2, NtDIR4, StDIR3) exhibited a robust response to pathogen infection, highlighting their essential roles in disease resistance.

CONCLUSION

The genome-wide identification, evolutionary analysis, and expression profiling of DIR genes in response to various pathogen infection in tobacco and potato have provided valuable insights into the roles of these genes under various stress conditions. Our results could provide a basis for further functional analysis of the DIR gene family under pathogen infection conditions.

Genome-wide identification and expression analysis of calmodulin and calmodulin-like genes in passion fruit (Passiflora edulis) and their involvement in flower and fruit development

BACKGROUND

The calmodulin (CaM) and calmodulin-like (CML) proteins play regulatory roles in plant growth and development, responses to biotic and abiotic stresses, and other biological processes. As a popular fruit and ornamental crop, it is important to explore the regulatory mechanism of flower and fruit development of passion fruit.

RESULTS

In this study, 32 PeCaM/PeCML genes were identified from passion fruit genome and were divided into 9 groups based on phylogenetic analysis. The structural analysis, including conserved motifs, gene structure and homologous modeling, illustrates that the PeCaM/PeCML in the same subgroup have relative conserved structural features. Collinearity analysis suggested that the expansion of the CaM/CML gene family likely took place mainly by segmental duplication, and the whole genome replication events were closely related with the rapid expansion of the gene group. PeCaM/PeCMLs were potentially required for different floral tissues development. Significantly, PeCML26 had extremely high expression levels during ovule and fruit development compared with other PeCML genes, suggesting that PeCML26 had potential functions involved in the development of passion fruit flowers and fruits. The co-presence of various cis-elements associated with growth and development, hormone responsiveness, and stress responsiveness in the promoter regions of these PeCaM/PeCMLs might contribute to their diverse regulatory roles. Furthermore, PeCaM/PeCMLs were also induced by various abiotic stresses. This work provides a comprehensive understanding of the CaM/CML gene family and valuable clues for future studies on the function and evolution of CaM/CML genes in passion fruit.

CONCLUSION

A total of 32 PeCaM/PeCML genes were divided into 9 groups. The PeCaM/PeCML genes showed differential expression patterns in floral tissues at different development stages. It is worth noting that PeCML26, which is highly homologous to AtCaM2, not only interacts with multiple BBR-BPC TFs, but also has high expression levels during ovule and fruit development, suggesting that PeCML26 had potential functions involved in the development of passion fruit flowers and fruits. This research lays the foundation for future investigations and validation of the potential function of PeCaM/PeCML genes in the growth and development of passion fruit.

Genome-wide identification and expression pattern analysis of MIKC-Type MADS-box genes in Chionanthus retusus, an androdioecy plant

BACKGROUND

The MADS-box gene family is widely distributed in the plant kingdom, and its members typically encoding transcription factors to regulate various aspects of plant growth and development. In particular, the MIKC-type MADS-box genes play a crucial role in the determination of floral organ development and identity recognition. As a type of androdioecy plant, Chionanthus retusus have unique gender differentiation. Manifested as male individuals with only male flowers and female individuals with only bisexual flowers. However, due to the lack of reference genome information, the characteristics of MIKC-type MADS-box genes in C. retusus and its role in gender differentiation of C. retusus remain largely unknown. Therefore, it is necessary to identify and characterize the MADS-box gene family within the genome of the C. retusus.

RESULTS

In this study, we performed a genome-wide identification and analysis of MIKC-type MADS-box genes in C. retusus (2n = 2x = 46), utilizing the latest reference genome, and studied its expression pattern in individuals of different genders. As a result, we identified a total of 61 MIKC-type MADS-box genes in C. retusus. 61 MIKC-type MADS-box genes can be divided into 12 subfamilies and distributed on 18 chromosomes. Genome collinearity analysis revealed their conservation in evolution, while gene structure, domains and motif analysis indicated their conservation in structure. Finally, based on their expression patterns in floral organs of different sexes, we have identified that CrMADS45 and CrMADS60 may potentially be involved in the gender differentiation of C. retusus.

CONCLUSIONS

Our studies have provided a general understanding of the conservation and characteristics of the MIKC-type MADS-box genes family in C. retusus. And it has been demonstrated that members of the AG subfamily, CrMADS45 and CrMADS60, may play important roles in the gender differentiation of C. retusus. This provides a reference for future breeding efforts to improve flower types in C. retusus and further investigate the role of MIKC-type MADS-box genes in gender differentiation.

Genome-wide identification and expression analysis of histone deacetylase and histone acetyltransferase genes in response to drought in poplars

BACKGROUND

Histone deacetylases (HDACs) and histone acetyltransferases (HATs) are involved in plant growth and development as well as in response to environmental changes, by dynamically regulating gene acetylation levels. Although there have been numerous reports on the identification and function of HDAC and HAT in herbaceous plants, there are fewer report related genes in woody plants under drought stress.

RESULTS

In this study, we performed a genome-wide analysis of the HDAC and HAT families in Populus trichocarpa, including phylogenetic analysis, gene structure, conserved domains, and expression analysis. A total of 16 PtrHDACs and 12 PtrHATs were identified in P. trichocarpa genome. Analysis of cis-elements in the promoters of PtrHDACs and PtrHATs revealed that both gene families could respond to a variety of environmental signals, including hormones and drought. Furthermore, real time quantitative PCR indicated that PtrHDA906 and PtrHAG3 were significantly responsive to drought. PtrHDA906, PtrHAC1, PtrHAC3, PtrHAG2, PtrHAG6 and PtrHAF1 consistently responded to abscisic acid, methyl jasmonate and salicylic acid under drought conditions.

CONCLUSIONS

Our study demonstrates that PtrHDACs and PtrHATs may respond to drought through hormone signaling pathways, which helps to reveal the hub of acetylation modification in hormone regulation of abiotic stress.

Crape myrtle LiGAoxs displaying activities of gibberellin oxidases respond to branching architecture

In the realm of ornamental horticulture, crape myrtle (Lagerstroemia indica) stands out for its aesthetic appeal, attributed largely to its vibrant flowers and distinctive branching architecture. This study embarked on a comprehensive exploration of the gibberellin oxidase (GAox) gene family in crape myrtle, illuminating its pivotal role in regulating GA levels, a key determinant of plant developmental processes. We identified and characterized 36 LiGAox genes, subdivided into GA2ox, GA3ox, GA20ox, and GAox-like subgroups, through genomic analyses. These genes' evolutionary trajectories were delineated, revealing significant gene expansions attributed to segmental duplication events. Functional analyses highlighted the divergent expression patterns of LiGAox genes across different crape myrtle varieties, associating them with variations in flower color and branching architecture. Enzymatic activity assays on selected LiGA2ox enzymes exhibited pronounced GA2 oxidase activity, suggesting a potential regulatory role in GA biosynthesis. Our findings offered a novel insight into the molecular underpinnings of GA-mediated growth and development in L. indica, providing a foundational framework for future genetic enhancements aimed at optimizing ornamental traits.

Genome-wide characterization and functional analysis of the melon TGA gene family in disease resistance through ectopic overexpression in Arabidopsis thaliana

TGA-binding (TGA) transcription factors, characterized by the basic region/leucine zipper motif (bZIP), have been recognized as pivotal regulators in plant growth, development, and stress responses through their binding to the as-1 element. In this study, the TGA gene families in melon, watermelon, cucumber, pumpkin, and zucchini were comprehensively characterized, encompassing analyses of gene/protein structures, phylogenetic relationships, gene duplication events, and cis-acting elements in gene promoters. Upon transient expression in Nicotiana benthamiana, the melon CmTGAs, with typical bZIP and DOG1 domains, were observed to localize within the nucleus. Biochemical investigation revealed specific interactions between CmTGA2/3/5/8/9 and CmNPR3 or CmNPR4. The CmTGA genes exhibited differential expression patterns in melon plants in response to different hormones like salicylic acid, methyl jasmonate, and ethylene, as well as a fungal pathogen, Stagonosporopsis cucurbitacearum that causes gummy stem blight in melon. The overexpression of CmTGA3, CmTGA8, and CmTGA9 in Arabidopsis plants resulted in the upregulation of AtPR1 and AtPR5 expression, thereby imparting enhanced resistance to Pseudomonas syringae pv. Tomato DC3000. In contrast, the overexpression of CmTGA7 or CmTGA9 resulted in a compromised resistance to Botrytis cinerea, coinciding with a concomitant reduction in the expression levels of AtPDF1.2 and AtMYC2 following infection with B. cinerea. These findings shed light on the important roles of specific CmTGA genes in plant immunity, suggesting that genetic manipulation of these genes could be a promising avenue for enhancing plant immune responses.

Detection of colinear blocks and synteny and evolutionary analyses based on utilization of MCScanX

As different taxa evolve, gene order often changes slowly enough that chromosomal 'blocks' with conserved gene orders (synteny) are discernible. The MCScanX toolkit ( https://github.com/wyp1125/MCScanX ) was published in 2012 as freely available software for the detection of such 'colinear blocks' and subsequent synteny and evolutionary analyses based on genome-wide gene location and protein sequence information. Owing to its simplicity and high efficiency for colinear block detection, MCScanX provides a powerful tool for conducting diverse synteny and evolutionary analyses. Moreover, the detection of colinear blocks has been embraced as an integral step for pangenome graph construction. Here, new application trends of MCScanX are explored, striving to better connect this increasingly used tool to other tools and accelerate insight generation from exponentially growing sequence data. We provide a detailed protocol that covers how to install MCScanX on diverse platforms, tune parameters, prepare input files from data from the National Center for Biotechnology Information, run MCScanX and its visualization and evolutionary analysis tools, and connect MCScanX with external tools, including MCScanX-transposed, Circos and SynVisio. This protocol is easily implemented by users with minimal computational background and is adaptable to new data of interest to them. The data and utility programs for this protocol can be obtained from http://bdx-consulting.com/mcscanx-protocol .

Allopolyploid subgenome identification and implications for evolutionary analysis

Whole-genome duplications (WGDs) are widespread genomic events in eukaryotes that are hypothesized to contribute to the evolutionary success of many lineages, including flowering plants, Saccharomyces yeast, and vertebrates. WGDs generally can be classified into autopolyploids (ploidy increase descended from one species) or allopolyploids (ploidy increase descended from multiple species). Assignment of allopolyploid progenitor species (called subgenomes in the polyploid) is important to understanding the biology and evolution of polyploids, including the asymmetric subgenome evolution following hybridization (biased fractionation). Here, I review the different methodologies used to identify the ancestors of allopolyploid subgenomes, discuss the advantages and disadvantages of these methods, and outline the implications of how these methods affect the subsequent evolutionary analysis of these genomes.

W chromosome sequences of two bombycid moths provide an insight into the origin of Fem

Fem is a W-linked gene that encodes a piRNA precursor, and its product, Fem piRNA, is a master factor of female determination in Bombyx mori. Fem has low similarity to any known sequences, and the origin of Fem remains unclear. So far, two hypotheses have been proposed for the origin of Fem: The first hypothesis is that Fem is an allele of Masc, which assumes that the W chromosome was originally a homologous chromosome of the Z chromosome. The second hypothesis is that Fem arose by the transposition of Masc to the W chromosome. To explore the origin of Fem, we determined the W chromosome sequences of B. mori and, as a comparison, a closely relative bombycid species of Trilocha varians with a Fem-independent sex determination system. To our surprise, although the sequences of W and Z chromosomes show no homology to each other, a few pairs of homologues are shared by W and Z chromosomes, indicating the W chromosome of both species originated from Z chromosome. In addition, the W chromosome of T. varians lacks Fem, while the W chromosome of B. mori has over 100 copies of Fem. The high-quality assembly of the W chromosome of B. mori arose the third hypothesis about the origin of Fem: Fem is a chimeric sequence of multiple transposons. More than half of one transcriptional unit of Fem shows a significant homology to RTE-BovB. Moreover, the Fem piRNA-producing region could correspond to the boundary of the two transposons, gypsy and satellite DNA.

Functional characterization of geranyl/farnesyl diphosphate synthase in Wurfbainia villosa and Wurfbainia longiligularis

Wurfbainia villosa and Wurfbainia longiligularis are the two primary plant sources of Fructus Amomi, a traditional Chinese medicine. Both plants are rich in volatile terpenoids, including monoterpenes and sesquiterpenes, which are the primary medicinal components of Fructus Amomi. The trans-isopentenyl diphosphate synthase (TIDS) gene family plays a key part in determining terpenoid diversity and accumulation. However, the TIDS gene family have not been identified in W. villosa and W. longiligularis. This study identified thirteen TIDS genes in W. villosa and eleven TIDS genes in W. longiligularis, which may have expanded through segmental replication events. Based on phylogenetic analysis and expression levels, eight candidate WvTIDSs and five WlTIDSs were selected for cloning. Functional characterization in vitro demonstrated that four homologous geranyl diphosphate synthases (GPPSs) (WvGPPS1, WvGPPS2, WlGPPS1, WlGPPS2) and two geranylgeranyl diphosphate synthases (GGPPSs) (WvGGPPS and WlGGPPS) were responsible for catalyzing the biosynthesis of geranyl diphosphate (GPP), whereas two farnesyl diphosphate synthases (FPPSs) (WvFPPS and WlFPPS) catalysed the biosynthesis of the farnesyl diphosphate (FPP). A comparison of six proteins with identified GPPS functions showed that WvGGPPS and WlGGPPS exhibited the highest activity levels. These findings indicate that homologous GPPS and GGPPS together promote the biosynthesis of GPP in W. villosa and W. longiligularis, thus providing sufficient precursors for the synthesis of monoterpenes and providing key genetic elements for Fructus Amomi variety improvement and molecular breeding.

Comprehensive identification and analysis of DUF640 genes associated with rice growth

Protein domains with conserved amino acid sequences and uncharacterized functions are called domains of unknown function (DUF). The DUF640 gene family plays a crucial role in plant growth, particularly in light regulation, floral organ development, and fruit development. However, there exists a lack of systematic understanding of the evolutionary relationships and functional differentiation of DUF640 within the Oryza genus. In this study, 61 DUF640 genes were identified in the Oryza genus. The expression of DUF640s is induced by multiple hormonal stressors including abscisic acid (ABA), cytokinin (CK), ethylene (ETH), and indole-3-acetic acid (IAA). Specifically, OiDUF640-10 expression significantly increased after ETH treatment. Transgenic experiments showed that overexpressing OiDUF640-10 lines were sensitive to ETH, and seedling length was obstructed. Evolutionary analysis revealed differentiation of the OiDUF640-10 gene in O. sativa ssp. indica and japonica varieties, likely driven by natural selection during the domestication of cultivated rice. These results indicate that OiDUF640-10 plays a vital role in the regulation of rice seedling length.

Haplotype-resolved genomes of octoploid species in Phyllanthaceae family reveal a critical role for polyploidization and hybridization in speciation

The Phyllanthaceae family comprises a diverse range of plants with medicinal, edible, and ornamental value, extensively cultivated worldwide. Polyploid species commonly occur in Phyllanthaceae. Due to the rather complex genomes and evolutionary histories, their speciation process has been still lacking in research. In this study, we generated chromosome-scale haplotype-resolved genomes of two octoploid species (Phyllanthus emblica and Sauropus spatulifolius) in Phyllanthaceae family. Combined with our previously reported one tetraploid (Sauropus androgynus) and one diploid species (Phyllanthus cochinchinensis) from the same family, we explored their speciation history. The three polyploid species were all identified as allopolyploids with subgenome A/B. Each of their two distinct subgenome groups from various species was uncovered to independently share a common diploid ancestor (Ancestor-AA and Ancestor-BB). Via different evolutionary routes, comprising various scenarios of bifurcating divergence, allopolyploidization (hybrid polyploidization), and autopolyploidization, they finally evolved to the current tetraploid S. androgynus, and octoploid S. spatulifolius and P. emblica, respectively. We further discuss the variations in copy number of alleles and the potential impacts within the two octoploids. In addition, we also investigated the fluctuation of metabolites with medical values and identified the key factor in its biosynthesis process in octoploids species. Our study reconstructed the evolutionary history of these Phyllanthaceae species, highlighting the critical roles of polyploidization and hybridization in their speciation processes. The high-quality genomes of the two octoploid species provide valuable genomic resources for further research of evolution and functional genomics.

Genome-wide identification of the LEA gene family in Panax ginseng: Evidence for the role of PgLEA2-50 in plant abiotic stress response

Ginseng frequently encounters environmental stress during its growth and development. Late Embryogenesis Abundant (LEA) proteins play a crucial role in combating adversity stress, particularly against abiotic challenges In this study, 107 LEA genes from ginseng, spanning eight subfamilies, were identified, demonstrating significant evolutionary conservation, with the LEA2 subfamily being most prominent. Gene duplication events, primarily segmental duplications, have played a major role in the expansion of the LEA gene family, which has undergone strong purifying selection. PgLEAs were unevenly distributed across 22 chromosomes, with each subfamily featuring unique structural domains and conserved motifs. PgLEAs were expressed in various tissues, exhibiting distinct variations in abundance and tissue specificity. Numerous regulatory cis-elements, related to abiotic stress and hormones, were identified in the promoter region. Additionally, PgLEAs were regulated by a diverse array of abiotic stress-related transcription factors. A total of 35 PgLEAs were differentially expressed following treatments with ABA, GA, and IAA. Twenty-three PgLEAs showed significant but varied responses to drought, extreme temperatures, and salinity stress. The transformation of tobacco with the key gene PgLEA2-50 enhanced osmoregulation and antioxidant levels in transgenic lines, improving their resistance to abiotic stress. This study offers insights into functional gene analysis, focusing on LEA proteins, and establishes a foundational framework for research on ginseng's resilience to abiotic stress.

Integrated analysis of genomics and transcriptomics revealed the genetic basis for goaty flavor formation in goat milk

Goat milk exhibits a robust and distinctive "goaty" flavor. However, the underlying genetic basis of goaty flavor remains elusive and requires further elucidation at the genomic level. Through comparative genomics analysis, we identified divergent signatures of certain proteins in goat, sheep, and cow. MMUT has undergone a goat-specific mutation in the B12 binding domain. We observed the goat FASN exhibits nonsynonymous mutations in the acyltransferase domain. Structural variations in these key proteins may enhance the capacity for synthesizing goaty flavor compounds in goat. Integrated omics analysis revealed the catabolism of branched-chain amino acids contributed to the goat milk flavor. Furthermore, we uncovered a regulatory mechanism in which the transcription factor ZNF281 suppresses the expression of the ECHDC1 gene may play a pivotal role in the accumulation of flavor substances in goat milk. These findings provide insights into the genetic basis underlying the formation of goaty flavor in goat milk. STATEMENT OF SIGNIFICANCE: Branched-chain fatty acids (BCFAs) play a crucial role in generating the distinctive "goaty" flavor of goat milk. Whether there is an underlying genetic basis associated with goaty flavor is unknown. To begin deciphering mechanisms of goat milk flavor development, we collected transcriptomic data from mammary tissue of goat, sheep, cow, and buffalo at peak lactation for cross-species transcriptome analysis and downloaded nine publicly available genomes for comparative genomic analysis. Our data indicate that the catabolic pathway of branched-chain amino acids (BCAAs) is under positive selection in the goat genome, and most genes involved in this pathway exhibit significantly higher expression levels in goat mammary tissue compared to other species, which contributes to the development of flavor in goat milk. Furthermore, we have elucidated the regulatory mechanism by which the transcription factor ZNF281 suppresses ECHDC1 gene expression, thereby exerting an important influence on the accumulation of flavor compounds in goat milk. These findings provide insights into the genetic mechanisms underlying flavor formation in goat milk and suggest further research to manipulate the flavor of animal products.

Genome-wide identification of walnut (Juglans regia) PME gene family members and expression analysis during infection with Cryptosphaeria pullmanensis pathogens

Walnuts exhibit a higher resistance to diseases, though they are not completely immune. This study focuses on the Pectin methylesterase (PME) gene family to investigate whether it is involved in disease resistance in walnuts. These 21 genes are distributed across 12 chromosomes, with four pairs demonstrating homology. Variations in conserved motifs and gene structures suggest diverse functions within the gene family. Phylogenetic and collinear gene pairs of the PME family indicate that the gene family has evolved in a relatively stable way. The cis-acting elements and gene ontology enrichment of these genes, underscores their potential role in bolstering walnuts' defense mechanisms. Transcriptomic analyses were conducted under conditions of Cryptosphaeria pullmanensis infestation and verified by RT-qPCR. The results showed that certain JrPME family genes were activated in response, leading to the hypothesis that some members may confer resistance to the disease.

Genome assembly and resequencing shed light on evolution, population selection, and sex identification in Vernicia montana

Vernicia montana is a dioecious plant widely cultivated for high-quality tung oil production and ornamental purposes in the Euphorbiaceae family. The lack of genomic information has severely hindered molecular breeding for genetic improvement and early sex identification in V. montana. Here, we present a chromosome-level reference genome of a male V. montana with a total size of 1.29 Gb and a contig N50 of 3.69 Mb. Genome analysis revealed that different repeat lineages drove the expansion of genome size. The model of chromosome evolution in the Euphorbiaceae family suggests that polyploidization-induced genomic structural variation reshaped the chromosome structure, giving rise to the diverse modern chromosomes. Based on whole-genome resequencing data and analyses of selective sweep and genetic diversity, several genes associated with stress resistance and flavonoid synthesis such as CYP450 genes and members of the LRR-RLK family, were identified and presumed to have been selected during the evolutionary process. Genome-wide association studies were conducted and a putative sex-linked insertion and deletion (InDel) (Chr 2: 102 799 917-102 799 933 bp) was identified and developed as a polymorphic molecular marker capable of effectively detecting the gender of V. montana. This InDel is located in the second intron of VmBASS4, suggesting a possible role of VmBASS4 in sex determination in V. montana. This study sheds light on the genome evolution and sex identification of V. montana, which will facilitate research on the development of agronomically important traits and genomics-assisted breeding.

Genome-scale identification, expression and evolution analysis of B-box members in Dendrobium huoshanense

B-box (BBX) proteins have been recognized as vital determinants in plant development, morphogenesis, and adaptive responses to a myriad of environmental stresses. These zinc-finger proteins play a pivotal role in various biological processes. Their influence spans photomorphogenesis, the regulation of flowering, and imparting resilience to a wide array of challenges, encompassing both biotic and abiotic factors. Chromosome localization, gene structure and conserved motifs, phylogenetic analysis, collinearity analysis, expression profiling, fluorescence quantitative analysis, and tobacco transient transformation methods were used for functional localization and expression pattern analysis of the DhBBX gene. A total of 23 DhBBX members were identified from Dendrobium huoshanense. Subsequent phylogenetic evaluations effectively segregated these genes into five discrete evolutionary subsets. The predictions of subcellular localizations revealed that all these proteins were localized in the nucleus. The genetic composition and patterns showed that the majority of these genes consisted of several exons, with a few variations that could be attributed to transposon insertion. A comprehensive analysis using qRT-PCR was conducted to unravel the expression patterns of these genes in D. huoshanense, with a specific concentration on their responses to various hormone treatments and cold stress. Subcellular localization reveals that DhBBX21 and DhBBX9 are located in the nucleus. Our results provide a deep comprehension of the complex regulatory mechanisms of BBXs in response to various environmental and hormonal stimuli. These discoveries encourage further detailed and focused investigations into the operational dynamics of the BBX gene family in a wider range of plant species.

Genome sequencing of Coryphaenoides yaquinae reveals convergent and lineage-specific molecular evolution in deep-sea adaptation

Abyssal (3501-6500 m) and hadal (>6500 m) fauna evolve under harsh abiotic stresses, characterized by high hydrostatic pressure, darkness and food shortage, providing unique opportunities to investigate mechanisms underlying environmental adaptation. Genomes of several hadal species have recently been reported. However, the genetic adaptation of deep sea species across a broad spectrum of ocean depths has yet to be thoroughly investigated, due to the challenges imposed by collecting the deep sea species. To elucidate the correlation between genetic innovation and vertical distribution, we generated a chromosome-level genome assembly of the macrourids Coryphaenoides yaquinae, which is widely distributed in the abyssal/hadal zone ranging from 3655 to 7259 m in depth. Genomic comparisons among shallow, abyssal and hadal-living species identified idiosyncratic and convergent genetic alterations underlying the extraordinary adaptations of deep-sea species including light perception, circadian regulation, hydrostatic pressure and hunger tolerance. The deep-sea fishes (Coryphaenoides Sp. and Pseudoliparis swirei) venturing into various ocean depths independently have undergone convergent amino acid substitutions in multiple proteins such as rhodopsin 1, pancreatic and duodenal homeobox 1 and melanocortin 4 receptor which are known or verified in zebrafish to be related with vision adaptation and energy expenditure. Convergent evolution events were also identified in heat shock protein 90 beta family member 1 and valosin-containing protein genes known to be related to hydrostatic pressure adaptation specifically in fishes found around the hadal range. The uncovering of the molecular convergence among the deep-sea species shed new light on the common genetic innovations required for deep-sea adaptation by the fishes.

Genome-wide identification, expression analysis of WRKY transcription factors in Citrus ichangensis and functional validation of CiWRKY31 in response to cold stress

BACKGROUND

Ichang papeda (Citrus ichangensis), a wild perennial plant of the Rutaceae family, is a cold-hardy plant. WRKY transcription factors are crucial regulators of plant growth and development as well as abiotic stress responses. However, the WRKY genes in C. ichangensis (CiWRKY) and their expression patterns under cold stress have not been thoroughly investigated, hindering our understanding of their role in cold tolerance.

RESULTS

In this study, a total of 52 CiWRKY genes identified in the genome of C. ichangensis were classified into three main groups and five subgroups based on phylogenetic analysis. Comprehensive analyses of motif features, conserved domains, and gene structures were performed. Segmental duplication plays a significant role in the CiWRKY gene family expansion. Cis-acting element analysis revealed the presence of various stress-responsive elements in the promoters of the majority of CiWRKYs. Gene ontology (GO) analysis and protein-protein interaction predictions indicate that the CiWRKYs exhibit crucial roles in regulation of both development and stress response. Expression profiling analysis demonstrates that 14 CiWRKYs were substantially induced under cold stress. Virus-induced gene silencing (VIGS) assay confirmed that CiWRKY31, one of the cold-induced WRKYs, functions positively in regulation of cold tolerance.

CONCLUSION

Sequence and protein properties of CiWRKYs were systematically analyzed. Among the 52 CiWRKY genes 14 members exhibited cold-responsive expression patterns, and CiWRKY31 was verified to be a positive regulator of cold tolerance. These findings pave way for future investigations to understand the molecular functions of CiWRKYs in cold tolerance and contribute to unravelling WRKYs that may be used for engineering cold tolerance in citrus.

Comprehensive identification and expression analysis of FAR1/FHY3 genes under drought stress in maize (Zea mays L.)

Background

FAR1/FHY3 transcription factors are derived from transposase, which play important roles in light signal transduction, growth and development, and response to stress by regulating downstream gene expression. Although many FAR1/FHY3 members have been identified in various species, the FAR1/FHY3 genes in maize are not well characterized and their function in drought are unknown.

Method

The FAR1/FHY3 family in the maize genome was identified using PlantTFDB, Pfam, Smart, and NCBI-CDD websites. In order to investigate the evolution and functions of FAR1 genes in maize, the information of protein sequences, chromosome localization, subcellular localization, conserved motifs, evolutionary relationships and tissue expression patterns were analyzed by bioinformatics, and the expression patterns under drought stress were detected by quantitative real-time polymerase chain reaction (qRT-PCR).

Results

A total of 24 ZmFAR members in maize genome, which can be divided into five subfamilies, with large differences in protein and gene structures among subfamilies. The promoter regions of ZmFARs contain abundant abiotic stress-responsive and hormone-respovensive cis-elements. Among them, drought-responsive cis-elements are quite abundant. ZmFARs were expressed in all tissues detected, but the expression level varies widely. The expression of ZmFARs were mostly down-regulated in primary roots, seminal roots, lateral roots, and mesocotyls under water deficit. Most ZmFARs were down-regulated in root after PEG-simulated drought stress.

Conclusions

We performed a genome-wide and systematic identification of FAR1/FHY3 genes in maize. And most ZmFARs were down-regulated in root after drought stress. These results indicate that FAR1/FHY3 transcription factors have important roles in drought stress response, which can lay a foundation for further analysis of the functions of ZmFARs in response to drought stress.

Chromosome-level genome assembly and comparative genomics shed light on Helicoverpa assulta ecology and pest management

BACKGROUND

The Oriental tobacco budworm, Helicoverpa assulta, a specialist herbivorous insect that exclusively feeds on plants of the Solanaceae family, causes considerable damage to crops, such as tobacco and hot pepper. The absence of a genome sequence for this species hinders further research on its pest management and ecological adaptation.

RESULTS

Here, we present a high-quality chromosome-level genome of a Korean strain of H. assulta (Pyeongchang strain, K18). The total assembly spans 424.4 Mb with an N50 length of 14.54 Mb and 37% GC content. The assembled genome (ASM2961881v1) comprises 31 chromosomes, similar to other congeneric generalist species including H. armigera and H. zea. In terms of genomic assembly quality, the complete BUSCOs and repeat content accounted for 98.3% and 33.01% of the genome, respectively. Based on this assembly, 19 485 protein-coding genes were predicted in the genome annotation. A comparative analysis was conducted using the identified number of protein-coding genes in H. armigera (24154) and H. zea (23696). Out of the 19 485 predicted genes, 137 genes in 15 orthogroups were found to have expanded significantly in H. assulta, while 149 genes in 95 orthogroups contracted rapidly.

CONCLUSION

This study revealed specific gene expansions and contractions in H. assulta compared to those in its close relatives, indicating potential adaptations related to its specialized feeding habits. Also, the comparative genome analysis provides valuable insights for the integrated pest management of H. assulta and other globally significant pests in the Heliothinae subfamily. © 2024 Society of Chemical Industry.

Genome-wide identification of Rosa roxburghii CML family genes identifies an RrCML13-RrGGP2 interaction involved in calcium-mediated regulation of ascorbate biosynthesis

Calmodulin-like proteins (CMLs) are an essential family of calcium sensors involved in multiple Ca2+-mediated cellular processes in plants. Rosa roxburghii Tratt, known for the abundance of L-ascorbic acid (AsA) in its fruits, is widely distributed in calcium-rich soil of the karst region in southwestern China. The aim of this study was to identify key CMLs that respond to exogenous Ca2+ levels and regulate AsA biosynthesis in R. roxburghii. A genome-wide scan revealed the presence of 41 RrCML genes with 1-4 EF-hand motif (s) unevenly distributed across the 7 chromosomes of R. roxburghii. qRT-PCR analysis revealed that RrCML13, RrCML10, and RrCML36 responded significantly to exogenous Ca2+ treatment, and RrCML13 was positively correlated with GDP-L-galactose phosphorylase encoding gene (RrGGP2) expression and AsA content in the developing fruit. Overexpression of RrCML13 in fruits and roots significantly promoted the transcription of RrGGP2 and the accumulation of AsA, while virus-induced silencing of RrCML13 reduced the transcription of RrGGP2 and the content of AsA. Furthermore, Moreover, the yeast two-hybrid and bimolecular fluorescence complementation (BiFC) analysis confirmed the interaction between RrCML13 and RrGGP2 proteins, indicating that RrCML13 plays a regulatory role in calcium-mediated AsA biosynthesis. This study enhances our understanding of R. roxburghii CMLs and sheds light on the calcium-mediated regulation of AsA biosynthesis.

Deciphering recursive polyploidization in Lamiales and reconstructing their chromosome evolutionary trajectories

Lamiales is an order of core eudicots with abundant diversity, and many Lamiales plants have important medicinal and ornamental values. Here, we comparatively reanalyzed 11 Lamiales species with well-assembled genome sequences and found evidence that Lamiales plants, in addition to a hexaploidization or whole-genome triplication (WGT) shared by core eudicots, experienced further polyploidization events, establishing new groups in the order. Notably, we identified a whole-genome duplication (WGD) occurred just before the split of Scrophulariaceae from the other Lamiales families, such as Acanthaceae, Bignoniaceae, and Lamiaceae, suggesting its likely being the causal reason for the establishment and fast divergence of these families. We also found that a WGT occurred ∼68 to 78 million years ago (Mya), near the split of Oleaceae from the other Lamiales families, implying that it may have caused their fast divergence and the establishment of the Oleaceae family. Then, by exploring and distinguishing intra- and intergenomic chromosomal homology due to recursive polyploidization and speciation, respectively, we inferred that the Lamiales ancestral cell karyotype had 11 proto-chromosomes. We reconstructed the evolutionary trajectories from these proto-chromosomes to form the extant chromosomes in each Lamiales plant under study. We must note that most of the inferred 11 proto-chromosomes, duplicated during a WGD thereafter, have been well preserved in jacaranda (Jacaranda mimosifolia) genome, showing the credibility of the present inference implementing a telomere-centric chromosome repatterning model. These efforts are important to understand genome repatterning after recursive polyploidization, especially shedding light on the origin of new plant groups and angiosperm cell karyotype evolution.

The characteristics of PtHSP40 gene family in Phaeodactylum tricornutum and its response to environmental stresses

Diatom has evolved response mechanisms to cope with multiple environmental stresses. Heat shock protein 40 (HSP40) plays a key role in these response mechanisms. HSP40 gene family in higher plants has been well-studied. However, the HSP40 gene family has not been systematically investigated in marine diatom. In this study, the bioinformatic characteristics, phylogenetic relationship, conserved motifs, gene structure, chromosome distribution and the transcriptional response of PtHSP40 to different environmental stresses were analyzed in the diatom Phaeodactylum tricornutum, and quantitative real-time PCR was conducted. Totally, 55 putative PtHSP40 genes are distributed to 21 chromosomes. All PtHSP40 proteins can be divided into four groups based on their evolutionary relationship, and 54 of them contain a conserved HPD (histidine-proline-aspartic acid tripeptide) motif. Additionally, six, eleven, ten and four PtHSP40 genes were significantly upregulated under the treatments of nitrogen starvation, phosphorus deprivation, 2,2',4,4'-tetrabrominated biphenyl ether (BDE-47) and ocean acidification, respectively. More interestingly, the expression level of 9 PtHSP40 genes was obviously upregulated in response to nickel stress, suggesting the sensitive to metal stress. The different expression models of PtHSP40 genes to environmental stresses imply the specificity of PtHSP40 proteins under different stresses. This study provides a systematic understanding of the PtHSP40 gene family in P. tricornutum and a comprehensive cognition in its functions and response mechanisms to environmental stresses.

YABBY transcription factor family in the octoploid Fragaria × ananassa and five diploid Fragaria species

YABBY genes encode specific TFs of seed plants involved in development and formation of leaves, flowers, and fruit. In the present work, genome-wide and expression analyses of the YABBY gene family were performed in six species of the Fragaria genus: Fragaria × ananassa, F. daltoniana, F. nilgerrensis, F. pentaphylla, F. viridis, and F. vesca. The chromosomal location, synteny pattern, gene structure, and phylogenetic analyses were carried out. By combining RNA-seq data and RT-qPCR analysis we explored specific expression of YABBYs in F. × ananassa and F. vesca. We also analysed the promoter regions of FaYABBYs and performed MeJA application to F. × ananassa fruit to observe effects on gene expression. We identified and characterized 25 YABBY genes in F. × ananassa and six in each of the other five species, which belong to FIL/YAB3 (YABBY1), YAB2 (YABBY2), YAB5 (YABBY5), CRC, and INO clades previously described. Division of the YABBY1 clade into YABBY1.1 and YABBY1.2 subclades is reported. We observed differential expression according to tissue, where some FaYABBYs are expressed mainly in leaves and flowers and to a minor extent during fruit development of F. × ananassa. Specifically, the FaINO genes contain jasmonate-responsive cis-acting elements in their promoters which may be functional since FaINOs are upregulated in F. × ananassa fruit under MeJA treatment. This study suggests that YABBY TFs play an important role in the development- and environment-associated responses of the Fragaria genus.

QTL mapping of downy mildew resistance in foxtail millet by SLAF‑seq and BSR-seq analysis

KEY MESSAGE

Key message Three major QTLs for resistance to downy mildew were located within an 0.78 Mb interval on chromosome 8 in foxtail millet. Downy mildew, a disease caused by Sclerospora graminicola, is a serious problem that jeopardizes the yield and quality of foxtail millet. Breeding resistant varieties represents one of the most economical and effective solutions, yet there is a lack of molecular markers related to the resistance. Here, a mapping population comprising of 158 F6:7 recombinant inbred lines (RILs) was constructed from the crossing of G1 and JG21. Based on the specific locus amplified fragment sequencing results, a high-density linkage map of foxtail millet with 1031 bin markers, spanning 1041.66 cM was constructed. Based on the high-density linkage map and the phenotype data in four environments, a total of nine quantitative trait loci (QTL) associated with resistance to downy mildew were identified. Further BSR-seq confirmed the genomic regions containing the potential candidate genes related to downy mildew resistance. Interestingly, a 0.78-Mb interval between C8M257 and C8M268 on chromosome 8 was highlighted because of its presence in three major QTL, qDM8_1, qDM8_2, and qDM8_4, which contains 10 NBS-LRR genes. Haplotype analysis in RILs and natural population suggest that 9 SNP loci on Seita8G.199800, Seita8G.195900, Seita8G.198300, and Seita.8G199300 genes were significantly correlated with disease resistance. Furthermore, we found that those genes were taxon-specific by collinearity analysis of pearl millet and foxtail millet genomes. The identification of these new resistance QTL and the prediction of resistance genes against downy mildew will be useful in breeding for resistant varieties and the study of genetic mechanisms of downy mildew disease resistance in foxtail millet.

Vole genomics links determinate and indeterminate growth of teeth

Continuously growing teeth are an important innovation in mammalian evolution, yet genetic regulation of continuous growth by stem cells remains incompletely understood. Dental stem cells responsible for tooth crown growth are lost at the onset of tooth root formation. Genetic signaling that initiates this loss is difficult to study with the ever-growing incisor and rooted molars of mice, the most common mammalian dental model species, because signals for root formation overlap with signals that pattern tooth size and shape (i.e., cusp patterns). Different species of voles (Cricetidae, Rodentia, Glires) have evolved rooted and unrooted molars that have similar size and shape, providing alternative models for studying roots. We assembled a de novo genome of Myodes glareolus, a vole with high-crowned, rooted molars, and performed genomic and transcriptomic analyses in a broad phylogenetic context of Glires (rodents and lagomorphs) to assess differential selection and evolution in tooth forming genes. We identified 15 dental genes with changing synteny relationships and six dental genes undergoing positive selection across Glires, two of which were undergoing positive selection in species with unrooted molars, Dspp and Aqp1. Decreased expression of both genes in prairie voles with unrooted molars compared to bank voles supports the presence of positive selection and may underlie differences in root formation. Bulk transcriptomics analyses of embryonic molar development in bank voles also demonstrated conserved patterns of dental gene expression compared to mice, with species-specific variation likely related to developmental timing and morphological differences between mouse and vole molars. Our results support ongoing evolution of dental genes across Glires, revealing the complex evolutionary background of convergent evolution for ever-growing molars.

Convergent evolution of desiccation tolerance in grasses

Desiccation tolerance has evolved repeatedly in plants as an adaptation to survive extreme environments. Plants use similar biophysical and cellular mechanisms to survive life without water, but convergence at the molecular, gene and regulatory levels remains to be tested. Here we explore the evolutionary mechanisms underlying the recurrent evolution of desiccation tolerance across grasses. We observed substantial convergence in gene duplication and expression patterns associated with desiccation. Syntenic genes of shared origin are activated across species, indicative of parallel evolution. In other cases, similar metabolic pathways are induced but using different gene sets, pointing towards phenotypic convergence. Species-specific mechanisms supplement these shared core mechanisms, underlining the complexity and diversity of evolutionary adaptations to drought. Our findings provide insight into the evolutionary processes driving desiccation tolerance and highlight the roles of parallel and convergent evolution in response to environmental challenges.

Genome-wide identification and expression analysis of the GAD family reveal their involved in embryogenesis and hormones responses in Dimocarpus longan Lour

Glutamate decarboxylase (GAD) is involved in GABA metabolism and plays an essential regulatory role in plant growth, abiotic stresses, and hormone response. This study investigated the expression mechanism of the GAD family during longan early somatic embryogenesis (SE) and identified 6 GAD genes based on the longan genome. Homology analysis indicated that DlGAD genes had a closer relationship with dicotyledonous plants. The analysis of cis-acting elements in the promoter region suggests that the GAD genes were associated with various stress responses and hormones. RNA sequencing (RNA-Seq) and the qRT-PCR data indicated that most DlGAD genes were highly expressed in the incomplete compact pro-embryogenic cultures (ICpEC) and upregulated in longan embryogenic callus (EC) after treatments with 2,4-D, high temperature (35 °C), IAA, and ABA. Moreover, the RNA-Seq analysis also revealed that DlGADs exhibit different expression patterns in various tissues and organs. The subcellular localization results showed that DlGAD5 was localized in the cytoplasm, suggesting that it played a role in the cytoplasm. Transient overexpression of DlGAD5 enhanced the expression levels of DlGADs and increased the activity of glutamate decarboxylase in longan embryogenic callus (EC), while the content of glutamic acid decreased. Thus, the DlGAD gene can play an important role in the early somatic embryogenesis of longan by responding to hormones such as IAA and ABA. DlGAD5 can affect the growth and development of longan by stimulating the expression of the DlGAD gene family, thereby increasing the GAD activity in the early SE of longan, participating in hormone synthesis and signaling pathways.

Transcriptome- and genome-wide systematic identification of expansin gene family and their expression in tuberous root development and stress responses in sweetpotato (Ipomoea batatas)

Introduction

Expansins (EXPs) are essential components of the plant cell wall that function as relaxation factors to directly promote turgor-driven expansion of the cell wall, thereby controlling plant growth and development and diverse environmental stress responses. EXPs genes have been identified and characterized in numerous plant species, but not in sweetpotato.

Results and methods

In the present study, a total of 59 EXP genes unevenly distributed across 14 of 15 chromosomes were identified in the sweetpotato genome, and segmental and tandem duplications were found to make a dominant contribution to the diversity of functions of the IbEXP family. Phylogenetic analysis showed that IbEXP members could be clustered into four subfamilies based on the EXPs from Arabidopsis and rice, and the regularity of protein motif, domain, and gene structures was consistent with this subfamily classification. Collinearity analysis between IbEXP genes and related homologous sequences in nine plants provided further phylogenetic insights into the EXP gene family. Cis-element analysis further revealed the potential roles of IbEXP genes in sweetpotato development and stress responses. RNA-seq and qRT-PCR analysis of eight selected IbEXPs genes provided evidence of their specificity in different tissues and showed that their transcripts were variously induced or suppressed under different hormone treatments (abscisic acid, salicylic acid, jasmonic acid, and 1-aminocyclopropane-1-carboxylic acid) and abiotic stresses (low and high temperature).

Discussion

These results provide a foundation for further comprehensive investigation of the functions of IbEXP genes and indicate that several members of this family have potential applications as regulators to control plant development and enhance stress resistance in plants.

Genome-Wide Identification of NDPK Family Genes and Expression Analysis under Abiotic Stress in Brassica napus

The NDPK gene family is an important group of genes in plants, playing a crucial role in regulating energy metabolism, growth, and differentiation, cell signal transduction, and response to abiotic stress. However, our understanding of the NDPK gene family in Brassica napus L. remains limited. This paper systematically analyzes the NDPK gene family in B. napus, particularly focusing on the evolutionary differences within the species. In this study, sixteen, nine, and eight NDPK genes were identified in B. napus and its diploid ancestors, respectively. These genes are not only homologous but also highly similar in their chromosomal locations. Phylogenetic analysis showed that the identified NDPK proteins were divided into four clades, each containing unique motif sequences, with most NDPKs experiencing a loss of introns/exons during evolution. Collinearity analysis revealed that the NDPK genes underwent whole-genome duplication (WGD) events, resulting in duplicate copies, and most of these duplicate genes were subjected to purifying selection. Cis-acting element analysis identified in the promoters of most NDPK genes elements related to a light response, methyl jasmonate response, and abscisic acid response, especially with an increased number of abscisic acid response elements in B. napus. RNA-Seq results indicated that NDPK genes in B. napus exhibited different expression patterns across various tissues. Further analysis through qRT-PCR revealed that BnNDPK genes responded significantly to stress conditions such as salt, drought, and methyl jasmonate. This study enhances our understanding of the NDPK gene family in B. napus, providing a preliminary theoretical basis for the functional study of NDPK genes and offering some references for further revealing the phenomenon of polyploidization in plants.

Genome-Wide Identification and Evolutionary and Expression Analyses of the Cyclin B Gene Family in Brassica napus

Cyclin B (CYCB) is a regulatory subunit of cyclin-dependent kinase (CDK), the concentration of which fluctuates to regulate cell cycle progression. Extensive studies have been performed on cyclins in numerous species, yet the evolutionary relationships and biological functions of the CYCB family genes in Brassica napus remain unclear. In this study, we identified 299 CYCB genes in 11 B. napus accessions. Phylogenetic analysis suggests that CYCB genes could be divided into three subfamilies in angiosperms and that the CYCB3 subfamily members may be a newer group that evolved in eudicots. The expansion of BnaCYCB genes underwent segmental duplication and purifying selection in genomes, and a number of drought-responsive and light-responsive cis-elements were found in their promoter regions. Additionally, expression analysis revealed that BnaCYCBs were strongly expressed in the developing seed and silique pericarp, as confirmed by the obviously reduced seed size of the mutant cycb3;1 in Arabidopsis thaliana compared with Col-0. This study provides a comprehensive evolutionary analysis of CYCB genes as well as insight into the biological function of CYCB genes in B. napus.