JGI Plant Gene Atlas: an updateable transcriptome resource to improve functional gene descriptions across the plant kingdom.

Gene functional descriptions offer a crucial line of evidence for candidate genes underlying trait variation. Conversely, plant responses to environmental cues represent important resources to decipher gene function and subsequently provide molecular targets for plant improvement through gene editing. However, biological roles of large proportions of genes across the plant phylogeny are poorly annotated. Here we describe the Joint Genome Institute (JGI) Plant Gene Atlas, an updateable data resource consisting of transcript abundance assays spanning 18 diverse species. To integrate across these diverse genotypes, we analyzed expression profiles, built gene clusters that exhibited tissue/condition specific expression, and tested for transcriptional response to environmental queues. We discovered extensive phylogenetically constrained and condition-specific expression profiles for genes without any previously documented functional annotation. Such conserved expression patterns and tightly co-expressed gene clusters let us assign expression derived additional biological information to 64 495 genes with otherwise unknown functions. The ever-expanding Gene Atlas resource is available at JGI Plant Gene Atlas (https://plantgeneatlas.jgi.doe.gov) and Phytozome (https://phytozome.jgi.doe.gov/), providing bulk access to data and user-specified queries of gene sets. Combined, these web interfaces let users access differentially expressed genes, track orthologs across the Gene Atlas plants, graphically represent co-expressed genes, and visualize gene ontology and pathway enrichments.

Bioenergy sorghum stem growth regulation: intercalary meristem localization, development, and gene regulatory network analysis.

Bioenergy sorghum is a highly productive drought tolerant C4 grass that accumulates 80% of its harvestable biomass in approximately 4 m length stems. Stem internode growth is regulated by development, shading, and hormones that modulate cell proliferation in intercalary meristems (IMs). In this study, sorghum stem IMs were localized above the pulvinus at the base of elongating internodes using magnetic resonance imaging, microscopy, and transcriptome analysis. A change in cell morphology/organization occurred at the junction between the pulvinus and internode where LATERAL ORGAN BOUNDARIES (SbLOB), a boundary layer gene, was expressed. Inactivation of an AGCVIII kinase in DDYM (dw2) resulted in decreased SbLOB expression, disrupted IM localization, and reduced internode cell proliferation. Transcriptome analysis identified approximately 1000 genes involved in cell proliferation, hormone signaling, and other functions selectively upregulated in the IM compared with a non-meristematic stem tissue. This cohort of genes is expressed in apical dome stem tissues before localization of the IM at the base of elongating internodes. Gene regulatory network analysis identified connections between genes involved in hormone signaling and cell proliferation. The results indicate that gibberellic acid induces accumulation of growth regulatory factors (GRFs) known to interact with ANGUSTIFOLIA (SbAN3), a master regulator of cell proliferation. GRF:AN3 was predicted to induce SbARF3/ETT expression and regulate SbAN3 expression in an auxin-dependent manner. GRFs and ARFs regulate genes involved in cytokinin and brassinosteroid signaling and cell proliferation. The results provide a molecular framework for understanding how hormone signaling regulates the expression of genes involved in cell proliferation in the stem IM.

Ruggedized, field-ready snapshot light-guide-based imaging spectrometer for environmental and remote sensing applications.

A field-ready, fiber-based high spatial sampling snapshot imaging spectrometer was developed for applications such as environmental monitoring and smart farming. The system achieves video rate frame transfer and exposure times down to a few hundred microseconds in typical daylight conditions with ∼63,000 spatial points and 32 spectral channels across the 470nm to 700nm wavelength range. We designed portable, ruggedized opto-mechanics to allow for imaging from an airborne platform. To ensure successful data collection prior to flight, imaging speed and signal-to-noise ratio was characterized for imaging a variety of land covers from the air. The system was validated by performing a series of observations including: Liriope Muscari plants under a range of water-stress conditions in a controlled laboratory experiment and field observations of sorghum plants in a variety of soil conditions. Finally, we collected data from a series of engineering flights and present reassembled images and spectral sampling of rural and urban landscapes.

Low-field magnetic resonance imaging of roots in intact clayey and silty soils.

The development of a robust method to non-invasively visualize root morphology in natural soils has been hampered by the opaque, physical, and structural properties of soils. In this work we describe a novel technology, low field magnetic resonance imaging (LF-MRI), for imaging energy sorghum (Sorghum bicolor (L.) Moench) root morphology and architecture in intact soils. The use of magnetic fields much weaker than those used with traditional MRI experiments reduces the distortion due to magnetic material naturally present in agricultural soils. A laboratory based LF-MRI operating at 47 mT magnetic field strength was evaluated using two sets of soil cores: 1) soil/root cores of Weswood silt loam (Udifluventic Haplustept) and a Belk clay (Entic Hapluderts) from a conventionally tilled field, and 2) soil/root cores from rhizotrons filled with either a Houston Black (Udic Haplusterts) clay or a sandy loam purchased from a turf company. The maximum soil water nuclear magnetic resonance (NMR) relaxation time T2 (4 ms) and the typical root water relaxation time T2 (100 ms) are far enough apart to provide a unique contrast mechanism such that the soil water signal has decayed to the point of no longer being detectable during the data collection time period. 2-D MRI projection images were produced of roots with a diameter range of 1.5-2.0 mm using an image acquisition time of 15 min with a pixel resolution of 1.74 mm in four soil types. Additionally, we demonstrate the use of a data-driven machine learning reconstruction approach, Automated Transform by Manifold Approximation (AUTOMAP) to reconstruct raw data and improve the quality of the final images. The application of AUTOMAP showed a SNR (Signal to Noise Ratio) improvement of two fold on average. The use of low field MRI presented here demonstrates the possibility of applying low field MRI through intact soils to root phenotyping and agronomy to aid in understanding of root morphology and the spatial arrangement of roots in situ.

