Institute of Crop Sciences

QTL IciMapping is freely available public software capable of building high-density linkage maps and mapping quantitative trait loci (QTL) in biparental populations. Eight functionalities are integrated in this software package: (1) BIN: binning of redundant markers; (2) MAP: construction of linkage maps in biparental populations; (3) CMP: consensus map construction from multiple linkage maps sharing common markers; (4) SDL: mapping of segregation distortion loci; (5) BIP: mapping of additive, dominant, and digenic epistasis genes; (6) MET: QTL-by-environment interaction analysis; (7) CSL: mapping of additive and digenic epistasis genes with chromosome segment substitution lines; and (8) NAM: QTL mapping in NAM populations. Input files can be arranged in plain text, MS Excel 2003, or MS Excel 2007 formats. Output files have the same prefix name as the input but with different extensions. As examples, there are two output files in BIN, one for summarizing the identified bin groups and deleted markers in each bin, and the other for using the MAP functionality. Eight output files are generated by MAP, including summary of the completed linkage maps, Mendelian ratio test of individual markers, estimates of recombination frequencies, LOD scores, and genetic distances, and the input files for using the BIP, SDL, and MET functionalities. More than 30 output files are generated by BIP, including results at all scanning positions, identified QTL, permutation tests, and detection powers for up to six mapping methods. Three supplementary tools have also been developed to display completed genetic linkage maps, to estimate recombination frequency between two loci, and to perform analysis of variance for multi-environmental trials.

Here we analyse genetic variation, population structure and diversity among 3,010 diverse Asian cultivated rice (Oryza sativa L.) genomes from the 3,000 Rice Genomes Project. Our results are consistent with the five major groups previously recognized, but also suggest several unreported subpopulations that correlate with geographic location. We identified 29 million single nucleotide polymorphisms, 2.4 million small indels and over 90,000 structural variations that contribute to within- and between-population variation. Using pan-genome analyses, we identified more than 10,000 novel full-length protein-coding genes and a high number of presence-absence variations. The complex patterns of introgression observed in domestication genes are consistent with multiple independent rice domestication events. The public availability of data from the 3,000 Rice Genomes Project provides a resource for rice genomics research and breeding.

Drought and salinity are major abiotic stresses to crop production. Here, we show that overexpression of stress responsive gene SNAC1 (STRESS-RESPONSIVE NAC 1) significantly enhances drought resistance in transgenic rice (22-34% higher seed setting than control) in the field under severe drought stress conditions at the reproductive stage while showing no phenotypic changes or yield penalty. The transgenic rice also shows significantly improved drought resistance and salt tolerance at the vegetative stage. Compared with WT, the transgenic rice are more sensitive to abscisic acid and lose water more slowly by closing more stomatal pores, yet display no significant difference in the rate of photosynthesis. SNAC1 is induced predominantly in guard cells by drought and encodes a NAM, ATAF, and CUC (NAC) transcription factor with transactivation activity. DNA chip analysis revealed that a large number of stress-related genes were up-regulated in the SNAC1-overexpressing rice plants. Our data suggest that SNAC1 holds promising utility in improving drought and salinity tolerance in rice.

Flowering time is a complex trait that controls adaptation of plants to their local environment in the outcrossing species Zea mays (maize). We dissected variation for flowering time with a set of 5000 recombinant inbred lines (maize Nested Association Mapping population, NAM). Nearly a million plants were assayed in eight environments but showed no evidence for any single large-effect quantitative trait loci (QTLs). Instead, we identified evidence for numerous small-effect QTLs shared among families; however, allelic effects differ across founder lines. We identified no individual QTLs at which allelic effects are determined by geographic origin or large effects for epistasis or environmental interactions. Thus, a simple additive model accurately predicts flowering time for maize, in contrast to the genetic architecture observed in the selfing plant species rice and Arabidopsis.

