Integrative Analysis of Oleosin Genes Provides Insights into Lineage-Specific Family Evolution in Brassicales

Oleosins (OLEs) are a class of small but abundant structural proteins that play essential roles in the formation and stabilization of lipid droplets (LDs) in seeds of oil crops. Despite the proposal of five oleosin clades (i.e., U, SL, SH, T, and M) in angiosperms, their evolution in eudicots has not been well-established. In this study, we employed Brassicales, an economically important order of flowering plants possessing the lineage-specific T clade, as an example to address this issue. Three to 10 members were identified from 10 species representing eight plant families, which include Caricaceae, Moringaceae, Akaniaceae, Capparaceae, and Cleomaceae. Evolutionary and reciprocal best hit-based homologous analyses assigned 98 oleosin genes into six clades (i.e., U, SL, SH, M, N, and T) and nine orthogroups (i.e., U1, U2, SL, SH1, SH2, SH3, M, N, and T). The newly identified N clade represents an ancient group that has already appeared in the basal angiosperm Amborella trichopoda, which are constitutively expressed in the tree fruit crop Carica papaya, including pulp and seeds of the fruit. Moreover, similar to Clade N, the previously defined M clade is actually not Lauraceae-specific but an ancient and widely distributed group that diverged before the radiation of angiosperm. Compared with A. trichopoda, lineage-specific expansion of the family in Brassicales was largely contributed by recent whole-genome duplications (WGDs) as well as the ancient γ event shared by all core eudicots. In contrast to the flower-preferential expression of Clade T, transcript profiling revealed an apparent seed/embryo/endosperm-predominant expression pattern of most oleosin genes in Arabidopsis thaliana and C. papaya. Moreover, the structure and expression divergence of paralogous pairs was frequently observed, and a good example is the lineage-specific gain of an intron. These findings provide insights into lineage-specific family evolution in Brassicales, which facilitates further functional studies in nonmodel plants such as C. papaya.

Gene Characterization of Nocturnin Paralogues in Goldfish: Full Coding Sequences, Structure, Phylogeny and Tissue Expression

The aim of this work is the full characterization of all the nocturnin (noc) paralogues expressed in a teleost, the goldfish. An in silico analysis of the evolutive origin of noc in Osteichthyes is performed, including the splicing variants and new paralogues appearing after teleostean 3R genomic duplication and the cyprinine 4Rc. After sequencing the full-length mRNA of goldfish, we obtained two isoforms for noc-a (noc-aa and noc-ab) with two splice variants (I and II), and only one for noc-b (noc-bb) with two transcripts (II and III). Using the splicing variant II, the prediction of the secondary and tertiary structures renders a well-conserved 3D distribution of four α-helices and nine β-sheets in the three noc isoforms. A synteny analysis based on the localization of noc genes in the patrilineal or matrilineal subgenomes and a phylogenetic tree of protein sequences were accomplished to stablish a classification and a long-lasting nomenclature of noc in goldfish, and valid to be extrapolated to allotetraploid Cyprininae. Finally, both goldfish and zebrafish showed a broad tissue expression of all the noc paralogues. Moreover, the enriched expression of specific paralogues in some tissues argues in favour of neo- or subfunctionalization.

Genome-wide identification of the C2H2 zinc finger gene family and expression analysis under salt stress in sweetpotato

Introduction

The higher plant transcription factor C2H2 zinc finger protein (C2H2-ZFP) is essential for plant growth, development, and stress response. There are limited studies on C2H2-ZFP genes in sweetpotato, despite a substantial number of C2H2-ZFP genes having been systematically found in plants.

Methods

In this work, 178 C2H2-ZFP genes were found in sweetpotato, distributed randomly on 15 chromosomes, and given new names according to where they were located. These members of the zinc finger gene family are separated into six branches, as shown by the phylogenetic tree. 24 tandem repeats of IbZFP genes and 46 fragment repeats were identified, and a homology study revealed that IbZFP genes linked more regions with wild relative species of sweetpotato as well as rhizome plants like potato and cassava. And we analyzed the expression patterns of IbZFP genes during the early development of sweetpotato storage roots (SRs) and salt stress using transcriptome data, and identified 44 IbZFP genes that exhibited differences in expression levels during the early expansion of sweetpotato SRs in different varieties, and 92 IbZFP genes that exhibited differences in expression levels under salt stress in salt tolerant and salt sensitive sweetpotato varieties. Additionally, we cloned six IbZFP genes in sweetpotato and analyzed their expression patterns in different tissues, their expression patterns under abiotic stress and hormone treatment, and subcellular localization.

Results and discussion

The results showed that the IbZFP genes had tissue specificity in sweetpotato and were induced to varying degrees by drought and salt stress. ABA and GA3 treatments also affected the expression of the IbZFP genes. We selected IbZFP105, which showed significant differences in expression levels under salt stress and ABA treatment, to be heterologously expressed in Arabidopsis thaliana. We found that IbZFP105 OE lines exhibited higher tolerance to salt stress and ABA stress. This indicates that IbZFP105 can enhance the salt tolerance of plants. These results systematically identified the evolution and expression patterns of members of the C2H2-ZFP gene family in sweetpotato, providing a theoretical basis for studying the role of IbZFP genes in the development of sweetpotato SRs and in resistance to stress.

Genome-Wide Identification and Analysis of the WRKY Gene Family in Asparagus officinalis

In recent years, the related research of the WRKY gene family has been gradually promoted, which is mainly reflected in the aspects of environmental stress and hormone response. However, to make the study of the WRKY gene family more complete, we also need to focus on the whole-genome analysis and identification of the family. In previous studies, the whole WRKY gene family of Arabidopsis, legumes and other plants has been thoroughly studied. However, since the publication of Asparagus officinalis genome-wide data, there has never been an analysis of the whole WRKY gene family. To understand more broadly the function of the WRKY gene family, the whole genome and salt stress transcriptome data of asparagus were used for comprehensive analysis in this study, including WRKY gene family identification, phylogenetic tree construction, analysis of conserved mods and gene domains, extraction of cis-acting elements, intron/exon analysis, species collinearity analysis, and WRKY expression analysis under salt stress. The results showed that a total of 70 genes were selected and randomly distributed on 10 chromosomes and one undefined chromosome. According to the functional classification of Arabidopsis thaliana, the WRKY family of asparagus was divided into 11 subgroups (C1-C9, U1, U2). It is worth considering that the distribution rules of gene-conserved motifs, gene domains and introns/exons in the same subfamily are similar, which suggests that genes in the same subfamily may regulate similar physiological processes. In this study, 11 cis-acting elements of WRKY family were selected, among which auxin, gibberellin, abscisic acid, salicylic acid and other hormone-regulated induction elements were involved. In addition, environmental stress (such as drought stress and low-temperature response) also accounted for a large proportion. Interestingly, we analyzed a total of two tandem duplicate genes and 13 segmental duplication genes, suggesting that this is related to the amplification of the WRKY gene family. Transcriptome data analysis showed that WRKY family genes could regulate plant growth and development by up-regulating and down-regulating gene expression under salt stress. Volcanic maps showed that 3 and 15 AoWRKY genes were significantly up-regulated or down-regulated in NI&NI+S and AMF&AMF+S, respectively. These results provide a new way to analyze the evolution and function of the WRKY gene family, and can provide a reference for the production and research of asparagus.

