Characterization of de novo transcriptome for waterhemp (Amaranthus tuberculatus) using GS-FLX 454 pyrosequencing and its application for studies of herbicide target-site genes

Weed biology and management in the multi-omics era: Progress and perspectives

Weeds pose a significant threat to crop production, resulting in substantial yield reduction. In addition, they possess robust weedy traits that enable them to survive in extreme environments and evade human control. In recent years, the application of multi-omics biotechnologies has helped to reveal the molecular mechanisms underlying these weedy traits. In this review, we systematically describe diverse applications of multi-omics platforms for characterizing key aspects of weed biology, including the origins of weed species, weed classification, and the underlying genetic and molecular bases of important weedy traits such as crop-weed interactions, adaptability to different environments, photoperiodic flowering responses, and herbicide resistance. In addition, we discuss limitations to the application of multi-omics techniques in weed science, particularly compared with their extensive use in model plants and crops. In this regard, we provide a forward-looking perspective on the future application of multi-omics technologies to weed science research. These powerful tools hold great promise for comprehensively and efficiently unraveling the intricate molecular genetic mechanisms that underlie weedy traits. The resulting advances will facilitate the development of sustainable and highly effective weed management strategies, promoting greener practices in agriculture.

Comparative analysis of dioecious Amaranthus plastomes and phylogenomic implications within Amaranthaceae s.s

BACKGROUND

The genus Amaranthus L. consists of 70-80 species distributed across temperate and tropical regions of the world. Nine species are dioecious and native to North America; two of which are agronomically important weeds of row crops. The genus has been described as taxonomically challenging and relationships among species including the dioecious ones are poorly understood. In this study, we investigated the phylogenetic relationships among the dioecious amaranths and sought to gain insights into plastid tree incongruence. A total of 19 Amaranthus species' complete plastomes were analyzed. Among these, seven dioecious Amaranthus plastomes were newly sequenced and assembled, an additional two were assembled from previously published short reads sequences and 10 other plastomes were obtained from a public repository (GenBank).

RESULTS

Comparative analysis of the dioecious Amaranthus species' plastomes revealed sizes ranged from 150,011 to 150,735 bp and consisted of 112 unique genes (78 protein-coding genes, 30 transfer RNAs and 4 ribosomal RNAs). Maximum likelihood trees, Bayesian inference trees and splits graphs support the monophyly of subgenera Acnida (7 dioecious species) and Amaranthus; however, the relationship of A. australis and A. cannabinus to the other dioecious species in Acnida could not be established, as it appears a chloroplast capture occurred from the lineage leading to the Acnida + Amaranthus clades. Our results also revealed intraplastome conflict at some tree branches that were in some cases alleviated with the use of whole chloroplast genome alignment, indicating non-coding regions contribute valuable phylogenetic signals toward shallow relationship resolution. Furthermore, we report a very low evolutionary distance between A. palmeri and A. watsonii, indicating that these two species are more genetically related than previously reported.

CONCLUSIONS

Our study provides valuable plastome resources as well as a framework for further evolutionary analyses of the entire Amaranthus genus as more species are sequenced.

Genomic profiling of dioecious Amaranthus species provides novel insights into species relatedness and sex genes

BACKGROUND

Amaranthus L. is a diverse genus consisting of domesticated, weedy, and non-invasive species distributed around the world. Nine species are dioecious, of which Amaranthus palmeri S. Watson and Amaranthus tuberculatus (Moq.) J.D. Sauer are troublesome weeds of agronomic crops in the USA and elsewhere. Shallow relationships among the dioecious Amaranthus species and the conservation of candidate genes within previously identified A. palmeri and A. tuberculatus male-specific regions of the Y (MSYs) in other dioecious species are poorly understood. In this study, seven genomes of dioecious amaranths were obtained by paired-end short-read sequencing and combined with short reads of seventeen species in the family Amaranthaceae from NCBI database. The species were phylogenomically analyzed to understand their relatedness. Genome characteristics for the dioecious species were evaluated and coverage analysis was used to investigate the conservation of sequences within the MSY regions.

RESULTS

We provide genome size, heterozygosity, and ploidy level inference for seven newly sequenced dioecious Amaranthus species and two additional dioecious species from the NCBI database. We report a pattern of transposable element proliferation in the species, in which seven species had more Ty3 elements than copia elements while A. palmeri and A. watsonii had more copia elements than Ty3 elements, similar to the TE pattern in some monoecious amaranths. Using a Mash-based phylogenomic analysis, we accurately recovered taxonomic relationships among the dioecious Amaranthus species that were previously identified based on comparative morphology. Coverage analysis revealed eleven candidate gene models within the A. palmeri MSY region with male-enriched coverages, as well as regions on scaffold 19 with female-enriched coverage, based on A. watsonii read alignments. A previously reported FLOWERING LOCUS T (FT) within A. tuberculatus MSY contig was also found to exhibit male-enriched coverages for three species closely related to A. tuberculatus but not for A. watsonii reads. Additional characterization of the A. palmeri MSY region revealed that 78% of the region is made of repetitive elements, typical of a sex determination region with reduced recombination.

CONCLUSIONS

The results of this study further increase our understanding of the relationships among the dioecious species of the Amaranthus genus as well as revealed genes with potential roles in sex function in the species.

Can allele-specific loop-mediated isothermal amplification be used for rapid detection of target-site herbicide resistance in Lolium spp.?

BACKGROUND

Herbicide resistance is one of the threats to modern agriculture and its early detection is one of the most effective components for sustainable resistance management strategies. Many techniques have been used for target-site-resistance detection. Allele-Specific Loop-Mediated Isothermal Amplification (AS-LAMP) was evaluated as a possible rapid diagnostic method for acetyl-CoA carboxylase (ACCase) and acetolactate synthase (ALS) inhibiting herbicides resistance in Lolium spp.

RESULTS

AS-LAMP protocols were set up for the most frequent mutations responsible for herbicide resistance to ALS (positions 197, 376 and 574) and ACCase (positions 1781, 2041 and 2078) inhibitors in previously characterized and genotyped Lolium spp.

POPULATIONS

A validation step on new putative resistant populations gave the overview of a possible use of this tool for herbicide resistance diagnosis in Lolium spp. Regarding the ACCase inhibitor pinoxaden, in more than 65% of the analysed plants, the LAMP assay and genotyping were in keeping, whereas the results were not consistent when ALS inhibitors resistance was considered. Limitations on the use of this technique for herbicide resistance detection in the allogamous Lolium spp. are discussed.

CONCLUSIONS

The LAMP method used for the detection of target-site resistance in weed species could be applicable with target genes that do not have high genetic variability, such as ACCase gene in Lolium spp.

Genetic variation and structure of complete chloroplast genome in alien monoecious and dioecious Amaranthus weeds

Amaranthus is a complex taxon with economic importance as well as harmful weeds. We studied the genetic variation and structure of the chloroplast genomes of 22 samples from 17 species of three subgenera. It was found that the length of the chloroplast genome of Amaranthus varied from 149,949 bp of A. polygonoides to 150,757 bp of A. albus. The frequencies of SNPs and InDels in chloroplast genomes were 1.79% and 2.86%, and the variation mainly occurred in the non-coding regions. The longest InDel was 387 bp, which occurred on ycf2, followed by 384 bp InDel on psbM-trnD. Two InDels in ndhE-I on the SSC make the three subgenera clearly distinguished. In LSC, SSC and IRs regions, there were four 30 bp forward and reverse repeats, and the repeats in SSC and LSC were in nearly opposite positions in circular genome structure, and almost divided the circular genome into symmetrical structures. In the topological tree constructed by chloroplast genome, species in subgen. Amaranthus and subgen. Acnida form monophyletic branches separately and cluster together. A. albus, A. blitoides and A. polygonoides were separated from subgen. Albersia, and the rest of subgen. Albersia were clustered into a monophyletic branch. The rpoC2, ycf1, ndhF-rpl32 were good at distinguishing most amaranths. The trnk-UUU-atpF, trnT-UGU-atpB, psbE-clpP, rpl14-rps19, and ndhF-D can distinguish several similar species. In general, the chloroplast genome is of certain value for the identification of the similar species of Amaranthus, which provides more evidence for clarifying the phylogenetic relationships within the genus.

