The Arabidopsis Information Resource (TAIR): improved gene annotation and new tools

Arabidopsis mRNA export factor MOS11: molecular interactions and role in abiotic stress responses

Transcription and export (TREX) is a multi-subunit complex that links synthesis, processing and export of mRNAs. It interacts with the RNA helicase UAP56 and export factors such as MOS11 and ALYs to facilitate nucleocytosolic transport of mRNAs. Plant MOS11 is a conserved, but sparsely researched RNA-binding export factor, related to yeast Tho1 and mammalian CIP29/SARNP. Using biochemical approaches, the domains of Arabidopsis thaliana MOS11 required for interaction with UAP56 and RNA-binding were identified. Further analyses revealed marked genetic interactions between MOS11 and ALY genes. Cell fractionation in combination with transcript profiling demonstrated that MOS11 is required for export of a subset of mRNAs that are shorter and more GC-rich than MOS11-independent transcripts. The central α-helical domain of MOS11 proved essential for physical interaction with UAP56 and for RNA-binding. MOS11 is involved in the nucleocytosolic transport of mRNAs that are upregulated under stress conditions and accordingly mos11 mutant plants turned out to be sensitive to elevated NaCl concentrations and heat stress. Collectively, our analyses identify functional interaction domains of MOS11. In addition, the results establish that mRNA export is critically involved in the plant response to stress conditions and that MOS11 plays a prominent role at this.

A data-driven genome annotation approach for cassava

Genome annotation files play a critical role in dictating the quality of downstream analyses by providing essential predictions for gene positions and structures. These files are pivotal in decoding the complex information encoded within DNA sequences. Here, we generated experimental data resolving RNA 5'- and 3'-ends as well as full-length RNAs for cassava TME12 sticklings in ambient temperature and cold. We used these data to generate genome annotation files using the TranscriptomeReconstructoR (TR) tool. A careful comparison to high-quality genome annotations suggests that our new TR genome annotations identified additional genes, resolved the transcript boundaries more accurately and identified additional RNA isoforms. We enhanced existing cassava genome annotation files with the information from TR that maintained the different transcript models as RNA isoforms. The resultant merged annotation was subsequently utilized for comprehensive analysis. To examine the effects of genome annotation files on gene expression studies, we compared the detection of differentially expressed genes during cold using the same RNA-seq data but alternative genome annotation files. We found that our merged genome annotation that included cold-specific TR gene models identified about twice as many cold-induced genes. These data indicate that environmentally induced genes may be missing in off-the-shelf genome annotation files. In conclusion, TR offers the opportunity to enhance crop genome annotations with implications for the discovery of differentially expressed candidate genes during plant-environment interactions.

Genome-wide association identifies a BAHD acyltransferase activity that assembles an ester of glucuronosylglycerol and phenylacetic acid

Genome-wide association studies (GWAS) are an effective approach to identify new specialized metabolites and the genes involved in their biosynthesis and regulation. In this study, GWAS of Arabidopsis thaliana soluble leaf and stem metabolites identified alleles of an uncharacterized BAHD-family acyltransferase (AT5G57840) associated with natural variation in three structurally related metabolites. These metabolites were esters of glucuronosylglycerol, with one metabolite containing phenylacetic acid as the acyl component of the ester. Knockout and overexpression of AT5G57840 in Arabidopsis and heterologous overexpression in Nicotiana benthamiana and Escherichia coli demonstrated that it is capable of utilizing phenylacetyl-CoA as an acyl donor and glucuronosylglycerol as an acyl acceptor. We, thus, named the protein Glucuronosylglycerol Ester Synthase (GGES). Additionally, phenylacetyl glucuronosylglycerol increased in Arabidopsis CYP79A2 mutants that overproduce phenylacetic acid and was lost in knockout mutants of UDP-sulfoquinovosyl: diacylglycerol sulfoquinovosyl transferase, an enzyme required for glucuronosylglycerol biosynthesis and associated with glycerolipid metabolism under phosphate-starvation stress. GGES is a member of a well-supported clade of BAHD family acyltransferases that arose by duplication and neofunctionalized during the evolution of the Brassicales within a larger clade that includes HCT as well as enzymes that synthesize other plant-specialized metabolites. Together, this work extends our understanding of the catalytic diversity of BAHD acyltransferases and uncovers a pathway that involves contributions from both phenylalanine and lipid metabolism.

Genome resources for three modern cotton lines guide future breeding efforts

Cotton (Gossypium hirsutum L.) is the key renewable fibre crop worldwide, yet its yield and fibre quality show high variability due to genotype-specific traits and complex interactions among cultivars, management practices and environmental factors. Modern breeding practices may limit future yield gains due to a narrow founding gene pool. Precision breeding and biotechnological approaches offer potential solutions, contingent on accurate cultivar-specific data. Here we address this need by generating high-quality reference genomes for three modern cotton cultivars ('UGA230', 'UA48' and 'CSX8308') and updating the 'TM-1' cotton genetic standard reference. Despite hypothesized genetic uniformity, considerable sequence and structural variation was observed among the four genomes, which overlap with ancient and ongoing genomic introgressions from 'Pima' cotton, gene regulatory mechanisms and phenotypic trait divergence. Differentially expressed genes across fibre development correlate with fibre production, potentially contributing to the distinctive fibre quality traits observed in modern cotton cultivars. These genomes and comparative analyses provide a valuable foundation for future genetic endeavours to enhance global cotton yield and sustainability.

A high-quality haplotype genome of Michelia alba DC reveals differences in methylation patterns and flower characteristics

Michelia alba DC is a highly valuable ornamental plant of the Magnoliaceae family. This evergreen tropical tree commonly grows in Southeast Asia and is adored for its delightful fragrance. Our study assembled the M. alba haplotype genome MC and MM by utilizing Nanopore ultralong reads, Pacbio Hifi long reads and parental second-generation data. Moreover, the first methylation map of Magnoliaceae was constructed based on the methylation site data obtained using Nanopore data. Metabolomic datasets were generated from the flowers of three different species to assess variations in pigment and volatile compound accumulation. Finally, transcriptome data were generated to link genomic, methylation, and morphological patterns to reveal the reasons underlying the differences between M. alba and its parental lines in petal color, flower shape, and fragrance. We found that the AP1 and AP2 genes are crucial in M. alba petal formation, while the 4CL, PAL, and C4H genes control petal color. The data generated in this study serve as a foundation for future physiological and biochemical research on M. alba, facilitate the targeted improvement of M. alba varieties, and offer a theoretical basis for molecular research on Michelia L.

Genome-wide identification of phasiRNAs in Arabidopsis thaliana, and insights into biogenesis, temperature sensitivity, and organ specificity

The knowledge of biogenesis and target regulation of the phased small interfering RNAs (phasiRNAs) needs continuous update, since the phasiRNA loci are dynamically evolved in plants. Here, hundreds of phasiRNA loci of Arabidopsis thaliana were identified in distinct tissues and under different temperature. In flowers, most of the 24-nt loci are RNA-dependent RNA polymerase 2 (RDR2)-dependent, while the 21-nt loci are RDR6-dependent. Among the RDR-dependent loci, a significant portion is Dicer-like 1-dependent, indicating the involvement of microRNAs in their expression. Besides, two TAS candidates were discovered. Some interesting features of the phasiRNA loci were observed, such as the strong strand bias of phasiRNA generation, and the capacity of one locus for producing phasiRNAs by different increments. Both organ specificity and temperature sensitivity were observed for phasiRNA expression. In leaves, the TAS genes are highly activated under low temperature. Several trans-acting siRNA-target pairs are also temperature-sensitive. In many cases, the phasiRNA expression patterns correlate well with those of the processing signals. Analysis of the rRNA-depleted degradome uncovered several phasiRNA loci to be RNA polymerase II-independent. Our results should advance the understanding on phasiRNA biogenesis and regulation in plants.

PSAURON: a tool for assessing protein annotation across a broad range of species