Identification of GA2ox Family Genes and Expression Analysis under Gibberellin Treatment in Pineapple (Ananas comosus (L.) Merr.)

Gibberellin (GAs) plays an important regulatory role in the development and growth of pineapple (Ananas comosus (L.) Merr.). Bioinformatics was used to confirm the differential expression of GA2 gibberellin oxidase gene AcGA2oxs in the pineapple genome, which laid the foundation for exploring its role in pineapple. In this study, 42 GA2ox genes (AcGA2oxs) were identified in the pineapple genome, named from AcGA2ox1 to AcGA2ox42, and divided into four groups according to phylogenetic analysis. We also analyzed the gene structure, conserved motifs and chromosome localization of AcGA2oxs. AcGA2oxs within the same group had similar gene structure and motifs composition. Collinear analysis and cis-element analysis provided the basis for understanding the evolution and function of GA2ox genes in pineapple. In addition, we selected different tissue parts to analyze the expression profile of AcGA2oxs, and the results show that 41 genes were expressed, except for AcGA2ox18. AcGA2ox18 may not be expressed in these sites or may be pseudogenes. qRT-PCR (real-time fluorescence quantitative PCR) was used to detect the relative expression levels of the GA2ox gene family under different concentrations of GA3 treatment, and it was found that AcGA2ox gene expression was upregulated in different degrees under GA3 treatment. These results provide useful information for further study on the evolution and function of the GA2ox family in pineapple.

Corrigendum: Genome-wide identification of genes encoding cystathionine beta synthase domain-containing proteins in wheat and its relationship with anther male sterility under heat stress

[This corrects the article DOI: 10.3389/fpls.2022.1061472.].

The lineage-specific evolution of the oleosin family in Theaceae

Oleosins play essential roles in stabilization of lipid droplets (LDs) and seed oil production. However, evolution of this gene family has not been reported in Theaceae, a large plant family that contains many important tea and oil tea species. In this study, a total of 65 oleosin genes were identified in nine genome-sequenced Theaceae species. Among these genomes, the gene number of oleosin showed significant difference, with Camellia sinensis var. sinensis cv. Shuchazao and Camellia lanceoleosa displayed more oleosin numbers than other species. Phylogenetic analyses revealed that Theaceae oleosin genes were classified into three clades (U, SL, SH) respectively. Proteins within the same clade had similar gene structure and motif composition. Segmental duplication was the primary driving force for the evolution of oleosin genes in Shuchazao (SCZ), Huangdan (HD), C.lanceoleosa (Cla), and wild tea (DASZ). Synteny analysis showed that most oleosin genes displayed inter-species synteny among tea and oil tea species. Expression analysis demonstrated that oleosin genes were specifically expressed in seed and kernel of Huangdan (HD) and C.lanceoleosa. Moreover, expression divergence was observed in paralogous pairs and ∼1-2 oleosin genes in each clade have become activate. This study leads to a comprehensive understanding of evolution of oleosin family in Theaceae, and provides a rich resource to further address the functions of oleosin in tea and oil tea species.

[Molecular characterization and transcriptional analysis of VrWOX genes in mungbean [Vigna radiate (L.) Wilczek]]

WUSCHEL-related homebox (WOX) gene family is a type of plant specific transcription factor, and belongs to the homeobox (HB) transcription factor superfamily. WOX genes play an important role in plant development, such as stem cell regulation and reproductive progress, and have been identified in many plant species. However, the information of mungbean VrWOX genes is limited. In this study, we identified 42 VrWOX genes in mungbean genome using Arabidopsis AtWOX genes as BLAST queries. VrWOX genes are unevenly distributed on 11 mungbean chromosomes, and chromosome 7 contains the most VrWOX genes. VrWOX genes are classified into three subgroups, the ancient group, the intermediate group and the modern/WUSCHEL group, which contains 19, 12 and 11 VrWOX members, respectively. Intraspecific synteny analysis revealed 12 VrWOX duplicated gene pairs in mungbean. Mungbean and Arabidopsis thaliana have 15 orthologous genes, and mungbean and Phaseolus vulgaris have 22 orthologous genes, respectively. The gene structure and conserved motif are different among VrWOX genes, indicating their functional diversity. The promoter regions of VrWOX genes contain different number and type of cis-acting elements, and VrWOX genes show distinct expression levels in eight mungbean tissues. Our study investigated the bioinformation and expression profiles of VrWOX genes, and provided essential information for further functional characterization of VrWOX genes.

Identification and Functional Characterization of WRKY, PHD and MYB Three Salt Stress Responsive Gene Families in Mungbean (Vigna radiata L.)

WRKY-, PHD-, and MYB-like proteins are three important types of transcription factors in mungbeans, and play an important role in development and stress resistance. The genes' structures and characteristics were clearly reported and were shown to contain the conservative WRKYGQK heptapeptide sequence, Cys4-His-cys3 zinc binding motif, and HTH (helix) tryptophan cluster W structure, respectively. Knowledge on the response of these genes to salt stress is largely unknown. To address this issue, 83 VrWRKYs, 47 VrPHDs, and 149 VrMYBs were identified by using comparative genomics, transcriptomics, and molecular biology methods in mungbeans. An intraspecific synteny analysis revealed that the three gene families had strong co-linearity and an interspecies synteny analysis showed that mungbean and Arabidopsis were relatively close in genetic relationship. Moreover, 20, 10, and 20 genes showed significantly different expression levels after 15 days of salt treatment (p < 0.05; Log2 FC > 0.5), respectively. Additionally, in the qRT-PCR analysis, VrPHD14 had varying degrees of response to NaCl and PEG treatments after 12 h. VrWRKY49 was upregulated by ABA treatment, especially in the beginning (within 24 h). VrMYB96 was significantly upregulated in the early stages of ABA, NaCl, and PEG stress treatments (during the first 4 h). VrWRKY38 was significantly upregulated by ABA and NaCl treatments, but downregulated by PEG treatment. We also constructed a gene network centered on the seven DEGs under NaCl treatment; the results showed that VrWRKY38 was in the center of the PPI network and most of the homologous Arabidopsis genes of the interacted genes were reported to have response to biological stress. Candidate genes identified in this study provide abundant gene resources for the study of salt tolerance in mungbeans.

