Navigating rice seedling cold resilience: QTL mapping in two inbred line populations and the search for genes.

Due to global climate change resulting in extreme temperature fluctuations, it becomes increasingly necessary to explore the natural genetic variation in model crops such as rice to facilitate the breeding of climate-resilient cultivars. To uncover genomic regions in rice involved in managing cold stress tolerance responses and to identify associated cold tolerance genes, two inbred line populations developed from crosses between cold-tolerant and cold-sensitive parents were used for quantitative trait locus (QTL) mapping of two traits: degree of membrane damage after 1 week of cold exposure quantified as percent electrolyte leakage (EL) and percent low-temperature seedling survivability (LTSS) after 1 week of recovery growth. This revealed four EL QTL and 12 LTSS QTL, all overlapping with larger QTL regions previously uncovered by genome-wide association study (GWAS) mapping approaches. Within the QTL regions, 25 cold-tolerant candidate genes were identified based on genomic differences between the cold-tolerant and cold-sensitive parents. Of those genes, 20% coded for receptor-like kinases potentially involved in signal transduction of cold tolerance responses; 16% coded for transcription factors or factors potentially involved in regulating cold tolerance response effector genes; and 64% coded for protein chaperons or enzymes potentially serving as cold tolerance effector proteins. Most of the 25 genes were cold temperature regulated and had deleterious nucleotide variants in the cold-sensitive parent, which might contribute to its cold-sensitive phenotype.

An Allelic Variant of the Broad-Spectrum Blast Resistance Gene Ptr in Weedy Rice Is Associated with Resistance to the Most Virulent Blast Race IB-33.

Rice resistance (R) genes have been effectively deployed to prevent blast disease caused by the fungal pathogen Magnaporthe oryzae, one of the most serious threats for stable rice production worldwide. Weedy rice competing with cultivated rice may carry novel or lost R genes. The quantitative trait locus qBR12.3b was previously mapped between two single nucleotide polymorphism markers at the 10,633,942-bp and 10,820,033-bp genomic positions in a black-hull-awned (BHA) weed strain using a weed-crop-mapping population under greenhouse conditions. In this study, we found a portion of the known resistance gene Ptr encoding a protein with four armadillo repeats and confers a broad spectrum of blast resistance. We then analyzed the sequences of the Ptr gene from weedy rice, PtrBHA, and identified a unique amino acid glutamine at protein position 874. Minor changes of protein conformation of the PtrBHA gene were predicted through structural analysis of PtrBHA, suggesting that the product of PtrBHA is involved in disease resistance. A gene-specific codominant marker HJ17-13 from PtrBHA was then developed to distinguish alleles in weeds and crops. The PtrBHA gene existed in 207 individuals of the same mapping population, where qBR12.3b was mapped using this gene-specific marker. Disease reactions of 207 individuals and their parents to IB-33 were evaluated. The resistant individuals had PtrBHA whereas the susceptible individuals did not, suggesting that HJ17-13 is reliable to predict qBR12.3b. Taken together, this newly developed marker, and weedy rice genotypes carrying qBR12.3b, are useful for blast improvement using marker assisted selection.

Mapping Blast Resistance Genes in Rice Varieties 'Minghui 63' and 'M-202'.

Rice blast caused by the fungus Magnaporthe oryzae (syn. Magnaporthe grisea) is one of the most lethal diseases for sustainable rice production worldwide. Blast resistance mediated by major resistance genes is often broken down after a short period of deployment, while minor blast resistance genes, each providing a small effect on disease reactions, are more durable. In the present study, we first evaluated disease reactions of two rice breeding parents 'Minghui 63' and 'M-202' with 11 blast races, IA45, IB1, IB45, IB49, IB54, IC1, IC17, ID1, IE1, IG1, and IH1, commonly present in the United States, under greenhouse conditions using a category disease rating resembling infection types under field conditions. 'Minghui 63' exhibited differential resistance responses in comparison with those of 'M-202' to the tested blast races. A recombinant inbred line (RIL) population of 275 lines from a cross between 'Minghui 63' and 'M-202' was also evaluated with the above-mentioned blast races. The population was genotyped with 156 simple sequence repeat (SSR) and insertion and deletion (Indel) markers. A linkage map with a genetic distance of 1,022.84 cM was constructed using inclusive composite interval mapping (ICIM) software. A total of 10 resistance QTLs, eight from 'Minghui 63' and two from 'M-202', were identified. One major QTL, qBLAST2 on chromosome 2, was identified by seven races/isolates. The remaining nine minor resistance QTLs were mapped on chromosomes 1, 3, 6, 9, 10, 11, and 12. These findings provide useful genetic markers and resources to tag minor blast resistance genes for marker-assisted selection in rice breeding program and for further studies of underlying genes.

