Agile Genetics: Single gene resolution without the fuss.

Gene discovery reveals new biology, expands the utility of marker-assisted selection, and enables targeted mutagenesis. Still, such discoveries can take over a decade. We present a general strategy, "Agile Genetics," that uses nested, structured populations to overcome common limits on gene resolution. Extensive simulation work on realistic genetic architectures shows that, at population sizes of >5000 samples, single gene-resolution can be achieved using bulk segregant pools. At this scale, read depth and technical replication become major drivers of resolution. Emerging enrichment methods to address coverage are on the horizon; we describe one possibility - iterative depth sequencing (ID-seq). In addition, graph-based pangenomics in experimental populations will continue to maximize accuracy and improve interpretation. Based on this merger of agronomic scale with molecular and bioinformatic innovation, we predict a new age of rapid gene discovery.

A FLOWERING LOCUS T ortholog is associated with photoperiod-insensitive flowering in hemp (Cannabis sativa L.).

Hemp (Cannabis sativa L.) is an extraordinarily versatile crop, with applications ranging from medicinal compounds to seed oil and fibre products. Cannabis sativa is a short-day plant, and its flowering is highly controlled by photoperiod. However, substantial genetic variation exists for photoperiod sensitivity in C. sativa, and photoperiod-insensitive ("autoflower") cultivars are available. Using a bi-parental mapping population and bulked segregant analysis, we identified Autoflower2, a 0.5 Mbp locus significantly associated with photoperiod-insensitive flowering in hemp. Autoflower2 contains an ortholog of the central flowering time regulator FLOWERING LOCUS T (FT) from Arabidopsis thaliana which we termed CsFT1. We identified extensive sequence divergence between alleles of CsFT1 from photoperiod-sensitive and insensitive cultivars of C. sativa, including a duplication of CsFT1 and sequence differences, especially in introns. Furthermore, we observed higher expression of one of the CsFT1 copies found in the photoperiod-insensitive cultivar. Genotyping of several mapping populations and a diversity panel confirmed a correlation between CsFT1 alleles and photoperiod response, affirming that at least two independent loci involved in the photoperiodic control of flowering, Autoflower1 and Autoflower2, exist in the C. sativa gene pool. This study reveals the multiple independent origins of photoperiod insensitivity in C. sativa, supporting the likelihood of a complex domestication history in this species. By integrating the genetic relaxation of photoperiod sensitivity into novel C. sativa cultivars, expansion to higher latitudes will be permitted, thus allowing the full potential of this versatile crop to be reached.

Identification of candidate genes associated with peanut pod length by combined analysis of QTL-seq and RNA-seq.

Pod length (PL) is one of the major traits determining pod size and yield of peanut. Discovering the quantitative trait loci (QTL) and identifying candidate genes associated with PL are essential for breeding high-yield peanut. In this study, quantitative trait loci sequencing (QTL-seq) was performed using the F2 population constructed by a short-pod variety Tifrunner (Tif) and a long-pod line Lps, and a 0.77 Mb genomic region on chromosome 07 was identified as the candidate region for PL. Then, the candidate region was narrowed to a 265.93 kb region by traditional QTL approach. RNA-seq analysis showed that there were four differentially expressed genes (DEGs) in the candidate region, among which Arahy.PF2L6F (AhCDC48) and Arahy.P4LK2T (AhTAA1) were speculated to be PL-related candidate genes. These results were informative for the elucidation of the underlying regulatory mechanism in peanut pod length and would facilitate further identification of valuable target genes.

Mapping and Candidate Gene Analysis of the Low-Temperature-Sensitive Albino Gene OsLTSA8 in Rice Seedlings.

