GWAS of grain color and tannin content in Chinese sorghum based on whole-genome sequencing.

KEY MESSAGE: Seventy-three QTL related to grain color and tannin content were identified in Chinese sorghum accessions, and a new recessive allelic variant of TAN2 gene was discovered. Sorghum is mainly used for brewing distilled liquors in China. Since grain tannins play an important role in liquor brewing, accurately understanding the relationship between grain color and tannin content can provide basis for selection standards of tannin sorghum. We resequenced a panel of 242 Chinese sorghum accessions and performed population structure and genome-wide association study (GWAS) to identify quantitative trait locus (QTL) affecting pericarp color, testa pigment, and tannin content. Phylogenetic analysis, principal component analysis (PCA), and admixture model were used to infer population structure. Two distinct genetic sub-populations were identified according to their corresponding northern and southern geographic origin. To investigate the genetic basis of natural variation in sorghum grain color, GWAS with 2,760,264 SNPs was conducted in four environments using multiple models (Blink, FarmCPU, GLM, and MLM). Seventy-three QTL were identified to be associated for the color of exocarp, mesocarp, testa, and tannin content on all chromosomes except chromosome 5, of which 47 might be novel QTL. Some important QTL were found to colocalize with orthologous genes in the flavonoid biosynthetic pathway from other plants, including orthologous of Arabidopsis (Arabidopsis thaliana) TT2, TT7, TT12, TT16 and AT5G41220 (GST), as well as orthologous of rice (Oryza sativa) MYB61 and OsbHLH025. Our investigation of the variation in grain color and tannin content in Chinese sorghum germplasm may help guide future sorghum breeding for liquor brewing.

Sorghum qTGW1a encodes a G-protein subunit and acts as a negative regulator of grain size.

Grain size is a major determinant of grain yield in sorghum and other cereals. Over 100 quantitative trait loci (QTLs) of grain size have been identified in sorghum. However, no gene underlying any grain size QTL has been cloned. Here, we describe the fine mapping and cloning of one grain size QTL. From an F8 recombinant inbred line population derived from a cross between inbred lines 654 and LTR108, we identified 44 grain size QTLs. One QTL, qTGW1a, was detected consistently on the long arm of chromosome 1 in the span of 4 years. Using the extreme recombinants from an F2:3 fine-mapping population, qTGW1a was delimited within a ~33 kb region containing three predicted genes. One of them, SORBI_3001G341700, predicted to encode a G-protein γ subunit and homologous to GS3 in rice, is likely to be the causative gene for qTGW1a. qTGW1a appears to act as a negative regulator of grain size in sorghum. The functional allele of the putatively causative gene of qTGW1a from inbred line 654 decreased grain size, plant height, and grain yield in transgenic rice. Identification of the gene underlying qTGW1a advances our understanding of the regulatory mechanisms of grain size in sorghum and provides a target to manipulate grain size through genome editing.

Efficient generation of marker-free transgenic rice plants using an improved transposon-mediated transgene reintegration strategy.

Marker-free transgenic plants can be developed through transposon-mediated transgene reintegration, which allows intact transgene insertion with defined boundaries and requires only a few primary transformants. In this study, we improved the selection strategy and validated that the maize (Zea mays) Activator/Dissociation (Ds) transposable element can be routinely used to generate marker-free transgenic plants. A Ds-based gene of interest was linked to green fluorescent protein in transfer DNA (T-DNA), and a green fluorescent protein-aided counterselection against T-DNA was used together with polymerase chain reaction (PCR)-based positive selection for the gene of interest to screen marker-free progeny. To test the efficacy of this strategy, we cloned the Bacillus thuringiensis (Bt) δ-endotoxin gene into the Ds elements and transformed transposon vectors into rice (Oryza sativa) cultivars via Agrobacterium tumefaciens. PCR assays of the transposon empty donor site exhibited transposition in somatic cells in 60.5% to 100% of the rice transformants. Marker-free (T-DNA-free) transgenic rice plants derived from unlinked germinal transposition were obtained from the T1 generation of 26.1% of the primary transformants. Individual marker-free transgenic rice lines were subjected to thermal asymmetric interlaced-PCR to determine Ds(Bt) reintegration positions, reverse transcription-PCR and enzyme-linked immunosorbent assay to detect Bt expression levels, and bioassays to confirm resistance against the striped stem borer Chilo suppressalis. Overall, we efficiently generated marker-free transgenic plants with optimized transgene insertion and expression. The transposon-mediated marker-free platform established in this study can be used in rice and possibly in other important crops.