Maize genetic diversity has been used to understand the molecular basis of phenotypic variation and to improve agricultural efficiency and sustainability. We crossed 25 diverse inbred maize lines to the B73 reference line, capturing a total of 136,000 recombination events. Variation for recombination frequencies was observed among families, influenced by local (cis) genetic variation. We identified evidence for numerous minor single-locus effects but little two-locus linkage disequilibrium or segregation distortion, which indicated a limited role for genes with large effects and epistatic interactions on fitness. We observed excess residual heterozygosity in pericentromeric regions, which suggested that selection in inbred lines has been less efficient in these regions because of reduced recombination frequency. This implies that pericentromeric regions may contribute disproportionally to heterosis.

Abstract Composite interval mapping (CIM) is the most commonly used method for mapping quantitative trait loci (QTL) with populations derived from biparental crosses. However, the algorithm implemented in the popular QTL Cartographer software may not completely ensure all its advantageous properties. In addition, different background marker selection methods may give very different mapping results, and the nature of the preferred method is not clear. A modified algorithm called inclusive composite interval mapping (ICIM) is proposed in this article. In ICIM, marker selection is conducted only once through stepwise regression by considering all marker information simultaneously, and the phenotypic values are then adjusted by all markers retained in the regression equation except the two markers flanking the current mapping interval. The adjusted phenotypic values are finally used in interval mapping (IM). The modified algorithm has a simpler form than that used in CIM, but a faster convergence speed. ICIM retains all advantages of CIM over IM and avoids the possible increase of sampling variance and the complicated background marker selection process in CIM. Extensive simulations using two genomes and various genetic models indicated that ICIM has increased detection power, a reduced false detection rate, and less biased estimates of QTL effects.

Rice blast is one of the most destructive diseases affecting rice worldwide. The adoption of host resistance has proven to be the most economical and effective approach to control rice blast. In recent years, sequence-specific nucleases (SSNs) have been demonstrated to be powerful tools for the improvement of crops via gene-specific genome editing, and CRISPR/Cas9 is thought to be the most effective SSN. Here, we report the improvement of rice blast resistance by engineering a CRISPR/Cas9 SSN (C-ERF922) targeting the OsERF922 gene in rice. Twenty-one C-ERF922-induced mutant plants (42.0%) were identified from 50 T0 transgenic plants. Sanger sequencing revealed that these plants harbored various insertion or deletion (InDel) mutations at the target site. We showed that all of the C-ERF922-induced allele mutations were transmitted to subsequent generations. Mutant plants harboring the desired gene modification but not containing the transferred DNA were obtained by segregation in the T1 and T2 generations. Six T2 homozygous mutant lines were further examined for a blast resistance phenotype and agronomic traits, such as plant height, flag leaf length and width, number of productive panicles, panicle length, number of grains per panicle, seed setting percentage and thousand seed weight. The results revealed that the number of blast lesions formed following pathogen infection was significantly decreased in all 6 mutant lines compared with wild-type plants at both the seedling and tillering stages. Furthermore, there were no significant differences between any of the 6 T2 mutant lines and the wild-type plants with regard to the agronomic traits tested. We also simultaneously targeted multiple sites within OsERF922 by using Cas9/Multi-target-sgRNAs (C-ERF922S1S2 and C-ERF922S1S2S3) to obtain plants harboring mutations at two or three sites. Our results indicate that gene modification via CRISPR/Cas9 is a useful approach for enhancing blast resistance in rice.

Jasmonic acid (JA) and its precursors and dervatives, referred as jasmonates (JAs) are important molecules in the regulation of many physiological processes in plant growth and development, and especially the mediation of plant responses to biotic and abiotic stresses. JAs biosynthesis, perception, transport, signal transduction and action have been extensively investigated. In this review, we will discuss the initiation of JA signaling with a focus on environmental signal perception and transduction, JA biosynthesis and metabolism, transport of signaling molecules (local transmission, vascular bundle transmission, and airborne transportation), and biological function (JA signal receptors, regulated transcription factors, and biological processes involved).