Characterization of Peptaibols Produced by a Marine Strain of the Fungus Trichoderma endophyticum via Mass Spectrometry, Genome Mining and Phylogeny-Based Prediction

Trichoderma is recognized as a prolific producer of nonribosomal peptides (NRPs) known as peptaibols, which have remarkable biological properties, such as antimicrobial and anticancer activities, as well as the ability to promote systemic resistance in plants against pathogens. In this study, the sequencing of 11-, 14- and 15-res peptaibols produced by a marine strain of Trichoderma isolated from the ascidian Botrylloides giganteus was performed via liquid chromatography coupled to high-resolution tandem mass spectrometry (LC-MS/MS). Identification, based on multilocus phylogeny, revealed that our isolate belongs to the species T. endophyticum, which has never been reported in marine environments. Through genome sequencing and genome mining, 53 biosynthetic gene clusters (BGCs) were identified as being related to bioactive natural products, including two NRP-synthetases: one responsible for the biosynthesis of 11- and 14-res peptaibols, and another for the biosynthesis of 15-res. Substrate prediction, based on phylogeny of the adenylation domains in combination with molecular networking, permitted extensive annotation of the mass spectra related to two new series of 15-res peptaibols, which are referred to herein as "endophytins". The analyses of synteny revealed that the origin of the 15-module peptaibol synthetase is related to 18, 19 and 20-module peptaibol synthetases, and suggests that the loss of modules may be a mechanism used by Trichoderma species for peptaibol diversification. This study demonstrates the importance of combining genome mining techniques, mass spectrometry analysis and molecular networks for the discovery of new natural products.

Genome-Wide Identification of GYF-Domain Encoding Genes in Three Brassica Species and Their Expression Responding to Sclerotinia sclerotiorum in Brassica napus

GYF (glycine-tyrosine-phenylalanine)-domain-containing proteins, which were reported to participate in many aspects of biological processes in yeast and animals, are highly conserved adaptor proteins existing in almost all eukaryotes. Our previous study revealed that GYF protein MUSE11/EXA1 is involved in nucleotide-binding leucine-rich repeat (NLR) receptor-mediated defense in Arabidopsis thaliana. However, the GYF-domain encoding homologous genes are still not clear in other plants. Here, we performed genome-wide identification of GYF-domain encoding genes (GYFs) from Brassica napus and its parental species, Brassica rapa and Brassica oleracea. As a result, 26 GYFs of B. napus (BnaGYFs), 11 GYFs of B. rapa (BraGYFs), and 14 GYFs of B. oleracea (BolGYFs) together with 10 A. thaliana (AtGYFs) were identified, respectively. We, then, conducted gene structure, motif, cis-acting elements, duplication, chromosome localization, and phylogenetic analysis of these genes. Gene structure analysis indicated the diversity of the exon numbers of these genes. We found that the defense and stress responsiveness element existed in 23 genes and also identified 10 motifs in these GYF proteins. Chromosome localization exhibited a similar distribution of BnaGYFs with BraGYFs or BolGYFs in their respective genomes. The phylogenetic and gene collinearity analysis showed the evolutionary conservation of GYFs among B. napus and its parental species as well as Arabidopsis. These 61 identified GYF domain proteins can be classified into seven groups according to their sequence similarity. Expression of BnaGYFs induced by Sclerotinia sclerotiorum provided five highly upregulated genes and five highly downregulated genes, which might be candidates for further research of plant-fungal interaction in B. napus.

Genome-wide identification of cystathionine beta synthase genes in wheat and its relationship with anther male sterility under heat stress

Cystathionine beta synthase (CBS) domains containing proteins (CDCPs) plays an important role in plant development through regulation of the thioredoxin system, as well as its ability to respond to biotic and abiotic stress conditions. Despite this, no systematic study has examined the wheat CBS gene family and its relation to high temperature-induced male sterility. In this study, 66 CBS family members were identified in the wheat genome, and their gene or protein sequences were used for subsequent analysis. The TaCBS gene family was found to be unevenly distributed on 21 chromosomes, and they were classified into four subgroups according to their gene structure and phylogeny. The results of collinearity analysis showed that there were 25 shared orthologous genes between wheat, rice and Brachypodium distachyon, and one shared orthologous gene between wheat, millet and barley. The cis-regulatory elements of the TaCBS were related to JA, IAA, MYB, etc. GO and KEGG pathway analysis identified these TaCBS genes to be associated with pollination, reproduction, and signaling and cellular processes, respectively. A heatmap of wheat plants based on transcriptome data showed that TaCBS genes were expressed to a higher extent in spikelets relative to other tissues. In addition, 29 putative tae-miRNAs were identified, targeting 41 TaCBS genes. Moreover, qRT-PCR validation of six TaCBS genes indicated their critical role in anther development, as five of them were expressed at lower levels in heat-stressed male sterile anthers than in Normal anthers. Together with anther phenotypes, paraffin sections, starch potassium iodide staining, and qRT-PCR data, we hypothesized that the TaCBS gene has a very important connection with the heat-stressed sterility process in wheat, and these data provide a basis for further insight into their relationship.

Genome-wide identification and comprehensive analyses of NAC transcription factor gene family and expression analysis under Fusarium kyushuense and drought stress conditions in Passiflora edulis

The NAC gene family is one of the largest plant transcription factors (TFs) families and plays important roles in plant growth, development, metabolism, and biotic and abiotic stresses. However, NAC gene family has not been reported in passion fruit (Passiflora edulis). In this study, a total of 105 NAC genes were identified in the passion fruit genome and were unevenly distributed across all nine-passion fruit chromomere, with a maximum of 48 PeNAC genes on chromosome one. The physicochemical features of all 105 PeNAC genes varied including 120 to 3,052 amino acids, 3 to 8 conserved motifs, and 1 to 3 introns. The PeNAC genes were named (PeNAC001-PeNAC105) according to their chromosomal locations and phylogenetically grouped into 15 clades (NAC-a to NAC-o). Most PeNAC proteins were predicted to be localized in the nucleus. The cis-element analysis indicated the possible roles of PeNAC genes in plant growth, development, light, hormones, and stress responsiveness. Moreover, the PeNAC gene duplications including tandem (11 gene pairs) and segmental (12 gene pairs) were identified and subjected to purifying selection. All PeNAC proteins exhibited similar 3D structures, and a protein-protein interaction network analysis with known Arabidopsis proteins was predicted. Furthermore, 17 putative ped-miRNAs were identified to target 25 PeNAC genes. Potential TFs including ERF, BBR-BPC, Dof, and bZIP were identified in promoter region of all 105 PeNAC genes and visualized in a TF regulatory network. GO and KEGG annotation analysis exposed that PeNAC genes were related to different biological, molecular, and cellular terms. The qRT-PCR expression analysis discovered that most of the PeNAC genes including PeNAC001, PeNAC003, PeNAC008, PeNAC028, PeNAC033, PeNAC058, PeNAC063, and PeNAC077 were significantly upregulated under Fusarium kyushuense and drought stress conditions compared to controls. In conclusion, these findings lay the foundation for further functional studies of PeNAC genes to facilitate the genetic improvement of plants to stress resistance.

