Bioinformatics analysis and function prediction of NBS-LRR gene family in Broussonetia papyrifera.

Most of the currently available disease resistance (R) genes have NBS (nucleotide-binding site) and LRR (leucine-rich-repeat) domain which belongs to the NBS-LRR gene family. The whole genome sequencing of Broussonetia papyrifera provides an important bioinformatics database for the study of the NBS-LRR gene family. In this study, 328 NBS-LRR family genes were identified and classified in B. papyrifera according to different classification schemes, where there are 92 N types, 47 CN type, 54 CNL type, 29 NL types, 55 TN type, and 51 TNL type. Subsequently, we conducted bioinformatics analysis of the NBS-LRR gene family. Classification, motif analysis of protein sequences, and phylogenetic tree studies of the NBS-LRR genes in B. papyrifera provide important basis for the functional study of NBS-LRR family genes. Additionally, we performed structural analysis of the chromosomal location, physicochemical properties, and sequences identified by genetic characterization. In addition, through the analysis of GO enrichment, it was found that NBS-LRR genes were involved in defense responses and were significantly enriched in biological stimulation, immune response, and abiotic stress. In addition, we found that Bp06g0955 was the most sensitive to low temperature and encoded the RPM1 protein by analyzing the low temperature transcriptome data of B. papyrifera. Quantitative results of gene expression after 48 h of Fusarium infection showed that Bp01g3293 increased 14 times after infection, which encodes RPM1 protein. The potential of NBS-LRR gene responsive to biotic and abiotic stresses can be exploited to improve the resistance of B. papyrifera.

Practical method for dynamic color holographic display.

A practical method for dynamic color holographic display by using a computer-generated hologram (CGH) with a high space-bandwidth product is proposed, and a dynamic color holographic display system is designed by a space-division method. First, three primary color CGHs of different frames from a color movie are fabricated on holographic recording material by a self-made CGH microfilming system. Secondly, the CGH is fixed on an X-Y moving stage, which is controlled by the system in order to bring the CGH to the appointed position. Thirdly, three primary color lasers are used to reconstruct the CGH. The switch of the lasers is controlled by the system synchronous with the X-Y moving stage. The color video with high quality can be obtained after filtering the three primary color reconstructed wavefronts. The experimental results demonstrate that the proposed dynamic color holographic display method is effective. It has practical application value in high-quality CGH display.

Identification and Investigation of the Genetic Variations and Candidate Genes Responsible for Seed Weight via GWAS in Paper Mulberry.

Seeds directly determine the survival and population size of woody plants, but the genetic basis of seed weight in woody plants remain poorly explored. To identify genetic variations and candidate genes responsible for seed weight in natural woody populations, we investigated the hundred-seed weight of 198 paper mulberry individuals from different areas. Our results showed that the hundred-seed weight of paper mulberry was significantly associated with the bioclimatic variables of sampling sites, which increased from south to north along the latitudinal-temperature gradient. Using 2,414,978 high-quality SNPs from re-sequencing data, the genome-wide association analysis of the hundred-seed weight was performed under three models, which identified 148, 19 and 12 associated genes, respectively. Among them, 25 candidate genes were directly hit by the significant SNPs, including the WRKY transcription factor, fatty acid desaturase, F-box protein, etc. Most importantly, we identified three crucial genetic variations in the coding regions of candidate genes (Bp02g2123, Bp01g3291 and Bp10g1642), and significant differences in the hundred-seed weight were detected among the individuals carrying different genotypes. Further analysis revealed that Bp02g2123 encoding a fatty acid desaturase (FAD) might be a key factor affecting the seed weight and local climate adaptation of woody plants. Furthermore, the genome-wide investigation and expression analysis of FAD genes were performed, and the results suggested that BpFADs widely expressed in various tissues and responded to multiple phytohormone and stress treatments. Overall, our study identifies valuable genetic variations and candidate genes, and provides a better understanding of the genetic basis of seed weight in woody plants.

Genome-Wide Analysis of BpYABs and Function Identification Involving in the Leaf and Silique Development in Transgenic Arabidopsis.

