Comparative transcriptome analysis of soybean response to bean pyralid larvae

The soybean plasma membrane GmDR1 protein conferring broad-spectrum disease and pest resistance regulates several receptor kinases and NLR proteins

Overexpression of Glycine max disease resistant 1 (GmDR1) exhibits broad-spectrum resistance against Fusarium virguliforme, Heterodera glycines (soybean cyst nematode), Tetranychus urticae (Koch) (spider mites), and Aphis glycines Matsumura (soybean aphids) in soybean. To understand the mechanisms of broad-spectrum immunity mediated by GmDR1, the transcriptomes of a strong and a weak GmDR1-overexpressor following treatment with chitin, a pathogen- and pest-associated molecular pattern (PAMP) common to these organisms, were investigated. The strong and weak GmDR1-overexpressors exhibited altered expression of 6098 and 992 genes, respectively, as compared to the nontransgenic control following chitin treatment. However, only 192 chitin- and 115 buffer-responsive genes exhibited over two-fold changes in expression levels in both strong and weak GmDR1-overexpressors as compared to the control. MapMan analysis of the 192 chitin-responsive genes revealed 64 biotic stress-related genes, of which 53 were induced and 11 repressed as compared to the control. The 53 chitin-induced genes include nine genes that encode receptor kinases, 13 encode nucleotide-binding leucine-rich repeat (NLR) receptor proteins, seven encode WRKY transcription factors, four ethylene response factors, and three MYB-like transcription factors. Investigation of a subset of these genes revealed three receptor protein kinases, seven NLR proteins, and one WRKY transcription factor genes that are induced following F. virguliforme and H. glycines infection. The integral plasma membrane GmDR1 protein most likely recognizes PAMPs including chitin and activates transcription of genes encoding receptor kinases, NLR proteins and defense-related genes. GmDR1 could be a pattern recognition receptor that regulates the expression of several NLRs for expression of PAMP-triggered immunity and/or priming the effector triggered immunity.

Transcriptome analysis reveals the potential mechanism of the response to scale insects in Camellia sasanqua Thunb

BACKGROUND

Camellia sasanqua Thunb. is an essential woody ornamental plant. Our continuous observation found that scale insects often infest C. sasanqua all year round in Kunming, China, resulting in poor growth. Scientifically preventing and controlling the infestation of scale insects should be paid attention to, and the mechanism of scale insects influencing C. sasanqua should be used as the research basis.

RESULTS

The scale insect was identified as Pseudaulacaspis sasakawai Takagi. We analyzed transcriptome sequencing data from leaves of C. sasanqua infested with scale insects. A total of 1320 genes were either up-regulated or down-regulated and differed significantly in response to scale insects. GO (Gene Ontology) annotation analysis showed that the pathway of catalytic activity, binding, membrane part, cell part, and cellular process were affected. KEGG (Kyoto Encyclopedia of Genes and Genomes) pathway analysis showed that most DEGs (differentially expressed genes) involved in plant hormone signal transduction, MAPK signaling pathway, flavonoid biosynthesis, tropane, piperidine and pyridine alkaloid biosynthesis. We also observed that the expression of galactose metabolism and carotenoid biosynthesis were significantly influenced. In addition, qRT-PCR (quantitative real-time PCR) validated the expression patterns of DEGs, which showed an excellent agreement with the transcriptome sequencing.

CONCLUSIONS

Our transcriptomic analysis revealed that the C. sasanqua had an intricate resistance strategy to cope with scale insect attacks. After sensing the attack signal of scale insects, C. sasanqua activated the early signal MAPK (mitogen-activated protein kinase) to activate further transcription factors and Auxin, ET, JA, ABA, and other plant hormone signaling pathways, ultimately leading to the accumulation of lignin, scopolin, flavonoids and other secondary metabolites, produces direct and indirect resistance to scale insects. Our results suggested that it provided some potential resources of defense genes that would benefit the following resistance breeding in C. sasanqua to scale insects.

