The extracellular EXO protein mediates cell expansion in Arabidopsis leaves

Wake up: the regulation of dormancy release and bud break in perennial plants

In order to survive harsh winter conditions, perennial trees in the temperate and frigid regions enter a dormant state and cease growth in late summer after vigorous growth in spring and summer. After experiencing prolonged cold temperature and short days in winter, trees release their dormancy, and they resume growth to produce new buds in the following spring, a process known as bud break. The establishment/release of bud dormancy and bud break are crucial for the adaptations of woody plants and their survival in the natural environment. Photoperiod and temperature are key regulators in the bud dormancy and break cycle. In recent years, significant progress has been made in understanding the molecular mechanism for how photoperiod and temperature regulate seasonal growth and dormancy. Here, we summarized the regulatory network and mechanisms underlying the seasonal growth of perennial woody plants in the temperate and frigid regions, focusing on several molecular modules including the photoperiod, circadian clock, EARLY BUD BREAK 1 (EBB1) - SHORT VEGETATIVE PHASE Like (SVL) - EARLY BUD BREAK 3 (EBB3) module and hormone regulation. Through these modules, we will summarize how perennial trees release dormancy and bud break in order to better understand their differences and connections. By elucidating the interactions among these factors, we also point out the questions and challenges need to be addressed in understanding the bud dormancy and break cycle of perennial plants.

Analysis of the CRK expressions in bottle gourd (Lagenaria siceraria) under Fusarium oxysporum f. sp. lagenariae stress by using genome-wide identification strategy

BACKGROUND

The cysteine-rich receptor-like kinases (CRKs) family in plants have been reported to perform multiple functions against various stresses. However, the CRK family in bottle gourd (Lagenaria siceraria) has not been well-explored. Herein, a comprehensive genome-wide identification and characterization of the CRK gene family has been carried out in bottle gourd under Genome-wide characterization of CRK genes in bottle gourds under Fusarium oxysporum f. sp. lagenariae infection.

RESULTS

A stringent set of bioinformatic analyses identified 18 LsCRKs in the bottle gourd genome. Chromosomal mapping of the identified LsCRKs revealed that the LsCRKs were distributed in 4 chromosomes in the bottle gourd genome. The phylogenetic analysis of LsCRKs divided them into two subgroups on the tree. The synteny and collinearity analysis of the LsCRKs among themselves and other plant CRKs provided insights into their conservancy and expansion. Gene ontology analysis of the identified LsCRKs suggested their possible roles in regulating different physiological processes and stress responses in bottle gourd. To assess the involvement of the LsCRKs under F. oxysporum f. sp. lagenariae infection, bottle gourd seedlings were transplanted into the pots with F. oxysporum-infected soil. The expression analysis revealed that multiple LsCRKs exhibited induced expression, suggesting their involvement in bottle gourd-F. oxysporum interactions. Additionally, the protein-protein interaction analysis suggested some important interacting partners of LsCRKs crucial to different physiological processes in bottle gourd.

CONCLUSIONS

The present work explored and analyzed the LsCRKs in bottle gourd. Functional predictions and interaction network analysis suggested the roles of LsCRKs in modulating multiple physiological processes in bottle gourd. The expression dynamics of LsCKRs under fungal pathogen infection suggest their involvement in stress response in bottle gourds. Overall, the results of the study provide basic information about the CRK family in bottle gourds and their involvement in fungal pathogen response.

Time-resolved oxidative signal convergence across the algae-embryophyte divide

The earliest land plants faced a significant challenge in adapting to environmental stressors. Stress on land is unique in its dynamics, entailing swift and drastic changes in light and temperature. While we know that land plants share with their closest streptophyte algal relatives key components of the genetic makeup for dynamic stress responses, their concerted action is little understood. Here, we combine time-course stress profiling using photophysiology, transcriptomics on 2.7 Tbp of data, and metabolite profiling analyses on 270 distinct samples, to study stress kinetics across three 600-million-year-divergent streptophytes. Through co-expression analysis and Granger causal inference we predict a gene regulatory network that retraces a web of ancient signal convergences at ethylene signaling components, osmosensors, and chains of major kinases. These kinase hubs already integrated diverse environmental inputs since before the dawn of plants on land.

Maize COMPACT PLANT 3 regulates plant architecture and facilitates high-density planting

Compact plant architecture allows more efficient light capture under higher planting density. Thus, it is a crucial strategy for improving crop yield, particularly in maize (Zea mays L.). Here, we isolated a maize gene, COMPACT PLANT 3 (CT3), regulating plant architecture, using map-based cloning. CT3, encoding a GRAS protein, interacts with an AP2 transcription factor (TF), DWARF AND IRREGULAR LEAF 1 (DIL1). The genetic analysis showed that CT3 and DIL1 regulate leaf angle and plant height via the same pathway, supporting the biological role of their interaction by forming a complex. Transcriptome and DNA profiling analyses revealed that these 2 TFs share many common target genes. We further observed that CT3 functions as a co-regulator to enhance the DNA affinity and transcriptional activity of DIL1. This finding was further supported by the direct binding of DIL1 to 2 cell wall-related genes, ZmEXO1 and ZmXTH14, which were downregulated in the ct3 mutant. Furthermore, ZmEXO1 regulated plant architecture in a manner similar to CT3- and DIL1-mediated regulation. Zmexo1, ct3, and dil1 mutants showed defective cell wall integrity and had reduced cell wall-related components. The introduction of the ct3 or dil1 mutant allele into elite maize hybrids led to a more compact architecture and increased yield under high planting density. Our findings reveal a regulatory pathway of maize plant architecture and provide targets to increase yield under high planting density.

Mixed DAMP/MAMP oligosaccharides promote both growth and defense against fungal pathogens of cucumber

Plants recognize a variety of environmental molecules, thereby triggering appropriate responses to biotic or abiotic stresses. Substances containing microbes-associated molecular patterns (MAMPs) and damage-associated molecular patterns (DAMPs) are representative inducers of pathogen resistance and damage repair, thus treatment of healthy plants with such substances can pre-activate plant immunity and cell repair functions. In this study, the effects of DAMP/MAMP oligosaccharides mixture (Oligo-Mix) derived from plant cell wall (cello-oligosaccharide and xylo-oligosaccharide), and fungal cell wall (chitin-oligosaccharide) were examined in cucumber. Treatment of cucumber with Oligo-Mix promoted root germination and plant growth, along with increased chlorophyll contents in the leaves. Oligo-Mix treatment also induced typical defense responses such as MAP kinase activation and callose deposition in leaves. Pretreatment of Oligo-Mix enhanced disease resistance of cucumber leaves against pathogenic fungi Podosphaera xanthii (powdery mildew) and Colletotrichum orbiculare (anthracnose). Oligo-Mix treatment increased the induction of hypersensitive cell death around the infection site of pathogens, which inhibited further infection and the conidial formation of pathogens on the cucumber leaves. RNA-seq analysis revealed that Oligo-Mix treatment upregulated genes associated with plant structural reinforcement, responses to abiotic stresses and plant defense. These results suggested that Oligo-Mix has beneficial effects on growth and disease resistance in cucumber, making it a promising biostimulant for agricultural application.

Extracellular vesicles of Norway spruce contain precursors and enzymes for lignin formation and salicylic acid

Lignin is a phenolic polymer in plants that rigidifies the cell walls of water-conducting tracheary elements and support-providing fibers and stone cells. Different mechanisms have been suggested for the transport of lignin precursors to the site of lignification in the cell wall. Extracellular vesicle (EV)-enriched samples isolated from a lignin-forming cell suspension culture of Norway spruce (Picea abies L. Karst.) contained both phenolic metabolites and enzymes related to lignin biosynthesis. Metabolomic analysis revealed mono-, di-, and oligolignols in the EV isolates, as well as carbohydrates and amino acids. In addition, salicylic acid (SA) and some proteins involved in SA signaling were detected in the EV-enriched samples. A proteomic analysis detected several laccases, peroxidases, β-glucosidases, putative dirigent proteins, and cell wall-modifying enzymes, such as glycosyl hydrolases, transglucosylase/hydrolases, and expansins in EVs. Our findings suggest that EVs are involved in transporting enzymes required for lignin polymerization in Norway spruce, and radical coupling of monolignols can occur in these vesicles.

GhEXL3 participates in brassinosteroids regulation of fiber elongation in Gossypium hirsutum