YABs play an important role in the leaf development of the paper mulberry (Broussonetia papyrifera) and of the heterophylly. Thus, we investigated the function of BpYABs. Gene cloning, phylogenetic analysis, motif identification, subcellular localization, transactivation activity assay, qRT-PCR, in situ hybridization, and ectopic expression were used in our study. Six BpYABs were isolated, and four of them had transcriptional activity. BpYAB1, BpYAB3, BpYAB4, and BpYAB5 were localized to the nucleus. BpYAB1 was only expressed in the flower, while BpYAB6 was not expressed in any detected tissues; the four remaining BpYABs were expressed in the bud, leaf and flower, and their expression level decreased with leaf development. Further in situ hybridization showed that BpYAB3 and BpYAB5 were expressed in the vascular tissues and lamina, but neither showed the adaxial-abaxial polarity distribution pattern in the mature leaf lamina. Ectopic expression of BpYAB2, BpYAB3, BpYAB4 and BpYAB5 induced increased expression of AtWOX1 and caused the leaf of Arabidopsis to become smaller and curl downwards. Ectopic expression also led to shorter siliques and smaller seeds, but not for BpYAB5. These results suggest that BpYABs have functional divergency and redundancy in regulating leaf and silique development.

Natural population re-sequencing detects the genetic basis of local adaptation to low temperature in a woody plant.

KEY MESSAGE: Total of 14 SNPs associated with overwintering-related traits and 75 selective regions were detected. Important candidate genes were identified and a possible network of cold-stress responses in woody plants was proposed. Local adaptation to low temperature is essential for woody plants to against changeable climate and safely survive the winter. To uncover the specific molecular mechanism of low temperature adaptation in woody plants, we sequenced 134 core individuals selected from 494 paper mulberry (Broussonetia papyrifera), which naturally distributed in different climate zones and latitudes. The population structure analysis, PCA analysis and neighbor-joining tree analysis indicated that the individuals were classified into three clusters, which showed forceful geographic distribution patterns because of the adaptation to local climate. Using two overwintering phenotypic data collected at high latitudes of 40°N and one bioclimatic variable, genome-phenotype and genome-environment associations, and genome-wide scans were performed. We detected 75 selective regions which possibly undergone temperature selection and identified 14 trait-associated SNPs that corresponded to 16 candidate genes (including LRR-RLK, PP2A, BCS1, etc.). Meanwhile, low temperature adaptation was also supported by other three trait-associated SNPs which exhibiting significant differences in overwintering traits between alleles within three geographic groups. To sum up, a possible network of cold signal perception and responses in woody plants were proposed, including important genes that have been confirmed in previous studies while others could be key potential candidates of woody plants. Overall, our results highlighted the specific and complex molecular mechanism of low temperature adaptation and overwintering of woody plants.

Identification of the PP2C gene family in paper mulberry (Broussonetia papyrifera) and its roles in the regulation mechanism of the response to cold stress.

OBJECTIVES: To study the possible roles of type-2C protein phosphatases (PP2Cs) which have been confirmed to play roles in the response to diverse abiotic stresses in paper mulberry, we launched a series of genomic and functional studies of BpPP2Cs. RESULTS: Sixty-three PP2C proteins in paper mulberry (Broussonetia papyrifera) were classified into 13 clades. Four BpPP2Cs with kinase domains were verified to be highly conserved in organisms ranging from algae to dicots. Seven pairs of BpPP2C genes were found to be expanding, and 18 BpPP2C genes had orthologues in Arabidopsis. BpPP2Cs showed broad expression in different tissues; the expression levels of 18 BpPP2Cs were changed and the phosphorylation levels of seven BpPP2C proteins increased at low temperature. Cold-response elements were found in the promoter region of 31 BpPP2Cs. Finally, Bp01g0320 was found to act as a hub protein and Bp01g0512 and Bp09g1278 played key roles in the ABA-signaling pathway and MAPK cascades, respectively. CONCLUSION: These results suggest that the PP2C gene family of paper mulberry is evolutionarily conserved and participates the regulation of the response to cold stress, which will play a vital role in further research on phosphatases in paper mulberry.