Multi-omics assisted breeding for biotic stress resistance in soybean

Biotic stress is a critical factor limiting soybean growth and development. Soybean responses to biotic stresses such as insects, nematodes, fungal, bacterial, and viral pathogens are governed by complex regulatory and defense mechanisms. Next-generation sequencing has availed research techniques and strategies in genomics and post-genomics. This review summarizes the available information on marker resources, quantitative trait loci, and marker-trait associations involved in regulating biotic stress responses in soybean. We discuss the differential expression of related genes and proteins reported in different transcriptomics and proteomics studies and the role of signaling pathways and metabolites reported in metabolomic studies. Recent advances in omics technologies offer opportunities to reshape and improve biotic stress resistance in soybean by altering gene regulation and/or other regulatory networks. We suggest using 'integrated omics' to precisely understand how soybean responds to different biotic stresses. We also discuss the potential challenges of integrating multi-omics for the functional analysis of genes and their regulatory networks and the development of biotic stress-resistant cultivars. This review will help direct soybean breeding programs to develop resistance against different biotic stresses.

MePAL6 regulates lignin accumulation to shape cassava resistance against two-spotted spider mite

INTRODUCTION

The two-spotted spider mite (TSSM) is a devastating pest of cassava production in China. Lignin is considered as an important defensive barrier against pests and diseases, several genes participate in lignin biosynthesis, however, how these genes modulate lignin accumulation in cassava and shape TSSM-resistance is largely unknown.

METHODS

To fill this knowledge gap, while under TSSM infestation, the cassava lignin biosynthesis related genes were subjected to expression pattern analysis followed by family identification, and genes with significant induction were used for further function exploration.

RESULTS

Most genes involved in lignin biosynthesis were up-regulated when the mite-resistant cassava cultivars were infested by TSSM, noticeably, the MePAL gene presented the most vigorous induction among these genes. Therefore, we paid more attention to dissect the function of MePAL gene during cassava-TSSM interaction. Gene family identification showed that there are 6 MePAL members identified in cassava genome, further phylogenetic analysis, gene duplication, cis-elements and conserved motif prediction speculated that these genes may probably contribute to biotic stress responses in cassava. The transcription profile of the 6 MePAL genes in TSSM-resistant cassava cultivar SC9 indicated a universal up-regulation pattern. To further elucidate the potential correlation between MePAL expression and TSSM-resistance, the most strongly induced gene MePAL6 were silenced using virus-induced gene silencing (VIGS) assay, we found that silencing of MePAL6 in SC9 not only simultaneously suppressed the expression of other lignin biosynthesis genes such as 4-coumarate--CoA ligase (4CL), hydroxycinnamoyltransferase (HCT) and cinnamoyl-CoA reductase (CCR), but also resulted in decrease of lignin content. Ultimately, the suppression of MePAL6 in SC9 can lead to significant deterioration of TSSM-resistance.

DISCUSSION

This study accurately identified MePAL6 as critical genes in conferring cassava resistance to TSSM, which could be considered as promising marker gene for evaluating cassava resistance to insect pest.

Structure, Function, and Applications of Soybean Calcium Transporters

Glycine max is a calcium-loving crop. The external application of calcium fertilizer is beneficial to the increase of soybean yield. Indeed, calcium is a vital nutrient in plant growth and development. As a core metal ion in signaling transduction, calcium content is maintained in dynamic balance under normal circumstances. Now, eight transporters were found to control the uptake and efflux of calcium. Though these calcium transporters have been identified through genome-wide analysis, only a few of them were functionally verified. Therefore, in this study, we summarized the current knowledge of soybean calcium transporters in structural features, expression characteristics, roles in stress response, and prospects. The above results will be helpful in understanding the function of cellular calcium transport and provide a theoretical basis for elevating soybean yield.

Methyl jasmonate-induced resistance to Delia platura (Diptera: Anthomyiidae) in Lupinus mutabilis

BACKGROUND

