Mapping quantitative trait loci in selected breeding populations: A segregation distortion approach

Mapping quantitative trait loci for seminal root angle in a selected durum wheat population

Seminal root angle (SRA) is an important root architectural trait associated with drought adaptation in cereal crops. To date, all attempts to dissect the genetic architecture of SRA in durum wheat (Triticum durum Desf.) have used large association panels or structured mapping populations. Identifying changes in allele frequency generated by selection provides an alternative genetic mapping approach that can increase the power and precision of QTL detection. This study aimed to map quantitative trait loci (QTL) for SRA by genotyping durum lines created through divergent selection using a combination of marker-assisted selection (MAS) for the major SRA QTL (qSRA-6A) and phenotypic selection for SRA over multiple generations. The created 11 lines (BC1F2:5) were genotyped with genome-wide single-nucleotide polymorphism (SNP) markers to map QTL by identifying markers that displayed segregation distortion significantly different from the Mendelian expectation. QTL regions were further assessed in an independent validation population to confirm their associations with SRA. The experiment revealed 14 genomic regions under selection, 12 of which have not previously been reported for SRA. Five regions, including qSRA-6A, were confirmed in the validation population. The genomic regions identified in this study indicate that the genetic control of SRA is more complex than previously anticipated. Our study demonstrates that selection mapping is a powerful approach to complement genome-wide association studies for QTL detection. Moreover, the verification of qSRA-6A in an elite genetic background highlights the potential for MAS, although it is necessary to combine additional QTL to develop new cultivars with extreme SRA phenotypes.

Secondary metabolite fingerprinting, anti-pathogenic activity, elite chemotype selection and conservation of Curcuma caesia- an ethnomedicinally underutilized species

Curcuma caesia Roxb. is an ethnomedicinally important, essential oil (EO) yielding aromatic plant. A total of twelve accessions of this plant rhizome were collected from six different agro-climatic zones of West Bengal, India and evaluated for their antimicrobial activities against eight disease-causing, multi-drug-resistant pathogenic strains of urinary-tract infection and respiratory-tract infection. The EO and extracts demonstrated antibacterial activity, with the highest inhibition zone of 18.00 ± 0.08 and 17.50 ± 0.14 mm against Klebsiella pneumoniae by accession 06, even where all the broad-spectrum antibiotics failed to respond. In this study, we employed high-performance thin-layer chromatography (HPTLC) to quantify curcumin, the primary secondary metabolite of C. caesia, and the highest 0.228 mg/gm of curcumin resulted from accession 06. Hence, on the basis of all aspects, accession 06 was identified as the elite chemotype among all twelve accessions. The chemical profiling of EO from accession 06 was done using gas chromatography-mass spectroscopy (GC-MS). Conceivably, about 13 medicinally significant compounds were detected. As this plant species is seasonal and has difficulties in conventional breeding due to dormancy, it must be conserved through in vitro tissue culture for a steady supply throughout the year in massive amounts for agricultural demand. A maximum number of 19.28 ± 0.37 shoots has been obtained in MS medium fortified with 6-Benzylaminopurine, Kinetin, and Naphthalene acetic acid. The genetic uniformity of the plants has been studied through Start Codon Targeted Polymorphism. Therefore, this study must help meet the need for essential phytoactive compounds through a simple, validated, and reproducible plant tissue culture protocol throughout the year.

Unveiling the Genetic Basis Underlying Rice Anther Culturability via Segregation Distortion Analysis in Doubled Haploid Population

Anther culture (AC) is a valuable technique in rice breeding. However, the genetic mechanisms underlying anther culturability remain elusive, which has hindered its widespread adoption in rice breeding programs. During AC, microspores carrying favorable alleles for AC are selectively regenerated, leading to segregation distortion (SD) of chromosomal regions linked to these alleles in the doubled haploid (DH) population. Using the AC method, a DH population was generated from the japonica hybrid rice Shenyou 26. A genetic map consisting of 470 SNPs was constructed using this DH population, and SD analysis was performed at both the single- and two-locus levels to dissect the genetic basis underlying anther culturability. Five segregation distortion loci (SDLs) potentially linked to anther culturability were identified. Among these, SDL5 exhibited an overrepresentation of alleles from the female parent, while SDL1.1, SDL1.2, SDL2, and SDL7 displayed an overrepresentation of alleles from the male parent. Furthermore, six pairs of epistatic interactions (EPIs) that influenced two-locus SDs in the DH population were discovered. A cluster of genetic loci, associated with EPI-1, EPI-3, EPI-4, and EPI-5, overlapped with SDL1.1, indicating that the SDL1.1 locus may play a role in regulating anther culturability via both additive and epistatic mechanisms. These findings provide valuable insights into the genetic control of anther culturability in rice and lay the foundation for future research focused on identifying the causal genes associated with anther culturability.