Characterization and molecular cloning of a serine hydroxymethyltransferase 1 (OsSHM1) in rice.

Serine hydroxymethyltransferase (SHMT) is important for one carbon metabolism and photorespiration in higher plants for its participation in plant growth and development, and resistance to biotic and abiotic stresses. A rice serine hydroxymethyltransferase gene, OsSHM1, an ortholog of Arabidopsis SHM1, was isolated using map-based cloning. The osshm1 mutant had chlorotic lesions and a considerably smaller, lethal phenotype under natural ambient CO2 concentrations, but could be restored to wild type with normal growth under elevated CO2 levels (0.5% CO2 ), showing a typical photorespiratory phenotype. The data from antioxidant enzymes activity measurement suggested that osshm1 was subjected to significant oxidative stress. Also, OsSHM1 was expressed in all organs tested (root, culm, leaf, and young panicle) but predominantly in leaves. OsSHM1 protein is localized to the mitochondria. Our study suggested that molecular function of the OsSHM1 gene is conserved in rice and Arabidopsis.

Young Leaf Chlorosis 2 encodes the stroma-localized heme oxygenase 2 which is required for normal tetrapyrrole biosynthesis in rice.

MAIN CONCLUSION: Rice heme oxygenase 2 (OsHO2) mutants are chlorophyll deficient with distinct tetrapyrrole metabolite and transcript profiles, suggesting a potential regulatory role of the stromal-localized OsHO2 in tetrapyrrole biosynthesis. In plants, heme oxygenases (HOs) are classified into the subfamilies HO1 and HO2. HO1 are highly conserved plastid enzymes required for synthesizing the chromophore in phytochromes which mediate a number of light-regulated responses. However, the physiological and biochemical functions of HO2, which are distantly related to HO1, are not well understood, especially in crop plants. From a population of (60)Coγ-irradiated rice mutants, we identified the ylc2 (young leaf chlorosis 2) mutant which displays a chlorosis phenotype in seedlings with substantially reduced chlorophyll content. Normal leaf pigmentation is gradually restored in older plants while newly emerged leaves remain yellow. Transmission electron microscopy further revealed defective chloroplast structures in the ylc2 seedlings. Map-based cloning located the OsYLC2 gene on chromosome 3 and it encodes the OsHO2 protein. The gene identification was confirmed by complementation and T-DNA mutant analyses. Subcellular localization and chloroplast fractionation experiments indicated that OsHO2 resides in the stroma. However, recombinant enzyme assay demonstrated that OsHO2 is not a functional HO enzyme. Analysis of tetrapyrrole metabolites revealed the reduced levels of most chlorophyll and phytochromobilin precursors in the ylc2 mutant. On the other hand, elevated accumulation of 5-aminolevulinic acid and Mg-protoporphyrin IX was observed. These unique metabolite changes are accompanied by consistent changes in the expression levels of the corresponding tetrapyrrole biosynthesis genes. Taken together, our work suggests that OsHO2 has a potential regulatory role for tetrapyrrole biosynthesis in rice.

Rice RNA-dependent RNA polymerase 6 acts in small RNA biogenesis and spikelet development.

Higher plants have evolved multiple RNA-dependent RNA polymerases (RDRs), which work with Dicer-like (DCL) proteins to produce different classes of small RNAs with specialized molecular functions. Here we report that OsRDR6, the rice (Oryza sativa L.) homolog of Arabidopsis RDR6, acts in the biogenesis of various types and sizes of small RNAs. We isolated a rice osrdr6-1 mutant, which was temperature sensitive and showed spikelet defects. This mutant displays reduced accumulation of tasiR-ARFs, the conserved trans-acting siRNAs (tasiRNAs) derived from the TAS3 locus, and ectopic expression of tasiR-ARF target genes, the Auxin Response Factors (including ARF2 and ARF3/ETTIN). The loss of tasiR-mediated repression of ARFs in osrdr6-1 can explain its morphological defects, as expression of two non-targeted ARF3 gene constructs (ARF3muts) in a wild-type background mimics the osrdr6 and osdcl4-1 mutant phenotypes. Small RNA high-throughput sequencing also reveals that besides tasiRNAs, 21-nucleotide (nt) phased small RNAs are also largely dependent on OsRDR6. Unexpectedly, we found that osrdr6-1 has a strong impact on the accumulation of 24-nt phased small RNAs, but not on unphased ones. Our work uncovers the key roles of OsRDR6 in small RNA biogenesis and directly illustrates the crucial functions of tasiR-ARFs in rice development.

