Genome-wide identification, classification, evolutionary analysis and gene expression patterns of the protein kinase gene family in wheat and Aegilops tauschii

Genome-wide evolutionary analysis of TKL_CTR1-DRK-2 gene family and functional characterization reveals that TaCTR1 positively regulates flowering time in wheat

BACKGROUND

Flowering time has an important effect on regional adaptation and yields for crops. The tyrosine kinase-like (TKL) gene family is widely existed and participates in many biological processes in plants. Furthermore, only few TKLs have been characterized functions in controlling flowering time in wheat.

RESULTS

Here, we report that TaCTR1, a tyrosine kinase-like (TKL) gene, regulates flowering time in wheat. Based on identification and evolutionary analysis of TKL_CTR1-DRK-2 subfamily in 15 plants, we proposed an evolutionary model for TaCTR1, suggesting that occurrence of some exon fusion events during evolution. The overexpression of TaCTR1 caused early flowering time in transgenic lines. Transcriptomics analysis enabled identification of mass differential expression genes including plant hormone (ET, ABA, IAA, BR) signaling, flavonoid biosynthesis, phenolamides and antioxidant, and flowering-related genes in TaCTR1 overexpression transgenic lines compared with WT plants. qRT-PCR results showed that the expression levels of ethylene (ET) signal-related genes (ETR, EIN, ERF) and flowering-related genes (FT, PPD1, CO, PRR, PHY) were altered in TaCTR1-overexpressing wheat compared with WT plants. Metabonomics analysis showed that flavonoid contents were altered.

CONCLUSIONS

Thus, the results show that TaCTR1 plays a positive role in controlling flowering time by activating various signaling pathways and regulating flowering-related genes, and will provide new insights on the mechanisms of wheat flowering regulation.

Phylogeny of the plant receptor-like kinase (RLK) gene family and expression analysis of wheat RLK genes in response to biotic and abiotic stresses

BACKGROUND

The receptor-like kinase (RLK) gene families in plants contains a large number of members. They are membrane proteins with an extracellular receptor domain and participate in biotic and abiotic stress responses.

RESULTS

In this study, we identified RLKs in 15 representative plant genomes, including wheat, and classified them into 64 subfamilies by using four types of phylogenetic trees and HMM models. Conserved exon‒intron structures with conserved exon phases in the kinase domain were found in many RLK subfamilies from Physcomitrella patens to Triticum aestivum. Domain distributions of RLKs were also diagrammed. Collinearity events and tandem gene clusters suggested that polyploidization and tandem duplication events contributed to the member expansions of T. aestivum RLKs. Global expression pattern analysis was performed by using public transcriptome data. These analyses were involved in T. aestivum, Aegilops tauschii and Brachypodium distachyon RLKs under biotic and abiotic stresses. We also selected 9 RLKs to validate the transcriptome prediction by using qRT‒PCR under drought treatment and with Fusarium graminearum infection. The expression trends of these 9 wheat RLKs from public transcriptome data were consistent with the results of qRT‒PCR, indicating that they might be stress response genes under drought or F. graminearum treatments.

CONCLUSION

In this study, we identified, classified, evolved, and expressed RLKs in wheat and related plants. Thus, our results will provide insights into the evolutionary history and molecular mechanisms of wheat RLKs.

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

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

The rubber tree kinome: Genome-wide characterization and insights into coexpression patterns associated with abiotic stress responses