Cotton fiber (Gossypium hirsutum) serves as an ideal model for investigating the molecular mechanisms of plant cell elongation at the single-cell level. Brassinosteroids (BRs) play a crucial role in regulating plant growth and development. However, the mechanism by which BR influences cotton fiber elongation remains incompletely understood. In this study, we identified EXORDIUM-like (GhEXL3) through transcriptome analysis of fibers from BR-deficient cotton mutant pagoda 1 (pag1) and BRI1-EMS-SUPPRESSOR 1 (GhBES1.4, encoding a central transcription factor of BR signaling) overexpression cotton lines. Knockout of GhEXL3 using CRISPR/Cas9 was found to impede cotton fiber elongation, while its overexpression promoted fiber elongation, suggesting a positive regulatory function for GhEXL3 in fiber elongation. Furthermore, in vitro ovule culture experiments revealed that the overexpression of GhEXL3 partially counteracted the inhibitory effects of brassinazole (BRZ) on cotton fiber elongation, providing additional evidence of GhEXL3 involvement in BR signaling pathways. Moreover, our findings demonstrate that GhBES1.4 directly binds to the E-box (CACGTG) motif in the GhEXL3 promoter region and enhances its transcription. RNA-seq analysis revealed that overexpression of GhEXL3 upregulated the expression of EXPs, XTHs, and other genes associated with fiber cell elongation. Overall, our study contributes to understanding the mechanism by which BR regulates the elongation of cotton fibers through the direct modulation of GhEXL3 expression by GhBES1.4.

Multilayer regulatory landscape and new regulators identification for bud dormancy release and bud break in Populus

For trees originating from boreal and temperate regions, the dormancy-to-active transition, also known as bud dormancy release and bud break, are crucial processes that allow trees to reactive growth in the spring. The molecular mechanisms underlying these two processes remain poorly understood. Here, through integrative multiomics analysis of the transcriptome, DNA methylome, and proteome, we gained insights into the reprogrammed cellular processes associated with bud dormancy release and bud break. Our findings revealed multilayer regulatory landscapes governing bud dormancy release and bud break regulation, providing a valuable reference framework for future functional studies. Based on the multiomics analysis, we have determined a novel long intergenic noncoding RNA named Phenology Responsive Intergenic lncRNA 1 (PRIR1) plays a role in the activation of bud break. that the molecular mechanism of PRIR1 has been preliminary explored, and it may partially promote bud break by activating its neighbouring gene, EXORDIUM LIKE 5 (PtEXL5), which has also been genetically confirmed as an activator for bud break. This study has revealed a lncRNA-mediated regulatory mechanism for the control of bud break in Populus, operating independently of known regulatory pathways.

Using RNA sequencing to unravel molecular changes underlying the defense response in chickpea induced by Phytophthora medicaginis

Phytophthora root rot (PRR), caused by Phytophthora medicaginis, is a major soil-borne disease of chickpea in Australia. Breeding for PRR resistance is an effective approach to avoid significant yield loss. Genetic resistance has been identified in cultivated chickpea (Cicer arietinum) and in the wild relative C. echinospermum, with previous studies identifying independent genetic loci associated with each of these sources. However, the molecular mechanisms associated with PRR resistance are not known. RNA sequencing analysis employed in this study identified changes in gene expression in roots of three chickpea genotypes grown hydroponically, early post-infection with P. medicaginis zoospores. Analyses of differentially expressed genes (DEG) identified the activation of a higher number of non-specific R-genes in a PRR-susceptible variety than in the resistant genotypes, suggesting a whole plant resistance response occurring in chickpea against the pathogen. Contrasting molecular changes in signaling profiles, proteolysis and transcription factor pathways were observed in the cultivated and wild Cicer-derived resistant genotypes. DEG patterns supported a hypothesis that increased root elongation and reduced adventitious root formation limit the pathogen entry points in the genotype containing the wild Cicer source of PRR resistance. Candidate resistance genes, including an aquaporin and a maltose transporter in the wild Cicer source and GDSL esterases/lipases in the cultivated source of resistance, were oppositely regulated. Increased knowledge of these genes and pathways will improve our understanding of molecular mechanisms controlling PRR resistance in chickpea, and support the development of elite chickpea varieties through molecular breeding approaches.

Comparative transcriptomic analysis of the nodulation-competent zone and inference of transcription regulatory network in silicon applied Glycine max [L.]-Merr. Roots

KEY MESSAGE

The study unveils Si's regulatory influence by regulating DEGs, TFs, and TRs. Further bHLH subfamily and auxin transporter pathway elucidates the mechanisms enhancing root development and nodulation. Soybean is a globally important crop serving as a primary source of vegetable protein for millions of individuals. The roots of these plants harbour essential nitrogen fixing structures called nodules. This study investigates the multifaceted impact of silicon (Si) application on soybean, with a focus on root development, and nodulation employing comprehensive transcriptomic analyses and gene regulatory network. RNA sequence analysis was utilised to examine the change in gene expression and identify the noteworthy differentially expressed genes (DEGs) linked to the enhancement of soybean root nodulation and root development. A set of 316 genes involved in diverse biological and molecular pathways are identified, with emphasis on transcription factors (TFs) and transcriptional regulators (TRs). The study uncovers TF and TR genes, categorized into 68 distinct families, highlighting the intricate regulatory landscape influenced by Si in soybeans. Upregulated most important bHLH subfamily and the involvement of the auxin transporter pathway underscore the molecular mechanisms contributing to enhanced root development and nodulation. The study bridges insights from other research, reinforcing Si's impact on stress-response pathways and phenylpropanoid biosynthesis crucial for nodulation. The study reveals significant alterations in gene expression patterns associated with cellular component functions, root development, and nodulation in response to Si.

PagEXPA1 combines with PagCDKB2;1 to regulate plant growth and the elongation of fibers in Populus alba × Populus glandulosa

Expansins are important plant cell wall proteins. They can loosen and soften the cell walls and lead to wall extension and cell expansion. To investigate their role in wood formation and fiber elongation, the PagEXPA1 that highly expressed in cell differentiation and expansion tissues was cloned from 84K poplar (Populus alba × P. glandulosa). The subcellular localization showed that PagEXPA1 located in the cell wall and it was highly expressed in primary stems and young leaves. Compared with non-transgenic 84K poplar, overexpression of PagEXPA1 can promote plant-growth, lignification, and fiber cell elongation, while PagEXPA1 Cas9-editing mutant lines exhibited the opposite phenotype. Transcriptome analysis revealed that DEGs were mainly enriched in some important processes, which are associated with cell wall formation and cellulose synthesis. The protein interaction prediction and expression analysis showed that PagCDKB2:1 and PagEXPA1 might have an interaction relationship. The luciferase complementary assay and bimolecular fluorescence complementary assay validated that PagEXPA1 can combined with PagCDKB2;1. So they promoted the expansion of xylem vascular tissues and the development of poplar though participating in the regulation of cell division and differentiation by programming the cell-cycle. It provides good foundation for molecular breeding of fast-growing and high-quality poplar varieties.

CsEXL3 regulate mechanical harvest-related droopy leaves under the transcriptional activation of CsBES1.2 in tea plant

Due to a labor shortage, the mechanical harvesting of tea plantations has become a focal point. However, mechanical harvest efficiency was hampered by droopy leaves, leading to a high rate of broken tea shoots and leaves. Here, we dissected the genetic structure of leaf droopiness in tea plants using genome-wide association studies (GWAS) on 146 accessions, combined with transcriptome from two accessions with contrasting droopy leaf phenotypes. A set of 16 quantitative trait loci (QTLs) containing 54 SNPs and 34 corresponding candidate genes associated with droopiness were then identified. Among these, CsEXL3 (EXORDIUM-LIKE 3) from Chromosome 1 emerged as a candidate gene. Further investigations revealed that silencing CsEXL3 in tea plants resulted in weaker vascular cell malformation and brassinosteroid-induced leaf droopiness. Additionally, brassinosteroid signal factor CsBES1.2 was proved to participate in CsEXL3-induced droopiness and vascular cell malformation via using the CsBES1.2-silencing tea plant. Notably, CsBES1.2 bound on the E-box of CsEXL3 promoter to transcriptionally activate CsEXL3 expression as CUT&TAG based ChIP-qPCR and ChIP-seq suggested in vivo as well as EMSA and Y1H indicated in vitro. Furthermore, CsEXL3 instead of CsBES1.2 decreased lignin content and the expressing levels of lignin biosynthesis genes. Overall, our findings suggest that CsEXL3 regulates droopy leaves, partially through the transcriptional activation of CsBES1.2, with the potential to improve mechanical harvest efficiency in tea plantations.

Genome-wide association study and network analysis of in vitro transformation in Populus trichocarpa support key roles of diverse phytohormone pathways and cross talk

Wide variation in amenability to transformation and regeneration (TR) among many plant species and genotypes presents a challenge to the use of genetic engineering in research and breeding. To help understand the causes of this variation, we performed association mapping and network analysis using a population of 1204 wild trees of Populus trichocarpa (black cottonwood). To enable precise and high-throughput phenotyping of callus and shoot TR, we developed a computer vision system that cross-referenced complementary red, green, and blue (RGB) and fluorescent-hyperspectral images. We performed association mapping using single-marker and combined variant methods, followed by statistical tests for epistasis and integration of published multi-omic datasets to identify likely regulatory hubs. We report 409 candidate genes implicated by associations within 5 kb of coding sequences, and epistasis tests implicated 81 of these candidate genes as regulators of one another. Gene ontology terms related to protein-protein interactions and transcriptional regulation are overrepresented, among others. In addition to auxin and cytokinin pathways long established as critical to TR, our results highlight the importance of stress and wounding pathways. Potential regulatory hubs of signaling within and across these pathways include GROWTH REGULATORY FACTOR 1 (GRF1), PHOSPHATIDYLINOSITOL 4-KINASE β1 (PI-4Kβ1), and OBF-BINDING PROTEIN 1 (OBP1).

