Rapid Identification by Resequencing-Based QTL Mapping of a Novel Allele RGA1-FH Decreasing Grain Length in a Rice Restorer Line 'Fuhui212'

Grain size is one of the most frequently selected traits during domestication and modern breeding. The continued discovery and characterization of new genes and alleles in controlling grain size are important in safeguarding the food supply for the world's growing population. Previously, a small grain size was observed in a rice restorer line 'Fuhui212', while the underlying genetic factors controlling this trait were unknown. In this study, by combining QTL mapping, variant effect prediction, and complementation experiments, we recovered a novel allele RGA1-FH that explains most of the phenotypic changes. The RGA1-FH allele contains an A-to-T splicing site variant that disrupts the normal function of RGA1. While population analysis suggests extremely strong artificial selection in maintaining a functional allele of RGA1, our study is the first, to the best of our knowledge, to prove that a dysfunctional RGA1 allele can also be beneficial in real agricultural production. Future breeding programs would benefit from paying more attention to the rational utilization of those overlooked 'unfavored' alleles.

High level of somatic mutations detected in a diploid banana wild relative Musa basjoo

Plants are thought to lack an early segregating germline and often retain both asexual and sexual reproduction, both of which may allow somatic mutations to enter the gametes or clonal progeny, and thereby impact plant evolution. It is yet unclear how often these somatic mutations occur during plant development and what proportion is transmitted to their sexual or cloned offspring. Asexual "seedless" propagation has contributed greatly to the breeding in many fruit crops, such as citrus, grapes and bananas. Whether plants in these lineages experience substantial somatic mutation accumulation is unknown. To estimate the somatic mutation accumulation and inheritance among a clonal population of plant, here we assess somatic mutation accumulation in Musa basjoo, a diploid banana wild relative, using 30 whole-genome resequenced samples collected from five structures, including leaves, sheaths, panicle, roots and underground rhizome connecting three clonal individuals. We observed 18.5 high proportion de novo somatic mutations on average between each two adjacent clonal suckers, equivalent to ~ 2.48 × 10^-8 per site per asexual generation, higher than the per site per sexual generation rates (< 1 × 10^-8) reported in Arabidopsis and peach. Interestingly, most of these inter-ramet somatic mutations were shared simultaneously in different tissues of the same individual with a high level of variant allele fractions, suggesting that these somatic mutations arise early in ramet development and that each individual may develop only from a few apical stem cells. These results thus suggest substantial mutation accumulation in a wild relative of banana. Our work reveals the significance of somatic mutation in Musa basjoo genetics variations and contribute to the trait improvement breeding of bananas and other asexual clonal crops.

In rice splice variants that restore the reading frame after frameshifting indel introduction are common, often induced by the indels and sometimes lead to organism-level rescue

The introduction of frameshifting non-3n indels enables the identification of gene-trait associations. However, it has been hypothesised that recovery of the original reading frame owing to usage of non-canonical splice forms could cause rescue. To date there is very little evidence for organism-level rescue by such a mechanism and it is unknown how commonly indels induce, or are otherwise associated with, frame-restoring splice forms. We perform CRISPR/Cas9 editing of randomly selected loci in rice to investigate these issues. We find that the majority of loci have a frame-restoring isoform. Importantly, three quarters of these isoforms are not seen in the absence of the indels, consistent with indels commonly inducing novel isoforms. This is supported by analysis in the context of NMD knockdowns. We consider in detail the two top rescue candidates, in wax deficient anther 1 (wda1) and brittle culm (bc10), finding that organismal-level rescue in both cases is strong but owing to different splice modification routes. More generally, however, as frame-restoring isoforms are low abundance and possibly too disruptive, such rescue we suggest to be the rare exception, not the rule. Nonetheless, assuming that indels commonly induce frame-restoring isoforms, these results emphasize the need to examine RNA level effects of non-3n indels and suggest that multiple non-3n indels in any given gene are advisable to probe a gene’s trait associations.