The protein kinase (PK) superfamily constitutes one of the largest and most conserved protein families in eukaryotic genomes, comprising core components of signaling pathways in cell regulation. Despite its remarkable relevance, only a few kinase families have been studied in Hevea brasiliensis. A comprehensive characterization and global expression analysis of the PK superfamily, however, is currently lacking. In this study, with the aim of providing novel inferences about the mechanisms associated with the stress response developed by PKs and retained throughout evolution, we identified and characterized the entire set of PKs, also known as the kinome, present in the Hevea genome. Different RNA-sequencing datasets were employed to identify tissue-specific expression patterns and potential correspondences between different rubber tree genotypes. In addition, coexpression networks under several abiotic stress conditions, such as cold, drought and latex overexploitation, were employed to elucidate associations between families and tissues/stresses. A total of 1,809 PK genes were identified using the current reference genome assembly at the scaffold level, and 1,379 PK genes were identified using the latest chromosome-level assembly and combined into a single set of 2,842 PKs. These proteins were further classified into 20 different groups and 122 families, exhibiting high compositional similarities among family members and with two phylogenetically close species Manihot esculenta and Ricinus communis. Through the joint investigation of tandemly duplicated kinases, transposable elements, gene expression patterns, and coexpression events, we provided insights into the understanding of the cell regulation mechanisms in response to several conditions, which can often lead to a significant reduction in rubber yield.

Genome-wide evolutionary analysis of AUX/IAA gene family in wheat identifies a novel gene TaIAA15-1A regulating flowering time by interacting with ARF

Flowering time is a critical agronomic trait that has strong effects on crop yields. Auxin signaling pathway plays an important role in various development processes, such as flowering, grain development. However, no Aux/IAA gene had been reported to have functions involving in wheat flowering time. Here, we systematically performed genome-wide identification, classification, domain distribution, exon-intron structure, chromosome locations and global expression pattern of Aux/IAA gene family in 14 plant genomes (including Triticum aestivum). A phylogenetic model was proposed to infer the Aux/IAA evolutionary history involving in a central exon-intron structure "2121" during evolution. Overexpression of TaIAA15-1A caused an early flowering time in Brachypodium. RNA-seq analysis showed that TaIAA15-1A overexpression alters various pathways including phytohormone signaling pathway, flowering-related pathway, and polyamine biosynthesis pathway. Screening of auxin response factor (ARF) genes identified BdARF16 that interacted with TaIAA15-1A. Exogenous polyamine (spermidine and spermine) treatments promoted early flowering and (putrescine and DCHA) delayed flowering time of WT plants. Our finding will provide insights on mechanisms of Aux/IAAs gene family and TaIAA15-1A, illustrating the potential during crop improvement programs.

Microspore embryogenesis induction by mannitol and TSA results in a complex regulation of epigenetic dynamics and gene expression in bread wheat

Reprogramming of microspores development towards embryogenesis mediated by stress treatment constitutes the basis of doubled haploid production. Recently, compounds that alter histone post-translational modifications (PTMs) have been reported to enhance microspore embryogenesis (ME), by altering histones acetylation or methylation. However, epigenetic mechanisms underlying ME induction efficiency are poorly understood. In this study, the epigenetic dynamics and the expression of genes associated with histone PTMs and ME induction were studied in two bread wheat cultivars with different ME response. Microspores isolated at 0, 3 and 5 days, treated with 0.7M mannitol (MAN) and 0.7M mannitol plus 0.4µM trichostatin A (TSA), which induced ME more efficiently, were analyzed. An additional control of gametophytic development was included. Microspores epigenetic state at the onset of ME induction was distinctive between cultivars by the ratio of H3 variants and their acetylated forms, the localization and percentage of labeled microspores with H3K9ac, H4K5ac, H4K16ac, H3K9me2 and H3K27me3, and the expression of genes related to pollen development. These results indicated that microspores of the high responding cultivar could be at a less advanced stage in pollen development. MAN and TSA resulted in a hyperacetylation of H3.2, with a greater effect of TSA. Histone PTMs were differentially affected by both treatments, with acetylation being most concerned. The effect of TSA was observed in the H4K5ac localization pattern at 3dT in the mid-low responding cultivar. Three gene networks linked to ME response were identified. TaHDT1, TaHAG2, TaYAO, TaNFD6-A, TabZIPF1 and TaAGO802-B, associated with pollen development, were down-regulated. TaHDA15, TaHAG3, TaHAM, TaYUC11D, Ta-2B-LBD16 TaMS1 and TaDRM3 constituted a network implicated in morphological changes by auxin signaling and cell wall modification up-regulated at 3dT. The last network included TaHDA18, TaHAC1, TaHAC4, TaABI5, TaATG18fD, TaSDG1a-7A and was related to ABA and ethylene hormone signaling pathways, DNA methylation and autophagy processes, reaching the highest expression at 5dT. The results indicated that TSA mainly modified the regulation of genes related to pollen and auxin signaling. This study represents a breakthrough in identifying the epigenetic dynamics and the molecular mechanisms governing ME induction efficiency, with relevance to recalcitrant wheat genotypes and other crops.

