Genomic Dissection and Expression Profiling Revealed Functional Divergence in Triticum aestivum Leucine Rich Repeat Receptor Like Kinases (TaLRRKs)

Comprehensive analysis of sesame LRR-RLKs: structure, evolution and dynamic expression profiles under Macrophomina phaseolina stress

Leucine-rich repeat receptor-like kinases (LRR-RLKs) can participate in the regulation of plant growth and development, immunity and signal transduction. Sesamum indicum, one of the most important oil crops, has a significant role in promoting human health. In this study, 175 SiLRR-RLK genes were identified in S. indicum, and they were subdivided into 12 subfamilies by phylogenetic analysis. Gene duplication analysis showed that the expansion of the SiLRR-RLK family members in the sesame was mainly due to segmental duplication. Moreover, the gene expansion of subfamilies IV and III contributed to the perception of stimuli under M. phaseolina stress in the sesame. The collinearity analysis with other plant species revealed that the duplication of SiLRR-RLK genes occurred after the differentiation of dicotyledons and monocotyledons. The expression profile analysis and functional annotation of SiLRR-RLK genes indicated that they play a vital role in biotic stress. Furthermore, the protein-protein interaction and coexpression networks suggested that SiLRR-RLKs contributed to sesame resistance to Macrophomina phaseolina by acting alone or as a polymer with other SiLRR-RLKs. In conclusion, the comprehensive analysis of the SiLRR-RLK gene family provided a framework for further functional studies on SiLRR-RLK genes.

In silico characterization of five novel disease-resistance proteins in Oryza sativa sp. japonica against bacterial leaf blight and rice blast diseases

In the current study, gene network analysis revealed five novel disease-resistance proteins against bacterial leaf blight (BB) and rice blast (RB) diseases caused by Xanthomonas oryzae pv. oryzae (Xoo) and Magnaporthe oryzae (M. oryzae), respectively. In silico modeling, refinement, and model quality assessment were performed to predict the best structures of these five proteins and submitted to ModelArchive for future use. An in-silico annotation indicated that the five proteins functioned in signal transduction pathways as kinases, phospholipases, transcription factors, and DNA-modifying enzymes. The proteins were localized in the nucleus and plasma membrane. Phylogenetic analysis showed the evolutionary relation of the five proteins with disease-resistance proteins (XA21, OsTRX1, PLD, and HKD-motif-containing proteins). This indicates similar disease-resistant properties between five unknown proteins and their evolutionary-related proteins. Furthermore, gene expression profiling of these proteins using public microarray data showed their differential expression under Xoo and M. oryzae infection. This study provides an insight into developing disease-resistant rice varieties by predicting novel candidate resistance proteins, which will assist rice breeders in improving crop yield to address future food security through molecular breeding and biotechnology.

Supplementary Information

The online version contains supplementary material available at 10.1007/s13205-023-03893-5.

An update on evolutionary, structural, and functional studies of receptor-like kinases in plants

All living organisms must develop mechanisms to cope with and adapt to new environments. The transition of plants from aquatic to terrestrial environment provided new opportunities for them to exploit additional resources but made them vulnerable to harsh and ever-changing conditions. As such, the transmembrane receptor-like kinases (RLKs) have been extensively duplicated and expanded in land plants, increasing the number of RLKs in the advanced angiosperms, thus becoming one of the largest protein families in eukaryotes. The basic structure of the RLKs consists of a variable extracellular domain (ECD), a transmembrane domain (TM), and a conserved kinase domain (KD). Their variable ECDs can perceive various kinds of ligands that activate the conserved KD through a series of auto- and trans-phosphorylation events, allowing the KDs to keep the conserved kinase activities as a molecular switch that stabilizes their intracellular signaling cascades, possibly maintaining cellular homeostasis as their advantages in different environmental conditions. The RLK signaling mechanisms may require a coreceptor and other interactors, which ultimately leads to the control of various functions of growth and development, fertilization, and immunity. Therefore, the identification of new signaling mechanisms might offer a unique insight into the regulatory mechanism of RLKs in plant development and adaptations. Here, we give an overview update of recent advances in RLKs and their signaling mechanisms.

Genome-wide identification of the LRR-RLK gene family in peanut and functional characterization of AhLRR-RLK265 in salt and drought stresses

Leucine-rich repeat receptor-like kinases (LRR-RLKs) play important roles in plant developmental regulations and various stress responses. Peanut (Arachis hypogaea L.) is a worldwide important oil crop; however, no systematic identification or analysis of the peanut LRR-RLK gene family has been reported. In present study, 495 LRR-RLK genes in peanut were identified and analyzed. The 495 AhLRR-RLK genes were classed into 14 groups and 10 subgroups together with their Arabidopsis homologs according to phylogenetic analyses, and 491 of 495 AhLRR-RLK genes unequally located on 20 chromosomes. Analyses of gene structure and protein motif organization revealed similarity in exon/intron and motif organization among members of the same subgroup, further supporting the phylogenetic results. Gene duplication events were found in peanut LRR-RLK gene family via syntenic analysis, which were important in LRR-RLK gene family expansion in peanut. We found that the expression of AhLRR-RLK genes was detected in different tissues using RNA-seq data, implying that AhLRR-RLK genes may differ in function. In addition, Arabidopsis plants overexpressing stress-induced AhLRR-RLK265 displayed lower seed germination rates and root lengths compared to wild-type under exogenous ABA treatment. Notably, overexpression of AhLRR-RLK265 enhanced tolerance to salt and drought stresses in transgenic Arabidopsis. Moreover, the AhLRR-RLK265-OE lines were found to have higher activities of superoxide dismutase (SOD), catalase (CAT), and peroxidase (POD) under salt and drought stress treatments. We believe these results may provide valuable information about the function of peanut LRR-RLK genes for further analysis.

Genome-Wide Association Study of Asian and European Common Wheat Accessions for Yield-Related Traits and Stripe Rust Resistance

Identifying novel loci of yield-related traits and resistance to stripe rust (caused by Puccinia striiformis f. sp. tritici) in wheat will help in breeding wheat that can meet projected demands in diverse environmental and agricultural practices. We performed a genome-wide association study with 24,767 single nucleotide polymorphisms (SNPs) in 180 wheat accessions that originated in 16 Asian or European countries between latitudes 30°N and 45°N. We detected seven accessions with desirable yield-related traits and 42 accessions that showed stable, high degrees of stripe rust resistance in multienvironment field assessments. A marker-trait association analysis of yield-related traits detected 18 quantitative trait loci (QTLs) in at least two test environments and two QTLs related to stripe rust resistance in at least three test environments. Five of these QTLs were identified as potentially novel QTLs by comparing their physical locations with those of known QTLs in the Chinese Spring (CS) reference genome RefSeq v1.1 published by the International Wheat Genome Sequencing Consortium; two were for spike length, one was for grain number per spike, one was for spike number, and one was for stripe rust resistance at the adult plant stage. We also identified 14 candidate genes associated with the five novel QTLs. These QTLs and candidate genes will provide breeders with new germplasm and can be used to conduct marker-assisted selection in breeding wheat with improved yield and stripe rust resistance.