The genetic architecture and population genomic signatures of glyphosate resistance in Amaranthus tuberculatus

Much of what we know about the genetic basis of herbicide resistance has come from detailed investigations of monogenic adaptation at known target-sites, despite the increasingly recognized importance of polygenic resistance. Little work has been done to characterize the broader genomic basis of herbicide resistance, including the number and distribution of genes involved, their effect sizes, allele frequencies and signatures of selection. In this work, we implemented genome-wide association (GWA) and population genomic approaches to examine the genetic architecture of glyphosate (Round-up) resistance in the problematic agricultural weed Amaranthus tuberculatus. A GWA was able to correctly identify the known target-gene but statistically controlling for two causal target-site mechanisms revealed an additional 250 genes across all 16 chromosomes associated with non-target-site resistance (NTSR). The encoded proteins had functions that have been linked to NTSR, the most significant of which is response to chemicals, but also showed pleiotropic roles in reproduction and growth. Compared to an empirical null that accounts for complex population structure, the architecture of NTSR was enriched for large effect sizes and low allele frequencies, suggesting the role of pleiotropic constraints on its evolution. The enrichment of rare alleles also suggested that the genetic architecture of NTSR may be population-specific and heterogeneous across the range. Despite their rarity, we found signals of recent positive selection on NTSR-alleles by both window- and haplotype-based statistics, and an enrichment of amino acid changing variants. In our samples, genome-wide single nucleotide polymorphisms explain a comparable amount of the total variation in glyphosate resistance to monogenic mechanisms, even in a collection of individuals where 80% of resistant individuals have large-effect TSR mutations, indicating an underappreciated polygenic contribution to the evolution of herbicide resistance in weed populations.

Herbicide resistance in Amaranthus tuberculatus†

Amaranthus tuberculatus is the major weed species in many midwestern US row-crop production fields, and it is among the most problematic weeds in the world in terms of its ability to evolve herbicide resistance. It has now evolved resistance to herbicides spanning seven unique sites of action, with populations and even individual plants often possessing resistance to several herbicides/herbicide groups. Historically, herbicide target-site changes accounted for most of the known resistance mechanisms in this weed; however, over the last few years, non-target-site mechanisms, particularly enhanced herbicide detoxification, have become extremely common in A. tuberculatus. Unravelling the genetics and molecular details of non-target-site resistance mechanisms, understanding the extent to which they confer cross resistance to other herbicides, and understanding how they evolve remain as critical research endeavors. Transcriptomic and genomics approaches are already facilitating such studies, the results of which hopefully will inform better resistance-mitigation strategies. The largely unprecedented level of herbicide resistance in A. tuberculatus is not only a fascinating example of evolution in action, but it is a serious and growing threat to the sustainability of midwestern US cropping systems. © 2020 Society of Chemical Industry.

Draft Genomes of Amaranthus tuberculatus, Amaranthus hybridus, and Amaranthus palmeri

Amaranthus tuberculatus, Amaranthus hybridus, and Amaranthus palmeri are agronomically important weed species. Here, we present the most contiguous draft assemblies of these three species to date. We utilized a combination of Pacific Biosciences long-read sequencing and chromatin contact mapping information to assemble and order sequences of A. palmeri to near-chromosome-level resolution, with scaffold N50 of 20.1 Mb. To resolve the issues of heterozygosity and coassembly of alleles in diploid species, we adapted the trio binning approach to produce haplotype assemblies of A. tuberculatus and A. hybridus. This approach resulted in an improved assembly of A. tuberculatus, and the first genome assembly for A. hybridus, with contig N50s of 2.58 and 2.26 Mb, respectively. Species-specific transcriptomes and information from related species were used to predict transcripts within each assembly. Syntenic comparisons of these species and Amaranthus hypochondriacus identified sites of genomic rearrangement, including duplication and translocation, whereas genetic map construction within A. tuberculatus highlighted the need for further ordering of the A. hybridus and A. tuberculatus contigs. These multiple reference genomes will accelerate genomic studies in these species to further our understanding of weedy evolution within Amaranthus.

Identification and expression of main genes involved in non-target site resistance mechanisms to fenoxaprop-p-ethyl in Beckmannia syzigachne

BACKGROUND

Non-target-site resistance (NTSR) to herbicides is a serious threat to global agriculture. Although metabolic resistance is the dominant mechanism of NTSR, the molecular mechanisms are not yet well-characterized. This study aimed to uncover the likely metabolism-related genes in Beckmannia syzigachne (American sloughgrass) resistant to fenoxaprop-p-ethyl.

RESULTS

Ultra-performance liquid chromatography - tandem mass spectrometry experiments showed that the resistant American sloughgrass biotype (R, SD-04-SS) showed enhanced degradation of this herbicide compared to the susceptible biotype (S, SD-12). R and S biotype were harvested at 24 h after fenoxaprop-p-ethyl treatment to conduct RNA sequencing (RNA-Seq) analysis to investigate the likely fenoxaprop-p-ethyl metabolic genes. The RNA-Seq libraries yield 417 969 980 clean reads. The de novo assembly generated 115 112 unigenes, of which 57 906 unigenes were annotated. Finally, we identified 273 cytochrome P450s, 178 oxidases, 47 glutathione S-transferases (GSTs), 166 glucosyltransferases (GTs) and 180 ABC transporter genes to determine the likely fenoxaprop-p-ethyl metabolism-related genes in R biotype. Twelve overlapping up-regulated genes in the R biotype (fenoxaprop-p-ethyl-treated R/non-treated R, fenoxaprop-p-ethyl-treated R/fenoxaprop-p-ethyl-treated S) were identified by RNA-Seq and the results were validated using qRT-PCR. Ten were identified as fenoxaprop-p-ethyl metabolism-related genes, including three P450s (homologous to CYP71D7, CYP99A2 and CYP71D10), one GST (homologous to GSTF1), two GTs (homologous to UGT90A1 and UGT83A1) and four oxidase genes.

CONCLUSION

This work demonstrates that the NTSR mechanism by means of enhanced detoxification of fenoxaprop-p-ethyl in American sloughgrass is very likely driven by herbicide metabolism related genes. The RNA-Seq data presented here provide a valuable resource for understanding the molecular mechanism of NTSR in American sloughgrass. © 2020 Society of Chemical Industry.

Species identification, phylogenetic analysis and detection of herbicide-resistant biotypes of Amaranthus based on ALS and ITS

The taxonomically challenging genus Amaranthus (Family Amaranthaceae) includes important agricultural weed species that are being spread globally as grain contaminants. We hypothesized that the ALS gene will help resolve these taxonomic challenges and identify potentially harmful resistant biotypes. We obtained 153 samples representing 26 species from three Amaranthus subgenera and included in that incorporated ITS, ALS (domains C, A and D) and ALS (domains B and E) sequences. Subgen. Albersia was well supported, but subgen. Amaranthus and subgen. Acnida were not. Amaranthus tuberculatus, A. palmeri and A. spinosus all showed different genetic structuring. Unique SNPs in ALS offered reliable diagnostics for most of the sampled Amaranthus species. Resistant ALS alleles were detected in sixteen A. tuberculatus samples (55.2%), eight A. palmeri (27.6%) and one A. arenicola (100%). These involved Ala₁₂₂Asn, Pro₁₉₇Ser/Thr/Ile, Trp₅₇₄Leu, and Ser₆₅₃Thr/Asn/Lys substitutions, with Ala₁₂₂Asn, Pro₁₉₇Thr/Ile and Ser₆₅₃Lys being reported in Amaranthus for the first time. Moreover, different resistant mutations were present in different A. tuberculatus populations. In conclusion, the ALS gene is important for species identification, investigating population genetic diversity and understanding resistant evolution within the genus Amaranthus.

Synthesis and herbicidal activities of aryloxyacetic acid derivatives as HPPD inhibitors