Silveira S, Brommonschenkel SH, Fontes EPB. RLPredictiOme, a Machine Learning-Derived Method for High-Throughput Prediction of Plant Receptor-like Proteins, Reveals Novel Classes of Transmembrane Receptors

Cell surface receptors play essential roles in perceiving and processing external and internal signals at the cell surface of plants and animals. The receptor-like protein kinases (RLK) and receptor-like proteins (RLPs), two major classes of proteins with membrane receptor configuration, play a crucial role in plant development and disease defense. Although RLPs and RLKs share a similar single-pass transmembrane configuration, RLPs harbor short divergent C-terminal regions instead of the conserved kinase domain of RLKs. This RLP receptor structural design precludes sequence comparison algorithms from being used for high-throughput predictions of the RLP family in plant genomes, as has been extensively performed for RLK superfamily predictions. Here, we developed the RLPredictiOme, implemented with machine learning models in combination with Bayesian inference, capable of predicting RLP subfamilies in plant genomes. The ML models were simultaneously trained using six types of features, along with three stages to distinguish RLPs from non-RLPs (NRLPs), RLPs from RLKs, and classify new subfamilies of RLPs in plants. The ML models achieved high accuracy, precision, sensitivity, and specificity for predicting RLPs with relatively high probability ranging from 0.79 to 0.99. The prediction of the method was assessed with three datasets, two of which contained leucine-rich repeats (LRR)-RLPs from Arabidopsis and rice, and the last one consisted of the complete set of previously described Arabidopsis RLPs. In these validation tests, more than 90% of known RLPs were correctly predicted via RLPredictiOme. In addition to predicting previously characterized RLPs, RLPredictiOme uncovered new RLP subfamilies in the Arabidopsis genome. These include probable lipid transfer (PLT)-RLP, plastocyanin-like-RLP, ring finger-RLP, glycosyl-hydrolase-RLP, and glycerophosphoryldiester phosphodiesterase (GDPD, GDPDL)-RLP subfamilies, yet to be characterized. Compared to the only Arabidopsis GDPDL-RLK, molecular evolution studies confirmed that the ectodomain of GDPDL-RLPs might have undergone a purifying selection with a predominance of synonymous substitutions. Expression analyses revealed that predicted GDPGL-RLPs display a basal expression level and respond to developmental and biotic signals. The results of these biological assays indicate that these subfamily members have maintained functional domains during evolution and may play relevant roles in development and plant defense. Therefore, RLPredictiOme provides a framework for genome-wide surveys of the RLP superfamily as a foundation to rationalize functional studies of surface receptors and their relationships with different biological processes.

Identification of key genes associated with secondary metabolites biosynthesis by system network analysis in Valeriana officinalis

