Footprints of parasitism in the genome of the parasitic flowering plant Cuscuta campestris

Horizontally transferred mitochondrial DNA tracts become circular by microhomology-mediated repair pathways

The holoparasitic plant Lophophytum mirabile exhibits remarkable levels of mitochondrial horizontal gene transfer (HGT). Gathering comparative data from other individuals and host plants can provide insights into the HGT process. We sequenced the mitochondrial genome (mtDNA) from individuals of two species of Lophophytum and from mimosoid hosts. We applied a stringent phylogenomic approach to elucidate the origin of the whole mtDNAs, estimate the timing of the transfers, and understand the molecular mechanisms involved. Ancestral and recent HGT events replaced and enlarged the multichromosomal mtDNA of Lophophytum spp., with the foreign DNA ascending to 74%. A total of 14 foreign mitochondrial chromosomes originated from continuous regions in the host mtDNA flanked by short direct repeats. These foreign tracts are circularized by microhomology-mediated repair pathways and replicate independently until they are lost or they eventually recombine with other chromosomes. The foreign noncoding chromosomes are variably present in the population and likely evolve by genetic drift. We present the 'circle-mediated HGT' model in which foreign mitochondrial DNA tracts become circular and are maintained as plasmid-like molecules. This model challenges the conventional belief that foreign DNA must be integrated into the recipient genome for successful HGT.

Chromosome-level genome assembly of Pedicularis kansuensis illuminates genome evolution of facultative parasitic plant

Parasitic plants have a heterotrophic lifestyle, in which they withdraw all or part of their nutrients from their host through the haustorium. Despite the release of many draft genomes of parasitic plants, the genome evolution related to the parasitism feature of facultative parasites remains largely unknown. In this study, we present a high-quality chromosomal-level genome assembly for the facultative parasite Pedicularis kansuensis (Orobanchaceae), which invades both legume and grass host species in degraded grasslands on the Qinghai-Tibet Plateau. This species has the largest genome size compared with other parasitic species, and expansions of long terminal repeat retrotransposons accounting for 62.37% of the assembly greatly contributed to the genome size expansion of this species. A total of 42,782 genes were annotated, and the patterns of gene loss in P. kansuensis differed from other parasitic species. We also found many mobile mRNAs between P. kansuensis and one of its host species, but these mobile mRNAs could not compensate for the functional losses of missing genes in P. kansuensis. In addition, we identified nine horizontal gene transfer (HGT) events from rosids and monocots, as well as one single-gene duplication events from HGT genes, which differ distinctly from that of other parasitic species. Furthermore, we found evidence for HGT through transferring genomic fragments from phylogenetically remote host species. Taken together, these findings provide genomic insights into the evolution of facultative parasites and broaden our understanding of the diversified genome evolution in parasitic plants and the molecular mechanisms of plant parasitism.

Horizontal gene transfer in eukaryotes: aligning theory with data

Horizontal gene transfer (HGT), or lateral gene transfer, is the non-sexual movement of genetic information between genomes. It has played a pronounced part in bacterial and archaeal evolution, but its role in eukaryotes is less clear. Behaviours unique to eukaryotic cells - phagocytosis and endosymbiosis - have been proposed to increase the frequency of HGT, but nuclear genomes encode fewer HGTs than bacteria and archaea. Here, I review the existing theory in the context of the growing body of data on HGT in eukaryotes, which suggests that any increased chance of acquiring new genes through phagocytosis and endosymbiosis is offset by a reduced need for these genes in eukaryotes, because selection in most eukaryotes operates on variation not readily generated by HGT.

Parasitic dodder expresses an arsenal of secreted cellulases with multi-substrate specificity during host invasion

Cuscuta campestris is a common and problematic parasitic plant which relies on haustoria to connect to and siphon nutrients from host plants. Glycoside hydrolase family 9 (GH9) cellulases (EC 3.2.1.4) play critical roles in plant cell wall biosynthesis and disassembly, but their roles during Cuscuta host invasion remains underexplored. In this study, we identified 22 full-length GH9 cellulase genes in C. campestris genome, which encoded fifteen secreted and seven membrane-anchored cellulases that showed distinct phylogenetic relationships. Expression profiles suggested that some of the genes are involved in biosynthesis and remodeling of the parasite's cell wall during haustoriogenesis, while other genes encoding secreted B- and C-type cellulases are tentatively associated with degrading host cell walls during invasion. Transcriptomic data in a host-free system and in the presence of susceptible or partially resistant tomato hosts, showed for especially GH9B7, GH9B11 and GH9B12 a shift in expression profiles in the presence of hosts, being more highly expressed during host attachment, indicating that Cuscuta can tune cellulase expression in response to a host. Functional analyses of recombinant B- and C-type cellulases showed endoglucanase activities over wide pH and temperature conditions, and activities towards multiple cellulose and hemicellulose substrates. These findings improve our understanding of host cell wall disassembly by Cuscuta, and cellulase activity towards broad substrate range potentially explain its wide host range. This is the first study to provide a broad biochemical insight into Cuscuta GH9 cellulases, which based on our study may have potential applications in industrial bioprocessing.

Variations and reduction of plastome are associated with the evolution of parasitism in Convolvulaceae

Parasitic lifestyle can often relax the constraint on the plastome, leading to gene pseudogenization and loss, and resulting in diverse genomic structures and rampant genome degradation. Although several plastomes of parasitic Cuscuta have  been reported, the evolution of parasitism in the family Convolvulaceae which is linked to structural variations and reduction of plastome has not been well investigated. In this study, we assembled and collected 40 plastid genomes belonging to 23 species representing four subgenera of Cuscuta and ten species of autotrophic Convolvulaceae. Our findings revealed nine types of structural variations and six types of inverted repeat (IR) boundary variations in the plastome of Convolvulaceae spp. These structural variations were associated with the shift of parasitic lifestyle, and IR boundary shift, as well as the abundance of long repeats. Overall, the degradation of Cuscuta plastome proceeded gradually, with one clade exhibiting an accelerated degradation rate. We observed five stages of gene loss in Cuscuta, including NAD(P)H complex → PEP complex → Photosynthesis-related → Ribosomal protein subunits → ATP synthase complex. Based on our results, we speculated that the shift of parasitic lifestyle in early divergent time promoted relaxed selection on plastomes, leading to the accumulation of microvariations, which ultimately resulted in the plastome reduction. This study provides new evidence towards a better understanding of plastomic evolution, variation, and reduction in the genus Cuscuta.

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

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

Comparative secretome analysis of Striga and Cuscuta species identifies candidate virulence factors for two evolutionarily independent parasitic plant lineages

BACKGROUND

Many parasitic plants of the genera Striga and Cuscuta inflict huge agricultural damage worldwide. To form and maintain a connection with a host plant, parasitic plants deploy virulence factors (VFs) that interact with host biology. They possess a secretome that represents the complement of proteins secreted from cells and like other plant parasites such as fungi, bacteria or nematodes, some secreted proteins represent VFs crucial to successful host colonisation. Understanding the genome-wide complement of putative secreted proteins from parasitic plants, and their expression during host invasion, will advance understanding of virulence mechanisms used by parasitic plants to suppress/evade host immune responses and to establish and maintain a parasite-host interaction.

RESULTS

