qPrimerDB: a thermodynamics-based gene-specific qPCR primer database for 147 organisms

Members of WRKY Group III Transcription Factors Are Important in Mite Infestation in Strawberry (Fragaria × ananassa Duch.)

Strawberry is frequently attacked by mites, which directly affects the yield and quality of this fruit species. The WRKY Group III transcription factors (TFs) play an important role in plant tolerance to biotic sources of stress, such as pathogens and insect pests. In this study, six Group III WRKY TFs (FaWRKY25, FaWRKY31, FaWRKY32, FaWRKY43, FaWRKY44, and FaWRKY45) were identified in strawberry. A phylogenetic analysis showed that the six WRKY III TFs were divided into two clades and all had a conserved WRKYGQK domain and the C-X7-C-X23-H-T-C zinc finger motif. An interaction network analysis revealed that FaWRKY44 was co-expressing with FaWRKY25 and FaWRKY45. The expression patterns showed that the WRKY Group III genes responded to plant hormones and mite infestation in strawberry. To further verify the role of FaWRKY25 in plant resistance to mites, we cloned the FaWRKY25 gene and overexpressed it in transgenic plants. An in vivo subcellular localization analysis indicated that the FaWRKY25 protein was localized in the nucleus. Fewer mites were also detected on the wild-type plants than on FaWRKY25-overexpressing transgenic plants, suggesting that FaWRKY25 negatively regulates the resistance of strawberry to mites. The present study advances our understanding on a potential target that mites use to manipulate host plant defenses.

Comprehensive Identification and Expression Analysis of the Multidrug and Toxic Compound Extrusion (MATE) Gene Family in Brachypodium distachyon

The Multidrug and Toxic Compound Extrusion (MATE) proteins serve as pivotal transporters responsible for the extrusion of metabolites, thereby playing a significant role in both plant development and the detoxification of toxins. The MATE gene family within the Brachypodium distachyon, which is an important model organism of the Poaceae family, remains largely unexplored. Here, a comprehensive identification and analysis of MATE genes that complement B. distachyon were conducted. The BdMATE genes were systematically categorized into five distinct groups, predicated on an assessment of their phylogenetic affinities and protein structure. Furthermore, our investigation revealed that dispersed duplication has significantly contributed to the expansion of the BdMATE genes, with tandem and segmental duplications showing important roles, suggesting that the MATE genes in Poaceae species have embarked on divergent evolutionary trajectories. Examination of ω values demonstrated that BdMATE genes underwent purifying selection throughout the evolutionary process. Furthermore, collinearity analysis has confirmed a high conservation of MATE genes between B. distachyon and rice. The cis-regulatory elements analysis within BdMATEs promoters, coupled with expression patterns, suggests that BdMATEs play important roles during plant development and in response to phytohormones. Collectively, the findings presented establish a foundational basis for the subsequent detailed characterization of the MATE gene family members in B. distachyon.

Analysis of the Mechanism of Wood Vinegar and Butyrolactone Promoting Rapeseed Growth and Improving Low-Temperature Stress Resistance Based on Transcriptome and Metabolomics

Rapeseed is an important oil crop in the world. Wood vinegar could increase the yield and abiotic resistance of rapeseed. However, little is known about the underlying mechanisms of wood vinegar or its valid chemical components on rapeseed. In the present study, wood vinegar and butyrolactone (γ-Butyrolactone, one of the main components of wood vinegar) were applied to rapeseed at the seedling stage, and the molecular mechanisms of wood vinegar that affect rapeseed were studied by combining transcriptome and metabolomic analyses. The results show that applying wood vinegar and butyrolactone increases the biomass of rapeseed by increasing the leaf area and the number of pods per plant, and enhances the tolerance of rapeseed under low temperature by reducing membrane lipid oxidation and improving the content of chlorophyll, proline, soluble sugar, and antioxidant enzymes. Compared to the control, 681 and 700 differentially expressed genes were in the transcriptional group treated with wood vinegar and butyrolactone, respectively, and 76 and 90 differentially expressed metabolites were in the metabolic group. The combination of transcriptome and metabolomic analyses revealed the key gene-metabolic networks related to various pathways. Our research shows that after wood vinegar and butyrolactone treatment, the amino acid biosynthesis pathway of rapeseed may be involved in mediating the increase in rapeseed biomass, the proline metabolism pathway of wood vinegar treatment may be involved in mediating rapeseed's resistance to low-temperature stress, and the sphingolipid metabolism pathway of butyrolactone treatment may be involved in mediating rapeseed's resistance to low-temperature stress. It is suggested that the use of wood vinegar or butyrolactone are new approaches to increasing rapeseed yield and low-temperature resistance.

qPrimerDB 2.0: an updated comprehensive gene-specific qPCR primer database for 1172 organisms

High-quality primer design is essential for the success of all polymerase chain reaction (PCR)-based experiments. We previously developed a thermodynamics-based gene-specific quantitative PCR (qPCR) primer database for 147 organisms, which has been used extensively in gene expression studies. However, the number of organisms and the imperfection of function in the database limits its potential applications. Here, we improved the functionality of qPrimerDB to create a more comprehensive primer resource. Specifically, we (i) developed an improved primer design tool, qPrimer, building upon the previous qPrimerDB pipeline, to enhance the efficiency and simplicity of genome-scale qPCR primer design; (ii) pre-computed qPCR primer resources from 1 308 genomes of 1172 organisms and (iii) introduced a complete system for identifying, designing, checking, marking, and submitting qPCR primers. qPrimerDB 2.0 is freely available at https://qprimerdb.biodb.org. The qPrimer source code is available at https://github.com/swu1019lab/qPrimer.

Genome-wide characterization of Related to ABI3/VP1 transcription factors among U's triangle Brassica species reveals a negative role for BnaA06.RAV3L in seed size

The transcription factors Related to ABI3/VP1 (RAV) are crucial for various plant processes and stress responses. Although the U's triangle Brassica species genomes have been released, the knowledge regarding the RAV family is still limited. In this study, we identified 123 putative RAV genes across the six U's triangle Brassica species (Brassica rapa, 14; Brassica oleracea, 14; Brassica nigra, 13; Brassica carinata, 27; Brassica juncea, 28; Brassica napus, 27). Phylogenetic analysis categorized them into three groups. The RAV genes exhibited diversity in both functional and structural aspects, particularly in gene structure and cis-acting elements within their promoters. The expression analysis revealed that BnaRAV genes in Group 1/2 exhibited diverse expression patterns across various tissues, while those in Group 3 did not show expression except for BnaRAV3L-2 and BnaRAV3L-6, which were exclusively expressed in seeds. Furthermore, the seed-specific expression of BnaA06. RAV3L (BnaRAV3L-2) was confirmed through promoter-GUS staining. Subcellular localization studies demonstrated that BnaA06.RAV3L is localized to the nucleus. The overexpression of BnaA06. RAV3L in Arabidopsis led to a remarkable inhibition of seed-specific traits such as seed width, seed length, seed area, and seed weight. This study provides insights into the functional evolution of the RAV gene family in U triangle Brassica species. It establishes a foundation for uncovering the molecular mechanisms underlying the negative role of RAV3L in seed development.

ZbMYB111 Expression Positively Regulates ZbUFGT-Mediated Anthocyanin Biosynthesis in Zanthoxylum bungeanum with the Involvement of ZbbHLH2

Anthocyanin (ACN)-derived pigmentation in the red Zanthoxylum bungeanum peel is an essential commercial trait. Therefore, exploring the metabolic regulatory networks involved in peel ACN levels in this species is crucial for improving its quality. However, its underlying transcriptional regulatory mechanisms are still unknown. This transcriptomic and bioinformatics study not only discovered a new TF (ZbMYB111) as a potential regulator for ACN biosynthesis in Z. bungeanum peel, but also deciphered the underlying molecular mechanisms of ACN biosynthesis. Overexpression of ZbMYB111 and flavonoid 3-O-glucosyltransferase (ZbUFGT) induced ACN accumulation in both Z. bungeanum peels and callus along with Arabidopsis thaliana and tobacco flowers, whereas their silencing impaired ACN biosynthesis. Therefore, the dual-luciferase reporter, yeast-one-hybrid, and electrophoretic mobility shift assays showed that ZbMYB111 directly interacted with the ZbUFGT promoter to activate its expression. This diverted the secondary metabolism toward the ACN pathway, thereby promoting ACN accumulation.

GhRCD1 promotes cotton tolerance to cadmium by regulating the GhbHLH12-GhMYB44-GhHMA1 transcriptional cascade