Andean lupin (Lupinus mutabilis Sweet) is an important leguminous crop from South America with a high protein content. In Ecuador, lupin yields are severely affected by the infestation of Delia platura larvae on germinating seeds. The application of elicitor molecules with activity against herbivorous insects to control D. platura infestation constitutes an unexplored and promising alternative for chemical insecticides. In this study, methyl jasmonate (MeJA), hexanoic acid, menadione sodium bisulfite, and DL-β-aminobutyric acid were evaluated for their ability to induce resistance against D. platura in three commercial lupin cultivars.

RESULTS

Only seeds pretreated with MeJA significantly impaired insect performance during choice and no-choice assays. Additionally, fitness indicators such as seed germination and growth were not affected by MeJA treatment. To investigate the molecular mechanisms behind the MeJA-mediated resistance, RT-qPCR assays were performed. First, RT-qPCR reference genes were validated, showing that LmUBC was the most stable reference gene. Next, expression analysis over time revealed that MeJA application up-regulated the activity of the jasmonic acid biosynthetic genes LmLOX2 and LmAOS, together with other jasmonate-related defense genes, such as LmTPS1, LmTPS4, LmPI2, LmMBL, LmL/ODC, LmCSD1, and LmPOD.

CONCLUSION

This study indicates that MeJA can be used as an environmentally friendly elicitor molecule to protect Andean lupin from D. platura attack without fitness cost. MeJA application induces plant defense responses to insects in Andean lupin that may be modulated by the onset of terpenoid biosynthesis, proteinase inhibitors, lectins, polyamines, and antioxidative enzymes. © 2021 Society of Chemical Industry.

Methylome and transcriptome analyses of soybean response to bean pyralid larvae

BACKGROUND

Bean pyralid is one of the major leaf-feeding insects that affect soybean crops. DNA methylation can control the networks of gene expressions, and it plays an important role in responses to biotic stress. However, at present the genome-wide DNA methylation profile of the soybean resistance to bean pyralid has not been reported so far.

RESULTS

Using whole-genome bisulfite sequencing (WGBS) and RNA-sequencing (RNA-seq), we analyzed the highly resistant material (Gantai-2-2, HRK) and highly susceptible material (Wan82-178, HSK), under bean pyralid larvae feeding 0 h and 48 h, to clarify the molecular mechanism of the soybean resistance and explore its insect-resistant genes. We identified 2194, 6872, 39,704 and 40,018 differentially methylated regions (DMRs), as well as 497, 1594, 9596 and 9554 differentially methylated genes (DMGs) in the HRK0/HRK48, HSK0/HSK48, HSK0/HRK0 and HSK48/HRK48 comparisons, respectively. Through the analysis of global methylation and transcription, 265 differentially expressed genes (DEGs) were negatively correlated with DMGs, there were 34, 49, 141 and 116 negatively correlated genes in the HRK0/HRK48, HSK0/HSK48, HSK0/HRK0 and HSK48/HRK48, respectively. The MapMan cluster analysis showed that 114 negatively correlated genes were clustered in 24 pathways, such as protein biosynthesis and modification; primary metabolism; secondary metabolism; cell cycle, cell structure and component; RNA biosynthesis and processing, and so on. Moreover, CRK40; CRK62; STK; MAPK9; L-type lectin-domain containing receptor kinase VIII.2; CesA; CSI1; fimbrin-1; KIN-14B; KIN-14 N; KIN-4A; cytochrome P450 81E8; BEE1; ERF; bHLH25; bHLH79; GATA26, were likely regulatory genes involved in the soybean responses to bean pyralid larvae. Finally, 5 DMRs were further validated that the genome-wide DNA data were reliable through PS-PCR and 5 DEGs were confirmed the relationship between DNA methylation and gene expression by qRT-PCR. The results showed an excellent agreement with deep sequencing.

CONCLUSIONS

