Selective modes determine evolutionary rates, gene compactness and expression patterns in Brassica

An evolutionary dynamics analysis of the plant DEK gene family reveals the role of BnaA02g08940D in drought tolerance

DEK is a chromatin protein that interacts with DNA to influence chromatin formation, thereby affecting plant growth, development, and stress response. This study investigates the molecular evolution of the DEK family in plants, with a particular focus on the Brassica species. A total of 127 DEK genes were identified in 34 plants and classified into seven groups based on the phylogenetic analysis. The distribution of motifs and gene structure is similar within each group, indicating a high degree of conservation. The results of the collinearity analysis indicated that the DEK protein has undergone a certain degree of evolutionary conservation. The expansion of the DEK family is primarily attributable to whole-genome duplication (WGD) or segmental duplication events. The DEK protein has undergone purification during its evolutionary history, and several positively selected sites have been identified. Moreover, the examination of cis-acting elements and expression patterns revealed that the BnDEKs play a significant role in plant growth and stress response. The protein-protein interaction network identified several noteworthy proteins that interact with DEK. These analyses enhance our comprehension of the DEK gene family and establish the foundation for additional validation of its function. Further research demonstrated that the overexpression of one DEK family member, BnaA02g08940D, enhanced the transgenic Arabidopsis tolerance to drought and osmosis. This indicates that the DEK family may respond when plants are subjected to drought stress, thereby strengthening the plant's resilience.

Contrasting levels of transcriptome-wide SNP diversity and adaptive molecular variation among conifers

Adaptive convergence can arise when response to natural selection involves shared molecular or functional mechanisms among multiple taxa. Conifers are archaic species of ancient origin with delayed sexual maturity related to their woody perennial nature. Thus, they represent a relevant plant group to assess if convergence from selection may have become disconnected between molecular and functional levels. In this purpose, transcriptome-wide SNP diversity was assessed in seven partially sympatric and reproductively isolated conifer species (118 individuals from 67 populations) populating the temperate and boreal forests of northeastern North America. SNP diversity was found highly heterogeneous among species, which would relate to variation in species-specific demography and history. Rapidly evolving genes with signatures of positive selection were identified, and their relative abundance among species reflected differences in transcriptome-wide SNP diversity. The analysis of sequence homology also revealed very limited convergence among taxa in spite of sampling same tissues at same age. However, convergence increased gradually at the levels of gene families and biological processes, which were largely related to stress response and regulatory mechanisms in all species. Given their multiple small to large gene families and long time since inception, conifers may have had sufficient gene network flexibility and gene functional redundancy for evolving alternative adaptive genes for similar metabolic responses to environmental selection pressures. Despite a long divergence time of ~350 Mya between conifers and Angiosperms, we also uncovered a set of 17 key genes presumably under positive selection in both lineages.

Genome-Wide Profiling of the ACTIN Gene Family and Its Implications for Agronomic Traits in Brassica napus: A Bioinformatics Study

ACTINs are key structural proteins in plants, which form the actin cytoskeleton and are engaged in numerous routine cellular processes. Meanwhile, ACTIN, recognized as a housekeeping gene, has not yet been thoroughly investigated in Brassica napus. The current research has led to the detection of 69 actin genes in B. napus, which were organized into six distinct subfamilies on the basis of phylogenetic relationships. Functional enrichment analysis, along with the construction of protein interaction networks, suggested that BnACTINs play roles in Preserving cell morphology and facilitating cytoplasmic movement, plant development, and adaptive responses to environmental stress. Moreover, the BnACTIN genes presented a wide range of expression levels among different tissues, whereas the majority experienced a substantial increase in expression when subjected to various abiotic stresses, demonstrating a pronounced sensitivity to abiotic factors. Furthermore, association mapping analysis indicated that some BnACTINs potentially affected certain key agronomic traits. Overall, our research deepens the knowledge of BnACTIN genes, promotes the cultivation of improved B. napus strains, and lays the groundwork for subsequent functional research.

Genome-wide identification of the E-class gene family in wheat: evolution, expression, and interaction

Introduction

Wheat (Triticum aestivum L.) is among themost important crop worldwide. Given a growing population and changing climate, enhancing wheat yield is of great importance. Yield is closely associated with flower and spike development, and E-class genes play important roles in the flower and kernel development of plants. Currently, the absence of systematic analysis on the E gene family hinders our comprehension of their roles in plant growth and development.

Methods

Identify E-class genes based on homologous sequence searches. Analyze the identified E-class genes through a series of gene family analyses. Determine the expression levels of wheat E-class genes by searching public databases. Validate the functions of these genes by transforming them into Arabidopsis. Finally, determine the interactions between the genes through yeast two-hybrid experiments.

Results

Fifteen E-class genes (TaEs) were identified in common wheat. Nine E-class genes were detected in five ancestral/closely related species, including one in Aegilops tauschii (AtE), one in T. Urartu (TuEs), two in T. turgidum (TtEs), two in T. dicoccoides (TdEs), and three in T. spelta (TsEs). The 24 E-class genes were classified into three subgroups using a phylogenetic approach. All genes were highly expressed in spikes, and most were only highly expressed at the floret meristem stage. The effects of TaSEP5-A on flowering and growth cycles were confirmed in homologous mutants and transgenic Arabidopsis thaliana. The E-class genes were able to regulate the growth cycle of Arabidopsis. Finally, we confirmed the interactions between TaSEP5-A and other wheat E-class genes based on yeast two-hybrid assays.

Discussion

Our findings provide information regarding the E-class genes in wheat and will potentially promote the application of these genes in wheat improvement.

Genetic factors explaining anthocyanin pigmentation differences

BACKGROUND

Anthocyanins are important contributors to coloration across a wide phylogenetic range of plants. Biological functions of anthocyanins span from reproduction to protection against biotic and abiotic stressors. Owing to a clearly visible phenotype of mutants, the anthocyanin biosynthesis and its sophisticated regulation have been studied in numerous plant species. Genes encoding the anthocyanin biosynthesis enzymes are regulated by a transcription factor complex comprising MYB, bHLH and WD40 proteins.

RESULTS

A systematic comparison of anthocyanin-pigmented vs. non-pigmented varieties was performed within numerous plant species covering the taxonomic diversity of flowering plants. The literature was screened for cases in which genetic factors causing anthocyanin loss were reported. Additionally, transcriptomic data sets from four previous studies were reanalyzed to determine the genes possibly responsible for color variation based on their expression pattern. The contribution of different structural and regulatory genes to the intraspecific pigmentation differences was quantified. Differences concerning transcription factors are by far the most frequent explanation for pigmentation differences observed between two varieties of the same species. Among the transcription factors in the analyzed cases, MYB genes are significantly more prone to account for pigmentation differences compared to bHLH or WD40 genes. Among the structural genes, DFR genes are most often associated with anthocyanin loss.

CONCLUSIONS

These findings support previous assumptions about the susceptibility of transcriptional regulation to evolutionary changes and its importance for the evolution of novel coloration phenotypes. Our findings underline the particular significance of MYBs and their apparent prevalent role in the specificity of the MBW complex.

Comparative genomics, microsynteny, ancestral state reconstruction and selection pressure analysis across distinctive genomes and sub-genomes of Brassicaceae for analysis of evolutionary history of VQ gene family