Chloroplasts are organelles responsible for photosynthesis in plants, providing energy for growth and development. However, the genetic regulatory mechanisms underlying early chloroplast development in rice remain incompletely understood. In this study, we identified a rice seedling thermosensitive chlorophyll-deficient mutant, osltsa8, and the genetic analysis of two F2 populations suggested that this trait may be controlled by more than one pair of alleles. Through reciprocal F2 populations and QTL-seq technology, OsLTSA8 was mapped to the interval of 24,280,402-25,920,942 bp on rice chromosome 8, representing a novel albino gene in rice. Within the candidate gene region of OsLTSA8, there were 258 predicted genes, among which LOC_Os08g39050, LOC_Os08g39130, and LOC_Os08g40870 encode pentatricopeptide repeat (PPR) proteins. RNA-seq identified 18 DEGs (differentially expressed genes) within the candidate interval, with LOC_Os08g39420 showing homology to the pigment biosynthesis-related genes Zm00001d017656 and Sb01g000470; LOC_Os08g39430 and LOC_Os08g39850 were implicated in chlorophyll precursor synthesis. RT-qPCR was employed to assess the expression levels of LOC_Os08g39050, LOC_Os08g39130, LOC_Os08g40870, LOC_Os08g39420, LOC_Os08g39430, and LOC_Os08g39850 in the wild-type and mutant plants. Among them, the differences in the expression levels of LOC_Os08g39050 and LOC_Os08g39430 were the most significant. This study will contribute to further elucidating the molecular mechanisms of rice chloroplast development.

QTL-seq identifies genomic region associated with the crown root development under Jasmonic acid response.

Rice root system plays a crucial role in plant adaptation under adverse conditions, particularly drought stress. However, the regulatory gene networks that govern rice root development during stress exposure remain largely unexplored. In this study, we applied a QTL sequencing method to identify QTL/gene controlling the crown root development under Jasmonic acid simulation using the Bulk-segregant analysis. Two rice cultivars with contrasting phenotypes from the Vietnamese traditional rice collection were used as parent pairs for crossing. The single-seed descent method was employed to generate an F2 population of progenies. This F2/3 population was further segregated based on root count under JA stress. Pooled DNA from the two extreme groups in this population was sequenced, and SNP indexes across all loci in these pools were calculated. We detected a significant genomic region on chromosome 10, spanned from 20.39-20.50 Mb, where two rice RLKs were located, OsPUB54 and OsPUB58. Receptor-like kinases (RLKs) are pivotal in regulating various aspects of root development in plants, and the U-box E3 ubiquitination ligase class was generally known for its degradation of some protein complexes. Notably, OsPUB54 was strongly induced by JA treatment, suggesting its involvement in the degradation of the Aux/IAA protein complex, thereby influencing crown root initiation. Besides, the Eukaryotic translation initiation of factor 3 subunit L (eIF3l) and the Mitogen-activated protein kinase kinase kinase 37 (MAPKKK 37) proteins identified from SNPs with high score index which suggests their significant roles in the translation initiation process and cellular signaling pathways, respectively. This information suggests several clues of how these candidates are involved in modifying the rice root system under stress conditions.

Homoeologous exchanges contribute to branch angle variations in rapeseed: Insights from transcriptome, QTL-seq and gene functional analysis.

Branch angle (BA) is a critical morphological trait that significantly influences planting density, light interception and ultimately yield in plants. Despite its importance, the regulatory mechanism governing BA in rapeseed remains poorly understood. In this study, we generated 109 transcriptome data sets for 37 rapeseed accessions with divergent BA phenotypes. Relative to adaxial branch segments, abaxial segments accumulated higher levels of auxin and exhibited lower expression of six TCP1 homologues and one GA20ox3. A co-expression network analysis identified two modules highly correlated with BA. The modules contained homologues to known BA control genes, such as FUL, YUCCA6, TCP1 and SGR3. Notably, a homoeologous exchange (HE), occurring at the telomeres of A09, was prevalent in large BA accessions, while an A02-C02 HE was common in small BA accessions. In their corresponding regions, these HEs explained the formation of hub gene hotspots in the two modules. QTL-seq analysis confirmed that the presence of a large A07-C06 HE (~8.1 Mb) was also associated with a small BA phenotype, and BnaA07.WRKY40.b within it was predicted as candidate gene. Overexpressing BnaA07.WRKY40.b in rapeseed increased BA by up to 20°, while RNAi- and CRISPR-mediated mutants (BnaA07.WRKY40.b and BnaC06.WRKY40.b) exhibited decreased BA by up to 11.4°. BnaA07.WRKY40.b was exclusively localized to the nucleus and exhibited strong expression correlations with many genes related to gravitropism and plant architecture. Taken together, our study highlights the influence of HEs on rapeseed plant architecture and confirms the role of WRKY40 homologues as novel regulators of BA.

