Genome-Wide Analysis of japonica Rice Performance under Limited Water and Permanent Flooding Conditions.

A rice GWAS panel of 281 accessions of japonica rice was phenotypically characterized for 26 traits related to phenology, plant and seed morphology, physiology and yield for 2 years in field conditions under permanent flooding (PF) and limited water (LW). A genome-wide analysis uncovered a total of 160 significant marker-trait associations (MTAs), of which 32 were LW-specific, 59 were PF-specific, and 69 were in common between the two water management systems. LW-specific associations were identified for several agronomic traits including days to maturation, days from flowering to maturation, leaf traits, plant height, panicle and seed traits, hundred grain weight, yield and tillering. Significant MTAs were detected across all the 12 rice chromosomes, while clusters of effects influencing different traits under LW or in both watering conditions were, respectively, observed on chromosomes 4, 8, and 12 and on chromosomes 1, 3, 4, 5, and 8. The analysis of genes annotated in the Nipponbare reference sequence and included in the regions associated to traits related to plant morphology, grain yield, and physiological parameters allowed the identification of genes that were demonstrated to affect the respective traits. Among these, three (OsOFP2, Dlf1, OsMADS56) and seven (SUI1, Sd1, OsCOL4, Nal1, OsphyB, GW5, Ehd1) candidate genes were, respectively, identified to co-localize with LW-specific associations and associations in common between the two water treatments. For several LW-specific MTAs, or in common among the two treatments, positional co-localizations with previously identified QTLs for rice adaptation to water shortages were observed, a result that further supports the role of the loci identified in this work in conferring adaptation to LW. The most robust associations identified here could represent suitable targets for genomic selection approaches to improve yield-related traits under LW.

Knockdown of MLO genes reduces susceptibility to powdery mildew in grapevine.

Erysiphe necator is the causal agent of powdery mildew (PM), one of the most destructive diseases of grapevine. PM is controlled by sulfur-based and synthetic fungicides, which every year are dispersed into the environment. This is why PM-resistant varieties should become a priority for sustainable grapevine and wine production. PM resistance can be achieved in other crops by knocking out susceptibility S-genes, such as those residing at genetic loci known as MLO (Mildew Locus O). All MLO S-genes of dicots belong to the phylogenetic clade V, including grapevine genes VvMLO7, 11 and 13, which are upregulated during PM infection, and VvMLO6, which is not upregulated. Before adopting a gene-editing approach to knockout candidate S-genes, the evidence that loss of function of MLO genes can reduce PM susceptibility is necessary. This paper reports the knockdown through RNA interference of VvMLO6, 7, 11 and 13. The knockdown of VvMLO6, 11 and 13 did not decrease PM severity, whereas the knockdown of VvMLO7 in combination with VvMLO6 and VvMLO11 reduced PM severity up to 77%. The knockdown of VvMLO7 and VvMLO6 seemed to be important for PM resistance, whereas a role for VvMLO11 does not seem likely. Cell wall appositions (papillae) were present in both resistant and susceptible lines in response to PM attack. Thirteen genes involved in defense were less upregulated in infected mlo plants, highlighting the early mlo-dependent disruption of PM invasion.

Genetic analysis of durable resistance to Magnaporthe oryzae in the rice accession Gigante Vercelli identified two blast resistance loci.