Phylogenetic origin and dispersal pattern of Taiwan weedy rice.

BACKGROUND: Knowledge of the genetic diversity and spatial structure of Taiwan weedy red rice (WRR) populations, which adapted in a transplanting system, will facilitate the design of effective methods to control this weed by tracing its origins and dispersal patterns in a given region. RESULTS: Taiwan WRR is genetically most similar to Taiwan indica cultivars and landraces according to genetic distance. The inbreeding coefficient of the Taiwan WRR population is greater than 0.8, which is similar to the inbred cultivars. The ancestry coefficients map suggests a dispersal pattern of long-distance and seed-mediated contamination across Taiwan, often from warmer, earlier-planted regions to cooler, later-planted regions. Parentage analysis of Taiwan WRR revealed that mostly early indica landraces and indica cultivars were present in the genetic pool; in rare cases temperate japonica was present. Based on the above results, the phylogenetic origin of most Taiwan weedy rice appears to be from hybrid progenies of old cultivated red rice accessions crossed with 'DGWG'. The inbreeding coefficient trend of the six TWR clusters suggests a temporal shift from 'old' indica landraces with red bran (high inbreeding coefficient) to modern indica varieties (low inbreeding coefficient). CONCLUSION: Although there were sustained efforts to remove these old red rice accessions from paddy fields before 1945, some farmers continued to use low purity seed. This practice, along with volunteer cultivation of these old varieties in the second cropping season, apparently has facilitated the long-distance, seed-mediated contamination of rice seed, and the increase in weedy rice seed in paddy soil. © 2019 Society of Chemical Industry.

Validation of Yield Component Traits Identified by Genome-Wide Association Mapping in a tropical japonica × tropical japonica Rice Biparental Mapping Population.

The Rice Diversity Panel 1 (RDP1) was developed for genome-wide association (GWA) studies to explore five rice ( L.) subpopulations (, , , , and ). The RDP1 was evaluated for over 30 traits, including agronomic, panicle architecture, seed, and disease traits and genotyped with 700,000 single nucleotide polymorphisms (SNPs). Most rice grown in the southern United States is and thus the diversity in this subpopulation is interesting to U.S. breeders. Among the RDP1 accessions, 'Estrela' and 'NSFTV199' are both phenotypically and genotypically diverse, thus making them excellent parents for a biparental mapping population. The objectives were to (i) ascertain the GWA QTLs from the RDP1 GWA studies that overlapped with the QTLs uncovered in an Estrela × NSFTV199 recombinant inbred line (RIL) population evaluated for 15 yield traits, and (ii) identify known or novel genes potentially controlling specific yield component traits. The 256 RILs were genotyped with 132 simple sequence repeat markers and 70 QTLs were found. Perl scripts were developed for automatic identification of the underlying candidate genes in the GWA QTL regions. Approximately 100 GWA QTLs overlapped with 41 Estrela × NSFTV199 QTL (RIL QTL) regions and 47 known genes were identified. Two seed trait RIL QTLs with overlapping GWA QTLs were not associated with a known gene. Segregating SNPs in the overlapping GWA QTLs for RIL QTLs with high values will be evaluated as potential DNA markers useful to breeding programs for the associated yield trait.

The rice blast resistance gene Ptr encodes an atypical protein required for broad-spectrum disease resistance.

Plant resistance genes typically encode proteins with nucleotide binding site-leucine rich repeat (NLR) domains. Here we show that Ptr is an atypical resistance gene encoding a protein with four Armadillo repeats. Ptr is required for broad-spectrum blast resistance mediated by the NLR R gene Pi-ta and by the associated R gene Pi-ta2. Ptr is expressed constitutively and encodes two isoforms that are mainly localized in the cytoplasm. A two base pair deletion within the Ptr coding region in the fast neutron-generated mutant line M2354 creates a truncated protein, resulting in susceptibility to M. oryzae. Targeted mutation of Ptr in a resistant cultivar using CRISPR/Cas9 leads to blast susceptibility, further confirming its resistance function. The cloning of Ptr may aid in the development of broad spectrum blast resistant rice.

