mrMLM v4.0.2: An R Platform for Multi-locus Genome-wide Association Studies

Genome-Wide Association Study of Brown Rot (Monilinia spp.) Tolerance in Peach

Brown rot, caused by Monilinia spp., is one of the most important diseases on stone fruit worldwide. Severe yield loss can be caused by pre- and post-harvest fruit decay. Although some degree of tolerance has been reported in peach and almond, the genetic resistance in peach cultivars is still lacking. To date, only few genomic regions associated with brown rot response in fruit skin and flesh have been detected in peach. Previous studies suggested brown rot tolerance in peach being a polygenic quantitative trait. More information is needed to uncover the genetics behind brown rot tolerance in peach. To identify the genomic regions in peach associated with this trait, 26 cultivars and progeny from 9 crosses with 'Bolinha' sources of tolerance, were phenotyped across two seasons (2015 and 2016) for brown rot disease severity index in wounded and non-wounded fruits and genotyped using a newly developed 9+9K peach SNP array. Genome wide association study using single- and multi-locus methods by GAPIT version 3, mrMLM 4.0, GAPIT and G Model, revealed 14 reliable SNPs significantly associated with brown rot infection responses in peach skin (10) and flesh (4) across whole genome except for chromosome 3. Candidate gene analysis within the haplotype regions of the detected markers identified 25 predicted genes associated with pathogen infection response/resistance. Results presented here facilitate further understanding of genetics behind brown rot tolerance in peach and provide an important foundation for DNA-assisted breeding.

Uncovering Novel Genomic Regions and Candidate Genes for Senescence-Related Traits by Genome-Wide Association Studies in Upland Cotton (Gossypium hirsutum L.)

Senescence in plants is a complex trait, which is controlled by both genetic and environmental factors and can affect the yield and quality of cotton. However, the genetic basis of cotton senescence remains relatively unknown. In this study, we reported genome-wide association studies (GWAS) based on 185 accessions of upland cotton and 26,999 high-quality single-nucleotide polymorphisms (SNPs) to reveal the genetic basis of cotton senescence. To determine cotton senescence, we evaluated eight traits/indices. Our results revealed a high positive correlation (r>0.5) among SPAD value 20 days after topping (SPAD20d), relative difference of SPAD (RSPAD), nodes above white flower on topping day (NAWF0d), nodes above white flower 7 days after topping (NAWF7d), and number of open bolls on the upper four branches (NB), and genetic analysis revealed that all traits had medium or high heritability ranging from 0.53 to 0.86. Based on a multi-locus method (FASTmrMLM), a total of 63 stable and significant quantitative trait nucleotides (QTNs) were detected, which represented 50 genomic regions (GWAS risk loci) associated with cotton senescence. We observed three reliable loci located on chromosomes A02 (A02_105891088_107196428), D03 (D03_37952328_38393621) and D13 (D13_59408561_60730103) because of their high repeatability. One candidate gene (Ghir_D03G011060) was found in the locus D03_37952328_38393621, and its Arabidopsis thaliana homologous gene (AT5G23040) encodes a cell growth defect factor-like protein (CDF1), which might be involved in chlorophyll synthesis and cell death. Moreover, qRT-PCR showed that the transcript level of Ghir_D03G011060 was down-regulated in old cotton leaves, and virus-induced gene silencing (VIGS) indicated that silencing of Ghir_D03G011060 resulted in leaf chlorosis and promoted leaf senescence. In addition, two candidate genes (Ghir_A02G017660 and Ghir_D13G021720) were identified in loci A02_105891088_107196428 and D13_59408561_60730103, respectively. These results provide new insights into the genetic basis of cotton senescence and will serve as an important reference for the development and implementation of strategies to prevent premature senescence in cotton breeding programs.

Genome-Wide Association Study (GWAS) for Resistance to Sclerotinia sclerotiorum in Common Bean

White mold (WM) is a devastating fungal disease affecting common bean (Phaseolus vulgaris L.). In this research, a genome-wide association study (GWAS) for WM resistance was conducted using 294 lines of the Spanish diversity panel. One single-locus method and six multi-locus methods were used in the GWAS. Response to this fungus showed a continuous distribution, and 28 lines were identified as potential resistance sources, including lines of Andean and Mesoamerican origin, as well as intermediate lines between the two gene pools. Twenty-two significant associations were identified, which were organized into 15 quantitative trait intervals (QTIs) located on chromosomes Pv01, Pv02, Pv03, Pv04, Pv08, and Pv09. Seven of these QTIs were identified for the first time, whereas eight corresponded to chromosome regions previously identified in the WM resistance. In all, 468 genes were annotated in these regions, 61 of which were proposed potential candidate genes for WM resistance, based on their function related to the three main defense stages on the host: recognition (22), signal transduction (8), and defense response (31). Results obtained from this work will contribute to a better understanding of the complex quantitative resistance to WM in common bean and reveal information of significance for future breeding programs.