Full-length transcriptome analysis of Phytolacca americana and its congener P. icosandra and gene expression normalization in three Phytolaccaceae species.

BACKGROUND: Phytolaccaceae species in China are not only ornamental plants but also perennial herbs that are closely related to human health. However, both large-scale full-length cDNA sequencing and reference gene validation of Phytolaccaceae members are still lacking. Therefore, single-molecule real-time sequencing technology was employed to generate full-length transcriptome in invasive Phytolacca americana and non-invasive exotic P. icosandra. Based on the transcriptome data, RT-qPCR was employed to evaluate the gene expression stability in the two plant species and another indigenous congener P. acinosa. RESULTS: Total of 19.96 Gb and 19.75 Gb clean reads of P. americana and P. icosandra were generated, including 200,857 and 208,865 full length non-chimeric (FLNC) reads, respectively. Transcript clustering analysis of FLNC reads identified 89,082 and 98,448 consensus isoforms, including 86,989 and 96,764 high-quality ones. After removing redundant reads, 46,369 and 50,220 transcripts were obtained. Based on structure analysis, total 1675 and 1908 alternative splicing variants, 25,641 and 31,800 simple sequence repeats (SSR) as well as 34,971 and 36,841 complete coding sequences were detected separately. Furthermore, 3574 and 3833 lncRNA were predicted and 41,676 and 45,050 transcripts were annotated respectively. Subsequently, seven reference genes in the two plant species and a native species P. acinosa were selected and evaluated by RT-qPCR for gene expression analysis. When tested in different tissues (leaves, stems, roots and flowers), 18S rRNA showed the highest stability in P. americana, whether infested by Spodoptera litura or not. EF2 had the most stable expression in P. icosandra, while EF1-α was the most appropriate one when attacked by S. litura. EF1-α showed the highest stability in P.acinosa, whereas GAPDH was recommended when infested by S. litura. Moreover, EF1-α was the most stable one among the three plant species whenever germinating seeds or flowers only were considered. CONCLUSION: Full-length transcriptome of P. americana and P. icosandra were produced individually. Based on the transcriptome data, the expression stability of seven candidate reference genes under different experimental conditions was evaluated. These results would facilitate further exploration of functional and comparative genomic studies in Phytolaccaceae and provide insights into invasion success of P. americana.

A comparison between constitutive and inducible transgenic expression of the PhRIP I gene for broad-spectrum resistance against phytopathogens in potato.

OBJECTIVES: To engineer broad spectrum resistance in potato using different expression strategies. RESULTS: The previously identified Ribosome-Inactivating Protein from Phytolacca heterotepala was expressed in potato under a constitutive or a wound-inducible promoter. Leaves and tubers of the plants constitutively expressing the transgene were resistant to Botrytis cinerea and Rhizoctonia solani, respectively. The wound-inducible promoter was useful in driving the expression upon wounding and fungal damage, and conferred resistance to B. cinerea. The observed differences between the expression strategies are discussed considering the benefits and features offered by the two systems. CONCLUSIONS: Evidence is provided of the possible impact of promoter sequences to engineer BSR in plants, highlighting that the selection of a suitable expression strategy has to balance specific needs and target species.

Enhanced cytotoxic activity of a bifunctional chimeric protein containing a type 1 ribosome-inactivating protein and a serine protease inhibitor.

Both ribosome-inactivating proteins (RIPs) and plant proteinase inhibitors, belong to protein families known to regulate cellular homeostasis and likely involved in plant defense. Nevertheless the interest in these protein classes is due to their potential use for the treatment of several important human diseases such as cancer. Thus, in the present study, type 1 ribosome-inactivating protein and wheat subtilisin/chymotrypsin inhibitor, were engineered into a chimeric protein with cytotoxic action selective for murine tumor cells, while lacking any appreciable toxicity on murine normal cells. This chimeric protein selectively sensitizes to apoptotic death cells derived from Simian-virus-40-transformed mouse fibroblasts (SVT2 cells). The cytotoxicity of this new recombinant product has been detected also on three different human malignant cells. Therefore action on tumor cells of this protein could represent a potentially very attractive novel tool for anticancer drug design.

The role of the glycan moiety on the structure-function relationships of PD-L1, type 1 ribosome-inactivating protein from P. dioica leaves.