Evaluating the accuracy of protein-coding sequences in genome annotations is a challenging problem for which there is no broadly applicable solution. In this manuscript we introduce PSAURON (Protein Sequence Assessment Using a Reference ORF Network), a novel software tool developed to assess the quality of protein-coding gene annotations. Utilizing a machine learning model trained on a diverse dataset from over 1000 plant and animal genomes, PSAURON assigns a score to coding DNA or protein sequence that reflects the likelihood that the sequence is a genuine protein coding region. PSAURON scores can be used for genome-wide protein annotation assessment as well as the rapid identification of potentially spurious annotated proteins. Validation against established benchmarks demonstrates PSAURON's effectiveness and correlation with recognized measures of protein quality, highlighting its potential use as a general-purpose method to evaluate gene annotation. PSAURON is open source and freely available at https://github.com/salzberg-lab/PSAURON .

One-Sentence Summary

PSAURON is a machine learning-based tool for rapid assessment of protein coding gene annotation.

Combining DNA and protein alignments to improve genome annotation with LiftOn

As the number and variety of assembled genomes continues to grow, the number of annotated genomes is falling behind, particularly for eukaryotes. DNA-based mapping tools help to address this challenge, but they are only able to transfer annotation between closely-related species. Here we introduce LiftOn, a homology-based software tool that integrates DNA and protein alignments to enhance the accuracy of genome-scale annotation and to allow mapping between relatively distant species. LiftOn's protein-centric algorithm considers both types of alignments, chooses optimal open reading frames, resolves overlapping gene loci, and finds additional gene copies where they exist. LiftOn can reliably transfer annotation between genomes representing members of the same species, as we demonstrate on human, mouse, honey bee, rice, and Arabidopsis thaliana. It can further map annotation effectively across species pairs as far apart as mouse and rat or Drosophila melanogaster and D. erecta.

Identification and expression analysis of the Xyloglucan transglycosylase/hydrolase (XTH) gene family under abiotic stress in oilseed (Brassica napus L.)

XTH genes are key genes that regulate the hydrolysis and recombination of XG components and plays role in the structure and composition of plant cell walls. Therefore, clarifying the changes that occur in XTHs during plant defense against abiotic stresses is informative for the study of the plant stress regulatory mechanism mediated by plant cell wall signals. XTH proteins in Arabidopsis thaliana was selected as the seed sequences in combination with its protein structural domains, 80 members of the BnXTH gene family were jointly identified from the whole genome of the Brassica napus ZS11, and analyzed for their encoded protein physicochemical properties, phylogenetic relationships, covariance relationships, and interoperating miRNAs. Based on the transcriptome data, the expression patterns of BnXTHs were analyzed in response to different abiotic stress treatments. The relative expression levels of some BnXTH genes under Al, alkali, salt, and drought treatments after 0, 6, 12 and 24 h were analyzed by using qRT-PCR to explore their roles in abiotic stress tolerance in B. napus. BnXTHs showed different expression patterns in response to different abiotic stress signals, indicating that the response mechanisms of oilseed rape against different abiotic stresses are also different. This paper provides a theoretical basis for clarifying the function and molecular genetic mechanism of the BnXTH gene family in abiotic stress tolerance in rapeseed.

The Arabidopsis Information Resource in 2024

Since 1999, The Arabidopsis Information Resource (www.arabidopsis.org) has been curating data about the Arabidopsis thaliana genome. Its primary focus is integrating experimental gene function information from the peer-reviewed literature and codifying it as controlled vocabulary annotations. Our goal is to produce a "gold standard" functional annotation set that reflects the current state of knowledge about the Arabidopsis genome. At the same time, the resource serves as a nexus for community-based collaborations aimed at improving data quality, access, and reuse. For the past decade, our work has been made possible by subscriptions from our global user base. This update covers our ongoing biocuration work, some of our modernization efforts that contribute to the first major infrastructure overhaul since 2011, the introduction of JBrowse2, and the resource's role in community activities such as organizing the structural reannotation of the genome. For gene function assessment, we used gene ontology annotations as a metric to evaluate: (1) what is currently known about Arabidopsis gene function and (2) the set of "unknown" genes. Currently, 74% of the proteome has been annotated to at least one gene ontology term. Of those loci, half have experimental support for at least one of the following aspects: molecular function, biological process, or cellular component. Our work sheds light on the genes for which we have not yet identified any published experimental data and have no functional annotation. Drawing attention to these unknown genes highlights knowledge gaps and potential sources of novel discoveries.

Assembly and analysis of the genome of Notholithocarpus densiflorus

Tanoak (Notholithocarpus densiflorus) is an evergreen tree in the Fagaceae family found in California and southern Oregon. Historically, tanoak acorns were an important food source for Native American tribes, and the bark was used extensively in the leather tanning process. Long considered a disjunct relictual element of the Asian stone oaks (Lithocarpus spp.), phylogenetic analysis has determined that the tanoak is an example of convergent evolution. Tanoaks are deeply divergent from oaks (Quercus) of the Pacific Northwest and comprise a new genus with a single species. These trees are highly susceptible to "sudden oak death" (SOD), a plant pathogen (Phytophthora ramorum) that has caused widespread deaths of tanoaks. In this study, we set out to assemble the genome and perform comparative studies among a number of individuals that demonstrated varying levels of susceptibility to SOD. First, we sequenced and de novo assembled a draft reference genome of N. densiflorus using cobarcoded library processing methods and an MGI DNBSEQ-G400 sequencer. To increase the contiguity of the final assembly, we also sequenced Oxford Nanopore long reads to 30× coverage. To our knowledge, the draft genome reported here is one of the more contiguous and complete genomes of a tree species published to date, with a contig N50 of ∼1.2 Mb, a scaffold N50 of ∼2.1 Mb, and a complete gene score of 95.5% through BUSCO analysis. In addition, we sequenced 11 genetically distinct individuals and mapped these onto the draft reference genome, enabling the discovery of almost 25 million single nucleotide polymorphisms and ∼4.4 million small insertions and deletions. Finally, using cobarcoded data, we were able to generate a complete haplotype coverage of all 11 genomes.

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.

A pan-genome of 69 Arabidopsis thaliana accessions reveals a conserved genome structure throughout the global species range

Although originally primarily a system for functional biology, Arabidopsis thaliana has, owing to its broad geographical distribution and adaptation to diverse environments, developed into a powerful model in population genomics. Here we present chromosome-level genome assemblies of 69 accessions from a global species range. We found that genomic colinearity is very conserved, even among geographically and genetically distant accessions. Along chromosome arms, megabase-scale rearrangements are rare and typically present only in a single accession. This indicates that the karyotype is quasi-fixed and that rearrangements in chromosome arms are counter-selected. Centromeric regions display higher structural dynamics, and divergences in core centromeres account for most of the genome size variations. Pan-genome analyses uncovered 32,986 distinct gene families, 60% being present in all accessions and 40% appearing to be dispensable, including 18% private to a single accession, indicating unexplored genic diversity. These 69 new Arabidopsis thaliana genome assemblies will empower future genetic research.

Complementing model species with model clades

Model species continue to underpin groundbreaking plant science research. At the same time, the phylogenetic resolution of the land plant tree of life continues to improve. The intersection of these 2 research paths creates a unique opportunity to further extend the usefulness of model species across larger taxonomic groups. Here we promote the utility of the Arabidopsis thaliana model species, especially the ability to connect its genetic and functional resources, to species across the entire Brassicales order. We focus on the utility of using genomics and phylogenomics to bridge the evolution and diversification of several traits across the Brassicales to the resources in Arabidopsis, thereby extending scope from a model species by establishing a "model clade." These Brassicales-wide traits are discussed in the context of both the model species Arabidopsis and the family Brassicaceae. We promote the utility of such a "model clade" and make suggestions for building global networks to support future studies in the model order Brassicales.

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.

Cloning and functional identification of apple LATERAL ORGAN BOUNDARY DOMAIN 3 (LBD3) transcription factor in the regulation of drought and salt stress

MAIN CONCLUSION

Overexpression of MdLBD3 in Arabidopsis reduced sensitivity to salt and drought stresses and was instrumental in promoting early flowering. Salt and drought stresses have serious effects on plant growth. LATERAL ORGAN BOUNDARY DOMAIN (LBD) proteins are a plant-specific transcription factors (TFs) family and play important roles in plants in resisting to abiotic stress. However, about the function of LBDs in apple and other woody plants is little known. In this study, protein sequences of the LBD family TFs in apples were identified which contained conserved LOB domains. The qRT-PCR analysis showed that the MdLBD3 gene was widely expressed in various tissues and organs. The subcellular localization assay showed that the MdLBD3 protein was localized in the nucleus. Ectopic expression of MdLBD3 in Arabidopsis positively regulated its salt and drought resistance, and promoted early flowering. Collectively, these results showed that MdLBD3 improved the abiotic stress resistance, plant growth and development. Overall, this study provided a new gene for breeding that can increase the abiotic stress tolerance in apple.

