Genomics of Compositae weeds: EST libraries, microarrays, and evidence of introgression

Transcriptomic analysis of Amaranthus retroflex resistant to PPO-inhibitory herbicides

Amaranthus retroflexus L. is one of the malignant weeds which can cause a reduction in the soybean yield. We found a population of A. retroflexus (R-Q) resistant to fomesafen through the initial screening of whole-plant dose response bioassay in the research. The resistance index of the population (R-Q) was 183 times of the sensitive population (S-N). The resistant and sensitive populations were used as experimental materials in the paper. Strand-specific RNA-Seq analyses of R‒Q and S‒N populations obtained from herbicide-treated and mock-treated leaf samples after treatment were conducted to generate a full-length transcriptome database. We analyzed differentially expressed genes (DEGs) among the R-Q and S‒N A. retroflexus populations treated with recommended dose and mock-treated on the 1st (24 h) and 3rd (72 h) days to identify genes involved in fomesafen resistance. All 82,287 unigenes were annotated by Blastx search with E-value < 0.00001 from 7 databases. A total of 94,815 DEGs among the three group comparisons were identified. Two nuclear genes encoding PPO (PPX1 and PPX2) and five unigenes belonging to the AP2-EREBP, GRAS, NAC, bHLH and bZIP families exhibited different expression patterns between individuals of S‒N and R-Q populations. The A. retroflexus transcriptome and specific transcription factor families which can respond to fomesafen in resistant and susceptible genotypes were reported in this paper. The PPX1 and PPX2 genes of the target enzyme were identified. The study establishes the foundation for future research and provides opportunities to manage resistant weeds better.

Wild sunflower goes viral: Citizen science and comparative genomics allow tracking the origin and establishment of invasive sunflower in the Levant

Globalization and intensified volume of trade and transport around the world are accelerating the rate of biological invasions. It is therefore increasingly important to understand the processes through which invasive species colonize new habitats, often to the detriment of native flora. The initial steps of an invasion are particularly critical, as the introduced species relies on limited genetic diversity to adapt to a new environment. However, our understanding of this critical stage of the invasion is currently limited. We used a citizen science approach and social media to survey the distribution of invasive sunflower in Israel. We then sampled and sequenced a representative collection and compared it with available genomic data sets of North American wild sunflower, landraces and cultivars. We show that invasive wild sunflower is rapidly establishing throughout Israel, probably from a single, recent introduction from Texas, while maintaining high genetic diversity through ongoing gene flow. Since its introduction, invasive sunflower has spread quickly to most regions, and differentiation was detected despite extensive gene flow between clusters. Our findings suggest that rapid spread followed by continuous gene flow between diverging populations can serve as an efficient mechanism for maintaining sufficient genetic diversity at the early stages of invasion, promoting rapid adaptation and establishment in the new territory.

Selective ancestral sorting and de novo evolution in the agricultural invasion of Amaranthus tuberculatus

The relative role of hybridization, de novo evolution, and standing variation in weed adaptation to agricultural environments is largely unknown. In Amaranthus tuberculatus, a widespread North American agricultural weed, adaptation is likely influenced by recent secondary contact and admixture of two previously isolated lineages. We characterized the extent of adaptation and phenotypic differentiation accompanying the spread of A. tuberculatus into agricultural environments and the contribution of ancestral divergence. We generated phenotypic and whole-genome sequence data from a manipulative common garden experiment, using paired samples from natural and agricultural populations. We found strong latitudinal, longitudinal, and sex differentiation in phenotypes, and subtle differences among agricultural and natural environments that were further resolved with ancestry inference. The transition into agricultural environments has favored southwestern var. rudis ancestry that leads to higher biomass and treatment-specific phenotypes: increased biomass and earlier flowering under reduced water availability, and reduced plasticity in fitness-related traits. We also detected de novo adaptation in individuals from agricultural habitats independent of ancestry effects, including marginally higher biomass, later flowering, and treatment-dependent divergence in time to germination. Therefore, the invasion of A. tuberculatus into agricultural environments has drawn on adaptive variation across multiple timescales-through both preadaptation via the preferential sorting of var. rudis ancestry and de novo local adaptation.

Pilot RNA-seq data from 24 species of vascular plants at Harvard Forest

PREMISE

Large-scale projects such as the National Ecological Observatory Network (NEON) collect ecological data on entire biomes to track climate change. NEON provides an opportunity to launch community transcriptomic projects that ask integrative questions in ecology and evolution. We conducted a pilot study to investigate the challenges of collecting RNA-seq data from diverse plant communities.

METHODS

We generated >650 Gbp of RNA-seq for 24 vascular plant species representing 12 genera and nine families at the Harvard Forest NEON site. Each species was sampled twice in 2016 (July and August). We assessed transcriptome quality and content with TransRate, BUSCO, and Gene Ontology annotations.

RESULTS

Only modest differences in assembly quality were observed across multiple k-mers. On average, transcriptomes contained hits to >70% of loci in the BUSCO database. We found no significant difference in the number of assembled and annotated transcripts between diploid and polyploid transcriptomes.

DISCUSSION

We provide new RNA-seq data sets for 24 species of vascular plants in Harvard Forest. Challenges associated with this type of study included recovery of high-quality RNA from diverse species and access to NEON sites for genomic sampling. Overcoming these challenges offers opportunities for large-scale studies at the intersection of ecology and genomics.

Genomic Analyses of Phenotypic Differences Between Native and Invasive Populations of Diffuse Knapweed (Centaurea diffusa)

Invasive species represent excellent opportunities to study the evolutionary potential of traits important to success in novel environments. Although some ecologically important traits have been identified in invasive species, little is typically known about the genetic mechanisms that underlie invasion success in non-model species. Here, we use a genome-wide association (GWAS) approach to identify the genetic basis of trait variation in the non-model, invasive, diffuse knapweed [Centaurea diffusa Lam. (Asteraceae)]. To assist with this analysis, we have assembled the first draft genome reference and fully annotated plastome assembly for this species, and one of the first from this large, weedy, genus, which is of major ecological and economic importance. We collected phenotype data from 372 individuals from four native and four invasive populations of C. diffusa grown in a common environment. Using these individuals, we produced reduced-representation genotype-by-sequencing (GBS) libraries and identified 7,058 SNPs. We identify two SNPs associated with leaf width in these populations, a trait which significantly varies between native and invasive populations. In this rosette forming species, increased leaf width is a major component of increased biomass, a common trait in invasive plants correlated with increased fitness. Finally, we use annotations from Arabidopsis thaliana to identify 98 candidate genes that are near the associated SNPs and highlight several good candidates for leaf width variation.

The Role of Hybridisation in the Making of the Species-Rich Arctic-Alpine Genus Saxifraga (Saxifragaceae)

