Systems Identification and Characterization of Cell Wall Reassembly and Degradation Related Genes in Glycine max (L.) Merill, a Bioenergy Legume

In Silico Comparison of WRKY Transcription Factors in Wild and Cultivated Soybean and Their Co-expression Network Arbitrating Disease Resistance

WRKY Transcription factors (TFs) play critical roles in plant defence mechanisms that are activated in response to biotic and abiotic stresses. However, information on the Glycine soja WRKYs (GsoWRKYs) is scarce. Owing to its importance in soybean breeding, here we identified putative WRKY TFs in wild soybean, and compared the results with Glycine max WRKYs (GmaWRKYs) by phylogenetic, conserved motif, and duplication analyses. Moreover, we explored the expression trends of WRKYs in G. max (oomycete, fungi, virus, bacteria, and soybean cyst nematode) and G. soja (soybean cyst nematode), and identified commonly expressed WRKYs and their co-expressed genes. We identified, 181 and 180 putative WRKYs in G. max and G. soja, respectively. Though the number of WRKYs in both studied species is almost the same, they differ in many ways, i.e., the number of WRKYs on corresponding chromosomes, conserved domain structures, WRKYGQK motif variants, and zinc-finger motifs. WRKYs in both species grouped in three major clads, i.e., I-III, where group-II had sub-clads IIa-IIe. We found that GsoWRKYs expanded mostly through segmental duplication. A large number of WRKYs were expressed in response to biotic stresses, i.e., Phakospora pachyrhizi, Phytoplasma, Heterodera glycines, Macrophomina phaseolina, and Soybean mosaic virus; 56 GmaWRKYs were commonly expressed in soybean plants infected with these diseases. Finally, 30 and 63 GmaWRKYs and GsoWRKYs co-expressed with 205 and 123 non-WRKY genes, respectively, indicating that WRKYs play essential roles in biotic stress tolerance in Glycine species.

Comparative Transcriptome Sequencing and Endogenous Phytohormone Content of Annual Grafted Branches of Zelkova schneideriana and Its Dwarf Variety HenTianGao

Zelkova schneideriana is a fast-growing tree species endemic to China. Recent surveys and reports have highlighted a continued decline in its natural populations; therefore, it is included in the Red List of Threatened Species by The International Union for Conservation of Nature. A new variety "HenTianGao" (H) has been developed with smaller plant height, slow growth, and lower branching points. In this study, we attempted to understand the differences in plant height of Z. schneideriana (J) and its dwarf variety H. We determined the endogenous hormone content in the annual grafted branches of both J and H. J exhibited higher gibberellic acid (GA)-19 and trans-Zeatin (tZ) levels, whereas H had higher levels of indole-3-acetic acid (IAA) catabolite 2-oxindole-3-acetic acid (OxIAA), IAA-Glu conjugate, and jasmonic acid (JA) (and its conjugate JA-Ile). The transcriptome comparison showed differential regulation of 20,944 genes enriched in growth and development, signaling, and metabolism-related pathways. The results show that the differential phytohormone level (IAA, JA, tZ, and GA) was consistent with the expression of the genes associated with their biosynthesis. The differences in relative OxIAA, IAA-Glu, GA19, trans-Zeatin, JA, and JA-Ile levels were linked to changes in respective signaling-related genes. We also observed significant differences in the expression of cell size, number, proliferation, cell wall biosynthesis, and remodeling-related genes in J and H. The differences in relative endogenous hormone levels, expression of biosynthesis, and signaling genes provide a theoretical basis for understanding the plant height differences in Z. schneideriana.

Integrated transcriptome and metabolome analysis unveil the response mechanism in wild rice (Zizania latifolia griseb.) against sheath rot infection

Sheath rot disease (SRD) is one of the most devastating diseases of Manchurian wild rice (MWR) (Zizania latifolia Griseb). Pilot experiments in our laboratory have shown that an MWR cultivar "Zhejiao NO.7"exhibits signs of SRD tolerance. To explore the responses of Zhejiao No. 7 to SRD infection, we used a combined transcriptome and metabolome analysis approach. A total of 136 differentially accumulated metabolites (DAMs, 114 up- and 22 down-accumulated in FA compared to CK) were detected. These up-accumulated metabolites were enriched in tryptophan metabolism, amino acid biosynthesis, flavonoids, and phytohormone signaling. Transcriptome sequencing results showed the differential expression of 11,280 genes (DEGs, 5,933 up-, and 5,347 downregulated in FA compared to CK). The genes expressed in tryptophan metabolism, amino acid biosynthesis, phytohormone biosynthesis and signaling, and reactive oxygen species homeostasis confirmed the metabolite results. In addition, genes related to the cell wall, carbohydrate metabolism, and plant-pathogen interaction (especially hypersensitive response) showed changes in expression in response to SRD infection. These results provide a basis for understanding the response mechanisms in MWR to FA attack that can be used for breeding SRD-tolerant MWR.

Twenty years of mining salt tolerance genes in soybean

