PLANT PROTEIN SERINE/THREONINE KINASES: Classification and Functions

RNA-seq analysis reveals transcriptome reprogramming and alternative splicing during early response to salt stress in tomato root

Salt stress is one of the dominant abiotic stress conditions that cause severe damage to plant growth and, in turn, limiting crop productivity. It is therefore crucial to understand the molecular mechanism underlying plant root responses to high salinity as such knowledge will aid in efforts to develop salt-tolerant crops. Alternative splicing (AS) of precursor RNA is one of the important RNA processing steps that regulate gene expression and proteome diversity, and, consequently, many physiological and biochemical processes in plants, including responses to abiotic stresses like salt stress. In the current study, we utilized high-throughput RNA-sequencing to analyze the changes in the transcriptome and characterize AS landscape during the early response of tomato root to salt stress. Under salt stress conditions, 10,588 genes were found to be differentially expressed, including those involved in hormone signaling transduction, amino acid metabolism, and cell cycle regulation. More than 700 transcription factors (TFs), including members of the MYB, bHLH, and WRKY families, potentially regulated tomato root response to salt stress. AS events were found to be greatly enhanced under salt stress, where exon skipping was the most prevalent event. There were 3709 genes identified as differentially alternatively spliced (DAS), the most prominent of which were serine/threonine protein kinase, pentatricopeptide repeat (PPR)-containing protein, E3 ubiquitin-protein ligase. More than 100 DEGs were implicated in splicing and spliceosome assembly, which may regulate salt-responsive AS events in tomato roots. This study uncovers the stimulation of AS during tomato root response to salt stress and provides a valuable resource of salt-responsive genes for future studies to improve tomato salt tolerance.

Characteristics of ZjCIPKs and ZjbHLH74-ZjCIPK5 regulated cold tolerance in jujube

CIPKs are kind of serine/threonine (Ser/Thr) protein kinases which play important roles in response to biotic and abiotic stresses, and in plant growth and development. However, CIPKs in jujube (Ziziphus jujuba Mill.) had limited information, especially regarding their response to cold stress. In the current study, a total of 18 ZjCIPKs were identified in jujube genome which unevenly distributed on seven chromosomes. Conserved motif and gene structural analysis depicted them with conserved DEGLSA and APE motifs and similar structures. Phylogenetic analysis indicated that CIPKs were classified into five subgroups (I-V). In addition, three pairs of ZjCIPKs exhibited tandem duplication while the segmental duplication of ZjCIPKs was not identified. Study on the cis-acting elements indicted that stress or hormone related cis-acting elements were distributed unevenly on ZjCIPKs promoters and most ZjCIPKs were down- or up-regulated by the cold stress. VIGS induced silencing of ZjCIPK5 decreased the cold tolerance of sour jujube. Subcellular location analysis showed ZjCIPK5 located in nucleus. Moreover, transcription factor ZjbHLH74 which was induced at 6 h under cold stress could interact with the promoter of ZjCIPK5 to regulate jujube cold tolerance. These findings provided insights to a molecular basis of CIPK5 in jujube cold tolerance breeding for future.

Genomic decoding of Theobroma grandiflorum (cupuassu) at chromosomal scale: evolutionary insights for horticultural innovation

BACKGROUND

Theobroma grandiflorum (Malvaceae), known as cupuassu, is a tree indigenous to the Amazon basin, valued for its large fruits and seed pulp, contributing notably to the Amazonian bioeconomy. The seed pulp is utilized in desserts and beverages, and its seed butter is used in cosmetics. Here, we present the sequenced telomere-to-telomere genome of cupuassu, disclosing its genomic structure, evolutionary features, and phylogenetic relationships within the Malvaceae family.

FINDINGS

The cupuassu genome spans 423 Mb, encodes 31,381 genes distributed in 10 chromosomes, and exhibits approximately 65% gene synteny with the Theobroma cacao genome, reflecting a conserved evolutionary history, albeit punctuated with unique genomic variations. The main changes are pronounced by bursts of long-terminal repeat retrotransposons at postspecies divergence, retrocopied and singleton genes, and gene families displaying distinctive patterns of expansion and contraction. Furthermore, positively selected genes are evident, particularly among retained and dispersed tandem and proximal duplicated genes associated with general fruit and seed traits and defense mechanisms, supporting the hypothesis of potential episodes of subfunctionalization and neofunctionalization following duplication, as well as impact from distinct domestication process. These genomic variations may underpin the differences observed in fruit and seed morphology, ripening, and disease resistance between cupuassu and the other Malvaceae species.

CONCLUSIONS

The cupuassu genome offers a foundational resource for both breeding improvement and conservation biology, yielding insights into the evolution and diversity within the genus Theobroma.

Characterization of PIF4 Phosphorylation by SPA1

The PHYTOCHROME INTERACTING FACTORs (PIFs) play pivotal roles in regulating thermo- and photo-morphogenesis in Arabidopsis. One of the main hubs in thermomorphogenesis is PIF4, which regulates plant development under high ambient temperature along with other PIFs. PIF4 enhances its own transcription and PIF4 protein is stabilized under high ambient temperature. However, the mechanisms of thermo-stabilization of PIF4 are less understood. Recently, it was shown that SUPPRESSOR OF PHYA-105 1 (SPA1) can function as a serine/threonine kinase to phosphorylate PIF4 in vitro, and the phosphorylated form of PIF4 is more stable under high ambient temperature conditions. In this chapter, we describe the in vitro kinase assay of PIF4 by SPA1. In principle, this protocol can be applied for other putative substrates and kinases.

Genome-wide association analysis identified molecular markers and candidate genes for flower traits in Chinese orchid (Cymbidium sinense)

The orchid, the champagne of flowers, brings luxury, elegance, and novelty to nature. Cymbidium sinense is a symbol of gigantic floral variability on account of wavering shapes and sizes of floral organs, although marker-trait association (MTA) has not been studied for its floral traits. We evaluated markers associated with 14 floral traits of C. sinense through a genome-wide association study (GWAS) of 195 accessions. A total of 65 318 522 single-nucleotide polymorphisms (SNPs) and 3 906 176 insertion/deletion (InDel) events were identified through genotyping-by-sequencing. Among these, 4694 potential SNPs and 477 InDels were identified as MTAs at -log10 P > 5. The genes related to these SNPs and InDels were largely associated with floral regulators, hormonal pathways, cell division, and metabolism, playing essential roles in tailoring floral morphology. Moreover, 20 candidate SNPs/InDels linked to 11 genes were verified, 8 of which were situated on exons, one was located in the 5'-UTR and two were positioned in introns. Here, the multitepal trait-related gene RABBIT EARS (RBE) was found to be the most crucial gene. We analyzed the role of CsRBE in the regulation of flower-related genes via efficient transient overexpression in C. sinense protoplasts, and found that the floral homeotic genes CsAP3 and CsPI, as well as organ boundary regulators, including CsCUC and CsTCP genes, were regulated by CsRBE. Thus, we obtained key gene loci for important ornamental traits of orchids using genome-wide association analysis of populations with natural variation. The findings of this study can do a great deal to expedite orchid breeding programs for shape variability.

Genomic and transcriptomic analyses reveal polygenic architecture for ecologically important traits in aspen (Populus tremuloides Michx.)

Intraspecific genetic variation in foundation species such as aspen (Populus tremuloides Michx.) shapes their impact on forest structure and function. Identifying genes underlying ecologically important traits is key to understanding that impact. Previous studies, using single-locus genome-wide association (GWA) analyses to identify candidate genes, have identified fewer genes than anticipated for highly heritable quantitative traits. Mounting evidence suggests that polygenic control of quantitative traits is largely responsible for this "missing heritability" phenomenon. Our research characterized the genetic architecture of 30 ecologically important traits using a common garden of aspen through genomic and transcriptomic analyses. A multilocus association model revealed that most traits displayed a highly polygenic architecture, with most variation explained by loci with small effects (likely below the detection levels of single-locus GWA methods). Consistent with a polygenic architecture, our single-locus GWA analyses found only 38 significant SNPs in 22 genes across 15 traits. Next, we used differential expression analysis on a subset of aspen genets with divergent concentrations of salicinoid phenolic glycosides (key defense traits). This complementary method to traditional GWA discovered 1243 differentially expressed genes for a polygenic trait. Soft clustering analysis revealed three gene clusters (241 candidate genes) involved in secondary metabolite biosynthesis and regulation. Our work reveals that ecologically important traits governing higher-order community- and ecosystem-level attributes of a foundation forest tree species have complex underlying genetic structures and will require methods beyond traditional GWA analyses to unravel.

Genome-wide association study of resistance to anthracnose in pepper (Capsicum chinense) germplasm

BACKGROUND

Anthracnose is a fungal disease caused by Colletotrichum spp. that has a significant impact on worldwide pepper production. Colletotrichum scovillei is the most common pathogenic anthracnose-causing species in the Republic of Korea.

