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

Genomic characterisation and ecological distribution of Mantoniella tinhauana: a novel Mamiellophycean green alga from the Western Pacific

Mamiellophyceae are dominant marine algae in much of the ocean, the most prevalent genera belonging to the order Mamiellales: Micromonas, Ostreococcus and Bathycoccus, whose genetics and global distributions have been extensively studied. Conversely, the genus Mantoniella, despite its potential ecological importance, remains relatively under-characterised. In this study, we isolated and characterised a novel species of Mamiellophyceae, Mantoniella tinhauana, from subtropical coastal waters in the South China Sea. Morphologically, it resembles other Mantoniella species; however, a comparative analysis of the 18S and ITS2 marker genes revealed its genetic distinctiveness. Furthermore, we sequenced and assembled the first genome of Mantoniella tinhauana, uncovering significant differences from previously studied Mamiellophyceae species. Notably, the genome lacked any detectable outlier chromosomes and exhibited numerous unique orthogroups. We explored gene groups associated with meiosis, scale and flagella formation, shedding light on species divergence, yet further investigation is warranted. To elucidate the biogeography of Mantoniella tinhauana, we conducted a comprehensive analysis using global metagenomic read mapping to the newly sequenced genome. Our findings indicate this species exhibits a cosmopolitan distribution with a low-level prevalence worldwide. Understanding the intricate dynamics between Mamiellophyceae and the environment is crucial for comprehending their impact on the ocean ecosystem and accurately predicting their response to forthcoming environmental changes.

Identification and characterization of the CRK gene family in the wheat genome and analysis of their expression profile in response to high temperature-induced male sterility

Cysteine-rich receptor-like kinases (CRKs) play many important roles during plant development, including defense responses under both biotic and abiotic stress, reactive oxygen species (ROS) homeostasis, callose deposition and programmed cell death (PCD). However, there are few studies on the involvement of the CRK family in male sterility due to heat stress in wheat (Triticum aestivum L.). In this study, a genome-wide characterization of the CRK family was performed to investigate the structural and functional attributes of the wheat CRKs in anther sterility caused by heat stress. A total of 95 CRK genes were unevenly distributed on 18 chromosomes, with the most genes distributed on chromosome 2B. Paralogous homologous genes with Ka/Ks ratios less than 1 may have undergone strong purifying selection during evolution and are more functionally conserved. The collinearity analysis results of CRK genes showed that wheat and Arabidopsis (A. thaliana), foxtail millet, Brachypodium distachyon (B. distachyon), and rice have three, 12, 15, and 11 pairs of orthologous genes, respectively. In addition, the results of the network interactions of genes and miRNAs showed that five miRNAs were in the hub of the interactions map, namely tae-miR9657b-5p, tae-miR9780, tae-miR9676-5p, tae-miR164, and tae-miR531. Furthermore, qRT-PCR validation of the six TaCRK genes showed that they play key roles in the development of the mononuclear stage anthers, as all six genes were expressed at highly significant levels in heat-stressed male sterile mononuclear stage anthers compared to normal anthers. We hypothesized that the TaCRK gene is significant in the process of high-temperature-induced sterility in wheat based on the combination of anther phenotypes, paraffin sections, and qRT-PCR data. These results improve our understanding of their relationship.

Genome-Wide Analysis of the SAUR Gene Family and Its Expression Profiles in Response to Salt Stress in Santalum album

The SAUR (small auxin-up RNA) family constitutes a category of genes that promptly respond to the hormone auxin and play a pivotal role in diverse biological processes encompassing plant growth and the response to abiotic stress. Santalum album L., a semi-parasitic evergreen tree, is renowned for its economically valuable essential oils, positioning it among the most prized tree species. In this study, a meticulous identification and comprehensive analysis of 43 SAUR genes was conducted within S. album. Based on phylogenetic relationships, the SaSAUR genes were systematically categorized into five groups. A collinearity analysis revealed intriguing insights, disclosing 14 segmental duplications and 9 tandem duplications within the SaSAUR genes, emphasizing the pivotal role of duplication in the expansion of this gene family. Noteworthy variations in the expression levels of SaSAUR genes were observed by delving into the SaSAUR transcriptome data from various tissues, including leaves, roots, and heartwood, as well as under salt-stress conditions. Notably, SaSAUR08 and SaSAUR13 were significantly upregulated in heartwood compared with roots and leaves, while SaSAUR18 was markedly more expressed in roots compared with heartwood and leaves. Furthermore, SaSAUR27 and SaSAUR28 were found to respond closely to salt stress, hinting at their potential involvement in the salt-stress response mechanism. This research offers a comprehensive investigation of SAUR genes in S. album and establishes a foundation for future exploration of the SAUR gene family, particularly its relation to growth and salt-stress responses.

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.

A chromosome-level haplotype-resolved genome assembly of oriental tobacco budworm (Helicoverpa assulta)

Oriental tobacco budworm (Helicoverpa assulta) and cotton bollworm (Helicoverpa armigera) are two closely related species within the genus Helicoverpa. They have similar appearances and consistent damage patterns, often leading to confusion. However, the cotton bollworm is a typical polyphagous insect, while the oriental tobacco budworm belongs to the oligophagous insects. In this study, we used Nanopore, PacBio, and Illumina platforms to sequence the genome of H. assulta and used Hifiasm to create a haplotype-resolved draft genome. The Hi-C technique helped anchor 33 primary contigs to 32 chromosomes, including two sex chromosomes, Z and W. The final primary haploid genome assembly was approximately 415.19 Mb in length. BUSCO analysis revealed a high degree of completeness, with 99.0% gene coverage in this genome assembly. The repeat sequences constituted 38.39% of the genome assembly, and we annotated 17093 protein-coding genes. The high-quality genome assembly of the oriental tobacco budworm serves as a valuable genetic resource that enhances our comprehension of how they select hosts in a complex odour environment. It will also aid in developing an effective control policy.

A high-quality chromosome-scale genome assembly of blood orange, an important pigmented sweet orange variety

Blood orange (BO) is a rare red-fleshed sweet orange (SWO) with a high anthocyanin content and is associated with numerous health-related benefits. Here, we reported a high-quality chromosome-scale genome assembly for Neixiu (NX) BO, reaching 336.63 Mb in length with contig and scaffold N50 values of 30.6 Mb. Furthermore, 96% of the assembled sequences were successfully anchored to 9 pseudo-chromosomes. The genome assembly also revealed the presence of 37.87% transposon elements and 7.64% tandem repeats, and the annotation of 30,395 protein-coding genes. A high level of genome synteny was observed between BO and SWO, further supporting their genetic similarity. The speciation event that gave rise to the Citrus species predated the duplication event found within them. The genome-wide variation between NX and SWO was also compared. This first high-quality BO genome will serve as a fundamental basis for future studies on functional genomics and genome evolution.

Genome-Wide Identification and Evolutionary Analysis of Receptor-like Kinase Family Genes Provides Insights into Anthracnose Resistance of Dioscorea alata

Dioscorea alata, commonly known as "greater yam", is a vital crop in tropical and subtropical regions of the world, yet it faces significant threats from anthracnose disease, mainly caused by Colletotrichum gloeosporioides. However, exploring disease resistance genes in this species has been challenging due to the difficulty of genetic mapping resulting from the loss of the flowering trait in many varieties. The receptor-like kinase (RLK) gene family represents essential immune receptors in plants. In this study, genomic analysis revealed 467 RLK genes in D. alata. The identified RLKs were distributed unevenly across chromosomes, likely due to tandem duplication events. However, a considerable number of ancient whole-genome or segmental duplications dating back over 100 million years contributed to the diversity of RLK genes. Phylogenetic analysis unveiled at least 356 ancient RLK lineages in the common ancestor of Dioscoreaceae, which differentially inherited and expanded to form the current RLK profiles of D. alata and its relatives. The analysis of cis-regulatory elements indicated the involvement of RLK genes in diverse stress responses. Transcriptome analysis identified RLKs that were up-regulated in response to C. gloeosporioides infection, suggesting their potential role in resisting anthracnose disease. These findings provide novel insights into the evolution of RLK genes in D. alata and their potential contribution to disease resistance.