Any unfavorable condition that affects the metabolism, growth, or development of plants is considered plant stress. The molecular response of plants towards abiotic stresses involves signaling to cellular components, repressing transcription factors, and subsequently induced metabolic changes. Most valine-glutamine (VQ) motif-containing genes in plants encode regulatory proteins that interact with transcription factors and modulate their activity as transcription regulators. Several VQ proteins regulate plant development and stress responses. In spite of the functional importance of VQs, there is relatively little information about their evolutionary history in Brassicaceae or beyond. Brassicaceae is characterized by paleoploidy, mesopolyploidy, and neopolyploidy, offering a resource for studying evolution and diversification. In current study we performed phylogeny of the VQ gene family along with comparative genomics, microsynteny and evolutionary rates analysis across seven species of Brassicaceae. Our findings revealed the following; (1) a large segmental duplication in the shared common ancestor of the family Brassicaceae, resulted in paralogies of VQ1-VQ10, VQ15-VQ24, VQ16-VQ23, VQ17-VQ25, VQ18-VQ26, VQ22-VQ27; (2) chromosomal mapping revealed diverse distributions of the gene family; (3) duplicated segments undergo varying degrees of retention and loss; and (4) Out of the 12 paralogous members, most of the genes are under purifying selection. However, VQ23 in Brassicaceae stands out as it is under positive selection, indicating the need for further investigation. Overall, our results clearly establish that the ancestral VQ1/VQ10, VQ15/VQ24, VQ16/VQ23, VQ17/VQ25, VQ18/VQ26, VQ22/VQ27 genes duplicated in shared common ancestor of Brassicaceae.

Supplementary Information

The online version contains supplementary material available at 10.1007/s12298-023-01347-z.

Review: WRKY transcription factors: Understanding the functional divergence

WRKY transcription factors (TFs) play crucial roles in the growth and development of plants and their response to environmental changes. WRKY TFs have been detected in sequenced plant genomes. The functions and regulatory networks of many WRKY TFs, especially from Arabidopsis thaliana (AtWRKY TFs), have been revealed, and the origin of WRKY TFs in plants is clear. Nonetheless, the relationship between WRKY TFs function and classification is unclear. Furthermore, the functional divergence of homologous WRKY TFs in plants is unclear. In this review, WRKY TFs were explored based on WRKY-related literature published from 1994 to 2022. WRKY TFs were identified in 234 species at the genome and transcriptome levels. The biological functions of ∼ 71 % of AtWRKY TFs were uncovered. Although functional divergence occurred in homologous WRKY TFs, different WRKY TF groups had no preferential function.

Comparative Analysis Reveals Different Evolutionary Fates and Biological Functions in Wheat Duplicated Genes (Triticum aestivum L.)

Wheat (Triticum aestivum L.) is a staple food crop that provides 20% of total human calorie consumption. Gene duplication has been considered to play an important role in evolution by providing new genetic resources. However, the evolutionary fates and biological functions of the duplicated genes in wheat remain to be elucidated. In this study, the resulting data showed that the duplicated genes evolved faster with shorter gene lengths, higher codon usage bias, lower expression levels, and higher tissue specificity when compared to non-duplicated genes. Our analysis further revealed functions of duplicated genes in various biological processes with significant enrichment to environmental stresses. In addition, duplicated genes derived from dispersed, proximal, tandem, transposed, and whole-genome duplication differed in abundance, evolutionary rate, gene compactness, expression pattern, and genetic diversity. Tandem and proximal duplicates experienced stronger selective pressure and showed a more compact gene structure with diverse expression profiles than other duplication modes. Moreover, genes derived from different duplication modes showed an asymmetrical evolutionary pattern for wheat A, B, and D subgenomes. Several candidate duplication hotspots associated with wheat domestication or polyploidization were characterized as potential targets for wheat molecular breeding. Our comprehensive analysis revealed the evolutionary trajectory of duplicated genes and laid the foundation for future functional studies on wheat.

Comparative analysis of Diospyros (Ebenaceae) plastomes: Insights into genomic features, mutational hotspots, and adaptive evolution

Diospyros (Ebenaceae) is a widely distributed genus of trees and shrubs from pantropical to temperate regions, with numerous species valued for their fruits (persimmons), timber, and medicinal values. However, information regarding their plastomes and chloroplast evolution is scarce. The present study performed comparative genomic and evolutionary analyses on plastomes of 45 accepted Diospyros species, including three newly sequenced ones. Our study showed a highly conserved genomic structure across the Diospyros species, with 135-136 encoding genes, including 89 protein-coding genes, 1-2 pseudogenes (Ψycf1 for all, Ψrps19 for a few), 37 tRNA genes and 8 rRNA genes. Comparative analysis of Diospyros identified three intergenic regions (ccsA-ndhD, rps16-psbK and petA-psbJ) and five genes (rpl33, rpl22, petL, psaC and rps15) as the mutational hotspots in these species. Phylogenomic analysis identified the phylogenetic position of three newly sequenced ones and well supported a monophylogenetic (sub)temperate taxa and four clades in the pantropical taxa. The analysis codon usage identified 30 codons with relative synonymous codon usage (RSCU) values >1 and 29 codons ending with A and U bases. A total of three codons (UUA, GCU, and AGA) with highest RSCU values were identified as the optimal codons. Effective number of codons (ENC)-plot indicated the significant role of mutational pressure in shaping codon usage, while most protein-coding genes in Diospyros experienced relaxed purifying selection (d N/d S < 1). Additionally, the psbH gene showed positive selection (d N/d S > 1) in the (sub)temperate species. Thus, the results provide a meaningful foundation for further elaborating Diospyros's genetic architecture and taxonomy, enriching genetic diversity and conserving genetic resources.

Genome-wide characterization of ubiquitin-conjugating enzyme gene family explores its genetic effects on the oil content and yield of Brassica napus

Ubiquitin-conjugating enzyme (UBC) is a critical part of the ubiquitin-proteasome pathway and plays crucial roles in growth, development and abiotic stress response in plants. Although UBC genes have been detected in several plant species, characterization of this gene family at the whole-genome level has not been conducted in Brassica napus. In the present study, 200 putative BnUBCs were identified in B. napus, which were clustered into 18 subgroups based on phylogenetic analysis. BnUBCs within each subgroup showed relatively conserved gene architectures and motifs. Moreover, the gene expression patterns in various tissues as well as the identification of cis-acting regulatory elements in BnUBC promoters suggested further investigation of their potential functions in plant growth and development. Furthermore, three BnUBCs were predicted as candidate genes for regulating agronomic traits related to oil content and yield through association mapping. In conclusion, this study provided a wealth of information on the UBC family in B. napus and revealed their effects on oil content and yield, which will aid future functional research and genetic breeding of B. napus.

The evolutionary patterns, expression profiles, and genetic diversity of expanded genes in barley

Gene duplication resulting from whole-genome duplication (WGD), small-scale duplication (SSD), or unequal hybridization plays an important role in the expansion of gene families. Gene family expansion can also mediate species formation and adaptive evolution. Barley (Hordeum vulgare) is the world's fourth largest cereal crop, and it contains valuable genetic resources due to its ability to tolerate various types of environmental stress. In this study, 27,438 orthogroups in the genomes of seven Poaceae were identified, and 214 of them were significantly expanded in barley. The evolutionary rates, gene properties, expression profiles, and nucleotide diversity between expanded and non-expanded genes were compared. Expanded genes evolved more rapidly and experienced lower negative selection. Expanded genes, including their exons and introns, were shorter, they had fewer exons, their GC content was lower, and their first exons were longer compared with non-expanded genes. Codon usage bias was also lower for expanded genes than for non-expanded genes; the expression levels of expanded genes were lower than those of non-expanded genes, and the expression of expanded genes showed higher tissue specificity than that of non-expanded genes. Several stress-response-related genes/gene families were identified, and these genes could be used to breed barley plants with greater resistance to environmental stress. Overall, our analysis revealed evolutionary, structural, and functional differences between expanded and non-expanded genes in barley. Additional studies are needed to clarify the functions of the candidate genes identified in our study and evaluate their utility for breeding barley plants with greater stress resistance.

