Polyploidization events shaped the transcription factor repertoires in legumes (Fabaceae)

Haplotype-resolved genome assembly provides insights into evolutionary history of the Actinidia arguta tetraploid

Actinidia arguta, known as hardy kiwifruit, is a widely cultivated species with distinct botanical characteristics such as small and smooth-fruited, rich in beneficial nutrients, rapid softening and tolerant to extremely low temperatures. It contains the most diverse ploidy types, including diploid, tetraploid, hexaploid, octoploid, and decaploid. Here we report a haplotype-resolved tetraploid genome (A. arguta cv. 'Longcheng No.2') containing four haplotypes, each with 40,859, 41,377, 39,833 and 39,222 protein-coding genes. We described the phased genome structure, synteny, and evolutionary analyses to identify and date possible WGD events. Ks calculations for both allelic and paralogous genes pairs throughout the assembled haplotypic individuals showed its tetraploidization is estimated to have formed ~ 1.03 Mya following Ad-α event occurred ~ 18.7 Mya. Detailed annotations of NBS-LRRs or CBFs highlight the importance of genetic variations coming about after polyploidization in underpinning ability of immune responses or environmental adaptability. WGCNA analysis of postharvest quality indicators in combination with transcriptome revealed several transcription factors were involved in regulating ripening kiwi berry texture. Taking together, the assembly of an A. arguta tetraploid genome provides valuable resources in deciphering complex genome structure and facilitating functional genomics studies and genetic improvement for kiwifruit and other crops.

Genome-wide identification of potassium channels in maize showed evolutionary patterns and variable functional responses to abiotic stresses

Potassium (K) channels are essential components of plant biology, mediating not only K ion (K+) homeostasis but also regulating several physiological processes and stress tolerance. In the current investigation, we identified 27 K+ channels in maize and deciphered the evolution and divergence pattern with four monocots and five dicot species. Chromosomal localization and expansion of K+ channel genes showed uneven distribution and were independent of genome size. The dispersed duplication is the major force in expanding K+ channels in the target genomes. The mean Ka/Ks ratio of <0.5 in paralogs and orthologs indicates horizontal and vertical expansions of K+ channel genes under strong purifying selection. The one-to-one K+ channel orthologs were prominent among the closely related species, with higher synteny between maize and the rest of the monocots. Comprehensive K+ channels promoter analysis revealed various cis-regulatory elements mediating stress tolerance with the predominance of MYB and STRE binding sites. The regulatory network showed AP2-EREBP TFs, miR164 and miR399 are prominent regulatory elements of K+ channels. The qRT-PCR analysis of K+ channels and regulatory miRNAs showed significant expressions in response to drought and waterlogging stresses. The present study expanded the knowledge on K+ channels in maize and will serve as a basis for an in-depth functional analysis.

Whole-genome Duplications and the Long-term Evolution of Gene Regulatory Networks in Angiosperms

Angiosperms have a complex history of whole-genome duplications (WGDs), with varying numbers and ages of WGD events across clades. These WGDs have greatly affected the composition of plant genomes due to the biased retention of genes belonging to certain functional categories following their duplication. In particular, regulatory genes and genes encoding proteins that act in multiprotein complexes have been retained in excess following WGD. Here, we inferred protein-protein interaction (PPI) networks and gene regulatory networks (GRNs) for seven well-characterized angiosperm species and explored the impact of both WGD and small-scale duplications (SSDs) in network topology by analyzing changes in frequency of network motifs. We found that PPI networks are enriched in WGD-derived genes associated with dosage-sensitive intricate systems, and strong selection pressures constrain the divergence of WGD-derived genes at the sequence and PPI levels. WGD-derived genes in network motifs are mostly associated with dosage-sensitive processes, such as regulation of transcription and cell cycle, translation, photosynthesis, and carbon metabolism, whereas SSD-derived genes in motifs are associated with response to biotic and abiotic stress. Recent polyploids have higher motif frequencies than ancient polyploids, whereas WGD-derived network motifs tend to be disrupted on the longer term. Our findings demonstrate that both WGD and SSD have contributed to the evolution of angiosperm GRNs, but in different ways, with WGD events likely having a more significant impact on the short-term evolution of polyploids.

Polyploidy and diploidization in soybean

Polyploidy is widespread and particularly common in angiosperms. The prevalence of polyploidy in the plant suggests it as a crucial driver of diversification and speciation. The paleopolyploid soybean (Glycine max) is one of the most important crops of plant protein and oil for humans and livestock. Soybean experienced two rounds of whole genome duplication around 13 and 59 million years ago. Due to the relatively slow process of post-polyploid diploidization, most genes are present in multiple copies across the soybean genome. Growing evidence suggests that polyploidization and diploidization could cause rapid and dramatic changes in genomic structure and epigenetic modifications, including gene loss, transposon amplification, and reorganization of chromatin architecture. This review is focused on recent progresses about genetic and epigenetic changes during polyploidization and diploidization of soybean and represents the challenges and potentials for application of polyploidy in soybean breeding.