Identification of QTLs associated with tissue culture response through sequencing-based genotyping of RILs derived from 93-11 × Nipponbare in rice (Oryza sativa).

KEY MESSAGE : The performance of callus induction and callus differentiation was evaluated by 9 indices for 140 RILs; 2 major QTLs associated with plant regeneration were identified. In order to investigate the genetic mechanisms of tissue culture response, 140 recombinant inbred lines (RILs) derived from 93-11 (Oryza sativa ssp. indica) × Nipponbare (Oryza sativa ssp. japonica) and a high quality genetic map based on the SNPs generated from deep sequencing of the RIL genomes, were used to identify the quantitative trait loci (QTLs) associated with in vitro tissue culture response (TCR) from mature seed in rice. The performance of callus induction was evaluated by indices of induced-callus color (ICC), induced-callus size (ICS), induced-callus friability (ICF) and callus induction rate (CIR), respectively, and the performance of callus differentiation was evaluated by indices of callus proliferation ability (CPA), callus browning tendency (CBT), callus greening ability (CGA), the average number of regenerated shoots per callus (NRS) and regeneration rate (%, RR), respectively. A total of 25 QTLs, 2 each for ICC, ICS, ICF, CIR and CBA, 3 for CPA, 4 each for CGA, NRS and RR, respectively, were detected and located on 8 rice chromosomes. Significant correlations were observed among the traits of CGA, NRS and RR, and QTLs identified for these three indices were co-located on chromosomes 3 and 7, and the additive effects came from both Nipponbare and 93-11, respectively. The results obtained from this study provide guidance for further fine mapping and gene cloning of the major QTL of TCR and the knowledge of the genes underlying the traits investigated would be very helpful for revealing the molecular bases of tissue culture response.

Identification of QTLs for eight agronomically important traits using an ultra-high-density map based on SNPs generated from high-throughput sequencing in sorghum under contrasting photoperiods.

The productivity of sorghum is mainly determined by agronomically important traits. The genetic bases of these traits have historically been dissected and analysed through quantitative trait locus (QTL) mapping based on linkage maps with low-throughput molecular markers, which is one of the factors that hinder precise and complete information about the numbers and locations of the genes or QTLs controlling the traits. In this study, an ultra-high-density linkage map based on high-quality single nucleotide polymorphisms (SNPs) generated from low-coverage sequences (~0.07 genome sequence) in a sorghum recombinant inbred line (RIL) population was constructed through new sequencing technology. This map consisted of 3418 bin markers and spanned 1591.4 cM of genome size with an average distance of 0.5 cM between adjacent bins. QTL analysis was performed and a total of 57 major QTLs were detected for eight agronomically important traits under two contrasting photoperiods. The phenotypic variation explained by individual QTLs varied from 3.40% to 33.82%. The high accuracy and quality of this map was evidenced by the finding that genes underlying two cloned QTLs, Dw3 for plant height (chromosome 7) and Ma1 for flowering time (chromosome 6), were localized to the correct genomic regions. The close associations between two genomic regions on chromosomes 6 and 7 with multiple traits suggested the existence of pleiotropy or tight linkage. Several major QTLs for heading date, plant height, numbers of nodes, stem diameter, panicle neck length, and flag leaf width were detected consistently under both photoperiods, providing useful information for understanding the genetic mechanisms of the agronomically important traits responsible for the change of photoperiod.

Selection Signatures in Chinese Sorghum Reveals Its Unique Liquor-Making Properties.