The Sorghum bicolor reference genome: improved assembly, gene annotations, a transcriptome atlas, and signatures of genome organization.

Sorghum bicolor is a drought tolerant C4 grass used for the production of grain, forage, sugar, and lignocellulosic biomass and a genetic model for C4 grasses due to its relatively small genome (approximately 800 Mbp), diploid genetics, diverse germplasm, and colinearity with other C4 grass genomes. In this study, deep sequencing, genetic linkage analysis, and transcriptome data were used to produce and annotate a high-quality reference genome sequence. Reference genome sequence order was improved, 29.6 Mbp of additional sequence was incorporated, the number of genes annotated increased 24% to 34 211, average gene length and N50 increased, and error frequency was reduced 10-fold to 1 per 100 kbp. Subtelomeric repeats with characteristics of Tandem Repeats in Miniature (TRIM) elements were identified at the termini of most chromosomes. Nucleosome occupancy predictions identified nucleosomes positioned immediately downstream of transcription start sites and at different densities across chromosomes. Alignment of more than 50 resequenced genomes from diverse sorghum genotypes to the reference genome identified approximately 7.4 M single nucleotide polymorphisms (SNPs) and 1.9 M indels. Large-scale variant features in euchromatin were identified with periodicities of approximately 25 kbp. A transcriptome atlas of gene expression was constructed from 47 RNA-seq profiles of growing and developed tissues of the major plant organs (roots, leaves, stems, panicles, and seed) collected during the juvenile, vegetative and reproductive phases. Analysis of the transcriptome data indicated that tissue type and protein kinase expression had large influences on transcriptional profile clustering. The updated assembly, annotation, and transcriptome data represent a resource for C4 grass research and crop improvement.

Sorghum Dw2 Encodes a Protein Kinase Regulator of Stem Internode Length.

