Genetic structure of Cannabis sativa var. indica cultivars based on genomic SSR (gSSR) markers: Implications for breeding and germplasm management

Abstract Cannabis sativa L. is cultivated for its fiber or seeds (var. sativa ; hemp), or for its high content in cannabinoids (var. indica ; marijuana). Knowledge of the genetic structure of C. sativa var. indica is important for selection and breeding of cultivars with medicinal interest. We used six genomic SSRs (gSSRs) for genotyping 154 individual plants of 20 cultivars of C. sativa var. indica , plus two cultivars of C. sativa var. sativa . A very high polymorphism was observed, with an average of 17 alleles and 23.8 genotypes per locus. Expected ( H e ) and observed ( H o ) heterozygosities were high, with average values of 0.753 and 0.429, respectively. In some cultivars H e and H o presented similar values, while in others H e was considerably higher than H o suggesting that consanguinity and fixation had taken place during its development. In addition, some cultivars had a reduced number of alleles per locus (in some cases only two) indicating that a genetic bottleneck had taken place during its development. Gene flow ( Nm ) between both botanical varieties was high, with Nm  = 1.736. The molecular analysis of variance (AMOVA) revealed that only 31.94% of the molecular variation observed was caused by differences among cultivars, while the variation among plants of the same cultivar was of 37.11%, and within individual variation, due to heterozygosity, was of 30.96%. This indicates that a large variation exists within cultivars, which can be exploited for selection, but also complicates germplasm management and regeneration. The population structure analysis identified 14 genetic clusters, with most individuals of a single cultivar clustering together. This analysis, together with UPGMA cluster analysis shows that the two C. sativa var. sativa cultivars studied are differentiated from C. sativa var. indica , and that some cultivars of C. sativa var. indica seem to represent different selections from a common original cultivar. Our results represent the first comprehensive study of intra-varietal diversity in C. sativa var. indica and provide information of relevance for selection, breeding, and germplasm conservation, as well as for forensic studies in this crop.

[1]  Assessment of genotypic variation among cultivated durum wheat based on EST-SSRS and genomic SSRS , 2001 .

[2]  L. Cattivelli,et al.  The sexual differentiation of Cannabis sativa L.: A morphological and molecular study , 2004, Euphytica.

[3]  J. Ours.  The Long and Winding Road to Cannabis Legalization , 2012 .

[4]  J. Prohens,et al.  Use of Embryos Extracted from Individual Cannabis sativa Seeds for Genetic Studies and Forensic Applications , 2016, Journal of forensic sciences.

[5]  Rod Peakall,et al.  GenAlEx 6.5: genetic analysis in Excel. Population genetic software for teaching and research—an update , 2012, Bioinform..

[6]  G. Delogu,et al.  Differentiation between fiber and drug types of hemp (Cannabis sativa L.) from a collection of wild and domesticated accessions , 2013, Genetic Resources and Crop Evolution.

[7]  E S Buckler,et al.  Structure of linkage disequilibrium and phenotypic associations in the maize genome , 2001, Proceedings of the National Academy of Sciences of the United States of America.

[8]  Bridgett M. vonHoldt,et al.  STRUCTURE HARVESTER: a website and program for visualizing STRUCTURE output and implementing the Evanno method , 2011, Conservation Genetics Resources.

[9]  E. Small.  Evolution and Classification of Cannabis sativa (Marijuana, Hemp) in Relation to Human Utilization , 2015, The Botanical Review.

[10]  M. Plazas,et al.  Diversity, relationships, and genetic fingerprinting of the Listada de Gandía eggplant landrace using genomic SSRs and EST-SSRs , 2011 .

[11]  R. Peakall,et al.  Developmental Validation of a Cannabis sativa STR Multiplex System for Forensic Analysis , 2008, Journal of forensic sciences.

[12]  J. Doebley,et al.  Population structure and genetic diversity of New World maize races assessed by DNA microsatellites. , 2008, American journal of botany.

[13]  S. Wright.  THE INTERPRETATION OF POPULATION STRUCTURE BY F‐STATISTICS WITH SPECIAL REGARD TO SYSTEMS OF MATING , 1965 .