Efficient construction of a linkage map and haplotypes for Mentha suaveolens using sequence capture

The sustainability of many crops is hindered by the lack of genomic resources and a poor understanding of natural genetic diversity. Particularly, application of modern breeding requires high-density linkage maps that are integrated into a highly contiguous reference genome. Here, we present a rapid method for deriving haplotypes and developing linkage maps, and its application to Mentha suaveolens, one of the diploid progenitors of cultivated mints. Using sequence-capture via DNA hybridization to target single nucleotide polymorphisms (SNPs), we successfully genotyped ∼5000 SNPs within the genome of >400 individuals derived from a self cross. After stringent quality control, and identification of nonredundant SNPs, 1919 informative SNPs were retained for linkage map construction. The resulting linkage map defined a total genetic space of 942.17 cM divided among 12 linkage groups, ranging from 56.32 to 122.61 cM in length. The linkage map is in good agreement with pseudomolecules from our preliminary genome assembly, proving this resource effective for the correction and validation of the reference genome. We discuss the advantages of this method for the rapid creation of linkage maps.

QTL Mapping and Favorable Allele Mining of Nitrogen Deficiency Tolerance Using an Interconnected Breeding Population in Rice

Nitrogen is one of the most important nutrients for rice growth and development. Breeding of nitrogen deficiency tolerance (NDT) variety is considered to be the most economic measure to solve the constrain of low nitrogen stress on grain yield in rice. An interconnected breeding (IB) population of 497 lines developed using Huanghuazhan (HHZ) as the recurrent parent and eight elite lines as the donor parents were tested for five traits including grain yield, biomass, harvest index, thousand grain weight, and spikelet fertility under two nitrogen treatments in three growing seasons. Association analysis using 7,388 bins generated by sequencing identified a total of 14, 14, and 12 QTLs for the five traits under low nitrogen (LN), normal nitrogen (NN), and LN/NN conditions, respectively, across three seasons. Favorable alleles were dissected for the 40 QTLs at the 10 NDT regions, and OM1723 was considered as the most important parent with the highest frequency of favorable alleles contributing to NDT-related traits. Six superior lines all showed significantly higher GY in LN environments and similar GY under NN environments except for H10. Substitution mapping using near-isogenic introgression lines delimited the qTGW2-1, which was identified on chromosome 2 under LN, NN, and LN/NN conditions into two QTLs, which were located in the two regions of about 200 and 350 kb with different favorable alleles. The bins 16, 1301, 1465, 1486, 3464, and 6249 harbored the QTLs for NDT detected in this study, and the QTLs/genes previously identified for NDT or nitrogen use efficiency (NUE) could be used for enhancing NDT and NUE by marker-assisted selection (MAS).

Mapping of QTL and identification of candidate genes conferring spontaneous haploid genome doubling in maize (Zea mays L.)

In vivo doubled haploid (DH) technology is widely used in commercial maize (Zea mays L.) breeding. Haploid genome doubling is a critical step in DH breeding. In this study, inbred lines GF1 (0.65), GF3(0.29), and GF5 (0) with high, moderate, and poor spontaneous haploid genome doubling (SHGD), respectively, were selected to develop mapping populations for SHGD. Three QTL, qshgd1, qshgd2, and qshgd3, related to SHGD were identified by selective genotyping. With the exception of qshgd3, the source of haploid genome doubling alleles were derived from GF1. Furthermore, RNA-Seq was conducted to identify putative candidate genes between GF1 and GF5 within the qshgd1 region. A differentially expressed formin-like protein 5 transcript was identified within the qshgd1 region.

Molecular mapping of major QTL conferring resistance to orange wheat blossom midge (Sitodiplosis mosellana) in Chinese wheat varieties with selective populations