Chinese sorghum (S. bicolor) has been a historically critical ingredient for brewing famous distilled liquors ever since Yuan Dynasty (749 ∼ 652 years BP). Incomplete understanding of the population genetics and domestication history limits its broad applications, especially that the lack of genetics knowledge underlying liquor-brewing properties makes it difficult to establish scientific standards for sorghum breeding. To unravel the domestic history of Chinese sorghum, we re-sequenced 244 Chinese sorghum lines selected from 16 provinces. We found that Chinese sorghums formed three distinct genetic sub-structures, referred as the Northern, the Southern, and the Chishui groups, following an obviously geographic pattern. These sorghum accessions were further characterized in liquor brewing traits and identified selection footprints associated with liquor brewing efficiency. An importantly selective sweep region identified includes several homologous genes involving in grain size, pericarp thickness, and architecture of inflorescence. Our result also demonstrated that pericarp strength rather than grain size determines the ability of the grains to resist repeated cooking during brewing process. New insight into the traits beneficial to the liquor-brewing process provides both a better understanding on Chinese sorghum domestication and a guidance on breeding sorghum as a multiple use crop in China.

Seed size is determined by the combinations of the genes controlling different seed characteristics in rice.

Rice seed size is an important agronomic trait in determining the yield potential, and four seed size related genes (GS3, GW2, qSW5/GW5 and GIF1) have been cloned in rice so far. However, the relationship among these four genes is still unclear, which will impede the process of gene pyramiding breeding program to some extent. To shade light on the relationship of above four genes, gene expression analysis was performed with GS3-RNAi, GW2-RNAi lines and CSSL of qSW5 at the transcriptional level. The results clearly showed that qSW5 and GW2 positively regulate the expression of GS3. Meanwhile, qSW5 can be down-regulated by repression of GW2 transcription. Additionally, GIF1 expression was found to be positively regulated by qSW5 but negatively by GW2 and GS3. Moreover, the allelic effects of qSW5 and GS3 were detailedly characterized based on a natural population consisting of 180 rice cultivars. It was indicated that mutual interactions exist between the two genes, in which, qSW5 affecting seed length is masked by GS3 alleles, and GS3 affecting seed width is masked by qSW5 alleles. These findings provide more insights into the molecular mechanisms underlying seed size development in rice and are likely to be useful for improving rice grain yield.

Molecular dissection of the pathogen-inducible 3-deoxyanthocyanidin biosynthesis pathway in sorghum.

3-Deoxyanthocyanidins are the unique phytoalexins synthesized by sorghum in response to fungal inoculation. They are structurally related to anthocyanins but the final steps of their pathogen-inducible biosynthesis are not fully understood. We have identified new flavonoid structural genes from the recently completed sorghum BTx623 genome sequence. The biochemical functions of the different expressed sorghum genes were established in planta by complementation in the appropriate Arabidopsis transparent testa mutants. There is a family of nine chalcone synthase genes which are all inducible by fungal inoculation in sorghum seedlings. Specific dihydroflavonol 4-reductase (DFR) genes responsive to conditions which stimulated anthocyanin accumulation (SbDFR1) or 3-deoxyanthocyanidin production (SbDFR3) were identified. Recombinant SbDFR1 and SbDFR3 were found to function as typical DFRs by accepting dihydroflavonol substrates. On the other hand, both DFRs showed substantially lower but detectable NADPH-dependent activities toward flavanones. Reduction of flavanones to flavan-4-ols is a reaction step required for 3-deoxyanthocyanidin production. Flavanone 3-hydroxylase (F3H) converts flavanones to dihydroflavonols for anthocyanin biosynthesis. In sorghum seedlings, expression of two F3H genes was either absent or strongly suppressed during the accumulation of 3-deoxyanthocyanidins. Under such conditions, most flavanones are expected to be reduced by the pathogen-induced SbDFR3 for the formation of flavan-4-ols. Our work also revealed that 3-deoxyanthocyanidin accumulation and SbDFR3 expression were induced by methyl jasmonate treatment in sorghum roots but the stimulation effects were antagonized by salicylic acid.

OsNOA1/RIF1 is a functional homolog of AtNOA1/RIF1: implication for a highly conserved plant cGTPase essential for chloroplast function.