Evolutionary processes fuelling rapid species diversification are not yet fully understood, although their major contribution to overall patterns of plant biodiversity is well established. Hybridisation is among the least understood of these processes, despite its multifaceted role in speciation processes being widely accepted. Species of the large arctic-alpine genus Saxifraga are notorious for their ability to hybridise; however, the overall role of hybridisation and polyploidisation for the diversification of this genus remains unknown. Here, we provide a comprehensive genus-wide review of hybridisation accounts and ploidy levels. We find that the sections of Saxifraga vary greatly in their propensity to hybridise. The majority of natural hybridisation accounts are from recent localised events (n = 71). Hybridisation hotspots were located in the Pyrenees and the European Alps, thus contrasting with the overall distribution of species richness in the genus. Hybrids or hybrid populations are often short-lived in Saxifraga due to a multitude of reproductive barriers, most commonly low F1 hybrid fertility. However, these barriers are not always fully effective, allowing for backcrossing and the formation of hybrid swarms. In addition, we find that the incidence of polyploidy varies widely across different sections of Saxifraga, with species-rich sections Porphyrion and Saxifraga showing divergent polyploidy proportions. Overall, we show that hybridisation and polyploidisation played differential roles in the diversification of this large genus. Nevertheless, a significant proportion of species are yet to be scrutinised, particularly among the Asian Saxifraga species, illustrating the need for systematic further study to fully unravel the role of hybridisation during the evolution of Saxifraga.

Integrated Analysis of Transcriptome and Metabolome in Cirsium japonicum Fisch ex DC

Cirsium japonicum Fisch ex DC belongs to the Compositae family and has been used as a folk remedy source in Asian countries because of its health-promoting properties. It is known that C. japonicum contains flavonoids, furans, long-chain alcohols, sterols, and volatile oils. Nevertheless, the molecular mechanism of secondary metabolite biosynthesis remains poorly understood. Therefore, transcriptome analysis and metabolic profiling were performed using different parts of C. japonicum to investigate phenylpropanoid metabolism. Based on the BLASTX search results, we identified 29 orthologs of enzymes responsible for phenylpropanoid biosynthesis. Additionally, 75 metabolites were identified in C. japonicum. Most of the flavonoid biosynthetic genes were significantly expressed ranging from 2.6- to 500-fold higher in the flowers than those in the leaves. Correspondently, the total content of flavonols was 21-fold higher in the flowers than in the roots. However, the total level of flavones showed 58-fold higher amounts in the leaves than in the flowers. Additionally, the total content of flavanols was 19-fold higher in the leaves than in the roots. The results of this study provide transcriptomic and metabolic information to elucidate the tissue-specific phenylpropanoid metabolism of C. japonicum.

De Novo Plant Transcriptome Assembly and Annotation Using Illumina RNA-Seq Reads

The ability to identify and quantify transcribed sequences from a multitude of organisms using high-throughput RNA sequencing has revolutionized our understanding of genetics and plant biology. However, a number of computational tools used in these analyses still require a reference genome sequence, something that is seldom available for non-model organisms. Computational tools employing de Bruijn graphs to reconstruct full-length transcripts from short sequence reads allow for de novo transcriptome assembly. Here we provide detailed methods for generating and annotating de novo transcriptome assembly from plant RNA-seq data.

North European invasion by common ragweed is associated with early flowering and dominant changes in FT/TFL1 expression

During the last two centuries, the North American common ragweed (Ambrosia artemisiifolia L.) invaded a large part of the globe. Local adaptation of this species was revealed by a common garden experiment, demonstrating that the distribution of the species in Europe could extend considerably to the North. Our study compares two populations of common ragweed (one from the native range and one from the invaded range) that differ in flowering time in the wild: the invasive population flowers earlier than the native population under non-inductive long-day photoperiods. Experiments conducted in controlled environments established that the two populations differ in their flowering time even under inductive short-day photoperiods, suggesting a change in autonomous flowering control. Genetic analysis revealed that early flowering is dominantly inherited and accompanied by the increased expression of the floral activator AaFTL1 and decreased expression of the floral repressor AaFTL2. Early flowering is also accompanied by reduced reproductive output, which is evolutionarily disadvantageous under long vegetation periods. In contrast, under short vegetation periods, only early-flowering plants can produce any viable seeds, making the higher seed set of late-flowering plants irrelevant. Thus, earlier flowering appears to be a specific adaptation to the higher latitudes of northern Europe.

Improved transcriptome sampling pinpoints 26 ancient and more recent polyploidy events in Caryophyllales, including two allopolyploidy events

Studies of the macroevolutionary legacy of polyploidy are limited by an incomplete sampling of these events across the tree of life. To better locate and understand these events, we need comprehensive taxonomic sampling as well as homology inference methods that accurately reconstruct the frequency and location of gene duplications. We assembled a data set of transcriptomes and genomes from 168 species in Caryophyllales, of which 43 transcriptomes were newly generated for this study, representing one of the most densely sampled genomic-scale data sets available. We carried out phylogenomic analyses using a modified phylome strategy to reconstruct the species tree. We mapped the phylogenetic distribution of polyploidy events by both tree-based and distance-based methods, and explicitly tested scenarios for allopolyploidy. We identified 26 ancient and more recent polyploidy events distributed throughout Caryophyllales. Two of these events were inferred to be allopolyploidy. Through dense phylogenomic sampling, we show the propensity of polyploidy throughout the evolutionary history of Caryophyllales. We also provide a framework for utilizing transcriptome data to detect allopolyploidy, which is important as it may have different macroevolutionary implications compared with autopolyploidy.

Gene flow in Argentinian sunflowers as revealed by genotyping-by-sequencing data

Gene flow can have several different applied consequences, ranging from extinction to the escape of transgenes to the evolution of weedy or invasive lineages. Here, we describe patterns of hybridization and gene flow involving domesticated and wild sunflowers in Argentina. To address the risks of introgression of variants from the cultivated sunflower into invasive wild Helianthus, we used genotyping‐by‐sequencing (GBS) to genotype 182 samples from 11 sites in Argentina, along with previously published data from samples from the native range (North America), to determine the native source populations of the Argentinian samples and to detect admixture. We unexpectedly discovered two distinctive forms of H. petiolaris in Argentina, one from H. petiolaris subsp. petiolaris as expected, but the other from an unknown source. Extensive admixture was observed among Argentinian sunflowers, largely confirming phenotypic predictions. While many hybrids are F1s, there were signals consistent with introgression from the domesticated sunflower into H. petiolaris. Whether this introgression is incidental or a causal driver of invasiveness is not yet clear, but it seems likely that genes found in the domesticated sunflower genome (whether engineered or not) will quickly find their way into wild Argentinian sunflower populations.

De novo transcriptome assembly reveals high transcriptional complexity in Pisum sativum axillary buds and shows rapid changes in expression of diurnally regulated genes

BACKGROUND

The decision for a bud to grow into a branch is a key regulatory process affecting plant architecture. In order to study molecular processes regulating axillary bud outgrowth in the model plant garden pea (Pisum sativum), we sequenced the axillary bud transcriptome and performed de novo transcriptome assembly.

RESULTS