Current combined challenges of rising food demand, climate change and farmland degradation exert enormous pressure on agricultural production. Worldwide soil salinization, in particular, necessitates the development of salt-tolerant crops. Soybean, being a globally important produce, has its genetic resources increasingly examined to facilitate crop improvement based on functional genomics. In response to the multifaceted physiological challenge that salt stress imposes, soybean has evolved an array of defences against salinity. These include maintaining cell homeostasis by ion transportation, osmoregulation, and restoring oxidative balance. Other adaptations include cell wall alterations, transcriptomic reprogramming, and efficient signal transduction for detecting and responding to salt stress. Here, we reviewed functionally verified genes that underly different salt tolerance mechanisms employed by soybean in the past two decades, and discussed the strategy in selecting salt tolerance genes for crop improvement. Future studies could adopt an integrated multi-omic approach in characterizing soybean salt tolerance adaptations and put our existing knowledge into practice via omic-assisted breeding and gene editing. This review serves as a guide and inspiration for crop developers in enhancing soybean tolerance against abiotic stresses, thereby fulfilling the role of science in solving real-life problems.

Supplementary Information

The online version contains supplementary material available at 10.1007/s11032-023-01383-3.

Transcriptome analysis reveals the key pathways and candidate genes involved in salt stress responses in Cymbidium ensifolium leaves

BACKGROUND

Cymbidium ensifolium L. is known for its ornamental value and is frequently used in cosmetics. Information about the salt stress response of C. ensifolium is scarce. In this study, we reported the physiological and transcriptomic responses of C. ensifolium leaves under the influence of 100 mM NaCl stress for 48 (T48) and 96 (T96) hours.

RESULTS

Leaf Na⁺ content, activities of the antioxidant enzymes i.e., superoxide dismutase, glutathione S-transferase, and ascorbate peroxidase, and malondialdehyde content were increased in salt-stressed leaves of C. ensifolium. Transcriptome analysis revealed that a relatively high number of genes were differentially expressed in CKvsT48 (17,249) compared to CKvsT96 (5,376). Several genes related to salt stress sensing (calcium signaling, stomata closure, cell-wall remodeling, and ROS scavenging), ion balance (Na⁺ and H⁺), ion homeostasis (Na⁺/K⁺ ratios), and phytohormone signaling (abscisic acid and brassinosteroid) were differentially expressed in CKvsT48, CKvsT96, and T48vsT96. In general, the expression of genes enriched in these pathways was increased in T48 compared to CK while reduced in T96 compared to T48. Transcription factors (TFs) belonging to more than 70 families were differentially expressed; the major families of differentially expressed TFs included bHLH, NAC, MYB, WRKY, MYB-related, and C3H. A Myb-like gene (CenREV3) was further characterized by overexpressing it in Arabidopsis thaliana. CenREV3's expression was decreased with the prolongation of salt stress. As a result, the CenREV3-overexpression lines showed reduced root length, germination %, and survival % suggesting that this TF is a negative regulator of salt stress tolerance.

CONCLUSION

These results provide the basis for future studies to explore the salt stress response-related pathways in C. ensifolium.

Overexpression of the GmEXPA1 gene reduces plant susceptibility to Meloidogyne incognita

The overexpression of the soybean GmEXPA1 gene reduces plant susceptibility to M. incognita by the increase of root lignification. Plant expansins are enzymes that act in a pH-dependent manner in the plant cell wall loosening and are associated with improved tolerance or resistance to abiotic or biotic stresses. Plant-parasitic nematodes (PPN) can alter the expression profile of several expansin genes in infected root cells. Studies have shown that overexpression or downregulation of particular expansin genes can reduce plant susceptibility to PPNs. Root-knot nematodes (RKN) are obligate sedentary endoparasites of the genus Meloidogyne spp. of which M. incognita is one of the most reported species. Herein, using a transcriptome dataset and real-time PCR assays were identified an expansin A gene (GmEXPA1; Glyma.02G109100) that is upregulated in the soybean nematode-resistant genotype PI595099 compared to the susceptible cultivar BRS133 during plant parasitism by M. incognita. To understand the role of the GmEXPA1 gene during the interaction between soybean plant and M. incognita were generated stable A. thaliana and N. tabacum transgenic lines. Remarkably, both A. thaliana and N. tabacum transgenic lines overexpressing the GmEXPA1 gene showed reduced susceptibility to M. incognita. Furthermore, plant growth, biomass accumulation, and seed yield were not affected in these transgenic lines. Interestingly, significant upregulation of the NtACC oxidase and NtEFE26 genes, involved in ethylene biosynthesis, and NtCCR and Nt4CL genes, involved in lignin biosynthesis, was observed in roots of the N. tabacum transgenic lines, which also showed higher lignin content. These data suggested a possible link between GmEXPA1 gene expression and increased lignification of the root cell wall. Therefore, these data support that engineering of the GmEXPA1 gene in soybean offers a powerful biotechnology tool to assist in RKN management.

Genome-Wide Identification and Expression Analysis of the Xyloglucan Endotransglucosylase/Hydrolase Gene Family in Sweet Potato [Ipomoea batatas (L.) Lam]