Identification of Photoperiod- and Phytohormone-Responsive DNA-Binding One Zinc Finger (Dof) Transcription Factors in Akebia trifoliata via Genome-Wide Expression Analysis

As a kind of plant-specific transcription factor (TF), DNA-Binding One Zinc Finger (Dof) is widely involved in the response to environmental change, and as an evolutionarily important perennial plant species, Akebia trifoliata is ideal for studying environmental adaptation. In this study, a total of 41 AktDofs were identified in the A. trifoliata genome. First, the characteristics, including the length, exon number, and chromosomal distribution of the AktDofs and the isoelectric point (PI), amino acid number, molecular weight (MW), and conserved motifs of their putative proteins, were reported. Second, we found that all AktDofs evolutionarily underwent strong purifying selection, and many (33, 80.5%) of them were generated by whole-genome duplication (WGD). Third, we outlined their expression profiles by the use of available transcriptomic data and RT-qPCR analysis. Finally, we identified four candidate genes (AktDof21, AktDof20, AktDof36, and AktDof17) and three other candidate genes (AktDof26, AktDof16, and AktDof12) that respond to long day (LD) and darkness, respectively, and that are closely associated with phytohormone-regulating pathways. Overall, this research is the first to identify and characterize the AktDofs family and is very helpful for further research on A. trifoliata adaptation to environmental factors, especially photoperiod changes.

Genome-wide characterization of the common bean kinome: Catalog and insights into expression patterns and genetic organization

The protein kinase (PK) superfamily is one of the largest superfamilies in plants and is the core regulator of cellular signaling. Even considering this substantial importance, the kinome of common bean (Phaseolus vulgaris) has not been profiled yet. Here, we identified and characterised the complete set of kinases of common bean, performing an in-depth investigation with phylogenetic analyses and measurements of gene distribution, structural organization, protein properties, and expression patterns over a large set of RNA-Sequencing data. Being composed of 1,203 PKs distributed across all P. vulgaris chromosomes, this set represents 3.25% of all predicted proteins for the species. These PKs could be classified into 20 groups and 119 subfamilies, with a more pronounced abundance of subfamilies belonging to the receptor-like kinase (RLK)-Pelle group. In addition to provide a vast and rich reservoir of data, our study supplied insights into the compositional similarities between PK subfamilies, their evolutionary divergences, highly variable functional profile, structural diversity, and expression patterns, modeled with coexpression networks for investigating putative interactions associated with stress response.

TVIR: a comprehensive vegetable information resource database for comparative and functional genomic studies