Sorghum is an important C4 grass crop grown for grain, forage, sugar, and bioenergy production. While tall, late flowering landraces are commonly grown in Africa, short early flowering varieties were selected in US grain sorghum breeding programs to reduce lodging and to facilitate machine harvesting. Four loci have been identified that affect stem length (Dw1-Dw4). Subsequent research showed that Dw3 encodes an ABCB1 auxin transporter and Dw1 encodes a highly conserved protein involved in the regulation of cell proliferation. In this study, Dw2 was identified by fine-mapping and further confirmed by sequencing the Dw2 alleles in Dwarf Yellow Milo and Double Dwarf Yellow Milo, the progenitor genotypes where the recessive allele of dw2 originated. The Dw2 locus was determined to correspond to Sobic.006G067700, a gene that encodes a protein kinase that is homologous to KIPK, a member of the AGCVIII subgroup of the AGC protein kinase family in Arabidopsis.

CONSTANS is a photoperiod regulated activator of flowering in sorghum.

BACKGROUND: Sorghum genotypes used for grain production in temperate regions are photoperiod insensitive and flower early avoiding adverse environments during the reproductive phase. In contrast, energy sorghum hybrids are highly photoperiod sensitive with extended vegetative phases in long days, resulting in enhanced biomass accumulation. SbPRR37 and SbGHD7 contribute to photoperiod sensitivity in sorghum by repressing expression of SbEHD1 and FT-like genes, thereby delaying flowering in long days with minimal influence in short days (PNAS_108:16469-16474, 2011; Plant Genome_in press, 2014). The GIGANTEA (GI)-CONSTANS (CO)-FLOWERING LOCUS T (FT) pathway regulates flowering time in Arabidopsis and the grasses (J Exp Bot_62:2453-2463, 2011). In long day flowering plants, such as Arabidopsis and barley, CONSTANS activates FT expression and flowering in long days. In rice, a short day flowering plant, Hd1, the ortholog of CONSTANS, activates flowering in short days and represses flowering in long days. RESULTS: Quantitative trait loci (QTL) that modify flowering time in sorghum were identified by screening Recombinant Inbred Lines (RILs) derived from BTx642 and Tx7000 in long days, short days, and under field conditions. Analysis of the flowering time QTL on SBI-10 revealed that BTx642 encodes a recessive CONSTANS allele containing a His106Tyr substitution in B-box 2 known to inactivate CONSTANS in Arabidopsis thaliana. Genetic analysis characterized sorghum CONSTANS as a floral activator that promotes flowering by inducing the expression of EARLY HEADING DATE 1 (SbEHD1) and sorghum orthologs of the maize FT genes ZCN8 (SbCN8) and ZCN12 (SbCN12). The floral repressor PSEUDORESPONSE REGULATOR PROTEIN 37 (PRR37) inhibits sorghum CONSTANS activity and flowering in long days. CONCLUSION: Sorghum CONSTANS is an activator of flowering that is repressed post-transcriptionally in long days by the floral inhibitor PRR37, contributing to photoperiod sensitive flowering in Sorghum bicolor, a short day plant.

Extensive variation in the density and distribution of DNA polymorphism in sorghum genomes.

Sorghum genotypes currently used for grain production in the United States were developed from African landraces that were imported starting in the mid-to-late 19(th) century. Farmers and plant breeders selected genotypes for grain production with reduced plant height, early flowering, increased grain yield, adaptation to drought, and improved resistance to lodging, diseases and pests. DNA polymorphisms that distinguish three historically important grain sorghum genotypes, BTx623, BTx642 and Tx7000, were characterized by genome sequencing, genotyping by sequencing, genetic mapping, and pedigree-based haplotype analysis. The distribution and density of DNA polymorphisms in the sequenced genomes varied widely, in part because the lines were derived through breeding and selection from diverse Kafir, Durra, and Caudatum race accessions. Genomic DNA spanning dw1 (SBI-09) and dw3 (SBI-07) had identical haplotypes due to selection for reduced height. Lower SNP density in genes located in pericentromeric regions compared with genes located in euchromatic regions is consistent with background selection in these regions of low recombination. SNP density was higher in euchromatic DNA and varied >100-fold in contiguous intervals that spanned up to 300 Kbp. The localized variation in DNA polymorphism density occurred throughout euchromatic regions where recombination is elevated, however, polymorphism density was not correlated with gene density or DNA methylation. Overall, sorghum chromosomes contain distal euchromatic regions characterized by extensive, localized variation in DNA polymorphism density, and large pericentromeric regions of low gene density, diversity, and recombination.