We conducted a comparative analysis of the secretomes of root (Striga spp.) and shoot (Cuscuta spp.) parasitic plants, to enable prediction of candidate VFs. Using orthogroup clustering and protein domain analyses we identified gene families/functional annotations common to both Striga and Cuscuta species that were not present in their closest non-parasitic relatives (e.g. strictosidine synthase like enzymes), or specific to either the Striga or Cuscuta secretomes. For example, Striga secretomes were strongly associated with 'PAR1' protein domains. These were rare in the Cuscuta secretomes but an abundance of 'GMC oxidoreductase' domains were found, that were not present in the Striga secretomes. We then conducted transcriptional profiling of genes encoding putatively secreted proteins for the most agriculturally damaging root parasitic weed of cereals, S. hermonthica. A significant portion of the Striga-specific secretome set was differentially expressed during parasitism, which we probed further to identify genes following a 'wave-like' expression pattern peaking in the early penetration stage of infection. We identified 39 genes encoding putative VFs with functions such as cell wall modification, immune suppression, protease, kinase, or peroxidase activities, that are excellent candidates for future functional studies.

CONCLUSIONS

Our study represents a comprehensive secretome analysis among parasitic plants and revealed both similarities and differences in candidate VFs between Striga and Cuscuta species. This knowledge is crucial for the development of new management strategies and delaying the evolution of virulence in parasitic weeds.

Rising from the shadows: Selective foraging in model shoot parasitic plants

Despite being sessile, plants nonetheless forage for resources by modulating their growth. Adaptative foraging in response to changes in resource availability and presence of neighbours has strong implications for performance and fitness. It is an even more pressing issue for parasitic plants, which draw resources directly from other plants. Indeed, parasitic plants were demonstrated over the years to direct their growth towards preferred hosts and invest resources in parasitism relative to host quality. In contrast to root parasites that rely mostly on chemical cues, some shoot parasites seem to profit from the ability to integrate different types of abiotic and biotic cues. While significant progress in this field has been made recently, there are still many open questions regarding the molecular perception and the integration of diverse signalling pathways under different ecological contexts. Addressing how different cues are integrated in parasitic plants will be important when unravelling variations in plant interaction pathways, and essential to predict the spread of parasites in natural and agricultural environments. In this review, we discuss this with a focus on Cuscuta species as an emerging parasitic model, and provide research perspectives based on the recent advances in the topic and plant-plant interactions in general.

Empirical evidence for metabolic drift in plant and algal lipid biosynthesis pathways

Metabolic pathway drift has been formulated as a general principle to help in the interpretation of comparative analyses between biosynthesis pathways. Indeed, such analyses often indicate substantial differences, even in widespread pathways that are sometimes believed to be conserved. Here, our purpose is to check how much this interpretation fits to empirical data gathered in the field of plant and algal biosynthesis pathways. After examining several examples representative of the diversity of lipid biosynthesis pathways, we explain why it is important to compare closely related species to gain a better understanding of this phenomenon. Furthermore, this comparative approach brings us to the question of how much biotic interactions are responsible for shaping this metabolic plasticity. We end up introducing some model systems that may be promising for further exploration of this question.

Cuscuta campestris fine-tunes gene expression during haustoriogenesis as an adaptation to different hosts

The Cuscuta genus comprises obligate parasitic plants that have an unusually wide host range. Whether Cuscuta uses different infection strategies for different hosts or whether the infection strategy is mechanistically and enzymatically conserved remains unknown. To address this, we investigated molecular events during the interaction between field dodder (Cuscuta campestris) and two host species of the Solanum genus that are known to react differently to parasitic infection. We found that host gene induction, particularly of cell wall fortifying genes, coincided with a differential induction of genes for cell wall degradation in the parasite in the cultivated tomato (Solanum lycopersicum) but not in a wild relative (Solanum pennellii). This indicates that the parasite can adjust its gene expression in response to its host. This idea was supported by the increased expression of C. campestris genes encoding an endo-β-1,4-mannanase in response to exposure of the parasite to purified mono- and polysaccharides in a host-independent infection system. Our results suggest multiple key roles of the host cell wall in determining the outcome of an infection attempt.

A roadmap of haustorium morphogenesis in parasitic plants

Parasitic plants invade their host through their invasive organ, the haustorium. This organ connects to the vasculature of the host roots and hijacks water and nutrients. Although parasitism has evolved independently in plants, haustoria formation follows a similar mechanism throughout different plant species, highlighting the developmental plasticity of plant tissues. Here, we compare three types of haustoria formed by the root and shoot in the plant parasites Striga and Cuscuta. We discuss mechanisms underlying the interactions with their hosts and how different approaches have contributed to major understanding of haustoria formation and host invasion. We also illustrate the role of auxin and cytokinin in controlling this process.

Invited Review Beyond parasitic convergence: unravelling the evolution of the organellar genomes in holoparasites

BACKGROUND

The molecular evolution of organellar genomes in angiosperms has been studied extensively, with some lineages, such as parasitic ones, displaying unique characteristics. Parasitism has emerged 12 times independently in angiosperm evolution. Holoparasitism is the most severe form of parasitism, and is found in ~10 % of parasitic angiosperms. Although a few holoparasitic species have been examined at the molecular level, most reports involve plastomes instead of mitogenomes. Parasitic plants establish vascular connections with their hosts through haustoria to obtain water and nutrients, which facilitates the exchange of genetic information, making them more susceptible to horizontal gene transfer (HGT). HGT is more prevalent in the mitochondria than in the chloroplast or nuclear compartments.

SCOPE

This review summarizes current knowledge on the plastid and mitochondrial genomes of holoparasitic angiosperms, compares the genomic features across the different lineages, and discusses their convergent evolutionary trajectories and distinctive features. We focused on Balanophoraceae (Santalales), which exhibits extraordinary traits in both their organelles.

CONCLUSIONS

Apart from morphological similarities, plastid genomes of holoparasitic plants also display other convergent features, such as rampant gene loss, biased nucleotide composition and accelerated evolutionary rates. In addition, the plastomes of Balanophoraceae have extremely low GC and gene content, and two unexpected changes in the genetic code. Limited data on the mitochondrial genomes of holoparasitic plants preclude thorough comparisons. Nonetheless, no obvious genomic features distinguish them from the mitochondria of free-living angiosperms, except for a higher incidence of HGT. HGT appears to be predominant in holoparasitic angiosperms with a long-lasting endophytic stage. Among the Balanophoraceae, mitochondrial genomes exhibit disparate evolutionary paths with notable levels of heteroplasmy in Rhopalocnemis and unprecedented levels of HGT in Lophophytum. Despite their differences, these Balanophoraceae share a multichromosomal mitogenome, a feature also found in a few free-living angiosperms.

Balanophora genomes display massively convergent evolution with other extreme holoparasites and provide novel insights into parasite-host interactions

Parasitic plants have evolved to be subtly or severely dependent on host plants to complete their life cycle. To provide new insights into the biology of parasitic plants in general, we assembled genomes for members of the sandalwood order Santalales, including a stem hemiparasite (Scurrula) and two highly modified root holoparasites (Balanophora) that possess chimaeric host-parasite tubers. Comprehensive genome comparisons reveal that hemiparasitic Scurrula has experienced a relatively minor degree of gene loss compared with autotrophic plants, consistent with its moderate degree of parasitism. Nonetheless, patterns of gene loss appear to be substantially divergent across distantly related lineages of hemiparasites. In contrast, Balanophora has experienced substantial gene loss for the same sets of genes as an independently evolved holoparasite lineage, the endoparasitic Sapria (Malpighiales), and the two holoparasite lineages experienced convergent contraction of large gene families through loss of paralogues. This unprecedented convergence supports the idea that despite their extreme and strikingly divergent life histories and morphology, the evolution of these and other holoparasitic lineages can be shaped by highly predictable modes of genome reduction. We observe substantial evidence of relaxed selection in retained genes for both hemi- and holoparasitic species. Transcriptome data also document unusual and novel interactions between Balanophora and host plants at the host-parasite tuber interface tissues, with evidence of mRNA exchange, substantial and active hormone exchange and immune responses in parasite and host.