Genome-Wide Identification and Structural Characterization of Growth-Regulating Factors (GRFs) in Actinida eriantha and Actinidia chinensis

Growth-regulating factors (GRFs) encode plant-specific transcription factors that play a vital role in regulation of plant growth, development, and stress response. Although GRFs have been identified in various plants, there is no reported work available in Actinidia (commonly known as kiwifruit) so far. In the present study, we identified 22 GRF genes on A. chinensis (hereafter A. chinensis is referred to as Ac, and GRF genes in A. chinensis are referred to as AcGRF) distributed on 17 chromosomes and one contig, and 26 GRF genes in A. eriantha (hereafter A. eriantha is referred to as Ae, and GRF genes in A. eriantha are referred to as AeGRF) distributed on 21 chromosomes. Phylogenetic analysis showed that kiwifruit GRF proteins were clustered into five distinct groups. Additionally, kiwifruit GRFs showed motif composition and gene structure similarities within the same group. Synteny analysis showed that whole-genome duplication played a key role in the expansion of the GRF family in kiwifruit. The higher expression levels of kiwifruit GRFs in young tissues and under stress conditions indicated their regulatory role in kiwifruit growth and development. We observed two genes in Ae (AeGRF6.1, AeGRF 6.2) and two genes in Ac (AcGRF 6.1, AeGRF 6.2) significantly upregulated in different RNA-seq datasets. The presence of conserved protein structures and cis-regulatory elements caused functional divergence in duplicated gene pairs. The subcellular localization indicated the presence of kiwifruit GRFs in the nucleus of the plant cell. Protein-protein interaction analysis predicted AtGIF protein orthologs for AcGRFs and AeGRFs. Taken together, we systematically analyzed the characterization of kiwifruit GRF family members for their potential role in kiwifruit development and Pseudomonas syringae pv. actinidiae (Psa.) invasion response. Further functional studies of kiwifruit GRFs in plant growth, development, and stress response will provide valuable insights for kiwifruit breeders.

Genomic Architecture and Evolution of the Cellulose synthase Gene Superfamily as Revealed by Phylogenomic Analysis

The Cellulose synthase superfamily synthesizes cellulose and different hemicellulosic polysaccharides in plant cell walls. While much has been discovered about the evolution and function of these genes, their genomic architecture and relationship with gene (sub-)functionalization and evolution remains unclear. By using 242 genomes covering plant evolution from green algae to eudicots, we performed a large-scale analysis of synteny, phylogenetic, and functional data of the CesA superfamily. Results revealed considerable gene copy number variation across species and gene families, and also two patterns - singletons vs. tandem arrays - in chromosomic gene arrangement. Synteny analysis revealed exceptional conservation of gene architecture across species, but also lineage-specific patterns across gene (sub-)families. Synteny patterns correlated with gene sub-functionalization into primary and secondary CesAs and distinct CslD functional isoforms. Furthermore, a genomic context shift of a group of cotton secondary CesAs was associated with peculiar properties of cotton fiber synthesis. Finally, phylogenetics suggested that primary CesA sequences appeared before the secondary CesAs, while phylogenomic analyses unveiled the genomic trace of the CslD duplication that initiated the CslF family. Our results describe in detail the genomic architecture of the CesA superfamily in plants, highlighting its crucial relevance for gene diversification and sub-functionalization, and for understanding their evolution.

Genome-Wide Identification of DNA Binding with One Finger (Dof) Gene Family in Tartary Buckwheat (Fagopyrum tataricum) and Analysis of Its Expression Pattern after Exogenous Hormone Stimulation

Simple Summary

A number of studies have demonstrated that DNA binding with one finger (Dof) proteins are involved in multiple biological processes. In the present study, Dof genes or proteins in Tartary buckwheat (FtDofs) were systematically analysed, including their physical properties, phylogenetic relationships, structure, motif composition, cis-acting elements present in promoter regions, chromosomal distribution, gene duplication events, syntenic relationships, expression patterns in different tissues and different fruit developmental stages and responses to exogenous hormone stimulation. The results indicated that the expansion of FtDofs was mainly due to segmental duplication. The tissue-specific expression patterns of FtDofs and their positive responses to exogenous hormone stimulation suggest that they play important roles in the growth and development of Tartary buckwheat as well as in the adaptation to environmental changes. Collectively, this study lays a foundation for further exploration of the function of FtDof genes in Tartary buckwheat.

Abstract

DNA binding with one finger (Dof) proteins have been proven to be involved in multiple biological processes. However, genome-wide identification of the Dof gene family has not been reported for Tartary buckwheat (Fagopyrum tataricum). In this study, 35 FtDof proteins were identified, and they could be divided into nine phylogenetic subgroups. Proteins within the same subgroup had similar gene structure and motif composition. Moreover, abundant cis-acting elements were present in the promoter regions of FtDof genes. Segmental duplication was the primary driving force for the evolution of the FtDof gene family. Synteny analysis indicated that Tartary buckwheat was closer to dicotyledons, and more orthologous Dof genes existed among them. The expression pattern of FtDofs in different tissues and at different fruit developmental stages varied. Different tissues contained several genes that were specifically expressed. FtDof expression was mainly upregulated under methyl jasmonate treatment and downregulated under other hormone treatments. Taken together, FtDofs may play important roles in the growth and development of Tartary buckwheat and in response to abiotic and biotic stresses. Therefore, the genome-wide identification and expression pattern analysis of the Tartary buckwheat Dof gene family lays a foundation for further exploration of the functional characteristics of FtDofs in the future.

Genome-Wide Identification of JRL Genes in Moso Bamboo and Their Expression Profiles in Response to Multiple Hormones and Abiotic Stresses

Jacalin-related lectins (JRLs) are a new subfamily of plant lectins that has recently been recognized and plays an important role in plant growth, development, and abiotic stress response. Although moso bamboo (Phyllostachys edulis) is an economically and industrially important bamboo worldwide, there has been no systematic identification of JRLs in this species. Here, we identified 25 JRL genes in moso bamboo, and these genes are unequally distributed among 10 genome scaffolds. Phylogenetic analysis showed that the moso bamboo JRLs were clustered into four JRL subgroups: I, II, V, and VII. Numerous stress-responsive and hormone-regulated cis-elements were detected in the upstream promoter regions of the JRLs. Genome collinearity analyses showed that the JRL genes of moso bamboo are more closely related to those of Brachypodium distachyon than to those of Oryza sativa and Zea mays. Sixty-four percent of the PeJRL genes are present as segmental and tandem duplicates. qRT-PCR expression analysis showed that JRL genes in the same subgroup were significantly downregulated in response to salicylic acid (SA), abscisic acid (ABA), and methyl jasmonate (MeJA) treatments and significantly upregulated under low temperature, drought, and salt stress; they also exhibited tissue-specific expression patterns. Subcellular localization experiments revealed that PeJRL04 and PeJRL13 were localized to the cell membrane, nucleus, and cytoplasm. Three dimensional structure prediction and yeast two-hybrid assays were used to verify that PeJRL13 exists as a self-interacting homodimer in vivo. These findings provide an important reference for understanding the functions of specific moso bamboo JRL genes and for the effective selection of stress-related genes.