GLK transcription factors accompany ELONGATED HYPOCOTYL5 to orchestrate light-induced seedling development in Arabidopsis

Light-induced de-etiolation is an important aspect of seedling photomorphogenesis. GOLDEN2 LIKE (GLK) transcriptional regulators are involved in chloroplast development, but to what extent they participate in photomorphogenesis is not clear. Here, we show that ELONGATED HYPOCOTYL5 (HY5) binds to GLK promoters to activate their expression, and also interacts with GLK proteins in Arabidopsis (Arabidopsis thaliana). The chlorophyll content in the de-etiolating Arabidopsis seedlings of the hy5 glk2 double mutants was lower than that in the hy5 single mutant. GLKs inhibited hypocotyl elongation, and the phenotype could superimpose on the hy5 phenotype. Correspondingly, GLK2 regulated the expression of photosynthesis and cell elongation genes partially independent of HY5. Before exposure to light, DE-ETIOLATED 1 (DET1) affected accumulation of GLK proteins. The enhanced etioplast development and photosystem gene expression observed in the det1 mutant were attenuated in the det1 glk2 double mutant. Our study reveals that GLKs act downstream of HY5, or additive to HY5, and are likely quantitatively adjusted by DET1, to orchestrate multiple developmental traits during the light-induced skotomorphogenesis-to-photomorphogenesis transition in Arabidopsis.

Ambient temperature regulates root circumnutation in rice through the ethylene pathway: transcriptome analysis reveals key genes involved

Plant roots are constantly prepared to adjust their growth trajectories to avoid unfavorable environments, and their ability to reorient is particularly crucial for survival. Under laboratory conditions, this continuous reorientation of the root tip is manifested as coiling or waving, which we refer to as root circumnutation. However, the effect of ambient temperature (AT) on root circumnutation remains unexplored. In this study, rice seedlings were employed to assess the impact of varying ATs on root circumnutation. The role of ethylene in mediating root circumnutation under elevated AT was examined using the ethylene biosynthesis inhibitor aminooxyacetic acid (AOA) and the ethylene perception antagonist silver thiosulfate (STS). Furthermore, transcriptome sequencing, weighted gene co-expression network analysis, and real-time quantitative PCR were utilized to analyze gene expressions in rice root tips under four distinct treatments: 25°C, 35°C, 35°C+STS, and 35°C+AOA. As a result, genes associated with ethylene synthesis and signaling (OsACOs and OsERFs), auxin synthesis and transport (OsYUCCA6, OsABCB15, and OsNPFs), cell elongation (OsEXPAs, OsXTHs, OsEGL1, and OsEXORDIUMs), as well as the inhibition of root curling (OsRMC) were identified. Notably, the expression levels of these genes increased with rising temperatures above 25°C. This study is the first to demonstrate that elevated AT can induce root circumnutation in rice via the ethylene pathway and proposes a potential molecular model through the identification of key genes. These findings offer valuable insights into the growth regulation mechanism of plant roots under elevated AT conditions.

Natural genetic variation in GLK1-mediated photosynthetic acclimation in response to light

BACKGROUND

GOLDEN-like (GLK) transcription factors are central regulators of chloroplast biogenesis in Arabidopsis and other species. Findings from Arabidopsis show that these factors also contribute to photosynthetic acclimation, e.g. to variation in light intensity, and are controlled by retrograde signals emanating from the chloroplast. However, the natural variation of GLK1-centered gene-regulatory networks in Arabidopsis is largely unexplored.

RESULTS

By evaluating the activities of GLK1 target genes and GLK1 itself in vegetative leaves of natural Arabidopsis accessions grown under standard conditions, we uncovered variation in the activity of GLK1 centered regulatory networks. This is linked with the ecogeographic origin of the accessions, and can be associated with a complex genetic variation across loci acting in different functional pathways, including photosynthesis, ROS and brassinosteroid pathways. Our results identify candidate upstream regulators that contribute to a basal level of GLK1 activity in rosette leaves, which can then impact the capacity to acclimate to different environmental conditions. Indeed, accessions with higher GLK1 activity, arising from habitats with a high monthly variation in solar radiation levels, may show lower levels of photoinhibition at higher light intensities.

CONCLUSIONS

Our results provide evidence for natural variation in GLK1 regulatory activities in vegetative leaves. This variation is associated with ecogeographic origin and can contribute to acclimation to high light conditions.

Flexure wood formation via growth reprogramming in hybrid aspen involves jasmonates and polyamines and transcriptional changes resembling tension wood development

Stem bending in trees induces flexure wood but its properties and development are poorly understood. Here, we investigated the effects of low-intensity multidirectional stem flexing on growth and wood properties of hybrid aspen, and on its transcriptomic and hormonal responses. Glasshouse-grown trees were either kept stationary or subjected to several daily shakes for 5 wk, after which the transcriptomes and hormones were analyzed in the cambial region and developing wood tissues, and the wood properties were analyzed by physical, chemical and microscopy techniques. Shaking increased primary and secondary growth and altered wood differentiation by stimulating gelatinous-fiber formation, reducing secondary wall thickness, changing matrix polysaccharides and increasing cellulose, G- and H-lignin contents, cell wall porosity and saccharification yields. Wood-forming tissues exhibited elevated jasmonate, polyamine, ethylene and brassinosteroids and reduced abscisic acid and gibberellin signaling. Transcriptional responses resembled those during tension wood formation but not opposite wood formation and revealed several thigmomorphogenesis-related genes as well as novel gene networks including FLA and XTH genes encoding plasma membrane-bound proteins. Low-intensity stem flexing stimulates growth and induces wood having improved biorefinery properties through molecular and hormonal pathways similar to thigmomorphogenesis in herbaceous plants and largely overlapping with the tension wood program of hardwoods.

Environmental gradients reveal stress hubs pre-dating plant terrestrialization

Plant terrestrialization brought forth the land plants (embryophytes). Embryophytes account for most of the biomass on land and evolved from streptophyte algae in a singular event. Recent advances have unravelled the first full genomes of the closest algal relatives of land plants; among the first such species was Mesotaenium endlicherianum. Here we used fine-combed RNA sequencing in tandem with a photophysiological assessment on Mesotaenium exposed to a continuous range of temperature and light cues. Our data establish a grid of 42 different conditions, resulting in 128 transcriptomes and ~1.5 Tbp (~9.9 billion reads) of data to study the combinatory effects of stress response using clustering along gradients. Mesotaenium shares with land plants major hubs in genetic networks underpinning stress response and acclimation. Our data suggest that lipid droplet formation and plastid and cell wall-derived signals have denominated molecular programmes since more than 600 million years of streptophyte evolution-before plants made their first steps on land.

Proteome-Level Investigation of Vitis amurensis Calli Transformed with a Constitutively Active, Ca2+-Independent Form of the Arabidopsis AtCPK1 Gene

Calcium-dependent protein kinases (CDPKs) are one of the main Ca2+ decoders in plants. Among them, Arabidopsis thaliana AtCPK1 is one of the most studied CDPK genes as a positive regulator of plant responses to biotic and abiotic stress. The mutated form of AtCPK1, in which the autoinhibitory domain is inactivated (AtCPK1-Ca), provides constitutive kinase activity by mimicking a stress-induced increase in the Ca2+ flux. In the present study, we performed a proteomic analysis of Vitis amurensis calli overexpressing the AtCPK1-Ca form using untransformed calli as a control. In our previous studies, we have shown that the overexpression of this mutant form leads to the activation of secondary metabolism in plant cell cultures, including an increase in resveratrol biosynthesis in V. amurensis cell cultures. We analyzed upregulated and downregulated proteins in control and transgenic callus cultures using two-dimensional gel electrophoresis, and Matrix-assisted laser desorption and ionization time-of-flight mass spectrometry (MALDI-TOF). In calli transformed with AtCPK1-Ca, an increased amounts of pathogenesis-related proteins were found. A quantitative real-time PCR analysis confirmed this result.

Integration of GWAS and RNA-Seq Analysis to Identify SNPs and Candidate Genes Associated with Alkali Stress Tolerance at the Germination Stage in Mung Bean