A series of aryloxyacetic acid derivatives were designed and synthesized as 4-hydoxyphenylpyruvate dioxygenase (HPPD) inhibitors. Preliminary bioassay results reveal that these derivatives are promising Arabidopsis thaliana HPPD (AtHPPD) inhibitors, in particular compounds I12 (K_i = 0.011 µM) and I23 (K_i = 0.012 µM), which exhibit similar activities to that of mesotrione, a commercial HPPD herbicide (K_i = 0.013 µM). Furthermore, the newly synthesized compounds show significant greenhouse herbicidal activities against tested weeds at dosages of 150 g ai/ha. In particular, II4 exhibited high herbicidal activity for pre-emergence treatment that was slightly better than that of mesotrione. In addition, compound II4 was safe for weed control in maize fields at a rate of 150 g ai/ha, and was identified as the most potent candidate for a novel HPPD inhibitor herbicide. The compounds described herein may provide useful guidance for the design of new HPPD inhibiting herbicides and their modification.

Multiple modes of convergent adaptation in the spread of glyphosate-resistant Amaranthus tuberculatus

While evolution has been thought of as playing out over millions of years, adaptation to new environments can occur very rapidly, presenting us with key opportunities to understand evolutionary dynamics. One of the most amazing examples of real-time evolution comes from agriculture, where due to the intense use of a few herbicides, many plant species have evolved herbicide resistance to become aggressive weeds. An important question has been whether herbicide resistance arises only rarely and then spreads quickly, or whether herbicide resistance arises all the time de novo. Our work with glyphosate resistance in US Midwestern and Canadian populations of Amaranthus tuberculatus reveals the answer to be, “it depends,” as we surprisingly find examples for both modes of evolution.

The selection pressure exerted by herbicides has led to the repeated evolution of herbicide resistance in weeds. The evolution of herbicide resistance on contemporary timescales in turn provides an outstanding opportunity to investigate key questions about the genetics of adaptation, in particular the relative importance of adaptation from new mutations, standing genetic variation, or geographic spread of adaptive alleles through gene flow. Glyphosate-resistant Amaranthus tuberculatus poses one of the most significant threats to crop yields in the Midwestern United States, with both agricultural populations and herbicide resistance only recently emerging in Canada. To understand the evolutionary mechanisms driving the spread of resistance, we sequenced and assembled the A. tuberculatus genome and investigated the origins and population genomics of 163 resequenced glyphosate-resistant and susceptible individuals from Canada and the United States. In Canada, we discovered multiple modes of convergent evolution: in one locality, resistance appears to have evolved through introductions of preadapted US genotypes, while in another, there is evidence for the independent evolution of resistance on genomic backgrounds that are historically nonagricultural. Moreover, resistance on these local, nonagricultural backgrounds appears to have occurred predominantly through the partial sweep of a single haplotype. In contrast, resistant haplotypes arising from the Midwestern United States show multiple amplification haplotypes segregating both between and within populations. Therefore, while the remarkable species-wide diversity of A. tuberculatus has facilitated geographic parallel adaptation of glyphosate resistance, more recently established agricultural populations are limited to adaptation in a more mutation-limited framework.

The Triple Amino Acid Substitution TAP-IVS in the EPSPS Gene Confers High Glyphosate Resistance to the Superweed Amaranthus hybridus

The introduction of glyphosate-resistant (GR) crops revolutionized weed management; however, the improper use of this technology has selected for a wide range of weeds resistant to glyphosate, referred to as superweeds. We characterized the high glyphosate resistance level of an Amaranthus hybridus population (GRH)-a superweed collected in a GR-soybean field from Cordoba, Argentina-as well as the resistance mechanisms that govern it in comparison to a susceptible population (GSH). The GRH population was 100.6 times more resistant than the GSH population. Reduced absorption and metabolism of glyphosate, as well as gene duplication of 5-enolpyruvylshikimate-3-phosphate synthase (EPSPS) or its overexpression did not contribute to this resistance. However, GSH plants translocated at least 10% more 14C-glyphosate to the rest of the plant and roots than GRH plants at 9 h after treatment. In addition, a novel triple amino acid substitution from TAP (wild type, GSH) to IVS (triple mutant, GRH) was identified in the EPSPS gene of the GRH. The nucleotide substitutions consisted of ATA102, GTC103 and TCA106 instead of ACA102, GCG103, and CCA106, respectively. The hydrogen bond distances between Gly-101 and Arg-105 positions increased from 2.89 Å (wild type) to 2.93 Å (triple-mutant) according to the EPSPS structural modeling. These results support that the high level of glyphosate resistance of the GRH A. hybridus population was mainly governed by the triple mutation TAP-IVS found of the EPSPS target site, but the impaired translocation of herbicide also contributed in this resistance.

Using RNA-seq to characterize responses to 4-hydroxyphenylpyruvate dioxygenase (HPPD) inhibitor herbicide resistance in waterhemp (Amaranthus tuberculatus)

BACKGROUND

Waterhemp (Amaranthus tuberculatus (Moq.) J.D. Sauer) is a problem weed commonly found in the Midwestern United States that can cause crippling yield losses for both maize (Zea mays L.) and soybean (Glycine max L. Merr). In 2011, 4-hydroxyphenylpyruvate-dioxygenase (HPPD, EC 1.13.11.27) inhibitor herbicide resistance was first reported in two waterhemp populations. Since the discovery of HPPD-herbicide resistance, studies have identified the mechanism of resistance and described the inheritance of the herbicide resistance. However, no studies have examined genome-wide gene expression changes in response to herbicide treatment in herbicide resistant and susceptible waterhemp.

RESULTS

We conducted RNA-sequencing (RNA-seq) analyses of two waterhemp populations (HPPD-herbicide resistant and susceptible), from herbicide-treated and mock-treated leaf samples at three, six, twelve, and twenty-four hours after treatment (HAT). We performed a de novo transcriptome assembly using all sample sequences. Following assessments of our assembly, individual samples were mapped to the de novo transcriptome allowing us to identify transcripts specific to a genotype, herbicide treatment, or time point. Our results indicate that the response of HPPD-herbicide resistant and susceptible waterhemp genotypes to HPPD-inhibiting herbicide is rapid, established as soon as 3 hours after herbicide treatment. Further, there was little overlap in gene expression between resistant and susceptible genotypes, highlighting dynamic differences in response to herbicide treatment. In addition, we used stringent analytical methods to identify candidate single nucleotide polymorphisms (SNPs) that distinguish the resistant and susceptible genotypes.

CONCLUSIONS

The waterhemp transcriptome, herbicide-responsive genes, and SNPs generated in this study provide valuable tools for future studies by numerous plant science communities. This collection of resources is essential to study and understand herbicide effects on gene expression in resistant and susceptible weeds. Understanding how herbicides impact gene expression could allow us to develop novel approaches for future herbicide development. Additionally, an increased understanding of the prolific traits intrinsic in weed success could lead to crop improvement.

Sequencing the Plastid Genome of Giant Ragweed (Ambrosia trifida, Asteraceae) From a Herbarium Specimen

We report the first plastome sequence of giant ragweed (Ambrosia trifida); with this new genome information, we assessed the phylogeny of Asteraceae and the transcriptional profiling against glyphosate resistance in giant ragweed. Assembly and genic features show a normal angiosperm quadripartite plastome structure with no signatures of deviation in gene directionality. Comparative analysis revealed large inversions across the plastome of giant ragweed and the previously sequenced members of the plant family. Asteraceae plastid genomes contain two inversions of 22.8 and 3.3 kb; the former is located between trnS-GCU and trnG-UCC genes, and the latter between trnE-UUC and trnT-GGU genes. The plastid genome sequences of A. trifida and the related species, Ambrosia artemisiifolia, are identical in gene content and arrangement, but they differ in length. The phylogeny is well-resolved and congruent with previous hypotheses about the phylogenetic relationship of Asteraceae. Transcriptomic analysis revealed divergence in the relative expressions at the exonic and intronic levels, providing hints toward the ecological adaptation of the genus. Giant ragweed shows various levels of glyphosate resistance, with introns displaying higher expression patterns at resistant time points after the assumed herbicide treatment.

Non-target site-based resistance to tribenuron-methyl and essential involved genes in Myosoton aquaticum (L.)

BACKGROUND

Water chickweed (Myosoton aquaticum (L.)) is a dicot broadleaf weed that is widespread in winter fields in China, and has evolved serious resistance to acetolactate synthase (ALS) inhibiting herbicides.