RESULTS

The resistances of 197 pepper (Capsicum chinense) accessions deposited in Korea's National Agrobiodiversity Center were evaluated for their response against the virulent pathogens Colletotrichum acutatum isolate 'KSCa-1' and C. scovillei isolate 'Hana') in the field and in vitro methods for three consecutive years (2018 to 2020). The severity of the disease was recorded and compared between inoculation methods. Six phenotypically resistant pepper accessions were selected based on three years of disease data. All of the selected resistant pepper accessions outperformed the control resistant pepper in terms of resistance (PI 594,137). A genome-wide association study (GWAS) was carried out to identify single nucleotide polymorphisms (SNPs) associated with anthracnose resistance. An association analysis was performed using 53,518 SNPs and the disease score of the 2020 field and in vitro experiment results. Both field and in vitro experiments revealed 25 and 32 significantly associated SNPs, respectively. These SNPs were found on all chromosomes except Ch06 and Ch07 in the field experiment, whereas in the in vitro experiment they were found on all chromosomes except Ch04 and Ch11.

CONCLUSION

In this study, six resistant C. chinense accessions were selected. Additionally, in this study, significantly associated SNPs were found in a gene that codes for a protein kinase receptor, such as serine/threonine-protein kinase, and other genes that are known to be involved in disease resistance. This may strengthen the role of these genes in the development of anthracnose resistance in Capsicum spp. As a result, the SNPs discovered to be strongly linked in this study can be used to identify a potential marker for selecting pepper material resistant to anthracnose, which will assist in the development of resistant varieties.

Signs of local adaptation by genetic selection and isolation promoted by extreme temperature and salinity in the Mediterranean seagrass Posidonia oceanica

Adaptation to local conditions is known to occur in seagrasses; however, knowledge of the genetic basis underlying this phenomenon remains scarce. Here, we analysed Posidonia oceanica from six sites within and around the Stagnone di Marsala, a semi-enclosed coastal lagoon where salinity and temperature exceed the generally described tolerance thresholds of the species. Sea surface temperatures (SSTs) were measured and plant samples were collected for the assessment of morphology, flowering rate and for screening genome-wide polymorphisms using double digest restriction-site-associated DNA sequencing. Results demonstrated more extreme SSTs and salinity levels inside the lagoon than the outer lagoon regions. Morphological results showed significantly fewer and shorter leaves and reduced rhizome growth of P. oceanica from the inner lagoon and past flowering events were recorded only for a meadow farthest away from the lagoon. Using an array of 51,329 single nucleotide polymorphisms, we revealed a clear genetic structure among the study sites and confirmed the genetic isolation and high clonality of the innermost site. In all, 14 outlier loci were identified and annotated with several proteins including those relate to plant stress response, protein transport and regulators of plant-specific developmental events. Especially, five outlier loci showed maximum allele frequency at the innermost site, likely reflecting adaptation to the extreme temperature and salinity regimes, possibly due to the selection of more resistant genotypes and the progressive restriction of gene flow. Overall, this study helps us to disentangle the genetic basis of seagrass adaptation to local environmental conditions and may support future works on assisted evolution in seagrasses.

Genome-Wide Analysis of Proline-Rich Extensin-Like Receptor Kinases (PERKs) Gene Family Reveals Their Roles in Plant Development and Stress Conditions in Oryza sativa L

Proline-rich extensin-like receptor kinases (PERKs) play a crucial role in a wide range of biological processes in plants. In model plants like Arabidopsis, the PERK gene family has been well investigated. Conversely, no information available on the PERK gene family and their biological functions largely remained unknown in rice. This study analyzed the basic physicochemical properties, phylogeny, gene structure, cis-acting elements, Gene ontology (GO) annotation and protein-protein interaction of OsPERK gene family members using various bioinformatics tools based on the whole-genome data of O. sativa. Thus, in this work, 8 PERK genes in rice were identified, and their roles in plant development, growth, and response to various stresses were studied. A phylogenetic study revealed that OsPERKs are grouped into seven classes. Chromosomal mapping also displayed that 8 PERK genes were unevenly distributed on 12 chromosomes. Further, the prediction of subcellular localization indicated that OsPERKs were mainly located at the endomembrane system. Gene structure analysis of OsPERKs has shown a distinctive evolutionary path. In addition, synteny analysis exhibited the 40 orthologous gene pairs in Arabidopsis thaliana, Triticum aestivum, Hordeum vulgare and Medicago truncatula. Furthermore, Ka to Ks proportion shows that most OsPERK genes experienced resilient purifying selection during evolutionary processes. The OsPERK promoters contained several cis-acting regulatory, which are crucial for plant development processes, phytohormone signaling, stress, and defense response. Moreover, the expression pattern of OsPERK family members showed differential expression patterns in different tissues and various stress conditions. Taken together, these results provide clear messages for a better understanding the roles of OsPERK genes in various development stages, tissues, and multifactorial stress as well as enriched the related research of OsPERK family members in rice.

Phase-specific transcriptional patterns of the oomycete pathogen Phytophthora sojae unravel genes essential for asexual development and pathogenic processes

Oomycetes are filamentous microorganisms easily mistaken as fungi but vastly differ in physiology, biochemistry, and genetics. This commonly-held misconception lead to a reduced effectiveness by using conventional fungicides to control oomycetes, thus it demands the identification of novel functional genes as target for precisely design oomycetes-specific microbicide. The present study initially analyzed the available transcriptome data of the model oomycete pathogen, Phytophthora sojae, and constructed an expression matrix of 10,953 genes across the stages of asexual development and host infection. Hierarchical clustering, specificity, and diversity analyses revealed a more pronounced transcriptional plasticity during the stages of asexual development than that in host infection, which drew our attention by particularly focusing on transcripts in asexual development stage to eventually clustered them into 6 phase-specific expression modules. Three of which respectively possessing a serine/threonine phosphatase (PP2C) expressed during the mycelial and sporangium stages, a histidine kinase (HK) expressed during the zoospore and cyst stages, and a bZIP transcription factor (bZIP32) exclusive to the cyst germination stage were selected for down-stream functional validation. In this way, we demonstrated that PP2C, HK, and bZIP32 play significant roles in P. sojae asexual development and virulence. Thus, these findings provide a foundation for further gene functional annotation in oomycetes and crop disease management.

β-Cyclocitral-Mediated Metabolic Changes Optimize Growth and Defense Responses in Solanum lycopersicum L

β-cyclocitral (βCC) is one of the significant oxidative products of β-carotene. It primes plants for multiple stress acclimation without compromising plant growth. Metabolic reorganization is necessary to maintain a balance between growth and defense. However, the βCC-mediated changes in a plant’s metabolic network are unknown. Here, we demonstrate how βCC-induced metabolic changes enable Solanum lycopersicum L. (tomato) plants to promote defense and maintain growth under stress. An analysis of early (0–240 min) and late (72 h) changes in the tomato metabolome after βCC-treatment using liquid chromatography and tandem mass spectrometry identified 57 compounds. A principal coordinate analysis suggested that βCC treatment significantly changes the metabolite profile. A variable importance in projection (VIP) analysis revealed 16 and 19 discriminant metabolites from early and late samples, respectively (VIP ≥ 1.0). Upregulated metabolites were mainly amino acids and phytophenols. Pathway enrichment analysis showed that βCC treatment influenced amino acid metabolism at early and later times; however, phenylpropanoid and isoquinoline biosynthesis were influenced only at the later time. A 66.6% similarity in the upregulated metabolites of βCC- and simulated-herbivory-treated plants confirmed βCC’s role against herbivores. We conclude that βCC steers a temporal separation in amino acids and defense metabolite accumulation that optimizes resource allocation to growth and defense.

MAPKKKs in Plants: Multidimensional Regulators of Plant Growth and Stress Responses

Mitogen-activated protein kinase kinase kinase (MAPKKK, MAP3K) is located upstream of the mitogen-activated protein kinase (MAPK) cascade pathway and is responsible for receiving and transmitting external signals to the downstream MAPKKs. Although a large number of MAP3K genes play important roles in plant growth and development, and response to abiotic and biotic stresses, only a few members' functions and cascade signaling pathways have been clarified, and the downstream MAPKKs and MAPKs of most MAP3Ks are still unknown. As more and more signaling pathways are discovered, the function and regulatory mechanism of MAP3K genes will become clearer. In this paper, the MAP3K genes in plants were classified and the members and basic characteristics of each subfamily of MAP3K were briefly described. Moreover, the roles of plant MAP3Ks in regulating plant growth and development and stress (abiotic and biotic) responses are described in detail. In addition, the roles of MAP3Ks involved in plant hormones signal transduction pathway were briefly introduced, and the future research focus was prospected.

Whole-genome re-sequencing, diversity analysis, and stress-resistance analysis of 77 grape rootstock genotypes