Comparison of the impact of two key fungal signalling pathways on Zymoseptoria tritici infection reveals divergent contribution to invasive growth through distinct regulation of infection-associated genes

The lifecycle of Zymoseptoria tritici requires a carefully regulated asymptomatic phase within the wheat leaf following penetration of the mesophyll via stomata. Here we compare the roles in this process of two key fungal signalling pathways, mutants of which were identified through forward genetics due to their avirulence on wheat. Whole-genome resequencing of avirulent Z. tritici T-DNA transformants identified disruptive mutations in ZtBCK1 from the kinase cascade of the cell wall integrity (CWI) pathway, and the adenylate cyclase gene ZtCYR1. Targeted deletion of these genes abolished the pathogenicity of the fungus and led to similar in vitro phenotypes to those associated with disruption of putative downstream kinases, both supporting previous studies and confirming the importance of these pathways in virulence. RNA sequencing was used to investigate the effect of ZtBCK1 and ZtCYR1 deletion on gene expression in both the pathogen and host during infection. ZtBCK1 was found to be required for the adaptation to the host environment, controlling expression of infection-associated secreted proteins, including known virulence factors. Meanwhile, ZtCYR1 is implicated in controlling the switch to necrotrophy, regulating expression of effectors associated with this transition. This represents the first study to compare the influence of CWI and cAMP signalling on in planta transcription of a fungal plant pathogen, providing insights into their differential regulation of candidate effectors during invasive growth.

Receptor-like Kinases (LRR-RLKs) in Response of Plants to Biotic and Abiotic Stresses

Plants live under different biotic and abiotic stress conditions, and, to cope with the adversity and severity, plants have well-developed resistance mechanisms. The mechanism starts with perception of the stimuli followed by molecular, biochemical, and physiological adaptive measures. The family of LRR-RLKs (leucine-rich repeat receptor-like kinases) is one such group that perceives biotic and abiotic stimuli and also plays important roles in different biological processes of development. This has been mostly studied in the model plant, Arabidopsis thaliana, and to some extent in other plants, such as Solanum lycopersicum, Nicotiana benthamiana, Brassica napus, Oryza sativa, Triticum aestivum, Hordeum vulgare, Brachypodium distachyon, Medicago truncatula, Gossypium barbadense, Phaseolus vulgaris, Solanum tuberosum, and Malus robusta. Most LRR-RLKs tend to form different combinations of LRR-RLKs-complexes (dimer, trimer, and tetramers), and some of them were observed as important receptors in immune responses, cell death, and plant development processes. However, less is known about the function(s) of LRR-RLKs in response to abiotic and biotic stresses. Here, we give recent updates about LRR-RLK receptors, specifically focusing on their involvement in biotic and abiotic stresses in the model plant, A. thaliana. Furthermore, the recent studies on LRR-RLKs that are homologous in other plants is also reviewed in relation to their role in triggering stress response processes against biotic and abiotic stimuli and/or in exploring their additional function(s). Furthermore, we present the interactions and combinations among LRR-RLK receptors that have been confirmed through experiments. Moreover, based on GENEINVESTIGATOR microarray database analysis, we predict some potential LRR-RLK genes involved in certain biotic and abiotic stresses whose function and mechanism may be explored.

Genome-Wide Identification and Characterization of Receptor-Like Protein Kinase 1 (RPK1) Gene Family in Triticum aestivum Under Drought Stress

Receptor-like protein kinase1 (RPK1) genes play crucial roles in plant growth and development processes, root architecture, and abiotic stress regulation. A comprehensive study of the RPK1 gene family has not been reported in bread wheat (Triticum aestivum). Here, we reported the genome-wide identification, characterization, and expression patterns of the RPK1 gene family in wheat. Results confirmed 15 TaRPK1 genes, classified mainly into three sub-clades based on a phylogenetic tree. The TaRPK1 genes were mapped on chromosomes 1–3 in the respective A, B, and D genomes. Gene structure, motif conservation, collinearity prediction, and synteny analysis were carried out systematically. A Gene ontology study revealed that TaRPK1 genes play a vital role during molecular and biological processes. We also identified 18 putative miRNAs targeting TaRPK1 genes, suggesting their roles in growth, development, and stress responses. Cis-Regulatory elements interpreted the presence of light-related elements, hormone responsiveness, and abiotic stress-related motifs in the promoter regions. The SWISS_MODEL predicted the successful models of TaRPK1 proteins with at least 30% identity to the template, a widely accepted threshold for successful modeling. In silico expression analysis in different tissues and stages suggested that TaRPK1 genes exhibited the highest expression in root tissues. Moreover, qRT-PCR further validated the higher expression of TaRPK1 genes in roots of drought-tolerant varieties compared to the drought-susceptible variety. Collectively, the present study renders valuable information on the functioning of TaRPK1 genes in wheat that will be useful in further functional validation of these genes in future studies.

Insight into the Roles of Proline-Rich Extensin-like Receptor Protein Kinases of Bread Wheat (Triticum aestivum L.)

Proline-rich extensin-like receptor protein kinases (PERKs) are known for their roles in the developmental processes and stress responses of many plants. We have identified 30 TaPERK genes in the genome of T. aestivum, exploring their evolutionary and syntenic relationship and analyzing their gene and protein structures, various cis-regulatory elements, expression profiling, and interacting miRNAs. The TaPERK genes formed 12 homeologous groups and clustered into four phylogenetic clades. All the proteins exhibited a typical domain organization of PERK and consisted of conserved proline residue repeats and serine-proline and proline-serine repeats. Further, the tyrosine-x-tyrosine (YXY) motif was also found conserved in thirteen TaPERKs. The cis-regulatory elements and expression profiling under tissue developmental stages suggested their role in plant growth processes. Further, the differential expression of certain TaPERK genes under biotic and abiotic stress conditions suggested their involvement in defense responses as well. The interaction of TaPERK genes with different miRNAs further strengthened evidence for their diverse biological roles. In this study, a comprehensive analysis of obtained TaPERK genes was performed, enriching our knowledge of TaPERK genes and providing a foundation for further possible functional analyses in future studies.

Evolutionary and expression dynamics of LRR-RLKs and functional establishment of KLAVIER homolog in shoot mediated regulation of AON in chickpea symbiosis