LysM-mediated signaling in Marchantia polymorpha highlights the conservation of pattern-triggered immunity in land plants

Pattern-recognition receptor (PRR)-triggered immunity (PTI) wards off a wide range of pathogenic microbes, playing a pivotal role in angiosperms. The model liverwort Marchantia polymorpha triggers defense-related gene expression upon sensing components of bacterial and fungal extracts, suggesting the existence of PTI in this plant model. However, the molecular components of the putative PTI in M. polymorpha and the significance of PTI in bryophytes have not yet been described. We here show that M. polymorpha has four lysin motif (LysM)-domain-containing receptor homologs, two of which, LysM-receptor-like kinase (LYK) MpLYK1 and LYK-related (LYR) MpLYR, are responsible for sensing chitin and peptidoglycan fragments, triggering a series of characteristic immune responses. Comprehensive phosphoproteomic analysis of M. polymorpha in response to chitin treatment identified regulatory proteins that potentially shape LysM-mediated PTI. The identified proteins included homologs of well-described PTI components in angiosperms as well as proteins whose roles in PTI are not yet determined, including the blue-light receptor phototropin MpPHOT. We revealed that MpPHOT is required for negative feedback of defense-related gene expression during PTI. Taken together, this study outlines the basic framework of LysM-mediated PTI in M. polymorpha and highlights conserved elements and new aspects of pattern-triggered immunity in land plants.

The Mitochondrial Genome of the Holoparasitic Plant Thonningia sanguinea Provides Insights into the Evolution of the Multichromosomal Structure

Multichromosomal mitochondrial genome (mitogenome) structures have repeatedly evolved in many lineages of angiosperms. However, the underlying mechanism remains unclear. The mitogenomes of three genera of Balanophoraceae, namely Lophophytum, Ombrophytum, and Rhopalocnemis, have already been sequenced and assembled, all showing a highly multichromosomal structure, albeit with different genome and chromosome sizes. It is expected that characterization of additional lineages of this family may expand the knowledge of mitogenome diversity and provide insights into the evolution of the plant mitogenome structure and size. Here, we assembled and characterized the mitogenome of Thonningia sanguinea, which, together with Balanophora, forms a clade sister to the clade comprising Lophophytum, Ombrophytum, and Rhopalocnemis. The mitogenome of T. sanguinea possesses a multichromosomal structure of 18 circular chromosomes of 8.7-19.2 kb, with a total size of 246,247 bp. There are very limited shared regions and poor chromosomal correspondence between T. sanguinea and other Balanophoraceae species, suggesting frequent rearrangements and rapid sequence turnover. Numerous medium- and small-sized repeats were identified in the T. sanguinea mitogenome; however, no repeat-mediated recombination was detected, which was verified by Illumina reads mapping and PCR experiments. Intraspecific mitogenome variations in T. sanguinea are mostly insertions and deletions, some of which can lead to degradation of perfect repeats in one or two accessions. Based on the mitogenome features of T. sanguinea, we propose a mechanism to explain the evolution of a multichromosomal mitogenome from a master circle, which involves mutation in organellar DNA replication, recombination and repair genes, decrease of recombination, and repeat degradation.

Anatomical observation and transcriptome analysis of branch-twisted mutations in Chinese jujube

BACKGROUND

Plant organs grow in a certain direction and organ twisted growth, a rare and distinctive trait, is associated with internal structure changes and special genes. The twisted branch mutant of Chinese jujube jujube, an important fruit tree native to China and introduced to nearly 50 countries, provides new typical materials for exploration of plant twisted growth.

RESULTS

In this study, the cytological characteristics and related genes of twisted branches in Chinese jujube were revealed by microscopy observation and transcriptome analysis. The unique coexistence of primary and secondary structures appeared in the twisted parts of branches, and special structures such as collateral bundle, cortical bundles, and internal phloem were formed. Ninety differentially expressed genes of 'Dongzao' and its twisted mutant were observed, in which ZjTBL43, ZjFLA11, ZjFLA12 and ZjIQD1 were selected as candidate genes. ZjTBL43 was homologous to AtTBL43 in Arabidopsis, which was involved in the synthesis and deposition of cellular secondary wall cellulose. The attbl43 mutant showed significant inflorescence stem bending growth. The transgenic lines of attbl43 with overexpression of ZjTBL43 were phenotypically normal.The branch twisted growth may be caused by mutations in ZjTBL43 in Chinese jujube. AtIQD10, AtFLA11 and AtFLA12 were homologous to ZjIQD1, ZjFLA11 and ZjFLA12, respectively. However, the phenotype of their function defect mutants was normal.

CONCLUSION

In summary, these findings will provide new insights into the plant organ twisted growth and a reference for investigation of controlling mechanisms of plant growth direction.

Host shift promotes divergent evolution between closely related holoparasitic species

Distinct hosts have been hypothesized to possess the potential for affecting species differentiation and genome evolution of parasitic organisms. However, what host shift history is experienced by the closely related parasites and whether disparate evolution of their genomes occur remain largely unknown. Here, we screened horizontal gene transfer (HGT) events in a pair of sister species of holoparasitic Boschniakia (Orobanchaceae) having obligate hosts from distinct families to recall the former host-parasite associations and performed a comparative analysis to investigate the difference of their organelle genomes. Except those from the current hosts (Ericaceae and Betulaceae), we identified a number of HGTs from Rosaceae supporting the occurrence of unexpected ancient host shifts. Different hosts transfer functional genes which changed nuclear genomes of this sister species. Likewise, different donors transferred sequences to their mitogenomes, which vary in size due to foreign and repetitive elements rather than other factors found in other parasites. The plastomes are both severely reduced, and the degree of difference in reduction syndrome reaches the intergeneric level. Our findings provide new insights into the genome evolution of parasites adapting to different hosts and extend the mechanism of host shift promoting species differentiation to parasitic plant lineages.

Trans-species microRNA loci in the parasitic plant Cuscuta campestris have a U6-like snRNA promoter

Small regulatory RNAs can move between organisms and regulate gene expression in the recipient. Whether the trans-species small RNAs being exported are distinguished from the normal endogenous small RNAs of the source organism is not known. The parasitic plant Cuscuta campestris (dodder) produces many microRNAs that specifically accumulate at the host-parasite interface, several of which have trans-species activity. We found that induction of C. campestris interface-induced microRNAs is similar regardless of host species and occurs in C. campestris haustoria produced in the absence of any host. The loci-encoding C. campestris interface-induced microRNAs are distinguished by a common cis-regulatory element. This element is identical to a conserved upstream sequence element (USE) used by plant small nuclear RNA loci. The properties of the interface-induced microRNA primary transcripts strongly suggest that they are produced via U6-like transcription by RNA polymerase III. The USE promotes accumulation of interface-induced miRNAs (IIMs) in a heterologous system. This promoter element distinguishes C. campestris IIM loci from other plant small RNAs. Our data suggest that C. campestris IIMs are produced in a manner distinct from canonical miRNAs. All confirmed C. campestris microRNAs with documented trans-species activity are interface-induced and possess these features. We speculate that RNA polymerase III transcription of IIMs may allow these miRNAs to be exported to hosts.

The plant vampire diaries: a historic perspective on Cuscuta research