Genome-Wide Analysis of Transcription Factor R2R3-MYB Gene Family and Gene Expression Profiles during Anthocyanin Synthesis in Common Walnut (Juglans regia L.)

The R2R3-MYB gene family, encoding plant transcriptional regulators, participates in many metabolic pathways of plant physiology and development, including flavonoid metabolism and anthocyanin synthesis. This study proceeded as follows: the JrR2R3-MYB gene family was analyzed genome-wide, and the family members were identified and characterized using the high-quality walnut reference genome "Chandler 2.0". All 204 JrR2R3-MYBs were established and categorized into 30 subgroups via phylogenetic analysis. JrR2R3-MYBs were unevenly distributed over 16 chromosomes. Most JrR2R3-MYBs had similar structures and conservative motifs. The cis-acting elements exhibit multiple functions of JrR2R3-MYBs such as light response, metabolite response, and stress response. We found that the expansion of JrR2R3-MYBs was mainly caused by WGD or segmental duplication events. Ka/Ks analysis indicated that these genes were in a state of negative purifying selection. Transcriptome results suggested that JrR2R3-MYBs were widely entangled in the process of walnut organ development and differentially expressed in different colored varieties of walnuts. Subsequently, we identified 17 differentially expressed JrR2R3-MYBs, 9 of which may regulate anthocyanin biosynthesis based on the results of a phylogenetic analysis. These genes were present in greater expression levels in 'Zijing' leaves than in 'Lvling' leaves, as revealed by the results of qRT-PCR experiments. These results contributed to the elucidation of the functions of JrR2R3-MYBs in walnut coloration. Collectively, this work provides a foundation for exploring the functional characteristics of the JrR2R3-MYBs in walnuts and improving the nutritional value and appearance quality of walnuts.

De novo assembly of the complete mitochondrial genome of pepino (Solanum muricatum) using PacBio HiFi sequencing: insights into structure, phylogenetic implications, and RNA editing

BACKGROUND

Solanum muricatum is an emerging horticultural fruit crop with rich nutritional and antioxidant properties. Although the chromosome-scale genome of this species has been sequenced, its mitochondrial genome sequence has not been reported to date.

RESULTS

PacBio HiFi sequencing was used to assemble the circular mitogenome of S. muricatum, which was 433,466 bp in length. In total, 38 protein-coding, 19 tRNA, and 3 rRNA genes were annotated. The reticulate mitochondrial conformations with multiple junctions were verified by polymerase chain reaction, and codon usage, sequence repeats, and gene migration from chloroplast to mitochondrial genome were determined. A collinearity analysis of eight Solanum mitogenomes revealed high structural variability. Overall, 585 RNA editing sites in protein coding genes were identified based on RNA-seq data. Among them, mttB was the most frequently edited (52 times), followed by ccmB (46 times). A phylogenetic analysis based on the S. muricatum mitogenome and those of 39 other taxa (including 25 Solanaceae species) revealed the evolutionary and taxonomic status of S. muricatum.

CONCLUSIONS

We provide the first report of the assembled and annotated S. muricatum mitogenome. This information will help to lay the groundwork for future research on the evolutionary biology of Solanaceae species. Furthermore, the results will assist the development of molecular breeding strategies for S. muricatum based on the most beneficial agronomic traits of this species.

Genome-wide characterization of the NBLRR gene family provides evolutionary and functional insights into blast resistance in pearl millet (Cenchrus americanus (L.) Morrone)

MAIN CONCLUSION

The investigation is the first report on genome-wide identification and characterization of NBLRR genes in pearl millet. We have shown the role of gene loss and purifying selection in the divergence of NBLRRs in Poaceae lineage and candidate CaNBLRR genes for resistance to Magnaporthe grisea infection. Plants have evolved multiple integral mechanisms to counteract the pathogens' infection, among which plant immunity through NBLRR (nucleotide-binding site, leucine-rich repeat) genes is at the forefront. The genome-wide mining in pearl millet (Cenchrus americanus (L.) Morrone) revealed 146 CaNBLRRs. The variation in the branch length of NBLRRs showed the dynamic nature of NBLRRs in response to evolving pathogen races. The orthology of NBLRRs showed a predominance of many-to-one orthologs, indicating the divergence of NBLRRs in the pearl millet lineage mainly through gene loss events followed by gene gain through single-copy duplications. Further, the purifying selection (Ka/Ks < 1) shaped the expansion of NBLRRs within the lineage of pear millet and other members of Poaceae. Presence of cis-acting elements, viz. TCA element, G-box, MYB, SARE, ABRE and conserved motifs annotated with P-loop, kinase 2, RNBS-A, RNBS-D, GLPL, MHD, Rx-CC and LRR suggests their putative role in disease resistance and stress regulation. The qRT-PCR analysis in pearl millet lines showing contrasting responses to Magnaporthe grisea infection identified CaNBLRR20, CaNBLRR33, CaNBLRR46 CaNBLRR51, CaNBLRR78 and CaNBLRR146 as putative candidates. Molecular docking showed the involvement of three and two amino acid residues of LRR domains forming hydrogen bonds (histidine, arginine and threonine) and salt bridges (arginine and lysine) with effectors. Whereas 14 and 20 amino acid residues of CaNBLRR78 and CaNBLRR20 showed hydrophobic interactions with 11 and 9 amino acid residues of effectors, Mg.00g064570.m01 and Mg.00g006570.m01, respectively. The present investigation gives a comprehensive overview of CaNBLRRs and paves the foundation for their utility in pearl millet resistance breeding through understanding of host-pathogen interactions.

Functional Divergence in Orthologous Transcription Factors: Insights from AtCBF2/3/1 and OsDREB1C

Despite traditional beliefs of orthologous genes maintaining similar functions across species, growing evidence points to their potential for functional divergence. C-repeat binding factors/dehydration-responsive element binding protein 1s (CBFs/DREB1s) are critical in cold acclimation, with their overexpression enhancing stress tolerance but often constraining plant growth. In contrast, a recent study unveiled a distinctive role of rice OsDREB1C in elevating nitrogen use efficiency (NUE), photosynthesis, and grain yield, implying functional divergence within the CBF/DREB1 orthologs across species. Here, we delve into divergent molecular mechanisms of OsDREB1C and AtCBF2/3/1 by exploring their evolutionary trajectories across rice and Arabidopsis genomes, regulatomes, and transcriptomes. Evolutionary scrutiny shows discrete clades for OsDREB1C and AtCBF2/3/1, with the Poaceae-specific DREB1C clade mediated by a transposon event. Genome-wide binding profiles highlight OsDREB1C's preference for GCCGAC compared to AtCBF2/3/1's preference for A/GCCGAC, a distinction determined by R12 in the OsDREB1C AP2/ERF domain. Cross-species multiomic analyses reveal shared gene orthogroups (OGs) and underscore numerous specific OGs uniquely bound and regulated by OsDREB1C, implicated in NUE, photosynthesis, and early flowering, or by AtCBF2/3/1, engaged in hormone and stress responses. This divergence arises from gene gains/losses (∼16.7% to 25.6%) and expression reprogramming (∼62.3% to 66.2%) of OsDREB1C- and AtCBF2/3/1-regulated OGs during the extensive evolution following the rice-Arabidopsis split. Our findings illustrate the regulatory evolution of OsDREB1C and AtCBF2/3/1 at a genomic scale, providing insights on the functional divergence of orthologous transcription factors following gene duplications across species.

Genome-wide identification and expression analysis of the SET domain-containing gene family in potato (Solanum tuberosum L.)

Genes containing the SET domain can catalyse histone lysine methylation, which in turn has the potential to cause changes to chromatin structure and regulation of the transcription of genes involved in diverse physiological and developmental processes. However, the functions of SET domain-containing (StSET) genes in potato still need to be studied. The objectives of our study can be summarized as in silico analysis to (i) identify StSET genes in the potato genome, (ii) systematically analyse gene structure, chromosomal distribution, gene duplication events, promoter sequences, and protein domains, (iii) perform phylogenetic analyses, (iv) compare the SET domain-containing genes of potato with other plant species with respect to protein domains and orthologous relationships, (v) analyse tissue-specific expression, and (vi) study the expression of StSET genes in response to drought and heat stresses. In this study, we identified 57 StSET genes in the potato genome, and the genes were physically mapped onto eleven chromosomes. The phylogenetic analysis grouped these StSET genes into six clades. We found that tandem duplication through sub-functionalisation has contributed only marginally to the expansion of the StSET gene family. The protein domain TDBD (PFAM ID: PF16135) was detected in StSET genes of potato while it was absent in all other previously studied species. This study described three pollen-specific StSET genes in the potato genome. Expression analysis of four StSET genes under heat and drought in three potato clones revealed that these genes might have non-overlapping roles under different abiotic stress conditions and durations. The present study provides a comprehensive analysis of StSET genes in potatoes, and it serves as a basis for further functional characterisation of StSET genes towards understanding their underpinning biological mechanisms in conferring stress tolerance.