We assembled a pea axillary bud transcriptome into 81,774 transcripts comprised of 194,067 isoforms. This new pea transcriptome resource is both comprehensive and representative, as shown by comparison to other available pea sequence resources. Over half of the transcriptome could be annotated based on sequence homology to Arabidopsis thaliana proteins, while almost one quarter of the isoforms were identified as putative long non-coding RNAs (lncRNAs). This transcriptome will be useful in studies of pea buds because it includes genes expressed specifically in buds which are not represented in other transcriptome studies. We also investigated the impact of a short time collection series on gene expression. Differential gene expression analysis identified 142 transcripts changing within the short 170 min time frame that the buds were harvested within. Thirty-three of these transcripts are implicated in diurnal fluctuations in other flowering plants, while the remaining transcripts include 31 putative lncRNA. Further investigation of the differentially expressed transcripts found an enrichment of genes involved in post-transcriptional regulation, including RNA processing and modification, as well as genes involved in fatty acid biosynthesis and oxidative phosphorylation.

CONCLUSIONS

We have sequenced and assembled a high quality pea bud transcriptome containing both coding and non-coding RNA transcripts that will be useful for further studies into axillary bud outgrowth. Over the short sample collection time frame of just 170 min, we identified differentially expressed coding and non-coding RNA, some of which are implicated in diurnal regulation, highlighting the utility of our transcriptome resource in identifying gene expression changes and informing future experimental designs.

De novo transcriptome assembly analysis of weed Apera spica-venti from seven tissues and growth stages

BACKGROUND

Loose silky bentgrass (Apera spica-venti) is an important weed in Europe with a recent increase in herbicide resistance cases. The lack of genetic information about this noxious weed limits its biological understanding such as growth, reproduction, genetic variation, molecular ecology and metabolic herbicide resistance. This study produced a reference transcriptome for A. spica-venti from different tissues (leaf, root, stem) and various growth stages (seed at phenological stages 05, 07, 08, 09). The de novo assembly was performed on individual and combined dataset followed by functional annotations. Individual transcripts and gene families involved in metabolic based herbicide resistance were identified.

RESULTS

Eight separate transcriptome assemblies were performed and compared. The combined transcriptome assembly consists of 83,349 contigs with an N50 and average contig length of 762 and 658 bp, respectively. This dataset contains 74,724 transcripts consisting of total 54,846,111 bp. Among them 94% had a homologue to UniProtKB, 73% retrieved a GO mapping, and 50% were functionally annotated. Compared with other grass species, A. spica-venti has 26% proteins in common to Brachypodium distachyon, and 41% to Lolium spp. Glycosyltransferases had the highest number of transcripts in each tissue followed by the cytochrome P450s. The GSTF1 and CYP89A2 transcripts were recovered from the majority of tissues and aligned at a maximum of 66 and 30% to proven herbicide resistant allele from Alopecurus myosuroides and Lolium rigidum, respectively.

CONCLUSIONS

De novo transcriptome assembly enabled the generation of the first reference transcriptome of A. spica-venti. This can serve as stepping stone for understanding the metabolic herbicide resistance as well as the general biology of this problematic weed. Furthermore, this large-scale sequence data is a valuable scientific resource for comparative transcriptome analysis for Poaceae grasses.

Chloroplast genome resources and molecular markers differentiate rubber dandelion species from weedy relatives

BACKGROUND

Rubber dandelion (Taraxacum kok-saghyz, TK) is being developed as a domestic source of natural rubber to meet increasing global demand. However, the domestication of TK is complicated by its colocation with two weedy dandelion species, Taraxacum brevicorniculatum (TB) and the common dandelion (Taraxacum officinale, TO). TB is often present as a seed contaminant within TK accessions, while TO is a pandemic weed, which may have the potential to hybridize with TK. To discriminate these species at the molecular level, and facilitate gene flow studies between the potential rubber crop, TK, and its weedy relatives, we generated genomic and marker resources for these three dandelion species.

RESULTS

Complete chloroplast genome sequences of TK (151,338 bp), TO (151,299 bp), and TB (151,282 bp) were obtained using the Illumina GAII and MiSeq platforms. Chloroplast sequences were analyzed and annotated for all the three species. Phylogenetic analysis within Asteraceae showed that TK has a closer genetic distance to TB than to TO and Taraxacum species were most closely related to lettuce (Lactuca sativa). By sequencing multiple genotypes for each species and testing variants using gel-based methods, four chloroplast Single Nucleotide Polymorphism (SNP) variants were found to be fixed between TK and TO in large populations, and between TB and TO. Additionally, Expressed Sequence Tag (EST) resources developed for TO and TK permitted the identification of five nuclear species-specific SNP markers.

CONCLUSIONS

The availability of chloroplast genomes of these three dandelion species, as well as chloroplast and nuclear molecular markers, will provide a powerful genetic resource for germplasm differentiation and purification, and the study of potential gene flow among Taraxacum species.

Identification and characterization of contrasting sunflower genotypes to early leaf senescence process combining molecular and physiological studies (Helianthus annuus L.)

Leaf senescence is a complex mechanism ruled by multiple genetic and environmental variables that affect crop yields. It is the last stage in leaf development, is characterized by an active decline in photosynthetic rate, nutrients recycling and cell death. The aim of this work was to identify contrasting sunflower inbred lines differing in leaf senescence and to deepen the study of this process in sunflower. Ten sunflower genotypes, previously selected by physiological analysis from 150 inbred genotypes, were evaluated under field conditions through physiological, cytological and molecular analysis. The physiological measurement allowed the identification of two contrasting senescence inbred lines, R453 and B481-6, with an increase in yield in the senescence delayed genotype. These findings were confirmed by cytological and molecular analysis using TUNEL, genomic DNA gel electrophoresis, flow sorting and gene expression analysis by qPCR. These results allowed the selection of the two most promising contrasting genotypes, which enables future studies and the identification of new biomarkers associated to early senescence in sunflower. In addition, they allowed the tuning of cytological techniques for a non-model species and its integration with molecular variables.

Most Compositae (Asteraceae) are descendants of a paleohexaploid and all share a paleotetraploid ancestor with the Calyceraceae

PREMISE OF THE STUDY

Like many other flowering plants, members of the Compositae (Asteraceae) have a polyploid ancestry. Previous analyses found evidence for an ancient duplication or possibly triplication in the early evolutionary history of the family. We sought to better place this paleopolyploidy in the phylogeny and assess its nature.

METHODS

We sequenced new transcriptomes for Barnadesia, the lineage sister to all other Compositae, and four representatives of closely related families. Using a recently developed algorithm, MAPS, we analyzed nuclear gene family phylogenies for evidence of paleopolyploidy.

KEY RESULTS

We found that the previously recognized Compositae paleopolyploidy is also in the ancestry of the Calyceraceae. Our phylogenomic analyses uncovered evidence for a successive second round of genome duplication among all sampled Compositae except Barnadesia.

CONCLUSIONS

Our analyses of new samples with new tools provide a revised view of paleopolyploidy in the Compositae. Together with results from a high density Lactuca linkage map, our results suggest that the Compositae and Calyceraceae have a common paleotetraploid ancestor and that most Compositae are descendants of a paleohexaploid. Although paleohexaploids have been previously identified, this is the first example where the paleotetraploid and paleohexaploid lineages have survived over tens of millions of years. The complex polyploidy in the ancestry of the Compositae and Calyceraceae represents a unique opportunity to study the long-term evolutionary fates and consequences of different ploidal levels.

Is hybridization driving the evolution of climatic niche in Alyssum montanum