The angiosperm genus Cuscuta lives as an almost achlorophyllous root- and leafless holoparasite and has therefore occupied scientists for more than a century. The 'evolution' of Cuscuta research started with early studies that established the phylogenetic framework for this unusual genus. It continued to produce groundbreaking cytological, morphological, and physiological insight throughout the second half of the 20th century and culminated in the last two decades in exciting discoveries regarding the molecular basis of Cuscuta parasitism that were facilitated by the modern 'omics' tools and traceable fluorescent marker technologies of the 21st century. This review will show how present activities are inspired by those past breakthroughs. It will describe significant milestones and recurring themes of Cuscuta research and connect these to the remaining as well as newly evolving questions and future directions in this research field that is expected to sustain its strong growth in the future.

Rethinking convergence in plant parasitism through the lens of molecular and population genetic processes

The autotrophic lifestyle of photosynthetic plants has profoundly shaped their body plan, physiology, and gene repertoire. Shifts to parasitism and heterotrophy have evolved at least 12 times in more than 4000 species, and this transition has consequently left major evolutionary footprints among these parasitic lineages. Features that are otherwise rare at the molecular level and beyond have evolved repetitively, including reduced vegetative bodies, carrion-mimicking during reproduction, and the incorporation of alien genetic material. Here, I propose an integrated conceptual model, referred to as the funnel model, to define the general evolutionary trajectory of parasitic plants and provide a mechanistic explanation for their convergent evolution. This model connects our empirical understanding of gene regulatory networks in flowering plants with classical theories of molecular and population genetics. It emphasizes that the cascading effects brought about by the loss of photosynthesis may be a major force constraining the physiological capacity of parasitic plants and shaping their genomic landscapes. Here I review recent studies on the anatomy, physiology, and genetics of parasitic plants that lend support to this photosynthesis-centered funnel model. Focusing on nonphotosynthetic holoparasites, I elucidate how they may inevitably reach an evolutionary terminal status (i.e., extinction) and highlight the utility of a general, explicitly described and falsifiable model for future studies of parasitic plants.

Weed genomics: yielding insights into the genetics of weedy traits for crop improvement

Weeds cause tremendous economic and ecological damage worldwide. The number of genomes established for weed species has sharply increased during the recent decade, with some 26 weed species having been sequenced and de novo genomes assembled. These genomes range from 270 Mb (Barbarea vulgaris) to almost 4.4 Gb (Aegilops tauschii). Importantly, chromosome-level assemblies are now available for 17 of these 26 species, and genomic investigations on weed populations have been conducted in at least 12 species. The resulting genomic data have greatly facilitated studies of weed management and biology, especially origin and evolution. Available weed genomes have indeed revealed valuable weed-derived genetic materials for crop improvement. In this review, we summarize the recent progress made in weed genomics and provide a perspective for further exploitation in this emerging field.

Genomic and Epigenomic Mechanisms of the Interaction between Parasitic and Host Plants

Parasitic plants extract nutrients from the other plants to finish their life cycle and reproduce. The control of parasitic weeds is notoriously difficult due to their tight physical association and their close biological relationship to their hosts. Parasitic plants differ in their susceptible host ranges, and the host species differ in their susceptibility to parasitic plants. Current data show that adaptations of parasitic plants to various hosts are largely genetically determined. However, multiple cases of rapid adaptation in genetically homogenous parasitic weed populations to new hosts strongly suggest the involvement of epigenetic mechanisms. Recent progress in genome-wide analyses of gene expression and epigenetic features revealed many new molecular details of the parasitic plants' interactions with their host plants. The experimental data obtained in the last several years show that multiple common features have independently evolved in different lines of the parasitic plants. In this review we discuss the most interesting new details in the interaction between parasitic and host plants.

Assembly, annotation, and comparative analysis of Ipomoea chloroplast genomes provide insights into the parasitic characteristics of Cuscuta species

In the Convolvulaceae family, around 1650 species belonging to 60 genera are widely distributed globally, mainly in the tropical and subtropical regions of America and Asia. Although a series of chloroplast genomes in Convolvulaceae were reported and investigated, the evolutionary and genetic relationships among the chloroplast genomes of the Convolvulaceae family have not been extensively elucidated till now. In this study, we first reported the complete chloroplast genome sequence of Ipomoea pes-caprae, a widely distributed coastal plant with medical values. The chloroplast genome of I. pes-caprae is 161667 bp in length, and the GC content is 37.56%. The chloroplastic DNA molecule of I. pes-caprae is a circular structure composed of LSC (large-single-copy), SSC (small-single-copy), and IR (inverted repeat) regions, with the size of the three regions being 88210 bp, 12117 bp, and 30670 bp, respectively. The chloroplast genome of I. pes-caprae contains 141 genes, and 35 SSRs are identified in the chloroplast genome. Our research results provide important genomic information for the molecular phylogeny of I. pes-caprae. The Phylogenetic analysis of 28 Convolvulaceae chloroplast genomes showed that the relationship of I. pes-caprae with I. involucrata or I. obscura was much closer than that with other Convolvulaccae species. Further comparative analyses between the Ipomoea species and Cuscuta species revealed the mechanism underlying the formation of parasitic characteristics of Cuscuta species from the perspective of the chloroplast genome.

Phylotranscriptomic Analyses of Mycoheterotrophic Monocots Show a Continuum of Convergent Evolutionary Changes in Expressed Nuclear Genes From Three Independent Nonphotosynthetic Lineages

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.

The Current Developments in Medicinal Plant Genomics Enabled the Diversification of Secondary Metabolites' Biosynthesis

Medicinal plants produce important substrates for their adaptation and defenses against environmental factors and, at the same time, are used for traditional medicine and industrial additives. Plants have relatively little in the way of secondary metabolites via biosynthesis. Recently, the whole-genome sequencing of medicinal plants and the identification of secondary metabolite production were revolutionized by the rapid development and cheap cost of sequencing technology. Advances in functional genomics, such as transcriptomics, proteomics, and metabolomics, pave the way for discoveries in secondary metabolites and related key genes. The multi-omics approaches can offer tremendous insight into the variety, distribution, and development of biosynthetic gene clusters (BGCs). Although many reviews have reported on the plant and medicinal plant genome, chemistry, and pharmacology, there is no review giving a comprehensive report about the medicinal plant genome and multi-omics approaches to study the biosynthesis pathway of secondary metabolites. Here, we introduce the medicinal plant genome and the application of multi-omics tools for identifying genes related to the biosynthesis pathway of secondary metabolites. Moreover, we explore comparative genomics and polyploidy for gene family analysis in medicinal plants. This study promotes medicinal plant genomics, which contributes to the biosynthesis and screening of plant substrates and plant-based drugs and prompts the research efficiency of traditional medicine.

Cuscuta species: Model organisms for haustorium development in stem holoparasitic plants

Parasitic plants are notorious for causing serious agricultural losses in many countries. Specialized intrusive organs, haustoria, confer on parasitic plants the ability to acquire water and nutrients from their host plants. Investigating the mechanism involved in haustorium development not only reveals the fascinating mystery of how autotrophic plants evolved parasitism but also provides the foundation for developing more effective methods to control the agricultural damage caused by parasitic plants. Cuscuta species, also known as dodders, are one of the most well-known and widely spread stem holoparasitic plants. Although progress has been made recently in understanding the evolution and development of haustoria in root parasitic plants, more and more studies indicate that the behaviors between root and stem haustorium formation are distinct, and the mechanisms involved in the formation of these organs remain largely unknown. Unlike most endoparasites and root holoparasitic plants, which have high host-specificity and self- or kin-recognition to avoid forming haustoria on themselves or closely related species, auto-parasitism and hyper-parasitism are commonly observed among Cuscuta species. In this review, we summarize the current understanding of haustorium development in dodders and the unique characteristics of their parasitizing behaviors. We also outline the advantages of using Cuscuta species as model organisms for haustorium development in stem holoparasitic plants, the current unknown mysteries and limitations in the Cuscuta system, and potential future research directions to overcome these challenges.