RESULTS

We identified a M. aquaticum population exhibiting moderate (6.15-fold) resistance to tribenuron-methyl (TM). Target-site ALS gene sequencing revealed no known resistance mutations in these plants, and the in vitro ALS activity assays showed no differences in enzyme sensitivity between susceptible and resistant populations; however, resistance was reversed by pretreatment with the cytochrome P450 (CYP) monooxygenase inhibitor malathion. An RNA sequencing transcriptome analysis was performed to identify candidate genes involved in metabolic resistance, and the unigenes obtained by de novo transcriptome assembly were annotated across seven databases. In total, 34 differentially expressed genes selected by digital gene expression analysis were validated by quantitative real-time (qRT)-PCR. Ten consistently overexpressed contigs, including four for CYP, four for ATP-binding cassette (ABC) transporter, and two for peroxidase were further validated by qRT-PCR using additional plants from resistant and susceptible populations. Three CYP genes (with homology to CYP734A1, CYP76C1, and CYP86B1) and one ABC transporter gene (with homology to ABCC10) were highly expressed in all resistant plants.

CONCLUSION

The mechanism of TM resistance in M. aquaticum is controlled by NTSR rather than TSR. Four genes, CYP734A1, CYP76C1, CYP86B1, and ABCC10 could play essential role in metabolic resistance to TM and justify further functional studies. To our knowledge, this is the first large-scale transcriptome analysis of genes associated with NTSR in M. aquaticum using the Illumina platform. Our data provide resource for M. aquaticum biology, and will facilitate the study of herbicide resistance mechanism at the molecular level in this species as well as in other weeds.

Reference genes for RT-qPCR normalisation in different tissues, developmental stages and stress conditions of amaranth

Studies of gene expression are very important for the identification of genes that participate in different biological processes. Currently, reverse transcription quantitative real-time PCR (RT-qPCR) is a high-throughput, sensitive and widely used method for gene expression analysis. Nevertheless, RT-qPCR requires precise normalisation of data to avoid the misinterpretation of experimental data. In this sense, the selection of reference genes is critical for gene expression analysis. At this time, several studies focus on the selection of reference genes in several species. However, the identification and validation of reference genes for the normalisation of RT-qPCR have not been described in amaranth. A set of seven housekeeping genes were analysed using RT-qPCR, to determine the most stable reference genes in amaranth for normalisation of gene expression analysis. Transcript stability and gene expression level of candidate reference genes were analysed in different tissues, at different developmental stages and under different types of stress. The data were compared using the geNorm, NormFinder and Bestkeeper statistical methods. The reference genes optimum for normalisation of data varied with respect to treatment. The results indicate that AhyMDH, AhyGAPDH, AhyEF-1α and AhyACT would be optimum for accurate normalisation of experimental data, when all treatment are analysed in the same experiment. This study presents the most stable reference genes for normalisation of gene expression analysis in amaranth, which will contribute significantly to future gene studies of this species.

Changes in the proteome of the problem weed blackgrass correlating with multiple-herbicide resistance

Herbicide resistance in grass weeds is now one of the greatest threats to sustainable cereal production in Northern Europe. Multiple-herbicide resistance (MHR), a poorly understood multigenic and quantitative trait, is particularly problematic as it provides tolerance to most classes of chemistries currently used for post-emergence weed control. Using a combination of transcriptomics and proteomics, the evolution of MHR in populations of the weed blackgrass (Alopecurus myosuroides) has been investigated. While over 4500 genes showed perturbation in their expression in MHR versus herbicide sensitive (HS) plants, only a small group of proteins showed >2-fold changes in abundance, with a mere eight proteins consistently associated with this class of resistance. Of the eight, orthologues of three of these proteins are also known to be associated with multiple drug resistance (MDR) in humans, suggesting a cross-phyla conservation in evolved tolerance to chemical agents. Proteomics revealed that MHR could be classified into three sub-types based on the association with resistance to herbicides with differing modes of action (MoA), being either global, specific to diverse chemistries acting on one MoA, or herbicide specific. Furthermore, the proteome of MHR plants were distinct from that of HS plants exposed to a range of biotic (insect feeding, plant-microbe interaction) and abiotic (N-limitation, osmotic, heat, herbicide safening) challenges commonly encountered in the field. It was concluded that MHR in blackgrass is a uniquely evolving trait(s), associated with changes in the proteome that are distinct from responses to conventional plant stresses, but sharing common features with MDR in humans.

Biochemical characterization of metabolism-based atrazine resistance in Amaranthus tuberculatus and identification of an expressed GST associated with resistance

Rapid detoxification of atrazine in naturally tolerant crops such as maize (Zea mays) and grain sorghum (Sorghum bicolor) results from glutathione S-transferase (GST) activity. In previous research, two atrazine-resistant waterhemp (Amaranthus tuberculatus) populations from Illinois, U.S.A. (designated ACR and MCR), displayed rapid formation of atrazine-glutathione (GSH) conjugates, implicating elevated rates of metabolism as the resistance mechanism. Our main objective was to utilize protein purification combined with qualitative proteomics to investigate the hypothesis that enhanced atrazine detoxification, catalysed by distinct GSTs, confers resistance in ACR and MCR. Additionally, candidate AtuGST expression was analysed in an F₂ population segregating for atrazine resistance. ACR and MCR showed higher specific activities towards atrazine in partially purified ammonium sulphate and GSH affinity-purified fractions compared to an atrazine-sensitive population (WCS). One-dimensional electrophoresis of these fractions displayed an approximate 26-kDa band, typical of GST subunits. Several phi- and tau-class GSTs were identified by LC-MS/MS from each population, based on peptide similarity with GSTs from Arabidopsis. Elevated constitutive expression of one phi-class GST, named AtuGSTF2, correlated strongly with atrazine resistance in ACR and MCR and segregating F₂ population. These results indicate that AtuGSTF2 may be linked to a metabolic mechanism that confers atrazine resistance in ACR and MCR.

Single-molecule sequencing and Hi-C-based proximity-guided assembly of amaranth (Amaranthus hypochondriacus) chromosomes provide insights into genome evolution

Background

Amaranth (Amaranthus hypochondriacus) was a food staple among the ancient civilizations of Central and South America that has recently received increased attention due to the high nutritional value of the seeds, with the potential to help alleviate malnutrition and food security concerns, particularly in arid and semiarid regions of the developing world. Here, we present a reference-quality assembly of the amaranth genome which will assist the agronomic development of the species.

Results

Utilizing single-molecule, real-time sequencing (Pacific Biosciences) and chromatin interaction mapping (Hi-C) to close assembly gaps and scaffold contigs, respectively, we improved our previously reported Illumina-based assembly to produce a chromosome-scale assembly with a scaffold N50 of 24.4 Mb. The 16 largest scaffolds contain 98% of the assembly and likely represent the haploid chromosomes (n = 16). To demonstrate the accuracy and utility of this approach, we produced physical and genetic maps and identified candidate genes for the betalain pigmentation pathway. The chromosome-scale assembly facilitated a genome-wide syntenic comparison of amaranth with other Amaranthaceae species, revealing chromosome loss and fusion events in amaranth that explain the reduction from the ancestral haploid chromosome number (n = 18) for a tetraploid member of the Amaranthaceae.

Conclusions

The assembly method reported here minimizes cost by relying primarily on short-read technology and is one of the first reported uses of in vivo Hi-C for assembly of a plant genome. Our analyses implicate chromosome loss and fusion as major evolutionary events in the 2n = 32 amaranths and clearly establish the homoeologous relationship among most of the subgenome chromosomes, which will facilitate future investigations of intragenomic changes that occurred post polyploidization.

Electronic supplementary material

The online version of this article (doi:10.1186/s12915-017-0412-4) contains supplementary material, which is available to authorized users.

Selection of relatively exact reference genes for gene expression studies in goosegrass (Eleusine indica) under herbicide stress