Sequencing and analysing the diploid genome and transcriptome of Aegilops tauschii provide new insights into the role of this genome in enabling the adaptation of bread wheat and are a step towards understanding the very large and complicated hexaploid genomes of wheat species. The hexaploid genome of bread wheat Triticum aestivum, designated AABBDD, evolved as a result of hybridization between three ancestral grasses. Two papers published in the issue of Nature present genome sequences and analysis of two of these wheat progenitors. First, the genome sequence of the diploid wild wheat T. urartu (ancestor of the A genome), which resembles cultivated wheat more strongly than either Aegilops speltoides (the B ancestor) or Ae. tauschii (the D donor). And second, the Ae. tauschii genome, together with an analysis of its transcriptome. These genomes and their analyses will be powerful tools for the study of complex, polyploid wheat genomes and a valuable resource for genetic improvement of wheat. About 8,000 years ago in the Fertile Crescent, a spontaneous hybridization of the wild diploid grass Aegilops tauschii (2n = 14; DD) with the cultivated tetraploid wheat Triticum turgidum (2n = 4x = 28; AABB) resulted in hexaploid wheat (T. aestivum; 2n = 6x = 42; AABBDD)1,2. Wheat has since become a primary staple crop worldwide as a result of its enhanced adaptability to a wide range of climates and improved grain quality for the production of baker’s flour2. Here we describe sequencing the Ae. tauschii genome and obtaining a roughly 90-fold depth of short reads from libraries with various insert sizes, to gain a better understanding of this genetically complex plant. The assembled scaffolds represented 83.4% of the genome, of which 65.9% comprised transposable elements. We generated comprehensive RNA-Seq data and used it to identify 43,150 protein-coding genes, of which 30,697 (71.1%) were uniquely anchored to chromosomes with an integrated high-density genetic map. Whole-genome analysis revealed gene family expansion in Ae. tauschii of agronomically relevant gene families that were associated with disease resistance, abiotic stress tolerance and grain quality. This draft genome sequence provides insight into the environmental adaptation of bread wheat and can aid in defining the large and complicated genomes of wheat species.

Wild relatives of crops are an important source of genetic diversity for agriculture, but their gene repertoire remains largely unexplored. We report the establishment and analysis of a pan-genome of Glycine soja, the wild relative of cultivated soybean Glycine max, by sequencing and de novo assembly of seven phylogenetically and geographically representative accessions. Intergenomic comparisons identified lineage-specific genes and genes with copy number variation or large-effect mutations, some of which show evidence of positive selection and may contribute to variation of agronomic traits such as biotic resistance, seed composition, flowering and maturity time, organ size and final biomass. Approximately 80% of the pan-genome was present in all seven accessions (core), whereas the rest was dispensable and exhibited greater variation than the core genome, perhaps reflecting a role in adaptation to diverse environments. This work will facilitate the harnessing of untapped genetic diversity from wild soybean for enhancement of elite cultivars.

Abstract Grassland diversity can support sustainable intensification of grassland production through increased yields, reduced inputs and limited weed invasion. We report the effects of diversity on weed suppression from 3 years of a 31‐site continental‐scale field experiment. At each site, 15 grassland communities comprising four monocultures and 11 four‐species mixtures based on a wide range of species' proportions were sown at two densities and managed by cutting. Forage species were selected according to two crossed functional traits, “method of nitrogen acquisition” and “pattern of temporal development”. Across sites, years and sown densities, annual weed biomass in mixtures and monocultures was 0.5 and 2.0 t DM ha −1 (7% and 33% of total biomass respectively). Over 95% of mixtures had weed biomass lower than the average of monocultures, and in two‐thirds of cases, lower than in the most suppressive monoculture (transgressive suppression). Suppression was significantly transgressive for 58% of site‐years. Transgressive suppression by mixtures was maintained across years, independent of site productivity. Based on models, average weed biomass in mixture over the whole experiment was 52% less (95% confidence interval: 30%–75%) than in the most suppressive monoculture. Transgressive suppression of weed biomass was significant at each year across all mixtures and for each mixture. Weed biomass was consistently low across all mixtures and years and was in some cases significantly but not largely different from that in the equiproportional mixture. The average variability (standard deviation) of annual weed biomass within a site was much lower for mixtures (0.42) than for monocultures (1.77). Synthesis and applications . Weed invasion can be diminished through a combination of forage species selected for complementarity and persistence traits in systems designed to reduce reliance on fertiliser nitrogen. In this study, effects of diversity on weed suppression were consistently strong across mixtures varying widely in species' proportions and over time. The level of weed biomass did not vary greatly across mixtures varying widely in proportions of sown species. These diversity benefits in intensively managed grasslands are relevant for the sustainable intensification of agriculture and, importantly, are achievable through practical farm‐scale actions.