Selective modes affect gene feature and function differentiation of tetraploid Brassica species in their evolution and domestication

The genus Brassica contains a diverse group of important vegetables and oilseed crops. Genome sequencing has been completed for the six species (B. rapa, B. oleracea, B. nigra, B. carinata, B. napus, and B. juncea) in U's triangle model. The purpose of the study is to investigate whether positively and negatively selected genes (PSGs and NSGs) affect gene feature and function differentiation of Brassica tetraploids in their evolution and domestication. A total of 9,701 PSGs were found in the A, B and C subgenomes of the three tetraploids, of which, a higher number of PSGs were identified in the C subgenome as comparing to the A and B subgenomes. The PSGs of the three tetraploids had more tandem duplicated genes, higher single copy, lower multi-copy, shorter exon length and fewer exon number than the NSGs, suggesting that the selective modes affected the gene feature of Brassica tetraploids. The PSGs of all the three tetraploids enriched in a few common KEGG pathways relating to environmental adaption (such as Phenylpropanoid biosynthesis, Riboflavin metabolism, Isoflavonoid biosynthesis, Plant-pathogen interaction and Tropane, piperidine and pyridine alkaloid biosynthesis) and reproduction (Homologous recombination). Whereas, the NSGs of the three tetraploids significantly enriched in dozens of biologic processes and pathways without clear relationships with evolution. Moreover, the PSGs of B. carinata were found specifically enriched in lipid biosynthesis and metabolism which possibly contributed to the domestication of B. carinata as an oil crop. Our data suggest that selective modes affected the gene feature of Brassica tetraploids, and PSGs contributed in not only the evolution but also the domestication of Brassica tetraploids.

Genome-wide identification, characterization, evolution, and expression pattern analyses of the typical thioredoxin gene family in wheat (Triticum aestivum L.)

Typical thioredoxin (TRX) plays an important role in maintaining redox balance in plants. However, the typical TRX genes in wheat still need to be comprehensively and deeply studied. In this research, a total of 48 typical TaTRX genes belonging to eight subtypes were identified via a genome-wide search in wheat, and the gene structures, protein conserved motifs, and protein 3D structures of the same subtype were very similar. Evolutionary analysis showed that there are two pairs of tandem duplication genes and 14 clusters of segmental duplication genes in typical TaTRX family members; TaTRX15, TaTRX36, and TaTRX42 had positive selection compared with the orthologs of their ancestral species; rice and maize have 11 and 13 orthologous typical TRXs with wheat, respectively. Gene Ontology (GO) analysis indicated that typical TaTRXs were involved in maintaining redox homeostasis in wheat cells. Estimation of ROS content, determination of antioxidant enzyme activity, and gene expression analysis in a line overexpressing one typical TaTRX confirmed that TRX plays an important role in maintaining redox balance in wheat. A predictive analysis of cis-acting elements in the promoter region showed that typical TaTRXs were extensively involved in various hormone metabolism and response processes to stress. The results predicted using public databases or verified using RT-qPCR show that typical TaTRXs were able to respond to biotic and abiotic stresses, and their expression in wheat was spatiotemporal. A total of 16 wheat proteins belonging to four different families interacting with typical TaTRXs were predicted. The above comprehensive analysis of typical TaTRX genes can enrich our understanding of this gene family in wheat and provide valuable insights for further gene function research.

Chloroplast gene expression level is negatively correlated with evolutionary rates and selective pressure while positively with codon usage bias in Ophioglossum vulgatum L

BACKGROUND

Characterization of the key factors determining gene expression level has been of significant interest. Previous studies on the relationship among evolutionary rates, codon usage bias, and expression level mostly focused on either nuclear genes or unicellular/multicellular organisms but few in chloroplast (cp) genes. Ophioglossum vulgatum is a unique fern and has important scientific and medicinal values. In this study, we sequenced its cp genome and transcriptome to estimate the evolutionary rates (dN and dS), selective pressure (dN/dS), gene expression level, codon usage bias, and their correlations.

RESULTS

The correlation coefficients between dN, dS, and dN/dS, and Transcripts Per Million (TPM) average values were -0.278 (P = 0.027 < 0.05), -0.331 (P = 0.008 < 0.05), and -0.311 (P = 0.013 < 0.05), respectively. The codon adaptation index (CAI) and tRNA adaptation index (tAI) were significantly positively correlated with TPM average values (P < 0.05).

CONCLUSIONS

Our results indicated that when the gene expression level was higher, the evolutionary rates and selective pressure were lower, but the codon usage bias was stronger. We provided evidence from cp gene data which supported the E-R (E stands for gene expression level and R stands for evolutionary rate) anti-correlation.

Spruce giga-genomes: structurally similar yet distinctive with differentially expanding gene families and rapidly evolving genes

Spruces (Picea spp.) are coniferous trees widespread in boreal and mountainous forests of the northern hemisphere, with large economic significance and enormous contributions to global carbon sequestration. Spruces harbor very large genomes with high repetitiveness, hampering their comparative analysis. Here, we present and compare the genomes of four different North American spruces: the genome assemblies for Engelmann spruce (Picea engelmannii) and Sitka spruce (Picea sitchensis) together with improved and more contiguous genome assemblies for white spruce (Picea glauca) and for a naturally occurring introgress of these three species known as interior spruce (P. engelmannii × glauca × sitchensis). The genomes were structurally similar, and a large part of scaffolds could be anchored to a genetic map. The composition of the interior spruce genome indicated asymmetric contributions from the three ancestral genomes. Phylogenetic analysis of the nuclear and organelle genomes revealed a topology indicative of ancient reticulation. Different patterns of expansion of gene families among genomes were observed and related with presumed diversifying ecological adaptations. We identified rapidly evolving genes that harbored high rates of non-synonymous polymorphisms relative to synonymous ones, indicative of positive selection and its hitchhiking effects. These gene sets were mostly distinct between the genomes of ecologically contrasted species, and signatures of convergent balancing selection were detected. Stress and stimulus response was identified as the most frequent function assigned to expanding gene families and rapidly evolving genes. These two aspects of genomic evolution were complementary in their contribution to divergent evolution of presumed adaptive nature. These more contiguous spruce giga-genome sequences should strengthen our understanding of conifer genome structure and evolution, as their comparison offers clues into the genetic basis of adaptation and ecology of conifers at the genomic level. They will also provide tools to better monitor natural genetic diversity and improve the management of conifer forests. The genomes of four closely related North American spruces indicate that their high similarity at the morphological level is paralleled by the high conservation of their physical genome structure. Yet, the evidence of divergent evolution is apparent in their rapidly evolving genomes, supported by differential expansion of key gene families and large sets of genes under positive selection, largely in relation to stimulus and environmental stress response.

Molecular evolutionary analysis of the SHI/STY gene family in land plants: A focus on the Brassica species