INTRODUCTION

Grape rootstocks play critical role in the development of the grape industry over the globe for their higher adaptability to various environments, and the evaluation of their genetic diversity among grape genotypes is necessary to the conservation and utility of genotypes.

METHODS

To analyze the genetic diversity of grape rootstocks for a better understanding multiple resistance traits, whole-genome re-sequencing of 77 common grape rootstock germplasms was conducted in the present study.

RESULTS

About 645 billion genome sequencing data were generated from the 77 grape rootstocks at an average depth of ~15.5×, based on which the phylogenic clusters were generated and the domestication of grapevine rootstocks was explored. The results indicated that the 77 rootstocks originated from five ancestral components. Through phylogenetic, principal components, and identity-by-descent (IBD) analyses, these 77 grape rootstocks were assembled into ten groups. It is noticed that the wild resources of V. amurensis and V. davidii, originating from China and being generally considered to have stronger resistance against biotic and abiotic stresses, were sub-divided from the other populations. Further analysis indicated that a high level of linkage disequilibrium was found among the 77 rootstock genotypes, and a total of 2,805,889 single nucleotide polymorphisms (SNPs) were excavated, GWAS analysis among the grape rootstocks located 631, 13, 9, 2, 810, and 44 SNP loci that were responsible to resistances to phylloxera, root-knot nematodes, salt, drought, cold and waterlogging traits.

DISCUSSION

This study generated a significant amount of genomic data from grape rootstocks, thus providing a theoretical basis for further research on the resistance mechanism of grape rootstocks and the breeding of resistant varieties. These findings also reveal that China originated V. amurensis and V. davidii could broaden the genetic background of grapevine rootstocks and be important germplasm used in breeding high stress-resistant grapevine rootstocks.

Genome wide association analysis for grain micronutrients and anti-nutritional traits in mungbean [Vigna radiata (L.) R. Wilczek] using SNP markers

Mungbean is an important food grain legume for human nutrition and nutritional food due to its nutrient-dense seed, liked palatability, and high digestibility. However, anti-nutritional factors pose a significant risk to improving nutritional quality for bio-fortification. In the present study, genetic architecture of grain micronutrients (grain iron and zinc concentration) and anti-nutritional factors (grain phytic acid and tannin content) in association mapping panel of 145 diverse mungbean were evaluated. Based on all four parameters genotypes PUSA 1333 and IPM 02-19 were observed as desired genotypes as they had high grain iron and zinc concentration but low grain phytic acid and tannin content. The next generation sequencing (NGS)-based genotyping by sequencing (GBS) identified 14,447 genome-wide SNPs in a diverse selected panel of 127 mungbean genotypes. Population admixture analysis revealed the presence of four different ancestries among the genotypes and LD decay of ∼57.6 kb kb physical distance was noted in mungbean chromosomes. Association mapping analysis revealed that a total of 20 significant SNPs were shared by both GLM and Blink models associated with grain micronutrient and anti-nutritional factor traits, with Blink model identifying 35 putative SNPs. Further, this study identified the 185 putative candidate genes. Including potential candidate genes Vradi07g30190, Vradi01g09630, and Vradi09g05450 were found to be associated with grain iron concentration, Vradi10g04830 with grain zinc concentration, Vradi08g09870 and Vradi01g11110 with grain phytic acid content and Vradi04g11580 and Vradi06g15090 with grain tannin content. Moreover, two genes Vradi07g15310 and Vradi09g05480 showed significant variation in protein structure between native and mutated versions. The identified SNPs and candidate genes are potential powerful tools to provide the essential information for genetic studies and marker-assisted breeding program for nutritional improvement in mungbean.

Review: Status and prospects of association mapping in grapevine

Thanks to current advances in sequencing technologies, novel bioinformatics tools, and efficient modeling solutions, association mapping has become a widely accepted approach to unravel the link between genotype and phenotype diversity in numerous crops. In grapevine, this strategy has been used in the last decades to understand the genetic basis of traits of agronomic interest (fruit quality, crop yield, biotic and abiotic resistance), of special relevance nowadays to improve crop resilience to cope with future climate scenarios. Genome-wide association studies have identified many putative causative loci for different traits, some of them overlapping well-known causal genes identified by conventional quantitative trait loci studies in biparental progenies, and/or validated by functional approaches. In addition, candidate-gene association studies have been useful to pinpoint the causal mutation underlying phenotypic variation for several traits of high interest in breeding programs (like berry color, seedlessness, and muscat flavor), information that has been used to develop highly informative and useful markers already in use in marker-assisted selection processes. Thus, association mapping has proved to represent a valuable step towards high quality and sustainable grape production. This review summarizes current applications of association mapping in grapevine research and discusses future prospects in view of current viticulture challenges.

Genome-wide identification and characterization of exapted transposable elements in the large genome of sunflower (Helianthus annuus L.)

Transposable elements (TEs) are an important source of genome variability, playing many roles in the evolution of eukaryotic species. Besides well-known phenomena, TEs may undergo the exaptation process and generate the so-called exapted transposable element genes (ETEs). Here we present a genome-wide survey of ETEs in the large genome of sunflower (Helianthus annuus L.), in which the massive amount of TEs, provides a significant source for exaptation. A library of sunflower TEs was used to build TE-specific Hidden Markov Model profiles, to search for all available sunflower gene products. In doing so, 20 016 putative ETEs were identified and further investigated for the characteristics that distinguish TEs from genes, leading to the validation of 3530 ETEs. The analysis of ETEs transcription patterns under different stress conditions showed a differential regulation triggered by treatments mimicking biotic and abiotic stress; furthermore, the distribution of functional domains of differentially regulated ETEs revealed a relevant presence of domains involved in many aspects of cellular functions. A comparative genomic investigation was performed including species representative of Asterids and appropriate outgroups: the bulk of ETEs that resulted were specific to the sunflower, while few ETEs presented orthologues in the genome of all analyzed species, making the hypothesis of a conserved function. This study highlights the crucial role played by exaptation, actively contributing to species evolution.

In silico and in vitro approaches allow the identification of the Prosystemin molecular network

Tomato Prosystemin (ProSys), the precursor of Systemin, a small peptidic hormone, is produced at very low concentration in unchallenged plants, while its expression greatly increases in response to several different stressors triggering an array of defence responses. The molecular mechanisms that underpin such a wide array of defence barriers are not fully understood and are likely correlated with the intrinsically disordered (ID) structure of the protein. ID proteins interact with different protein partners forming complexes involved in the modulation of different biological mechanisms. Here we describe the ProSys-protein network that shed light on the molecular mechanisms underpinning ProSys associated defence responses. Three different approaches were used. In silico prediction resulted in 98 direct interactors, most clustering in phytohormone biosynthesis, transcription factors and signal transduction gene classes. The network shows the central role of ProSys during defence responses, that reflects its role as central hub. In vitro ProSys interactors, identified by Affinity Purification-Mass Spectrometry (AP-MS), revealed over three hundred protein partners, while Bimolecular Fluorescent Complementation (BiFC) experiments validated in vivo some interactors predicted in silico and in vitro. Our results demonstrate that ProSys interacts with several proteins and reveal new key molecular events in the ProSys-dependent defence response of tomato plant.

The Gain-of-Function Mutation, OsSpl26, Positively Regulates Plant Immunity in Rice

Rice spotted-leaf mutants are ideal materials to study the molecular mechanism underlying programmed cell death and disease resistance in plants. LOC_Os07g04820 has previously been identified as the candidate gene responsible for the spotted-leaf phenotype in rice Spotted-leaf 26 (Spl26) mutant. Here, we cloned and validated that LOC_Os07g04820 is the locus controlling the spotted-leaf phenotype of Spl26 by reverse functional complementation and CRISPR/Cas9-mediated knockout of the mutant allele. The recessive wild-type spl26 allele (Oryza sativa spotted-leaf 26, Osspl26) is highly conservative in grass species and encodes a putative G-type lectin S-receptor-like serine/threonine protein kinase with 444 amino acid residuals. OsSPL26 localizes to the plasma membrane and can be detected constitutively in roots, stems, leaves, sheaths and panicles. The single base substitution from T to A at position 293 leads to phenylalanine/tyrosine replacement at position 98 in the encoded protein in the mutant and induces excessive accumulation of H₂O₂, leading to oxidative damage to cells, and finally, formation of the spotted-leaf phenotype in Spl26. The formation of lesions not only affects the growth and development of the plants but also activates the defense response and enhances the resistance to the bacterial blight pathogen, Xanthomonas oryzae pv. oryzae. Our results indicate that the gain-of-function by the mutant allele OsSpl26 positively regulates cell death and immunity in rice.

Development of SSR Markers Linked to Stress Responsive Genes along Tomato Chromosome 3 (Solanum lycopersicum L.)

