Miki Okada

Complete Switchgrass Genetic Maps Reveal Subgenome Collinearity, Preferential Pairing and Multilocus Interactions

Polyploidy is an important aspect of the evolution of flowering plants. The potential of gene copies to diverge and evolve new functions is influenced by meiotic behavior of chromosomes leading to segregation as a single locus or duplicated loci. Switchgrass (Panicum virgatum) linkage maps were constructed using a full-sib population of 238 plants and SSR and STS markers to access the degree of preferential pairing and the structure of the tetraploid genome and as a step toward identification of loci underlying biomass feedstock quality and yield. The male and female framework map lengths were 1645 and 1376 cM with 97% of the genome estimated to be within 10 cM of a mapped marker in both maps. Each map coalesced into 18 linkage groups arranged into nine homeologous pairs. Comparative analysis of each homology group to the diploid sorghum genome identified clear syntenic relationships and collinear tracts. The number of markers with PCR amplicons that mapped across subgenomes was significantly fewer than expected, suggesting substantial subgenome divergence, while both the ratio of coupling to repulsion phase linkages and pattern of marker segregation indicated complete or near complete disomic inheritance. The proportion of transmission ratio distorted markers was relatively low, but the male map was more extensively affected by distorted transmission ratios and multilocus interactions, associated with spurious linkages.

De Novo Genome Assembly of the Economically Important Weed Horseweed Using Integrated Data from Multiple Sequencing Platforms

De novo genome assembly and genomic resources of horseweed will be useful to understand the genetic and molecular bases of weediness. Horseweed (Conyza canadensis), a member of the Compositae (Asteraceae) family, was the first broadleaf weed to evolve resistance to glyphosate. Horseweed, one of the most problematic weeds in the world, is a true diploid (2n = 2x = 18), with the smallest genome of any known agricultural weed (335 Mb). Thus, it is an appropriate candidate to help us understand the genetic and genomic bases of weediness. We undertook a draft de novo genome assembly of horseweed by combining data from multiple sequencing platforms (454 GS-FLX, Illumina HiSeq 2000, and PacBio RS) using various libraries with different insertion sizes (approximately 350 bp, 600 bp, 3 kb, and 10 kb) of a Tennessee-accessed, glyphosate-resistant horseweed biotype. From 116.3 Gb (approximately 350× coverage) of data, the genome was assembled into 13,966 scaffolds with 50% of the assembly = 33,561 bp. The assembly covered 92.3% of the genome, including the complete chloroplast genome (approximately 153 kb) and a nearly complete mitochondrial genome (approximately 450 kb in 120 scaffolds). The nuclear genome is composed of 44,592 protein-coding genes. Genome resequencing of seven additional horseweed biotypes was performed. These sequence data were assembled and used to analyze genome variation. Simple sequence repeat and single-nucleotide polymorphisms were surveyed. Genomic patterns were detected that associated with glyphosate-resistant or -susceptible biotypes. The draft genome will be useful to better understand weediness and the evolution of herbicide resistance and to devise new management strategies. The genome will also be useful as another reference genome in the Compositae. To our knowledge, this article represents the first published draft genome of an agricultural weed.

Evolution and spread of glyphosate resistance in Conyza canadensis in California.

Recent increases in glyphosate use in perennial crops of California, USA, are hypothesized to have led to an increase in selection and evolution of resistance to the herbicide in Conyza canadensis populations. To gain insight into the evolutionary origins and spread of resistance and to inform glyphosate resistance management strategies, we investigated the geographical distribution of glyphosate resistance in C. canadensis across and surrounding the Central Valley, its spatial relationship to groundwater protection areas (GWPA), and the genetic diversity and population structure and history using microsatellite markers. Frequencies of resistant individuals in 42 sampled populations were positively correlated with the size of GWPA within counties. Analyses of population genetic structure also supported spread of resistance in these areas. Bayesian clustering and approximate Bayesian computation (ABC) analyses revealed multiple independent origins of resistance within the Central Valley. Based on parameter estimation in the ABC analyses, resistant genotypes underwent expansion after glyphosate use began in agriculture, but many years before it was detected. Thus, diversity in weed control practices prior to herbicide regulation in GWPA probably kept resistance frequencies low. Regionally coordinated efforts to reduce seed dispersal and selection pressure are needed to manage glyphosate resistance in C. canadensis.

De novo genome assembly of the economically-important weed Conyza canadensis using integrated data from multiple sequencing platforms