Heavy metal pollution poses a significant risk to human health and wreaks havoc on agricultural productivity. Phytoremediation, a plant-based, environmentally benign, and cost-effective method, is employed to remove heavy metals from contaminated soil, particularly in agricultural or heavy metal-sensitive lands. However, the phytoremediation capacity of various plant species and germplasm resources display significant genetic diversity, and the mechanisms underlying these differences remain hitherto obscure. Given its potential benefits, genetic improvement of plants is essential for enhancing their uptake of heavy metals, tolerance to harmful levels, as well as overall growth and development in contaminated soil. In this study, we uncover a molecular cascade that regulates cadmium (Cd2+) tolerance in cotton, involving GhRCD1, GhbHLH12, GhMYB44, and GhHMA1. We identified a Cd2+-sensitive cotton T-DNA insertion mutant with disrupted GhRCD1 expression. Genetic knockout of GhRCD1 by CRISPR/Cas9 technology resulted in reduced Cd2+ tolerance in cotton seedlings, while GhRCD1 overexpression enhanced Cd2+ tolerance. Through molecular interaction studies, we demonstrated that, in response to Cd2+ presence, GhRCD1 directly interacts with GhbHLH12. This interaction activates GhMYB44, which subsequently activates a heavy metal transporter, GhHMA1, by directly binding to a G-box cis-element in its promoter. These findings provide critical insights into a novel GhRCD1-GhbHLH12-GhMYB44-GhHMA1 regulatory module responsible for Cd2+ tolerance in cotton. Furthermore, our study paves the way for the development of elite Cd2+-tolerant cultivars by elucidating the molecular mechanisms governing the genetic control of Cd2+ tolerance in cotton.

Exploring silique number in Brassica napus L.: Genetic and molecular advances for improving yield

Silique number is a crucial yield-related trait for the genetic enhancement of rapeseed (Brassica napus L.). The intricate molecular process governing the regulation of silique number involves various factors. Despite advancements in understanding the mechanisms regulating silique number in Arabidopsis (Arabidopsis thaliana) and rice (Oryza sativa), the molecular processes involved in controlling silique number in rapeseed remain largely unexplored. In this review, we identify candidate genes and review the roles of genes and environmental factors in regulating rapeseed silique number. We use genetic regulatory networks for silique number in Arabidopsis and grain number in rice to uncover possible regulatory pathways and molecular mechanisms involved in regulating genes associated with rapeseed silique number. A better understanding of the genetic network regulating silique number in rapeseed will provide a theoretical basis for the genetic improvement of this trait and genetic resources for the molecular breeding of high-yielding rapeseed.

Transcription factor NtNAC56 regulates jasmonic acid-induced leaf senescence in tobacco

Leaf senescence is a vital aspect of plant physiology and stress responses and is induced by endogenous factors and environmental cues. The plant-specific NAC (NAM, ATAF1/2, CUC2) transcription factor family influences growth, development, and stress responses in Arabidopsis (Arabidopsis thaliana) and other species. However, the roles of NACs in tobacco (Nicotiana tabacum) leaf senescence are still unclear. Here, we report that NtNAC56 regulates leaf senescence in tobacco. Transgenic plants overexpressing NtNAC56 (NtNAC56-OE) showed induction of senescence-related genes and exhibited early senescence and lower chlorophyll content compared to wild-type (WT) plants and the Ntnac56-19 mutant. In addition, root development and seed germination were inhibited in the NtNAC56-OE lines. Transmission electron microscopy observations accompanied by physiological and biochemical assays revealed that NtNAC56 overexpression triggers chloroplast degradation and reactive oxygen species accumulation in tobacco leaves. Transcriptome analysis demonstrated that NtNAC56 activates leaf senescence-related genes and jasmonic acid (JA) biosynthesis pathway genes. In addition, the JA content of NtNAC56-OE plants was higher than in WT plants, and JA treatment induced NtNAC56 expression. We performed DNA affinity purification sequencing to identify direct targets of NtNAC56, among which we focused on LIPOXYGENASE 5 (NtLOX5), a key gene in JA biosynthesis. A dual-luciferase reporter assay and a yeast one-hybrid assay confirmed that NtNAC56 directly binds to the TTTCTT motif in the NtLOX5 promoter. Our results reveal a mechanism whereby NtNAC56 regulates JA-induced leaf senescence in tobacco and provide a strategy for genetically manipulating leaf senescence and plant growth.

Stomatal density suppressor PagSDD1 is a "generalist" gene that promotes plant growth and improves water use efficiency

The serine protease SDD1 regulates stomatal density, but its potential impact on plant vegetative growth is unclear. Our study reveals a substantial upregulation of SDD1 in triploid poplar apical buds and leaves, suggesting its possible role in their growth regulation. We cloned PagSDD1 from poplar 84 K (Populus alba × P. glandulosa) and found that overexpression in poplar, soybean, and lettuce led to decreased leaf stomatal density. Furthermore, PagSDD1 represses PagEPF1, PagEPF2, PagEPFL9, PagSPCH, PagMUTE, and PagFAMA expression. In contrast, PagSDD1 promotes the expression of its receptors, PagTMM and PagERECTA. PagSDD1-OE poplars showed stronger drought tolerance than wild-type poplars. Simultaneously, PagSDD1-OE poplar, soybean, and lettuce had vegetative growth advantages. RNA sequencing revealed a significant upregulation of genes PagLHCB2.1 and PagGRF5, correlating positively with photosynthetic rate, and PagCYCA3;4 and PagEXPA8 linked to cell division and differentiation in PagSDD1-OE poplars. This increase promoted leaf photosynthesis, boosted auxin and cytokinin accumulation, and enhanced vegetative growth. SDD1 overexpression can increase the biomass of poplar, soybean, and lettuce by approximately 70, 176, and 155 %, respectively, and increase the water use efficiency of poplar leaves by over 52 %, which is of great value for the molecular design and breeding of plants with growth and water-saving target traits.

Heterologous overexpression of heat shock protein 20 genes of different species of yellow Camellia in Arabidopsis thaliana reveals their roles in high calcium resistance

Yellow Camellia (Camellia sect. chrysantha) is a rare ornamental plant and an important germplasm resource globally. Camellia nitidissima thrives in normal acidic soils, while Camellia limonia can adapt to the calcareous soils found in karst areas. Our previous study on the karst adaptation of yellow camellias revealed that the expression levels of heat shock protein 20(HSP20) were higher in Camellia limonia than in Camellia nitidissima. However, the functions of the HSP20 gene of Camellia limonia remain unclear to data. In this study, the HSP20 genes of Camellia limonia (ClHSP20-OE lines) and Camellia. nitidissima (CnHSP20-OE lines) were cloned and overexpressed heterologously in Arabidopsis thaliana. Additionally, we overexpressed the HSP20 gene of Arabidopsis (AtHSP20-OE lines) was also overexpressed, and the T-DNA inserted mutants (athspmutant lines) were also used to determine the functions of HSP20 genes. Under high calcium stress, the chlorophyll, nitrogen, water content and humidity of leaves were increased in ClHSP20-OE lines, while those of other lines were declined. The size of the stomatal apertures, stomatal conductance, and the photosynthetic efficiency of ClHSP20-OE lines were higher than those of the other lines. However, the accumulation of H2O2 and O2- in the leaves of ClHSP20-OE lines was the lowest among all the lines. Energy spectrum scanning revealed that the percentage of calcium on the surfaces of the leaves of ClHSP20-OE lines was relatively low, while that of athspmutant lines was the highest. The ClHSP20 gene can also affected soil humidity and the contents of soil nitrogen, phosphorus, and potassium. Transcriptome analysis revealed that the expressions of FBA5 and AT5G10770 in ClHSP20-OE lines was significantly up-regulated compared to that of CnHSP20-OE lines. Compared to that of athspmutant lines, the expressions of DREB1A and AT3G30460 was significantly upregulated in AtHSP20-OE lines, and the expression of POL was down-regulated. Our findings suggest that the HSP20 gene plays a crucial role in maintained photosynthetic rate and normal metabolism by regulating the expression of key genes under high-calcium stress. This study elucidates the mechanisms underlying the karst adaptation in Camellia. limonia and provides novel insights for future research on karst plants.

Regulation of micro- and small-exon retention and other splicing processes by GRP20 for flower development

Pre-mRNA splicing is crucial for gene expression and depends on the spliceosome and splicing factors. Plant exons have an average size of ~180 nucleotides and typically contain motifs for interactions with spliceosome and splicing factors. Micro exons (<51 nucleotides) are found widely in eukaryotes and in genes for plant development and environmental responses. However, little is known about transcript-specific regulation of splicing in plants and about the regulators for micro exon splicing. Here we report that glycine-rich protein 20 (GRP20) is an RNA-binding protein and required for splicing of ~2,100 genes including those functioning in flower development and/or environmental responses. Specifically, GRP20 is required for micro-exon retention in transcripts of floral homeotic genes; these micro exons are conserved across angiosperms. GRP20 is also important for small-exon (51-100 nucleotides) splicing. In addition, GRP20 is required for flower development. Furthermore, GRP20 binds to poly-purine motifs in micro and small exons and a spliceosome component; both RNA binding and spliceosome interaction are important for flower development and micro-exon retention. Our results provide new insights into the mechanisms of micro-exon retention in flower development.

Loss function of NtGA3ox1 delays flowering through impairing gibberellins metabolite synthesis in Nicotiana tabacum