A new member of the AP2/ERF transcription factor family, GmERF3, was isolated from soybean. Sequence analysis showed that GmERF3 contained an AP2/ERF domain of 58 amino acids and two putative nuclear localization signal (NLS) domains. It belonged to a group IV protein in the ERF (ethylene response factor) subfamily as typified by a conserved N-terminal motif [MCGGAI(I/L)]. Expression of GmERF3 was induced by treatments with high salinity, drought, abscisic acid (ABA), salicylic acid (SA), jasmonic acid (JA), ethylene (ET), and soybean mosaic virus (SMV), whereas there was no significant GmERF3 mRNA accumulation under cold stress treatment. GmERF3 could bind to the GCC box and DRE/CRT element, and was targeted to the nucleus when transiently expressed in onion epidermal cells. The GmERF3 protein fused to the GAL4 DNA-binding domain to activate transcription of reporter genes in yeast. Ectopic expression of the GmERF3 gene in transgenic tobacco plants induced the expression of some PR genes and enhanced resistance against infection by Ralstonia solanacearum, Alternaria alternata, and tobacco mosaic virus (TMV), and gave tolerance to high salinity and dehydration stresses. Furthermore, overexpression of GmERF3 in transgenic tobacco led to higher levels of free proline and soluble carbohydrates compared to wild-type plants under drought conditions. The overall results suggested that GmERF3 as an AP2/ERF transcription factor may play dual roles in response to biotic and abiotic stresses in plants.

Rice blast, caused by the fungal pathogen Magnaporthe grisea, is one of the most devastating diseases in rice worldwide. The dominant resistance gene, Pi-d2 [previously named Pi-d(t)2], present in the rice variety Digu, confers gene-for-gene resistance to the Chinese blast strain, ZB15. Pi-d2 was previously mapped close to the centromere of chromosome 6. In this study, the Pi-d2 gene was isolated by a map-based cloning strategy. Pi-d2 encodes a receptor-like kinase protein with a predicted extracellular domain of a bulb-type mannose specific binding lectin (B-lectin) and an intracellular serine-threonine kinase domain. Pi-d2 is a single-copy gene that is constitutively expressed in the rice variety Digu. Transgenic plants carrying the Pi-d2 transgene confer race-specific resistance to the M. grisea strain, ZB15. The Pi-d2 protein is plasma membrane localized. A single amino acid difference at position 441 of Pi-d2 distinguishes resistant and susceptible alleles of rice blast resistance gene Pi-d2. Because of its novel extracellular domain, Pi-d2 represents a new class of plant resistance genes.

Abstract Genetic diversity is key to crop improvement. Owing to pervasive genomic structural variation, a single reference genome assembly cannot capture the full complement of sequence diversity of a crop species (known as the ‘pan-genome’ 1 ). Multiple high-quality sequence assemblies are an indispensable component of a pan-genome infrastructure. Barley ( Hordeum vulgare L.) is an important cereal crop with a long history of cultivation that is adapted to a wide range of agro-climatic conditions 2 . Here we report the construction of chromosome-scale sequence assemblies for the genotypes of 20 varieties of barley—comprising landraces, cultivars and a wild barley—that were selected as representatives of global barley diversity. We catalogued genomic presence/absence variants and explored the use of structural variants for quantitative genetic analysis through whole-genome shotgun sequencing of 300 gene bank accessions. We discovered abundant large inversion polymorphisms and analysed in detail two inversions that are frequently found in current elite barley germplasm; one is probably the product of mutation breeding and the other is tightly linked to a locus that is involved in the expansion of geographical range. This first-generation barley pan-genome makes previously hidden genetic variation accessible to genetic studies and breeding.