De novo genome assembly and genomic resources of horseweed will be useful to understand the genetic and molecular bases of weediness. Horseweed (Conyza canadensis), a member of the Compositae (Asteraceae) family, was the first broadleaf weed to evolve resistance to glyphosate. Horseweed, one of the most problematic weeds in the world, is a true diploid (2n = 2x = 18), with the smallest genome of any known agricultural weed (335 Mb). Thus, it is an appropriate candidate to help us understand the genetic and genomic bases of weediness. We undertook a draft de novo genome assembly of horseweed by combining data from multiple sequencing platforms (454 GS-FLX, Illumina HiSeq 2000, and PacBio RS) using various libraries with different insertion sizes (approximately 350 bp, 600 bp, 3 kb, and 10 kb) of a Tennessee-accessed, glyphosate-resistant horseweed biotype. From 116.3 Gb (approximately 350× coverage) of data, the genome was assembled into 13,966 scaffolds with 50% of the assembly = 33,561 bp. The assembly covered 92.3% of the genome, including the complete chloroplast genome (approximately 153 kb) and a nearly complete mitochondrial genome (approximately 450 kb in 120 scaffolds). The nuclear genome is composed of 44,592 protein-coding genes. Genome resequencing of seven additional horseweed biotypes was performed. These sequence data were assembled and used to analyze genome variation. Simple sequence repeat and single-nucleotide polymorphisms were surveyed. Genomic patterns were detected that associated with glyphosate-resistant or -susceptible biotypes. The draft genome will be useful to better understand weediness and the evolution of herbicide resistance and to devise new management strategies. The genome will also be useful as another reference genome in the Compositae. To our knowledge, this article represents the first published draft genome of an agricultural weed.

Microsatellite variation points to local landscape plantings as sources of invasive pampas grass (Cortaderia selloana) in California

International trade in horticultural plants is a major pathway of introduction of invasive species. Pampas grass (Cortaderia selloana) is an invasive species of horticultural origin that is native to South America but cultivated as an ornamental in regions with Mediterranean climates worldwide. To gain insight into the introduction history of invasive populations in California, we analysed microsatellite marker variation in cultivated and invasive C. selloana. We sampled 275 cultivated plants from diverse sources and 698 invasive plants from 33 populations in four geographical regions of California. A model‐based Bayesian clustering analysis identified seven distinct gene pools in cultivated C. selloana. Probabilities of assignment of invasive individuals to cultivated gene pools indicated that two gene pools accounted for the genomic origin of 78% of the invasive C. selloana sampled. Extensive admixture between cultivated source gene pools was detected within invasive individuals. Sources of admixed invasive individuals are probably landscape plantings. Consistent with the Bayesian assignment results indicating that multiple cultivated gene pools and landscape plantings are probable sources of invasive populations, FST and neighbour‐joining clustering analyses indicated multiple escapes from shared sources in each geographical region. No isolation by distance or geographical trend in reduction of genetic diversity was evident. Furthermore, a generally random and discontinuous distribution of proportional assignments of invasive populations to cultivated gene pools suggests that introductions occurred recurrently within each geographical region. Our results strongly suggest that dispersal through local landscape plantings has contributed to the range expansion of invasive C. selloana in California.

Genetic Diversity and Spread of Thiobencarb Resistant Early Watergrass (Echinochloa oryzoides) in California

Resistance to the thiocarbamates has been selected in early watergrass populations within the rice-growing region of California. To elucidate the processes contributing to the spread of resistance among rice fields, we characterized the genetic diversity and differentiation of thiobencarb-resistant (R) and thiobencarb-susceptible (S) populations across the Central Valley using microsatellite markers. A total of 406 individuals from 22 populations were genotyped using seven nuclear microsatellite primer pairs. Three analytical approaches (unshared allele, Shannon–Weaver, and allelic-phenotype statistics) were used to assess genetic diversity and differentiation in the allohexaploid species. Low levels of genetic variation were detected within populations, consistent with other highly selfing species, with S populations tending to be more diverse than R populations. FST values indicated that populations were genetically differentiated and that genetic differentiation was greater among S populations than R populations. Principal coordinate analysis generated two orthogonal axes that explained 88% of the genetic variance among early watergrass populations and differentiated populations by geographical region, which was associated with resistance phenotype. A Mantel test revealed that genetic distances between R populations were positively correlated with the geographical distances separating populations. Taken together, our results suggest that both short- and long-distance seed dispersal, and multiple local and independent evolutionary events, are involved in the spread of thiobencarb-resistant early watergrass across rice fields in the Sacramento Valley. In contrast, resistance was not detected in early watergrass populations in the San Joaquin Valley. Nomenclature: Early watergrass; Echinochloa oryzoides (Ard.) Fritsch ECHOR

