Image_1_Detection of Breeding-Relevant Fruit Cracking and Fruit Firmness Quantitative Trait Loci in Sweet Cherry via Pedigree-Based and Genome-Wide Association Approaches.JPEG

Breeding for decreased fruit cracking incidence and increased fruit firmness in sweet cherry creates an attractive alternative to variable results from cultural management practices. DNA-informed breeding increases its efficiency, yet upstream research is needed to identify the genomic regions associated with the trait variation of a breeding-relevant magnitude, as well as to identify the parental sources of favorable alleles. The objectives of this research were to identify the quantitative trait loci (QTLs) associated with fruit cracking incidence and firmness, estimate the effects of single nucleotide polymorphism (SNP) haplotypes at the detected QTLs, and identify the ancestral source(s) of functional haplotypes. Fruit cracking incidence and firmness were evaluated for multiple years on 259 unselected seedlings representing 22 important breeding parents. Phenotypic data, in conjunction with genome-wide genotypic data from the RosBREED cherry 6K SNP array, were used in the QTL analysis performed via Pedigree-Based Analysis using the FlexQTL™ software, supplemented by a Genome-Wide Association Study using the BLINK software. Haplotype analysis was conducted on the QTLs to identify the functional SNP haplotypes and estimate their phenotypic effects, and the haplotypes were tracked through the pedigree. Four QTLs (two per trait) were consistent across the years and/or both analysis methods and validated the previously reported QTLs. qCrack-LG1.1m (the label given to a consistent QTL for cracking incidence on chromosome 1) explained 2–15.1% of the phenotypic variance, while qCrack-LG5.1m, qFirm-LG1.2m, and qFirm-LG3.2m explained 7.6–13.8, 8.8–21.8, and 1.7–10.1% of the phenotypic variance, respectively. At each QTL, at least two SNP haplotypes had significant effects and were considered putative functional SNP haplotypes. Putative low-cracking SNP haplotypes were tracked to an unnamed parent of ‘Emperor Francis’ and ‘Schmidt’ and unnamed parents of ‘Napoleon’ and ‘Hedelfingen,’ among others, and putative high-firmness haplotypes were tracked to an unnamed parent of ‘Emperor Francis’ and ‘Schmidt,’ an unnamed grandparent of ‘Black Republican,’ ‘Rube,’ and an unknown parent of ‘Napoleon.’ These four stable QTLs can now be targeted for DNA test development, with the goal of translating information discovered here into usable tools to aid in breeding decisions.</p

Image_3_Detection of Breeding-Relevant Fruit Cracking and Fruit Firmness Quantitative Trait Loci in Sweet Cherry via Pedigree-Based and Genome-Wide Association Approaches.JPEG

Breeding for decreased fruit cracking incidence and increased fruit firmness in sweet cherry creates an attractive alternative to variable results from cultural management practices. DNA-informed breeding increases its efficiency, yet upstream research is needed to identify the genomic regions associated with the trait variation of a breeding-relevant magnitude, as well as to identify the parental sources of favorable alleles. The objectives of this research were to identify the quantitative trait loci (QTLs) associated with fruit cracking incidence and firmness, estimate the effects of single nucleotide polymorphism (SNP) haplotypes at the detected QTLs, and identify the ancestral source(s) of functional haplotypes. Fruit cracking incidence and firmness were evaluated for multiple years on 259 unselected seedlings representing 22 important breeding parents. Phenotypic data, in conjunction with genome-wide genotypic data from the RosBREED cherry 6K SNP array, were used in the QTL analysis performed via Pedigree-Based Analysis using the FlexQTL™ software, supplemented by a Genome-Wide Association Study using the BLINK software. Haplotype analysis was conducted on the QTLs to identify the functional SNP haplotypes and estimate their phenotypic effects, and the haplotypes were tracked through the pedigree. Four QTLs (two per trait) were consistent across the years and/or both analysis methods and validated the previously reported QTLs. qCrack-LG1.1m (the label given to a consistent QTL for cracking incidence on chromosome 1) explained 2–15.1% of the phenotypic variance, while qCrack-LG5.1m, qFirm-LG1.2m, and qFirm-LG3.2m explained 7.6–13.8, 8.8–21.8, and 1.7–10.1% of the phenotypic variance, respectively. At each QTL, at least two SNP haplotypes had significant effects and were considered putative functional SNP haplotypes. Putative low-cracking SNP haplotypes were tracked to an unnamed parent of ‘Emperor Francis’ and ‘Schmidt’ and unnamed parents of ‘Napoleon’ and ‘Hedelfingen,’ among others, and putative high-firmness haplotypes were tracked to an unnamed parent of ‘Emperor Francis’ and ‘Schmidt,’ an unnamed grandparent of ‘Black Republican,’ ‘Rube,’ and an unknown parent of ‘Napoleon.’ These four stable QTLs can now be targeted for DNA test development, with the goal of translating information discovered here into usable tools to aid in breeding decisions.</p