Chickpea shoot exogenously treated with cytokinin showed stunted phenotype of root, shoot and significantly reduced nodule numbers. Genome-wide identification of LRR-RLKs in chickpea and Medicago resulted in 200 and 371 genes respectively. Gene duplication analysis revealed that LRR-RLKs family expanded through segmental duplications in chickpea and tandem duplications in Medicago. Expression profiling of LRR-RLKs revealed their involvement in cytokinin signaling and plant organ development. Overexpression of KLAVIER ortholog of chickpea, Ca_LRR-RLK147, in roots revealed its localization in the membrane but showed no effect on root nodulation despite increased cle peptide levels. Two findings (i) drastic effect on nodule number by exogenous cytokinin treatment to only shoot and restoration to normal nodulation by treatment to both root and shoot tissue and (ii) no effect on nodule number by overexpression of Ca_LRR-RLK147 establishes the fact that despite presence of cle peptides in root, the function of Ca_LRR-RLK147 was shoot mediated during AON.

Insilco identification and characterization of superoxide dismutase gene family in Brassica rapa

Superoxide Dismutase SODs are defense associated proteins that detoxify ROS and primarily serve as scavengers. They have been described in numerous plant species, but their in-depth characterization in Brassica rapa has not been reported. Therefore, the present investigation on genome wide study of SOD gene family was conducted to identify BrSOD genes, their domain-based organization, gene structure analysis, phylogenetic analysis, intron-exon structure of genes and expression analysis. The sequence characterization of Super oxide dismutase gene family in Brassica rapa, their syntenic associateship of conserved motifs and phylogenetic correlationship, prediction of cis-elements and determing the expression analysis in distinct tissues namely plant callus, root, stem, leaf, flower, and silique under abiotic conditions have been analysed using different software’s. The study on SOD gene family identified 17 BrSOD genes which were grouped into eight BrCu-ZnSODs and nine BrFe-MnSODs domain-based organization. Furthermore, the conserved character of BrSODs were confirmed by intron-exon organisation, motif arrangements and domain architectural investigations. Expression analysis using RNA Sequence data of different developmental stages proclaimed that genes were manifested in all six tissues with an exception of BrCu-ZnSOD3, which was not manifested in roots; however, whose transcript was detected in all other tested tissues. The study has genome wide insight into the occurrence and functional specifications of BrSOD gene family in Brassica rapa that can be potentially utilized in breeding program for resilience to climate change and abiotic stresses tolerance Brassica variety.

Genome-wide identification and evolutionary analysis of RLKs involved in the response to aluminium stress in peanut

BACKGROUND

As an important cash crop, the yield of peanut is influenced by soil acidification and pathogen infection. Receptor-like protein kinases play important roles in plant growth, development and stress responses. However, little is known about the number, location, structure, molecular phylogeny, and expression of RLKs in peanut, and no comprehensive analysis of RLKs in the Al stress response in peanuts have been reported.

RESULTS

A total of 1311 AhRLKs were identified from the peanut genome. The AhLRR-RLKs and AhLecRLKs were further divided into 24 and 35 subfamilies, respectively. The AhRLKs were randomly distributed across all 20 chromosomes in the peanut. Among these AhRLKs, 9.53% and 61.78% originated from tandem duplications and segmental duplications, respectively. The ka/ks ratios of 96.97% (96/99) of tandem duplication gene pairs and 98.78% (646/654) of segmental duplication gene pairs were less than 1. Among the tested tandem duplication clusters, there were 28 gene conversion events. Moreover, all total of 90 Al-responsive AhRLKs were identified by mining transcriptome data, and they were divided into 7 groups. Most of the Al-responsive AhRLKs that clustered together had similar motifs and evolutionarily conserved structures. The gene expression patterns of these genes in different tissues were further analysed, and tissue-specifically expressed genes, including 14 root-specific Al-responsive AhRLKs were found. In addition, all 90 Al-responsive AhRLKs which were distributed unevenly in the subfamilies of AhRLKs, showed different expression patterns between the two peanut varieties (Al-sensitive and Al-tolerant) under Al stress.

CONCLUSIONS

In this study, we analysed the RLK gene family in the peanut genome. Segmental duplication events were the main driving force for AhRLK evolution, and most AhRLKs subject to purifying selection. A total of 90 genes were identified as Al-responsive AhRLKs, and the classification, conserved motifs, structures, tissue expression patterns and predicted functions of Al-responsive AhRLKs were further analysed and discussed, revealing their putative roles. This study provides a better understanding of the structures and functions of AhRLKs and Al-responsive AhRLKs.

Comprehensive analysis of structural, functional, and evolutionary dynamics of Leucine Rich Repeats-RLKs in Thinopyrum elongatum

Leucine Rich Repeats-receptor-like protein kinases (LRR-RLKs) regulate several critical biological processes ranging from growth and development to stress response. Thinopyrum elongatum harbours many desirable traits such as biotic and abiotic stress resistance and therefore commonly used by wheat breeders. In the present investigation, in-silico analysis of LRR-RLKs yielded 589 genes of which 431 were membrane surface RLKs and 158 were receptor like cytoplasmic kinases. An insight into the gene and protein structure revealed quite a conserved nature of these proteins within subgroups. A large expansion in LRR-RLKs was due to tandem and segmental duplication event. Maximum number of tandem and segmentally duplicated pairs was observed in LRR-VI and LRR-XII subfamily, respectively. Furthermore, syntenic analyses revealed that chromosome 6 harboured more (48) tandem duplicated genes while chromosome 7 possessed more (47) segmentally duplicated genes. A detailed analysis about the gene duplication events coupled with expression profiles during Fusarium graminearum infection and water deficiency unravelled the expansion of the gene family with sub functionalization and neofunctionalization. Interaction network analysis showed that LRR-RLKs can heterodimerize upon ligand binding to perform various plant functional attributes.

Identification of Candidate Genes for Root Traits Using Genotype-Phenotype Association Analysis of Near-Isogenic Lines in Hexaploid Wheat (Triticum aestivum L.)

Global wheat (Triticum aestivum L.) production is constrained by different biotic and abiotic stresses, which are increasing with climate change. An improved root system is essential for adaptability and sustainable wheat production. In this study, 10 pairs of near-isogenic lines (NILs)-targeting four genomic regions (GRs) on chromosome arms 4BS, 4BL, 4AS, and 7AL of hexaploid wheat-were used to phenotype root traits in a semi-hydroponic system. Seven of the 10 NIL pairs significantly differed between their isolines for 11 root traits. The NIL pairs targeting qDSI.4B.1 GR varied the most, followed by the NIL pair targeting qDT.4A.1 and QHtscc.ksu-7A GRs. For pairs 5-7 targeting qDT.4A.1 GR, pair 6 significantly differed in the most root traits. Of the 4 NIL pairs targeting qDSI.4B.1 GR, pairs 2 and 4 significantly differed in 3 and 4 root traits, respectively. Pairs 9 and 10 targeting QHtscc.ksu-7A GR significantly differed in 1 and 4 root traits, respectively. Using the wheat 90K Illumina iSelect array, we identified 15 putative candidate genes associated with different root traits in the contrasting isolines, in which two UDP-glycosyltransferase (UGT)-encoding genes, TraesCS4A02G185300 and TraesCS4A02G442700, and a leucine-rich repeat receptor-like protein kinase (LRR-RLK)-encoding gene, TraesCS4A02G330900, also showed important functions for root trait control in other crops. This study characterized, for the first time, that these GRs control root traits in wheat, and identified candidate genes, although the candidate genes will need further confirmation and validation for marker-assisted wheat breeding.

