Transcript profiling reveals expression differences in wild-type and glabrous soybean lines

Transcriptomic and photosynthetic responses to grafting of the Nod1 gene in nodulated and non-nodulated soybeans

Legume plants form symbiotic relationships with rhizobia to convert N2 into ammonia, and the nodulation status can affect plant development including photosynthesis. However, the relationship between nitrogen fixation and photosynthesis during carbon and nitrogen metabolism remains unclear. This study was undertaken to unravel regulation of nodulation and photosynthesis using a spontaneous nonnodulated soybean mutant by grafting. The results of inheritance and gene mapping showed that the nonnodulated mutant was controlled by a recessive gene overlapped with the reported rj1 locus, and might be a new rj1 allele with 1 bp deletion in the fourth exon in comparison to the sequence of normal nodulation plants. According to grafting results, soybean nodulation is obviously determined by the roots, not the seedlings. Moreover, nitrogen content along with related metabolic enzyme activity, and photosynthetic capacity were enhanced by nonnodulated scions grafted with nodulated roots. Contrary results were obtained for nodulated scions grafted with nonnodulated roots. A total of 853 differentially expressed genes (DEGs) in the leaves and 1874 in the roots were identified by transcriptome analyses of the grafting treatments. We identified 285 differential gene ontology (GO) terms and 57 differential pathway terms identified in the leaves, while 856 differential GO terms and 207 differential pathway terms in the roots. Twenty DEGs interacting at translation level were selected, and the results of transcriptome analyses were verified by q-PCR. These findings indicated that the nodulation-related Nod allelic gene increases the nitrogen content of nonnodulated plants, which affects the enzymes involved in nitrogen metabolism, leading to changes in hormone levels and further regulation of photosynthesis and carbon metabolism.

Non-glandular trichomes of sunflower are important in the absorption and translocation of foliar-applied Zn

Trichomes are potentially important for absorption of foliar fertilizers. A study has shown that the non-glandular trichromes (NGTs) of sunflower (Helianthus annuus) accumulated high concentrations of foliar-applied zinc (Zn); however, the mechanisms of Zn accumulation in the NGTs and the fate of this Zn are unclear. Here we investigated how foliar-applied Zn accumulates in the NGTs and the subsequent translocation of this Zn. Time-resolved synchrotron-based X-ray fluorescence microscopy and transcriptional analyses were used to probe the movement of Zn in the NGTs, with the cuticle composition of the NGTs examined using confocal Raman microscopy. The accumulation of Zn in the NGTs is both an initial preferential absorption process and a subsequent translocation process. This preferred absorption is likely because the NGT base has a higher hydrophilicity, whilst the subsequent translocation is due to the presence of plasmodesmata, Zn-chelating ligands, and Zn transporters in the NGTs. Furthermore, the Zn sequestered in the NGTs was eventually translocated out of the trichome once the leaf Zn concentration had decreased, suggesting that the NGTs are also important in maintaining leaf Zn homeostasis. This study demonstrates for the first time that trichomes have a key structural and functional role in the absorption and translocation of foliar-applied Zn.

Multi-omic dissection of the drought resistance traits of soybean landrace LX

With diverse genetic backgrounds, soybean landraces are valuable resource for breeding programs. Herein, we apply multi-omic approaches to extensively characterize the molecular basis of drought tolerance in the soybean landrace LX. Initial screens established that LX performed better with PEG6000 treatment than control cultivars. LX germinated better than William 82 under drought conditions and accumulated more anthocyanin and flavonoids. Untargeted mass spectrometry in combination with transcriptomic analyses revealed the chemical diversity and genetic basis underlying the overall performance of LX landrace. Under control and drought conditions, significant differences in the expression of a suite of secondary metabolism genes, particularly those involved in the general phenylpropanoid pathway and flavonoid but not lignin biosynthesis, were seen in LX and William 82. The expression of these genes correlated with the corresponding metabolites in LX plants. Further correlation analysis between metabolites and transcripts identified pathway structural genes and transcription factors likely are responsible for the LX agronomic traits. The activities of some key biosynthetic genes or regulators were confirmed through heterologous expression in transgenic Arabidopsis and hairy root transformation in soybean. We propose a regulatory mechanism based on flavonoid secondary metabolism and adaptive traits of this landrace which could be of relevance to cultivated soybean.