Goosegrass (Eleusine indica) is one of the most serious annual grassy weeds worldwide, and its evolved herbicide-resistant populations are more difficult to control. Quantitative real-time PCR (qPCR) is a common technique for investigating the resistance mechanism; however, there is as yet no report on the systematic selection of stable reference genes for goosegrass. This study proposed to test the expression stability of 9 candidate reference genes in goosegrass in different tissues and developmental stages and under stress from three types of herbicide. The results show that for different developmental stages and organs (control), eukaryotic initiation factor 4 A (eIF-4) is the most stable reference gene. Chloroplast acetolactate synthase (ALS) is the most stable reference gene under glyphosate stress. Under glufosinate stress, eIF-4 is the best reference gene. Ubiquitin-conjugating enzyme (UCE) is the most stable reference gene under quizalofop-p-ethyl stress. The gene eIF-4 is the recommended reference gene for goosegrass under the stress of all three herbicides. Moreover, pairwise analysis showed that seven reference genes were sufficient to normalize the gene expression data under three herbicides treatment. This study provides a list of reliable reference genes for transcript normalization in goosegrass, which will facilitate resistance mechanism studies in this weed species.

Transcriptional Responses of Chilean Quinoa (Chenopodium quinoa Willd.) Under Water Deficit Conditions Uncovers ABA-Independent Expression Patterns

HIGHLIGHTS R49 genotype displayed best performance on selected physiological parameters and highest tolerance to drought.R49 drought over-represented transcripts has exhibited 19% of genes (306 contigs) that presented no homology to published databases.Expression pattern for canonical responses to drought such as ABA biosynthesis and other genes induced in response to drought were assessed by qPCR. Global freshwater shortage is one of the biggest challenges of our time, often associated to misuse, increased consumption demands and the effects of climate change, paralleled with the desertification of vast areas. Chenopodium quinoa (Willd.) represents a very promising species, due to both nutritional content and cultivation under water constraint. We characterized drought tolerance of three Chilean genotypes and selected Genotype R49 (Salares ecotype) based upon Relative Water Content (RWC), Electrolyte Leakage (EL) and maximum efficiency of photosystem II (Fv/Fm) after drought treatment, when compared to another two genotypes. Exploratory RNA-Seq of R49 was generated by Illumina paired-ends method comparing drought and control irrigation conditions. We obtained 104.8 million reads, with 54 million reads for control condition and 51 million reads for drought condition. Reads were assembled in 150,952 contigs, were 31,523 contigs have a reading frame of at least 300 nucleotides (100 aminoacids). BLAST2GO annotation showed a 15% of genes without homology to NCBI proteins, but increased to 19% (306 contigs) when focused into drought-induced genes. Expression pattern for canonical drought responses such as ABA biosynthesis and other genes induced were assessed by qPCR, suggesting novelty of R49 drought responses.

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.

Investigating the mechanisms of glyphosate resistance in goosegrass (Eleusine indica (L.) Gaertn.) by RNA sequencing technology

Glyphosate is an important non-selective herbicide that is in common use worldwide. However, evolved glyphosate-resistant (GR) weeds significantly affect crop yields. Unfortunately, the mechanisms underlying resistance in GR weeds, such as goosegrass (Eleusine indica (L.) Gaertn.), an annual weed found worldwide, have not been fully elucidated. In this study, transcriptome analysis was conducted to further assess the potential mechanisms of glyphosate resistance in goosegrass. The RNA sequencing libraries generated 24 597 462 clean reads. De novo assembly analysis produced 48 852 UniGenes with an average length of 847 bp. All UniGenes were annotated using seven databases. Sixteen candidate differentially expressed genes selected by digital gene expression analysis were validated by quantitative real-time PCR (qRT-PCR). Among these UniGenes, the EPSPS and PFK genes were constitutively up-regulated in resistant (R) individuals and showed a higher copy number than that in susceptible (S) individuals. The expressions of four UniGenes relevant to photosynthesis were inhibited by glyphosate in S individuals, and this toxic response was confirmed by gas exchange analysis. Two UniGenes annotated as glutathione transferase (GST) were constitutively up-regulated in R individuals, and were induced by glyphosate both in R and S. In addition, the GST activities in R individuals were higher than in S. Our research confirmed that two UniGenes (PFK, EPSPS) were strongly associated with target resistance, and two GST-annotated UniGenes may play a role in metabolic glyphosate resistance in goosegrass.

Transcriptome Profiling to Identify Genes Involved in Mesosulfuron-Methyl Resistance in Alopecurus aequalis

Non-target-site resistance (NTSR) to herbicides is a worldwide concern for weed control. However, as the dominant NTSR mechanism in weeds, metabolic resistance is not yet well-characterized at the genetic level. For this study, we have identified a shortawn foxtail (Alopecurus aequalis Sobol.) population displaying both TSR and NTSR to mesosulfuron-methyl and fenoxaprop-P-ethyl, yet the molecular basis for this NTSR remains unclear. To investigate the mechanisms of metabolic resistance, an RNA-Seq transcriptome analysis was used to find candidate genes that may confer metabolic resistance to the herbicide mesosulfuron-methyl in this plant population. The RNA-Seq libraries generated 831,846,736 clean reads. The de novo transcriptome assembly yielded 95,479 unigenes (averaging 944 bp in length) that were assigned putative annotations. Among these, a total of 29,889 unigenes were assigned to 67 GO terms that contained three main categories, and 14,246 unigenes assigned to 32 predicted KEGG metabolic pathways. Global gene expression was measured using the reads generated from the untreated control (CK), water-only control (WCK), and mesosulfuron-methyl treatment (T) of R and susceptible (S). Contigs that showed expression differences between mesosulfuron-methyl-treated R and S biotypes, and between mesosulfuron-methyl-treated, water-treated and untreated R plants were selected for further quantitative real-time PCR (qRT-PCR) validation analyses. Seventeen contigs were consistently highly expressed in the resistant A. aequalis plants, including four cytochrome P450 monooxygenase (CytP450) genes, two glutathione S-transferase (GST) genes, two glucosyltransferase (GT) genes, two ATP-binding cassette (ABC) transporter genes, and seven additional contigs with functional annotations related to oxidation, hydrolysis, and plant stress physiology. These 17 contigs could serve as major candidate genes for contributing to metabolic mesosulfuron-methyl resistance; hence they deserve further functional study. This is the first large-scale transcriptome-sequencing study to identify NTSR genes in A. aequalis that uses the Illumina platform. This work demonstrates that NTSR is likely driven by the differences in the expression patterns of a set of genes. The assembled transcriptome data presented here provide a valuable resource for A. aequalis biology, and should facilitate the study of herbicide resistance at the molecular level in this and other weed species.

454 Pyrosequencing of Olive (Olea europaea L.) Transcriptome in Response to Salinity

Olive (Olea europaea L.) is one of the most important crops in the Mediterranean region. The expansion of cultivation in areas irrigated with low quality and saline water has negative effects on growth and productivity however the investigation of the molecular basis of salt tolerance in olive trees has been only recently initiated. To this end, we investigated the molecular response of cultivar Kalamon to salinity stress using next-generation sequencing technology to explore the transcriptome profile of olive leaves and roots and identify differentially expressed genes that are related to salt tolerance response. Out of 291,958 obtained trimmed reads, 28,270 unique transcripts were identified of which 35% are annotated, a percentage that is comparable to similar reports on non-model plants. Among the 1,624 clusters in roots that comprise more than one read, 24 were differentially expressed comprising 9 down- and 15 up-regulated genes. Respectively, inleaves, among the 2,642 clusters, 70 were identified as differentially expressed, with 14 down- and 56 up-regulated genes. Using next-generation sequencing technology we were able to identify salt-response-related transcripts. Furthermore we provide an annotated transcriptome of olive as well as expression data, which are both significant tools for further molecular studies in olive.

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.

RNA-Seq transcriptome analysis of maize inbred carrying nicosulfuron-tolerant and nicosulfuron-susceptible alleles

