Genome-Wide Identification of Soybean ABC Transporters Relate to Aluminum Toxicity

ATP-binding cassette (ABC) transporter proteins are a gene super-family in plants and play vital roles in growth, development, and response to abiotic and biotic stresses. The ABC transporters have been identified in crop plants such as rice and buckwheat, but little is known about them in soybean. Soybean is an important oil crop and is one of the five major crops in the world. In this study, 255 ABC genes that putatively encode ABC transporters were identified from soybean through bioinformatics and then categorized into eight subfamilies, including 7 ABCAs, 52 ABCBs, 48 ABCCs, 5 ABCDs, 1 ABCEs, 10 ABCFs, 111 ABCGs, and 21 ABCIs. Their phylogenetic relationships, gene structure, and gene expression profiles were characterized. Segmental duplication was the main reason for the expansion of the GmABC genes. Ka/Ks analysis suggested that intense purifying selection was accompanied by the evolution of GmABC genes. The genome-wide collinearity of soybean with other species showed that GmABCs were relatively conserved and that collinear ABCs between species may have originated from the same ancestor. Gene expression analysis of GmABCs revealed the distinct expression pattern in different tissues and diverse developmental stages. The candidate genes GmABCB23, GmABCB25, GmABCB48, GmABCB52, GmABCI1, GmABCI5, and GmABCI13 were responsive to Al toxicity. This work on the GmABC gene family provides useful information for future studies on ABC transporters in soybean and potential targets for the cultivation of new germplasm resources of aluminum-tolerant soybean.     

Genome-wide Identification and Characterization of FCS-Like Zinc Finger (FLZ) Family Genes in Maize (Zea mays) and Functional Analysis of ZmFLZ25 in Plant Abscisic Acid Response

FCS-like zinc finger family proteins (FLZs), a class of plant-specific scaffold of SnRK1 complex, are involved in the regulation of various aspects of plant growth and stress responses. Most information of FLZ family genes was obtained from the studies in Arabidopsis thaliana, whereas little is known about the potential functions of FLZs in crop plants. In this study, 37 maize FLZ (ZmFLZ) genes were identified to be asymmetrically distributed on 10 chromosomes and can be divided into three subfamilies. Protein interaction and subcellular localization assays demonstrated that eight typical ZmFLZs interacted and partially co-localized with ZmKIN10, the catalytic α-subunit of the SnRK1 complex in maize leaf mesophyll cells. Expression profile analysis revealed that several ZmFLZs were differentially expressed across various tissues and actively responded to diverse abiotic stresses. In addition, ectopic overexpression of ZmFLZ25 in Arabidopsis conferred hypersensitivity to exogenous abscisic acid (ABA) and triggered higher expression of ABA-induced genes, pointing to the positive regulatory role of ZmFLZ25 in plant ABA signaling, a scenario further evidenced by the interactions between ZmFLZ25 and ABA receptors. In summary, these data provide the most comprehensive information on FLZ family genes in maize, and shed light on the biological function of ZmFLZ25 in plant ABA signaling.     

Evolution of RGF/GLV/CLEL Peptide Hormones and Their Roles in Land Plant Growth and Regulation

Rooting is a key innovation during plant terrestrialization. RGFs/GLVs/CLELs are a family of secreted peptides, playing key roles in root stem cell niche maintenance and pattern formation. The origin of this peptide family is not well characterized. RGFs and their receptor genes, RGIs, were investigated comprehensively using phylogenetic and genetic analyses. We identified 203 RGF genes from 24 plant species, representing a variety of land plant lineages. We found that the RGF genes originate from land plants and expand via multiple duplication events. The lineage-specific RGF duplicates are retained due to their regulatory divergence, while a majority of RGFs experienced strong purifying selection in most land plants. Functional analysis indicated that RGFs and their receptor genes, RGIs, isolated from liverwort, tomato, and maize possess similar biological functions with their counterparts from Arabidopsis in root development. RGFs and RGIs are likely coevolved in land plants. Our studies shed light on the origin and functional conservation of this important peptide family in plant root development.     

Distinct Evolutionary Profiles and Functions of microRNA156 and microRNA529 in Land Plants