Table_1_Detection of Breeding-Relevant Fruit Cracking and Fruit Firmness Quantitative Trait Loci in Sweet Cherry via Pedigree-Based and Genome-Wide Association Approaches.XLSX

Breeding for decreased fruit cracking incidence and increased fruit firmness in sweet cherry creates an attractive alternative to variable results from cultural management practices. DNA-informed breeding increases its efficiency, yet upstream research is needed to identify the genomic regions associated with the trait variation of a breeding-relevant magnitude, as well as to identify the parental sources of favorable alleles. The objectives of this research were to identify the quantitative trait loci (QTLs) associated with fruit cracking incidence and firmness, estimate the effects of single nucleotide polymorphism (SNP) haplotypes at the detected QTLs, and identify the ancestral source(s) of functional haplotypes. Fruit cracking incidence and firmness were evaluated for multiple years on 259 unselected seedlings representing 22 important breeding parents. Phenotypic data, in conjunction with genome-wide genotypic data from the RosBREED cherry 6K SNP array, were used in the QTL analysis performed via Pedigree-Based Analysis using the FlexQTL™ software, supplemented by a Genome-Wide Association Study using the BLINK software. Haplotype analysis was conducted on the QTLs to identify the functional SNP haplotypes and estimate their phenotypic effects, and the haplotypes were tracked through the pedigree. Four QTLs (two per trait) were consistent across the years and/or both analysis methods and validated the previously reported QTLs. qCrack-LG1.1m (the label given to a consistent QTL for cracking incidence on chromosome 1) explained 2–15.1% of the phenotypic variance, while qCrack-LG5.1m, qFirm-LG1.2m, and qFirm-LG3.2m explained 7.6–13.8, 8.8–21.8, and 1.7–10.1% of the phenotypic variance, respectively. At each QTL, at least two SNP haplotypes had significant effects and were considered putative functional SNP haplotypes. Putative low-cracking SNP haplotypes were tracked to an unnamed parent of ‘Emperor Francis’ and ‘Schmidt’ and unnamed parents of ‘Napoleon’ and ‘Hedelfingen,’ among others, and putative high-firmness haplotypes were tracked to an unnamed parent of ‘Emperor Francis’ and ‘Schmidt,’ an unnamed grandparent of ‘Black Republican,’ ‘Rube,’ and an unknown parent of ‘Napoleon.’ These four stable QTLs can now be targeted for DNA test development, with the goal of translating information discovered here into usable tools to aid in breeding decisions.</p

Image_2_Detection of Breeding-Relevant Fruit Cracking and Fruit Firmness Quantitative Trait Loci in Sweet Cherry via Pedigree-Based and Genome-Wide Association Approaches.JPEG