Cooperation between Broussonetia papyrifera and Its Symbiotic Fungal Community To Improve Local Adaptation of the Host.

The genetic basis of plant local adaptation has been extensively studied, yet the interplay between local adaptation, plant genetic divergence, and the microbial community remains unclear. Our study used the restriction-site associated DNA sequencing (RAD-seq) approach to explore genetic divergence in Broussonetia papyrifera and used internal transcribed spacers (ITS) to characterize fungal community. RAD-seq results show that B. papyrifera individuals could be divided into three genotypes; this genotyping result was consistent with the classification of climate type at the sample site. Most of the 101 highly differentiated genes were related to stress resistance and the microbiome. Moreover, ß-diversity results indicated that genetic divergence had a significant effect on fungal community across all compartments (P < 0.01). At genus and operational taxonomic unit (OTU) level, Mortierella, Hannaella oryzae, OTU81578 (Mortierella), and OTU1665209 (H. oryzae) were found to be the major OTUs that contribute to differences in fungal community. The properties of cooccurrence networks vary greatly among three genotypes. The results of redundancy analysis (RDA) indicated that B. papyrifera-associated fungal community was significantly related to its local adaptability. Our findings suggest that genetic divergence of B. papyrifera is closely related to local adaptation, with significant effects on the associated fungal community, which in turn would enhance host local adaptability. This improves present understanding about the coevolution of microbial communities and the host plant.IMPORTANCE The coevolution of plants with the associated fungal community and its effect on plant adaptability are not clear, especially for native trees. This study focuses on the genetic basis of local adaptation in plants and the effect of genetic divergence of Broussonetia papyrifera on the associated fungal community. We identified genes related to the microbiome that are important for local adaptation of the host. Our results show that genetic divergence in B. papyrifera significantly affects the fungal community, which has a close connection with local adaptation. This helps us to understand the relationship between local adaptation, genetic divergence, and associated fungal communities. This study highlights the effect of plant genetic divergence on associated fungal community for native trees and establishes a close connection between this effect and local adaptability in the host. In addition, these observations lay a foundation for the research of coevolution of plants and their symbiotic microbiome through genome-wide association study (GWAS).

The effect of plant compartments on the Broussonetia papyrifera-associated fungal and bacterial communities.

Plants associate with numerous microbes, but little is known about how microbiome components, especially fungi, adapt to specific plant compartments. The adaptability of microbial function to the plant compartment is also not clear especially for woody species. Here, we characterized the bacterial and fungal communities in root endosphere, stems, and rhizospheres of 33 Broussonetia papyrifera seedlings, based on amplification of 16S and ITS rRNA. Results showed that the α-diversity indexes of the bacterial community were significantly different in different plant compartments and they significantly increased from stem to root endosphere to the rhizosphere, whereas those of the fungal community were similar (p > 0.05). However, the result of constrained PCoA (CPCoA) and analysis of similarity (ANOSIM) showed that both bacterial and fungal compositions were significantly affected by plant compartments (p < 0.01). In detail, the operational taxonomic units (OTUs) distribution of the bacterial community was significantly different, but 249 of 252 fungal OTUs were shared in different plant compartments. Both the bacterial and fungal compositions were significantly influenced by plant compartments, based on the result on phyla, core OTUs, and indicator OTUs level. Further, 40 of 42 enriched KEGG pathways involving the bacteria also differed significantly among plant compartments (p < 0.01). This study provides an understanding of the influence of plant compartments on the microbiome and confirms that the disperse limitation of fungal OTUs across different plant compartments is smaller. This study sheds light on how the microbial community adapts to and thrives in different plant compartments.

The effect of environment on the microbiome associated with the roots of a native woody plant under different climate types in China.