Soil salt-alkalization seriously impacts crop growth and productivity worldwide. Breeding and applying tolerant varieties is the most economical and effective way to address soil alkalization. However, genetic resources for breeders to improve alkali tolerance are limited in mung bean. Here, a genome-wide association study (GWAS) was performed to detect alkali-tolerant genetic loci and candidate genes in 277 mung bean accessions during germination. Using the relative values of two germination traits, 19 QTLs containing 32 SNPs significantly associated with alkali tolerance on nine chromosomes were identified, and they explained 3.6 to 14.6% of the phenotypic variance. Moreover, 691 candidate genes were mined within the LD intervals containing significant trait-associated SNPs. Transcriptome sequencing of alkali-tolerant accession 132-346 under alkali and control conditions after 24 h of treatment was conducted, and 2565 DEGs were identified. An integrated analysis of the GWAS and DEGs revealed six hub genes involved in alkali tolerance responses. Moreover, the expression of hub genes was further validated by qRT-PCR. These findings improve our understanding of the molecular mechanism of alkali stress tolerance and provide potential resources (SNPs and genes) for the genetic improvement of alkali tolerance in mung bean.

Transcriptomic analysis of Chinese yam (Dioscorea polystachya Turcz.) variants indicates brassinosteroid involvement in tuber development

Dioscorea is an important but underutilized genus of flowering plants that grows predominantly in tropical and subtropical regions. Several species, known as yam, develop large underground tubers and aerial bulbils that are used as food. The Chinese yam (D. polystachya Turcz.) is one of the few Dioscorea species that grows well in temperate regions and has been proposed as a climate-resilient crop to enhance food security in Europe. However, the fragile, club-like tubers are unsuitable for mechanical harvesting, which is facilitated by shorter and thicker storage organs. Brassinosteroids (BRs) play a key role in plant cell division, cell elongation and proliferation, as well as in the gravitropic response. We collected RNA-Seq data from the head, middle and tip of two tuber shape variants: F60 (long, thin) and F2000 (short, thick). Comparative transcriptome analysis of F60 vs. F2000 revealed 30,229 differentially expressed genes (DEGs), 1,393 of which were differentially expressed in the growing tip. Several DEGs are involved in steroid/BR biosynthesis or signaling, or may be regulated by BRs. The quantification of endogenous BRs revealed higher levels of castasterone (CS), 28-norCS, 28-homoCS and brassinolide in F2000 compared to F60 tubers. The highest BR levels were detected in the growing tip, and CS was the most abundant (439.6 ± 196.41 pmol/g in F2000 and 365.6 ± 112.78 pmol/g in F60). Exogenous 24-epi-brassinolide (epi-BL) treatment (20 nM) in an aeroponic system significantly increased the width-to-length ratio (0.045 ± 0.002) compared to the mock-treated plants (0.03 ± 0.002) after 7 weeks, indicating that exogenous epi-BL produces shorter and thicker tubers. In this study we demonstrate the role of BRs in D. polystachya tuber shape, providing insight into the role of plant hormones in yam storage organ development. We found that BRs can influence tuber shape in Chinese yam by regulating the expression of genes involved cell expansion. Our data can help to improve the efficiency of Chinese yam cultivation, which could provide an alternative food source and thus contribute to future food security in Europe.

A high-confidence Physcomitrium patens plasmodesmata proteome by iterative scoring and validation reveals diversification of cell wall proteins during evolution

Plasmodesmata (PD) facilitate movement of molecules between plant cells. Regulation of this movement is still not understood. Plasmodesmata are hard to study, being deeply embedded within cell walls and incorporating several membrane types. Thus, structure and protein composition of PD remain enigmatic. Previous studies of PD protein composition identified protein lists with few validations, making functional conclusions difficult. We developed a PD scoring approach in iteration with large-scale systematic localization, defining a high-confidence PD proteome of Physcomitrium patens (HC300). HC300, together with bona fide PD proteins from literature, were placed in Pddb. About 65% of proteins in HC300 were not previously PD-localized. Callose-degrading glycolyl hydrolase family 17 (GHL17) is an abundant protein family with representatives across evolutionary scale. Among GHL17s, we exclusively found members of one phylogenetic clade with PD localization and orthologs occur only in species with developed PD. Phylogenetic comparison was expanded to xyloglucan endotransglucosylases/hydrolases and Exordium-like proteins, which also diversified into PD-localized and non-PD-localized members on distinct phylogenetic clades. Our high-confidence PD proteome HC300 provides insights into diversification of large protein families. Iterative and systematic large-scale localization across plant species strengthens the reliability of HC300 as basis for exploring structure, function, and evolution of this important organelle.

Genetic Adaptation of Siberian Larch (Larix sibirica Ledeb.) to High Altitudes

Forest trees growing in high altitude conditions offer a convenient model for studying adaptation processes. They are subject to a whole range of adverse factors that are likely to cause local adaptation and related genetic changes. Siberian larch (Larix sibirica Ledeb.), whose distribution covers different altitudes, makes it possible to directly compare lowland with highland populations. This paper presents for the first time the results of studying the genetic differentiation of Siberian larch populations, presumably associated with adaptation to the altitudinal gradient of climatic conditions, based on a joint analysis of altitude and six other bioclimatic variables, together with a large number of genetic markers, single nucleotide polymorphisms (SNPs), obtained from double digest restriction-site-associated DNA sequencing (ddRADseq). In total, 25,143 SNPs were genotyped in 231 trees. In addition, a dataset of 761 supposedly selectively neutral SNPs was assembled by selecting SNPs located outside coding regions in the Siberian larch genome and mapped to different contigs. The analysis using four different methods (PCAdapt, LFMM, BayeScEnv and RDA) revealed 550 outlier SNPs, including 207 SNPs whose variation was significantly correlated with the variation of some of environmental factors and presumably associated with local adaptation, including 67 SNPs that correlated with altitude based on either LFMM or BayeScEnv and 23 SNPs based on both of them. Twenty SNPs were found in the coding regions of genes, and 16 of them represented non-synonymous nucleotide substitutions. They are located in genes involved in the processes of macromolecular cell metabolism and organic biosynthesis associated with reproduction and development, as well as organismal response to stress. Among these 20 SNPs, nine were possibly associated with altitude, but only one of them was identified as associated with altitude by all four methods used in the study, a nonsynonymous SNP in scaffold_31130 in position 28092, a gene encoding a cell membrane protein with uncertain function. Among the studied populations, at least two main groups (clusters), the Altai populations and all others, were significantly genetically different according to the admixture analysis based on any of the three SNP datasets as follows: 761 supposedly selectively neutral SNPs, all 25,143 SNPs and 550 adaptive SNPs. In general, according to the AMOVA results, genetic differentiation between transects or regions or between population samples was relatively low, although statistically significant, based on 761 neutral SNPs (F_ST = 0.036) and all 25,143 SNPs (F_ST = 0.017). Meanwhile, the differentiation based on 550 adaptive SNPs was much higher (F_ST = 0.218). The data showed a relatively weak but highly significant linear correlation between genetic and geographic distances (r = 0.206, p = 0.001).

Nectary development in Cleome violacea

Nectaries are a promising frontier for plant evo-devo research, and are particularly fascinating given their diversity in form, position, and secretion methods across angiosperms. Emerging model systems permit investigations of the molecular basis for nectary development and nectar secretion across a range of taxa, which addresses fundamental questions about underlying parallelisms and convergence. Herein, we explore nectary development and nectar secretion in the emerging model taxa, Cleome violacea (Cleomaceae), which exhibits a prominent adaxial nectary. First, we characterized nectary anatomy and quantified nectar secretion to establish a foundation for quantitative and functional gene experiments. Next, we leveraged RNA-seq to establish gene expression profiles of nectaries across three key stages of development: pre-anthesis, anthesis, and post-fertilization. We then performed functional studies on five genes that were putatively involved in nectary and nectar formation: CvCRABSCLAW (CvCRC), CvAGAMOUS (CvAG), CvSHATTERPROOF (CvSHP), CvSWEET9, and a highly expressed but uncharacterized transcript. These experiments revealed a high degree of functional convergence to homologues from other core Eudicots, especially Arabidopsis. CvCRC, redundantly with CvAG and CvSHP, are required for nectary initiation. Concordantly, CvSWEET9 is essential for nectar formation and secretion, which indicates that the process is eccrine based in C. violacea. While demonstration of conservation is informative to our understanding of nectary evolution, questions remain. For example, it is unknown which genes are downstream of the developmental initiators CvCRC, CvAG, and CvSHP, or what role the TCP gene family plays in nectary initiation in this family. Further to this, we have initiated a characterization of associations between nectaries, yeast, and bacteria, but more research is required beyond establishing their presence. Cleome violacea is an excellent model for continued research into nectary development because of its conspicuous nectaries, short generation time, and close taxonomic distance to Arabidopsis.

GhBES1 mediates brassinosteroid regulation of leaf size by activating expression of GhEXO2 in cotton (Gossypium hirsutum)