Breeding for decreased fruit cracking incidence and increased fruit firmness in sweet cherry creates an attractive alternative to variable results from cultural management practices. DNA-informed breeding increases its efficiency, yet upstream research is needed to identify the genomic regions associated with the trait variation of a breeding-relevant magnitude, as well as to identify the parental sources of favorable alleles. The objectives of this research were to identify the quantitative trait loci (QTLs) associated with fruit cracking incidence and firmness, estimate the effects of single nucleotide polymorphism (SNP) haplotypes at the detected QTLs, and identify the ancestral source(s) of functional haplotypes. Fruit cracking incidence and firmness were evaluated for multiple years on 259 unselected seedlings representing 22 important breeding parents. Phenotypic data, in conjunction with genome-wide genotypic data from the RosBREED cherry 6K SNP array, were used in the QTL analysis performed via Pedigree-Based Analysis using the FlexQTL™ software, supplemented by a Genome-Wide Association Study using the BLINK software. Haplotype analysis was conducted on the QTLs to identify the functional SNP haplotypes and estimate their phenotypic effects, and the haplotypes were tracked through the pedigree. Four QTLs (two per trait) were consistent across the years and/or both analysis methods and validated the previously reported QTLs. qCrack-LG1.1m (the label given to a consistent QTL for cracking incidence on chromosome 1) explained 2–15.1% of the phenotypic variance, while qCrack-LG5.1m, qFirm-LG1.2m, and qFirm-LG3.2m explained 7.6–13.8, 8.8–21.8, and 1.7–10.1% of the phenotypic variance, respectively. At each QTL, at least two SNP haplotypes had significant effects and were considered putative functional SNP haplotypes. Putative low-cracking SNP haplotypes were tracked to an unnamed parent of ‘Emperor Francis’ and ‘Schmidt’ and unnamed parents of ‘Napoleon’ and ‘Hedelfingen,’ among others, and putative high-firmness haplotypes were tracked to an unnamed parent of ‘Emperor Francis’ and ‘Schmidt,’ an unnamed grandparent of ‘Black Republican,’ ‘Rube,’ and an unknown parent of ‘Napoleon.’ These four stable QTLs can now be targeted for DNA test development, with the goal of translating information discovered here into usable tools to aid in breeding decisions.</p

Data_Sheet_2_Comparison of Single-Trait and Multi-Trait Genome-Wide Association Models and Inclusion of Correlated Traits in the Dissection of the Genetic Architecture of a Complex Trait in a Breeding Program.xlsx

Unknown genetic architecture makes it difficult to characterize the genetic basis of traits and associated molecular markers because of the complexity of small effect quantitative trait loci (QTLs), environmental effects, and difficulty in phenotyping. Seedling emergence of wheat (Triticum aestivum L.) from deep planting, has a poorly understood genetic architecture, is a vital factor affecting stand establishment and grain yield, and is historically correlated with coleoptile length. This study aimed to dissect the genetic architecture of seedling emergence while accounting for correlated traits using one multi-trait genome-wide association study (MT-GWAS) model and three single-trait GWAS (ST-GWAS) models. The ST-GWAS models included one single-locus model [mixed-linear model (MLM)] and two multi-locus models [fixed and random model circulating probability unification (FarmCPU) and Bayesian information and linkage-disequilibrium iteratively nested keyway (BLINK)]. We conducted GWAS using two populations. The first population consisted of 473 varieties from a diverse association mapping panel phenotyped from 2015 to 2019. The second population consisted of 279 breeding lines phenotyped in 2015 in Lind, WA, with 40,368 markers. We also compared the inclusion of coleoptile length and markers associated with reduced height as covariates in our ST-GWAS models. ST-GWAS found 107 significant markers across 19 chromosomes, while MT-GWAS found 82 significant markers across 14 chromosomes. The FarmCPU and BLINK models, including covariates, were able to identify many small effect markers while identifying large effect markers on chromosome 5A. By using multi-locus model breeding, programs can uncover the complex nature of traits to help identify candidate genes and the underlying architecture of a trait, such as seedling emergence.</p

Data_Sheet_3_Comparison of Single-Trait and Multi-Trait Genome-Wide Association Models and Inclusion of Correlated Traits in the Dissection of the Genetic Architecture of a Complex Trait in a Breeding Program.xlsx

Unknown genetic architecture makes it difficult to characterize the genetic basis of traits and associated molecular markers because of the complexity of small effect quantitative trait loci (QTLs), environmental effects, and difficulty in phenotyping. Seedling emergence of wheat (Triticum aestivum L.) from deep planting, has a poorly understood genetic architecture, is a vital factor affecting stand establishment and grain yield, and is historically correlated with coleoptile length. This study aimed to dissect the genetic architecture of seedling emergence while accounting for correlated traits using one multi-trait genome-wide association study (MT-GWAS) model and three single-trait GWAS (ST-GWAS) models. The ST-GWAS models included one single-locus model [mixed-linear model (MLM)] and two multi-locus models [fixed and random model circulating probability unification (FarmCPU) and Bayesian information and linkage-disequilibrium iteratively nested keyway (BLINK)]. We conducted GWAS using two populations. The first population consisted of 473 varieties from a diverse association mapping panel phenotyped from 2015 to 2019. The second population consisted of 279 breeding lines phenotyped in 2015 in Lind, WA, with 40,368 markers. We also compared the inclusion of coleoptile length and markers associated with reduced height as covariates in our ST-GWAS models. ST-GWAS found 107 significant markers across 19 chromosomes, while MT-GWAS found 82 significant markers across 14 chromosomes. The FarmCPU and BLINK models, including covariates, were able to identify many small effect markers while identifying large effect markers on chromosome 5A. By using multi-locus model breeding, programs can uncover the complex nature of traits to help identify candidate genes and the underlying architecture of a trait, such as seedling emergence.</p