Proton-pumping pyrophosphatase homeolog expression is a dynamic trait in bread wheat (Triticum aestivum)

Proton-pumping pyrophosphatases (H+-PPases) have been shown to enhance biomass and yield. However, to date, there has been little work towards identify genes encoding H+-PPases in bread wheat (Triticum aestivum) (TaVPs) and limited knowledge on how the expression of these genes varies across different growth stages and tissue types. In this study, the IWGSC database was used to identify two novel TaVP genes, TaVP4 and TaVP5, and elucidate the complete homeolog sequences of the three known TaVP genes, bringing the total number of bread wheat TaVPs from 9 to 15. Gene expression levels of each TaVP homeolog were assessed using quantitative real-time PCR (qRT-PCR) in four diverse wheat varieties in terms of phenotypic traits related to high vacuolar pyrophosphatase expression. Homeolog expression was analyzed across multiple tissue types and developmental stages. Expression levels of the TaVP homeologs were found to vary significantly between varieties, tissues and plant developmental stages. During early development (Z10 and Z13), expressions of TaVP1 and TaVP2 homeologs were higher in shoot tissue than root tissue, with both shoot and root expression increasing in later developmental stages (Z22). TaVP2-D was expressed in all varieties and tissue types and was the most highly expressed homeolog at all developmental stages. Expression of the TaVP3 homeologs was restricted to developing grain (Z75), while TaVP4 homeolog expression was higher at Z22 than earlier developmental stages. Variation in TaVP4B was detected among varieties at Z22 and Z75, with Buck Atlantico (high biomass) and Scout (elite Australian cultivar) having the highest levels of expression. These findings offer a comprehensive overview of the bread wheat H+-PPase family and identify variation in TaVP homeolog expression that will be of use to improve the growth, yield, and abiotic stress tolerance of bread wheat.

The Insights of Genomic Synteny and Codon Usage Preference on Genera Demarcation of Iridoviridae Family

The members of the family Iridoviridae are large, double-stranded DNA viruses that infect various hosts, including both vertebrates and invertebrates. Although great progress has been made in genomic and phylogenetic analyses, the adequacy of the existing criteria for classification within the Iridoviridae family remains unknown. In this study, we redetermined 23 Iridoviridae core genes by re-annotation, core-pan analysis and local BLASTN search. The phylogenetic tree based on the 23 re-annotated core genes (Maximum Likelihood, ML-Tree) and amino acid sequences (composition vector, CV-Tree) were found to be consistent with previous reports. Furthermore, the information provided by synteny analysis and codon usage preference (relative synonymous codon usage, correspondence analysis, ENC-plot and Neutrality plot) also supports the phylogenetic relationship. Collectively, our results will be conducive to understanding the genera demarcation within the Iridoviridae family based on genomic synteny and component (codon usage preference) and contribute to the existing taxonomy methods for the Iridoviridae family.

Genome-Wide Analysis of Dof Genes and Their Response to Abiotic Stress in Rose (Rosa chinensis)

Dof (DNA binding with one finger) proteins play important roles in plant development and defense regulatory networks. In the present study, we report a genome-wide analysis of rose Dof genes (RchDof), including phylogenetic inferences, gene structures, chromosomal locations, gene duplications, and expression diversity. A total of 24 full-length RchDof genes were identified in Rosa chinensis, which were assigned to nine distinct subgroups. These RchDof genes were unevenly distributed on rose chromosomes. The genome-scale analysis of synteny indicated that segmental duplication events may have played a major role in the evolution of the RchDof gene family. Analysis of cis-acting elements revealed putative functions of Dofs in rose during development as well as under numerous biotic and abiotic stress conditions. Moreover, the expression profiles derived from qRT-PCR experiments demonstrated distinct expression patterns in various tissues, and gene expression divergence existed among the duplicated RchDof genes, suggesting a fundamentally functional divergence of the duplicated Dof paralogs in rose. The gene expression analysis of RchDofs under drought and salt stress conditions was also performed. The present study offered novel insights into the evolution of RchDofs and can aid in the further functional characterization of its candidate genes.

Genomic Organization, Phylogenetic Comparison, and Differential Expression of the Nuclear Factor-Y Gene Family in Apple (Malus Domestica)

The nuclear factor Y (NF-Y) as a transcription factor plays an important role in plants growth and development, and response to stress. However, few genome-wide analyzes and functional research of the NF-Y family has been undertaken in apple (Malus domestica Borkh.) so far. In this study, we comprehensively identified the 43 MdNF-Y genes in apple, which dispersedly distributed among the three subgroups based on their sequence alignment analysis, including 11 MdNF-YAs, 22 MdNF-YBs and 10 MdNF-YCs. The members in the same subgroups had similar evolution relationships, gene structures, and conserved motifs. The gene duplication analysis suggested that all the genes were dispersed followed by 27 segmental duplication. Moreover, based on synteny analysis of MdNF-Ys with eight plant species results suggested that some ortholog genes were preserved during the evolution of these species. Cis-element analysis showed potential functions of MdNF-Ys in apple growth and development and responded to abiotic stress. Furthermore, the interaction among MdNF-Ys protein were investigated in yeast two-hybrid assays. The expression patterns of MdNF-Ys in tissue-specific response reveled divergence and might play important role in apple growth and development. Subsequently, whole MdNF-Y genes family was carried out for RT-PCR in response to five abiotic stress (ABA, drought, heat, cold, and salinity) to identify their expression patterns. Taken together, our study will provide a foundation for the further study to the molecular mechanism of apple in growing development and response to abiotic stresses.

Comparative Genomics Studies on the dmrt Gene Family in Fish

