BACKGROUND
horizontal gene transfer  is the nonsexual transmission of genetic material across species boundaries  <cit> .
hgt is well known in bacteria, where hgt often results in adaptive gains of novel genes and traits  <cit> .
there are fewer well-documented cases of hgt among eukaryotes  <cit>  and the large majority of these cases appear to result in short-lived, nonfunctional sequences  <cit> .
consequently, the long-term evolutionary impact of hgt in multicellular eukaryotes remains largely unknown.
several cases of hgt are known or suspected in plants  <cit> , most involving mitochondrial sequences, and/or parasitic plants  <cit> .
parasitic plants form direct haustorial connections with their host plants and are capable of obtaining a wide range of macromolecules from their hosts, including viruses  <cit> , gene silencing signals  <cit> , and messenger rnas  <cit> .
consequently, parasites may have many opportunities for hgt events and an increased likelihood that some of these result in functional, and potentially adaptive, gene transfers.
two recent reports by yoshida et al  <cit>  and xi et al  <cit>  were the first indications that nuclear protein coding sequences, likely obtained from their respective host species, could be integrated into the genomes of parasitic plants by hgt.
these were important advances, but they provided few clues as to the long term impact of hgt, how the transgenes evolve, and how they may function.
we hypothesized that systematic analysis of genome-scale datasets from parasitic plants could lead to evidence for acquisition and long-term maintenance of functional gene sequences in plants that had been acquired via hgt.
albumin  <dig> genes are known only from a subset of species in the legume family  of angiosperms where they encode seed storage proteins and insect toxins  <cit> .
the albumin  <dig> proteins in legumes are  <dig> to  <dig> amino acids in length and rich in cysteine residues.
they form a unique protein structure known as a knottin, which has three disulfide bonds and is characterized by a “disulfide through disulfide knot”  <cit> .
the knottins are famous for their intruguing “disulfide through disulfide knot” structure and have been extensively studied in various fields, most of which are related with potentials in drug design  <cit> .
albumin  <dig> genes may have originated early in the diversification of papilionoid legumes  <cit> , but multiple homologous gene copies have been found only in species that are members of the more derived “millettioid s.l.” and “hologalegina” clades  <cit> .
orobanche s. l., often known by the common name “broomrape,” includes 150- <dig> obligate parasitic plant species in the family orobanchaceae.
growing evidence supports the segregation of broomrapes into four genera  <cit> : aphyllon , myzorrhiza , phelipanche , orobanche s. str.
.
most broomrape species have a narrow host spectrum and grow exclusively on perennial eudicot host plants  <cit> , with members of the leguminosae, solanaceae, and asteraceae among the more common hosts  <cit> .
as a member of order lamiales, orobanchaceae is phylogenetically well-separated from host members in these lineages, particularly legume hosts in the rosid order fabales .
a few broomrape species  have become devastating pests of important crop plants, affecting their growth and resource allocation and imparting significant losses in yield  <cit> .
p. aegyptiaca, the focal species in this study, has a broad host range that includes members of the eudicot families apiaceae, asteraceae, brassicaceae, cucurbitaceae, leguminosae, and solanaceae.
here we show that a gene encoding albumin  <dig> knottin-like protein, closely related to the albumin  <dig> genes, only known from papilionoid legumes, serving dual roles in food storage and as insect toxins, was found in phelipanche aegyptiaca and related parasitic species of family orobanchaceae, and was likely acquired by a phelipanche ancestor via hgt from a legume host based on phylogenetic analyses.
according to genomic sequences from nine related parasite species, 3d protein structure simulation tests, and evolutionary constraint analyses, the broomrape xenogene we identified here retains the intron structure, six highly conserved cysteine residues necessary to form a knottin protein, and displays levels of purifying selection like those seen in legumes.
the albumin  <dig> xenogene has evolved through > <dig> speciation events over ca.
<dig> million years, forming a small family of differentially expressed genes that may confer novel functions in the parasites.
RESULTS
the albumin  <dig> transcript was first identified as a hgt candidate in the transcriptome of p. aegyptiaca  using a blast-based  <cit>  bioinformatic screen .
albumin  <dig> transcripts were then searched further, using blastx, against the ncbi nr database and the plantgdb database  <cit> .
top hits were seen  to medicago truncatula albumin  <dig> sequences, with expected values of 5e- <dig> and 1e- <dig>  additional blast, including hidden markov model -based psi-blast searches with the sequence from p. aegyptiaca were performed to attempt to detect homologs in three other members of orobanchaceae with large transcriptome datasets  .
several large public databases, including phytozome  <cit> , plantgdb, and sol genomics network  <cit> , were also searched.
after searching  <dig> sequenced genomes and transcriptomes of  <dig> additional plant species, albumin  <dig> homologs were detected only in legumes and the transcriptome libraries of p. aegyptiaca.
having identified the albumin  <dig> sequence in the p. aegyptiaca transcriptome, genomic sequences encoding albumin  <dig> were then obtained from p. aegyptiaca and eight additional broomrape species, including p. schultzii, p. ramosa, p. mutelli, p. nana, and orobanche hederae, o. minor, o. cernua and o. ballotae.
the nucleotide sequence and inferred gene structures of the albumin  <dig> genes in broomrape species  are closely comparable, with inferred protein alignments  <dig> - <dig> % identical and  <dig> %- <dig> % similar  in ungapped regions between the legume and parasite proteins.
two albumin  <dig> genes were identified in phelipanche species, and are identified here as copy_ <dig> and copy_ <dig>  or albumin1- <dig> and albumin1- <dig>  respectively.
an intron disrupts the coding region at the same position in both genes and the intron sequences are similar but contain a number of insertion and deletion mutations.
only one albumin  <dig> gene was detected from orobanche species.
although the intron length in albumin  <dig> genes of phelipanche and legume species is not well conserved, several critical intron features are shared .
first, the starting position of the intron in both the p. aegyptiaca and m. truncatula sequences are the same, and the first nine base pairs are identical.
second, the introns have characteristic splicing sites at their 5’ and 3’ ends; 5’ ends often have gt/gu and 3’ ends often have ag, and these motifs are found in both m. truncatula and phelipanche albumin  <dig> introns .
albumin <dig> gene sequences from phelipanche were also searched with blastn against the ncbi nt database in order to search for high frequency repeats and mobile elements, but no such features were identified.
phylogenetic analysis  of all known plant albumin  <dig> sequences showed a strongly supported clade containing all of the albumin  <dig> sequences from broomrapes  boostrap  <dig>  bayesian inference  posterior probabilities   <dig> ) nested deeply within the irlc  of papilionoid legumes  <cit> .
among legumes, the next most closely related sequences  are from onobrychis argentea and onobrychis viciifolia.
because the node supporting the position of the broomrape clade  within the papilionoid legumes is relatively weakly supported, we also tested the hypothesis that the broomrape clade of albumin  <dig> sequences falls outside the larger clade of legumes represented in this analysis .
this hypothesis was rejected  relative to the maximum likelihood position as indicated in this tree.
two albumin  <dig> genes are resolved as sister clades in phelipanche species, which are in turn resolved as sister to the single gene obtained from orobanche species.
gene structures supported a similar conclusion .
the amino acid sequence alignments of albumin  <dig> from legumes to p. aegyptiaca show conservation of all cysteine residues essential for disulfide bond formation in albumin  <dig> proteins .
we investigated whether the predicted albumin  <dig> proteins from parasites maintain the characteristic knottin structures found in the legume albumin  <dig> proteins using knoter1d  <cit> .
simulated 3d structures show that the phelipanche albumin  <dig> proteins form a characteristic knottin structure with three-disulfide bonds and a “disulfide through disulfide knot”.
knottin protein structures are also predicted in all of the other full-length albumin  <dig> genes in phelipanche species.
knoter1d assigned scores greater than  <dig> to each phelipanche albumin  <dig> sequence; a score greater than  <dig> in this analysis passes the knoter1d criteria for identification as an albumin  <dig> structure.
the predicted 3d structures for p_aegyptiaca_albumin1- <dig>  and p_aegyptiaca_albumin1- <dig>  are very similar to the insect toxic albumin  <dig> protein from m. truncatula.
albumin  <dig>  a non-knottin legume protein, has no discernable homology with the albumin  <dig> protein in legumes .
having found that the horizontally acquired albumin <dig> genes were present in related species of broomrapes we then asked if the genes are evolving under purifying selection indicative of a functional protein coding sequence.
dn , ds  and dn/ds were calculated for all three lineages of the broomrape albumin  <dig> clade  and for the albumin  <dig> sequences from three closely related legume sequences; astragalus monspessulanus, onobrychis argentea and onobrychis viciifolia.
synonymous substitutions in the albumin  <dig> genes  outnumber non-synonymous substitutions  by at least 3: <dig> in most lineages , and dn/ds, reflecting the level of purifying selection, is similar in broomrapes to the value estimated for closely related albumin <dig> sequences from legumes.
all cysteine residues were also identified as evolving under purifying selection, suggesting that the horizontally acquired albumin  <dig> genes in broomrapes are functional
having observed evidence for selection for structural conservation, we investigated whether these genes exhibit transcription profiles that suggest a new or unique pattern of expression in parasites.
normalized expression levels of both albumin  <dig> genes in p. aegyptiaca were estimated as reads per kilobase per million reads  for eight libraries representing major stages of belowground and aboveground parasite development .
both genes displayed lowest expression levels at stage  <dig>  and highest at stage  <dig> .
transcripts were particularly abundant at stage  <dig>  , more than 1000x higher than the haustorial stage.
discussion
biogeographic overlap and common feeding interactions between diverse broomrapes and temperate papilionoid legumes increase the likelihood that the hgt event occurred in a common ancestor of the parasites that was in direct contact with legume host plants.
an alternative  explanation is that another organism or virus that co-occurred in the same habitats as the ancestral lineages served as a “stepping stone” for a two- or more step transfer.
however, this is not supported in our searches of the sequence databases.
based on fossil-calibrated age estimates of legume lineages  <cit> , we estimate that this horizontal acquisition occurred in an ancestral broomrape that lived in the miocene epoch, about  <dig> mya.
both parasite and their legume host groups have northern temperate distributions, and their lineages likely overlapped in the past as they do now, providing a minimal requirement for a horizontal gene transfer to occur.
another possibility, however unlikely, is that albumin <dig> was a more recent acquisition that underwent strong convergence at the protein level with this legume lineage.
however the branch lengths we observed in the phylogeny  were not unusually long in our dna-based phylogeny, and given the large collection of related sequences we obtained from other broomrape species, we have reduced any tendency the orobanche/phelipanche lineage may have had to connect by chance to a deep branch.
thus, the convergence hypothesis is not supported.
because the breadth of phelipanche and orobanche species we have sampled spans the deepest branches of broomrape diversity  <cit> , the albumin gene can be inferred to have survived through an extended evolutionary radiation of at least  <dig> species  <cit>  or more, if the number of now-extinct broomrape species could be estimated.
because the introns of phelipanche albumin  <dig> xenogenes maintain critical splicing sites and share the same starting positions and first nine base pairs with the known m. truncatula albumin  <dig> intron, it is likely that the hgt event in broomrapes involved transfer of a genomic sequence rather than a separate cdna.
following the transfer, albumin  <dig> genes in broomrape species have evolved under purifying selection consistent with what is observed in related legume albumin  <dig> genes.
this observation, as well as the stage-specific transcription patterns, conserved cysteine residues and predicted 3d knottin protein structures, strongly suggest that albumin  <dig> genes encode functional proteins in broomrape species, and could potentially serve a function similar to its role in legumes, providing a large pool of sulfur storage and exhibiting toxicity to insect herbivores in certain legumes  <cit> .
a recent report involves panicoid grass species with c <dig> or c <dig> photosynthetic pathways.
evidence was presented that nuclear genes were horizontally transferred between panicoid species and were subsequently adapted into the existing pathways with the effect of advancing the extent of c <dig> photosynthesis in some lineages  <cit> .
these results indicate that hgt may promote the sharing of adaptive traits among related species.
in comparison, the albumin example described here shows how a completely novel and highly specialized trait has been acquired at an ancestral stage from a distantly related donor species and maintained by the recipient lineage throughout an extended period of evolutionary history.
the albumin  <dig> genes in p. aegyptiaca are highly transcribed in most of the developmental stages we examined.
transcripts are more abundant in reproductive tissue, and lowest in the young haustorium , which represents the earliest point in our tissue sampling where the parasite is in direct contact with the host plant.
this suggests that the novel gene in p. aegyptiaca is probably not encoding a protein that is playing a direct role in the process of haustorial formation, and that albumin  <dig> expression is down-regulated as the parasite devotes energy to the essential process of establishing host vascular connections.
it is also possible that the low expression in the haustorial stage could help the parasite avoid detection or minimize a negative impact on the health of the host plant during early stages of parasite contact and feeding.
several other parasitic lineages, including members of cuscuta , cassytha , apodanthaceae, hydnoraceae, and the order santalales, regularly feed upon legumes  <cit>  and therefore might also have had opportunities to acquire albumin  <dig> sequences through hgt.
large transcriptome datasets are currently available for only two of these, the generalist parasite cuscuta pentagona  and for the legume specialist feeder pilostyles thurberi  <cit> .
both of these parasites, and other species in these genera, feed widely on legumes.
no homolog of albumin  <dig> was detected in blast searches of the pilostyles transcriptome in the 1kp dataset  <cit> .
however, albumin  <dig> sequences were detected in the same dataset and in two additional transcriptome libraries from cuscuta pentagona .
phylogenetic analysis nests the cuscuta sequences well within leguminosae, but on an independent branch from the broomrape sequences , suggesting that these transcripts in cuscuta represent a different hgt event into cuscuta from a lineage of papilionoid legumes that was different from the source of the broomrape albumin  <dig> xenogene.
the putative cuscuta albumin  <dig> similarly encodes a protein predicted to have knottin structure .
no other albumin  <dig> sequences were identified elsewhere in searches of refseq or publicly available plant transcriptome datasets.
CONCLUSIONS
because of their extensive, intimate contacts with host plant tissues, and the wide range of materials that are commonly transmitted across haustorial connections  <cit> , parasitic plants play an important role as recipients and donors for hgt in plants  <cit>  <cit> .
as parasitic plants increasingly become the targets for genome-scale analyses, it should become possible to estimate the frequency and likely mechanisms of hgt events between parasites and hosts involving albumin  <dig> and other genes, the likelihood of more complex stepping-stone models, and how often hgt leads to long-term maintenance of new genes and novel traits.

 horizontal gene transfer is the nonsexual transmission of genetic material across species boundaries. hgt is well known in bacteria, where hgt often results in adaptive gains of novel genes and traits. there are fewer well-documented cases of hgt among eukaryotes and the large majority of these cases appear to result in short-lived, nonfunctional sequences. consequently, the long-term evolutionary impact of hgt in multicellular eukaryotes remains largely unknown. several cases of hgt are known or suspected in plants, most involving mitochondrial sequences, and/or parasitic plants. parasitic plants form direct haustorial connections with their host plants and are capable of obtaining a wide range of macromolecules from their hosts, including viruses, gene silencing signals, and messenger rnas. consequently, parasites may have many opportunities for hgt events and an increased likelihood that some of these result in functional, and potentially adaptive, gene transfers. two recent reports by yoshida et al and xi et al were the first indications that nuclear protein coding sequences, likely obtained from their respective host species, could be integrated into the genomes of parasitic plants by hgt. these were important advances, but they provided few clues as to the long term impact of hgt, how the transgenes evolve, and how they may function. we hypothesized that systematic analysis of genome-scale datasets from parasitic plants could lead to evidence for acquisition and long-term maintenance of functional gene sequences in plants that had been acquired via hgt. albumin genes are known only from a subset of species in the legume family of angiosperms where they encode seed storage proteins and insect toxins. the albumin proteins in legumes are to amino acids in length and rich in cysteine residues. they form a unique protein structure known as a knottin, which has three disulfide bonds and is characterized by a “disulfide through disulfide knot”. the knottins are famous for their intruguing “disulfide through disulfide knot” structure and have been extensively studied in various fields, most of which are related with potentials in drug design. albumin genes may have originated early in the diversification of papilionoid legumes, but multiple homologous gene copies have been found only in species that are members of the more derived “millettioid s.l.” and “hologalegina” clades. orobanche s. l., often known by the common name “broomrape,” includes 150- obligate parasitic plant species in the family orobanchaceae. growing evidence supports the segregation of broomrapes into four genera: aphyllon, myzorrhiza, phelipanche, orobanche s. str.. most broomrape species have a narrow host spectrum and grow exclusively on perennial eudicot host plants, with members of the leguminosae, solanaceae, and asteraceae among the more common hosts. as a member of order lamiales, orobanchaceae is phylogenetically well-separated from host members in these lineages, particularly legume hosts in the rosid order fabales. a few broomrape species have become devastating pests of important crop plants, affecting their growth and resource allocation and imparting significant losses in yield. p. aegyptiaca, the focal species in this study, has a broad host range that includes members of the eudicot families apiaceae, asteraceae, brassicaceae, cucurbitaceae, leguminosae, and solanaceae. here we show that a gene encoding albumin knottin-like protein, closely related to the albumin genes, only known from papilionoid legumes, serving dual roles in food storage and as insect toxins, was found in phelipanche aegyptiaca and related parasitic species of family orobanchaceae, and was likely acquired by a phelipanche ancestor via hgt from a legume host based on phylogenetic analyses. according to genomic sequences from nine related parasite species, 3d protein structure simulation tests, and evolutionary constraint analyses, the broomrape xenogene we identified here retains the intron structure, six highly conserved cysteine residues necessary to form a knottin protein, and displays levels of purifying selection like those seen in legumes. the albumin xenogene has evolved through > speciation events over ca. million years, forming a small family of differentially expressed genes that may confer novel functions in the parasites.  the albumin transcript was first identified as a hgt candidate in the transcriptome of p. aegyptiaca using a blast-based bioinformatic screen. albumin transcripts were then searched further, using blastx, against the ncbi nr database and the plantgdb database. top hits were seen to medicago truncatula albumin sequences, with expected values of 5e- and 1e- additional blast, including hidden markov model -based psi-blast searches with the sequence from p. aegyptiaca were performed to attempt to detect homologs in three other members of orobanchaceae with large transcriptome datasets. several large public databases, including phytozome, plantgdb, and sol genomics network, were also searched. after searching sequenced genomes and transcriptomes of additional plant species, albumin homologs were detected only in legumes and the transcriptome libraries of p. aegyptiaca. having identified the albumin sequence in the p. aegyptiaca transcriptome, genomic sequences encoding albumin were then obtained from p. aegyptiaca and eight additional broomrape species, including p. schultzii, p. ramosa, p. mutelli, p. nana, and orobanche hederae, o. minor, o. cernua and o. ballotae. the nucleotide sequence and inferred gene structures of the albumin genes in broomrape species are closely comparable, with inferred protein alignments - % identical and %- % similar in ungapped regions between the legume and parasite proteins. two albumin genes were identified in phelipanche species, and are identified here as copy_ and copy_ or albumin1- and albumin1- respectively. an intron disrupts the coding region at the same position in both genes and the intron sequences are similar but contain a number of insertion and deletion mutations. only one albumin gene was detected from orobanche species. although the intron length in albumin genes of phelipanche and legume species is not well conserved, several critical intron features are shared. first, the starting position of the intron in both the p. aegyptiaca and m. truncatula sequences are the same, and the first nine base pairs are identical. second, the introns have characteristic splicing sites at their 5’ and 3’ ends; 5’ ends often have gt/gu and 3’ ends often have ag, and these motifs are found in both m. truncatula and phelipanche albumin introns. albumin gene sequences from phelipanche were also searched with blastn against the ncbi nt database in order to search for high frequency repeats and mobile elements, but no such features were identified. phylogenetic analysis of all known plant albumin sequences showed a strongly supported clade containing all of the albumin sequences from broomrapes boostrap bayesian inference posterior probabilities ) nested deeply within the irlc of papilionoid legumes. among legumes, the next most closely related sequences are from onobrychis argentea and onobrychis viciifolia. because the node supporting the position of the broomrape clade within the papilionoid legumes is relatively weakly supported, we also tested the hypothesis that the broomrape clade of albumin sequences falls outside the larger clade of legumes represented in this analysis. this hypothesis was rejected relative to the maximum likelihood position as indicated in this tree. two albumin genes are resolved as sister clades in phelipanche species, which are in turn resolved as sister to the single gene obtained from orobanche species. gene structures supported a similar conclusion. the amino acid sequence alignments of albumin from legumes to p. aegyptiaca show conservation of all cysteine residues essential for disulfide bond formation in albumin proteins. we investigated whether the predicted albumin proteins from parasites maintain the characteristic knottin structures found in the legume albumin proteins using knoter1d. simulated 3d structures show that the phelipanche albumin proteins form a characteristic knottin structure with three-disulfide bonds and a “disulfide through disulfide knot”. knottin protein structures are also predicted in all of the other full-length albumin genes in phelipanche species. knoter1d assigned scores greater than to each phelipanche albumin sequence; a score greater than in this analysis passes the knoter1d criteria for identification as an albumin structure. the predicted 3d structures for p_aegyptiaca_albumin1- and p_aegyptiaca_albumin1- are very similar to the insect toxic albumin protein from m. truncatula. albumin a non-knottin legume protein, has no discernable homology with the albumin protein in legumes. having found that the horizontally acquired albumin genes were present in related species of broomrapes we then asked if the genes are evolving under purifying selection indicative of a functional protein coding sequence. dn, ds and dn/ds were calculated for all three lineages of the broomrape albumin clade and for the albumin sequences from three closely related legume sequences; astragalus monspessulanus, onobrychis argentea and onobrychis viciifolia. synonymous substitutions in the albumin genes outnumber non-synonymous substitutions by at least 3: in most lineages, and dn/ds, reflecting the level of purifying selection, is similar in broomrapes to the value estimated for closely related albumin sequences from legumes. all cysteine residues were also identified as evolving under purifying selection, suggesting that the horizontally acquired albumin genes in broomrapes are functional having observed evidence for selection for structural conservation, we investigated whether these genes exhibit transcription profiles that suggest a new or unique pattern of expression in parasites. normalized expression levels of both albumin genes in p. aegyptiaca were estimated as reads per kilobase per million reads for eight libraries representing major stages of belowground and aboveground parasite development. both genes displayed lowest expression levels at stage and highest at stage. transcripts were particularly abundant at stage, more than 1000x higher than the haustorial stage. discussion biogeographic overlap and common feeding interactions between diverse broomrapes and temperate papilionoid legumes increase the likelihood that the hgt event occurred in a common ancestor of the parasites that was in direct contact with legume host plants. an alternative explanation is that another organism or virus that co-occurred in the same habitats as the ancestral lineages served as a “stepping stone” for a two- or more step transfer. however, this is not supported in our searches of the sequence databases. based on fossil-calibrated age estimates of legume lineages, we estimate that this horizontal acquisition occurred in an ancestral broomrape that lived in the miocene epoch, about mya. both parasite and their legume host groups have northern temperate distributions, and their lineages likely overlapped in the past as they do now, providing a minimal requirement for a horizontal gene transfer to occur. another possibility, however unlikely, is that albumin was a more recent acquisition that underwent strong convergence at the protein level with this legume lineage. however the branch lengths we observed in the phylogeny were not unusually long in our dna-based phylogeny, and given the large collection of related sequences we obtained from other broomrape species, we have reduced any tendency the orobanche/phelipanche lineage may have had to connect by chance to a deep branch. thus, the convergence hypothesis is not supported. because the breadth of phelipanche and orobanche species we have sampled spans the deepest branches of broomrape diversity, the albumin gene can be inferred to have survived through an extended evolutionary radiation of at least species or more, if the number of now-extinct broomrape species could be estimated. because the introns of phelipanche albumin xenogenes maintain critical splicing sites and share the same starting positions and first nine base pairs with the known m. truncatula albumin intron, it is likely that the hgt event in broomrapes involved transfer of a genomic sequence rather than a separate cdna. following the transfer, albumin genes in broomrape species have evolved under purifying selection consistent with what is observed in related legume albumin genes. this observation, as well as the stage-specific transcription patterns, conserved cysteine residues and predicted 3d knottin protein structures, strongly suggest that albumin genes encode functional proteins in broomrape species, and could potentially serve a function similar to its role in legumes, providing a large pool of sulfur storage and exhibiting toxicity to insect herbivores in certain legumes. a recent report involves panicoid grass species with c or c photosynthetic pathways. evidence was presented that nuclear genes were horizontally transferred between panicoid species and were subsequently adapted into the existing pathways with the effect of advancing the extent of c photosynthesis in some lineages. these results indicate that hgt may promote the sharing of adaptive traits among related species. in comparison, the albumin example described here shows how a completely novel and highly specialized trait has been acquired at an ancestral stage from a distantly related donor species and maintained by the recipient lineage throughout an extended period of evolutionary history. the albumin genes in p. aegyptiaca are highly transcribed in most of the developmental stages we examined. transcripts are more abundant in reproductive tissue, and lowest in the young haustorium, which represents the earliest point in our tissue sampling where the parasite is in direct contact with the host plant. this suggests that the novel gene in p. aegyptiaca is probably not encoding a protein that is playing a direct role in the process of haustorial formation, and that albumin expression is down-regulated as the parasite devotes energy to the essential process of establishing host vascular connections. it is also possible that the low expression in the haustorial stage could help the parasite avoid detection or minimize a negative impact on the health of the host plant during early stages of parasite contact and feeding. several other parasitic lineages, including members of cuscuta, cassytha, apodanthaceae, hydnoraceae, and the order santalales, regularly feed upon legumes and therefore might also have had opportunities to acquire albumin sequences through hgt. large transcriptome datasets are currently available for only two of these, the generalist parasite cuscuta pentagona and for the legume specialist feeder pilostyles thurberi. both of these parasites, and other species in these genera, feed widely on legumes. no homolog of albumin was detected in blast searches of the pilostyles transcriptome in the 1kp dataset. however, albumin sequences were detected in the same dataset and in two additional transcriptome libraries from cuscuta pentagona. phylogenetic analysis nests the cuscuta sequences well within leguminosae, but on an independent branch from the broomrape sequences, suggesting that these transcripts in cuscuta represent a different hgt event into cuscuta from a lineage of papilionoid legumes that was different from the source of the broomrape albumin xenogene. the putative cuscuta albumin similarly encodes a protein predicted to have knottin structure. no other albumin sequences were identified elsewhere in searches of refseq or publicly available plant transcriptome datasets.  because of their extensive, intimate contacts with host plant tissues, and the wide range of materials that are commonly transmitted across haustorial connections, parasitic plants play an important role as recipients and donors for hgt in plants. as parasitic plants increasingly become the targets for genome-scale analyses, it should become possible to estimate the frequency and likely mechanisms of hgt events between parasites and hosts involving albumin and other genes, the likelihood of more complex stepping-stone models, and how often hgt leads to long-term maintenance of new genes and novel traits. 