Genome mining reveals abiotic stress resistance genes in plant genomes acquired from microbes via HGT

Colonization by beneficial microbes can enhance plant tolerance to abiotic stresses. However, there are still many unknown fields regarding the beneficial plant-microbe interactions. In this study, we have assessed the amount or impact of horizontal gene transfer (HGT)-derived genes in plants that have potentials to confer abiotic stress resistance. We have identified a total of 235 gene entries in fourteen high-quality plant genomes belonging to phyla Chlorophyta and Streptophyta that confer resistance against a wide range of abiotic pressures acquired from microbes through independent HGTs. These genes encode proteins contributed to toxic metal resistance (e.g., ChrA, CopA, CorA), osmotic and drought stress resistance (e.g., Na+/proline symporter, potassium/proton antiporter), acid resistance (e.g., PcxA, ArcA, YhdG), heat and cold stress resistance (e.g., DnaJ, Hsp20, CspA), oxidative stress resistance (e.g., GST, PoxA, glutaredoxin), DNA damage resistance (e.g., Rad25, Rad51, UvrD), and organic pollutant resistance (e.g., CytP450, laccase, CbbY). Phylogenetic analyses have supported the HGT inferences as the plant lineages are all clustering closely with distant microbial lineages. Deep-learning-based protein structure prediction and analyses, in combination with expression assessment based on codon adaption index (CAI) further corroborated the functionality and expressivity of the HGT genes in plant genomes. A case-study applying fold comparison and molecular dynamics (MD) of the HGT-driven CytP450 gave a more detailed illustration on the resemblance and evolutionary linkage between the plant recipient and microbial donor sequences. Together, the microbe-originated HGT genes identified in plant genomes and their participation in abiotic pressures resistance indicate a more profound impact of HGT on the adaptive evolution of plants.

Genome-enabled discovery of candidate virulence loci in Striga hermonthica, a devastating parasite of African cereal crops

Parasites have evolved proteins, virulence factors (VFs), that facilitate plant colonisation, however VFs mediating parasitic plant-host interactions are poorly understood. Striga hermonthica is an obligate, root-parasitic plant of cereal hosts in sub-Saharan Africa, causing devastating yield losses. Understanding the molecular nature and allelic variation of VFs in S. hermonthica is essential for breeding resistance and delaying the evolution of parasite virulence. We assembled the S. hermonthica genome and identified secreted proteins using in silico prediction. Pooled sequencing of parasites growing on a susceptible and a strongly resistant rice host allowed us to scan for loci where selection imposed by the resistant host had elevated the frequency of alleles contributing to successful colonisation. Thirty-eight putatively secreted VFs had very different allele frequencies with functions including host cell wall modification, protease or protease inhibitor and kinase activities. These candidate loci had significantly higher Tajima's D than the genomic background, consistent with balancing selection. Our results reveal diverse strategies used by S. hermonthica to overcome different layers of host resistance. Understanding the maintenance of variation at virulence loci by balancing selection will be critical to managing the evolution of virulence as part of a sustainable control strategy.

Parasitic plant small RNA analyses unveil parasite-specific signatures of microRNA retention, loss, and gain

Parasitism is a successful life strategy that has evolved independently in several families of vascular plants. The genera Cuscuta and Orobanche represent examples of the two profoundly different groups of parasites: one parasitizing host shoots and the other infecting host roots. In this study, we sequenced and described the overall repertoire of small RNAs from Cuscuta campestris and Orobanche aegyptiaca. We showed that C. campestris contains a number of novel microRNAs (miRNAs) in addition to a conspicuous retention of miRNAs that are typically lacking in other Solanales, while several typically conserved miRNAs seem to have become obsolete in the parasite. One new miRNA appears to be derived from a horizontal gene transfer event. The exploratory analysis of the miRNA population (exploratory due to the absence of a full genomic sequence for reference) from the root parasitic O. aegyptiaca also revealed a loss of a number of miRNAs compared to photosynthetic species from the same order. In summary, our study shows partly similar evolutionary signatures in the RNA silencing machinery in both parasites. Our data bear proof for the dynamism of this regulatory mechanism in parasitic plants.

Biological Features and In Planta Transcriptomic Analyses of a Microviridae Phage (CLasMV1) in "Candidatus Liberibacter asiaticus"

“Candidatus Liberibacter asiaticus” (CLas) is the causal agent of citrus Huanglongbing (HLB, also called citrus greening disease), a highly destructive disease threatening citrus production worldwide. A novel Microviridae phage (named CLasMV1) has been found to infect CLas, providing a potential therapeutic strategy for CLas/HLB control. However, little is known about the CLasMV1 biology. In this study, we analyzed the population dynamics of CLasMV1 between the insect vector of CLas, the Asian citrus psyllid (ACP, Diaphorina citri Kuwayama) and the holoparasitic dodder plant (Cuscuta campestris Yunck.); both acquired CLasMV1-infected CLas from an HLB citrus. All CLas-positive dodder samples were CLasMV1-positive, whereas only 32% of CLas-positive ACP samples were identified as CLasMV1-positive. Quantitative analyses showed a similar distribution pattern of CLasMV1 phage and CLas among eight citrus cultivars by presenting at highest abundance in the fruit pith and/or the center axis of the fruit. Transcriptome analyses revealed the possible lytic activity of CLasMV1 on CLas in fruit pith as evidenced by high-level expressions of CLasMV1 genes, and CLas genes related to cell wall biogenesis and remodeling to maintain the CLas cell envelope integrity. The up-regulation of CLas genes were involved in restriction–modification system that could involve possible phage resistance for CLas during CLasMV1 infection. In addition, the regulation of CLas genes involved in cell surface components and Sec pathway by CLasMV1 phage could be beneficial for phage infection. This study expanded our knowledge of CLasMV1 phage that will benefit further CLas phage research and HLB control.

A chromosome-level reference genome of a Convolvulaceae species Ipomoea cairica

Ipomoea cairica is a perennial creeper that has been widely introduced as a garden ornamental across tropical, subtropical, and temperate regions. Because it grows extremely fast and spreads easily, it has been listed as an invasive species in many countries. Here, we constructed the chromosome-level reference genome of Ipomoea cairica by Pacific Biosciences HiFi and Hi-C sequencing, with the assembly size of 733.0 Mb, the contig N50 of 43.8 Mb, the scaffold N50 of 45.7 Mb, and the Benchmarking Universal Single-Copy Orthologs complete rate of 98.0%. Hi-C scaffolding assigned 97.9% of the contigs to 15 pseudo-chromosomes. Telomeric repeat analysis reveals that 7 of the 15 pseudo-chromosomes are gapless and telomere to telomere. The transposable element content of Ipomoea cairica is 73.4%, obviously higher than that of other Ipomoea species. A total of 38,115 protein-coding genes were predicted, with the Benchmarking Universal Single-Copy Orthologs complete rate of 98.5%, comparable to that of the genome assembly, and 92.6% of genes were functional annotated. In addition, we identified 3,039 tRNA genes and 2,403 rRNA genes in the assembled genome. Phylogenetic analysis showed that Ipomoea cairica formed a clade with Ipomoea aquatica, and they diverged from each other 8.1 million years ago. Through comparative genome analysis, we reconfirmed that a whole genome triplication event occurred specific to Convolvulaceae family and in the ancestor of the genus Ipomoea and Cuscuta. This high-quality reference genome of Ipomoea cairica will greatly facilitate the studies on the molecular mechanisms of its rapid growth and invasiveness.

