Transcriptome Assembly and Comparison of an Allotetraploid Weed Species, Annual Bluegrass, with its Two Diploid Progenitor Species, Schrad and Kunth

Recent allopolyploidization and transcriptomic asymmetry in the mangrove shrub Acanthus tetraploideus

BACKGROUND

Mangrove species are vital to the ecosystems of tropical and subtropical coastlines worldwide. Despite the underexplored role of polyploidization in these species, deciphering its impact on gene expression is essential for understanding the connection between polyploidization and species diversification. Our initial investigation, integrating multiple nuclear loci with morphological and cytological data, indicates that the tetraploid Acanthus tetraploideus likely originated from allopolyploidization events involving the diploid species A. ilicifolius and A. ebracteatus. Expanding on these insights, this study utilizes genome-wide evidence to confirm the divergence patterns among extant Acanthus mangrove diploids and to investigate the origin and transcriptome asymmetry of the tetraploid A. tetraploideus.

RESULTS

Phylogenetic analysis and molecular dating revealed a closer evolutionary relationship between A. ebracteatus and A. volubilis than between A. ebracteatus and A. ilicifolius, diverged approximately 6.92 Mya and 9.59 Mya, respectively. Analysis of individual whole transcriptomes revealed that homeologous sequences in A. tetraploideus were preferentially clustered with A. ilicifolius and A. ebracteatus, rather than A. volubilis, in a roughly 1:1 ratio. The high similarity in nucleotide sequences and homologous polymorphisms between the tetraploid A. tetraploideus and its two parental diploids, A. ebracteatus and A. ilicifolius, supports the hypothesis of a recent allopolyploid origin for A. tetraploideus. Estimation of homeolog expression revealed a general attenuation of homeolog expression divergence in A. tetraploideus compared to the in silico parental mix, with 22.87% and 67.66% of genes exhibiting biased homeolog expression, respectively. Further investigation identified remarkable retention of parental expression dominance in the tetraploid, suggesting that parental genetic legacy substantially influences the reconfiguration of homeolog expression in the derived tetraploid. Meanwhile, the observation of numerous novel expression patterns between the two homeolog sets suggests that the transcriptome shock (i.e., the transcriptomic changes induced by interspecific hybridization) associated with allopolyploidization and subsequent post-polyploid evolutionary processes also significantly impact transcriptome asymmetry in A. tetraploideus. While no strong evidence directly links transcriptomic changes to specific adaptive traits, the patterns in unbiased and novelly biased genes in A. tetraploideus suggest adaptations to stable polyploidy. Unbiased genes involved in fundamental cellular processes and novelly biased genes related to chromosome dynamics and cell cycle regulation may stabilize polyploid genomes, supporting the species' establishment and long-term success. These findings underscore the role of transcriptomic stability in polyploid adaptation.

CONCLUSIONS

Our study sheds light on the evolutionary origins and the intricate transcriptional reconfiguration of the tetraploid A. tetraploideus. These insights significantly enhance our comprehension of the pivotal role that polyploidization plays in speciation and adaptative evolution of mangrove species.

Homoeologous evolution of the allotetraploid genome of Poa annua L

BACKGROUND

Poa annua (annual bluegrass) is an allotetraploid turfgrass, an agronomically significant weed, and one of the most widely dispersed plant species on earth. Here, we report the chromosome-scale genome assemblies of P. annua's diploid progenitors, P. infirma and P. supina, and use multi-omic analyses spanning all three species to better understand P. annua's evolutionary novelty.

RESULTS

We find that the diploids diverged from their common ancestor 5.5 - 6.3 million years ago and hybridized to form P. annua ≤ 50,000 years ago. The diploid genomes are similar in chromosome structure and most notably distinguished by the divergent evolutionary histories of their transposable elements, leading to a 1.7 × difference in genome size. In allotetraploid P. annua, we find biased movement of retrotransposons from the larger (A) subgenome to the smaller (B) subgenome. We show that P. annua's B subgenome is preferentially accumulating genes and that its genes are more highly expressed. Whole-genome resequencing of several additional P. annua accessions revealed large-scale chromosomal rearrangements characterized by extensive TE-downsizing and evidence to support the Genome Balance Hypothesis.