Valeriana officinalis is a medicinal plant, a source of bioactive chemical compounds and secondary metabolites which are applied in pharmaceutical industries. The advent of ethnomedicine has provided alternatives for disease treatment and has increased demands for natural products and bioactive compounds. A set of preliminary steps to answers for such demands can include integrative omics for systems metabolic engineering, as an approach that contributes to the understanding of cellular metabolic status. There is a growing trend of this approach for genetically engineering metabolic pathways in plant systems, by which natural and synthetic compounds can be produced. As in the case of most medicinal plants, there are no sufficient information about molecular mechanisms involved in the regulation of metabolic pathways in V. officinalis. In this research, systems biology was performed on the RNA-seq transcriptome and metabolome data to find key genes that contribute to the synthesis of major secondary metabolites in V. officinalis. The R Package Weighted Gene Co-Expression Network Analysis (WGCNA) was employed to analyze the data. Based on the results, some major modules and hub genes were identified to be associated with the valuable secondary metabolites. In addition, some TF-encoding genes, including AP2/ERF-ERF, WRKY and NAC TF families, as well as some regulatory factors including protein kinases and transporters were identified. The results showed that several novel hub genes, such as PCMP-H24, RPS24B, ANX1 and PXL1, may play crucial roles in metabolic pathways. The current findings provide an overall insight into the metabolic pathways of V. officinalis and can expand the potential for engineering genome-scale pathways and systems metabolic engineering to increase the production of bioactive compounds by plants.

The Wild Sugarcane and Sorghum Kinomes: Insights Into Expansion, Diversification, and Expression Patterns

The protein kinase (PK) superfamily is one of the largest superfamilies in plants and the core regulator of cellular signaling. Despite this substantial importance, the kinomes of sugarcane and sorghum have not been profiled. Here, we identified and profiled the complete kinomes of the polyploid Saccharum spontaneum (Ssp) and Sorghum bicolor (Sbi), a close diploid relative. The Sbi kinome was composed of 1,210 PKs; for Ssp, we identified 2,919 PKs when disregarding duplications and allelic copies, and these were related to 1,345 representative gene models. The Ssp and Sbi PKs were grouped into 20 groups and 120 subfamilies and exhibited high compositional similarities and evolutionary divergences. By utilizing the collinearity between the species, this study offers insights into Sbi and Ssp speciation, PK differentiation and selection. We assessed the PK subfamily expression profiles via RNA-Seq and identified significant similarities between Sbi and Ssp. Moreover, coexpression networks allowed inference of a core structure of kinase interactions with specific key elements. This study provides the first categorization of the allelic specificity of a kinome and offers a wide reservoir of molecular and genetic information, thereby enhancing the understanding of Sbi and Ssp PK evolutionary history.

Genome-wide identification and expression analysis of YTH domain-containing RNA-binding protein family in common wheat

BACKGROUND

N6-Methyladenosine (m6A) is the most widespread RNA modification that plays roles in the regulation of genes and genome stability. YT521-B homology (YTH) domain-containing RNA-binding proteins are important RNA binding proteins that affect the fate of m6A-containing RNA by binding m6A. Little is known about the YTH genes in common wheat (Triticum aestivum L.), one of the most important crops for humans.

RESULTS

A total of 39 TaYTH genes were identified in common wheat, which are comprised of 13 homologous triads, and could be mapped in 18 out of the 21 chromosomes. A phylogenetic analysis revealed that the TaYTHs could be divided into two groups: YTHDF (TaDF) and YTHDC (TaDC). The TaYTHs in the same group share similar motif distributions and domain organizations, which indicates functional similarity between the closely related TaYTHs. The TaDF proteins share only one domain, which is the YTH domain. In contrast, the TaDCs possess three C3H1-type zinc finger repeats at their N-termini in addition to their central YTH domain. In TaDFs, the predicated aromatic cage pocket that binds the methylysine residue of m6A is composed of tryptophan, tryptophan, and tryptophan (WWW). In contrast, the aromatic cage pocket in the TaDCs is composed of tryptophan, tryptophan, and tyrosine (WWY). In addition to the general aspartic acid or asparagine residue used to form a hydrogen bond with N¹ of m6A, histidine might be utilized in some TaDFb proteins. An analysis of the expression using both online RNA-Seq data and quantitative real-time PCR verification revealed that the TaDFa and TaDFb genes are highly expressed in various tissues/organs compared with that of TaDFcs and TaDCs. In addition, the expression of the TaYTH genes is changed in response to various abiotic stresses.

CONCLUSIONS

