MOROKOSHI: transcriptome database in Sorghum bicolor

Evolution and expression analysis of the sorghum ubiquitin-conjugating enzyme family

Ubiquitin-conjugating enzymes (UBCs), which catalyse the transfer of ubiquitin to substrate or E3 ligases, are key enzymes in ubiquitination modifications of target proteins. Current knowledge regarding the sorghum (Sorghum bicolor (L.) Moench) ubiquitin-conjugating enzyme (SbUBC) family remains very limited. We identified 53 UBC-encoding genes in the sorghum genome and divided these into 18 groups according to their phylogenetic relationship with Arabidopsis thaliana (L.) Heynh., which was further supported by conserved motif and gene structure analyses. Different expression levels under a variety of abiotic stresses suggested that these might participate in distinct signalling pathways and that they underwent functional divergence during evolution. Furthermore, several SbUBC genes responded to single treatments, and individual SbUBC genes responded to multiple treatments, suggesting that sorghum UBCs may mediate crosstalk among different signalling pathways. Overall, the results provide valuable information for better understanding the classification and putative functions of sorghum UBC-encoding genes.

Identification of O-acetylserine(thiol)lyase (OASTL) genes in sorghum (Sorghum bicolor) and gene expression analysis under cadmium stress

Cysteine (Cys) is the first identified molecule in plant metabolism which includes both sulfur and nitrogen. It can be synthesized in three cellular compartments, containing chloroplast, cytoplasm and mitochondrion. The final step of cysteine biosynthesis is catalyzed by the O-acetylserine(thiol)lyase enzyme (OASTL, E.C. 4.2.99). In the present study, seven members of the OASTL gene family in the sorghum (Sorghum bicolor) genome were identified at a genome-wide scale and comparative bioinformatics analyses were performed between sorghum and Arabidopsis OASTLs. In all OASTL proteins, a pyridoxal-phosphate dependent domain structure (PALP, PF00291) was identified. The gene ontology annotations also revealed that all sorghum OASTL genes have KOG1252 (Cystathionine beta-synthase and related enzyme) and K01738 (cysteine synthase A) activities. In promotor sequences of OASTL genes, diverse cis-acting elements were found, including hormone and light responsiveness, abiotic stress responsiveness, and tissue-specific ones (meristem and endosperm). Sorghum OASTL genes demonstrated medium or high level expressions in anatomical parts and developmental stages based on the digital expression data. Expression of OASTL genes were also analyzed under cadmium (Cd) stress in sorghum by Real Time-quantitative PCR (RT-qPCR). The results exclusively showed that OASTL A1-2 gene was 1.12 fold up-regulated in roots, whereas cysteine synthase 26 was 2.25 fold down-regulated in leaves. The predicted 3D structure of OASTLs indicated some structural diversities as well as variations in the secondary structures.

Overexpression of the Sorghum bicolor SbCCoAOMT alters cell wall associated hydroxycinnamoyl groups

Sorghum (Sorghum bicolor) is a drought tolerant crop, which is being developed as a bioenergy feedstock. The monolignol biosynthesis pathway is a major focus for altering the abundance and composition of lignin. Caffeoyl coenzyme-A O-methyltransferase (CCoAOMT) is an S-adenosyl methionine (SAM)-dependent O-methyltransferase that methylates caffeoyl-CoA to generate feruloyl-CoA, an intermediate required for the biosynthesis of both G- and S-lignin. SbCCoAOMT was overexpressed to assess the impact of increasing the amount of this enzyme on biomass composition. SbCCoAOMT overexpression increased both soluble and cell wall-bound (esterified) ferulic and sinapic acids, however lignin concentration and its composition (S/G ratio) remained unaffected. This increased deposition of hydroxycinnamic acids in these lines led to an increase in total energy content of the stover. In stalk and leaf midribs, the increased histochemical staining and autofluorescence in the cell walls of the SbCCoAOMT overexpression lines also indicate increased phenolic deposition within cell walls, which is consistent with the chemical analyses of soluble and wall-bound hydroxycinnamic acids. The growth and development of overexpression lines were similar to wild-type plants. Likewise, RNA-seq and metabolite profiling showed that global gene expression and metabolite levels in overexpression lines were also relatively similar to wild-type plants. Our results demonstrate that SbCCoAOMT overexpression significantly altered cell wall composition through increases in cell wall associated hydroxycinnamic acids without altering lignin concentration or affecting plant growth and development.

DREB2 (dehydration-responsive element-binding protein 2) type transcription factor in sorghum (Sorghum bicolor): genome-wide identification, characterization and expression profiles under cadmium and salt stresses

Biotic and abiotic stresses negatively affect fitness, biomass production, and crop yield in plants. The dehydration-responsive element-binding proteins (DREB) are important transcription factors (TFs), and are induced by abiotic and biotic stresses. In this study, genome-wide identification, in silico sequence, and phylogenetic analyses and expression analyses of DREB2 genes under cadmium (Cd) and salt (NaCl) stresses in sorghum (Sorghum bicolor, Sb) were performed. Six putative SbDREB2 genes were identified in sorghum genome and all contained AP2 domain (PF00847). Nucleotide diversities in SbDREB2 genes were calculated as π: 0.53 and θ: 0.39, respectively. While exon numbers of them were either one or two, length of SbDREB2 proteins ranged from 238 to 388 amino acid residues. Fifty-six cis-acting regulatory elements, which are tissue specific, light, hormone, and stress responsive, were identified in the promotor regions of SbDREB2 genes. Analyses on digital expression data indicated that SbDREB2A and SbDREB2B are more expressed genes than other SbDREB genes in sorghum. Under Cd and NaCl stresses, expressions of SbDREB2 genes were induced at different levels. All SbDREB2 genes in root were up-regulated under salt stress. In case of Cd stress, SbDREB2D gene was particularly up-regulated in leaves and roots. Co-expression analyses revealed four of TFs in co-expression network, indicating that they have roles in transcriptional cascade. Furthermore, five miRNA target regions were identified for four SbDREB2 genes, indicating their roles in post-transcriptional regulation. The predicted 3D structure of SbDREB2 proteins showed some structural divergences and structure overlap between rice and sorghum varied at between 26.58 and 50%. Finally, obtained data could be used in breeding of stress-tolerant plants, particularly genetically engineered DREB2 expressing plants. Findings in this study would also contribute to the understanding of DREB2 genes in plants, especially in sorghum.