As protein coding genes are read in units of three (codons), insertions or deletions (indels) that are not a multiple of three long (non 3n) are expected to be especially harmful. Whether they are is important both for interpreting the results of non-3n indel experiments to probe a gene’s functional importance and for diagnostics. Particularly enigmatic are incidences where some non-3n changes in a gene compromise phenotypes while other seemingly comparable ones do not. One explanation for the latter is that a non-3n indel might be rescued via a frame-restoring splice form. Here we examine this hypothesis by inducing non-3n indels in many genes in rice and find that many non-3n indels are associated with a splice form that restores the reading frame. In the majority of these cases the indel appears to induce the potential rescuing splice form. We examine two top hit cases in detail and show functional rescue by splice modification. More generally, the frame-restoring forms are, however, low abundance and probably result in compromised proteins. We conclude then that splice mediated rescue is possible, but probably uncommon. Nonetheless it should not be overlooked in experimental design and interpretation.

Ancient rapid functional differentiation and fixation of the duplicated members in rice Dof genes after whole genome duplication

Whole genome duplication (WGD) in plants is typically followed by genomic downsizing, where large portions of the new genome are lost. Whether this downsizing is accompanied by increased or decreased evolutionary rates of the remaining genes is poorly known, not least because homeolog pairings are often obscured by chromosomal rearrangement. Here, we use the newly published genome from a sedge, namely Kobresia littledalei, and CRISPR/Cas-9 editing to investigate how the Rho WGD event 70 million years ago (MYA) affected transcription factor evolutionary rates, fates, and function in rice (Oryza sativa) and sorghum (Sorghum bicolor). We focus on the 30-member DNA-binding with one zinc finger (Dof) transcription factor family in both crops due to their agronomic importance. Using the known speciation dates of rice from Kobresia (97 MYA) and sorghum (50 MYA), we find that rates of amino acid substitution in the critical Dof domain region were over twofold higher during the 20-million-year period following the WGD than before or afterward. Through comparison of synteny blocks, we report that at least 11% of Dof genes were purged from 70 to 50 MYA, while only 6% have been lost in the most recent 50-million-year interval. CRISPR/Cas9 editing revealed widespread fitness-related defects in flowering and lack of redundancy of paired members, as well as significant differences in expression between gene pairs. Together these findings demonstrate the strength of Dof genes as a model for deep evolutionary study and offer one of the most detailed portraits yet of the Rho WGD impact on a gene lineage.

Mitotic gene conversion can be as important as meiotic conversion in driving genetic variability in plants and other species without early germline segregation

In contrast to common meiotic gene conversion, mitotic gene conversion, because it is so rare, is often ignored as a process influencing allelic diversity. We show that if there is a large enough number of premeiotic cell divisions, as seen in many organisms without early germline sequestration, such as plants, this is an unsafe position. From examination of 1.1 million rice plants, we determined that the rate of mitotic gene conversion events, per mitosis, is 2 orders of magnitude lower than the meiotic rate. However, owing to the large number of mitoses between zygote and gamete and because of long mitotic tract lengths, meiotic and mitotic gene conversion can be of approximately equivalent importance in terms of numbers of markers converted from zygote to gamete. This holds even if we assume a low number of premeiotic cell divisions (approximately 40) as witnessed in Arabidopsis. A low mitotic rate associated with long tracts is also seen in yeast, suggesting generality of results. For species with many mitoses between each meiotic event, mitotic gene conversion should not be overlooked.

Gene conversion associated with meiosis has long been a focus of attention in population genomics, but mitotic conversion has been relatively overlooked as it was thought to be rare. Analysis in plants suggests that this could be a mistake; long tract lengths and multiple mitoses in species lacking germline sequestration suggest that mitotic conversion, although rare per mitosis, should not be ignored.

Comparative population genomic analysis provides insights into breeding of modern indica rice in China

Comparative genomic analysis within Asian cultivated rice (Oryza sativa L.) populations has greatly enriched our knowledge regarding rice domestication and the divergence of the indica and japonica subspecies, while study on genomic regions associated with improvement within the indica subspecies is still limited. Here, through combined investigation of 2,429 indica cultivar genomes from public sequencing projects, we depict the improvement of modern indica rice in China. We identify three subgroups within indica populations: two geographically distinct, historical subgroups indica I (Ind_I) and indica III (Ind_III) and a modern subgroup indica II (Ind_II). The modern indica subgroup Ind_II shows admixture of the other two subgroups and enrichment of alleles that had been low-frequency in the other two subgroups. The Chinese indica cultivars exhibit a strong subgroup component change from Ind_I to Ind_II in the 1980s. Through haplotype-based comparative analysis, we detect 187 regions associated with separation of Ind_II compared to Ind_I or Ind_III. Within those regions we find strong representation of beneficial agricultural production-related alleles in Ind_II and a positive correlation between grain yield and number of differentiated haplotypes. Phenotypic features of long and slender grain, small tiller angle and decreased flowering time were detected for Ind_II. Through haplotype-based comparative analysis between rice subpopulations and subspecies, we find differentiated haplotypes not only from indica itself but also from japonica and aus, suggesting that introgression from other rice sub-populations has substantially contributed to modern indica rice breeds. These results help clarify the evolutionary landscape of modern indica rice in China and provide useful targets for future improvement.