The xyloglucan endotransglucosylase/hydrolase (XET/XEH, also named XTH) family is a multigene family, the function of which plays a significant role in cell-wall rebuilding and stress tolerance in plants. However, the specific traits of the XTH gene family members and their expression pattern in different tissues and under stress have not been carried out in sweet potato. Thirty-six XTH genes were identified in I. batatas, all of which had conserved structures (Glyco_hydro_16). Based on Neighbor-Joining phylogenetic analysis the IbXTHs can be divided into three subfamilies-the I/II, IIIA, and IIIB subfamilies, which were unevenly distributed on 13 chromosomes, with the exception of Chr9 and Chr15. Multiple cis-acting regions related to growth and development, as well as stress responses, may be found in the IbXTH gene promoters. The segmental duplication occurrences greatly aided the evolution of IbXTHs. The results of a collinearity analysis showed that the XTH genes of sweet potato shared evolutionary history with three additional species, including A. thaliana, G. max, and O. sativa. Additionally, based on the transcriptome sequencing data, the results revealed that the IbXTHs have different expression patterns in leaves, stems, the root body (RB), the distal end (DE), the root stock (RS), the proximal end (PE), the initiative storage root (ISR), and the fibrous root (FR), and many of them are well expressed in the roots. Differentially expressed gene (DEG) analysis of FRs after hormone treatment of the roots indicated that IbXTH28 and IbXTH30 are up-regulated under salicylic acid (SA) treatment but down-regulated under methyl jasmonate (MeJA) treatment. Attentionally, there were only two genes showing down-regulation under the cold and drought treatment. Collectively, all of the findings suggested that genes from the XTH family are crucial for root specificity. This study could provide a theoretical basis for further research on the molecular function of sweet potato XTH genes.

Identification of the Xyloglucan Endotransglycosylase/Hydrolase (XTH) Gene Family Members Expressed in Boehmeria nivea in Response to Cadmium Stress

Xyloglucan endotransglycosylase/hydrolase (XTH) genes play an important role in plant resistance to abiotic stress. However, systematic studies of the response of Boehmeria nivea (ramie) XTH genes (BnXTHs) to cadmium (Cd) stress are lacking. We sought to identify the BnXTH-family genes in ramie through bioinformatics analyses and to investigate their responses to Cd stress. We identified 19 members of the BnXTH gene family from the ramie genome, referred to as BnXTH1-19, among which BnXTH18 and BnXTH19 were located on no chromosomes and the remaining genes were unevenly distributed across 11 chromosomes. The 19 members were divided into four groups, Groups I/II/IIIA/IIIB, according to their phylogenetic relationships, and these groups were supported by analyses of intron-exon structure and conserved motif composition. A highly conserved catalytic site (HDEIDFEFLG) was observed in all BnXTH proteins. Additionally, three gene pairs (BnXTH6-BnXTH16, BnXTH8-BnXTH9, and BnXTH17-BnXTH18) were obtained with a fragment and tandem-repeat event analysis of the ramie genome. An analysis of cisregulatory elements revealed that BnXTH expression might be regulated by multiple hormones and abiotic and biotic stress responses. In particular, 17 cisregulatory elements related to abiotic and biotic stress responses and 11 cisregulatory elements related to hormone responses were identified. We also found that most BnXTH genes responded to Cd stress, and BnXTH1, BnXTH3, BnXTH6, and BnXTH15 were most likely to contribute to the Cd tolerance of ramie, as evidenced by the substantial increases in expression under Cd treatment. Heterologous expression of BnXTH1, BnXTH6, and BnXTH15 significantly enhanced the Cd tolerance of transgenic yeast cells. These results suggest that the BnXTH gene family is involved in Cd stress responses, laying a theoretical foundation for functional studies of BnXTH genes and the innovative breeding of Cd-tolerant ramie.

Leaf-transcriptome profiles of phoebe bournei provide insights into temporal drought stress responses

Phoebe bournei (Hemsl.) Yang is used as a commercial wood in China and is enlisted as a near-threatened species. Prolonged droughts pose a serious threat to young seedlings (1-2 years old). A transcriptome sequencing approach, together with the measurement of growth parameters and biochemical analyses were used to understand P. bournei's drought responses on 15d, 30d, and 45d of drought stress treatment. The stem and root dry weights decreased significantly with drought stress duration. Activities of antioxidative enzymes i.e., peroxidase (POD), superoxide dismutase (SOD), and catalase (CAT) increased significantly with the increase in drought stress duration. A total of 13,274, 15,648, and 9,949 genes were differentially expressed in CKvs15d, CKvs30d, and CKvs45d, respectively. The differential expression analyses showed that photosystem I and II underwent structural changes, chlorophyll biosynthesis, and photosynthesis were reduced. The genes annotated as POD, SOD, and CAT were upregulated in drought-treated leaves as compared to control. Additionally, plant-hormone signal transduction, MAPK signaling-plant, phenylpropanoid biosynthesis, flavonoid biosynthesis, and starch and sucrose metabolism pathways showed large-scale expression changes in major genes. We also found that members of 25 transcription factor families were differentially expressed. Our study presents and discusses these transcriptome signatures. Overall, our findings represent key data for breeding towards drought stress tolerance in P. bournei.