Cross-species multiple environmental stress responses: An integrated approach to identify candidate genes for multiple stress tolerance in sorghum (Sorghum bicolor (L.) Moench) and related model species

BACKGROUND

Crop response to the changing climate and unpredictable effects of global warming with adverse conditions such as drought stress has brought concerns about food security to the fore; crop yield loss is a major cause of concern in this regard. Identification of genes with multiple responses across environmental stresses is the genetic foundation that leads to crop adaptation to environmental perturbations.

METHODS

In this paper, we introduce an integrated approach to assess candidate genes for multiple stress responses across-species. The approach combines ontology based semantic data integration with expression profiling, comparative genomics, phylogenomics, functional gene enrichment and gene enrichment network analysis to identify genes associated with plant stress phenotypes. Five different ontologies, viz., Gene Ontology (GO), Trait Ontology (TO), Plant Ontology (PO), Growth Ontology (GRO) and Environment Ontology (EO) were used to semantically integrate drought related information.

RESULTS

Target genes linked to Quantitative Trait Loci (QTLs) controlling yield and stress tolerance in sorghum (Sorghum bicolor (L.) Moench) and closely related species were identified. Based on the enriched GO terms of the biological processes, 1116 sorghum genes with potential responses to 5 different stresses, such as drought (18%), salt (32%), cold (20%), heat (8%) and oxidative stress (25%) were identified to be over-expressed. Out of 169 sorghum drought responsive QTLs associated genes that were identified based on expression datasets, 56% were shown to have multiple stress responses. On the other hand, out of 168 additional genes that have been evaluated for orthologous pairs, 90% were conserved across species for drought tolerance. Over 50% of identified maize and rice genes were responsive to drought and salt stresses and were co-located within multifunctional QTLs. Among the total identified multi-stress responsive genes, 272 targets were shown to be co-localized within QTLs associated with different traits that are responsive to multiple stresses. Ontology mapping was used to validate the identified genes, while reconstruction of the phylogenetic tree was instrumental to infer the evolutionary relationship of the sorghum orthologs. The results also show specific genes responsible for various interrelated components of drought response mechanism such as drought tolerance, drought avoidance and drought escape.

CONCLUSIONS

We submit that this approach is novel and to our knowledge, has not been used previously in any other research; it enables us to perform cross-species queries for genes that are likely to be associated with multiple stress tolerance, as a means to identify novel targets for engineering stress resistance in sorghum and possibly, in other crop species.

Construction of Pará rubber tree genome and multi-transcriptome database accelerates rubber researches

BACKGROUND

Natural rubber is an economically important material. Currently the Pará rubber tree, Hevea brasiliensis is the main commercial source. Little is known about rubber biosynthesis at the molecular level. Next-generation sequencing (NGS) technologies brought draft genomes of three rubber cultivars and a variety of RNA sequencing (RNA-seq) data. However, no current genome or transcriptome databases (DB) are organized by gene.

RESULTS

A gene-oriented database is a valuable support for rubber research. Based on our original draft genome sequence of H. brasiliensis RRIM600, we constructed a rubber tree genome and transcriptome DB. Our DB provides genome information including gene functional annotations and multi-transcriptome data of RNA-seq, full-length cDNAs including PacBio Isoform sequencing (Iso-Seq), ESTs and genome wide transcription start sites (TSSs) derived from CAGE technology. Using our original and publically available RNA-seq data, we calculated co-expressed genes for identifying functionally related gene sets and/or genes regulated by the same transcription factor (TF). Users can access multi-transcriptome data through both a gene-oriented web page and a genome browser. For the gene searching system, we provide keyword search, sequence homology search and gene expression search; users can also select their expression threshold easily.

CONCLUSION

The rubber genome and transcriptome DB provides rubber tree genome sequence and multi-transcriptomics data. This DB is useful for comprehensive understanding of the rubber transcriptome. This will assist both industrial and academic researchers for rubber and economically important close relatives such as R. communis, M. esculenta and J. curcas. The Rubber Transcriptome DB release 2017.03 is accessible at http://matsui-lab.riken.jp/rubber/ .

Expression patterns of the native Shrunken-2 promoter in Sorghum bicolor visualised through use of the GFP reporter gene