Phylogenomic profiles of whole-genome duplications in Poaceae and landscape of differential duplicate retention and losses among major Poaceae lineages

Poaceae members shared a whole-genome duplication called rho. However, little is known about the evolutionary pattern of the rho-derived duplicates among Poaceae lineages and implications in adaptive evolution. Here we present phylogenomic/phylotranscriptomic analyses of 363 grasses covering all 12 subfamilies and report nine previously unknown whole-genome duplications. Furthermore, duplications from a single whole-genome duplication were mapped to multiple nodes on the species phylogeny; a whole-genome duplication was likely shared by woody bamboos with possible gene flow from herbaceous bamboos; and recent paralogues of a tetraploid Oryza are implicated in tolerance of seawater submergence. Moreover, rho duplicates showing differential retention among subfamilies include those with functions in environmental adaptations or morphogenesis, including ACOT for aquatic environments (Oryzoideae), CK2β for cold responses (Pooideae), SPIRAL1 for rapid cell elongation (Bambusoideae), and PAI1 for drought/cold responses (Panicoideae). This study presents a Poaceae whole-genome duplication profile with evidence for multiple evolutionary mechanisms that contribute to gene retention and losses.

Fossil-calibrated molecular clock data enable reconstruction of steps leading to differentiated multicellularity and anisogamy in the Volvocine algae

BACKGROUND

Throughout its nearly four-billion-year history, life has undergone evolutionary transitions in which simpler subunits have become integrated to form a more complex whole. Many of these transitions opened the door to innovations that resulted in increased biodiversity and/or organismal efficiency. The evolution of multicellularity from unicellular forms represents one such transition, one that paved the way for cellular differentiation, including differentiation of male and female gametes. A useful model for studying the evolution of multicellularity and cellular differentiation is the volvocine algae, a clade of freshwater green algae whose members range from unicellular to colonial, from undifferentiated to completely differentiated, and whose gamete types can be isogamous, anisogamous, or oogamous. To better understand how multicellularity, differentiation, and gametes evolved in this group, we used comparative genomics and fossil data to establish a geologically calibrated roadmap of when these innovations occurred.

RESULTS

Our ancestral-state reconstructions, show that multicellularity arose independently twice in the volvocine algae. Our chronograms indicate multicellularity evolved during the Carboniferous-Triassic periods in Goniaceae + Volvocaceae, and possibly as early as the Cretaceous in Tetrabaenaceae. Using divergence time estimates we inferred when, and in what order, specific developmental changes occurred that led to differentiated multicellularity and oogamy. We find that in the volvocine algae the temporal sequence of developmental changes leading to differentiated multicellularity is much as proposed by David Kirk, and that multicellularity is correlated with the acquisition of anisogamy and oogamy. Lastly, morphological, molecular, and divergence time data suggest the possibility of cryptic species in Tetrabaenaceae.

CONCLUSIONS

Large molecular datasets and robust phylogenetic methods are bringing the evolutionary history of the volvocine algae more sharply into focus. Mounting evidence suggests that extant species in this group are the result of two independent origins of multicellularity and multiple independent origins of cell differentiation. Also, the origin of the Tetrabaenaceae-Goniaceae-Volvocaceae clade may be much older than previously thought. Finally, the possibility of cryptic species in the Tetrabaenaceae provides an exciting opportunity to study the recent divergence of lineages adapted to live in very different thermal environments.

Halophytes and heavy metals: A multi-omics approach to understand the role of gene and genome duplication in the abiotic stress tolerance of Cakile maritima

PREMISE

The origin of diversity is a fundamental biological question. Gene duplications are one mechanism that provides raw material for the emergence of novel traits, but evolutionary outcomes depend on which genes are retained and how they become functionalized. Yet, following different duplication types (polyploidy and tandem duplication), the events driving gene retention and functionalization remain poorly understood. Here we used Cakile maritima, a species that is tolerant to salt and heavy metals and shares an ancient whole-genome triplication with closely related salt-sensitive mustard crops (Brassica), as a model to explore the evolution of abiotic stress tolerance following polyploidy.

METHODS

Using a combination of ionomics, free amino acid profiling, and comparative genomics, we characterize aspects of salt stress response in C. maritima and identify retained duplicate genes that have likely enabled adaptation to salt and mild levels of cadmium.

RESULTS

Cakile maritima is tolerant to both cadmium and salt treatments through uptake of cadmium in the roots. Proline constitutes greater than 30% of the free amino acid pool in C. maritima and likely contributes to abiotic stress tolerance. We find duplicated gene families are enriched in metabolic and transport processes and identify key transport genes that may be involved in C. maritima abiotic stress tolerance.

CONCLUSIONS

These findings identify pathways and genes that could be used to enhance plant resilience and provide a putative understanding of the roles of duplication types and retention on the evolution of abiotic stress response.

A comprehensive benchmark of graph-based genetic variant genotyping algorithms on plant genomes for creating an accurate ensemble pipeline

BACKGROUND

Although sequencing technologies have boosted the measurement of the genomic diversity of plant crops, it remains challenging to accurately genotype millions of genetic variants, especially structural variations, with only short reads. In recent years, many graph-based variation genotyping methods have been developed to address this issue and tested for human genomes. However, their performance in plant genomes remains largely elusive. Furthermore, pipelines integrating the advantages of current genotyping methods might be required, considering the different complexity of plant genomes.

RESULTS

Here we comprehensively evaluate eight such genotypers in different scenarios in terms of variant type and size, sequencing parameters, genomic context, and complexity, as well as graph size, using both simulated and real data sets from representative plant genomes. Our evaluation reveals that there are still great challenges to applying existing methods to plants, such as excessive repeats and variants or high resource consumption. Therefore, we propose a pipeline called Ensemble Variant Genotyper (EVG) that can achieve better genotyping performance in almost all experimental scenarios and comparably higher genotyping recall and precision even using 5× reads. Furthermore, we demonstrate that EVG is more robust with an increasing number of graphed genomes, especially for insertions and deletions.

CONCLUSIONS

Our study will provide new insights into the development and application of graph-based genotyping algorithms. We conclude that EVG provides an accurate, unbiased, and cost-effective way for genotyping both small and large variations and will be potentially used in population-scale genotyping for large, repetitive, and heterozygous plant genomes.

Genome-Wide Identification and Analysis of the EIN3/EIL Transcription Factor Gene Family in Doubled Haploid (DH) Poplar

Ethylene (ET) is an important phytohormone that regulates plant growth, development and stress responses. The ethylene-insensitive3/ethylene-insensitive3-like (EIN3/EIL) transcription factor family, as a key regulator of the ET signal transduction pathway, plays an important role in regulating the expression of ET-responsive genes. Although studies of EIN3/EIL family members have been completed in many species, their role in doubled haploid (DH) poplar derived from another culture of diploid Populus simonii × P. nigra (donor tree, DT) remains ambiguous. In this study, a total of seven EIN3/EIL gene family members in the DH poplar genome were identified. Basic physical and chemical property analyses of these genes were performed, and these proteins were predicted to be localized to the nucleus. According to the phylogenetic relationship, EIN3/EIL genes were divided into two groups, and the genes in the same group had a similar gene structure and conserved motifs. The expression patterns of EIN3/EIL genes in the apical buds of different DH poplar plants were analyzed based on transcriptome data. At the same time, the expression patterns of PsnEIL1, PsnEIN3, PsnEIL4 and PsnEIL5 genes in different tissues of different DH plants were detected via RT-qPCR, including the apical buds, young leaves, functional leaves, xylem, cambium and roots. The findings presented above indicate notable variations in the expression levels of PsnEIL genes across various tissues of distinct DH plants. Finally, the PsnEIL1 gene was overexpressed in DT, and the transgenic plants showed a dwarf phenotype, indicating that the PsnEIL1 gene was involved in regulating the growth and development of poplar. In this study, the EIN3/EIL gene family of DH poplar was analyzed and functionally characterized, which provides a theoretical basis for the future exploration of the EIN3/EIL gene function.