AmDHN4, a winter accumulated SKn-type dehydrin from Ammopiptanthus mongolicus, and regulated by AmWRKY45, enhances the tolerance of Arabidopsis to low temperature and osmotic stress

Ammopiptanthus mongolicus, a rare temperate evergreen broadleaf shrub, exhibits remarkable tolerance to low temperature and drought stress in winter. Late embryogenesis abundant (LEA) proteins, a kind of hydrophilic protein with a protective function, play significant roles in enhancing plant tolerance to abiotic stress. In this present study, we analyzed the evolution and expression of LEA genes in A. mongolicus, and investigated the function and regulatory mechanism of dehydrin under abiotic stresses. Evolutionary analysis revealed that 14 AmLEA genes underwent tandem duplication events, and 36 AmLEA genes underwent segmental duplication events Notably, an expansion in SKn-type dehydrins was observed. Expression analysis showed that AmDHN4, a SKn-type dehydrin, was up-regulated in winter and under low temperature and osmotic stresses. Functional analysis showcased that the heterologous expression of the AmDHN4 enhanced the tolerance of yeast and tobacco to low temperature stress. Additionally, the overexpression of AmDHN4 significantly improved the tolerance of transgenic Arabidopsis to low temperature, drought, and osmotic stress. Further investigations identified AmWRKY45, a downstream transcription factor in the jasmonic acid signaling pathway, binding to the AmDHN4 promoter and positively regulating its expression. In summary, these findings contribute to a deeper understanding of the functional and regulatory mechanisms of dehydrin.

Multi-omics analyses provide insights into the evolutionary history and the synthesis of medicinal components of the Chinese wingnut

Chinese wingnut (Pterocarya stenoptera) is a medicinally and economically important tree species within the family Juglandaceae. However, the lack of high-quality reference genome has hindered its in-depth research. In this study, we successfully assembled its chromosome-level genome and performed multi-omics analyses to address its evolutionary history and synthesis of medicinal components. A thorough examination of genomes has uncovered a significant expansion in the Lateral Organ Boundaries Domain gene family among the winged group in Juglandaceae. This notable increase may be attributed to their frequent exposure to flood-prone environments. After further differentiation between Chinese wingnut and Cyclocarya paliurus, significant positive selection occurred on the genes of NADH dehydrogenase related to mitochondrial aerobic respiration in Chinese wingnut, enhancing its ability to cope with waterlogging stress. Comparative genomic analysis revealed Chinese wingnut evolved more unique genes related to arginine synthesis, potentially endowing it with a higher capacity to purify nutrient-rich water bodies. Expansion of terpene synthase families enables the production of increased quantities of terpenoid volatiles, potentially serving as an evolved defense mechanism against herbivorous insects. Through combined transcriptomic and metabolomic analysis, we identified the candidate genes involved in the synthesis of terpenoid volatiles. Our study offers essential genetic resources for Chinese wingnut, unveiling its evolutionary history and identifying key genes linked to the production of terpenoid volatiles.

NnSUS1 encodes a sucrose synthase involved in sugar accumulation in lotus seed cotyledons

Fresh lotus seeds are gaining favor with consumers for their crunchy texture and natural sweetness. However, the intricacies of sugar accumulation in lotus seeds remain elusive, which greatly hinders the quality improvement of fresh lotus seeds. This study endeavors to elucidate this mechanism by identifying and characterizing the sucrose synthase (SUS) gene family in lotus. Comprising five distinct members, namely NnSUS1 to NnSUS5, each gene within this family features a C-terminal glycosyl transferase1 (GT1) domain. Among them, NnSUS1 is the predominately expressed gene, showing high transcript abundance in the floral organs and cotyledons. NnSUS1 was continuously up-regulated from 6 to 18 days after pollination (DAP) in lotus cotyledons. Furthermore, NnSUS1 demonstrates co-expression relationships with numerous genes involved in starch and sucrose metabolism. To investigate the function of NnSUS1, a transient overexpression system was established in lotus cotyledons, which confirmed the gene's contribution to sugar accumulation. Specifically, transient overexpression of NnSUS1 in seed cotyledons leads to a significant increase in the levels of total soluble sugar, including sucrose and fructose. These findings provide valuable theoretical insights for improving sugar content in lotus seeds through molecular breeding methods.

The asymmetric expression of plasma membrane H+-ATPase family genes in response to pulvinus-driven leaf phototropism movement in Vitis vinifera

Phototropism movement is crucial for plants to adapt to various environmental changes. Plant P-type H+-ATPase (HA) plays diverse roles in signal transduction during cell expansion, regulation of cellular osmotic potential and stomatal opening, and circadian movement. Despite numerous studies on the genome-wide analysis of Vitis vinifera, no research has been done on the P-type H+-ATPase family genes, especially concerning pulvinus-driven leaf movement. In this study, 55 VvHAs were identified and classified into nine distinct subgroups (1 to 9). Gene members within the same subgroups exhibit similar features in motif, intron/exon, and protein tertiary structures. Furthermore, four pairs of genes were derived by segmental duplication in grapes. Cis-acting element analysis identified numerous light/circadian-related elements in the promoters of VvHAs. qRT-PCR analysis showed that several genes of subgroup 7 were highly expressed in leaves and pulvinus during leaf movement, especially VvHA14, VvHA15, VvHA16, VvHA19, VvHA51, VvHA52, and VvHA54. Additionally, we also found that the VvHAs genes were asymmetrically expressed on both sides of the extensor and flexor cell of the motor organ, the pulvinus. The expression of VvHAs family genes in extensor cells was significantly higher than that in flexor cells. Overall, this study serves as a foundation for further investigations into the functions of VvHAs and contributes to the complex mechanisms underlying grapevine pulvinus growth and development.

De novo assembly and analysis of Sonneratia ovata genome and population analysis

Mangroves are an important part of coastal and estuarine ecosystems where they serve as nurseries for marine species and prevent coastal erosion. Here we report the genome of Sonneratia ovata, which is a true mangrove that grows in estuarine environments and can tolerate moderate salt exposure. We sequenced the S. ovata genome and assembled it into chromosome-level scaffolds through the use of Hi-C. The genome is 212.3 Mb and contains 12 chromosomes that range in size from 12.2 to 23.2 Mb. Annotation identified 29,829 genes with a BUSCO completeness of 95.9%. We identified salt genes and found copy number expansion of salt genes such as ADP-ribosylation factor 1, and elongation factor 1-alpha. Population analysis identified a low level of genetic variation and a lack of population structure within S. ovata.

Chromosome-scale genomes of Quercus sichourensis and Quercus rex provide insights into the evolution and adaptation of Fagaceae

The Fagaceae, a plant family with a wide distribution and diverse adaptability, has garnered significant interest as a subject of study in plant speciation and adaptation. Meanwhile, certain Fagaceae species are regarded as highly valuable wood resources due to the exceptional quality of their wood. In this study, we present two high-quality, chromosome-scale genome sequences for Quercus sichourensis (848.75 Mb) and Quercus rex (883.46 Mb). Comparative genomics analysis reveals that the difference in the number of plant disease resistance genes and the nonsynonymous and synonymous substitution ratio (Ka/Ks) of protein-coding genes among Fagaceae species are related to different environmental adaptations. Interestingly, most genes related to starch synthesis in the investigated Quercoideae species are located on a single chromosome, as compared to the outgroup species, Fagus sylvatica. Furthermore, resequencing and population analysis of Q. sichourensis and Q. rex reveal that Q. sichourensis has lower genetic diversity and higher deleterious mutations compared to Q. rex. The high-quality, chromosome-level genomes and the population genomic analysis of the critically endangered Q. sichourensis and Q. rex will provide an invaluable resource as well as insights for future study in these two species, even the genus Quercus, to facilitate their conservation.