Genome-wide DNA methylation profile of soybean response to bean pyralid was obtained for the first time. Several specific DMGs which participated in protein kinase, cell and organelle, flavonoid biosynthesis and transcription factor were further identified to be likely associated with soybean response to bean pyralid. Our data will provide better understanding of DNA methylation alteration and their potential role in soybean insect resistance.

Disclosure of salicylic acid and jasmonic acid-responsive genes provides a molecular tool for deciphering stress responses in soybean

Hormones orchestrate the physiology of organisms. Measuring the activity of defense hormone-responsive genes can help understanding immune signaling and facilitate breeding for plant health. However, different from model species like Arabidopsis, genes that respond to defense hormones salicylic acid (SA) and jasmonic acid (JA) have not been disclosed in the soybean crop. We performed global transcriptome analyses to fill this knowledge gap. Upon exogenous application, endogenous levels of SA and JA increased in leaves. SA predominantly activated genes linked to systemic acquired resistance and defense signaling whereas JA mainly activated wound response-associated genes. In general, SA-responsive genes were activated earlier than those responding to JA. Consistent with the paradigm of biotrophic pathogens predominantly activating SA responses, free SA and here identified most robust SA marker genes GmNIMIN1, GmNIMIN1.2 and GmWRK40 were induced upon inoculation with Phakopsora pachyrhizi, whereas JA marker genes did not respond to infection with the biotrophic fungus. Spodoptera exigua larvae caused a strong accumulation of JA-Ile and JA-specific mRNA transcripts of GmBPI1, GmKTI1 and GmAAT whereas neither free SA nor SA-marker gene transcripts accumulated upon insect feeding. Our study provides molecular tools for monitoring the dynamic accumulation of SA and JA, e.g. in a given stress condition.

Overexpression of a plasma membrane protein generated broad-spectrum immunity in soybean

Plants fight-off pathogens and pests by manifesting an array of defence responses using their innate immunity mechanisms. Here we report the identification of a novel soybean gene encoding a plasma membrane protein, transcription of which is suppressed following infection with the fungal pathogen, Fusarium virguliforme. Overexpression of the protein led to enhanced resistance against not only against F. virguliforme, but also against spider mites (Tetranychus urticae, Koch), soybean aphids (Aphis glycines, Matsumura) and soybean cyst nematode (Heterodera glycines). We, therefore, name this protein as Glycine max disease resistance 1 (GmDR1; Glyma.10g094800). The homologues of GmDR1 have been detected only in legumes, cocoa, jute and cotton. The deduced GmDR1 protein contains 73 amino acids. GmDR1 is predicted to contain an ecto- and two transmembrane domains. Transient expression of the green fluorescent protein fused GmDR1 protein in soybean leaves showed that it is a plasma membrane protein. We investigated if chitin, a pathogen-associated molecular pattern (PAMP), common to all pathogen and pests considered in this study, can significantly enhance defence pathways among the GmDR1-overexpressed transgenic soybean lines. Chitin induces marker genes of the salicylic- and jasmonic acid-mediated defence pathways, but suppresses the defence pathway regulated by ethylene. Chitin induced SA- and JA-regulated defence pathways may be one of the mechanisms involved in generating broad-spectrum resistance among the GmDR1-overexpressed transgenic soybean lines against two serious pathogens and two pests including spider mites, against which no known resistance genes have been identified in soybean and among the most other crop species.

Influence of Perfluorobutanoic Acid (PFBA) on the Developmental Cycle and Damage Potential of the Beet Armyworm Spodoptera exigua (Hübner) (Insecta: Lepidoptera: Noctuidae)

