Comparative evolutionary analysis of chalcone synthase and alcohol dehydrogenase loci in Arabidopsis, Arabis, and related genera (Brassicaceae)

Genome-wide characterization and expression profiling of E2F/DP gene family members in response to abiotic stress in tomato (Solanum lycopersicum L.)

BACKGROUND

E2F/DP (Eukaryotic 2 transcription factor/dimerization partner) family proteins play an essential function in the cell cycle development of higher organisms. E2F/DP family genes have been reported only in a few plant species. However, comprehensive genome-wide characterization analysis of the E2F/DP gene family of Solanum lycopersicum has not been reported so far.

RESULTS

This study identified eight nonredundant SlE2F/DP genes that were classified into seven groups in the phylogenetic analysis. All eight genes had a single E2F-TDP domain and few genes had additional domains. Two segmental duplication gene pairs were observed within tomato, in addition to cis-regulatory elements, miRNA target sites and phosphorylation sites which play an important role in plant development and stress response in tomato. To explore the three-dimensional (3D) models and gene ontology (GO) annotations of SlE2F/DP proteins, we pointed to their putative transporter activity and their interaction with several putative ligands. The localization of SlE2F/DP-GFP fused proteins in the nucleus and endoplasmic reticulum suggested that they may act in other biological functions. Expression studies revealed the differential expression pattern of most of the SlE2F/DP genes in various organs. Moreover, the expression of E2F/DP genes against abiotic stress, particularly SlE2F/DP2 and/or SlE2F/DP7, was upregulated in response to heat, salt, cold and ABA treatment. Furthermore, the co-expression analysis of SlE2F/DP genes with multiple metabolic pathways was co-expressed with defence genes, transcription factors and so on, suggested their crucial role in various biological processes.

CONCLUSIONS

Overall, our findings provide a way to understand the structure and function of SlE2F/DP genes; it might be helpful to improve fruit development and tolerance against abiotic stress through marker-assisted selection or transgenic approaches.

A new chromosome-scale genome of wild Brassica oleracea provides insights into the domestication of Brassica crops

The cultivated diploid Brassica oleracea is an important vegetable crop, but the genetic basis of its domestication remains largely unclear in the absence of high-quality reference genomes of wild B. oleracea. Here, we report the first chromosome-level assembly of the wild Brassica oleracea L. W03 genome (total genome size, 630.7 Mb; scaffold N50, 64.6 Mb). Using the newly assembled W03 genome, we constructed a gene-based B. oleracea pangenome and identified 29 744 core genes, 23 306 dispensable genes, and 1896 private genes. We re-sequenced 53 accessions, representing six potential wild B. oleracea progenitor species. The results of the population genomic analysis showed that the wild B. oleracea populations had the highest level of diversity and represents the most closely related population to modern-day horticultural B. oleracea. In addition, the WUSCHEL gene was found to play a decisive role in domestication and to be involved in cauliflower and broccoli curd formation. We also illustrate the loss of disease-resistance genes during selection for domestication. Our results provide new insights into the domestication of B. oleracea and will facilitate the future genetic improvement of Brassica crops.

Upgraded durian genome reveals the role of chromosome reshuffling during ancestral karyotype evolution, lignin biosynthesis regulation, and stress tolerance

Durian (Durio zibethinus) is a tropical fruit that has a unique flavor and aroma. It occupies a significant phylogenetic position within the Malvaceae family. Extant core-eudicot plants are reported to share seven ancestral karyotypes that have undergone reshuffling, resulting in an abundant genomic diversity. However, the ancestral karyotypes of the Malvaceae family, as well as the evolution trajectory leading to the 28 chromosomes in durian, remain poorly understood. Here, we report the high-quality assembly of the durian genome with comprehensive comparative genomic analyses. By analyzing the collinear blocks between cacao and durian, we inferred 11 Malvaceae ancestral karyotypes. These blocks were present in a single-copy form in cacao and mainly in triplicates in durian, possibly resulting from a recent whole genome triplication (WGT) event that led to hexaploidization of the durian genome around 20 (17-24) million years ago. A large proportion of the duplicated genes in durian, such as those involved in the lignin biosynthesis module for phenylpropane biosynthesis, are derived directly from whole genome duplication, which makes it an important force in reshaping its genomic architecture. Transcriptome studies have revealed that genes involved in feruloyl-CoA formations were highly preferentially expressed in fruit peels, indicating that the thorns produced on durian fruit may comprise guaiacyl and syringyl lignins. Among all the analyzed transcription factors (TFs), members of the heat shock factor family (HSF) were the most significantly upregulated under heat stress. All subfamilies of genes encoding heat shock proteins (HSPs) in the durian genome appear to have undergone expansion. The potential interactions between HSF Dzi05.397 and HSPs were examined and experimentally verified. Our study provides a high-quality durian genome and reveals the reshuffling mechanism of ancestral Malvaceae chromosomes to produce the durian genome. We also provide insights into the mechanism underlying lignin biosynthesis and heat stress tolerance.

Identification and Functional Verification of the Glycosyltransferase Gene Family Involved in Flavonoid Synthesis in Rubus chingii Hu

Glycosylation is catalyzed by UDP-glycosyltransferase (UGT) and plays an important role in enriching the diversity of flavonoids. Rubus plants contain a lot of natural flavonoid glycosides, which are important plants with a homology of medicine and food. However, information about the Rubus UGT gene family is very limited. In this study, we carried out genome-wide analysis and identified the 172, 121, 130, 121 UGT genes in R. chingii, R. corchorifolius, R. idaeus, and R. occidentalis, respectively, and divided them into 18 groups. The analysis of the protein motif and gene structure showed that there were structural and functional conservations in the same group, but there were differences among different groups. Gene replication analysis showed that raspberry and dicotyledons had a higher homology. The expansion of the UGTs gene family was mainly driven by tandem replication events, and experienced purified selection during the long evolution of the raspberry. Cis-acting element analysis showed that they were related to plant growth and development, hormone regulation, and stress response. In addition, according to a comprehensive analysis of the co-expression network constructed by transcriptome data and phylogenetic homology, RchUGT169 was identified as a flavonoid glucosyltransferase. Through the transient expression in tobacco, it was verified that RchUGT169 could catalyze the conversion of kaempferol and quercetin to the corresponding flavonoid glycosides. In conclusion, this research enriched the understanding of the diversity of UGTs in Rubus and determined that RcUGT169 can catalyze flavonoids.

Enhancing rutin accumulation in Tartary buckwheat through a novel flavonoid transporter protein FtABCC2

Tartary buckwheat (Fagopyrum tataricum) is an annual coarse cereal from the Polygonaceae family, known for its high content of flavonoid compounds, particularly rutin. But so far, the mechanisms of the flavonoid transport and storage in Tartary buckwheat (TB) remain largely unexplored. This study focuses on ATP-binding cassette transporters subfamily C (ABCC) members, which are crucial for the biosynthesis and transport of flavonoids in plants. The evolutionary and expression pattern analyses of the ABCC genes in TB identified an ABCC protein gene, FtABCC2, that is highly correlated with rutin synthesis. Subcellular localization analysis revealed that FtABCC2 protein is specifically localized to the vacuole membrane. Heterologous expression of FtABCC2 in Saccharomyces cerevisiae confirmed that its transport ability of flavonoid glycosides such as rutin and isoquercetin, but not the aglycones such as quercetin and dihydroquercetin. Overexpression of FtABCC2 in TB hairy root lines resulted in a significant increase in total flavonoid and rutin content (P < 0.01). Analysis of the FtABCC2 promoter revealed potential cis-acting elements responsive to hormones, cold stress, mechanical injury and light stress. Overall, this study demonstrates that FtABCC2 can efficiently facilitate the transport of rutin into vacuoles, thereby enhancing flavonoids accumulation. These findings suggest that FtABCC2 is a promising candidate for molecular-assisted breeding aimed at developing high-flavonoid TB varieties.

The super-pangenome of Populus unveils genomic facets for its adaptation and diversification in widespread forest trees

Understanding the underlying mechanisms and links between genome evolution and adaptive innovations stands as a key goal in evolutionary studies. Poplars, among the world's most widely distributed and cultivated trees, exhibit extensive phenotypic diversity and environmental adaptability. In this study, we present a genus-level super-pangenome comprising 19 Populus genomes, revealing the likely pivotal role of private genes in facilitating local environmental and climate adaptation. Through the integration of pangenomes with transcriptomes, methylomes, and chromatin accessibility mapping, we unveil that the evolutionary trajectories of pangenes and duplicated genes are closely linked to local genomic landscapes of regulatory and epigenetic architectures, notably CG methylation in gene-body regions. Further comparative genomic analyses have enabled the identification of 142 202 structural variants across species that intersect with a significant number of genes and contribute substantially to both phenotypic and adaptive divergence. We have experimentally validated a ∼180-bp presence/absence variant affecting the expression of the CUC2 gene, crucial for leaf serration formation. Finally, we developed a user-friendly web-based tool encompassing the multi-omics resources associated with the Populus super-pangenome (http://www.populus-superpangenome.com). Together, the present pioneering super-pangenome resource in forest trees not only aids in the advancement of breeding efforts of this globally important tree genus but also offers valuable insights into potential avenues for comprehending tree biology.

Genome-wide identification, expression analysis, and abiotic stress response of the CBL and CIPK gene families in Artocarpus nanchuanensis