Two novel midge resistance QTL were mapped to a 4.9-Mb interval on chromosome arm 4AL based on the genetic maps constructed with SNP markers. Orange wheat blossom midge (OWBM) is a devastating insect pest affecting wheat production. In order to detect OWBM resistance genes and quantitative trait loci (QTL) for wheat breeding, two recombinant inbred line (RIL) populations were established and used for molecular mapping. A total of seven QTL were detected on chromosomes 2D, 4A, 4D and 7D, respectively, of which positive alleles were all from the resistant parents except for the QTL on 7D. Two stable QTL (QSm.hbau-4A.2-1 and QSm.hbau-4A.2-2) were detected in both populations with the LOD scores ranging from 5.58 to 29.22 under all three environments, and they explained a combined phenotypic variation of 24.4-44.8%. These two novel QTL were mapped to a 4.9-Mb physical interval. The single-nucleotide polymorphism (SNP) markers AX-109543456, AX-108942696 and AX-110928325 were closely linked to the QTL and could be used for marker-assisted selection (MAS) for OWBM resistance in wheat breeding programs.

Seed protein content and its relationships with agronomic traits in pigeonpea is controlled by both main and epistatic effects QTLs

The genetic architecture of seed protein content (SPC) and its relationships to agronomic traits in pigeonpea is poorly understood. Accordingly, five F₂ populations segregating for SPC and four agronomic traits (seed weight (SW), seed yield (SY), growth habit (GH) and days to first flowering (DFF)) were phenotyped and genotyped using genotyping-by-sequencing approach. Five high-density population-specific genetic maps were constructed with an average inter-marker distance of 1.6 to 3.5 cM, and subsequently, integrated into a consensus map with average marker spacing of 1.6 cM. Based on analysis of phenotyping data and genotyping data, 192 main effect QTLs (M-QTLs) with phenotypic variation explained (PVE) of 0.7 to 91.3% were detected for the five traits across the five populations. Major effect (PVE ≥ 10%) M-QTLs included 14 M-QTLs for SPC, 16 M-QTLs for SW, 17 M-QTLs for SY, 19 M-QTLs for GH and 24 M-QTLs for DFF. Also, 573 epistatic QTLs (E-QTLs) were detected with PVE ranging from 6.3 to 99.4% across traits and populations. Colocalization of M-QTLs and E-QTLs explained the genetic basis of the significant (P < 0.05) correlations of SPC with SW, SY, DFF and GH. The nature of genetic architecture of SPC and its relationship with agronomic traits suggest that genomics-assisted breeding targeting genome-wide variations would be effective for the simultaneous improvement of SPC and other important traits.

Transmission ratio distortion is frequent in Arabidopsis thaliana controlled crosses

The equal probability of transmission of alleles from either parent during sexual reproduction is a central tenet of genetics and evolutionary biology. Yet, there are many cases where this rule is violated. The preferential transmission of alleles or genotypes is termed transmission ratio distortion (TRD). Examples of TRD have been identified in many species, implying that they are universal, but the resolution of species-wide studies of TRD are limited. We have performed a species-wide screen for TRD in over 500 segregating F2 populations of Arabidopsis thaliana using pooled reduced-representation genome sequencing. TRD was evident in up to a quarter of surveyed populations. Most populations exhibited distortion at only one genomic region, with some regions being repeatedly affected in multiple populations. Our results begin to elucidate the species-level architecture of biased transmission of genetic material in A. thaliana, and serve as a springboard for future studies into the biological basis of TRD in this species.

QTL and QTL networks for cold tolerance at the reproductive stage detected using selective introgression in rice