Flowering time, plays a crucial role in tobacco ecological adaptation besides its substantial influence on tobacco production and leaf quality. Meanwhile, it is sensitive to biotic or abiotic challenges. The plant hormones Gibberellins (GAs), controlling a number of metabolic processes, govern plants growth and development. In this study, we created a late flowering mutant HG14 through knocking out NtGA3ox1 by CRISPR/Cas9. It took around 13.0 and 12.1 days longer to budding and flowering compared to wild type Honghuadajinyuan. Nearly all of the evaluated agronomic characters deteriorated in HG14, showing slower growth and noticeably shorter and narrower leaves. We found that NtGA3ox was more prevalent in flowers through quantitative reverse transcription PCR analysis. Transcriptome profiling detected 4449, 2147, and 4567 differently expressed genes at the budding, flowering, and mature stages, respectively. The KEGG pathway enrichment analysis identified the plant-pathogen interaction, plant hormone signal transduction pathway, and MAPK signaling pathway are the major clusters controlled by NtGA3ox1 throughout the budding and flowering stages. Together with the abovementioned signaling pathway, biosynthesis of monobactam, metabolism of carbon, pentose, starch, and sucrose were enriched at the mature stage. Interestingly, 108 up- and 73 down- regulated DEGs, impairing sugar metabolism, diterpenoid biosynthesis, linoleic and alpha-linolenic acid metabolism pathway, were continuously detected accompanied with the development of HG14. This was further evidenced by the decreasing content of GA metabolites such as GA4 and GA7, routine chemicals, alkaloids, amino acids, and organic acids Therefore, we discovered a novel tobacco flowering time gene NtGA3ox1 and resolved its regulatory network, which will be beneficial to the improvement of tobacco varieties.

Comprehensive identification and functional characterization of GhpPLA gene family in reproductive organ development

BACKGROUND

Phospholipases As (PLAs) are acyl hydrolases that catalyze the release of free fatty acids in phospholipids and play multiple functions in plant growth and development. The three families of PLAs are: PLA1, PLA2 (sPLA), and patatin-related PLA (pPLA). The diverse functions that pPLAs play in the growth and development of a broad range of plants have been demonstrated by prior studies.

METHODS

Genome-wide analysis of the pPLA gene family and screening of genes for expression verification and gene silencing verification were conducted. Additionally, pollen vitality testing, analysis of the pollen expression pattern, and the detection of POD, SOD, CAT, MDA, and H2O2 were performed.

RESULT

In this study, 294 pPLAs were identified from 13 plant species, including 46 GhpPLAs that were divided into three subfamilies (I-III). Expression patterns showed that the majority of GhpPLAs were preferentially expressed in the petal, pistil, anther, and ovule, among other reproductive organs. Particularly, GhpPLA23 and GhpPLA44, were found to be potentially important for the reproductive development of G. hirsutum. Functional validation was demonstrated by VIGS which showed that reduced expression levels of GhpPLA23 and GhpPLA44 in the silenced plants were associated with a decrease in pollen activity. Moreover, a substantial shift in ROS and ROS scavengers and a considerable increase in POD, CAT, SOD, and other physiological parameters was found out in these silenced plants. Our results provide plausibility to the hypothesis that GhpPLA23 and GhpPLA44 had a major developmental impact on cotton reproductive systems. These results also suggest that pPLAs are important for G. hirsutum's reproductive development and suggest that they could be employed as potential genes for haploid induction.

CONCLUSIONS

The findings of the present research indicate that pPLA genes are essential for the development of floral organs and sperm cells in cotton. Consequently, this family might be important for the reproductive development of cotton and possibly for inducing the plant develop haploid progeny.

LESION MIMIC MUTANT 1 confers basal resistance to Sclerotinia sclerotiorum in rapeseed via a salicylic acid-dependent pathway

Rapeseed (Brassica napus) is a major edible oilseed crop consumed worldwide. However, its yield is seriously affected by infection from the broad-spectrum non-obligate pathogen Sclerotinia sclerotiorum due to a lack of highly resistant germplasm. Here, we identified a Sclerotinia-resistant and light-dependent lesion mimic mutant from an ethyl methanesulfonate-mutagenized population of the rapeseed inbred Zhongshuang 11 (ZS11) named lesion mimic mutant 1 (lmm1). The phenotype of lmm1 is controlled by a single recessive gene, named LESION MIMIC MUTANT 1 (LMM1), which mapped onto chromosome C04 by bulked segregant analysis within a 2.71-Mb interval. Histochemical analysis indicated that H2O2 strongly accumulated and cell death occurred around the lesion mimic spots. Among 877 differentially expressed genes (DEGs) between ZS11 and lmm1 leaves, 188 DEGs were enriched in the defense response, including 95 DEGs involved in systemic acquired resistance, which is consistent with the higher salicylic acid levels in lmm1. Combining bulked segregant analysis and transcriptome analysis, we identified a significantly up-regulated gene, BnaC4.PR2, which encodes β-1,3-glucanase, as the candidate gene for LMM1. Overexpression of BnaC4.PR2 may induce a reactive oxygen species burst to trigger partial cell death and systemic acquired resistance. Our study provides a new genetic resource for S. sclerotiorum resistance as well as new insights into disease resistance breeding in B. napus.

Comparative transcriptional and co-expression network analysis of two upland cotton accessions with extreme phenotypic differences reveals molecular mechanisms of fiber development

Introduction

Upland cotton (Gossypium hirsutum) is the main source of natural fiber in the global textile industry, and thus its fiber quality and yield are important parameters. In this study, comparative transcriptomics was used to analyze differentially expressed genes (DEGs) due to its ability to effectively screen candidate genes during the developmental stages of cotton fiber. However, research using this method is limited, particularly on fiber development. The aim of this study was to uncover the molecular mechanisms underlying the whole period of fiber development and the differences in transcriptional levels.

Methods

Comparative transcriptomes are used to analyze transcriptome data and to screen for differentially expressed genes. STEM and WGCNA were used to screen for key genes involved in fiber development. qRT-PCR was performed to verify gene expression of selected DEGs and hub genes.

Results

Two accessions of upland cotton with extreme phenotypic differences, namely EZ60 and ZR014121, were used to carry out RNA sequencing (RNA-seq) on fiber samples from different fiber development stages. The results identified 704, 376, 141, 269, 761, and 586 genes that were upregulated, and 1,052, 476, 355, 259, 702, and 847 genes that were downregulated at 0, 5, 10, 15, 20, and 25 days post anthesis, respectively. Similar expression patterns of DEGs were monitored using short time-series expression miner (STEM) analysis, and associated pathways of DEGs within profiles were investigated. In addition, weighted gene co-expression network analysis (WGCNA) identified five key modules in fiber development and screened 20 hub genes involved in the development of fibers.

Discussion

Through the annotation of the genes, it was found that the excessive expression of resistance-related genes in the early fiber development stages affects the fiber yield, whereas the sustained expression of cell elongation-related genes is critical for long fibers. This study provides new information that can be used to improve fibers in newly developed upland cotton genotypes.

Brassica rapa Nitrate Transporter 2 (BrNRT2) Family Genes, Identification, and Their Potential Functions in Abiotic Stress Tolerance

Nitrate transporter 2 (NRT2) proteins play vital roles in both nitrate (NO3-) uptake and translocation as well as abiotic stress responses in plants. However, little is known about the NRT2 gene family in Brassica rapa. In this study, 14 NRT2s were identified in the B. rapa genome. The BrNRT2 family members contain the PLN00028 and MATE_like superfamily domains. Cis-element analysis indicated that regulatory elements related to stress responses are abundant in the promoter sequences of BrNRT2 genes. BrNRT2.3 expression was increased after drought stress, and BrNRT2.1 and BrNRT2.8 expression were significantly upregulated after salt stress. Furthermore, protein interaction predictions suggested that homologs of BrNRT2.3, BrNRT2.1, and BrNRT2.8 in Arabidopsis thaliana may interact with the known stress-regulating proteins AtNRT1.1, AtNRT1.5, and AtNRT1.8. In conclusion, we suggest that BrNRT2.1, BrNRT2.3, and BrNRT2.8 have the greatest potential for inducing abiotic stress tolerance. Our findings will aid future studies of the biological functions of BrNRT2 family genes.

Real-Time Quantitative PCR: Primer Design, Reference Gene Selection, Calculations and Statistics

Real-time quantitative PCR is a technique that can measure the content of the target nucleic acid sequence of interest in a given sample. It is mainly divided into absolute and relative quantitative methods. The relative quantification is mainly used in gene expressions for functional genomic and transcriptome studies. However, to use this technology accurately, there are some key points to master. First, specific primers need to be designed to ensure amplification of the gene of interest (GOI). Second, the appropriate reference gene or reference gene combination has to be selected. Finally, scientific gene expression level calculations and statistics are required to obtain accurate results. Therefore, this work proposes a workflow for relative quantitative PCR and introduces the relevant points so that beginners can better understand and use this technology.

Identification of a Rice Leaf Width Gene Narrow Leaf 22 (NAL22) through Genome-Wide Association Study and Gene Editing Technology