CONCLUSIONS

The divergent evolutions of the diploid progenitors played a central role in conferring onto P. annua its remarkable phenotypic plasticity. We find that plant genes (guided by selection and drift) and transposable elements (mostly guided by host immunity) each respond to polyploidy in unique ways and that P. annua uses whole-genome duplication to purge highly parasitized heterochromatic sequences. The findings and genomic resources presented here will enable the development of homoeolog-specific markers for accelerated weed science and turfgrass breeding.

Expanding the Repertoire of the Plant-Infecting Ophioviruses through Metatranscriptomics Data

Ophioviruses (genus Ophiovirus, family Aspiviridae) are plant-infecting viruses with non-enveloped, filamentous, naked nucleocapsid virions. Members of the genus Ophiovirus have a segmented single-stranded negative-sense RNA genome (ca. 11.3–12.5 kb), encompassing three or four linear segments. In total, these segments encode four to seven proteins in the sense and antisense orientation, both in the viral and complementary strands. The genus Ophiovirus includes seven species with viruses infecting both monocots and dicots, mostly trees, shrubs and some ornamentals. From a genomic perspective, as of today, there are complete genomes available for only four species. Here, by exploring large publicly available metatranscriptomics datasets, we report the identification and molecular characterization of 33 novel viruses with genetic and evolutionary cues of ophioviruses. Genetic distance and evolutionary insights suggest that all the detected viruses could correspond to members of novel species, which expand the current diversity of ophioviruses ca. 4.5-fold. The detected viruses increase the tentative host range of ophioviruses for the first time to mosses, liverwort and ferns. In addition, the viruses were linked to several Asteraceae, Orchidaceae and Poaceae crops/ornamental plants. Phylogenetic analyses showed a novel clade of mosses, liverworts and fern ophioviruses, characterized by long branches, suggesting that there is still plenty of unsampled hidden diversity within the genus. This study represents a significant expansion of the genomics of ophioviruses, opening the door to future works on the molecular and evolutionary peculiarity of this virus genus.

Chromosome-Scale Genome Assembly and Annotation of Allotetraploid Annual Bluegrass (Poa annua L.)

Poa annua L. is a globally distributed grass with economic and horticultural significance as a weed and as a turfgrass. This dual significance, and its phenotypic plasticity and ecological adaptation, have made P. annua an intriguing plant for genetic and evolutionary studies. Because of the lack of genomic resources and its allotetraploid (2n = 4x = 28) nature, a reference genome sequence would be a valuable asset to better understand the significance and polyploid origin of P. annua. Here we report a genome assembly with scaffolds representing the 14 haploid chromosomes that are 1.78 Gb in length with an N50 of 112 Mb and 96.7% of BUSCO orthologs. Seventy percent of the genome was identified as repetitive elements, 91.0% of which were Copia- or Gypsy-like long-terminal repeats. The genome was annotated with 76,420 genes spanning 13.3% of the 14 chromosomes. The two subgenomes originating from Poa infirma (Knuth) and Poa supina (Schrad) were sufficiently divergent to be distinguishable but syntenic in sequence and annotation with repetitive elements contributing to the expansion of the P. infirma subgenome.

Crystal structure of Arabidopsis thaliana HPPK/DHPS, a bifunctional enzyme and target of the herbicide asulam