Molecular characterization revealed the role of catalases under abiotic and arsenic stress in bread wheat (Triticum aestivum L.)

Catalases are crucial antioxidant enzymes that reduce the excessive level of H₂O₂ caused by various environmental stresses and metal toxicity and hence protect the plant cells. In this study, a total of ten TaCAT genes, forming three homeologous groups, were identified in the genome of bread wheat (Triticum aestivum L.) and named as per the wheat gene symbolization guidelines. The identified catalases were characterized for various structural and physicochemical features. The proximal active-site (F(D/A)RERIPERVVHAKGASA) and heme-ligand (R(I/V)F(S/A)Y(A/S)DTQ) signature motifs, catalytic residues and peroxisomal targeting peptides were found conserved. Phylogenetic analysis clustered TaCATs into three classes, which showed conserved functional specialization based on their tissue specific expression. Modulated spatio-temporal expression of various TaCAT genes and alteration in total catalase enzyme activity during heat, drought, salt and arsenic (As^III and As^V) treatment suggested their roles in abiotic stress response and arsenic tolerance. Molecular cloning and overexpression of TaCAT3-B gene in Escherichia coli showed tolerance against heat, drought, salt and varied concentrations of As^III and As^V treatments. The results further confirmed their role in stress tolerance and recommended that these genes can be used in future stress management strategies for the development of abiotic and arsenic stress resistant transgenic crops.

Genome-Wide Identification and Evolution of Receptor-Like Kinases (RLKs) and Receptor like Proteins (RLPs) in Brassica juncea

Brassica juncea, an allotetraploid species, is an important germplasm resource for canola improvement, due to its many beneficial agronomic traits, such as heat and drought tolerance and blackleg resistance. Receptor-like kinase (RLK) and receptor-like protein (RLP) genes are two types of resistance gene analogues (RGA) that play important roles in plant innate immunity, stress response and various development processes. In this study, genome wide analysis of RLKs and RLPs is performed in B. juncea. In total, 493 RLKs (LysM-RLKs and LRR-RLKs) and 228 RLPs (LysM-RLPs and LRR-RLPs) are identified in the genome of B. juncea, using RGAugury. Only 13.54% RLKs and 11.79% RLPs are observed to be grouped within gene clusters. The majority of RLKs (90.17%) and RLPs (52.83%) are identified as duplicates, indicating that gene duplications significantly contribute to the expansion of RLK and RLP families. Comparative analysis between B. juncea and its progenitor species, B. rapa and B. nigra, indicate that 83.62% RLKs and 41.98% RLPs are conserved in B. juncea, and RLPs are likely to have a faster evolution than RLKs. This study provides a valuable resource for the identification and characterisation of candidate RLK and RLP genes.

Genome-wide characterization and expression analysis suggested diverse functions of the mechanosensitive channel of small conductance-like (MSL) genes in cereal crops

Mechanosensitive ion channels are pore-forming transmembrane proteins that allow ions to move down their electrochemical gradient in response to mechanical stimuli. They participate in many plant developmental processes including the maintenance of plastid shape, pollen tube growth, etc. Herein, a total of 11, 10, 6, 30, 9, and 8 MSL genes were identified in Aegilops tauschii, Hordeum vulgare, Sorghum bicolor, Triticum aestivum, Triticum urartu, and Zea mays, respectively. These genes were located on various chromosomes of their respective cereal, while MSLs of T. urartu were found on scaffolds. The phylogenetic analysis, subcellular localization, and sequence homology suggested clustering of MSLs into two classes. These genes consisted of cis-regulatory elements related to growth and development, responsive to light, hormone, and stress. Differential expression of various MSL genes in tissue developmental stages and stress conditions revealed their precise role in development and stress responses. Altered expression during CaCl2 stress suggested their role in Ca2+ homeostasis and signaling. The co-expression analysis suggested their interactions with other genes involved in growth, defense responses etc. A comparative expression profiling of paralogous genes revealed either retention of function or pseudo-functionalization. The present study unfolded various characteristics of MSLs in cereals, which will facilitate their in-depth functional characterization in future studies.

Resistance Gene Analogs in the Brassicaceae: Identification, Characterization, Distribution, and Evolution

The Brassicaceae consists of a wide range of species, including important Brassica crop species and the model plant Arabidopsis (Arabidopsis thaliana). Brassica spp. crop diseases impose significant yield losses annually. A major way to reduce susceptibility to disease is the selection in breeding for resistance gene analogs (RGAs). Nucleotide binding site-leucine rich repeats (NLRs), receptor-like kinases (RLKs), and receptor-like proteins (RLPs) are the main types of RGAs; they contain conserved domains and motifs and play specific roles in resistance to pathogens. Here, all classes of RGAs have been identified using annotation and assembly-based pipelines in all available genome annotations from the Brassicaceae, including multiple genome assemblies of the same species where available (total of 32 genomes). The number of RGAs, based on genome annotations, varies within and between species. In total 34,065 RGAs were identified, with the majority being RLKs (21,691), then NLRs (8,588) and RLPs (3,786). Analysis of the RGA protein sequences revealed a high level of sequence identity, whereby 99.43% of RGAs fell into several orthogroups. This study establishes a resource for the identification and characterization of RGAs in the Brassicaceae and provides a framework for further studies of RGAs for an ultimate goal of assisting breeders in improving resistance to plant disease.

Thaumatin-like protein kinases: Molecular characterization and transcriptional profiling in five cereal crops