Rice leaf width (RLW) is a crucial determinant of photosynthetic area. Despite the discovery of several genes controlling RLW, the underlying genetic architecture remains unclear. In order to better understand RLW, this study conducted a genome-wide association study (GWAS) on 351 accessions from the rice diversity population II (RDP-II). The results revealed 12 loci associated with leaf width (LALW). In LALW4, we identified one gene, Narrow Leaf 22 (NAL22), whose polymorphisms and expression levels were associated with RLW variation. Knocking out this gene in Zhonghua11, using CRISPR/Cas9 gene editing technology, resulted in a short and narrow leaf phenotype. However, seed width remained unchanged. Additionally, we discovered that the vein width and expression levels of genes associated with cell division were suppressed in nal22 mutants. Gibberellin (GA) was also found to negatively regulate NAL22 expression and impact RLW. In summary, we dissected the genetic architecture of RLW and identified a gene, NAL22, which provides new loci for further RLW studies and a target gene for leaf shape design in modern rice breeding.

Effect of a suitable treatment period on the genetic transformation efficiency of the plant leaf disc method

BACKGROUND

Agrobacterium tumefaciens-mediated leaf disc genetic transformation is an important way to achieve transgenics or gene editing. Ensuring stable and efficient genetic transformation is still an important problem in modern biology. It is assumed that the difference in the development status of genetic transformation cells of receptor materials is the main reason for the difference and instability of genetic transformation efficiency; the stable and efficient genetic transformation rate can be obtained by defining the appropriate treatment period of the receptor material and applying genetic transformation in a timely manner.

RESULTS

Based on these assumptions, we studied and established an efficient and stable Agrobacterium-mediated plant transformation system with hybrid poplar (Populus alba × Populus glandulosa, 84 K) leaves, stem segments and tobacco leaves as the research objects. There were differences in the development process of leaf bud primordial cells from different explants, and the genetic transformation efficiency was significantly related to the cell development stage of the in vitro cultured materials. Among them, the genetic transformation rate of poplar and tobacco leaves was the highest on the 3rd and 2nd day of culture, reaching 86.6% and 57.3%, respectively. The genetic transformation rate of poplar stem segments was the highest on the 4th day of culture, reaching 77.8%. The best treatment period was from the development of leaf bud primordial cells to the S phase of the cell cycle. The number of cells detected using flow cytometry and 5-ethynyl-2'-deoxyuridine (EdU) staining, the expression of cell cycle-related protein CDKB1; 2, CDKD1; 1, CYCA3; 4, CYCD1; 1, CYCD3; 2, CYCD6; 1, and CYCH; 1 of explants, and morphological changes of explants can be used as indicators to determine the appropriate treatment period for genetic transformation.

CONCLUSIONS

Our study provides a new and universal set of methods and characteristics to identify the S phase of the cell cycle and apply genetic transformation treatments at the appropriate time. Our results are of great significance for improving the efficiency and stability of plant leaf disc genetic transformation.

Comparative transcriptome and co-expression network analysis revealed the genes associated with senescence and polygalacturonase activity involved in pod shattering of rapeseed

BACKGROUND

The pod shattering (PS) trait negatively affects the crop yield in rapeseed especially under dry conditions. To better understand the trait and cultivate higher resistance varieties, it's necessary to identify key genes and unravel the PS mechanism thoroughly.

RESULTS

In this study, we conducted a comparative transcriptome analysis between two materials significantly different in silique shatter resistance lignin deposition and polygalacturonase (PG) activity. Here, we identified 10,973 differentially expressed genes at six pod developmental stages. We found that the late pod development stages might be crucial in preparing the pods for upcoming shattering events. GO enrichment results from K-means clustering and weighed gene correlation network analysis (WGCNA) both revealed senescence-associated genes play an important role in PS. Two hub genes Bna.A05ABI5 and Bna.C03ERF/AP2-3 were selected from the MEyellow module, which possibly regulate the PS through senescence-related mechanisms. Further investigation found that senescence-associated transcription factor Bna.A05ABI5 upregulated the expression of SAG2 and ERF/AP2 to control the shattering process. In addition, the upregulation of Bna.C03ERF/AP2-3 is possibly involved in the transcription of downstream SHP1/2 and LEA proteins to trigger the shattering mechanism. We also analyzed the PS marker genes and found Bna.C07SHP1/2 and Bna.PG1/2 were significantly upregulated in susceptible accession. Furthermore, the role of auxin transport by Bna.WAG2 was also observed, which could reduce the PG activity to enhance the PS resistance through the cell wall loosening process.

CONCLUSION

Based on comparative transcriptome evaluation, this study delivers insights into the regulatory mechanism primarily underlying the variation of PS in rapeseed. Taken together, these results provide a better understanding to increase the yield of rapeseed by reducing the PS through better engineered crops.

Genome-wide identification and characterization of polycomb repressive complex 2 core components in upland cotton (Gossypium hirsutum L.)

BACKGROUND

The evolutionarily conserved Polycomb Repressive Complex 2 (PRC2) plays a vital role in epigenetic gene repression by depositing tri-methylation on lysine residue K27 of histone H3 (H3K27me3) at the target loci, thus participating in diverse biological processes. However, few reports about PRC2 are available in plant species with large and complicated genomes, like cotton.

RESULTS

Here, we performed a genome-wide identification and comprehensive analysis of cotton PRC2 core components, especially in upland cotton (Gossypium hirsutum). Firstly, a total of 8 and 16 PRC2 core components were identified in diploid and tetraploid cotton species, respectively. These components were classified into four groups, E(z), Su(z)12, ESC and p55, and the members in the same group displayed good collinearity, similar gene structure and domain organization. Next, we cloned G. hirsutum PRC2 (GhPRC2) core components, and found that most of GhPRC2 proteins were localized in the nucleus, and interacted with each other to form multi-subunit complexes. Moreover, we analyzed the expression profile of GhPRC2 genes. The transcriptome data and quantitative real-time PCR (qRT-PCR) assays indicated that GhPRC2 genes were ubiquitously but differentially expressed in various tissues, with high expression levels in reproductive organs like petals, stamens and pistils. And the expressions of several GhPRC2 genes, especially E(z) group genes, were responsive to various abiotic and biotic stresses, including drought, salinity, extreme temperature, and Verticillium dahliae (Vd) infection.

CONCLUSION

We identified PRC2 core components in upland cotton, and systematically investigated their classifications, phylogenetic and synteny relationships, gene structures, domain organizations, subcellular localizations, protein interactions, tissue-specific and stresses-responsive expression patterns. Our results will provide insights into the evolution and composition of cotton PRC2, and lay the foundation for further investigation of their biological functions and regulatory mechanisms.

Novel TLR7/8 agonists promote activation of HIV-1 latent reservoirs and human T and NK cells

Antiretroviral therapy can successfully suppress HIV-1 replication to undetectable levels but fails to eliminate latent and persistent HIV-1 reservoirs. Recent studies have focused on the immunomodulatory agents such as Toll-like receptor 7 and 8 (TLR7 and TLR8) capable of activating, thereby rendering the reservoir susceptible to antiretroviral inhibition and immune recognition and elimination. In this context, this study focused on generating a diverse repertoire of TLR7/8 agonists to identify more potent candidates for activating latent HIV-1 and immune cells' response. Through combinational strategies of computer-aided design and biological characterization, 159 pyrido [3,2-d] pyrimidine and pyridine-2-amine-based derivatives were synthesized. Of which, two TLR7/8 dual and one TLR8-specific agonists with exceptionally high potency in activating HIV-1 latent reservoirs in cell lines and PBMCs of patients with persistent and durable virologic controls were identified. Particularly, these agonists appeared to enhance NK and T cells activity, which were correlated with the degree of surface activation markers. The outcome of this study highlights the remarkable potential of TLR7/8 agonists in simultaneously activating HIV-1 from the latently infected cells and augmenting immune effector cells.

Stepwise Optimization of Real-Time RT-PCR Analysis

Quantitative real-time reverse transcription PCR (qRT-PCR) analysis has been used routinely to quantify gene expression levels. Primer design and the optimization of qRT-PCR parameters are critical for the accuracy and reproducibility of qRT-PCR analysis. Computational tool-assisted primer design often overlooks the presence of homologous sequences of the gene of interest and the sequence similarities between homologous genes in a plant genome. This sometimes results in skipping the optimization of qRT-PCR parameters due to the false confidence in the quality of the designed primers. Here we present a stepwise optimization protocol for single nucleotide polymorphisms (SNPs)-based sequence-specific primer design and sequential optimization of primer sequences, annealing temperatures, primer concentrations, and cDNA concentration range for each reference and target gene. The goal of this optimization protocol is to achieve a standard cDNA concentration curve with an R² ≥ 0.9999 and efficiency (E) = 100 ± 5% for the best primer pair of each gene, which serves as the prerequisite for using the 2^-ΔΔCT method for data analysis.

Enrichment of health-promoting lutein and zeaxanthin in tomato fruit through metabolic engineering