Postemergence applications of nicosulfuron can cause great damage to certain maize inbred lines and hybrids. Variation among different responses to nicosulfuron may be attributed to differential rates of herbicide metabolism. We employed RNA-Seq analysis to compare transcriptome responses between nicosulfuron-treated and untreated in both tolerant and susceptible maize plants. A total of 71.8 million paired end Illumina RNA-Seq reads were generated, representing the transcription of around 40,441 unique reads. About 345,171 gene ontology (GO) term assignments were conducted for the annotation in terms of biological process, cellular component and molecular function categories, and 6413 sequences with 108 enzyme commission numbers were assigned to 134 predicted Kyoto Encyclopedia of Genes and Genomes (KEGG) metabolic pathways. Digital gene expression profile (DGE) analysis using Solexa sequencing was performed within the susceptible and tolerant maize between the nicosulfuron-treated and untreated conditions, 13 genes were selected as the candidates most likely involved in herbicide metabolism, and quantitative RT-PCR validated the RNA-Seq results for eight genes. This transcriptome data may provide opportunities for the study of sulfonylurea herbicides susceptibility emergence of Zea mays.

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.

Glyphosate resistance: state of knowledge

Studies of mechanisms of resistance to glyphosate have increased current understanding of herbicide resistance mechanisms. Thus far, single-codon non-synonymous mutations of EPSPS (5-enolypyruvylshikimate-3-phosphate synthase) have been rare and, relative to other herbicide mode of action target-site mutations, unconventionally weak in magnitude for resistance to glyphosate. However, it is possible that weeds will emerge with non-synonymous mutations of two codons of EPSPS to produce an enzyme endowing greater resistance to glyphosate. Today, target-gene duplication is a common glyphosate resistance mechanism and could become a fundamental process for developing any resistance trait. Based on competition and substrate selectivity studies in several species, rapid vacuole sequestration of glyphosate occurs via a transporter mechanism. Conversely, as the chloroplast requires transporters for uptake of important metabolites, transporters associated with the two plastid membranes may separately, or together, successfully block glyphosate delivery. A model based on finite glyphosate dose and limiting time required for chloroplast loading sets the stage for understanding how uniquely different mechanisms can contribute to overall glyphosate resistance.

The draft genome and transcriptome of Amaranthus hypochondriacus: a C4 dicot producing high-lysine edible pseudo-cereal

Grain amaranths, edible C4 dicots, produce pseudo-cereals high in lysine. Lysine being one of the most limiting essential amino acids in cereals and C4 photosynthesis being one of the most sought-after phenotypes in protein-rich legume crops, the genome of one of the grain amaranths is likely to play a critical role in crop research. We have sequenced the genome and transcriptome of Amaranthus hypochondriacus, a diploid (2n = 32) belonging to the order Caryophyllales with an estimated genome size of 466 Mb. Of the 411 linkage single-nucleotide polymorphisms (SNPs) reported for grain amaranths, 355 SNPs (86%) are represented in the scaffolds and 74% of the 8.6 billion bases of the sequenced transcriptome map to the genomic scaffolds. The genome of A. hypochondriacus, codes for at least 24,829 proteins, shares the paleohexaploidy event with species under the superorders Rosids and Asterids, harbours 1 SNP in 1,000 bases, and contains 13.76% of repeat elements. Annotation of all the genes in the lysine biosynthetic pathway using comparative genomics and expression analysis offers insights into the high-lysine phenotype. As the first grain species under Caryophyllales and the first C4 dicot genome reported, the work presented here will be beneficial in improving crops and in expanding our understanding of angiosperm evolution.

Transcriptome analysis by Illumina high-throughout paired-end sequencing reveals the complexity of differential gene expression during in vitro plantlet growth and flowering in Amaranthus tricolor L

Amaranthus tricolor L. is a C4 plant, which is consumed as a major leafy vegetable in some tropical countries. Under conditions of high temperature and short daylight, Am. tricolor readily bolts and blooms, degrading leaf quality. A preliminary in vitro flowering study demonstrated that the flowering control pathway in Am. tricolor may differ from that of Arabidopsis. Nevertheless, no transcriptome analysis of the flowering process in Amaranthus has been conducted. To study Am. tricolor floral regulatory mechanisms, we conducted a large-scale transcriptome analysis--based on Illumina HiSeq sequencing of cDNA libraries generated from Am. tricolor at young seedling (YSS), adult seedling (ASS), flower bud (FBS), and flowering (FS) stages. A total of 99,312 unigenes were obtained. Using BLASTX, 43,088 unigenes (43.39%) were found to have significant similarity with accessions in Nr, Nt, and Swiss-Prot databases. Of these unigenes, 11,291 were mapped to 266 KEGG pathways. Further analysis of the four digital transcriptomes revealed that 735, 17,184, 274, and 206 unigenes were specifically expressed during YSS, ASS, FBS, and FS, respectively, with 59,517 unigenes expressed throughout the four stages. These unigenes were involved in many metabolic pathways related to in vitro flowering. Among these pathways, 259 unigenes were associated with ubiquitin-mediated proteolysis, indicating its importance for in vitro flowering in Am. tricolor. Other pathways, such as circadian rhythm and cell cycle, also had important roles. Finally, 26 unigenes were validated by qRT-PCR in samples from Am. tricolor at YSS, ASS, FBS, and FS; their differential expressions at the various stages indicate their possible roles in Am. tricolor growth and development, but the results were somewhat similar to Arabidopsis. Because unigenes involved in many metabolic pathways or of unknown function were revealed to regulate in vitro plantlet growth and flowering in Am. tricolor, the process appears to be highly complex in this species.

RNA-Seq transcriptome analysis to identify genes involved in metabolism-based diclofop resistance in Lolium rigidum

Weed control failures due to herbicide resistance are an increasing and worldwide problem that significantly affect crop yields. Metabolism-based herbicide resistance (referred to as metabolic resistance) in weeds is not well characterized at the genetic level. An RNA-Seq transcriptome analysis was used to find candidate genes that conferred metabolic resistance to the herbicide diclofop in a diclofop-resistant population (R) of the major global weed Lolium rigidum. A reference cDNA transcriptome (19 623 contigs) was assembled and assigned putative annotations. Global gene expression was measured using Illumina reads from untreated control, adjuvant-only control, and diclofop treatment of R and susceptible (S). Contigs that showed constitutive expression differences between untreated R and untreated S were selected for further validation analysis, including 11 contigs putatively annotated as cytochrome P450 (CytP450), glutathione transferase (GST), or glucosyltransferase (GT), and 17 additional contigs with annotations related to metabolism or signal transduction. In a forward genetics validation experiment, nine contigs had constitutive up-regulation in R individuals from a segregating F2 population, including three CytP450, one nitronate monooxygenase (NMO), three GST, and one GT. Principal component analysis using these nine contigs differentiated F2 -R from F2 -S individuals. In a physiological validation experiment in which 2,4-D pre-treatment induced diclofop protection in S individuals due to increased metabolism, seven of the nine genetically validated contigs were induced significantly. Four contigs (two CytP450, NMO, and GT) were consistently highly expressed in nine field-evolved metabolic resistant L. rigidum populations. These four contigs were strongly associated with the resistance phenotype and are major candidates for contributing to metabolic diclofop resistance.

Distinct detoxification mechanisms confer resistance to mesotrione and atrazine in a population of waterhemp

Previous research reported the first case of resistance to mesotrione and other 4-hydroxyphenylpyruvate dioxygenase (HPPD) herbicides in a waterhemp (Amaranthus tuberculatus) population designated MCR (for McLean County mesotrione- and atrazine-resistant). Herein, experiments were conducted to determine if target site or nontarget site mechanisms confer mesotrione resistance in MCR. Additionally, the basis for atrazine resistance was investigated in MCR and an atrazine-resistant but mesotrione-sensitive population (ACR for Adams County mesotrione-sensitive but atrazine-resistant). A standard sensitive population (WCS for Wayne County herbicide-sensitive) was also used for comparison. Mesotrione resistance was not due to an alteration in HPPD sequence, HPPD expression, or reduced herbicide absorption. Metabolism studies using whole plants and excised leaves revealed that the time for 50% of absorbed mesotrione to degrade in MCR was significantly shorter than in ACR and WCS, which correlated with previous phenotypic responses to mesotrione and the quantity of the metabolite 4-hydroxy-mesotrione in excised leaves. The cytochrome P450 monooxygenase inhibitors malathion and tetcyclacis significantly reduced mesotrione metabolism in MCR and corn (Zea mays) excised leaves but not in ACR. Furthermore, malathion increased mesotrione activity in MCR seedlings in greenhouse studies. These results indicate that enhanced oxidative metabolism contributes significantly to mesotrione resistance in MCR. Sequence analysis of atrazine-resistant (MCR and ACR) and atrazine-sensitive (WCS) waterhemp populations detected no differences in the psbA gene. The times for 50% of absorbed atrazine to degrade in corn, MCR, and ACR leaves were shorter than in WCS, and a polar metabolite of atrazine was detected in corn, MCR, and ACR that cochromatographed with a synthetic atrazine-glutathione conjugate. Thus, elevated rates of metabolism via distinct detoxification mechanisms contribute to mesotrione and atrazine resistance within the MCR population.