Artocarpus nanchuanensis, the northernmost species in the jackfruit genus, has great economic and horticultural value due to its nutritious fruit and beautiful tree shape. Calcineurin B-like proteins (CBLs) act as plant-specific Ca2+ sensors and participate in regulating plant responses to various abiotic stresses by interacting with CBL-interacting protein kinases (CIPKs). However, the characteristics and functions of the CBL and CIPK genes in A. nanchuanensis are still unclear. Here, we identified 14 CBL and 33 CIPK genes from the A. nanchuanensis genome, and based on phylogenetic analysis, they were divided into 4 and 7 clades, respectively. Gene structure and motif analysis indicated that the AnCBL and AnCIPK genes were relatively conserved. Colinear analysis showed that segmental duplication contributed to the expansion of the AnCBL and AnCIPK gene families. Expression analysis showed that AnCBL and AnCIPK genes were widely expressed in various tissues of A. nanchuanensis and exhibited tissue-specific expression. In addition, three genes (AnCBL6, AnCIPK7/8) may play important roles in response to salt, cold, and drought stresses. In summary, this study lays an important foundation for the improvement of stress resistance in A. nanchuanensis and provides new insight for the functional research on CBL and CIPK gene families.

Genome-wide analysis and expression profile of TCP gene family under drought and salinity stress condition in cowpea (Vigna unguiculata (L.) Walp.)

TCP transcription factors are known to regulate abiotic stress condition, but their role in V. unguiculata remains unexplored. So, in silico analysis and expression profile of the TCP gene family were performed in V. unguiculata to understand its role in response to heat and drought stress. A genome-wide search detected 28 TCPs (designated as VuTCPs) that were grouped into three subclasses by phylogenetic analysis. Gene structure, synteny, and phylogeny analyses of VuTCPs have shown a typical evolutionary path. One tandem and eight segmental duplication events were identified. Furthermore, identified duplicated, and orthologous VuTCP genes were under strong purifying selection pressure. A total of 15 SSRs were identified in the 12 VuTCPs, while 10 VuTCP genes were regulated by different miRNAs having a major role in abiotic stress tolerance. Analysed physicochemical properties, cis-acting elements, and gene ontology suggested that VuTCPs play various roles, including salinity and drought stress tolerance. qRT-PCR analysis showed that 11 and 15 VuTCPs were upregulated under drought and salinity stress conditions, respectively. Our findings provide comprehensive insights into the genomic characterization of the VuTCPs gene family in V. unguiculata, offering a foundation for understanding their structure, evolution, and role in abiotic stress tolerance.

Supplementary Information

The online version contains supplementary material available at 10.1007/s13205-024-03976-x.

Genome-Wide Identification of Glutathione S-Transferase Genes in Eggplant (Solanum melongena L.) Reveals Their Potential Role in Anthocyanin Accumulation on the Fruit Peel

Anthocyanins are ubiquitous pigments derived from the phenylpropanoid compound conferring red, purple and blue pigmentations to various organs of horticultural crops. The metabolism of flavonoids in the cytoplasm leads to the biosynthesis of anthocyanin, which is then conveyed to the vacuoles for storage by plant glutathione S-transferases (GST). Although GST is important for transporting anthocyanin in plants, its identification and characterization in eggplant (Solanum melongena L.) remains obscure. In this study, a total of 40 GST genes were obtained in the eggplant genome and classified into seven distinct chief groups based on the evolutionary relationship with Arabidopsis thaliana GST genes. The seven subgroups of eggplant GST genes (SmGST) comprise: dehydroascorbate reductase (DHAR), elongation factor 1Bγ (EF1Bγ), Zeta (Z), Theta(T), Phi(F), Tau(U) and tetra-chlorohydroquinone dehalogenase TCHQD. The 40 GST genes were unevenly distributed throughout the 10 eggplant chromosomes and were predominantly located in the cytoplasm. Structural gene analysis showed similarity in exons and introns within a GST subgroup. Six pairs of both tandem and segmental duplications have been identified, making them the primary factors contributing to the evolution of the SmGST. Light-related cis-regulatory elements were dominant, followed by stress-related and hormone-responsive elements. The syntenic analysis of orthologous genes indicated that eggplant, Arabidopsis and tomato (Solanum lycopersicum L.) counterpart genes seemed to be derived from a common ancestry. RNA-seq data analyses showed high expression of 13 SmGST genes with SmGSTF1 being glaringly upregulated on the peel of purple eggplant but showed no or low expression on eggplant varieties with green or white peel. Subsequently, SmGSTF1 had a strong positive correlation with anthocyanin content and with anthocyanin structural genes like SmUFGT (r = 0.9), SmANS (r = 0.85), SmF3H (r = 0.82) and SmCHI2 (r = 0.7). The suppression of SmGSTF1 through virus-induced gene silencing (VIGs) resulted in a decrease in anthocyanin on the infiltrated fruit surface. In a nutshell, results from this study established that SmGSTF1 has the potential of anthocyanin accumulation in eggplant peel and offers viable candidate genes for the improvement of purple eggplant. The comprehensive studies of the SmGST family genes provide the foundation for deciphering molecular investigations into the functional analysis of SmGST genes in eggplant.

Origin, evolution, and diversification of inositol 1,4,5-trisphosphate 3-kinases in plants and animals

BACKGROUND

In Eukaryotes, inositol polyphosphates (InsPs) represent a large family of secondary messengers and play crucial roes in various cellular processes. InsPs are synthesized through a series of pohophorylation reactions catalyzed by various InsP kinases in a sequential manner. Inositol 1,4,5-trisphosphate 3-kinase (IP3 3-kinase/IP3K), one member of InsP kinase, plays important regulation roles in InsPs metabolism by specifically phosphorylating inositol 1,4,5-trisphosphate (IP3) to inositol 1,3,4,5-tetrakisphosphate (IP4) in animal cells. IP3Ks were widespread in fungi, plants and animals. However, its evolutionary history and patterns have not been examined systematically.

RESULTS

A total of 104 and 31 IP3K orthologues were identified across 57 plant genomes and 13 animal genomes, respectively. Phylogenetic analyses indicate that IP3K originated in the common ancestor before the divergence of fungi, plants and animals. In most plants and animals, IP3K maintained low-copy numbers suggesting functional conservation during plant and animal evolution. In Brassicaceae and vertebrate, IP3K underwent one and two duplication events, respectively, resulting in multiple gene copies. Whole-genome duplication (WGD) was the main mechanism for IP3K duplications, and the IP3K duplicates have experienced functional divergence. Finally, a hypothetical evolutionary model for the IP3K proteins is proposed based on phylogenetic theory.

CONCLUSION

Our study reveals the evolutionary history of IP3K proteins and guides the future functions of animal, plant, and fungal IP3K proteins.

Genome-Wide Investigation of Class III Peroxidase Genes in Brassica napus Reveals Their Responsiveness to Abiotic Stresses

Brassica napus (B. napus) is susceptible to multiple abiotic stresses that can affect plant growth and development, ultimately reducing crop yields. In the past, many genes that provide tolerance to abiotic stresses have been identified and characterized. Peroxidase (POD) proteins, members of the oxidoreductase enzyme family, play a critical role in protecting plants against abiotic stresses. This study demonstrated a comprehensive investigation of the POD gene family in B. napus. As a result, a total of 109 POD genes were identified across the 19 chromosomes and classified into five distinct subgroups. Further, 44 duplicate events were identified; of these, two gene pairs were tandem and 42 were segmental. Synteny analysis revealed that segmental duplication was more prominent than tandem duplication among POD genes. Expression pattern analysis based on the RNA-seq data of B. napus indicated that BnPOD genes were expressed differently in various tissues; most of them were expressed in roots rather than in other tissues. To validate these findings, we performed RT-qPCR analysis on ten genes; these genes showed various expression levels under abiotic stresses. Our findings not only furnish valuable insights into the evolutionary dynamics of the BnPOD gene family but also serve as a foundation for subsequent investigations into the functional roles of POD genes in B. napus.

Molecular insights into self-incompatibility systems: From evolution to breeding

Plants have evolved diverse self-incompatibility (SI) systems for outcrossing. Since Darwin's time, considerable progress has been made toward elucidating this unrivaled reproductive innovation. Recent advances in interdisciplinary studies and applications of biotechnology have given rise to major breakthroughs in understanding the molecular pathways that lead to SI, particularly the strikingly different SI mechanisms that operate in Solanaceae, Papaveraceae, Brassicaceae, and Primulaceae. These best-understood SI systems, together with discoveries in other "nonmodel" SI taxa such as Poaceae, suggest a complex evolutionary trajectory of SI, with multiple independent origins and frequent and irreversible losses. Extensive exploration of self-/nonself-discrimination signaling cascades has revealed a comprehensive catalog of male and female identity genes and modifier factors that control SI. These findings also enable the characterization, validation, and manipulation of SI-related factors for crop improvement, helping to address the challenges associated with development of inbred lines. Here, we review current knowledge about the evolution of SI systems, summarize key achievements in the molecular basis of pollen‒pistil interactions, discuss potential prospects for breeding of SI crops, and raise several unresolved questions that require further investigation.

Comprehensive Analysis of NAC Transcription Factors Reveals Their Evolution in Malvales and Functional Characterization of AsNAC019 and AsNAC098 in Aquilaria sinensis