Physiological and Transcriptomic Responses of Growth in Neolamarckia cadamba Stimulated by Exogenous Gibberellins

(1) The phytohormones gibberellins (GAs) play a crucial role in plant growth and development, such as seed germination, flowering, fruiting, and stem elongation. Although many biological roles of GAs have been studied intensively, the molecular mechanisms of GAs in woody plants are still unclear. (2) In this study, we investigated the effects of exogenous application of GAs on Neolamarckia cadamba. (3) The height and biomass of N. cadamba increased after 7 days of GA treatment, especially on the second internode. Transcriptome analysis showed that although the majority of genes involved in the GA signaling pathway were up-regulated, the expression of GA20 oxidase (GA20ox) and GA3 oxidase (GA3ox) was down-regulated in the 3 days GA-treated group compared to the CK group. The expression of the cell elongation-related basic helix-loop-helix genes bHLH74 and bHLH49 was up-regulated in the GA-treated group compared with the CK group. Transcriptional expression levels of transcription factors involved in hormone signaling were changed, mainly including bHLH, ethylene response factor (ERF), and WRKY families. In addition, the transcriptional expression level of the key enzymes engaged in the phenylalanine pathway was downregulated after GA treatment. (4) In brief, our findings reveal the physiological and molecular mechanisms of exogenous GA treatment stimulation in N. cadamba.

Combined transcriptome and metabolome analysis of Nerium indicum L. elaborates the key pathways that are activated in response to witches' broom disease

BACKGROUND

Nerium indicum Mill. is an ornamental plant that is found in parks, riversides, lakesides, and scenic areas in China and other parts of the world. Our recent survey indicated the prevalence of witches' broom disease (WBD) in Guangdong, China. To find out the possible defense strategies against WBD, we performed a MiSeq based ITS sequencing to identify the possible casual organism, then did a de novo transcriptome sequencing and metabolome profiling in the phloem and stem tip of N. indicum plants suffering from WBD compared to healthy ones.

RESULTS

The survey showed that Wengyuen county and Zengcheng district had the highest disease incidence rates. The most prevalent microbial species in the diseased tissues was Cophinforma mamane. The transcriptome sequencing resulted in the identification of 191,224 unigenes of which 142,396 could be annotated. There were 19,031 and 13,284 differentially expressed genes (DEGs) between diseased phloem (NOWP) and healthy phloem (NOHP), and diseased stem (NOWS) and healthy stem (NOHS), respectively. The DEGs were enriched in MAPK-signaling (plant), plant-pathogen interaction, plant-hormone signal transduction, phenylpropanoid and flavonoid biosynthesis, linoleic acid and α-linoleic acid metabolism pathways. Particularly, we found that N. indicum plants activated the phytohormone signaling, MAPK-signaling cascade, defense related proteins, and the biosynthesis of phenylpropanoids and flavonoids as defense responses to the pathogenic infection. The metabolome profiling identified 586 metabolites of which 386 and 324 metabolites were differentially accumulated in NOHP vs NOWP and NOHS and NOWS, respectively. The differential accumulation of metabolites related to phytohormone signaling, linoleic acid metabolism, phenylpropanoid and flavonoid biosynthesis, nicotinate and nicotinamide metabolism, and citrate cycle was observed, indicating the role of these pathways in defense responses against the pathogenic infection.

CONCLUSION

Our results showed that Guangdong province has a high incidence of WBD in most of the surveyed areas. C. mamane is suspected to be the causing pathogen of WBD in N. indicum. N. indicum initiated the MAPK-signaling cascade and phytohormone signaling, leading to the activation of pathogen-associated molecular patterns and hypersensitive response. Furthermore, N. indicum accumulated high concentrations of phenolic acids, coumarins and lignans, and flavonoids under WBD. These results provide scientific tools for the formulation of control strategies of WBD in N. indicum.

Investigating the Mechanism of Unilateral Cross Incompatibility in Longan (Dimocarpus longan Lour.) Cultivars (Yiduo × Shixia)

Longan (Dimocarpus longan Lour.) is an important subtropical fruit tree in China. Nearly 90% of longan fruit imports from Thailand are from the cultivar Yiduo. However, we have observed that there exists a unilateral cross incompatibility (UCI) when Yiduo is used as a female parent and Shixia (a famous Chinese cultivar) as a male parent. Here, we performed a comparative transcriptome analysis coupled with microscopy of pistils from two reciprocal pollination combinations [Shixia♂ × Yiduo♀(SY) and Yiduo♀ × Shixia♂(YS)] 4, 8, 12, and 24 h after pollination. We also explored endogenous jasmonic acid (JA) and jasmonyl isoleucine (JA-Ile) levels in pistils of the crosses. The microscopic observations showed that the UCI was sporophytic. The endogenous JA and JA-Ile levels were higher in YS than in SY at the studied time points. We found 7,251 differentially expressed genes from the transcriptome analysis. Our results highlighted that genes associated with JA biosynthesis and signaling, pollen tube growth, cell wall modification, starch and sucrose biosynthesis, and protein processing in endoplasmic reticulum pathways were differentially regulated between SY and YS. We discussed transcriptomic changes in the above-mentioned pathways regarding the observed microscopic and/or endogenous hormone levels. This is the first report on the elaboration of transcriptomic changes in longan reciprocal pollination combination showing UCI. The results presented here will enable the longan breeding community to better understand the mechanisms of UCI.

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

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