Rice cultivars exhibiting durable resistance to blast, the most important rice fungal disease provoking up to 30 % of rice losses, are very rare and searching for sources of such a resistance represents a priority for rice-breeding programs. To this aim we analyzed Gigante Vercelli (GV) and Vialone Nano (VN), two temperate japonica rice cultivars in Italy displaying contrasting response to blast, with GV showing a durable and broad-spectrum resistance, whereas VN being highly susceptible. An SSR-based genetic map developed using a GV × VN population segregating for blast resistance identified two blast resistance loci, localized to the long arm of chromosomes 1 and 4 explaining more than 78 % of the observed phenotypic variation for blast resistance. The pyramiding of two blast resistance QTLs was therefore involved in the observed durable resistance in GV. Mapping data were integrated with information obtained from RNA-seq expression profiling of all classes of resistance protein genes (resistance gene analogs, RGAs) and with the map position of known cloned or mapped blast resistance genes to search candidates for the GV resistant response. A co-localization of RGAs with the LOD peak or the marker interval of the chromosome 1 QTL was highlighted and a valuable tool for selecting the resistance gene during breeding programs was developed. Comparative analysis with known blast resistance genes revealed co-positional relationships between the chromosome 1 QTL with the Pi35/Pish blast resistance alleles and showed that the chromosome 4 QTL represents a newly identified blast resistance gene. The present genetic analysis has therefore allowed the identification of two blast resistance loci in the durable blast-resistant rice cultivar GV and tools for molecular selection of these resistance genes.

Deep sequencing transcriptional fingerprinting of rice kernels for dissecting grain quality traits.

BACKGROUND: Rice represents one the most important foods all over the world. In Europe, Italy is the first rice producer and Italian production is driven by tradition and quality. All main rice grain quality traits, like cooking properties, texture, gelatinization temperature, chalkiness and yield, are related to the content and composition of starch and seed-storage proteins in the endosperm and to grain shape. In addition, a number of nutraceutical compounds and allergens are known to have a significant effect on grain quality determination. To investigate the genetic bases underlying the qualitative differences that characterize traditional Italian rice cultivars, a comparative RNA-Seq-based transcriptomic analysis of developing caryopsis was conducted at 14 days after flowering on six popular Italian varieties (Carnaroli, Arborio, Balilla, Vialone Nano, Gigante Vercelli and Volano) phenotypically differing for qualitative grain-related traits. RESULTS: Co-regulation analyses of differentially expressed genes showing the same expression patterns in the six genotypes highlighted clusters of loci up or down-regulated in specific varieties, with respect to the others. Among them, we detected loci involved in cell wall biosynthesis, protein metabolism and redox homeostasis, classes of genes affecting in chalkiness determination. Moreover, loci encoding for seed-storage proteins, allergens or involved in the biosynthesis of specific nutraceutical compounds were also present and specifically regulated in the different clusters. A wider investigation of all the DEGs detected in pair-wise comparisons revealed transcriptional variation, among the six genotypes, for quality-related loci involved in starch biosynthesis (e.g. GBSSI, starch synthases and AGPase), genes encoding for transcription factors, additional seed storage proteins, allergens or belonging to additional nutraceutical compounds biosynthetic pathways and loci affecting grain size. Putative functional SNPs associated to amylose content in starch, gelatinization temperature and grain size were also identified. CONCLUSIONS: The present work represents a more extended phenotypic characterization of a set of rice accessions that present a wider genetic variability than described nowadays in literature. The results provide the first transcriptional picture for several of the grain quality differences observed among the Italian rice varieties analyzed and reveal that each variety is characterized by the over-expression of a peculiar set of loci affecting grain appearance and quality. A list of candidates and SNPs affecting specific grain properties has been identified offering a starting point for further works aimed to characterize genes and molecular markers for breeding programs.

Improvement of marker-based predictability of Apparent Amylose Content in japonica rice through GBSSI allele mining.