We proposed a working model of BR to promote leaf size through cell expansion. In the BR signaling pathway, GhBES1 affects cotton leaf size by binding to and activating the expression of the E-box element in the GhEXO2 promoter region. Brassinosteroid (BR) is an essential phytohormone that controls plant growth. However, the mechanisms of BR regulation of leaf size remain to be determined. Here, we found that the BR deficient cotton mutant pagoda1 (pag1) had a smaller leaf size than wild-type CRI24. The expression of EXORDIUM (GhEXO2) gene, was significantly downregulated in pag1. Silencing of BRI1-EMS-SUPPRESSOR 1 (GhBES1), inhibited leaf cell expansion and reduced leaf size. Overexpression of GhBES1.4 promoted leaf cell expansion and enlarged leaf size. Expression analysis showed GhEXO2 expression positively correlated with GhBES1 expression. In plants, altered expression of GhEXO2 promoted leaf cell expansion affecting leaf size. Furthermore, GhBES1.4 specifically binds to the E-box elements in the GhEXO2 promoter, inducing its expression. RNA-seq data revealed many down-regulated genes related to cell expansion in GhEXO2 silenced plants. In summary, we discovered a novel mechanism of BR regulation of leaf size through GhBES1 directly activating the expression of GhEXO2.

Gene expression and genetic divergence in oak species highlight adaptive genes to soil water constraints

Drought and waterlogging impede tree growth and may even lead to tree death. Oaks, an emblematic group of tree species, have evolved a range of adaptations to cope with these constraints. The two most widely distributed European species, pedunculate (PO; Quercus robur L.) and sessile oak (SO; Quercus petraea Matt. Lieb), have overlapping ranges, but their respective distribution are highly constrained by local soil conditions. These contrasting ecological preferences between two closely related and frequently hybridizing species constitute a powerful model to explore the functional bases of the adaptive responses in oak. We exposed oak seedlings to waterlogging and drought, conditions typically encountered by the two species in their respective habitats, and studied changes in gene expression in roots using RNA-seq. We identified genes that change in expression between treatments differentially depending on species. These "species × environment"-responsive genes revealed adaptive molecular strategies involving adventitious and lateral root formation, aerenchyma formation in PO, and osmoregulation and ABA regulation in SO. With this experimental design, we also identified genes with different expression between species independently of water conditions imposed. Surprisingly, this category included genes with functions consistent with a role in intrinsic reproductive barriers. Finally, we compared our findings with those for a genome scan of species divergence and found that the expressional candidate genes included numerous highly differentiated genetic markers between the two species. By combining transcriptomic analysis, gene annotation, pathway analyses, as well as genome scan for genetic differentiation among species, we were able to highlight loci likely involved in adaptation of the two species to their respective ecological niches.

GhTCE1-GhTCEE1 dimers regulate transcriptional reprogramming during wound-induced callus formation in cotton

Wounded plant cells can form callus to seal the wound site. Alternatively, wounding can cause adventitious organogenesis or somatic embryogenesis. These distinct developmental pathways require specific cell fate decisions. Here, we identify GhTCE1, a basic helix-loop-helix family transcription factor, and its interacting partners as a central regulatory module of early cell fate transition during in vitro dedifferentiation of cotton (Gossypium hirsutum). RNAi- or CRISPR/Cas9-mediated loss of GhTCE1 function resulted in excessive accumulation of reactive oxygen species (ROS), arrested callus cell elongation, and increased adventitious organogenesis. In contrast, GhTCE1-overexpressing tissues underwent callus cell growth, but organogenesis was repressed. Transcriptome analysis revealed that several pathways depend on proper regulation of GhTCE1 expression, including lipid transfer pathway components, ROS homeostasis, and cell expansion. GhTCE1 bound to the promoters of the target genes GhLTP2 and GhLTP3, activating their expression synergistically, and the heterodimer TCE1-TCEE1 enhances this activity. GhLTP2- and GhLTP3-deficient tissues accumulated ROS and had arrested callus cell elongation, which was restored by ROS scavengers. These results reveal a unique regulatory network involving ROS and lipid transfer proteins, which act as potential ROS scavengers. This network acts as a switch between unorganized callus growth and organized development during in vitro dedifferentiation of cotton cells.

Insight into the mechanism of indium toxicity in rice

Indium is widely used in the technology industry and is an emerging form of environmental pollution. The presence of indium in soil and groundwater inhibits shoot and root growth in crops, thus reducing yields. However, the underlying mechanisms are unknown, making it difficult to design effective countermeasures. We explored the spatiotemporal effects of excess indium on the morphological, physiological and biochemical properties of rice (Oryza sativa L.). Indium accumulated mainly in the roots, severely restricting their growth and causing the acute perturbation of phosphorus, magnesium and iron homeostasis. Other effects included leaf necrosis and anatomical changes in the roots (thinned sclerenchyma and enlarged epidermal and exodermal layers). Whole-transcriptome sequencing revealed that rice immediately responded to indium stress by activating genes involved in heavy metal tolerance and phosphate starvation responses, including the expression of genes encoding phosphate-regulated transcription factors and transporters in the roots. Direct indium toxicity rather than phosphate deficiency was identified as the major factor affecting the growth of rice plants, resulting in the profound phenotypic changes we observed. The application of exogenous phosphate alleviated indium toxicity by reducing indium uptake. Our results suggest that indium immobilization could be used to prevent indium toxicity in the field.

WHIRLY1 functions in the nucleus to regulate barley leaf development and associated metabolite profiles

The WHIRLY (WHY) DNA/RNA binding proteins fulfil multiple but poorly characterised functions in leaf development. Here, we show that WHY1 transcript levels were highest in the bases of 7-day old barley leaves. Immunogold labelling revealed that the WHY1 protein was more abundant in the nuclei than the proplastids of the leaf bases. To identify transcripts associated with leaf development we conducted hierarchical clustering of differentially abundant transcripts along the developmental gradient of wild-type leaves. Similarly, metabolite profiling was employed to identify metabolites exhibiting a developmental gradient. A comparative analysis of transcripts and metabolites in barley lines (W1–1 and W1–7) lacking WHY1, which show delayed greening compared with the wild type revealed that the transcript profile of leaf development was largely unchanged in W1–1 and W1–7 leaves. However, there were differences in levels of several transcripts encoding transcription factors associated with chloroplast development. These include a barley homologue of the Arabidopsis GATA transcription factor that regulates stomatal development, greening and chloroplast development, NAC1; two transcripts with similarity to Arabidopsis GLK1 and two transcripts encoding ARF transcriptions factors with functions in leaf morphogenesis and development. Chloroplast proteins were less abundant in the W1–1 and W1–7 leaves than the wild type. The levels of tricarboxylic acid cycle metabolites and GABA were significantly lower in WHY1 knockdown leaves than the wild type. This study provides evidence that WHY1 is localised in the nuclei of leaf bases, contributing the regulation of nuclear-encoded transcripts that regulate chloroplast development.

Genome Wide Association Study Identifies Candidate Genes Related to the Earlywood Tracheid Properties in Picea crassifolia Kom

Picea crassifolia Kom. is one of the timber and ecological conifers in China and its wood tracheid traits directly affect wood formation and adaptability under harsh environment. Molecular studies on P. crassifolia remain inadequate because relatively few genes have been associated with these traits. To identify markers and candidate genes that can potentially be used for genetic improvement of wood tracheid traits, we examined 106 clones of P. crassifolia, and investigated phenotypic data for 14 wood tracheid traits before specific-locus amplified fragment sequencing (SLAF-seq) was employed to perform a genome wide association study (GWAS). Subsequently, the results were used to screen single nucleotide polymorphism (SNP) loci and candidate genes that exhibited a significant correlation with the studied traits. We developed 4,058,883 SLAF-tags and 12,275,765 SNP loci, and our analyses identified a total of 96 SNP loci that showed significant correlations with three earlywood tracheid traits using a mixed linear model (MLM). Next, candidate genes were screened in the 100 kb zone (50 kb upstream, 50 kb downstream) of each of the SNP loci, whereby 67 candidate genes were obtained in earlywood tracheid traits, including 34 genes of known function and 33 genes of unknown function. We provide the most significant SNP for each trait-locus combination and candidate genes occurring within the GWAS hits. These resources provide a foundation for the development of markers that could be used in wood traits improvement and candidate genes for the development of earlywood tracheid in P. crassifolia.

The genome and preliminary single-nuclei transcriptome of Lemna minuta reveals mechanisms of invasiveness

The ability to trace every cell in some model organisms has led to the fundamental understanding of development and cellular function. However, in plants the complexity of cell number, organ size, and developmental time makes this a challenge even in the diminutive model plant Arabidopsis (Arabidopsis thaliana). Duckweed, basal nongrass aquatic monocots, provide an opportunity to follow every cell of an entire plant due to their small size, reduced body plan, and fast clonal growth habit. Here we present a chromosome-resolved genome for the highly invasive Lesser Duckweed (Lemna minuta) and generate a preliminary cell atlas leveraging low cell coverage single-nuclei sequencing. We resolved the 360 megabase genome into 21 chromosomes, revealing a core nonredundant gene set with only the ancient tau whole-genome duplication shared with all monocots, and paralog expansion as a result of tandem duplications related to phytoremediation. Leveraging SMARTseq2 single-nuclei sequencing, which provided higher gene coverage yet lower cell count, we profiled 269 nuclei covering 36.9% (8,457) of the L. minuta transcriptome. Since molecular validation was not possible in this nonmodel plant, we leveraged gene orthology with model organism single-cell expression datasets, gene ontology, and cell trajectory analysis to define putative cell types. We found that the tissue that we computationally defined as mesophyll expressed high levels of elemental transport genes consistent with this tissue playing a role in L. minuta wastewater detoxification. The L. minuta genome and preliminary cell map provide a paradigm to decipher developmental genes and pathways for an entire plant.