The AGPase large subunit (shrunken-2) promoter was demonstrated to be active in the placentochalaza and endosperm of developing grain as well as the root tips in transgenic sorghum. The temporal and spatial expression patterns of the Sorghum bicolor Shrunken-2 (Sh2) promoter were evaluated using the green fluorescence protein reporter gene (gfp) in transgenic sorghum, within the context of upregulating starch biosynthesis in the developing grain. GFP fluorescence was analysed throughout development in various tissue types using confocal laser scanning microscopy techniques. Sh2 promoter activity was first detected in the placentochalaza region of the developing caryopsis and apoplasm adjacent to the nucellar epidermis at 7 days post anthesis (dpa) where fluorescence remained relatively constant until 17 dpa. Fluorescence in this region weakened by 20 dpa and disappeared by 25 dpa. Expression was also detected in the developing endosperm, but not until 12 dpa, continuing until 25 dpa. Whilst the endosperm expression was expected, the fluorescence detected in the placentochalaza was completely unexpected. Although transcript presence does not mean the resulting biochemistry is also present, these preliminary findings may suggest alternate spatial activity of ADP-glucose pyrophosphorylase prior to uptake by the developing grain. Sh2 promoter activity was also unexpectedly detected in the root tips at all developmental time points. Sh2 promoter activity was not detected in any reproductive floral tissue (both pre and post anthesis) or in pollen. Similarly, no expression was detected in leaf tissue at any stage.

Large-scale collection of full-length cDNA and transcriptome analysis in Hevea brasiliensis

Natural rubber has unique physical properties that cannot be replaced by products from other latex-producing plants or petrochemically produced synthetic rubbers. Rubber from Hevea brasiliensis is the main commercial source for this natural rubber that has a cis-polyisoprene configuration. For sustainable production of enough rubber to meet demand elucidation of the molecular mechanisms involved in the production of latex is vital. To this end, we firstly constructed rubber full-length cDNA libraries of RRIM 600 cultivar and sequenced around 20,000 clones by the Sanger method and over 15,000 contigs by Illumina sequencer. With these data, we updated around 5,500 gene structures and newly annotated around 9,500 transcription start sites. Second, to elucidate the rubber biosynthetic pathways and their transcriptional regulation, we carried out tissue- and cultivar-specific RNA-Seq analysis. By using our recently published genome sequence, we confirmed the expression patterns of the rubber biosynthetic genes. Our data suggest that the cytoplasmic mevalonate (MVA) pathway is the main route for isoprenoid biosynthesis in latex production. In addition to the well-studied polymerization factors, we suggest that rubber elongation factor 8 (REF8) is a candidate factor in cis-polyisoprene biosynthesis. We have also identified 39 transcription factors that may be key regulators in latex production. Expression profile analysis using two additional cultivars, RRIM 901 and PB 350, via an RNA-Seq approach revealed possible expression differences between a high latex-yielding cultivar and a disease-resistant cultivar.

Genome-Wide Association Study of Developing Leaves' Heat Tolerance during Vegetative Growth Stages in a Sorghum Association Panel

Heat stress reduces grain yield and quality worldwide. Enhancing heat tolerance of crops at all developmental stages is one of the essential strategies required for sustaining agricultural production especially as frequency of temperature extremes escalates in response to climate change. Although heat tolerance mechanisms have been studied extensively in model plant species, little is known about the genetic control underlying heat stress responses of crop plants at the vegetative stage under field conditions. To dissect the genetic basis of heat tolerance in sorghum [ (L.) Moench], we performed a genome-wide association study (GWAS) for traits responsive to heat stress at the vegetative stage in an association panel. Natural variation in leaf firing (LF) and leaf blotching (LB) were evaluated separately for 3 yr in experimental fields at three locations where sporadic heat waves occurred throughout the sorghum growing season. We identified nine single-nucleotide polymorphisms (SNPs) that were significantly associated with LF and five SNPs that were associated with LB. Candidate genes near the SNPs were investigated and 14 were directly linked to biological pathways involved in plant stress responses including heat stress response. The findings of this study provide new knowledge on the genetic control of leaf traits responsive to heat stress in sorghum, which could aid in elucidating the genetic and molecular mechanisms of vegetative stage heat tolerance in crops. The results also provide candidate markers for molecular breeding of enhanced heat tolerance in cereal and bioenergy crops.

Sweet sorghum as biofuel feedstock: recent advances and available resources

Sweet sorghum is a promising target for biofuel production. It is a C4 crop with low input requirements and accumulates high levels of sugars in its stalks. However, large-scale planting on marginal lands would require improved varieties with optimized biofuel-related traits and tolerance to biotic and abiotic stresses. Considering this, many studies have been carried out to generate genetic and genomic resources for sweet sorghum. In this review, we discuss various attributes of sweet sorghum that make it an ideal candidate for biofuel feedstock, and provide an overview of genetic diversity, tools, and resources available for engineering and/or marker-assisting breeding of sweet sorghum. Finally, the progress made so far, in identification of genes/quantitative trait loci (QTLs) important for agronomic traits and ongoing molecular breeding efforts to generate improved varieties, has been discussed.

Characterization and Regulation of Aquaporin Genes of Sorghum [Sorghum bicolor (L.) Moench] in Response to Waterlogging Stress

Waterlogging is a significant environmental constraint to crop production, and a better understanding of plant responses is critical for the improvement of crop tolerance to waterlogged soils. Aquaporins (AQPs) are a class of channel-forming proteins that play an important role in water transport in plants. This study aimed to examine the regulation of AQP genes under waterlogging stress and to characterize the genetic variability of AQP genes in sorghum (Sorghum bicolor). Transcriptional profiling of AQP genes in response to waterlogging stress in nodal root tips and nodal root basal regions of two tolerant and two sensitive sorghum genotypes at 18 and 96 h after waterlogging stress imposition revealed significant gene-specific pattern with regard to genotype, root tissue sample, and time point. For some tissue sample and time point combinations, PIP2-6, PIP2-7, TIP2-2, TIP4-4, and TIP5-1 expression was differentially regulated in tolerant compared to sensitive genotypes. The differential response of these AQP genes suggests that they may play a tissue specific role in mitigating waterlogging stress. Genetic analysis of sorghum revealed that AQP genes were clustered into the same four subfamilies as in maize (Zea mays) and rice (Oryza sativa) and that residues determining the AQP channel specificity were largely conserved across species. Single nucleotide polymorphism (SNP) data from 50 sorghum accessions were used to build an AQP gene-based phylogeny of the haplotypes. Phylogenetic analysis based on single nucleotide polymorphisms of sorghum AQP genes placed the tolerant and sensitive genotypes used for the expression study in distinct groups. Expression analyses suggested that selected AQPs may play a pivotal role in sorghum tolerance to water logging stress. Further experimentation is needed to verify their role and to leverage phylogenetic analyses and AQP expression data to improve waterlogging tolerance in sorghum.