Transcriptomic responses to drought stress in Polygonatum kingianum tuber

BACKGROUND

Polygonatum kingianum Coll. et Hemsl. is an important plant in Traditional Chinese Medicine. The extracts from its tubers are rich in polysaccharides and other metabolites such as saponins. It is a well-known concept that growing medicinal plants in semi-arid (or drought stress) increases their natural compounds concentrations. This study was conducted to explore the morpho-physiological responses of P. kingianum plants and transcriptomic signatures of P. kingianum tubers exposed to mild, moderate, and severe drought and rewatering.

RESULTS

The stress effects on the morpho-physiological parameters were dependent on the intensity of the drought stress. The leaf area, relative water content, chlorophyll content, and shoot fresh weight decreased whereas electrolyte leakage increased with increase in drought stress intensity. A total of 53,081 unigenes were obtained; 59% of which were annotated. We observed that 1352 and 350 core genes were differentially expressed in drought and rewatering, respectively. Drought stress driven differentially expressed genes (DEGs) were enriched in phenylpropanoid biosynthesis, flavonoid biosynthesis, starch and sucrose metabolism, and stilbenoid diarylheptanoid and gingerol biosynthesis, and carotenoid biosynthesis pathways. Pathways such as plant-pathogen interaction and galactose metabolism were differentially regulated between severe drought and rewatering. Drought reduced the expression of lignin, gingerol, and flavonoid biosynthesis related genes and rewatering recovered the tubers from stress by increasing the expression of the genes. Increased expression of carotenoid biosynthesis pathway related genes under drought suggested their important role in stress endurance. An increase in starch and sucrose biosynthesis was evident from transcriptomic changes under drought stress. Rewatering recovered the drought affected tubers as evident from the contrasting expression profiles of genes related to these pathways. P. kingianum tuber experiences an increased biosynthesis of sucrose, starch, and carotenoid under drought stress. Drought decreases the flavonoids, phenylpropanoids, gingerol, and lignin biosynthesis. These changes can be reversed by rewatering the P. kingianum plants.

CONCLUSIONS

These results provide a transcriptome resource for P. kingianum and expands the knowledge on the effect of drought and rewatering on important pathways. This study also provides a large number of candidate genes that could be manipulated for drought stress tolerance and managing the polysaccharide and secondary metabolites' contents in P. kingianum.

Transcriptome and Metabolome Comparison of Smooth and Rough Citrus limon L. Peels Grown on Same Trees and Harvested in Different Seasons

Background: Farmers harvest two batches fruits of Lemons (Citrus limon L. Burm. f.) i.e., spring flowering fruit and autumn flowering fruit in dry-hot valley in Yunnan, China. Regular lemons harvested in autumn have smooth skin. However, lemons harvested in spring have rough skin, which makes them less attractive to customers. Furthermore, the rough skin causes a reduction in commodity value and economical losses to farmers. This is a preliminary study that investigates the key transcriptomic and metabolomic differences in peels of lemon fruits (variety Yuning no. 1) harvested 30, 60, 90, 120, and 150 days after flowering from the same trees in different seasons. Results: We identified 5,792, 4,001, 3,148, and 5,287 differentially expressed genes (DEGs) between smooth peel (C) and rough peel (D) 60, 90, 120, and 150 days after flowering, respectively. A total of 1,193 metabolites differentially accumulated (DAM) between D and C. The DEGs and DAMs were enriched in the mitogen-activated protein kinase (MAPK) and plant hormone signaling, terpenoid biosynthesis, flavonoid, and phenylalanine biosynthesis, and ribosome pathways. Predominantly, in the early stages, phytohormonal regulation and signaling were the main driving force for changes in peel surface. Changes in the expression of genes associated with asymmetric cell division were also an important observation. The biosynthesis of terpenoids was possibly reduced in rough peels, while the exclusive expression of cell wall synthesis-related genes could be a possible reason for the thick peel of the rough-skinned lemons. Additionally, cell division, cell number, hypocotyl growth, accumulation of fatty acids, lignans and coumarins- related gene expression, and metabolite accumulation changes were major observations. Conclusion: The rough peels fruit (autumn flowering fruit) and smooth peels fruit (spring flowering fruit) matured on the same trees are possibly due to the differential regulation of asymmetric cell division, cell number regulation, and randomization of hypocotyl growth related genes and the accumulation of terpenoids, flavonoids, fatty acids, lignans, and coumarins. The preliminary results of this study are important for increasing the understanding of peel roughness in lemon and other citrus species.