Evolutionary divergence of subgenomes in common carp provides insights into speciation and allopolyploid success

Hybridization and polyploidization have made great contributions to speciation, heterosis, and agricultural production within plants, but there is still limited understanding and utilization in animals. Subgenome structure and expression reorganization and cooperation post hybridization and polyploidization are essential for speciation and allopolyploid success. However, the mechanisms have not yet been comprehensively assessed in animals. Here, we produced a high-fidelity reference genome sequence for common carp, a typical allotetraploid fish species cultured worldwide. This genome enabled in-depth analysis of the evolution of subgenome architecture and expression responses. Most genes were expressed with subgenome biases, with a trend of transition from the expression of subgenome A during the early stages to that of subgenome B during the late stages of embryonic development. While subgenome A evolved more rapidly, subgenome B contributed to a greater level of expression during development and under stressful conditions. Stable dominant patterns for homoeologous gene pairs both during development and under thermal stress suggest a potential fixed heterosis in the allotetraploid genome. Preferentially expressing either copy of a homoeologous gene at higher levels to confer development and response to stress indicates the dominant effect of heterosis. The plasticity of subgenomes and their shifting of dominant expression during early development, and in response to stressful conditions, provide novel insights into the molecular basis of the successful speciation, evolution, and heterosis of the allotetraploid common carp.

Genome-wide characterization of cytochrome P450 genes reveals the potential roles in fruit ripening and response to cold stress in tomato

Plant cytochrome P450 (CYP) superfamily, the largest enzyme metabolism family, has been identified in many species and plays a vital role in plant development and stress response via secondary metabolite biosynthesis. A comprehensive identification and functional investigation of CYPs in tomato plants would contribute to deeper understanding of their biological significance. In this study, 268 tomato CYP genes were identified and found to be unevenly located on 12 chromosomes. Based on the phylogenetic analysis, these 268 SlCYPs were classed into two distinct clades (A-type and non-A-type) and nine clans, including 48 families. Moreover, 67 tandem and 22 WGD (whole genome duplication)/segmental duplication events were detected, of which 12 SlCYP genes experienced both WGD/segmental and tandem duplication events, indicating that tandem duplication plays a major role in the expansion of the SlCYP family. Besides, 48 pairs containing 41 SlCYP and 44 AtCYP genes were orthologous, while 216 orthologous pairs were obtained between tomato and potato. The expression level of all SlCYP genes in tomato tissues at different development stages was analyzed, and most expressed SlCYPs showed a tissue-specific pattern. Meanwhile, 143 differentially expressed SlCYPs were identified under cold stress. Furthermore, the RT-qPCR results indicated that SlCYPs may be involved in fruit ripening and cold tolerance in tomato seedlings. These findings provide valuable insights into the evolutionary relationships and functional characteristics of SlCYPs, which can be utilized for further investigation of fruit metabolic pathways and cold tolerance in tomato.

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.

A chromosome-scale genome of Rhus chinensis Mill. provides new insights into plant-insect interaction and gallotannins biosynthesis

Rhus chinensis Mill., an economically valuable Anacardiaceae species, is parasitized by the galling aphid Schlechtendalia chinensis, resulting in the formation of the Chinese gallnut (CG). Here, we report a chromosomal-level genome assembly of R. chinensis, with a total size of 389.40 Mb and scaffold N50 of 23.02 Mb. Comparative genomic and transcriptome analysis revealed that the enhanced structure of CG and nutritional metabolism contribute to improving the adaptability of R. chinensis to S. chinensis by supporting CG and galling aphid growth. CG was observed to be abundant in hydrolysable tannins (HT), particularly gallotannin and its isomers. Tandem repeat clusters of dehydroquinate dehydratase/shikimate dehydrogenase (DQD/SDH) and serine carboxypeptidase-like (SCPL) and their homologs involved in HT production were determined as specific to HT-rich species. The functional differentiation of DQD/SDH tandem duplicate genes and the significant contraction in the phenylalanine ammonia-lyase (PAL) gene family contributed to the accumulation of gallic acid and HT while minimizing the production of shikimic acid, flavonoids, and condensed tannins in CG. Furthermore, we identified one UDP glucosyltransferase (UGT84A), three carboxylesterase (CXE), and six SCPL genes from conserved tandem repeat clusters that are involved in gallotannin biosynthesis and hydrolysis in CG. We then constructed a regulatory network of these genes based on co-expression and transcription factor motif analysis. Our findings provide a genomic resource for the exploration of the underlying mechanisms of plant-galling insect interaction and highlight the importance of the functional divergence of tandem duplicate genes in the accumulation of secondary metabolites.

The NCBI Comparative Genome Viewer (CGV) is an interactive visualization tool for the analysis of whole-genome eukaryotic alignments