Structural and Biochemical Characterization of Cinnamoyl-CoA Reductases

Cinnamoyl-coenzyme A reductase (CCR) catalyzes the reduction of hydroxycinnamoyl-coenzyme A (CoA) esters using NADPH to produce hydroxycinnamyl aldehyde precursors in lignin synthesis. The catalytic mechanism and substrate specificity of cinnamoyl-CoA reductases from sorghum (Sorghum bicolor), a strategic plant for bioenergy production, were deduced from crystal structures, site-directed mutagenesis, and kinetic and thermodynamic analyses. Although SbCCR1 displayed higher affinity for caffeoyl-CoA or p-coumaroyl-CoA than for feruloyl-CoA, the enzyme showed significantly higher activity for the latter substrate. Through molecular docking and comparisons between the crystal structures of the Vitis vinifera dihydroflavonol reductase and SbCCR1, residues threonine-154 and tyrosine-310 were pinpointed as being involved in binding CoA-conjugated phenylpropanoids. Threonine-154 of SbCCR1 and other CCRs likely confers strong substrate specificity for feruloyl-CoA over other cinnamoyl-CoA thioesters, and the T154Y mutation in SbCCR1 led to broader substrate specificity and faster turnover. Through data mining using our structural and biochemical information, four additional putative CCR genes were discovered from sorghum genomic data. One of these, SbCCR2, displayed greater activity toward p-coumaroyl-CoA than did SbCCR1, which could imply a role in the synthesis of defense-related lignin. Taken together, these findings provide knowledge about critical residues and substrate preference among CCRs and provide, to our knowledge, the first three-dimensional structure information for a CCR from a monocot species.

Genome-wide association analysis of seedling traits in diverse Sorghum germplasm under thermal stress

BACKGROUND

Climate variability due to fluctuation in temperature is a worldwide concern that imperils crop production. The need to understand how the germplasm variation in major crops can be utilized to aid in discovering and developing breeding lines that can withstand and adapt to temperature fluctuations is more necessary than ever. Here, we analyzed the genetic variation associated with responses to thermal stresses in a sorghum association panel (SAP) representing major races and working groups to identify single nucleotide polymorphisms (SNPs) that are associated with resilience to temperature stress in a major cereal crop.

RESULTS

The SAP exhibited extensive variation for seedling traits under cold and heat stress. Genome-wide analyses identified 30 SNPs that were strongly associated with traits measured at seedling stage under cold stress and tagged genes that act as regulators of anthocyanin expression and soluble carbohydrate metabolism. Meanwhile, 12 SNPs were significantly associated with seedling traits under heat stress and these SNPs tagged genes that function in sugar metabolism, and ion transport pathways. Evaluation of co-expression networks for genes near the significantly associated SNPs indicated complex gene interactions for cold and heat stresses in sorghum. We focused and validated the expression of four genes in the network of Sb06g025040, a basic-helix-loop-helix (bHLH) transcription factor that was proposed to be involved in purple color pigmentation of leaf, and observed that genes in this network were upregulated during cold stress in a moderately tolerant line as compared to the more sensitive line.

CONCLUSION

This study facilitated the tagging of genome regions associated with variation in seedling traits of sorghum under cold and heat stress. These findings show the potential of genotype information for development of temperature resilient sorghum cultivars and further characterization of genes and their networks responsible for adaptation to thermal stresses. Knowledge on the gene networks from this research can be extended to the other cereal crops to better understand the genetic basis of resilience to temperature fluctuations during plant developmental stages.

Plant Carbonic Anhydrases: Structures, Locations, Evolution, and Physiological Roles

Carbonic anhydrases (CAs) are zinc metalloenzymes that catalyze the interconversion of CO2 and HCO3- and are ubiquitous in nature. Higher plants contain three evolutionarily distinct CA families, αCAs, βCAs, and γCAs, where each family is represented by multiple isoforms in all species. Alternative splicing of CA transcripts appears common; consequently, the number of functional CA isoforms in a species may exceed the number of genes. CAs are expressed in numerous plant tissues and in different cellular locations. The most prevalent CAs are those in the chloroplast, cytosol, and mitochondria. This diversity in location is paralleled in the many physiological and biochemical roles that CAs play in plants. In this review, the number and types of CAs in C3, C4, and crassulacean acid metabolism (CAM) plants are considered, and the roles of the α and γCAs are briefly discussed. The remainder of the review focuses on plant βCAs and includes the identification of homologs between species using phylogenetic approaches, a consideration of the inter- and intracellular localization of the proteins, along with the evidence for alternative splice forms. Current understanding of βCA tissue-specific expression patterns and what controls them are reviewed, and the physiological roles for which βCAs have been implicated are presented.

Efficient Identification of Causal Mutations through Sequencing of Bulked F 2 from Two Allelic Bloomless Mutants of Sorghum bicolor