Different knockout genotypes of OsIAA23 in rice using CRISPR/Cas9 generating different phenotypes

We have isolated several Osiaa23 rice mutants with different knockout genotypes, resulting in different phenotypes, which suggested that different genetic backgrounds or mutation types influence gene function. The Auxin/Indole-3-Acetic Acid (Aux/IAA) gene family performs critical roles in auxin signal transduction in plants. In rice, the gene OsIAA23 (Os06t0597000) is known to affect development of roots and shoots, but previous knockouts in OsIAA23 have been sterile and difficult for research continuously. Here, we isolate new Osiaa23 mutants using the CRISPR/Cas9 system in japonica (Wuyunjing24) and indica (Kasalath) rice, with extensive genome re-sequencing to confirm the absence of off-target effects. In Kasalath, mutants with a 13-amino acid deletion showed profoundly greater dwarfing, lateral root developmental disorder, and fertility deficiency, relative to mutants with a single amino acid deletion, demonstrating that those 13 amino acids in Kasalath are essential to gene function. In Wuyunjing24, we predicted that mutants with a single base-pair frameshift insertion would experience premature termination and strong phenotypic defects, but instead these lines exhibited negligible phenotypic difference and normal fertility. Through RNA-seq, we show here that new mosaic transcripts of OsIAA23 were produced de novo, which circumvented the premature termination and thereby preserved the wild-type phenotype. This finding is a notable demonstration in plants that mutants can mask loss of function CRISPR/Cas9 editing of the target gene through de novo changes in alternative splicing.

High-resolution insight into recombination events at the SD1 locus in rice

Recombination during meiosis plays an important role in genome evolution by reshuffling existing genetic variations into fresh combinations with the possibility of recovery of lost ancestral genotypes. While crossover (CO) events have been well studied, gene conversion events (GCs), which represent non-reciprocal information transfer between chromosomes, are poorly documented and difficult to detect due to their relatively small converted tract size. Here, we document these GC events and their phenotypic effects at an important locus in rice containing the SD1 gene, where multiple defective alleles contributed to the semi-dwarf phenotype of rice in the 'Green Revolution' of the 1960s. Here, physical separation of two defects allows recombination to generate the wild-type SD1 gene, for which plant height can then be used as a reporter. By screening 18 000 F₂ progeny from a cross between two semi-dwarf cultivars that carry these different defective alleles, we detected 24 GC events, indicating a conversion rate of ~3.3 × 10^-4 per marker per generation in a single meiotic cycle in rice. Furthermore, our data show that indels and single-nucleotide polymorphisms (SNPs) do not differ significantly in GC rates, at least at the SD1 locus. Our results provide strong evidence that GC by itself can regain an ancestral phenotype that was lost through mutation. This GC detection approach is likely to be broadly applicable to natural or artificial alleles of other phenotype-related functional genes, which are abundant in other plant genomes.

Salicylic acid inhibits jasmonic acid-induced resistance of Arabidopsis thaliana to Spodoptera exigua