Vegetables are an indispensable part of the daily diet of humans. Therefore, it is vital to systematically study the genomic data of vegetables and build a platform for data sharing and analysis. In this study, a comprehensive platform for vegetables with a user-friendly Web interface-The Vegetable Information Resource (TVIR, http://tvir.bio2db.com)-was built based on the genomes of 59 vegetables. TVIR database contains numerous important functional genes, including 5215 auxin genes, 2437 anthocyanin genes, 15 002 flowering genes, 79 830 resistance genes, and 2639 glucosinolate genes of 59 vegetables. In addition, 2597 N6-methyladenosine (m6A) genes were identified, including 513 writers, 1058 erasers, and 1026 readers. A total of 2 101 501 specific clustered regularly interspaced short palindromic repeat (CRISPR) guide sequences and 17 377 miRNAs were detected and deposited in TVIR database. Information on gene synteny, duplication, and orthologs is also provided for 59 vegetable species. TVIR database contains 2 346 850 gene annotations by the Swiss-Prot, TrEMBL, Gene Ontology (GO), Pfam, and Non-redundant (Nr) databases. Synteny, Primer Design, Blast, and JBrowse tools are provided to facilitate users in conducting comparative genomic analyses. This is the first large-scale collection of vegetable genomic data and bioinformatic analysis. All genome and gene sequences, annotations, and bioinformatic results can be easily downloaded from TVIR. Furthermore, transcriptome data of 98 vegetables have been collected and collated, and can be searched by species, tissues, or different growth stages. TVIR is expected to become a key hub for vegetable research globally. The database will be updated with newly assembled vegetable genomes and comparative genomic studies in the future.

Genome-wide analysis of the homeodomain-leucine zipper family in Lotus japonicus and the overexpression of LjHDZ7 in Arabidopsis for salt tolerance

The homeodomain-leucine zipper (HD-Zip) family participates in plant growth, development, and stress responses. Here, 40 HD-Zip transcription factors of Lotus japonicus were identified and gave an overview of the phylogeny and gene structures. The expression pattern of these candidate genes was determined in different organs and their response to abiotic stresses, including cold, heat, polyethylene glycol and salinity. The expression of the LjHDZ7 was strongly induced by abiotic stress, especially salt stress. Subsequently, LjHDZ7 gene was overexpressed in Arabidopsis. The transgenic plants grew obviously better than Col-0 plants under salt stress. Furthermore, LjHDZ7 transgenic lines accumulated higher proline contents and showed lower electrolyte leakage and MDA contents than Col-0 plants under salt stress. Antioxidant activities of the LjHDZ7 overexpression lines leaf were significantly higher than those of the Col-0 plants under salt stress. The concentration of Na⁺ ion in LjHDZ7 overexpression lines was significantly lower than that of Col-0 in leaf and root parts. The concentration of K⁺ ion in LjHDZ7 overexpression lines was significantly higher than that of Col-0 in the leaf parts. Therefore, these results showed that overexpression of LjHDZ7 increased resistance to salt stress in transgenic Arabidopsis plants, and certain genes of this family can be used as valuable tools for improving abiotic stresses.

Integrating omics approaches to discover and prioritize candidate genes involved in oil biosynthesis in soybean

Soybean is a major source of edible protein and oil. Oil content is a quantitative trait that is significantly determined by genetic and environmental factors. Over the past 30 years, a large volume of soybean genetic, genomic, and transcriptomic data have been accumulated. Nevertheless, integrative analyses of such data remain scarce, in spite of their importance for crop improvement. We hypothesized that the co-occurrence of genomic regions for oil-related traits in different studies may reveal more stable regions encompassing important genetic determinants of oil content and quality in soybean. We integrated publicly available data, obtained with distinct techniques, to discover and prioritize candidate genes involved in oil biosynthesis and regulation in soybean. We detected key fatty acid biosynthesis genes (e.g., BCCP2 and ACCase, FADs, KAS family proteins) and several transcription factors, which are likely regulators of oil biosynthesis. In addition, we identified new candidates for seed oil accumulation and quality, such as Glyma.03G213300 and Glyma.19G160700, which encode a translocator protein homolog and a histone acetyltransferase, respectively. Further, oil and protein genomic hotspots are strongly associated with breeding and not with domestication, suggesting that soybean domestication prioritized other traits. The genes identified here are promising targets for breeding programs and for the development of soybean lines with increased oil content and quality.

Exploring the complexity of soybean (Glycine max) transcriptional regulation using global gene co-expression networks

MAIN CONCLUSION

We report a soybean gene co-expression network built with data from 1284 RNA-Seq experiments, which was used to identify important regulators, modules and to elucidate the fates of gene duplicates. Soybean (Glycine max (L.) Merr.) is one of the most important crops worldwide, constituting a major source of protein and edible oil. Gene co-expression networks (GCN) have been extensively used to study transcriptional regulation and evolution of genes and genomes. Here, we report a soybean GCN using 1284 publicly available RNA-Seq samples from 15 distinct tissues. We found modules that are differentially regulated in specific tissues, comprising processes such as photosynthesis, gluconeogenesis, lignin metabolism, and response to biotic stress. We identified transcription factors among intramodular hubs, which probably integrate different pathways and shape the transcriptional landscape in different conditions. The top hubs for each module tend to encode proteins with critical roles, such as succinate dehydrogenase and RNA polymerase subunits. Importantly, gene essentiality was strongly correlated with degree centrality and essential hubs were enriched in genes involved in nucleic acids metabolism and regulation of cell replication. Using a guilt-by-association approach, we predicted functions for 93 of 106 hubs without functional description in soybean. Most of the duplicated genes had different transcriptional profiles, supporting their functional divergence, although paralogs originating from whole-genome duplications (WGD) are more often preserved in the same module than those from other mechanisms. Together, our results highlight the importance of GCN analysis in unraveling key functional aspects of the soybean genome, in particular those associated with hub genes and WGD events.

Systematic analysis of 1298 RNA-Seq samples and construction of a comprehensive soybean (Glycine max) expression atlas

Soybean (Glycine max [L.] Merr.) is a major crop in animal feed and human nutrition, mainly for its rich protein and oil contents. The remarkable rise in soybean transcriptome studies over the past 5 years generated an enormous amount of RNA-seq data, encompassing various tissues, developmental conditions and genotypes. In this study, we have collected data from 1298 publicly available soybean transcriptome samples, processed the raw sequencing reads and mapped them to the soybean reference genome in a systematic fashion. We found that 94% of the annotated genes (52 737/56 044) had detectable expression in at least one sample. Unsupervised clustering revealed three major groups, comprising samples from aerial, underground and seed/seed-related parts. We found 452 genes with uniform and constant expression levels, supporting their roles as housekeeping genes. On the other hand, 1349 genes showed heavily biased expression patterns towards particular tissues. A transcript-level analysis revealed that 95% (70 963 of 74 490) of the assembled transcripts have intron chains exactly matching those from known transcripts, whereas 3256 assembled transcripts represent potentially novel splicing isoforms. The dataset compiled here constitute a new resource for the community, which can be downloaded or accessed through a user-friendly web interface at http://venanciogroup.uenf.br/resources/. This comprehensive transcriptome atlas will likely accelerate research on soybean genetics and genomics.