Marc Libault

Soybean (Glycine max L.) is a major crop providing an important source of protein and oil, which can also be converted into biodiesel. A major milestone in soybean research was the recent sequencing of its genome. The sequence predicts 69,145 putative soybean genes, with 46,430 predicted with high confidence. In order to examine the expression of these genes, we utilized the Illumina Solexa platform to sequence cDNA derived from 14 conditions (tissues). The result is a searchable soybean gene expression atlas accessible through a browser (http://digbio.missouri.edu/soybean_atlas). The data provide experimental support for the transcription of 55,616 annotated genes and also demonstrate that 13,529 annotated soybean genes are putative pseudogenes, and 1736 currently unannotated sequences are transcribed. An analysis of this atlas reveals strong differences in gene expression patterns between different tissues, especially between root and aerial organs, but also reveals similarities between gene expression in other tissues, such as flower and leaf organs. In order to demonstrate the full utility of the atlas, we investigated the expression patterns of genes implicated in nodulation, and also transcription factors, using both the Solexa sequence data and large-scale qRT-PCR. The availability of the soybean gene expression atlas allowed a comparison with gene expression documented in the two model legume species, Medicago truncatula and Lotus japonicus, as well as data available for Arabidopsis thaliana, facilitating both basic and applied aspects of soybean research.

Gene expression analysis requires the use of reference genes constitutively expressed independently of tissues or environmental conditions. Housekeeping genes (e.g., actin, tubulin, ribosomal, polyubiquitin, and elongation factor 1‐α) are commonly used as reference genes with the assumption that they are uniformly expressed. In many cases, however, this assumption was shown to be incorrect. To provide reliable reference genes in soybean [ Glycine max (L.)], we surveyed a set of genes that showed little variation in a nodulation study across a series of soybean microarray experiments. More than 200 putative reference genes were identified. We focused on 18 for further analysis using additional cDNA and Affymetrix arrays and quantitative reverse‐transcription polymerase chain reactions. Taken together, these experiments allowed us to test the expression stability of these genes in 130 different conditions, confirming four soybean genes as new reference genes (annotated as ATP‐binding cassette [ABC] transporter, F‐box protein family, metalloprotease, and CDPK‐related protein kinase). These genes should be useful for normalization of gene expression studies in soybean, an important crop plant.

Polyploidy is generally not tolerated in animals, but is widespread in plant genomes and may result in extensive genetic redundancy. The fate of duplicated genes is poorly understood, both functionally and evolutionarily. Soybean (Glycine max L.) has undergone two separate polyploidy events (13 and 59 million years ago) that have resulted in 75% of its genes being present in multiple copies. It therefore constitutes a good model to study the impact of whole-genome duplication on gene expression. Using RNA-seq, we tested the functional fate of a set of approximately 18 000 duplicated genes. Across seven tissues tested, approximately 50% of paralogs were differentially expressed and thus had undergone expression sub-functionalization. Based on gene ontology and expression data, our analysis also revealed that only a small proportion of the duplicated genes have been neo-functionalized or non-functionalized. In addition, duplicated genes were often found in collinear blocks, and several blocks of duplicated genes were co-regulated, suggesting some type of epigenetic or positional regulation. We also found that transcription factors and ribosomal protein genes were differentially expressed in many tissues, suggesting that the main consequence of polyploidy in soybean may be at the regulatory level.

Similar to other complex organisms, plants consist of diverse and specialized cell types. The gain of unique biological functions of these different cell types is the consequence of the establishment of cell-type-specific transcriptional programs. As a necessary step in gaining a deeper understanding of the regulatory mechanisms controlling plant gene expression, we report the use of single-nucleus RNA sequencing (sNucRNA-seq) and single-nucleus assay for transposase accessible chromatin sequencing (sNucATAC-seq) technologies on Arabidopsis roots. The comparison of our single-nucleus transcriptomes to the published protoplast transcriptomes validated the use of nuclei as biological entities to establish plant cell-type-specific transcriptomes. Furthermore, our sNucRNA-seq results uncovered the transcriptomes of additional cell subtypes not identified by single-cell RNA-seq. Similar to our transcriptomic approach, the sNucATAC-seq approach led to the distribution of the Arabidopsis nuclei into distinct clusters, suggesting the differential accessibility of chromatin between groups of cells according to their identity. To reveal the impact of chromatin accessibility on gene expression, we integrated sNucRNA-seq and sNucATAC-seq data and demonstrated that cell-type-specific marker genes display cell-type-specific patterns of chromatin accessibility. Our data suggest that the differential chromatin accessibility is a critical mechanism to regulate gene activity at the cell-type level.