The plant-specific SHORT INTERNODES/STYLISH (SHI/STY) proteins belong to a family of transcription factors that are involved in the formation and development of early lateral roots. However, the molecular evolution of this family is rarely reported. Here, a total of 195 SHI/STY genes were identified in 21 terrestrial plants, and the Brassica species is the focus of our research. Their physicochemical properties, chromosome location and duplication, motif distribution, exon-intron structures, genetic evolution, and expression patterns were systematically analyzed. These genes are divided into four clades (Clade 1/2/3/4) based on phylogenetic analysis. Motif distribution and gene structure are similar in each clade. SHI/STY proteins are localized in the nucleus by the prediction of subcellular localization. Collinearity analysis indicates that the SHI/STYs are relatively conserved in evolution. Whole-genome duplication is the main factor for their expansion. SHI/STYs have undergone intense purifying selection, but several positive selection sites are also identified. Most promoters of SHI/STY genes contain different types of cis-elements, such as light, stress, and hormone-responsive elements, suggesting that they may be involved in many biological processes. Protein-protein interaction predicted some important SHI/STY interacting proteins, such as LPAT4, MBOATs, PPR, and UBQ3. In addition, the RNA-seq and qRT-PCR analysis were studied in detail in rape. As a result, SHI/STYs are highly expressed in root and bud, and can be affected by Sclerotinia sclerotiorum, drought, cold, and heat stresses. Moreover, quantitative real-time PCR (qRT-PCR) analyses indicates that expression levels of BnSHI/STYs are significantly altered in different treatments (cold, salt, drought, IAA, auxin; ABA, abscisic acid; 6-BA, cytokinin). It provides a new understanding of the evolution and expansion of the SHI/STY family in land plants and lays a foundation for further research on their functions.

Genomic insights into positive selection during barley domestication

BACKGROUND

Cultivated barley (Hordeum vulgare) is widely used in animal feed, beverages, and foods and has become a model crop for molecular evolutionary studies. Few studies have examined the evolutionary fates of different types of genes in barley during the domestication process.

RESULTS

The rates of nonsynonymous substitution (Ka) to synonymous substitution (Ks) were calculated by comparing orthologous genes in different barley groups (wild vs. landrace and landrace vs. improved cultivar). The rates of evolution, properties, expression patterns, and diversity of positively selected genes (PSGs) and negatively selected genes (NSGs) were compared. PSGs evolved more rapidly, possessed fewer exons, and had lower GC content than NSGs; they were also shorter and had shorter intron, exon, and first exon lengths. Expression levels were lower, the tissue specificity of expression was higher, and codon usage bias was weaker for PSGs than for NSGs. Nucleotide diversity analysis revealed that PSGs have undergone a more severe genetic bottleneck than NSGs. Several candidate PSGs were involved in plant growth and development, which might make them as excellent targets for the molecular breeding of barley.

CONCLUSIONS

Our comprehensive analysis of the evolutionary, structural, and functional divergence between PSGs and NSGs in barley provides new insight into the evolutionary trajectory of barley during domestication. Our findings also aid future functional studies of PSGs in barley.

Triplication is the main evolutionary driving force of NLP transcription factor family in Chinese cabbage and related species

The NODULE-INCEPTION-like protein (NLP) is a plant-specific transcription factor (TF) family that plays an important role in both signal transduction and nitrate assimilation. However, the NLP gene family in Chinese cabbage (Brassica rapa) has yet to be studied. Here we identified 17, 16, and 32 NLP genes in Chinese cabbage, Brassica oleracea, and Brassica napus, respectively. We found that duplication of those NLP genes almost always originated from genome-wide duplication events. Further analysis (using Arabidopsis as a reference) revealed that the NLP family in Chinese cabbage and B. oleracea was characterized by direct expansion caused by whole-genome duplication. By contrast, indirect expansion characterized B. napus, which arose from hybridization and fusion of the two species. In addition, phylogenetic and homology analyses showed that the Brassica NLP gene family has been highly conserved in evolution. Finally, we also identified optimal codons for four studied species. Altogether, through comparative genome analysis methods, we presented compelling evidence that triplication is the main driving force for the NLP TF family's evolution in Chinese cabbage and related Brassica plants, a process evidently highly conserved. This work will help in better understanding the impact of genome-wide duplication on gene families of plants.

The Landscape of Autophagy-Related (ATG) Genes and Functional Characterization of TaVAMP727 to Autophagy in Wheat

Autophagy is an indispensable biological process and plays crucial roles in plant growth and plant responses to both biotic and abiotic stresses. This study systematically identified autophagy-related proteins (ATGs) in wheat and its diploid and tetraploid progenitors and investigated their genomic organization, structure characteristics, expression patterns, genetic variation, and regulation network. We identified a total of 77, 51, 29, and 30 ATGs in wheat, wild emmer, T. urartu and A. tauschii, respectively, and grouped them into 19 subfamilies. We found that these autophagy-related genes (ATGs) suffered various degrees of selection during the wheat’s domestication and breeding processes. The genetic variations in the promoter region of Ta2A_ATG8a were associated with differences in seed size, which might be artificially selected for during the domestication process of tetraploid wheat. Overexpression of TaVAMP727 improved the cold, drought, and salt stresses resistance of the transgenic Arabidopsis and wheat. It also promoted wheat heading by regulating the expression of most ATGs. Our findings demonstrate how ATGs regulate wheat plant development and improve abiotic stress resistance. The results presented here provide the basis for wheat breeding programs for selecting varieties of higher yield which are capable of growing in colder, drier, and saltier areas.

Distribution, expression and methylation analysis of positively selected genes provides insights into the evolution in Brassica rapa

It is believed that positive selection is one of the major evolutionary forces underlying organism phenotypic diversification. Nevertheless, the characteristics of positively selected genes (PSGs), have not been well investigated. In this study, we performed a genome-wide analysis of orthologous genes between Brassica rapa (B. rapa) and Brassica oleracea (B. oleracea), and identified 468 putative PSGs. Our data show that, (1) PSGs are enriched in plant hormone signal transduction pathway and the transcription factor family; (2) PSGs are significantly lower expressed than randomly selected non-PSGs; (3) PSGs with tissue specificity are significantly higher expressed in the callus and reproductive tissues (flower and silique) than in vegetable tissues (root, stem and leaf); (4) the proportion of PSGs is positively correlated with the number of retained triplication gene copies, but the expression level of PSGs decay with the increasing of triplication gene copies; (5) the CG and CHG methylation levels of PSGs are significantly higher in introns and UTRs than in the promoter and exon regions; (6) the percent of transposable element is in proportion to the methylation level, and DNA methylation (especially in the CG content) has the tendency to reduce the expression of PSGs. This study provides insights into the characteristics, evolution, function, expression and methylation of PSGs in B. rapa.

The evolution of the expansin gene family in Brassica species

Expansin gene (EXP) family plays important roles in plant growth and crop improvement. However, it has not been well studied in the Brassica genus that includes several important agricultural and horticultural crops. To get insight to the evolution and expansion of EXP family in Brassica, Brassica EXPs which are homologues of 35 known AtEXPs of Arabidopsis were comprehensively and systematically analyzed in the present study. In total, 340 Brassica EXPs were clustered into four groups that corresponded multiple alignment to four subfamilies of AtEXPs, with divergent conserved motifs and cis-acting elements among groups. To understand the expansion of EXP family, an integrated genomic block system was constructed among Arabidopsis and Brassica species based on 24 known ancestral karyotype blocks. Obvious gene loss, segmental duplication, tandem duplication and DNA sequence repeat events were found during the expansion of Brassica EXPs, of which the segmental duplication was possibly the major driving force. The divergence time was estimated in 1109 orthologs pairs of EXPs, revealing the divergence of Brassica EXPs from AtEXPs during ~30 MYA, and the divergence of EXPs among Brassica species during 13.50-17.94 MYA. Selective mode analysis revealed that the purifying selection was the major contributor to expansion of Brassica EXPs. This study provides new insights into the evolution and expansion of the EXP family in Brassica genus.

Rates and patterns of molecular evolution in bryophyte genomes, with focus on complex thalloid liverworts, Marchantiopsida