MicroRNA156 (miR156) and miR529 have high sequence similarity and recognize overlapping sites in the same target genes, SQUAMOSA promoter binding protein-like (SPL or SBP box) genes, making it difficult to accurately distinguish their roles in regulatory networks that affect numerous biological functions. Here, we collected data about miR156 and miR529 family members from representative land plants and performed sequence comparisons, phylogenetic analysis, small RNA sequencing, and parallel analysis of RNA ends (PARE) analysis to dissect their evolutionary and functional differences. Although miR156 and miR529 are highly similar, there are differences in their mismatch-sensitive regions, which are essential for target recognition. In land plants, miR156 precursors are conserved mainly within the hairpin region, whereas miR529 precursors are conserved outside the hairpin region, including both the 5’ and 3’ arms. Phylogenetic analysis showed that MIR156 and MIR529 evolved independently, through divergent evolutionary patterns. The two genes also exhibit different expression patterns, with MIR529 preferentially expressed in reproductive tissues and MIR156 in other tissues. PARE analysis revealed that miR156 and miR529 possess specific targets in addition to common targets in maize, pointing to functional differences between them. Based on our findings, we developed a method for the rapid identification of miR529 and miR156 family members and uncovered the evolutionary divergence of these families, providing insights into their different regulatory roles in plant growth and development.     

Rice SUT and SWEET Transporters

Sugar transporters play important or even indispensable roles in sugar translocation among adjacent cells in the plant. They are mainly composed of sucrose–proton symporter SUT family members and SWEET family members. In rice, 5 and 21 members are identified in these transporter families, and some of their physiological functions have been characterized on the basis of gene knockout or knockdown strategies. Existing evidence shows that most SUT members play indispensable roles, while many SWEET members are seemingly not so critical in plant growth and development regarding whether their mutants display an aberrant phenotype or not. Generally, the expressions of SUT and SWEET genes focus on the leaf, stem, and grain that represent the source, transport, and sink organs where carbohydrate production, allocation, and storage take place. Rice SUT and SWEET also play roles in both biotic and abiotic stress responses in addition to plant growth and development. At present, these sugar transporter gene regulation mechanisms are largely unclear. In this review, we compare the expressional profiles of these sugar transporter genes on the basis of chip data and elaborate their research advances. Some suggestions concerning future investigation are also proposed.     

Genome-Wide Analysis of the Cyclin Gene Family and Their Expression Profile in Medicago truncatula

Cyclins, together with highly conserved cyclin-dependent kinases (CDKs), play an important role in the process of cell cycle in plants, but less is known about the functions of cyclins in legume plants, especially Medicago truncatula. Our genome-wide analysis identified 58, 103, and 51 cyclin members in the M. truncatula, Glycine max, and Phaseolus vulgaris genomes. Phylogenetic analysis suggested that these cyclins could be classified into 10 types, and the CycB-like types (CycBL1-BL8) were the specific subgroups in M. truncatula, which was one reason for the expansion of the B-type in M. truncatula. All putative cyclin genes were mapped onto their own chromosomes of each genome, and 9 segmental duplication gene pairs involving 20 genes were identified in M. truncatula cyclins. Determined by quantitative real-time PCR, the expression profiling suggested that 57 cyclins in M. truncatula were differentially expressed in 9 different tissues, while a few genes were expressed in some specific tissues. Using the publicly available RNAseq data, the expression of Mtcyclins in the wild-type strain A17 and three nodule mutants during rhizobial infection showed that 23 cyclins were highly upregulated in the nodulation (Nod) factor-hypersensitive mutant sickle (skl) mutant after 12 h of rhizobium inoculation. Among these cyclins, six cyclin genes were also specifically expressed in roots and nodules, which might play specific roles in the various phases of Nod factor-mediated cell cycle activation and nodule development. Our results provide information about the cyclin gene family in legume plants, serving as a guide for further functional research on plant cyclins.     

Molecular Cloning and Functional Analysis of GmLACS2-3 Reveals Its Involvement in Cutin and Suberin Biosynthesis along with Abiotic Stress Tolerance