PREMISE OF THE STUDY

After decades of interest, the contribution of hybridization to ecological diversification remains unclear. Hybridization is a potent source of novelty, but nascent hybrid lineages must overcome reproductive and ecological competition from their parental species. Here, we assess whether hybrid speciation is advantageous over alternative modes of speciation, by comparing the geographical and ecological ranges and climatic niche evolutionary rates of stabilized allopolyploid vs. autopolyploids in the Alyssum montanum species complex.

METHODS

We combined an extensive review of studies addressing the systematics and genetic diversity of A. montanum s.l., with flow cytometry and cloning of nuclear markers, to establish the ploidy level and putative hybrid nature of 205 populations. The respective geographic distribution and climatic niche evolution dynamics of the allo- and autopolyploids were investigated using multivariate analyses and comparative phylogenetic approaches.

KEY RESULTS

As expected by theory, allopolyploids occur mainly along contact zones and are generally spatially overlapping with their diploid counterparts. However, they demonstrate higher rates of niche evolution and expand into different climatic conditions than those of their diploid congeners. In contrast, autopolyploids show lower rates of niche evolution, occupy ecological niches similar to their ancestors and are restricted to less competitive and peripheral geographic areas.

CONCLUSIONS

Hybridization thus seems advantageous by promoting ecological niche evolution and more readily allowing escape from competitive exclusion.

The population genomic basis of geographic differentiation in North American common ragweed (Ambrosia artemisiifolia L.)

Common ragweed (Ambrosia artemisiifolia L.) is an invasive, wind-pollinated plant nearly ubiquitous in disturbed sites in its eastern North American native range and present across growing portions of Europe, Africa, Asia, and Australia. Phenotypic divergence between European and native-range populations has been described as rapid evolution. However, a recent study demonstrated major human-mediated shifts in ragweed genetic structure before introduction to Europe and suggested that native-range genetic structure and local adaptation might fully explain accelerated growth and other invasive characteristics of introduced populations. Genomic differentiation that potentially influenced this structure has not yet been investigated, and it remains unclear whether substantial admixture during historical disturbance of the native range contributed to the development of invasiveness in introduced European ragweed populations. To investigate fine-scale population genetic structure across the species' native range, we characterized diallelic SNP loci via a reduced-representation genotyping-by-sequencing (GBS) approach. We corroborate phylogeographic domains previously discovered using traditional sequencing methods, while demonstrating increased power to resolve weak genetic structure in this highly admixed plant species. By identifying exome polymorphisms underlying genetic differentiation, we suggest that geographic differentiation of this important invasive species has occurred more often within pathways that regulate growth and response to defense and stress, which may be associated with survival in North America's diverse climatic regions.

Genomics of crop wild relatives: expanding the gene pool for crop improvement

Plant breeders require access to new genetic diversity to satisfy the demands of a growing human population for more food that can be produced in a variable or changing climate and to deliver the high-quality food with nutritional and health benefits demanded by consumers. The close relatives of domesticated plants, crop wild relatives (CWRs), represent a practical gene pool for use by plant breeders. Genomics of CWR generates data that support the use of CWR to expand the genetic diversity of crop plants. Advances in DNA sequencing technology are enabling the efficient sequencing of CWR and their increased use in crop improvement. As the sequencing of genomes of major crop species is completed, attention has shifted to analysis of the wider gene pool of major crops including CWR. A combination of de novo sequencing and resequencing is required to efficiently explore useful genetic variation in CWR. Analysis of the nuclear genome, transcriptome and maternal (chloroplast and mitochondrial) genome of CWR is facilitating their use in crop improvement. Genome analysis results in discovery of useful alleles in CWR and identification of regions of the genome in which diversity has been lost in domestication bottlenecks. Targeting of high priority CWR for sequencing will maximize the contribution of genome sequencing of CWR. Coordination of global efforts to apply genomics has the potential to accelerate access to and conservation of the biodiversity essential to the sustainability of agriculture and food production.

Selecting Superior De Novo Transcriptome Assemblies: Lessons Learned by Leveraging the Best Plant Genome

Whereas de novo assemblies of RNA-Seq data are being published for a growing number of species across the tree of life, there are currently no broadly accepted methods for evaluating such assemblies. Here we present a detailed comparison of 99 transcriptome assemblies, generated with 6 de novo assemblers including CLC, Trinity, SOAP, Oases, ABySS and NextGENe. Controlled analyses of de novo assemblies for Arabidopsis thaliana and Oryza sativa transcriptomes provide new insights into the strengths and limitations of transcriptome assembly strategies. We find that the leading assemblers generate reassuringly accurate assemblies for the majority of transcripts. At the same time, we find a propensity for assemblers to fail to fully assemble highly expressed genes. Surprisingly, the instance of true chimeric assemblies is very low for all assemblers. Normalized libraries are reduced in highly abundant transcripts, but they also lack 1000s of low abundance transcripts. We conclude that the quality of de novo transcriptome assemblies is best assessed through consideration of a combination of metrics: 1) proportion of reads mapping to an assembly 2) recovery of conserved, widely expressed genes, 3) N50 length statistics, and 4) the total number of unigenes. We provide benchmark Illumina transcriptome data and introduce SCERNA, a broadly applicable modular protocol for de novo assembly improvement. Finally, our de novo assembly of the Arabidopsis leaf transcriptome revealed ~20 putative Arabidopsis genes lacking in the current annotation.

Utilizing next-generation sequencing to study homeologous polymorphisms and herbicide-resistance-endowing mutations in Poa annua acetolactate synthase genes

BACKGROUND

Detection of single nucleotide polymorphisms (SNPs) related to herbicide resistance in non-model polyploid weed species is fraught with difficulty owing to the gene duplication and lack of reference sequences. Our research seeks to overcome these obstacles by Illumina HiSeq read mapping, SNP calling and allele frequency determinations. Our focus is on the acetolactate synthase (ALS) gene, the target site of ALS-inhibiting herbicides, in Poa annua, an allotetraploid weed species originating from two diploid parents, P. supina and P. infirma.

RESULTS

ALS contigs with complete coding regions of P. supina, P. infirma and P. annua were assembled and compared with ALS genes from other plant species. The ALS infirma-homeolog of P. annua showed higher levels of nucleotide sequence variability than the supina-homeolog. Comparisons of read mappings of P. annua and a simulated P. supina × P. infirma hybrid showed high resemblance. Two homeolog-specific primer pairs were designed and used to amplify a 1860 bp region covering all resistance-conferring codons in the ALS gene. Four P. annua populations, GN, RB, GW and LG, showed high resistance to two ALS inhibitors, bispyribac-sodium and foramsulfuron, and two populations, HD and RS, showed lower resistance in the rate-response trial. Mutations conferring Trp-574-Leu substitution were observed in the infirma-homeolog of GN and RB and in the supina-homeolog of GW and LG, but no resistance-conferring mutation was observed in the two populations of lower resistance, HD and RS.

CONCLUSION