Doublesex and mab-3-related transcription factor (dmrt) genes are widely distributed across various biological groups and play critical roles in sex determination and neural development. Here, we applied bioinformatics methods to exam cross-species changes in the dmrt family members and evolutionary relationships of the dmrt genes based on genomes of 17 fish species. All the examined fish species have dmrt1–5 while only five species contained dmrt6. Most fish harbored two dmrt2 paralogs (dmrt2a and dmrt2b), with dmrt2b being unique to fish. In the phylogenetic tree, 147 DMRT are categorized into eight groups (DMRT1–DMRT8) and then clustered in three main groups. Selective evolutionary pressure analysis indicated purifying selections on dmrt1–3 genes and the dmrt1–3–2(2a) gene cluster. Similar genomic conservation patterns of the dmrt1–dmrt3–dmrt2(2a) gene cluster with 20-kb upstream/downstream regions in fish with various sex-determination systems were observed except for three regions with remarkable diversity. Synteny analysis revealed that dmrt1, dmrt2a, dmrt2b, and dmrt3–5 were relatively conserved in fish during the evolutionary process. While dmrt6 was lost in most species during evolution. The high conservation of the dmrt1–dmrt3–dmrt2(2a) gene cluster in various fish genomes suggests their crucial biological functions while various dmrt family members and sequences across fish species suggest different biological roles during evolution. This study provides a molecular basis for fish dmrt functional analysis and may serve as a reference for in-depth phylogenomics.

Genome-Wide Identification of the ABA Receptors Genes and Their Response to Abiotic Stress in Apple

The pyrabactin resistance (PYR)/PYR1-like (PYL)/regulatory components of ABA receptor (RCAR) (known as PYLs for short) have been identified and characterized as the ABA receptors in some plants. However, little is known about the details regarding PYL family genes in the apple (Malusdomestica). In this study, we identified 13 apple PYLs, termed MdPYL1-13, which could be classified into four groups according to structural features of the amino acid sequence. The gene structures and conserved motifs analysis found that the majority of MdPYLs had a similar number of exons and similar conserved motif profile in the same group. In addition, 11 gene pairs were identified to exhibit synteny by synteny analysis between the apple and Arabidopsis. Furthermore, we investigated MdPYLs transcript level in various organs of the red-fleshed apple (Malussieversii f. Neidzwetzkyana (Dieck) Langenf) ‘Xinjiang No.1’. The results suggested all MdPYLs within group I were expressed at relatively higher levels in all of the organs tested. However, the genes of group IV had little or no variation. Additionally, we found various hormone and stress-related cis-elements in the promoters of MdPYLs by analyzing cis-elements. Therefore, the expression levels of all MdPYLs were further detected under ABA, PEG, salt, and cold stresses in ‘Xinjiang No.1’ seedlings. We found that all MdPYLs except for MdPYL11 were upregulated by ABA treatment, 10 genes were upregulated by PEG treatment, 12 genes were upregulated by NaCl treatment, and six genes were upregulated by cold treatment (4 °C) while seven genes were downregulated. Thus, these MdPYLs might be involved in the defense against abiotic stresses. In addition, the interaction between 13 MdPYLs and two 2C protein phosphatases in the apple (MdPP2C65 and MdPP2C72) was investigated in yeast two-hybrid assays. These results suggested that MdPYLs may bind to MdPP2C65 and MdPP2C72 in different manners and with different intensity. Our studies provide useful information for further investigating and researching the regulatory mechanisms of PYL family genes in response to abiotic stresses in the apple.

Cinnamic Acid and Sorbic acid Conversion Are Mediated by the Same Transcriptional Regulator in Aspergillus niger

Cinnamic acid is an aromatic compound commonly found in plants and functions as a central intermediate in lignin synthesis. Filamentous fungi are able to degrade cinnamic acid through multiple metabolic pathways. One of the best studied pathways is the non-oxidative decarboxylation of cinnamic acid to styrene. In Aspergillus niger, the enzymes cinnamic acid decarboxylase (CdcA, formally ferulic acid decarboxylase) and the flavin prenyltransferase (PadA) catalyze together the non-oxidative decarboxylation of cinnamic acid and sorbic acid. The corresponding genes, cdcA and padA, are clustered in the genome together with a putative transcription factor previously named sorbic acid decarboxylase regulator (SdrA). While SdrA was predicted to be involved in the regulation of the non-oxidative decarboxylation of cinnamic acid and sorbic acid, this was never functionally analyzed. In this study, A. niger deletion mutants of sdrA, cdcA, and padA were made to further investigate the role of SdrA in cinnamic acid metabolism. Phenotypic analysis revealed that cdcA, sdrA and padA are exclusively involved in the degradation of cinnamic acid and sorbic acid and not required for other related aromatic compounds. Whole genome transcriptome analysis of ΔsdrA grown on different cinnamic acid related compounds, revealed additional target genes, which were also clustered with cdcA, sdrA, and padA in the A. niger genome. Synteny analysis using 30 Aspergillus genomes demonstrated a conserved cinnamic acid decarboxylation gene cluster in most Aspergilli of the Nigri clade. Aspergilli lacking certain genes in the cluster were unable to grow on cinnamic acid, but could still grow on related aromatic compounds, confirming the specific role of these three genes for cinnamic acid metabolism of A. niger.

Identification and characterization of BoPUB3: a novel interaction protein with S-locus receptor kinase in Brassica oleracea L

Armadillo repeat containing 1 (ARC1) is phosphorylated by S-locus receptor kinase (SRK) and functions as a positive regulator in self-incompatibility response of Brassica. However, ARC1 only causes partial breakdown of the self-incompatibility response, and other SRK downstream factors may also participate in the self-incompatibility signaling pathway. In the present study, to search for SRK downstream targets, a plant U-box protein 3 (BoPUB3) was identified from the stigma of Brassica oleracea L. BoPUB3 was highly expressed in the stigma, and its expression was increased with the stigma development and reached to the highest level in the mature-stage stigma. BoPUB3, a 76.8-kDa protein with 697 amino acids, is a member of the PUB-ARM family and contains three domain characteristics of BoARC1, including a U-box N-terminal domain, a U-box motif, and a C-terminal arm repeat domain. The phylogenic tree showed that BoPUB3 was close to BoARC1. The synteny analysis revealed that B. oleracea chromosomal region containing BoPUB3 had high synteny with the Arabidopsis thaliana chromosomal region containing AtPUB3 (At3G54790). In addition, the subcellular localization analysis showed that BoPUB3 primarily localized in the plasma membrane and also in the cytoplasm. The combination of the yeast two-hybrid and in vitro binding assay showed that both BoPUB3 and BoARC1 could interact with SRK kinase domain, and SRK showed much higher level of β-galactosidase activity in its interaction with BoPUB3 than with BoARC1. These results implied that BoPUB3 is a novel interactor with SRK, which lays a basis for further research on whether PUB3 participates in the self-incompatibility signaling pathway.

Identification, Molecular Cloning, and Functional Characterization of a Wheat UDP-Glucosyltransferase Involved in Resistance to Fusarium Head Blight and to Mycotoxin Accumulation