Sorghum bicolor's transcriptome response to dehydration, high salinity and ABA.

Genome wide changes in gene expression were monitored in the drought tolerant C4 cereal Sorghum bicolor, following exposure of seedlings to high salinity (150 mM NaCl), osmotic stress (20% polyethylene glycol) or abscisic acid (125 microM ABA). A sorghum cDNA microarray providing data on 12,982 unique gene clusters was used to examine gene expression in roots and shoots at 3- and 27-h post-treatment. Expression of approximately 2200 genes, including 174 genes with currently unknown functions, of which a subset appear unique to monocots and/or sorghum, was altered in response to dehydration, high salinity or ABA. The modulated sorghum genes had homology to proteins involved in regulation, growth, transport, membrane/protein turnover/repair, metabolism, dehydration protection, reactive oxygen scavenging, and plant defense. Real-time PCR was used to quantify changes in relative mRNA abundance for 333 genes that responded to ABA, NaCl or osmotic stress. Osmotic stress inducible sorghum genes identified for the first time included a beta-expansin expressed in shoots, actin depolymerization factor, inositol-3-phosphate synthase, a non-C4 NADP-malic enzyme, oleosin, and three genes homologous to 9-cis-epoxycarotenoid dioxygenase that may be involved in ABA biosynthesis. Analysis of response profiles demonstrated the existence of a complex gene regulatory network that differentially modulates gene expression in a tissue- and kinetic-specific manner in response to ABA, high salinity and water deficit. Modulation of genes involved in signal transduction, chromatin structure, transcription, translation and RNA metabolism contributes to sorghum's overlapping but nonetheless distinct responses to ABA, high salinity, and osmotic stress. Overall, this study provides a foundation of information on sorghum's osmotic stress responsive gene complement that will accelerate follow up biochemical, QTL and comparative studies.

Intrinsic curvature in the VP1 gene of SV40: comparison of solution and gel results.

DNA restriction fragments that are stably curved are usually identified by polyacrylamide gel electrophoresis because curved fragments migrate more slowly than normal fragments containing the same number of basepairs. In free solution, curved DNA molecules can be identified by transient electric birefringence (TEB) because they exhibit rotational relaxation times that are faster than those of normal fragments of the same size. In this article, the results observed in free solution and in polyacrylamide gels are compared for a highly curved 199-basepair (bp) restriction fragment taken from the VP1 gene in Simian Virus 40 (SV40) and various sequence mutants and insertion derivatives. The TEB method of overlapping fragments was used to show that the 199-bp fragment has an apparent bend angle of 46 +/- 2 degrees centered at sequence position 1922 +/- 2 bp. Four unphased A- and T-tracts and a mixed A3T4-tract occur within a span of approximately 60 bp surrounding the apparent bend center; for brevity, this 60-bp sequence element is called a curvature module. Modifying any of the A- or T-tracts in the curvature module by site-directed mutagenesis decreases the curvature of the fragment; replacing all five A- and T-tracts by random-sequence DNA causes the 199-bp mutant to adopt a normal conformation, with normal electrophoretic mobilities and birefringence relaxation times. Hence, stable curvature in this region of the VP1 gene is due to the five unphased A- and T- tracts surrounding the apparent bend center. Discordant solution and gel results are observed when long inverted repeats are inserted within the curvature module. These insertion derivatives migrate anomalously slowly in polyacrylamide gels but have normal, highly flexible conformations in free solution. Discordant solution and gel results are not observed if the insert does not contain a long inverted repeat or if the long inverted repeat is added to the 199-bp fragment outside the curvature module. The results suggest that long inverted repeats can form hairpins or cruciforms when they are located within a region of the helix backbone that is intrinsically curved, leading to large mobility anomalies in polyacrylamide gels. Hairpin/cruciform formation is not observed in free solution, presumably because of rapid conformational exchange. Hence, DNA restriction fragments that migrate anomalously slowly in polyacrylamide gels are not necessarily stably curved in free solution.