Functional Interactions Between Major Rice Blast Resistance Genes, Pi-ta and Pi-b, and Minor Blast Resistance Quantitative Trait Loci.

Major blast resistance (R) genes confer resistance in a gene-for-gene manner. However, little information is available on interactions between R genes. In this study, interactions between two rice blast R genes, Pi-ta and Pi-b, and other minor blast resistance quantitative trait loci (QTLs) were investigated in a recombinant inbred line (RIL) population comprising 243 RILs from a Cybonnet (CYBT) × Saber (SB) cross. CYBT has the R gene Pi-ta and SB has Pi-b. Ten differential isolates of four Magnaporthe oryzae races (IB-1, IB-17, IB-49, and IE-1K) were used to evaluate disease reactions of the 243 RILs under greenhouse conditions. Five resistance QTLs were mapped on chromosomes 2, 3, 8, 9, and 12 with a linkage map of 179 single nucleotide polymorphism markers. Among them, qBR12 (Q1), was mapped at the Pi-ta locus and accounted for 45.41% of phenotypic variation while qBR2 (Q2) was located at the Pi-b locus and accounted for 24.81% of disease reactions. The additive-by-additive epistatic interaction between Q1 (Pi-ta) and Q2 (Pi-b) was detected; they can enhance the disease resistance by an additive 0.93 using the 0 to 9 standard phenotyping method. These results suggest that Pi-ta interacts synergistically with Pi-b.

Dynamic Changes in the Rice Blast Population in the United States Over Six Decades.

Rice blast disease caused by Magnaporthe oryzae is one of the most destructive diseases of rice. Field isolates of M. oryzae rapidly adapt to their hosts and climate. Tracking the genetic and pathogenic variability of field isolates is essential to understand how M. oryzae interacts with hosts and environments. In this study, a total of 1,022 United States field isolates collected from 1959 to 2015 were analyzed for pathogenicity toward eight international rice differentials. A subset of 457 isolates was genotyped with 10 polymorphic simple sequence repeat (SSR) markers. The average polymorphism information content value of markers was 0.55, suggesting that the SSR markers were highly informative to capture the population variances. Six genetic clusters were identified by both STRUCTURE and discriminant analysis of principal components methods. Overall, Nei's diversity of M. oryzae in the United States was 0.53, which is higher than previously reported in a world rice blast collection (0.19). The observed subdivision was associated with collection time periods but not with geographic origin of the isolates. Races such as IC-17, IE-1, and IB-49 have been identified across almost all collection periods and all clusters; races such as IA-1, IB-17, and IH-1 have a much higher frequency in certain periods and clusters. Both genomic and pathogenicity changes of United States blast isolates were associated with collection year, suggesting that hosts are a driving force for the genomic variability of rice blast fungus.

Genome-Wide Association of Rice Blast Disease Resistance and Yield-Related Components of Rice.

Robust disease resistance may require an expenditure of energy that may limit crop yield potential. In the present study, a subset of a United States Department of Agriculture rice core collection consisting of 151 accessions was selected using a major blast resistance (R) gene, Pi-ta, marker and was genotyped with 156 simple sequence repeat (SSR) markers. Disease reactions to Magnaporthe oryzae, the causal agent of rice blast disease, were evaluated under greenhouse and field conditions, and heading date, plant height, paddy and brown seed weight in two field environments were analyzed, using an association mapping approach. A total of 21 SSR markers distributed among rice chromosomes 2 to 12 were associated with blast resistance, and 16 SSR markers were associated with seed weight, heading date, and plant height. Most noticeably, shorter plants were significantly correlated with resistance to blast, rice genomes with Pi-ta were associated with lighter seed weights, and the susceptible alleles of RM171 and RM6544 were associated with heavier seed weight. These findings unraveled a complex relationship between disease resistance and yield-related components.

Identification of candidate genes associated with positive and negative heterosis in rice.