Plants commonly referred to as "bryophytes" belong to three major lineages of non-vascular plants: the liverworts, the hornworts and the mosses. They are unique among land plants in having a dominant haploid generation and a short-lived diploid sporophytic generation. The dynamics of selection acting on a haploid genome differs from those acting on a diploid genome: new mutations are directly exposed to selection. The general aim of this paper is to investigate the diversification rateof bryophytes - measured as silent site substitution rate representing neutral evolution (mutation rate) and the nonsynonymous to synonymous substitution rate ratio (dN/dS) representing selective evolution - and compare it with earlier studies on vascular plants. Results show that the silent site substitution rate is lower for liverworts as compared to angiosperms, but not as low as for gymnosperms. The selection pressure, measured as dN/dS, isnot remarkably lower for bryophytes as compared to other diploid dominant plants as would be expected by the masking hypothesis, indicating that other factors are more important than ploidy.

Chloroplast genome features of Moricandia arvensis (Brassicaceae), a C3-C4 intermediate photosynthetic species

Moricandia arvensis, a plant species originating from the Mediterranean, has been classified as a rare C3-C4 intermediate species, and it is a possible bridge during the evolutionary process from C3 to C4 plant photosynthesis in the family Brassicaceae. Understanding the genomic structure, gene order, and gene content of chloroplasts (cp) of such species can provide a glimpse into the evolution of photosynthesis. In the present study, we obtained a well-annotated cp genome of M. arvensis using long PacBio and short Illumina reads with a de novo assembly strategy. The M. arvensis cp genome was a quadripartite circular molecule with the length of 153,312 bp, including two inverted repeats (IR) regions of 26,196 bp, divided by a small single copy (SSC) region of 17,786 bp and a large single copy (LSC) region of 83,134 bp. We detected 112 unigenes in this genome, comprising 79 protein-coding genes, 29 tRNAs, and four rRNAs. Forty-nine long repeat sequences and 51 simple sequence repeat (SSR) loci of 15 repeat types were identified. The analysis of Ks (synonymous) and Ka (non-synonymous) substitution rates indicated that the genes associated with “subunits of ATP synthase” (atpB), “subunits of NADH-dehydrogenase” (ndhG and ndhE), and “self-replication” (rps12 and rpl16) showed relatively higher Ka/Ks values than those of the other genes. The gene content, gene order, and LSC/IR/SSC boundaries and adjacent genes of the M. arvensis cp genome were highly conserved compared to those in related C3 species. Our phylogenetic analysis demonstrated that M. arvensis was clustered into a subclade with cultivated Brassica species and Raphanus sativus, indicating that M. arvensis was not involved in an independent evolutionary origin event. These results will open the way for further studies on the evolutionary process from C3 to C4 photosynthesis and hopefully provide guidance for utilizing M. arvensis as a resource for improvinng photosynthesis efficiency in cultivated Brassica species.

A prescient evolutionary model for genesis, duplication and differentiation of MIR160 homologs in Brassicaceae

MicroRNA160 is a class of nitrogen-starvation responsive genes which governs establishment of root system architecture by down-regulating AUXIN RESPONSE FACTOR genes (ARF10, ARF16 and ARF17) in plants. The high copy number of MIR160 variants discovered by us from land plants, especially polyploid crop Brassicas, posed questions regarding genesis, duplication, evolution and function. Absence of studies on impact of whole genome and segmental duplication on retention and evolution of MIR160 homologs in descendent plant lineages prompted us to undertake the current study. Herein, we describe ancestry and fate of MIR160 homologs in Brassicaceae in context of polyploidy driven genome re-organization, copy number and differentiation. Paralogy amongst Brassicaceae MIR160a, MIR160b and MIR160c was inferred using phylogenetic analysis of 468 MIR160 homologs from land plants. The evolutionarily distinct MIR160a was found to represent ancestral form and progenitor of MIR160b and MIR160c. Chronology of evolutionary events resulting in origin and diversification of genomic loci containing MIR160 homologs was delineated using derivatives of comparative synteny. A prescient model for causality of segmental duplications in establishment of paralogy in Brassicaceae MIR160, with whole genome duplication accentuating the copy number increase, is being posited in which post-segmental duplication events viz. differential gene fractionation, gene duplications and inversions are shown to drive divergence of chromosome segments. While mutations caused the diversification of MIR160a, MIR160b and MIR160c, duplicated segments containing these diversified genes suffered gene rearrangements via gene loss, duplications and inversions. Yet the topology of phylogenetic and phenetic trees were found congruent suggesting similar evolutionary trajectory. Over 80% of Brassicaceae genomes and subgenomes showed a preferential retention of single copy each of MIR160a, MIR160b and MIR160c suggesting functional relevance. Thus, our study provides a blue-print for reconstructing ancestry and phylogeny of MIRNA gene families at genomics level and analyzing the impact of polyploidy on organismal complexity. Such studies are critical for understanding the molecular basis of agronomic traits and deploying appropriate candidates for crop improvement.

Asymmetric epigenome maps of subgenomes reveal imbalanced transcription and distinct evolutionary trends in Brassica napus

The complexity of the epigenome landscape and transcriptional regulation is significantly increased during plant polyploidization, which drives genome evolution and contributes to the increased adaptability to diverse environments. However, a comprehensive epigenome map of Brassica napus is still unavailable. In this study, we performed integrative analysis of five histone modifications, RNA polymerase II occupancy, DNA methylation, and transcriptomes in two B. napus lines (2063A and B409), and established global maps of regulatory elements, chromatin states, and their dynamics for the whole genome (including the An and Cn subgenomes) in four tissue types (young leaf, flower bud, silique, and root) of these two lines. Approximately 65.8% of the genome was annotated with different epigenomic signals. Compared with the Cn subgenome, the An subgenome possesses a higher level of active epigenetic marks and lower level of repressive epigenetic marks. Genes from subgenome-unique regions contribute to the major differences between the An and Cn subgenomes. Asymmetric histone modifications between homeologous gene pairs reflect their biased expression patterns. We identified a novel bivalent chromatin state (with H3K4me1 and H3K27me3) in B. napus that is associated with tissue-specific gene expression. Furthermore, we observed that different types of duplicated genes have discrepant patterns of histone modification and DNA methylation levels. Collectively, our findings provide a valuable epigenetic resource for allopolyploid plants.

Chloroplast genomes of two Mediterranean Bupleurum species and the phylogenetic relationship inferred from combined analysis with East Asian species

The chloroplast genomes of Mediterranean Bupleurum species are reported for the first time. Phylogenetic analysis supports the species as a basal clade of Bupleurum with divergence time at 35.40 Ma. Bupleurum is one of the most species-rich genus with high medicinal value in Apiaceae. Although infrageneric classifications of Bupleurum have been the subject of numerous studies, it still remains controversial. Chloroplast genome information will prove essential in advancing our understanding on phylogenetic study. Here we report cp genomes of two woody Bupleurum species (Bupleurum gibraltaricum and B. fruticosum) endemic to Mediterranean. The complete cp genomes of the two species were 157,303 and 157,391 bp in size, respectively. They encoded 114 unique genes including 30 tRNA genes, 4 rRNA genes and 80 protein coding genes. Genome structure, distributions of SDRs and SSRs, gene content exhibited similarities among Bupleurum species. High variable hotspots were detected in eight intergenic spacers and four genes. Most of genes were under purifying selection with two exceptions: atpF and clpP. The phylogenetic analysis based on 80 coding genes revealed that the genus was divided into 2 distinct clades corresponding to the 2 subgenera (subg. Penninervia, subg. Bupleurum) with divergence time at the end of collision of India with Eurasia. Most species diversified mainly during the later period of uplift of Qinghai-Tibetan Plateau. The cp genomes of the two Bupleurum species can be significant complementary to insights into the cp genome characteristics of this genus. The comparative chloroplast genomes and phylogenetic analysis advances our understanding of the evolution of cp genomes and phylogeny in Bupleurum.

Genome-wide identification of PYL gene family in wheat: Evolution, expression and 3D structure analysis

