Standing genetic variation fuels rapid evolution of herbicide resistance in blackgrass

Repeated herbicide applications exert enormous selection on blackgrass (Alopecurus myosuroides), a major weed in cereal crops of the temperate climate zone including Europe. This inadvertent large-scale experiment gives us the opportunity to look into the underlying genetic mechanisms and evolutionary processes of rapid adaptation, which can occur both through mutations in the direct targets of herbicides and through changes in other, often metabolic, pathways, known as non-target-site resistance. How much either type of adaptation relies on de novo mutations versus pre-existing standing variation is important for developing strategies to manage herbicide resistance. We generated a chromosome-level reference genome for A. myosuroides for population genomic studies of herbicide resistance and genome-wide diversity across Europe in this species. Bulked-segregant analysis evidenced that non-target-site resistance has a complex genetic architecture. Through empirical data and simulations, we showed that, despite its simple genetics, target-site resistance mainly results from standing genetic variation, with only a minor role for de novo mutations.

Identification and Characterization of an Early Leaf Senescence Gene ELS1 in Soybean

Early leaf senescence phenotype in soybean could be helpful to shorten the maturation period and prevent green stem disorder. From a high-density mutation library, we identified two early leaf senescence soybean mutant lines, els1-1 (early leaf senescence 1) and els1-2. The chlorophyll contents of both els1-1 and els1-2 were low in pre-senescent leaves. They degraded rapidly in senescent leaves, revealing that ELS1 is involved in chlorophyll biosynthesis during leaf development and chlorophyll degradation during leaf senescence. The causal mutations in els1 were identified by next-generation sequencing-based bulked segregant analysis. ELS1 encodes the ortholog of the Arabidopsis CaaX-like protease BCM1, which is localized in chloroplasts. Soybean ELS1 was highly expressed in green tissue, especially in mature leaves. The accumulation of photosystem I core proteins and light-harvesting proteins in els1 was low even in pre-senescent leaves, and their degradation was accelerated during leaf senescence. These results suggest that soybean ELS1 is involved in both chlorophyll synthesis and degradation, consistent with the findings in Arabidopsis BCM1. The gene els1, characterized by early leaf senescence and subsequent early maturation, does not affect the flowering time. Hence, the early leaf senescence trait regulated by els1 helps shorten the harvesting period because of early maturation characteristics. The els1-1 allele with weakly impaired function of ELS1 has only a small effect on agricultural traits and could contribute to practical breeding.

Next-generation sequencing-based bulked segregant analysis without sequencing the parental genomes

Genomic regions that control traits of interest can be rapidly identified using BSA-Seq, a technology in which next-generation sequencing (NGS) is applied to bulked segregant analysis (BSA). We recently developed the significant structural variant method for BSA-Seq data analysis that exhibits higher detection power than standard BSA-Seq analysis methods. Our original algorithm was developed to analyze BSA-Seq data in which genome sequences of one parent served as the reference sequences in genotype calling, and thus required the availability of high-quality assembled parental genome sequences. Here we modified the original script to effectively detect the genomic region-trait associations using only bulk genome sequences. We analyzed two public BSA-Seq datasets using our modified method and the standard allele frequency and G-statistic methods with and without the aid of the parental genome sequences. Our results demonstrate that the genomic region(s) associated with the trait of interest could be reliably identified via the significant structural variant method without using the parental genome sequences.

Exome Sequencing from Bulked Segregant Analysis Identifies a Gene for All Stage Resistance to Stripe Rust on Chromosome 1AL in Chinese Wheat Landrace Xiaohemai

Stripe rust caused by Puccinia striiformis f. sp. tritici (Pst) is one of the most destructive diseases of wheat. Identifying novel resistance genes applicable for developing disease resistant cultivars is important for the sustainable control of wheat stripe rust. Chinese wheat landrace Xiaohemai (XHM) is an elite germplasm line with all-stage resistance (ASR) effective against predominant Chinese Pst races. In this study, we performed a bulked segregant analysis coupled with exome capture sequencing (BSE-seq) to identify a candidate genomic region strongly associated with stripe rust resistance on chromosome 1AL in 173 F2:3 lines derived from cross XHM × Avocet S. The gene, designated as YrXH-1AL, was validated by a conventional quantitative trait locus analysis using newly developed Kompetitive allele-specific PCR (KASP) markers, explaining up to 48.50% of the phenotypic variance. By testing a secondary mapping population comprising 144 lines from the same cross at the seedling stage with prevalent Pst race CYR34, YrXH-1AL was identified as a single Mendelian factor in a 1.5 cM interval flanked by KASP markers KP1A_484.33 and KP1A_490.09. This region corresponded to a 5.76 Mb genomic interval on Chinese Spring chromosome 1AL. Furthermore, two co-segregating KASP markers showed high polymorphisms among 130 Chinese wheat cultivars and could be used for marker-assisted selection. Because no other Yr genes for ASR that originated from common wheat have been detected on chromosome 1AL, YrXH-1AL is likely a novel gene that can be incorporated into modern breeding materials to develop wheat cultivars with enhanced stripe rust resistance.