To identify the genes responsible for yield related traits, and heterosis, massively parallel signature sequencing (MPSS) libraries were constructed from leaves, roots and meristem tissues from the two parents, 'Nipponbare' and '93-11', and their F1 hybrid. From the MPSS libraries, 1-3 million signatures were obtained. Using cluster analysis, commonly and specifically expressed genes in the parents and their F1 hybrid were identified. To understand heterosis in the F1 hybrid, the differentially expressed genes in the F1 hybrid were mapped to yield related quantitative trait loci (QTL) regions using a linkage map constructed from 131 polymorphic simple sequence repeat markers with 266 recombinant inbred lines derived from a cross between Nipponbare and 93-11. QTLs were identified for yield related traits including days to heading, plant height, plant type, number of tillers, main panicle length, number of primary branches per main panicle, number of kernels per main panicle, total kernel weight per main panicle, 1000 grain weight and total grain yield per plant. Seventy one QTLs for these traits were mapped, of which 3 QTLs were novel. Many highly expressed chromatin-related genes in the F1 hybrid encoding histone demethylases, histone deacetylases, argonaute-like proteins and polycomb proteins were located in these yield QTL regions. A total of 336 highly expressed transcription factor (TF) genes belonging to 50 TF families were identified in the yield QTL intervals. These findings provide the starting genomic materials to elucidate the molecular basis of yield related traits and heterosis in rice.

Similar traits, different genes? Examining convergent evolution in related weedy rice populations.

Convergent phenotypic evolution may or may not be associated with convergent genotypic evolution. Agricultural weeds have repeatedly been selected for weed-adaptive traits such as rapid growth, increased seed dispersal and dormancy, thus providing an ideal system for the study of convergent evolution. Here, we identify QTL underlying weedy traits and compare their genetic architecture to assess the potential for convergent genetic evolution in two distinct populations of weedy rice. F(2) offspring from crosses between an indica cultivar and two individuals from genetically differentiated U.S. weedy rice populations were used to map QTL for four quantitative (heading date, seed shattering, plant height and growth rate) and two qualitative traits. We identified QTL on nine of the twelve rice chromosomes, yet most QTL locations do not overlap between the two populations. Shared QTL among weed groups were only seen for heading date, a trait for which weedy groups have diverged from their cultivated ancestors and from each other. Sharing of some QTL with wild rice also suggests a possible role in weed evolution for genes under selection during domestication. The lack of overlapping QTL for the remaining traits suggests that, despite a close evolutionary relationship, weedy rice groups have adapted to the same agricultural environment through different genetic mechanisms.

Indica and japonica crosses resulting in linkage block and recombination suppression on rice chromosome 12.

Understanding linkage block size and molecular mechanisms of recombination suppression is important for plant breeding. Previously large linkage blocks ranging from 14 megabases to 27 megabases were observed around the rice blast resistance gene Pi-ta in rice cultivars and backcross progeny involving an indica and japonica cross. In the present study, the same linkage block was further examined in 456 random recombinant individuals of rice involving 5 crosses ranging from F(2) to F(10) generation, with and without Pi-ta containing genomic indica regions with both indica and japonica germplasm. Simple sequence repeat markers spanning the entire chromosome 12 were used to detect recombination break points and to delimit physical size of linkage blocks. Large linkage blocks ranging from 4.1 megabases to 10 megabases were predicted from recombinant individuals involving genomic regions of indica and japonica. However, a significantly reduced block from less than 800 kb to 2.1megabases was identified from crosses of indica with indica rice regardless of the existence of Pi-ta. These findings suggest that crosses of indica and japonica rice have significant recombination suppression near the centromere on chromosome 12.

Development of genetic markers linked to straighthead resistance through fine mapping in rice (Oryza sativa L.).

Straighthead, a physiological disorder characterized by sterile florets and distorted spikelets, causes significant yield losses in rice, and occurs in many countries. The current control method of draining paddies early in the season stresses plants, is costly, and wastes water. Development of resistant cultivar is regarded as the most efficient way for its control. We mapped a QTL for straighthead resistance using two recombinant inbred line (RIL) F(9) populations that were phenotyped over two years using monosodium methanearsonate (MSMA) to induce the symptoms. One population of 170 RILs was genotyped with 136 SSRs and the other population of 91 RILs was genotyped with 159 SSRs. A major QTL qSH-8 was identified in an overlapping region in both populations, and explained 46% of total variation in one and 67% in another population for straighthead resistance. qSH-8 was fine mapped from 1.0 Mbp to 340 kb using 7 SSR markers and further mapped to 290 kb in a population between RM22573 and InDel 27 using 4 InDel markers. SSR AP3858-1 and InDel 11 were within the fine mapped region, and co-segregated with straighthead resistance in both RIL populations, as well as in a collection of diverse global accessions. These results demonstrate that AP3858-1 and InDel 11 can be used for marker-assisted selection (MAS) for straighthead resistant cultivars, which is especially important because there is no effective way to directly evaluate straighthead resistance.