Herbicides are vital for modern agriculture, but their utility is threatened by genetic or metabolic resistance in weeds, as well as regulatory barriers. Of the known herbicide modes of action, 7,8-dihydropterin synthase (DHPS), which is involved in folate biosynthesis, is targeted by just one commercial herbicide, asulam. A mimic of the substrate para-aminobenzoic acid, asulam is chemically similar to sulfonamide antibiotics, and although it is still in widespread use, asulam has faced regulatory scrutiny. With an entire mode of action represented by just one commercial agrochemical, we sought to improve the understanding of its plant target. Here we solve a 2.3 Å resolution crystal structure for Arabidopsis thaliana DHPS that is conjoined to 6-hydroxymethyl-7,8-dihydropterin pyrophosphokinase (HPPK), and we reveal a strong structural conservation with bacterial counterparts at the sulfonamide-binding pocket of DHPS. We demonstrate that asulam and the antibiotic sulfamethoxazole have herbicidal as well as antibacterial activity, and we explore the structural basis of their potency by modeling these compounds in mitochondrial HPPK/DHPS. Our findings suggest limited opportunity for the rational design of plant selectivity from asulam and indicate that pharmacokinetic or delivery differences between plants and microbes might be the best ways to safeguard this mode of action.

Differential expression of genes associated with non-target site resistance in Poa annua with target site resistance to acetolactate synthase inhibitors

BACKGROUND

Poa annua is a pervasive grassy, self-pollinating, weed that has evolved resistance to 10 different herbicide modes-of-action, third most of all weed species. We investigated constitutive overexpression of genes associated with non-target site resistance (NTSR) in POAAN-R3 and the response of those genes when treated with trifloxysulfuron despite the biotype having a known target site mutation in acetolactate synthase (ALS).

RESULTS

Despite having an ALS target site mutation, POAAN-R3 still had a transcriptomic response to herbicide application that differed from a susceptible biotype. We observed differential expression of genes associated with transmembrane transport and oxidation-reduction activities, with differences being most pronounced prior to herbicide treatment.

CONCLUSIONS

In the P. annua biotype we studied with confirmed target site resistance to ALS inhibitors, we also observed constitutive expression of genes regulating transmembrane transport, as well as differential expression of genes associated with oxidative stress after treatment with trifloxysulfuron. This accumulation of mechanisms, in addition to the manifestation of target site resistance, could potentially increase the chance of survival when plants are challenged by different modes of action.

Detection of subgenome bias using an anchored syntenic approach in Eleusine coracana (finger millet)

Background

Finger millet (Eleusine coracana 2n = 4x = 36) is a hardy, nutraceutical, climate change tolerant, orphan crop that is consumed throughout eastern Africa and India. Its genome has been sequenced multiple times, but A and B subgenomes could not be separated because no published genome for E. indica existed. The classification of A and B subgenomes is important for understanding the evolution of this crop and provide a means to improve current and future breeding programs.

Results

We produced subgenome calls for 704 syntenic blocks and inferred A or B subgenomic identity for 59,377 genes 81% of the annotated genes. Phylogenetic analysis of a super matrix containing 455 genes shows high support for A and B divergence within the Eleusine genus. Synonymous substitution rates between A and B genes support A and B calls. The repetitive content on highly supported B contigs is higher than that on similar A contigs. Analysis of syntenic singletons showed evidence of biased fractionation showed a pattern of A genome dominance, with 61% A, 37% B and 1% unassigned, and was further supported by the pattern of loss observed among cyto-nuclear interacting genes.

Conclusion

The evidence of individual gene calls within each syntenic block, provides a powerful tool for inference for subgenome classification. Our results show the utility of a draft genome in resolving A and B subgenomes calls, primarily it allows for the proper polarization of A and B syntenic blocks. There have been multiple calls for the use of phylogenetic inference in subgenome classification, our use of synteny is a practical application in a system that has only one parental genome available.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07447-y.

Mechanisms of evolved herbicide resistance