Sorghum (Sorghum bicolor Moench, L.) plant accumulates copious layers of epi-cuticular wax (EW) on its aerial surfaces, to a greater extent than most other crops. EW provides a vapor barrier that reduces water loss, and is therefore considered to be a major determinant of sorghum's drought tolerance. However, little is known about the genes responsible for wax accumulation in sorghum. We isolated two allelic mutants, bloomless40-1 (bm40-1) and bm40-2, from a mutant library constructed from ethyl methane sulfonate (EMS) treated seeds of an inbred, BTx623. Both mutants were nearly devoid of the EW layer. Each bm mutant was crossed to the un-mutated BTx623 to generated F₂ populations that segregated for the bm phenotype. Genomic DNA from 20 bm F₂ plants from each population was bulked for whole genome sequencing. A single gene, Sobic.001G228100, encoding a GDSL-like lipase/acylhydrolase, had unique homozygous mutations in each bulked F₂ population. Mutant bm40-1 harbored a missense mutation in the gene, whereas bm40-2 had a splice donor site mutation. Our findings thus provide strong evidence that mutation in this GDSL-like lipase gene causes the bm phenotype, and further demonstrate that this approach of sequencing two independent allelic mutant populations is an efficient method for identifying causal mutations. Combined with allelic mutants, MutMap provides powerful method to identify all causal genes for the large collection of bm mutants in sorghum, which will provide insight into how sorghum plants accumulate such abundant EW on their aerial surface. This knowledge may facilitate the development of tools for engineering drought-tolerant crops with reduced water loss.

Comparative phylogenomic analysis provides insights into TCP gene functions in Sorghum

Sorghum is a highly efficient C4 crop with potential to mitigate challenges associated with food, feed and fuel. TCP proteins are of particular interest for crop improvement programs due to their well-demonstrated roles in crop domestication and shaping plant architecture thereby, affecting agronomic traits. We identified 20 TCP genes from Sorghum. Except SbTCP8, all are either intronless or contain introns in the untranslated regions. Comparative phylogenetic analysis of Arabidopsis, rice, Brachypodium and Sorghum TCP proteins revealed two distinct classes categorized into ten sub-clades. Sub-clade F is dicot-specific, whereas A2, G1 and I1 groups only contained genes from grasses. Sub-clade B was missing in Sorghum, whereas group A1 was missing in rice indicating species-specific divergence of TCP proteins. TCP proteins of Sorghum are enriched in disorder promoting residues with class I containing higher percent disorder than class II proteins. Seven pairs of paralogous TCP genes were identified from Sorghum, five of which seem to predate Rice-Sorghum divergence. All of them have diverged in their expression. Based on the expression and orthology analysis, five Sorghum genes have been shortlisted for further investigation for their roles in regulating plant morphology. Whereas, three genes have been identified as candidates for engineering abiotic stress tolerance.

Domestication and the storage starch biosynthesis pathway: signatures of selection from a whole sorghum genome sequencing strategy

Next-generation sequencing of complete genomes has given researchers unprecedented levels of information to study the multifaceted evolutionary changes that have shaped elite plant germplasm. In conjunction with population genetic analytical techniques and detailed online databases, we can more accurately capture the effects of domestication on entire biological pathways of agronomic importance. In this study, we explore the genetic diversity and signatures of selection in all predicted gene models of the storage starch synthesis pathway of Sorghum bicolor, utilizing a diversity panel containing lines categorized as either 'Landraces' or 'Wild and Weedy' genotypes. Amongst a total of 114 genes involved in starch synthesis, 71 had at least a single signal of purifying selection and 62 a signal of balancing selection and others a mix of both. This included key genes such as STARCH PHOSPHORYLASE 2 (SbPHO2, under balancing selection), PULLULANASE (SbPUL, under balancing selection) and ADP-glucose pyrophosphorylases (SHRUNKEN2, SbSH2 under purifying selection). Effectively, many genes within the primary starch synthesis pathway had a clear reduction in nucleotide diversity between the Landraces and wild and weedy lines indicating that the ancestral effects of domestication are still clearly identifiable. There was evidence of the positional rate variation within the well-characterized primary starch synthesis pathway of sorghum, particularly in the Landraces, whereby low evolutionary rates upstream and high rates downstream in the metabolic pathway were expected. This observation did not extend to the wild and weedy lines or the minor starch synthesis pathways.

The biosynthetic gene cluster for the cyanogenic glucoside dhurrin in Sorghum bicolor contains its co-expressed vacuolar MATE transporter

Genomic gene clusters for the biosynthesis of chemical defence compounds are increasingly identified in plant genomes. We previously reported the independent evolution of biosynthetic gene clusters for cyanogenic glucoside biosynthesis in three plant lineages. Here we report that the gene cluster for the cyanogenic glucoside dhurrin in Sorghum bicolor additionally contains a gene, SbMATE2, encoding a transporter of the multidrug and toxic compound extrusion (MATE) family, which is co-expressed with the biosynthetic genes. The predicted localisation of SbMATE2 to the vacuolar membrane was demonstrated experimentally by transient expression of a SbMATE2-YFP fusion protein and confocal microscopy. Transport studies in Xenopus laevis oocytes demonstrate that SbMATE2 is able to transport dhurrin. In addition, SbMATE2 was able to transport non-endogenous cyanogenic glucosides, but not the anthocyanin cyanidin 3-O-glucoside or the glucosinolate indol-3-yl-methyl glucosinolate. The genomic co-localisation of a transporter gene with the biosynthetic genes producing the transported compound is discussed in relation to the role self-toxicity of chemical defence compounds may play in the formation of gene clusters.