This study aimed to develop novel SSR markers in tomato. Several BAC clones along chromosome 3 in tomato were selected based on their content. The criteria was the availability of genes, either directly or indirectly related to stress response (drought, salinity, and heat) in tomato. A total of 20 novel in silico SSR markers were developed and 96 important nearby genes were identified. The identified nearby genes represent different tomato genes involved in plant growth and development and biotic and abiotic stress tolerance. The developed SSR markers were assessed using tomato landraces. A total of 29 determinate and semi-determinate local tomato landraces collected from diverse environments were utilized. A total of 33 alleles with mean of 1.65 alleles per locus were scored, showing 100% polymorphic patterns, with a mean of 0.18 polymorphism information content (PIC) values. The mean of observed and expected heterozygosity were 0.19 and 0.24, respectively. The mean value of the Jaccard similarity index was used for clustering the landraces. The developed microsatellite markers showed potential to assess genetic variability among tomato landraces. The genetic distance information reported in this study can be used by breeders in future genetic improvement of tomato for tolerance against diverse stresses.

Identification of Candidate Genes for a Major Quantitative Disease Resistance Locus From Soybean PI 427105B for Resistance to Phytophthora sojae

Phytophthora root and stem rot is a yield-limiting soybean disease caused by the soil-borne oomycete Phytophthora sojae. Although multiple quantitative disease resistance loci (QDRL) have been identified, most explain <10% of the phenotypic variation (PV). The major QDRL explaining up to 45% of the PV were previously identified on chromosome 18 and represent a valuable source of resistance for soybean breeding programs. Resistance alleles from plant introductions 427105B and 427106 significantly increase yield in disease-prone fields and result in no significant yield difference in fields with less to no disease pressure. In this study, high-resolution mapping reduced the QDRL interval to 3.1 cm, and RNA-seq analysis of near-isogenic lines (NILs) varying at QDRL-18 pinpointed a single gene of interest which was downregulated in inoculated NILs carrying the resistant allele compared to inoculated NILs with the susceptible allele. This gene of interest putatively encodes a serine-threonine kinase (STK) related to the AtCR4 family and may be acting as a susceptibility factor, based on the specific increase of jasmonic acid concentration in inoculated NILs. This work facilitates further functional analyses and marker-assisted breeding efforts by prioritizing candidate genes and narrowing the targeted region for introgression.

Physical mapping of QTL associated with agronomic and end-use quality traits in spring wheat under conventional and organic management systems

Using phenotypic data of four biparental spring wheat populations evaluated at multiple environments under two management systems, we discovered 152 QTL and 22 QTL hotspots, of which two QTL accounted for up to 37% and 58% of the phenotypic variance, consistently detected in all environments, and fell within genomic regions harboring known genes. Identification of the physical positions of quantitative trait loci (QTL) would be highly useful for developing functional markers and comparing QTL results across multiple independent studies. The objectives of the present study were to map and characterize QTL associated with nine agronomic and end-use quality traits (tillering ability, plant height, lodging, grain yield, grain protein content, thousand kernel weight, test weight, sedimentation volume, and falling number) in hard red spring wheat recombinant inbred lines (RILs) using the International Wheat Genome Sequencing Consortium (IWGSC) RefSeq v2.0 physical map. We evaluated a total of 698 RILs from four populations derived from crosses involving seven parents at 3-8 conventionally (high N) and organically (low N) managed field environments. Using the phenotypic data combined across all environments per management, and the physical map between 1058 and 6526 markers per population, we identified 152 QTL associated with the nine traits, of which 29 had moderate and 2 with major effects. Forty-nine of the 152 QTL mapped across 22 QTL hotspot regions with each region coincident to 2-6 traits. Some of the QTL hotspots were physically located close to known genes. QSv.dms-1A and QPht.dms-4B.1 individually explained up to 37% and 58% of the variation in sedimentation volume and plant height, respectively, and had very large LOD scores that varied from 19.0 to 35.7 and from 16.7 to 55.9, respectively. We consistently detected both QTL in the combined and all individual environments, laying solid ground for further characterization and possibly for cloning.

Targeted designing functional markers revealed the role of retrotransposon derived miRNAs as mobile epigenetic regulators in adaptation responses of pistachio

We developed novel miRNA-based markers based on salt responsive miRNA sequences to detect polymorphisms in miRNA sequences and locations. The validation of 76 combined miRNA + miRNA and miRNA + ISSR markers in the three extreme pistachio populations led to the identification of three selected markers that could link salt tolerance phenotype to genotype and divided pistachio genotypes and Pistacia species into three clusters. This novel functional marker system, in addition to more efficient performance, has higher polymorphisms than previous miRNA-based marker systems. The functional importance of the target gene of five miRNAs in the structure of the three selected markers in regulation of different genes such as ECA2, ALA10, PFK, PHT1;4, PTR3, KUP2, GRAS, TCP, bHLH, PHD finger, PLATZ and genes involved in developmental, signaling and biosynthetic processes shows that the polymorphism associated with these selected miRNAs can make a significant phenotypic difference between salt sensitive and tolerant pistachio genotypes. The sequencing results of selected bands showed the presence of conserved miRNAs in the structure of the mitochondrial genome. Further notable findings of this study are that the sequences of PCR products of two selected markers were annotated as Gypsy and Copia retrotransposable elements. The transposition of retrotransposons with related miRNAs by increasing the number of miRNA copies and changing their location between nuclear and organellar genomes can affect the regulatory activity of these molecules. These findings show the crucial role of retrotransposon-derived miRNAs as mobile epigenetic regulators between intracellular genomes in regulating salt stress responses as well as creating new and tolerant phenotypes for adaptation to environmental conditions.

A Method for Assaying of Protein Kinase Activity In Vivo and Its Use in Studies of Signal Transduction in Strawberry Fruit Ripening

Strawberry (Fragaria × ananassa) fruit ripening is regulated by a complex of cellular signal transduction networks, in which protein kinases are key components. Here, we report a relatively simple method for assaying protein kinase activity in vivo and specifically its application to study the kinase, FaMPK6, signaling in strawberry fruit. Green fluorescent protein (GFP)-tagged FaMPK6 was transiently expressed in strawberry fruit and after stimuli were applied to the fruit it was precipitated using an anti-GFP antibody. The precipitated kinase activity was measured in vitro using ³²P-ATP and myelin basic protein (MBP) as substrates. We also report that FaMPK6 is not involved in the abscisic acid (ABA) signaling cascade, which is closely associated with FaMPK6 signaling in other plant species. However, methyl jasmonate (MeJA), low temperature, and high salt treatments were all found to activate FaMPK6. Transient manipulation of FaMPK6 expression was observed to cause significant changes in the expression patterns of 2749 genes, of which 264 were associated with MeJA signaling. The data also suggest a role for FaMPK6 in modulating cell wall metabolism during fruit ripening. Taken together, the presented method is powerful and its use will contribute to a profound exploration to the signaling mechanism of strawberry fruit ripening.

Expression studies of stress responsive genes in cotton Gossypium hirsutum L

BACKGROUND

Cotton is the world's richest source of natural fiber. Meanwhile cotton plant is top ranked stress sensitive plant thereby affecting its yield and fiber quality. But, in climate change scenario, fiber yield and quality are being affected due to environmental stresses, especially heat, drought and salinity. Present study is aimed to identify cotton genotype harboring prominently expressed stress responsive genes.

METHODS

Four cotton genotypes (IUB-13, IUB-222, IUB-09 and MM-58) were evaluated under drought and salinity stress for yield traits and expression of different stress responsive genes (GhWRKY3, GhDREB2 and GhRDR6).

RESULTS

Pronounced expression of GhWRKY3, GhDREB2and GhRDR6 was observed in cotton variety IUB-13 in stress condition (drought and salinity) as compared to control followed by IUB-222 which revealed that these genotypes might possess substantial potential to cope with environmental hazards encountered in growing season CONCLUSION: Utilization of cotton genotypes i.e., IUB-13 and IUB-222 in cotton breeding program can be very much fruitful for developing cotton genotypes adoptable to climate change.

Genome Wide Association Mapping of Root Traits in the Andean Genepool of Common Bean (Phaseolus vulgaris L.) Grown With and Without Aluminum Toxicity