N-glycosylation is one of the major naturally occurring covalent co-translational modifications of proteins in plants, being involved in proteins structure, folding, stability and biological activity. In the present work the influence of carbohydrate moieties on the structure-function relationships of type 1 ribosome-inactivating proteins (RIPs) was investigated. To this aim, PD-Ls, RIPs isolated from Phytolacca dioica L. leaves, differing for their glycosylation degree, were used as an experimental system. In particular, comparative structural and biological analyses were performed using native and unglycosylated recombinant PD-L1, the most glycosylated P. dioica RIP isoform. The glycans influence on protein synthesis inhibition and adenine polynucleotide glycosidase activity was investigated. The interaction with adenine, the product of the de-adenylation reaction, was also investigated for native and recombinant PD-L1 by fluorescence spectroscopy. Furthermore, the crystal structure of PD-L1 in complex with adenine was determined. Our data confirm that the absence of glycan moieties did not affect the biological activity in terms of protein synthesis inhibition. However, the removal of carbohydrate chains significantly increased the deadenylation capability, likely as a consequence of the increased accessibility of substrates to the active site pocket. Furthermore, preliminary data on cellular uptake showed that all PD-L isoforms were internalized and, for the first time, that the vesicular distribution within cells could be influenced by the protein glycosylation degree.

Structural characterization and comparative modeling of PD-Ls 1-3, type 1 ribosome-inactivating proteins from summer leaves of Phytolacca dioica L.

The amino acid sequence and glycan structure of PD-L1, PD-L2 and PD-L3, type 1 ribosome-inactivating proteins isolated from Phytolacca dioica L. leaves, were determined using a combined approach based on peptide mapping, Edman degradation and ESI-Q-TOF MS in precursor ion discovery mode. The comparative analysis of the 261 amino acid residue sequences showed that PD-L1 and PD-L2 have identical primary structure, as it is the case of PD-L3 and PD-L4. Furthermore, the primary structure of PD-Ls 1-2 and PD-Ls 3-4 have 81.6% identity (85.1% similarity). The ESI-Q-TOF MS analysis confirmed that PD-Ls 1-3 were glycosylated at different sites. In particular, PD-L1 contained three glycidic chains with the well known paucidomannosidic structure (Man)(3) (GlcNAc)(2) (Fuc)(1) (Xyl)(1) linked to Asn10, Asn43 and Asn255. PD-L2 was glycosylated at Asn10 and Asn43, and PD-L3 was glycosylated only at Asn10. PD-L4 was confirmed to be not glycosylated. Despite an overall high structural similarity, the comparative modeling of PD-L1, PD-L2, PD-L3 and PD-L4 has shown potential influences of the glycidic chains on their adenine polynucleotide glycosylase activity on different substrates.

Type 1 ribosome-inactivating proteins from Phytolacca dioica L. leaves: differential seasonal and age expression, and cellular localization.

The expression of type 1 ribosome-inactivating proteins (RIPs) in Phytolacca dioica L. leaves was investigated. Fully expanded leaves of young P. dioica plants (up to 3 years old) expressed two novel RIPs, dioicin 1 and dioicin 2. The former was also found in developing leaves from adult P. dioica within about two and a half weeks after leaf development, and the latter continuously synthesized, with no seasonal or ontogenetic constraint. Fully expanded leaves from adult P. dioica expressed four RIPs (PD-Ls1-4) exhibiting seasonal variation. RIPs were localized in the extracellular space, in the vacuole and in the Golgi apparatus of mesophyll cells. Dioicin 1 and dioicin 2 showed rRNA N-beta-glycosidase activity and displayed the following properties, respectively: (1) Mr values of 30,047.00 and 29,910.00, (2) pIs of 8.74 and 9.37, and (3) IC(50) values of 19.74 (0.658 nM) and 6.85 ng/mL (0.229 nM). Furthermore, they showed adenine polynucleotide glycosylase activity and nicked pBR322 dsDNA. The amino acid sequence of dioicin 2 had 266 amino acid residues, and the highest percentage identity (81.6%) and similarity (84.6%) with PAP-II from Phytolacca americana, while its identity with other RIPs from Phytolaccaceae was around 40%.