Identification, Analysis, and Confirmation of Seed Storability-Related Loci in Dongxiang Wild Rice (Oryza rufipogon Griff.)

Dongxiang wild rice (Oryza rufipogon Griff.) (DXWR) has strong seed storability and identifying its elite gene resources may facilitate genetic improvements in rice seed storability. In this study, we developed two backcross inbred lines (BILs) populations, with DXWR as a common donor parent and two rice varieties (F6 and R974) as recipient parents. Bulked segregant analysis via whole genome sequencing (BSA-seq) was used to identify seed storability-related loci in the DXWR and F6 population. Two main genomic regions containing 18,550,000–20,870,000 bp on chromosome 4 and 7,860,000–9,780,000 bp on chromosome 9 were identified as candidate loci of DXWR seed storability; these overlapped partially with seed storability-related quantitative trait loci (QTLs) discovered in previous studies, suggesting that these loci may provide important regions for isolating the responsible genes. In total, 448 annotated genes were predicted within the identified regions, of which 274 and 82 had nonsynonymous and frameshift mutations, respectively. We detected extensive metabolic activities and cellular processes during seed storability and confirmed the effects of the seed storability-related candidate loci using four BILs from DXWR and R974. These results may facilitate the cloning of DXWR seed storability-related genes, thereby elucidating rice seed storability and its improvement potential.

Identification, analysis, and confirmation of seed storability-related loci in Dongxiang wild rice (Oryza rufipogon Griff.)

Dongxiang wild rice (Oryza rufipogon Griff.) (DXWR) has strong seed storability and identifying its elite gene resources may facilitate genetic improvements in rice seed storability. In this study, we developed two backcross inbred lines (BILs) populations, with DXWR as a common donor parent and two rice varieties (F6 and R974) as recipient parents. Bulked segregant analysis via whole genome sequencing (BSAseq) was used to identify seed storability related loci in the DXWR and F6 population. Two main genomic regions containing 18,550,000 to 20,870,000 bp on chromosome 4 and 7,860,000 to 9,780,000 bp on chromosome 9 were identified as candidate loci of DXWR seed storability; these overlapped partially with seed storability related quantitative trait loci (QTLs) discovered in previous studies, suggesting that these loci may provide important regions for isolating the responsible genes. In total, 448 annotated genes were predicted within the identified regions, of which 274 and 82 had nonsynonymous and frameshift mutations, respectively. We detected extensive metabolic activities and cellular processes during seed storability and confirmed the effects of the seed storability related candidate loci using four BILs from DXWR and R974. These results may facilitate the cloning of DXWR seed storability-related genes, thereby elucidating rice seed storability and its improvement potential.

Identification of a Novel Salt Tolerance-Related Locus in Wild Soybean (Glycine soja Sieb. & Zucc.)

Salinity is an important abiotic stress factor that affects growth and yield of soybean. NY36-87 is a wild soybean germplasm with high salt tolerance. In this study, two F2:3 mapping populations derived from NY36-87 and two salt-sensitive soybean cultivars, Zhonghuang39 and Peking, were used to map salt tolerance-related genes. The two populations segregated as 1 (tolerant):2 (heterozygous):1 (sensitive), indicating a Mendelian segregation model. Using simple sequence repeat (SSR) markers together with the bulked segregant analysis (BSA) mapping strategy, we mapped a salt tolerance locus on chromosome 03 in F2:3 population Zhonghuang39×NY36-87 to a 98-kb interval, in which the known gene GmSALT3 co-segregated with the salt tolerance locus. In the F2:3 population of Peking×NY36-87, the dominant salt tolerance-associated gene was detected and mapped on chromosome 18. We named this gene GmSALT18 and fine mapped it to a 241-kb region. Time course analysis and a grafting experiment confirmed that Peking accumulated more Na+ in the shoot via a root-based mechanism. These findings reveal that the tolerant wild soybean line NY36-87 contains salt tolerance-related genes GmSALT3 and GmSALT18, providing genetic material and a novel locus for breeding salt-tolerant soybean.