Thaumatin-like protein kinases (TLPKs) are defense related proteins having antimicrobial property. Herein, we identified two TLPKs in the genome of Brachypodium distachyon and Oryza sativa, four in Hordeum vulgare and Sorghum bicolor, and 16 in Triticum aestivum. All the TLPKs were located at only one chromosome in each plant except T. aestivum, where they were located on chromosome 2 and chromosome 3. Paralogous analysis suggested the occurrence of one duplication event (DE) in each B. distachyon and O. sativa, two in H. vulgare while four DEs in T. aestivum genome during the evolution of TLPKs. The majority of TLPKs were intron less, while a few contains one or two introns. The introns were found in each 0, 1 and 2 phase. Protein structure analysis suggested the occurrence of a thaumatin and a kinase domain with a transmembrane (TM) helix in each TLPK. Further, a thaumatin family signature motif "GX[GF]XCXT[GA]DCX(1,2)GX(2,3)C", a "REDDD" motif and 16 cysteine residues were found conserved in the majority of TLPKs. Expression analysis indicated variable expression of TLPKs in various tissues of different cereal crops. They were high expressing in reproductive tissues in B. distachyon, while in leaves in T. aestivum. Modulated expression of TaTLPKs in the presence of fungal pathogen, and heat, drought and salt stress in T. aestivum suggested their roles in stress response. Co-expression analysis showed interaction of TLPKs with various development and stress related genes. The results indicated diverse roles of TLPKs, which can be utilized for the development of eco-friendly pest resistant crops in future.

Genome-Wide Association Study for Multiple Biotic Stress Resistance in Synthetic Hexaploid Wheat

Genetic resistance against biotic stress is a major goal in many wheat breeding programs. However, modern wheat cultivars have a limited genetic variation for disease and pest resistance and there is always a possibility of the evolution of new diseases and pests to overcome previously identified resistance genes. A total of 125 synthetic hexaploid wheats (SHWs; 2n = 6x = 42, AABBDD, Triticum aestivum L.) were characterized for resistance to fungal pathogens that cause wheat rusts (leaf; Puccinia triticina, stem; P. graminis f.sp. tritici, and stripe; P. striiformis f.sp. tritici) and crown rot (Fusarium spp.); cereal cyst nematode (Heterodera spp.); and Hessian fly (Mayetiola destructor). A wide range of genetic variation was observed among SHWs for multiple (two to five) biotic stresses and 17 SHWs that were resistant to more than two stresses. The genomic regions and potential candidate genes conferring resistance to these biotic stresses were identified from a genome-wide association study (GWAS). This GWAS study identified 124 significant marker-trait associations (MTAs) for multiple biotic stresses and 33 of these were found within genes. Furthermore, 16 of the 33 MTAs present within genes had annotations suggesting their potential role in disease resistance. These results will be valuable for pyramiding novel genes/genomic regions conferring resistance to multiple biotic stresses from SHWs into elite bread wheat cultivars and providing further insights on a wide range of stress resistance in wheat.

Hidden genomic evolution in a morphospecies-The landscape of rapidly evolving genes in Tetrahymena

A morphospecies is defined as a taxonomic species based wholly on morphology, but often morphospecies consist of clusters of cryptic species that can be identified genetically or molecularly. The nature of the evolutionary novelty that accompanies speciation in a morphospecies is an intriguing question. Morphospecies are particularly common among ciliates, a group of unicellular eukaryotes that separates 2 kinds of nuclei—the silenced germline nucleus (micronucleus [MIC]) and the actively expressed somatic nucleus (macronucleus [MAC])—within a common cytoplasm. Because of their very similar morphologies, members of the Tetrahymena genus are considered a morphospecies. We explored the hidden genomic evolution within this genus by performing a comprehensive comparative analysis of the somatic genomes of 10 species and the germline genomes of 2 species of Tetrahymena. These species show high genetic divergence; phylogenomic analysis suggests that the genus originated about 300 million years ago (Mya). Seven universal protein domains are preferentially included among the species-specific (i.e., the youngest) Tetrahymena genes. In particular, leucine-rich repeat (LRR) genes make the largest contribution to the high level of genome divergence of the 10 species. LRR genes can be sorted into 3 different age groups. Parallel evolutionary trajectories have independently occurred among LRR genes in the different Tetrahymena species. Thousands of young LRR genes contain tandem arrays of exactly 90-bp exons. The introns separating these exons show a unique, extreme phase 2 bias, suggesting a clonal origin and successive expansions of 90-bp–exon LRR genes. Identifying LRR gene age groups allowed us to document a Tetrahymena intron length cycle. The youngest 90-bp exon LRR genes in T. thermophila are concentrated in pericentromeric and subtelomeric regions of the 5 micronuclear chromosomes, suggesting that these regions act as genome innovation centers. Copies of a Tetrahymena Long interspersed element (LINE)-like retrotransposon are very frequently found physically adjacent to 90-bp exon/intron repeat units of the youngest LRR genes. We propose that Tetrahymena species have used a massive exon-shuffling mechanism, involving unequal crossing over possibly in concert with retrotransposition, to create the unique 90-bp exon array LRR genes.

Genomic comparison of ten morphologically very similar species of ciliate from the genus Tetrahymena reveals how parallel microevolutionary processes have shaped their genomes and created unique genes through retrotransposition.

Genomic dissection and transcriptional profiling of Cysteine-rich receptor-like kinases in five cereals and functional characterization of TaCRK68-A

Cysteine-rich receptor-like kinases (CRK) constitute one of the largest subfamily of receptor-like kinases, which play crucial roles in plant development and stress response. In total, 43, 37, 36, 38 and 170 CRK genes including duplicated genes were identified in the genome of Brachypodium distachyon, Hordeum vulgare, Oryza sativa, Sorghum bicolor and Triticum aestivum, respectively. These CRK proteins were tightly clustered into four phylogenetic groups and exhibited close syntenic relationship among orthologous genes. Majority of CRK proteins contain a transmembrane domain for plasma membrane localization. The organization of exon/intron, domains and motifs were variably conserved. Tissue-specific expression suggested the involvement of certain CRK genes in plant development. Modulated expression revealed their specific stress-responsive functions. Co-expression and interaction analysis indicated their role in signaling. Ks value and divergence time analysis suggested duplication of TaCRK genes before the hybridization of T. aestivum sub-genomes. Expression comparison of duplicated TaCRK genes revealed functional retention, neofunctionalization or pseudo-functionalization. Recombinant expression of a stress-responsive gene TaCRK68-A in Escherichia coli and Saccharomyces cerevisiae displayed enhanced tolerance against heat, drought, cold and salinity stresses. The study suggested vital functions of CRKs during development and stresses, and provides the basis for functional characterization of each gene in future studies.

Genome-wide association study of resistance to stripe rust (Puccinia striiformis f. sp. tritici) in Sichuan wheat

BACKGROUND

Stripe rust (also called yellow rust) is a common and serious fungal disease of wheat (Triticum aestivum L.) caused by Puccinia striiformis f. sp. tritici. The narrow genetic basis of modern wheat cultivars and rapid evolution of the rust pathogen have been responsible for periodic and devastating epidemics of wheat rust diseases. In this study, we conducted a genome-wide association study with 44,059 single nucleotide polymorphism markers to identify loci associated with resistance to stripe rust in 244 Sichuan wheat accessions, including 79 landraces and 165 cultivars, in six environments.