In this study we have demonstrated the use of NGS data to study homeologous polymorphisms, parentage and herbicide resistance in an allotetraploid weed species, P. annua. Complete coding sequences of the ALS gene were assembled for P. infirma, P. supina, infirma-homeolog and supina-homeolog in P. annua. A pipeline consisting of read mapping, SNP calling and allele frequency calculation was developed to study the parentage of P. annua, which provided a new perspective on this topic besides the views of morphology, karyotype and phylogeny. Our two homeolog-specific primer pairs can be utilized in future research to separate the homeologs of the ALS gene in P. annua and cover all the codons that have been reported to confer herbicide resistance.

Adaptive plasticity and niche expansion in an invasive thistle

Phenotypic differentiation in size and fecundity between native and invasive populations of a species has been suggested as a causal driver of invasion in plants. Local adaptation to novel environmental conditions through a micro-evolutionary response to natural selection may lead to phenotypic differentiation and fitness advantages in the invaded range. Local adaptation may occur along a stress tolerance trade-off, favoring individuals that, in benign conditions, shift resource allocation from stress tolerance to increased vigor and fecundity and, therefore, invasiveness. Alternately, the typically disturbed invaded range may select for a plastic, generalist strategy, making phenotypic plasticity the main driver of invasion success. To distinguish between these hypotheses, we performed a field common garden and tested for genetically based phenotypic differentiation, resource allocation shifts in response to water limitation, and local adaptation to the environmental gradient which describes the source locations for native and invasive populations of diffuse knapweed (Centaurea diffusa). Plants were grown in an experimental field in France (naturalized range) under water addition and limitation conditions. After accounting for phenotypic variation arising from environmental differences among collection locations, we found evidence of genetic variation between the invasive and native populations for most morphological and life-history traits under study. Invasive C. diffusa populations produced larger, later maturing, and therefore potentially fitter individuals than native populations. Evidence for local adaptation along a resource allocation trade-off for water limitation tolerance is equivocal. However, native populations do show evidence of local adaptation to an environmental gradient, a relationship which is typically not observed in the invaded range. Broader analysis of the climatic niche inhabited by the species in both ranges suggests that the physiological tolerances of C. diffusa may have expanded in the invaded range. This observation could be due to selection for plastic, "general-purpose" genotypes with broad environmental tolerances.

RNA-Seq analysis identifies genes associated with differential reproductive success under drought-stress in accessions of wild barley Hordeum spontaneum

BACKGROUND

The evolutionary basis of reproductive success in different environments is of major interest in the study of plant adaptation. Since the reproductive stage is particularly sensitive to drought, genes affecting reproductive success during this stage are key players in the evolution of adaptive mechanisms. We used an ecological genomics approach to investigate the reproductive response of drought-tolerant and sensitive wild barley accessions originating from different habitats in the Levant.

RESULTS

We sequenced mRNA extracted from spikelets at the flowering stage in drought-treated and control plants. The barley genome was used for a reference-guided assembly and differential expression analysis. Our approach enabled to detect biological processes affecting grain production under drought stress. We detected novel candidate genes and differentially expressed alleles associated with drought tolerance. Drought associated genes were shown to be more conserved than non-associated genes, and drought-tolerance genes were found to evolve more rapidly than other drought associated genes.

CONCLUSIONS

We show that reproductive success under drought stress is not a habitat-specific trait but a shared physiological adaptation that appeared to evolve recently in the evolutionary history of wild barley. Exploring the genomic basis of reproductive success under stress in crop wild progenitors is expected to have considerable ecological and economical applications.

Single Nucleotide Polymorphism Identification in Polyploids: A Review, Example, and Recommendations

Understanding the relationship between genotype and phenotype is a major biological question and being able to predict phenotypes based on molecular genotypes is integral to molecular breeding. Whole-genome duplications have shaped the history of all flowering plants and present challenges to elucidating the relationship between genotype and phenotype, especially in neopolyploid species. Although single nucleotide polymorphisms (SNPs) have become popular tools for genetic mapping, discovery and application of SNPs in polyploids has been difficult. Here, we summarize common experimental approaches to SNP calling, highlighting recent polyploid successes. To examine the impact of software choice on these analyses, we called SNPs among five peanut genotypes using different alignment programs (BWA-mem and Bowtie 2) and variant callers (SAMtools, GATK, and Freebayes). Alignments produced by Bowtie 2 and BWA-mem and analyzed in SAMtools shared 24.5% concordant SNPs, and SAMtools, GATK, and Freebayes shared 1.4% concordant SNPs. A subsequent analysis of simulated Brassica napus chromosome 1A and 1C genotypes demonstrated that, of the three software programs, SAMtools performed with the highest sensitivity and specificity on Bowtie 2 alignments. These results, however, are likely to vary among species, and we therefore propose a series of best practices for SNP calling in polyploids.

Genome scans reveal candidate domestication and improvement genes in cultivated sunflower, as well as post-domestication introgression with wild relatives

The development of modern crops typically involves both selection and hybridization, but to date most studies have focused on the former. In the present study, we explore how both processes, and their interactions, have molded the genome of the cultivated sunflower (Helianthus annuus), a globally important oilseed. To identify genes targeted by selection during the domestication and improvement of sunflower, and to detect post-domestication hybridization with wild species, we analyzed transcriptome sequences of 80 genotypes, including wild, landrace, and modern lines of H. annuus, as well as two cross-compatible wild relatives, Helianthus argophyllus and Helianthus petiolaris. Outlier analyses identified 122 and 15 candidate genes associated with domestication and improvement, respectively. As in several previous studies, genes putatively involved in oil biosynthesis were the most extreme outliers. Additionally, several promising associations were observed with previously mapped quantitative trait loci (QTLs), such as branching. Admixture analyses revealed that all the modern cultivar genomes we examined contained one or more introgressions from wild populations, with every chromosome having evidence of introgression in at least one modern line. Cumulatively, introgressions cover c. 10% of the cultivated sunflower genome. Surprisingly, introgressions do not avoid candidate domestication genes, probably because of the reintroduction of branching.

Seed fates in crop-wild hybrid sunflower: crop allele and maternal effects

Domestication has resulted in selection upon seed traits found in wild populations, yet crop-wild hybrids retain some aspects of both parental phenotypes. Seed fates of germination, dormancy, and mortality can influence the success of crop allele introgression in crop-wild hybrid zones, especially if crop alleles or crop-imparted seed coverings result in out-of-season germination. We performed a seed burial experiment using crop, wild, and diverse hybrid sunflower (Helianthus annuus) cross types to test how a cross type's maternal parent and nuclear genetic composition might affect its fate under field conditions. We observed higher maladaptive fall germination in the crop- and F₁- produced seeds than wild-produced seeds and, due to an interaction with percent crop alleles, fall germination was higher for cross types with more crop-like nuclear genetics. By spring, crop-produced cross types had the highest overwintering mortality, primarily due to higher fall germination. Early spring germination was identical across maternal types, but germination continued for F₁-produced seeds. In conclusion, the more wild-like the maternal parent or the less proportion of the cross type's genome contributed by the crop, the greater likelihood a seed will remain ungerminated than die. Wild-like dormancy may facilitate introgression through future recruitment from the soil seed bank.

Comparative genomics in the Asteraceae reveals little evidence for parallel evolutionary change in invasive taxa