NAC is a class of plant-specific transcription factors that are widely involved in the growth, development and (a)biotic stress response of plants. However, their molecular evolution has not been extensively studied in Malvales, especially in Aquilaria sinensis, a commercial and horticultural crop that produces an aromatic resin named agarwood. In this study, 1502 members of the NAC gene family were identified from the genomes of nine species from Malvales and three model plants. The macroevolutionary analysis revealed that whole genome duplication (WGD) and dispersed duplication (DSD) have shaped the current architectural structure of NAC gene families in Malvales plants. Then, 111 NAC genes were systemically characterized in A. sinensis. The phylogenetic analysis suggests that NAC genes in A. sinensis can be classified into 16 known clusters and four new subfamilies, with each subfamily presenting similar gene structures and conserved motifs. RNA-seq analysis showed that AsNACs presents a broad transcriptional response to the agarwood inducer. The expression patterns of 15 AsNACs in A. sinensis after injury treatment indicated that AsNAC019 and AsNAC098 were positively correlated with the expression patterns of four polyketide synthase (PKS) genes. Additionally, AsNAC019 and AsNAC098 were also found to bind with the AsPKS07 promoter and activate its transcription. This comprehensive analysis provides valuable insights into the molecular evolution of the NAC gene family in Malvales plants and highlights the potential mechanisms of AsNACs for regulating secondary metabolite biosynthesis in A. sinensis, especially for the biosynthesis of 2-(2-phenyl) chromones in agarwood.

The Population Divergence and Genetic Basis of Local Adaptation of Wild Soybean (Glycine soja) in China

Glycine soja is the wild relative species of cultivated soybean. In this study, we investigated the population divergence and genetic basis of the local adaptation of wild soybean in China using genome-wide single-nucleotide polymorphisms (SNPs) of a population of 72 G. soja accessions. Using phylogenetic analysis, we observed that G. soja accessions clustered into three distinct groups, each corresponding to a specific geographic region, the northeastern region (NER), central region (CR), and southern region (SR), consistent with previous studies. Notably, we found a significant positive correlation between genetic and geographic distances. Further population structure analysis revealed each group was associated with an ancestral population and a specific geographic area. By utilizing the genome sequencing data of accessions from 16 different locations, we inferred the population history of these wild soybean groups. Our results indicate that the three groups diverged ~25,000 years ago, coinciding with the time of the last glacial maximum. The effective population size of the SR group expanded first, and subsequently, the NER and CR groups expanded approximately 5000 and 2500 years ago, respectively. Moreover, 83, 104, and 101 significant associated loci (SALs) were identified using genome-wide association analysis for annual mean temperature, annual precipitation, and latitude, respectively. Functional analysis of genes located in SALs highlighted candidate genes related to local adaptation. This study highlights the significant role of geographic isolation and environmental factors in shaping the genetic structure and adaptability of wild soybean populations. Furthermore, it emphasizes the value of wild soybean as a crucial genetic resource for enhancing the adaptability of cultivated soybeans, which have experienced a loss of genetic diversity due to domestication and intensive breeding practices. The insights gained from our research provide valuable information for the protection, conservation, and utilization of this important genetic resource.

Genome-wide characterization of regulator of chromosome condensation 1 (RCC1) gene family in Artemisia annua L. revealed a conservation evolutionary pattern

BACKGROUND

Artemisia annua is the major source for artemisinin production. The artemisinin content in A. annua is affected by different types of light especially the UV light. UVR8, a member of RCC1 gene family was found to be the UV-B receptor in plants. The gene structures, evolutionary history and expression profile of UVR8 or RCC1 genes remain undiscovered in A. annua.

RESULTS

Twenty-two RCC1 genes (AaRCC1) were identified in each haplotype genome of two diploid strains of A. annua, LQ-9 and HAN1. Varied gene structures and sequences among paralogs were observed. The divergence of most RCC1 genes occurred at 46.7 - 51 MYA which overlapped with species divergence of core Asteraceae during the Eocene, while no recent novel RCC1 members were found in A. annua genome. The number of RCC1 genes remained stable among eudicots and RCC1 genes underwent purifying selection. The expression profile of AaRCC1 is analogous to that of Arabidopsis thaliana (AtRCC1) when responding to environmental stress.

CONCLUSIONS

This study provided a comprehensive characterization of the AaRCC1 gene family and suggested that RCC1 genes were conserved in gene number, structures, constitution of amino acids and expression profiles among eudicots.

KIPEs3: Automatic annotation of biosynthesis pathways

Flavonoids and carotenoids are pigments involved in stress mitigation and numerous other processes. Both pigment classes can contribute to flower and fruit coloration. Flavonoid aglycones and carotenoids are produced by a pathway that is largely conserved across land plants. Glycosylations, acylations, and methylations of the flavonoid aglycones can be species-specific and lead to a plethora of biochemically diverse flavonoids. We previously developed KIPEs for the automatic annotation of biosynthesis pathways and presented an application on the flavonoid aglycone biosynthesis. KIPEs3 is an improved version with additional features and the potential to identify not just the core biosynthesis players, but also candidates involved in the decoration steps and in the transport of flavonoids. Functionality of KIPEs3 is demonstrated through the analysis of the flavonoid biosynthesis in Arabidopsis thaliana Nd-1, Capsella grandiflora, and Dioscorea dumetorum. We demonstrate the applicability of KIPEs to other pathways by adding the carotenoid biosynthesis to the repertoire. As a technical proof of concept, the carotenoid biosynthesis was analyzed in the same species and Daucus carota. KIPEs3 is available as an online service to enable access without prior bioinformatics experience. KIPEs3 facilitates the automatic annotation and analysis of biosynthesis pathways with a consistent and high quality in a large number of plant species. Numerous genome sequencing projects are generating a huge amount of data sets that can be analyzed to identify evolutionary patterns and promising candidate genes for biotechnological and breeding applications.

Genome resequencing reveals the evolutionary history of garlic reproduction traits

The propagation of cultivated garlic relies on vegetative cloves, thus flowers become non-essential for reproduction in this species, driving the evolution of reproductive feature-derived traits. To obtain insights into the evolutionary alteration of reproductive traits in the clonally propagated garlic, the evolutionary histories of two main reproduction-related traits, bolting and flower differentiation, were explored by genome analyses using 134 accessions displaying wide diversity in these two traits. Resequencing identified 272.8 million variations in the garlic genome, 198.0 million of which represent novel variants. Population analysis identified five garlic groups that have evolved into two clades. Gene expression, single-cell transcriptome sequencing, and genome-wide trait association analyses have identified numerous candidates that correlate with reproductive transition and flower development, some of which display distinct selection signatures. Selective forces acting on the B-box zinc finger protein-encoding Asa2G00291.1, the global transcription factor group E protein-encoding Asa5G01527.1, and VERNALIZATION INSENSITIVE 3-like Asa3G03399.1 appear to be representative of the evolution of garlic bolting. Plenty of novel genomic variations and trait-related candidates represent valuable resources for biological studies of garlic. Numerous selective signatures from genes associated with the two chosen reproductive traits provide important insights into the evolutionary history of reproduction in this clonally propagated crop.

Global Brassicaceae phylogeny based on filtering of 1,000-gene dataset

The mustard family (Brassicaceae) is a scientifically and economically important family, containing the model plant Arabidopsis thaliana and numerous crop species that feed billions worldwide. Despite its relevance, most phylogenetic trees of the family are incompletely sampled and often contain poorly supported branches. Here, we present the most complete Brassicaceae genus-level family phylogenies to date (Brassicaceae Tree of Life or BrassiToL) based on nuclear (1,081 genes, 319 of the 349 genera; 57 of the 58 tribes) and plastome (60 genes, 265 genera; all tribes) data. We found cytonuclear discordance between the two, which is likely a result of rampant hybridization among closely and more distantly related lineages. To evaluate the impact of such hybridization on the nuclear phylogeny reconstruction, we performed five different gene sampling routines, which increasingly removed putatively paralog genes. Our cleaned subset of 297 genes revealed high support for the tribes, whereas support for the main lineages (supertribes) was moderate. Calibration based on the 20 most clock-like nuclear genes suggests a late Eocene to late Oligocene origin of the family. Finally, our results strongly support a recently published new family classification, dividing the family into two subfamilies (one with five supertribes), together representing 58 tribes. This includes five recently described or re-established tribes, including Arabidopsideae, a monogeneric tribe accommodating Arabidopsis without any close relatives. With a worldwide community of thousands of researchers working on Brassicaceae and its diverse members, our new genus-level family phylogeny will be an indispensable tool for studies on biodiversity and plant biology.

Population genetics of Camellia granthamiana, an endangered plant species with extremely small populations in China

Introduction: Camellia, the largest genus of Theaceae, is well-known for having high economic values. Camellia granthamiana demonstrates large beautiful flowers with some primitive characters, such as multiple large and persistent bracteoles and sepals, was listed as Vulnerable species on the IUCN Red List. Methods: In this study, we investigated all possible records of the species, and sampled four natural populations and five cultivated individuals. By applying shallow-genome sequencing for nine individuals and RAD-seq sequencing for all the sampled 77 individuals, we investigated population genetic diversity and population structure of the species. Results and discussion: The results showed that the population sampled from Fengkai, previously identified as C. albogigias, possessed different plastid genome from other species possibly due to plastid capture; the species possesses strong population structure possibly due to the effect of isolation by distance, habitat fragmentation, and self-crossing tendency of the species, whose effective population size declined quickly in the past 4,000 years. Nevertheless, C. granthamiana maintains a medium level of genetic diversity within population, and significant differentiation was observed among the four investigated populations, it is anticipated that more populations are expected to be found and all these extant populations should be taken into instant protection.

The Late Embryogenesis Abundant Proteins in Soybean: Identification, Expression Analysis, and the Roles of GmLEA4_19 in Drought Stress