Perfluorobutanoic acid (PFBA), one of the short-chain replacement perfluoroalkyl substances, has been shown to accumulate in plants. The potential of PFBA to modulate the developmental cycle of the beet armyworm, Spodoptera exigua, a polyphagous pest, was investigated. Second-instar larvae were fed with PFBA-spiked artificial diets and leaves from soybean plants grown with PFBA-spiked irrigation water. Spiked PFBA concentrations were 200 μg/kg for the artificial diet, whereas 405 to 15,190 ng/kg accumulated in the soybean leaves. The larvae fed with the PFBA-spiked diet showed a significant increase in weight gain compared with the controls over a 7-day exposure period. A similar weight gain trend was observed with larvae fed with the PFBA-containing soybean leaves, with the dose-response data fitting into a Brain-Cousens hormesis model with a 57% stimulation over controls. The artificial diet treatments showed 66.7% metamorphosed larva to pupa at 9 days after exposure (dpe) compared with 33.3% of the controls. The adult emergence at 16-dpe followed a similar trend with 57.7% and 33.3%, respectively, for the exposed and control groups. The duration of transition from larvae to adults was more symmetrical and 0.5 day faster for the exposed groups over controls. The beet armyworm caused more damage on leaves from the PFBA exposed plants in a nonmonotonic dose-response manner. The results suggest PFBA may have a stimulatory impact on some hormonal signaling pathways at low doses.

Common bean resistance to Xanthomonas is associated with upregulation of the salicylic acid pathway and downregulation of photosynthesis

BACKGROUND

Common bacterial blight (CBB) caused by Xanthomonas phaseoli pv. phaseoli and Xanthomonas citri pv. fuscans is one of the major threats to common bean crops (Phaseolus vulgaris L.). Resistance to CBB is particularly complex as 26 quantitative resistance loci to CBB have been described so far. To date, transcriptomic studies after CBB infection have been very scarce and the molecular mechanisms underlying susceptibility or resistance are largely unknown.

RESULTS

We sequenced and annotated the genomes of two common bean genotypes being either resistant (BAT93) or susceptible (JaloEEP558) to CBB. Reciprocal BLASTp analysis led to a list of 20,787 homologs between these genotypes and the common bean reference genome (G19833), which provides a solid dataset for further comparative analyses. RNA-Seq after inoculation with X. phaseoli pv. phaseoli showed that the susceptible genotype initiated a more intense and diverse biological response than the resistant genotype. Resistance was linked to upregulation of the salicylic acid pathway and downregulation of photosynthesis and sugar metabolism, while susceptibility was linked to downregulation of resistance genes and upregulation of the ethylene pathway and of genes involved in cell wall modification.

CONCLUSIONS

This study helps better understanding the mechanisms occurring during the early colonization phase of common bean by Xanthomonas and unveils new actors potentially important for resistance and susceptibility to CBB. We discuss the potential link between the pathways induced during bean colonization and genes induced by transcription activator-like effectors (TALEs), as illustrated in other Xanthomonas pathovars.

Common bean resistance to Xanthomonas is associated with upregulation of the salicylic acid pathway and downregulation of photosynthesis

Background

Common bacterial blight (CBB) caused by Xanthomonas phaseoli pv. phaseoli and Xanthomonas citri pv. fuscans is one of the major threats to common bean crops (Phaseolus vulgaris L.). Resistance to CBB is particularly complex as 26 quantitative resistance loci to CBB have been described so far. To date, transcriptomic studies after CBB infection have been very scarce and the molecular mechanisms underlying susceptibility or resistance are largely unknown.

Results

We sequenced and annotated the genomes of two common bean genotypes being either resistant (BAT93) or susceptible (JaloEEP558) to CBB. Reciprocal BLASTp analysis led to a list of 20,787 homologs between these genotypes and the common bean reference genome (G19833), which provides a solid dataset for further comparative analyses. RNA-Seq after inoculation with X. phaseoli pv. phaseoli showed that the susceptible genotype initiated a more intense and diverse biological response than the resistant genotype. Resistance was linked to upregulation of the salicylic acid pathway and downregulation of photosynthesis and sugar metabolism, while susceptibility was linked to downregulation of resistance genes and upregulation of the ethylene pathway and of genes involved in cell wall modification.

Conclusions