Few studies have investigated the effect of environment on the root-associated microbiome, especially for woody plants in their native environment. The roots and rhizosphere soils of a native woody species (Broussonetia papyrifera) sampled across four different climate types in China were used to elucidate the influence of environment on the root-associated microbiome. Our results showed that the B. papyrifera root-associated microbiome contained abundant Proteobacteria and Basidiomycota, especially Pseudomonas and Rhizobium. The root-associated microbiomes were found to be significantly different under different climate types except for the bacterial community in the rhizosphere, and the proportion of bacterial operational taxonomic units (OTUs) shared among different climate types was lower than that of fungi. More than 50% of the total variance between microbiomes could be explained by 15 environmental factors, six of which, especially soil concentration phosphate and nitrate, had a significant effect. This study provided a comprehensive understanding of the root-associated microbiome of B. papyrifera and further confirmed the effect of environment on the root-associated microbiome of B. papyrifera under different climate types, with some exceptions in the rhizobacterial community and fungal OTUs. Our findings advanced knowledge of the effect of environment through an exploration of environmental factors and found that the nitrogen and phosphorus content represented the key factors.

Global Phosphoproteomic Analysis Reveals the Defense and Response Mechanisms of Jatropha Curcas Seedling under Chilling Stress.

As a promising energy plant for biodiesel, Jatropha curcas is a tropical and subtropical shrub and its growth is affected by one of major abiotic stress, chilling. Therefore, we adopt the phosphoproteomic analysis, physiological measurement and ultrastructure observation to illustrate the responsive mechanism of J. curcas seedling under chilling (4 °C) stress. After chilling for 6 h, 308 significantly changed phosphoproteins were detected. Prolonged the chilling treatment for 24 h, obvious physiological injury can be observed and a total of 332 phosphoproteins were examined to be significantly changed. After recovery (28 °C) for 24 h, 291 phosphoproteins were varied at the phosphorylation level. GO analysis showed that significantly changed phosphoproteins were mainly responsible for cellular protein modification process, transport, cellular component organization and signal transduction at the chilling and recovery periods. On the basis of protein-protein interaction network analysis, phosphorylation of several protein kinases, such as SnRK2, MEKK1, EDR1, CDPK, EIN2, EIN4, PI4K and 14-3-3 were possibly responsible for cross-talk between ABA, Ca2+, ethylene and phosphoinositide mediated signaling pathways. We also highlighted the phosphorylation of HOS1, APX and PIP2 might be associated with response to chilling stress in J. curcas seedling. These results will be valuable for further study from the molecular breeding perspective.

Rhizobia promote the growth of rice shoots by targeting cell signaling, division and expansion.

KEY MESSAGE: The growth-promotion of rice seedling following inoculation with Sinorhizobium meliloti 1021 was a cumulative outcome of elevated expression of genes that function in accelerating cell division and enhancing cell expansion. Various endophytic rhizobacteria promote the growth of cereal crops. To achieve a better understanding of the cellular and molecular bases of beneficial cereal-rhizobia interactions, we performed computer-assisted microscopy and transcriptomic analyses of rice seedling shoots (Oryza sativa) during early stages of endophytic colonization by the plant growth-promoting Sinorhizobium meliloti 1021. Phenotypic analyses revealed that plants inoculated with live rhizobia had increased shoot height and dry weight compared to control plants inoculated with heat-killed cells of the same microbe. At 6 days after inoculation (DAI) with live cells, the fourth-leaf sheaths showed significant cytological differences including their enlargement of parenchyma cells and reduction in shape complexity. Transcriptomic analysis of shoots identified 2,414 differentially-expressed genes (DEGs) at 1, 2, 5 and 8 DAI: 195, 1390, 1025 and 533, respectively. Among these, 46 DEGs encoding cell-cycle functions were up-regulated at least 3 days before the rhizobia ascended from the roots to the shoots, suggesting that rhizobia are engaged in long-distance signaling events during early stages of this plant-microbe interaction. DEGs involved in phytohormone production, photosynthetic efficiency, carbohydrate metabolism, cell division and wall expansion were significantly elevated at 5 and 8 DAI, consistent with the observed phenotypic changes in rice cell morphology and shoot growth-promotion. Correlation analysis identified 104 height-related DEGs and 120 dry-weight-related DEGs that represent known quantitative-trait loci for seedling vigor and increased plant height. These findings provide multiple evidences of plant-microbe interplay that give insight into the growth-promotion processes associated with this rhizobia-rice beneficial association.