Late embryogenesis abundant (LEA) proteins play important roles in regulating plant growth and responses to various abiotic stresses. In this research, a genome-wide survey was conducted to recognize the LEA genes in Glycine max. A total of 74 GmLEA was identified and classified into nine subfamilies based on their conserved domains and the phylogenetic analysis. Subcellular localization, the duplication of genes, gene structure, the conserved motif, and the prediction of cis-regulatory elements and tissue expression pattern were then conducted to characterize GmLEAs. The expression profile analysis indicated that the expression of several GmLEAs was a response to drought and salt stress. The co-expression-based gene network analysis suggested that soybean LEA proteins may exert regulatory effects through the metabolic pathways. We further explored GnLEA4_19 function in Arabidopsis and the results suggests that overexpressed GmLEA4_19 in Arabidopsis increased plant height under mild or serious drought stress. Moreover, the overexpressed GmLEA4_19 soybean also showed a drought tolerance phenotype. These results indicated that GmLEA4_19 plays an important role in the tolerance to drought and will contribute to the development of the soybean transgenic with enhanced drought tolerance and better yield. Taken together, this study provided insight for better understanding the biological roles of LEA genes in soybean.

A Hyb-Seq phylogeny of Boechera and related genera using a combination of Angiosperms353 and Brassicaceae-specific bait sets

PREMISE

Although Boechera (Boechereae, Brassicaceae) has become a plant model system for both ecological genomics and evolutionary biology, all previous phylogenetic studies have had limited success in resolving species relationships within the genus. The recent effective application of sequence data from target enrichment approaches to resolve the evolutionary relationships of several other challenging plant groups prompted us to investigate their usefulness in Boechera and Boechereae.

METHODS

To resolve the phylogeny of Boechera and closely related genera, we utilized the Hybpiper pipeline to analyze two combined bait sets: Angiosperms353, with broad applicability across flowering plants; and a Brassicaceae-specific bait set designed for use in the mustard family. Relationships for 101 samples representing 81 currently recognized species were inferred from a total of 1114 low-copy nuclear genes using both supermatrix and species coalescence methods.

RESULTS

Our analyses resulted in a well-resolved and highly supported phylogeny of the tribe Boechereae. Boechereae is divided into two major clades, one comprising all western North American species of Boechera, the other encompassing the eight other genera of the tribe. Our understanding of relationships within Boechera is enhanced by the recognition of three core clades that are further subdivided into robust regional species complexes.

CONCLUSIONS

This study presents the first broadly sampled, well-resolved phylogeny for most known sexual diploid Boechera. This effort provides the foundation for a new phylogenetically informed taxonomy of Boechera that is crucial for its continued use as a model system.

Strong postmating reproductive isolation in Mimulus section Eunanus

Postmating reproductive isolation can help maintain species boundaries when premating barriers to reproduction are incomplete. The strength and identity of postmating reproductive barriers are highly variable among diverging species, leading to questions about their genetic basis and evolutionary drivers. These questions have been tackled in model systems but are less often addressed with broader phylogenetic resolution. In this study we analyse patterns of genetic divergence alongside direct measures of postmating reproductive barriers in an overlooked group of sympatric species within the model monkeyflower genus, Mimulus. Within this Mimulus brevipes species group, we find substantial divergence among species, including a cryptic genetic lineage. However, rampant gene discordance and ancient signals of introgression suggest a complex history of divergence. In addition, we find multiple strong postmating barriers, including postmating prezygotic isolation, hybrid seed inviability and hybrid male sterility. M. brevipes and M. fremontii have substantial but incomplete postmating isolation. For all other tested species pairs, we find essentially complete postmating isolation. Hybrid seed inviability appears linked to differences in seed size, providing a window into possible developmental mechanisms underlying this reproductive barrier. While geographic proximity and incomplete mating isolation may have allowed gene flow within this group in the distant past, strong postmating reproductive barriers today have likely played a key role in preventing ongoing introgression. By producing foundational information about reproductive isolation and genomic divergence in this understudied group, we add new diversity and phylogenetic resolution to our understanding of the mechanisms of plant speciation.

Whole-genome identification and multiple abiotic stresses expression pattern profiling analysis of PLATZ transcription factor family members in Pecan (Carya illinoensis)

Plant AT-rich sequence and zinc-binding (PLATZ), as a plant-specific transcription factor, have been identified and studied in a variety of plants. However, there are no reports about PLATZ proteins in Carya illinoensis (pecan). Here, 24 C. illinoensis CiPLATZs have been identified and divided into 4 groups. Gene structure, motif composition, conserved domain and cis-acting elements analysis indicated that the PLATZ gene family was highly conserved. Transcriptome data combined with qRT-PCR analysis revealed that CiPLATZ6, CiPLATZ12, CiPLATZ13, CiPLATZ14 and CiPLATZ23 were highly expressed in multiple tissues of C. illinoensis and strongly responded to drought, salt and heat stress. Among them, CiPLATZ6, CiPLATZ12 and CiPLATZ23 were all located in the nucleus and had no transcriptional autoactivation ability in yeast cells, and acted as transcriptional suppressors in plants. In addition, the CiPLATZ23-overexpressing transgenic Arabidopsis thaliana showed enhanced tolerance to drought. Measurements of physiological indicators and analysis of stress-related genes expression levels in transgenic A. thaliana were used to support this conclusion. The results of this study are helpful to understand the structural feature and function of CiPLATZs, and provide candidate genes for molecular breeding of drought tolerance of C. illinoensis.

Genome-wide identification and characterization of ADH gene family and the expression under different abiotic stresses in tomato (Solanum lycopersicum L.)

The SlADH gene plays a key role in environmental stress response. However, limited studies exist regarding the tomato SlADH gene. In this study, we identified 35 SlADH genes in tomato by genome-wide identification. Among the 12 chromosomes of tomato, SlADH gene is distributed on 10 chromosomes, among which the 7th and 10th chromosomes have no family members, while the 11th chromosome has the most members with 8 family members. Members of this gene family are characterized by long coding sequences, few amino acids, and introns that make up a large proportion of the genetic structure of most members of this family. Moreover, the molecular weight of the proteins of the family members was similar, and the basic proteins were mostly, and the overall distribution was relatively close to neutral (pI = 7). This may indicate that proteins in this family have a more conserved function. In addition, a total of four classes of cis-acting elements were detected in all 35 SlADH promoter regions, most of which were associated with biotic and abiotic stresses. The results indicate that SlADH gene had a certain response to cold stress, salt stress, ABA treatment and PEG stress. This study provides a new candidate gene for improving tomato stress resistance.

Systematic evolution of bZIP transcription factors in Malvales and functional exploration of AsbZIP14 and AsbZIP41 in Aquilaria sinensis

Introduction

Agarwood, the dark-brown resin produced by Aquilaria trees, has been widely used as incense, spice, perfume or traditional medicine and 2-(2-phenethyl) chromones (PECs) are the key markers responsible for agarwood formation. But the biosynthesis and regulatory mechanism of PECs were still not illuminated. The transcription factor of basic leucine zipper (bZIP) presented the pivotal regulatory roles in various secondary metabolites biosynthesis in plants, which might also contribute to regulate PECs biosynthesis. However, molecular evolution and function of bZIP are rarely reported in Malvales plants, especially in Aquilaria trees.

Methods and results

Here, 1,150 bZIPs were comprehensively identified from twelve Malvales and model species genomes and the evolutionary process were subsequently analyzed. Duplication types and collinearity indicated that bZIP is an ancient or conserved TF family and recent whole genome duplication drove its evolution. Interesting is that fewer bZIPs in A. sinensis than that species also experienced two genome duplication events in Malvales. 62 AsbZIPs were divided into 13 subfamilies and gene structures, conservative domains, motifs, cis-elements, and nearby genes of AsbZIPs were further characterized. Seven AsbZIPs in subfamily D were significantly regulated by ethylene and agarwood inducer. As the typical representation of subfamily D, AsbZIP14 and AsbZIP41 were localized in nuclear and potentially regulated PECs biosynthesis by activating or suppressing type III polyketide synthases (PKSs) genes expression via interaction with the AsPKS promoters.

Discussion

Our results provide a basis for molecular evolution of bZIP gene family in Malvales and facilitate the understanding the potential functions of AsbZIP in regulating 2-(2-phenethyl) chromone biosynthesis and agarwood formation.

Heat shock factor binding protein BrHSBP1 regulates seed and pod development in Brassica rapa

Plant heat shock factor binding proteins (HSBPs) are well known for their implication in the negative regulation of heat stress response (HSR) pathways. Herein, we report on the hitherto unknown functions of HSBP1 in Brassica rapa (BrHSBP1). BrHBSP1 was found to be predominant in flower buds and young leaves, while its segmental duplicate, BrHSBP1-like, was abundant in green siliques. Exposure to abiotic stress conditions, such as heat, drought, cold, and H2O2, and to phytohormones was found to differentially regulate BrHSBP1. The activity of BrHSBP1-GFP fusion proteins revealed their cellular localization in nuclei and cytosols. Transgenic overexpression of BrHSBP1 (BrHSBP1OX) improved pod and seed sizes, while CRISPR-Cas BrHSBP1 knock-out mutants (Brhsbp1_KO) were associated with aborted seed and pod development. The transcriptomic signatures of BrHSBP1OX and Brhsbp1_KO lines revealed that 360 and 2381 genes, respectively, were differentially expressed (Log2FC≥2, padj<0.05) expressed relative to control lines. In particular, developmental processes, including plant reproductive structure development (RSD)-related genes, were relatively downregulated in Brhsbp1_KO. Furthermore, yeast two-hybrid assays confirmed that BrHSBP1 can physically bind to RSD and other genes. Taking the findings together, it is clear that BrHSBP1 is involved in seed development via the modulation of RSD genes. Our findings represent the addition of a new regulatory player in seed and pod development in B. rapa.