RESULTS

In all the field assessments, 24 accessions displayed stable high resistance to stripe rust. Significant correlations among environments were observed for both infection (IT) and disease severity (DS), and high heritability levels were found for both IT and DS. Using mixed linear models, 12 quantitative trait loci (QTLs) significantly associated with IT and/or DS were identified. Two QTLs were mapped on chromosomes 5AS and 5AL and were distant from previously identified stripe rust resistance genes or QTL regions, indicating that they may be novel resistance loci.

CONCLUSIONS

Our results revealed that resistance alleles to stripe rust were accumulated in Sichuan wheat germplasm, implying direct or indirect selection for improved stripe rust resistance in elite wheat breeding programs. The identified stable QTLs or favorable alleles could be important chromosome regions in Sichuan wheat that controlled the resistance to stripe rust. These markers can be used molecular marker-assisted breeding of Sichuan wheat cultivars, and will be useful in the ongoing effort to develop new wheat cultivars with strong resistance to stripe rust.

Organization and evolution of the chalcone synthase gene family in bread wheat and relative species

Background

Flavonoid compounds are secondary plant metabolites, having a functional importance in plant development, protection from pathogens and unfavorable environmental factors. Chalcone synthase (CHS) is a key enzyme in the biosynthesis of flavonoids; it is involved in biosynthesis of all classes of flavonoid compounds. Nevertheless, the Chs gene family in bread wheat (Triticum aestivum L.) has been not characterized yet. The aim of the current study was to investigate structural and functional organization of the Chs genes and evolution of this gene family in bread wheat and relative species.

Results

The nucleotide sequences of the eight Chs copies in T. aestivum were identified. Among them, two homoeologous sets of the Chs genes were located on the short (Chs-A1, −B1, −D1) and the long (Chs-A4, −B4, −D4) arms of homoeologous group 2 chromosomes. Two paralogous gene copies in the B-genome (Chs-B2, −B3) were located in the distal regions of 2BS chromosome. To clarify the origin of Chs duplications in the B-genome the phylogenetic analysis with the Chs sequences of Triticum and Aegilops species carrying ancestral genomes was conducted. It was estimated that the first duplication event occurred in the genome of the common ancestor of Triticum and Aegilops genera about 10–12 million years ago (MYA), then another copy was formed in the ancestor of the B-genome about 6–7 MYA. A homology modeling revealed high sequence similarity of bread wheat CHS enzymes. A number of short deletions in coding regions of some Chs sequences are not expected to have any significant functional effects. Estimation of transcriptional activity of the Chs copies along with a comparative analysis of their promoters structure suggested their functional specialization, which likely contributed to the maintaining of the duplicated Chs genes in wheat genome.

Conclusions

From possible ten Chs copies in bread wheat genome, eight members of this family retained their intact structure and activity, while two copies appear to be lost at the level of diploid and tetraploid ancestors. Transcriptional assay along with a comparative analysis of the cis-regulatory elements revealed their functional diversification. The multiple functions supported by the Chs family are assumed to be a driving force for duplications of the Chs gene and their retention in plant genome.

Electronic supplementary material

The online version of this article (10.1186/s12863-019-0727-y) contains supplementary material, which is available to authorized users.

Identification, characterization and expression profiling of cation-proton antiporter superfamily in Triticum aestivum L. and functional analysis of TaNHX4-B

The monovalent cation proton antiporter (CPA) superfamily comprises Na⁺/H⁺ exchanger (NHX), K⁺ efflux antiporter (KEA), and cation/H⁺ exchanger (CHX) family proteins, which play vital functions in plants. A total of 107 TaCPA proteins were identified in Triticum aestivum, and phylogenetically classified into 35 TaNHX, 24 TaKEA and 48 TaCHX proteins. These families had representatives derived from all three sub-genomes. TaKEA genes consisted of higher number of exons, followed by TaNHXs and TaCHXs. The occurrence of about 10 transmembrane regions and higher composition of helices and coils support their membrane-bound and hydrophobic nature. Diverse expression in various tissues and modulated expression under stress conditions suggested their role in development and in response to stress. Co-expression analyses revealed their complex interaction networks. Expression of TaNHX4-B.1 and TaNHX4-B.4 facilitated differential abiotic stress tolerance to Escherichia coli. Our study provides comprehensive information about CPA genes, which would be useful in their future functional characterization.

Two-State Co-Expression Network Analysis to Identify Genes Related to Salt Tolerance in Thai rice

Khao Dawk Mali 105 (KDML105) rice is one of the most important crops of Thailand. It is a challenging task to identify the genes responding to salinity in KDML105 rice. The analysis of the gene co-expression network has been widely performed to prioritize significant genes, in order to select the key genes in a specific condition. In this work, we analyzed the two-state co-expression networks of KDML105 rice under salt-stress and normal grown conditions. The clustering coefficient was applied to both networks and exhibited significantly different structures between the salt-stress state network and the original (normal-grown) network. With higher clustering coefficients, the genes that responded to the salt stress formed a dense cluster. To prioritize and select the genes responding to the salinity, we investigated genes with small partners under normal conditions that were highly expressed and were co-working with many more partners under salt-stress conditions. The results showed that the genes responding to the abiotic stimulus and relating to the generation of the precursor metabolites and energy were the great candidates, as salt tolerant marker genes. In conclusion, in the case of the complexity of the environmental conditions, gaining more information in order to deal with the co-expression network provides better candidates for further analysis.

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.

Genome-wide identification and characterization of LRR-RLKs reveal functional conservation of the SIF subfamily in cotton (Gossypium hirsutum)

BACKGROUND

As one of the largest subfamilies of the receptor-like protein kinases (RLKs) in plants, Leucine Rich Repeats-RLKs (LRR-RLKs) are involved in many critical biological processes including growth, development and stress responses in addition to various physiological roles. Arabidopsis contains 234 LRR-RLKs, and four members of Stress Induced Factor (SIF) subfamily (AtSIF1-AtSIF4) which are involved in abiotic and biotic stress responses. Herein, we aimed at identification and functional characterization of SIF subfamily in cultivated tetraploid cotton Gossypium hirsutum.

RESULTS

Genome-wide analysis of cotton LRR-RLK gene family identified 543 members and phylogenetic analysis led to the identification of 6 cotton LRR-RLKs with high homology to Arabidopsis SIFs. Of the six SIF homologs, GhSIF1 is highly conserved exhibiting 46-47% of homology with AtSIF subfamily in amino acid sequence. The GhSIF1 was transiently silenced using Virus-Induced Gene Silencing system specifically targeting the 3' Untranslated Region. The transiently silenced cotton seedlings showed enhanced salt tolerance compared to the control plants. Further, the transiently silenced plants showed better growth, lower electrolyte leakage, and higher chlorophyll and biomass contents.