BACKGROUND: Apparent Amylose Content (AAC), regulated by the Waxy gene, represents the key determinant of rice cooking properties. In occidental countries high AAC rice represents the most requested market class but the availability of molecular markers allowing specific selection of high AAC varieties is limited. RESULTS: In this study, the effectiveness of available molecular markers in predicting AAC was evaluated in a collection of 127 rice accessions (125 japonica ssp. and 2 indica ssp.) characterized by AAC values from glutinous to 26%. The analyses highlighted the presence of several different allelic patterns identifiable by a few molecular markers, and two of them, i.e., the SNPs at intron1 and exon 6, were able to explain a maximum of 79.5% of AAC variation. However, the available molecular markers haplotypes did not provide tools for predicting accessions with AAC higher than 24.5%. To identify additional polymorphisms, the re-sequencing of the Waxy gene and 1kbp of the putative upstream regulatory region was performed in 21 genotypes representing all the AAC classes identified. Several previously un-characterized SNPs were identified and four of them were used to develop dCAPS markers. CONCLUSIONS: The addition of the SNPs newly identified slightly increased the AAC explained variation and allowed the identification of a haplotype almost unequivocally associated to AAC higher than 24.5%. Haplotypes at the waxy locus were also associated to grain length and length/width (L/W) ratio. In particular, the SNP at the first intron, which identifies the Wxa and Wxb alleles, was associated with differences in the width of the grain, the L/W ratio and the length of the kernel, most likely as a result of human selection.

Haplotype variability and identification of new functional alleles at the Rdg2a leaf stripe resistance gene locus.

The barley Rdg2a locus confers resistance to the leaf stripe pathogen Pyrenophora graminea and, in the barley genotype Thibaut, it is composed of a gene family with three highly similar paralogs. Only one member of the gene family (called as Rdg2a) encoding for a CC-NB-LRR protein is able to confer resistance to the leaf stripe isolate Dg2. To study the genome evolution and diversity at the Rdg2a locus, sequences spanning the Rdg2a gene were compared in two barley cultivars, Thibaut and Morex, respectively, resistant and susceptible to leaf stripe. An overall high level of sequence conservation interrupted by several rearrangements that included three main deletions was observed in the Morex contig. The main deletion of 13,692 bp was most likely derived from unequal crossing over between Rdg2a paralogs leading to the generation of a chimeric Morex rdg2a gene which was not associated to detectable level of resistance toward leaf stripe. PCR-based analyses of genic and intergenic regions at the Rdg2a locus in 29 H. vulgare lines and one H. vulgare ssp. spontaneum accession indicated large haplotype variability in the cultivated barley gene pool suggesting rapid and recent divergence at this locus. Barley genotypes showing the same haplotype as Thibaut at the Rdg2a locus were selected for a Rdg2a allele mining through allele re-sequencing and two lines with polymorphic nucleotides leading to amino acid changes in the CC-NB and LRR encoding domains, respectively, were identified. Analysis of nucleotide diversity of the Rdg2a alleles revealed that the polymorphic sites were subjected to positive selection. Moreover, strong positively selected sites were located in the LRR encoding domain suggesting that both positive selection and divergence at homologous loci are possibly representing the molecular mechanism for the generation of high diversity at the Rdg2a locus in the barley gene pool.

Comparative transcriptome profiling of the early response to Magnaporthe oryzae in durable resistant vs susceptible rice (Oryza sativa L.) genotypes.

Durable resistance to blast, the most significant fungal disease of rice, represents an agronomically relevant character. Gigante Vercelli (GV) and Vialone Nano (VN) are two old temperate japonica Italian rice cultivars with contrasting response to blast infection: GV displays durable and broad resistance while VN is highly susceptible. RNA-seq was used to dissect the early molecular processes deployed during the resistance response of GV at 24 h after blast inoculation. Differential gene expression analysis identified 1,070 and 1,484 modulated genes, of which 726 and 699 were up regulated in response to infection in GV and VN, respectively. Gene ontology (GO) enrichment analyses revealed a set of GO terms enriched in both varieties but, despite this commonality, the gene sets contributing to common GO enriched terms were dissimilar. The expression patterns of genes grouped in GV-specific enriched GO terms were examined in detail including at the transcript isoform level. GV exhibited a dramatic up-regulation of genes encoding diterpene phytoalexin biosynthetic enzymes, flavin-containing monooxygenase, class I chitinase and glycosyl hydrolase 17. The sensitivity and high dynamic range of RNA-seq allowed the identification of genes critically involved in conferring GV resistance during the early steps of defence perception-signalling. These included chitin oligosaccharides sensing factors, wall associated kinases, MAPK cascades and WRKY transcription factors. Candidate genes with expression patterns consistent with a potential role as GV-specific functional resistance (R) gene(s) were also identified. This first application of RNA-seq to dissect durable blast resistance supports a crucial role of the prompt induction of a battery of responses including defence-related genes as well as members of gene families involved in signalling and pathogen-related gene expression regulation.