Whole Genome Sequencing Reveals Potential New Targets for Improving Nitrogen Uptake and Utilization in Sorghum bicolor

Nitrogen (N) fertilizers are a major agricultural input where more than 100 million tons are supplied annually. Cereals are particularly inefficient at soil N uptake, where the unrecovered nitrogen causes serious environmental damage. Sorghum bicolor (sorghum) is an important cereal crop, particularly in resource-poor semi-arid regions, and is known to have a high NUE in comparison to other major cereals under limited N conditions. This study provides the first assessment of genetic diversity and signatures of selection across 230 fully sequenced genes putatively involved in the uptake and utilization of N from a diverse panel of sorghum lines. This comprehensive analysis reveals an overall reduction in diversity as a result of domestication and a total of 128 genes displaying signatures of purifying selection, thereby revealing possible gene targets to improve NUE in sorghum and cereals alike. A number of key genes appear to have been involved in selective sweeps, reducing their sequence diversity. The ammonium transporter (AMT) genes generally had low allelic diversity, whereas a substantial number of nitrate/peptide transporter 1 (NRT1/PTR) genes had higher nucleotide diversity in domesticated germplasm. Interestingly, members of the distinct race Guinea margaritiferum contained a number of unique alleles, and along with the wild sorghum species, represent a rich resource of new variation for plant improvement of NUE in sorghum.

The Structure and Catalytic Mechanism of Sorghum bicolor Caffeoyl-CoA O-Methyltransferase

Caffeoyl-coenzyme A 3-O-methyltransferase (CCoAOMT) is an S-adenosyl methionine (SAM)-dependent O-methyltransferase responsible for methylation of the meta-hydroxyl group of caffeoyl-coenzyme A (CoA) on the pathway to monolignols, with their ring methoxylation status characteristic of guaiacyl or syringyl units in lignin. In order to better understand the unique class of type 2 O-methyltransferases from monocots, we have characterized CCoAOMT from sorghum (Sorghum bicolor; SbCCoAOMT), including the SAM binary complex crystal structure and steady-state enzyme kinetics. Key amino acid residues were validated with site-directed mutagenesis. Isothermal titration calorimetry data indicated a sequential binding mechanism for SbCCoAOMT, wherein SAM binds prior to caffeoyl-CoA, and the enzyme showed allosteric behavior with respect to it. 5-Hydroxyferuloyl-CoA was not a substrate for SbCCoAOMT. We propose a catalytic mechanism in which lysine-180 acts as a catalytic base and deprotonates the reactive hydroxyl group of caffeoyl-CoA. This deprotonation is facilitated by the coordination of the reactive hydroxyl group by Ca(2+) in the active site, lowering the pKa of the 3'-OH group. Collectively, these data give a new perspective on the catalytic mechanism of CCoAOMTs and provide a basis for the functional diversity exhibited by type 2 plant OMTs that contain a unique insertion loop (residues 208-231) conferring affinity for phenylpropanoid-CoA thioesters. The structural model of SbCCoAOMT can serve as the basis for protein engineering approaches to enhance the nutritional, agronomic, and industrially relevant properties of sorghum.

Genome-Wide Association Studies of Grain Yield Components in Diverse Sorghum Germplasm

Grain yield and its primary determinants, grain number and weight, are important traits in cereal crops that have been well studied; however, the genetic basis of and interactions between these traits remain poorly understood. Characterization of grain yield per primary panicle (YPP), grain number per primary panicle (GNP), and 1000-grain weight (TGW) in sorghum [ (L.) Moench], a hardy C cereal with a genome size of ∼730 Mb, was implemented in a diversity panel containing 390 accessions. These accessions were genotyped to obtain 268,830 single-nucleotide polymorphisms (SNPs). Genome-wide association studies (GWAS) were performed to identify loci associated with each grain yield component and understand the genetic interactions between these traits. Genome-wide association studies identified associations across the genome with YPP, GNP, and TGW that were located within previously mapped sorghum QTL for panicle weight, grain yield, and seed size, respectively. There were no significant associations between GNP and TGW that were within 100 kb, much greater than the average linkage disequilibrium (LD) in sorghum. The identification of nonoverlapping loci for grain number and weight suggests these traits may be manipulated independently to increase the grain yield of sorghum. Following GWAS, genomic regions surrounding each associated SNP were mined for candidate genes. Previously published expression data indicated several TGW candidate genes, including an ethylene receptor homolog, were primarily expressed within developing seed tissues to support GWAS. Furthermore, maize ( L.) homologs of identified TGW candidates were differentially expressed within the seed between small- and large-kernel lines from a segregating maize population.

SorghumFDB: sorghum functional genomics database with multidimensional network analysis

Sorghum (Sorghum bicolor [L.] Moench) has excellent agronomic traits and biological properties, such as heat and drought-tolerance. It is a C₄ grass and potential bioenergy-producing plant, which makes it an important crop worldwide. With the sorghum genome sequence released, it is essential to establish a sorghum functional genomics data mining platform. We collected genomic data and some functional annotations to construct a sorghum functional genomics database (SorghumFDB). SorghumFDB integrated knowledge of sorghum gene family classifications (transcription regulators/factors, carbohydrate-active enzymes, protein kinases, ubiquitins, cytochrome P450, monolignol biosynthesis related enzymes, R-genes and organelle-genes), detailed gene annotations, miRNA and target gene information, orthologous pairs in the model plants Arabidopsis, rice and maize, gene loci conversions and a genome browser. We further constructed a dynamic network of multidimensional biological relationships, comprised of the co-expression data, protein–protein interactions and miRNA-target pairs. We took effective measures to combine the network, gene set enrichment and motif analyses to determine the key regulators that participate in related metabolic pathways, such as the lignin pathway, which is a major biological process in bioenergy-producing plants.

