Polymorphism and structure of style–specific arabinogalactan proteins as determinants of pollen tube growth in Nicotiana
© The Author(s). 2017
Received: 1 February 2017
Accepted: 3 July 2017
Published: 10 August 2017
Pollen tube growth and fertilization are key processes in angiosperm sexual reproduction. The transmitting tract (TT) of Nicotiana tabacum controls pollen tube growth in part by secreting pistil extensin-like protein III (PELPIII), transmitting-tract-specific (TTS) protein and 120 kDa glycoprotein (120 K) into the stylar extracellular matrix. The three arabinogalactan proteins (AGP) are referred to as stylar AGPs and are the focus of this research. The transmitting tract regulates pollen tube growth, promoting fertilization or rejecting pollen tubes.
The N-terminal domain (NTD) of the stylar AGPs is proline rich and polymorphic among Nicotiana spp. The NTD was predicted to be mainly an intrinsically disordered region (IDR), making it a candidate for protein-protein interactions. The NTD is also the location for the majority of the predicted O-glycosylation sites that were variable among Nicotiana spp. The C-terminal domain (CTD) contains an Ole e 1-like domain, that was predicted to form beta-sheets that are similar in position and length among Nicotiana spp. and among stylar AGPs. The TTS protein had the greatest amino acid and predicted O-glycosylation conservation among Nicotiana spp. relative to the PELPIII and 120 K. The PELPIII, TTS and 120 K genes undergo negative selection, with dn/ds ratios of 0.59, 0.29 and 0.38 respectively. The dn/ds ratio for individual species ranged from 0.4 to 0.9 and from 0.1 to 0.8, for PELPIII and TTS genes, respectively. These data indicate that PELPIII and TTS genes are under different selective pressures. A newly discovered AGP gene, Nicotiana tabacum Proline Rich Protein (NtPRP), was found with a similar intron-exon configuration and protein structure resembling other stylar AGPs, particularly TTS.
Further studies of the NtPRP gene are necessary to elucidate its biological role. Due to its high similarity to the TTS gene, NtPRP may be involved in pollen tube guidance and growth. In contrast to TTS, both PELPIII and 120 K genes are more diverse indicating a possible role in speciation or mating preference of Nicotiana spp. We hypothesize that the stylar AGPs and NtPRP share a common origin from a single gene that duplicated and diversified into four distinct genes involved in pollen-style interactions.
Pollen-pistil interactions are dynamic, complex and spatially differentiated. The pollen tube delivers the male gamete to the female gametophyte, beginning with pollen grain hydration and germination on the stigma. The transmitting tract regulates pollen tube growth, promoting fertilization (plant compatibility) or rejecting pollen tubes (incompatibility). The stigma, style and TT have a role in genetic isolation of plant populations and consequently species evolution . Plants evolved multiple prezygotic mechanisms to control fertilization. Self-incompatibility (SI) is a barrier that helps maintain species genetic variation  and interspecific incompatibility (II) prevents gene flow among species, preserving species genetic integrity.
The highly differentiated TT evolved with enclosed ovules of angiosperms and is a pathway for pollen tube growth from the stigma to the ovules [9, 42]. Pollen tubes grow rapidly  and the fastest growing pollen tubes reach the ovules first giving rise to progeny [4, 16], making pollen tube growth a key step where natural selection may act . The initial rate of pollen tube growth through the style is slower but increases as it grows . This is associated with the transition from autotrophic (nutrients obtained from the pollen grain) to heterotrophic growth (nutrients obtained from transmitting tract; [15, 37, 38, 50]). The final step of pollen tube growth is pollen tube-synergid attraction [31, 33] and finally fertilization [1, 46].
Arabinogalactan proteins are found in the plasma membrane, the cell wall, as well as the apoplastic space of the pollen tube [20, 52, 61] and are involved in many diverse processes . Stylar AGPs are very abundant in the TT extracellular matrix and heterogeneous due to post-translational modifications [23, 56, 69]. The AGPs belong to a family of structurally related glycoproteins/proteoglycans, the Pro/Hyp-rich glycoproteins with attached peripheral sugars that produce large protein diversity . Stylar AGPs have a hydroxyproline rich, highly O-glycosylated NTD and a cysteine-rich CTD [2, 86]. Three Nicotiana spp. AGPs, class III pistil extension-like protein (PELPIII), transmitting tissue-specific proteins (TTS) and 120 kDa protein (120 K), accumulate in the extra cellular matrix, interact with growing pollen tubes, and are developmentally regulated and involved in regulation of pollen tube growth [26, 88]. de Graaf  showed that the N. tabacum PELPIII (pMG15) CTD, in particular the cysteine pattern was highly similar to this of N. alata 120 K, N. alata PELPIII, N. alata GaRSGP, Phaseolus vulgaris PvPRP1, and N. tabacum TTS-1. It was suspected that the PELPIII gene has two exons, but the CTD of the gene was not fully described previously . Current genomic resources of N. tabacum, both ancestral species N. sylvestris and N. tomentosiformis [73, 74] provide the possibility to fully describe intron-exon configuration of stylar AGPs.
The PELPIII protein is incorporated into the pollen tube wall of both compatible and incompatible pollen tubes [10, 13, 26]. Gardner et al.,  produced a transmitting tract ablated line (TT-ablated) of N. tabacum that does not have a mature TT and has greatly reduced accumulation of the stylar AGPs. The TT-ablated line was used as a female in controlled pollinations with several species of Nicotiana. Nicotiana tabacum pollen tube growth occurred, albeit at a slightly reduced rate, suggesting the TT and AGPs are not essential for self-pollen tube growth. However, TT-ablation in N. tabacum did alter interspecific pollen tube growth and was essential for II with N. obtusifolia and N. repanda .  showed that PELPIII was not essential for self N. tabacum pollen tube growth or seed set, but was essential for inhibition of N. obtusifolia and N. repanda pollen tube growth. Eberle et al.,  found that N. obtusifolia and N. repanda pollen tubes grew significantly longer in N. tabacum styles where expression of PELPIII was suppressed. The TTS protein promotes self N. tabacum pollen tube growth in vivo and in vitro and acts as a chemical attractant for N. tabacum pollen tubes [7, 85]. During growth through the style, pollen tubes walls incorporate TTS and de-glycosylate it and possibly use the freed arabinogalactan as a source of energy . The 120 K protein is localized to the lumen and vacuolar membranes in N. alata pollen tubes  and was shown to be required for S-specific pollen rejection . Plants with no detectable PELPIII or 120 K and greatly reduced TTS all set self-seed, although pollen tube growth was reduced in plants with lower TTS accumulation [7, 17, 30]. In higher plants, the S-RNase is the female component of SI and the S-locus F-box is the male interactor . While the 120 K is essential for S-specific pollen rejection, two other stylar AGPs (PELPIII and TTS) were found as S-RNase binding proteins . Despite their abundance and regulatory functions, little is known regarding the specific mechanism of stylar AGP action in relation to pollen tube-style interactions [20, 56].
Arabinogalactan proteins undergo extensive O-glycosylation at hydroxyproline and serine residues [59, 80], which leads to the different molecular weights of stylar AGPs [2, 17, 30, 87]. A difference between PELPIII and TTS is the presence of repeating units of P3–6 in the NTD of PELPIII that are absent in TTS. Those repeats were predicted to be sites of post-translational modifications [2, 13] such as O- and N-glycosylation and may be important to AGP function in regulating pollen tube growth. However, Bosch et al.,  found that PELPIII was not N-glycosylated. The level of O-glycosylation of TTS is higher at the top of the style than at the bottom . The 120 K glycosylation patterns among closely related Nicotiana spp. showed differences in protein molecular weights that may result from differences in protein sequence rather than being the result of differential glycosylation . de Graaf et al.,  concluded that glycoproteins with homologous amino acid sequences may have different functions based on their distinct post-translational O-glycosylation patterns, which can be developmentally and spatially regulated. Algorithms for predicting plant-specific O-glycosylation are not fully developed; however, a consensus amino acid motif of [ASTV]-P(1,4)-X(0,10)-[ASTV]-P(1,4) was proposed by Gomord et al.  and can be useful in predicting O-glycosylation patterns of the stylar AGPs. There is evidence that O-glycosylation, phosphorylation and acetylation (but not N-glycosylation) occurs predominantly in the intrinsically disordered regions (IDRs) of plant proteins, making them hot spots for post-translational modifications [40, 41, 58]. Since O-glycosylation is common and associated with AGP functions [70, 71] it is important to make use of predictions for O-glycosylation sites to better understand AGP-protein interactions and their function in pollen tube-style interactions.
A feature of the stylar AGPs is the conserved CTD and less conserved proline-rich NTD [13, 26]. The cysteine-rich CTD shows high similarity to the conserved Ole e 1 domain [55, 68] that was first identified as the main allergen from olive pollen as well as growing pollen tubes [44, 83]. Pollen Ole e 1 was localized in extracellular space in close proximity of the pollen tube wall . Muschietti et al.,  hypothesized that Ole e 1 proteins participate in pollen tube emergence and guidance based on sequence similarities between Ole e 1 protein from olive and its homolog in tomato, the LAT52 gene. de Dios et al.,  found a significant increase of Ole e 1 protein during and after pollen tube germination. The N. tabacum PELPIII, 120 K and TTS each has a conserved Ole e 1-like domain [26, 68]. Similar to the stylar AGPs, the pollen Ole e 1 protein is glycosylated, resulting in multiple glycosylation variants . The petunia PhPRP1 protein has high similarity to N. alata NaTTS (83%) and N. tabacum TTS-1 (81%) and has six conserved cysteine residues of an Ole e 1-like domain, further confirming conservation of this domain outside of Oleaceae family, and its common presence in Solanaceae.
Reproductive proteins that mediate sexual reproduction by taking part in gamete recognition diverge rapidly due to adaptive evolution [49, 78]. Speciation genes prevent gene flow among populations that can result in the divergence of populations, creation of new species and prevent inbreeding depression. Signatures of natural selection  are identified by comparing orthologous genes and provide insights into adaptation and the processes of speciation . Rapid gene evolution (gene sequence divergence) would be indicative of natural selection acting on a gene. The dn/ds ratio (where: dn = rate of nonsynonymous substitution; ds = rate of synonymous substitution) is a method to quantify how amino acid changes accumulate during the course of evolution . A high dn/ds ratio (above 1) suggests that adaptive evolution has been frequent with a high rate of functional protein divergence arising from positive selection . The gametophytic SI locus, the S locus in Solanaceae encodes a RNase and has a signature of positive selection with a dn/ds ratio greater then 1 . Interspecific incompatibility and SI in Nicotiana spp. act as prezygotic isolation mechanisms . The stylar AGPs, particularly PELPIII and 120 K take part in II and SI, respectively [17, 30] and are excellent candidates to test whether they have undergone positive selection. In contrast, TTS is known to take part in regulation of pollen-tube growth and it could have a distinct dn/ds ratio when compared to PELPIII and 120 K, as proteins involved in II and SI, respectively.
To better understand the role that stylar AGPs play in sexual reproduction, the regulation of pollen tube growth and the mechanisms of reproductive barriers, PELPIII (12 species) and TTS (10 species) cDNAs from phylogenetically diverse Nicotiana spp. were sequenced and analyzed. Newly discovered NtPRP gene was also included to fully describe relationship among stylar AGPs. Nicotiana tabacum ancestral species were added to describe the diversification of these genes that occurred post-hybridization . Due to the overlapping components between II and SI, 120 K sequences were also added to our analysis .
Plant material and sequence source
The species used for coding sequences of the PELPIII included: N. kawakamii (PI# 459106; NCBI sequence accession number: MF278946), N. otophora (PI# 555542; MF278952), N. paniculata (PI# 266380; MF278948), N. repanda (PI# 555551; MF278947), N. rustica (PI# 499174; MF278951), N. setchellii (PI# 555557; MF278950), N. tomentosa (PI# 574525; MF278949), N. stocktonii (PI# 555538; MF278953), TTS: N. clevelandii (PI# 555491; MF278937), N. debenyi (PI# 503320; MF278943) N. kawakamii (PI# 459106; MF278936), N. miersii (PI# 555537; MF278940; partial cDNA without signal sequence region), N. occidentalis (PI# 555541; MF402942), N. rependa (PI# 555551; MF40941), N. setchellii (PI# 555557; MF278942), N. tomentosa (PI# 574525; MF278938), N. velutina (PI# 244630; MF278945), N. paniculata (PI# 266380; MF278939), N. otophora (PI# 555542; MF278944), N. alata  and N. obtusifolia (PI# 555543). Plant material sources were previously described . The above species represent distinct phylogenetic clades of Nicotiana spp. .
Sequences available from the NCBI database for PELPIII are N. alata [U45958.1], N. sylvestris [XM_009798359.1], N. tomentosiformis [XM_009619111.1] and N. tabacum [Z14019.1]; for 120 K: N. alata [U88587.1], N. sylvestris [XM_009798358.1], N. tomentosiformis [XM_009619116.1], N. bonariensis [AY886518.1], N. forgetiana [AY886517.1], N. langsdorffii [AY886516.1.], N. longiflora [AY886513.1], N. plumbaginifolia [AY886512.1] and N. tabacum [AY886511.1]; TTS: N. alata [X70441.1], N. sylvestris [XM_009760521.1], N. tomentosiformis [XM_009604038], N. tabacum [Z16403.1] and [Z16404.1].
The alignment of multiple EST sequences of TTS mRNA from expression analysis studies (NCBI dbEST) showed that Z16403.1 (TTS-1) has an additional cytosine (C) at position 687 bp from the start codon that creates a frame shift. Deletion of C687 restored the open reading frame making Z16403.1 highly similar to the TTS-2 gene [Z16404.1].
RNA isolation and cDNA sequencing
Styles were collected from mature flowers, from plants grown in a temperature controlled greenhouse (average temperature of 24.4 °C) with a photoperiod of 14 h day (supplemental light from metal halide lamps) , soil mix LC8 (Sun Gro® Horticulture), 20 cm nursery pots and stored at −80 °C. Total RNA was extracted from 3 to 5 styles with the RNeasy Plant Mini Kit (Qiagen, Frankfurt, Gremany). The mRNA was eluted with nuclease-free water and stored at −80 °C until used. First-strand cDNA synthesis was performed based on Pinto and Lindblad  using the “combined” method except for a primer (CDS) which was for initiating reverse-transcription instead of a gene specific primer. The cDNA was stored at −20 °C until use. A complete list of primers is shown in Additional file 2: Table S2.
All PCR reactions were performed using Q5 High-Fidelity 2X Master Mix (NEB, Frankfurt, Germany) or Phusion High-Fidelity DNA Polymerase (NEB, Frankfurt, Germany) in a Bio-Rad iCycler Thermal Cycler. Each reaction was prepared according to the user manual for the enzyme or master mix being used. The cycling conditions began with 2 min at 98 °C, followed by 35 cycles of 98 °C for 10 s, 65–68 °C for 30 s plus 72 °C for 30 s and a final extension step at 72 °C for 2 min. A PCR with one gene specific primer and one universal primer at either the 3′-end or 5′-end of the cDNA (3′/5′-RACE) was performed to obtain the full coding sequence information of a gene. A second-round of amplification (nested-PCR) was done using another species- and gene-specific primer.
All PCR products, including 5′ and 3′ RACE were separated on 1% agarose gels and stained with ethidium bromide. The cDNA product was purified using Zymoclean Gel DNA Recovery Kit (Zymo Research Corp, Irvine, USA) and sequenced. Molecular cloning was performed with the use of a NEB PCR Cloning Kit (NEB, Frankfurt, Germany), but only on purified cDNAs that did not produce sequences by direct PCR product sequencing. Sequencing was accomplished by the dideoxy chain termination method by the Genomics Center at the University of Minnesota. The final gene sequence was assembled from multiple (at least three independent) overlapping fragments.
The VSL2B predictor (part of DisProt database of protein disorder) was used to detect protein intrinsically disordered regions [60, 72]. For details on disorder prediction please see Additional file 3: Table S3. Secondary protein structure was predicted using Phyre2 (Protein Homology/analogy Recognition Engine V 2.0; ) with default settings using the intensive modeling mode. Signal sequence prediction was performed using the SignalP 4.1 Server  with the default settings. O-glycosylation patterns were detected using ScanProsite [11, 75] with the amino acid motif [ASTV]-P(1,4)-X(0,10)-[ASTV]-P(1,4) with greedy and no overlap options . Sequence alignment, editing, annotation and manipulation were done using Geneious 8.1.3 software , BioEdit v7.2.5  and MEGA7.0 software . Sequencing files were edited and quality checked using 4Peaks software (http://nucleobytes.com/4peaks/). INDEL diversity and average INDEL length was calculated by DnaSP 5.10 software . Gene names for N. tabacum include their ancestral donor with S and T indicating similarity to N. sylvestris or N. tomentosiformis, respectively. Evolutionary analysis was performed using MEGA 7.0 software and MUSCLE algorithm (default settings) with minor manual adjustments. The evolutionary history was inferred using the Neighbor-Joining method  and the bootstrap test was performed for each tree (500 replicates; ). The trees are drawn to scale, with branch lengths in the same units as those of the evolutionary distances used to infer the phylogenetic tree. The evolutionary distances were computed using the Tamura-Nei method  and are in the units of the number of base substitutions per site. The rate of variation among sites was modeled with a gamma distribution (shape parameter = 2). The analysis used 12 nucleotide sequences of PELPIII and TTS, and 10 nucleotide sequences of 120 K. Codon positions included were 1st + 2nd + 3rd + noncoding. To calculate the dN/dS ratio for PELPIII, TTS and 120 K genes Nei-Gojobori method was used. In both methods, all positions containing gaps and missing data were eliminated. Trees were drawn to scale, with branch lengths in the same units as those of the evolutionary distances used to infer the phylogenetic trees.
Identification of a novel NtPRP gene similar to TTS
Percent identity among NtPRP, PELPIII, 120 K and TTS in TN90 variety of N. tabacum
Stylar AGPs and NtPRP have similar intron-exon configuration
To identify the ancestral origin of each stylar AGP gene and NtPRP, genomic sequences of both N. tabacum ancestral species were analyzed . As expected, two copies of the PELPIII, TTS and NtPRP genes were identified in each of the three N. tabacum genomes. A single copy of the 120 K gene was found, which was most similar to the ancestral N. tomentosiformis. An NCBI EST database search showed a single gene sequence of the 120 K EST, further confirming the existence of only one copy of this gene in the N. tabacum genome. These results suggest that the 120 K-S gene was most likely lost post-hybridization. A single copy of the stylar AGPs and NtPRP genes were found in the N. sylvestris and N. tomentosiformis sequenced genomes.
Stylar AGPs and NtPRP have intrinsically disordered and single globular region
Each stylar AGP and NtPRP contains two predicted intrinsically disordered regions (VSL2B predictor; ). The IDR1 is at the N-terminal, while IDR2 with a single globular domain is located at the C-terminal. The IDR1 is similar in length among all PELPIII and 120 K proteins, but is significantly shorter in TTS and NtPRP proteins. IDR2 is similar in length among all stylar AGPs and NtPRP.
Stylar AGPs and NtPRP are polymorphic among N. tabacum genotypes, N. sylvestris and N. tomentosiformis and Nicotiana spp.
Summary of stylar AGP and NtPRP polymorphisms among N. sylvestris and N. tomentosiformis and BX, TN90, K326 genotypes
G/C45 (L/F), T/C94, INDEL190–201, C/T434 (P/L), G/T516 (P/S), C/T695 (P/L), G/A733 (A/T), A/C790 (T/P), A/C842 (Q/P), G/A978, G/A982 (A/T), G/C998 (G/A), T/G1189(V/F)
T/C816, T/A1097 (I/K)
G/A15, G/A161, G/C190 (A/P), T/C245, C/A264 (N/K), A/G267, G/T288, G/T696, A/T70, A/G732
Cysteine residues are conserved among all PELPIII proteins, one additional cysteine is found in N. otophora at position 119. Predicted sequence of N. clevelandii obtained from genomic DNA aligns with other PELPIII sequences. The N. clevelandii genomic sequence of PELPIII diverges significantly from other PELPIII sequences at amino acid 154. cDNA sequencing produced truncated transcript that was shorter when compared to PELPIII of N. tabacum. However, insertion of two nucleotides restores an amino acid sequence that is similar to other PELPIII. Alignment of PELPIII sequence from N. clevelandii with the additional nucleotides allows extension of the amino acid sequence and restoration of the cysteine residue at a similar position to PELPIII from N. tabacum (Additional file 4: Figure S1). In TTS of N. repanda KPPTKPPTYSPSKPPAKSP sequence is duplicated, additionally with KPPT sequence found in three places near region of duplication. Similarly to PELPIII, INDELs are present in IDR1 and IDR2 of TTS (Additional file 5: Figure S2). Analogous features can be seen in 120 K (; Fig. 2), where multiple INDELs were described. Signal peptides from stylar AGPs were obtained from multiple Nicotiana spp., and are relatively conserved (with minor amino acid polymorphism) among each stylar AGP (Additional files 4 and 5: Figure S1 and S2). However, each signal sequence is characteristic for each stylar AGP.
Multiple INDELs were found in PELPIII among Nicotiana spp., mainly in IDR1 region, with a single amino acid INDEL in the globular region and IDR2 (Additional file 4 Figure S1). In PELPIII of N. setchellii and N. tomentosa the sequence PPPVKAPSPSPAKQP is repeated and the sequence PAKQP was found in three positions in close proximity. A second repeated sequence PSPAKQSPPPP is found twice in N. otophora, but only once in other Nicotiana spp. Similarly, there was a short amino acid sequence PSPA found in three positions near each other in N. otophora and twice in other species.
INDEL polymorphism of PELPIII, TTS and 120 K among Nicotiana spp.
Ole e 1-like
Average INDEL length (bp)
INDEL diversity π(i)
Average INDEL length (bp)
INDEL diversity π(i)
120 K (Hancock et al., 2004)
Average INDEL length (bp)
INDEL diversity π(i)
Stylar AGPs and NtPRP have variable predicted O-glycosylation patterns
Stylar AGPs are very heterogenous proteins due to the high degree of variable post- translational modifications, in particular O-glycosylation. The O-glycosylation is thought to be important in their role as regulators of pollen tube growth [2, 30, 87]. The amino acid motif [ASTV]-P(1,4)-X(0,10)-[ASTV]-P(1,4) was used to predict O-glycosylation sites . O-glycosylation predictions showed variation in IDRs, with IDR2 being much more uniform among the AGPs among Nicotiana spp. No O-glycosylation sites were found in the Ole e 1-like domain in any of the AGPs. The IDR2 region of the TTS proteins were unique, lacking predicted O-glycosylation sites. When comparing predicted O-glycosylation patterns among the same species it is apparent that TTS has a more conserved pattern of predicted O-glycosylation sites than the other AGPs. The predicted O-glycosylation pattern of 120 K resembles that of the PELPIII. The predicted O-glycosylation of NtPRP was conserved and most similar to the TTS gene, however this conclusion is made based on limited sequence data for NtPRP. Relative conservation of TTS and NtPRP O-glycosylation patterns among Nicotiana spp., when compared to PELPIII and 120 K may relate to the known function of TTS as a regulator of pollen tube growth in general.
Stylar AGPs and NtPRP are under negative selection
The reproductive AGPs act during pollen tube growth through the style and can participate in prezygotic barriers that maintain species [3, 17, 65, 76]. Multiple methods exist to estimate signatures of selection that could indicate a possible role of the reproductive AGPs in species diversification. The Nei and Gojobori, 1986 algorithm was used to estimate the signatures of selection in stylar AGPs. The dn/ds ratio analysis provided evidence that stylar AGPs are under negative selection with dn/ds ratios lower than 1. PELPIII and TTS had distinct dn/ds ratios on each branch using both a branch and branch-site analysis. Overall dn/ds ratio for PELPIII (0.59) is higher than that for TTS (0.29) indicating that negative selection acts differently on TTS within the same set of Nicotiana spp. The lowest dn/ds value for PELPIII was N. alata (0.42) and for TTS was N. paniculata (0.09). The highest dn/ds values for PELPIII were N. tabacum-T (0.89) and N. tomentosiformis (0.89) and for TTS was N. tomentosiformis (0.81).  selected multiple 120 K genes from SI and SC Nicotiana spp., based on this selection, dN/dS ratio was calculated using Nei and Gojobori algorithm. The overall dN/dS ratio was 0.38, which indicates that negative selection also took place for 120 K gene (Additional file 6: Figure S3). In summary, AGP gene dn/ds ratio analysis indicate that there has been no positive selection acting on PELPIII, TTS and 120 K genes.
NtPRP and stylar AGPs intron-exon configuration
Stylar AGPs have been studied because of their role in regulating pollen tube growth [2, 5, 7, 14, 48]. The PELPIII is a specific inhibitor of N. obtusifolia and N. repanda pollen tube growth , the TTS protein promotes pollen tube growth in vivo and in vitro  and the 120 K protein is required for N. alata S-specific pollen rejection . Despite the progress in functional analysis of stylar AGPs, little is known about the mechanisms of AGP regulation of pollen tube growth and the relationship among stylar AGPs. Availability of genomic sequence  of N. tabacum allowed discovery of a fourth AGP (named NtPRP) that is similar to the stylar AGPs, and contributes to further understanding of pollen tube-style interactions. The N. tabacum PELPIII, 120 K, TTS and the newly discovered NtPRP genes share a very similar intron-exon configuration with two exons, separated by a variable length intron (Fig. 1). Exon 1 protein sequence lacks conservation and may be important in discriminating the distinct functions of the AGPs with specific pollen tube genotypes. Exon 2 protein sequence, which contains the Ole e 1 like domain, is highly conserved among stylar AGPs and NtPRP and may therefore play an important role in pollen – style biology. The NtPRP may have a similar function to TTS serving as redundancy in promoting pollen tube growth or it may have additional functions given its mRNA accumulation can occur outside the mature style. The mRNA accumulation of NtPRP suggests a potential role of this gene in leaf, seedling and root tissues. Identification of NtPRP provides a new opportunity to investigate its role in pollen tube growth regulation.
Stylar AGPs as interactors
It is possible that variation of the O-glycosylation patterns may be important in the regulation of interspecific incompatibility (PELPIII) or self-incompatibility (120 K). Two species N. repanda and N. alata showed very distinct O-glycosylation patterns of PELPIII. PELPIII inhibits pollen tube growth of N. repanda and N. obtusifolia when grown in N. tabacum styles , indicating PELPIII has a role in II. When compared among species, the predicted O-glycosylation patterns for N. repanda PELPIII and TTS and N. alata 120 K have the most distinct patterns of O-glycosylation across protein, and most of amino acid polymorphism within each glycosylation site (Fig. 4).
The O-glycosylation sites of PELPIII and 120 K among Nicotiana spp. have unique patterns and may be integral to the regulation of pollen tube growth determining whether a pollen-pistil interaction is compatible or incompatible thus influencing II (PELPIII) and SI (120 K). In contrast to these results, the O-glycosylation pattern of TTS is more conserved among Nicotiana spp. The TTS gene from N. alata, N. miersii and N. tomentosa, N. otophora share a similar O-glycosylation pattern. Additionally, there is conservation of O-glycosylation sites among TTS in divergent species. Considering that TTS is known to facilitate pollen tube growth and has no known function in II or SI, the low O-glycosylation variation, suggest the important role of TTS in Nicotiana spp. NtPRP appears to have a level of conservation similar to that of TTS, when compared among N. otophora, N. sylvestris and N. tomentosiformis, but are unique in PELPIII, 120 K and TTS.
Small O-glycosylation differences were found among N. tabacum genotypes. When PELPIII amino acid sequence from N. tabacum genotypes were compared, there was a difference in K326-T predicted O-glycosylation in the first glycosylation site, when compared to the other N. tabacum genotypes originating from the same ancestral donor. Lack of the fifth O-glycosylation site was recognized in BX-S PELPIII, when compared to the TN90 and K326. There is a difference in amino acid sequence polymorphism in case of the first glycosylated site in N. sylvestris, when compared to N. tabacum PELPIII-S proteins. TN90-S and K326-S have small changes in amino acid sequence within predicted O-glycosylation sites, in addition to a deletion in the K326-S that causes a glycosylation site shift in relation to other N. tabacum genotypes. Similarly, INDELs in K326-S protein create difference in the distance of O-glycosylation sites, but not in number of glycosylation sites. Those differences were observed among Nicotiana genotypes that were changed among genotypes due to selection during breeding, that could be consider genetically highly conserved.
Another feature that can be important for protein-protein interactions is the INDEL diversity found in both IDRs of stylar AGPs and possibly NtPRP among Nicotiana spp. (Table 3). The longest and most diverse INDELs are present in 120 K and PELPIII in the IDR1 region among Nicotiana spp. TTS has fewer and shorter INDELs overall and the lowest INDEL diversity when compared to PELPIII and 120 K. The presence of INDELS as long as ~26 aa (PELPIII; on average), changes the distance between O-glycosylation positions in the IDR regions, and may affect the interaction with other proteins.
The Ole e 1-like domain is present in all stylar AGPs and NtPRP and is conserved among other AGPs from many taxa. An Ole e I-like domain is found in the stylar AGP homolog LAT52 protein of Solanum lycopersicum  and was originally discovered in Olea europaea pollen (olive; ). Preservation of the Ole e I-like domain among stylar AGPs and NtPRP suggests a conserved and important biological function and can be a key region of protein-protein interaction .
Hancock et al.  found that N. plumbaginifolia (SC) and N. longiflora (SC) have a 10-amino acid deletion relative to other Nicotiana spp. in the IDR2 (part of CTD) of the 120 K protein and concluded that the deletion would not be expected to inhibit protein folding or function. The predicted pattern of 120 K O-glycosylation (Fig. 4) suggest that the deletion in N. plumbaginifolia and N. longiflora would not change the O-glycosylation pattern relative to the diversity of other 120 K proteins at IDR2 (Fig. 4). However, the 10 amino acids may affect interactions with other proteins that may occur in this region due to the shift O-glycosylation sites.
Stylar AGPs are under negative selection
Stylar AGPs as paralogs
Stylar AGPs and the newly discovered NtPRP share similar intron-exon configuration and secondary structure. The location of the single intron, the presence of two intrinsically disordered regions and the conserved Ole e 1-like domain strongly suggest that the stylar AGPs have a common evolutionary origin. These findings create the possibility to manipulate the composition of stylar AGPs and NtPRP as closely related proteins. Additionally, dn/ds ratios calculated with use of Nei-Gojobori method for PELPIII, TTS and 120 K provide evidence that these genes are under negative selection. Future studies that swap domains among PELPIII, 120 K, TTS and NtPRP will show how the divergent and conserved domains of the AGPs could influence regulation pollen tube growth and II among Nicotiana spp.
The authors wish to thank Jamie Knutson, Marie Sorensen and Jessie Rydeen for their support throughout the research.
This work was funded by the Minnesota Agricultural Experiment Station.
Availability of data and materials
The cDNA sequences produced will be deposited and available in NCBI database Accession numbers for sequences used are located in Material and Methods.
AKN and AGS designed the research; YL and AKN performed the research, AKN analyzed data, KT analyzed data in regard to dn/ds analysis, AKN and AGS drafted and edited the manuscript, AGS managed the project. All authors read and approved the final draft of the manuscript.
Ethics approval and consent to participate
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
- Beale KM, Leydon AR, Johnson MA. Gamete fusion is required to block multiple pollen tubes from entering an Arabidopsis Ovule. Curr Biol. 2012;22:1090–4.View ArticlePubMedPubMed CentralGoogle Scholar
- Bosch M, Knudsen JS, Derksen J, Mariani C. Class III pistil-specific extensin-like proteins from tobacco have characteristics of arabinogalactan proteins. Plant Physiol. 2001;125:2180–8.View ArticlePubMedPubMed CentralGoogle Scholar
- Busch JW, Schoen DJ. The evolution of self-incompatibility when mates are limiting. Trends Plant Sci. 2008;13(3):128–36. doi:10.1016/j.tplants.2008.01.002.View ArticlePubMedGoogle Scholar
- Carlson AL, Gong H, Toomajian C, Swanson RJ. Parental genetic distance and patterns in nonrandom mating and seed yield in predominately selfing Arabidopsis thaliana. Plant Reprod. 2013;26(4):317–28.View ArticlePubMedPubMed CentralGoogle Scholar
- Chen C-G, Mau S-L, Clarke AE. Nucleotide sequence and style-specific expression of a novel proline-rich protein gene from Nicotiana alata. Plant Mol Biol. 1993;21:391–5.View ArticlePubMedGoogle Scholar
- Cheung AY, May B, Kawata EE, Gu Q, Wu HM. Characterization of cDNAs for stylar transmitting tissue-specific proline-rich proteins in tobacco. Plant J. 1993;3(1):151–60.Google Scholar
- Cheung AY, Wang H, Wu H-M. A floral transmitting tissue-specific glycoprotein attracts pollen tubes and stimulates their growth. Cell. 1995;82:383–93.View ArticlePubMedGoogle Scholar
- Clarkson JJ, Knapp S, Garcia VF, Olmstead RG, Leitch AR, Chase MW. Phylogenetic relationships in Nicotiana (Solanaceae) inferred from multiple plastid DNA regions. Mol Phylogenet Evol. 2004;33(1):75–90.View ArticlePubMedGoogle Scholar
- Crawford B, Ditta G, Yanofsky M. The NTT transmitting tract gene is required for transmitting-tract development in carpels of Arabidopsis thaliana. Curr Biol. 2007;17:1101–8.View ArticlePubMedGoogle Scholar
- Cruz-Garcia F, Hancock NC, Kim D, McClure B. Stylar glycoproteins bind to S-RNase in vitro. Plant J. 2005;42(3):295–304.View ArticlePubMedGoogle Scholar
- de Castro E, Sigrist CJA, Gattiker A, Bulliard V, Langendijk-Genevaux PS, Gasteiger E, Bairoch A, Hulo N. ScanProsite: detection of PROSITE signature matches and ProRule-associated functional and structural residues in proteins. Nucleic Acids Res. 2006;34(Web Server issue):W362–5.View ArticlePubMedPubMed CentralGoogle Scholar
- de Dios AJ, M’rani-Alaoui M, Jesús Castro A, Rodríguez-García MI. Ole e 1, the major allergen from olive (Olea europaea L.) pollen, increases its expression and is released to the culture medium during in vitro germination. Plant Cell Physiol. 2004;45(9):1149–57.View ArticleGoogle Scholar
- de Graaf BHJ, Knuiman BA, Derksen J, Mariani C. Characterization and localization of the transmitting tissue-specific PELPIII proteins of Nicotiana tabacum. J Exp Bot. 2003;54(380):55–63.View ArticlePubMedGoogle Scholar
- de Graaf BHJ. Pistil proline-rich proteins in Nicotiana tabacum. Their involvement in pollen-pistil interaction, PhD thesis. Nijmegen: Catholic University of Nijmegen; 1999.Google Scholar
- de Nettancourt D. Incompatibility and incongruity in wild and cultivated plants. eBook. Berlin: Springer; 2010.Google Scholar
- Delph LF, Weinig C, Sullivan K. Why fast-growing pollen tubes give rise to vigorous progeny: the test of a new mechanism. Proc R Soc Lond. 1998;265:935–9.View ArticleGoogle Scholar
- Eberle CA, Anderson NO, Clasen BM, Hegeman AD, Smith AG. PELPIII: the class III pistil-specific extensin-like Nicotiana tabacum proteins are essential for interspecific incompatibility. Plant J. 2013;74(4):805–14.View ArticlePubMedGoogle Scholar
- Eberle CA, Clasen BM, Anderson NO, Smith AG. A novel pollen tube growth assay utilizing a transmitting tract-ablated Nicotiana tabacum style. Sex Plant Reprod. 2012;25:27–37.View ArticlePubMedGoogle Scholar
- Ellegren H. Comparative genomics and the study of evolution by natural selection. Mol Ecol. 2008;17(21):4586–96.View ArticlePubMedGoogle Scholar
- Ellis M, Egelund J, Schultz CJ, Bacic A. Arabinogalactan-proteins: key regulators at the cell surface? Plant Physiol. 2010;153(2):403–19.View ArticlePubMedPubMed CentralGoogle Scholar
- Felsenstein J. Phylogenies and the comparative method. Am Nat. 1985;125(1):1–15.View ArticleGoogle Scholar
- Ferris PJ, Goodenough UW. Mating type in Chlamydomonas is specified by mid, the minus-dominance gene. Genetics. 1997;146:859–69.PubMedPubMed CentralGoogle Scholar
- Fincher GB, Stone BA, Clarke AE. Arabinogalactan-proteins: structure, biosynthesis, and function. Annu Rev Plant Physiol. 1983;34:47–70.View ArticleGoogle Scholar
- Finn RD, Bateman A, Clements J, Coggill P, Eberhardt RY, Eddy SR, Heger A, Hetherington K, Holm L, Mistry J, Sonnhammer ELL, Tate J, Punta M. Pfam: the protein families database. Nucleic Acids Res. 2014;42(D1):D222–30.View ArticlePubMedGoogle Scholar
- Gardner N, Felsheim R, Smith AG. Production of male- and female-sterile plants through reproductive tissue ablation. J Plant Physiol. 2009;166(8):871–81.View ArticlePubMedGoogle Scholar
- Goldman MHS, Pezzotti M, Seurinck J. Mariani C. Developmental expression of tobacco pistil-specific genes encoding novel extensin-like proteins. Plant Cell 1992;4:1041–1051.Google Scholar
- Goldraij A, Kondo K, Lee CB, Hancock CN, Sivaguru M, Vazquez-Santana S, Kim S, Phillips TE, Cruz-Garcia F, McClure B. Compartmentalization of S-RNase and HT-B degradation in self-incompatible Nicotiana. Nature. 2006;439:805–10.View ArticlePubMedGoogle Scholar
- Gomord V, Fitchette AC, Menu-Bouaouiche L, Saint-Jore-Dupas C, Plasson C, Michaud D, Faye L. Plant-specific glycosylation patterns in the context of therapeutic protein production. Plant Biotechnol J. 2010;8(5):564–87.View ArticlePubMedGoogle Scholar
- Hall TA. BioEdit: a user-friendly biological sequence alignment program for windows 95/98/NT. Nucleic Acids Symposium. 1999;41:95–8.Google Scholar
- Hancock CN, Kent L, McClure BA. The stylar 120 kDa glycoprotein is required for S-specific pollen rejection in Nicotiana. Plant J. 2005;43(5):716–23.View ArticlePubMedGoogle Scholar
- Higashiyama T, Yabe S, Sasaki N, Nishimura Y, Miyagishima S. Pollen tube attraction by the synergid cell. Science. 2001;293:1480–3.View ArticlePubMedGoogle Scholar
- Hsu WL, Oldfield C, Meng J, Huang F, Xue B, Uversky VN, Romero P, Dunker AK. Intrinsic protein disorder and protein-protein interactions. Pac Symp Biocomput. 2012:116–27.Google Scholar
- Kanaoka MM, Higashiyama T. Peptide signaling in pollen tube guidance. Curr Opin Plant Biol. 2015;28:127–36.View ArticlePubMedGoogle Scholar
- Kearse M, Moir R, Wilson A, Stones-Havas S, Cheung M, Sturrock S, Buxton S, Cooper A, Markowitz S, Duran C, Thierer T, Ashton B, Meintjes P, Drummond A. Geneious basic. An integrated and extendable desktop software platform for the organization and analysis of sequence data. Bioinformatics. 2012;28(12):1647–9.View ArticlePubMedPubMed CentralGoogle Scholar
- Kelley LA, Sternberg MJ. Protein structure prediction on the web: a case study using the Phyre server. Nat Protoc. 2009;4(3):363–71.View ArticlePubMedGoogle Scholar
- Kimura M. Preponderance of synonymous changes as evidence for the neutral theory of molecular evolution. Nature. 1977;267(5608):275–6.View ArticlePubMedGoogle Scholar
- Knox RB. Pollen–pistil interactions. Encycl Plant Physiol. 1984;17:506–608.Google Scholar
- Kroh M, Miki-Hirosige H, Rosen W, Loewus F. Distribution and utilization of label from myoinositol-U 14C and −2-3H by detached flowers and pistils of Lilium longiflorum. Plant Physiol. 1970;45:86–91.View ArticlePubMedPubMed CentralGoogle Scholar
- Kumar S, Stecher G, Tamura K. MEGA7. Molecular evolutionary genetics analysis version 7.0 for bigger datasets. Mol Biol Evol. 2016;33(7):1870–4.View ArticlePubMedGoogle Scholar
- Kurotani A, Sakurai T. In Silico analysis of correlations between protein disorder and post-translational modifications in algae. Int J Mol Sci. 2015;16:19812–35.Google Scholar
- Kurotani A, Tokmakov AA, Kuroda Y, Fukami Y, Shinozaki K, Sakurai T. Correlations between predicted protein disorder and post-translational modifications in plants. Bioinformatics. 2014;30:1095–103.Google Scholar
- Labarca C, Kroh M, Loewus F. The composition of stigmatic exudate from Lilium longiflorum. Labelling studies with myo-inositol, D-glucose and L-proline. Plant Physiol. 1970;46:150–6.View ArticlePubMedPubMed CentralGoogle Scholar
- Lankinen A, Maad J, Ambruster WS. Pollen-tube growth rates in Collinsia heterophylla (Plantaginaceae): one-donor crosses reveal heritability but no effect on sporophytic-offspring fitness. 2009;103(6):941–950.Google Scholar
- Lauzurica P, Maruri N, Galocha B, Gonzalez J, Diaz R, Palomino P, Hernandez D, Garcia R, Lahoz C. Olive (Olea europea) pollen allergens--II. Isolation and characterization of two major antigens. Mol Immunol. 1988;25(4):337–44.View ArticlePubMedGoogle Scholar
- Leitch IJ, Hanson L, Lim KY, Kovarik A, Chase MW, Clarkson JJ, Leitch AR. The ups and downs of genome size evolution in polyploid species of Nicotiana (Solanaceae). Ann Bot. 2008;101(6):805–14.View ArticlePubMedPubMed CentralGoogle Scholar
- Leydon AR, Tsukamoto T, Dunatunga D, Qin Y, Johnson MA, Palanivelu R. Pollen tube discharge completes the process of synergid degeneration that is initiated by pollen tube-synergid interaction in Arabidopsis. Plant Physiol. 2015;169:485–96.View ArticlePubMedPubMed CentralGoogle Scholar
- Librado P, Rozas J. DnaSP v5. A software for comprehensive analysis of DNA polymorphism data. Bioinformatics. 2009;25:1451–2.View ArticlePubMedGoogle Scholar
- Lind JL, Bacic A, Clarke AE, Anderson MA. A style specific hydroxyproline-rich glycoprotein with properties of both extensins and arabinogalactan proteins. Plant J. 1994;6:491–502.View ArticlePubMedGoogle Scholar
- Lipinska AP, Van Damme EJM, Clerck O. Molecular evolution of candidate male reproductive genes in the brown algal model Ectocarpus. BMC Evol Biol. 2016;16:5.View ArticlePubMedPubMed CentralGoogle Scholar
- Lord EM, Sanders LC. Roles for the extracellular matrix in plant development and pollination: a special case of cell movement in plants. Dev Biol. 1992;153:16–28.View ArticlePubMedGoogle Scholar
- Lowry DB, Modliszewski JL, Wright KM, Wu CA, Willis JH. The strength and genetic basis of reproductive isolating barriers in flowering plants. Philos Trans R Soc Lond Ser B Biol Sci. 2008;363(1506):3009–21.View ArticleGoogle Scholar
- Majewska-Sawka A, Nothnagel EA. The multiple roles of arabinogalactan proteins in plant development. Plant Physiol. 2000;122(1):3–10.View ArticlePubMedPubMed CentralGoogle Scholar
- Mascarenhas JP. Molecular mechanisms of pollen tube growth and differentiation. Plant Cell. 1993;5:1303–14.View ArticlePubMedPubMed CentralGoogle Scholar
- Murfett J, Atherton TL, Mou B, Gasser CS, McClure BA. S-RNase expressed in transgenic Nicotiana causes S-allele-specific pollen rejection. Nature. 1994;367:563–6.View ArticlePubMedGoogle Scholar
- Muschietti J, Dircks L, Vancanneyt G, McCormick S. LAT52 protein is essential for tomato pollen development: pollen expressing antisense LAT52 RNA hydrates and germinates abnormally and cannot achieve fertilization. Plant J. 1994;6:321–38.View ArticlePubMedGoogle Scholar
- Nguema-Ona E, Coimbra S, Vicré-Gibouin M, Mollet JC, Driouich A. Arabinogalactan proteins in root and pollen-tube cells: distribution and functional aspects. Ann Bot. 2012;110(2):3383–404.View ArticleGoogle Scholar
- Nielsen R. Molecular signatures of natural selection. Annu Rev Genet. 2005;39(1):651.View ArticleGoogle Scholar
- Nishikawa I, Nakajima Y, Ito M, Fukuchi S, Homma K, Nishikawa K. Computational prediction of O-linked glycosylation sites that preferentially map on intrinsically disordered regions of extracellular proteins. Int J Mol Sci. 2010;11:4992–5009.View ArticleGoogle Scholar
- Nothnagel EA. Proteoglycans and related components in plant cells. Int Rev Cytol. 1997;174:195–291.View ArticlePubMedGoogle Scholar
- Obradovic Z, Peng K, Vucetic S, Radivojac P, Dunker AK. Exploiting heterogeneous sequence properties improves prediction of protein disorder. Proteins. 2005;61(7):176–82.View ArticlePubMedGoogle Scholar
- Pereira AM, Pereira LG, Coimbra S. 2015. Arabinogalactan proteins: rising attention from plant biologists. Plant Reprod. 2015;28(1):1–15.View ArticlePubMedGoogle Scholar
- Petersen TN, Brunak S, von Heijne G, Nielsen H. SignalP 4.0. Discriminating signal peptides from transmembrane regions. Nat Methods. 2011;8(10):785–6.View ArticlePubMedGoogle Scholar
- Pinto FL, Lindblad P. A guide for in-house design of template-switch-based 5′ rapid amplification of cDNA ends systems. Anal Biochem. 2010;397(2):227–32.View ArticlePubMedGoogle Scholar
- Richman AD, Kohn JR. Evolutionary genetics of self-incompatibility in the Solanaceae. Plant Mol Biol. 2000;42:169–79.View ArticlePubMedGoogle Scholar
- Rieseberg LH, Willis JH. Plant speciation. Science (New York, NY). 2007;317(5840):910–4.View ArticleGoogle Scholar
- Rodríguez R, Villalba M, Batanero E, González EM, Monsalve RI, Huecas S, Tejera ML, Ledesma A. Allergenic diversity of the olive pollen. Allergy. 2002;57(Suppl. 71):6–16.View ArticlePubMedGoogle Scholar
- Saitou N, Nei M. The neighbor-joining method: a new method for reconstructing phylogenetic trees. Mol Biol Evol. 1987;4(4):406–25.PubMedGoogle Scholar
- Schultz CJ, Hauser K, Lind JL. Molecular characterization of a cDNA sequence encoding the backbone of a style-specific 120-kDa glycoprotein which has features of both extensins and arabinogalactan proteins. Plant Mol Biol. 1997;35:833–45.View ArticlePubMedGoogle Scholar
- Schultz CJ, Johnsonb KL, Currieb G, Bacica A. The classical arabinogalactan protein gene family of Arabidopsis. Plant Cell. 2000;12(9):1751–67.View ArticlePubMedPubMed CentralGoogle Scholar
- Seifert GJ, Roberts K. The biology of arabinogalactan proteins. Annu Rev Plant Biol. 2007;58:137–61.View ArticlePubMedGoogle Scholar
- Showalter AM. Arabinogalactan-proteins: structure, expression and function. Cell Mol Life Sci. 2001;58(10):1399–417.View ArticlePubMedGoogle Scholar
- Sickmeier M, Hamilton JA, LeGall T, Vacic V, Cortese MS, Tantos A, Szabo B, Tompa P, Chen J, Uversky VN, Obradovic Z, Dunker AK. DisProt: the database of disordered proteins. Nucleic Acids Res. 2007;35(Database issue):D786–93. Epub 2006 Dec 1View ArticlePubMedGoogle Scholar
- Sierro N, Battey JND, Ouadi S, Bakaher N, Bovet L, Willig A, Goepfert S, Peitsch MC, Ivanov NV. The tobacco genome sequences and its comparison with those of tomato and potato. Nat Commun. 2014;5:3833.View ArticlePubMedPubMed CentralGoogle Scholar
- Sierro N, Battey JND, Ouadi S, Bovet L, Goepfert S, Bakaher N, Peitsch MN, Ivanov VI. Reference genomes and transcriptomes of Nicotiana sylvestris and Nicotiana tomentosiformis. Genome Biol. 2013;14(6):R60.View ArticlePubMedPubMed CentralGoogle Scholar
- Sigrist CJA, de Castro E, Cerutti L, Cuche BA, Hulo N, Bridge A, Bougueleret L, Xenarios I. New and continuing developments at PROSITE. Nucleic Acids Res. 2013;41:D344–7.View ArticlePubMedGoogle Scholar
- Smith AG, Eberle CA, Anderson NO, Clasen BM, Hegeman AD. The transmitting tissue of Nicotiana tabacum is not essential to pollen tube growth, and its ablation can reverse prezygotic interspecific barriers. Plant Reprod. 2013;26:339–50.View ArticlePubMedGoogle Scholar
- Stephenson AG, Travers SE, Mena-Ali JA, Winsor JA. Pollen performance before and during the autotrophic–heterotrophic transition of pollen tube growth. Phil Trans R Soc London. 2003;358:1009–18.View ArticleGoogle Scholar
- Swanson WJ, Vacquier VS. The rapid evolution of reproductive proteins. Nat Rev Genet. 2002;3(2):137–44.View ArticlePubMedGoogle Scholar
- Tamura K, Nei M. Estimation of the number of nucleotide substitutions in the control region of mitochondrial DNA in humans and chimpanzees. Mol Biol Evol. 1993;10(3):512–26.PubMedGoogle Scholar
- Tan L, Leykam JF, Kieliszewski MJ. Glycosylation motifs that that direct arabinogalactan addition to arabinogalactan-proteins. Plant Physiol. 2003;132(3):1362–9.View ArticlePubMedPubMed CentralGoogle Scholar
- Tompa P, Schad E, Tantos A, Kalmar L. Intrinsically disordered proteins: emerging interaction specialists. Curr Opin Struct Biol. 2015;35:49–59.View ArticlePubMedGoogle Scholar
- Uversky VN. Intrinsic disorder-based protein interactions and their modulators. Curr Pharm Des. 2013;19(23):4191–213.View ArticlePubMedGoogle Scholar
- Villalba M, Batanero E, López-Otín C, Sánchez LM, Monsalve RI. González de la Peña, MA, Lahoz C. Rodríguez R. The amino acid sequence of ole e I, the major allergen from olive tree (Olea europaea) pollen. Eur J Biochem. 1993;216:863–9.View ArticlePubMedGoogle Scholar
- Wolf JB, Lindell J, Backström N. Speciation genetics: current status and evolving approaches. Philos Trans R Soc B. 2010;365:1716–33.Google Scholar
- Wu H, Wang H, Cheung AY. A pollen tube growth stimulatory glycoprotein is deglycosylated by pollen tubes and displays a glycosylation gradient in the flower. Cell. 1995;83:395–403.View ArticleGoogle Scholar
- Wu HM, de Graaf B, Mariani C, Cheung, AY. Hydroxyproline-rich glycoproteins in plant reproductive tissues: structure, functions and regulation. Cell Mol Life Sci. 2001;58:1418–29.Google Scholar
- Wu HM, Wong E, Ogdahl J, Cheung AY. A pollen tube growth-promoting arabinogalactan protein from Nicotiana alata is similar to the tobacco TTS protein. Plant J. 2000;22(2):165–76.View ArticlePubMedGoogle Scholar
- Zhang XL, Ma HL, Qi HD, Zhao J. Roles of hydroxyproline-rich glycoproteins in the pollen tube and style cell growth of tobacco (Nicotiana tabacum L.). J. Plant Physiol. 2014;171(12):1036–45.View ArticleGoogle Scholar