This study helps better understanding the mechanisms occurring during the early colonization phase of common bean by Xanthomonas and unveils new actors potentially important for resistance and susceptibility to CBB. We discuss the potential link between the pathways induced during bean colonization and genes induced by transcription activator-like effectors (TALEs), as illustrated in other Xanthomonas pathovars.

GenFam: A web application and database for gene family-based classification and functional enrichment analysis

Genome-scale studies using high-throughput sequencing (HTS) technologies generate substantial lists of differentially expressed genes under different experimental conditions. These gene lists need to be further mined to narrow down biologically relevant genes and associated functions in order to guide downstream functional genetic analyses. A popular approach is to determine statistically overrepresented genes in a user-defined list through enrichment analysis tools, which rely on functional annotations of genes based on Gene Ontology (GO) terms. Here, we propose a new computational approach, GenFam, which allows annotation, classification, and enrichment of genes based on their gene family, thus simplifying identification of candidate gene families and associated genes that may be relevant to the query. GenFam and its integrated database comprises of three hundred and eighty-four unique gene families and supports gene family analyses for sixty plant genomes. Four comparative case studies with plant species belonging to different clades and families were performed using GenFam which demonstrated its robustness and comprehensiveness over preexisting functional enrichment tools. To make it readily accessible for plant biologists, GenFam is available as a web-based application where users can input gene IDs and export enrichment results in both tabular and graphical formats. Users can also customize analysis parameters by choosing from the various statistical enrichment tests and multiple testing correction methods. Additionally, the web-based application, source code, and database are freely available to use and download. Website: http://mandadilab.webfactional.com/home/. Source code and database: http://mandadilab.webfactional.com/home/dload/.

Determination of the MiRNAs Related to Bean Pyralid Larvae Resistance in Soybean Using Small RNA and Transcriptome Sequencing

Soybean is one of the most important oil crops in the world. Bean pyralid is a major leaf-feeding insect of soybean. In order to screen out the functional genes and regulatory pathways related to the resistance for bean pyralid larvae, the small RNA and transcriptome sequencing were performed based on the highly resistant material (Gantai-2-2) and highly susceptible material (Wan 82-178) of soybean. The results showed that, when comparing 48 h feeding with 0 h feeding, 55 differentially expressed miRNAs were identified in Gantai-2-2 and 58 differentially expressed miRNAs were identified in Wan82-178. When comparing Gantai-2-2 with Wan82-178, 77 differentially expressed miRNAs were identified at 0 h feeding, and 70 differentially expressed miRNAs were identified at 48 h feeding. The pathway analysis of the predicted target genes revealed that the plant hormone signal transduction, RNA transport, protein processing in the endoplasmic reticulum, zeatin biosynthesis, ubiquinone and other terpenoid-quinone biosynthesis, and isoquinoline alkaloid biosynthesis may play important roles in soybean's defense against the stress caused by bean pyralid larvae. According to conjoint analysis of the miRNA/mRNA, a total of 20 differentially expressed miRNAs were negatively correlated with 26 differentially expressed target genes. The qRT-PCR analysis verified that the small RNA sequencing results were credible. According to the analyses of the differentially expressed miRNAs, we speculated that miRNAs are more likely to play key roles in the resistance to insects. Gma-miR156q, Gma-miR166u, Gma-miR166b, Gma-miR166j-3p, Gma-miR319d, Gma-miR394a-3p, Gma-miR396e, and so on-as well as their negatively regulated differentially expressed target genes-may be involved in the regulation of soybean resistance to bean pyralid larvae. These results laid a foundation for further in-depth research regarding the action mechanisms of insect resistance.

Transcriptome Profile Analysis of Winter Rapeseed (Brassica napus L.) in Response to Freezing Stress, Reveal Potentially Connected Events to Freezing Stress