Emerging knowledge from genome sequencing of crop species.

Extensive insights into the genome composition, organization, and evolution have been gained from the plant genome sequencing and annotation ongoing projects. The analysis of crop genomes provided surprising evidences with important implications in plant origin and evolution: genome duplication, ancestral re-arrangements and unexpected polyploidization events opened new doors to address fundamental questions related to species proliferation, adaptation, and functional modulations. Detailed paleogenomic analysis led to many speculation on how chromosomes have been shaped over time in terms of gene content and order. The completion of the genome sequences of several major crops, prompted to a detailed identification and annotation of transposable elements: new hypothesis related to their composition, chromosomal distribution, insertion models, amplification rate, and evolution patterns are coming up. Availability of full genome sequence of several crop species as well as from many accessions within species is providing new keys for biodiversity exploitation and interpretation. Re-sequencing is enabling high-throughput genotyping to identify a wealth of SNP and afterward to produce haplotype maps necessary to accurately associate molecular variation to phenotype. Conservation genomics is emerging as a powerful tool to explain adaptation, genetic drift, natural selection, hybridization and to estimate genetic variation, fitness and population's viability.

Identification and mapping of the leaf stripe resistance gene Rdg1a in Hordeum spontaneum.

Leaf stripe of barley, caused by Pyrenophora graminea, is an important seed-borne disease in organically grown as well as in conventionally grown Nordic and Mediterranean barley districts. Two barley segregating populations represented by 103 recombinant inbred lines (RILs) of the cross L94 (susceptible) x Vada (resistant) and 194 RILs of the cross Arta (susceptible) x Hordeum spontaneum 41-1 (resistant) were analysed with two highly virulent leaf stripe isolates, Dg2 and Dg5, to identify loci for P. graminea resistance. A major gene with its positive allele contributed by Vada and H. spontaneum 41-1 was detected in both populations and for both pathogen isolates on chromosome 2HL explaining 44.1 and 91.8% R (2), respectively for Dg2 and Dg5 in L94 x Vada and 97.8 and 96.1% R (2), respectively for Dg2 and Dg5 in Arta x H. spontaneum 41-1. Common markers in the gene region of the two populations enabled map comparison and highlighted an overlapping for the region of the resistance locus. Since the map position of the resistance locus identified in this report is the same as that for the leaf stripe resistance gene Rdg1a, mapped earlier in Alf and derived from the 'botanical' barley line H. laevigatum, we propose that leaf stripe resistance in Vada and H. spontaneum 41-1 is governed by the same gene, namely by Rdg1a, and that Rdg1a resistance could be traced back to H. spontaneum, the progenitor of cultivated barley. PCR-based molecular markers that can be used for marker-assisted selection (MAS) of Rdg1a were identified. An Rdg1a syntenic interval with the rice chromosome arm 4L was identified on the basis of rice orthologs of EST-based barley markers. Analysis of the rice genes annotated into the syntenic interval did not reveal sequences strictly belonging to the major class (nucleotide-binding site plus leucine-rich repeat) of the resistance genes. Nonetheless, four genes coding for domains that are present in the major disease-resistance genes, namely receptor-like protein kinase and ATP/GTP-binding proteins, were identified together with a homolog of the barley powdery mildew resistance gene mlo. Three (out of five) homologs of these genes were mapped in the Rdg1a region in barley and the mlo homolog map position was tightly associated with the LOD score peak in both populations.