Common bean is one of the most important grain legumes for human diets but is produced on marginal lands with unfavorable soil conditions; among which Aluminum (Al) toxicity is a serious and widespread problem. Under low pH, stable forms of Al dissolve into the soil solution and as phytotoxic ions inhibit the growth and function of roots through injury to the root apex. This results in a smaller root system that detrimentally effects yield. The goal of this study was to evaluate 227 genotypes from an Andean diversity panel (ADP) of common bean and determine the level of Al toxicity tolerance and candidate genes for this abiotic stress tolerance through root trait analysis and marker association studies. Plants were grown as seedlings in hydroponic tanks at a pH of 4.5 with a treatment of high Al concentration (50 μM) compared to a control (0 μM). The roots were harvested and scanned to determine average root diameter, root volume, root surface area, number of root links, number of root tips, and total root length. Percent reduction or increase was calculated for each trait by comparing treatments. Genome wide association study (GWAS) was conducted by testing phenotypic data against single nucleotide polymorphism (SNP) marker genotyping data for the panel. Principal components and a kinship matrix were included in the mixed linear model to correct for population structure. Analyses of variance indicated the presence of significant difference between genotypes. The heritability of traits ranged from 0.67 to 0.92 in Al-treated and reached similar values in non-treated plants. GWAS revealed significant associations between root traits and genetic markers on chromosomes Pv01, Pv04, Pv05, Pv06, and Pv11 with some SNPs contributing to more than one trait. Candidate genes near these loci were analyzed to explain the detected association and included an Al activated malate transporter gene and a multidrug and toxic compound extrusion gene. This study showed that polygenic inheritance was critical to aluminum toxicity tolerance in common beans roots. Candidate genes found suggested that exudation of malate and citrate as organic acids would be important for Al tolerance. Possible cross-talk between mechanisms of aluminum tolerance and resistance to other abiotic stresses are discussed.

Identification and characterization of long non-coding RNA and their response against citrus bark cracking viroid infection in Humulus lupulus

Long non-coding RNAs (lncRNAs) are a highly heterogeneous class of non-protein-encoding transcripts that play an essential regulatory role in diverse biological processes, including stress responses. The severe stunting disease caused by Citrus bark cracking viroid (CBCVd) poses a major threat to the production of Humulus lupulus (hop) plants. In this study, we systematically investigate the characteristics of the lncRNAs in hop and their role in CBCVd-infection using RNA-sequencing data. Following a stringent filtration criterion, a total of 3598 putative lncRNAs were identified with a high degree of certainty, of which 19% (684) of the lncRNAs were significantly differentially expressed (DE) in CBCVd-infected hop, which were predicted to be mainly involved in plant-pathogen interactions, kinase cascades, secondary metabolism and phytohormone signal transduction. Besides, several lncRNAs and CBCVd-responsive lncRNAs were identified as the precursor of microRNAs and predicted as endogenous target mimics (eTMs) for hop microRNAs involved in CBCVd-infection.

Molecular Programming of Drought-Challenged Trichoderma harzianum-Bioprimed Rice (Oryza sativa L.)

Trichoderma biopriming enhances rice growth in drought-stressed soils by triggering various plant metabolic pathways related to antioxidative defense, secondary metabolites, and hormonal upregulation. In the present study, transcriptomic analysis of rice cultivar IR64 bioprimed with Trichoderma harzianum under drought stress was carried out in comparison with drought-stressed samples using next-generation sequencing techniques. Out of the 2,506 significant (p < 0.05) differentially expressed genes (DEGs), 337 (15%) were exclusively expressed in drought-stressed plants, 382 (15%) were expressed in T. harzianum-treated drought-stressed plants, and 1,787 (70%) were commonly expressed. Furthermore, comparative analysis of upregulated and downregulated genes under stressed conditions showed that 1,053 genes (42%) were upregulated and 733 genes (29%) were downregulated in T. harzianum-treated drought-stressed rice plants. The genes exclusively expressed in T. harzianum-treated drought-stressed plants were mostly photosynthetic and antioxidative such as plastocyanin, small chain of Rubisco, PSI subunit Q, PSII subunit PSBY, osmoproteins, proline-rich protein, aquaporins, stress-enhanced proteins, and chaperonins. The Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis states that the most enriched pathways were metabolic (38%) followed by pathways involved in the synthesis of secondary metabolites (25%), carbon metabolism (6%), phenyl propanoid (7%), and glutathione metabolism (3%). Some of the genes were selected for validation using real-time PCR which showed consistent expression as RNA-Seq data. Furthermore, to establish host-T. harzianum interaction, transcriptome analysis of Trichoderma was also carried out. The Gene Ontology (GO) analysis of T. harzianum transcriptome suggested that the annotated genes are functionally related to carbohydrate binding module, glycoside hydrolase, GMC oxidoreductase, and trehalase and were mainly upregulated, playing an important role in establishing the mycelia colonization of rice roots and its growth. Overall, it can be concluded that T. harzianum biopriming delays drought stress in rice cultivars by a multitude of molecular programming.

Peach PpSnRK1α interacts with bZIP11 and maintains trehalose balance in plants

The nonreducing disaccharide trehalose is widespread in nature. It plays a very important role in plant growth and development. In plants, trehalose is present in trace amounts. High concentration of trehalose disrupts energy balance and inhibits normal growth and development. Studies have shown that high levels of trehalose and trehalose-6-phosphate (T6P), the metabolic precursor of trehalose, inhibit sucrose non-fermenting-1-related protein kinase1 (SnRK1) activity, which affect plant growth and development. However, the role of SnRK1, the energy balance center, in the regulation of trehalose metabolism in plants is unknown. In this study, exogenous trehalose at higher concentrations inhibited the expression of SnRK1 genes, especially PpSnRK1α in peach (Prunus persica) seedlings. This change in gene expression was dependent on trehalose concentration. Furthermore, overexpression of peach PpSnRK1α in Arabidopsis thaliana significantly promoted trehalase activity, reduced T6P content, and suppressed the trehalose synthesis related genes (TPSs, TPPB) expression, promoted the trehalose metabolism of gene expression (TRE1), in addition the transgenic plants alleviated photosynthetic product distribution imbalance (aboveground and underground parts), and enhanced root growth. Yeast two-hybrid and bimolecular fluorescence assays revealed the interaction between PpSnRK1α and peach basic domain leucine zipper transcription factor 11 (PpbZIP11), a key transcription factor of trehalose metabolism, in the nucleus. To summarize, PpSnRK1α overexpression improved bZIP11 transcriptional activity and regulated trehalose metabolism to protect the plants against trehalose-induced damage. This study preliminarily explained the mechanism of SnRK1 regulating trehalose metabolism balance in plants, which laid a foundation for further understanding of energy metabolism and function of SnRK1 in plants.

Transcriptome profiling and weighted gene co-expression network analysis of early floral development in Aquilegia coerulea

The earliest phases of floral development include a number of crucial processes that lay the foundation for the subsequent morphogenesis of floral organs and success in reproduction. Currently, key transcriptional changes during this developmental window have been characterized in the model species Arabidopsis thaliana, but little is known about how transcriptional dynamics change over the course of these developmental processes in other plant systems. Here, we have conducted the first in-depth transcriptome profiling of early floral development in Aquilegia at four finely dissected developmental stages, with eight biological replicates per stage. Using differential gene expression analysis and weighted gene co-expression network analysis, we identified both crucial genes whose expression changes mark the transitions between developmental stages and hub genes in co-expression modules. Our results support the potential functional conservation of key genes in early floral development that have been identified in other systems, but also reveal a number of previously unknown or overlooked loci that are worthy of further investigation. In addition, our results highlight not only the dynamics of transcriptional regulation during early floral development, but also the potential involvement of the complex, essential networks of small RNA and post-translational regulation to these developmental stages.

Genome-wide survey of sucrose non-fermenting 1-related protein kinase 2 in Rosaceae and expression analysis of PbrSnRK2 in response to ABA stress

Background

The members of the sucrose non-fermenting 1-related protein kinase 2 (SnRK2) family are specific serine/threonine protein kinases in plants that play important roles in stress signal transduction and adaptation. Because of their positive regulatory roles in response to adverse conditions, the genes encoding thes proteins are considered potential candidates for breeding of plants for disease resistance and genetic improvement. However, there is far less information about this kinase family, and the function of these genes has not been explored in Rosaceae.

Results

A genome-wide survey and analysis of the genes encoding members of the SnRK2 family were performed in pear (Pyrus bretschneideri) and seven other Rosaceae species. A total of 71 SnRK2 genes were identified from the eight Rosaceae species and classified into three subgroups based on phylogenetic analysis and structural characteristics. Purifying selection played a crucial role in the evolution of SnRK2 genes, and whole-genome duplication and dispersed duplication were the primary forces underlying the characteristics of the SnRK2 gene family in Rosaceae. Transcriptome data and qRT-PCR assay results revealed that the distribution of PbrSnRK2s was very extensive, including across the roots, leaves, pollen, styles, and flowers, although most of them were mainly expressed in leaves. In addition, under stress conditions, the transcript levels of some of the genes were upregulated in leaves in response to ABA treatment.

Conclusions

This study provides useful information and a theoretical introduction for the study of the evolution, expression, and functions of the SnRK2 gene family in plants.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-020-07201-w.

MAPK Enzymes: a ROS Activated Signaling Sensors Involved in Modulating Heat Stress Response, Tolerance and Grain Stability of Wheat under Heat Stress