Carotenoids constitute a large group of natural pigments widely distributed in nature. These compounds not only provide fruits and flowers with distinctive colors, but also have significant health benefits for humans. Lutein and zeaxanthin, both oxygen-containing carotenoids, are considered to play vital roles in promoting ocular development and maintaining eye health. However, humans and mammals cannot synthesize these carotenoid derivatives, which can only be taken from certain fruits or vegetables. Here, by introducing four endogenous synthetic genes, SlLCYE, SlLCYB, SlHYDB, and SlHYDE under fruit-specific promoters, we report the metabolic engineering of lutein/zeaxanthin biosynthesis in tomato fruit. Transgenic lines overexpression of one (SlLCYE), two (SlLCYE and SlLCYB; SlLCYB and SlHYDB), and all these four synthetic genes re-established the lutein/zeaxanthin biosynthetic pathways in the ripe tomato fruit and thus resulted in various types of carotenoid riched lines. Metabolic analyses of these engineered tomato fruits showed the strategy involved expression of SlLCYE tends to produce α-carotene and lutein, as well as a higher content of β-carotene and zeaxanthin was detected in lines overexpressing SlLCYB. In addition, the different combinations of engineered tomatoes with riched carotenoids showed higher antioxidant capacity and were associated with a significantly extended shelf life during postharvest storage. This work provides a successful example of accurate metabolic engineering in tomato fruit, suggesting the potential utility for synthetic biology to improve agronomic traits in crops. These biofortified tomato fruits could be also exploited as new research subjects for studying the health benefits of carotenoid derivatives.

Genome-Wide Identification and Expression Analysis of BraGLRs Reveal Their Potential Roles in Abiotic Stress Tolerance and Sexual Reproduction

Glutamate receptors (GLRs) are involved in multiple functions during the plant life cycle through affecting the Ca²⁺ concentration. However, GLRs in Brassica species have not yet been reported. In this study, 16 glutamate receptor-like channels (GLR) belonged to two groups were identified in the Brassica rapa (B. rapa) genome by bioinformatic analysis. Most members contain domains of ANF_receptor, Peripla_BP_6, Lig_chan, SBP_bac_3, and Lig_chan_Glu_bd that are closely related to glutamate receptor channels. This gene family contains many elements associated with drought stress, low temperature stress, methyl jasmonate (MeJA), salicylic acid (SA), and other stress resistance. Gene expression profiles showed that BraGLR genes were expressed in roots, stems, leaves, flowers, and siliques. BraGLR5 expression was elevated after drought stress in drought-sensitive plants. BraGLR1, BraGLR8, and BraGLR11 expression were significantly upregulated after salt stress. BraGLR3 expression is higher in the female sterile-line mutants than in the wild type. The expression levels of BraGLR6, BraGLR9, BraGLR12, and BraGLR13 were significantly higher in the male sterile-line mutants than in the wild type. The expression of most BraGLRs increased after self-pollination, with BraGLR9 exhibiting the greatest increase. These results suggest that BraGLRs play an important role in abiotic stress tolerance and sexual reproduction.

Genome-wide association study reveals a GLYCOGEN SYNTHASE KINASE 3 gene regulating plant height in Brassica napus

Rapeseed (Brassica napus) is an allotetraploid crop that is the main source of edible oils and feed proteins in the world. The ideal plant architecture breeding is a major objective of rapeseed breeding and determining the appropriate plant height is a key element of the ideal plant architecture. Therefore, this study aims to improve the understanding of the genetic controls underlying plant height. The plant heights of 230 rapeseed accessions collected worldwide were investigated in field experiments over two consecutive years in Wuhan, China. Whole-genome resequencing of these accessions yielded a total of 1,707,194 informative single nucleotide polymorphisms (SNPs) that were used for genome-wide association analysis (GWAS). GWAS and haplotype analysis showed that BnaA01g09530D, which encodes BRASSINOSTEROID-INSENSITIVE 2 and belongs to the GLYCOGEN SYNTHASE KINASE 3 (GSK3) family, was significantly associated with plant height in B. napus. Moreover, a total of 31 BnGSK3s with complete domains were identified from B. napus genome and clustered into four groups according to phylogenetic analysis, gene structure, and motif distribution. The expression patterns showed that BnGSK3s exhibited significant differences in 13 developmental tissues in B. napus, suggesting that BnGSK3s may be involved in tissue-specific development. Sixteen BnGSK3 genes were highly expressed the in shoot apical meristem, which may be related to plant height or architecture development. These results are important for providing new haplotypes of plant height in B. napus and for extending valuable genetic information for rapeseed genetic improvement of plant architecture.

Titanium nanoparticles activate a transcriptional response in Arabidopsis that enhances tolerance to low phosphate, osmotic stress and pathogen infection

Titanium is a ubiquitous element with a wide variety of beneficial effects in plants, including enhanced nutrient uptake and resistance to pathogens and abiotic stresses. While there is numerous evidence supporting the beneficial effects that Ti fertilization give to plants, there is little information on which genetic signaling pathways the Ti application activate in plant tissues. In this study, we utilize RNA-seq and ionomics technologies to unravel the molecular signals that Arabidopsis plants unleash when treated with Ti. RNA-seq analysis showed that Ti activates abscisic acid and salicylic acid signaling pathways and the expression of NUCLEOTIDE BINDING SITE-LEUCINE RICH REPEAT receptors likely by acting as a chemical priming molecule. This activation results in enhanced resistance to drought, high salinity, and infection with Botrytis cinerea in Arabidopsis. Ti also grants an enhanced nutritional state, even at suboptimal phosphate concentrations by upregulating the expression of multiple nutrient and membrane transporters and by modifying or increasing the production root exudates. Our results suggest that Ti might act similarly to the beneficial element Silicon in other plant species.

Genome-wide association studies reveal novel QTLs, QTL-by-environment interactions and their candidate genes for tocopherol content in soybean seed

Genome-wide association studies (GWAS) is an efficient method to detect quantitative trait locus (QTL), and has dissected many complex traits in soybean [Glycine max (L.) Merr.]. Although these results have undoubtedly played a far-reaching role in the study of soybean biology, environmental interactions for complex traits in traditional GWAS models are frequently overlooked. Recently, a new GWAS model, 3VmrMLM, was established to identify QTLs and QTL-by-environment interactions (QEIs) for complex traits. In this study, the GLM, MLM, CMLM, FarmCPU, BLINK, and 3VmrMLM models were used to identify QTLs and QEIs for tocopherol (Toc) content in soybean seed, including δ-Tocotrienol (δ-Toc) content, γ-Tocotrienol (γ-Toc) content, α-Tocopherol (α-Toc) content, and total Tocopherol (T-Toc) content. As a result, 101 QTLs were detected by the above methods in single-environment analysis, and 57 QTLs and 13 QEIs were detected by 3VmrMLM in multi-environment analysis. Among these QTLs, some QTLs (Group I) were repeatedly detected three times or by at least two models, and some QTLs (Group II) were repeatedly detected only by 3VmrMLM. In the two Groups, 3VmrMLM was able to correctly detect all known QTLs in group I, while good results were achieved in Group II, for example, 8 novel QTLs were detected in Group II. In addition, comparative genomic analysis revealed that the proportion of Glyma_max specific genes near QEIs was higher, in other words, these QEIs nearby genes are more susceptible to environmental influences. Finally, around the 8 novel QTLs, 11 important candidate genes were identified using haplotype, and validated by RNA-Seq data and qRT-PCR analysis. In summary, we used phenotypic data of Toc content in soybean, and tested the accuracy and reliability of 3VmrMLM, and then revealed novel QTLs, QEIs and candidate genes for these traits. Hence, the 3VmrMLM model has broad prospects and potential for analyzing the genetic structure of complex quantitative traits in soybean.

Transcriptomic analysis provides insight into defensive strategies in response to continuous cropping in strawberry (Fragaria × ananassa Duch.) plants

BACKGROUND

Strawberries are an important economic fruit crop world-wide. In strawberry cultivation, continuous cropping (CC) can seriously threaten yield and quality. However, our understanding of the gene expression changes in response to CC and during subsequent defense processes is limited. In this study, we analyzed the impact of CC on the transcriptome of strawberry roots using RNA-Seq technology to elucidate the effect of CC and the subsequent molecular changes.

RESULTS

We found that CC significantly affects the growth of strawberry plants. The transcriptome analysis identified 136 differentially expressed genes (DEGs), including 49 up-regulated and 87 down-regulated DEGs. A Gene Ontology (GO) analysis indicated that the up-regulated DEGs were mainly assigned to defense-related GO terms, and most down-regulated DEGs were assigned to nutrient-related GO terms. Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis revealed that the responsive DEGs were classified in a large number of important biological pathways, such as phenylalanine metabolism, starch and sucrose metabolism, phenylpropanoid biosynthesis, glutathione metabolism and plant-pathogen interaction. We also found that four WRKY transcription factors and three peroxidase genes involved in plant defense pathways were up-regulated in the roots of strawberry plants subjected to CC.

CONCLUSION

Several unigenes involved in plant defense processes, such as CNGCs, WRKY transcription factors, PR1, and peroxidase genes with highly variable expression levels between non-CC and CC treatments may be involved in the regulation of CC in strawberry. These results indicate that strawberry roots reallocate development resources to defense mechanisms in response to CC. This study will further deepen our understanding of the fundamental regulatory mechanisms of strawberry resource reallocation in response to CC.