The role of salicylic acid (SA) in plant responses to pathogens has been well documented, but its direct and indirect effects on plant responses to insects are not so well understood. We examined the effects of SA, alone and in combination with jasmonic acid (JA), on the performance of the generalist herbivore, Spodoptera exigua, in wild-type and mutant Arabidopsis thaliana genotypes that varied genetically in their ability to mount SA- and JA-mediated defence responses. In one experiment, growth of S. exigua larvae was highest on the Wassilewskija wild-type, intermediate on the Columbia wild-type and the JA-deficient fad mutant, and lowest on the nim1-1 and jar1-mutants, which are defective in the SA and JA pathways, respectively. Activity of guaiacol peroxidase, polyphenoloxidase, n-acetylglucosaminidase, and trypsin inhibitor varied by genotype but did not correlate with insect performance. SA treatment increased growth of S. exigua larvae by approximately 35% over all genotypes, but had no discernable effect on activities of the four defence proteins. In a second experiment, growth of S. exigua was highest across treatments on the cep1 mutant, a constitutive expressor of high SA levels and systemic acquired resistance, and lowest on the fad mutant, which is JA-deficient. JA treatment generally increased activity of all four defence proteins, increased total glucosinolate levels and reduced insect growth by approximately 25% over all genotypes. SA generally inhibited expression of JA-induced resistance to S. exigua when both hormones were applied simultaneously. Across genotypes and treatments, larval mass was negatively correlated with the activity of trypsin inhibitor and polyphenoloxidase and with total glucosinolate levels, and insect damage was negatively correlated with the activity of polyphenoloxidase. SA had little effect on the induction of defence protein activity by JA. However, SA attenuated the induction of glucosinolates by JA and therefore may explain better the interactive effects of SA and JA on insect performance. This study illustrates that direct and indirect cross-effects of SA on resistance to S. exigua can occur in A. thaliana. Effects of SA may be mediated through effects on plant defence chemistry or other aspects of the suitability of foliage for insect feeding and growth.

Negative cross-talk between salicylate- and jasmonate-mediated pathways in the Wassilewskija ecotype of Arabidopsis thaliana

Plants often respond to attack by insect herbivores and necrotrophic pathogens with induction of jasmonate-dependent resistance traits, but respond to attack by biotrophic pathogens with induction of salicylate-dependent resistance traits. To assess the degree to which the jasmonate- and salicylate-dependent pathways interact, we compared pathogenesis-related protein activity and bacterial performance in four mutant Arabidopsis thaliana lines relative to their wild-type backgrounds. We found that two salicylate-dependent pathway mutants (cep1, nim1-1) exhibited strong effects on the growth of the generalist biotrophic pathogen, Pseudomonas syringae pv. tomato, whereas two jasmonate-dependent pathway mutants (fad3-2fad7-2fad8, jar1-1) did not. Leaf peroxidase and exochitinase activity were negatively correlated with bacterial growth, whereas leaf polyphenol oxidase activity and trypsin inhibitor concentration were not. Interestingly, leaf total glucosinolate concentration was positively correlated with bacterial growth. In the same experiment, we also found that application of jasmonic acid generally increased leaf peroxidase activity and trypsin inhibitor concentration in the mutant lines. However, the cep1 mutant, shown previously to overexpress salicylic acid, exhibited no detectable biological or chemical responses to jasmonic acid, suggesting that high levels of salicylic acid may have inhibited a plant response. In a second experiment, we compared the effect of jasmonic acid and/or salicylic acid on two ecotypes of A. thaliana. Application of salicylic acid to the Wassilewskija ecotype decreased bacterial growth. However, this effect was not observed when both salicylic acid and jasmonic acid were applied, suggesting that jasmonic acid negated the beneficial effect of salicylic acid. Collectively, our results confirm that the salicylate-dependent pathway is more important than the jasmonate-dependent pathway in determining growth of P. syringae pv. tomato in A. thaliana, and suggest important negative interactions between these two major defensive pathways in the Wassilewskija ecotype. In contrast, the Columbia ecotype exhibited little evidence of negative interactions between the two pathways, suggesting intraspecific variability in how these pathways interact in A. thaliana.

Fitness costs of R-gene-mediated resistance in Arabidopsis thaliana

Resistance genes (R-genes) act as an immune system in plants by recognizing pathogens and inducing defensive pathways. Many R-gene loci are present in plant genomes, presumably reflecting the need to maintain a large repertoire of resistance alleles. These loci also often segregate for resistance and susceptibility alleles that natural selection has maintained as polymorphisms within a species for millions of years. Given the obvious advantage to an individual of being disease resistant, what prevents these resistance alleles from being driven to fixation by natural selection? A cost of resistance is one potential explanation; most models require a lower fitness of resistant individuals in the absence of pathogens for long-term persistence of susceptibility alleles. Here we test for the presence of a cost of resistance at the RPM1 locus of Arabidopsis thaliana. Results of a field experiment comparing the fitness of isogenic strains that differ in the presence or absence of RPM1 and its natural promoter reveal a large cost of RPM1, providing the first evidence that costs contribute to the maintenance of an ancient R-gene polymorphism.