Isolation and characterization of heterotepalins, type 1 ribosome-inactivating proteins from Phytolacca heterotepala leaves.

Leaves from Phytolacca heterotepala H. Walter (Mexican pokeweed) contain at least 10 type 1 RIP isoforms, named heterotepalins. Their Mr values are included in the range 28,000-36,000, as shown by SDS-PAGE performed under reduced conditions and the pI values in the pH range 8.50-9.50. Some heterotepalins are glycosylated. ESI-QTOF mass spectrometry provides the accurate Mr of heterotepalin 4 (29,326.00) and heterotepalin 5b (30,477.00), two isoforms purified to homogeneity by conventional chromatographic techniques. The N-terminal sequences up to residue 35, show that heterotepalins exhibit an high percentage identity with other type 1 RIPs isolated from Phytolaccaceae. Some heterotepalins cross-react with antisera raised against RIPs isolated from Phytolacca dioica leaves. The complete amino acid sequence of heterotepalin 4 matches that of Phytolacca heterotepala anti-viral protein PAP (RIP1), deduced from the cDNA sequence of PhRIP1 gene (AC: AY327475), with one exception concerning residue 245 which, in the native protein, is Ile instead of Met. This substitution, found by mass spectrometry mapping, has been directly confirmed by Edman degradation sequencing of the C-terminal tryptic peptide 242-262. The results show the high potential of mass spectrometry and Edman degradation to verify and to uncover possible amino acid substitutions between native proteins and their cDNA deduced sequences.

Pokeweed antiviral protein inhibits brome mosaic virus replication in plant cells.

Pokeweed antiviral protein (PAP) is a ribosome-inactivating protein isolated from the pokeweed plant (Phytolacca americana) that inhibits the proliferation of several plant and animal viruses. We have shown previously that PAP and nontoxic mutants of PAP can directly depurinate brome mosaic virus (BMV) RNA in vitro, resulting in reduced viral protein translation. Here we expand on these initial studies and, using a barley protoplast system, demonstrate that recombinant PAP and nontoxic mutants isolated from E. coli are able to reduce the accumulation of BMV RNAs in vivo. Pretreatment of only BMV RNA3 with PAP prior to transfection of barley protoplasts reduced the accumulation of all BMV RNAs, with a more severe effect on subgenomic RNA4 levels. Using in vitro RNA synthesis assays, we show that a depurinated template causes the BMV replicase to stall at the template nucleotide adjacent to the missing base. These results provide new insight into the antiviral mechanism of PAP, namely that PAP depurination of BMV RNA impedes both RNA replication and subgenomic RNA transcription. These novel activities are distinct from the PAP-induced reduction of viral RNA translation and represent new targets for the inhibition of viral infection.

Distinction between wild and cultivated enset (Ensete ventricosum) gene pools in Ethiopia using RAPD markers.