Integrated analysis of transcriptome and metabolome reveals molecular mechanisms of salt tolerance in seedlings of upland rice landrace 17SM-19

Salt stress is a major abiotic stress that threatens global rice production. It is particularly important to improve salt tolerance in upland rice because of its growth environment. Upland rice landrace 17SM-19 with high salt tolerance was obtained from a previous study. In this study, an integrated analysis of transcriptome and metabolome was performed to determine the responses of the rice seedling to salt stress. When treated with 100 mm NaCl, the rice seedling growth was significantly inhibited at 5 d, with inhibition first observed in shoot dry weight (SDW). Changes in potassium (K⁺) content were associated with changes in SDW. In omics analyses, 1,900 differentially expressed genes (DEGs) and 659 differentially abundant metabolites (DAMs) were identified at 3 d after salt stress (DAS), and 1,738 DEGs and 657 DAMs were identified at 5 DAS. Correlation analyses between DEGs and DAMs were also conducted. The results collectively indicate that salt tolerance of upland rice landrace 17SM-19 seedlings involves many molecular mechanisms, such as those involved with osmotic regulation, ion balance, and scavenging of reactive oxygen species.

Overexpression of PagSTOMAGEN, a Positive Regulator of Stomatal Density, Promotes Vegetative Growth in Poplar

Poplar is an important fast-growing tree, and its photosynthetic capacity directly affects its vegetative growth. Stomatal density is closely related to photosynthetic capacity and growth characteristics in plants. Here, we isolated PagSTOMAGEN from the hybrid poplar (Populus alba × Populus glandulosa) clone 84K and investigated its biological function in vegetative growth. PagSTOMAGEN was expressed predominantly in young tissues and localized in the plasma membrane. Compared with wild-type 84K poplars, PagSTOMAGEN-overexpressing plants displayed an increased plant height, leaf area, internode number, basal diameter, biomass, IAA content, IPR content, and stomatal density. Higher stomatal density improved the net photosynthetic rate, stomatal conductance, intercellular CO₂ concentration, and transpiration rate in transgenic poplar. The differential expression of genes related to stomatal development showed a diverged influence of PagSTOMAGEN at different stages of stomatal development. Finally, transcriptomic analysis showed that PagSTOMAGEN affected vegetative growth by affecting the expression of photosynthesis and plant hormone-related genes (such as SAUR75, PQL2, PSBX, ERF1, GNC, GRF5, and ARF11). Taken together, our data indicate that PagSTOMAGEN could positively regulate stomatal density and increase the photosynthetic rate and plant hormone content, thereby promoting vegetative growth in poplar. Our study is of great significance for understanding the relationship between stoma, photosynthesis, and yield breeding in poplar.

Isolation and Characterization of an LBD Transcription Factor CsLBD39 from Tea Plant (Camellia sinensis) and Its Roles in Modulating Nitrate Content by Regulating Nitrate-Metabolism-Related Genes

Nitrate nitrogen is an important nitrogen source for tea plants’ growth and development. LBD transcription factors play important roles in response to the presence of nitrate in plants. The functional study of LBD transcription factors in tea plants remains limited. In this study, the LBD family gene CsLBD39 was isolated and characterized from tea plants. Sequence analysis indicated that CsLBD39 contained a highly conserved CX₂CX₆CX₃CX domain. The phylogenetic tree assay showed that CsLBD39 belonged to class II subfamily of the LBD family. CsLBD39 was highly expressed in flowers and root; we determined that its expression could be induced by nitrate treatment. The CsLBD39 protein was located in the nucleus and has transcriptional activation activity in yeast. Compared with the wild type, overexpression of CsLBD39 gene in Arabidopsis resulted in smaller rosettes, shorter main roots, reduced lateral roots and lower plant weights. The nitrate content and the expression levels of genes related to nitrate transport and regulation were decreased in transgenic Arabidopsis hosting CsLBD39 gene. Compared with the wild type, CsLBD39 overexpression in transgenic Arabidopsis had smaller cell structure of leaves, shorter diameter of stem cross section, and slender and compact cell of stem longitudinal section. Under KNO₃ treatment, the contents of nitrate, anthocyanins, and chlorophyll in leaves, and the content of nitrate in roots of Arabidopsis overexpressing CsLBD39 were reduced, the expression levels of nitrate transport and regulation related genes were decreased. The results revealed that CsLBD39 may be involved in nitrate signal transduction in tea plants as a negative regulator and laid the groundwork for future studies into the mechanism of nitrate response.

Overexpression of the Arabidopsis MACPF Protein AtMACP2 Promotes Pathogen Resistance by Activating SA Signaling

Immune response in plants is tightly regulated by the coordination of the cell surface and intracellular receptors. In animals, the membrane attack complex/perforin-like (MACPF) protein superfamily creates oligomeric pore structures on the cell surface during pathogen infection. However, the function and molecular mechanism of MACPF proteins in plant pathogen responses remain largely unclear. In this study, we identified an Arabidopsis MACP2 and investigated the responsiveness of this protein during both bacterial and fungal pathogens. We suggest that MACP2 induces programmed cell death, bacterial pathogen resistance, and necrotrophic fungal pathogen sensitivity by activating the biosynthesis of tryptophan-derived indole glucosinolates and the salicylic acid signaling pathway dependent on the activity of enhanced disease susceptibility 1 (EDS1). Moreover, the response of MACP2 mRNA isoforms upon pathogen attack is differentially regulated by a posttranscriptional mechanism: alternative splicing. In comparison to previously reported MACPFs in Arabidopsis, MACP2 shares a redundant but nonoverlapping role in plant immunity. Thus, our findings provide novel insights and genetic tools for the MACPF family in maintaining SA accumulation in response to pathogens in Arabidopsis.

Effect of Metabolic Changes in Aphis gossypii-Damaged Cotton Plants on Oviposition Preference and Larval Development of Subsequent Helicoverpa armigera

Aphis gossypii and Helicoverpa armigera are two important agricultural pests in cotton plants. However, whether early colonization of A. gossypii affects subsequent H. armigera is unknown. We implemented ecological experiments to reveal that A. gossypii-damaged cotton plants [Bacillus thuringiensis (Bt) and non-Bt] had a significant avoidance effect on the oviposition preference of H. armigera adults. However, A. gossypii-damaged cotton plants (non-Bt) increased the weight and pupation rate and reduced the mortality of H. armigera larvae. Transcriptomic and metabolomic analyses showed that 13 and 9 genes were significantly upregulated to be involved in salicylic acid (SA) and indole acetic acid (IAA) biosynthesis, and SA and IAA contents were significantly increased, respectively. However, 15 genes involved in jasmonic acid (JA) biosynthesis were significantly downregulated as a result of the antagonism of SA and JA. Moreover, there was significant upregulation in multiple genes involved in the biosynthesis of l-histidine, fructose, maltotetraose, melezitose, lecithin, stearidonic acid, and mannitol, in which metabolites were confirmed to promote the growth and development of H. armigera. Our study is a reference for investigating the evolutionary relationships and provides insights into implementing effective insect biocontrol between H. armigera and A. gossypii.

Ethylene Plays a Dual Role during Infection by Plasmodiophora brassicae of Arabidopsis thaliana

Plasmodiophora brassicae infection leads to hypertrophy of host roots and subsequent formation of galls, causing huge economic losses to agricultural producers of Cruciferae plants. Ethylene (ET) has been reported to play a vital role against necrotrophic pathogens in the classic immunity system. More clues suggested that the defense to pathogens in roots may be different from the acrial. The ET pathway may play a positive role in the infection of P. brassicae, as shown by recent transcriptome profiling. However, the molecular basis of ET remains poorly understood. In this study, we investigated the potential role of ethylene against P. brassicae infection in an ein3/eil1 double-mutant of Arabidopsis thaliana (A. thaliana). After infection, ein3/eil1 (Disease Index/DI: 93) showed more susceptibility compared with wild type (DI: 75). Then, we inoculated A. thaliana Columbia-0 (Col-0) with P. brassicae by 1-aminocyclopropane-1-carboxylic acid (ACC) and pyrazinamide (PZA), respectively. It was found that the symptoms of infected roots with ACC were more serious than those with PZA at 20 dpi (day post infection). However, the DI were almost the same in different treatments at 30 dpi. WRKY75 can be directly regulated by ET and was upregulated at 7 dpi with ACC, as shown by qRT-PCR. The wrky75-c mutant of A. thaliana (DI: 93.75) was more susceptible than the wild type in Arabidopsis. Thus, our work reveals the dual roles of ET in infection of P. brassicae and provides evidence of ET in root defense against pathogens.

Helicoverpa armigera herbivory negatively impacts Aphis gossypii populations via inducible metabolic changes

BACKGROUND