Phosphoproteomic Analysis of Paper Mulberry Reveals Phosphorylation Functions in Chilling Tolerance.

Paper mulberry is a valuable woody species with a good chilling tolerance. In this study, phosphoproteomic analysis, physiological measurement, and mRNA quantification were employed to explore the molecular mechanism of chilling (4 °C) tolerance in paper mulberry. After chilling for 6 h, 427 significantly changed phosphoproteins were detected in paper mulberry seedlings without obvious physiological injury. When obvious physiological injury occurred after chilling for 48 h, a total of 611 phosphoproteins were found to be significantly changed at the phosphorylation level. Several protein kinases, especially CKII, were possibly responsible for these changes according to conserved sequence analysis. The results of Gene Ontology analysis showed that phosphoproteins were mainly responsible for signal transduction, protein modification, and translation during chilling. Additionally, transport and cellular component organization were enriched after chilling for 6 and 48 h, respectively. On the basis of the protein-protein interaction network analysis, a protein kinase and phosphatases hub protein (P1959) were found to be involved in cross-talk between Ca2+, BR, ABA, and ethylene-mediated signaling pathways. We also highlighted the phosphorylation of BpSIZ1 and BpICE1 possibly impacted on the CBF/DREB-responsive pathway. From these results, we developed a schematic for the chilling tolerance mechanism at phosphorylation level.

Identification and Functional Divergence Analysis of WOX Gene Family in Paper Mulberry.

The WOX (WUSCHEL-related homeobox) is a plant-specific transcription factor involved in plant development and stress response. However, few studies have been reported on the WOX gene in woody plants. In this study, 10 BpWOX genes were isolated from paper mulberry by RACE-PCR and categorized into three clades through phylogenetic analysis, ancient, intermediate and WUS clade. Among them, five members had the transcriptional activity detected by yeast one-hybrid and seven were uniquely localized to the nucleus through green fluorescent protein (GFP) observation. The expression patterns of BpWOX genes in different tissues and under diverse treatments were quantified by the qRT-PCR method. Results showed that BpWUS was expressed in the apical bud, stem and root, BpWOX5 and BpWOX7 functioned only in the root tip, and three BpWOXs regulated leaf development redundantly. BpWOX9 and BpWOX10 were induced by indole-3-acetic acid (IAA) or jasmonic acid (JA), while BpWOX2 was repressed by five phytohormones. Interestingly, most BpWOX genes were responsive to the abiotic stress stimuli of drought, salt, cold, and cadmium (CdCl2). Together, our study revealed that BpWOXs were functionally divergent during paper mulberry development and environmental adaptation, which might be related to their evolutionary relationships. Our work will benefit the systematic understanding of the precise function of WOX in plant development and environmental stress responses.

Genome-wide identification of the Jatropha curcas MYB family and functional analysis of the abiotic stress responsive gene JcMYB2.