Winter rapeseed is not only an important oilseed crop, but also a winter cover crop in Northern China, where its production was severely limited by freezing stress. As an overwinter crop, the production is severely limited by freezing stress. Therefore, understanding the physiological and molecular mechanism of winter rapeseed (Brassica napus L.) in freezing stress responses becomes essential for the improvement and development of freezing-tolerant varieties of Brassica napus. In this study, morphological, physiological, ultrastructure and transcriptome changes in the Brassica napus line "2016TS(G)10" (freezing-tolerance line) that was exposed to -2 °C for 0 h, 1 h, 3 h and 24 h were characterized. The results showed that freezing stress caused seedling dehydration, and chloroplast dilation and degradation. The content of malondialdehyde (MDA), proline, soluble protein and soluble sugars were increased, as well as the relative electrolyte leakage (REL) which was significantly increased at frozen 24 h. Subsequently, RNA-seq analysis revealed a total of 98,672 UniGenes that were annotated in Brassica napus and 3905 UniGenes were identified as differentially expressed genes after being exposed to freezing stress. Among these genes, 2312 (59.21%) were up-regulated and 1593 (40.79%) were down-regulated. Most of these DEGs were significantly annotated in the carbohydrates and energy metabolism, signal transduction, amino acid metabolism and translation. Most of the up-regulated DEGs were especially enriched in plant hormone signal transduction, starch and sucrose metabolism pathways. Transcription factor enrichment analysis showed that the AP2/ERF, WRKY and MYB families were also significantly changed. Furthermore, 20 DEGs were selected to validate the transcriptome profiles via quantitative real-time PCR (qRT-PCR). In conclusion, the results provide an overall view of the dynamic changes in physiology and insights into the molecular regulation mechanisms of winter Brassica napus in response to freezing treatment, expanding our understanding on the complex molecular mechanism in plant response to freezing stress.

RNA-Seq analysis reveals transcript diversity and active genes after common cutworm (Spodoptera litura Fabricius) attack in resistant and susceptible wild soybean lines

BACKGROUND

Common cutworm (CCW) is highly responsible for destabilizing soybean productivity. Wild soybean is a resource used by breeders to discover elite defensive genes.

RESULTS

The transcriptomes of two wild accessions (W11 and W99) with different resistance to CCW were analyzed at early- and late-induction time points. After induction, the susceptible accession W11 differentially expressed 1268 and 508 genes at the early and late time points, respectively. Compared with W11, the resistant accession W99 differentially expressed 1270 genes at the early time point and many more genes (2308) at the late time point. In total, 3836 non-redundant genes were identified in both lines. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analyses revealed that the differentially expressed genes (DEGs) in W99 at the late time point were mostly associated with specific processes and pathways. Among the non-redundant genes, 146 genes were commonly up-regulated in the treatment condition compared with the control condition at the early- and late-induction time points in both accessions used in this experiment. Approximately 40% of the common DEGs were related to secondary metabolism, disease resistance, and signal transduction based on their putative function. Excluding the common DEGs, W99 expressed more unique DEGs than W11. Further analysis of the 3836 DEGs revealed that the induction of CCW not only up-regulated defense-related genes, including 37 jasmonic acid (JA)-related genes, 171 plant-pathogen-related genes, and 17 genes encoding protease inhibitors, but also down-regulated growth-related genes, including 35 photosynthesis-related genes, 48 nutrition metabolism genes, and 28 auxin metabolism genes. Therefore, representative defense-related and growth-related genes were chosen for binding site prediction via co-expression of transcription factors (TFs) and spatial expression pattern analyses. In total, 53 binding sites of 28 TFs were identified based on 3 defense-related genes and 3 growth-related genes. Phosphate transporter PT1, which is a representative growth-related gene, was transformed into soybean, and the transgenic soybean plants were susceptible to CCW.

CONCLUSIONS

In summary, we described transcriptome reprograming after herbivore induction in wild soybean, identified the susceptibility of growth-related genes, and provided new resources for the breeding of herbivore-resistant cultivated soybeans.

Transcriptome Analysis of Chinese Chestnut (Castanea mollissima Blume) in Response to Dryocosmus kuriphilus Yasumatsu Infestation