Comparative analysis of differential gene expression indicates divergence in ontogenetic strategies of leaves in two conifer genera

In land plants, heteroblasty broadly refers to a drastic change in morphology during growth through ontogeny. Juniperus flaccida and Pinus cembroides are conifers of independent lineages known to exhibit leaf heteroblasty between the juvenile and adult life stage of development. Juvenile leaves of P. cembroides develop spirally on the main stem and appear decurrent, flattened, and needle-like; whereas adult photosynthetic leaves are triangular or semi-circular needle-like, and grow in whorls on secondary or tertiary compact dwarf shoots. By comparison, J. flaccida juvenile leaves are decurrent and needle-like, and adult leaves are compact, short, and scale-like. Comparative analyses were performed to evaluate differences in anatomy and gene expression patterns between developmental phases in both species. RNA from 12 samples was sequenced and analyzed with available software. They were assembled de novo from the RNA-Seq reads. Following assembly, 63,741 high-quality transcripts were functionally annotated in P. cembroides and 69,448 in J. flaccida. Evaluation of the orthologous groups yielded 4140 shared gene families among the four references (adult and juvenile from each species). Activities related to cell division and development were more abundant in juveniles than adults in P. cembroides, and more abundant in adults than juveniles in J. flaccida. Overall, there were 509 up-regulated and 81 down-regulated genes in the juvenile condition of P. cembroides and 14 up-regulated and 22 down-regulated genes in J. flaccida. Gene interaction network analysis showed evidence of co-expression and co-localization of up-regulated genes involved in cell wall and cuticle formation, development, and phenylpropanoid pathway, in juvenile P. cembroides leaves. Whereas in J. flaccida, differential expression and gene interaction patterns were detected in genes involved in photosynthesis and chloroplast biogenesis. Although J. flaccida and P. cembroides both exhibit leaf heteroblastic development, little overlap was detected, and unique genes and pathways were highlighted in this study.

Transcriptome analysis at mid-stage seed development in litchi with contrasting seed size

Litchi is a sub-tropical fruit crop with genotypes that bear fruits with variable seed size. Small seed size is a desirable trait in litchi, as it improves consumers' preference and facilitates fruit processing. Seed specific transcriptome analysis was performed in two litchi genotypes with contrasting seed size to identify the genes associated with seed development. The transcriptomic sequence data from seeds at mid-development stages (16-28 days after anthesis) were de-novo assembled into 1,39,608 Trinity transcripts. Out of these, 6325 transcripts expressed differentially between the two contrasting genotypes. Several putative genes for salicylic acid, jasmonic acid and brassinosteriod pathways were down-regulated in seeds of the small-seeded litchi. The putative regulators of seed maturation and seed storage were down-regulated in the small-seeded genotype. Embryogenesis, cell expansion, seed size and stress related Trinity transcripts exhibited differential expression. Further studies on gene characterization will reveal the early regulators of seed size in litchi.

SUPPLEMENTARY INFORMATION

The online version contains supplementary material available at 10.1007/s13205-021-03098-8.

Integrative analysis of the metabolome and transcriptome reveal the phosphate deficiency response pathways of alfalfa

Understanding the mechanisms underlying the responses to inorganic phosphate (Pi) deficiency in alfalfa will help enhance Pi acquisition efficiency and the sustainable use of phosphorous resources. Integrated global metabolomic and transcriptomic analyses of mid-vegetative alfalfa seedlings under 12-day Pi deficiency were conducted. Limited seedling growth were found, including 13.24%, 16.85% and 33.36% decreases in height, root length and photosynthesis, and a 24.10% increase in root-to-shoot ratio on day 12. A total of 322 and 448 differentially abundant metabolites and 1199 and 1061 differentially expressed genes were identified in roots and shoots. Increased (>3.68-fold) inorganic phosphate transporter 1;4 and SPX proteins levels in the roots (>2.15-fold) and shoots (>2.50-fold) were related to Pi absorption and translocation. The levels of phospholipids and Pi-binding carbohydrates and nucleosides were decreased, while those of phosphatases and pyrophosphatases in whole seedlings were induced under reduced Pi. In addition, nitrogen assimilation was affected by inhibiting high-affinity nitrate transporters (NRT2.1 and NRT3.1), and nitrate reductase. Increased delphinidin-3-glucoside might contribute to the gray-green leaves induced by Pi limitation. Stress-induced MYB, WRKY and ERF transcription factors were identified. The responses of alfalfa to Pi deficiency were summarized as local systemic signaling pathways, including root growth, stress-related responses consisting of enzymatic and nonenzymatic systems, and hormone signaling and systemic signaling pathways including Pi recycling and Pi sensing in the whole plant, as well as Pi recovery, and nitrate and metal absorption in the roots. This study provides important information on the molecular mechanism of the response to Pi deficiency in alfalfa.

The XTH Gene Family in Schima superba: Genome-Wide Identification, Expression Profiles, and Functional Interaction Network Analysis

Xyloglucan endotransglucosylase/hydrolase (XTH), belonging to glycoside hydrolase family 16, is one of the key enzymes in plant cell wall remodeling. Schima superba is an important timber and fireproof tree species in southern China. However, little is known about XTHs in S. superba. In the present study, a total of 34 SsuXTHs were obtained, which were classified into three subfamilies based on the phylogenetic relationship and unevenly distributed on 18 chromosomes. Furthermore, the intron-exon structure and conserved motif composition of them supported the classification and the members belonging to the same subfamily shared similar gene structures. Segmental and tandem duplication events did not lead to SsuXTH gene family expansion, and strong purifying selection pressures during evolution led to similar structure and function of SsuXTH gene family. The interaction network and cis-acting regulatory elements analysis revealed the SsuXTH expression might be regulated by multiple hormones, abiotic stresses and transcription factors. Finally, expression profiles and GO enrichment analysis showed most of the tandem repeat genes were mainly expressed in the phloem and xylem and they mainly participated in glycoside metabolic processes through the transfer and hydrolysis of xyloglucan in the cell wall and then regulated fiber elongation.

Genetic Dissection of Mature Root Characteristics by Genome-Wide Association Studies in Rapeseed (Brassica napus L.)

Roots are complicated quantitative characteristics that play an essential role in absorbing water and nutrients. To uncover the genetic variations for root-related traits in rapeseed, twelve mature root traits of a Brassica napus association panel were investigated in the field within three environments. All traits showed significant phenotypic variation among genotypes, with heritabilities ranging from 55.18% to 79.68%. Genome-wide association studies (GWAS) using 20,131 SNPs discovered 172 marker-trait associations, including 103 significant SNPs (-log10 (p) > 4.30) that explained 5.24-20.31% of the phenotypic variance. With the linkage disequilibrium r² > 0.2, these significant associations were binned into 40 quantitative trait loci (QTL) clusters. Among them, 14 important QTL clusters were discovered in two environments and/or with phenotypic contributions greater than 10%. By analyzing the genomic regions within 100 kb upstream and downstream of the peak SNPs within the 14 loci, 334 annotated genes were found. Among these, 32 genes were potentially associated with root development according to their expression analysis. Furthermore, the protein interaction network using the 334 annotated genes gave nine genes involved in a substantial number of interactions, including a key gene associated with root development, BnaC09g36350D. This research provides the groundwork for deciphering B. napus' genetic variations and improving its root system architecture.

Comparing Early Transcriptomic Responses of 18 Soybean (Glycine max) Genotypes to Iron Stress

Iron deficiency chlorosis (IDC) is an abiotic stress that negatively affects soybean (Glycine max [L.] Merr.) production. Much of our knowledge of IDC stress responses is derived from model plant species. Gene expression, quantitative trait loci (QTL) mapping, and genome-wide association studies (GWAS) performed in soybean suggest that stress response differences exist between model and crop species. Our current understanding of the molecular response to IDC in soybeans is largely derived from gene expression studies using near-isogenic lines differing in iron efficiency. To improve iron efficiency in soybeans and other crops, we need to expand gene expression studies to include the diversity present in germplasm collections. Therefore, we collected 216 purified RNA samples (18 genotypes, two tissue types [leaves and roots], two iron treatments [sufficient and deficient], three replicates) and used RNA sequencing to examine the expression differences of 18 diverse soybean genotypes in response to iron deficiency. We found a rapid response to iron deficiency across genotypes, most responding within 60 min of stress. There was little evidence of an overlap of specific differentially expressed genes, and comparisons of gene ontology terms and transcription factor families suggest the utilization of different pathways in the stress response. These initial findings suggest an untapped genetic potential within the soybean germplasm collection that could be used for the continued improvement of iron efficiency in soybean.