BACKGROUND: The MYB family is one of the most abundant transcription factor families in plants. MYB proteins are involved in plant development, abiotic stress tolerance, hormone signal transduction and disease resistance. Here we perform genome-wide identification of MYB family transcription factors in an energy plant J. curcas, including determining family composition, phylogenetic evolution and functional prediction analysis. In addition, we further elucidate the function of the JcMYB2 gene. METHODS: The phylogenetic trees were constructed by using the neighbor-joining method in MEGA 5.2. The biological functions of some JcMYBs were predicted according to orthology. The full length cDNA of JcMYB2 was cloned by using the RACE method. GUS histochemical staining was used to test the activity of the JcMYB2 promoter. Expression patterns of JcMYB2 were detected by using qPCR Transcriptional activity JcMYB2 were confirmed through yeast one hybrid. Subcellular Localization of JcMYB2 Protein were demonstrated by transient expression in the tobacco leaf. The function of JcMYB2 in salt and freezing tolerance were detected in transgenic plants. RESULTS: A genome-wide analysis identified 128 MYB genes, including 123 R2R3-MYBs, 4 R1R2R3-MYBs and 1 4R-MYB. All of the R2R3-MYBs are further classified into 19 groups which indicated functional conservation among previously identified groups of R2R3-MYB proteins. Among of these newly identified MYBs, the JcMYB2 belongs to group G11 and its expression is induced obviously by cold, salt and MeJA (Methyl Jasmonate) and slightly by ABA (abscisic acid). JcMYB2 is localized to the nucleus and has transcriptional activity. JcMYB2 overexpressing plants are more tolerant to salt and cold stress than wild type plants. Tissue specific expression profiles showed that the JcMYB2 gene was expressed ubiquitously throughout the plant, with higher expression levels observed in the root. CONCLUSION: A comprehensive genome-wide analysis and phylogenetic relationship of R2R3-MYB subfamily in J. curcas present the global identification and functional prediction of JcR2R3-MYBs. Additionally, JcMYB2 regulates the stress response signaling networks by interacting with MeJA and ABA signaling pathway and functions in the root development of J. curcas.

The cold responsive mechanism of the paper mulberry: decreased photosynthesis capacity and increased starch accumulation.

BACKGROUND: Most studies on the paper mulberry are mainly focused on the medicated and pharmacology, fiber quality, leaves feed development, little is known about its mechanism of adaptability to abiotic stress. Physiological measurement, transcriptomics and proteomic analysis were employed to understand its response to cold stress in this study. METHODS: The second to fourth fully expanded leaves from up to down were harvested at different stress time points forthe transmission electron microscope (TEM) observation. Physiological characteristics measurement included the relative electrolyte leakage (REL), SOD activity assay, soluble sugar content, and Chlorophyll fluorescence parameter measurement. For screening of differentially expressed genes, the expression level of every transcript in each sample was calculated by quantifying the number of Illumina reads. To identify the differentially expressed protein, leaves of plants under 0, 6, 12, 24, 48 and 72 h cold stress wereharvested for proteomic analysis. Finally, real time PCR was used to verify the DEG results of the RNA-seq and the proteomics data. RESULTS: Results showed that at the beginning of cold stress, respiratory metabolism was decreased and the transportation and hydrolysis of photosynthetic products was inhibited, leading to an accumulation of starch in the chloroplasts. Total of 5800 unigenes and 38 proteins were affected, including the repressed expression of photosynthesis and the enhanced expression in signal transduction, stress defense pathway as well as secondary metabolism. Although the transcriptional level of a large number of genes has been restored after 12 h, sustained cold stress brought more serious injury to the leaf cells, including the sharp rise of the relative electrolyte leakage, the declined Fv/Fm value, swelled chloroplast and the disintegrated membrane system. CONCLUSION: The starch accumulation and the photoinhibition might be the main adaptive mechanism of the paper mulberry responded to cold stress. Most of important, enhancing the transport and hydrolysis of photosynthetic products could be the potential targets for improving the cold tolerance of the paper mulberry.

Transcriptional regulation of the paper mulberry under cold stress as revealed by a comprehensive analysis of transcription factors.