Genome-wide identification, characterization, evolution and expression analysis of the DIR gene family in potato (Solanum tuberosum)

The dirigent (DIR) gene is a key player in environmental stress response and has been identified in many multidimensional tube plant species. However, there are few studies on the StDIR gene in potato. In this study, we used genome-wide identification to identify 31 StDIR genes in potato. Among the 12 potato chromosomes, the StDIR gene was distributed on 11 chromosomes, among which the third chromosome did not have a family member, while the tenth chromosome had the most members with 11 members. 22 of the 31 StDIRs had a classical DIR gene structure, with one exon and no intron. The conserved DIR domain accounts for most of the proteins in the 27 StDIRs. The structure of the StDIR gene was analyzed and ten different motifs were detected. The StDIR gene was divided into three groups according to its phylogenetic relationship, and 22 duplicate genes were identified. In addition, four kinds of cis-acting elements were detected in all 31 StDIR promoter regions, most of which were associated with biotic and abiotic stress. The findings demonstrated that the StDIR gene exhibited specific responses to cold stress, salt stress, ABA, and drought stress. This study provides new candidate genes for improving potato's resistance to stress.

Identification and Expression Analysis of the Alfin-like Gene Family in Tomato and the Role of SlAL3 in Salt and Drought Stresses

Alfin-like (AL) transcription factors are a family of plant-specific genes with a PHD-finger-like structural domain at the C-terminus and a DUF3594 structural domain at the N-terminus that play important roles in plant development and stress response. In the present study, genome-wide identification and analysis were performed of the AL protein family in cultivated tomato (Solanum lycopersicum) and three wild relatives (S. pennellii, S. pimpinellifolium, and S. lycopersicoides) to evaluate their response to different abiotic stresses. A total of 39 ALs were identified and classified into four groups and based on phylogenetic tree and evolutionary analysis were shown to have formed prior to the differentiation of monocotyledons and dicots. Moreover, cis-acting element analysis revealed that various phytohormone response and abiotic stress response elements were highly existed in tomato. In addition, further analysis of the SlAL3 gene revealed that its expression was induced by drought and salt stresses and localized to the nucleus. In conclusion, our findings concerning AL genes provide useful information for further studies on their functions and regulatory mechanisms and provide theoretical references for studying AL gene response to abiotic stresses in plants.

High-Quality Genome Assembly and Genome-Wide Association Study of Male Sterility Provide Resources for Flax Improvement

Flax is an economic crop with a long history. It is grown worldwide and is mainly used for edible oil, industry, and textiles. Here, we reported a high-quality genome assembly for "Neiya No. 9", a popular variety widely grown in China. Combining PacBio long reads, Hi-C sequencing, and a genetic map reported previously, a genome assembly of 473.55 Mb was constructed, which covers ~94.7% of the flax genome. These sequences were anchored onto 15 chromosomes. The N50 lengths of the contig and scaffold were 0.91 Mb and 31.72 Mb, respectively. A total of 32,786 protein-coding genes were annotated, and 95.9% of complete BUSCOs were found. Through morphological and cytological observation, the male sterility of flax was considered dominant nuclear sterility. Through GWAS analysis, the gene LUSG00017705 (cysteine synthase gene) was found to be closest to the most significant SNP, and the expression level of this gene was significantly lower in male sterile plants than in fertile plants. Among the significant SNPs identified in the GWAS analysis, only two were located in the coding region, and these two SNPs caused changes in the protein encoded by LUSG00017565 (cysteine protease gene). It was speculated that these two genes may be related to male sterility in flax. This is the first time the molecular mechanism of male sterility in flax has been reported. The high-quality genome assembly and the male sterility genes revealed, provided a solid foundation for flax breeding.

Genome-wide identification and characterization of glutathione S-transferase gene family in quinoa (Chenopodium quinoa Willd.)

The present investigation was envisaged for large scale in-silico genome wide identification and characterization of glutathione S-transferases (GSTs) in Chenopodium quinoa. In this study, a total of 120 GST genes (CqGSTs) were identified and divided into 11 classes of which tau and phi were highest in numbers. The average protein length of protein was found to be 279.06 with their corresponding average molecular weight of 31,819.4 kDa. The subcellular localization analysis results showed that proteins were centrally localized in the cytoplasm followed by chloroplast, mitochondria and plastids. Structural analysis revealed the presence of 2 -14 exons in CqGST genes. Most of the proteins possessed two exon one intron organization. MEME analysis identified 15 significantly conserved motifs with a width of 6-50 amino acids. Motifs 1, 3, 2, 5, 6, 8, 9 and 13 were found specifically in tau class family; motifs 3, 4, 5, 6, 7 and 9 were found in phi class gene family, while motifs 3, 4, 13 and 14 were found in metaxin class. Multiple sequence alignment revealed highly conserved N-terminus with active site serine (Ser; S) or cysteine (Cys; C) residue for the activation of GSH binding and GST catalytic activity. The gene loci were found to be unevenly distributed across 18 different chromosomes with a maximum of 17 genes located on chromosome number 7. Dominance of alpha helix was followed by coil, extended strand and beta turns. Gene duplication analysis revealed that segmental duplication and purifying type selection were highest in number and found to be main source of expansion of GST gene family. Cis acting regulatory elements analysis showed the presence of 21 different elements involved in stress, hormone and light response and cellular development. The evolutionary relationship of CqGST proteins carried out using maximum likelihood method revealed that all the tau and phi class GSTs were closely associated with those of G. max, O. sativa and A. thaliana. Molecular docking of GST molecules with the fungicide metalaxyl showed that the CqGSTF1 had the lowest binding energy. The comprehensive study of CqGST gene family in quinoa provides groundwork for further functional analysis of CqGST genes in the species at molecular level and has potential applications in plant breeding.

A High-Continuity Genome Assembly of Chinese Flowering Cabbage (Brassica rapa var. parachinensis) Provides New Insights into Brassica Genome Structure Evolution

Chinese flowering cabbage (Brassica rapa var. parachinensis) is a popular and widely cultivated leaf vegetable crop in Asia. Here, we performed a high quality de novo assembly of the 384 Mb genome of 10 chromosomes of a typical cultivar of Chinese flowering cabbage with an integrated approach using PacBio, Illumina, and Hi-C technology. We modeled 47,598 protein-coding genes in this analysis and annotated 52% (205.9/384) of its genome as repetitive sequences including 17% in DNA transposons and 22% in long terminal retrotransposons (LTRs). Phylogenetic analysis reveals the genome of the Chinese flowering cabbage has a closer evolutionary relationship with the AA diploid progenitor of the allotetraploid species, Brassica juncea. Comparative genomic analysis of Brassica species with different subgenome types (A, B and C) reveals that the pericentromeric regions on chromosome 5 and 6 of the AA genome have been significantly expanded compared to the orthologous genomic regions in the BB and CC genomes, largely driven by LTR-retrotransposon amplification. Furthermore, we identified a large number of structural variations (SVs) within the B. rapa lines that could impact coding genes, suggesting the functional significance of SVs on Brassica genome evolution. Overall, our high-quality genome assembly of the Chinese flowering cabbage provides a valuable genetic resource for deciphering the genome evolution of Brassica species and it can potentially serve as the reference genome guiding the molecular breeding practice of B. rapa crops.

Molecular evolution, diversification, and expression assessment of MADS gene family in Setaria italica, Setaria viridis, and Panicum virgatum

KEY MESSAGE

This paper sheds light on the evolution and expression patterns of MADS genes in Setaria and Panicum virgatum. SiMADS51 and SiMADS64 maybe involved in the ABA-dependent pathway of drought response. The MADS gene family is a key regulatory factor family that controls growth, reproduction, and response to abiotic stress in plants. However, the molecular evolution of this family is rarely reported. Here, a total of 265 MADS genes were identified in Setaria italica (foxtail millet), Setaria viridis (green millet), and Panicum virgatum (switchgrass) and analyzed by bioinformatics, including physicochemical characteristics, subcellular localization, chromosomal position and duplicate, motif distribution, genetic structure, genetic evolvement, and expression patterns. Phylogenetic analysis was used to categorize these genes into M and MIKC types. The distribution of motifs and gene structure were similar for the corresponding types. According to a collinearity study, the MADS genes have been mostly conserved during evolution. The principal cause of their expansion is segmental duplication. However, the MADS gene family tends to shrink in foxtail millet, green millet, and switchgrass. The MADS genes were subjected to purifying selection, but several positive selection sites were also identified in three species. And most of the promoters of MADS genes contain cis-elements related to stress and hormonal response. RNA-seq and quantitative Real-time PCR (qRT-PCR) analysis also were examined. SiMADS genes expression levels are considerably changed in reaction to various treatments, following qRT-PCR analysis. This sheds fresh light on the evolution and expansion of the MADS family in foxtail millet, green millet, and switchgrass, and lays the foundation for further research on their functions.

Evolution of the 14-3-3 gene family in monocotyledons and dicotyledons and validation of MdGRF13 function in transgenic Arabidopsis thaliana

KEY MESSAGE