The three most important agronomic traits of rice (Oryza sativa), yield, plant height, and flowering time, are controlled by many quantitative trait loci (QTLs). In this study, a newly identified QTL, DTH8 (QTL for days to heading on chromosome 8), was found to regulate these three traits in rice. Map-based cloning reveals that DTH8 encodes a putative HAP3 subunit of the CCAAT-box-binding transcription factor and the complementary experiment increased significantly days to heading, plant height, and number of grains per panicle in CSSL61 (a chromosome segment substitution line that carries the nonfunctional DTH8 allele) with the Asominori functional DTH8 allele under long-day conditions. DTH8 is expressed in most tissues and its protein is localized to the nucleus exclusively. The quantitative real-time PCR assay revealed that DTH8 could down-regulate the transcriptions of Ehd1 (for Early heading date1) and Hd3a (for Heading date3a; a rice ortholog of FLOWERING LOCUS T) under long-day conditions. Ehd1 and Hd3a can also be down-regulated by the photoperiodic flowering genes Ghd7 and Hd1 (a rice ortholog of CONSTANS). Meanwhile, the transcription of DTH8 has been proved to be independent of Ghd7 and Hd1, and the natural mutation of this gene caused weak photoperiod sensitivity and shorter plant height. Taken together, these data indicate that DTH8 probably plays an important role in the signal network of photoperiodic flowering as a novel suppressor as well as in the regulation of plant height and yield potential.

Chlorophyll (Chl) synthase catalyzes esterification of chlorophyllide to complete the last step of Chl biosynthesis. Although the Chl synthases and the corresponding genes from various organisms have been well characterized, Chl synthase mutants have not yet been reported in higher plants. In this study, a rice (Oryza Sativa) Chl-deficient mutant, yellow-green leaf1 (ygl1), was isolated, which showed yellow-green leaves in young plants with decreased Chl synthesis, increased level of tetrapyrrole intermediates, and delayed chloroplast development. Genetic analysis demonstrated that the phenotype of ygl1 was caused by a recessive mutation in a nuclear gene. The ygl1 locus was mapped to chromosome 5 and isolated by map-based cloning. Sequence analysis revealed that it encodes the Chl synthase and its identity was verified by transgenic complementation. A missense mutation was found in a highly conserved residue of YGL1 in the ygl1 mutant, resulting in reduction of the enzymatic activity. YGL1 is constitutively expressed in all tissues, and its expression is not significantly affected in the ygl1 mutant. Interestingly, the mRNA expression of the cab1R gene encoding the Chl a/b-binding protein was severely suppressed in the ygl1 mutant. Moreover, the expression of some nuclear genes associated with Chl biosynthesis or chloroplast development was also affected in ygl1 seedlings. These results indicate that the expression of nuclear genes encoding various chloroplast proteins might be feedback regulated by the level of Chl or Chl precursors.

Cereals high in amylose content (AC) and resistant starch (RS) offer potential health benefits. Previous studies using chemical mutagenesis or RNA interference have demonstrated that starch branching enzyme (SBE) plays a major role in determining the fine structure and physical properties of starch. However, it remains a challenge to control starch branching in commercial lines. Here, we use CRISPR/Cas9 technology to generate targeted mutagenesis in SBEI and SBEIIb in rice. The frequencies of obtained homozygous or bi-allelic mutant lines with indels in SBEI and SBEIIb in T0 generation were from 26.7 to 40%. Mutations in the homozygous T0 lines stably transmitted to the T1 generation and those in the bi-allelic lines segregated in a Mendelian fashion. Transgene-free plants carrying only the frame-shifted mutagenesis were recovered in T1 generation following segregation. Whereas no obvious differences were observed between the sbeI mutants and wild type, sbeII mutants showed higher proportion of long chains presented in debranched amylopectin, significantly increased AC and RS content to as higher as 25.0% and 9.8%, respectively, and thus altered fine structure and nutritional properties of starch. Taken together, our results demonstrated for the first time the feasibility to create high-amylose rice through CRISPR/Cas9-mediated editing of SBEIIb.