Low temperature stress is one of the major abiotic stresses limiting the productivity of Geng (japonica) rice grown the temperate regions as well as in tropical high lands worldwide. To develop rice varieties with improved cold tolerance (CT) at the reproductive stage, 84 BC2 CT introgression lines (ILs) were developed from five populations through backcross breeding. These CT ILs plus 310 random ILs from the same BC populations were used for dissecting genetic networks underlying CT in rice by detecting QTLs and functional genetic units (FGUs) contributing to CT. Seventeen major QTLs for CT were identified using five selective introgression populations and the method of segregation distortion. Of them, three QTLs were confirmed using the random populations and seven others locate in the regions with previously reported CT QTLs/genes. Using multi-locus probability tests and linkage disequilibrium (LD) analyses, 46 functional genetic units (FGUs) (37 single loci and 9 association groups or AGs) distributed in 37 bins (~20%) across the rice genome for CT were detected. Together, each of the CT loci (bins) was detected in 1.7 populations, including 18 loci detected in two or more populations. Putative genetic networks (multi-locus structures) underlying CT were constructed based on strong non-random associations between or among donor alleles at the unlinked CT loci/FGUs identified in the CT ILs, suggesting the presence of strong epistasis among the detected CT loci. Our results demonstrated the power and usefulness of using selective introgression for simultaneous improvement and genetic dissection of complex traits such as CT in rice.

Genetic insights into underground responses to Fusarium graminearum infection in wheat

The ongoing global intensification of wheat production will likely be accompanied by a rising pressure of Fusarium diseases. While utmost attention was given to Fusarium head blight (FHB) belowground plant infections of the pathogen have largely been ignored. The current knowledge about the impact of soil borne Fusarium infection on plant performance and the underlying genetic mechanisms for resistance remain very limited. Here, we present the first large-scale investigation of Fusarium root rot (FRR) resistance using a diverse panel of 215 international wheat lines. We obtained data for a total of 21 resistance-related traits, including large-scale Real-time PCR experiments to quantify fungal spread. Association mapping and subsequent haplotype analyses discovered a number of highly conserved genomic regions associated with resistance, and revealed a significant effect of allele stacking on the stembase discoloration. Resistance alleles were accumulated in European winter wheat germplasm, implying indirect prior selection for improved FRR resistance in elite breeding programs. Our results give first insights into the genetic basis of FRR resistance in wheat and demonstrate how molecular parameters can successfully be explored in genomic prediction. Ongoing work will help to further improve our understanding of the complex interactions of genetic factors influencing FRR resistance.

Simultaneous Improvement and Genetic Dissection of Drought Tolerance Using Selected Breeding Populations of Rice

Drought is the most important factor limiting rice yield in most rainfed areas of Asia and Africa. Four large BC₂F₂ populations consisted of 3,200 individuals, which were derived from crosses between an elite Geng variety, Jigeng88, and four donors from three different countries, were screened and progeny tested under severe drought stress, resulting in the development of 72 introgression lines (ILs) with significantly improved yield compared to the recurrent parent Jigeng88. These DT ILs plus four random populations (without drought selection population) from the same crosses were evaluated in replicated trials under both drought stress and non-stress conditions in two environments, and characterized with simple sequence repeat (SSR) markers to understand how directional selection was operating on the genetic variation of DT of rice. Thirteen DT QTLs of large effect were identified based on the significant allelic and genotypic frequency shits in the DT ILs by using the joint segregation distortion method. The 13 QTLs were validated by the genotypic differences at individual QTL in the random populations. Putative genetic networks consisting of 30 loci in 29 functional genetic units underlying DT were detected by X² tests and non-random associations between or among DT loci in DT ILs from the four populations. Most large-effect DT QTLs were previously reported and located in the upstream of the genetic networks as putative regulators, and were either mapped to important regulatory genes for DT or drought responsiveness reported previously. In our study, five promising ILs with significantly improved yield were selected under both drought and normal irrigated conditions. The QTLs and their genetic networks underlying DT detected provided useful genetic information for further improving DT and yield using designed QTL pyramiding.

Genetic Dissection and Simultaneous Improvement of Drought and Low Nitrogen Tolerances by Designed QTL Pyramiding in Rice