Molecular evolution of the rice blast resistance gene Pi-ta in invasive weedy rice in the USA.

The Pi-ta gene in rice has been effectively used to control rice blast disease caused by Magnaporthe oryzae worldwide. Despite a number of studies that reported the Pi-ta gene in domesticated rice and wild species, little is known about how the Pi-ta gene has evolved in US weedy rice, a major weed of rice. To investigate the genome organization of the Pi-ta gene in weedy rice and its relationship to gene flow between cultivated and weedy rice in the US, we analyzed nucleotide sequence variation at the Pi-ta gene and its surrounding 2 Mb region in 156 weedy, domesticated and wild rice relatives. We found that the region at and around the Pi-ta gene shows very low genetic diversity in US weedy rice. The patterns of molecular diversity in weeds are more similar to cultivated rice (indica and aus), which have never been cultivated in the US, rather than the wild rice species, Oryza rufipogon. In addition, the resistant Pi-ta allele (Pi-ta) found in the majority of US weedy rice belongs to the weedy group strawhull awnless (SH), suggesting a single source of origin for Pi-ta. Weeds with Pi-ta were resistant to two M. oryzae races, IC17 and IB49, except for three accessions, suggesting that component(s) required for the Pi-ta mediated resistance may be missing in these accessions. Signatures of flanking sequences of the Pi-ta gene and SSR markers on chromosome 12 suggest that the susceptible pi-ta allele (pi-ta), not Pi-ta, has been introgressed from cultivated to weedy rice by out-crossing.

Genotypic and phenotypic characterization of genetic differentiation and diversity in the USDA rice mini-core collection.

A rice mini-core collection consisting of 217 accessions has been developed to represent the USDA core and whole collections that include 1,794 and 18,709 accessions, respectively. To improve the efficiency of mining valuable genes and broadening the genetic diversity in breeding, genetic structure and diversity were analyzed using both genotypic (128 molecular markers) and phenotypic (14 numerical traits) data. This mini-core had 13.5 alleles per locus, which is the most among the reported germplasm collections of rice. Similarly, polymorphic information content (PIC) value was 0.71 in the mini-core which is the highest with one exception. The high genetic diversity in the mini-core suggests there is a good possibility of mining genes of interest and selecting parents which will improve food production and quality. A model-based clustering analysis resulted in lowland rice including three groups, aus (39 accessions), indica (71) and their admixtures (5), upland rice including temperate japonica (32), tropical japonica (40), aromatic (6) and their admixtures (12) and wild rice (12) including glaberrima and four other species of Oryza. Group differentiation was analyzed using both genotypic distance Fst from 128 molecular markers and phenotypic (Mahalanobis) distance D(2) from 14 traits. Both dendrograms built by Fst and D(2) reached similar-differentiative relationship among these genetic groups, and the correlation coefficient showed high value 0.85 between Fst matrix and D(2) matrix. The information of genetic and phenotypic differentiation could be helpful for the association mapping of genes of interest. Analysis of genotypic and phenotypic diversity based on genetic structure would facilitate parent selection for broadening genetic base of modern rice cultivars via breeding effort.

RL-SAGE and microarray analysis of the rice transcriptome after Rhizoctonia solani infection.

Sheath blight caused by the fungal pathogen Rhizoctonia solani is an emerging problem in rice production worldwide. To elucidate the molecular basis of rice defense to the pathogen, RNA isolated from R. solani-infected leaves of Jasmine 85 was used for both RL-SAGE library construction and microarray hybridization. RL-SAGE sequence analysis identified 20,233 and 24,049 distinct tags from the control and inoculated libraries, respectively. Nearly half of the significant tags (> or =2 copies) from both libraries matched TIGR annotated genes and KOME full-length cDNAs. Among them, 42% represented sense and 7% antisense transcripts, respectively. Interestingly, 60% of the library-specific (> or =10 copies) and differentially expressed (>4.0-fold change) tags were novel transcripts matching genomic sequence but not annotated genes. About 70% of the genes identified in the SAGE libraries showed similar expression patterns (up or down-regulated) in the microarray data obtained from three biological replications. Some candidate RL-SAGE tags and microarray genes were located in known sheath blight QTL regions. The expression of ten differentially expressed RL-SAGE tags was confirmed with RT-PCR. The defense genes associated with resistance to R. solani identified in this study are useful genomic materials for further elucidation of the molecular basis of the defense response to R. solani and fine mapping of target sheath blight QTLs.