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.

DNA restriction fragments ranging from 79 to 789 base pairs in length have been characterized by transient electric birefringence (TEB) measurements at various temperatures between 4 and 43 degrees C. The DNA fragments do not contain runs of four or more adenine residues in a row and migrate with normal electrophoretic mobilities in polyacrylamide gels, indicating that they are not intrinsically curved or bent. The low ionic strength buffers used for the measurements contained 1 mM Tris Cl, pH 8.0, EDTA, and variable concentrations of Na(+) or Mg(2+) ions. The rotational relaxation times were obtained by fitting the TEB field-free decay signals with a nonlinear least-squared fitting program; the decay of the birefringence was monoexponential for fragments < or = 241 base pair (bp) in length and multiexponential for larger fragments. The terminal relaxation times, characteristic of the end-over-end rotation of the DNA molecules, were then used to determine the persistence length (p) and hydrodynamic radius (r) of DNA as a function of temperature and ionic strength, using several different hydrodynamic models. The specific values obtained for p and r are model dependent. The wormlike chain model of P. J. Hagerman and B. H. Zimm (Biopolymers 1981, Vol. 20, pp. 1481-1502) combined with the revised Broersma equation (J. Newman et al., Journal of Mol Biol 1997, Vol. 116, pp. 593-606) appears to be the most suitable for describing the flexibility of DNA in low ionic strength solutions. The values of p and r obtained from the global least squares fitting of this equation are independent of DNA length, and the deviations of the individual values from the average are reasonably small. The consensus r value calculated for DNA in various low ionic strength solutions containing 1 mM Tris buffer is 14.7 +/- 0.4 A at 20 degrees C. The consensus p values decrease from 814 approximately 564 A in solutions containing 1 mM Tris buffer plus 0.2-1 mM NaCl and decrease still further to 440 A in solutions containing 0.2 mM Mg(2+) ions. The persistence length exhibits a shallow maximum at 20 degrees C and decreases slowly upon either increasing or decreasing the temperature, regardless of the model used to fit the data. By contrast, the consensus values of the hydrodynamic radius are independent of temperature. The calculated persistence lengths and hydrodynamic radii are compared with other data in the literature.

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.

BackgroundSorghum 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.ResultsQuantitative 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.ConclusionSorghum 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.

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.

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 19th 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.

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.

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.