An insight into transcriptome of Cyathus bulleri for lignocellulase expression on wheat bran

To identify enzymes that can be effectively used for hydrolysis of lignocellulosic biomass, an attractive carbon source in biorefineries, transcriptome analysis was carried out of wheat bran grown fungus, Cyathus bulleri. A comprehensive set of transcripts, encoding carbohydrate active enzymes, were identified. These belonged to 55, 32, 12, 11 and 7 different families of the enzyme classes of Glycoside Hydrolases (GHs), Glycosyl Transferases (GTs), Auxiliary Activities (AAs), Carbohydrate Esterases (CEs) and Polysaccharide Lyases (PLs) respectively. Higher levels of transcripts were obtained for proteins encoding cellulose and hemicellulose degrading activities (of the GH class) with the highest diversity found in the transcripts encoding the hemicellulases. Several transcripts encoding pectin degrading activity were also identified indicating close association of the pectin with the cellulose/hemicellulose in the cell wall of this fungus. Transcripts encoding ligninases were categorized into Cu radical oxidase, Glucose-Methanol-Choline oxidoreductase (with 37 different transcripts in the AA3 sub-family), Laccase and Manganese peroxidases. Temporal gene expression profile for laccase isoforms was studied to understand their role in lignin degradation. To our knowledge, this is the first analysis of the transcriptome of a member belonging to the family Nidulariaceae.

Transcriptome and metabolome profiling provide insights into molecular mechanism of pseudostem elongation in banana

BACKGROUND

Banana plant height is an important trait for horticultural practices and semi-dwarf cultivars show better resistance to damages by wind and rain. However, the molecular mechanisms controlling the pseudostem height remain poorly understood. Herein, we studied the molecular changes in the pseudostem of a semi-dwarf banana mutant Aifen No. 1 (Musa spp. Pisang Awak sub-group ABB) as compared to its wild-type dwarf cultivar using a combined transcriptome and metabolome approach.

RESULTS

A total of 127 differentially expressed genes and 48 differentially accumulated metabolites were detected between the mutant and its wild type. Metabolites belonging to amino acid and its derivatives, flavonoids, lignans, coumarins, organic acids, and phenolic acids were up-regulated in the mutant. The transcriptome analysis showed the differential regulation of genes related to the gibberellin pathway, auxin transport, cell elongation, and cell wall modification. Based on the regulation of gibberellin and associated pathway-related genes, we discussed the involvement of gibberellins in pseudostem elongation in the mutant banana. Genes and metabolites associated with cell wall were explored and their involvement in cell extension is discussed.

CONCLUSIONS

The results suggest that gibberellins and associated pathways are possibly developing the observed semi-dwarf pseudostem phenotype together with cell elongation and cell wall modification. The findings increase the understanding of the mechanisms underlying banana stem height and provide new clues for further dissection of specific gene functions.

Transcriptional response of Asarum heterotropoides Fr. Schmidt var. mandshuricum (Maxim.) Kitag. leaves grown under full and partial daylight conditions

BACKGROUND

Asarum heterotropides Fr. Schmidt var. mandshuricum (Maxim.) Kitag. is an important medicinal and industrial plant, which is used in the treatment of various diseases. The main bioactive ingredient is the volatile oil having more than 82 identified components of which methyleugenol, safrole, myristicin, and toluene account for about 70% of the total volume. As a sciophyte plant, the amount of light it absorbs through leaves is an important factor for growth and metabolism.

RESULTS

We grew Asarum plants under full, 50, 28, and 12% sunlight conditions to investigate the effect of different light irradiances on the four major volatile oil components. We employed de novo transcriptome sequencing to understand the transcriptional behavior of Asarum leaves regarding the biosynthetic pathways of the four volatile oil components, photosynthesis and biomass accumulation, and hormone signaling. Our results demonstrated that the increasing light conditions promoted higher percent of the four components. Under full sunlight conditions, cinnamyl alcohol dehydrogenase and cytochrome p450719As were upregulated and led the increased methyleugenol, safrole, and myristicin. The transcriptomic data also showed that Asarum leaves, under full sunlight conditions, adjust their photosynthesis-antenna proteins as a photoprotective response with the help of carotenoids. Plant hormone-signaling related genes were also differentially expressed between full sunlight and low light conditions.

CONCLUSIONS

High light induces accumulation of major bioactive ingredients A. heterotropides volatile oil and this is ascribed to upregulation of key genes such as cinnamyl alcohol dehydrogenase and cytochrome p450719As. The transcriptome data presented here lays the foundation of further understanding of light responses in sciophytes and provides guidance for increasing bioactive molecules in Asarum.

Genome-wide identification, and phylogenetic and expression profiling analyses, of XTH gene families in Brassica rapa L. and Brassica oleracea L

BACKGROUND