*The bacterial protein YqeH is a circularly permuted GTPase with homologs encoded by plant nuclear genomes. The rice homolog OsNOA1/RIF1 is encoded by the single-copy gene Os02g01440. OsNOA1/RIF1 is expressed in different tissues and is light-inducible. The OsNOA1/RIF1-EYFP fusion protein was targeted to chloroplasts in transgenic Arabidopsis plants. In addition, the rice homolog was able to rescue most of the growth phenotypes in an Arabidopsis rif1 mutant. *Rice (Oryza sativa) OsNOA1/RIF1 RNAi mutant seedlings were chlorotic with reduced pigment contents and lower photosystem II (PSII) efficiency. However, the expressions of the chloroplast-encoded genes rbcL, atpB, psaA and psbA were not affected. By contrast, reduced abundance of the chloroplast 16S rRNA was observed in the mutant. *Quantitative iTRAQ-LC-MS/MS proteomics investigations revealed proteome changes in the rice mutant consistent with the expected functional role of OsNOA1/RIF1 in chloroplast translation. The RNAi mutant showed significantly decreased expression levels of chloroplast-encoded proteins as well as nuclear-encoded components of chloroplast enzyme complexes. Conversely, upregulation of some classes of nonchloroplastic proteins, such as glycolytic and phenylpropanoid pathway enzymes, was detected. *Our work provides independent indications that a highly conserved nuclear-encoded cGTPase of likely prokaryotic origin is essential for proper chloroplast ribosome assembly and/or translation in plants.

Exogenous gamma-aminobutyric acid alleviates oxidative damage caused by aluminium and proton stresses on barley seedlings.

BACKGROUND: Proton (H(+)) and aluminium (Al(3+)) toxicities are major factors limiting crop production on acid soils, while gamma-aminobutyric acid (GABA) is a non-protein amino acid involved in various stress tolerances in plants. In this study, to determine whether exogenous GABA is functional in alleviating oxidative stress induced by H(+) and Al(3+) toxicities, the antioxidant defence response regulated by GABA was investigated in barley (Hordeum vulgare L.). RESULTS: After 24 h treatments of seedlings under H(+), Al(3+) and combined stresses with and without GABA, morphological and biochemical assays were conducted. It was observed that the inhibition of seedling root elongation caused by Al(3+) and H(+) toxicities was significantly mitigated by GABA. The amount of carbonylated proteins with molecular masses of 14.4-97 kDa was decreased. The activities of antioxidant enzymes were enhanced, the content of malondialdehyde was reduced and the accumulation of reactive oxygen species (ROS), as shown by staining roots with nitroblue tetrazolium, declined in GABA-treated seedlings. CONCLUSION: GABA can alleviate oxidative damage caused by H(+) and Al(3+) toxicities in barley seedlings by activating antioxidant defence responses and reducing the elevated levels of carbonylated proteins caused by ROS.

Genome-wide analysis of the auxin response factors (ARF) gene family in rice (Oryza sativa).

Auxin response factors (ARFs) are transcription factors that bind with specificity to TGTCTC-containing auxin response elements (AuxREs) found in promoters of primary/early auxin response genes and mediate responses to the plant hormone auxin. The ARF genes are represented by a large multigene family in plants. A comprehensive genome-wide analysis was carried out in this study to find all ARFs in Arabidopsis (Arabidopsis thaliana) and rice (Oryza sativa subsp. japonica), 23 and 25 ARF genes, named as AtARFs and OsARFs, were identified, respectively. Chromosomal locations of all OsARFs were presented and it was found that the duplication of OsARFs was associated with only the chromosomal block duplications but not local tandem duplications. A phylogenetic tree was generated from alignments of the full-length protein sequences of 25 OsARFs and 23 AtARFs to examine the phylogenetic relationships of rice and Arabidopsis ARF proteins. All 48 members of ARF gene families fell into three major classes, a total of 13 sister pairs, including 9 OsARF-OsARF, 2 AtARF-AtARF and 2 AtARF-OsARF sister pairs were formed, showing different orthologous relationships between AtARFs and OsARFs. EST analysis and RT-PCR assays demonstrated that 24 of all 25 OsARF genes were active and the transcript abundance of some OsARF genes was affected by auxin treatment or light- and dark-grown conditions. The outcome of the present study provides basic genomic information for the rice ARF gene family and will pave the way for elucidating the precise role of OsARFs in plant growth and development in the future.

Molecular cloning and expression analysis of duplicated polyphenol oxidase genes reveal their functional differentiations in sorghum.