Molecular Analysis Uncovers the Mechanism of Fertility Restoration in Temperature-Sensitive Polima Cytoplasmic Male-Sterile Brassica napus

Temperature-sensitive male sterility is a heritable agronomic trait affected by genotype-environment interactions. In rapeseed (Brassica napus), Polima (pol) temperature-sensitive cytoplasmic male sterility (TCMS) is commonly used for two-line breeding, as the fertility of pol TCMS lines can be partially restored at certain temperatures. However, little is known about the underlying molecular mechanism that controls fertility restoration. Therefore, we aimed to investigate the fertility conversion mechanism of the pol TCMS line at two different ambient temperatures (16 °C and 25 °C). Our results showed that the anthers developed and produced vigorous pollen at 16 °C but not at 25 °C. In addition, we identified a novel co-transcript of orf224-atp6 in the mitochondria that might lead to fertility conversion of the pol TCMS line. RNA-seq analysis showed that 1637 genes were significantly differentially expressed in the fertile flowers of 596-L when compared to the sterile flower of 1318 and 596-H. Detailed analysis revealed that differentially expressed genes were involved in temperature response, ROS accumulation, anther development, and mitochondrial function. Single-molecule long-read isoform sequencing combined with RNA sequencing revealed numerous genes produce alternative splicing transcripts at high temperatures. Here, we also found that alternative oxidase, type II NAD(P)H dehydrogenases, and transcription factor Hsfs might play a crucial role in male fertility under the low-temperature condition. RNA sequencing and bulked segregant analysis coupled with whole-genome sequencing identified the candidate genes involved in the post-transcriptional modification of orf224. Overall, our study described a putative mechanism of fertility restoration in a pol TCMS line controlled by ambient temperature that might help utilise TCMS in the two-line breeding of Brassica crops.

Genetic Mapping and Identification of the Candidate Genes for Mottled Rind in Cucumis melo L.

The rind appearance of melon is one of the most vital commercial quality traits which determines the preferences and behavior of consumers toward the consumption of melon. In this study, we constructed an F2 population derived from SC (mottled rind) and MG (non-mottled rind) lines for mapping the mottled rind gene(s) in melon. Genetic analysis showed that there were two dominant genes (CmMt1 and CmMt2) with evidence of epistasis controlling the mottled rind. Meanwhile, the phenotypic segregation ratio implied that the immature rind color had an epistatic effect on the mottled rind, which was regulated by CmAPRR2. A Kompetitive Allele-Specific PCR (KASP) DNA marker (CmAPRR2SNP(G/T)) was developed and shown to co-segregate with rind color, confirming that CmAPRR2 was CmMt1. Using bulked segregant analysis sequencing and KASP assays, CmMt2 was fine-mapped to an interval of 40.6 kb with six predicted genes. Functional annotation, expression analysis, and sequence variation analyses confirmed that AtCPSFL1 homolog, MELO3C026282, was the most likely candidate gene for CmMt2. Moreover, pigment content measurement and transmission electron microscopy analysis demonstrated that CmMt2 might participate in the development of chloroplast, which, in turn, decreases the accumulation of chlorophyll. These results provide insight into the molecular mechanism underlying rind appearance and reveal valuable information for marker-assisted selection breeding in melon.

Optimization of BSA-seq experiment for QTL mapping

Deep sequencing-based bulked segregant analysis (BSA-seq) has become a popular approach for quantitative trait loci (QTL) mapping in recent years. Effective statistical methods for BSA-seq have been developed, but how to design a suitable experiment for BSA-seq remains unclear. In this paper, we show in theory how the major experimental factors (including population size, pool proportion, pool balance, and generation) and the intrinsic factors of a QTL (including heritability and degree of dominance) affect the power of QTL detection and the precision of QTL mapping in BSA-seq. Increasing population size can improve the power and precision, depending on the QTL heritability. The best proportion of each pool in the population is around 0.25. So, 0.25 is generally applicable in BSA-seq. Small pool proportion can greatly reduce the power and precision. Imbalance of pool pair in size also causes decrease of the power and precision. Additive effect is more important than dominance effect for QTL mapping. Increasing the generation of filial population produced by selfing can significantly increase the power and precision, especially from F2 to F3. These findings enable researchers to optimize the experimental design for BSA-seq. A web-based program named BSA-seq Design Tool is available at http://124.71.74.135/BSA-seqDesignTool/ and https://github.com/huanglikun/BSA-seqDesignTool.