In this study, we identified 39 TaYTH genes from common wheat. The phylogenetic structure, chromosome distribution, and patterns of expression of these genes and their protein structures were analyzed. Our results provide a foundation for the functional analysis of TaYTHs in the future.

Marker-Trait Associations for Enhancing Agronomic Performance, Disease Resistance, and Grain Quality in Synthetic and Bread Wheat Accessions in Western Siberia

Exploiting genetically diverse lines to identify genes for improving crop performance is needed to ensure global food security. A genome-wide association study (GWAS) was conducted using 46,268 SNP markers on a diverse panel of 143 hexaploid bread and synthetic wheat to identify potential genes/genomic regions controlling agronomic performance (yield and 26 yield-related traits), disease resistance, and grain quality traits. From phenotypic evaluation, we found large genetic variation among the 35 traits and recommended five lines having a high yield, better quality, and multiple disease resistance for direct use in a breeding program. From a GWAS, we identified a total of 243 significant marker-trait associations (MTAs) for 35 traits that explained up to 25% of the phenotypic variance. Of these, 120 MTAs have not been reported in the literature and are potentially novel MTAs. In silico gene annotation analysis identified 116 MTAs within genes and of which, 21 MTAs were annotated as a missense variant. Furthermore, we were able to identify 23 co-located multi-trait MTAs that were also phenotypically correlated to each other, showing the possibility of simultaneous improvement of these traits. Additionally, most of the co-located MTAs were within genes. We have provided genomic fingerprinting for significant markers with favorable and unfavorable alleles in the diverse set of lines for developing elite breeding lines from useful trait-integration. The results from this study provided a further understanding of genetically complex traits and would facilitate the use of diverse wheat accessions for improving multiple traits in an elite wheat breeding program.

Genome-wide and evolutionary analysis of the class III peroxidase gene family in wheat and Aegilops tauschii reveals that some members are involved in stress responses

Background

The class III peroxidase (PRX) gene family is a plant-specific member of the PRX superfamily that is closely related to various physiological processes, such as cell wall loosening, lignification, and abiotic and biotic stress responses. However, its classification, evolutionary history and gene expression patterns are unclear in wheat and Aegilops tauschii.

Results

Here, we identified 374, 159 and 169 PRXs in Triticum aestivum, Triticum urartu and Ae. tauschii, respectively. Together with PRXs detected from eight other plants, they were classified into 18 subfamilies. Among subfamilies V to XVIII, a conserved exon-intron structure within the “001” exon phases was detected in the PRX domain. Based on the analysis, we proposed a phylogenetic model to infer the evolutionary history of the exon-intron structures of PRX subfamilies. A comparative genomics analysis showed that subfamily VII could be the ancient subfamily that originated from green algae (Chlamydomonas reinhardtii). Further integrated analysis of chromosome locations and collinearity events of PRX genes suggested that both whole genome duplication (WGD) and tandem duplication (TD) events contributed to the expansion of T. aestivum PRXs (TaePRXs) during wheat evolution. To validate functions of these genes in the regulation of various physiological processes, the expression patterns of PRXs in different tissues and under various stresses were studied using public microarray datasets. The results suggested that there were distinct expression patterns among different tissues and PRXs could be involved in biotic and abiotic responses in wheat. qRT-PCR was performed on samples exposed to drought, phytohormone treatments and Fusarium graminearum infection to validate the microarray predictions. The predicted subcellular localizations of some TaePRXs were consistent with the confocal microscopy results. We predicted that some TaePRXs had hormone-responsive cis-elements in their promoter regions and validated these predicted cis-acting elements by sequencing promoters.

Conclusion

In this study, identification, classification, evolution, and expression patterns of PRXs in wheat and relative plants were performed. Our results will provide information for further studies on the evolution and molecular mechanisms of wheat PRXs.

Electronic supplementary material

The online version of this article (10.1186/s12864-019-6006-5) contains supplementary material, which is available to authorized users.

