Genome‐Wide Association Study of Resistance to Cassava Green Mite Pest and Related Traits in Cassava

Cassava green mite [CGM, Mononychellus tanajoa (Bondar)] is a dry-season pest that usually feeds on the underside of young leaves causing leaf chlorosis, stunted growth, and root yield reduction by 80%. Since cassava (Manihot esculenta Crantz) leaves and roots serve as a primary staple food source, a decline in cassava yield can lead to household food, nutrition, and income insecurity. To evaluate the existence of CGM resistance alleles in the available germplasm, a diversity panel of 845 advanced breeding lines obtained from IITA, CIAT, and the National Root Crops Research Institute (NRCRI) were evaluated for CGM severity (CGMS), leaf pubescence (LP), leaf retention (LR), stay green, shoot tip compactness, and shoot tip size. A genome-wide association mapping detected 35 single-nucleotide polymorphisms (SNPs) markers significantly associated with CGMS, LP, and LR on chromosome 8. Colocalization of the most significant SNP associated with CGMS, LP, and LR on chromosome 8 is possibly an indication of pleiotropy or the presence of closely linked genes that regulate these traits. Seventeen candidate genes were found to be associated to CGM resistance. These candidate genes were subdivided into seven categories according to their protein structure namely, Zn finger, pentatricopeptide, MYB, MADS, homeodomain, trichome birefringence-related protein, and ethylene-responsive transcription factor genes. This study revealed significant loci associated with CGM, not previously reported, which together represent potential sources for the ongoing effort to develop multiple pestand disease-resistant cassava cultivars. L. Ezenwaka and C. Egesi, National Root Crops Research Institute (NRCRI), P.M.B, 7006 Umudike, Nigeria; L. Ezenwaka, E. Danquah, I. Asante, A. Danquah, and E. Blay, West Africa Center for Crop Improvement (WACCI), Univ. of Ghana, Legon, Accra, Ghana; D.P.D. Carpio, J.-L. Jannink, and C. Egesi, Dep. of Plant Breeding and Genetics, Cornell Univ., Ithaca, NY 14850; J.-L. Jannink, USDAARS, Ithaca, NY 14850; I. Rabbi and C. Egesi, IITA, Ibadan, Nigeria; E. Danquah, A. Danquah, and E. Blay, Dep. of Crop Science, Univ. of Ghana, Legon, Accra, Ghana; I. Asante, Dep. of Botany, Univ. of Ghana, Legon, Accra, Ghana. Received 11 Jan. 2018. Accepted 25 May 2018. *Corresponding author (cne22@cornell.edu). Assigned to Associate Editor Yiqun Weng. Abbreviations: BLUP, best linear unbiased prediction; CGM, cassava green mite; CGMS, cassava green mite severity; GBS, genotypingby-sequencing; GLM, general linear model; GWAS, genome-wide association study; H2, broad-sense heritability; LP, leaf pubescence; LR, leaf retention; MAF, minor allele frequency; MLM, mixed linear model; NRCRI, National Roots Crop Research Institute; PC, principal component; PCA, principal component analysis; QTL, quantitative trait locus; SG, stay green; SNP, single-nucleotide polymorphism; SSR, simple sequence repeat; STC, shoot tip compactness; STS, shoot tip size. Published in Crop Sci. 58:1907–1918 (2018). doi: 10.2135/cropsci2018.01.0024 © Crop Science Society of America | 5585 Guilford Rd., Madison, WI 53711 USA This is an open access article distributed under the CC BY license (https:// creativecommons.org/licenses/by/4.0/). Published online August 23, 2018

[1]  D. Connor,et al.  Cassava , 2019, Crop Physiology Case Histories for Major Crops.

[2]  D. Rokhsar,et al.  QTL Mapping for Pest and Disease Resistance in Cassava and Coincidence of Some QTL with Introgression Regions Derived from Manihot glaziovii , 2017, Front. Plant Sci..

[3]  D. P. Carpio,et al.  Prospects for Genomic Selection in Cassava Breeding , 2017, bioRxiv.

[4]  M. Gedil,et al.  Genome-Wide Association Mapping of Correlated Traits in Cassava: Dry Matter and Total Carotenoid Content , 2017, bioRxiv.

[5]  D. P. Carpio,et al.  Genome‐Wide Association and Prediction Reveals Genetic Architecture of Cassava Mosaic Disease Resistance and Prospects for Rapid Genetic Improvement , 2016, The plant genome.

[6]  M. Labuschagne,et al.  Stability and genotype by environment interaction of provitamin A carotenoid and dry matter content in cassava in Uganda , 2016, Breeding science.

[7]  J. Grimwood,et al.  Sequencing wild and cultivated cassava and related species reveals extensive interspecific hybridization and genetic diversity , 2016, Nature Biotechnology.

[8]  Guoqing Chen,et al.  Genome-wide association study of the husk number and weight in maize (Zea mays L.) , 2016, Euphytica.

[9]  R. Snowdon,et al.  Association mapping of seed quality traits in Brassica napus L. using GWAS and candidate QTL approaches , 2015, Molecular Breeding.

[10]  R. Cai,et al.  Transcriptome profiling of trichome-less reveals genes associated with multicellular trichome development in Cucumis sativus , 2015, Molecular Genetics and Genomics.

[11]  R. Melis,et al.  Inheritance of resistance to cassava green mite and other useful agronomic traits in cassava grown in Zambia , 2015, Euphytica.

[12]  I. Rabbi,et al.  The Effects of Restriction-Enzyme Choice on Properties of Genotyping-by-Sequencing Libraries: A Study in Cassava ( Manihot esculenta ) , 2014 .

[13]  D. Bates,et al.  Fitting Linear Mixed-Effects Models Using lme4 , 2014, 1406.5823.

[14]  Stephen D. Turner,et al.  qqman: an R package for visualizing GWAS results using Q-Q and manhattan plots , 2014, bioRxiv.

[15]  Robert J. Elshire,et al.  TASSEL-GBS: A High Capacity Genotyping by Sequencing Analysis Pipeline , 2014, PloS one.

[16]  M. Resende,et al.  Genome-wide selection in cassava , 2012, Euphytica.

[17]  Martin A. Fregene,et al.  Identificación de microsatélites en yuca asociados con la resistencia al ácaro Mononychellus tanajoa (Acari: Tetranychidae) , 2012, Revista Colombiana de Entomología.

[18]  Maoni Chao,et al.  Genome-wide association analysis detecting significant single nucleotide polymorphisms for chlorophyll and chlorophyll fluorescence parameters in soybean (Glycine max) landraces , 2012, Euphytica.

[19]  F. V. van Eeuwijk,et al.  Genome-wide association studies for agronomical traits in a world wide spring barley collection , 2012, BMC Plant Biology.

[20]  R. Hanna,et al.  Potential of the predatory mite, Amblyseius swirskii to suppress the Broad Mite, Polyphagotarsonemus latus on the Gboma Eggplant, Solanum macrocarpon , 2012, Journal of insect science.

[21]  H. Ceballos,et al.  Cassava Breeding: Current Status, Bottlenecks and the Potential of Biotechnology Tools , 2012, Tropical Plant Biology.

[22]  David M. Goodstein,et al.  Phytozome: a comparative platform for green plant genomics , 2011, Nucleic Acids Res..

[23]  Robert J. Elshire,et al.  A Robust, Simple Genotyping-by-Sequencing (GBS) Approach for High Diversity Species , 2011, PloS one.

[24]  J. Selbig,et al.  Involvement of TBL/DUF231 proteins into cell wall biology , 2010, Plant signaling & behavior.

[25]  B. Browning,et al.  A unified approach to genotype imputation and haplotype-phase inference for large data sets of trios and unrelated individuals. , 2009, American journal of human genetics.

[26]  M. Daly,et al.  Genetic Mapping in Human Disease , 2008, Science.

[27]  M. Calus,et al.  Accuracy of Genomic Selection Using Different Methods to Define Haplotypes , 2008, Genetics.

[28]  Edward S. Buckler,et al.  TASSEL: software for association mapping of complex traits in diverse samples , 2007, Bioinform..

[29]  F. Ariel,et al.  The true story of the HD-Zip family. , 2007, Trends in plant science.

[30]  P. Broun,et al.  Integration of cytokinin and gibberellin signalling by Arabidopsis transcription factors GIS, ZFP8 and GIS2 in the regulation of epidermal cell fate , 2007, Development.

[31]  M. McMullen,et al.  A unified mixed-model method for association mapping that accounts for multiple levels of relatedness , 2006, Nature Genetics.

[32]  G. Wasteneys,et al.  Transient Exposure to Ethylene Stimulates Cell Division and Alters the Fate and Polarity of Hypocotyl Epidermal Cells1 , 2004, Plant Physiology.

[33]  R. Hanna,et al.  Cassava green mite in Africa: a unique example of successful classical biological control of a mite pest on a continental scale. , 2003 .

[34]  D. Arnold,et al.  Heterologous myb genes distinct from GL1 enhance trichome production when overexpressed in Nicotiana tabacum. , 1999, Development.

[35]  A. Bellotti,et al.  The cassava mites , 1983 .

[36]  Cedric A. B. Smith,et al.  Introduction to Quantitative Genetics , 1960 .

[37]  R Core Team,et al.  R: A language and environment for statistical computing. , 2014 .

[38]  Ira M. Hall,et al.  BEDTools: a flexible suite of utilities for comparing genomic features , 2010, Bioinform..

[39]  A. Hallauer,et al.  Testers and Combining Ability , 2010 .

[40]  K. Leuschner Pest control for cassava and sweet potato. , 1982 .