Impairment in a member of AP2/ERF and F-box family protein enhances complete panicle exsertion, an economically important trait in rice.

Complete panicle exsertion (CPE) is an economically important quantitative trait that contributes to grain yield in rice. We deployed an integrated approach for understanding the molecular mechanism of CPE using a stable EMS mutant line, CPE-109 of Samba Mahsuri (SM) exhibiting CPE. Two consistent genomic regions have been identified for CPE through QTL mapping [qCPE-4 (28.24-31.22 Mb) and qCPE-12 (2.30-3.18 Mb)] and QTL-sequencing [Chr-4 (31.21-33.69 Mb) and Chr-12 (0.12-3.15 Mb)]. Two non-synonymous SNPs, viz; KASP 12-12 (T→C; Chr12:1269983) in Os12g0126300; AP2/ERF transcription factor and KASP 12-16 (G→A; Chr12:1515198) in Os12g0131400; F-box domain-containing protein explained 81.05 and 59.61% phenotypic variance respectively and exhibited strong co-segregation with CPE in F2 mapping populations, advanced generation lines and CPE exhibiting SM mutants through KASP assays. The downregulation of these genes in CPE-109 compared to SM was observed in transcriptome sequencing of flag leaves which was validated through qRT-PCR. We propose that the abrogation of Os12g0126300 and Os12g0131400 in CPE-109 combinatorially influences the downregulation of ethylene biosynthetic genes viz. ACC synthase, ethylene-responsive factor-2, and up-regulation of gibberellic acid synthetic genes viz. ent-kaurene synthase and two cytokinin biosynthesis genes viz. cytokinin-O-glucosyltransferase 2, carboxy-lyase which result in complete panicle exsertion.

Combined QTL-sequencing, linkage mapping, and RNA-sequencing identify candidate genes and KASP markers for low-temperature germination in Oryza sativa L. ssp. Japonica.

MAIN CONCLUSION: Through QTL-seq, QTL mapping and RNA-seq, six candidate genes of qLTG9 can be used as targets for cold tolerance functional characterization, and six KASP markers can be used for marker-assisted breeding to improve the germination ability of japonica rice at low temperature. The development of direct-seeded rice at high latitudes and altitudes depends on the seed germination ability of rice under a low-temperature environment. However, the lack of regulatory genes for low-temperature germination has severely limited the application of genetics in improving the breeds. Here, we used cultivars DN430 and DF104 with significantly different low-temperature germination (LTG) and 460 F2:3 progeny derived from them to identify LTG regulators by combining QTL-sequencing, linkage mapping, and RNA-sequencing. The QTL-sequencing mapped qLTG9 within a physical interval of 3.4 Mb. In addition, we used 10 Kompetitive allele-specific PCR (KASP) markers provided by the two parents, and qLTG9 was optimized from 3.4 Mb to a physical interval of 397.9 kb and accounted for 20.4% of the phenotypic variation. RNA-sequencing identified qLTG9 as eight candidate genes with significantly different expression within the 397.9 kb interval, six of which possessed SNPs on the promoter and coding regions. Quantitative reverse transcription-polymerase chain reaction (qRT-PCR) completely validated the results of these six genes in RNA-sequencing. Subsequently, six non-synonymous SNPs were designed using variants in the coding region of these six candidates. Genotypic analysis of these SNPs in 60 individuals with extreme phenotypes indicated these SNPs determined the differences in cold tolerance between parents. The six candidate genes of qLTG9 and the six KASP markers could be used together for marker-assisted breeding to improve LTG.