Asteraceae, the largest family of flowering plants, has given rise to many notorious invasive species. Using publicly available transcriptome assemblies from 35 Asteraceae, including six major invasive species, we examined evidence for micro- and macro-evolutionary genomic changes associated with invasion. To detect episodes of positive selection repeated across multiple introductions, we conducted comparisons between native and introduced genotypes from six focal species and identified genes with elevated rates of amino acid change (dN/dS). We then looked for evidence of positive selection at a broader phylogenetic scale across all taxa. As invasive species may experience founder events during colonization and spread, we also looked for evidence of increased genetic load in introduced genotypes. We rarely found evidence for parallel changes in orthologous genes in the intraspecific comparisons, but in some cases we identified changes in members of the same gene family. Using among-species comparisons, we detected positive selection in 0.003-0.69% and 2.4-7.8% of the genes using site and stochastic branch-site models, respectively. These genes had diverse putative functions, including defence response, stress response and herbicide resistance, although there was no clear pattern in the GO terms. There was no indication that introduced genotypes have a higher proportion of deleterious alleles than native genotypes in the six focal species, suggesting multiple introductions and admixture mitigated the impact of drift. Our findings provide little evidence for common genomic responses in invasive taxa of the Asteraceae and hence suggest that multiple evolutionary pathways may lead to adaptation during introduction and spread in these species.

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

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.

Identification of candidate genes associated with leaf senescence in cultivated sunflower (Helianthus annuus L.)

Cultivated sunflower (Helianthus annuus L.), an important source of edible vegetable oil, shows rapid onset of senescence, which limits production by reducing photosynthetic capacity under specific growing conditions. Carbon for grain filling depends strongly on light interception by green leaf area, which diminishes during grain filling due to leaf senescence. Transcription factors (TFs) regulate the progression of leaf senescence in plants and have been well explored in model systems, but information for many agronomic crops remains limited. Here, we characterize the expression profiles of a set of putative senescence associated genes (SAGs) identified by a candidate gene approach and sunflower microarray expression studies. We examined a time course of sunflower leaves undergoing natural senescence and used quantitative PCR (qPCR) to measure the expression of 11 candidate genes representing the NAC, WRKY, MYB and NF-Y TF families. In addition, we measured physiological parameters such as chlorophyll, total soluble sugars and nitrogen content. The expression of Ha-NAC01, Ha-NAC03, Ha-NAC04, Ha-NAC05 and Ha-MYB01 TFs increased before the remobilization rate increased and therefore, before the appearance of the first physiological symptoms of senescence, whereas Ha-NAC02 expression decreased. In addition, we also examined the trifurcate feed-forward pathway (involving ORE1, miR164, and ethylene insensitive 2) previously reported for Arabidopsis. We measured transcription of Ha-NAC01 (the sunflower homolog of ORE1) and Ha-EIN2, along with the levels of miR164, in two leaves from different stem positions, and identified differences in transcription between basal and upper leaves. Interestingly, Ha-NAC01 and Ha-EIN2 transcription profiles showed an earlier up-regulation in upper leaves of plants close to maturity, compared with basal leaves of plants at pre-anthesis stages. These results suggest that the H. annuus TFs characterized in this work could play important roles as potential triggers of leaf senescence and thus can be considered putative candidate genes for senescence in sunflower.

Conservation and divergence of gene expression plasticity following c. 140 million years of evolution in lodgepole pine (Pinus contorta) and interior spruce (Picea glauca×Picea engelmannii)

Species respond to environmental stress through a combination of genetic adaptation and phenotypic plasticity, both of which may be important for survival in the face of climatic change. By characterizing the molecular basis of plastic responses and comparing patterns among species, it is possible to identify how such traits evolve. Here, we used de novo transcriptome assembly and RNAseq to explore how patterns of gene expression differ in response to temperature, moisture, and light regime treatments in lodgepole pine (Pinus contorta) and interior spruce (a natural hybrid population of Picea glauca and Picea engelmannii). We found wide evidence for an effect of treatment on expression within each species, with 6413 and 11,658 differentially expressed genes identified in spruce and pine, respectively. Comparing patterns of expression among these species, we found that 74% of all orthologs with differential expression had a pattern that was conserved in both species, despite 140 million yr of evolution. We also found that the specific treatments driving expression patterns differed between genes with conserved versus diverged patterns of expression. We conclude that natural selection has probably played a role in shaping plastic responses to environment in these species.

Selection of candidate reference genes for real-time PCR studies in lettuce under abiotic stresses

The process of selection and validation of reference genes is the first step in studies of gene expression by real-time quantitative polymerase chain reaction (RT-qPCR). The genome of lettuce, the most popular leaf vegetable cultivated worldwide, has recently been sequenced; therefore, suitable reference genes for reliable results in RT-qPCR analyses are required. In the present study, 17 candidate reference genes were selected, and their expression stability in lettuce leaves under drought, salt, heavy metal, and UV-C irradiation conditions and under the application of abscisic acid (ABA) was evaluated using geNorm and NormFinder software. The candidate reference genes included protein-coding traditional and novel reference genes and microRNAs (miRNAs). The results indicate that the expression stability is dependent on the experimental conditions. The novel protein-coding reference genes were more suitable than the traditional reference genes under drought, UV-C irradiation, and heavy metal conditions and under the application of ABA. Only under salinity conditions were the traditional protein-coding reference genes more stable than the novel genes. In addition, the miRNAs, mainly MIR169, MIR171/170 and MIR172, were stably expressed under the abiotic stresses evaluated, representing a suitable alternative approach for gene expression data normalization. The expression of phenylalanine ammonia lyase (PAL) and 4-hydroxyphenylpyruvate dioxygenase (HPPD) was used to further confirm the validated protein-coding reference genes, and the expression of MIR172 and MIR398 was used to confirm the validated miRNA genes, showing that the use of an inappropriate reference gene induces erroneous results. This work is the first survey of the stability of reference genes in lettuce and provides guidelines to obtain more accurate RT-qPCR results in lettuce studies.

Rapid evolution of an invasive weed

Trade-offs between performance and the ability to tolerate abiotic and biotic stress have been suggested to explain both the success of invasive species and phenotypic differentiation between native and invasive populations. It is critical to sample broadly across both ranges and to account for latitudinal clines and maternal effects when testing this premise. Wild-collected Centaurea diffusa seeds were grown in benign and stressful conditions (drought, flooding, nutrient stress and simulated herbivory), to evaluate whether native and introduced individuals differ in performance or life history phenotypes. A second experiment used glasshouse-grown seeds to evaluate whether patterns remain comparable when maternal environment is consistent. Many traits differed between ranges, and in all cases but one, invasive individuals grew larger, performed better, or matured later. No trade-off in performance with herbivore defense was evident. Invasive populations may have been released from a trade-off between growth and drought tolerance apparent in the native range. Larger individuals with delayed maturity and greater reproductive potential have evolved in invasive populations, a pattern evident across broad population sampling, and after latitude and maternal environment were considered. Release from abiotic stress tolerance trade-offs may be important for the invasion success of Centaurea diffusa.

Maintenance of strong morphological differentiation despite ongoing natural hybridization between sympatric species of Lomatia (Proteaceae)

BACKGROUND AND AIMS