Database URL:http://structuralbiology.cau.edu.cn/sorghum/index.html.

Molecular Approaches to Understand Nutritional Potential of Coarse Cereals

Coarse grains are important group of crops that constitutes staple food for large population residing primarily in the arid and semi-arid regions of the world. Coarse grains are designated as nutri-cereals as they are rich in essential amino acids, minerals and vitamins. In spite of having several nutritional virtues in coarse grain as mentioned above, there is still scope for improvement in quality parameters such as cooking qualities, modulation of nutritional constituents and reduction or elimination of anti-nutritional factors. Besides its use in traditional cooking, coarse grains have been used mainly in the weaning food preparation and other malted food production. Improvement in quality parameters will certainly increase consumer's preference for coarse grains and increase their demand. The overall genetic gain in quality traits of economic importance in the cultivated varieties will enhance their industrial value and simultaneously increase income of farmers growing these varieties. The urgent step for improvement of quality traits in coarse grains requires a detailed understanding of molecular mechanisms responsible for varied level of different nutritional contents in different genotypes of these crops. In this review we have discussed the progresses made in understanding of coarse grain biology with various omics tool coupled with modern breeding approaches and the current status with regard to our effort towards dissecting traits related to improvement of quality and nutritional constituents of grains.

Mapping QTLs and association of differentially expressed gene transcripts for multiple agronomic traits under different nitrogen levels in sorghum

BACKGROUND

Sorghum is an important C4 crop which relies on applied Nitrogen fertilizers (N) for optimal yields, of which substantial amounts are lost into the atmosphere. Understanding the genetic variation of sorghum in response to limited nitrogen supply is important for elucidating the underlying genetic mechanisms of nitrogen utilization.

RESULTS

A bi-parental mapping population consisting of 131 recombinant inbred lines (RILs) was used to map quantitative trait loci (QTLs) influencing different agronomic traits evaluated under normal N (100 kg.ha(-1) fertilizer) and low N (0 kg.ha(-1) fertilizer) conditions. A linkage map spanning 1614 cM was developed using 642 polymorphic single nucleotide polymorphisms (SNPs) detected in the population using Genotyping-By-Sequencing (GBS) technology. Composite interval mapping detected a total of 38 QTLs for 11 agronomic traits tested under different nitrogen levels. The phenotypic variation explained by individual QTL ranged from 6.2 to 50.8%. Illumina RNA sequencing data generated on seedling root tissues revealed 726 differentially expressed gene (DEG) transcripts between parents, of which 108 were mapped close to the QTL regions.

CONCLUSIONS

Co-localized regions affecting multiple traits were detected on chromosomes 1, 5, 6, 7 and 9. These potentially pleiotropic regions were coincident with the genomic regions of cloned QTLs, including genes associated with flowering time, Ma3 on chromosome 1 and Ma1 on chromosome 6, gene associated with plant height, Dw2 on chromosome 6. In these regions, RNA sequencing data showed differential expression of transcripts related to nitrogen metabolism (Ferredoxin-nitrate reductase), glycolysis (Phosphofructo-2-kinase), seed storage proteins, plant hormone metabolism and membrane transport. The differentially expressed transcripts underlying the pleiotropic QTL regions could be potential targets for improving sorghum performance under limited N fertilizer through marker assisted selection.

Tonoplast Sugar Transporters (SbTSTs) putatively control sucrose accumulation in sweet sorghum stems

Carbohydrates are differentially partitioned in sweet versus grain sorghums. While the latter preferentially accumulate starch in the grain, the former primarily store large amounts of sucrose in the stem. Previous work determined that neither sucrose metabolizing enzymes nor changes in Sucrose transporter (SUT) gene expression accounted for the carbohydrate partitioning differences. Recently, 2 additional classes of sucrose transport proteins, Tonoplast Sugar Transporters (TSTs) and SWEETs, were identified; thus, we examined whether their expression tracked sucrose accumulation in sweet sorghum stems. We determined 2 TSTs were differentially expressed in sweet vs. grain sorghum stems, likely underlying the massive difference in sucrose accumulation. A model illustrating potential roles for different classes of sugar transport proteins in sorghum sugar partitioning is discussed.

The sorghum SWEET gene family: stem sucrose accumulation as revealed through transcriptome profiling

BACKGROUND

SWEET is a newly identified family of sugar transporters. Although SWEET transporters have been characterized by using Arabidopsis and rice, very little knowledge of sucrose accumulation in the stem region is available, as these model plants accumulate little sucrose in their stems. To elucidate the expression of key SWEET genes involved in sucrose accumulation of sorghum, we performed transcriptome profiling by RNA-seq, categorization using phylogenetic trees, analysis of chromosomal synteny, and comparison of amino acid sequences between SIL-05 (a sweet sorghum) and BTx623 (a grain sorghum).

RESULTS