Automated trichome counting in soybean using advanced image-processing techniques

PREMISE

Trichomes are hair-like appendages extending from the plant epidermis. They serve many important biotic roles, including interference with herbivore movement. Characterizing the number, density, and distribution of trichomes can provide valuable insights on plant response to insect infestation and define the extent of plant defense capability. Automated trichome counting would speed up this research but poses several challenges, primarily because of the variability in coloration and the high occlusion of the trichomes.

METHODS AND RESULTS

We developed a simplified method for image processing for automated and semi-automated trichome counting. We illustrate this process using 30 leaves from 10 genotypes of soybean (Glycine max) differing in trichome abundance. We explored various heuristic image-processing methods including thresholding and graph-based algorithms to facilitate trichome counting. Of the two automated and two semi-automated methods for trichome counting tested and with the help of regression analysis, the semi-automated manually annotated trichome intersection curve method performed best, with an accuracy of close to 90% compared with the manually counted data.

CONCLUSIONS

We address trichome counting challenges including occlusion by combining image processing with human intervention to propose a semi-automated method for trichome quantification. This provides new opportunities for the rapid and automated identification and quantification of trichomes, which has applications in a wide variety of disciplines.

Transcription factors and glyoxylate cycle genes prominent in the transition of soybean cotyledons to the first functional leaves of the seedling

During early seedling growth, the cotyledons transition from a storage tissue to become effectively the first leaf-like structures of the plant. In this programmed developmental process, they likely undergo a massive change in gene expression to redirect their metabolism and physiological processes. To define the developmental shifts in gene expression and begin to understand the gene regulatory networks that set this transition in motion, we carried out high-throughput RNA sequencing of cotyledons from seven developmental stages of soybean seedlings. We identified 154 gene models with high expression exclusively in the early seedling stages. A significant number (about 25 %) of those genes with known annotations were involved in carbohydrate metabolism. A detailed examination of glyoxylate cycle genes revealed the upregulation of their expression in the early stages of development. A total of approximately 50 % of the highly expressed genes whose expression peaked in the mid-developmental stages encoded ribosomal family proteins. Our analysis also identified 219 gene models with high expression at late developmental stages. The majority of these genes are involved in photosynthesis, including photosystem I- and II-associated genes. Additionally, the advantage of RNA-Seq to detect genes expressed at low levels revealed approximately 460 transcription factors with notable expression in at least one stage of the developing soybean seedling. Relatively over-represented transcription factor genes encode AP2, zinc finger, NAC, WRKY, and MYB families. These transcription factor genes may lead to the transcriptional reprogramming during the transition of seedling cotyledons from storage tissue to metabolically active organs that serve as the first functional leaves of the plant.

Bayesian phylogeny of sucrose transporters: ancient origins, differential expansion and convergent evolution in monocots and dicots