A novel glucagon-like peptide 1 peptide identified from Ophisaurus harti

Glucagon-like peptide 1 receptor (GLP1R) is a promising target for the treatment of type 2 diabetes. Because of the short half-life of endogenous GLP1 peptide, other GLP1R agonists are considered to be appealing therapeutic candidates. A high-throughput assay has been established to screen for GLP1R agonists in a 60 000-well natural product compound library fractionated from 670 different herbs/materials widely used in traditional Chinese medicines (TCMs). The screening is based on primary screen of GLP1R⁺ reporter gene assay with the counter screen in GLP1R⁻ cell line. An active fraction, A089-147, was identified from the screening. Fraction A089-147 was isolated from dried Ophisaurus harti, and the fact that its GLP1R agonist activity was sensitive to trypsin treatment indicates its peptidic nature. The active ingredient of A089-147 was later identified as O. harti GLP1 through transcriptome analysis. Chemically synthesized O. harti GLP1 showed GLP1R agonist activity and sensitivity to dipeptidase IV digestion. This study illustrated a comprehensive screening strategy to identify novel GLP1R agonists from TCMs libraries and at the same time underlined the difficulty of identifying a non-peptidic GLP1R agonist. The novel O. harti GLP1 peptide yielded from this study confirmed broader application of TCMs libraries in active peptide identification.

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.

Unravelling the genetic bases of non-target-site-based resistance (NTSR) to herbicides: a major challenge for weed science in the forthcoming decade

Non-target-site-based resistance (NTSR) can confer unpredictable cross-resistance to herbicides. However, the genetic determinants of NTSR remain poorly known. The current, urgent challenge for weed scientists is thus to elucidate the bases of NTSR so that detection tools are developed, the evolution of NTSR is understood, the efficacy of the shrinking herbicide portfolio is maintained and integrated weed management strategies, including fully effective herbicide applications, are designed and implemented. In this paper, the importance of NTSR in resistance to herbicides is underlined. The most likely way in which NTSR evolves-by accumulation of different mechanisms within individual plants-is described. The NTSR mechanisms, which can interfere with herbicide penetration, translocation and accumulation at the target site, and/or protect the plant against the consequences of herbicide action, are then reviewed. NTSR is a part of the plant stress response. As such, NTSR is a dynamic process unrolling over time that involves 'protectors' directly interfering with herbicide action, and also regulators controlling 'protector' expression. NTSR is thus a quantitative trait. On this basis, a three-step procedure is proposed, based on the use of the 'omics' (genomics, transcriptomics, proteomics or metabolomics), to unravel the genetic bases of NTSR.

A catalogue of putative unique transcripts from Douglas-fir (Pseudotsuga menziesii) based on 454 transcriptome sequencing of genetically diverse, drought stressed seedlings

Background

Douglas-fir (Pseudotsuga menziesii) extends over a wide range of contrasting environmental conditions, reflecting substantial local adaptation. For this reason, it is an interesting model species to study plant adaptation and the effects of global climate change such as increased temperatures and significant periods of drought on individual trees and the forest landscape in general. However, genomic data and tools for studying genetic variation in natural populations to understand the genetic and physiological mechanisms of adaptation are currently missing for Douglas-fir. This study represents a first step towards characterizing the Douglas-fir transcriptome based on 454 sequencing of twelve cDNA libraries. The libraries were constructed from needle and wood tissue of coastal and interior provenances subjected to drought stress experiments.

Results

The 454 sequencing of twelve normalized cDNA libraries resulted in 3.6 million reads from which a set of 170,859 putative unique transcripts (PUTs) was assembled. Functional annotation by BLAST searches and Gene Ontology mapping showed that the composition of functional classes is very similar to other plant transcriptomes and demonstrated that a large fraction of the Douglas-fir transcriptome is tagged by the PUTs. Based on evolutionary conservation, we identified about 1,000 candidate genes related to drought stress. A total number of 187,653 single nucleotide polymorphisms (SNPs) were detected by three SNP detection tools. However, only 27,688 SNPs were identified by all three methods, indicating that SNP detection depends on the particular method used. The two alleles of about 60% of the 27,688 SNPs are segregating simultaneously in both coastal and interior provenances, which indicates a high proportion of ancestral shared polymorphisms or a high level of gene flow between these two ecologically and phenotypically different varieties.

Conclusions

We established a catalogue of PUTs and large SNP database for Douglas-fir. Both will serve as a useful resource for the further characterization of the genome and transcriptome of Douglas-fir and for the analysis of genetic variation using genotyping or resequencing methods.

De novo sequencing and analysis of the Ulva linza transcriptome to discover putative mechanisms associated with its successful colonization of coastal ecosystems

BACKGROUND

The green algal genus Ulva Linnaeus (Ulvaceae, Ulvales, Chlorophyta) is well known for its wide distribution in marine, freshwater, and brackish environments throughout the world. The Ulva species are also highly tolerant of variations in salinity, temperature, and irradiance and are the main cause of green tides, which can have deleterious ecological effects. However, limited genomic information is currently available in this non-model and ecologically important species. Ulva linza is a species that inhabits bedrock in the mid to low intertidal zone, and it is a major contributor to biofouling. Here, we presented the global characterization of the U. linza transcriptome using the Roche GS FLX Titanium platform, with the aim of uncovering the genomic mechanisms underlying rapid and successful colonization of the coastal ecosystems.

RESULTS

De novo assembly of 382,884 reads generated 13,426 contigs with an average length of 1,000 bases. Contiguous sequences were further assembled into 10,784 isotigs with an average length of 1,515 bases. A total of 304,101 reads were nominally identified by BLAST; 4,368 isotigs were functionally annotated with 13,550 GO terms, and 2,404 isotigs having enzyme commission (EC) numbers were assigned to 262 KEGG pathways. When compared with four other full sequenced green algae, 3,457 unique isotigs were found in U. linza and 18 conserved in land plants. In addition, a specific photoprotective mechanism based on both LhcSR and PsbS proteins and a C4-like carbon-concentrating mechanism were found, which may help U. linza survive stress conditions. At least 19 transporters for essential inorganic nutrients (i.e., nitrogen, phosphorus, and sulphur) were responsible for its ability to take up inorganic nutrients, and at least 25 eukaryotic cytochrome P450s, which is a higher number than that found in other algae, may be related to their strong allelopathy. Multi-origination of the stress related proteins, such as glutamate dehydrogenase, superoxide dismutases, ascorbate peroxidase, catalase and heat-shock proteins, may also contribute to colonization of U. linza under stress conditions.

CONCLUSIONS

The transcriptome of U. linza uncovers some potential genomic mechanisms that might explain its ability to rapidly and successfully colonize coastal ecosystems, including the land-specific genes; special photoprotective mechanism based on both LhcSR and PsbS; development of C4-like carbon-concentrating mechanisms; muti-origin transporters for essential inorganic nutrients; multiple and complex P450s; and glutamate dehydrogenase, superoxide dismutases, ascorbate peroxidase, catalase, and heat-shock proteins that are related to stress resistance.

Evaluating characteristics of de novo assembly software on 454 transcriptome data: a simulation approach

Background

The quantity of transcriptome data is rapidly increasing for non-model organisms. As sequencing technology advances, focus shifts towards solving bioinformatic challenges, of which sequence read assembly is the first task. Recent studies have compared the performance of different software to establish a best practice for transcriptome assembly. Here, we adapted a simulation approach to evaluate specific features of assembly programs on 454 data. The novelty of our study is that the simulation allows us to calculate a model assembly as reference point for comparison.

Findings

