Genome-wide association study reveals novel SNPs and genes in Gossypium hirsutum underlying Aphis gossypii resistance

A. gossypii resistance showed great variability in G. hirsutum varieties. One hundred and seventy-six SNPs associated with A. gossypii resistance were identified using GWAS. Four candidate resistance genes were functionally validated. Aphis gossypii is an economically important sap-feeding pest and is widely distributed in the world's cotton-producing regions. Identification of cotton genotypes and developing cultivars with improved A. gossypii resistance (AGR) is essential and desirable for sustainable agriculture. In the present study, A. gossypii was offered no choice but to propagate on 200 Gossypium hirsutum accessions. A relative aphid reproduction index (RARI) was used to evaluate the AGR, which showed large variability in cotton accessions and was classified into 6 grades. A significantly positive correlation was found between AGR and Verticillium wilt resistance. A total of 176 SNPs significantly associated with the RARI were identified using GWAS. Of these, 21 SNPs could be repeatedly detected in three replicates. Cleaved amplified polymorphic sequence, a restriction digestion-based genotyping assay, was developed using SNP1 with the highest observed -log₁₀(P-value). Four genes within the 650 kb region of SNP1 were further identified, including GhRem (remorin-like), GhLAF1 (long after far-red light 1), GhCFIm25 (pre-mRNA cleavage factor Im 25 kDa subunit) and GhPMEI (plant invertase/pectin methylesterase inhibitor superfamily protein). The aphid infection could induce their expression and showed a significant difference between resistant and susceptible cotton varieties. Silencing of GhRem, GhLAF1 or GhCFIm25 could significantly increase aphid reproduction on cotton seedlings. Silencing of GhRem significantly reduced callose deposition, which is reasonably believed to be the cause for the higher AGR. Our results provide insights into understanding the genetic regulation of AGR in cotton and suggest candidate germplasms, SNPs and genes for developing cultivars with improved AGR.