Genome-wide association studies identifies genetic loci related to fatty acid and branched-chain amino acid metabolism and histone modifications under varying nitrogen treatments in safflower (Carthamus tinctorius)

Effective identification and usage of genetic variation are prerequisites for developing nutrient-efficient cultivars. A collection of 94 safflower (Carthamus tinctorius ) genotypes (G) was investigated for important morphological and photosynthetic traits at four nitrogen (N) treatments. We found significant variation for all the studied traits except chlorophyll b (chl b ) among safflower genotypes, nitrogen treatments and G×N interaction. The examined traits showed a 2.82-50.00% increase in response to N application. Biological yield (BY) reflected a significantly positive correlation with fresh shoot weight (FSW), root length (RL), fresh root weight (FRW) and number of leaves (NOL), while a significantly positive correlation was also observed among carotenoids (C), chlorophyll a (chl a ), chl b and total chlorophyll content (CT) under all treatments. Superior genotypes with respect to plant height (PH), FSW, NOL, RL, FRW and BY were clustered into Group 3, while genotypes with better mean performance regarding chl a , chl b C and CT were clustered into Group 2 as observed in principal component analysis. The identified eight best-performing genotypes could be useful to develop improved nitrogen efficient cultivars. Genome-wide association analysis resulted in 32 marker-trait associations (MTAs) under four treatments. Markers namely DArT-45481731 , DArT-17812864 , DArT-15670279 and DArT-45482737 were found consistent. Protein-protein interaction networks of loci associated with MTAs were related to fatty acid and branched-chain amino acid metabolism and histone modifications.

Mapping protein-RNA binding in plants with individual-nucleotide-resolution UV cross-linking and immunoprecipitation (plant iCLIP2)