In southwest Ethiopia, the cultivation area of Ensete ventricosum (enset) overlaps with the natural distribution area of this species. Analyses of genetic diversity were undertaken using RAPD to provide information for conservation strategies as well as evidence of possible gene flow between the different gene pools, which can be of interest for future improvement of cultivated enset. The extent of RAPD variation in wild enset was investigated in 5 populations in the Bonga area (Kefficho administrative region) and 9 cultivated clones. Comparisons were also made with some Musa samples of potential relevance for crop improvement. Nine oligonucleotide primers amplified 72 polymorphic loci. Population differentiation was estimated with the Shannon index (G'(ST)=0.10), Nei's G(ST) (0.12) and AMOVA (Phi(ST)=0.12), and appears to be relatively low when compared with outbreeding, perennial species in general. Cluster analysis (UPGMA) and principal component analysis (PCA) similarly indicated low population differentiation, and also demonstrated that cultivated clones essentially clustered distinctly from wild enset samples, suggesting that the present-day cultivated enset clones have been introduced to domestication from a limited number of wild progenitors. In addition, subsequent gene flow between wild and cultivated enset may have been prohibited by differences between modes of propagation and harvesting time; cultivated enset is propagated vegetatively through sucker production and the plant is generally harvested before maturity or flower set, thereby hindering pollination by wild enset or vice versa. A significant correlation was not found between genetic and geographical distances. The relatively high total RAPD diversity suggests that wild enset populations in the Bonga area harbour genetic variability which could potentially act as a source for useful or rare genes in the improvement of cultivated enset. As expected, E. ventricosum was clearly differentiated from the analysed Musa samples, that clustered in accordance with the present morphology- and molecular marker-based taxonomy of the genus.

[Induction and culture of hairy boots in Phytolacca esculenta and its saponin production].

Hairy roots appeared in vitro 20 days after inoculation of Phytolacca esculenta leaf explants with the strain of Agrobacterium rhizogenes R1601. The frequency of leaf explants transformed by R1601 was up to 70%. Hairy roots can be incited directly from the veins of explants or via callus. Hairy roots induced by R1601 grew rapidly on medium MS without hormone and were 85.6% higher in respiration rate than control roots. Transformation was confirmed by opine detection and the amplification of rol B and rol C genes from the hairy roots of P. esculenta. The total saponin content in hairy roots of P. esculenta was about 1.54 times as much as in natural roots whereas polysaccharides content was about 70% times as much as in the natural roots.

Genetic relationships among ten endod types as revealed by a combination of morphological, RAPD and AFLP markers.

The genetic relationships among ten types of endod (Phytolacca dodecandra) cultivated by the Institute of Pathobiology of the Addis Ababa University to combat the disease bilharzia in Ethiopia were studied using morphology and molecular markers. A total of 18 morphological characters, 194 amplified fragment length polymorphism (AFLP) and 42 random amplified polymorphic DNA (RAPD) markers were used to determine genetic proximity between types. Genetic distance and cluster analysis of the AFLP data revealed the lack of genetic difference between E47 and E48 but relatively wider genetic difference among the other endod types. Cluster and principal component analyses performed on the AFLP and RAPD markers demonstrated the presence of distinct separation of E56 but not that of E44 from the others. The AFLP and RAPD data, thcrefore, did not support the hypothesis that the superiority of E44 in agronomic traits and molluscicidal potency is linked to its distinct genetic difference from the other endod types. Matrices correspondence tests demonstrated the presence of greater correspondence between AFLP and RAPD data (r = 0.842) but not between the morphology and that of AFLP and RAPD. This indicates the correspondence more between the two DNA markers systems than either of them with morphological traits. The cophenetic correlation coefficients also revealed poor fit for morphology (r = 0.716), good fit for RAPD (r = 0.872) and very good fit for AFLP (r = 0.975), reflecting the hyper-variability and higher resolving power of AFLP.

Analysis of genetic diversity and structure in Ethiopian populations of Phytolacca dodecandra using RAPD.

The genetic diversity and structure in 17 wild populations (249 individuals) of Phytolacca dodecandra (endod) sampled along altitudinal gradients of 1600-3000 meters above sea level (m.a.s.l.) in Ethiopia was studied using random amplified polymorphic DNA (RAPD). A total of 70 polymorphic loci (P) scored from 12 RAPD primers were used to calculate different diversity indices within and between populations, habitats, geographical regions, climatic zones and altitude groups. The number of polymorphic loci and overall Shannon information measure (H) in the populations varied from 30 to 55 and from 0.228 to 0.418, respectively. In general, differences in population variability were found significantly correlated to effective population size. Both P and H were significantly higher in an undisturbed than in a disturbed habitat, and in the lowland and central-highland than in the highland altitude group. However, for both parameters the differences were not statistically significant between regions and climatic zones. Genetic distance between populations varied from 0.301 to 0.628. Cluster analysis performed using the genetic distance matrix revealed a clear separation of the highland populations (2501-3000 m.a.s.l.) from those of the lowland/central-highlands (1600-2500 m.a.s.l.) irrespective of their geographical regions and climatic zones. Analysis of molecular variance (AMOVA) indicated that differences in habitat, geographical regions and climatic zones explained 4.6%, 2.5% and 4.6%, respectively. But none of these differences were significant. Altitude explained 17.2% of the total variance and was highly significant. The data, therefore, clearly indicated the association of genetic structure in endod with altitude. The proportion of RAPD variation found among populations (21.2-35.0%) was somewhat intermediate between values reported for selfing and outcrossing species. The fixation index (FST) values (0.350 to 0.384) indicated very high genetic differentiation among populations.