Many plants increase in freezing tolerance in response to low temperature, a process known as cold acclimation. In Arabidopsis, cold acclimation involves action of the CBF cold response pathway. Key components of the pathway include rapid cold-induced expression of three homologous genes encoding transcriptional activators, CBF1, 2 and 3 (also known as DREB1b, c and a, respectively), followed by expression of CBF-targeted genes, the CBF regulon, that increase freezing tolerance. Unlike Arabidopsis, tomato cannot cold acclimate raising the question of whether it has a functional CBF cold response pathway. Here we show that tomato, like Arabidopsis, encodes three CBF homologs, LeCBF1-3 (Lycopersicon esculentum CBF1-3), that are present in tandem array in the genome. Only the tomato LeCBF1 gene, however, was found to be cold-inducible. As is the case for Arabidopsis CBF1-3, transcripts for LeCBF1-3 did accumulate in response to mechanical agitation, but not in response to drought, ABA or high salinity. Constitutive overexpression of LeCBF1 in transgenic Arabidopsis plants induced expression of CBF-targeted genes and increased freezing tolerance indicating that LeCBF1 encodes a functional homolog of the Arabidopsis CBF1-3 proteins. However, constitutive overexpression of either LeCBF1 or AtCBF3 in transgenic tomato plants did not increase freezing tolerance. Gene expression studies, including the use of a cDNA microarray representing approximately 8000 tomato genes, identified only four genes that were induced 2.5-fold or more in the LeCBF1 or AtCBF3 overexpressing plants, three of which were putative members of the tomato CBF regulon as they were also upregulated in response to low temperature. Additional experiments indicated that of eight tomato genes that were likely orthologs of Arabidopsis CBF regulon genes, none were responsive to CBF overexpression in tomato. From these results, we conclude that tomato has a complete CBF cold response pathway, but that the tomato CBF regulon differs from that of Arabidopsis and appears to be considerably smaller and less diverse in function.

Plants have acquired sophisticated stress response systems to adapt to changing environments. It is important to understand plants' stress response mechanisms in the effort to improve crop productivity under stressful conditions. The AP2/ERF transcription factors are known to regulate diverse processes of plant development and stress responses. In this study, the molecular characteristics and biological functions of AP2/ERFs in a variety of plant species were analyzed. AP2/ERFs, especially those in DREB and ERF subfamilies, are ideal candidates for crop improvement because their overexpression enhances tolerances to drought, salt, freezing, as well as resistances to multiple diseases in the transgenic plants. The comprehensive analysis of physiological functions is useful in elucidating the biological roles of AP2/ERF family genes in gene interaction, pathway regulation, and defense response under stress environments, which should provide new opportunities for the crop tolerance engineering.

BACKGROUND: Real-time quantitative reverse transcription PCR (RT-qPCR) data needs to be normalized for its proper interpretation. Housekeeping genes are routinely employed for this purpose, but their expression level cannot be assumed to remain constant under all possible experimental conditions. Thus, a systematic validation of reference genes is required to ensure proper normalization. For soybean, only a small number of validated reference genes are available to date. RESULTS: A systematic comparison of 14 potential reference genes for soybean is presented. These included seven commonly used (ACT2, ACT11, TUB4, TUA5, CYP, UBQ10, EF1b) and seven new candidates (SKIP16, MTP, PEPKR1, HDC, TIP41, UKN1, UKN2). Expression stability was examined by RT-qPCR across 116 biological samples, representing tissues at various developmental stages, varied photoperiodic treatments, and a range of soybean cultivars. Expression of all 14 genes was variable to some extent, but that of SKIP16, UKN1 and UKN2 was overall the most stable. A combination of ACT11, UKN1 and UKN2 would be appropriate as a reference panel for normalizing gene expression data among different tissues, whereas the combination SKIP16, UKN1 and MTP was most suitable for developmental stages. ACT11, TUA5 and TIP41 were the most stably expressed when the photoperiod was altered, and TIP41, UKN1 and UKN2 when the light quality was changed. For six different cultivars in long day (LD) and short day (SD), their expression stability did not vary significantly with ACT11, UKN2 and TUB4 being the most stable genes. The relative gene expression level of GmFTL3, an ortholog of Arabidopsis FT (FLOWERING LOCUS T) was detected to validate the reference genes selected in this study. CONCLUSION: None of the candidate reference genes was uniformly expressed across all experimental conditions, and the most suitable reference genes are conditional-, tissue-specific-, developmental-, and cultivar-dependent. Most of the new reference genes performed better than the conventional housekeeping genes. These results should guide the selection of reference genes for gene expression studies in soybean.