The 14-3-3 family is more highly conserved among monocotyledons, and overexpression of MdGRF13 improved drought and salt tolerance in transgenic Arabidopsis thaliana. The 14-3-3 are highly conserved regulatory proteins found in eukaryotes and play an essential role in plant growth, development and stress response. However, the 14-3-3 gene family evolution in monocotyledons and dicotyledons and the biological functions of the MdGRF13 under abiotic stress remain unknown. In our study, 195 members of the 14-3-3 family were identified from 12 species and divided into ε group and the Non-ε group. Synteny analysis within the 14-3-3 family indicated that segmental duplication events contributed to the expansion of the family. Selective pressure analysis indicated that purifying selection was a vital force in the 14-3-3 genes evolution, and monocotyledons had a lower million years ago (Mya) mean values than dicotyledons. Meanwhile, the codon adaptation index (CAI) and frequency of optical codons (FOP) are higher and the effective number of codons (Nc) is lower in monocotyledons 14-3-3 genes compared to dicotyledons. Moreover, the yeast two-hybrid (Y2H) demonstrated that MdGRF13 interacts with MdRD22, MdLHP1a and MdMORF1. Significantly, the malondialdehyde (MDA) content and relative conductivity were decreased, while the superoxide dismutase (SOD), peroxidase (POD), and catalase (CAT) activities were increased in transgenic Arabidopsis compared to the wild type (WT) under drought and salt stress. These results suggest that overexpression of MdGRF13 significantly improved the tolerance to drought and salt stress in transgenic Arabidopsis. Thus, our results provide a theoretical basis for exploring the evolution and function of the 14-3-3 gene family in monocotyledons and dicotyledons.

Comprehensive Analysis of the INDETERMINATE DOMAIN (IDD) Gene Family and Their Response to Abiotic Stress in Zea mays

Transcription factors (TFs) are important regulators of numerous gene expressions due to their ability to recognize and combine cis-elements in the promoters of target genes. The INDETERMINATE DOMAIN (IDD) gene family belongs to a subfamily of C2H2 zinc finger proteins and has been identified only in terrestrial plants. Nevertheless, little study has been reported concerning the genome-wide analysis of the IDD gene family in maize. In total, 22 ZmIDD genes were identified, which can be distributed on 8 chromosomes in maize. On the basis of evolutionary relationships and conserved motif analysis, ZmIDDs were categorized into three clades (1, 2, and 3), each owning 4, 6, and 12 genes, respectively. We analyzed the characteristics of gene structure and found that 3 of the 22 ZmIDD genes do not contain an intron. Cis-element analysis of the ZmIDD promoter showed that most ZmIDD genes possessed at least one ABRE or MBS cis-element, and some ZmIDD genes owned the AuxRR-core, TCA-element, TC-rich repeats, and LTR cis-element. The Ka:Ks ratio of eight segmentally duplicated gene pairs demonstrated that the ZmIDD gene families had undergone a purifying selection. Then, the transcription levels of ZmIDDs were analyzed, and they showed great differences in diverse tissues as well as abiotic stresses. Furthermore, regulatory networks were constructed through the prediction of ZmIDD-targeted genes and miRNAs, which can inhibit the transcription of ZmIDDs. In total, 6 ZmIDDs and 22 miRNAs were discovered, which can target 180 genes and depress the expression of 9 ZmIDDs, respectively. Taken together, the results give us valuable information for studying the function of ZmIDDs involved in plant development and climate resilience in maize.

Genome-wide characterization of ubiquitin-conjugating enzyme gene family explores its genetic effects on the oil content and yield of Brassica napus

Ubiquitin-conjugating enzyme (UBC) is a critical part of the ubiquitin-proteasome pathway and plays crucial roles in growth, development and abiotic stress response in plants. Although UBC genes have been detected in several plant species, characterization of this gene family at the whole-genome level has not been conducted in Brassica napus. In the present study, 200 putative BnUBCs were identified in B. napus, which were clustered into 18 subgroups based on phylogenetic analysis. BnUBCs within each subgroup showed relatively conserved gene architectures and motifs. Moreover, the gene expression patterns in various tissues as well as the identification of cis-acting regulatory elements in BnUBC promoters suggested further investigation of their potential functions in plant growth and development. Furthermore, three BnUBCs were predicted as candidate genes for regulating agronomic traits related to oil content and yield through association mapping. In conclusion, this study provided a wealth of information on the UBC family in B. napus and revealed their effects on oil content and yield, which will aid future functional research and genetic breeding of B. napus.

PlantLTRdb: An interactive database for 195 plant species LTR-retrotransposons

LTR-retrotransposons (LTR-RTs) are a large group of transposable elements that replicate through an RNA intermediate and alter genome structure. The activities of LTR-RTs in plant genomes provide helpful information about genome evolution and gene function. LTR-RTs near or within genes can directly alter gene function. This work introduces PlantLTRdb, an intact LTR-RT database for 195 plant species. Using homology- and de novo structure-based methods, a total of 150.18 Gbp representing 3,079,469 pseudomolecules/scaffolds were analyzed to identify, characterize, annotate LTR-RTs, estimate insertion ages, detect LTR-RT-gene chimeras, and determine nearby genes. Accordingly, 520,194 intact LTR-RTs were discovered, including 29,462 autonomous and 490,732 nonautonomous LTR-RTs. The autonomous LTR-RTs included 10,286 Gypsy and 19,176 Copia, while the nonautonomous were divided into 224,906 Gypsy, 218,414 Copia, 1,768 BARE-2, 3,147 TR-GAG and 4,2497 unknown. Analysis of the identified LTR-RTs located within genes showed that a total of 36,236 LTR-RTs were LTR-RT-gene chimeras and 11,619 LTR-RTs were within pseudo-genes. In addition, 50,026 genes are within 1 kbp of LTR-RTs, and 250,587 had a distance of 1 to 10 kbp from LTR-RTs. PlantLTRdb allows researchers to search, visualize, BLAST and analyze plant LTR-RTs. PlantLTRdb can contribute to the understanding of structural variations, genome organization, functional genomics, and the development of LTR-RT target markers for molecular plant breeding. PlantLTRdb is available at https://bioinformatics.um6p.ma/PlantLTRdb.

Enrichment of intrinsically disordered residues in ohnologs facilitates abiotic stress resilience in Brassica rapa

Arabidopsis thaliana and Brassica rapa are in the same evolutionary lineage, although the latter experienced an additional whole genome triplication event. Therefore, it would be intriguing to investigate the traits that gene duplication imposes to mediate plant stress tolerance. Here, we noticed that B. rapa abiotic stress resistance (ASR) genes which code at least one stress responsive domain have a significantly higher number of paralogs than A. thaliana. Analysing the disordered content of the ASR genes in both species, we found that intrinsically disordered residues (IDR) are specifically enriched in whole genome duplication (WGD) derived paralogs. Subsequently, domain similarity analysis between WGD pairs of both species has revealed that majority of WGD pairs in B. rapa did not share domains with each other. Furthermore, domain enrichment analysis has shown that B. rapa paralogs contain 36 distinct stress responsive enriched domains, significantly higher than A. thaliana paralogs. Next, we performed MSA to investigate the domain conservation between orthologs and ohnologs pairs, we found that 80.13% of B. rapa ohnologs acquire new domains, depicting the fact that ohnologs play a significant role in stress-related behaviours. The average IDR content of the ohnologs enriching new domains after gene duplication in B. rapa (0.19), is also significantly higher than A. thaliana (0.04). Interestingly, we also found that all of these attributes i.e., exhibiting higher number of WGD paralogs and enhancement of IDR in ASR genes of B. rapa compared to A. thaliana is exclusive for ASR genes only. No such significant differences were observed in randomly selected non-ASR genes between the two species. Together these results provide strong support for the hypothesis that augmentation of IDR content followed by a whole genome duplication event imposes the stress resistance potentiality in B. rapa. This research will shed light on the mechanism of how B. rapa is able to successfully adapt to stress over the evolutionary timescale.

Gnawing pressure led to the expansion of JAZ genes in angiosperms

A long-standing problem in evolutionary biology is why some populations differentiate into many species while the majority do not. Angiosperms is an excellent group for investigating this problem because their diversity is unevenly distributed in space and phylogeny. Plant hormone participates in growth, development and defense. However, jasmonic acid (JA) was the only hormone response to bites. We first searched jasmonate ZIM-domain (JAZ), AUXIN/INDOLE ACETIC ACID (IAA / aux), PYR/PYL/RCAR (PYL), DELLA, and SUPPRESSOR OF MAX2 1-like (SMAX) in 272 plant species. We found the gene number change trends were consistent with origination rates and species numbers of angiosperms. So, 26 representative species were selected as an example for further analysis. The results showed JAZ had experienced two lineage-specific gene expansion events in angiosperms, which coincided with increases in mammalian body size and dental diversity. The proliferation of large herbivores as a results of mammalian prosperity after dinosaur extinction may be related to angiosperm evolution and bursting. The proliferation of large herbivores as the result of mammalian prosperity after the extinction of the dinosaurs was related to angiosperm evolution and bursting. Overall, our study uncovered a previously unknown co-evolution mechanism in terrestrial plants exposed to extreme environmental conditions.

Identification of HubHLH family and key role of HubHLH159 in betalain biosynthesis by activating the transcription of HuADH1, HuCYP76AD1-1, and HuDODA1 in pitaya

