The origin of novel traits is recognized as an important process underlying many major evolutionary radiations. We studied the genetic basis for the evolution of haustoria, the novel feeding organs of parasitic flowering plants, using comparative transcriptome sequencing in three species of Orobanchaceae. Around 180 genes are upregulated during haustorial development following host attachment in at least two species, and these are enriched in proteases, cell wall modifying enzymes, and extracellular secretion proteins. Additionally, about 100 shared genes are upregulated in response to haustorium inducing factors prior to host attachment. Collectively, we refer to these newly identified genes as putative “parasitism genes.” Most of these parasitism genes are derived from gene duplications in a common ancestor of Orobanchaceae and Mimulus guttatus, a related nonparasitic plant. Additionally, the signature of relaxed purifying selection and/or adaptive evolution at specific sites was detected in many haustorial genes, and may play an important role in parasite evolution. Comparative analysis of gene expression patterns in parasitic and nonparasitic angiosperms suggests that parasitism genes are derived primarily from root and floral tissues, but with some genes co-opted from other tissues. Gene duplication, often taking place in a nonparasitic ancestor of Orobanchaceae, followed by regulatory neofunctionalization, was an important process in the origin of parasitic haustoria.

Highlights d The Striga genome reflects a three-phase model of parasitic plant genome evolution d A family of strigolactone receptors has undergone a striking expansion in Striga d Genes in lateral root development are coordinately induced in a parasitic organ d Host genes and retrotransposons are horizontally transferred into Striga

Plastid genomes of photosynthetic flowering plants are usually highly conserved in both structure and gene content. However, the plastomes of parasitic and mycoheterotrophic plants may be released from selective constraint due to the reduction or loss of photosynthetic ability. Here we present the greatly reduced and highly divergent, yet functional, plastome of the nonphotosynthetic holoparasite Hydnora visseri (Hydnoraceae, Piperales). The plastome is 27 kb in length, with 24 genes encoding ribosomal proteins, ribosomal RNAs, tRNAs, and a few nonbioenergetic genes, but no genes related to photosynthesis. The inverted repeat and the small single copy region are only approximately 1.5 kb, and intergenic regions have been drastically reduced. Despite extreme reduction, gene order and orientation are highly similar to the plastome of Piper cenocladum, a related photosynthetic plant in Piperales. Gene sequences in Hydnora are highly divergent and several complementary approaches using the highest possible sensitivity were required for identification and annotation of this plastome. Active transcription is detected for all of the protein-coding genes in the plastid genome, and one of two introns is appropriately spliced out of rps12 transcripts. The whole-genome shotgun read depth is 1,400× coverage for the plastome, whereas the mitochondrial genome is covered at 40× and the nuclear genome at 2×. Despite the extreme reduction of the genome and high sequence divergence, the presence of syntenic, long transcriptionally active open-reading frames with distant similarity to other plastid genomes and a high plastome stoichiometry relative to the mitochondrial and nuclear genomes suggests that the plastome remains functional in H. visseri. A four-stage model of gene reduction, including the potential for complete plastome loss, is proposed to account for the range of plastid genomes in nonphotosynthetic plants.

The Parasitic Plant Genome Project has sequenced transcripts from three parasitic species and a nonparasitic relative in the Orobanchaceae with the goal of understanding genetic changes associated with parasitism. The species studied span the trophic spectrum from free-living nonparasite to obligate holoparasite. Parasitic species used wereTriphysaria versicolor, a photosynthetically competent species that opportunistically parasitizes roots of neighboring plants;Striga hermonthica, a hemiparasite that has an obligate need for a host; andOrobanche aegyptiaca, a holoparasite with absolute nutritional dependence on a host.Lindenbergia philippensisrepresents the closest nonparasite sister group to the parasitic Orobanchaceae and was included for comparative purposes. Tissues for transcriptome sequencing from each plant were gathered to identify expressed genes for key life stages from seed conditioning through anthesis. Two of the species studied,S. hermonthicaandO. aegyptiaca, are economically important weeds and the data generated by this project are expected to aid in research and control of these species and their relatives. The sequences generated through this project will provide an abundant resource of molecular markers for understanding population dynamics, as well as provide insight into the biology of parasitism and advance progress toward understanding parasite virulence and host resistance mechanisms. In addition, the sequences provide important information on target sites for herbicide action or other novel control strategies such as trans-specific gene silencing.