Integrated genomic approaches delineate a novel role of ROP1 ENHANCER1 in controlling seed protein content of chickpea.

Utilizing a combinatorial approach of quantitative trait locus (QTL)-Seq and candidate gene-based association mapping, the QTLs and genes responsible for seed protein content (SPC), a major quality trait in chickpea, were identified. Whole genome re-sequencing based QTL-Seq analysis of bulked recombinant inbred lines from a mapping population contrasting for SPC led to the identification of two QTLs [0.94 Mb on Linkage Group (LG)5 and 1.16 Mb on LG6] encompassing three SNPs, displaying the highest ΔSNP index. These highly significant SNPs and their associated genes were validated in 211 chickpea mini-core accessions varying in SPC, revealing a tightly associated marker affecting CaREN1 (ROP1 ENHANCER1) and explaining a phenotypic variation of 23%. This SNP was subsequently converted into a cost effective allele-specific PCR-based marker that could be utilized for rapid screening of SPC during marker assisted breeding. Furthermore, in planta functional validation via knockdown of CaREN1 transcripts led to significant reduction in SPC of chickpea. This decrease in seed protein is likely due to disruption in the formation of CaREN1 protein complexes comprising chaperones, phosphopeptide-binding proteins, and GTPases that mediate folding, transport and accumulation of seed storage proteins, as indicated through affinity purification-mass spectrometry. Taken together, our data will expedite tailoring of chickpea cultivars with augmented SPC.

Deploying QTL-seq rapid identification and separation of the major QTLs of tassel branch number for fine-mapping in advanced maize populations.

The tassel competes with the ear for nutrients and shields the upper leaves, thereby reducing the yield of grain. The tassel branch number (TBN) is a pivotal determinant of tassel size, wherein the reduced TBN has the potential to enhance the transmission of light and reduce the consumption of nutrients, which should ultimately result in increased yield. Consequently, the TBN has emerged as a vital target trait in contemporary breeding programs that focus on compact maize varieties. In this study, QTL-seq technology and advanced population mapping were used to rapidly identify and dissect the major effects of the TBN on QTL. Advanced mapping populations (BC4F2 and BC4F3) were derived from the inbred lines 18-599 (8-11 TBN) and 3237 (0-1 TBN) through phenotypic recurrent selection. First, 13 genomic regions associated with the TBN were detected using quantitative trait locus (QTL)-seq and were located on chromosomes 2 and 5. Subsequently, validated loci within these regions were identified by QTL-seq. Three QTLs for TBN were identified in the BC4F2 populations by traditional QTL mapping, with each QTL explaining the phenotypic variation of 6.13-18.17%. In addition, for the major QTL (qTBN2-2 and qTBN5-1), residual heterozygous lines (RHLs) were developed from the BC4F2 population. These two major QTLs were verified in the RHLs by QTL mapping, with the phenotypic variation explained (PVE) of 21.57% and 30.75%, respectively. Near-isogenic lines (NILs) of qTBN2-2 and qTBN5-1 were constructed. There were significant differences between the NILs in TBN. These results will enhance our understanding of the genetic basis of TBN and provide a solid foundation for the fine-mapping of TBN. Supplementary Information: The online version contains supplementary material available at 10.1007/s11032-023-01431-y.

Mapping of qChalk1 controlling grain chalkiness in japonica rice.