Helicoverpa armigera and Aphis gossypii are two important insect species that feed on cotton plants. These insects have distinct abilities to induce plant resistance and tolerate plant toxins, which results in interspecific competition imbalance that may be fatal to the low-tolerance A. gossypii and force these insects to develop avoidance behaviors and subsequently separate from their niche. We implemented ecological experiments to test the effects of H. armigera-induced plant resistance and behavioral avoidance in A. gossypii, and employed transcriptomics and metabolomics analyses to reveal changes in resistance genes and metabolites in plants.

RESULTS

Our results demonstrate that cotton plants induced by H. armigera cause significant inhibitory and avoidance effects on A. gossypii insect populations. Electrical penetration graph (EPG) analysis showed changes in plant resistance induced by H. armigera leading to a decreased feeding efficiency of A. gossypii. In addition, genes associated with jasmonic acid and ethylene signaling pathways were significantly up-regulated in cotton plants after H. armigera induction, which led to a significant up-regulation of metabolites inducing plant resistance. These observations were corroborated by bioactivity analysis on metabolites, which showed that jasmonic acid, gossypol and tannins have significant inhibitory effects on A. gossypii populations. In contrast, methylparaben is associated with avoidance behaviors on A. gossypii populations.

CONCLUSION

Our study suggests that the differences in the ability to induce plant resistance and tolerance between two non-predatory insects were lethal to low-tolerance A. gossypii insects, which might be a major factor determining their niche differentiation. This was further demonstrated by screening anti-insect and bio-hormonal metabolites. Our study provides a reference for investigating the evolutionary relationship between non-predatory insects and insights to implement effective insect biocontrol. © 2022 Society of Chemical Industry.

Transcriptome and Small RNA Sequencing Reveal the Mechanisms Regulating Harvest Index in Brassica napus

Harvest index (HI), the ratio of harvested seed weight to total aboveground biomass weight, is an economically critical value reflecting the convergence of complex agronomic traits. HI values in rapeseed (Brassica napus) remain much lower than in other major crops, and the underlying regulatory network is largely unknown. In this study, we performed mRNA and small RNA sequencing to reveal the mechanisms shaping HI in B. napus during the seed-filling stage. A total of 8,410 differentially expressed genes (DEGs) between high-HI and low-HI accessions in four tissues (silique pericarp, seed, leaves, and stem) were identified. Combining with co-expression network, 72 gene modules were identified, and a key gene BnaSTY46 was found to participate in retarded establishment of photosynthetic capacity to influence HI. Further research found that the genes involved in circadian rhythms and response to stimulus may play important roles in HI and that their transcript levels were modulated by differentially expressed microRNAs (DEMs), and we identified 903 microRNAs (miRNAs), including 46 known miRNAs and 857 novel miRNAs. Furthermore, transporter activity-related genes were critical to enhancing HI in good cultivation environments. Of 903 miRNAs, we found that the bna-miR396-Bna.A06SRp34a/Bna.A01EMB3119 pair may control the seed development and the accumulation of storage compounds, thus contributing to higher HI. Our findings uncovered the underlying complex regulatory network behind HI and offer potential approaches to rapeseed improvement.

SlWRKY35 positively regulates carotenoid biosynthesis by activating the MEP pathway in tomato fruit

Carotenoids are vital phytonutrients widely recognised for their health benefits. Therefore, it is vital to thoroughly investigate the metabolic regulatory network underlying carotenoid biosynthesis and accumulation to open new leads towards improving their contents in vegetables and crops. The outcome of our study defines SlWRKY35 as a positive regulator of carotenoid biosynthesis in tomato. SlWRKY35 can directly activate the expression of the 1-deoxy-d-xylulose 5-phosphate synthase (SlDXS1) gene to reprogramme metabolism towards the 2-C-methyl-d-erythritol 4-phosphate (MEP) pathway, leading to enhanced carotenoid accumulation. We also show that the master regulator SlRIN directly regulates the expression of SlWRKY35 during tomato fruit ripening. Compared with the SlLCYE overexpression lines, coexpression of SlWRKY35 and SlLCYE can further enhance lutein production in transgenic tomato fruit, indicating that SlWRKY35 represents a potential target towards designing innovative metabolic engineering strategies for carotenoid derivatives. In addition to providing new insights into the metabolic regulatory network associated with tomato fruit ripening, our data define a new tool for improving fruit content in specific carotenoid compounds.

Genome-Wide Analysis of the Protein Phosphatase 2C Genes in Tomato

The plant protein phosphatase 2C (PP2C) plays an irreplaceable role in phytohormone signaling, developmental processes, and manifold stresses. However, information about the PP2C gene family in tomato (Solanum lycopersicum) is relatively restricted. In this study, a genome-wide investigation of the SlPP2C gene family was performed. A total of 92 SlPP2C genes were identified, they were distributed on 11 chromosomes, and all the SlPP2C proteins have the type 2C phosphatase domains. Based on phylogenetic analysis of PP2C genes in Arabidopsis, rice, and tomato, SlPP2C genes were divided into eight groups, designated A–H, which is also supported by the analyses of gene structures and protein motifs. Gene duplication analysis revealed that the duplication of whole genome and chromosome segments was the main cause of SLPP2Cs expansion. A total of 26 cis-elements related to stress, hormones, and development were identified in the 3 kb upstream region of these SlPP2C genes. Expression profile analysis revealed that the SlPP2C genes display diverse expression patterns in various tomato tissues. Furthermore, we investigated the expression patterns of SlPP2C genes in response to Ralstonia solanacearum infection. RNA-seq and qRT-PCR data reveal that nine SlPP2Cs are correlated with R. solanacearum. The above evidence hinted that SlPP2C genes play multiple roles in tomato and may contribute to tomato resistance to bacterial wilt. This study obtained here will give an impetus to the understanding of the potential function of SlPP2Cs and lay a solid foundation for tomato breeding and transgenic resistance to plant pathogens.

Activation of the VQ Motif-Containing Protein Gene VQ28 Compromised Nonhost Resistance of Arabidopsis thaliana to Phytophthora Pathogens

Nonhost resistance refers to resistance of a plant species to all genetic variants of a non-adapted pathogen. Such resistance has the potential to become broad-spectrum and durable crop disease resistance. We previously employed Arabidopsis thaliana and a forward genetics approach to identify plant mutants susceptible to the nonhost pathogen Phytophthora sojae, which resulted in identification of the T-DNA insertion mutant esp1 (enhanced susceptibility to Phytophthora). In this study, we report the identification of VQ motif-containing protein 28 (VQ28), whose expression was highly up-regulated in the mutant esp1. Stable transgenic A. thaliana plants constitutively overexpressing VQ28 compromised nonhost resistance (NHR) against P. sojae and P. infestans, and supported increased infection of P. parasitica. Transcriptomic analysis showed that overexpression of VQ28 resulted in six differentially expressed genes (DEGs) that are involved in the response to abscisic acid (ABA). High performance liquid chromatography-mass spectrometry (HPLC-MS) detection showed that the contents of endogenous ABA, salicylic acid (SA), and jasmonate (JA) were enriched in VQ28 overexpression lines. These findings suggest that overexpression of VQ28 may lead to an imbalance in plant hormone homeostasis. Furthermore, transient overexpression of VQ28 in Nicotiana benthamiana rendered plants more susceptible to Phytophthora pathogens. Deletion mutant analysis showed that the C-terminus and VQ-motif were essential for plant susceptibility. Taken together, our results suggest that VQ28 negatively regulates plant NHR to Phytophthora pathogens.

Nitrogen application and differences in leaf number retained after topping affect the tobacco (Nicotiana tabacum) transcriptome and metabolome

BACKGROUND

Agronomic treatments such as the application of nitrogen fertilizer and topping (removal of the inflorescence and top leaves) cause substantial changes in plant metabolism. To explore these changes, we conducted comparative transcriptomic and metabolomic analyses of leaves collected from four positions along the stem on plants exposed to two nitrogen doses and with different numbers of leaves retained after topping in tobacco (Nicotiana tabacum).

RESULTS

We identified 13,330 unique differentially expressed genes and 32 differentially abundant metabolites. Through RNA-seq and WGCNA analyze, we constructed 2 co-expression networks (green and blue) highly correlation to N application and leaf number retained, predicted a hub gene NtGER3 may play an important role in N metabolism related to amino acid (cysteine) through CK pathway in tobacco leaves, NtARFs may participated in modulating the auxin signal and N in bottom leaves and NtRAP2.12 as key gene involved in N regulation by ethylene pathway. What's more, our data prove C/N transformation and balance affect the "source - flow - sink" redistribution and remobilization in tobacco during growth and development process.

CONCLUSIONS

Overall, this comparative transcriptomics study provides novel insight into the complex molecular mechanisms underlying plant responses to different levels of nitrogen application and the number of leaves remaining after topping in plants.

qPrimerDB: A Powerful and User-Friendly Database for qPCR Primer Design

Real-time quantitative polymerase chain reaction (qPCR) is a powerful tool for analyzing and quantifying gene expression, and its primer design is the first and most important step. In order to improve the efficiency and effectiveness of primer design, we designed a database qPrimerDB, based on thermodynamics gene-specific for multispecies qPCR primers design.In this chapter, we explained the working principle of the database and detailed the step-by-step practical steps with examples. The valuable and time-saving qPrimerDB database is publicly accessible at http://biodb.swu.edu.cn/qprimerdb and will be routinely updated.