Genome-Wide Association Study Reveals Novel Genomic Regions for Grain Yield and Yield-Related Traits in Drought-Stressed Synthetic Hexaploid Wheat

Synthetic hexaploid wheat (SHW; 2n = 6x = 42, AABBDD, Triticum aestivum L.) is produced from an interspecific cross between durum wheat (2n = 4x = 28, AABB, T. turgidum L.) and goat grass (2n = 2x = 14, DD, Aegilops tauschii Coss.) and is reported to have significant novel alleles-controlling biotic and abiotic stresses resistance. A genome-wide association study (GWAS) was conducted to unravel these loci [marker⁻trait associations (MTAs)] using 35,648 genotyping-by-sequencing-derived single nucleotide polymorphisms in 123 SHWs. We identified 90 novel MTAs (45, 11, and 34 on the A, B, and D genomes, respectively) and haplotype blocks associated with grain yield and yield-related traits including root traits under drought stress. The phenotypic variance explained by the MTAs ranged from 1.1% to 32.3%. Most of the MTAs (120 out of 194) identified were found in genes, and of these 45 MTAs were in genes annotated as having a potential role in drought stress. This result provides further evidence for the reliability of MTAs identified. The large number of MTAs (53) identified especially on the D-genome demonstrate the potential of SHWs for elucidating the genetic architecture of complex traits and provide an opportunity for further improvement of wheat under rapidly changing climatic conditions.

Variation and diversity of the breakpoint sequences on 4AL for the 4AL/5AL translocation in Triticum

The translocation of 4AL/5AL in Triticum, which occurred before the differentiation of T. urartu and einkorn, is an important chromosomal rearrangement. Recently, the first identification of breakpoint sequence on 4AL for this translocation provides the opportunity to analyze the variation and diversity of breakpoints in Triticum. In this study, the breakpoint regions of 52 accessions from 21 species were isolated and further characterized. The sequences were divided into 12 types based on their lengths, which ranged from 2009 to 2552 bp. Cluster analysis showed that they were further divided into three groups. Interesting evolutionary relationships among a few of the species were observed and discussed. Multiple sequence alignment of the 52 sequences made it possible to detect 13 insertion and deletion length polymorphisms (InDels) and 101 single nucleotide polymorphisms (SNPs). Furthermore, several species- or accession-specific SNPs or InDels were also identified. Based on BLAST analysis of the conserved sequences, the breakpoint was narrowed down to a 125 bp fragment. Taken together, the results obtained in this study enrich our understanding of chromosomal breakpoints and will be useful for the identification of other breakpoints in wheat.

Genome-wide classification, evolutionary analysis and gene expression patterns of the kinome in Gossypium

The protein kinase (PK, kinome) family is one of the largest families in plants and regulates almost all aspects of plant processes, including plant development and stress responses. Despite their important functions, comprehensive functional classification, evolutionary analysis and expression patterns of the cotton PK gene family has yet to be performed on PK genes. In this study, we identified the cotton kinomes in the Gossypium raimondii, Gossypium arboretum, Gossypium hirsutum and Gossypium barbadense genomes and classified them into 7 groups and 122-24 subfamilies using software HMMER v3.0 scanning and neighbor-joining (NJ) phylogenetic analysis. Some conserved exon-intron structures were identified not only in cotton species but also in primitive plants, ferns and moss, suggesting the significant function and ancient origination of these PK genes. Collinearity analysis revealed that 16.6 million years ago (Mya) cotton-specific whole genome duplication (WGD) events may have played a partial role in the expansion of the cotton kinomes, whereas tandem duplication (TD) events mainly contributed to the expansion of the cotton RLK group. Synteny analysis revealed that tetraploidization of G. hirsutum and G. barbadense contributed to the expansion of G. hirsutum and G. barbadense PKs. Global expression analysis of cotton PKs revealed stress-specific and fiber development-related expression patterns, suggesting that many cotton PKs might be involved in the regulation of the stress response and fiber development processes. This study provides foundational information for further studies on the evolution and molecular function of cotton PKs.