Search

The mitochondrial genome of Liriodendron tulipifera exhibits many ancestral angiosperm features and a remarkably slow evolutionary rate, while mitochondrial genomes of other magnoliids remain yet to be characterized. We assembled nine new mitochondrial genomes, representing all genera of perianth-bearing Piperales, as well as for a member of the sister clade: three complete or nearly complete mitochondrial genomes from Aristolochiaceae and six additional draft assemblies including Thottea, Asaraceae, Lactoridaceae, and Hydnoraceae. For comparative purpose, a complete mitochondrial genome was assembled for Saururus, a member of the perianth-less Piperales. The average number of short repeats (50–99 bp) was much larger in genus Aristolochia than in other angiosperm mitochondrial genomes, and approximately 30% of repeats ( T:A substitutions of all other investigated angiosperm groups. Our study reports the first mitochondrial genomes for Piperales and uses this new information for a better understanding of the evolutionary patterns of magnoliids and angiosperms in general.

Mycoheterotrophy is an alternative nutritional strategy whereby plants obtain sugars and other nutrients from soil fungi. Mycoheterotrophy and associated loss of photosynthesis have evolved repeatedly in plants, particularly in monocots. Although reductive evolution of plastomes in mycoheterotrophs is well documented, the dynamics of nuclear genome evolution remains largely unknown. Transcriptome datasets were generated from four mycoheterotrophs in three families (Orchidaceae, Burmanniaceae, Triuridaceae) and related green plants and used for phylogenomic analyses to resolve relationships among the mycoheterotrophs, their relatives, and representatives across the monocots. Phylogenetic trees based on 602 genes were mostly congruent with plastome phylogenies, except for an Asparagales + Liliales clade inferred in the nuclear trees. Reduction and loss of chlorophyll synthesis and photosynthetic gene expression and relaxation of purifying selection on retained genes were progressive, with greater loss in older nonphotosynthetic lineages. One hundred seventy-four of 1375 plant benchmark universally conserved orthologous genes were undetected in any mycoheterotroph transcriptome or the genome of the mycoheterotrophic orchid Gastrodia but were expressed in green relatives, providing evidence for massively convergent gene loss in nonphotosynthetic lineages. We designate this set of deleted or undetected genes Missing in Mycoheterotrophs (MIM). MIM genes encode not only mainly photosynthetic or plastid membrane proteins but also a diverse set of plastid processes, genes of unknown function, mitochondrial, and cellular processes. Transcription of a photosystem II gene (psb29) in all lineages implies a nonphotosynthetic function for this and other genes retained in mycoheterotrophs. Nonphotosynthetic plants enable novel insights into gene function as well as gene expression shifts, gene loss, and convergence in nuclear genomes.

Aristolochia, a genus in the magnoliid order Piperales, has been famous for centuries for its highly specialized flowers and wide medicinal applications. Here, we present a new, high-quality genome sequence of Aristolochia fimbriata, a species that, similar to Amborella trichopoda, lacks further whole-genome duplications since the origin of extant angiosperms. As such, the A. fimbriata genome is an excellent reference for inferences of angiosperm genome evolution, enabling detection of two novel whole-genome duplications in Piperales and dating of previously reported whole-genome duplications in other magnoliids. Genomic comparisons between A. fimbriata and other angiosperms facilitated the identification of ancient genomic rearrangements suggesting the placement of magnoliids as sister to monocots, whereas phylogenetic inferences based on sequence data we compiled yielded ambiguous relationships. By identifying associated homologues and investigating their evolutionary histories and expression patterns, we revealed highly conserved floral developmental genes and their distinct downstream regulatory network that may contribute to the complex flower morphology in A. fimbriata. Finally, we elucidated the genetic basis underlying the biosynthesis of terpenoids and aristolochic acids in A. fimbriata.

Fitness responses to environment can shape species distributions, though opposing eco-evolutionary processes can obscure environmental effects. For example, host specificity influences parasite dynamics, but is unclear how adaptation of parasites to local host communities may scale up to continental distributions. Here, we develop a macroecological framework to determine how host community structure affects the distribution of specialist and generalist populations of Striga hermonthica, an African parasitic plant of cereal crops. Combining data from global crop production and parasite experimental trials, we find that parasites perform best on the host species that is most common in their location of origin. Moreover, niche model contrasts predict parasite specialization on two hosts that evolved alongside Striga during domestication (pearl millet and sorghum), indicating that specialist parasites may be most likely to occur where host niches differ most in multivariate environmental space. Our study demonstrates that patterns of parasite local adaptation to host communities can emerge at continental scales and that differential environmental tolerances of hosts indirectly shape the distribution of specialist and generalist parasites. By predicting spatial dynamics of parasite specialization versus generalization directly from environmental data, our approach may help inform current and future management of pests and disease.

PremiseParasitic plants and their hosts are model systems for studying genetic variation in species interactions across environments. The parasitic plant Striga hermonthica (witchweed) attacks a range of cereal crop hosts in Africa and exhibits substantial variation in performance on different host species. Some of this variation is due to local adaptation, but the genetic basis of specialization on certain hosts is unknown.MethodsTo identify genomic regions that are strongly differentiated between parasites attacking different host species, we present an alignment-free analysis of S. hermonthica population diversity using whole genome sequencing (WGS) data for 68 individuals from western Kenya. We validate our findings with germination experiments and analyses based on a de novo assembled draft genome.ResultsReference-free and reference-based analyses suggest that only a small portion of the S. hermonthica genome is strongly differentiated by host species in populations from western Kenya. Analysis of host-associated k-mers implicated genes involved in development of the parasite haustorium (a specialized structure used to establish vascular connections with host roots) and a potential role of chemocyanins in molecular host-parasitic plant interactions. Conversely, no phenotypic or genomic evidence was observed suggesting host-specific selection on parasite response to strigolactones, hormones exuded by host roots and required for parasite germination.ConclusionsThis study demonstrates the utility of WGS for plant species with large, complex genomes and no available reference. Contrasting with theory emphasizing the role of early recognition loci for host specificity, our findings support host-specific selection on later interaction stages, recurring each generation after homogenizing gene flow.

The rapid development of sequencing technologies has led to a deeper understanding of horticultural plant genomes. However, experimental evidence connecting genes to important agronomic traits is still lacking in most non-model organisms. For instance, the genetic mechanisms underlying plant architecture are poorly understood in pome fruit trees, creating a major hurdle in developing new cultivars with desirable architecture, such as dwarfing rootstocks in European pear (Pyrus communis). Further, the quality and content of genomes vary widely. Therefore, it can be challenging to curate a list of genes with high-confidence gene models across reference genomes. This is often an important first step towards identifying key genetic factors for important traits. Here we present a draft genome of P. communis ‘d’Anjou’ and an improved assembly of the latest P. communis ‘Bartlett’ genome. To study gene families involved in tree architecture in European pear and other rosaceous species, we developed a workflow using a collection of bioinformatic tools towards curation of gene families of interest across genomes. This lays the groundwork for future functional studies in pear tree architecture. Importantly, our workflow can be easily adopted for other plant genomes and gene families of interest.

Genomic structural variants (SVs) can play important roles in adaptation and speciation. Yet the overall fitness effects of SVs are poorly understood, partly because accurate population-level identification of SVs requires multiple high-quality genome assemblies. Here, we use 31 chromosome-scale, haplotype-resolved genome assemblies of Theobroma cacao—an outcrossing, long-lived tree species that is the source of chocolate—to investigate the fitness consequences of SVs in natural populations. Among the 31 accessions, we find over 160,000 SVs, which together cover eight times more of the genome than single-nucleotide polymorphisms and short indels (125 versus 15 Mb). Our results indicate that a vast majority of these SVs are deleterious: they segregate at low frequencies and are depleted from functional regions of the genome. We show that SVs influence gene expression, which likely impairs gene function and contributes to the detrimental effects of SVs. We also provide empirical support for a theoretical prediction that SVs, particularly inversions, increase genetic load through the accumulation of deleterious nucleotide variants as a result of suppressed recombination. Despite the overall detrimental effects, we identify individual SVs bearing signatures of local adaptation, several of which are associated with genes differentially expressed between populations. Genes involved in pathogen resistance are strongly enriched among these candidates, highlighting the contribution of SVs to this important local adaptation trait. Beyond revealing empirical evidence for the evolutionary importance of SVs, these 31 de novo assemblies provide a valuable resource for genetic and breeding studies in T. cacao.

Horizontal gene transfer (HGT), the movement and genomic integration of DNA across species boundaries, is commonly associated with bacteria and other microorganisms, but functional HGT (fHGT) is increasingly being recognized in heterotrophic parasitic plants that obtain their nutrients and water from their host plants through direct haustorial feeding. Here, in the holoparasitic stem parasite Cuscuta, we identify 108 transcribed and probably functional HGT events in Cuscuta campestris and related species, plus 42 additional regions with host-derived transposon, pseudogene and non-coding sequences. Surprisingly, 18 Cuscuta fHGTs were acquired from the same gene families by independent HGT events in Orobanchaceae parasites, and the majority are highly expressed in the haustorial feeding structures in both lineages. Convergent retention and expression of HGT sequences suggests an adaptive role for specific additional genes in parasite biology. Between 16 and 20 of the transcribed HGT events are inferred as ancestral in Cuscuta based on transcriptome sequences from species across the phylogenetic range of the genus, implicating fHGT in the successful radiation of Cuscuta parasites. Genome sequencing of C. campestris supports transfer of genomic DNA-rather than retroprocessed RNA-as the mechanism of fHGT. Many of the C. campestris genes horizontally acquired are also frequent sources of 24-nucleotide small RNAs that are typically associated with RNA-directed DNA methylation. One HGT encoding a leucine-rich repeat protein kinase overlaps with a microRNA that has been shown to regulate host gene expression, suggesting that HGT-derived parasite small RNAs may function in the parasite-host interaction. This study enriches our understanding of HGT by describing a parasite-host system with unprecedented gene exchange that points to convergent evolution of HGT events and the functional importance of horizontally transferred coding and non-coding sequences.

Host-specific interactions can maintain genetic and phenotypic diversity in parasites that attack multiple host species. Host diversity, in turn, may promote parasite diversity by selection for genetic divergence or plastic responses to host type. The parasitic weed purple witchweed [Striga hermonthica (Delile) Benth.] causes devastating crop losses in sub-Saharan Africa and is capable of infesting a wide range of grass hosts. Despite some evidence for host adaptation and host-by-Striga genotype interactions, little is known about intraspecific Striga genomic diversity. Here we present a study of transcriptomic diversity in populations of S. hermonthica growing on different hosts (maize [Zea mays L.] vs. grain sorghum [Sorghum bicolor (L.) Moench]). We examined gene expression variation and differences in allelic frequency in expressed genes of aboveground tissues from populations in western Nigeria parasitizing each host. Despite low levels of host-based genome-wide differentiation, we identified a set of parasite transcripts specifically associated with each host. Parasite genes in several different functional categories implicated as important in host–parasite interactions differed in expression level and allele on different hosts, including genes involved in nutrient transport, defense and pathogenesis, and plant hormone response. Overall, we provide a set of candidate transcripts that demonstrate host-specific interactions in vegetative tissues of the emerged parasite S. hermonthica. Our study shows how signals of host-specific processes can be detected aboveground, expanding the focus of host–parasite interactions beyond the haustorial connection.

Temperature conditioning is an effective approach for enhancing some aspects of postharvest apple fruit quality, particularly avoidance of chilling injury. The molecular mechanisms of warming regimes that alter apple fruit quality are poorly understood, as is how warming contributes to avoiding development of the peel disorder superficial scald. Here we report transcriptional responses of ‘Granny Smith’ fruit peel during the early phases of long-term cold storage in response to intermittent warming, an ostensibly organic compliant strategy that effectively reduced scald incidence. We observed two temporally distinct classes of gene expression, which were discovered by co-expression network analyses. One profile is largely concordant with recovery from chilling stress, whereas the other reveals transient shifts marked by hormone signaling and transcription and translation machinery. Altogether, our analyses point to novel aspects of superficial scald etiology and circumscribes a list of candidate genes that may be useful to uncover molecular processes that promote, as well as mitigate, the peel disorder.

International audience; The genetic basis and evolution of sex determination in dioecious plants is emerging as an active area of research with exciting advances in genome sequencing and analysis technologies. As the sole species within the sister lineage to all other extant flowering plants, Amborella trichopoda is an important model for understanding the evolution and development of flowers. Plants typically produce only male or female flowers, but sex determination mechanisms are unknown for the species. Sequence data derived from plants of natural origin and an F1 mapping population were used to identify sex-linked genes and the nonrecombining region. Amborella trichopoda has a ZW sex determination system. Analysis of genes in a 4 Mb nonrecombining sex-determination region reveals recent divergence of Z and W gametologs, and few Z- and W-specific genes. The sex chromosomes of A. trichopoda evolved less than 16.5 Myr ago, long after the divergence of the extant angiosperms.

Aristolochia, a genus in the magnoliid order Piperales, has been famous for centuries for its highly specialized flowers and wide medicinal applications. Here, we present a new, high-quality genome sequence of Aristolochia fimbriata, a species that, similar to Amborella trichopoda, lacks further whole-genome duplications since the origin of extant angiosperms. As such, the A. fimbriata genome is an excellent reference for inferences of angiosperm genome evolution, enabling detection of two novel whole-genome duplications in Piperales and dating of previously reported whole-genome duplications in other magnoliids. Genomic comparisons between A. fimbriata and other angiosperms facilitated the identification of ancient genomic rearrangements suggesting the placement of magnoliids as sister to monocots, whereas phylogenetic inferences based on sequence data we compiled yielded ambiguous relationships. By identifying associated homologues and investigating their evolutionary histories and expression patterns, we revealed highly conserved floral developmental genes and their distinct downstream regulatory network that may contribute to the complex flower morphology in A. fimbriata. Finally, we elucidated the genetic basis underlying the biosynthesis of terpenoids and aristolochic acids in A. fimbriata.

Sex determination is poorly understood in plants. Amborella trichopoda is a well-known plant model for evo-devo studies, which is also dioecious (has male and female individuals), with an unknown sex determination mechanism. A. trichopoda is a “sex switcher”, which points to possible environmental factors that act on sex, but populations grown from seed under greenhouse conditions exhibit a 50:50 sex ratio, which indicates the operation of genetic factors. Here, we use a new method (SDpop) to identify sex-linked genes from genotyping data of male and female individuals sampled in the field, and find that A. trichopoda has a ZW sex-chromosome system. The sex-linked genes map to a 4 Mb sex-determining region on chromosome 9. The low extent of ZW divergence suggests these sex chromosomes are of recent origin, which is consistent with dioecy being derived character in the A. trichopoda lineage. Our work has uncovered clearly formed sex chromosomes in a species in which both genetic and environmental factors can influence sex.One Sentence SummaryAmborella trichopoda, a dioecious species in which both genetics and the environment influence sex, possesses a pair of quite recently evolved ZW chromosomes.

Fitness responses to environment can shape species distributions, though opposing eco-evolutionary processes can obscure environmental effects. For example, host specificity influences parasite dynamics, but is unclear how specialization of individual parasites or populations scales up to continental distributions. Here, we develop a macroecological framework to determine how host community structure affects continent-scale specialization in Striga hermonthica, an African parasitic plant of cereal crops. We find regional abundance of hosts in cultivated cereal communities is associated with parasite specialization observed in experiments. Moreover, abiotic environment at location of origin predicts parasite performance on pearl millet and sorghum but not maize, possibly due to the shorter coevolutionary history for maize and Striga. Our study demonstrates that patterns of parasite local adaptation to host communities can emerge at continental scales and that differential environmental tolerances of hosts indirectly shape the distribution of specialist and generalist parasites.

Host-parasite coevolution can maintain high levels of genetic diversity in traits involved in species interactions. In many systems, host traits exploited by parasites are constrained by use in other functions, leading to complex selective pressures across space and time. Here, we study genome-wide variation in the staple crop Sorghum bicolor (L.) Moench and its association with the parasitic weed Striga hermonthica (Delile) Benth., a major constraint to food security in Africa. We hypothesize that geographic selection mosaics across gradients of parasite occurrence maintain genetic diversity in sorghum landrace resistance. Suggesting a role in local adaptation to parasite pressure, multiple independent loss-of-function alleles at sorghum LOW GERMINATION STIMULANT 1 (LGS1) are broadly distributed among African landraces and geographically associated with S. hermonthica occurrence. However, low frequency of these alleles within S. hermonthica-prone regions and their absence elsewhere implicate potential tradeoffs restricting their fixation. LGS1 is thought to cause resistance by changing stereochemistry of strigolactones, hormones that control plant architecture and belowground signaling to mycorrhizae and are required to stimulate parasite germination. Consistent with tradeoffs, we find signatures of balancing selection surrounding LGS1 and other candidates from analysis of genome-wide associations with parasite distribution. Experiments with CRISPR-Cas9 edited sorghum further indicate the benefit of LGS1-mediated resistance strongly depends on parasite genotype and abiotic environment and comes at the cost of reduced photosystem gene expression. Our study demonstrates long-term maintenance of diversity in host resistance genes across smallholder agroecosystems, providing a valuable comparison to both industrial farming systems and natural communities.SIGNIFICANCE STATEMENTUnderstanding co-evolution in crop-parasite systems is critical to management of myriad pests and pathogens confronting modern agriculture. In contrast to wild plant communities, parasites in agricultural ecosystems are usually expected to gain the upper hand in co-evolutionary ‘arms races’ due to limited genetic diversity of host crops in cultivation. Here, we develop a framework to characterize associations between genome variants in global landraces (traditional varieties) of the staple crop sorghum with the distribution of the devastating parasitic weed Striga hermonthica. We find long-term maintenance of diversity in genes related to parasite resistance, highlighting an important role of host adaptation for co-evolutionary dynamics in smallholder agroecosystems.

The tropical tree Theobroma cacao is the source of chocolate, and its seeds are a major export from many producing countries in Central and South America, Africa, and Asia. Every year, 30–40% of pre-harvest yield is lost due to disease damage. Host plant resistance is the most efficient and environmentally friendly approach for disease management. Historically, cacao germplasm resources have been underutilized in efforts to introduce novel sources of disease tolerance into breeding programs. Maintenance of cacao germplasm also relies on clonally propagated live collections, as cacao seeds do not exhibit dormancy and cannot be stored for more than a few weeks. In this study, we use a 90 SNP array to verify genetic identity of a set of clones in the International Cocoa Collection at CATIE, Costa Rica, and assign the clones into known genetic groups. We also used a detached leaf inoculation technique to measure the susceptibility of 60 genotypes to Phytophthora palmivora, a major cacao pathogen with global importance. We identified 24 genotypes with disease tolerance statistically similar to a standard tolerant variety (SCA6) and another 24 which performed similarly to a standard susceptible variety (ICS1). Our results indicate that each of the four included genetic groups show variability for quantitative resistance to P. palmivora. These results provide a foundation for future genomic and transcriptomic analysis of disease tolerance and susceptibility in the field at CATIE and provide guidelines for breeders searching for novel sources of tolerance that can be introduced into breeding programs.

Trans-species small regulatory RNAs (sRNAs) are delivered to host plants from diverse pathogens and parasites and can target host mRNAs. How trans-species sRNAs can be effective on diverse hosts has been unclear. Multiple species of the parasitic plant Cuscuta produce trans-species sRNAs that collectively target many host mRNAs. Confirmed target sites are nearly always in highly conserved, protein-coding regions of host mRNAs. Cuscuta trans-species sRNAs can be grouped into superfamilies that have variation in a three-nucleotide period. These variants compensate for synonymous-site variation in host mRNAs. By targeting host mRNAs at highly conserved protein-coding sites, and simultaneously expressing multiple variants to cover synonymous-site variation, Cuscuta trans-species sRNAs may be able to successfully target homologous mRNAs from diverse hosts.One Sentence SummaryThe parasitic plant Cuscuta produces a diverse set of sRNAs that compensate for sequence variation in mRNA targets in diverse hosts.

In tropical forests, rarer species show increased sensitivity to species-specific soil pathogens and more negative effects of conspecific density on seedling survival (NDD). These patterns suggest a connection between ecology and immunity, perhaps because small population size disproportionately reduces genetic diversity of hyperdiverse loci such as immunity genes. In an experiment examining seedling roots from six species in one tropical tree community, we found that smaller populations have reduced amino acid diversity in pathogen resistance (R) genes but not the transcriptome in general. Normalized R gene amino acid diversity varied with local abundance and prior measures of differences in sensitivity to conspecific soil and NDD. After exposure to live soil, species with lower R gene diversity had reduced defence gene induction, more cosusceptibility of maternal cohorts to colonization by potentially pathogenic fungi, reduced root growth arrest (an R gene-mediated response) and their root-associated fungi showed lower induction of self-defence (antioxidants). Local abundance was not related to the ability to induce immune responses when pathogen recognition was bypassed by application of salicylic acid, a phytohormone that activates defence responses downstream of R gene signalling. These initial results support the hypothesis that smaller local tree populations have reduced R gene diversity and recognition-dependent immune responses, along with greater cosusceptibility to species-specific pathogens that may facilitate disease transmission and NDD. Locally rare species may be less able to increase their equilibrium abundance without genetic boosts to defence via immigration of novel R gene alleles from a larger and more diverse regional population.

Separating footprints of adaptation from demography is challenging. When selection has acted on a single locus with major effect, this issue can be alleviated through signatures left by selective sweeps. However, as adaptation is often driven by small allele frequency shifts at many loci, studies focusing on single genes are able to identify only a small portion of genomic variants responsible for adaptation. In face of this challenge, we utilize coexpression information to search for signals of polygenetic adaptation in Theobroma cacao, a tropical tree species that is the source of chocolate. Using transcriptomics and a weighted correlation network analysis, we group genes with similar expression patterns into functional modules. We then ask whether modules enriched for specific biological processes exhibit cumulative effects of differential selection in the form of high FST and dXY between populations. Indeed, modules putatively involved in protein modification, flowering, and water transport show signs of polygenic adaptation even though individual genes that are members of those groups do not bear strong signatures of selection. Modeling of demography, background selection, and the effects of genomic features reveal that these patterns are unlikely to arise by chance. We also find that specific modules are enriched for signals of strong or relaxed purifying selection, with one module bearing signs of adaptive differentiation and an excess of deleterious mutations. Our results provide insight into polygenic adaptation and contribute to understanding of population structure, demographic history, and genome evolution in T. cacao.

Background Parasitic plants engage in a complex molecular dialog with potential host plants to identify a host and overcome host defenses to initiate development of the parasitic feeding organ, the haustorium, invade host tissues, and withdraw water and nutrients. While one of two critical signaling events in the parasitic plant life cycle (germination via stimulant chemicals) has been relatively well-studied, the signaling event that triggers haustorium formation remains elusive. Elucidation of this poorly understood molecular dialogue will shed light on plant-plant communication, parasitic plant physiology, and the evolution of parasitism in plants. Results Here we present an experimental framework that develops easily quantifiable contrasts for the facultative generalist parasitic plant, Triphysaria, as it feeds across a broad range of diverse flowering plants. The contrasts, including variable parasite growth form and mortality when grown with different hosts, suggest a dynamic and host-dependent molecular dialogue between the parasite and host. Finally, by comparing transcriptome datasets from attached versus unattached parasites we gain insight into some of the physiological processes that are altered during parasitic behavior including shifts in photosynthesis-related and stress response genes. Conclusions This work sheds light on Triphysaria’s parasitic life habit and is an important step towards understanding the mechanisms of haustorium initiation factor perception, a unique form of plant-plant communication. Electronic supplementary material The online version of this article (10.1186/s12870-019-1856-1) contains supplementary material, which is available to authorized users.

This protocol is part of a collection of eighteen protocols used to isolate total RNA from plant tissue. (RNA Isolation from Plant Tissue Collection: https://www.protocols.io/view/rna-isolation-from-plant-tissue-439gyr6) and was originally published as part of Appendix S1 of "Evaluating Methods for Isolating Total RNA and Predicting the Success of Sequencing Phylogenetically Diverse Plant Transcriptomes" Marc T. J. Johnson et al. PLOS ONE, November 21, 2012. https://doi.org/10.1371/journal.pone.0050226

Parasitic weeds of the family Orobanchaceae attach to the roots of host plants via haustoria capable of drawing nutrients from host vascular tissue. The connection of the haustorium to the host marks a shift in parasite metabolism from autotrophy to at least partial heterotrophy, depending on the level of parasite dependence. Species within the family Orobanchaceae span the spectrum of host nutrient dependency, yet the diversity of parasitic plant metabolism remains poorly understood, particularly during the key metabolic shift surrounding haustorial attachment. Comparative profiling of major metabolites in the obligate holoparasite Phelipanche aegyptiaca and the facultative hemiparasite Triphysaria versicolor before and after attachment to the hosts revealed several metabolic shifts implicating remodeling of energy and amino acid metabolism. After attachment, both parasites showed metabolite profiles that were different from their respective hosts. In P. aegyptiaca, prominent changes in metabolite profiles were also associated with transitioning between different tissue types before and after attachment, with aspartate levels increasing significantly after the attachment. Based on the results from 15N labeling experiments, asparagine and/or aspartate-rich proteins were enriched in host-derived nitrogen in T. versicolor. These results point to the importance of aspartate and/or asparagine in the early stages of attachment in these plant parasites and provide a rationale for targeting aspartate-family amino acid biosynthesis for disrupting the growth of parasitic weeds. Published version

Background Root parasitic weeds are a major constraint to crop production worldwide causing significant yearly losses in yield and economic value. These parasites cause their destruction by attaching to their hosts with a unique organ, the haustorium, that allows them to obtain the nutrients (sugars, amino acids, etc.) needed to complete their lifecycle. Parasitic weeds differ in their nutritional requirements and degree of host dependency and the differential expression of sugar transporters is likely to be a critical component in the parasite’s post-attachment survival. Results We identified gene families encoding monosaccharide transporters (MSTs), sucrose transporters (SUTs), and SWEETs (Sugars Will Eventually be Exported Transporters) in three root-parasitic weeds differing in host dependency: Triphysaria versicolor (facultative hemiparasite), Phelipanche aegyptiaca (holoparasite), and Striga hermonthica (obligate hemiparasite). The phylogenetic relationship and differential expression profiles of these genes throughout parasite development were examined to uncover differences existing among parasites with different levels of host dependence. Differences in estimated gene numbers are found among the three parasites, and orthologs within the different sugar transporter gene families are found to be either conserved among the parasites in their expression profiles throughout development, or to display parasite-specific differences in developmentally-timed expression. For example, MST genes in the pGLT clade express most highly before host connection in Striga and Triphysaria but not Phelipanche, whereas genes in the MST ERD6-like clade are highly expressed in the post-connection growth stages of Phelipanche but highest in the germination and reproduction stages in Striga. Whether such differences reflect changes resulting from differential host dependence levels is not known. Conclusions While it is tempting to speculate that differences in estimated gene numbers and expression profiles among members of MST, SUT and SWEET gene families in Phelipanche, Striga and Triphysaria reflect the parasites’ levels of host dependence, additional evidence that altered transporter gene expression is causative versus consequential is needed. Our findings identify potential targets for directed manipulation that will allow for a better understanding of the nutrient transport process and perhaps a means for controlling the devastating effects of these parasites on crop productivity. Electronic supplementary material The online version of this article (10.1186/s12870-019-1786-y) contains supplementary material, which is available to authorized users.

Plastid genomes (plastomes) vary enormously in size and gene content among the many lineages of nonphotosynthetic plants, but key lineages remain unexplored. We therefore investigated plastome sequence and expression in the holoparasitic and morphologically bizarre Balanophoraceae. The two Balanophora plastomes examined are remarkable, exhibiting features rarely if ever seen before in plastomes or in any other genomes. At 15.5 kb in size and with only 19 genes, they are among the most reduced plastomes known. They have no tRNA genes for protein synthesis, a trait found in only three other plastid lineages, and thus Balanophora plastids must import all tRNAs needed for translation. Balanophora plastomes are exceptionally compact, with numerous overlapping genes, highly reduced spacers, loss of all cis -spliced introns, and shrunken protein genes. With A+T contents of 87.8% and 88.4%, the Balanophora genomes are the most AT-rich genomes known save for a single mitochondrial genome that is merely bloated with AT-rich spacer DNA. Most plastid protein genes in Balanophora consist of ≥90% AT, with several between 95% and 98% AT, resulting in the most biased codon usage in any genome described to date. A potential consequence of its radical compositional evolution is the novel genetic code used by Balanophora plastids, in which TAG has been reassigned from stop to tryptophan. Despite its many exceptional properties, the Balanophora plastome must be functional because all examined genes are transcribed, its only intron is correctly trans -spliced, and its protein genes, although highly divergent, are evolving under various degrees of selective constraint.

Parasitic plants in the genus Striga, commonly known as witchweeds, cause major crop losses in sub-Saharan Africa and pose a threat to agriculture worldwide. An understanding of Striga parasite biology, which could lead to agricultural solutions, has been hampered by the lack of genome information. Here, we report the draft genome sequence of Striga asiatica with 34,577 predicted protein-coding genes, which reflects gene family contractions and expansions that are consistent with a three-phase model of parasitic plant genome evolution. Striga seeds germinate in response to host-derived strigolactones (SLs) and then develop a specialized penetration structure, the haustorium, to invade the host root. A family of SL receptors has undergone a striking expansion, suggesting a molecular basis for the evolution of broad host range among Striga spp. We found that genes involved in lateral root development in non-parasitic model species are coordinately induced during haustorium development in Striga, suggesting a pathway that was partly co-opted during the evolution of the haustorium. In addition, we found evidence for horizontal transfer of host genes as well as retrotransposons, indicating gene flow to S. asiatica from hosts. Our results provide valuable insights into the evolution of parasitism and a key resource for the future development of Striga control strategies.

Because novel environmental conditions alter the selection pressure on genes or entire subgenomes, adaptive and nonadaptive changes will leave a measurable signature in the genomes, shaping their molecular evolution. We present herein a model of the trajectory of plastid genome evolution under progressively relaxed functional constraints during the transition from autotrophy to a nonphotosynthetic parasitic lifestyle. We show that relaxed purifying selection in all plastid genes is linked to obligate parasitism, characterized by the parasite’s dependence on a host to fulfill its life cycle, rather than the loss of photosynthesis. Evolutionary rates and selection pressure coevolve with macrostructural and microstructural changes, the extent of functional reduction, and the establishment of the obligate parasitic lifestyle. Inferred bursts of gene losses coincide with periods of relaxed selection, which are followed by phases of intensified selection and rate deceleration in the retained functional complexes. Our findings suggest that the transition to obligate parasitism relaxes functional constraints on plastid genes in a stepwise manner. During the functional reduction process, the elevation of evolutionary rates reaches several new rate equilibria, possibly relating to the modified protein turnover rates in heterotrophic plastids.

Parasitic plants in the Orobanchaceae cause serious agricultural problems worldwide. Parasitic plants develop a multicellular infectious organ called a haustorium after recognition of host-released signals. To understand the molecular events associated with host signal perception and haustorium development, we identified differentially regulated genes expressed during early haustorium development in the facultative parasite Phtheirospermum japonicum using a de novo assembled transcriptome and a customized microarray. Among the genes that were upregulated during early haustorium development, we identified YUC3, which encodes a functional YUCCA (YUC) flavin monooxygenase involved in auxin biosynthesis. YUC3 was specifically expressed in the epidermal cells around the host contact site at an early time point in haustorium formation. The spatio-temporal expression patterns of YUC3 coincided with those of the auxin response marker DR5, suggesting generation of auxin response maxima at the haustorium apex. Roots transformed with YUC3 knockdown constructs formed haustoria less frequently than nontransgenic roots. Moreover, ectopic expression of YUC3 at the root epidermal cells induced the formation of haustorium-like structures in transgenic P. japonicum roots. Our results suggest that expression of the auxin biosynthesis gene YUC3 at the epidermal cells near the contact site plays a pivotal role in haustorium formation in the root parasitic plant P. japonicum.

In this study we isolated and performed in silico analysis of a putative coclaurine N-methyltransferase (CNMT) from the basal angiosperm Aristolochia fimbriata. The Aristolochiaceae plant family produces alkaloids similar to the Papavaraceae family, and CNMTs are central enzymes in biosynthesis pathways producing compounds of ethnopharmacological interest. We used bioinformatics and computational tools to predict a three-dimensional homology model and to investigate the putative function of the protein and its mechanism for methylation. The putative CNMT is a unique (S)-adenosyl-L-methionine (SAM)-dependent N-methyltransferase, catalyzing transfer of a methyl group from SAM to the amino group of coclaurine. The model revealed a mixed α/β structure comprising seven twisted β-strands surrounded by twelve α-helices. Sequence comparisons and the model indicate an N-terminal catalytic Core domain and a C-terminal domain, of which the latter forms a pocket for coclaurine. An additional binding pocket for SAM is connected to the coclaurine binding pocket by a small opening. CNMT activity is proposed to follow an SN2-type mechanism as observed for a similarly conformed enzyme. Residues predicted for the methyl transfer reaction are Tyr79 and Glu96, which are conserved in the sequence from A. fimbriata and in homologous N-methyltransferases. The isolated CNMT is the first to be investigated from any basal angiosperm.

PREMISE OF THE STUDY We present the first plastome phylogeny encompassing all 77 monocot families, estimate branch support, and infer monocot-wide divergence times and rates of species diversification. METHODS We conducted maximum likelihood analyses of phylogeny and BAMM studies of diversification rates based on 77 plastid genes across 545 monocots and 22 outgroups. We quantified how branch support and ascertainment vary with gene number, branch length, and branch depth. KEY RESULTS Phylogenomic analyses shift the placement of 16 families in relation to earlier studies based on four plastid genes, add seven families, date the divergence between monocots and eudicots+Ceratophyllum at 136 Mya, successfully place all mycoheterotrophic taxa examined, and support recognizing Taccaceae and Thismiaceae as separate families and Arecales and Dasypogonales as separate orders. Only 45% of interfamilial divergences occurred after the Cretaceous. Net species diversification underwent four large-scale accelerations in PACMAD-BOP Poaceae, Asparagales sister to Doryanthaceae, Orchidoideae-Epidendroideae, and Araceae sister to Lemnoideae, each associated with specific ecological/morphological shifts. Branch ascertainment and support across monocots increase with gene number and branch length, and decrease with relative branch depth. Analysis of entire plastomes in Zingiberales quantifies the importance of non-coding regions in identifying and supporting short, deep branches. CONCLUSIONS We provide the first resolved, well-supported monocot phylogeny and timeline spanning all families, and quantify the significant contribution of plastome-scale data to resolving short, deep branches. We outline a new functional model for the evolution of monocots and their diagnostic morphological traits from submersed aquatic ancestors, supported by convergent evolution of many of these traits in aquatic Hydatellaceae (Nymphaeales).

GetOrganelle is a state-of-the-art toolkit to assemble accurate organelle genomes from NGS data. This toolkit recruit organelle-associated reads using a modified “baiting and iterative mapping” approach, conducts de novo assembly, filters and disentangles assembly graph, and produces all possible configurations of circular organelle genomes. For 50 published samples, we reassembled the circular plastome in 47 samples using GetOrganelle, but only in 12 samples using NOVOPlasty. In comparison with published/NOVOPlasty plastomes, we demonstrated that GetOrganelle assemblies are more accurate. Moreover, we assembled complete mitogenomes of fungi and animals using GetOrganelle. GetOrganelle is freely released under a GPL-3 license (https://github.com/Kinggerm/GetOrganelle).

3 PREMISE OF THE STUDY: Understanding fern (monilophyte) phylogeny and its evolutionary timescale is critical for broad investigations of the evolution of land plants, and for providing the point of comparison necessary for studying the evolution of the fern sister group, seed plants. Molecular phylogenetic investigations have revolutionized our understanding of fern phylogeny, however, to date, these studies have relied almost exclusively on plastid data. METHODS: Here we take a curated phylogenomics approach to infer the fibroad fern phylogeny from multiple nuclear loci, by combining broad taxon sampling (73 ferns and 12 outgroup species) with focused character sampling (25 loci comprising 35 877 bp), along with rigorous alignment, orthology inference and model selection. KEY RESULTS: Our phylogeny corroborates some earlier inferences and provides novel insights; in particular, we fi nd strong support for Equisetales as sister to the rest of ferns, Marattiales as sister to leptosporangiate ferns, and Dennstaedtiaceae as sister to the eupolypods. Our divergence-time analyses reveal that divergences among the extant fern orders all occurred prior to ~200 MYA. Finally, our species-tree inferences are congruent with analyses of concat- enated data, but generally with lower support. Those cases where species-tree support values are higher than expected involve relationships that have been supported by smaller plastid datasets, suggesting that deep coalescence may be reducing support from the concatenated nuclear data. CONCLUSIONS: Our study demonstrates the utility of a curated phylogenomics approach to inferring fern phylogeny, and highlights the need to consider underlying data characteristics, along with data quantity, in phylogenetic studies.

Coevolutionary interactions are thought to have spurred the evolution of key innovations and driven the diversification of much of life on Earth. However, the genetic and evolutionary basis of the innovations that facilitate such interactions remains poorly understood. We examined the coevolutionary interactions between plants (Brassicales) and butterflies (Pieridae), and uncovered evidence for an escalating evolutionary arms-race. Although gradual changes in trait complexity appear to have been facilitated by allelic turnover, key innovations are associated with gene and genome duplications. Furthermore, we show that the origins of both chemical defenses and of molecular counter adaptations were associated with shifts in diversification rates during the arms-race. These findings provide an important connection between the origins of biodiversity, coevolution, and the role of gene and genome duplications as a substrate for novel traits.

Strigolactones are plant hormones with multiple functions, including regulating various aspects of plant architecture such as shoot branching, facilitating the colonization of plant roots by arbuscular mycorrhizal fungi, and acting as seed germination stimulants for certain parasitic plants of the family Orobanchaceae. The obligate parasitic species Phelipanche aegyptiaca and Striga hermonthica require strigolactones for germination, while the facultative parasite Triphysaria versicolor does not. It has been hypothesized that P. aegyptiaca and S. hermonthica would have undergone evolutionary loss of strigolactone biosynthesis as a part of their mechanism to enable specific detection of exogenous strigolactones. We analyzed the transcriptomes of P. aegyptiaca, S. hermonthica and T. versicolor and identified genes known to act in strigolactone synthesis (D27, CCD7, CCD8, and MAX1), perception (MAX2 and D14) and transport (PDR12). These genes were then analyzed to assess likelihood of function. Transcripts of all strigolactone-related genes were found in P. aegyptiaca and S. hermonthica, and evidence points to their encoding functional proteins. Gene open reading frames were consistent with homologs from Arabidopsis and other strigolactone-producing plants, and all genes were expressed in parasite tissues. In general, the genes related to strigolactone synthesis and perception appeared to be evolving under codon-based selective constraints in strigolactone-dependent species. Bioassays of S. hermonthica root extracts indicated the presence of strigolactone class stimulants on germination of P. aegyptiaca seeds. Taken together, these results indicate that Phelipanche aegyptiaca and S. hermonthica have retained functional genes involved in strigolactone biosynthesis, suggesting that the parasites use both endogenous and exogenous strigolactones and have mechanisms to differentiate the two.