Systematics and Plastome Evolution in Schizaeaceae

While the family Schizaeaceae (Schizaeales) represents only about 0.4% of the extant fern species diversity, it differs from other ferns greatly in gross morphologies, niche preferences, and life histories. One of the most notable features in this family is its mycoheterotrophic life style in the gametophytic stage, which appears to be associated with extensive losses of plastid genes. However, the limited number of sequenced plastomes, and the lack of a well-resolved phylogenetic framework of Schizaeaceae, makes it difficult to gain any further insight. Here, with a comprehensive sampling of ~77% of the species diversity of this family, we first inferred a plastid phylogeny of Schizaeaceae using three DNA regions. To resolve the deep relationships within this family, we then reconstructed a plastome-based phylogeny focusing on a selection of representatives that covered all the major clades. From this phylogenomic backbone, we traced the evolutionary histories of plastid genes and examined whether gene losses were associated with the evolution of gametophytic mycoheterotrophy. Our results reveal that extant Schizaeaceae is comprised of four major clades-Microschizaea, Actinostachys, Schizaea, and Schizaea pusilla. The loss of all plastid NADH-like dehydrogenase (ndh) genes was confirmed to have occurred in the ancestor of extant Schizaeaceae, which coincides with the evolution of mycoheterotrophy in this family. For chlorophyll biosynthesis genes (chl), the losses were interpreted as convergent in Schizaeaceae, and found not only in Actinostachys, a clade producing achlorophyllous gametophytes, but also in S. pusilla with chlorophyllous gametophytes. In addition, we discovered a previously undescribed but phylogenetically distinct species hidden in the Schizaea dichotoma complex and provided a taxonomic treatment and morphological diagnostics for this new species-Schizaea medusa. Finally, our phylogenetic results suggest that the current PPG I circumscription of Schizaea is non-monophyletic, and we therefore proposed a three-genus classification moving a subset of Schizaea species sensu PPG I to a third genus-Microschizaea.

Mitochondrial phylogenomics of Cuscuta (Convolvulaceae) reveals a potentially functional horizontal gene transfer from the host

Horizontal gene transfers (HGTs) from host or other organisms have been reported in mitochondrial genomes of parasitic plants. Genes transferred in this fashion have usually been found non-functional. Several examples of HGT from the mitochondrial genome of parasitic Cuscuta (Convolvulaceae) to its hosts have been reported, but not vice versa. Here we used 31 protein-coding mitochondrial genes to infer the phylogeny of Cuscuta, and compared it with previous nuclear and plastid phylogenetic estimates. We also investigated the presence of HGTs within these lineages. Unlike in plastid genomes, we did not find extensive gene loss in their mitochondrial counterparts. Our results reveal the first example of organellar HGT from host to Cuscuta. Mitochondrial atp1 genes of South African subgenus Pachystigma were inferred to be transferred from Lamiales, with high support. Moreover, the horizontally transferred atp1 gene has functionally replaced the native, vertically transmitted copy, has an intact open reading frame, and is under strong purifying selection, all of which suggests that this xenolog remains functional. The mitochondrial phylogeny of Cuscuta is generally consistent with previous plastid and nuclear phylogenies, except for the misplacement of Pachystigma when atp1 is included. This incongruence may be caused by the HGT mentioned above. No example of HGT was found within mitochondrial genes of three other, independently evolved parasitic lineages we sampled: Cassytha/Laurales, Krameria/Zygophyllales, and Lennooideae/Boraginales.

Transcriptional Basis for Haustorium Formation and Host Establishment in Hemiparasitic Psittacanthus schiedeanus Mistletoes

The mistletoe Psittacanthus schiedeanus, a keystone species in interaction networks between plants, pollinators, and seed dispersers, infects a wide range of native and non-native tree species of commercial interest. Here, using RNA-seq methodology we assembled the whole circularized quadripartite structure of P. schiedeanus chloroplast genome and described changes in the gene expression of the nuclear genomes across time of experimentally inoculated seeds. Of the 140,467 assembled and annotated uniGenes, 2,000 were identified as differentially expressed (DEGs) and were classified in six distinct clusters according to their expression profiles. DEGs were also classified in enriched functional categories related to synthesis, signaling, homoeostasis, and response to auxin and jasmonic acid. Since many orthologs are involved in lateral or adventitious root formation in other plant species, we propose that in P. schiedeanus (and perhaps in other rootless mistletoe species), these genes participate in haustorium formation by complex regulatory networks here described. Lastly, and according to the structural similarities of P. schiedeanus enzymes with those that are involved in host cell wall degradation in fungi, we suggest that a similar enzymatic arsenal is secreted extracellularly and used by mistletoes species to easily parasitize and break through tissues of the host.

An artificial host system enables the obligate parasite Cuscuta campestris to grow and reproduce in vitro

Cuscuta campestris is an obligate parasitic plant that requires a host to complete its life cycle. Parasite-host connections occur via a haustorium, a unique organ that acts as a bridge for the uptake of water, nutrients, and macromolecules. Research on Cuscuta is often complicated by host influences, but comparable systems for growing the parasite in the absence of a host do not exist. We developed an axenic method to grow C. campestris on an artificial host system (AHS). We evaluated the effects of nutrients and phytohormones on parasite haustoria development and growth. Haustorium morphology and gene expression were also characterized. The AHS consists of an inert, fibrous stick that mimics a host stem, wicking water and nutrients to the parasite. It enables C. campestris to exhibit a parasitic habit and develop through all stages of its life cycle, including production of new shoots and viable seeds. The phytohormones 1-naphthaleneacetic acid and 6-benzylaminopurine affect haustoria morphology and increase parasite fresh weight and biomass. Unigene expression in AHS haustoria reflects processes similar to those in haustoria on living host plants. The AHS is a methodological improvement for studying Cuscuta biology by avoiding specific host effects on the parasite and giving researchers full control of the parasite environment.

An artificial host system enables the obligate parasite Cuscuta campestris to grow and reproduce in vitro

Cuscuta campestris is an obligate parasitic plant that requires a host to complete its life cycle. Parasite-host connections occur via a haustorium, a unique organ that acts as a bridge for the uptake of water, nutrients, and macromolecules. Research on Cuscuta is often complicated by host influences, but comparable systems for growing the parasite in the absence of a host do not exist. We developed an axenic method to grow C. campestris on an artificial host system (AHS). We evaluated the effects of nutrients and phytohormones on parasite haustoria development and growth. Haustorium morphology and gene expression were also characterized. The AHS consists of an inert, fibrous stick that mimics a host stem, wicking water and nutrients to the parasite. It enables C. campestris to exhibit a parasitic habit and develop through all stages of its life cycle, including production of new shoots and viable seeds. The phytohormones 1-naphthaleneacetic acid and 6-benzylaminopurine affect haustoria morphology and increase parasite fresh weight and biomass. Unigene expression in AHS haustoria reflects processes similar to those in haustoria on living host plants. The AHS is a methodological improvement for studying Cuscuta biology by avoiding specific host effects on the parasite and giving researchers full control of the parasite environment.

A host-free transcriptome for haustoriogenesis in Cuscuta campestris: Signature gene expression identifies markers of successive development stages

The development of the infection organ of the parasitic angiosperm genus Cuscuta is a dynamic process that is normally obscured from view as it happens endophytically in its host. We artificially induced haustoriogenesis in Cuscuta campestris by far-red light to define specific morphologically different stages and analyze their transcriptional patterns. This information enabled us to extract sets of high-confidence housekeeping and marker genes for the different stages, validated in a natural infection setting on a compatible host. This study provides a framework for more reproducible investigations of haustoriogenesis and the processes governing host-parasite interactions in shoot parasites, with C. campestris as a model species.