Xyloglucan endotransglucosylase/hydrolase genes (XTHs) are a multigene family and play key roles in regulating cell wall extensibility in plant growth and development. Brassica rapa and Brassica oleracea contain XTHs, but detailed identification and characterization of the XTH family in these species, and analysis of their tissue expression profiles, have not previously been carried out.

RESULTS

In this study, 53 and 38 XTH genes were identified in B. rapa and B. oleracea respectively, which contained some novel members not observed in previous studies. All XTHs of B. rapa, B. oleracea and Arabidopsis thaliana could be classified into three groups, Group I/II, III and the Early diverging group, based on phylogenetic relationships. Gene structures and motif patterns were similar within each group. All XTHs in this study contained two characteristic conserved domains (Glyco_hydro and XET_C). XTHs are located mainly in the cell wall but some are also located in the cytoplasm. Analyses of the mechanisms of gene family expansion revealed that whole-genome triplication (WGT) events and tandem duplication (TD) may have been the major mechanisms accounting for the expansion of the XTH gene family. Interestingly, TD genes all belonged to Group I/II, suggesting that TD was the main reason for the largest number of genes being in these groups. B. oleracea had lost more of the XTH genes, the conserved domain XET_C and the conserved active-site motif EXDXE compared with B. rapa, consistent with asymmetrical evolution between the two Brassica genomes. A majority of XTH genes exhibited different tissue-specific expression patterns based on RNA-seq data analyses. Moreover, there was differential expression of duplicated XTH genes in the two species, indicating that their functional differentiation occurred after B. rapa and B. oleracea diverged from a common ancestor.

CONCLUSIONS

We carried out the first systematic analysis of XTH gene families in B. rapa and B. oleracea. The results of this investigation can be used for reference in further studies on the functions of XTH genes and the evolution of this multigene family.

Overexpression of an ethylene-forming ACC oxidase (ACO) gene precedes the Minute Hilum seed coat phenotype in Glycine max

BACKGROUND

To elucidate features of seed development, we investigated the transcriptome of a soybean isoline from the germplasm collection that contained an introgressed allele known as minute hilum (mi) which confers a smaller hilum region where the seed attaches to the pod and also results in seed coat cracking surrounding the hilum region.

RESULTS

RNAs were extracted from immature seed from an extended hilum region (i.e., the hilum and a small ring of tissue surrounding the hilum in which the cracks form) at three different developmental stages:10-25, 25-50 and 50-100 mg seed fresh weight in two independent replicates for each stage. The transcriptomes of these samples from both the Clark isoline containing the mi allele (PI 547628, UC413, ii R t mi G), and its recurrent Clark 63 parent isoline (PI 548532, UC7, ii R T Mi g), which was used for six generations of backcrossing, were compared for differential expression of 88,648 Glyma models of the soybean genome Wm82.a2. The RNA sequence data obtained from the 12 cDNA libraries were subjected to padj value < 0.05 and at least two-fold expression differences to select with confidence genes differentially expressed in the hilum-containing tissue of the seed coat between the two lines. Glyma.09G008400 annotated as encoding an ethylene forming enzyme, ACC oxidase (ACO), was found to be highly overexpressed in the mi hilum region at 165 RPKMs (reads per kilobase per million mapped reads) compared to the standard line at just 0.03 RPKMs. Evidence of changes in expression of genes downstream of the ethylene pathway included those involved in auxin and gibberellin hormone action and extensive differences in expression of cell wall protein genes. These changes are postulated to determine the restricted hilum size and cracking phenotypes.

CONCLUSIONS

We present transcriptome and phenotypic evidence that substantially higher expression of an ethylene-forming ACO gene likely shifts hormone balance and sets in motion downstream changes resulting in a smaller hilum phenotype and the cracks observed in the minute hilum (mi) isoline as compared to its recurrent parent.

Transcriptomic and Metabolomic Changes Triggered by Fusarium solani in Common Bean (Phaseolus vulgaris L.)

Common bean (Phaseolus vulgaris L.) is a major legume and is frequently attacked by fungal pathogens, including Fusarium solani f. sp. phaseoli (FSP), which cause Fusarium root rot. FSP substantially reduces common bean yields across the world, including China, but little is known about how common bean plants defend themselves against this fungal pathogen. In the current study, we combined next-generation RNA sequencing and metabolomics techniques to investigate the changes in gene expression and metabolomic processes in common bean infected with FSP. There were 29,722 differentially regulated genes and 300 differentially regulated metabolites between control and infected plants. The combined omics approach revealed that FSP is perceived by PAMP-triggered immunity and effector-triggered immunity. Infected seedlings showed that common bean responded by cell wall modification, ROS generation, and a synergistic hormone-driven defense response. Further analysis showed that FSP induced energy metabolism, nitrogen mobilization, accumulation of sugars, and arginine and proline metabolism. Importantly, metabolic pathways were most significantly enriched, which resulted in increased levels of metabolites that were involved in the plant defense response. A correspondence between the transcript pattern and metabolite profile was observed in the discussed pathways. The combined omics approach enhances our understanding of the less explored pathosystem and will provide clues for the development of common bean cultivars' resistant to FSP.