BACKGROUND: Several studies have focused on cold tolerance in multiple regulated levels. However, a genome-scale molecular analysis of the regulated network under the control of transcription factors (TFs) is still lacking, especially for trees. To comprehensively identify the TFs that regulate cold stress response in the paper mulberry and understand their regulatory interactions, transcriptomic data was used to assess changes in gene expression induced by exposure to cold. RESULTS: Results indicated that 794 TFs, belonging to 47 families and comprising more than 59% of the total TFs of this plant, were involved in the cold stress response. They were clustered into three groups, namely early, intermediate and late responsive groups which contained 95, 550 and 149 TFs, respectively. Among of these differentially expressed TFs, one bHLH, two ERFs and three CAMTAs were considered to be the key TFs functioning in the primary signal transduction. After that, at the intermediate stage of cold stress, there were mainly two biological processes that were regulated by TFs, namely cold stress resistance (including 5 bHLH, 14 ERFs, one HSF, 4 MYBs, 3 NACs, 11 WRKYs and so on) and growth and development of lateral organ or apical meristem (including ARR-B, B3, 5 bHLHs, 2 C2H2, 4 CO-like, 2 ERF, 3 HD-ZIP, 3 YABBYs, G2-like, GATA, GRAS and TCP). In late responsive group, 3 ARR-B, C3H, 6 CO-like, 2 G2-like, 2 HSFs, 2 NACs and TCP. Most of them presented the up-regulated expression at 12 or 24 hours after cold stress implied their important roles for the new growth homeostasis under cold stress. CONCLUSIONS: Our study identified the key TFs that function in the regulatory cascades mediating the activation of downstream genes during cold tress tolerance in the paper mulberry. Based on the analysis, we found that the AP2/ERF, bHLH, MYB, NAC and WRKY families might play the central and significant roles during cold stress response in the paper mulberry just as in other species. Meanwhile, many other TF families previously reported as involving in regulation of growth and development, including ARF, DBB, G2-like, GRF, GRAS, LBD, WOX and YAABY exhibited their important potential function in growth regulation under cold stress.

Bovine serum albumin in saliva mediates grazing response in Leymus chinensis revealed by RNA sequencing.

BACKGROUND: Sheepgrass (Leymus chinensis) is an important perennial forage grass across the Eurasian Steppe and is adaptable to various environmental conditions, but little is known about its molecular mechanism responding to grazing and BSA deposition. Because it has a large genome, RNA sequencing is expensive and impractical except for the next-generation sequencing (NGS) technology. RESULTS: In this study, NGS technology was employed to characterize de novo the transcriptome of sheepgrass after defoliation and grazing treatments and to identify differentially expressed genes (DEGs) responding to grazing and BSA deposition. We assembled more than 47 M high-quality reads into 120,426 contigs from seven sequenced libraries. Based on the assembled transcriptome, we detected 2,002 DEGs responding to BSA deposition during grazing. Enrichment analysis of Gene ontology (GO), EuKaryotic Orthologous Groups (KOG) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways revealed that the effects of grazing and BSA deposition involved more apoptosis and cell oxidative changes compared to defoliation. Analysis of DNA fragments, cell oxidative factors and the lengths of leaf scars after grazing provided physiological and morphological evidence that BSA deposition during grazing alters the oxidative and apoptotic status of cells. CONCLUSIONS: This research greatly enriches sheepgrass transcriptome resources and grazing-stress-related genes, helping us to better understand the molecular mechanism of grazing in sheepgrass. The grazing-stress-related genes and pathways will be a valuable resource for further gene-phenotype studies.

The large-scale investigation of gene expression in Leymus chinensis stigmas provides a valuable resource for understanding the mechanisms of poaceae self-incompatibility.

BACKGROUND: Many Poaceae species show a gametophytic self-incompatibility (GSI) system, which is controlled by at least two independent and multiallelic loci, S and Z. Until currently, the gene products for S and Z were unknown. Grass SI plant stigmas discriminate between pollen grains that land on its surface and support compatible pollen tube growth and penetration into the stigma, whereas recognizing incompatible pollen and thus inhibiting pollination behaviors. Leymus chinensis (Trin.) Tzvel. (sheepgrass) is a Poaceae SI species. A comprehensive analysis of sheepgrass stigma transcriptome may provide valuable information for understanding the mechanism of pollen-stigma interactions and grass SI. RESULTS: The transcript abundance profiles of mature stigmas, mature ovaries and leaves were examined using high-throughput next generation sequencing technology. A comparative transcriptomic analysis of these tissues identified 1,025 specifically or preferentially expressed genes in sheepgrass stigmas. These genes contained a significant proportion of genes predicted to function in cell-cell communication and signal transduction. We identified 111 putative transcription factors (TFs) genes and the most abundant groups were MYB, C2H2, C3H, FAR1, MADS. Comparative analysis of the sheepgrass, rice and Arabidopsis stigma-specific or preferential datasets showed broad similarities and some differences in the proportion of genes in the Gene Ontology (GO) functional categories. Potential SI candidate genes identified in other grasses were also detected in the sheepgrass stigma-specific or preferential dataset. Quantitative real-time PCR experiments validated the expression pattern of stigma preferential genes including homologous grass SI candidate genes. CONCLUSIONS: This study represents the first large-scale investigation of gene expression in the stigmas of an SI grass species. We uncovered many notable genes that are potentially involved in pollen-stigma interactions and SI mechanisms, including genes encoding receptor-like protein kinases (RLK), CBL (calcineurin B-like proteins) interacting protein kinases, calcium-dependent protein kinase, expansins, pectinesterase, peroxidases and various transcription factors. The availability of a pool of stigma-specific or preferential genes for L. chinensis offers an opportunity to elucidate the mechanisms of SI in Poaceae.