Here, 38 wheat PYL genes (TaPYLs) belonging to 13 homoeologous groups were identified using the genome-search method, with 26 and 12 PYL genes identified in Triticum dicoccoides and Aegilops tauschii, respectively. Phylogenetic relationship, conserved domain and molecular evolution analysis revealed that PYL genes showed highly conservative between wheat and theprogenitors. Interaction network and miRNA target prediction found that TaPYLs could interact with the important components of ABA signaling pathway and Tae-miR966b-3p might be a hub regulator mediating wheat ABA signal network. Furthermore, the tissue-specific and stress-responsive TaPYLs were detected through RNA-seq analysis. Expressions of 10 TaPYLs were validated by QPCR analysis and the homoeologous genes showed significantly differential expression, suggesting subfunctionalization of them has occurred. Finally, 3D structures of the TaPYL proteins were predicted by homology modeling. This study lays the foundation for further functional study of PYL genes for development and stress tolerance improvement in wheat and beyond.

Chloroplast genome features of an important medicinal and edible plant: Houttuynia cordata (Saururaceae)

Houttuynia cordata (Saururaceae), an ancient and relic species, has been used as an important medicinal and edible plant in most parts of Asia. However, because of the lack of genome information and reliable molecular markers, studies on its population structure, or phylogenetic relationships with other related species are still rare. Here, we de novo assembled the complete chloroplast (cp) genome of H. cordata using the integration of the long PacBio and short Illumina reads. The cp genome of H. cordata showed a typical quadripartite cycle of 160,226 bp. This included a pair of inverted repeats (IRa and IRb) of 26,853 bp, separated by a large single-copy (LSC) region of 88,180 bp and a small single-copy (SSC) region of 18,340 bp. A total of 112 unique genes, including 79 protein-coding genes, 29 tRNA genes, and four rRNA genes, were identified in this cp genome. Eighty-one genes were located on the LSC region, 13 genes were located on the SSC region, and 17 two-copy genes were located on the IR region. Additionally, 48 repeat sequences and 86 SSR loci, which can be used as genomic markers for population structure analysis, were also detected. Phylogenetic analysis using 21 cp genomes of the Piperales family demonstrated that H. cordata had a close relationship with the species within the Aristolochia genus. Moreover, the results of mVISTA analysis and comparisons of IR regions demonstrated that the cp genome of H. cordata was conserved with that of the Aristolochia species. Our results provide valuable information for analyzing the genetic diversity and population structure of H. cordata, which can contribute to further its genetic improvement and breeding.

Comprehensive identification and analysis of DELLA genes throughout the plant kingdom

BACKGROUND

DELLAs play key roles in plant gibberellin signaling pathways and are generally important in plant development and growth. However, DELLAs in many plant taxa have not yet been systematically analyzed.

RESULTS

In our study, we searched for DELLA genes across 58 green plant genomes and found 181 DELLAs. Structure analysis showed some DELLA domains do not contain "D-E-L-L-A" sequences and instead contain similar domains, including DGLLA and DSLLH domains. "VHYNP" motifs in plant DELLAs comprise 23 types of sequences, while some DELLAs did not contain GRAS domains. In grape, we found that the DELLA protein GSVIVT01015465001 contains an F-box domain, while apple DELLA proteins MDP0000220512 and MDP0000403162 contain a WW domain and a BCIP domain, respectively. These DELLAs can be divided into 22 homologous groups and 17 orthologous groups, and 35 paralogous genes were identified. In total, 35 positively selected genes (PSGs) and 121 negatively selected genes (NSGs) were found among DELLAs based on selective pressure analysis, with an average K_s of NSGs that was significantly higher than that of PSGs (P < 0.05). Among the paralogous groups, CBI and Fop were significantly positively correlated with GC, GC1, GC2, GC12, and GC3, while CAI was significantly positively correlated with GC, GC1, GC12, and GC. The paralogous groups with ω values exceeding 1 had significantly higher K_a values. We also found some paralogous groups with ω values exceeding 1 that differed in their motifs.

CONCLUSIONS

This study provides helpful insights into the evolution of DELLA genes and offers exciting opportunities for the investigation of DELLA functions in different plants.

Identification, evolution and expression analyses of whole genome-wide TLP gene family in Brassica napus

BACKGROUND

Brassica is a very important genus of Brassicaceae, including many important oils, vegetables, forage crops, and ornamental horticultural plants. TLP family genes play important regulatory roles in the growth and development of plants. Therefore, this study used a bioinformatics approach to conduct the systematic comparative genomics analysis of TLP gene family in B. napus and other three important Brassicaceae crops.

RESULTS

Here, we identified a total of 29 TLP genes from B. napus genome, and they distributed on 16 chromosomes of B. napus. The evolutionary relationship showed that these genes could be divided into six groups from Group A to F. We found that the gene corresponding to Arabidopsis thaliana AT1G43640 was completely lost in B. rapa, B. oleracea and B. napus after whole genome triplication. The gene corresponding to AT1G25280 was retained in all the three species we analysed, belonging to 1:3:6 ratios. Our analyses suggested that there was a selective loss of some genes that might be redundant after genome duplication. This study proposed that the TLP genes in B. napus did not directly expansion compared with its diploid parents B. rapa, and B. oleracea. Instead, an indirect expansion of TLP gene family occurred in its two diploid parents. In addition, the study further utilized RNA-seq to detect the expression pattern of TLP genes between different tissues and two subgenomes.

CONCLUSIONS

This study systematically conducted the comparative analyses of TLP gene family in B. napus, discussed the loss and expansion of genes after genome duplication. It provided rich gene resources for exploring the molecular mechanism of TLP gene family. Meanwhile, it provided guidance and reference for the research of other gene families in B. napus.

The complete chloroplast genome sequence of yellow mustard (Sinapis alba L.) and its phylogenetic relationship to other Brassicaceae species

As a member of the large Brassicaceae family, yellow mustard (Sinapis alba L.) has been used as an important gene pool for the genetic improvement of cash crops in Brassicaceae. Understanding the phylogenetic relationship between Sinapis alba (S. alba) and other Brassicaceae crops can provide guidance on the introgression of its favorable alleles into related species. The chloroplast (cp) genome is an ideal model for assessing genome evolution and the phylogenetic relationships of complex angiosperm families. Herein, we de novo assembled the complete cp genome of S. alba by integrating the PacBio and Illumina sequencing platforms. A 153,760 bp quadripartite cycle without any gap was obtained, including a pair of inverted repeats (IRa and IRb) of 26,221 bp, separated by a large single copy (LSC) region of 83,506 bp and a small single copy (SSC) region of 17,821 bp. A total of 78 protein-coding genes, 30 tRNA genes, and four rRNA genes were identified in this cp genome, as were 89 simple sequence repeat (SSR) loci of 18 types. The codon usage analysis revealed a preferential use of the Leu codon with the A/U ending. The phylogenetic analysis using 82 Brassicaceae species demonstrated that S. alba had a close relationship with important Brassica and Raphanus species; moreover, it likely originated from a separate evolutionary pathway compared with the congeneric Sinapis arvensis. The synonymous (Ks) and non-synonymous (Ks) substitution rate analysis showed that genes encoding "Subunits of cytochrome b/f complex" were under the lowest purifying selection pressure, whereas those associated with "Maturase", "Subunit of acetyl-CoA", and "Subunits of NADH-dehydrogenase" underwent relatively higher purifying selection pressures. Our results provide valuable information for fully utilizing the S. alba cp genome as a potential genetic resource for the genetic improvement of Brassica and Raphanus species.

Systematic Comparisons of Positively Selected Genes between Gossypium arboreum and Gossypium raimondii Genomes

Background:

Studies of Positively Selected Genes (PSGs) in microorganisms and mammals have provided insights into the dynamics of genome evolution and the genetic basis of differences between species by using whole genome-wide scans. Systematic investigations and comparisons of PSGs in plants, however, are still limited.

Objective:

A systematic comparison of PSGs between the genomes of two cotton species, Gossypium arboreum (G. arboreum) and G. raimondii, will give the key answer for revealing molecular evolutionary differences in plants.

Methods:

Genome sequences of G. arboreum and G. raimondii were compared, including Whole Genome Duplication (WGD) events and genomic features such as gene number, gene length, codon bias index, evolutionary rate, number of expressed genes, and retention of duplicated copies.

Results:

Unlike the PSGs in G. raimondii, G. arboreum comprised more PSGs, smaller gene size and fewer expressed gene. In addition, the PSGs evolved at a higher rate of synonymous substitutions, but were subjected to lower selection pressure. The PSGs in G. arboreum were also retained with a lower number of duplicate gene copies than G. raimondii after a single WGD event involving Gossypium.

Conclusion:

These data indicate that PSGs in G. arboreum and G. raimondii differ not only in Ka/Ks, but also in their evolutionary, structural, and expression properties, indicating that divergence of G. arboreum and G. raimondii was associated with differences in PSGs in terms of evolutionary rates, gene length, expression patterns, and WGD retention in Gossypium.

Evolutionary Rate Heterogeneity and Functional Divergence of Orthologous Genes in Pyrus

Negatively selected genes (NSGs) and positively selected genes (PSGs) are the two types of most nuclear protein-coding genes in organisms. However, the evolutionary rates and characteristics of different types of genes have been rarely understood. In the present study, we investigate the rates of synonymous substitution (Ks) and the rates of non-synonymous substitution (Ka) by comparing the orthologous genes of two sequenced Pyrus species, Pyrus bretschneideri and Pyrus communis. Subsequently, we compared the evolutionary rates, gene structures, and expression profiles during different fruit development between PSGs and NSGs. Compared with the NSGs, the PSGs have fewer exons, shorter gene length, lower synonymous substitution rates and have higher evolutionary rates. Remarkably, gene expression patterns between two Pyrus species fruit indicated functional divergence for most of the orthologous genes derived from a common ancestor, and subfunctionalization for some of them. Overall, the present study shows that PSGs differs from NSGs not only under environmental selective pressure (Ka/Ks), but also in their structural, functional, and evolutionary properties. Additionally, our resulting data provides important insights for the evolution and highlights the diversification of orthologous genes in two Pyrus species.

Genome-Wide Identification and Expression Analysis of GA2ox, GA3ox, and GA20ox Are Related to Gibberellin Oxidase Genes in Grape (Vitis Vinifera L.)

Gibberellin (GAs) plays the important role in the regulation of grape developmental and growth processes. The bioinformatics analysis confirmed the differential expression of GA2, GA3, and GA20 gibberellin oxidase genes (VvGA2oxs, VvGA3oxs, and VvGA20oxs) in the grape genome, and laid a theoretical basis for exploring its role in grape. Based on the Arabidopsis GA2oxs, GA3oxs, and GA20oxs genes already reported, the VvGA2oxs, VvGA3oxs, and VvGA20oxs genes in the grape genome were identified using the BLAST software in the grape genome database. Bioinformatics analysis was performed using software such as DNAMAN v.5.0, Clustalx, MapGene2Chrom, MEME, GSDS v.2.0, ExPASy, DNAsp v.5.0, and MEGA v.7.0. Chip expression profiles were generated using grape Affymetrix GeneChip 16K and Grape eFP Browser gene chip data in PLEXdb. The expression of VvGA2oxs, VvGA3oxs, and VvGA20oxs gene families in stress was examined by qRT-PCR (Quantitative real-time-PCR). There are 24 GAoxs genes identified with the grape genome that can be classified into seven subgroups based on a phylogenetic tree, gene structures, and conserved Motifs in our research. The gene family has higher codon preference, while selectivity is negative selection of codon bias and selective stress was analyzed. The expression profiles indicated that the most of VvGAox genes were highly expressed under different time lengths of ABA (Abscisic Acid) treatment, NaCl, PEG and 5 °C. Tissue expression analysis showed that the expression levels of VvGA2oxs and VvGA20oxs in different tissues at different developmental stages of grapes were relatively higher than that of VvGA3oxs. Last but not least, qRT-PCR (Real-time fluorescent quantitative PCR) was used to determine the relative expression of the GAoxs gene family under the treatment of GA3 (gibberellin 3) and uniconazole, which can find that some VvGA2oxs was upregulated under GA3 treatment. Simultaneously, some VvGA3oxs and VvGA20oxs were upregulated under uniconazole treatment. In a nutshell, the GA2ox gene mainly functions to inactivate biologically active GAs, while GA20ox mainly degrades C20 gibberellins, and GA3ox is mainly composed of biologically active GAs. The comprehensive analysis of the three classes of VvGAoxs would provide a basis for understanding the evolution and function of the VvGAox gene family in a grape plant.

The complete chloroplast genome sequence of watercress (Nasturtium officinale R. Br.): Genome organization, adaptive evolution and phylogenetic relationships in Cardamineae

Watercress (Nasturtium officinale R. Br.), an aquatic leafy vegetable of the Brassicaceae family, is known as a nutritional powerhouse. Here, we de novo sequenced and assembled the complete chloroplast (cp) genome of watercress based on combined PacBio and Illumina data. The cp genome is 155,106 bp in length, exhibiting a typical quadripartite structure including a pair of inverted repeats (IRA and IRB) of 26,505 bp separated by a large single copy (LSC) region of 84,265 bp and a small single copy (SSC) region of 17,831 bp. The genome contained 113 unique genes, including 79 protein-coding genes, 30 tRNAs and 4 rRNAs, with 20 duplicate in the IRs. Compared with the prior cp genome of watercress deposited in GenBank, 21 single nucleotide polymorphisms (SNPs) and 27 indels were identified, mainly located in noncoding sequences. A total of 49 repeat structures and 71 simple sequence repeats (SSRs) were detected. Codon usage showed a bias for A/T-ending codons in the cp genome of watercress. Moreover, 45 RNA editing sites were predicted in 16 genes, all for C-to-U transitions. A comparative plastome study with Cardamineae species revealed a conserved gene order and high similarity of protein-coding sequences. Analysis of the Ka/Ks ratios of Cardamineae suggested positive selection exerted on the ycf2 gene in watercress, which might reflect specific adaptations of watercress to its particular living environment. Phylogenetic analyses based on complete cp genomes and common protein-coding genes from 56 species showed that the genus Nasturtium was a sister to Cardamine in the Cardamineae tribe. Our study provides valuable resources for future evolution, population genetics and molecular biology studies of watercress.

Genome-Wide Identification of Flowering-Time Genes in Brassica Species and Reveals a Correlation between Selective Pressure and Expression Patterns of Vernalization-Pathway Genes in Brassica napus

Flowering time is a key agronomic trait, directly influencing crop yield and quality. Many flowering-time genes have been identified and characterized in the model plant Arabidopsis thaliana; however, these genes remain uncharacterized in many agronomically important Brassica crops. In this study, we identified 1064, 510, and 524 putative orthologs of A. thaliana flowering-time genes from Brassica napus, Brassica rapa, and Brassica oleracea, respectively, and found that genes involved in the aging and ambient temperature pathways were fewer than those in other flowering pathways. Flowering-time genes were distributed mostly on chromosome C03 in B. napus and B. oleracea, and on chromosome A09 in B. rapa. Calculation of non-synonymous (Ka)/synonymous substitution (Ks) ratios suggested that flowering-time genes in vernalization pathways experienced higher selection pressure than those in other pathways. Expression analysis showed that most vernalization-pathway genes were expressed in flowering organs. Approximately 40% of these genes were highly expressed in the anther, whereas flowering-time integrator genes were expressed in a highly organ-specific manner. Evolutionary selection pressures were negatively correlated with the breadth and expression levels of vernalization-pathway genes. These findings provide an integrated framework of flowering-time genes in these three Brassica crops and provide a foundation for deciphering the relationship between gene expression patterns and their evolutionary selection pressures in Brassica napus.