The widely successful use of synthetic herbicides over the past 70 years has imposed strong and widespread selection pressure, leading to the evolution of herbicide resistance in hundreds of weed species. Both target-site resistance (TSR) and nontarget-site resistance (NTSR) mechanisms have evolved to most herbicide classes. TSR often involves mutations in genes encoding the protein targets of herbicides, affecting the binding of the herbicide either at or near catalytic domains or in regions affecting access to them. Most of these mutations are nonsynonymous SNPs, but polymorphisms in more than one codon or entire codon deletions have also evolved. Some herbicides bind multiple proteins, making the evolution of TSR mechanisms more difficult. Increased amounts of protein target, by increased gene expression or by gene duplication, are an important, albeit less common, TSR mechanism. NTSR mechanisms include reduced absorption or translocation and increased sequestration or metabolic degradation. The mechanisms that can contribute to NTSR are complex and often involve genes that are members of large gene families. For example, enzymes involved in herbicide metabolism-based resistances include cytochromes P450, GSH S-transferases, glucosyl and other transferases, aryl acylamidase, and others. Both TSR and NTSR mechanisms can combine at the individual level to produce higher resistance levels. The vast array of herbicide-resistance mechanisms for generalist (NTSR) and specialist (TSR and some NTSR) adaptations that have evolved over a few decades illustrate the evolutionary resilience of weed populations to extreme selection pressures. These evolutionary processes drive herbicide and herbicide-resistant crop development and resistance management strategies.

Transcriptional Sequencing and Gene Expression Analysis of Various Genes in Fruit Development of Three Different Black Pepper (Piper nigrum L.) Varieties

Black pepper (Piper nigrum) is a vital spice crop with uses ranging from culinary to pharmacological applications. However, limited genetic information has constrained the understanding of the molecular regulation of flower and fruit development in black pepper. In this study, a comparison among three different black pepper varieties, Semengok Aman (SA), Kuching (KC), and Semengok 1 (S1), with varying fruit characteristics was used to provide insight on the genetic regulation of flower and fruit development. Next-generation sequencing (NGS) technology was used to determine the flower and fruit transcriptomes by sequencing on an Illumina HiSeq 2500 platform followed by de novo assembly using SOAPdenovo-Trans. The high-quality assembly of 66,906 of unigenes included 64.4% of gene sequences (43,115) with similarity to one or more protein sequences from the GenBank database. Annotation with Blast2Go assigned 37,377 genes to one or more Gene Ontology terms. Of these genes, 5,874 genes were further associated with the biological pathways recorded in the KEGG database. Comparison of flower and fruit transcriptome data from the three different black pepper varieties revealed a large number of DEGs between flower and fruit of the SA variety. Gene Ontology (GO) enrichment analysis further supports functions of DEGs between flower and fruit in the categories of carbohydrate metabolic processes, embryo development, and DNA metabolic processes while the DEGs in fruit relate to biosynthetic process, secondary metabolic process, and catabolic process. The enrichment of DEGs in KEGG pathways was also investigated, and a large number of genes were found to belong to the nucleotide metabolism and carbohydrate metabolism categories. Gene expression profiling of flower formation-related genes reveals that other than regulating the flowering in black pepper, the flowering genes might also be implicated in the fruit development process. Transcriptional analysis of sugar transporter and carbohydrate metabolism genes in different fruit varieties suggested that the carbohydrate metabolism in black pepper fruit is developmentally regulated, and some genes might serve as potential genes for future crop quality improvement. Study on the piperine-related gene expression analysis suggested that lysine-derived products might present in all stages of fruit development, but the transportation was only active at the early stage of fruit development. These results indicate several candidate genes related to the development of flower and fruit in black pepper and provide a resource for future functional analysis and potentially for future crop improvement.

Transcriptome Analysis Reveals Unique Relationships Among Eleusine Species and Heritage of Eleusine coracana