When species cohesion is maintained despite ongoing natural hybridization, many questions are raised about the evolutionary processes operating in the species complex. This study examined the extensive natural hybridization between the Australian native shrubs Lomatia myricoides and L. silaifolia (Proteaceae). These species exhibit striking differences in morphology and ecological preferences, exceeding those found in most studies of hybridization to date.

METHODS

Nuclear microsatellite markers (nSSRs), genotyping methods and morphometric analyses were used to uncover patterns of hybridization and the role of gene flow in morphological differentiation between sympatric species.

KEY RESULTS

The complexity of hybridization patterns differed markedly between sites, however, signals of introgression were present at all sites. One site provided evidence of a large hybrid swarm and the likely presence of multiple hybrid generations and backcrosses, another site a handful of early generational hybrids and a third site only traces of admixture from a past hybridization event. The presence of cryptic hybrids and a pattern of morphological bimodality amongst hybrids often disguised the extent of underlying genetic admixture.

CONCLUSIONS

Distinct parental habitats and phenotypes are expected to form barriers that contribute to the rapid reversion of hybrid populations to their parental character state, due to limited opportunities for hybrid/intermediate advantage. Furthermore, strong genomic filters may facilitate continued gene flow between species without the danger of assimilation. Stochastic fire events facilitate temporal phenological isolation between species and may partly explain the bi-directional and site-specific patterns of hybridization observed. Furthermore, the findings suggest that F1 hybrids are rare, and backcrosses may occur rapidly following these initial hybridization events.

Genetic analysis of safflower domestication

BACKGROUND

Safflower (Carthamus tinctorius L.) is an oilseed crop in the Compositae (a.k.a. Asteraceae) that is valued for its oils rich in unsaturated fatty acids. Here, we present an analysis of the genetic architecture of safflower domestication and compare our findings to those from sunflower (Helianthus annuus L.), an independently domesticated oilseed crop within the same family.We mapped quantitative trait loci (QTL) underlying 24 domestication-related traits in progeny from a cross between safflower and its wild progenitor, Carthamus palaestinus Eig. Also, we compared QTL positions in safflower against those that have been previously identified in cultivated x wild sunflower crosses to identify instances of colocalization.

RESULTS

We mapped 61 QTL, the vast majority of which (59) exhibited minor or moderate phenotypic effects. The two large-effect QTL corresponded to one each for flower color and leaf spininess. A total of 14 safflower QTL colocalized with previously reported sunflower QTL for the same traits. Of these, QTL for three traits (days to flower, achene length, and number of selfed seed) had cultivar alleles that conferred effects in the same direction in both species.

CONCLUSIONS

As has been observed in sunflower, and unlike many other crops, our results suggest that the genetics of safflower domestication is quite complex. Moreover, our comparative mapping results indicate that safflower and sunflower exhibit numerous instances of QTL colocalization, suggesting that parallel trait transitions during domestication may have been driven, at least in part, by parallel genotypic evolution at some of the same underlying genes.

Genomic islands of divergence are not affected by geography of speciation in sunflowers

Genomic studies of speciation often report the presence of highly differentiated genomic regions interspersed within a milieu of weakly diverged loci. The formation of these speciation islands is generally attributed to reduced inter-population gene flow near loci under divergent selection, but few studies have critically evaluated this hypothesis. Here, we report on transcriptome scans among four recently diverged pairs of sunflower (Helianthus) species that vary in the geographical context of speciation. We find that genetic divergence is lower in sympatric and parapatric comparisons, consistent with a role for gene flow in eroding neutral differences. However, genomic islands of divergence are numerous and small in all comparisons, and contrary to expectations, island number and size are not significantly affected by levels of interspecific gene flow. Rather, island formation is strongly associated with reduced recombination rates. Overall, our results indicate that the functional architecture of genomes plays a larger role in shaping genomic divergence than does the geography of speciation.

Detecting and locating whole genome duplications on a phylogeny: a probabilistic approach

Whole genome duplications (WGDs) followed by massive gene loss occurred in the evolutionary history of many groups. WGDs are usually inferred from the age distribution of paralogs (Ks-based methods) or from gene collinearity data (synteny). However, Ks-based methods are restricted to detect the recent WGDs due to saturation effects and the difficulty to date old duplicates, and synteny is difficult to reconstruct for distantly related species. Recently, Jiao et al. (Jiao Y, Wickett N, Ayyampalayam S, Chanderbali AS, Landherr L, Ralph PE, Tomsho LP, Hu Y, Liang H, Soltis PS, et al. 2011. Ancestral polyploidy in seed plants and angiosperms. Nature 473:97-100) introduced an empirical method that aims to detect a peak in duplication ages among nodes selected from a previous phylogenetic analysis. In this context, we present here two rigorous methods based on data from multiple gene families and on a new probabilistic model. Our model assumes that all gene lineages are instantaneously duplicated at the WGD event with a possible almost-immediate loss of some extra copies. Our reconciliation method relies on aligned molecular sequences, whereas our gene count method relies only on gene count data across species. We show, using extensive simulations, that both methods have a good detection power. Surprisingly, the gene count method enjoys no loss of power compared with the reconciliation method, despite the fact that sequence information is not used. We finally illustrate the performance of our methods on a benchmark yeast data set. Both methods are able to detect the well-known WGD in the Saccharomyces cerevisiae clade and agree on a small retention rate at the WGD, as established by synteny-based methods.

De novo assembly and characterization of the Barnyardgrass (Echinochloa crus-galli) transcriptome using next-generation pyrosequencing

Background

Barnyardgrass (Echinochloa crus-galli) is an important weed that is a menace to rice cultivation and production. Rapid evolution of herbicide resistance in this weed makes it one of the most difficult to manage using herbicides. Since genome-wide sequence data for barnyardgrass is limited, we sequenced the transcriptomes of susceptible and resistant barnyardgrass biotypes using the 454 GS-FLX platform.

Results

454 pyrosequencing generated 371,281 raw reads with an average length of 341.8 bp, which made a total length of 126.89 Mb (SRX160526). De novo assembly produced 10,142 contigs (∼5.92 Mb) with an average length of 583 bp and 68,940 singletons (∼22.13 Mb) with an average length of 321 bp. About 244,653 GO term assignments to the biological process, cellular component and molecular function categories were obtained. A total of 6,092 contigs and singletons with 2,515 enzyme commission numbers were assigned to 151 predicted KEGG metabolic pathways. Digital abundance analysis using Illumina sequencing identified 78,124 transcripts among susceptible, resistant, herbicide-treated susceptible and herbicide-treated resistant barnyardgrass biotypes. From these analyses, eight herbicide target-site gene groups and four non-target-site gene groups were identified in the resistant biotype. These could be potential candidate genes involved in the herbicide resistance of barnyardgrass and could be used for further functional genomics research. C₄ photosynthesis genes including RbcS, RbcL, NADP-me and MDH with complete CDS were identified using PCR and RACE technology.

Conclusions

This is the first large-scale transcriptome sequencing of E. crus-galli performed using the 454 GS-FLX platform. Potential candidate genes involved in the evolution of herbicide resistance were identified from the assembled sequences. This transcriptome data may serve as a reference for further gene expression and functional genomics studies, and will facilitate the study of herbicide resistance at the molecular level in this species as well as other weeds.