A Non-Shedding Fruit Elaeis oleifera Palm Reveals Perturbations to Hormone Signaling, ROS Homeostasis, and Hemicellulose Metabolism

The developmentally programmed loss of a plant organ is called abscission. This process is characterized by the ultimate separation of adjacent cells in the abscission zone (AZ). The discovery of an American oil palm (Elaeis oleifera) variant that does not shed its has allowed for the study of the mechanisms of ripe fruit abscission in this species. A comparative transcriptome analysis was performed to compare the fruit AZs of the non-shedding E. oleifera variant to an individual of the same progeny that sheds its ripe fruit normally. The study provides evidence for widespread perturbation to gene expression in the AZ of the non-shedding variant, compared to the normal fruit-shedding control, and offers insight into abscission-related functions. Beyond the genes with known or suspected roles during organ abscission or indehiscence that were identified, a list of genes with hormone-related functions, including ethylene, jasmonic acid, abscisic acid, cytokinin and salicylic acid, in addition to reactive oxygen species (ROS) metabolism, transcriptional responses and signaling pathways, was compiled. The results also allowed a comparison between the ripe fruit abscission processes of the African and American oil palm species at the molecular level and revealed commonalities with environmental stress pathways.

Transcriptomic study of the night break in Chenopodium rubrum reveals possible upstream regulators of the floral activator CrFTL1

The transition from vegetative to reproductive phases is the most fundamental and tightly controlled switch in the life of flowering plants. The short-day plant Chenopodium rubrum is a fast cycling annual plant lacking a juvenile phase. It can be induced to flowering at the seedling stage by exposure to a single period of darkness. This floral induction may then be cancelled by a short pulse of red light at midnight called night break (NB), which also inhibits the floral activator FLOWERING LOCUS T LIKE 1 (CrFTL1). We performed a comparative transcriptomic study between C. rubrum seedlings treated by NB and ones growing through uninterrupted night, and found about six hundred differentially expressed genes, including the B-BOX DOMAIN (BBX) genes. We focused on the CrBBX19 and BOLTING TIME CONTROL 1 (BTC1) genes, homologous to the upstream regulators of the BvFT2, a floral inducer in sugar beet. The transcription patterns of the two genes were compatible with their putative role as a sensor of the dark period length optimal for flowering (CrBBX19), and a signal of lights-on (CrBTC1), but the participation of other genes cannot be excluded. The expression profiles of CrBBX19 and the homolog of the core endogenous clock gene LATE ELONGATED HYPOCOTYL (LHY) were highly similar, which suggested their co-regulation.

Identification and Analysis of GhEXO Gene Family Indicated That GhEXO7_At Promotes Plant Growth and Development Through Brassinosteroid Signaling in Cotton (Gossypium hirsutum L.)

Brassinosteroids (BRs), an efficient plant endogenous hormone, significantly promotes plant nutrient growth adapting to biological and abiotic adversities. BRs mainly promote plant cell elongation by regulating gene expression patterns. EXORDIUM (EXO) genes have been characterized as the indicators of BR response genes. Cotton, an ancient crop, is of great economic value and its fibers can be made into all kinds of fabrics. However, EXO gene family genes have not been full identified in cotton. 175 EXO genes were identified in nine plant species, of which 39 GhEXO genes in Gossypium hirsutum in our study. A phylogenetic analysis grouped all of the proteins encoded by the EXO genes into five major clades. Sequence identification of conserved amino acid residues among monocotyledonous and dicotyledonous species showed a high level of conservation across the N and C terminal regions. Only 25% the GhEXO genes contain introns besides conserved gene structure and protein motifs distribution. The 39 GhEXO genes were unevenly distributed on the 18 At and Dt sub-genome chromosomes. Most of the GhEXO genes were derived from gene duplication events, while only three genes showed evidence of tandem duplication. Homologous locus relationships showed that 15 GhEXO genes are located on collinear blocks and that all orthologous/paralogous gene pairs had Ka > Ks values, indicating purifying selection pressure. The GhEXO genes showed ubiquitous expression in all eight tested cotton tissues and following exposure to three phytohormones, IAA, GA, and BL. Furthermore, GhEXO7_At was mainly expressed in response to BL treatment, and was predominantly expressed in the fibers. GhEXO7_At was found to be a plasma membrane protein, and its ectopic expression in Arabidopsis mediated BR-regulated plant growth and development with altered expression of DWF4, CPD, KCS1, and EXP5. Additionally, the functions of GhEXO7_At were confirmed by virus-induced gene silencing (VIGS) in cotton. This study will provide important genetic resources for future cotton breeding programs.

CsLOB1 regulates susceptibility to citrus canker through promoting cell proliferation in citrus

Citrus sinensis lateral organ boundary 1 (CsLOB1) was previously identified as a critical disease susceptibility gene for citrus bacterial canker, which is caused by Xanthomonas citri subsp. citri (Xcc). However, the molecular mechanisms of CsLOB1 in citrus response to Xcc are still elusive. Here, we constructed transgenic plants overexpressing and RNAi-silencing of CsLOB1 using the canker-disease susceptible 'wanjincheng' orange (C. sinensis Osbeck) as explants. CsLOB1-overexpressing plants exhibited dwarf phenotypes with smaller and thicker leaf, increased branches and adventitious buds clustered on stems. These phenotypes were followed by a process of pustule- and canker-like development that exhibited enhanced cell proliferation. Pectin depolymerization and expansin accumulation were enhanced by CsLOB1 overexpression, while cellulose and hemicellulose synthesis were increased by CsLOB1 silence. Whilst overexpression of CsLOB1 increased susceptibility, RNAi-silencing of CsLOB1 enhanced resistance to canker disease without impairing pathogen entry. Transcriptome analysis revealed that CsLOB1 positively regulated cell wall degradation and modification processes, cytokinin metabolism, and cell division. Additionally, 565 CsLOB1-targeted genes were identified in chromatin immunoprecipitation-sequencing (ChIP-seq) experiments. Motif discovery analysis revealed that the most highly overrepresented binding sites had a conserved 6-bp 'GCGGCG' consensus DNA motif. RNA-seq and ChIP-seq data suggested that CsLOB1 directly activates the expression of four genes involved in cell wall remodeling, and three genes that participate in cytokinin and brassinosteroid hormone pathways. Our findings indicate that CsLOB1 promotes cell proliferation by mechanisms depending on cell wall remodeling and phytohormone signaling, which may be critical to citrus canker development and bacterial growth in citrus.

Analysis of the chromosomal clustering of Fusarium-responsive wheat genes uncovers new players in the defence against head blight disease

There is increasing evidence that some functionally related, co-expressed genes cluster within eukaryotic genomes. We present a novel pipeline that delineates such eukaryotic gene clusters. Using this tool for bread wheat, we uncovered 44 clusters of genes that are responsive to the fungal pathogen Fusarium graminearum. As expected, these Fusarium-responsive gene clusters (FRGCs) included metabolic gene clusters, many of which are associated with disease resistance, but hitherto not described for wheat. However, the majority of the FRGCs are non-metabolic, many of which contain clusters of paralogues, including those implicated in plant disease responses, such as glutathione transferases, MAP kinases, and germin-like proteins. 20 of the FRGCs encode nonhomologous, non-metabolic genes (including defence-related genes). One of these clusters includes the characterised Fusarium resistance orphan gene, TaFROG. Eight of the FRGCs map within 6 FHB resistance loci. One small QTL on chromosome 7D (4.7 Mb) encodes eight Fusarium-responsive genes, five of which are within a FRGC. This study provides a new tool to identify genomic regions enriched in genes responsive to specific traits of interest and applied herein it highlighted gene families, genetic loci and biological pathways of importance in the response of wheat to disease.

Genetic mapping and transcriptomic characterization of a new fuzzless-tufted cottonseed mutant

Fiber mutants are unique and valuable resources for understanding the genetic and molecular mechanisms controlling initiation and development of cotton fibers that are extremely elongated single epidermal cells protruding from the seed coat of cottonseeds. In this study, we reported a new fuzzless-tufted cotton mutant (Gossypium hirsutum) and showed that fuzzless-tufted near-isogenic lines (NILs) had similar agronomic traits and a higher ginning efficiency compared to their recurrent parents with normal fuzzy seeds. Genetic analysis revealed that the mutant phenotype is determined by a single incomplete dominant locus, designated N5. The mutation was fine mapped to an approximately 250-kb interval containing 33 annotated genes using a combination of bulked segregant sequencing, SNP chip genotyping, and fine mapping. Comparative transcriptomic analysis using 0-6 days post-anthesis (dpa) ovules from NILs segregating for the phenotypes of fuzzless-tufted (mutant) and normal fuzzy cottonseeds (wild-type) uncovered candidate genes responsible for the mutant phenotype. It also revealed that the flanking region of the N5 locus is enriched with differentially expressed genes (DEGs) between the mutant and wild-type. Several of those DEGs are members of the gene families with demonstrated roles in cell initiation and elongation, such as calcium-dependent protein kinase and expansin. The transcriptome landscape of the mutant was significantly reprogrammed in the 6 dpa ovules and, to a less extent, in the 0 dpa ovules, but not in the 2 and 4 dpa ovules. At both 0 and 6 dpa, the reprogrammed mutant transcriptome was mainly associated with cell wall modifications and transmembrane transportation, while transcription factor activity was significantly altered in the 6 dpa mutant ovules. These results imply a similar molecular basis for initiation of lint and fuzz fibers despite certain differences.