Sucrose transporters (SUTs) are essential for the export and efficient movement of sucrose from source leaves to sink organs in plants. The angiosperm SUT family was previously classified into three or four distinct groups, Types I, II (subgroup IIB), and III, with dicot-specific Type I and monocot-specific Type IIB functioning in phloem loading. To shed light on the underlying drivers of SUT evolution, Bayesian phylogenetic inference was undertaken using 41 sequenced plant genomes, including seven basal lineages at key evolutionary junctures. Our analysis supports four phylogenetically and structurally distinct SUT subfamilies, originating from two ancient groups (AG1 and AG2) that diverged early during terrestrial colonization. In both AG1 and AG2, multiple intron acquisition events in the progenitor vascular plant established the gene structures of modern SUTs. Tonoplastic Type III and plasmalemmal Type II represent evolutionarily conserved descendants of AG1 and AG2, respectively. Type I and Type IIB were previously thought to evolve after the dicot-monocot split. We show, however, that divergence of Type I from Type III SUT predated basal angiosperms, likely associated with evolution of vascular cambium and phloem transport. Type I SUT was subsequently lost in monocots along with vascular cambium, and independent evolution of Type IIB coincided with modified monocot vasculature. Both Type I and Type IIB underwent lineage-specific expansion. In multiple unrelated taxa, the newly-derived SUTs exhibit biased expression in reproductive tissues, suggesting a functional link between phloem loading and reproductive fitness. Convergent evolution of Type I and Type IIB for SUT function in phloem loading and reproductive organs supports the idea that differential vascular development in dicots and monocots is a strong driver for SUT family evolution in angiosperms.

Study of seed hair growth in Populus tomentosa, an important character of female floral bud development

Background

Poplar seed hair is an environmental annoyance in northern China due to its abundance and widespread airborne distribution after maturation. The morphogenesis and molecular mechanisms of its development are not well understood, and little attention has been focused on the dynamics of its development. To better understand the mechanism of poplar seed hair development, paraffin sections were used to examine the initiation and elongation of poplar seed hairs. RNA-seq technology was also employed to provide a comprehensive overview of transcriptional changes that occur during seed hair development.

Results

The placenta at the base of ovary, was identified as the origin of seed hair development, which is in sharp contrast to cotton fibers that originate from epidermal cells of the seed coat. An enlarged cell nucleus in seed hair cells was also observed, which was supported by our gene ontology enrichment analysis. The significant enriched GO term of “endoreduplication” indicated that cycles of endoreduplication, bypassing normal mitosis, is the underlying mechanisms for the maintenance of the uni-cellular structure of seed hairs. By analyzing global changes in the transcriptome, many genes regulating cell cycle, cell elongation, cell well modification were identified. Additionally, in an analysis of differential expression, cellulose synthesis and cell wall biosynthesis-related biological processes were enriched, indicating that this component of fiber structure in poplar seed hairs is consistent with what is found in cotton fibers. Differentially expressed transcription factors exhibited a stage-specific up-regulation. A dramatic down-regulation was also revealed during the mid-to-late stage of poplar seed hair development, which may point to novel mechanisms regulating cell fate determination and cell elongation.

Conclusions

This study revealed the initiation site of poplar seed hairs and also provided a comprehensive overview of transcriptome dynamics during the process of seed hair development. The high level of resolution on dynamic changes in the transcriptome provided in this study may serve as a valuable resource for developing a more complete understanding of this important biological process.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-475) contains supplementary material, which is available to authorized users.

RNA-Seq profiling of a defective seed coat mutation in Glycine max reveals differential expression of proline-rich and other cell wall protein transcripts

The plant cell wall performs a number of essential functions including providing shape to many different cell types and serving as a defense against potential pathogens. The net pattern mutation creates breaks in the seed coat of soybean (Glycine max) because of ruptured cell walls. Using RNA-Seq, we examined the seed coat transcriptome from three stages of immature seed development in two pairs of isolines with normal or defective seed coat phenotypes due to the net pattern. The genome-wide comparative study of the transcript profiles of these isolines revealed 364 differentially expressed genes in common between the two varieties that were further divided into different broad functional categories. Genes related to cell wall processes accounted for 19% of the differentially expressed genes in the middle developmental stage of 100-200 mg seed weight. Within this class, the cell wall proline-rich and glycine-rich protein genes were highly differentially expressed in both genetic backgrounds. Other genes that showed significant expression changes in each of the isoline pairs at the 100-200 mg seed weight stage were xylem serine proteinase, fasciclin-related genes, auxin and stress response related genes, TRANSPARENT TESTA 1 (TT1) and other transcription factors. The mutant appears to shift the timing of either the increase or decrease in the levels of some of the transcripts. The analysis of these data sets reveals the physiological changes that the seed coat undergoes during the formation of the breaks in the cell wall.