Cutin and wax are the main precursors of the cuticle that covers the aerial parts of plants and provide protection against biotic and abiotic stresses. Long-chain acyl-CoA synthetases (LACSs) play diversified roles in the synthesis of cutin, wax, and triacylglycerol (TAG). Most of the information concerned with LACS functions is obtained from model plants, whereas the roles of LACS genes in Glycine max are less known. Here, we have identified 19 LACS genes in Glycine max, an important crop plant, and further focused our attention on 4 LACS2 genes (named as GmLACS2-1, 2, 3, 4, respectively). These GmLACS2 genes display different expression patterns in various organs and also show different responses to abiotic stresses, implying that these genes might play diversified functions during plant growth and against stresses. To further identify the role of GmLACS2-3, greatly induced by abiotic stresses, we transformed a construct containing its full length of coding sequence into Arabidopsis. The expression of GmLACS2-3 in an Arabidopsis atlacs2 mutant greatly suppressed its phenotype, suggesting it plays conserved roles with that of AtLACS2. The overexpression of GmLACS2-3 in wild-type plants significantly increased the amounts of cutin and suberin but had little effect on wax amounts, indicating the specific role of GmLACS2-3 in the synthesis of cutin and suberin. In addition, these GmLACS2-3 overexpressing plants showed enhanced drought tolerance. Taken together, our study deepens our understanding of the functions of LACS genes in different plants and also provides a clue for cultivating crops with strong drought resistance.     

Characterization of the WAK Gene Family Reveals Genes for FHB Resistance in Bread Wheat (Triticum aestivum L.)

Wall-associated kinases (WAKs) are important receptor-like proteins that play major roles in plant defense against pathogens. Fusarium head blight (FHB), one of the most widespread and devastating crop diseases, reduces wheat yield and leads to quality deterioration. Although WAK gene families have been studied in many plants, systematic research on bread wheat (Triticum aestivum) and its role in FHB resistance, in particular, is lacking. In this study, we identified and characterized 320 genes of the TaWAK family in wheat distributed across all chromosomes except 4B and divided them into three phylogenetic groups. Duplication and synteny analyses provided valuable information on the evolutionary characteristics of the TaWAK genes. The gene expression pattern analysis suggested that TaWAK genes play diverse roles in plant biological processes and that at least 30 genes may be involved in the response to Fusarium infection in wheat spikes, with most of the genes contributing to pectin- and chitin-induced defense pathways. Furthermore, 45 TaWAK genes were identified within 17 hcmQTLs that are related to wheat FHB resistance. Our findings provide potential candidate genes for improving FHB resistance and insights into the future functional analysis of TaWAK genes in wheat.     

Molecular Evolution and Functional Divergence of Stress-Responsive Cu/Zn Superoxide Dismutases in Plants

Superoxide dismutases (SODs), a family of antioxidant enzymes, are the first line of defense against oxidative damage and are ubiquitous in every cell of all plant types. The Cu/Zn SOD, one of three types of SODs present in plant species, is involved in many of the biological functions of plants in response to abiotic and biotic stresses. Here, we carried out a comprehensive analysis of the Cu/Zn SOD gene family in different plant species, ranging from lower plants to higher plants, and further investigated their organization, sequence features, and expression patterns in response to biotic and abiotic stresses. Our results show that plant Cu/Zn SODs can be divided into two subfamilies (group I and group II). Group II appeared to be conserved only as single- or low-copy genes in all lineages, whereas group I genes underwent at least two duplication events, resulting in multiple gene copies and forming three different subgroups (group Ia, group Ib, and group Ic). We also found that, among these genes, two important events—the loss of introns and the loss of and variation in signal peptides—occurred over the long course of their evolution, indicating that they were involved in shifts in subcellular localization from the chloroplast to cytosol or peroxisome and underwent functional divergence. In addition, expression patterns of Cu/Zn SOD genes from Arabidopsis thaliana and Solanum lycopersicum were tested in different tissues/organs and developmental stages and under different abiotic stresses. The results indicate that the Cu/Zn SOD gene family possesses potential functional divergence and may play vital roles in ROS scavenging in response to various stresses in plants. This study will help establish a foundation for further understanding these genes’ function during stress responses.     

Genome-Wide Comprehensive Analysis of the GASA Gene Family in Populus