Integrative transcriptomics reveals genotypic impact on sugar beet storability

An integrative comparative transcriptomic approach on six sugar beet varieties showing different amount of sucrose loss during storage revealed genotype-specific main driver genes and pathways characterizing storability. Sugar beet is next to sugar cane one of the most important sugar crops accounting for about 15% of the sucrose produced worldwide. Since its processing is increasingly centralized, storage of beet roots over an extended time has become necessary. Sucrose loss during storage is a major concern for the sugar industry because the accumulation of invert sugar and byproducts severely affect sucrose manufacturing. This loss is mainly due to ongoing respiration, but changes in cell wall composition and pathogen infestation also contribute. While some varieties can cope better during storage, the underlying molecular mechanisms are currently undiscovered. We applied integrative transcriptomics on six varieties exhibiting different levels of sucrose loss during storage. Already prior to storage, well storable varieties were characterized by a higher number of parenchyma cells, a smaller cell area, and a thinner periderm. Supporting these findings, transcriptomics identified changes in genes involved in cell wall modifications. After 13 weeks of storage, over 900 differentially expressed genes were detected between well and badly storable varieties, mainly in the category of defense response but also in carbohydrate metabolism and the phenylpropanoid pathway. These findings were confirmed by gene co-expression network analysis where hub genes were identified as main drivers of invert sugar accumulation and sucrose loss. Our data provide insight into transcriptional changes in sugar beet roots during storage resulting in the characterization of key pathways and hub genes that might be further used as markers to improve pathogen resistance and storage properties.

EgPHI-1, a PHOSPHATE-INDUCED-1 gene from Eucalyptus globulus, is involved in shoot growth, xylem fiber length and secondary cell wall properties

MAIN CONCLUSION

EgPHI-1 is a member of PHI-1/EXO/EXL protein family. Its overexpression in tobacco resulted in changes in biomass partitioning, xylem fiber length, secondary cell wall thickening and composition, and lignification. Here, we report the functional characterization of a PHOSPHATE-INDUCED PROTEIN 1 homologue showing differential expression in xylem cells from Eucalyptus species of contrasting phenotypes for wood quality and growth traits. Our results indicated that this gene is a member of the PHI-1/EXO/EXL family. Analysis of the promoter cis-acting regulatory elements and expression responses to different treatments revealed that the Eucalyptus globulus PHI-1 (EgPHI-1) is transcriptionally regulated by auxin, cytokinin, wounding and drought. EgPHI-1 overexpression in transgenic tobacco changed the partitioning of biomass, favoring its allocation to shoots in detriment of roots. The stem of the transgenic plants showed longer xylem fibers and reduced cellulose content, while the leaf xylem had enhanced secondary cell wall thickness. UV microspectrophotometry of individual cell wall layers of fibers and vessels has shown that the transgenic plants exhibit differences in the lignification of S2 layer in both cell types. Taken together, the results suggest that EgPHI-1 mediates the elongation of secondary xylem fibers, secondary cell wall thickening and composition, and lignification, making it an attractive target for biotechnological applications in forestry and biofuel crops.

MiMIF-2 Effector of Meloidogyne incognita Exhibited Enzyme Activities and Potential Roles in Plant Salicylic Acid Synthesis

Plant-parasitic nematodes secrete a series of effectors to promote parasitism by modulating host immunity, but the detailed molecular mechanism is ambiguous. Animal parasites secrete macrophage migration inhibitory factor (MIF)-like proteins for evasion of host immune systems, in which their biochemical activities play essential roles. Previous research demonstrated that MiMIF-2 effector was secreted by Meloidogyne incognita and modulated host immunity by interacting with annexins. In this study, we show that MiMIF-2 had tautomerase activity and protected nematodes against H₂O₂ damage. MiMIF-2 expression not only decreased the amount of H₂O₂ generation during nematode infection in Arabidopsis, but also suppressed Bax-induced cell death by inhibiting reactive oxygen species burst in Nicotiana benthamiana. Further, RNA-seq transcriptome analysis and RT-qPCR showed that the expression of some heat-shock proteins was down regulated in MiMIF-2 transgenic Arabidopsis. After treatment with flg22, RNA-seq transcriptome analysis indicated that the differentially expressed genes in MiMIF-2 expressing Arabidopsis were pointed to plant hormone signal transduction, compound metabolism and plant defense. RT-qPCR and metabolomic results confirmed that salicylic acid (SA) related marker genes and SA content were significantly decreased. Our results provide a comprehensive understanding of how MiMIF-2 modulates plant immunity and broaden knowledge of the intricate relationship between M. incognita and host plants.

Transcriptomic Response to Water Deficit Reveals a Crucial Role of Phosphate Acquisition in a Drought-Tolerant Common Bean Landrace

Drought is one of the most critical factors limiting legume crop productivity. Understanding the molecular mechanisms of drought tolerance in the common bean is required to improve the yields of this important crop under adverse conditions. In this work, RNA-seq analysis was performed to compare the transcriptome profiles of drought-stressed and well-irrigated plants of a previously characterized drought-tolerant common bean landrace. The analysis revealed responses related with the abscisic acid signaling, including downregulation of a phosphatase 2C (PP2C) and an abscisic acid-8' hydroxylase, and upregulation of several key transcription factors and genes involved in cell wall remodeling, synthesis of osmoprotectants, protection of photosynthetic apparatus, and downregulation of genes involved in cell expansion. The results also highlighted a significant proportion of differentially expressed genes related to phosphate starvation response. In addition, the moderate detrimental effects of drought in the biomass of these tolerant plants were abolished by the addition of phosphate, thus indicating that, besides the ABA-mediated response, acquisition of phosphate could be crucial for the drought tolerance of this common bean genotype. These results provided information about the mechanisms involved in drought response of common bean response that could be useful for enhancing the drought tolerance of this important crop legume.

Overexpression of native Musa-miR397 enhances plant biomass without compromising abiotic stress tolerance in banana

Plant micro RNAs (miRNAs) control growth, development and stress tolerance but are comparatively unexplored in banana, whose cultivation is threatened by abiotic stress and nutrient deficiencies. In this study, a native Musa-miR397 precursor harboring 11 copper-responsive GTAC motifs in its promoter element was identified from banana genome. Musa-miR397 was significantly upregulated (8-10) fold in banana roots and leaves under copper deficiency, correlating with expression of root copper deficiency marker genes such as Musa-COPT and Musa-FRO2. Correspondingly, target laccases were significantly downregulated (>-2 fold), indicating miRNA-mediated silencing for Cu salvaging. No significant expression changes in the miR397-laccase module were observed under iron stress. Musa-miR397 was also significantly upregulated (>2 fold) under ABA, MV and heat treatments but downregulated under NaCl stress, indicating universal stress-responsiveness. Further, Musa-miR397 overexpression in banana significantly increased plant growth by 2-3 fold compared with wild-type but did not compromise tolerance towards Cu deficiency and NaCl stress. RNA-seq of transgenic and wild type plants revealed modulation in expression of 71 genes related to diverse aspects of growth and development, collectively promoting enhanced biomass. Summing up, our results not only portray Musa-miR397 as a candidate for enhancing plant biomass but also highlight it at the crossroads of growth-defense trade-offs.

A role for SPEECHLESS in the integration of leaf stomatal patterning with the growth vs disease trade-off in poplar

Occurrence of stomata on both leaf surfaces (amphistomaty) promotes higher stomatal conductance and photosynthesis while simultaneously increasing exposure to potential disease agents in black cottonwood (Populus trichocarpa). A genome-wide association study (GWAS) with 2.2M single nucleotide polymorphisms generated through whole-genome sequencing found 280 loci associated with variation in adaxial stomatal traits, implicating genes regulating stomatal development and behavior. Strikingly, numerous loci regulating plant growth and response to biotic and abiotic stresses were also identified. The most significant locus was a poplar homologue of SPEECHLESS (PtSPCH1). Individuals possessing PtSPCH1 alleles associated with greater adaxial stomatal density originated primarily from environments with shorter growing seasons (e.g. northern latitudes, high elevations) or with less precipitation. PtSPCH1 was expressed in developing leaves but not developing stem xylem. In developing leaves, RNA sequencing showed patterns of coordinated expression between PtSPCH1 and other GWAS-identified genes. The breadth of our GWAS results suggests that the evolution of amphistomaty is part of a larger, complex response in plants. Suites of genes underpin this response, retrieved through genetic association to adaxial stomata, and show coordinated expression during development. We propose that the occurrence of amphistomaty in P. trichocarpa involves PtSPCH1 and reflects selection for supporting rapid growth over investment in immunity.