Drought and low nitrogen are the most common abiotic stresses limiting rice productivity in the rainfed rice areas of Asia and Africa. Development and adoption of green super rice (GSR) varieties with greatly improved drought tolerance (DT) and low nitrogen tolerance (LNT) are the most efficient way to resolve this problem. In this study, using three sets of trait-specific introgression lines (ILs) in a Xian (indica) variety Huanghuazhan (HHZ) background, we identified nine DT-QTL and seven LNT-QTL by a segregation distortion approach and a genome-wide association study, respectively. Based on performances of DT and LNT and genotypes at the detected QTL, two ILs M79 and M387 with DT and LNT were selected for cross-making to validate the identified QTL and to develop DT and LNT rice lines by pyramiding two DT-QTL (qDT3.9 and qDT6.3) and two LNT-QTL (qGY1 and qSF8). Using four pairs of kompetitive allele specific PCR (KASP) SNP markers, we selected 66 F2 individuals with different combinations of the target DT- and LNT-QTL favorable alleles and they showed expected improvement in DT and/or LNT, which were further validated by the significant improvement in DT and/or LNT of their F3 progeny testing. Based on evaluation of pyramiding lines in F3 lines under drought, low nitrogen (LN) and normal conditions, four promising pyramiding lines having different QTL favorable alleles were selected, which showed significantly improved tolerances to drought and/or LN than HHZ and their IL parents. Our results demonstrated that trait-specific ILs could effectively connect QTL mapping and QTL pyramiding breeding, and designed QTL pyramiding (DQP) using ILs could be more effective in molecular rice breeding for complex quantitative traits.

The Application of Multi-Locus GWAS for the Detection of Salt-Tolerance Loci in Rice

Improving the salt-tolerance of direct-seeding rice at the seed germination stage is a major goal of breeders. Efficiently identifying salt tolerance loci will help researchers develop effective rice breeding strategies. In this study, six multi-locus genome-wide association studies (GWAS) methods (mrMLM, FASTmrMLM, FASTmrEMMA, pLARmEB, pKWmEB, and ISIS EM-BLASSO) were applied to identify quantitative trait nucleotides (QTNs) for the salt tolerance traits of 478 rice accessions with 162,529 SNPs at the seed germination stage. Among the 371 QTNs detected by the six methods, 56 were identified by at least three methods. Among these 56 QTNs, 12, 6, 7, 4, 13, 12, and 12 were found to be associated with SSI-GI, SSI-VI, SSI-MGT, SSI-IR-24h, SSI-IR-48h, SSI-GR-5d, and SSI-GR-10d, respectively. Additionally, 66 candidate genes were identified in the vicinity of the 56 QTNs, and two of these genes (LOC_Os01g45760 and LOC_Os10g04860) are involved in auxin biosynthesis according to the enriched GO terms and KEGG pathways. This information will be useful for identifying the genes responsible for rice salt tolerance. A comparison of the six methods revealed that ISIS EM-BLASSO identified the most co-detected QTNs and performed best, with the smallest residual errors and highest computing speed, followed by FASTmrMLM, pLARmEB, mrMLM, pKWmEB, and FASTmrEMMA. Although multi-locus GWAS methods are superior to single-locus GWAS methods, their utility for identifying QTNs may be enhanced by adding a bin analysis to the models or by developing a hybrid method that merges the results from different methods.

A Mixed Model Approach to Genome-Wide Association Studies for Selection Signatures, with Application to Mice Bred for Voluntary Exercise Behavior

Selection experiments and experimental evolution provide unique opportunities to study the genetics of adaptation because the target and intensity of selection are known relatively precisely. In contrast to natural selection, where populations are never strictly "replicated," experimental evolution routinely includes replicate lines so that selection signatures-genomic regions showing excessive differentiation between treatments-can be separated from possible founder effects, genetic drift, and multiple adaptive solutions. We developed a mouse model with four lines within a high running (HR) selection treatment and four nonselected controls (C). At generation 61, we sampled 10 mice of each line and used the Mega Mouse Universal Genotyping Array to obtain single nucleotide polymorphism (SNP) data for 25,318 SNPs for each individual. Using an advanced mixed model procedure developed in this study, we identified 152 markers that were significantly different in frequency between the two selection treatments. They occurred on all chromosomes except 1, 2, 8, 13, and 19, and showed a variety of patterns in terms of fixation (or the lack thereof) in the four HR and four C lines. Importantly, none were fixed for alternative alleles between the two selection treatments. The current state-of-the-art regularized F test applied after pooling DNA samples for each line failed to detect any markers. We conclude that when SNP or sequence data are available from individuals, the mixed model methodology is recommended for selection signature detection. As sequencing at the individual level becomes increasingly feasible, the new methodology may be routinely applied for detection of selection.

QTL mapping for haploid male fertility by a segregation distortion method and fine mapping of a key QTL qhmf4 in maize