Plant uridine diphosphate (UDP)-glucosyltransferases (UGT) catalyze the glucosylation of xenobiotic, endogenous substrates and phytotoxic agents produced by pathogens such as mycotoxins. The Bradi5g03300 UGT-encoding gene from the model plant Brachypodium distachyon was previously shown to confer tolerance to the mycotoxin deoxynivalenol (DON) through glucosylation into DON 3-O-glucose (D3G). This gene was shown to be involved in early establishment of quantitative resistance to Fusarium Head Blight, a major disease of small-grain cereals. In the present work, using a translational biology approach, we identified and characterized a wheat candidate gene, Traes_2BS_14CA35D5D, orthologous to Bradi5g03300 on the short arm of chromosome 2B of bread wheat (Triticum aestivum L.). We showed that this UGT-encoding gene was highly inducible upon infection by a DON-producing Fusarium graminearum strain while not induced upon infection by a strain unable to produce DON. Transformation of this wheat UGT-encoding gene into B. distachyon revealed its ability to confer FHB resistance and root tolerance to DON as well as to potentially conjugate DON into D3G in planta and its impact on total DON reduction. In conclusion, we provide a UGT-encoding candidate gene to include in selection process for FHB resistance.

A Comprehensive Analysis of Transcript-Supported De Novo Genes in Saccharomyces sensu stricto Yeasts

Novel genes arising from random DNA sequences (de novo genes) have been suggested to be widespread in the genomes of different organisms. However, our knowledge about the origin and evolution of de novo genes is still limited. To systematically understand the general features of de novo genes, we established a robust pipeline to analyze >20,000 transcript-supported coding sequences (CDSs) from the budding yeast Saccharomyces cerevisiae. Our analysis pipeline combined phylogeny, synteny, and sequence alignment information to identify possible orthologs across 20 Saccharomycetaceae yeasts and discovered 4,340 S. cerevisiae-specific de novo genes and 8,871 S. sensu stricto-specific de novo genes. We further combine information on CDS positions and transcript structures to show that >65% of de novo genes arose from transcript isoforms of ancient genes, especially in the upstream and internal regions of ancient genes. Fourteen identified de novo genes with high transcript levels were chosen to verify their protein expressions. Ten of them, including eight transcript isoform-associated CDSs, showed translation signals and five proteins exhibited specific cytosolic localizations. Our results suggest that de novo genes frequently arise in the S. sensu stricto complex and have the potential to be quickly integrated into ancient cellular network.

Genome-Wide Characterization of Heat-Shock Protein 70s from Chenopodium quinoa and Expression Analyses of Cqhsp70s in Response to Drought Stress

Heat-shock proteins (HSPs) are ubiquitous proteins with important roles in response to biotic and abiotic stress. The 70-kDa heat-shock genes (Hsp70s) encode a group of conserved chaperone proteins that play central roles in cellular networks of molecular chaperones and folding catalysts across all the studied organisms including bacteria, plants and animals. Several Hsp70s involved in drought tolerance have been well characterized in various plants, whereas no research on Chenopodium quinoa HSPs has been completed. Here, we analyzed the genome of C. quinoa and identified sixteen Hsp70 members in quinoa genome. Phylogenetic analysis revealed the independent origination of those Hsp70 members, with eight paralogous pairs comprising the Hsp70 family in quinoa. While the gene structure and motif analysis showed high conservation of those paralogous pairs, the synteny analysis of those paralogous pairs provided evidence for expansion coming from the polyploidy event. With several subcellular localization signals detected in CqHSP70 protein paralogous pairs, some of the paralogous proteins lost the localization information, indicating the diversity of both subcellular localizations and potential functionalities of those HSP70s. Further gene expression analyses revealed by quantitative polymerase chain reaction (qPCR) analysis illustrated the significant variations of Cqhsp70s in response to drought stress. In conclusion, the sixteen Cqhsp70s undergo lineage-specific expansions and might play important and varied roles in response to drought stress.

Constructing a 'Chromonome' of Yellowtail (Seriola quinqueradiata) for Comparative Analysis of Chromosomal Rearrangements

To investigate chromosome evolution in fish species, we newly mapped 181 markers that allowed us to construct a yellowtail (Seriola quinqueradiata) radiation hybrid (RH) physical map with 1,713 DNA markers, which was far denser than a previous map, and we anchored the de novo assembled sequences onto the RH physical map. Finally, we mapped a total of 13,977 expressed sequence tags (ESTs) on a genome sequence assembly aligned with the physical map. Using the high-density physical map and anchored genome sequences, we accurately compared the yellowtail genome structure with the genome structures of five model fishes to identify characteristics of the yellowtail genome. Between yellowtail and Japanese medaka (Oryzias latipes), almost all regions of the chromosomes were conserved and some blocks comprising several markers were translocated. Using the genome information of the spotted gar (Lepisosteus oculatus) as a reference, we further documented syntenic relationships and chromosomal rearrangements that occurred during evolution in four other acanthopterygian species (Japanese medaka, zebrafish, spotted green pufferfish and three-spined stickleback). The evolutionary chromosome translocation frequency was 1.5-2-times higher in yellowtail than in medaka, pufferfish, and stickleback.

Comprehensive Analysis of the CDPK-SnRK Superfamily Genes in Chinese Cabbage and Its Evolutionary Implications in Plants

The CDPK-SnRK (calcium-dependent protein kinase/Snf1-related protein kinase) gene superfamily plays important roles in signaling pathways for disease resistance and various stress responses, as indicated by emerging evidence. In this study, we constructed comparative analyses of gene structure, retention, expansion, whole-genome duplication (WGD) and expression patterns of CDPK-SnRK genes in Brassica rapa and their evolution in plants. A total of 49 BrCPKs, 14 BrCRKs, 3 BrPPCKs, 5 BrPEPRKs, and 56 BrSnRKs were identified in B. rapa. All BrCDPK-SnRK proteins had highly conserved kinase domains. By statistical analysis of the number of CDPK-SnRK genes in each species, we found that the expansion of the CDPK-SnRK gene family started from angiosperms. Segmental duplication played a predominant role in CDPK-SnRK gene expansion. The analysis showed that PEPRK was more preferentially retained than other subfamilies and that CPK was retained similarly to SnRK. Among the CPKs and SnRKs, CPKIII and SnRK1 genes were more preferentially retained than other groups. CRK was closest to CPK, which may share a common evolutionary origin. In addition, we identified 196 CPK genes and 252 SnRK genes in 6 species, and their different expansion and evolution types were discovered. Furthermore, the expression of BrCDPK-SnRK genes is dynamic in different tissues as well as in response to abiotic stresses, demonstrating their important roles in development in B. rapa. In summary, this study provides genome-wide insight into the evolutionary history and mechanisms of CDPK-SnRK genes following whole-genome triplication in B. rapa.

Evolutionary and Expression Analyses of the Apple Basic Leucine Zipper Transcription Factor Family