Combined De Novo Transcriptome and Metabolome Analysis of Common Bean Response to Fusarium oxysporum f. sp. phaseoli Infection

Molecular changes elicited by common bean (Phaseolus vulgaris L.) in response to Fusarium oxysproum f. sp. Phaseoli (FOP) remain elusive. We studied the changes in root metabolism during common bean–FOP interactions using a combined de novo transcriptome and metabolome approach. Our results demonstrated alterations of transcript levels and metabolite concentrations in common bean roots 24 h post infection as compared to control. The transcriptome and metabolome responses in common bean roots revealed significant changes in structural defense i.e., cell-wall loosening and weakening characterized by hyper accumulation of cell-wall loosening and degradation related transcripts. The levels of pathogenesis related genes were significantly higher upon FOP inoculation. Interestingly, we found the involvement of glycosylphosphatidylinositol- anchored proteins (GPI-APs) in signal transduction in response to FOP infection. Our results confirmed that hormones have strong role in signaling pathways i.e., salicylic acid, jasmonate, and ethylene pathways. FOP induced energy metabolism and nitrogen mobilization in infected common bean roots as compared to control. Importantly, the flavonoid biosynthesis pathway was the most significantly enriched pathway in response to FOP infection as revealed by the combined transcriptome and metabolome analysis. Overall, the observed modulations in the transcriptome and metabolome flux as outcome of several orchestrated molecular events are determinant of host’s role in common bean–FOP interactions.

Phytolith Formation in Plants: From Soil to Cell

Silica is deposited extra- and intracellularly in plants in solid form, as phytoliths. Phytoliths have emerged as accepted taxonomic tools and proxies for reconstructing ancient flora, agricultural economies, environment, and climate. The discovery of silicon transporter genes has aided in the understanding of the mechanism of silicon transport and deposition within the plant body and reconstructing plant phylogeny that is based on the ability of plants to accumulate silica. However, a precise understanding of the process of silica deposition and the formation of phytoliths is still an enigma and the information regarding the proteins that are involved in plant biosilicification is still scarce. With the observation of various shapes and morphologies of phytoliths, it is essential to understand which factors control this mechanism. During the last two decades, significant research has been done in this regard and silicon research has expanded as an Earth-life science superdiscipline. We review and integrate the recent knowledge and concepts on the uptake and transport of silica and its deposition as phytoliths in plants. We also discuss how different factors define the shape, size, and chemistry of the phytoliths and how biosilicification evolved in plants. The role of channel-type and efflux silicon transporters, proline-rich proteins, and siliplant1 protein in transport and deposition of silica is presented. The role of phytoliths against biotic and abiotic stress, as mechanical barriers, and their use as taxonomic tools and proxies, is highlighted.

Neonicotinoid Insecticides Alter the Transcriptome of Soybean and Decrease Plant Resistance

Neonicotinoids are widely used systemic insecticides that have been associated with spider mite outbreaks on diverse plants. These insecticides have complex effects on plant physiology, which have been speculated to drive enhanced performance of spider mites. We used RNA-Seq to explore how neonicotinoids modify gene expression in soybean thereby lowering plant resistance. We exposed soybean (Glycine max L.) to two neonicotinoid insecticides, thiamethoxam applied to seeds and imidacloprid applied as a soil drench, and we exposed a subset of these plants to spider mites (Tetranychus cinnabarinus). Applications of both insecticides downregulated genes involved in plant-pathogen interactions, phytohormone pathways, phenylpropanoid pathway, and cell wall biosynthesis. These effects were especially pronounced in plants exposed to thiamethoxam. Introduction of spider mites restored induction of genes in these pathways in plants treated with imidacloprid, while expression of genes involved in phenylpropanoid synthesis, in particular, remained downregulated in thiamethoxam-treated plants. Our outcomes indicate that both insecticides suppress genes in pathways relevant to plant⁻arthropod interactions, and suppression of genes involved in cell wall synthesis may explain lower plant resistance to spider mites, cell-content feeders. These effects appear to be particularly significant when plants are exposed to neonicotinoids applied to soybean seeds.

Characterization of Cellulose Synthase A (CESA) Gene Family in Eudicots

Cellulose synthase A (CESA) is a key enzyme involved in the complex process of plant cell wall biosynthesis, and it remains a productive subject for research. We employed systems biology approaches to explore structural diversity of eudicot CESAs by exon-intron organization, mode of duplication, synteny, and splice site analyses. Using a combined phylogenetics and comparative genomics approach coupled with co-expression networks we reconciled the evolution of cellulose synthase gene family in eudicots and found that the basic forms of CESA proteins are retained in angiosperms. Duplications have played an important role in expansion of CESA gene family members in eudicots. Co-expression networks showed that primary and secondary cell wall modules are duplicated in eudicots. We also identified 230 simple sequence repeat markers in 103 eudicot CESAs. The 13 identified conserved motifs in eudicots will provide a basis for gene identification and functional characterization in other plants. Furthermore, we characterized (in silico) eudicot CESAs against senescence and found that expression levels of CESAs decreased during leaf senescence.