Four QTL related to haploid male fertility were detected by a segregation distortion method and the key QTL qhmf4 was fine mapped to an interval of ~800 kb. Doubled haploid (DH) technology enables rapid development of homozygous lines in maize breeding programs. However, haploid genome doubling is a bottleneck for the commercialization of DH technology and is limited by haploid male fertility (HMF). This is the first study reporting the quantitative trait locus (QTL) analysis of HMF in maize. Four QTL, qhmf1, qhmf2, qhmf3, and qhmf4, controlling HMF have been identified by segregation distortion (SD) loci detection in the selected haploid population derived from 'Yu87-1/Zheng58'. Three loci, qhmf1, qhmf2, and qhmf4, were also detected in the selected haploid population derived from '4F1/Zheng58'. The QTL qhmf4 showed the strongest SD in both haploid populations. Based on the sequence information of 'Yu87-1' and 'Zheng58', thirteen markers being polymorphic between the two lines were developed to saturate the qhmf4 region. A total of 8168 H₁BC₂ (haploid backcross generation) plants produced from 'Yu87-1' and 'Zheng58' were screened for recombinants. All the 48 recombinants were backcrossed to 'Zheng58' to develop H₁BC₃ progeny. The heterozygous H₁BC₃ individuals were crossed with CAU5 to induce haploids. In each H₁BC₃ progeny, haploids were genotyped and evaluated for anther emergence score (AES). Significant (or no significant) difference (P < 0.05) between haploids with or without 'Yu87-1' donor segment indicated presence or absence of qhmf4 in the donor segment. The analysis of the 48 recombinants narrowed the qhmf4 locus down to an ~800 kb interval flanked by markers IND166 and IND1668.

Exploring the genetic characteristics of 93-11 and Nipponbare recombination inbred lines based on the GoldenGate SNP assay

Understanding genetic characteristics in rice populations will facilitate exploring evolutionary mechanisms and gene cloning. Numerous molecular markers have been utilized in linkage map construction and quantitative trait locus (QTL) mappings. However, segregation-distorted markers were rarely considered, which prevented understanding genetic characteristics in many populations. In this study, we designed a 384-marker GoldenGate SNP array to genotype 283 recombination inbred lines (RILs) derived from 93-11 and Nipponbare Oryza sativa crosses. Using 294 markers that were highly polymorphic between parents, a linkage map with a total genetic distance of 1,583.2 cM was constructed, including 231 segregation-distorted markers. This linkage map was consistent with maps generated by other methods in previous studies. In total, 85 significant quantitative trait loci (QTLs) with phenotypic variation explained (PVE) values≥5% were identified. Among them, 34 QTLs were overlapped with reported genes/QTLs relevant to corresponding traits, and 17 QTLs were overlapped with reported sterility-related genes/QTLs. Our study provides evidence that segregation-distorted markers can be used in linkage map construction and QTL mapping. Moreover, genetic information resulting from this study will help us to understand recombination events and segregation distortion. Furthermore, this study will facilitate gene cloning and understanding mechanism of inter-subspecies hybrid sterility and correlations with important agronomic traits in rice.

Genetic dissection of heterosis using epistatic association mapping in a partial NCII mating design

Heterosis refers to the phenomenon in which an F₁ hybrid exhibits enhanced growth or agronomic performance. However, previous theoretical studies on heterosis have been based on bi-parental segregating populations instead of F₁ hybrids. To understand the genetic basis of heterosis, here we used a subset of F₁ hybrids, named a partial North Carolina II design, to perform association mapping for dependent variables: original trait value, general combining ability (GCA), specific combining ability (SCA) and mid-parental heterosis (MPH). Our models jointly fitted all the additive, dominance and epistatic effects. The analyses resulted in several important findings: 1) Main components are additive and additive-by-additive effects for GCA and dominance-related effects for SCA and MPH, and additive-by-dominant effect for MPH was partly identified as additive effect; 2) the ranking of factors affecting heterosis was dominance > dominance-by-dominance > over-dominance > complete dominance; and 3) increasing the proportion of F₁ hybrids in the population could significantly increase the power to detect dominance-related effects, and slightly reduce the power to detect additive and additive-by-additive effects. Analyses of cotton and rapeseed datasets showed that more additive-by-additive QTL were detected from GCA than from trait phenotype, and fewer QTL were from MPH than from other dependent variables.