Mitogen-activated protein kinase (MAPK) signaling cascade is highly conserved across the species triggering the self-adjustment of the cells by transmitting the external signals to the nucleus. The cascade consists of MAPK kinase kinases (MAPKKKs), MAPK kinases (MAPKKs) and MAPKs. These kinases are functionally interrelated through activation by sequential phosphorylation. MAPK cascade is involved in modulating the tolerance and regulating the growth and developmental processes in plants through transcriptional programming. The cascade has been well characterized in Arabidopsis, Tobacco and rice, but limited information is available in wheat due to complexity of genome. MAPK-based sensors have been reported to be highly specific for the external or intracellular stimuli activating specific TF, stress-associated genes (SAGs) and stress-associated proteins (SAPs) linked with heat-stress tolerance and other biological functions especially size, number and quality of grains. Even, MAPKs have been reported to influence the activity of ATP-binding cassette (ABC) transporter superfamily involved in stabilizing the quality of the grains under adverse conditions. Wheat has also diverse network of MAPKs involved in transcriptional reprogramming upon sensing the terminal HS and in turn protect the plants. Current review mainly focuses on the role of MAPKs as signaling sensor and modulator of defense mechanism for mitigating the effect of heat on plants with focus on wheat. It also indirectly protects the nutrient depletion from the grains under heat stress. MAPKs, lying at pivotal positions, can be utilized for manipulating the heat-stress response (HSR) of wheat to develop plant for future (P4F).

Origin and Diversity of Plant Receptor-Like Kinases

Because of their high level of diversity and complex evolutionary histories, most studies on plant receptor-like kinase subfamilies have focused on their kinase domains. With the large amount of genome sequence data available today, particularly on basal land plants and Charophyta, more attention should be paid to primary events that shaped the diversity of the RLK gene family. We thus focus on the motifs and domains found in association with kinase domains to illustrate their origin, organization, and evolutionary dynamics. We discuss when these different domain associations first occurred and how they evolved, based on a literature review complemented by some of our unpublished results.

Genome-wide analysis of small RNAs from Odontoglossum ringspot virus and Cymbidium mosaic virus synergistically infecting Phalaenopsis

Cymbidium mosaic virus (CymMV) and Odontoglossum ringspot virus (ORSV) are the two most prevalent viruses infecting orchids and causing economic losses worldwide. Mixed infection of CymMV and ORSV could induce intensified symptoms as early at 10 days post-inoculation in inoculated Phalaenopsis amabilis, where CymMV pathogenesis was unilaterally enhanced by ORSV. To reveal the antiviral RNA silencing activity in orchids, we characterized the viral small-interfering RNAs (vsiRNAs) from CymMV and ORSV singly or synergistically infecting P. amabilis. We also temporally classified the inoculated leaf-tip tissues and noninoculated adjacent tissues as late and early stages of infection, respectively. Regardless of early or late stage with single or double infection, CymMV and ORSV vsiRNAs were predominant in 21- and 22-nt sizes, with excess positive polarity and under-represented 5'-guanine. While CymMV vsiRNAs mainly derived from RNA-dependent RNA polymerase-coding regions, ORSV vsiRNAs encompassed the coat protein gene and 3'-untranslated region, with a specific hotspot residing in the 3'-terminal pseudoknot. With double infection, CymMV vsiRNAs increased more than 5-fold in number with increasing virus titres. Most vsiRNA features remained unchanged with double inoculation, but additional ORSV vsiRNA hotspot peaks were prominent. The potential vsiRNA-mediated regulation of the novel targets in double-infected tissues thereby provides a different view of CymMV and ORSV synergism. Hence, temporally profiled vsiRNAs from taxonomically distinct CymMV and ORSV illustrate active antiviral RNA silencing in their natural host, Phalaenopsis, during both early and late stages of infection. Our findings provide insights into offence-defence interactions among CymMV, ORSV and orchids.

Genome-wide expression of low temperature response genes in Rosa hybrida L

Plants respond to low temperature stress during cold acclimation, a complex process involving changes in physiological and biochemical modifications. The rose serves as a good model to investigate low temperature responses in perennial ornamentals. In this study, a heterologous apple microarray is used to investigate genome-wide expression profiles in Rosa hybrida subjected to low temperature dark treatment. Transcriptome profiles are determined in floral buds at 0h, 2h, and 12h of low temperature treatment (4 °C). It is observed that a total of 134 transcripts are up-regulated and 169 transcripts are down-regulated in response to low temperature. Interestingly, a total of eight up-regulated genes, including those coding for two cytochrome P450 proteins, two ankyrin repeat family proteins, two metal ion binding proteins, and two zinc finger protein-related transcription factors, along with a single down-regulated gene, coding for a dynamin-like protein, are detected. Transcript profiles of 12 genes known to be involved in cold stress response are also validated using qRT-PCR. Furthermore, expression patterns of the AP2/ERF gene family of transcription factors are investigated in both floral buds and leaves. Overall, AP2/ERFs genes are more rapidly induced in leaves than in floral buds. Moreover, differential expression of several AP2/ERF genes are detected earlier in vegetative rather than in reproductive tissues. These findings highlight important roles of various low temperature response genes in mediating cold acclimation, thereby allowing roses to adapt to low temperatures, but without adversely affecting flower bud development and subsequent flowering, while vegetative tissues undergo early adaptation to low temperatures.

Phosphorylation of ADP-Glucose Pyrophosphorylase During Wheat Seeds Development

Starch is the dominant reserve polysaccharide accumulated in the seed of grasses (like wheat). It is the most common carbohydrate in the human diet and a material applied to the bioplastics and biofuels industry. Hence, the complete understanding of starch metabolism is critical to design rational strategies to improve its allocation in plant reserve tissues. ADP-glucose pyrophosphorylase (ADP-Glc PPase) catalyzes the key (regulated) step in the synthetic starch pathway. The enzyme comprises a small (S) and a large (L) subunit forming an S₂L₂ heterotetramer, which is allosterically regulated by orthophosphate, fructose-6P, and 3P-glycerate. ADP-Glc PPase was found in a phosphorylated state in extracts from wheat seeds. The amount of the phosphorylated protein increased along with the development of the seed and correlated with relative increases of the enzyme activity and starch content. Conversely, this post-translational modification was absent in seeds from Ricinus communis. In vitro, the recombinant ADP-Glc PPase from wheat endosperm was phosphorylated by wheat seed extracts as well as by recombinant Ca²⁺-dependent plant protein kinases. Further analysis showed that the preferential phosphorylation takes place on the L subunit. Results suggest that the ADP-Glc PPase is a phosphorylation target in seeds from grasses but not from oleaginous plants. Accompanying seed maturation and starch accumulation, a combined regulation of ADP-Glc PPase by metabolites and phosphorylation may provide an enzyme with stable levels of activity. Such concerted modulation would drive carbon skeletons to the synthesis of starch for its long-term storage, which later support seed germination.

De novo transcriptome assembly for four species of crustose coralline algae and analysis of unique orthologous genes

Crustose coralline algae (CCA) are calcifying red macroalgae that reef build in their own right and perform essential ecosystem functions on coral reefs worldwide. Despite their importance, limited genetic information exists for this algal group. De novo transcriptomes were compiled for four species of common tropical CCA using RNA-seq. Sequencing generated between 66 and 87 million raw reads. Transcriptomes were assembled, redundant contigs removed, and remaining contigs were annotated using Trinotate. Protein orthology analysis was conducted between CCA species and two noncalcifying red algae species from NCBI that have published genomes and transcriptomes, and 978 orthologous protein groups were found to be uniquely shared amongst CCA. Functional enrichment analysis of these 'CCA-specific' proteins showed a higher than expected number of sequences from categories relating to regulation of biological and cellular processes, such as actin related proteins, heat shock proteins, and adhesion proteins. Some proteins found within these enriched categories, i.e. actin and GH18, have been implicated in calcification in other taxa, and are thus candidates for involvement in CCA calcification. This study provides the first comprehensive investigation of gene content in these species, offering insights not only into the evolution of coralline algae but also of the Rhodophyta more broadly.

Whole-genome resequencing of 472 Vitis accessions for grapevine diversity and demographic history analyses

Understanding the Vitis species at the genomic level is important for cultivar improvement of grapevine. Here we report whole-genome genetic variation at single-base resolution of 472 Vitis accessions, which cover 48 out of 60 extant Vitis species from a wide geographic distribution. The variation helps to identify a recent dramatic expansion and contraction of effective population size in the domesticated grapevines and that cultivars from the pan-Black Sea region have a unique demographic history in comparison to the other domesticated cultivars. We also find selective sweeps for berry edibility and stress resistance improvement. Furthermore, we find associations between candidate genes and important agronomic traits, such as berry shape and aromatic compounds. These results demonstrate resource value of the resequencing data for illuminating the evolutionary biology of Vitis species and providing targets for grapevine genetic improvement.