SummaryAmino acids represent the major form of reduced nitrogen that is transported in plants. Amino acid transporters in plants often show tissue-specific expression patterns and are used by plants to transport these metabolites from source to sink during development and under changing environmental conditions. We identified one amino acid transporter, AtCAT6, which is expressed in sink tissues such as lateral root primordia, flowers and seeds. Additionally AtCAT6 was induced during infestation of roots by the plantparasitic root-knot nematode, Meloidogyne incognita. Quantitative reverse-transcriptase PCR revealed nematode inducibility throughout the duration of nematode infestation and in nematode-induced feeding sites. Promoter analyses confirmed expression in endogenous sink tissues and nematode-induced feeding sites. In Xenopus oocytes, AtCAT6 mediated electrogenic transport of proteinogenic as well as nonproteinogenic amino acids with moderate affinity. AtCAT6 transported large, neutral and cationic amino acids in preference to other amino acids. Knockout mutants of this transporter failed to grow on medium containing L-glutamine as the sole nitrogen source. Our data suggest that AtCAT6 plays a role in supplying amino acids to sink tissues of plants and nematode-induced feeding structures.

Horizontal gene transfer (HGT) is the transfer of genetic material across species boundaries and has been a driving force in prokaryotic evolution. HGT involving eukaryotes appears to be much less frequent, and the functional implications of HGT in eukaryotes are poorly understood. We test the hypothesis that parasitic plants, because of their intimate feeding contacts with host plant tissues, are especially prone to horizontal gene acquisition. We sought evidence of HGTs in transcriptomes of three parasitic members of Orobanchaceae, a plant family containing species spanning the full spectrum of parasitic capabilities, plus the free-living Lindenbergia. Following initial phylogenetic detection and an extensive validation procedure, 52 highconfidence horizontal transfer events were detected, often from lineages of known host plants and with an increasing number of HGT events in species with the greatest parasitic dependence. Analyses of intron sequences in putative donor and recipient lineages provide evidence for integration of genomic fragments far more often than retro-processed RNA sequences. Purifying selection predominates in functionally transferred sequences, with a small fraction of adaptively evolving sites. HGT-acquired genes are preferentially expressed in the haustorium-the organ of parasitic plants-and are strongly biased in predicted gene functions, suggesting that expression products of horizontally acquired genes are contributing to the unique adaptive feeding structure of parasitic plants.

Parasitism in flowering plants has evolved at least 11 times [1]. Only one family, Orobanchaceae, comprises all major nutritional types of parasites: facultative, hemiparasitic (partially photosynthetic), and holoparasitic (nonphotosynthetic) [2]. Additionally, the family includes Lindenbergia, a nonparasitic genus sister to all parasitic Orobanchaceae [3-6]. Parasitic Orobanchaceae include species with severe economic impacts: Striga (witchweed), for example, affects over 50 million hectares of crops in sub-Saharan Africa, causing more than $3 billion in damage annually [7]. Although gene losses and increased substitution rates have been characterized for parasitic plant plastid genomes [5, 8-11], the nuclear genome and transcriptome remain largely unexplored. The Parasitic Plant Genome Project (PPGP; http://ppgp.huck.psu.edu/) [2] is leveraging the natural variation in Orobanchaceae to explore the evolution and genomic consequences of parasitism in plants through a massive transcriptome and gene discovery project involving Triphysaria versicolor (facultative hemiparasite), Striga hermonthica (obligate hemiparasite), and Phelipanche aegyptiaca (Orobanche [12]; holoparasite). Here we present the first set of large-scale genomic resources for parasitic plant comparative biology. Transcriptomes of above-ground tissues reveal that, in addition to the predictable loss of photosynthesis-related gene expression in P. aegyptiaca, the nonphotosynthetic parasite retains an intact, expressed, and selectively constrained chlorophyll synthesis pathway.

Whereas de novo assemblies of RNA-Seq data are being published for a growing number of species across the tree of life, there are currently no broadly accepted methods for evaluating such assemblies. Here we present a detailed comparison of 99 transcriptome assemblies, generated with 6 de novo assemblers including CLC, Trinity, SOAP, Oases, ABySS and NextGENe. Controlled analyses of de novo assemblies for Arabidopsis thaliana and Oryza sativa transcriptomes provide new insights into the strengths and limitations of transcriptome assembly strategies. We find that the leading assemblers generate reassuringly accurate assemblies for the majority of transcripts. At the same time, we find a propensity for assemblers to fail to fully assemble highly expressed genes. Surprisingly, the instance of true chimeric assemblies is very low for all assemblers. Normalized libraries are reduced in highly abundant transcripts, but they also lack 1000s of low abundance transcripts. We conclude that the quality of de novo transcriptome assemblies is best assessed through consideration of a combination of metrics: 1) proportion of reads mapping to an assembly 2) recovery of conserved, widely expressed genes, 3) N50 length statistics, and 4) the total number of unigenes. We provide benchmark Illumina transcriptome data and introduce SCERNA, a broadly applicable modular protocol for de novo assembly improvement. Finally, our de novo assembly of the Arabidopsis leaf transcriptome revealed ~20 putative Arabidopsis genes lacking in the current annotation.