Identification of candidate flowering and sex genes in white Guinea yam (D. rotundata Poir.) by SuperSAGE transcriptome profiling.

Dioecy (distinct male and female individuals) and scarce to non-flowering are common features of cultivated yam (Dioscorea spp.). However, the molecular mechanisms underlying flowering and sex determination in Dioscorea are largely unknown. We conducted SuperSAGE transcriptome profiling of male, female and monoecious individuals to identify flowering and sex-related genes in white Guinea yam (D. rotundata), generating 20,236 unique tags. Of these, 13,901 were represented by a minimum of 10 tags. A total 88 tags were significantly differentially expressed in male, female and monoecious plants, of which 18 corresponded to genes previously implicated in flower development and sex determination in multiple plant species. We validated the SuperSAGE data with quantitative real-time PCR (qRT-PCR)-based analysis of the expression of three candidate genes. We further investigated the flowering patterns of 1938 D. rotundata accessions representing diverse geographical origins over two consecutive years. Over 85% of accessions were either male or non-flowering, less than 15% were female, while monoecious plants were rare. Intensity of flowering varied between male and female plants, with the former flowering more abundantly than the latter. Candidate genes identified in this study can be targeted for further validation and to induce regular flowering in poor to non-flowering cultivars. Findings of the study provide important inputs for further studies aiming to overcome the challenge of flowering in yams and to improve efficiency of yam breeding.

Next-generation sequencing-based bulked segregant analysis for QTL mapping in the heterozygous species Brassica rapa.

KEY MESSAGE: An improved protocol of QTL-seq, an NGS-based method for bulked segregant analysis we previously developed in rice, allowed successful mapping of QTLs of interest in the highly heterozygous genome of B. rapa, demonstrating the power of this elegant method for genetic analyses in heterozygous species of economic importance. Recent advances in next-generation sequencing (NGS) and the various NGS-based methods developed for rapidly identifying candidate genes of interest have accelerated genetic analysis mainly in the model plants rice and Arabidopsis. Brassica rapa includes several economically important crops such as Chinese cabbage, turnip and various leafy vegetables. The application of NGS-based approaches for the analysis of B. rapa has been limited mainly due to its highly heterozygous genome and poor quality of the reference genome sequence currently available for this species. In this study, we have improved QTL-seq, a method for NGS-based bulked segregant analysis we previously developed in rice, extending its applicability for accelerating the genetic analysis and molecular breeding of B. rapa. Addition of new filters to the original QTL-seq pipeline allowed removal of spurious single-nucleotide polymorphisms caused by alignment/sequencing errors and variability between parents, significantly improving accuracy of the analysis. As proof of principle, we successfully applied the new approach to identify candidate genomic regions controlling flowering and trichome formation using segregating F2 progeny obtained from crosses made between cultivars of B. rapa showing contrasting phenotypes for these traits. We strongly believe that the improved QTL-seq method reported here will extend the applicability of NGS-based genetic analysis not only to B. rapa but also to other plant species of economic importance with heterozygous genomes.

Strigolactone perception and deactivation by a hydrolase receptor DWARF14.

The perception mechanism for the strigolactone (SL) class of plant hormones has been a subject of debate because their receptor, DWARF14 (D14), is an α/ß-hydrolase that can cleave SLs. Here we show via time-course analyses of SL binding and hydrolysis by Arabidopsis thaliana D14, that the level of uncleaved SL strongly correlates with the induction of the active signaling state. In addition, we show that an AtD14D218A catalytic mutant that lacks enzymatic activity is still able to complement the atd14 mutant phenotype in an SL-dependent manner. We conclude that the intact SL molecules trigger the D14 active signaling state, and we also describe that D14 deactivates bioactive SLs by the hydrolytic degradation after signal transmission. Together, these results reveal that D14 is a dual-functional receptor, responsible for both the perception and deactivation of bioactive SLs.

AgriSeqDB: an online RNA-Seq database for functional studies of agriculturally relevant plant species.