An optimized protocol for stepwise optimization of real-time RT-PCR analysis

Computational tool-assisted primer design for real-time reverse transcription (RT) PCR (qPCR) analysis largely ignores the sequence similarities between sequences of homologous genes in a plant genome. It can lead to false confidence in the quality of the designed primers, which sometimes results in skipping the optimization steps for qPCR. However, the optimization of qPCR parameters plays an essential role in the efficiency, specificity, and sensitivity of each gene's primers. Here, we proposed an optimized approach to sequentially optimizing primer sequences, annealing temperatures, primer concentrations, and cDNA concentration range for each reference (and target) gene. Our approach started with a sequence-specific primer design that should be based on the single-nucleotide polymorphisms (SNPs) present in all the homologous sequences for each of the reference (and target) genes under study. By combining the efficiency calibrated and standard curve methods with the 2^-ΔΔCt method, the standard cDNA concentration curve with a logarithmic scale was obtained for each primer pair for each gene. As a result, an R² ≥ 0.9999 and the efficiency (E) = 100 ± 5% should be achieved for the best primer pair of each gene, which serve as the prerequisite for using the 2^-ΔΔCt method for data analysis. We applied our newly developed approach to identify the best reference genes in different tissues and at various inflorescence developmental stages of Tripidium ravennae, an ornamental and biomass grass, and validated their utility under varying abiotic stress conditions. We also applied this approach to test the expression stability of six reference genes in soybean under biotic stress treatment with Xanthomonas axonopodis pv. glycines (Xag). Thus, these case studies demonstrated the effectiveness of our optimized protocol for qPCR analysis.

Co-expression network and comparative transcriptome analysis for fiber initiation and elongation reveal genetic differences in two lines from upland cotton CCRI70 RIL population

Upland cotton is the most widely planted for natural fiber around the world, and either lint percentage (LP) or fiber length (FL) is the crucial component tremendously affecting cotton yield and fiber quality, respectively. In this study, two lines MBZ70-053 and MBZ70-236 derived from G. hirsutum CCRI70 recombinant inbred line (RIL) population presenting different phenotypes in LP and FL traits were chosen to conduct RNA sequencing on ovule and fiber samples, aiming at exploring the differences of molecular and genetic mechanisms during cotton fiber initiation and elongation stages. As a result, 249/128, 369/206, 4296/1198 and 3547/2129 up-/down- regulated differentially expressed genes (DGEs) in L2 were obtained at -3, 0, 5 and 10 days post-anthesis (DPA), respectively. Seven gene expression profiles were discriminated using Short Time-series Expression Miner (STEM) analysis; seven modules and hub genes were identified using weighted gene co-expression network analysis. The DEGs were mainly enriched into energetic metabolism and accumulating as well as auxin signaling pathway in initiation and elongation stages, respectively. Meanwhile, 29 hub genes were identified as 14-3-3ω , TBL35, GhACS, PME3, GAMMA-TIP, PUM-7, etc., where the DEGs and hub genes revealed the genetic and molecular mechanisms and differences during cotton fiber development.

Metabolite Profiling and Transcriptome Analysis Provide Insight into Seed Coat Color in Brassica juncea

The allotetraploid species Brassica juncea (mustard) is grown worldwide as oilseed and vegetable crops; the yellow seed-color trait is particularly important for oilseed crops. Here, to examine the factors affecting seed coat color, we performed a metabolic and transcriptomic analysis of yellow- and dark-seeded B. juncea seeds. In this study, we identified 236 compounds, including 31 phenolic acids, 47 flavonoids, 17 glucosinolates, 38 lipids, 69 other hydroxycinnamic acid compounds, and 34 novel unknown compounds. Of these, 36 compounds (especially epicatechin and its derivatives) accumulated significantly different levels during the development of yellow- and dark-seeded B. juncea. In addition, the transcript levels of BjuDFR, BjuANS,BjuBAN, BjuTT8, and BjuTT19 were closely associated with changes to epicatechin and its derivatives during seed development, implicating this pathway in the seed coat color determinant in B. juncea. Furthermore, we found numerous variations of sequences in the TT8A genes that may be associated with the stability of seed coat color in B. rapa, B. napus, and B. juncea, which might have undergone functional differentiation during polyploidization in the Brassica species. The results provide valuable information for understanding the accumulation of metabolites in the seed coat color of B. juncea and lay a foundation for exploring the underlying mechanism.

Resolving diurnal dynamics of the chloroplastic glutathione redox state in Arabidopsis reveals its photosynthetically derived oxidation

Plants are subjected to fluctuations in light intensity, and this might cause unbalanced photosynthetic electron fluxes and overproduction of reactive oxygen species (ROS). Electrons needed for ROS detoxification are drawn, at least partially, from the cellular glutathione (GSH) pool via the ascorbate-glutathione cycle. Here, we explore the dynamics of the chloroplastic glutathione redox potential (chl-EGSH) using high-temporal-resolution monitoring of Arabidopsis (Arabidopsis thaliana) lines expressing the reduction-oxidation sensitive green fluorescent protein 2 (roGFP2) in chloroplasts. This was carried out over several days under dynamic environmental conditions and in correlation with PSII operating efficiency. Peaks in chl-EGSH oxidation during dark-to-light and light-to-dark transitions were observed. Increasing light intensities triggered a binary oxidation response, with a threshold around the light saturating point, suggesting two regulated oxidative states of the chl-EGSH. These patterns were not affected in npq1 plants, which are impaired in non-photochemical quenching. Oscillations between the two oxidation states were observed under fluctuating light in WT and npq1 plants, but not in pgr5 plants, suggesting a role for PSI photoinhibition in regulating the chl-EGSH dynamics. Remarkably, pgr5 plants showed an increase in chl-EGSH oxidation during the nights following light stresses, linking daytime photoinhibition and nighttime GSH metabolism. This work provides a systematic view of the dynamics of the in vivo chloroplastic glutathione redox state during varying light conditions.

Identification and candidate gene screening of qCIR9.1, a novel QTL associated with anther culturability in rice (Oryza sativa L.)

A novel QTL, qCIR9.1, that controls callus induction rate in anther culture was identified on chromosome 9 in rice, and based on RNA-seq data, Os09g0551600 was the most promising candidate gene. Anther culture, a doubled haploid (DH) technique, has become an important technology in many plant-breeding programmes. Although anther culturability is the key factor in this technique, its genetic mechanisms in rice remain poorly understood. In this study, we mapped quantitative trait loci (QTLs) responsible for anther culturability by using 192 recombinant inbred lines (RILs) derived from YZX (Oryza sativa ssp. indica) × 02428 (Oryza sativa ssp. japonica) and a high-density bin map. A total of eight QTLs for anther culturability were detected in three environments. Among these QTLs, a novel major QTL for callus induction rate (CIR) named qCIR9.1 was repeatedly mapped to a ~ 100 kb genomic interval on chromosome 9 and explained 8.39-14.14% of the phenotypic variation. Additionally, RNA sequencing (RNA-seq) was performed for the parents (YZX and 02428), low- (L-Pool) and high-CIR RILs (H-Pool) after 16 and 26 days of culture. By using the RNA of the bulked RILs for background normalization, the number of differentially expressed genes (DEGs) both between the parents and between the bulked RILs after 26 days of culture was drastically reduced to only 78. Among these DEGs, only one gene, Os09g0551600, encoding a high-mobility group (HMG) protein, was located in the candidate region of qCIR9.1. qRT-PCR analysis of Os09g0551600 showed the same results as RNA-seq, and the expression of this gene was decreased in the low-callus-induction parent (YZX) and L-Pool. Our results provide a foundational step for further cloning of qCIR9.1 and will be very useful for improving anther culturability in rice.

Genome-Wide Analysis of the YABBY Transcription Factor Family in Rapeseed (Brassica napus L.)

The YABBY family of plant-specific transcription factors play important regulatory roles during the development of leaves and floral organs, but their functions in Brassica species are incompletely understood. Here, we identified 79 YABBY genes from Arabidopsis thaliana and five Brassica species (B. rapa, B. nigra, B. oleracea, B. juncea, and B. napus). A phylogenetic analysis of YABBY proteins separated them into five clusters (YAB1–YAB5) with representatives from all five Brassica species, suggesting a high degree of conservation and similar functions within each subfamily. We determined the gene structure, chromosomal location, and expression patterns of the 21 BnaYAB genes identified, revealing extensive duplication events and gene loss following polyploidization. Changes in exon–intron structure during evolution may have driven differentiation in expression patterns and functions, combined with purifying selection, as evidenced by K_a/K_s values below 1. Based on transcriptome sequencing data, we selected nine genes with high expression at the flowering stage. qRT-PCR analysis further indicated that most BnaYAB family members are tissue-specific and exhibit different expression patterns in various tissues and organs of B. napus. This preliminary study of the characteristics of the YABBY gene family in the Brassica napus genome provides theoretical support and reference for the later functional identification of the family genes.