BACKGROUND: Rice grain chalkiness is an undesirable characteristic that affects grain quality. The aim of this study was to map QTLs controlling grain chalkiness in japonica rice. METHODS AND RESULTS: In this study, two japonica rice cultivars with similar grain shapes but different grain chalkiness rates were crossed and the F2 and BC1F2 populations were subjected to QTL-seq analysis to map the QTLs controlling the grain chalkiness rate. QTL-seq analysis revealed SNP index differences on chromosome 1 in both of the segregating populations. Using polymorphic markers between the two parents, QTL mapping was conducted on 213 individual plants in the BC1F2 population. QTL mapping confined a QTL controlling grain chalkiness, qChalk1, to a 1.1 Mb genomic region on chromosome 1. qChalk1 explained 19.7% of the phenotypic variation. CONCLUSION: A QTL controlling grain chalkiness qChalk1 was detected in both F2 and BC1F2 segregating populations by QTL-Seq and QTL mapping methods. This result would be helpful for further cloning of the genes controlling grain chalkiness in japonica rice.

Identification of a Major-Effect Quantitative Trait Loci Associated with Begomovirus Resistance in Cucurbita moschata.

Begomoviruses, viz. squash leaf curl China virus and tomato leaf curl New Delhi virus causative diseases are major concerns of quantitative and qualitative losses in pumpkin (Cucurbita moschata) worldwide. Punjab Agricultural University (PAU) in India has identified a resistant source (PVR-1343) against mixed infection (MI-Sq/To) of these begomoviruses. Introgression of resistance in diverse genetic backgrounds requires the identification of quantitative trait loci (QTLs) associated with MI-Sq/To resistance. Phenotyping of 229 F2:3 progenies derived from the PVR-1343 × P-135 cross revealed digenic recessive inheritance against MI-Sq/To resistance in PVR-1343. To identify the genomic region, resistant and susceptible bulks were subjected to whole-genome resequencing along with their parents. The whole-genome resequence analysis of parents and bulks using QTLseq/QTLseqr approaches identified an overlapping 1.52 Mb region on chromosome 7 (qMI-Sq/To7.1), while chromosomal region spanning 0.87 Mb on chromosome17 (qMI-Sq/To17.1) was additionally identified by QTLseqr. However, the highest peak value on chromosome 7 with three algorithms {G', ∆(SNP-index) and -log10 (P value)} highlighted the major contribution of qMI-Sq/To7.1 in MI-Sq/To resistance. Nine polymorphic SNPs identified within the highly significant qMI-Sq/To7.1 region were converted into KASP markers. KASP genotyping of F2 individuals narrowed down the qMI-Sq/To7.1 interval to 103 kb region flanked by two markers, Cmo3914729 and Cmo4018182, which contained 16 annotated genes and accounted for 59.84% of phenotypic variation. The Cmo4018182 KASP marker accurately predicted disease reaction in 91% of diverse Cucurbita genotypes and showed nonsynonym substitutions in the coding region of putative candidate SYNTAXIN-121 gene. These findings pave the way for marker-assisted breeding and elucidating the underlying mechanism of begomovirus resistance in C. moschata.

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 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 P. striiformis f. sp. tritici 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 the 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 P. striiformis f. sp. tritici 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 cosegregating 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.

QTL and candidate gene identification of the node of the first fruiting branch (NFFB) by QTL-seq in upland cotton (Gossypium hirsutum L.).

BACKGROUND: The node of the first fruiting branch (NFFB) is an important precocious trait in cotton. Many studies have been conducted on the localization of quantitative trait loci (QTLs) and genes related to fiber quality and yield, but there has been little attention to traits related to early maturity, especially the NFFB, in cotton. RESULTS: To identify the QTL associated with the NFFB in cotton, a BC4F2 population comprising 278 individual plants was constructed. The parents and two DNA bulks for high and low NFFB were whole genome sequenced, and 243.8 Gb of clean nucleotide data were generated. A total of 449,302 polymorphic SNPs and 135,353 Indels between two bulks were identified for QTL-seq. Seventeen QTLs were detected and localized on 11 chromosomes in the cotton genome, among which two QTLs (qNFFB-Dt2-1 and qNFFB-Dt3-3) were located in hotspots. Two candidate genes (GhAPL and GhHDA5) related to the NFFB were identified using quantitative real-time PCR (qRT-PCR) and virus-induced gene silencing (VIGS) experiments in this study. Both genes exhibited higher expression levels in the early-maturing cotton material RIL182 during flower bud differentiation, and the silencing of GhAPL and GhHDA5 delayed the flowering time and increased the NFFB compared to those of VA plants in cotton. CONCLUSIONS: Our study preliminarily found that GhAPL and GhHDA5 are related to the early maturity in cotton. The findings provide a basis for the further functional verification of candidate genes related to the NFFB and contribute to the study of early maturity in cotton.

Identification of a major-effect QTL associated with pre-harvest sprouting in cucumber (Cucumis sativus L.) using the QTL-seq method.

BACKGROUND: Cucumber (Cucumis sativus L.) is cultivated worldwide, and it is essential to produce enough high-quality seeds to meet demand. Pre-harvest sprouting (PHS) in cucumber is a critical problem and causes serious damage to seed production and quality. Nevertheless, the genetic basis and molecular mechanisms underlying cucumber PHS remain unclear. QTL-seq is an efficient approach for rapid quantitative trait loci (QTL) identification that simultaneously takes advantage of bulked-segregant analysis (BSA) and whole-genome resequencing. In the present research, QTL-seq analysis was performed to identify QTLs associated with PHS in cucumber using an F2 segregating population. RESULTS: Two QTLs that spanned 7.3 Mb on Chromosome 4 and 0.15 Mb on Chromosome 5 were identified by QTL-seq and named qPHS4.1 and qPHS5.1, respectively. Subsequently, SNP and InDel markers selected from the candidate regions were used to refine the intervals using the extended F2 populations grown in the 2016 and 2017 seasons. Finally, qPHS4.1 was narrowed to 0.53 Mb on chromosome 4 flanked by the markers SNP-16 and SNP-24 and was found to explain 19-22% of the phenotypic variation in cucumber PHS. These results reveal that qPHS4.1 is a major-effect QTL associated with PHS in cucumber. Based on gene annotations and qRT-PCR expression analyses, Csa4G622760 and Csa4G622800 were proposed as the candidate genes. CONCLUSIONS: These results provide novel insights into the genetic mechanism controlling PHS in cucumber and highlight the potential for marker-assisted selection of PHS resistance breeding.

Combining QTL-seq and linkage mapping to fine map a candidate gene in qCTS6 for cold tolerance at the seedling stage in rice.

BACKGROUND: Cold stress caused by low temperatures is an important factor restricting rice production. Identification of cold-tolerance genes that can stably express in cold environments is crucial for molecular rice breeding. RESULTS: In this study, we employed high-throughput quantitative trait locus sequencing (QTL-seq) analyses in a 460-individual F2:3 mapping population to identify major QTL genomic regions governing cold tolerance at the seedling stage in rice. A novel major QTL (qCTS6) controlling the survival rate (SR) under low-temperature conditions of 9°C/10 days was mapped on the 2.60-Mb interval on chromosome 6. Twenty-seven single-nucleotide polymorphism (SNP) markers were designed for the qCST6 region based on re-sequencing data, and local QTL mapping was conducted using traditional linkage analysis. Eventually, we mapped qCTS6 to a 96.6-kb region containing 13 annotated genes, of which seven predicted genes contained 13 non-synonymous SNP loci. Quantitative reverse transcription PCR analysis revealed that only Os06g0719500, an OsbZIP54 transcription factor, was strongly induced by cold stress. Haplotype analysis confirmed that +376 bp (T>A) in the OsbZIP54 coding region played a key role in regulating cold tolerance in rice. CONCLUSION: We identified OsbZIP54 as a novel regulatory gene associated with rice cold-responsive traits, with its Dongfu-104 allele showing specific cold-induction expression serving as an important molecular variation for rice improvement. This result is expected to further exploration of the genetic mechanism of rice cold tolerance at the seedling stage and improve cold tolerance in rice varieties by marker-assisted selection.

Genetic mapping high protein content QTL from soybean 'Nanxiadou 25' and candidate gene analysis.

BACKGROUND: Soybean is a globally important legume crop that provides a primary source of high-quality vegetable protein and oil. Seed protein content (SPC) is a valuable quality trait controlled by multiple genes in soybean. RESULTS: In this study, we performed quantitative trait loci (QTL) mapping, QTL-seq, and RNA sequencing (RNA-seq) to reveal the genes controlling protein content in the soybean by using the high protein content variety Nanxiadou 25. A total of 50 QTL for SPC distributed on 14 chromosomes except chromosomes 4, 12, 14, 17, 18, and 19 were identified by QTL mapping using 178 recombinant inbred lines (RILs). Among these QTL, the major QTL qSPC_20-1 and qSPC_20-2 on chromosome 20 were repeatedly detected across six tested environments, corresponding to the location of the major QTL detected using whole-genome sequencing-based QTL-seq. 329 candidate DEGs were obtained within the QTL region of qSPC_20-1 and qSPC_20-2 via gene expression profile analysis. Nine of which were associated with SPC, potentially representing candidate genes. Clone sequencing results showed that different single nucleotide polymorphisms (SNPs) and indels between high and low protein genotypes in Glyma.20G088000 and Glyma.16G066600 may be the cause of changes in this trait. CONCLUSIONS: These results provide the basis for research on candidate genes and marker-assisted selection (MAS) in soybean breeding for seed protein content.

Combining QTL-seq and linkage mapping to uncover the genetic basis of single vs. paired spikelets in the advanced populations of two-ranked maize×teosinte.

BACKGROUND: Teosinte ear bears single spikelet, whereas maize ear bears paired spikelets, doubling the number of grains in each cupulate during maize domestication. In the past 20 years, genetic analysis of single vs. paired spikelets (PEDS) has been stagnant. A better understanding of genetic basis of PEDS could help fine mapping of quantitative trait loci (QTL) and cloning of genes. RESULTS: In this study, the advanced mapping populations (BC3F2 and BC4F2) of maize × teosinte were developed by phenotypic recurrent selection. Four genomic regions associated with PEDS were detected using QTL-seq, located on 194.64-299.52 Mb, 0-162.80 Mb, 12.82-97.17 Mb, and 125.06-157.01 Mb of chromosomes 1, 3, 6, and 8, respectively. Five QTL for PEDS were identified in the regions of QTL-seq using traditional QTL mapping. Each QTL explained 1.12-38.05% of the phenotypic variance (PVE); notably, QTL qPEDS3.1 with the average PVE of 35.29% was identified in all tests. Moreover, 14 epistatic QTL were detected, with the total PVE of 47.57-66.81% in each test. The QTL qPEDS3.1 overlapped with, or was close to, one locus of 7 epistatic QTL. Near-isogenic lines (NILs) of QTL qPEDS1.1, qPEDS3.1, qPEDS6.1, and qPEDS8.1 were constructed. All individuals of NIL-qPEDS6.1(MT1) and NIL-qPEDS8.1(MT1) showed paired spikelets (PEDS = 0), but the flowering time was 7 days shorter in the NIL-qPEDS8.1(MT1). The ratio of plants with PEDS > 0 was low (1/18 to 3/18) in the NIL-qPEDS1.1(MT1) and NIL-qPEDS3.1(MT1), maybe due to the epistatic effect. CONCLUSION: Our results suggested that major QTL, minor QTL, epistasis and photoperiod were associated with the variation of PEDS, which help us better understand the genetic basis of PEDS and provide a genetic resource for fine mapping of QTL.

Targeted amplicon sequencing + next-generation sequencing-based bulked segregant analysis identified genetic loci associated with preharvest sprouting tolerance in common buckwheat (Fagopyrum esculentum).

BACKGROUND: Common buckwheat (2n = 2x = 16) is an outcrossing pseudocereal whose seeds contain abundant nutrients and potential antioxidants. As these beneficial compounds are damaged by preharvest sprouting (PHS) and PHS is likely to increase with global warming, it is important to find efficient ways to develop new PHS-tolerant lines. However, genetic loci and selection markers associated with PHS in buckwheat have not been reported. RESULTS: By next-generation sequencing (NGS) of whole-genome of parental lines, we developed a genome-wide set of 300 markers. By NGS- based bulked segregant analysis (NGS-BSA), we developed 100 markers linked to PHS tolerance. To confirm the effectiveness of marker development from NGS-BSA data, we developed 100 markers linked to the self-compatibility (SC) trait from previous NGS-BSA data. Using these markers, we developed genetic maps with AmpliSeq technology, which can quickly detect polymorphisms by amplicon-based multiplex targeted NGS, and performed quantitative trait locus (QTL) analysis for PHS tolerance in combination with NGS-BSA. QTL analysis detected two major and two minor QTLs for PHS tolerance in a segregating population developed from a cross between the PHS-tolerant 'Kyukei 29' and the self-compatible susceptible 'Kyukei SC7'. We found different major and minor QTLs in other segregating populations developed from the PHS-tolerant lines 'Kyukei 28' and 'NARO-FE-1'. Candidate markers linked to PHS developed by NGS-BSA were located near these QTL regions. We also investigated the effectiveness of markers linked to these QTLs for selection of PHS-tolerant lines among other segregating populations. CONCLUSIONS: We efficiently developed genetic maps using a method combined with AmpliSeq technology and NGS-BSA, and detected QTLs associated with preharvest sprouting tolerance in common buckwheat. This is the first report to identify QTLs for PHS tolerance in buckwheat. Our marker development system will accelerate genetic research and breeding in common buckwheat.

Polyploid QTL-seq towards rapid development of tightly linked DNA markers for potato and sweetpotato breeding through whole-genome resequencing.

Potato (Solanum tuberosum L.) and sweetpotato (Ipomoea batatas L.), which are nutritionally and commercially important tuberous crops, possess a perplexing heredity because of their autopolyploid genomes. To reduce cross-breeding efforts for selecting superior cultivars from progenies with innumerable combinations of traits, DNA markers tightly linked to agronomical traits are required. To develop DNA markers, we developed a method for quantitative trait loci (QTL) mapping using whole-genome next-generation sequencing (NGS) in autopolyploid crops. To apply the NGS-based bulked segregant method, QTL-seq was modified. (1) Single parent-specific simplex (unique for one homologous chromosome) single-nucleotide polymorphisms (SNPs), which present a simple segregation ratio in the progenies, were exploited by filtering SNPs by SNP index (allele frequency). (2) Clusters of SNPs, which were inherited unevenly between bulked progenies with opposite phenotypes, especially those with an SNP index of 0 for the bulk that did not display the phenotypes of interest, were explored. These modifications allowed for separate tracking of alleles located on each of the multiple homologous chromosomes. By applying this method, clusters of SNPs linked to the potato cyst nematode resistance H1 gene and storage root anthocyanin (AN) content were identified in tetraploid potato and hexaploid sweetpotato, respectively, and completely linked DNA markers were developed at the site of the presented SNPs. Thus, polyploid QTL-seq is a versatile method that is free from specialized manipulation for sequencing and construction of elaborate linkage maps and facilitates rapid development of tightly linked DNA markers in autopolyploid crops, such as potato and sweetpotato.