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.

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.

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.

[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.

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.