BACKGROUND: The genome-wide expression profile of genes in different tissues/cell types and developmental stages is a vital component of many functional genomic studies. Transcriptome data obtained by RNA-sequencing (RNA-Seq) is often deposited in public databases that are made available via data portals. Data visualization is one of the first steps in assessment and hypothesis generation. However, these databases do not typically include visualization tools and establishing one is not trivial for users who are not computational experts. This, as well as the various formats in which data is commonly deposited, makes the processes of data access, sharing and utility more difficult. Our goal was to provide a simple and user-friendly repository that meets these needs for data-sets from major agricultural crops. DESCRIPTION: AgriSeqDB ( https://expression.latrobe.edu.au/agriseqdb ) is a database for viewing, analysing and interpreting developmental and tissue/cell-specific transcriptome data from several species, including major agricultural crops such as wheat, rice, maize, barley and tomato. The disparate manner in which public transcriptome data is often warehoused and the challenge of visualizing raw data are both major hurdles to data reuse. The popular eFP browser does an excellent job of presenting transcriptome data in an easily interpretable view, but previous implementation has been mostly on a case-by-case basis. Here we present an integrated visualisation database of transcriptome data-sets from six species that did not previously have public-facing visualisations. We combine the eFP browser, for gene-by-gene investigation, with the Degust browser, which enables visualisation of all transcripts across multiple samples. The two visualisation interfaces launch from the same point, enabling users to easily switch between analysis modes. The tools allow users, even those without bioinformatics expertise, to mine into data-sets and understand the behaviour of transcripts of interest across samples and time. We have also incorporated an additional graphic download option to simplify incorporation into presentations or publications. CONCLUSION: Powered by eFP and Degust browsers, AgriSeqDB is a quick and easy-to-use platform for data analysis and visualization in five crops and Arabidopsis. Furthermore, it provides a tool that makes it easy for researchers to share their data-sets, promoting research collaborations and data-set reuse.

Regulation of genome-wide DNA methylation by mobile small RNAs.

Contents Summary 540 I. Introduction 540 II. There are different types of sRNA mobility 541 III. Mechanisms of sRNA movement 541 IV. Long-distance, shoot-root, mobile siRNAs influence DNA methylation in recipient tissues 541 V. Classes of interactions between shoot-root mobile siRNAs and DNA methylation 542 VI. Loci targeted directly and indirectly by shoot-root mobile siRNAs are associated with different histone modifications 543 VII. Is mobile siRNA-regulated DNA methylation important in specific tissues or under specific conditions? 543 VIII. Mobile sRNAs can be used to modify plant traits 544 IX. Conclusions 544 Acknowledgements 544 References 544 SUMMARY: RNA-directed DNA methylation (RdDM) at cytosine residues regulates gene expression, silences transposable elements and influences genome stability. The mechanisms responsible for RdDM are guided to target loci by small RNAs (sRNAs) that can move within plants cell to cell and long distance. Here we discuss recent advances in the understanding of interactions between mobile sRNAs and DNA methylation. We describe the mechanisms of sRNA movement, the differences between known classes of mobile sRNA-DNA methylation interactions and the limits of current knowledge. Finally, we discuss potential applications of mobile sRNAs in modifying plant traits.

Genome sequencing of the staple food crop white Guinea yam enables the development of a molecular marker for sex determination.

BACKGROUND: Root and tuber crops are a major food source in tropical Africa. Among these crops are several species in the monocotyledonous genus Dioscorea collectively known as yam, a staple tuber crop that contributes enormously to the subsistence and socio-cultural lives of millions of people, principally in West and Central Africa. Yam cultivation is constrained by several factors, and yam can be considered a neglected "orphan" crop that would benefit from crop improvement efforts. However, the lack of genetic and genomic tools has impeded the improvement of this staple crop. RESULTS: To accelerate marker-assisted breeding of yam, we performed genome analysis of white Guinea yam (Dioscorea rotundata) and assembled a 594-Mb genome, 76.4% of which was distributed among 21 linkage groups. In total, we predicted 26,198 genes. Phylogenetic analyses with 2381 conserved genes revealed that Dioscorea is a unique lineage of monocotyledons distinct from the Poales (rice), Arecales (palm), and Zingiberales (banana). The entire Dioscorea genus is characterized by the occurrence of separate male and female plants (dioecy), a feature that has limited efficient yam breeding. To infer the genetics of sex determination, we performed whole-genome resequencing of bulked segregants (quantitative trait locus sequencing [QTL-seq]) in F1 progeny segregating for male and female plants and identified a genomic region associated with female heterogametic (male = ZZ, female = ZW) sex determination. We further delineated the W locus and used it to develop a molecular marker for sex identification of Guinea yam plants at the seedling stage. CONCLUSIONS: Guinea yam belongs to a unique and highly differentiated clade of monocotyledons. The genome analyses and sex-linked marker development performed in this study should greatly accelerate marker-assisted breeding of Guinea yam. In addition, our QTL-seq approach can be utilized in genetic studies of other outcrossing crops and organisms with highly heterozygous genomes. Genomic analysis of orphan crops such as yam promotes efforts to improve food security and the sustainability of tropical agriculture.

A chloroplast-localized protein LESION AND LAMINA BENDING affects defence and growth responses in rice.

Understanding how plants allocate their resources to growth or defence is of long-term importance to the development of new and improved varieties of different crops. Using molecular genetics, plant physiology, hormone analysis and Next-Generation Sequencing (NGS)-based transcript profiling, we have isolated and characterized the rice (Oryza sativa) LESION AND LAMINA BENDING (LLB) gene that encodes a chloroplast-targeted putative leucine carboxyl methyltransferase. Loss of LLB function results in reduced growth and yield, hypersensitive response (HR)-like lesions, accumulation of the antimicrobial compounds momilactones and phytocassanes, and constitutive expression of pathogenesis-related genes. Consistent with these defence-associated responses, llb shows enhanced resistance to rice blast (Magnaporthe oryzae) and bacterial blight (Xanthomonas oryzae pv. oryzae). The lesion and resistance phenotypes are likely to be caused by the over-accumulation of jasmonates (JAs) in the llb mutant including the JA precursor 12-oxo-phytodienoic acid. Additionally, llb shows an increased lamina inclination and enhanced early seedling growth due to elevated brassinosteroid (BR) synthesis and/or signalling. These findings show that LLB functions in the chloroplast to either directly or indirectly repress both JA- and BR-mediated responses, revealing a possible mechanism for controlling how plants allocate resources for defence and growth.

WIPPER: an accurate and efficient field phenotyping platform for large-scale applications.

More accurate, rapid, and easy phenotyping tools are required to match the recent advances in high-throughput genotyping for accelerating breeding and genetic analysis. The conventional data recording in field notebooks and then inputting data to computers for further analysis is inefficient, time-consuming, laborious, and prone to human error. Here, we report WIPPER (for Wireless Plant Phenotyper), a new phenotyping platform that combines field phenotyping and data recording with the aid of Bluetooth communication, thus saving time and labor not only for field data recoding but also for inputting data to computers. Additionally, it eliminates the risk of human error associated with phenotyping and inputting data. We applied WIPPER to 100 individuals of a rice recombinant inbred line (RIL) for measuring leaf width and relative chlorophyll content (SPAD value), and were able to record an accurate data in a significantly reduced time compared with the conventional method of data collection. We are currently using WIPPER for routine management of rice germplasm including recording and documenting information on phenotypic data, seeds, and DNA for their accelerated utilization in crop breeding.

A cytochrome P450, OsDSS1, is involved in growth and drought stress responses in rice (Oryza sativa L.).

Cytochrome P450s are among the largest protein coding gene families in plant genomes. However, majority of the genes remain uncharacterized. Here, we report the characterization of dss1, a rice mutant showing dwarfism and reduced grain size. The dss1 phenotype is caused by a non-synonymous point mutation we identified in DSS1, which is member of a P450 gene cluster located on rice chromosome 3 and corresponds to the previously reported CYP96B4/SD37 gene. Phenotypes of several dwarf mutants characterized in rice are associated with defects in the biosynthesis or perception of the phytohormones gibberellins (GAs) and brassinosteroids (BRs). However, both GA and BR failed to rescue the dss1 phenotype. Hormone profiling revealed the accumulation of abscisic acid (ABA) and ABA metabolites, as well as significant reductions in GA19 and GA53 levels, precursors of the bioactive GA1, in the mutant. The dss1 contents of cytokinin and auxins were not significantly different from wild-type plants. Consistent with the accumulation of ABA and metabolites, germination and early growth was delayed in dss1, which also exhibited an enhanced tolerance to drought. Additionally, expressions of members of the DSS1/CYP96B gene cluster were regulated by drought stress and exogenous ABA. RNA-seq-based transcriptome profiling revealed, among others, that cell wall-related genes and genes involved in lipid metabolism were up- and down-regulated in dss1, respectively. Taken together, these findings suggest that DSS1 mediates growth and stress responses in rice by fine-tuning GA-to-ABA balance, and might as well play a role in lipid metabolism.

The rice (Oryza sativa L.) LESION MIMIC RESEMBLING, which encodes an AAA-type ATPase, is implicated in defense response.

Lesion mimic mutants (LMMs) provide a useful tool to study defense-related programmed cell death (PCD) in plants. Although a number of LMMs have been identified in multiple species, most of the candidate genes are yet to be isolated. Here, we report the identification and characterization of a novel rice (Oryza sativa L.) lesion mimic resembling (lmr) mutant, and cloning of the corresponding LMR gene. The LMR locus was initially delineated to 1.2 Mb region on chromosome 6, which was further narrowed down to 155-kb using insertions/deletions (INDELs) and cleavage amplified polymorphic sequence markers developed in this study. We sequenced the open reading frames predicted within the candidate genomic region, and identified a G-A base substitution causing a premature translation termination in a gene that encodes an ATPase associated with various cellular activities type (AAA-type) protein. RNA interference transgenic lines with reduced LMR transcripts exhibited the lesion mimic phenotype similar to that of lmr plants. Furthermore, expression of the wild-type LMR in the mutant background complemented the lesion phenotype, confirming that the mutation identified in LMR is responsible for the mutant phenotype. The pathogenesis-related (PR) genes PBZ1 and PR1 were induced in lmr, which also showed enhanced resistance to rice blast (Magnaporthe oryzae) and bacterial blight (Xanthomonas oryzae pv. oryzae), suggesting LMR is a negative regulator of cell death in rice. The identification of lmr and cloning of the corresponding LMR gene provide an additional resource for the study of PCD in plants.

Protoplast cell death assay to study Magnaporthe oryzae AVR gene function in rice.

We describe a protocol for transient gene expression in rice protoplasts and its application to the study of Magnaporthe oryzae avirulence (AVR) gene function. In this assay the gene encoding the firefly luciferase protein is transfected into rice protoplasts by electroporation together with the candidate AVR genes. The luminescence can then be used to assess the viability of rice protoplasts. The hypersensitive response (HR) caused by the interaction between M. oryzae AVR and rice R genes can subsequently be monitored by recording the decrease in luminescence from the transfected cells.

The tillering phenotype of the rice plastid terminal oxidase (PTOX) loss-of-function mutant is associated with strigolactone deficiency.

The significance of plastid terminal oxidase (PTOX) in phytoene desaturation and chloroplast function has been demonstrated using PTOX-deficient mutants, particularly in Arabidopsis. However, studies on its role in monocots are lacking. Here, we report cloning and characterization of the rice (Oryza sativa) PTOX1 gene. Using Ecotype Targeting Induced Local Lesions IN Genomes (EcoTILLING) and TILLING as forward genetic tools, we identified the causative mutation of an EMS mutant characterized by excessive tillering, semi-dwarfism and leaf variegation that corresponded to the PTOX1 gene. The tillering and semi-dwarf phenotypes of the ptox1 mutant are similar to phenotypes of known strigolactone (SL)-related rice mutants, and both phenotypic traits could be rescued by application of the synthetic SL GR24. The ptox1 mutant accumulated phytoene in white leaf sectors with a corresponding deficiency in ß-carotene, consistent with the expected function of PTOX1 in promoting phytoene desaturase activity. There was also no accumulation of the carotenoid-derived SL ent-2'-epi-5-deoxystrigol in root exudates. Elevated concentrations of auxin were detected in the mutant, supporting previous observations that SL interaction with auxin is important in shoot branching control. Our results demonstrate that PTOX1 is required for both carotenoid and SL synthesis resulting in SL-deficient phenotypes in rice.

MutMap+: genetic mapping and mutant identification without crossing in rice.

Advances in genome sequencing technologies have enabled researchers and breeders to rapidly associate phenotypic variation to genome sequence differences. We recently took advantage of next-generation sequencing technology to develop MutMap, a method that allows rapid identification of causal nucleotide changes of rice mutants by whole genome resequencing of pooled DNA of mutant F2 progeny derived from crosses made between candidate mutants and the parental line. Here we describe MutMap+, a versatile extension of MutMap, that identifies causal mutations by comparing SNP frequencies of bulked DNA of mutant and wild-type progeny of M3 generation derived from selfing of an M2 heterozygous individual. Notably, MutMap+ does not necessitate artificial crossing between mutants and the wild-type parental line. This method is therefore suitable for identifying mutations that cause early development lethality, sterility, or generally hamper crossing. Furthermore, MutMap+ is potentially useful for gene isolation in crops that are recalcitrant to artificial crosses.

MutMap-Gap: whole-genome resequencing of mutant F2 progeny bulk combined with de novo assembly of gap regions identifies the rice blast resistance gene Pii.

Next-generation sequencing allows the identification of mutations responsible for mutant phenotypes by whole-genome resequencing and alignment to a reference genome. However, when the resequenced cultivar/line displays significant structural variation from the reference genome, mutations in the genome regions missing from the reference (gaps) cannot be identified by simple alignment. Here we report on a method called 'MutMap-Gap', which involves delineating a candidate region harboring a mutation of interest using the recently reported MutMap method, followed by de novo assembly, alignment, and identification of the mutation within genome gaps. We applied MutMap-Gap to isolate the blast resistant gene Pii from the rice cv Hitomebore using mutant lines that have lost Pii function. MutMap-Gap should prove useful for cloning genes that exhibit significant structural variations such as disease resistance genes of the nucleotide-binding site-leucine rich repeat (NBS-LRR) class.

QTL-seq: rapid mapping of quantitative trait loci in rice by whole genome resequencing of DNA from two bulked populations.

The majority of agronomically important crop traits are quantitative, meaning that they are controlled by multiple genes each with a small effect (quantitative trait loci, QTLs). Mapping and isolation of QTLs is important for efficient crop breeding by marker-assisted selection (MAS) and for a better understanding of the molecular mechanisms underlying the traits. However, since it requires the development and selection of DNA markers for linkage analysis, QTL analysis has been time-consuming and labor-intensive. Here we report the rapid identification of plant QTLs by whole-genome resequencing of DNAs from two populations each composed of 20-50 individuals showing extreme opposite trait values for a given phenotype in a segregating progeny. We propose to name this approach QTL-seq as applied to plant species. We applied QTL-seq to rice recombinant inbred lines and F2 populations and successfully identified QTLs for important agronomic traits, such as partial resistance to the fungal rice blast disease and seedling vigor. Simulation study showed that QTL-seq is able to detect QTLs over wide ranges of experimental variables, and the method can be generally applied in population genomics studies to rapidly identify genomic regions that underwent artificial or natural selective sweeps.

Mutation in OsLMS, a gene encoding a protein with two double-stranded RNA binding motifs, causes lesion mimic phenotype and early senescence in rice (Oryza sativa L.).

The rice (Oryza sativa L.) lesion mimic and senescence (lms) EMS-mutant, identified in a japonica cultivar Hitomebore, is characterized by a spontaneous lesion mimic phenotype during its vegetative growth, an accelerated senescence after flowering, and enhanced resistance to rice blast (Magnaporthe oryzae). To isolate the OsLMS gene, we crossed the lms mutant to Kasalath (indica), and used mutant F(2) plants to initially map the candidate region to about 322-kb on the long arm of chromosome 2. Illumina whole-genome re-sequencing of the mutant and aligning the reads to Hitomebore reference sequence within the candidate region delineated by linkage analysis identified a G to A nucleotide substitution. The mutation corresponded to the exon-intron splicing junction of a novel gene that encodes a carboxyl-terminal domain (CTD) phosphatase domain and two double stranded RNA binding motifs (dsRBM) containing protein. By PCR amplification, we confirmed that the mutation causes splicing error that is predicted to introduce a premature stop codon. RNA interference (RNAi) transgenic lines with suppressed expression of LMS gene exhibited the lesion mimic phenotype, confirming that the mutation identified in LMS is responsible for the mutant phenotype. OsLMS shares a moderate amino-acid similarity to the Arabidopsis FIERY2/CPL1 gene, which is known to control many plant processes such as stress response and development. Consistence with this similarity, the lms mutant shows sensitivity to cold stress at the early growth stage, suggesting that LMS is a negative regulator of stress response in rice.

Genome sequencing reveals agronomically important loci in rice using MutMap.

The majority of agronomic traits are controlled by multiple genes that cause minor phenotypic effects, making the identification of these genes difficult. Here we introduce MutMap, a method based on whole-genome resequencing of pooled DNA from a segregating population of plants that show a useful phenotype. In MutMap, a mutant is crossed directly to the original wild-type line and then selfed, allowing unequivocal segregation in second filial generation (F(2)) progeny of subtle phenotypic differences. This approach is particularly amenable to crop species because it minimizes the number of genetic crosses (n = 1 or 0) and mutant F(2) progeny that are required. We applied MutMap to seven mutants of a Japanese elite rice cultivar and identified the unique genomic positions most probable to harbor mutations causing pale green leaves and semidwarfism, an agronomically relevant trait. These results show that MutMap can accelerate the genetic improvement of rice and other crop plants.

A multifaceted genomics approach allows the isolation of the rice Pia-blast resistance gene consisting of two adjacent NBS-LRR protein genes.

The Oryza sativa (rice) resistance gene Pia confers resistance to the blast fungus Magnaporthe oryzae carrying the AVR-Pia avirulence gene. To clone Pia, we employed a multifaceted genomics approach. First, we selected 12 R-gene analog (RGA) genes encoding nucleotide binding site-leucine rich repeats (NBS-LRRs) proteins from a region on chromosome 11 that shows linkage to Pia. By using seven rice accessions, we examined the association between Pia phenotypes and DNA polymorphisms in the 10 genes, which revealed three genes (Os11gRGA3-Os11gRGA5) exhibiting a perfect association with the Pia phenotypes. We also screened ethyl methane sulfonate (EMS)-treated mutant lines of the rice cultivar 'Sasanishiki' harboring Pia, and isolated two mutants that lost the Pia phenotype. DNA sequencing of Os11gRGA3-Os11gRGA5 from the two mutant lines identified independent mutations of major effects in Os11gRGA4. The wild-type 'Sasanishiki' allele of Os11gRGA4 (SasRGA4) complemented Pia function in both mutants, suggesting that SasRGA4 is necessary for Pia function. However, when the rice cultivar 'Himenomochi' lacking Pia was transfected with SasRGA4, the Pia phenotype was not recovered. An additional complementation study revealed that the two NBS-LRR-type R genes, SasRGA4 and SasRGA5, that are located next to each other and oriented in the opposite direction are necessary for Pia function. A population genetics analysis of SasRGA4 and SasRGA5 suggests that the two genes are under long-term balancing selection.