Inheritance of Glyphosate Resistance in Hairy Fleabane (Conyza bonariensis) from California

Abstract Inheritance of glyphosate resistance was investigated in hairy fleabane populations from California as part of providing the information needed to predict and manage resistance and to gain insight into resistance mechanism (or mechanisms) present in the populations. Three glyphosate-resistant individuals grown from seed collected from distinct sites near Fresno, CA, were crossed to individuals from the same susceptible population to create reciprocal F1 populations. A single individual from each of the F1 populations was used to create a backcross population with a susceptible maternal parent, and an F2 population. Based on dose response analyses, reciprocal F1 populations were not statistically different from each other, more similar to the resistant parent, and statistically different from the susceptible parent, consistent with nuclear control of the trait and dominance to incomplete dominance of resistance over susceptibility in all three crosses. Glyphosate resistance in two of the three crosses segregated in the backcross and the F2 populations as a single-locus trait. In the remaining cross, the resistant parent had approximately half the resistance level as the other two resistant parents, and the segregation of glyphosate resistance in backcross and F2 populations conformed to a two-locus model with resistance alleles acting additively and at least two copies of the allele required for expression of resistance. This two-locus model of the segregation of glyphosate resistance has not been reported previously. Variation in the pattern of inheritance and the level of resistance indicate that multiple resistance mechanisms may be present in hairy fleabane populations in California. Nomenclature: Glyphosate; hairy fleabane, Conyza bonariensis (L.) Cronq.

High-Density Single Nucleotide Polymorphism Linkage Maps of Lowland Switchgrass using Genotyping-by-Sequencing

Switchgrass (Panicum virgatum L.) is a warm‐season perennial grass with promising potential as a bioenergy crop in the United States. However, the lack of genomic resources has slowed the development of plant lines with optimal characteristics for sustainable feedstock production. We generated high‐density single nucleotide polymorphism (SNP) linkage maps using a reduced‐representation sequencing approach by genotyping 231 F1 progeny of a cross between two parents of lowland ecotype from the cultivars Kanlow and Alamo. Over 350 million reads were generated and aligned, which enabled identification and ordering of 4611 high‐quality SNPs. The total lengths of the resulting framework maps were 1770 cM for the Kanlow parent and 2059 cM for the Alamo parent. These maps show collinearity with maps generated with polymerase chain reaction (PCR)‐based simple‐sequence repeat (SSR) markers, and new SNP markers were identified in previously unpopulated regions of the genome. Transmission segregation distortion affected all linkage groups (LGs) to differing degrees, and ordering of distorted markers highlighted several regions of unequal inheritance. Framework maps were adversely affected by the addition of distorted markers with varying severity, but distorted maps were of higher marker density and provided additional information for analysis. Alignment of these linkage maps with a draft version of the switchgrass genome assembly demonstrated high levels of collinearity and provides greater confidence in the validity of both resources. This methodology has proven to be a rapid and cost‐effective way to generate high‐quality linkage maps of an outcrossing species.

Molecular evidence for hybridization in Colias (Lepidoptera: Pieridae): are Colias hybrids really hybrids?

Gene flow and hybridization among species dramatically affect our understanding of the species as a biological unit, species relationships, and species adaptations. In North American Colias eurytheme and Colias eriphyle, there has been historical debate over the extent of hybridization occurring and the identity of phenotypically intermediate individuals as genetic hybrids. This study assesses the population structure of these two species to measure the extent of hybridization and the genetic identity of phenotypic intermediates as hybrids. Amplified fragment length polymorphism (AFLP) marker analysis was performed on 378 specimens collected from northern California and Nevada. Population structure was inferred using a Bayesian/Markov chain Monte Carlo method, which probabilistically assigns individuals to genetic clusters. Three genetic clusters provided the best fit for the data. C. eurytheme individuals were primarily assigned to two closely related clusters, and C. eriphyle individuals were mostly assigned to a third, more distantly related cluster. There appeared to be significant hybridization between the two species. Individuals of intermediate phenotype (putative hybrids) were found to be genetically indistinguishable from C. eriphyle, indicating that previous work based on the assumption that these intermediate forms are hybrids may warrant reconsideration.