We report a new visualization tool for analysis of whole-genome assembly-assembly alignments, the Comparative Genome Viewer (CGV) (https://ncbi.nlm.nih.gov/genome/cgv/). CGV visualizes pairwise same-species and cross-species alignments provided by National Center for Biotechnology Information (NCBI) using assembly alignment algorithms developed by us and others. Researchers can examine large structural differences spanning chromosomes, such as inversions or translocations. Users can also navigate to regions of interest, where they can detect and analyze smaller-scale deletions and rearrangements within specific chromosome or gene regions. RefSeq or user-provided gene annotation is displayed where available. CGV currently provides approximately 800 alignments from over 350 animal, plant, and fungal species. CGV and related NCBI viewers are undergoing active development to further meet needs of the research community in comparative genome visualization.

Nudix hydrolase WvNUDX24 is involved in borneol biosynthesis in Wurfbainia villosa

Borneol, camphor, and bornyl acetate are highly promising monoterpenoids widely used in medicine, flavor, food, and chemical applications. Bornyl diphosphate (BPP) serves as a common precursor for the biosynthesis of these monoterpenoids. Although bornyl diphosphate synthase (BPPS) that catalyzes the cyclization of geranyl diphosphate (GPP) to BPP has been identified in multiple plants, the enzyme responsible for the hydrolysis of BPP to produce borneol has not been reported. Here, we conducted in vitro and in vivo functional characterization to identify the Nudix hydrolase WvNUDX24 from W. villosa, which specifically catalyzes the hydrolysis of BPP to generate bornyl phosphate (BP), and then BP forms borneol under the action of phosphatase. Subcellular localization experiments indicated that the hydrolysis of BPP likely occurs in the cytoplasm. Furthermore, site-directed mutagenesis experiments revealed that four critical residues (R84, S96, P98, and G99) for the hydrolysis activity of WvNUDX24. Additionally, the functional identification of phosphatidic acid phosphatase (PAP) demonstrated that WvPAP5 and WvPAP10 were able to hydrolyze geranylgeranyl diphosphate (GGPP) and farnesyl diphosphate (FPP) to generate geranylgeranyl phosphate (GGP) and farnesyl phosphate (FP), respectively, but could not hydrolyze BPP, GPP, and neryl diphosphate (NPP) to produce corresponding monophosphate products. These findings highlight the essential role of WvNUDX24 in the first step of BPP hydrolysis to produce borneol and provide genetic elements for the production of BPP-related terpenoids through plant metabolic engineering and synthetic biology.

Convergent and/or parallel evolution of RNA-binding proteins in angiosperms after polyploidization

Increasing studies suggest that the biased retention of stress-related transcription factors (TFs) after whole-genome duplications (WGDs) could rewire gene transcriptional networks, facilitating plant adaptation to challenging environments. However, the role of posttranscriptional factors (e.g. RNA-binding proteins, RBPs) following WGDs has been largely ignored. Uncovering thousands of RBPs in 21 representative angiosperm species, we integrate genomic, transcriptomic, regulatomic, and paleotemperature datasets to unravel their evolutionary trajectories and roles in adapting to challenging environments. We reveal functional enrichments of RBP genes in stress responses and identify their convergent retention across diverse angiosperms from independent WGDs, coinciding with global cooling periods. Numerous RBP duplicates derived from WGDs are then identified as cold-induced. A significant overlap of 29 orthogroups between WGD-derived and cold-induced RBP genes across diverse angiosperms highlights a correlation between WGD and cold stress. Notably, we unveil an orthogroup (Glycine-rich RNA-binding Proteins 7/8, GRP7/8) and relevant TF duplicates (CCA1/LHY, RVE4/8, CBF2/4, etc.), co-retained in different angiosperms post-WGDs. Finally, we illustrate their roles in rewiring circadian and cold-regulatory networks at both transcriptional and posttranscriptional levels during global cooling. Altogether, we underline the adaptive evolution of RBPs in angiosperms after WGDs during global cooling, improving our understanding of plants surviving periods of environmental turmoil.

The promise and pitfalls of synteny in phylogenomics

Reconstructing the tree of life remains a central goal in biology. Early methods, which relied on small numbers of morphological or genetic characters, often yielded conflicting evolutionary histories, undermining confidence in the results. Investigations based on phylogenomics, which use hundreds to thousands of loci for phylogenetic inquiry, have provided a clearer picture of life's history, but certain branches remain problematic. To resolve difficult nodes on the tree of life, 2 recent studies tested the utility of synteny, the conserved collinearity of orthologous genetic loci in 2 or more organisms, for phylogenetics. Synteny exhibits compelling phylogenomic potential while also raising new challenges. This Essay identifies and discusses specific opportunities and challenges that bear on the value of synteny data and other rare genomic changes for phylogenomic studies. Synteny-based analyses of highly contiguous genome assemblies mark a new chapter in the phylogenomic era and the quest to reconstruct the tree of life.

Decoding the genome of bloodsucking midge Forcipomyia taiwana (Diptera: Ceratopogonidae): Insights into odorant receptor expansion

Biting midges, notably those within the Ceratopogonidae family, have long been recognized for their epidemiological significance, both as nuisances and vectors for disease transmission in vertebrates. Despite their impact, genomic insights into these insects, particularly beyond the Culicoides genus, remain limited. In this study, we assembled the Forcipomyia taiwana (Shiraki) genome, comprising 113 scaffolds covering 130.4 Mbps-with the longest scaffold reaching 7.6 Mbps and an N50 value of 2.6 Mbps-marking a pivotal advancement in understanding the genetic architecture of ceratopogonid biting midges. Phylogenomic analyses reveal a shared ancestry between F. taiwana and Culicoides sonorensis Wirth & Jones, dating back approximately 124 million years, and highlight a dynamic history of gene family expansions and contractions within the Ceratopogonidae family. Notably, a substantial expansion of the odorant receptor (OR) gene family was observed, which is crucial for the chemosensory capabilities that govern biting midges' interactions with their environment, including host seeking and oviposition behaviors. The distribution of OR genes across the F. taiwana genome displays notable clusters on scaffolds, indicating localized tandem gene duplication events. Additionally, several collinear regions were identified, hinting at segmental duplications, inversions, and translocations, contributing to the olfactory system's evolutionary complexity. Among the 156 ORs identified in F. taiwana, 134 are biting midge-specific ORs, distributed across three distinct clades, each exhibiting unique motif features that distinguish them from the others. Through weighted gene co-expression network analysis, we correlated distinct gene modules with sex and reproductive status, laying the groundwork for future investigations into the interplay between gene expression and adaptive behaviors in F. taiwana. In conclusion, our study not only highlights the unique olfactory repertoire of ceratopogonid biting midges but also sets the stage for future studies into the genetic underpinnings of their unique biological traits and ecological strategies.

The high-quality genome of Cryptotaenia japonica and comparative genomics analysis reveals anthocyanin biosynthesis in Apiaceae

Cryptotaenia japonica, a traditional medicinal and edible vegetable crops, is well-known for its attractive flavors and health care functions. As a member of the Apiaceae family, the evolutionary trajectory and biological properties of C. japonica are not clearly understood. Here, we first reported a high-quality genome of C. japonica with a total length of 427 Mb and N50 length 50.76 Mb, was anchored into 10 chromosomes, which confirmed by chromosome (cytogenetic) analysis. Comparative genomic analysis revealed C. japonica exhibited low genetic redundancy, contained a higher percentage of single-cope gene families. The homoeologous blocks, Ks, and collinearity were analyzed among Apiaceae species contributed to the evidence that C. japonica lacked recent species-specific WGD. Through comparative genomic and transcriptomic analyses of Apiaceae species, we revealed the genetic basis of the production of anthocyanins. Several structural genes encoding enzymes and transcription factor genes of the anthocyanin biosynthesis pathway in different species were also identified. The CjANSa, CjDFRb, and CjF3H gene might be the target of Cjaponica_2.2062 (bHLH) and Cjaponica_1.3743 (MYB). Our findings provided a high-quality reference genome of C. japonica and offered new insights into Apiaceae evolution and biology.

Functional analysis of fasciclin-like arabinogalactan in carotenoid synthesis during tomato fruit ripening

Carotenoids are important pigmented nutrients synthesized by tomato fruits during ripening. To reveal the molecular mechanism underlying carotenoid synthesis during tomato fruit ripening, we analyzed carotenoid metabolites and transcriptomes in six development stages of tomato fruits. A total of thirty different carotenoids were detected and quantified in tomato fruits from 10 to 60 DPA. Based on differential gene expression profiles and WGCNA, we explored several genes that were highly significant and negatively correlated with lycopene, all of which encode fasciclin-like arabinogalactan proteins (FLAs). The FLAs are involved in plant signal transduction, however the functional role of these proteins has not been studied in tomato. Genome-wide analysis revealed that cultivated and wild tomato species contained 18 to 22 FLA family members, clustered into four groups, and mainly evolved by means of segmental duplication. The functional characterization of FLAs showed that silencing of SlFLA1, 5, and 13 were found to contribute to the early coloration of tomato fruits, and the expression of carotenoid synthesis-related genes was up-regulated in fruits that changed phenotypically, especially in SlFLA13-silenced plants. Furthermore, the content of multiple carotenoids (including (E/Z)-phytoene, lycopene, γ-carotene, and α-carotene) was significantly increased in SlFLA13-silenced fruits, suggesting that SlFLA13 has a potential inhibitory function in regulating carotenoid synthesis in tomato fruits. The results of the present study broaden the idea of analyzing the biological functions of tomato FLAs and preliminary evidence for the inhibitory role of SlFLA13 in carotenoid synthesis in fruit, providing the theoretical basis and a candidate for improving tomato fruit quality.

New insights into evolution and functional diversification of Camellia sinensis LRR-RLKs

Leucine-rich repeat receptor-like kinases (LRR-RLKs) represent the largest subgroup of receptor-like kinases (RLKs) in plants. While some LRR-RLK members play a role in regulating various plant growth processes related to morphogenesis, disease resistance, and stress response, the functions of most LRR-RLK genes remain unclear. In this study, we identified 397 LRR-RLK genes from the genome of Camellia sinensis and categorized them into 16 subfamilies. Approximately 62% of CsLRR-RLK genes are situated in regions resulting from segmental duplications, suggesting that the expansion of CsLRR-RLK genes is due to segmental duplications. Analysis of gene expression patterns revealed differential expression of CsLRR-RLK genes across different tissues and in response to stress. Furthermore, we demonstrated that CssEMS1 localizes to the cell membrane and can complement Arabidopsis ems1 mutant. This study is the initial in-depth evolutionary examination of LRR-RLKs in tea and provides a basis for future investigations into their functionality.

Supplementary Information

The online version contains supplementary material available at 10.1007/s12298-024-01458-1.

Systematic Analysis of Zinc Finger-Homeodomain Transcription Factors (ZF-HDs) in Barley (Hordeum vulgare L.)

Zinc finger-homeodomain transcription factors (ZF-HDs) are pivotal in regulating plant growth, development, and diverse stress responses. In this study, we found 8 ZF-HD genes in barley genome. Theses eight HvZF-HD genes were located on five chromosomes, and classified into ZHD and MIF subfamily. The collinearity, gene structure, conserved motif, and cis-elements of HvZF-HD genes were also analyzed. Real-time PCR results suggested that the expression of HvZF-HD4, HvZF-HD6, HvZF-HD7 and HvZF-HD8 were up-regulated after hormones (ABA, GA3 and MeJA) or PEG treatments, especially HvZF-HD6 was significantly induced. These results provide useful information of ZF-HD genes to future study aimed at barley breeding.

Characterisation of the MLP genes in peach postharvest cold storage and the regulatory role of PpMLP10 in the chilling stress response

The major latex proteins/ripening-related proteins are a subfamily of the Bet v 1 protein superfamily and are commonly involved in plant development and responses to various stresses. However, the functions of MLPs in the postharvest cold storage of fruits remain uninvestigated. Herein, we identified 30 MLP genes in the peach (Prunus persica) genome that were clustered into three subgroups. Chromosomal location analysis revealed that the PpMLP genes were unevenly distributed on five of the eight peach chromosomes. Synteny analysis of the MLP genes between peach and seven other plant species (five dicotyledons and two monocotyledons) explored their evolutionary characteristics. Furthermore, the PpMLP promoters contained cis-elements for multiple hormones and stress responses. Gene expression analysis revealed that PpMLPs participated in chilling stress responses. Ectopic expression of PpMLP10 in Arabidopsis improved chilling stress tolerance by decreasing membrane damage and maintaining membrane stability. Additional research confirmed that PpWRKY2 participates in PpMLP10-mediated chilling stress by binding to its promoter. Collectively, these results suggest the role of PpMLP10 in enhancing chilling stress tolerance, which is significant for decreasing chilling injury during the postharvest cold storage of peaches.

Genome-wide identification of JAZ gene family in sugarcane and function analysis of ScJAZ1/2 in drought stress response and flowering regulation

The JASMONATE ZIM DOMAIN (JAZ) proteins are a key inhibitors of the jasmonic acid (JA) signaling pathway that play an important role in the regulation of plant growth and development and environmental stress responses. However, there is no systematic identification and functional analysis of JAZ gene family members in sugarcane. In this study, a total of 49 SsJAZ genes were identified from the wild sugarcane species Saccharum spontaneum genome that were unevenly distributed on 13 chromosomes. Phylogenetic analysis showed that all SsJAZ members can be divided into six groups, and most of the SsJAZ genes contained photoreactive and ABA-responsive elements. RNA-seq analysis revealed that SsJAZ1-1/2/3/4 and SsJAZ7-1 were significantly upregulated under drought stress. The transcript level of ScJAZ1 which is the homologous gene of SsJAZ1 in modern sugarcane cultivars was upregulated by JA, PEG, and abscisic acid (ABA). Moreover, ScJAZ1 can interact with three other JAZ proteins to form heterodimers. The spatial and temporal expression analysis showed that SsJAZ2-1/2/3/4 were highly expressed in different tissues and growth stages and during the day-night rhythm between 10:00 and 18:00. Overexpression of ScJAZ2 in Arabidopsis accelerated flowering through activating the expression of AtSOC1, AtFT, and AtLFY. Moreover, the transcription level of ScJAZ2 was about 30-fold in the early-flowering sugarcane variety than that of the non-flowering variety, indicating ScJAZ2 positively regulated flowering. This first systematic analysis of the JAZ gene family and function analysis of ScJAZ1/2 in sugarcane provide key candidate genes and lay the foundation for sugarcane breeding.

De novo genome assembly and population genomics of a shrub tree Barthea barthei (Hance) krass provide insights into the adaptive color variations

Flower color is a classic example of an ecologically important trait under selection in plants. Understanding the genetic mechanisms underlying shifts in flower color can provide key insights into ecological speciation. In this study, we investigated the genetic basis of flower color divergence in Barthea barthei, a shrub tree species exhibiting natural variation in flower color. We assembled a high-quality genome assembly for B. barthei with a contig N50 of 2.39 Mb and a scaffold N50 of 16.21 Mb. The assembly was annotated with 46,430 protein-coding genes and 1,560 non-coding RNAs. Genome synteny analysis revealed two recent tetraploidization events in B. barthei, estimated to have occurred at approximately 17 and 63 million years ago. These tetraploidization events resulted in massive duplicated gene content, with over 70% of genes retained in collinear blocks. Gene family members of the core regulators of the MBW complex were significantly expanded in B. barthei compared to Arabidopsis, suggesting that these duplications may have provided raw genetic material for the evolution of novel regulatory interactions and the diversification of anthocyanin pigmentation. Transcriptome profiling of B. barthei flowers revealed differential expression of 9 transcription factors related to anthocyanin biosynthesis between the two ecotypes. Six of these differentially expressed transcription factors were identified as high-confidence candidates for adaptive evolution based on positive selection signals. This study provides insights into the genetic basis of flower color divergence and the evolutionary mechanisms underlying ecological adaptation in plants.

Over-expression of Medicago Acyl-CoA-binding 2 genes enhance salt and drought tolerance in Arabidopsis

Acyl-CoA-binding proteins (ACBPs) are mainly involved in acyl-CoA ester binding and trafficking in eukaryotic cells, and they function in lipid metabolism, membrane biosynthesis, cellular signaling, stress response, disease resistance, and other biological activities in plants. However, the roles of ACBP family members in Medicago remain unclear. In this study, a total of eight ACBP genes were identified in the genome of Medicago truncatula and Medicago sativa, and they were clustered into four sub-families (Class I-IV). Many cis-acting elements related to abiotic response were identified in the promoter region of these ACBP genes, in particular light-responsive elements. These ACBP genes exhibited distinct expression pattern in various tissues, and the expression level of MtACBP1/MsACBP1 and MtACBP2/MsACBP2 gene pairs were significantly increased under NaCl treatment. Subcellular localization analysis showed that MtACBP1/MsACBP1 and MtACBP2/MsACBP2 were localized in the endoplasmic reticulum of tobacco epidermal cells. Arabidopsis seedlings over-expressing MtACBP2/MsACBP2 displayed increased root length than the wild type under short light, Cu2+, ABA, PEG, and NaCl treatments. Over-expression of MtACBP2/MsACBP2 also significantly enhanced Arabidopsis tolerance under NaCl and PEG treatments in mature plants. Collectively, our study identified salt and drought responsive ACBP genes in Medicago and verified their functions in increasing resistance against salt and drought stresses.

Genomes of multicellular algal sisters to land plants illuminate signaling network evolution

Zygnematophyceae are the algal sisters of land plants. Here we sequenced four genomes of filamentous Zygnematophyceae, including chromosome-scale assemblies for three strains of Zygnema circumcarinatum. We inferred traits in the ancestor of Zygnematophyceae and land plants that might have ushered in the conquest of land by plants: expanded genes for signaling cascades, environmental response, and multicellular growth. Zygnematophyceae and land plants share all the major enzymes for cell wall synthesis and remodifications, and gene gains shaped this toolkit. Co-expression network analyses uncover gene cohorts that unite environmental signaling with multicellular developmental programs. Our data shed light on a molecular chassis that balances environmental response and growth modulation across more than 600 million years of streptophyte evolution.

Identification and expression profiling of GAPDH family genes involved in response to Sclerotinia sclerotiorum infection and phytohormones in Brassica napus

Glyceraldehyde 3-phosphate dehydrogenase (GAPDH) is a crucial enzyme in glycolysis, an essential metabolic pathway for carbohydrate metabolism across all living organisms. Recent research indicates that phosphorylating GAPDH exhibits various moonlighting functions, contributing to plant growth and development, autophagy, drought tolerance, salt tolerance, and bacterial/viral diseases resistance. However, in rapeseed (Brassica napus), the role of GAPDHs in plant immune responses to fungal pathogens remains unexplored. In this study, 28 genes encoding GAPDH proteins were revealed in B. napus and classified into three distinct subclasses based on their protein structural and phylogenetic relationships. Whole-genome duplication plays a major role in the evolution of BnaGAPDHs. Synteny analyses revealed orthologous relationships, identifying 23, 26, and 26 BnaGAPDH genes with counterparts in Arabidopsis, Brassica rapa, and Brassica oleracea, respectively. The promoter regions of 12 BnaGAPDHs uncovered a spectrum of responsive elements to biotic and abiotic stresses, indicating their crucial role in plant stress resistance. Transcriptome analysis characterized the expression profiles of different BnaGAPDH genes during Sclerotinia sclerotiorum infection and hormonal treatment. Notably, BnaGAPDH17, BnaGAPDH20, BnaGAPDH21, and BnaGAPDH22 exhibited sensitivity to S. sclerotiorum infection, oxalic acid, hormone signals. Intriguingly, under standard physiological conditions, BnaGAPDH17, BnaGAPDH20, and BnaGAPDH22 are primarily localized in the cytoplasm and plasma membrane, with BnaGAPDH21 also detectable in the nucleus. Furthermore, the nuclear translocation of BnaGAPDH20 was observed under H2O2 treatment and S. sclerotiorum infection. These findings might provide a theoretical foundation for elucidating the functions of phosphorylating GAPDH.

Exploring the evolution of CHS gene family in plants

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

Identification and characterization of PAL genes involved in the regulation of stem development in Saccharum spontaneum L

BACKGROUND

Saccharum spontaneum L. is a closely related species of sugarcane and has become an important genetic component of modern sugarcane cultivars. Stem development is one of the important factors for affecting the yield, while the molecular mechanism of stem development remains poorly understanding in S. spontaneum. Phenylalanine ammonia-lyase (PAL) is a vital component of both primary and secondary metabolism, contributing significantly to plant growth, development and stress defense. However, the current knowledge about PAL genes in S. spontaneum is still limited. Thus, identification and characterization of the PAL genes by transcriptome analysis will provide a theoretical basis for further investigation of the function of PAL gene in sugarcane.

RESULTS

In this study, 42 of PAL genes were identified, including 26 SsPAL genes from S. spontaneum, 8 ShPAL genes from sugarcane cultivar R570, and 8 SbPAL genes from sorghum. Phylogenetic analysis showed that SsPAL genes were divided into three groups, potentially influenced by long-term natural selection. Notably, 20 SsPAL genes were existed on chromosomes 4 and 5, indicating that they are highly conserved in S. spontaneum. This conservation is likely a result of the prevalence of whole-genome replications within this gene family. The upstream sequence of PAL genes were found to contain conserved cis-acting elements such as G-box and SP1, GT1-motif and CAT-box, which collectively regulate the growth and development of S. spontaneum. Furthermore, quantitative reverse transcription polymerase chain reaction (qRT-PCR) analysis showed that SsPAL genes of stem had a significantly upregulated than that of leaves, suggesting that they may promote the stem growth and development, particularly in the + 6 stem (The sixth cane stalk from the top to down) during the growth stage.

CONCLUSIONS

The results of this study revealed the molecular characteristics of SsPAL genes and indicated that they may play a vital role in stem growth and development of S. spontaneum. Altogether, our findings will promote the understanding of the molecular mechanism of S. spontaneum stem development, and also contribute to the sugarcane genetic improving.

Identification and characterization of ACR gene family in maize for salt stress tolerance

Background

Members of the ACR gene family are commonly involved in various physiological processes, including amino acid metabolism and stress responses. In recent decades, significant progress has been made in the study of ACR genes in plants. However, little is known about their characteristics and function in maize.

Methods

In this study, ACR genes were identified from the maize genome, and their molecular characteristics, gene structure, gene evolution, gene collinearity analysis, cis-acting elements were analyzed. qRT-PCR technology was used to verify the expression patterns of the ZmACR gene family in different tissues under salt stress. In addition, Ectopic expression technique of ZmACR5 in Arabidopsis thaliana was utilized to identify its role in response to salt stress.

Results

A total of 28 ZmACR genes were identified, and their molecular characteristics were extensively described. Two gene pairs arising from segmented replication events were detected in maize, and 18 collinear gene pairs were detected between maize and 3 other species. Through phylogenetic analysis, three subgroups were revealed, demonstrating distinct divergence between monocotyledonous and dicotyledonous plants. Analysis of ZmACR cis-acting elements revealed the optional involvement of ZmACR genes in light response, hormone response and stress resistance. Expression analysis of 8 ZmACR genes under salt treatment clearly revealed their role in the response to salt stress. Ectopic overexpression of ZmACR5 in Arabidopsis notably reduced salt tolerance compared to that of the wild type under salt treatment, suggesting that ZmACR5 has a negative role in the response to salt stress.

Conclusion

Taken together, these findings confirmed the involvement of ZmACR genes in regulating salt stress and contributed significantly to our understanding of the molecular function of ACR genes in maize, facilitating further research in this field.

High-quality genome assembly enables prediction of allele-specific gene expression in hybrid poplar

Poplar (Populus) is a well-established model system for tree genomics and molecular breeding, and hybrid poplar is widely used in forest plantations. However, distinguishing its diploid homologous chromosomes is difficult, complicating advanced functional studies on specific alleles. In this study, we applied a trio-binning design and PacBio high-fidelity long-read sequencing to obtain haplotype-phased telomere-to-telomere genome assemblies for the 2 parents of the well-studied F1 hybrid "84K" (Populus alba × Populus tremula var. glandulosa). Almost all chromosomes, including the telomeres and centromeres, were completely assembled for each haplotype subgenome apart from 2 small gaps on one chromosome. By incorporating information from these haplotype assemblies and extensive RNA-seq data, we analyzed gene expression patterns between the 2 subgenomes and alleles. Transcription bias at the subgenome level was not uncovered, but extensive-expression differences were detected between alleles. We developed machine-learning (ML) models to predict allele-specific expression (ASE) with high accuracy and identified underlying genome features most highly influencing ASE. One of our models with 15 predictor variables achieved 77% accuracy on the training set and 74% accuracy on the testing set. ML models identified gene body CHG methylation, sequence divergence, and transposon occupancy both upstream and downstream of alleles as important factors for ASE. Our haplotype-phased genome assemblies and ML strategy highlight an avenue for functional studies in Populus and provide additional tools for studying ASE and heterosis in hybrids.

Genome-wide analysis of bZIP transcription factors and their expression patterns in response to methyl jasmonate and low-temperature stresses in Platycodon grandiflorus

Background

Platycodon grandiflorus belongs to the genus Platycodon and has many pharmacological effects, such as expectorant, antitussive, and anti-tumor properties. Among transcription factor families peculiar to eukaryotes, the basic leucine zipper (bZIP) family is one of the most important, which exists widely in plants and participates in many biological processes, such as plant growth, development, and stress responses. However, genomic analysis of the bZIP gene family and related stress response genes has not yet been reported in P. grandiflorus.

Methods

P. grandiflorus bZIP (PgbZIP) genes were first identified here, and the phylogenetic relationships and conserved motifs in the PgbZIPs were also performed. Meanwhile, gene structures, conserved domains, and the possible protein subcellular localizations of these PgbZIPs were characterized. Most importantly, the cis-regulatory elements and expression patterns of selected genes exposed to two different stresses were analyzed to provide further information on PgbZIPs potential biological roles in P. grandiflorus upon exposure to environmental stresses.

Conclusions

Forty-six PgbZIPs were identified in P. grandiflorus and divided into nine groups, as displayed in the phylogenetic tree. The results of the chromosomal location and the collinearity analysis showed that forty-six PgbZIP genes were distributed on eight chromosomes, with one tandem duplication event and eleven segmental duplication events identified. Most PgbZIPs in the same phylogenetic group have similar conserved motifs, domains, and gene structures. There are cis-regulatory elements related to the methyl jasmonate (MeJA) response, low-temperature response, abscisic acid response, auxin response, and gibberellin response. Ten PgbZIP genes were selected to study their expression patterns upon exposure to low-temperature and MeJA treatments, and all ten genes responded to these stresses. The real-time quantitative polymerase chain reaction (RT-qPCR) results suggest that the expression levels of most PgbZIPs decreased significantly within 6 h and then gradually increased to normal or above normal levels over the 90 h following MeJA treatment. The expression levels of all PgbZIPs were significantly reduced after 3 h of the low-temperature treatment. These results reveal the characteristics of the PgbZIP family genes and provide valuable information for improving P. grandiflorus's ability to cope with environmental stresses during growth and development.

Differential expression and global analysis of miR156/SQUAMOSA promoter binding-like proteins (SPL) module in oat

SQUAMOSA promoter binding-like proteins (SPLs) are important transcription factors that influence growth phase transition and reproduction in plants. SPLs are targeted by miR156 but the SPL/miR156 module is completely unknown in oat. We identified 28 oat SPL genes (AsSPLs) distributed across all 21 oat chromosomes except for 4C and 6D. The oat- SPL gene family represented six of eight SPL phylogenetic groups, with no AsSPLs in groups 3 and 7. A novel oat miR156 (AsmiR156) family with 21 precursors divided into 7 groups was characterized. A total of 16 AsSPLs were found to be targeted by AsmiR156. Intriguingly, AsSPL3s showed high transcript abundance during early inflorescence (GS-54), as compared to the lower abundance of AsmiR156, indicating their role in reproductive development. Unravelling the SPL/miR156 regulatory hub and alterations in expression patterns of AsSPLs could provide an essential toolbox for genetic improvement in the cultivated oat.

Genomic analysis reveals phylogeny of Zygophyllales and mechanism for water retention of a succulent xerophyte

Revealing the genetic basis for stress-resistant traits in extremophile plants will yield important information for crop improvement. Zygophyllum xanthoxylum, an extant species of the ancient Mediterranean, is a succulent xerophyte that can maintain a favorable water status under desert habitats; however, the genetic basis of this adaptive trait is poorly understood. Furthermore, the phylogenetic position of Zygophyllales, to which Z. xanthoxylum belongs, remains controversial. In this study, we sequenced and assembled the chromosome-level genome of Z. xanthoxylum. Phylogenetic analysis showed that Zygophyllales and Myrtales form a separated taxon as a sister to the clade comprising fabids and malvids, clarifying the phylogenetic position of Zygophyllales at whole-genome scale. Analysis of genomic and transcriptomic data revealed multiple critical mechanisms underlying the efficient osmotic adjustment using Na+ and K+ as "cheap" osmolytes that Z. xanthoxylum has evolved through the expansion and synchronized expression of genes encoding key transporters/channels and their regulators involved in Na+/K+ uptake, transport, and compartmentation. It is worth noting that ZxCNGC1;1 (cyclic nucleotide-gated channels) and ZxCNGC1;2 constituted a previously undiscovered energy-saving pathway for Na+ uptake. Meanwhile, the core genes involved in biosynthesis of cuticular wax also featured an expansion and upregulated expression, contributing to the water retention capacity of Z. xanthoxylum under desert environments. Overall, these findings boost the understanding of evolutionary relationships of eudicots, illustrate the unique water retention mechanism in the succulent xerophyte that is distinct from glycophyte, and thus provide valuable genetic resources for the improvement of stress tolerance in crops and insights into the remediation of sodic lands.

MYB transcription factors encoded by diversified tandem gene clusters cause varied Morella rubra fruit color

Chinese bayberry (Morella rubra) is a fruit tree with a remarkable variation in fruit color, ranging from white to dark red as determined by anthocyanin content. In dark red "Biqi" (BQ), red "Dongkui" (DK), pink "Fenhong" (FH), and white "Shuijing" (SJ), we identified an anthocyanin-related MYB transcription factor-encoding gene cluster of four members, i.e. MrMYB1.1, MrMYB1.2, MrMYB1.3, and MrMYB2. Collinear analysis revealed that the MYB tandem cluster may have occurred in a highly conserved region of many eudicot genomes. Two alleles of MrMYB1.1 were observed; MrMYB1.1-1 (MrMYB1.1n) was a full-length allele and homozygous in "BQ", MrMYB1.1-2 (MrMYB1.1d) was a nonfunctional allele with a single base deletion and homozygous in "SJ", and MrMYB1.1n/MrMYB1.1d were heterozygous in "DK" and "FH". In these four cultivars, expression of MrMYB1.1, MrMYB1.2, and MrMYB2 was enhanced during ripening. Both alleles were equally expressed in MrMYB1.1n/MrMYB1.1d heterozygous cultivars as revealed by a cleaved amplified polymorphic sequence marker. Expression of MrMYB1.3 was restricted to some dark red cultivars only. Functional characterization revealed that MrMYB1.1n and MrMYB1.3 can induce anthocyanin accumulation while MrMYB1.1d, MrMYB1.2, and MrMYB2 cannot. DNA-protein interaction assays indicated that MrMYB1.1n and MrMYB1.3 can directly bind to and activate the promoters of anthocyanin-related genes via interaction with a MYC-like basic helix-loop-helix protein MrbHLH1. We concluded that the specific genotype of MrMYB1.1 alleles, as well as the exclusive expression of MrMYB1.3 in some dark red cultivars, contributes to fruit color variation. The study provides insights into the mechanisms for regulation of plant anthocyanin accumulation by MYB tandem clusters.

Genome-wide identification of Aux/IAA gene family in white clover (Trifolium repens L.) and functional verification of TrIAA18 under different abiotic stress

BACKGROUND

White clover (Trifolium repens L.) is an excellent leguminous cool-season forage with a high protein content and strong nitrogen-fixing ability. Despite these advantages, its growth and development are markedly sensitive to environmental factors. Indole-3-acetic acid (IAA) is the major growth hormone in plants, regulating plant growth, development, and response to adversity. Nevertheless, the specific regulatory functions of Aux/IAA genes in response to abiotic stresses in white clover remain largely unexplored.

RESULTS

In this study, we identified 47 Aux/IAA genes in the white clover genome, which were categorized into five groups based on phylogenetic analysis. The TrIAAs promoter region co-existed with different cis-regulatory elements involved in developmental and hormonal regulation, and stress responses, which may be closely related to their diverse regulatory roles. Collinearity analysis showed that the amplification of the TrIAA gene family was mainly carried out by segmental duplication. White clover Aux/IAA genes showed different expression patterns in different tissues and under different stress treatments. In addition, we performed a yeast two-hybrid analysis to investigate the interaction between white clover Aux/IAA and ARF proteins. Heterologous expression indicated that TrIAA18 could enhance stress tolerance in both yeast and transgenic Arabidopsis thaliana.

CONCLUSION

These findings provide new scientific insights into the molecular mechanisms of growth hormone signaling in white clover and its functional characteristics in response to environmental stress.

The Spartina alterniflora genome sequence provides insights into the salt-tolerance mechanisms of exo-recretohalophytes

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

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

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

The genomic history and global migration of a windborne pest

Many insect pests, including the brown planthopper (BPH), undergo windborne migration that is challenging to observe and track. It remains controversial about their migration patterns and largely unknown regarding the underlying genetic basis. By analyzing 360 whole genomes from around the globe, we clarify the genetic sources of worldwide BPHs and illuminate a landscape of BPH migration showing that East Asian populations perform closed-circuit journeys between Indochina and the Far East, while populations of Malay Archipelago and South Asia undergo one-way migration to Indochina. We further find round-trip migration accelerates population differentiation, with highly diverged regions enriching in a gene desert chromosome that is simultaneously the speciation hotspot between BPH and related species. This study not only shows the power of applying genomic approaches to demystify the migration in windborne migrants but also enhances our understanding of how seasonal movements affect speciation and evolution in insects.

Tomato root specialized metabolites evolved through gene duplication and regulatory divergence within a biosynthetic gene cluster

Tremendous plant metabolic diversity arises from phylogenetically restricted specialized metabolic pathways. Specialized metabolites are synthesized in dedicated cells or tissues, with pathway genes sometimes colocalizing in biosynthetic gene clusters (BGCs). However, the mechanisms by which spatial expression patterns arise and the role of BGCs in pathway evolution remain underappreciated. In this study, we investigated the mechanisms driving acylsugar evolution in the Solanaceae. Previously thought to be restricted to glandular trichomes, acylsugars were recently found in cultivated tomato roots. We demonstrated that acylsugars in cultivated tomato roots and trichomes have different sugar cores, identified root-enriched paralogs of trichome acylsugar pathway genes, and characterized a key paralog required for root acylsugar biosynthesis, SlASAT1-LIKE (SlASAT1-L), which is nested within a previously reported trichome acylsugar BGC. Last, we provided evidence that ASAT1-L arose through duplication of its paralog, ASAT1, and was trichome-expressed before acquiring root-specific expression in the Solanum genus. Our results illuminate the genomic context and molecular mechanisms underpinning metabolic diversity in plants.