Early transcriptional response of soybean contrasting accessions to root dehydration

Drought is a significant constraint to yield increase in soybean. The early perception of water deprivation is critical for recruitment of genes that promote plant tolerance. DeepSuperSAGE libraries, including one control and a bulk of six stress times imposed (from 25 to 150 min of root dehydration) for drought-tolerant and sensitive soybean accessions, allowed to identify new molecular targets for drought tolerance. The survey uncovered 120,770 unique transcripts expressed by the contrasting accessions. Of these, 57,610 aligned with known cDNA sequences, allowing the annotation of 32,373 unitags. A total of 1,127 unitags were up-regulated only in the tolerant accession, whereas 1,557 were up-regulated in both as compared to their controls. An expression profile concerning the most representative Gene Ontology (GO) categories for the tolerant accession revealed the expression "protein binding" as the most represented for "Molecular Function", whereas CDPK and CBL were the most up-regulated protein families in this category. Furthermore, particular genes expressed different isoforms according to the accession, showing the potential to operate in the distinction of physiological behaviors. Besides, heat maps comprising GO categories related to abiotic stress response and the unitags regulation observed in the expression contrasts covering tolerant and sensitive accessions, revealed the unitags potential for plant breeding. Candidate genes related to "hormone response" (LOX, ERF1b, XET), "water response" (PUB, BMY), "salt stress response" (WRKY, MYB) and "oxidative stress response" (PER) figured among the most promising molecular targets. Additionally, nine transcripts (HMGR, XET, WRKY20, RAP2-4, EREBP, NAC3, PER, GPX5 and BMY) validated by RT-qPCR (four different time points) confirmed their differential expression and pointed that already after 25 minutes a transcriptional reorganization started in response to the new condition, with important differences between both accessions.

Functional genomics to study stress responses in crop legumes: progress and prospects

Legumes are important food crops worldwide, contributing to more than 33% of human dietary protein. The production of crop legumes is frequently impacted by abiotic and biotic stresses. It is therefore important to identify genes conferring resistance to biotic stresses and tolerance to abiotic stresses that can be used to both understand molecular mechanisms of plant response to the environment and to accelerate crop improvement. Recent advances in genomics offer a range of approaches such as the sequencing of genomes and transcriptomes, gene expression microarray as well as RNA-seq based gene expression profiling, and map-based cloning for the identification and isolation of biotic and abiotic stress-responsive genes in several crop legumes. These candidate stress associated genes should provide insights into the molecular mechanisms of stress tolerance and ultimately help to develop legume varieties with improved stress tolerance and productivity under adverse conditions. This review provides an overview on recent advances in the functional genomics of crop legumes that includes the discovery as well as validation of candidate genes.

Genome-wide expression profiles of seasonal bud dormancy at four critical stages in Japanese apricot

Dormancy is one of the most important adaptive mechanisms developed by perennial plants. To reveal the comprehensive mechanism of seasonal bud dormancy at four critical stages in Japanese apricot (Prunus persica), we applied Illumina sequencing to study differentially expressed genes (DEGs) at the transcriptional level. As a result, 19,759, 16,375, 19,749 and 20,800 tag-mapped genes were sequenced from libraries of paradormancy (R1), endodormancy (R2), ecodormancy (R3) and dormancy release (R4) stages based on the P. persica genome. Moreover, 6,199, 5,539, and 5,317 genes were differentially expressed in R1 versus R2, R2 versus R3, and R3 versus R4, respectively. Gene Ontology analysis of dormancy-related genes showed that these were mainly related to the cytoplasm, cytoplasmic part metabolism, intracellular metabolism and membrane-bound organelle metabolism. Pathway-enrichment annotation revealed that highly ranked genes were involved in ribosome pathways and protein processing in the endoplasmic reticulum. The results demonstrated that hormone response genes such as auxin, abscisic acid, ethylene and jasmonic acid, as well as zinc finger family protein genes are possibly involved in seasonal bud dormancy in Japanese apricot. The expression patterns of DEGs were verified using real-time quantitative RT-PCR. These results contribute to further understanding of the mechanism of bud dormancy in Japanese apricot.