Basic helix-loop-helix (bHLH) proteins are dimeric transcription factors (TFs) involved in various plant physiological and biological processes. Despite this, little is known about the molecular properties and roles of bHLH TFs in pitaya betalain biosynthesis. Here we report the identification of 165 HubHLH genes in H. undantus genome, their chromosomal distribution, physiochemical characteristics, conserved motifs, gene structure, phylogeny and synteny of HubHLH genes. Based on phylogenetic relationship analysis, the 165 HubHLHs were divided into 26 subfamilies and unequally distributed on the 11 chromosomes of pitaya. Based on the pitaya transcriptome data, a candidate gene HubHLH159 was obtained, and the real-time quantitative PCR analysis confirmed that HubHLH159 showed a high expression level in 'Guanhuahong' pitaya (red-pulp) at mature stage, indicating its role in betalain biosynthesis. HubHLH159 is a Group II protein and contains a bHLH domain. It is a nuclear protein with transcriptional activation activity. Dual luciferase reporter assays and virus-induced gene silencing (VIGS) experiments showed that HubHLH159 promotes betalain biosynthesis by activating the expression of HuADH1, HuCYP76AD1-1, and HuDODA1. The results of the present study lay a new theoretical reference for the regulation of pitaya betalain biosynthesis and also provides as essential basis for the future analysis of the functions of HubHLH gene family.

Unravelling structural, functional, evolutionary and genetic basis of SWEET transporters regulating abiotic stress tolerance in maize

Sugars Will Eventually be Exported Transporters (SWEETs) are the novel sugar transporters widely distributed among living systems. SWEETs play a crucial role in various bio-physiological processes, viz., plant developmental, nectar secretion, pollen development, and regulation of biotic and abiotic stresses, in addition to their prime sugar-transporting activity. Thus, in-depth structural, evolutionary, and functional characterization of maize SWEET transporters was performed for their utility in maize improvement. The mining of SWEET genes in the latest maize genome release (v.5) showed an uneven distribution of 20 ZmSWEETs. The comprehensive structural analyses and docking of ZmSWEETs with four sugars, viz., fructose, galactose, glucose, and sucrose, revealed frequent amino acid residues forming hydrogen (asparagine, valine, serine) and hydrophobic (tryptophan, glycine, and phenylalanine) interactions. Evolutionary analyses of SWEETs showed a mixed lineage with 50-100 % commonality of ortho-groups and -sequences evolved under strong purifying selection (Ka/Ks < 0.5). The duplication analysis showed non-functionalization (ZmSWEET18 in B73) and neo- and sub-functionalization (ZmSWEET3, ZmSWEET6, ZmSWEET9, ZmSWEET19, and ZmSWEET20) events in maize. Functional analyses of ZmSWEET genes through co-expression, in silico expression and qRT-PCR assays showed the relevance of ZmSWEETs expression in regulating drought, heat, and waterlogging stress tolerances in maize. The first ever ZmSWEET-regulatory network revealed 286 direct (ZmSWEET-TF: 140 ZmSWEET-miRNA: 146) and 1226 indirect (TF-TF: 597; TF-miRNA: 629) edges. The present investigation has given new insights into the complex transcriptional and post-transcriptional regulation and the regulatory and functional relevance of ZmSWEETs in assigning stress tolerance in maize.

An improved assembly of the "Cascade" hop (Humulus lupulus) genome uncovers signatures of molecular evolution and refines time of divergence estimates for the Cannabaceae family

We present a chromosome-level assembly of the Cascade hop (Humulus lupulus L. var. lupulus) genome. The hop genome is large (2.8 Gb) and complex, and early attempts at assembly were fragmented. Recent advances have made assembly of the hop genome more tractable, transforming the extent of investigation that can occur. The chromosome-level assembly of Cascade was developed by scaffolding the previously reported Cascade assembly generated with PacBio long-read sequencing and polishing with Illumina short-read DNA sequencing. We developed gene models and repeat annotations and used a controlled bi-parental mapping population to identify significant sex-associated markers. We assessed molecular evolution in gene sequences, gene family expansion and contraction, and time of divergence from Cannabis sativa and other closely related plant species using Bayesian inference. We identified the putative sex chromosome in the female genome based on significant sex-associated markers from the bi-parental mapping population. While the estimate of repeat content (~64%) is similar to the estimate for the hemp genome, syntenic blocks in hop contain a greater percentage of LTRs. Hop is enriched for disease resistance-associated genes in syntenic gene blocks and expanded gene families. The Cascade chromosome-level assembly will inform cultivation strategies and serve to deepen our understanding of the hop genomic landscape, benefiting hop researchers and the Cannabaceae genomics community.

Genome-wide analysis of the Populus trichocarpa laccase gene family and functional identification of PtrLAC23

INTRODUCTION

Biofuel is a kind of sustainable, renewable and environment friendly energy. Lignocellulose from the stems of woody plants is the main raw material for "second generation biofuels". Lignin content limits fermentation yield and is therefore a major obstacle in biofuel production. Plant laccase plays an important role in the final step of lignin formation, which provides a new strategy for us to obtain ideal biofuels by regulating the expression of laccase genes to directly gain the desired lignin content or change the composition of lignin.

METHODS

Multiple sequence alignment and phylogenetic analysis were used to classify PtrLAC genes; sequence features of PtrLACs were revealed by gene structure and motif composition analysis; gene duplication, interspecific collinearity and Ka/Ks analysis were conducted to identify ancient PtrLACs; expression levels of PtrLAC genes were measured by RNA-Seq data and qRT-PCR; domain analysis combine with cis-acting elements prediction together showed the potential function of PtrLACs. Furthermore, Alphafold2 was used to simulate laccase 3D structures, proLAC23::LAC23-eGFP transgenic Populus stem transects were applied to fluorescence observation.

RESULTS

A comprehensive analysis of the P. trichocarpa laccase gene (PtLAC) family was performed. Some ancient PtrLAC genes such as PtrLAC25, PtrLAC19 and PtrLAC41 were identified. Gene structure and distribution of conserved motifs clearly showed sequence characteristics of each PtrLAC. Combining published RNA-Seq data and qRT-PCR analysis, we revealed the expression pattern of PtrLAC gene family. Prediction results of cis-acting elements show that PtrLAC gene regulation was closely related to light. Through above analyses, we selected 5 laccases and used Alphafold2 to simulate protein 3D structures, results showed that PtrLAC23 may be closely related to the lignification. Fluorescence observation of proLAC23::LAC23-eGFP transgenic Populus stem transects and qRT-PCR results confirmed our hypothesis again.

DISCUSSION

In this study, we fully analyzed the Populus trichocarpa laccase gene family and identified key laccase genes related to lignification. These findings not only provide new insights into the characteristics and functions of Populus laccase, but also give a new understanding of the broad prospects of plant laccase in lignocellulosic biofuel production.

Genome-wide identification and expression analysis of two-component system genes in sweet potato (Ipomoea batatas L.)

Two-component system (TCS), which comprises histidine kinases (HKs), histidine phosphotransfer proteins (HPs), and response regulators (RRs), plays essential roles in regulating plant growth, development, and response to various environmental stimuli. TCS genes have been comprehensively identified in various plants, while studies on the genome-wide identification and analysis of TCS in sweet potato were still not reported. Therefore, in this study, a total of 90 TCS members consisting of 20 HK(L)s, 11 HPs, and 59 RRs were identified in the genome of Ipomoea batatas. Furthermore, their gene structures, conserved domains, and phylogenetic relationships were analyzed in detail. Additionally, the gene expression profiles in various organs were analyzed, and response patterns to adverse environmental stresses were investigated. The results showed that these 90 TCS genes were mapped on 15 chromosomes with a notably uneven distribution, and the expansion of TCS genes in sweet potato was attributed to both segmental and tandem duplications. The majority of the TCS genes showed distinct organ-specific expression profiles, especially in three types of roots (stem roots, fibrous roots, tuberous roots). Moreover, most of the TCS genes were either induced or suppressed upon treatment with abiotic stresses (drought, salinity, cold, heat) and exogenous phytohormone abscisic acid (ABA). In addition, the yeast-two hybrid system was used to reveal the HK-HP-RR protein-protein interactions. IbHP1, IbHP2, IbHP4, and IbHP5 could interact with three HKs (IbHK1a, IbHK1b, and IbHK5), and also interact with majority of the type-B RRs (IbRR20-IbRR28), while no interaction affinity was detected for IbHP3. Our systematic analyses could provide insights into the characterization of the TCS genes, and further the development of functional studies in sweet potato.

MeGATAs, functional generalists in interactions between cassava growth and development, and abiotic stresses

The proteins with DNA-binding preference to the consensus DNA sequence (A/T) GATA (A/G) belong to a GATA transcription factor family, with a wide array of biological processes in plants. Cassava (Manihot esculenta) is an important food crop with high production of starch in storage roots. Little was however known about cassava GATA domain-containing genes (MeGATAs). Thirty-six MeGATAs, MeGATA1 to MeGATA36, were found in this study. Some MeGATAs showed a collinear relationship with orthologous genes of Arabidopsis, poplar and potato, rice, maize and sorghum. Eight MeGATA-encoded proteins (MeGATAs) analysed were all localized in the nucleus. Some MeGATAs had potentials of binding ligands and/or enzyme activity. One pair of tandem-duplicated MeGATA17-MeGATA18 and 30 pairs of whole genome-duplicated MeGATAs were found. Fourteen MeGATAs showed low or no expression in the tissues. Nine analysed MeGATAs showed expression responses to abiotic stresses and exogenous phytohormones. Three groups of MeGATA protein interactions were found. Fifty-three miRNAs which can target 18 MeGATAs were identified. Eight MeGATAs were found to target other 292 cassava genes, which were directed to radial pattern formation and phyllome development by gene ontology enrichment, and autophagy by Kyoto Encyclopaedia of Genes and Genomes enrichment. These data suggest that MeGATAs are functional generalists in interactions between cassava growth and development, abiotic stresses and starch metabolism.