We identified 23 SWEET genes in the sorghum genome. In the leaf, SbSWEET8-1 was highly expressed and was grouped in the same clade as AtSWEET11 and AtSWEET12 that play a role in the efflux of photosynthesized sucrose. The key genes in sucrose synthesis (SPS3) and that in another step of sugar transport (SbSUT1 and SbSUT2) were also highly expressed, suggesting that sucrose is newly synthesized and actively exported from the leaf. In the stem, SbSWEET4-3 was uniquely highly expressed. SbSWEET4-1, SbSWEET4-2, and SbSWEET4-3 were categorized into the same clade, but their tissue specificities were different, suggesting that SbSWEET4-3 is a sugar transporter with specific roles in the stem. We found a putative SWEET4-3 ortholog in the corresponding region of the maize chromosome, but not the rice chromosome, suggesting that SbSWEET4-3 was copied after the branching of sorghum and maize from rice. In the panicle from the heading through to 36 days afterward, SbSWEET2-1 and SbSWEET7-1 were expressed and grouped in the same clade as rice OsSWEET11/Xa13 that is essential for seed development. SbSWEET9-3 was highly expressed in the panicle only just after heading and was grouped into the same clade as AtSWEET8/RPG1 that is essential for pollen viability. Five of 23 SWEET genes had SNPs that caused nonsynonymous amino acid substitutions between SIL-05 and BTx623.

CONCLUSIONS

We determined the key SWEET genes for technological improvement of sorghum in the production of biofuels: SbSWEET8-1 for efflux of sucrose from the leaf; SbSWEET4-3 for unloading sucrose from the phloem in the stem; SbSWEET2-1 and SbSWEET7-1 for seed development; SbSWEET9-3 for pollen nutrition.

Functional and expression analyses of transcripts based on full-length cDNAs of Sorghum bicolor

Sorghum bicolor is one of the most important crops for food and bioethanol production. Its small diploid genome and resistance to environmental stress make sorghum an attractive model for studying the functional genomics of the Saccharinae and other C₄ grasses. We analyzed the domain-based functional annotation of the cDNAs using the gene ontology (GO) categories for molecular function to characterize all the genes cloned in the full-length cDNA library of sorghum. The sorghum cDNA library successfully captured a wide range of cDNA-encoded proteins with various functions. To characterize the protein function of newly identified cDNAs, a search of their deduced domains and comparative analyses in the Oryza sativa and Zea mays genomes were carried out. Furthermore, genes on the sense strand corresponding to antisense transcripts were classified based on the GO of molecular function. To add more information about these genes, we have analyzed the expression profiles using RNA-Seq of three tissues (spikelet, seed and stem) during the starch-filling phase. We performed functional analysis of tissue-specific genes and expression analysis of genes of starch biosynthesis enzymes. This functional analysis of sorghum full-length cDNAs and the transcriptome information will facilitate further analysis of the Saccharinae and grass families.

Inheritance and molecular mapping of Rf6 locus with pollen fertility restoration ability on A1 and A2 cytoplasms in sorghum

Of the several male sterility cytoplasms available as an alternative to the widely exploited A1 (milo) cytoplasm in sorghum, A2 is more suitable for commercial exploitation. Diversification of genetic and cytoplasmic base of hybrids involving A2 cytoplasm necessitates mapping of fertility restorer (Rf) genes for use in marker-assisted restorer development. We mapped a major male fertility restoration locus on sorghum chromosome 4 tightly linked with SSR markers, SB2387 and SB2388. This new fertility locus, Rf6, was able to restore male fertility on both A1 and A2 cytoplasms. Analysis of the genomic region around the Rf6 locus identified six genes including a pentatricopeptide repeat (PPR) gene, Sobic.004G004100. With its similar restoration ability to Rf1, Rf2 and Rf5 loci in sorghum, it is most likely that the Rf6 is a member of the PPR gene family, and the PPR gene Sobic.004G004100 could be a candidate for fertility restoration on A1 and A2 cytoplasms.

Sucrose accumulation in sweet sorghum stems occurs by apoplasmic phloem unloading and does not involve differential Sucrose transporter expression

BACKGROUND

Sorghum (Sorghum bicolor L. Moench) cultivars store non-structural carbohydrates predominantly as either starch in seeds (grain sorghums) or sugars in stems (sweet sorghums). Previous research determined that sucrose accumulation in sweet sorghum stems was not correlated with the activities of enzymes functioning in sucrose metabolism, and that an apoplasmic transport step may be involved in stem sucrose accumulation. However, the sucrose unloading pathway from stem phloem to storage parenchyma cells remains unelucidated. Sucrose transporters (SUTs) transport sucrose across membranes, and have been proposed to function in sucrose partitioning differences between sweet and grain sorghums. The purpose of this study was to characterize the key differences in carbohydrate accumulation between a sweet and a grain sorghum, to define the path sucrose may follow for accumulation in sorghum stems, and to determine the roles played by sorghum SUTs in stem sucrose accumulation.

RESULTS

Dye tracer studies to determine the sucrose transport route revealed that, for both the sweet sorghum cultivar Wray and grain sorghum cultivar Macia, the phloem in the stem veins was symplasmically isolated from surrounding cells, suggesting sucrose was apoplasmically unloaded. Once in the phloem apoplasm, a soluble tracer diffused from the vein to stem parenchyma cell walls, indicating the lignified mestome sheath encompassing the vein did not prevent apoplasmic flux outside of the vein. To characterize carbohydrate partitioning differences between Wray and Macia, we compared the growth, stem juice volume, solute contents, SbSUTs gene expression, and additional traits. Contrary to previous findings, we detected no significant differences in SbSUTs gene expression within stem tissues.

CONCLUSIONS

Phloem sieve tubes within sweet and grain sorghum stems are symplasmically isolated from surrounding cells; hence, unloading from the phloem likely occurs apoplasmically, thereby defining the location of the previously postulated step for sucrose transport. Additionally, no changes in SbSUTs gene expression were detected in sweet vs. grain sorghum stems, suggesting alterations in SbSUT transcript levels do not account for the carbohydrate partitioning differences between cultivars. A model illustrating sucrose phloem unloading and movement to stem storage parenchyma, and highlighting roles for sucrose transport proteins in sorghum stems is discussed.