Regulation of pyrimidine metabolism in plants.

Pyrimidine nucleotides represent one of the most fundamental of cellular components. They are the building blocks for the direct synthesis of DNA and RNA that function in information storage and retrieval within the cell, but they also participate in the metabolism of a large number of other cellular components from sugar interconversion to cellular polysaccharides to glycoproteins and phospholipids. Thus, the metabolism of pyrimidine nucleotides and their intracellular pool sizes influence vast areas of normal cellular metabolism. The first pyrimidine, UMP, is synthesized by a de novo pathway that appears to be mechanistically invariant in all organisms. UMP is then further modified to form other pyrimidines. Breakdown of deoxyribo- and ribonucleic acids, the main sink for pyrimidine nucleotides, allows pyrimidines to be reutilized for resynthesis of these important cellular components. Pyrimidines are salvaged by converting the modified components into the free base, uracil for reutilization. Finally, pyrimidines are degraded into simple cellular metabolites permitting reutilization of nitrogen and carbon from pyrimidine ring systems into cellular metabolic pools. The regulation of pyrimidine metabolism is tightly controlled in plants. Additionally, plants produce toxic secondary metabolites derived from pyrimidines for use as defense compounds.

Analysis of DNA bending by transient electric birefringence.

Transient electric birefringence has been used to analyze DNA bending in six restriction fragments containing 171, 174, 207, 263, 289, and 471 bp in three different low ionic strength buffers. The target fragments contain sequences corresponding to the apparent bend centers in pUC19 and Litmus 28, previously identified by the circular permutation assay (Strutz, K.; Stellwagen, N. C. Electrophoresis 1996, 17, 989-995). The target fragments migrate anomalously slowly in polyacrylamide gels and exhibit birefringence relaxation times that are shorter than those of restriction fragments of the same size, taken from nonbent regions of the same plasmids. Apparent bend angles ranging from 30 degrees to 41 degrees were calculated for the target fragments by tau-ratio method. The bend angles of four of the target fragments were independent of temperature from 4 degrees C to 20 degrees C, but decreased when the temperature was increased to 37 degrees C. The bend angles of the other two target fragments were independent of temperature over the entire range examined, 4 degrees -37 degrees C. Hence, the thermal stability of sequence-dependent bends in random-sequence DNA is variable. The bend angles of five of the six target fragments were independent of the presence or absence of Mg2+ ions in the solution, indicating most of the target fragments were stably bent or curved, rather than anisometrically flexible. Restriction fragments containing 219 and 224 bp, with sequences somewhat offset from the sequence of the 207 bp fragment, were also studied. Comparison of the tau-ratios of these overlapping fragments allowed both the bend angle and bend position to be independently determined. These methods should be useful for analyzing sequence-dependent bending in other random-sequence DNAs.

Analysis of the intrinsic bend in the M13 origin of replication by atomic force microscopy.

Atomic force microscopy (AFM) has been used to image a 471-bp bent DNA restriction fragment derived from the M13 origin of replication in plasmid LITMUS 28, and a 476-bp normal, unbent fragment from plasmid pUC19. The most probable angle of curvature of the 471-bp DNA fragment is 40-50 degrees, in reasonably good agreement with the bend angle determined by transient electric birefringence, 38 degrees +/- 7 degrees. The normal 476-bp DNA fragment exhibited a Gaussian distribution of bend angles centered at 0 degrees, indicating that this fragment does not contain an intrinsic bend. The persistence length, P, was estimated to be 60 +/- 8 and 62 +/- 8 nm for the 471- and 476-bp fragments, respectively, from the observed mean-square end-to-end distances in the AFM images. Since the P-values of the normal and bent fragments are close to each other, the overall flexibility of DNA fragments of this size is only marginally affected by the presence of a stable bend. The close agreement of AFM and transient electric birefringence results validates the suitability of both methods for characterizing DNA bending and flexibility.