Transcription factors (TFs) play essential roles in the regulatory networks controlling many developmental processes in plants. Members of the basic leucine (Leu) zipper (bZIP) TF family, which is unique to eukaryotes, are involved in regulating diverse processes, including flower and vascular development, seed maturation, stress signaling, and defense responses to pathogens. The bZIP proteins have a characteristic bZIP domain composed of a DNA-binding basic region and a Leu zipper dimerization region. In this study, we identified 112 apple (Malus domestica Borkh) bZIP TF-encoding genes, termed MdbZIP genes. Synteny analysis indicated that segmental and tandem duplication events, as well as whole genome duplication, have contributed to the expansion of the apple bZIP family. The family could be divided into 11 groups based on structural features of the encoded proteins, as well as on the phylogenetic relationship of the apple bZIP proteins to those of the model plant Arabidopsis thaliana (AtbZIP genes). Synteny analysis revealed that several paired MdbZIP genes and AtbZIP gene homologs were located in syntenic genomic regions. Furthermore, expression analyses of group A MdbZIP genes showed distinct expression levels in 10 different organs. Moreover, changes in these expression profiles in response to abiotic stress conditions and various hormone treatments identified MdbZIP genes that were responsive to high salinity and drought, as well as to different phytohormones.

Genome-wide Identification and Expression Analysis of the CDPK Gene Family in Grape, Vitis spp

BACKGROUND

Calcium-dependent protein kinases (CDPKs) play vital roles in plant growth and development, biotic and abiotic stress responses, and hormone signaling. Little is known about the CDPK gene family in grapevine.

RESULTS

In this study, we performed a genome-wide analysis of the 12X grape genome (Vitis vinifera) and identified nineteen CDPK genes. Comparison of the structures of grape CDPK genes allowed us to examine their functional conservation and differentiation. Segmentally duplicated grape CDPK genes showed high structural conservation and contributed to gene family expansion. Additional comparisons between grape and Arabidopsis thaliana demonstrated that several grape CDPK genes occured in the corresponding syntenic blocks of Arabidopsis, suggesting that these genes arose before the divergence of grapevine and Arabidopsis. Phylogenetic analysis divided the grape CDPK genes into four groups. Furthermore, we examined the expression of the corresponding nineteen homologous CDPK genes in the Chinese wild grape (Vitis pseudoreticulata) under various conditions, including biotic stress, abiotic stress, and hormone treatments. The expression profiles derived from reverse transcription and quantitative PCR suggested that a large number of VpCDPKs responded to various stimuli on the transcriptional level, indicating their versatile roles in the responses to biotic and abiotic stresses. Moreover, we examined the subcellular localization of VpCDPKs by transiently expressing six VpCDPK-GFP fusion proteins in Arabidopsis mesophyll protoplasts; this revealed high variability consistent with potential functional differences.

CONCLUSIONS

Taken as a whole, our data provide significant insights into the evolution and function of grape CDPKs and a framework for future investigation of grape CDPK genes.

Patterns of evolutionary conservation of ascorbic acid-related genes following whole-genome triplication in Brassica rapa

Ascorbic acid (AsA) is an important antioxidant in plants and an essential vitamin for humans. Extending the study of AsA-related genes from Arabidopsis thaliana to Brassica rapa could shed light on the evolution of AsA in plants and inform crop breeding. In this study, we conducted whole-genome annotation, molecular-evolution and gene-expression analyses of all known AsA-related genes in B. rapa. The nucleobase-ascorbate transporter (NAT) gene family and AsA l-galactose pathway genes were also compared among plant species. Four important insights gained are that: 1) 102 AsA-related gene were identified in B. rapa and they mainly diverged 12-18 Ma accompanied by the Brassica-specific genome triplication event; 2) during their evolution, these AsA-related genes were preferentially retained, consistent with the gene dosage hypothesis; 3) the putative proteins were highly conserved, but their expression patterns varied; and 4) although the number of AsA-related genes is higher in B. rapa than in A. thaliana, the AsA contents and the numbers of expressed genes in leaves of both species are similar, the genes that are not generally expressed may serve as substitutes during emergencies. In summary, this study provides genome-wide insights into evolutionary history and mechanisms of AsA-related genes following whole-genome triplication in B. rapa.

Evolution and expression analysis of the grape (Vitis vinifera L.) WRKY gene family

WRKY proteins comprise a large family of transcription factors that play important roles in plant defence regulatory networks, including responses to various biotic and abiotic stresses. To date, no large-scale study of WRKY genes has been undertaken in grape (Vitis vinifera L.). In this study, a total of 59 putative grape WRKY genes (VvWRKY) were identified and renamed on the basis of their respective chromosome distribution. A multiple sequence alignment analysis using all predicted grape WRKY genes coding sequences, together with those from Arabidopsis thaliana and tomato (Solanum lycopersicum), indicated that the 59 VvWRKY genes can be classified into three main groups (I-III). An evaluation of the duplication events suggested that several WRKY genes arose before the divergence of the grape and Arabidopsis lineages. Moreover, expression profiles derived from semiquantitative PCR and real-time quantitative PCR analyses showed distinct expression patterns in various tissues and in response to different treatments. Four VvWRKY genes showed a significantly higher expression in roots or leaves, 55 responded to varying degrees to at least one abiotic stress treatment, and the expression of 38 were altered following powdery mildew (Erysiphe necator) infection. Most VvWRKY genes were downregulated in response to abscisic acid or salicylic acid treatments, while the expression of a subset was upregulated by methyl jasmonate or ethylene treatments.

Genome-wide analysis of respiratory burst oxidase homologs in grape (Vitis vinifera L.)

Plant respiratory burst oxidase homolog (rboh) genes appear to play crucial roles in plant development, defense reactions and hormone signaling. In this study, a total of seven rboh genes from grape were identified and characterized. Genomic structure and predicted protein sequence analysis indicated that the sequences of plant rboh genes are highly conserved. Synteny analysis demonstrated that several Vvrboh genes were found in corresponding syntenic blocks of Arabidopsis, suggesting that these genes arose before the divergence of the respective lineages. The expression pattern of Vvrboh genes in different tissues was assessed by qRT-PCR and two were constitutively expressed in all tissues tested. The expression profiles were similarly analyzed following exposure to various stresses and hormone treatments. It was shown that the expression levels of VvrbohA, VvrbohB and VvrbohC1 were significantly increased by salt and drought treatments. VvrbohB, VvrbohC2, and VvrbohD exhibited a dramatic up-regulation after powdery mildew (Uncinula necator (Schw.) Burr.) inoculation, while VvrbohH was down-regulated. Finally, salicylic acid treatment strongly stimulated the expression of VvrbohD and VvrbohH, while abscisic acid treatment induced the expression of VvrbohB and VvrbohH. These results demonstrate that the expression patterns of grape rboh genes exhibit diverse and complex stress-response expression signatures.