Polyphenol oxidase (PPO) is believed to play a role in plant growth, reproduction, and resistance to pathogens and pests. PPO causes browning of grains in cereals. In this study, genetic mapping of sorghum grain for phenol color reaction (PHR) was performed using a recombinant inbred line population. Only one locus was detected between SSR markers SM06072 and Xtxp176 on chromosome 6. Two linked orthologous genes (Sb06PPO1 and Sb06PPO2) within the mapped region were discovered and cloned. Transformation experiments using Nipponbare (a PHR negative rice cultivar) showed that Sb06PPO1 from LTR108 and two Sb06PPO2 alleles from both varieties could complement Nipponbare, whereas Sb06PPO1 from 654 could not. Subsequent quantitative real-time PCR (qPCR) experiments showed that Sb06PPO1 and Sb06PPO2 functioned diversely, Sb06PPO1 was mainly expressed in young panicles before flowering. Sb06PPO2 was strongly expressed in flowering panicles, especially in hulls and branches at filling stage. Moreover, the expression of Sb06PPO1 was found to be significantly up-regulated by exogenous ABA and salt, whereas Sb06PPO2 was not changed significantly, further demonstrating functional differentiation between the two genes.

Application of TILLING in plant improvement.

TILLING (Targeting induced local lesions in genomes) is a general reverse-genetic strategy that is used to locate an allelic series of induced point mutations in genes of interest. High-throughput TILLING allows the rapid and cost-effective detection of induced point mutations in populations of chemically mutagenized individuals. The technique can be applied not only to model organisms but also to economically important organisms in plants. Owing to its full of advantages such as simple procedure, high sensitivity, and high efficiency, TILLING provides a powerful approach for gene discovery, DNA polymorphism assessment, and plant improvement. Coupled with other genomic resources, TILLING and EcoTILLING can be used immediately as a haplotyping tool in plant breeding for identifying allelic variation in genes exhibiting expression correlating with phenotypes and establishing an allelic series at genetic loci for the traits of interest in germplasm or induced mutants.

OspTAC2 encodes a pentatricopeptide repeat protein and regulates rice chloroplast development.

Functional chloroplast generation depends on the precise coordination of gene expression between the plastid and the nucleus and is essential for plant growth and development. In this study, a rice (Oryza sativa) mutant that exhibited albino and seedling-lethal phenotypes was isolated from a60Co-irradiated rice population. The mutant gene was identified as an ortholog of the Arabidopsis plastid transcriptionally active chromosome protein 2 (pTAC2) gene, and the mutant strain was designated osptac2. Sequence and transcription analyses showed that OspTAC2 encodes a putative chloroplast protein consisting of 10 pentratricopeptide repeat (PPR) domains and a C-terminal small MutS-related (SMR) domain. Cytological observations via microscopy showed that the OspTAC2-green fluorescent fusion protein is localized in the chloroplasts. Transmission electron microscopy revealed that the chloroplast of the osptac2 mutant lacks an organized thylakoid membrane. The transcript levels of all investigated PEP (plastid-encoded RNA polymerase)-dependent genes were dramatically reduced in the osptac2 mutant, whereas the transcript levels of NEP (nuclear-encoded polymerase)-dependent genes were increased. These results suggest that OspTAC2 plays a critical role in chloroplast development and indicate that the molecular function of the OspTAC2 gene is conserved in rice and Arabidopsis.

Assessment of genetic diversity by simple sequence repeat markers among forty elite varieties in the germplasm for malting barley breeding.

The genetic diversity and relationship among 40 elite barley varieties were analyzed based on simple sequence repeat (SSR) genotyping data. The amplified fragments from SSR primers were highly polymorphic in the barley accessions investigated. A total of 85 alleles were detected at 35 SSR loci, and allelic variations existed at 29 SSR loci. The allele number per locus ranged from 1 to 5 with an average of 2.4 alleles per locus detected from the 40 barley accessions. A cluster analysis based on the genetic similarity coefficients was conducted and the 40 varieties were classified into two groups. Seven malting barley varieties from China fell into the same subgroup. It was found that the genetic diversity within the Chinese malting barley varieties was narrower than that in other barley germplasm sources, suggesting the importance and feasibility of introducing elite genotypes from different origins for malting barley breeding in China.