Relationships in the genus Eleusine were obtained through transcriptome analysis. Eleusine coracana (E. coracana ssp. coracana), also known as finger millet, is an allotetraploid minor crop primarily grown in East Africa and India. Domesticated E. coracana evolved from wild E. africana (E. coracana ssp. africana) with the maternal genome donor largely supported to be E. indica; however, the paternal genome donor remains elusive. We developed transcriptomes for six Eleusine species from fully developed seedlings using Illumina technology and three de novo assemblers (Trinity, Velvet, and SOAPdenovo2) with the redundancy-reducing EvidentialGene pipeline. Mapping E. coracana reads to the chloroplast genes of all Eleusine species detected fewer variants between E. coracana and E. indica compared to all other species. Phylogenetic analysis further supports E. indica as the maternal parent of E. coracana and E. africana, in addition to a close relationship between E. indica and E. tristachya, and between E. floccifolia and E. multiflora, and E. intermedia as a separate group. A close relationship between E. floccifolia and E. multiflora was unexpected considering they are reported to have distinct nuclear genomes, BB and CC, respectively. Further, it was expected that E. intermedia and E. floccifolia would have a closer relationship considering they have similar nuclear genomes, AB and BB, respectively. A rethinking of the labeling of ancestral genomes of E. floccifolia, E. multiflora, and E. intermedia is maybe needed based on this data.

EPSPS duplication and mutation involved in glyphosate resistance in the allotetraploid weed species Poa annua L

BACKGROUND

Poa annua is a widespread winter annual weed species in California. Recently, poor control of this species with glyphosate was reported by growers in an almond orchard in California with a history of repetitive glyphosate use. The objectives of this research were to evaluate the level of glyphosate resistance in a developed S4 P. annua line (R) and identify the mechanisms of resistance involved.

RESULTS

Whole-plant dose-response experiments confirmed glyphosate resistance in R, which required 18-fold more glyphosate to achieve a 50% growth reduction compared with a susceptible line (S), results that were supported by the lower shikimate accumulation observed in R compared with S. No differences in glyphosate absorption, translocation, or metabolism were observed, suggesting that non-target-site mechanisms of resistance are not involved in the resistance phenotype. A missense single nucleotide polymorphism was observed in EPSPS coding position 106 in R, resulting in a leucine to proline substitution. This polymorphism was observed exclusively in P. supina EPSPS homeologs. A seven-fold increase in the number of copies of EPSPS alleles was observed in R compared with S.

CONCLUSIONS

We report the first case of glyphosate resistance associated with both EPSPS duplication and target-site mutation at position 106, leading to high levels of glyphosate resistance in the allotetraploid weed species Poa annua L. Data obtained in this research will be useful for the development of diagnostic tools for rapid glyphosate resistance identification, monitoring and containment. © 2018 Society of Chemical Industry.

Transcriptome-derived investigation of biosynthesis of quinolizidine alkaloids in narrow-leafed lupin (Lupinus angustifolius L.) highlights candidate genes linked to iucundus locus

Unravelling the biosynthetic pathway of quinolizidine alkaloids (QAs), regarded as antinutritional compounds of narrow-leafed lupin (NLL) seeds, is fundamental to best exploit NLL as food or feed. We investigated 12 candidate genes connected to QA biosynthesis, selecting them by transcriptomic and genomic approaches, from the landscape of genes differentially expressed in leaves of the high- and low-alkaloid NLL accessions. Linkage analysis enabled the assessment of the location of the candidate genes in relation to iucundus, a major locus of unknown identity, that confers reduced QA content in seeds. The key finding was the identification of APETALA2/ethylene response transcription factor, RAP2-7, cosegregating with the iucundus locus and located within a region with highly significant QTLs that affect QA composition. We additionally identified a 4-hydroxy-tetrahydrodipicolinate synthase (DHDPS) gene involved in L-lysine biosynthesis as being closely linked to iucundus. The distributed location of other remaining candidates (including previously known QA genes) across different linkage groups, also indirectly supports the transcription factor as a possible regulator of lupin alkaloid biosynthesis. Our findings provide crucial insight into QA biosynthesis in NLL. Additionally, we evaluated and selected appropriate reference genes for qRT-PCRs to analyse the expression levels of QA genes in NLL.

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.