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.

Plant genome-scale resources are being generated at an increasing rate as sequencing technologies continue to improve and raw data costs continue to fall; however, the cost of downstream analyses remains large. This has resulted in a considerable range of genome assembly and annotation qualities across plant genomes due to their varying sizes, complexity, and the technology used for the assembly and annotation. To effectively work across genomes, researchers increasingly rely on comparative genomic approaches that integrate across plant community resources and data types. Such efforts have aided the genome annotation process and yielded novel insights into the evolutionary history of genomes and gene families, including complex non-model organisms. The essential tools to achieve these insights rely on gene family analysis at a genome-scale, but they are not well integrated for rapid analysis of new data, and the learning curve can be steep. Here we present PlantTribes2, a scalable, easily accessible, highly customizable, and broadly applicable gene family analysis framework with multiple entry points including user provided data. It uses objective classifications of annotated protein sequences from existing, high-quality plant genomes for comparative and evolutionary studies. PlantTribes2 can improve transcript models and then sort them, either genome-scale annotations or individual gene coding sequences, into pre-computed orthologous gene family clusters with rich functional annotation information. Then, for gene families of interest, PlantTribes2 performs downstream analyses and customizable visualizations including, (1) multiple sequence alignment, (2) gene family phylogeny, (3) estimation of synonymous and non-synonymous substitution rates among homologous sequences, and (4) inference of large-scale duplication events. We give examples of PlantTribes2 applications in functional genomic studies of economically important plant families, namely transcriptomics in the weedy Orobanchaceae and a core orthogroup analysis (CROG) in Rosaceae. PlantTribes2 is freely available for use within the main public Galaxy instance and can be downloaded from GitHub or Bioconda. Importantly, PlantTribes2 can be readily adapted for use with genomic and transcriptomic data from any kind of organism.

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.

The rapid development of sequencing technologies has led to a deeper understanding of plant genomes. However, direct experimental evidence connecting genes to important agronomic traits is still lacking in most non-model plants. 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). An efficient way to identify genetic factors for important traits in non-model organisms can be to transfer knowledge across genomes. However, major obstacles exist, including complex evolutionary histories and variable quality and content of publicly available plant genomes. As researchers aim to link genes to traits of interest, these challenges can impede the transfer of experimental evidence across plant species, namely in the curation of high-quality, high-confidence gene models in an evolutionary context. Here we present a workflow using a collection of bioinformatic tools for the curation of deeply conserved gene families of interest across plant genomes. To study gene families involved in tree architecture in European pear and other rosaceous species, we used our workflow, plus a draft genome assembly and high-quality annotation of a second P. communis cultivar, ‘d’Anjou.’ Our comparative gene family approach revealed significant issues with the most recent ‘Bartlett’ genome - primarily thousands of missing genes due to methodological bias. After correcting assembly errors on a global scale in the ‘Bartlett’ genome, we used our workflow for targeted improvement of our genes of interest in both P. communis genomes, thus laying the groundwork for future functional studies in pear tree architecture. Further, our global gene family classification of 15 genomes across 6 genera provides a valuable and previously unavailable resource for the Rosaceae research community. With it, orthologs and other gene family members can be easily identified across any of the classified genomes. Importantly, our workflow can be easily adopted for any other plant genomes and gene families of interest.

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.

Disentangling the phylogenetic relationships of taxonomically complex plant groups is often mired by challenges associated with recent speciation, hybridization, complex mating systems, and polyploidy. Here, we perform a phylogenomic analysis of eyebrights (Euphrasia), a group renowned for taxonomic complexity, with the aim of documenting the extent of phylogenetic discordance at both deep and at shallow phylogenetic scales. We generate whole-genome sequencing data and integrate this with prior genomic data to perform a comprehensive analysis of nuclear genomic, nuclear ribosomal (nrDNA), and complete plastid genomes from 57 individuals representing 36 Euphrasia species. The species tree analysis of 3,454 conserved nuclear scaffolds (46 Mb) reveals that at shallow phylogenetic scales postglacial colonization of North Western Europe occurred in multiple waves from discrete source populations, with most species not being monophyletic, and instead combining genomic variants from across clades. At a deeper phylogenetic scale, the Euphrasia phylogeny is structured by geography and ploidy, and partially by taxonomy. Comparative analyses show Southern Hemisphere tetraploids include a distinct subgenome indicative of independent polyploidy events from Northern Hemisphere taxa. In contrast to the nuclear genome analyses, the plastid genome phylogeny reveals limited geographic structure, while the nrDNA phylogeny is informative of some geographic and taxonomic affinities but more thorough phylogenetic inference is impeded by the retention of ancestral polymorphisms in the polyploids. Overall our results reveal extensive phylogenetic discordance at both deeper and shallower nodes, with broad-scale geographic structure of genomic variation but a lack of definitive taxonomic signal. This suggests that Euphrasia species either have polytopic origins or are maintained by narrow genomic regions in the face of extensive homogenizing gene flow. Moreover, these results suggest genome skimming will not be an effective extended barcode to identify species in groups such as Euphrasia, or many other postglacial species groups.

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.

Disentangling the phylogenetic relationships of taxonomically complex plant groups is often mired by challenges associated with recent speciation, hybridization, complex mating systems, and polyploidy. Here, we perform a phylogenomic analysis of eyebrights (

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.

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.

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.

Estimating maturity in pome fruits is a critical task that directs virtually all postharvest supply chain decisions. This is especially important for European pear (Pyrus communis)cultivars because losses due to spoilage and senescence must be minimized while ensuring proper ripening capacity is achieved (in part by satisfying a fruit chilling requirement). Reliable methods are lacking for accurate estimation of pear fruit maturity, and because ripening is maturity dependent it makes predicting ripening capacity a challenge. In this study of the European pear cultivar ‘d’Anjou’, we sorted fruit at harvest based upon on-tree fruit position to build contrasts of maturity. Our sorting scheme showed clear contrasts of maturity between canopy positions, yet there was substantial overlap in the distribution of values for the index of absorbance difference (IAD), a non-destructive spectroscopic measurement that has been used as a proxy for pome fruit maturity. This presented an opportunity to explore a contrast of maturity that was more subtle thanIADcould differentiate, and thus guided our subsequent transcriptome analysis of tissue samples taken at harvest and during storage. Using a novel approach that tests for condition-specific differences of co-expressed genes, we discovered genes with a phased character that mirrored our sorting scheme. The expression patterns of these genes are associated with fruit quality and ripening differences across the experiment. Functional profiles of these co-expressed genes are concordant with previous findings, and also offer new clues, and thus hypotheses, about genes involved in pear fruit quality, maturity, and ripening. This work may lead to new tools for enhanced postharvest management based on activity of gene co-expression modules, rather than individual genes. Further, our results indicate that modules may have utility within specific windows of time during postharvest management of ‘d’Anjou’ pear.

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.

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.

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.

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.

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.