Genome-wide identification and characterization of the bZIP gene family and their function in starch accumulation in Chinese chestnut (Castanea mollissima Blume)

The transcription factors of basic leucine zipper (bZIP) family genes play significant roles in stress response as well as growth and development in plants. However, little is known about the bZIP gene family in Chinese chestnut (Castanea mollissima Blume). To better understand the characteristics of bZIPs in chestnut and their function in starch accumulation, a series of analyses were performed including phylogenetic, synteny, co-expression and yeast one-hybrid analyses. Totally, we identified 59 bZIP genes that were unevenly distributed in the chestnut genome and named them CmbZIP01 to CmbZIP59. These CmbZIPs were clustered into 13 clades with clade-specific motifs and structures. A synteny analysis revealed that segmental duplication was the major driving force of expansion of the CmbZIP gene family. A total of 41 CmbZIP genes had syntenic relationships with four other species. The results from the co-expression analyses indicated that seven CmbZIPs in three key modules may be important in regulating starch accumulation in chestnut seeds. Yeast one-hybrid assays showed that transcription factors CmbZIP13 and CmbZIP35 might participate in starch accumulation in the chestnut seed by binding to the promoters of CmISA2 and CmSBE1_2, respectively. Our study provided basic information on CmbZIP genes, which can be utilized in future functional analysis and breeding studies.

Genome-wide characterization of the xyloglucan endotransglucosylase/hydrolase gene family in Solanum lycopersicum L. and gene expression analysis in response to arbuscular mycorrhizal symbiosis

Xyloglucan endotransglucosylase/hydrolases (XTHs) are a glycoside hydrolase protein family involved in the biosynthesis of xyloglucans, with essential roles in the regulation of plant cell wall extensibility. By taking advantage of the whole genome sequence in Solanum lycopersicum, 37 SlXTHs were identified in the present work. SlXTHs were classified into four subfamilies (ancestral, I/II, III-A, III-B) when aligned to XTHs of other plant species. Gene structure and conserved motifs showed similar compositions in each subfamily. Segmental duplication was the primary mechanism accounting for the expansion of SlXTH genes. In silico expression analysis showed that SlXTH genes exhibited differential expression in several tissues. GO analysis and 3D protein structure indicated that all 37 SlXTHs participate in cell wall biogenesis and xyloglucan metabolism. Promoter analysis revealed that some SlXTHs have MeJA- and stress-responsive elements. qRT-PCR expression analysis of nine SlXTHs in leaves and roots of mycorrhizal colonized vs. non-colonized plants showed that eight of these genes were differentially expressed in leaves and four in roots, suggesting that SlXTHs might play roles in plant defense induced by arbuscular mycorrhiza. Our results provide valuable insight into the function of XTHs in S. lycopersicum, in addition to the response of plants to mycorrhizal colonization.

Identification, characterization of Apyrase (APY) gene family in rice (Oryza sativa) and analysis of the expression pattern under various stress conditions

Apyrase (APY) is a nucleoside triphosphate (NTP) diphosphohydrolase (NTPDase) which is a member of the superfamily of guanosine diphosphatase 1 (GDA1)-cluster of differentiation 39 (CD39) nucleoside phosphatase. Under various circumstances like stress, cell growth, the extracellular adenosine triphosphate (eATP) level increases, causing a detrimental influence on cells such as cell growth retardation, ROS production, NO burst, and apoptosis. Apyrase hydrolyses eATP accumulated in the extracellular membrane during stress, wounds, into adenosine diphosphate (ADP) and adenosine monophosphate (AMP) and regulates the stress-responsive pathway in plants. This study was designed for the identification, characterization, and for analysis of APY gene expression in Oryza sativa. This investigation discovered nine APYs in rice, including both endo- and ecto-apyrase. According to duplication event analysis, in the evolution of OsAPYs, a significant role is performed by segmental duplication. Their role in stress control, hormonal responsiveness, and the development of cells is supported by the corresponding cis-elements present in their promoter regions. According to expression profiling by RNA-seq data, the genes were expressed in various tissues. Upon exposure to a variety of biotic as well as abiotic stimuli, including anoxia, drought, submergence, alkali, heat, dehydration, salt, and cold, they showed a differential expression pattern. The expression analysis from the RT-qPCR data also showed expression under various abiotic stress conditions, comprising cold, salinity, cadmium, drought, submergence, and especially heat stress. This finding will pave the way for future in-vivo analysis, unveil the molecular mechanisms of APY genes in stress response, and contribute to the development of stress-tolerant rice varieties.

Genome-Scale Investigation of GARP Family Genes Reveals Their Pivotal Roles in Nutrient Stress Resistance in Allotetraploid Rapeseed

The GARP genes are plant-specific transcription factors (TFs) and play key roles in regulating plant development and abiotic stress resistance. However, few systematic analyses of GARPs have been reported in allotetraploid rapeseed (Brassica napus L.) yet. In the present study, a total of 146 BnaGARP members were identified from the rapeseed genome based on the sequence signature. The BnaGARP TFs were divided into five subfamilies: ARR, GLK, NIGT1/HRS1/HHO, KAN, and PHL subfamilies, and the members within the same subfamilies shared similar exon-intron structures and conserved motif configuration. Analyses of the Ka/Ks ratios indicated that the GARP family principally underwent purifying selection. Several cis-acting regulatory elements, essential for plant growth and diverse biotic and abiotic stresses, were identified in the promoter regions of BnaGARPs. Further, 29 putative miRNAs were identified to be targeting BnaGARPs. Differential expression of BnaGARPs under low nitrate, ammonium toxicity, limited phosphate, deficient boron, salt stress, and cadmium toxicity conditions indicated their potential involvement in diverse nutrient stress responses. Notably, BnaA9.HHO1 and BnaA1.HHO5 were simultaneously transcriptionally responsive to these nutrient stresses in both hoots and roots, which indicated that BnaA9.HHO1 and BnaA1.HHO5 might play a core role in regulating rapeseed resistance to nutrient stresses. Therefore, this study would enrich our understanding of molecular characteristics of the rapeseed GARPs and will provide valuable candidate genes for further in-depth study of the GARP-mediated nutrient stress resistance in rapeseed.

Genome-Wide Identification, Characterization, and Expression Analysis of CHS Gene Family Members in Chrysanthemum nankingense

The chalcone synthase (CHS) gene family catalyzes the first committed step in the biosynthesis of flavonoids and plays key roles in various biological processes in plants. However, systematic studies of the CHS gene family in chrysanthemum remain unknown to date. In this study, 16 CnCHS genes were identified by searching the complete genome sequence of Chrysanthemum nankingense. Most contained two exons and one intron with Chal-sti-synt_N and Chal-sti-synt_C domains. A phylogenetic tree of CnCHSs indicated divergence into three major groups, including I, II, and III. Analyses of the genes and promoters of these genes indicated that there are many cis-acting elements that respond to light, phytohormones, stress, and developmental stages. The CnCHS genes have extensive patterns of expression in various tissues and stages of flower development. Tandemly repeated and segmental repeat genes were expressed at higher levels in different tissues than most of the CnCHS genes that have been identified. CnCHS10 is expressed at higher levels in various flower organs than in vegetative tissues, particularly in disc floret petals and pistils. Our study provides valuable information for the systematic analysis of the CnCHS gene family, which also contributes to further research on flavonoid synthesis and petal colors of chrysanthemum.

Comprehensive Genome-Wide Analysis and Expression Pattern Profiling of PLATZ Gene Family Members in Solanum Lycopersicum L. under Multiple Abiotic Stresses

PLATZ (plant AT-rich sequence and zinc-binding) family proteins with two conserved zinc-dependent DNA-binding motifs are transcription factors specific to the plant kingdom. The functions of PLATZ proteins in growth, development, and adaptation to multiple abiotic stresses have been investigated in various plant species, but their role in tomato has not been explored yet. In the present work, 20 non-redundant Solanum lycopersicum PLATZ (SlPLATZ) genes with three segmentally duplicated gene pairs and four tandemly duplicated gene pairs were identified on eight tomato chromosomes. The comparative modeling and gene ontology (GO) annotations of tomato PLATZ proteins indicated their probable roles in defense response, transcriptional regulation, and protein metabolic processes as well as their binding affinity for various ligands, including nucleic acids, peptides, and zinc. SlPLATZ10 and SlPLATZ17 were only expressed in 1 cm fruits and flowers, respectively, indicating their preferential involvement in the development of these organs. The expression of SlPLATZ1, SlPLATZ12, and SlPLATZ19 was up- or down-regulated following exposure to various abiotic stresses, whereas that of SlPLATZ11 was induced under temperature stresses (i.e., cold and heat stress), revealing their probable function in the abiotic stress tolerance of tomato. Weighted gene co-expression network analysis corroborated the aforementioned findings by spotlighting the co-expression of several stress-associated genes with SlPLATZ genes. Confocal fluorescence microscopy revealed the localization of SlPLATZ−GFP fusion proteins in the nucleus, hinting at their functions as transcription factors. These findings provide a foundation for a better understanding of the structure and function of PLATZ genes and should assist in the selection of potential candidate genes involved in the development and abiotic stress adaptation in tomato.