Using RNA-Seq to profile soybean seed development from fertilization to maturity

To understand gene expression networks leading to functional properties and compositional traits of the soybean seed, we have undertaken a detailed examination of soybean seed development from a few days post-fertilization to the mature seed using Illumina high-throughput transcriptome sequencing (RNA-Seq). RNA was sequenced from seven different stages of seed development, yielding between 12 million and 78 million sequenced transcripts. These have been aligned to the 79,000 gene models predicted from the soybean genome recently sequenced by the Department of Energy Joint Genome Institute. Over one hundred gene models were identified with high expression exclusively in young seed stages, starting at just four days after fertilization. These were annotated as being related to many basic components and processes such as histones and proline-rich proteins. Genes encoding storage proteins such as glycinin and beta-conglycinin had their highest expression levels at the stages of largest fresh weight, confirming previous knowledge that these storage products are being rapidly accumulated before the seed begins the desiccation process. Other gene models showed high expression in the dry, mature seeds, perhaps indicating the preparation of pathways needed later, in the early stages of imbibition. Many highly-expressed gene models at the dry seed stage are, as expected, annotated as hydrophilic proteins associated with low water conditions, such as late embryogenesis abundant (LEA) proteins and dehydrins, which help preserve the cellular structures and nutrients within the seed during desiccation. More significantly, the power of RNA-Seq to detect genes expressed at low levels revealed hundreds of transcription factors with notable expression in at least one stage of seed development. Results from a second biological replicate demonstrate high reproducibility of these data revealing a comprehensive view of the transciptome of seed development in the cultivar Williams, the reference cultivar for the first soybean genome sequence.

Identification of differentially-expressed genes associated with pistil abortion in Japanese apricot by genome-wide transcriptional analysis

The phenomenon of pistil abortion widely occurs in Japanese apricot, and imperfect flowers with pistil abortion seriously decrease the yield in production. Although transcriptome analyses have been extensively studied in the past, a systematic study of differential gene expression has not been performed in Japanese apricot. To investigate genes related to the pistil development of Japanese apricot, high-throughput sequencing technology (Illumina) was employed to survey gene expression profiles from perfect and imperfect Japanese apricot flower buds. 3,476,249 and 3,580,677 tags were sequenced from two libraries constructed from perfect and imperfect flower buds of Japanese apricot, respectively. There were 689 significant differentially-expressed genes between the two libraries. GO annotation revealed that highly ranked genes were those implicated in small molecule metabolism, cellular component organisation or biogenesis at the cellular level and fatty acid metabolism. According to the results, we assumed that late embryogenesis abundant protein (LEA), Dicer-like 3 (DCL3) Xyloglucan endotransglucosylase/hydrolase 2 (XTH2), Pectin lyase-like superfamily protein (PPME1), Lipid transfer protein 3 (LTP3), Fatty acid biosynthesis 1 (FAB1) and Fatty acid desaturase 5 (FAD5) might have relationships with the pistil abortion in Japanese apricot. The expression patterns of 36 differentially expressed genes were confirmed by real-time (RT)-PCR. This is the first report of the Illumina RNA-seq technique being used for the analysis of differentially-expressed gene profiles related to pistil abortion that both computationally and experimentally provides valuable information for the further functional characterisation of genes associated with pistil development in woody plants.