as a result of their direct plant-plant connections with their host plant, parasitic plants have special opportunities for horizontal gene transfer, the nonsexual transmission of genetic material across species boundaries. there is increasing evidence that parasitic plants have served as recipients and donors of horizontal gene transfer, but the long-term impacts of eukaryotic hgt in parasitic plants are largely unknown. parasitic plants, represented by several thousand species of angiosperms, use modified structures known as haustoria to tap into photosynthetic host plants and extract nutrients and water. moreover, further data show that a distantly related parasitic plant, cuscuta, obtained two copies of albumin knottin-like genes from legumes through a separate hgt event, suggesting that legume knottin structures have been repeatedly co-opted by parasitic plants. here we show that a gene encoding albumin knottin-like protein, closely related to the albumin genes only known from papilionoid legumes, where they serve dual roles as food storage and insect toxin, was found in phelipanche aegyptiaca and related parasitic species of family orobanchaceae, and was likely acquired by a phelipanche ancestor via hgt from a legume host based on phylogenetic analyses. genomic sequences from nine related parasite species were obtained, and 3d protein structure simulation tests and evolutionary constraint analyses were performed. the parasite gene we identified here retains the intron structure, six highly conserved cysteine residues necessary to form a knottin protein, and displays levels of purifying selection like those seen in legumes. million years, forming a small family of differentially expressed genes that may confer novel functions in the parasites.
 parasitic plants, represented by several thousand species of angiosperms, use modified structures known as haustoria to tap into photosynthetic host plants and extract nutrients and water. as a result of their direct plant-plant connections with their host plant, parasitic plants have special opportunities for horizontal gene transfer, the nonsexual transmission of genetic material across species boundaries. there is increasing evidence that parasitic plants have served as recipients and donors of horizontal gene transfer, but the long-term impacts of eukaryotic hgt in parasitic plants are largely unknown.  here we show that a gene encoding albumin knottin-like protein, closely related to the albumin genes only known from papilionoid legumes, where they serve dual roles as food storage and insect toxin, was found in phelipanche aegyptiaca and related parasitic species of family orobanchaceae, and was likely acquired by a phelipanche ancestor via hgt from a legume host based on phylogenetic analyses. the knottins are well known for their unique “disulfide through disulfide knot” structure and have been extensively studied in various contexts, including drug design. genomic sequences from nine related parasite species were obtained, and 3d protein structure simulation tests and evolutionary constraint analyses were performed. the parasite gene we identified here retains the intron structure, six highly conserved cysteine residues necessary to form a knottin protein, and displays levels of purifying selection like those seen in legumes. the albumin xenogene has evolved through > speciation events over ca. million years, forming a small family of differentially expressed genes that may confer novel functions in the parasites. moreover, further data show that a distantly related parasitic plant, cuscuta, obtained two copies of albumin knottin-like genes from legumes through a separate hgt event, suggesting that legume knottin structures have been repeatedly co-opted by parasitic plants.  the hgt-derived albumins in phelipanche represent a novel example of how plants can acquire genes from other plants via hgt that then go on to duplicate, evolve, and retain the specialized features required to perform a unique host-derived function. parasitic plantshorizontal gene transferphelipancheorobanchelegumeknottinalbumin 1evolution 