Despite the importance of grapevine cultivation in human history and the economic values of cultivar improvement, large-scale genomic variation data are lacking. Here the authors resequence 472 Vitis accessions and use the identified genetic variations for domestication history, demography, and GWAS analyses.

On a Cold Night: Transcriptomics of Grapevine Flower Unveils Signal Transduction and Impacted Metabolism

Low temperature is a critical environmental factor limiting plant productivity, especially in northern vineyards. To clarify the impact of this stress on grapevine flower, we used the Vitis array based on Roche-NimbleGen technology to investigate the gene expression of flowers submitted to a cold night. Our objectives were to identify modifications in the transcript levels after stress and during recovery. Consequently, our results confirmed some mechanisms known in grapes or other plants in response to cold stress, notably, (1) the pivotal role of calcium/calmodulin-mediated signaling; (2) the over-expression of sugar transporters and some genes involved in plant defense (especially in carbon metabolism), and (3) the down-regulation of genes encoding galactinol synthase (GOLS), pectate lyases, or polygalacturonases. We also identified some mechanisms not yet known to be involved in the response to cold stress, i.e., (1) the up-regulation of genes encoding G-type lectin S-receptor-like serine threonine-protein kinase, pathogen recognition receptor (PRR5), or heat-shock factors among others; (2) the down-regulation of Myeloblastosis (MYB)-related transcription factors and the Constans-like zinc finger family; and (3) the down-regulation of some genes encoding Pathogen-Related (PR)-proteins. Taken together, our results revealed interesting features and potentially valuable traits associated with stress responses in the grapevine flower. From a long-term perspective, our study provides useful starting points for future investigation.

Circadian Regulation of Alternative Splicing of Drought-Associated CIPK Genes in Dendrobium catenatum (Orchidaceae)

Dendrobium catenatum, an epiphytic and lithophytic species, suffers frequently from perennial shortage of water in the wild. The molecular mechanisms of this orchid's tolerance to abiotic stress, especially drought, remain largely unknown. It is well-known that CBL-interacting protein kinase (CIPKs) proteins play important roles in plant developmental processes, signal transduction, and responses to abiotic stress. To study the CIPKs' functions for D. catenatum, we first identified 24 CIPK genes from it. We divided them into three subgroups, with varying intron numbers and protein motifs, based on phylogeny analysis. Expression patterns of CIPK family genes in different tissues and in response to either drought or cold stresses suggested DcaCIPK11 may be associated with signal transduction and energy metabolism. DcaCIPK9, -14, and -16 are predicted to play critical roles during drought treatment specifically. Furthermore, transcript expression abundances of DcaCIPK16 showed polar opposites during day and night. Whether under drought treatment or not, DcaCIPK16 tended to emphatically express transcript1 during the day and transcript3 at night. This implied that expression of the transcripts might be regulated by circadian rhythm. qRT-PCR analysis also indicated that DcaCIPK3, -8, and -20 were strongly influenced by circadian rhythmicity. In contrast with previous studies, for the first time to our knowledge, our study revealed that the major CIPK gene transcript expressed was not always the same and was affected by the biological clock, providing a different perspective on alternative splicing preference.

Comparison of gene co-networks analysis provide a systems view of rice (Oryza sativa L.) response to Tilletia horrida infection

The biotrophic soil-borne fungus Tilletia horrida causes rice kernel smut, an important disease affecting the production of rice male sterile lines in most hybrid rice growing regions of the world. There are no successful ways of controlling this disease and there has been little study of mechanisms of resistance to T. horrida. Based on transcriptional data of different infection time points, we found 23, 782 and 23, 718 differentially expressed genes (fragments per kilobase of transcript sequence per million, FPKM >1) in Jiangcheng 3A (resistant to T. horrida) and 9311A (susceptible to T. horrida), respectively. In order to illuminate the differential responses of the two rice male sterile lines to T. horrida, we identified gene co-expression modules using the method of weighted gene co-expression network analysis (WGCNA) and compared the different biological functions of gene co-expression networks in key modules at different infection time points. The results indicated that gene co-expression networks in the two rice genotypes were different and that genes contained in some modules of the two groups may play important roles in resistance to T. horrida, such as DTH8 and OsHop/Sti1a. Furthermore, these results provide a global view of the responses of two different phenotypes to T. horrida, and assist our understanding of the regulation of expression changes after T. horrida infection.

Peanut Stunt Virus and Its Satellite RNA Trigger Changes in Phosphorylation in N. benthamiana Infected Plants at the Early Stage of the Infection

Signaling in host plants is an integral part of a successful infection by pathogenic RNA viruses. Therefore, identifying early signaling events in host plants that play an important role in establishing the infection process will help our understanding of the disease process. In this context, phosphorylation constitutes one of the most important post-translational protein modifications, regulating many cellular signaling processes. In this study, we aimed to identify the processes affected by infection with Peanut stunt virus (PSV) and its satellite RNA (satRNA) in Nicotiana benthamiana at the early stage of pathogenesis. To achieve this, we performed proteome and phosphoproteome analyses on plants treated with PSV and its satRNA. The analysis of the number of differentially phosphorylated proteins showed strong down-regulation in phosphorylation in virus-treated plants (without satRNA). Moreover, proteome analysis revealed more down-regulated proteins in PSV and satRNA-treated plants, which indicated a complex dependence between proteins and their modifications. Apart from changes in photosynthesis and carbon metabolism, which are usually observed in virus-infected plants, alterations in proteins involved in RNA synthesis, transport, and turnover were observed. As a whole, this is the first community (phospho)proteome resource upon infection of N. benthamiana with a cucumovirus and its satRNA and this resource constitutes a valuable data set for future studies.

Quantitative Proteomics Analysis of Lettuce (Lactuca sativa L.) Reveals Molecular Basis-Associated Auxin and Photosynthesis with Bolting Induced by High Temperature

Bolting is a key process in the growth and development of lettuce (Lactuca sativa L.). A high temperature can induce early bolting, which decreases both the quality and production of lettuce. However, knowledge of underlying lettuce bolting is still lacking. To better understand the molecular basis of bolting, a comparative proteomics analysis was conducted on lettuce stems, during the bolting period induced by a high temperature (33 °C) and a control temperature (20 °C) using iTRAQ-based proteomics, phenotypic measures, and biological verifications using qRT-PCR and Western blot. The high temperature induced lettuce bolting, while the control temperature did not. Of the 5454 identified proteins, 619 proteins presented differential abundance induced by high-temperature relative to the control group, of which 345 had an increased abundance and 274 had a decreased abundance. Proteins with an abundance level change were mainly enriched in pathways associated with photosynthesis and tryptophan metabolism involved in auxin (IAA) biosynthesis. Moreover, among the proteins with differential abundance, proteins associated with photosynthesis and tryptophan metabolism were increased. These findings indicate that a high temperature enhances the function of photosynthesis and IAA biosynthesis to promote the process of bolting, which is in line with the physiology and transcription level of IAA metabolism. Our data provide a first comprehensive dataset for gaining novel understanding of the molecular basis underlying lettuce bolting induced by high temperature. It is potentially important for further functional analysis and genetic manipulation for molecular breeding to breed new cultivars of lettuce to restrain early bolting, which is vital for improving vegetable quality.

The kinome of pineapple: catalog and insights into functions in crassulacean acid metabolism plants

BACKGROUND

Crassulacean acid metabolism (CAM) plants use water 20-80% more efficiently by shifting stomata opening and primary CO₂ uptake and fixation to the nighttime. Protein kinases (PKs) play pivotal roles in this biological process. However, few PKs have been functionally analyzed precisely due to their abundance and potential functional redundancy (caused by numerous gene duplications).

RESULTS

In this study, we systematically identified a total of 758 predicted PK genes in the genome of a CAM plant, pineapple (Ananas comosus). The pineapple kinome was classified into 20 groups and 116 families based on the kinase domain sequences. The RLK was the largest group, containing 480 members, and over half of them were predicted to locate at the plasma membrane. Both segmental and tandem duplications make important contributions to the expansion of pineapple kinome based on the synteny analysis. Ka/Ks ratios showed all of the duplication events were under purifying selection. The global expression analysis revealed that pineapple PKs exhibit different tissue-specific and diurnal expression patterns. Forty PK genes in a cluster performed higher expression levels in green leaf tip than in white leaf base, and fourteen of them had strong differential expression patterns between the photosynthetic green leaf tip and the non-photosynthetic white leaf base tissues.

CONCLUSIONS

Our findings provide insights into the evolution and biological function of pineapple PKs and a foundation for further functional analysis of PKs in CAM plants. The gene duplication, expression, and coexpression analysis helped us to rapidly identify the key candidates in pineapple kinome, which may play roles in the carbon fixation process in pineapple and help engineering CAM pathway into C3 crops for improved drought tolerance.