Chinese chestnut (Castanea mollissima Blume) can be infested by Dryocosmus kuriphilus Yasumatsu, resulting in gall formation and yield losses. Research on the control of gall wasps using genomics approaches is rarely reported. We used RNA-seq to investigate the dynamic changes in the genes of a chestnut species (C. mollissima B.) during four gall-formation stages caused by D. kuriphilus. A total of 21,306 genes were annotated by BLAST in databases. Transcriptome comparison between different gall-formation stages revealed many genes that were differentially expressed compared to the control. Among these, 2410, 7373, 6294, and 9412 genes were differentially expressed in four gall-formation stages: initiation stage (A), early growth stage (B), late growth stage (C), and maturation stage (D), respectively. Annotation analysis indicated that many metabolic processes (e.g., phenylpropanoid biosynthesis, secondary metabolism, plant⁻pathogen interaction) were affected. Interesting genes encoding putative components of signal transduction, stress response, and transcription factors were also differentially regulated. These genes might play important roles in response to D. kuriphilus gall formation. These new data on the mechanism by which D. kuriphilus infests chestnuts could help improve chestnut resistance.

Comprehensive RNA sequencing and co-expression network analysis to complete the biosynthetic pathway of coumestrol, a phytoestrogen

Coumestrol (CMS), a coumestan isoflavone, plays key roles in nodulation through communication with rhizobia, and has been used as phytoestrogens for hormone replacement therapy in humans. Because CMS content is controlled by multiple genetic factors, the genetic basis of CMS biosynthesis has remained unclear. We identified soybean genotypes with consistently high (Daewonkong) or low (SS0903-2B-21-1-2) CMS content over 2 years. We performed RNA sequencing of leaf samples from both genotypes at developmental stage R7, when CMS levels are highest. Within the phenylpropanoid biosynthetic pathway, 41 genes were tightly connected in a functional co-expression gene network; seven of these genes were differentially expressed between two genotypes. We identified 14 candidate genes involved in CMS biosynthesis. Among them, seven were annotated as encoding oxidoreductases that may catalyze the transfer of electrons from daidzein, a precursor of CMS. Two of the other genes, annotated as encoding a MYB domain protein and a MLP-like protein, may increase CMS accumulation in response to stress conditions. Our results will help to complete our understanding of the CMS biosynthetic pathway, and should facilitate development of soybean cultivars with high CMS content that could be used to promote the fitness of plants and human beings.

A mongolian pine specific endoplasmic reticulum localized CALMODULIN-LIKE calcium binding protein enhances arabidopsis growth

Stress-adapted wild plants are natural sources of novel genes for molecular breeding. Here, we conducted a transcriptional analysis of Pinus sylvestris var. mongolica Litv, an evergreen pine in northeastern China, to identify a novel CALMODULIN-LIKE protein-encoding gene, PsCML1, no significant homologs found in other plant species. PsCML1 encodes a protein predicted to have a single trans-membrane domain at its N-terminal. Four EF-hand motifs (calcium [Ca]-binding structures) are located at its C-terminal and showed Ca²⁺-specific affinity in isothermal titration calorimetric analysis. Transient expression of PsCML1 in Nicotiana benthamiana showed that the PsCML1 localizes to the endoplasmic reticulum (ER). Heterologous expression of PsCML1 in Arabidopsis significantly promoted seedling growth, and increased resistance to stress from NaCl and Ca²⁺ deficiency. The roots of the transgenic seedlings had higher contents of cellulose and pectin, but less hemicellulose than those of the wild type (WT). The biosynthesis of cell wall components is linked with protein glycosylation in the ER and reactive oxygen species (ROS) homeostasis. No significant difference was found in the extent of protein glycosylation between the transgenic and WT plants. However, the transgenic roots had higher steady-state levels of ROS, NADPH oxidase activity, and endo-membrane dynamics than those of the WT. A working model was proposed to delineate the interaction among Ca²⁺, ROS homeostasis, and cell wall loosening-dependent cell division.