[14]  M. Divashuk,et al.  Molecular cytogenetic analysis of monoecious hemp (Cannabis sativa L.) cultivars reveals its karyotype variations and sex chromosomes constitution , 2015, Protoplasma.

[15]  H. M. Coyle,et al.  Field Testing of Collection Cards for Cannabis sativa Samples With a Single Hexanucleotide DNA Marker *,† , 2011, Journal of forensic sciences.

[16]  P. Smouse,et al.  genalex 6: genetic analysis in Excel. Population genetic software for teaching and research , 2006 .

[17]  David J Potter.  A review of the cultivation and processing of cannabis (Cannabis sativa L.) for production of prescription medicines in the UK. , 2014, Drug testing and analysis.

[18]  James Robertson,et al.  Short tandem repeat (STR) DNA markers are hypervariable and informative in Cannabis sativa: implications for forensic investigations. , 2003, Forensic science international.

[19]  Theunis Piersma,et al.  The interplay between habitat availability and population differentiation , 2012 .

[20]  M. Nei.  Analysis of gene diversity in subdivided populations. , 1973, Proceedings of the National Academy of Sciences of the United States of America.

[21]  Yoshio Tateno,et al.  Accuracy of estimated phylogenetic trees from molecular data , 1983, Journal of Molecular Evolution.

[22]  H. M. Coyle,et al.  Analysis of the NMI01 Marker for a Population Database of Cannabis Seeds , 2013, Journal of forensic sciences.

[23]  Atina G. Coté,et al.  The draft genome and transcriptome of Cannabis sativa , 2011, Genome Biology.

[24]  S. Myles,et al.  The Genetic Structure of Marijuana and Hemp , 2015, PLoS ONE.

[25]  Kejun Liu,et al.  PowerMarker: an integrated analysis environment for genetic marker analysis , 2005, Bioinform..

[26]  B. Ulrich,et al.  The structure of linkage , 1987 .

[27]  S. Harbison,et al.  Differentiation of drug and non-drug Cannabis using a single nucleotide polymorphism (SNP) assay. , 2011, Forensic science international.

[28]  Lei Liu,et al.  A Belated Green Revolution for Cannabis: Virtual Genetic Resources to Fast-Track Cultivar Development , 2016, Front. Plant Sci..

[29]  R. Highton,et al.  The relationship between the number of loci and the statistical support for the topology of UPGMA trees obtained from genetic distance data. , 1993, Molecular phylogenetics and evolution.

[30]  G. Evanno,et al.  Detecting the number of clusters of individuals using the software structure: a simulation study , 2005, Molecular ecology.

[31]  L. Keizer,et al.  Patterns of diversity in Cannabis , 1996, Genetic Resources and Crop Evolution.

[32]  C. Richards,et al.  Genetic diversity and population structure in Malus sieversii, a wild progenitor species of domesticated apple , 2009, Tree Genetics & Genomes.

[33]  P. Donnelly,et al.  Inference of population structure using multilocus genotype data. , 2000, Genetics.

[34]  I. Khan,et al.  Propagation through alginate encapsulation of axillary buds of Cannabis sativa L. — an important medicinal plant , 2009, Physiology and Molecular Biology of Plants.

[35]  G. Weiblen,et al.  Genetic Variation in Hemp and Marijuana (Cannabis sativa L.) According to Amplified Fragment Length Polymorphisms * , 2006, Journal of forensic sciences.

[36]  M. Rai,et al.  Microsatellite markers: an overview of the recent progress in plants , 2011, Euphytica.

[37]  S WRIGHT,et al.  Genetical Structure of Populations , 1950, British medical journal.

[38]  Jianbin Hu,et al.  Microsatellite Diversity, Population Structure, and Core Collection Formation in Melon Germplasm , 2015, Plant Molecular Biology Reporter.

[39]  Stefano Amaducci,et al.  New developments in fiber hemp (Cannabis sativa L.) breeding , 2015 .

[40]  Ping Chen,et al.  Diversity Analysis in Cannabis sativa Based on Large-Scale Development of Expressed Sequence Tag-Derived Simple Sequence Repeat Markers , 2014, PloS one.

[41]  D. Botstein,et al.  Construction of a genetic linkage map in man using restriction fragment length polymorphisms. , 1980, American journal of human genetics.