CONCLUSIONS

Overall, 543 LRR-RLK genes were identified using genome-wide analysis in cultivated tetraploid cotton G. hirsutum. The present investigation also demonstrated the conserved salt tolerance function of SIF family member in cotton. The GhSIF1 gene can be knocked out using genome editing technologies to improve salt tolerance in cotton.

Molecular characterization and differential expression suggested diverse functions of P-type II Ca2+ATPases in Triticum aestivum L

Background

Plant P-type II Ca²⁺ATPases are formed by two distinct groups of proteins (ACAs and ECAs) that perform pumping of Ca²⁺ outside the cytoplasm during homeostasis, and play vital functions during development and stress management. In the present study, we have performed identification and characterisation of P-type II Ca²⁺ATPase gene family in an important crop plant Triticum aestivum.

Results

Herein, a total of 33 TaACA and 9 TaECA proteins were identified from the various chromosomes and sub-genomes of Triticum aestivum. Phylogenetic analysis revealed clustering of the homoeologous TaACA and TaECA proteins into 11 and 3 distinct groups that exhibited high sequence homology and comparable structural organization as well. Both TaACA and TaECA group proteins consisted of eight to ten transmembrane regions, and their respective domains and motifs. Prediction of sub-cellular localization was found variable for most of the proteins; moreover, it was consistent with the evolutionarily related proteins from rice and Arabidopsis in certain cases. The occurrence of assorted sets of cis-regulatory elements indicated their diverse functions. The differential expression of various TaACA and TaECA genes during developmental stages suggested their roles in growth and development. The modulated expression during heat, drought, salt and biotic stresses along with the occurrence of various stress specific cis-regulatory elements suggested their association with stress response. Interaction of these genes with numerous development and stress related genes indicated their decisive role in various biological processes and signaling.

Conclusion

T. aestivum genome consisted of a maximum of 42 P-type II Ca²⁺ATPase genes, derived from each A, B and D sub-genome. These genes may play diverse functions during plant growth and development. They may also be involved in signalling during abiotic and biotic stresses. The present study provides a comprehensive insight into the role of P-type II Ca²⁺ATPase genes in T. aestivum. However, the specific function of each gene needs to be established, which could be utilized in future crop improvement programs.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-4792-9) contains supplementary material, which is available to authorized users.

Duplication and diversification of lectin receptor-like kinases (LecRLK) genes in soybean

Lectin receptor-like kinases (LecRLKs) play important roles in plant development and stress responses. Although genome-wide studies of LecRLKs have been performed in several species, a comprehensive analysis including evolutionary, structural and functional analysis has not been carried out in soybean (Glycine max). In this study, we identified 185 putative LecRLK genes in the soybean genome, including 123 G-type, 60 L-type and 2 C-type LecRLK genes. Tandem duplication and segmental duplication appear to be the main mechanisms of gene expansion in the soybean LecRLK (GmLecRLK) gene family. According to our phylogenetic analysis, G-type and L-type GmLecRLK genes can be organized into fourteen and eight subfamilies, respectively. The subfamilies within the G-type GmLecRLKs differ from each other in gene structure and/or protein domains and motifs, which indicates that the subfamilies have diverged. The evolution of L-type GmLecRLKs has been more conservative: most genes retain the same gene structures and nearly the same protein domain and motif architectures. Furthermore, the expression profiles of G-type and L-type GmLecRLK genes show evidence of functional redundancy and divergence within each group. Our results contribute to a better understanding of the evolution and function of soybean LecRLKs and provide a framework for further functional investigation of them.

Gene architecture and expression analyses provide insights into the role of glutathione peroxidases (GPXs) in bread wheat (Triticum aestivum L.)

Glutathione peroxidases (GPXs) are redox sensor proteins that maintain a steady-state of H₂O₂ in plant cells. They exhibit distinct sub-cellular localization and have diverse functionality in response to different stimuli. In this study, a total of 14 TaGPX genes and three splice variants were identified in the genome of Triticum aestivum and evaluated for various physicochemical properties. The TaGPX genes were scattered on the various chromosomes of the A, B, and D sub-genomes and clustered into five homeologous groups based on high sequence homology. The majority of genes were derived from the B sub-genome and localized on chromosome 2. The intron-exon organization, motif and domain architecture, and phylogenetic analyses revealed the conserved nature of TaGPXs. The occurrence of both development-related and stress-responsive cis-acting elements in the promoter region, the differential expression of these genes during various developmental stages, and the modulation of expression in the presence of biotic and abiotic stresses suggested their diverse role in T. aestivum. The majority of TaGPX genes showed higher expression in various leaf developmental stages. However, TaGPX1-A1 was upregulated in the presence of each abiotic stress treatment. A co-expression analysis revealed the interaction of TaGPXs with numerous development and stress-related genes, which indicated their vital role in numerous biological processes. Our study revealed the opportunities for further characterization of individual TaGPX proteins, which might be useful in designing future crop improvement strategies.

Mapping and characterization of wheat stem rust resistance genes SrTm5 and Sr60 from Triticum monococcum

KEY MESSAGE

The new stem rust resistance gene Sr60 was fine-mapped to the distal region of chromosome arm 5AmS, and the TTKSK-effective gene SrTm5 could be a new allele of Sr22. The emergence and spread of new virulent races of the wheat stem rust pathogen (Puccinia graminis f. sp. tritici; Pgt), including the Ug99 race group, is a serious threat to global wheat production. In this study, we mapped and characterized two stem rust resistance genes from diploid wheat Triticum monococcum accession PI 306540. We mapped SrTm5, a previously postulated gene effective to Ug99, on chromosome arm 7AmL, completely linked to Sr22. SrTm5 displayed a different race specificity compared to Sr22 indicating that they are distinct. Sequencing of the Sr22 homolog in PI 306540 revealed a novel haplotype. Characterization of the segregating populations with Pgt race QFCSC revealed an additional resistance gene on chromosome arm 5AmS that was assigned the official name Sr60. This gene was also effective against races QTHJC and SCCSC but not against TTKSK (a Ug99 group race). Using two large mapping populations (4046 gametes), we mapped Sr60 within a 0.44 cM interval flanked by sequenced-based markers GH724575 and CJ942731. These two markers delimit a 54.6-kb region in Brachypodium distachyon chromosome 4 and a 430-kb region in the Chinese Spring reference genome. Both regions include a leucine-rich repeat protein kinase (LRRK123.1) that represents a potential candidate gene. Three CC-NBS-LRR genes were found in the colinear Brachypodium region but not in the wheat genome. We are currently developing a Bacterial Artificial Chromosome library of PI 306540 to determine which of these candidate genes are present in the T. monococcum genome and to complete the cloning of Sr60.