Data_Sheet_1_Comparison of Single-Trait and Multi-Trait Genome-Wide Association Models and Inclusion of Correlated Traits in the Dissection of the Genetic Architecture of a Complex Trait in a Breeding Program.docx

Unknown genetic architecture makes it difficult to characterize the genetic basis of traits and associated molecular markers because of the complexity of small effect quantitative trait loci (QTLs), environmental effects, and difficulty in phenotyping. Seedling emergence of wheat (Triticum aestivum L.) from deep planting, has a poorly understood genetic architecture, is a vital factor affecting stand establishment and grain yield, and is historically correlated with coleoptile length. This study aimed to dissect the genetic architecture of seedling emergence while accounting for correlated traits using one multi-trait genome-wide association study (MT-GWAS) model and three single-trait GWAS (ST-GWAS) models. The ST-GWAS models included one single-locus model [mixed-linear model (MLM)] and two multi-locus models [fixed and random model circulating probability unification (FarmCPU) and Bayesian information and linkage-disequilibrium iteratively nested keyway (BLINK)]. We conducted GWAS using two populations. The first population consisted of 473 varieties from a diverse association mapping panel phenotyped from 2015 to 2019. The second population consisted of 279 breeding lines phenotyped in 2015 in Lind, WA, with 40,368 markers. We also compared the inclusion of coleoptile length and markers associated with reduced height as covariates in our ST-GWAS models. ST-GWAS found 107 significant markers across 19 chromosomes, while MT-GWAS found 82 significant markers across 14 chromosomes. The FarmCPU and BLINK models, including covariates, were able to identify many small effect markers while identifying large effect markers on chromosome 5A. By using multi-locus model breeding, programs can uncover the complex nature of traits to help identify candidate genes and the underlying architecture of a trait, such as seedling emergence.</p

Statistical power under different ranges of type 1 error and heritability.

<p><b>A</b>) Statistical power was examined on a trait simulated from 27 causative mutations (QTNs) sampled from SNPs on chromosomes 1 to 9 in maize data. The SNPs on the last chromosome (10) were used to derive null distribution. Power was defined as the proportion of detected QTNs under type I error of 5%. A total of 100 replications was conducted for each method. The heritability of the trait was set to 75%. Four methods were examined: 1) SUPER; 2) EMMAX; 3) FaST-LMM-Select; and 4) General linear model (GLM). <b>B</b>) Statistical power of four methods under different heritability levels. The four methods are SUPER, LMM-Selected, EMMAX and GLM.</p

SNPs found to be significant by SUPER and other three methods for AIL mouse data.

<p>P values that reached the Bonferroni correction threshold (1.6E-5) are shown in bold.</p><p>SNPs found to be significant by SUPER and other three methods for AIL mouse data.</p

Computing time and memory usage of five software packages.

<p>Three statistical models were performed by the five packages: 1) GLM by PLINK; 2) MLMs by EMMAX, GenABEL, and MLMM; and 3) FarmCPU by FarmCPU. Computing time <b>(a)</b> and memory usage <b>(b)</b> in response to sample size are displayed. The analyses were performed on a laptop (Asus A53S) running a Linux system (Ubuntu 12.10, 64 bit) with a 4.0 Gb of Random-Access Memory (RAM) and an Inter duo Core i3-2310M processor at 2.1 GHz. One core was used for this test. All datasets had 60,000 markers, but response was measured as a function of sample size. The last data point indicates the maximum sample size each software package could process without freezing the computer, except for PLINK and FarmCPU. The limitations for these two software packages were not reached with the maximum sample size examined.</p