Unprecedented organelle genomic variations in morning glories reveal independent evolutionary scenarios of parasitic plants and the diversification of plant mitochondrial complexes

Background

The morning glories (Convolvulaceae) are distributed worldwide and produce economically important crops, medicinal herbs, and ornamentals. Members of this family are diverse in morphological characteristics and trophic modes, including the leafless parasitic Cuscuta (dodders). Organelle genomes were generally used for studying plant phylogeny and genomic variations. Notably, plastomes in parasitic plants always show non-canonical features, such as reduced size and accelerated rates. However, few organelle genomes of this group have been sequenced, hindering our understanding of their evolution, and dodder mitogenome in particular.

Results

We assembled 22 new mitogenomes and 12 new plastomes in Convolvulaceae. Alongside previously known ones, we totally analyzed organelle genomes of 23 species in the family. Our sampling includes 16 leafy autotrophic species and 7 leafless parasitic dodders, covering 8 of the 12 tribes. Both the plastid and mitochondrial genomes of these plants have encountered variations that were rarely observed in other angiosperms. All of the plastomes possessed atypical IR boundaries. Besides the gene and IR losses in dodders, some leafy species also showed gene and intron losses, duplications, structural variations, and insertions of foreign DNAs. The phylogeny reconstructed by plastid protein coding sequences confirmed the previous relationship of the tribes. However, the monophyly of ‘Merremieae’ and the sister group of Cuscuta remained uncertain. The mitogenome was significantly inflated in Cuscuta japonica, which has exceeded over 800 kb and integrated massive DNAs from other species. In other dodders, mitogenomes were maintained in small size, revealing divergent evolutionary strategies. Mutations unique to plants were detected in the mitochondrial gene ccmFc, which has broken into three fragments through gene fission and splicing shift. The unusual changes likely initially happened to the common ancestor of the family and were caused by a foreign insertion from rosids followed by double-strand breaks and imprecise DNA repairs. The coding regions of ccmFc expanded at both sides after the fission, which may have altered the protein structure.

Conclusions

Our family-scale analyses uncovered unusual scenarios for both organelle genomes in Convolvulaceae, especially in parasitic plants. The data provided valuable genetic resources for studying the evolution of Convolvulaceae and plant parasitism.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12915-022-01250-1.

Investigating Host and Parasitic Plant Interaction by Tissue-Specific Gene Analyses on Tomato and Cuscuta campestris Interface at Three Haustorial Developmental Stages

Parasitic weeds cause billions of dollars in agricultural losses each year worldwide. Cuscuta campestris (C. campestris), one of the most widespread and destructive parasitic plants in the United States, severely reduces yield in tomato plants. Reducing the spread of parasitic weeds requires understanding the interaction between parasites and hosts. Several studies have identified factors needed for parasitic plant germination and haustorium induction, and genes involved in host defense responses. However, knowledge of the mechanisms underlying the interactions between host and parasitic plants, specifically at the interface between the two organisms, is relatively limited. A detailed investigation of the crosstalk between the host and parasite at the tissue-specific level would enable development of effective parasite control strategies. To focus on the haustorial interface, we used laser-capture microdissection (LCM) with RNA-seq on early, intermediate and mature haustorial stages. In addition, the tomato host tissue that immediately surround the haustoria was collected to obtain tissue- resolution RNA-Seq profiles for C. campestris and tomato at the parasitism interface. After conducting RNA-Seq analysis and constructing gene coexpression networks (GCNs), we identified CcHB7, CcPMEI, and CcERF1 as putative key regulators involved in C. campestris haustorium organogenesis, and three potential regulators, SlPR1, SlCuRe1-like, and SlNLR, in tomatoes that are involved in perceiving signals from the parasite. We used host-induced gene silencing (HIGS) transgenic tomatoes to knock-down the candidate genes in C. campestris and produced CRISPR transgenic tomatoes to knock out candidate genes in tomatoes. The interactions of C. campestris with these transgenic lines were tested and compared with that in wild-type tomatoes. The results of this study reveal the tissue-resolution gene regulatory mechanisms at the parasitic plant-host interface and provide the potential of developing a parasite-resistant system in tomatoes.

Karyotype asymmetry in Cuscuta L. subgenus Pachystigma reflects its repeat DNA composition

Cuscuta is a cytogenetically diverse genus, with karyotypes varying 18-fold in chromosome number and 127-fold in genome size. Each of its four subgenera also presents particular chromosomal features, such as bimodal karyotypes in Pachystigma. We used low coverage sequencing of the Cuscuta nitida genome (subgenus Pachystigma), as well as chromosome banding and molecular cytogenetics of three subgenus representatives, to understand the origin of bimodal karyotypes. All three species, C. nitida, C. africana (2n = 28) and C. angulata (2n = 30), showed heterochromatic bands mainly in the largest chromosome pairs. Eighteen satellite DNAs were identified in C. nitida genome, two showing similarity to mobile elements. The most abundant were present at the largest pairs, as well as the highly abundant ribosomal DNAs. The most abundant Ty1/Copia and Ty3/Gypsy elements were also highly enriched in the largest pairs, except for the Ty3/Gypsy CRM, which also labelled the pericentromeric regions of the smallest chromosomes. This accumulation of repetitive DNA in the larger pairs indicates that these sequences are largely responsible for the formation of bimodal karyotypes in the subgenus Pachystigma. The repetitive DNA fraction is directly linked to karyotype evolution in Cuscuta.

Bidirectional lncRNA Transfer between Cuscuta Parasites and Their Host Plant

Dodder species (Cuscuta spp.) are holoparasites that have extensive material exchange with their host plants through vascular connections. Recent studies on cross-species transfer have provided breakthrough insights, but little is known about the interaction mechanisms of the inter-plant mobile substances in parasitic systems. We sequenced the transcriptomes of dodder growing on soybean hosts to characterize the long non-coding RNA (lncRNA) transfer between the two species, and found that lncRNAs can move in high numbers (365 dodder lncRNAs and 14 soybean lncRNAs) in a bidirectional manner. Reverse transcription-polymerase chain reaction further confirmed that individual lncRNAs were trafficked in the dodder-soybean parasitic system. To reveal the potential functions of mobile transcripts, the Gene Ontology terms of mobile lncRNA target genes were predicted, and mobile dodder target genes were found to be mainly enriched in "metabolic process", "catalytic activity", "signaling", and "response to stimulus" categories, whereas mobile soybean target genes were enriched in organelle-related categories, indicating that specific mobile lncRNAs may be important in regulating dodder parasitism. Our findings reveal that lncRNAs are transferred between dodder and its host soybean plants, which may act as critical regulators to coordinate the host-dodder interaction at the whole parasitic level.

The Evolution of Cytogenetic Traits in Cuscuta (Convolvulaceae), the Genus With the Most Diverse Chromosomes in Angiosperms