The simulation approach allows us to compare basic metrics of assemblies computed by different software applications (CAP3, MIRA, Newbler, and Oases) to a known optimal solution. We found MIRA and CAP3 are conservative in merging reads. This resulted in comparably high number of short contigs. In contrast, Newbler more readily merged reads into longer contigs, while Oases produced the overall shortest assembly. Due to the simulation approach, reads could be traced back to their correct placement within the transcriptome. Together with mapping reads onto the assembled contigs, we were able to evaluate ambiguity in the assemblies. This analysis further supported the conservative nature of MIRA and CAP3, which resulted in low proportions of chimeric contigs, but high redundancy. Newbler produced less redundancy, but the proportion of chimeric contigs was higher.

Conclusion

Our evaluation of four assemblers suggested that MIRA and Newbler slightly outperformed the other programs, while showing contrasting characteristics. Oases did not perform very well on the 454 reads. Our evaluation indicated that the software was either conservative (MIRA) or liberal (Newbler) about merging reads into contigs. This suggested that in choosing an assembly program researchers should carefully consider their follow up analysis and consequences of the chosen approach to gain an assembly.

[Application of next generation sequencing techniques in plant transcriptome]

With the development of DNA sequencing techniques, the next-generation sequencing (NGS) techniques with the characteristics of high-throughput and low cost have become the first choice for more and more researchers to carry out the biological researches. Among the next-generation sequencing techniques, the 454 sequencing platform is the first commercially available and relatively mature one and widely used in various fields of biological research. Taking 454 sequencing platform as an example, we illustrate the advantages and disadvantages of NGS technical principles, review their applications in plant transcriptome, and outlook their future development and applications in plant research field.

Designing a transcriptome next-generation sequencing project for a nonmodel plant species

The application of next-generation sequencing (NGS) to transcriptomics, commonly called RNA-seq, allows the nearly complete characterization of transcriptomic events occurring in a specific tissue. It has proven particularly useful in nonmodel species, which often lack the resources available for sequenced organisms. Mainly, RNA-seq does not require a reference genome to gain useful transcriptomic information. In this review, the application of RNA-seq to nonmodel plant species will be addressed. Important experimental considerations from presequencing issues to postsequencing analysis, including sample and platform selection, and useful bioinformatics tools for assembly and data analysis, are covered. Methods of assembling RNA-seq data and analyses commonly performed with RNA-seq data, including single nucleotide polymorphism detection and analysis of differential expression, are explored. In addition, studies that have used RNA-seq to elucidate nonmodel plant transcriptomics are highlighted.

RNA-Seq Assembly - Are We There Yet?

Transcriptomic sequence resources represent invaluable assets for research, in particular for non-model species without a sequenced genome. To date, the Next Generation Sequencing technologies 454/Roche and Illumina have been used to generate transcriptome sequence databases by mRNA-Seq for more than fifty different plant species. While some of the databases were successfully used for downstream applications, such as proteomics, the assembly parameters indicate that the assemblies do not yet accurately reflect the actual plant transcriptomes. Two different assembly strategies have been used, overlap consensus based assemblers for long reads and Eulerian path/de Bruijn graph assembler for short reads. In this review, we discuss the challenges and solutions to the transcriptome assembly problem. A list of quality control parameters and the necessary scripts to produce them are provided.

Interspecific hybridization transfers a previously unknown glyphosate resistance mechanism in Amaranthus species

A previously unknown glyphosate resistance mechanism, amplification of the 5-enolpyruvyl shikimate-3-phosphate synthase gene, was recently reported in Amaranthus palmeri. This evolved mechanism could introgress to other weedy Amaranthus species through interspecific hybridization, representing an avenue for acquisition of a novel adaptive trait. The objective of this study was to evaluate the potential for this glyphosate resistance trait to transfer via pollen from A. palmeri to five other weedy Amaranthus species (Amaranthus hybridus, Amaranthus powellii, Amaranthus retroflexus, Amaranthus spinosus, and Amaranthus tuberculatus). Field and greenhouse crosses were conducted using glyphosate-resistant male A. palmeri as pollen donors and the other Amaranthus species as pollen recipients. Hybridization between A. palmeri and A. spinosus occurred with frequencies in the field studies ranging from <0.01% to 0.4%, and 1.4% in greenhouse crosses. A majority of the A. spinosus × A. palmeri hybrids grown to flowering were monoecious and produced viable seed. Hybridization occurred in the field study between A. palmeri and A. tuberculatus (<0.2%), and between A. palmeri and A. hybridus (<0.01%). This is the first documentation of hybridization between A. palmeri and both A. spinosus and A. hybridus.

Transcriptomic analysis of grain amaranth (Amaranthus hypochondriacus) using 454 pyrosequencing: comparison with A. tuberculatus, expression profiling in stems and in response to biotic and abiotic stress

BACKGROUND

Amaranthus hypochondriacus, a grain amaranth, is a C4 plant noted by its ability to tolerate stressful conditions and produce highly nutritious seeds. These possess an optimal amino acid balance and constitute a rich source of health-promoting peptides. Although several recent studies, mostly involving subtractive hybridization strategies, have contributed to increase the relatively low number of grain amaranth expressed sequence tags (ESTs), transcriptomic information of this species remains limited, particularly regarding tissue-specific and biotic stress-related genes. Thus, a large scale transcriptome analysis was performed to generate stem- and (a)biotic stress-responsive gene expression profiles in grain amaranth.

RESULTS

A total of 2,700,168 raw reads were obtained from six 454 pyrosequencing runs, which were assembled into 21,207 high quality sequences (20,408 isotigs + 799 contigs). The average sequence length was 1,064 bp and 930 bp for isotigs and contigs, respectively. Only 5,113 singletons were recovered after quality control. Contigs/isotigs were further incorporated into 15,667 isogroups. All unique sequences were queried against the nr, TAIR, UniRef100, UniRef50 and Amaranthaceae EST databases for annotation. Functional GO annotation was performed with all contigs/isotigs that produced significant hits with the TAIR database. Only 8,260 sequences were found to be homologous when the transcriptomes of A. tuberculatus and A. hypochondriacus were compared, most of which were associated with basic house-keeping processes. Digital expression analysis identified 1,971 differentially expressed genes in response to at least one of four stress treatments tested. These included several multiple-stress-inducible genes that could represent potential candidates for use in the engineering of stress-resistant plants. The transcriptomic data generated from pigmented stems shared similarity with findings reported in developing stems of Arabidopsis and black cottonwood (Populus trichocarpa).

CONCLUSIONS

This study represents the first large-scale transcriptomic analysis of A. hypochondriacus, considered to be a highly nutritious and stress-tolerant crop. Numerous genes were found to be induced in response to (a)biotic stress, many of which could further the understanding of the mechanisms that contribute to multiple stress-resistance in plants, a trait that has potential biotechnological applications in agriculture.

Selection and validation of reference genes for quantitative real-time PCR in buckwheat (Fagopyrum esculentum) based on transcriptome sequence data

Quantitative reverse transcription PCR (qRT-PCR) is one of the most precise and widely used methods of gene expression analysis. A necessary prerequisite of exact and reliable data is the accurate choice of reference genes. We studied the expression stability of potential reference genes in common buckwheat (Fagopyrum esculentum) in order to find the optimal reference for gene expression analysis in this economically important crop. Recently sequenced buckwheat floral transcriptome was used as source of sequence information. Expression stability of eight candidate reference genes was assessed in different plant structures (leaves and inflorescences at two stages of development and fruits). These genes are the orthologs of Arabidopsis genes identified as stable in a genome-wide survey gene of expression stability and a traditionally used housekeeping gene GAPDH. Three software applications – geNorm, NormFinder and BestKeeper - were used to estimate expression stability and provided congruent results. The orthologs of AT4G33380 (expressed protein of unknown function, Expressed1), AT2G28390 (SAND family protein, SAND) and AT5G46630 (clathrin adapter complex subunit family protein, CACS) are revealed as the most stable. We recommend using the combination of Expressed1, SAND and CACS for the normalization of gene expression data in studies on buckwheat using qRT-PCR. These genes are listed among five the most stably expressed in Arabidopsis that emphasizes utility of the studies on model plants as a framework for other species.