Gibberellic acid-stimulated Arabidopsis (GASA) proteins, as cysteine-rich peptides (CRPs), play roles in development and reproduction and biotic and abiotic stresses. Although the GASA gene family has been identified in plants, the knowledge about GASAs in Populus euphratica, the woody model plant for studying abiotic stress, remains limited. Here, we referenced the well-sequenced Populus trichocarpa genome, and identified the GASAs in the whole genome of P. euphratica and P. trichocarpa. 21 candidate genes in P. trichocarpa and 19 candidate genes in P. euphratica were identified and categorized into three subfamilies by phylogenetic analysis. Most GASAs with signal peptides were located extracellularly. The GASA genes in Populus have experienced multiple gene duplication events, especially in the subfamily A. The evolution of the subfamily A, with the largest number of members, can be attributed to whole-genome duplication (WGD) and tandem duplication (TD). Collinearity analysis showed that WGD genes played a leading role in the evolution of GASA genes subfamily B. The expression patterns of P. trichocarpa and P. euphratica were investigated using the PlantGenIE database and the real-time quantitative PCR (qRT-PCR), respectively. GASA genes in P. trichocarpa and P. euphratica were mainly expressed in young tissues and organs, and almost rarely expressed in mature leaves. GASA genes in P. euphratica leaves were also widely involved in hormone responses and drought stress responses. GUS activity assay showed that PeuGASA15 was widely present in various organs of the plant, especially in vascular bundles, and was induced by auxin and inhibited by mannitol dramatically. In summary, this present study provides a theoretical foundation for further research on the function of GASA genes in P. euphratica.     

Genome-Wide Identification,Characterization and Expression Analysis of the Chalcone Synthase Family in Maize

Members of the chalcone synthase (CHS) family participate in the synthesis of a series of secondary metabolites in plants, fungi and bacteria. The metabolites play important roles in protecting land plants against various environmental stresses during the evolutionary process. Our research was conducted on comprehensive investigation of CHS genes in maize (Zea mays L.), including their phylogenetic relationships, gene structures, chromosomal locations and expression analysis. Fourteen CHS genes (ZmCHS01–14) were identified in the genome of maize, representing one of the largest numbers of CHS family members identified in one organism to date. The gene family was classified into four major classes (classes I–IV) based on their phylogenetic relationships. Most of them contained two exons and one intron. The 14 genes were unevenly located on six chromosomes. Two segmental duplication events were identified, which might contribute to the expansion of the maize CHS gene family to some extent. In addition, quantitative real-time PCR and microarray data analyses suggested that ZmCHS genes exhibited various expression patterns, indicating functional diversification of the ZmCHS genes. Our results will contribute to future studies of the complexity of the CHS gene family in maize and provide valuable information for the systematic analysis of the functions of the CHS gene family.     

HDACs Gene Family Analysis of Eight Rosaceae Genomes Reveals the Genomic Marker of Cold Stress in Prunus mume

Histone deacetylases (HDACs) play important roles in plant growth, development, and stress response. However, the pattern of how they are expressed in response to cold stress in the ornamental woody plant Prunus mume is poorly understood. Here, we identify 121 RoHDACs from eight Rosaceae plants of which 13 PmHDACs genes are from P. mume. A phylogenetic analysis suggests that the RoHDACs family is classified into three subfamilies, HDA1/RPD3, HD2, and SIR2. We identify 11 segmental duplication gene pairs of RoHDACs and find, via a sequence alignment, that the HDACs gene family, especially the plant-specific HD2 family, has experienced gene expansion and contraction at a recent genome evolution history. Each of the three HDACs subfamilies has its own conserved domains. The expression of PmHDACs in mei is found to be tissue-specific or tissue-wide. RNA-seq data and qRT-PCR experiments in cold treatments suggest that almost all PmHDACs genes—especially PmHDA1/6/14, PmHDT1, and PmSRT1/2—significantly respond to cold stress. Our analysis provides a fundamental insight into the phylogenetic relationship of the HDACs family in Rosaceae plants. Expression profiles of PmHDACs in response to cold stress could provide an important clue to improve the cold hardiness of mei.