Transcriptome divergence between introduced and native populations of Canada thistle, Cirsium arvense

Introduced plants may quickly evolve new adaptive traits upon their introduction. Canada thistle (Cirsium arvense - Cardueae, Asteraceae) is one of the worst invasive weeds worldwide. The goal of this study is to compare gene expression profiles of native (European) and introduced (North American) populations of this species, to elucidate the genetic mechanisms that may underlie such rapid adaptation. We explored the transcriptome of ten populations (five per range) of C. arvense in response to three treatments (control, nutrient deficiency and shading) using a customized microarray chip containing 63 690 expressed sequence tags (ESTs), and verified the expression level of 13 loci through real-time quantitative PCR. Only 2116 ESTs (3.5%) were found to be differentially expressed between the ranges, and 4458 ESTs (7.1%) exhibited a significant treatment-by-range effect. Among them was an overrepresentation of loci involved in stimulus and stress responses. Cirsium arvense has evolved different life history strategies on each continent. The two ranges notably differ with regard to R-protein mediated defence, sensitivity to abiotic stresses, and developmental timing. The fact that genotypes from the Midwest exhibit different expression kinetics than remaining North American samples further corroborates the hypothesis that the New World has been colonized twice, independently.

Reprint of: using nuclear gene data for plant phylogenetics: progress and prospects

The paper reviews the current state of low and single copy nuclear markers that have been applied successfully in plant phylogenetics to date, and discusses case studies highlighting the potential of massively parallel high throughput or next-generation sequencing (NGS) approaches for molecular phylogenetic and evolutionary investigations. The current state, prospects and challenges of specific single- or low-copy plant nuclear markers as well as phylogenomic case studies are presented and evaluated.

Allele identification for transcriptome-based population genomics in the invasive plant Centaurea solstitialis

Transcriptome sequences are becoming more broadly available for multiple individuals of the same species, providing opportunities to derive population genomic information from these datasets. Using the 454 Life Science Genome Sequencer FLX and FLX-Titanium next-generation platforms, we generated 11−430 Mbp of sequence for normalized cDNA for 40 wild genotypes of the invasive plant Centaurea solstitialis, yellow starthistle, from across its worldwide distribution. We examined the impact of sequencing effort on transcriptome recovery and overlap among individuals. To do this, we developed two novel publicly available software pipelines: SnoWhite for read cleaning before assembly, and AllelePipe for clustering of loci and allele identification in assembled datasets with or without a reference genome. AllelePipe is designed specifically for cases in which read depth information is not appropriate or available to assist with disentangling closely related paralogs from allelic variation, as in transcriptome or previously assembled libraries. We find that modest applications of sequencing effort recover most of the novel sequences present in the transcriptome of this species, including single-copy loci and a representative distribution of functional groups. In contrast, the coverage of variable sites, observation of heterozygosity, and overlap among different libraries are all highly dependent on sequencing effort. Nevertheless, the information gained from overlapping regions was informative regarding coarse population structure and variation across our small number of population samples, providing the first genetic evidence in support of hypothesized invasion scenarios.

The molecular basis of invasiveness: differences in gene expression of native and introduced common ragweed (Ambrosia artemisiifolia) in stressful and benign environments

Although the evolutionary and ecological processes that contribute to plant invasion have been the focus of much research, investigation into the molecular basis of invasion is just beginning. Common ragweed (Ambrosia artemisiifolia) is an annual weed native to North America and has been introduced to Europe where it has become invasive. Using a custom-designed NimbleGen oligoarray, we examined differences in gene expression between five native and six introduced populations of common ragweed in three different environments (control, light stress and nutrient stress), as well as two different time points. We identified candidate genes that may contribute to invasiveness in common ragweed based on differences in expression between native and introduced populations from Europe. Specifically, we found 180 genes where range explained a significant proportion of the variation in gene expression and a further 103 genes with a significant range by treatment interaction. Several of these genes are potentially involved in the metabolism of secondary compounds, stress response and the detoxification of xenobiotics. Previously, we found more rapid growth and greater reproductive success in introduced populations, particularly in benign and competitive (light stress) environments, and many of these candidate genes potentially underlie these growth differences. We also found expression differences among populations within each range, reflecting either local adaptation or neutral processes, although no associations with climate or latitude were identified. These data provide a first step in identifying genes that are involved with introduction success in an aggressive annual weed.

RNA-seq analysis of allele-specific expression, hybrid effects, and regulatory divergence in hybrids compared with their parents from natural populations

Hybridization is a prominent process among natural plant populations that can result in phenotypic novelty, heterosis, and changes in gene expression. The effects of intraspecific hybridization on F1 hybrid gene expression were investigated using parents from divergent, natural populations of Cirsium arvense, an invasive Compositae weed. Using an RNA-seq approach, the expression of 68,746 unigenes was quantified in parents and hybrids. The expression levels of 51% of transcripts differed between parents, a majority of which had less than 1.25× fold-changes. More unigenes had higher expression in the invasive parent (P1) than the noninvasive parent (P2). Of those that were divergently expressed between parents, 10% showed additive and 81% showed nonadditive (transgressive or dominant) modes of gene action in the hybrids. A majority of the dominant cases had P2-like expression patterns in the hybrids. Comparisons of allele-specific expression also enabled a survey of cis- and trans-regulatory effects. Cis- and trans-regulatory divergence was found at 70% and 68% of 62,281 informative single-nucleotide polymorphism sites, respectively. Of the 17% of sites exhibiting both cis- and trans-effects, a majority (70%) had antagonistic regulatory interactions (cis x trans); trans-divergence tended to drive higher expression of the P1 allele, whereas cis-divergence tended to increase P2 transcript abundance. Trans-effects correlated more highly than cis with parental expression divergence and accounted for a greater proportion of the regulatory divergence at sites with additive compared with nonadditive inheritance patterns. This study explores the nature of, and types of mechanisms underlying, expression changes that occur in upon intraspecific hybridization in natural populations.

The red queen in the corn: agricultural weeds as models of rapid adaptive evolution

Weeds are among the greatest pests of agriculture, causing billions of dollars in crop losses each year. As crop field management practices have changed over the past 12 000 years, weeds have adapted in turn to evade human removal. This evolutionary change can be startlingly rapid, making weeds an appealing system to study evolutionary processes that occur over short periods of time. An understanding of how weeds originate and adapt is needed for successful management; however, relatively little emphasis has been placed on genetically characterizing these systems. Here, we review the current literature on agricultural weed origins and their mechanisms of adaptation. Where possible, we have included examples that have been genetically well characterized. Evidence for three possible, non-mutually exclusive weed origins (from wild species, crop-wild hybrids or directly from crops) is discussed with respect to what is known about the microevolutionary signatures that result from these processes. We also discuss what is known about the genetic basis of adaptive traits in weeds and the range of genetic mechanisms that are responsible. With a better understanding of genetic mechanisms underlying adaptation in weedy species, we can address the more general process of adaptive evolution and what can be expected as we continue to apply selective pressures in agroecosystems around the world.