Transcriptional analysis identifies major pathways as response components to Sporisorium scitamineum stress in sugarcane

BACKGROUND

Sugarcane smut, which is caused by Sporisorium scitamineum, is a severe fungal disease affecting sugarcane. However, the major pathways involved in the interaction between sugarcane and S. scitamineum remains unclear.

RESULTS

In the present study, suppression subtractive hybridization (SSH) library construction, together with reverse northern blotting, was conducted on the most prevalent sugarcane genotype ROC22 challenged with S. scitamineum. After alignment and homologous expressed sequence tag (EST) assembly, a total of 155 differentially expressed unigenes were identified from SSH libraries. Totally, 26 of 155 differentially expressed unigenes were analyzed by qRT-PCR in sugarcane smut-resistant genotype YC05-179 and susceptible genotype ROC22. Genes encoded two unknown protein (Q1 and Q11), serine/threonine kinase (Q2), fiber protein (Q3), eukaryotic translation initiation factor 5A (Q23), and Sc14-3-3-like protein (Q24) were induced in sugarcane smut-resistant genotype YC05-179 but inhibited in susceptible genotype ROC22. Based on the differential expression data achieved from SSH libraries and qRT-PCR, we found that, serine/threonine kinases, Ca²⁺ sensors, mitogen-activated protein genes and some NBS-LRR genes may involve in the signal recognition and transduction of smut fungus infection in sugarcane. While in the plant hormone signaling pathways, the genes related to auxin, abscisic acid, salicylic acid and ethylene were more apparently in response to smut fungus invasion. The hypersensitive response, protein metabolism, polyamine synthesis, and cell wall formation may play an important role in sugarcane defense against smut fungus colonization. Additionally, the Sc14-3-3 might serve as a molecular modulator in sugarcane being immune to smut disease by interacting with proteins like ScGAPN (Q10), which have been further verified by BiFC assay.

CONCLUSIONS

The findings of the present study could provide a general view about gene pathways involving in sugarcane defense against smut disease and facilitate a better understanding of the molecular mechanism underlying sugarcane-S. scitamineum interaction.

Global analysis of threonine metabolism genes unravel key players in rice to improve the abiotic stress tolerance

The diversity in plant metabolites with improved phytonutrients is essential to achieve global food security and sustainable crop yield. Our study using computational metabolomics genome wide association study (cmGWAS) reports on a comprehensive profiling of threonine (Thr) metabolite in rice. Sixteen abiotic stress responsive (AbSR) – Thr metabolite producing genes (ThrMPG), modulate metabolite levels and play a significant role determining both physiological and nutritional importance of rice. These AbSR-ThrMPG were computationally analysed for their protein properties using OryzaCyc through plant metabolic network analyser. A total of 1373 and 1028 SNPs were involved in complex traits and genomic variations. Comparative mapping of AbSR-ThrMPG revealed the chromosomal colinearity with C4 grass species. Further, computational expression pattern of these genes predicted a differential expression profiling in diverse developmental tissues. Protein interaction of protein coding gene sequences revealed that the abiotic stresses (AbS) are multigenic in nature. In silico expression of AbSR-ThrMPG determined the putative involvement in response to individual AbS. This is the first comprehensive genome wide study reporting on AbSR –ThrMPG analysis in rice. The results of this study provide a pivotal resource for further functional investigation of these key genes in the vital areas of manipulating AbS signaling in rice improvement.

A gene co-expression network model identifies yield-related vicinity networks in Jatropha curcas shoot system

The plant shoot system consists of reproductive organs such as inflorescences, buds and fruits, and the vegetative leaves and stems. In this study, the reproductive part of the Jatropha curcas shoot system, which includes the aerial shoots, shoots bearing the inflorescence and inflorescence were investigated in regard to gene-to-gene interactions underpinning yield-related biological processes. An RNA-seq based sequencing of shoot tissues performed on an Illumina HiSeq. 2500 platform generated 18 transcriptomes. Using the reference genome-based mapping approach, a total of 64 361 genes was identified in all samples and the data was annotated against the non-redundant database by the BLAST2GO Pro. Suite. After removing the outlier genes and samples, a total of 12 734 genes across 17 samples were subjected to gene co-expression network construction using petal, an R library. A gene co-expression network model built with scale-free and small-world properties extracted four vicinity networks (VNs) with putative involvement in yield-related biological processes as follow; heat stress tolerance, floral and shoot meristem differentiation, biosynthesis of chlorophyll molecules and laticifers, cell wall metabolism and epigenetic regulations. Our VNs revealed putative key players that could be adapted in breeding strategies for J. curcas shoot system improvements.

Fine Mapping and Identification of a Novel Phytophthora Root Rot Resistance Locus RpsZS18 on Chromosome 2 in Soybean

Phytophthora root rot (PRR) caused by Phytophthora sojae is a major soybean disease that causes severe economic losses worldwide. Using soybean cultivars carrying a Rps resistance gene is the most effective strategy for controlling this disease. We previously detected a novel Phytophthora resistance gene, RpsZS18, on chromosome 2 of the soybean cultivar Zaoshu18. The aim of the present study was to identify and finely map RpsZS18. We used 232 F2:3 families generated from a cross between Zaoshu18 (resistant) and Williams (susceptible) as the mapping population. Simple sequence repeat (SSR) markers distributed on chromosome 2 were used to map RpsZS18. First, 12 SSR markers linked with RpsZS18 were identified by linkage analyses, including two newly developed SSR markers, ZCSSR33 and ZCSSR46, that flanked the gene at distances of 0.9 and 0.5 cM, respectively. Second, PCR-based InDel markers were developed based on sequence differences between the two parents and used to further narrow down the mapping region of RpsZS18 to 71.3 kb. Third, haplotype analyses were carried out in the RpsZS18 region using 14 soybean genotypes with whole-genome resequencing. We detected six genes with unique haplotype sequences in Zaoshu18. Finally, quantitative real-time PCR assays of the six genes revealed an EF-hand calcium-binding domain containing protein encoding gene (Glyma.02g245700), a pfkB carbohydrate kinase encoding gene (Glyma.02g245800), and a gene with no functional annotation (Glyma.02g246300), are putative candidate PRR resistance genes. This study provides useful information for breeding P. sojae-resistant soybean cultivars.

Genome-Wide Identification and Functional Analyses of the CRK Gene Family in Cotton Reveals GbCRK18 Confers Verticillium Wilt Resistance in Gossypium barbadense

Cysteine-rich receptor-like kinases (CRKs) are a large subfamily of plant receptor-like kinases that play a critical role in disease resistance in plants. However, knowledge about the CRK gene family in cotton and its function against Verticillium wilt (VW), a destructive disease caused by Verticillium dahliae that significantly reduces cotton yields is lacking. In this study, we identified a total of 30 typical CRKs in a Gossypium barbadense genome (GbCRKs). Eleven of these (>30%) are located on the A06 and D06 chromosomes, and 18 consisted of 9 paralogous pairs encoded in the A and D subgenomes. Phylogenetic analysis showed that the GbCRKs could be classified into four broad groups, the expansion of which has probably been driven by tandem duplication. Gene expression profiling of the GbCRKs in resistant and susceptible cotton cultivars revealed that a phylogenetic cluster of nine of the GbCRK genes were up-regulated in response to V. dahliae infection. Virus-induced gene silencing of each of these nine GbCRKs independently revealed that the silencing of GbCRK18 was sufficient to compromise VW resistance in G. barbadense. GbCRK18 expression could be induced by V. dahliae infection or jasmonic acid, and displayed plasma membrane localization. Therefore, our expression analyses indicated that the CRK gene family is differentially regulated in response to Verticillium infection, while gene silencing experiments revealed that GbCRK18 in particular confers VW resistance in G. barbadense.

Competitive Ability of Maize Pollen Grains Requires Paralogous Serine Threonine Protein Kinases STK1 and STK2

serine threonine kinase1 (stk1) and serine threonine kinase2 (stk2) are closely related maize paralogous genes predicted to encode serine/threonine protein kinases. Pollen mutated in stk1 or stk2 competes poorly with normal pollen, pointing to a defect in pollen tube germination or growth. Both genes are expressed in pollen, but not in most other tissues. In germination media, STK1 and STK2 fluorescent fusion proteins localize to the plasma membrane of the vegetative cell. RNA-seq experiments identified 534 differentially expressed genes in stk1 mutant pollen relative to wild type. Gene ontology (GO) molecular functional analysis uncovered several differentially expressed genes with putative ribosome initiation and elongation functions, suggesting that stk1 might affect ribosome function. Of the two paralogs, stk1 may play a more important role in pollen development than stk2, as stk2 mutations have a smaller pollen transmission effect. However, stk2 does act as an enhancer of stk1 because the double mutant combination is only infrequently pollen-transmitted in double heterozygotes. We conclude that the stk paralogs play an essential role in pollen development.