Despite crucial roles of RNA-binding proteins (RBPs) in plant physiology and development, methods for determining their transcriptome-wide binding landscape are less developed than those used in other model organisms. Cross-linking and immunoprecipitation (CLIP) methods (based on UV-mediated generation of covalent bonds between RNAs and cognate RBPs in vivo, purification of the cross-linked complexes and identification of the co-purified RNAs by high-throughput sequencing) have been applied mainly in mammalian cells growing in monolayers or in translucent tissue. We have developed plant iCLIP2, an efficient protocol for performing individual-nucleotide-resolution CLIP (iCLIP) in plants, tailored to overcome the experimental hurdles posed by plant tissue. We optimized the UV dosage to efficiently cross-link RNA and proteins in plants and expressed epitope-tagged RBPs under the control of their native promoters in loss-of-function mutants. We select epitopes for which nanobodies are available, allowing stringent conditions for immunopurification of the RNA-protein complexes to be established. To overcome the inherently high RNase content of plant cells, RNase inhibitors are added and the limited RNA fragmentation step is modified. We combine the optimized isolation of RBP-bound RNAs with iCLIP2, a streamlined protocol that greatly enhances the efficiency of library preparation for high-throughput sequencing. Plant researchers with experience in molecular biology and handling of RNA can complete this iCLIP2 protocol in ~5 d. Finally, we describe a bioinformatics workflow to determine targets of Arabidopsis RBPs from iCLIP data, covering all steps from downloading sequencing reads to identifying cross-linking events ( https://github.com/malewins/Plant-iCLIPseq ), and present the R/Bioconductor package BindingSiteFinder to extract reproducible binding sites ( https://bioconductor.org/packages/release/bioc/html/BindingSiteFinder.html ).

Genomic analysis of Nypa fruticans elucidates its intertidal adaptations and early palm evolution

Nypa fruticans (Wurmb), a mangrove palm species with origins dating back to the Late Cretaceous period, is a unique species for investigating long-term adaptation strategies to intertidal environments and the early evolution of palms. Here, we present a chromosome-level genome sequence and assembly for N. fruticans. We integrated the genomes of N. fruticans and other palm family members for a comparative genomic analysis, which confirmed that the common ancestor of all palms experienced a whole-genome duplication event around 89 million years ago, shaping the distinctive characteristics observed in this clade. We also inferred a low mutation rate for the N. fruticans genome, which underwent strong purifying selection and evolved slowly, thus contributing to its stability over a long evolutionary period. Moreover, ancient duplicates were preferentially retained, with critical genes having experienced positive selection, enhancing waterlogging tolerance in N. fruticans. Furthermore, we discovered that the pseudogenization of Early Methionine-labelled 1 (EM1) and EM6 in N. fruticans underly its crypto-vivipary characteristics, reflecting its intertidal adaptation. Our study provides valuable genomic insights into the evolutionary history, genome stability, and adaptive evolution of the mangrove palm. Our results also shed light on the long-term adaptation of this species and contribute to our understanding of the evolutionary dynamics in the palm family.

GeneMiner: A tool for extracting phylogenetic markers from next-generation sequencing data

The advancement of next-generation sequencing (NGS) technologies has been revolutionary for the field of evolutionary biology. This technology has led to an abundance of available genomes and transcriptomes for researchers to mine. Specifically, researchers can mine for various types of molecular markers that are vital for phylogenetic, evolutionary and ecological studies. Numerous tools have been developed to extract these molecular markers from NGS data. However, due to an insufficient number of well-annotated reference genomes for non-model organisms, it remains challenging to obtain these markers accurately and efficiently. Here, we present GeneMiner, an improved and expanded version of our previous tool, Easy353. GeneMiner combines the reference-guided de Bruijn graph assembly with seed self-discovery and greedy extension. Additionally, it includes a verification step using a parameter-bootstrap method to reduce the pitfalls associated with using a relatively distant reference. Our results, using both experimental and simulation data, showed GeneMiner can accurately acquire phylogenetic molecular markers for plants using transcriptomic, genomic and other NGS data. GeneMiner is designed to be user-friendly, fast and memory-efficient. Further, it is compatible with Linux, Windows and macOS. All source codes are publicly available on GitHub (https://github.com/sculab/GeneMiner) and Gitee (https://gitee.com/sculab/GeneMiner) for easy accessibility and transparency.

The miR166 targets CsHDZ3 genes to negatively regulate drought tolerance in tea plant (Camellia sinensis)

Drought is the stressor with a significant adverse impact on the yield stability of tea plants. HD-ZIP III transcription factors (TFs) play important regulatory roles in plant growth, development, and stress responses. However, whether and how HD-ZIP III TFs are involved in drought response and tolerance in tea plants remains unclear. Here, we identified seven HD-ZIP III genes (CsHDZ3-1 to CsHDZ3-7) in tea plant genome. The evolutionary analysis demonstrated that CsHDZ3 members were subjected to purify selection. Subcellular localization analysis revealed that all seven CsHDZ3s located in the nucleus. Yeast self-activation and dual-luciferase reporter assays demonstrated that CsHDZ3-1 to CsHDZ3-4 have trans-activation ability whereas CsHDZ3-5 to CsHDZ3-7 served as transcriptional inhibitors. The qRT-PCR assay showed that all seven CsHDZ3 genes could respond to simulated natural drought stress and polyethylene glycol treatment. Further assays verified that all CsHDZ3 genes can be cleaved by csn-miR166. Overexpression of csn-miR166 inhibited the expression of seven CsHDZ3 genes and weakened drought tolerance of tea leaves. In contrast, suppression of csn-miR166 promoted the expression of seven CsHDZ3 genes and enhanced drought tolerance of tea leaves. These findings established the foundation for further understanding the mechanism of CsHDZ3-miR166 modules' participation in drought responses and tolerance.

The Madagascar palm genome provides new insights on the evolution of Apocynaceae specialized metabolism

Specialized metabolites possess diverse interesting biological activities and some cardenolides- and monoterpene indole alkaloids- (MIAs) derived pharmaceuticals are currently used to treat human diseases such as cancers or hypertension. While these two families of biocompounds are produced by specific subfamilies of Apocynaceae, one member of this medicinal plant family, the succulent tree Pachypodium lamerei Drake (also known as Madagascar palm), does not produce such specialized metabolites. To explore the evolutionary paths that have led to the emergence and loss of cardenolide and MIA biosynthesis in Apocynaceae, we sequenced and assembled the P. lamerei genome by combining Oxford Nanopore Technologies long-reads and Illumina short-reads. Phylogenomics revealed that, among the Apocynaceae whose genomes have been sequenced, the Madagascar palm is so far the species closest to the common ancestor between MIA producers/non-MIA producers. Transposable elements, constituting 72.48% of the genome, emerge as potential key players in shaping genomic architecture and influencing specialized metabolic pathways. The absence of crucial MIA biosynthetic genes such as strictosidine synthase in P. lamerei and non-Rauvolfioideae species hints at a transposon-mediated mechanism behind gene loss. Phylogenetic analysis not only showcases the evolutionary divergence of specialized metabolite biosynthesis within Apocynaceae but also underscores the role of transposable elements in this intricate process. Moreover, we shed light on the low conservation of enzymes involved in the final stages of MIA biosynthesis in the distinct MIA-producing plant families, inferring independent gains of these specialized enzymes along the evolution of these medicinal plant clades. Overall, this study marks a leap forward in understanding the genomic dynamics underpinning the evolution of specialized metabolites biosynthesis in the Apocynaceae family, with transposons emerging as potential architects of genomics restructuring and gene loss.

An Escherichia coli-Based Phosphorylation System for Efficient Screening of Kinase Substrates

Posttranslational modifications (PTMs), particularly phosphorylation, play a pivotal role in expanding the complexity of the proteome and regulating diverse cellular processes. In this study, we present an efficient Escherichia coli phosphorylation system designed to streamline the evaluation of potential substrates for Arabidopsis thaliana plant kinases, although the technology is amenable to any. The methodology involves the use of IPTG-inducible vectors for co-expressing kinases and substrates, eliminating the need for radioactive isotopes and prior protein purification. We validated the system's efficacy by assessing the phosphorylation of well-established substrates of the plant kinase SnRK1, including the rat ACETYL-COA CARBOXYLASE 1 (ACC1) and FYVE1/FREE1 proteins. The results demonstrated the specificity and reliability of the system in studying kinase-substrate interactions. Furthermore, we applied the system to investigate the phosphorylation cascade involving the A. thaliana MKK3-MPK2 kinase module. The activation of MPK2 by MKK3 was demonstrated to phosphorylate the Myelin Basic Protein (MBP), confirming the system's ability to unravel sequential enzymatic steps in phosphorylation cascades. Overall, this E. coli phosphorylation system offers a rapid, cost-effective, and reliable approach for screening potential kinase substrates, presenting a valuable tool to complement the current portfolio of molecular techniques for advancing our understanding of kinase functions and their roles in cellular signaling pathways.

Rice kinase OsMRLK63 contributes to drought tolerance by regulating reactive oxygen species production

Drought is a major adverse environmental factor that plants face in nature but the molecular mechanism by which plants transduce stress signals and further endow themselves with tolerance remains unclear. Malectin/malectin-like domains containing receptor-like kinases (MRLKs) have been proposed to act as receptors in multiple biological signaling pathways, but limited studies show their roles in drought-stress signaling and tolerance. In this study, we demonstrate OsMRLK63 in rice (Oryza sativa L.) functions in drought tolerance by acting as the receptor of 2 rapid alkalization factors, OsRALF45 and OsRALF46. We show OsMRLK63 is a typical receptor-like kinase that positively regulates drought tolerance and reactive oxygen species (ROS) production. OsMRLK63 interacts with and phosphorylates several nicotinamide adenine dinucleotide phosphate (NADPH) oxidases with the primarily phosphorylated site at Ser26 in the N-terminal of RESPIRATORY BURST OXIDASE HOMOLOGUE A (OsRbohA). The application of the 2 small signal peptides (OsRALF45/46) on rice can greatly alleviate the dehydration of plants induced by mimic drought. This function depends on the existence of OsMRLK63 and the NADPH oxidase-dependent ROS production. The 2 RALFs interact with OsMRLK63 by binding to its extracellular domain, suggesting they may act as drought/dehydration signal sensors for the OsMRLK63-mediated process. Our study reveals a OsRALF45/46-OsMRLK63-OsRbohs module which contributes to drought-stress signaling and tolerance in rice.

A transcriptome-wide identification of ATP-binding cassette (ABC) transporters revealed participation of ABCB subfamily in abiotic stress management of Glycyrrhiza glabra L

Transcriptome-wide survey divulged a total of 181 ABC transporters in G. glabra which were phylogenetically classified into six subfamilies. Protein-Protein interactions revealed nine putative GgABCBs (-B6, -B14, -B15, -B25, -B26, -B31, -B40, -B42 &-B44) corresponding to five AtABCs orthologs (-B1, -B4, -B11, -B19, &-B21). Significant transcript accumulation of ABCB6 (31.8 folds), -B14 (147.5 folds), -B15 (17 folds), -B25 (19.7 folds), -B26 (18.31 folds), -B31 (61.89 folds), -B40 (1273 folds) and -B42 (51 folds) was observed under the influence of auxin. Auxin transport-specific inhibitor, N-1-naphthylphthalamic acid, showed its effectiveness only at higher (10 µM) concentration where it down regulated the expression of ABCBs, PINs (PIN FORMED) and TWD1 (TWISTED DWARF 1) genes in shoot tissues, while their expression was seen to enhance in the root tissues. Further, qRT-PCR analysis under various growth conditions (in-vitro, field and growth chamber), and subjected to abiotic stresses revealed differential expression implicating role of ABCBs in stress management. Seven of the nine genes were shown to be involved in the stress physiology of the plant. GgABCB6, 15, 25 and ABCB31 were induced in multiple stresses, while GgABCB26, 40 & 42 were exclusively triggered under drought stress. No study pertaining to the ABC transporters from G. glabra is available till date. The present investigation will give an insight to auxin transportation which has been found to be associated with plant growth architecture; the knowledge will help to understand the association between auxin transportation and plant responses under the influence of various conditions.

Identification and expression analyses of B3 genes reveal lineage-specific evolution and potential roles of REM genes in pepper

BACKGROUND

The B3 gene family, one of the largest plant-specific transcription factors, plays important roles in plant growth, seed development, and hormones. However, the B3 gene family, especially the REM subfamily, has not been systematically and functionally studied.

RESULTS

In this study, we performed genome-wide re-annotation of B3 genes in five Solanaceae plants, Arabidopsis thaliana, and Oryza sativa, and finally predicted 1,039 B3 genes, including 231 (22.2%) newly annotated genes. We found a striking abundance of REM genes in pepper species (Capsicum annuum, Capsicum baccatum, and Capsicum chinense). Comparative motif analysis revealed that REM and other subfamilies (ABI3/VP1, ARF, RAV, and HSI) consist of different amino acids. We verified that the large number of REM genes in pepper were included in the specific subgroup (G8) through the phylogenetic analysis. Chromosome location and evolutionary analyses suggested that the G8 subgroup genes evolved mainly via a pepper-specific recent tandem duplication on chromosomes 1 and 3 after speciation between pepper and other Solanaceae. RNA-seq analyses suggested the potential functions of REM genes under salt, heat, cold, and mannitol stress conditions in pepper (C. annuum).

CONCLUSIONS

Our study provides evolutionary and functional insights into the REM gene family in pepper.

Genome-wide identification and in silico characterization of major RNAi gene families in date palm (Phoenix dactylifera)

BACKGROUND

Dates contain various minerals that are essential for good health. The major RNA interference (RNAi) gene families play a vital role in plant growth and development by controlling the expression of protein-coding genes against different biotic and abiotic stresses. However, these gene families for date palm are not yet studied. Therefore, this study has explored major RNAi genes and their characteristics in date palm.

RESULTS

We have identified 4 PdDCLs, 7 PdAGOs, and 3 PdRDRs as RNAi proteins from the date palm genome by using AtRNAi genes as query sequences in BLASTp search. Domain analysis of predicted RNAi genes has revealed the Helicase_C, Dicer_dimer, PAZ, RNase III, and Piwi domains that are associated with the gene silencing mechanisms. Most PdRNAi proteins have been found in the nucleus and cytosol associated with the gene silencing actions. The gene ontology (GO) enrichment analysis has revealed some important GO terms including RNA interference, dsRNA fragmentation, and ribonuclease_III activity that are related to the protein-coding gene silencing mechanisms. Gene regulatory network (GRN) analysis has identified PAZ and SNF2 as the transcriptional regulators of PdRNAi genes. Top-ranked 10 microRNAs including Pda-miR156b, Pda-miR396a, Pda-miR166a, Pda-miR167d, and Pda-miR529a have been identified as the key post-transcriptional regulators of PdRNAi genes that are associated with different biotic/abiotic stresses. The cis-acting regulatory element analysis of PdRNAi genes has detected some vital cis-acting elements including ABRE, MBS, MYB, MYC, Box-4, G-box, I-box, and STRE that are linked with different abiotic stresses.

CONCLUSION

The results of this study might be valuable resources for the improvement of different characteristics in date palm by further studies in wet-lab.

Genome-wide identification and expression analysis of NF-Y gene family in tobacco (Nicotiana tabacum L.)

Nuclear factor Y (NF-Y) gene family is an important transcription factor composed of three subfamilies of NF-YA, NF-YB and NF-YC, which is involved in plant growth, development and stress response. In this study, 63 tobacco NF-Y genes (NtNF-Ys) were identified in Nicotiana tabacum L., including 17 NtNF-YAs, 30 NtNF-YBs and 16 NtNF-YCs. Phylogenetic analysis revealed ten pairs of orthologues from tomato and tobacco and 25 pairs of paralogues from tobacco. The gene structure of NtNF-YAs exhibited similarities, whereas the gene structure of NtNF-YBs and NtNF-YCs displayed significant differences. The NtNF-Ys of the same subfamily exhibited a consistent distribution of motifs and protein 3D structure. The protein interaction network revealed that NtNF-YC12 and NtNF-YC5 exhibited the highest connectivity. Many cis-acting elements related to light, stress and hormone response were found in the promoter of NtNF-Ys. Transcriptome analysis showed that more than half of the NtNF-Y genes were expressed in all tissues, and NtNF-YB9/B14/B15/B16/B17/B29 were specifically expressed in roots. A total of 15, 12, 5, and 6 NtNF-Y genes were found to respond to cold, drought, salt, and alkali stresses, respectively. The results of this study will lay a foundation for further study of NF-Y genes in tobacco and other Solanaceae plants.

Plant genomic resources at National Genomics Data Center: assisting in data-driven breeding applications

Genomic data serve as an invaluable resource for unraveling the intricacies of the higher plant systems, including the constituent elements within and among species. Through various efforts in genomic data archiving, integrative analysis and value-added curation, the National Genomics Data Center (NGDC), which is a part of the China National Center for Bioinformation (CNCB), has successfully established and currently maintains a vast amount of database resources. This dedicated initiative of the NGDC facilitates a data-rich ecosystem that greatly strengthens and supports genomic research efforts. Here, we present a comprehensive overview of central repositories dedicated to archiving, presenting, and sharing plant omics data, introduce knowledgebases focused on variants or gene-based functional insights, highlight species-specific multiple omics database resources, and briefly review the online application tools. We intend that this review can be used as a guide map for plant researchers wishing to select effective data resources from the NGDC for their specific areas of study.

Supplementary Information

The online version contains supplementary material available at 10.1007/s42994-023-00134-4.

Transposable elements contribute to the establishment of the glycine shuttle in Brassicaceae species

C3 -C4 intermediate photosynthesis has evolved at least five times convergently in the Brassicaceae, despite this family lacking bona fide C4 species. The establishment of this carbon concentrating mechanism is known to require a complex suite of ultrastructural modifications, as well as changes in spatial expression patterns, which are both thought to be underpinned by a reconfiguration of existing gene-regulatory networks. However, to date, the mechanisms which underpin the reconfiguration of these gene networks are largely unknown. In this study, we used a pan-genomic association approach to identify genomic features that could confer differential gene expression towards the C3 -C4 intermediate state by analysing eight C3 species and seven C3 -C4 species from five independent origins in the Brassicaceae. We found a strong correlation between transposable element (TE) insertions in cis-regulatory regions and C3 -C4 intermediacy. Specifically, our study revealed 113 gene models in which the presence of a TE within a gene correlates with C3 -C4 intermediate photosynthesis. In this set, genes involved in the photorespiratory glycine shuttle are enriched, including the glycine decarboxylase P-protein whose expression domain undergoes a spatial shift during the transition to C3 -C4 photosynthesis. When further interrogating this gene, we discovered independent TE insertions in its upstream region which we conclude to be responsible for causing the spatial shift in GLDP1 gene expression. Our findings hint at a pivotal role of TEs in the evolution of C3 -C4 intermediacy, especially in mediating differential spatial gene expression.

Strong prevalence of light regime-specific QTL in Arabidopsis detected using automated high-throughput phenotyping in fluctuating or constant light

Plants have evolved and adapted under dynamic environmental conditions, particularly to fluctuating light, but plant research has often focused on constant growth conditions. To quantitatively asses the adaptation to fluctuating light, a panel of 384 natural Arabidopsis thaliana accessions was analyzed in two parallel independent experiments under fluctuating and constant light conditions in an automated high-throughput phenotyping system upgraded with supplemental LEDs. While the integrated daily photosynthetically active radiation was the same under both light regimes, plants in fluctuating light conditions accumulated significantly less biomass and had lower leaf area during their measured vegetative growth than plants in constant light. A total of 282 image-derived architectural and/or color-related traits at six common time points, and 77 photosynthesis-related traits from one common time point were used to assess their associations with genome-wide natural variation for both light regimes. Out of the 3000 significant marker-trait associations (MTAs) detected, only 183 (6.1%) were common for fluctuating and constant light conditions. The prevalence of light regime-specific QTL indicates a complex adaptation. Genes in linkage disequilibrium with fluctuating light-specific MTAs with an adjusted repeatability value >0.5 were filtered for gene ontology terms containing "photo" or "light", yielding 15 selected candidates. The candidate genes are involved in photoprotection, PSII maintenance and repair, maintenance of linear electron flow, photorespiration, phytochrome signaling, and cell wall expansion, providing a promising starting point for further investigations into the response of Arabidopsis thaliana to fluctuating light conditions.

Sequence similarity networks bear out hierarchical relationships of green cytochrome P450

Land plants have diversified enzyme families. One of the most prominent is the cytochrome P450 (CYP or CYP450) family. With over 443,000 CYP proteins sequenced across the tree of life, CYPs are ubiquitous in archaea, bacteria, and eukaryotes. Here, we focused on land plants and algae to study the role of CYP diversification. CYPs, acting as monooxygenases, catalyze hydroxylation reactions crucial for specialized plant metabolic pathways, including detoxification and phytohormone production; the CYPome consists of one enormous superfamily that is divided into clans and families. Their evolutionary history speaks of high substrate promiscuity; radiation and functional diversification have yielded numerous CYP families. To understand the evolutionary relationships within the CYPs, we employed sequence similarity network analyses. We recovered distinct clusters representing different CYP families, reflecting their diversified sequences that we link to the prediction of functionalities. Hierarchical clustering and phylogenetic analysis further elucidated relationships between CYP clans, uncovering their shared deep evolutionary history. We explored the distribution and diversification of CYP subfamilies across plant and algal lineages, uncovering novel candidates and providing insights into the evolution of these enzyme families. This identified unexpected relationships between CYP families, such as the link between CYP82 and CYP74, shedding light on their roles in plant defense signaling pathways. Our approach provides a methodology that brings insights into the emergence of new functions within the CYP450 family, contributing to the evolutionary history of plants and algae. These insights can be further validated and implemented via experimental setups under various external conditions.

Proteome plasticity during Physcomitrium patens spore germination - from the desiccated phase to heterotrophic growth and reconstitution of photoautotrophy

The establishment of moss spores is considered a milestone in plant evolution. They harbor protein networks underpinning desiccation tolerance and accumulation of storage compounds that can be found already in algae and that are also utilized in seeds and pollen. Furthermore, germinating spores must produce proteins that drive the transition through heterotrophic growth to the autotrophic plant. To get insight into the plasticity of this proteome, we investigated it at five timepoints of moss (Physcomitrium patens) spore germination and in protonemata and gametophores. The comparison to previously published Arabidopsis proteome data of seedling establishment showed that not only the proteomes of spores and seeds are functionally related, but also the proteomes of germinating spores and young seedlings. We observed similarities with regard to desiccation tolerance, lipid droplet proteome composition, control of dormancy, and β-oxidation and the glyoxylate cycle. However, there were also striking differences. For example, spores lacked any obvious storage proteins. Furthermore, we did not detect homologs to the main triacylglycerol lipase in Arabidopsis seeds, SUGAR DEPENDENT1. Instead, we discovered a triacylglycerol lipase of the oil body lipase family and a lipoxygenase as being the overall most abundant proteins in spores. This finding indicates an alternative pathway for triacylglycerol degradation via oxylipin intermediates in the moss. The comparison of spores to Nicotiana tabacum pollen indicated similarities for example in regards to resistance to desiccation and hypoxia, but the overall developmental pattern did not align as in the case of seedling establishment and spore germination.

MicroRNAs constitute an additional layer in plant response to simultaneous bio- and abiotic stresses as exemplified by UV-B radiation and flg22-treatment on Arabidopsis thaliana

Plants are confronted with various environmental stresses and develop sophisticated adaptive mechanisms. Our previous work demonstrated that the crosstalk of flg22 and ultraviolet (UV)-B-induced signalling cascades reprograms the expression of flavonol pathway genes (FPGs), benefiting plant defence responses. Although several transcription factors have been identified to be involved in this crosstalk, the underlying mechanism is largely unclear. Here, we analyzed microRNAs (miRNAs) and identified 126, 129 and 113 miRNAs with altered abundances compared to untreated control in flg22-, UV-B- and flg22/UV-B-treated seedlings, respectively. Two distinct modules were identified: The first consists of 10 miRNAs repressed by UV-B but up-regulated by flg22, and the second with five miRNAs repressed by flg22 but up-regulated by UV-B. In Arabidopsis, the knockdown of miR858a, a representative of module I, increased the abundance of CHS (a marker gene for FPGs), whereas its overexpression reduced CHS. Conversely, knockout of miR164b from module II decreased CHS and its overexpression increased CHS transcript levels. These data suggest a decisive role of miRNAs in the crosstalk. In the next, we described the interaction between miR858a and its target MYB111 (a positive regulator of FPGs) from module I in detail. We showed that MYB111 was profoundly post-transcriptionally regulated by miR858a during the crosstalk, whose expression was specifically but antagonistically controlled by UVR8- and FLS2-mediated signallings. Moreover, transcriptional monitoring using the GUS reporter gene demonstrates that miRNA-mediated posttranscriptional regulation is the main driving force in reprogramming the expression of FPGs and regulates plant adaptation to multiple concurrent environmental stresses.

Innovative computational tools provide new insights into the polyploid wheat genome

Bread wheat (Triticum aestivum) is an important crop and serves as a significant source of protein and calories for humans, worldwide. Nevertheless, its large and allopolyploid genome poses constraints on genetic improvement. The complex reticulate evolutionary history and the intricacy of genomic resources make the deciphering of the functional genome considerably more challenging. Recently, we have developed a comprehensive list of versatile computational tools with the integration of statistical models for dissecting the polyploid wheat genome. Here, we summarize the methodological innovations and applications of these tools and databases. A series of step-by-step examples illustrates how these tools can be utilized for dissecting wheat germplasm resources and unveiling functional genes associated with important agronomic traits. Furthermore, we outline future perspectives on new advanced tools and databases, taking into consideration the unique features of bread wheat, to accelerate genomic-assisted wheat breeding.

Genome-wide identification and characterization of FAD family genes in barley

Fatty acid desaturases (FADs) play pivotal roles in determining plant stress tolerance. Barley is the most salt-tolerant cereal crop. In this study, we performed genome-wide identification and characterization analysis of the FAD gene family in barley (Hordeum vulgare). A total of 24 HvFADs were identified and divided into four subfamilies based on their amino acid sequence similarity. HvFADs unevenly distributed on six of seven barley chromosomes, and three clusters of HvFADs mainly occurred on the chromosome 2, 3 and 6. Segmental duplication events were found to be a main cause for the HvFAD gene family expansion. The same HvFAD subfamily showed the relatively consistent exon-intron composition and conserved motifs of HvFADs. Cis-element analysis in HvFAD promoters indicated that the expression of HvFADs may be subject to complex regulation, especially stress-responsive elements that may involve in saline-alkaline stress response. Combined transcriptomic data with quantitative experiments, at least five HvFADs highly expressed in roots under salt or alkali treatment, suggesting they may participate in saline or alkaline tolerance in barley. This study provides novel and valuable insights for underlying salt/alkali-tolerant mechanisms in barley.

GmGLU1 and GmRR4 contribute to iron deficiency tolerance in soybean

Iron deficiency chlorosis (IDC) is a form of abiotic stress that negatively impacts soybean yield. In a previous study, we demonstrated that the historical IDC quantitative trait locus (QTL) on soybean chromosome Gm03 was composed of four distinct linkage blocks, each containing candidate genes for IDC tolerance. Here, we take advantage of virus-induced gene silencing (VIGS) to validate the function of three high-priority candidate genes, each corresponding to a different linkage block in the Gm03 IDC QTL. We built three single-gene constructs to target GmGLU1 (GLUTAMATE SYNTHASE 1, Glyma.03G128300), GmRR4 (RESPONSE REGULATOR 4, Glyma.03G130000), and GmbHLH38 (beta Helix Loop Helix 38, Glyma.03G130400 and Glyma.03G130600). Given the polygenic nature of the iron stress tolerance trait, we also silenced the genes in combination. We built two constructs targeting GmRR4+GmGLU1 and GmbHLH38+GmGLU1. All constructs were tested on the iron-efficient soybean genotype Clark grown in iron-sufficient conditions. We observed significant decreases in soil plant analysis development (SPAD) measurements using the GmGLU1 construct and both double constructs, with potential additive effects in the GmRR4+GmGLU1 construct. Whole genome expression analyses (RNA-seq) revealed a wide range of affected processes including known iron stress responses, defense and hormone signaling, photosynthesis, and cell wall structure. These findings highlight the importance of GmGLU1 in soybean iron stress responses and provide evidence that IDC is truly a polygenic trait, with multiple genes within the QTL contributing to IDC tolerance. Finally, we conducted BLAST analyses to demonstrate that the Gm03 IDC QTL is syntenic across a broad range of plant species.

"Genome-wide identification of bZIP gene family in Pearl millet and transcriptional profiling under abiotic stress, phytohormonal treatments; and functional characterization of PgbZIP9"

Abiotic stresses are major constraints in crop production, and are accountable for more than half of the total crop loss. Plants overcome these environmental stresses using coordinated activities of transcription factors and phytohormones. Pearl millet an important C4 cereal plant having high nutritional value and climate resilient features is grown in marginal lands of Africa and South-East Asia including India. Among several transcription factors, the basic leucine zipper (bZIP) is an important TF family associated with diverse biological functions in plants. In this study, we have identified 98 bZIP family members (PgbZIP) in pearl millet. Phylogenetic analysis divided these PgbZIP genes into twelve groups (A-I, S, U and X). Motif analysis has shown that all the PgbZIP proteins possess conserved bZIP domains and the exon-intron organization revealed conserved structural features among the identified genes. Cis-element analysis, RNA-seq data analysis, and real-time expression analysis of PgbZIP genes suggested the potential role of selected PgbZIP genes in growth/development and abiotic stress responses in pearl millet. Expression profiling of selected PgbZIPs under various phytohormones (ABA, SA and MeJA) treatment showed differential expression patterns of PgbZIP genes. Further, PgbZIP9, a homolog of AtABI5 was found to localize in the nucleus and modulate gene expression in pearl millet under stresses. Our present findings provide a better understanding of bZIP genes in pearl millet and lay a good foundation for the further functional characterization of multi-stress tolerant PgbZIP genes, which could become efficient tools for crop improvement.

Interplay between coding and non-coding regulation drives the Arabidopsis seed-to-seedling transition

Translation of seed stored mRNAs is essential to trigger germination. However, when RNAPII re-engages RNA synthesis during the seed-to-seedling transition has remained in question. Combining csRNA-seq, ATAC-seq and smFISH in Arabidopsis thaliana we demonstrate that active transcription initiation is detectable during the entire germination process. Features of non-coding regulation such as dynamic changes in chromatin accessible regions, antisense transcription, as well as bidirectional non-coding promoters are widespread throughout the Arabidopsis genome. We show that sensitivity to exogenous ABSCISIC ACID (ABA) during germination depends on proximal promoter accessibility at ABA-responsive genes. Moreover, we provide genetic validation of the existence of divergent transcription in plants. Our results reveal that active enhancer elements are transcribed producing non-coding enhancer RNAs (eRNAs) as widely documented in metazoans. In sum, this study defining the extent and role of coding and non-coding transcription during key stages of germination expands our understanding of transcriptional mechanisms underlying plant developmental transitions.

Genomic evidence for rediploidization and adaptive evolution following the whole-genome triplication

Whole-genome duplication (WGD), or polyploidy, events are widespread and significant in the evolutionary history of angiosperms. However, empirical evidence for rediploidization, the major process where polyploids give rise to diploid descendants, is still lacking at the genomic level. Here we present chromosome-scale genomes of the mangrove tree Sonneratia alba and the related inland plant Lagerstroemia speciosa. Their common ancestor has experienced a whole-genome triplication (WGT) approximately 64 million years ago coinciding with a period of dramatic global climate change. Sonneratia, adapting mangrove habitats, experienced extensive chromosome rearrangements post-WGT. We observe the WGT retentions display sequence and expression divergence, suggesting potential neo- and sub-functionalization. Strong selection acting on three-copy retentions indicates adaptive value in response to new environments. To elucidate the role of ploidy changes in genome evolution, we improve a model of the polyploidization-rediploidization process based on genomic evidence, contributing to the understanding of adaptive evolution during climate change.

Crop-GPA: an integrated platform of crop gene-phenotype associations

With the increasing availability of large-scale biology data in crop plants, there is an urgent demand for a versatile platform that fully mines and utilizes the data for modern molecular breeding. We present Crop-GPA ( https://crop-gpa.aielab.net ), a comprehensive and functional open-source platform for crop gene-phenotype association data. The current Crop-GPA provides well-curated information on genes, phenotypes, and their associations (GPAs) to researchers through an intuitive interface, dynamic graphical visualizations, and efficient online tools. Two computational tools, GPA-BERT and GPA-GCN, are specifically developed and integrated into Crop-GPA, facilitating the automatic extraction of gene-phenotype associations from bio-crop literature and predicting unknown relations based on known associations. Through usage examples, we demonstrate how our platform enables the exploration of complex correlations between genes and phenotypes in crop plants. In summary, Crop-GPA serves as a valuable multi-functional resource, empowering the crop research community to gain deeper insights into the biological mechanisms of interest.

Characterization of soybean chitinase genes induced by rhizobacteria involved in the defense against Fusarium oxysporum

Rhizobacteria are capable of inducing defense responses via the expression of pathogenesis-related proteins (PR-proteins) such as chitinases, and many studies have validated the functions of plant chitinases in defense responses. Soybean (Glycine max) is an economically important crop worldwide, but the functional validation of soybean chitinase in defense responses remains limited. In this study, genome-wide characterization of soybean chitinases was conducted, and the defense contribution of three chitinases (GmChi01, GmChi02, or GmChi16) was validated in Arabidopsis transgenic lines against the soil-borne pathogen Fusarium oxysporum. Compared to the Arabidopsis Col-0 and empty vector controls, the transgenic lines with GmChi02 or GmChi16 exhibited fewer chlorosis symptoms and wilting. While GmChi02 and GmChi16 enhanced defense to F. oxysporum, GmChi02 was the only one significantly induced by Burkholderia ambifaria. The observation indicated that plant chitinases may be induced by different rhizobacteria for defense responses. The survey of 37 soybean chitinase gene expressions in response to six rhizobacteria observed diverse inducibility, where only 10 genes were significantly upregulated by at least one rhizobacterium and 9 genes did not respond to any of the rhizobacteria. Motif analysis on soybean promoters further identified not only consensus but also rhizobacterium-specific transcription factor-binding sites for the inducible chitinase genes. Collectively, these results confirmed the involvement of GmChi02 and GmChi16 in defense enhancement and highlighted the diverse inducibility of 37 soybean chitinases encountering F. oxysporum and six rhizobacteria.

Generation and characterisation of an Arabidopsis thaliana f3h/fls1/ans triple mutant that accumulates eriodictyol derivatives

BACKGROUND

Flavonoids are plant specialised metabolites, which derive from phenylalanine and acetate metabolism. They possess a variety of beneficial characteristics for plants and humans. Several modification steps in the synthesis of tricyclic flavonoids cause for the amazing diversity of flavonoids in plants. The 2-oxoglutarate-dependent dioxygenases (2-ODDs) flavanone 3-hydroxylase (F3H, synonym FHT), flavonol synthase (FLS) and anthocyanidin synthase (ANS, synonym leucoanthocyanidin dioxygenase (LDOX)), catalyse oxidative modifications to the central C ring. They are highly similar and have been shown to catalyse, at least in part, each other's reactions. FLS and ANS have been identified as bifunctional enzymes in many species, including Arabidopsis thaliana, stressing the capability of plants to bypass missing or mutated reaction steps on the way to flavonoid production. However, little is known about such bypass reactions and the flavonoid composition of plants lacking all three central flavonoid 2-ODDs.

RESULTS

To address this issue, we generated a f3h/fls1/ans mutant, as well as the corresponding double mutants and investigated the flavonoid composition of this mutant collection. The f3h/fls1/ans mutant was further characterised at the genomic level by analysis of a nanopore DNA sequencing generated genome sequence assembly and at the transcriptomic level by RNA-Seq analysis. The mutant collection established, including the novel double mutants f3h/fls1 and f3h/ans, was used to validate and analyse the multifunctionalities of F3H, FLS1, and ANS in planta. Metabolite analyses revealed the accumulation of eriodictyol and additional glycosylated derivatives in mutants carrying the f3h mutant allele, resulting from the conversion of naringenin to eriodictyol by flavonoid 3'-hydroxylase (F3'H) activity.

CONCLUSIONS

We describe the in planta multifunctionality of the three central flavonoid 2-ODDs from A. thaliana and identify a bypass in the f3h/fls1/ans triple mutant that leads to the formation of eriodictyol derivatives. As (homo-)eriodictyols are known as bitter taste maskers, the annotated eriodictyol (derivatives) and in particular the observations made on their in planta production, could provide valuable insights for the creation of novel food supplements.

Light controls mesophyll-specific post-transcriptional splicing of photoregulatory genes by AtPRMT5

Light plays a central role in plant growth and development, providing an energy source and governing various aspects of plant morphology. Previous study showed that many polyadenylated full-length RNA molecules within the nucleus contain unspliced introns (post-transcriptionally spliced introns, PTS introns), which may play a role in rapidly responding to changes in environmental signals. However, the mechanism underlying post-transcriptional regulation during initial light exposure of young, etiolated seedlings remains elusive. In this study, we used FLEP-seq2, a Nanopore-based sequencing technique, to analyze nuclear RNAs in Arabidopsis (Arabidopsis thaliana) seedlings under different light conditions and found numerous light-responsive PTS introns. We also used single-nucleus RNA sequencing (snRNA-seq) to profile transcripts in single nucleus and investigate the distribution of light-responsive PTS introns across distinct cell types. We established that light-induced PTS introns are predominant in mesophyll cells during seedling de-etiolation following exposure of etiolated seedlings to light. We further demonstrated the involvement of the splicing-related factor A. thaliana PROTEIN ARGININE METHYLTRANSFERASE 5 (AtPRMT5), working in concert with the E3 ubiquitin ligase CONSTITUTIVE PHOTOMORPHOGENIC 1 (COP1), a critical repressor of light signaling pathways. We showed that these two proteins orchestrate light-induced PTS events in mesophyll cells and facilitate chloroplast development, photosynthesis, and morphogenesis in response to ever-changing light conditions. These findings provide crucial insights into the intricate mechanisms underlying plant acclimation to light at the cell-type level.