Karyotypes are characterized by traits such as chromosome number, which can change through whole-genome duplication and dysploidy. In the parasitic plant genus Cuscuta (Convolvulaceae), chromosome numbers vary more than 18-fold. In addition, species of this group show the highest diversity in terms of genome size among angiosperms, as well as a wide variation in the number and distribution of 5S and 35S ribosomal DNA (rDNA) sites. To understand its karyotypic evolution, ancestral character state reconstructions were performed for chromosome number, genome size, and position of 5S and 35S rDNA sites. Previous cytogenetic data were reviewed and complemented with original chromosome counts, genome size estimates, and rDNA distribution assessed via fluorescence in situ hybridization (FISH), for two, seven, and 10 species, respectively. Starting from an ancestral chromosome number of x = 15, duplications were inferred as the prevalent evolutionary process. However, in holocentric clade (subgenus Cuscuta), dysploidy was identified as the main evolutionary mechanism, typical of holocentric karyotypes. The ancestral genome size of Cuscuta was inferred as approximately 1C = 12 Gbp, with an average genome size of 1C = 2.8 Gbp. This indicates an expansion of the genome size relative to other Convolvulaceae, which may be linked to the parasitic lifestyle of Cuscuta. Finally, the position of rDNA sites varied mostly in species with multiple sites in the same karyotype. This feature may be related to the amplification of rDNA sites in association to other repeats present in the heterochromatin. The data suggest that different mechanisms acted in different subgenera, generating the exceptional diversity of karyotypes in Cuscuta.

Regulatory Modules Involved in the Degradation and Modification of Host Cell Walls During Cuscuta campestris Invasion

Haustoria of parasitic plants have evolved sophisticated traits to successfully infect host plants. The degradation and modification of host cell walls enable the haustorium to effectively invade host tissues. This study focused on two APETALA2/ETHYLENE RESPONSE FACTOR (ERF) genes and a set of the cell wall enzyme genes principally expressed during the haustorial invasion of Cuscuta campestris Yuncker. The orthogroups of the TF and cell wall enzyme genes have been implicated in the cell wall degradation and modification activities in the abscission of tomatoes, which are currently the phylogenetically closest non-parasitic model species of Cuscuta species. Although haustoria are generally thought to originate from root tissues, our results suggest that haustoria have further optimized invasion potential by recruiting regulatory modules from other biological processes.

The minicircular and extremely heteroplasmic mitogenome of the holoparasitic plant Rhopalocnemis phalloides

The plastid and nuclear genomes of parasitic plants exhibit deeply altered architectures,^1-13 whereas the few examined mitogenomes range from deeply altered to conventional.^14-20 To provide further insight on mitogenome evolution in parasitic plants, we report the highly modified mitogenome of Rhopalocnemis phalloides, a holoparasite in Balanophoraceae. Its mitogenome is uniquely arranged in 21 minicircular chromosomes that vary in size from 4,949 to 7,861 bp, with a total length of only 130,713 bp. All chromosomes share an identical 896 bp conserved region, with a large stem-loop that acts as the origin of replication, flanked on each side by hypervariable and semi-conserved regions. Similar minicircular structures with shared and unique regions have been observed in parasitic animals and free-living protists,^21-24 suggesting convergent structural evolution. Southern blots confirm both the minicircular structure and the replication origin of the mitochondrial chromosomes. PacBio reads provide evidence for chromosome recombination and rolling-circle replication for the R. phalloides mitogenome. Despite its small size, the mitogenome harbors a typical set of genes and introns within the unique regions of each chromosome, yet introns are the smallest among seed plants and ferns. The mitogenome also exhibits extreme heteroplasmy, predominantly involving short indels and more complex variants, many of which cause potential loss-of-function mutations for some gene copies. All heteroplasmic variants are transcribed, and functional and nonfunctional protein-coding variants are spliced and RNA edited. Our findings offer a unique perspective into how mitogenomes of parasitic plants can be deeply altered and shed light on plant mitogenome replication.

An angiosperm NLR Atlas reveals that NLR gene reduction is associated with ecological specialization and signal transduction component deletion

Nucleotide-binding leucine-rich-repeat (NLR) genes comprise the largest family of plant disease-resistance genes. Angiosperm NLR genes are phylogenetically divided into the TNL, CNL, and RNL subclasses. NLR copy numbers and subclass composition vary tremendously across angiosperm genomes. However, the evolutionary associations between genomic NLR content and ecological adaptation, or between NLR content and signal transduction components, are poorly characterized because of limited genome availability. In this study, we established an angiosperm NLR atlas (ANNA, https://biobigdata.nju.edu.cn/ANNA/) that includes NLR genes from over 300 angiosperm genomes. Using ANNA, we revealed that NLR copy numbers differ up to 66-fold among closely related species owing to rapid gene loss and gain. Interestingly, NLR contraction was associated with adaptations to aquatic, parasitic, and carnivorous lifestyles. The convergent NLR reduction in aquatic plants resembles the lack of NLR expansion during the long-term evolution of green algae before the colonization of land. A co-evolutionary pattern between NLR subclasses and plant immune pathway components was also identified, suggesting that immune pathway deficiencies may drive TNL loss. Finally, we identified a conserved TNL lineage that may function independently of the EDS1-SAG101-NRG1 module. Collectively, these findings provide new insights into the evolution of NLR genes in the context of ecological adaptation and genome content variation.

A chromosome-scale Gastrodia elata genome and large-scale comparative genomic analysis indicate convergent evolution by gene loss in mycoheterotrophic and parasitic plants

Mycoheterotrophic and parasitic plants are heterotrophic and parasitize on fungi and plants, respectively, to obtain nutrients. Large-scale comparative genomics analysis has not been conducted in mycoheterotrophic or parasitic plants or between these two groups of parasites. We assembled a chromosome-level genome of the fully mycoheterotrophic plant Gastrodia elata (Orchidaceae) and performed comparative genomic analyses on the genomes of G. elata and four orchids (initial mycoheterotrophs), three parasitic plants (Cuscuta australis, Striga asiatica, and Sapria himalayana), and 36 autotrophs from various angiosperm lineages. It was found that while in the hemiparasite S. asiatica and initial mycoheterotrophic orchids, approximately 4-5% of the conserved orthogroups were lost, the fully heterotrophic G. elata and C. australis both lost approximately 10% of the conserved orthogroups, indicating that increased heterotrophy is positively associated with gene loss. Importantly, many genes that are essential for autotrophs, including those involved in photosynthesis, the circadian clock, flowering time regulation, immunity, nutrient uptake, and root and leaf development, were convergently lost in both G. elata and C. australis. The high-quality genome of G. elata will facilitate future studies on the physiology, ecology, and evolution of mycoheterotrophic plants, and our findings highlight the critical role of gene loss in the evolution of plants with heterotrophic lifestyles.

Mitochondrial genomes of two parasitic Cuscuta species lack clear evidence of horizontal gene transfer and retain unusually fragmented ccmFC genes

Background

The intimate association between parasitic plants and their hosts favours the exchange of genetic material, potentially leading to horizontal gene transfer (HGT) between plants. With the recent publication of several parasitic plant nuclear genomes, there has been considerable focus on such non-sexual exchange of genes. To enhance the picture on HGT events in a widely distributed parasitic genus, Cuscuta (dodders), we assembled and analyzed the organellar genomes of two recently sequenced species, C. australis and C. campestris, making this the first account of complete mitochondrial genomes (mitogenomes) for this genus.

Results

The mitogenomes are 265,696 and 275,898 bp in length and contain a typical set of mitochondrial genes, with 10 missing or pseudogenized genes often lost from angiosperm mitogenomes. Each mitogenome also possesses a structurally unusual ccmF_C gene, which exhibits splitting of one exon and a shift to trans-splicing of its intron. Based on phylogenetic analysis of mitochondrial genes from across angiosperms and similarity-based searches, there is little to no indication of HGT into the Cuscuta mitogenomes. A few candidate regions for plastome-to-mitogenome transfer were identified, with one suggestive of possible HGT.

Conclusions

The lack of HGT is surprising given examples from the nuclear genomes, and may be due in part to the relatively small size of the Cuscuta mitogenomes, limiting the capacity to integrate foreign sequences.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-08105-z.