A Pathogen-Responsive Leucine Rich Receptor Like Kinase Contributes to Fusarium Resistance in Cereals

Receptor-like kinases form the largest family of receptors in plants and play an important role in recognizing pathogen-associated molecular patterns and modulating the plant immune responses to invasive fungi, including cereal defenses against fungal diseases. But hitherto, none have been shown to modulate the wheat response to the economically important Fusarium head blight (FHB) disease of small-grain cereals. Homologous genes were identified on barley chromosome 6H (HvLRRK-6H) and wheat chromosome 6DL (TaLRRK-6D), which encode the characteristic domains of surface-localized receptor like kinases. Gene expression studies validated that the wheat TaLRRK-6D is highly induced in heads as an early response to both the causal pathogen of FHB disease, Fusarium graminearum, and its' mycotoxic virulence factor deoxynivalenol. The transcription of other wheat homeologs of this gene, located on chromosomes 6A and 6B, was also up-regulated in response to F. graminearum. Virus-induced gene silencing (VIGS) of the barley HvLRRK-6H compromised leaf defense against F. graminearum. VIGS of TaLRRK-6D in two wheat cultivars, CM82036 (resistant to FHB disease) and cv. Remus (susceptible to FHB), confirmed that TaLRRK-6D contributes to basal resistance to FHB disease in both genotypes. Although the effect of VIGS did not generally reduce grain losses due to FHB, this experiment did reveal that TaLRRK-6D positively contributes to grain development. Further gene expression studies in wheat cv. Remus indicated that VIGS of TaLRRK-6D suppressed the expression of genes involved in salicylic acid signaling, which is a key hormonal pathway involved in defense. Thus, this study provides the first evidence of receptor like kinases as an important component of cereal defense against Fusarium and highlights this gene as a target for enhancing cereal resistance to FHB disease.

Survey of High Throughput RNA-Seq Data Reveals Potential Roles for lncRNAs during Development and Stress Response in Bread Wheat

Long non-coding RNAs (lncRNAs) are a family of regulatory RNAs that play essential role in the various developmental processes and stress responses. Recent advances in sequencing technology and computational methods enabled identification and characterization of lncRNAs in certain plant species, but they are less known in Triticum aestivum (bread wheat). Herein, we analyzed 52 RNA seq data (>30 billion reads) and identified 44,698 lncRNAs in T. aestivum genome, which were characterized in comparison to the coding sequences (mRNAs). Similar to the mRNAs, lncRNAs were also derived from each sub-genome and chromosome, and showed tissue developmental stage specific and differential expression, as well. The modulated expression of lncRNAs during abiotic stresses like heat, drought, and salt indicated their putative role in stress response. The co-expression of lncRNAs with vital mRNAs including various transcription factors and enzymes involved in Abscisic acid (ABA) biosynthesis, and gene ontology mapping inferred their regulatory roles in numerous biological processes. A few lncRNAs were predicted as precursor (19 lncRNAs), while some as target mimics (1,047 lncRNAs) of known miRNAs involved in various regulatory functions. The results suggested numerous functions of lncRNAs in T. aestivum, and unfolded the opportunities for functional characterization of individual lncRNA in future studies.

Ca2+/Cation Antiporters (CaCA): Identification, Characterization and Expression Profiling in Bread Wheat (Triticum aestivum L.)

The Ca2+/cation antiporters (CaCA) superfamily proteins play vital function in Ca2+ ion homeostasis, which is an important event during development and defense response. Molecular characterization of these proteins has been performed in certain plants, but they are still not characterized in Triticum aestivum (bread wheat). Herein, we identified 34 TaCaCA superfamily proteins, which were classified into TaCAX, TaCCX, TaNCL, and TaMHX protein families based on their structural organization and evolutionary relation with earlier reported proteins. Since the T. aestivum comprises an allohexaploid genome, TaCaCA genes were derived from each A, B, and D subgenome and homeologous chromosome (HC), except chromosome-group 1. Majority of genes were derived from more than one HCs in each family that were considered as homeologous genes (HGs) due to their high similarity with each other. These HGs showed comparable gene and protein structures in terms of exon/intron organization and domain architecture. Majority of TaCaCA proteins comprised two Na_Ca_ex domains. However, TaNCLs consisted of an additional EF-hand domain with calcium binding motifs. Each TaCaCA protein family consisted of about 10 transmembrane and two α-repeat regions with specifically conserved signature motifs except TaNCL, which had single α-repeat. Variable expression of most of the TaCaCA genes during various developmental stages suggested their specified role in development. However, constitutively high expression of a few genes like TaCAX1-A and TaNCL1-B indicated their role throughout the plant growth and development. The modulated expression of certain genes during biotic (fungal infections) and abiotic stresses (heat, drought, salt) suggested their role in stress response. Majority of TaCCX and TaNCL family genes were found highly affected during various abiotic stresses. However, the role of individual gene needs to be established. The present study unfolded the opportunity for detail functional characterization of TaCaCA proteins and their utilization in future crop improvement programs.

Duplication and Divergence of Leucine-Rich Repeat Receptor-Like Protein Kinase (LRR-RLK) Genes in Basal Angiosperm Amborella trichopoda

Leucine-rich repeat receptor-like protein kinases (LRR-RLKs) are the largest group of receptor-like kinases, which are one of the largest protein superfamilies in plants, and play crucial roles in development and stress responses. Although the evolution of LRR-RLK families has been investigated in some eudicot and monocot plants, no comprehensive evolutionary studies have been performed for these genes in basal angiosperms like Amborella trichopoda. In this study, we identified 94 LRR-RLK genes in the genome of A. trichopoda. The number of LRR-RLK genes in the genome of A. trichopoda is only 17-50% of that of several eudicot and monocot species. Tandem duplication and whole-genome duplication have made limited contributions to the expansion of LRR-RLK genes in A. trichopoda. According to the phylogenetic analysis, all A. trichopoda LRR-RLK genes can be organized into 18 subfamilies, which roughly correspond to the LRR-RLK subfamilies defined in Arabidopsis thaliana. Most LRR-RLK subfamilies are characterized by highly conserved protein structures, motif compositions, and gene structures. The unique gene structure, protein structures, and protein motif compositions of each subfamily provide evidence for functional divergence among LRR-RLK subfamilies. Moreover, the expression data of LRR-RLK genes provided further evidence for the functional diversification of them. In addition, selection analyses showed that most LRR-RLK protein sites are subject to purifying selection. Our results contribute to a better understanding of the evolution of LRR-RLK gene family in angiosperm and provide a framework for further functional investigation on A. trichopoda LRR-RLKs.