Comparative analysis of selection mode reveals different evolutionary rate and expression pattern in Arachis duranensis and Arachis ipaënsis duplicated genes

Our results reveal that K_s is a determining factor affecting selective pressure and different evolution and expression patterns are detected between PSGs and NSGs in wild Arachis duplicates. Selective pressure, including purifying (negative) and positive selection, can be detected in organisms. However, studies on comparative evolutionary rates, gene expression patterns and gene features between negatively selected genes (NSGs) and positively selected genes (PSGs) are lagging in paralogs of plants. Arachis duranensis and Arachis ipaënsis are ancestors of the cultivated peanut, an important oil and protein crop. Here, we carried out a series of systematic analyses, comparing NSG and PSG in paralogs, using genome sequences and transcriptome datasets in A. duranensis and A. ipaënsis. We found that synonymous substitution rate (K_s) is a determining factor affecting selective pressure in A. duranensis and A. ipaënsis duplicated genes. Lower expression level, lower gene expression breadth, higher codon bias and shorter polypeptide length were found in PSGs and not in NSGs. The correlation analyses showed that gene expression breadth was positively correlated with polypeptide length and GC content at the first codon site (GC1) in PSGs and NSGs, respectively. There was a negative correlation between expression level and polypeptide length in PSGs. In NSGs, the K_s was positively correlated with expression level, gene expression breadth, GC1, and GC content at the third codon site (GC3), but selective pressure was negatively correlated with expression level, gene expression breadth, polypeptide length, GC1, and GC3 content. The function of most duplicated gene pairs was divergent under drought and nematode stress. Taken together, our results show that different evolution and expression patterns occur between PSGs and NSGs in paralogs of two wild Arachis species.

WRKY transcription factors in legumes

BACKGROUND

WRKY transcription factors, so named because of the WRKYGQK heptapeptide at the N-terminal end, are widely distributed in plants and play an important role in physiological changes and response to biotic and abiotic stressors. Many previous studies have focused on the evolution of WRKY transcription factors in a given plant; however, little is known about WRKY evolution in legumes. The gene expression pattern of duplicated WRKY transcription factors remains unclear.

RESULTS

We first identified the WRKY proteins in 12 legumes. We found that the WRKYGQK heptapeptide tended to mutate into WRKYGKK. The Q site in WRKYGQK preferentially mutated, while W, K, and Y were conserved. The phylogenetic tree shows that the WRKY proteins in legumes have multiple origins, especially group IIc. For example, WRKY64 from Lupinus angustifolius (LaWRKY64) contains three WRKY domains, of which the first two clustered together in the N-terminal WRKY domain of the group I WRKY protein, and the third WRKY domain grouped in the C-terminal WRKY domain of the group I WRKY protein. Orthologous WRKY genes have a faster evolutionary rate and are subject to constrained selective pressure, unlike paralogous WRKY genes. Different gene features were observed between duplicated WRKY genes and singleton WRKY genes. Duplicated Glycine max WRKY genes with similar gene features have gene expression divergence.

CONCLUSIONS

We analyzed the WRKY number and type in 12 legumes, concluding that the WRKY proteins have multiple origins. A novel WRKY protein, LaWRKY64, was found in L. angustifolius. The first two WRKY domains of LaWRKY64 have the same origin. The orthologous and paralogous WRKY proteins have different evolutionary rates. Duplicated WRKY genes have gene expression divergence under normal growth conditions in G. max. These results provide insight into understanding WRKY evolution and expression.

A segmental duplication in the common ancestor of Brassicaceae is responsible for the origin of the paralogs KCS6-KCS5, which are not shared with other angiosperms

Novel morphological structures allowed adaptation to dry conditions in early land plants. The cuticle, one such novelty, plays diverse roles in tolerance to abiotic and biotic stresses and plant development. Cuticular waxes represent a major constituent of the cuticle and are comprised of an assortment of chemicals that include, among others, very long chain fatty acids (VLCFAs). Members of the β-ketoacyl coenzyme A synthases (KCS) gene family code for enzymes that are essential for fatty acid biosynthesis. The gene KCS6 (CUT1) is known to be a key player in the production of VLCFA precursors essential for the synthesis of cuticular waxes in the model plant Arabidopsis thaliana (Brassicaceae). Despite its functional importance, relatively little is known about the evolutionary history of KCS6 or its paralog KCS5 in Brassicaceae or beyond. This lacuna becomes important when we extrapolate understanding of mechanisms gained from the model plant to its containing clades Brassicaceae, flowering plants, or beyond. The Brassicaceae, with several sequenced genomes and a known history of paleoploidy, mesopolyploidy and neopolyploidy, offer a system in which to study the evolution and diversification of the KCS6-KCS5 paralogy. Our phylogenetic analyses across green plants, combined with comparative genomic, microsynteny and evolutionary rates analyses across nine genomes of Brassicaceae, reveal that (1) the KCS6-KCS5 paralogy arose as the result of a large segmental duplication in the ancestral Brassicaceae, (2) the KCS6-KCS5 lineage is represented by a single copy in other flowering plant lineages, (3) the duplicated segments undergo different degrees of retention and loss, and (4) most of the genes in the KCS6 and KCS5 gene blocks (including KCS6 and KCS5 themselves) are under purifying selection. The last also true for most members of the KCS gene family in Brassicaceae, except for KCS8, KCS9 and KCS17, which are under positive selection and may be undergoing functional evolution, meriting further investigation. Overall, our results clearly establish that the ancestral KCS6/5 gene duplicated in the Brassicaceae lineage. It is possible that any specialized functions of KCS5 found in Brassicaceae are either part of a set of KCS6/5 gene functions in the rest of the flowering plants, or unique to Brassicaceae.

Genome-Wide Characterization of bHLH Genes in Grape and Analysis of their Potential Relevance to Abiotic Stress Tolerance and Secondary Metabolite Biosynthesis

Basic helix-loop-helix (bHLH) transcription factors are involved in many abiotic stress responses as well as flavonol and anthocyanin biosynthesis. In grapes (Vitis vinifera L.), flavonols including anthocyanins and condensed tannins are most abundant in the skins of the berries. Flavonols are important phytochemicals for viticulture and enology, but grape bHLH genes have rarely been examined. We identified 94 grape bHLH genes in a genome-wide analysis and performed Nr and GO function analyses for these genes. Phylogenetic analyses placed the genes into 15 clades, with some remaining orphans. 41 duplicate gene pairs were found in the grape bHLH gene family, and all of these duplicate gene pairs underwent purifying selection. Nine triplicate gene groups were found in the grape bHLH gene family and all of these triplicate gene groups underwent purifying selection. Twenty-two grape bHLH genes could be induced by PEG treatment and 17 grape bHLH genes could be induced by cold stress treatment including a homologous form of MYC2, VvbHLH007. Based on the GO or Nr function annotations, we found three other genes that are potentially related to anthocyanin or flavonol biosynthesis: VvbHLH003, VvbHLH007, and VvbHLH010. We also performed a cis-acting regulatory element analysis on some genes involved in flavonoid or anthocyanin biosynthesis and our results showed that most of these gene promoters contained G-box or E-box elements that could be recognized by bHLH family members.