Integrating genome-wide association and transcriptome analysis to provide molecular insights into heterophylly and eco-adaptability in woody plants.

Heterophylly is regard as an important adaptive mechanism in response to different environments within plants. However, the genetic mechanisms responsible for heterophylly in woody plants are still poorly understood. Herein, the divergence of heterophyllous leaves was investigated at morphogenesis and using microdissection and physiological indexes in paper mulberry, and the genetic basis of heterophylly was further revealed combined with genome-wide association study (GWAS), transcriptome analysis and weighted gene coexpression network analysis (WGCNA). Our results revealed that the flavonoid content and antioxidant activity increased gradually from the entire leaf to the palmatisect leaf, while the hormone content and net photosynthetic rate decreased. Through GWAS and transcriptome analysis, a total of 98 candidate genes and 2338 differentially expressed genes associated with heterophylly were identified. Importantly, we uncovered critical variations in the candidate genes Bp07g0981 (WOX) and Bp07g0920 (HHO), along with significant differences in haplotypes and expression levels among heterophyllous leaves. Our results also suggested that the genes involved in hormone signaling pathways, antioxidant activity, and flavonoid metabolism might be closely related to the heterophylly of paper mulberry, which could account for the physiological data. Indeed, CR-wox mutant lines showed significant changes in leaf phenotypes, and differential expression profile analysis also highlighted the expression of genes related to phytohormones and transcription factors. Together, the genetic variations and candidate genes detected in this study provide novel insights into the genetic mechanism of heterophylly, and would improve the understanding of eco-adaptability in heterophyllous woody plants.

Identification of leucine-rich repeat receptor-like protein kinase (LRR-RLK) genes in paper mulberry and their potential roles in response to cold stress.

Leucine-rich repeat receptor-like protein kinases (LRR-RLKs) represent the largest group of receptor-like kinases in plants, which have been previously reported to play vital roles in plant growth, development, stress adaptation and signal transduction. However, there is lack of comprehensive analysis of this family in paper mulberry (Broussonetia papyrifera). In the present investigation, a genome-wide scan revealed the presence of 236 LRR-RLK genes in paper mulberry, which were classified into 21 subgroups based on the maximum-likelihood phylogenetic tree. Gene structure and conserved motif analyses suggested genes in the same subgroup had highly consistent motif composition and intron/exon arrangement, but were divergent among subgroups. Total of 223 BpLRR-RLK genes were unevenly distributed across all 13 chromosomes, while the remaining 13 genes were localized to the unassembled scaffolds. Tandem and segmental duplications were confirmed to contribute to the expansion of BpLRR-RLK family. Further Ka/Ks showed that the duplicated BpLRR-RLKs had experienced strong purifying selection. The global promoter composition, transcriptome and phosphorylation analysis indicated that many of BpLRR-RLKs were associated with plant development, biotic and abiotic stress response, especially for cold stress. Furthermore, protein-protein interaction network was constructed for the 127 and 14 BpLRR-RLKs that responded to cold stress at the transcriptomics and phosphorylation level, respectively. All these findings will facilitate the studies on the evolutionary history of the LRR-RLK gene family in paper mulberry, also establish a solid foundation to further explore the potential functions of LRR-RLK genes in higher plants, particularly with regards to cold resistance.