- Research article
- Open Access
Evolutionary characterization and transcript profiling of β-tubulin genes in flax (Linum usitatissimum L.) during plant development
- Floriana Gavazzi†1,
- Gaia Pigna†1,
- Luca Braglia1,
- Silvia Gianì1,
- Diego Breviario1 and
- Laura Morello1Email authorView ORCID ID profile
© The Author(s). 2017
- Received: 4 April 2017
- Accepted: 29 November 2017
- Published: 8 December 2017
Microtubules, polymerized from alpha and beta-tubulin monomers, play a fundamental role in plant morphogenesis, determining the cell division plane, the direction of cell expansion and the deposition of cell wall material. During polarized pollen tube elongation, microtubules serve as tracks for vesicular transport and deposition of proteins/lipids at the tip membrane. Such functions are controlled by cortical microtubule arrays. Aim of this study was to first characterize the flax β-tubulin family by sequence and phylogenetic analysis and to investigate differential expression of β-tubulin genes possibly related to fibre elongation and to flower development.
We report the cloning and characterization of the complete flax β-tubulin gene family: exon-intron organization, duplicated gene comparison, phylogenetic analysis and expression pattern during stem and hypocotyl elongation and during flower development. Sequence analysis of the fourteen expressed β-tubulin genes revealed that the recent whole genome duplication of the flax genome was followed by massive retention of duplicated tubulin genes. Expression analysis showed that β-tubulin mRNA profiles gradually changed along with phloem fibre development in both the stem and hypocotyl. In flowers, changes in relative tubulin transcript levels took place at anthesis in anthers, but not in carpels.
Phylogenetic analysis supports the origin of extant plant β-tubulin genes from four ancestral genes pre-dating angiosperm separation. Expression analysis suggests that particular tubulin subpopulations are more suitable to sustain different microtubule functions such as cell elongation, cell wall thickening or pollen tube growth. Tubulin genes possibly related to different microtubule functions were identified as candidate for more detailed studies.
- Bast fibres
- Expression analysis
- Phylogenetic analysis
- Gene family
- Duplicated genes
- Pollen tube
Flax is an important crop, cultivated for both oil and fibre production. Bast fibres are phloem fibres characterized by an impressive cell elongation, up to several centimetres, and the deposition of gelatinous-type cell walls rich in cellulose, responsible for the particular tensile strength and flexibility of textile fibres. According to detailed studies performed by Gorshkova and co-workers [1, 2], different stem regions, corresponding to different fibre developmental stages, follow one another in fast growing flax plants. In the uppermost section, easily distinguished by smaller distance between leaves, most fibre cells elongate by coordinated growth, synchronously with the surrounding tissues and only the primary cell wall is present. Further downward, phloem fibres undergo intrusive growth by cell surface elongation. Finally, after elongation of bast fibres is completed, cell wall thickening takes place: the transition between cell elongation and cell wall thickening is marked by the snap point (SP), a region identified by the change in the mechanical properties of the stem, which becomes rigid and can’t anymore be easily broken by hand . Downstream the SP, the galactan-enriched matrix (Gn-layer) of the secondary cell wall is gradually modified into a mature cellulosic gelatinous-layer (G-layer), showing the typical crystalline morphology of thickened cell wall, characterized by microfibril axial orientation, xylan absence, and a low content of lignin . This change requires transcription of specific β-galactosidase genes (LuBGAL1 and LuBGAL2)  and β-galactosidase activity, involved in remodelling of the cell wall matrix. Toward the base of the stem all phloem vascular bundles have completed cell wall thickening. Phloem fibre differentiation proceeds similarly during hypocotyl elongation, where the onset of secondary cell wall thickening follows cell elongation and is marked by an increase in β-galactosidase activity and up-regulation of corresponding genes along time . For this reason, flax stems and hypocotyls have been used as a model to identify, through microarray studies, transcripts specifically expressed during phloem fibre differentiation [5, 6], and to study the expression of genes involved in cell wall development, such as glycosyl hydrolases (GH35) , cellulose synthase (CESA) [8, 9] and chitinase-like (CTL) genes .
Arrays of cortical microtubules (MT) control both the direction of cell elongation and cellulose deposition, fundamental processes related to fibre development. During bast fibre differentiation, as typical for diffusely elongating cells, MT orientation gradually changes from transversal, in the elongation stage, to longitudinal, re-orienting cellulose deposition at the onset of secondary cell wall thickening . Cortical microtubules have a role in guiding CESA complexes from the trans-Golgi network to the plasma membrane and to drive microfibril orientation through the direct interaction with CESA complex . Such interaction is mediated by cellulose-synthase interacting proteins (CSI) [11, 12]. It has also been reported that cellulose synthesis, defective in the AtCesa2 mutant of Arabidopsis, can in turn affect microtubule reorientation . MT also play a role in non-cellulosic polysaccharide assembly during cell wall biogenesis by recruiting vesicles carrying pectin and other cell wall matrix components to the cell cortex [14–17]. Finally, many signals triggering microtubule reorientation, such as auxin and light  were also known to differentially regulate transcription of α- and β-tubulin genes [19, 20]. This complex interplay between cell wall synthesis and cortical microtubule arrays might also involve specific tubulin isotypes. Tissue-specific expression of tubulin genes in plants has been widely reported  and up-regulation of particular tubulin isotypes has been related to specific cell functions, although direct evidence is limited. Overexpression of EgTub1has been reported to affect cellulose microfibril orientation in xylem fibres in Eucalyptus . In Populus, genes encoding for xylem-abundant tubulins are also up-regulated in tension wood, which is highly active in cellulose deposition . A subset of tubulin genes is preferentially expressed in elongating cotton fibres [24, 25]. One of these genes induces longitudinal growth when expressed in fission yeast . Another plant structure characterized by polarized fast growth is the pollen tube. MTs play a role in secretion and endocytosis, acting as tracks for vesicular transport and deposition of proteins/lipids at the tip membrane . Interestingly, overexpression of a Picea wilsonii pollen-specific alpha tubulin gene increases pollen tube growth rate in Arabidopsis . Although such studies are not direct evidence for specific functions of particular tubulin isotypes, we can hypothesize that MT assembled from different tubulin pools, with slightly different aminoacid sequences, may show different kinetic or binding properties.
The peculiarity of flax fibres coordinated development, with a spatial separation of elongation and thickening along the stem axis and a temporal scansion during hypocotyls development, offers an opportunity to identify candidate tubulin genes involved in specific cell growth stages. In this manuscript, the large β-tubulin gene family of flax was phylogenetically characterized and the expression pattern of each member was investigated during flower development and along with stem and hypocotyl growth. Tubulin genes possibly related to cell elongation or cell wall thickening were identified as candidate for more detailed studies.
Plant growth conditions and tissue collection
Flax seeds (Linum usitatissimum, cultivar Valoal) were purchased from Semfor SrL. Plants were grown in pots with soil, in a growth chamber, under a 14 h light/10 h dark cycle, at 26 °C day/20 °C night, with 50% humidity. Stem tissues were collected during the fast growth stage, one month after sowing, from plants with similar height, 22 cm on average. The following sections were collected, after leaflet removal: apex (A) the uppermost 0,5 cm section; top section (TS) from 0,5 cm to 1,5 cm below the apex, above the snap point; medium section (MS) from 3,5 cm to 5,5 cm below the apex, encompassing the SP (located approximately at 4,5 cm from the apex), and basal section (BS), 2 cm long, about 2 cm upstream to the cotyledonary leaves. Male and female reproductive organs were collected from mature plants the day before anthesis (still unopened flowers with unelongated anther filaments) and early in the anthesis day (open flowers with elongated anther filaments and mature pollen). For hypocotyl collection, sterilized seeds were sown in magenta boxes with solid MS medium, kept in the dark for the first three days and then exposed to a 16 h light/8 h dark cycle, at 25 °C. Hypocotyls of similar length were sectioned at 3, 5, 7 and 12 days after sowing (DAS). For all the tissues, three independent samples of 100 mg each were immediately frozen in liquid nitrogen and stored at −80 °C.
DNA extraction and capillary electrophoresis-TBP analysis (CE-TBP)
Genomic DNA was extracted from 100 mg of flax seedlings using the Qiacube® robotic workstation (Qiagen) with the DNeasy Plant Mini Kit (Qiagen). Extracted DNA was eluted in 150 μL Tris-EDTA buffer. DNA quality and quantity were assessed by NanoDrop-2000C spectrophotometer (Thermo Fisher Scientific). To generate the CE-TBP profiles, approximately 30 ng of genomic DNA was used as a template in TBP 1st and 2nd intron amplifications. The PCR reactions were performed in 30 μl according to . Negative PCR controls (no template) were always included in each analysis. 4 μL of each PCR product were loaded on a 2% agarose gel to define the proper dilution rate for CE analysis. Capillary Electrophoresis was performed on the 3500 Genetic Analyser (Thermo Fisher Scientific) as described . The data were collected using the Data Collection Software v. 2.0 (Thermo Fisher Scientific) and analysed by the GeneMapper Software v. 5.0 tool (Thermo Fisher Scientific).
RACE amplification and cloning of the β-tubulin coding sequences
Gene expression by qRT-PCR
Gene expression was studied in stem, hypocotyl and flower samples collected as described in the previous section. Total RNA was extracted from frozen samples as described for RACE. cDNA was obtained using the Revert Aid H Minus First Strand cDNA Synthesis Kit (Thermo Fisher Scientific), according to the manufacturer’s protocol, starting from 1 microgram of total RNA. Samples were finally treated with RNase H (Thermo Fisher Scientific) 0,1 U/ μL, at 37 °C for 20 min and were stored at −20 °C.
Isotype-specific primer pairs, suitable for qRT-PCR, were manually designed in the 3’UTR regions and evaluated with the help of the Primer Express software (Invitrogen). Primer specificity was evaluated by BLASTn search in the flax WGS database and confirmed in real-time PCR by re-association kinetics. LuBGAL1 (GenBank: HQ902252.1) primers were from . Specific primers for reference genes Elongation Factor 1-α (EF1α) and Glyceraldehyde-phosphate dehydrogenase (GAPDH) were from . All qPCR primer sequences are listed in Additional file 1.
A five-point standard curve of a 4-fold dilution series (1:4 to 1:1024) of a cDNA mix from various tissues, was used to calculate the PCR efficiency of each primer pair and the optimal cDNA dilution for quantitative analysis (Ct comprised between 20 and 28). For the expression analysis, 1,5 μL of a 1:16 dilution of each cDNA sample was used. qRT-PCR reactions were performed in 96 well plates, in a CFX Connect Real Time System (BioRad), in a final volume of 20 μl containing 1,5 μL of cDNA, 0,4 μM of each primer and 10 μL of 2× Sso Advanced TM Universal SYBR® Green Supermix (BioRad), under the following conditions: 95 °C for 3 min, 40 cycles of 15 s at 95 °C and 60 s at 62 °C. A 65–95 °C melting curve was performed to confirm the specificity of the products.
Each experiment was repeated twice on two biological replicates (RNA extraction and cDNA synthesis), with three technical repetitions each. Data analysis was performed using the BioRad CFX Manager Software. Threshold was manually set at the same level to compare results from different plates. Correction for specific amplification efficiency and inter-run variation were applied in Gene Studies. Ct values were normalized using the geometrical mean of EF1α and GAPDH as reference genes. Stability acceptance values (CV < 0.5; M < 1) for both reference genes were respected . The ΔΔCt method with Pfaffl correction for amplification efficiency [33, 34] was used to compare transcript fold increase, relative to zero (mean expression of each target) or to a control sample. Heat maps were generated using the same BioRad CFX Manager Software. Each data point represents the expression level of one target sequence in one sample, relative to the average of the reference genes. Data points are scaled by target and clustered by the degree of similarity of expression.
Cloned cDNA sequences were assembled and aligned using the Contig Express and Align X tools of the NTI Vector v.11 suite and manually refined. Phylogenetic analysis was performed using MEGA5 . Deduced amino acid sequences were aligned by MUSCLE with default settings. The phylogenetic analysis was performed choosing the “best substitution model” tool provided by the software. The chosen model (Jones Taylors-Thornton, JTT by Gamma distribution rate with Invariant sites, G + I) was applied to infer a phylogenetic tree by the Maximum Likelihood (ML) estimation method. The initial tree was inferred by Neighbour Joining analysis, following the Nearest-Neighbour-Interchange heuristic method. The statistical confidence of a particular group of accessions in the tree was evaluated by bootstrap test with 1000 replicates. The tree was rooted using Chlamydomonas reinhardtii β-tubulin as the outgroup. Accession numbers of all aminoacid sequences are provided in Additional file 2. Schematic drawing of the exon-intron structure of tubulin genes was obtained by alignment of the cDNA and genomic sequences in the Gene Structure Display Server website tool, v 2.0 (GSDS, http://gsds.cbi.pku.edu.cn, accessed 20 October 2016).
PCR amplification and cloning of the β-tubulin genes
To isolate simultaneously the coding region of all members of the β-tubulin gene family, we designed an amplification strategy based on degenerated primers effective on all plant β-tubulin genes. This is made possible because of the strong evolutionary conservation of tubulin coding sequences, due to the extraordinary structural constraints imposed to tubulin by its fundamental functions. Short nucleotide stretches, showing 70–80% overall identity, are useful to design slightly degenerated universal primers. Internal primer pairs Fex1-Rex1 and Fin2-Rin2, annealing to the ORF at the first and second intron borders respectively, have been abundantly tested because they are at the base of the genomic profiling method called tubulin-based polymorphism (TBP), exploiting on the different length of β-tubulin introns [29, 36]. A further, degenerated forward primer, named MREI, was designed at the very 5′ end of the β-tubulin open reading frame, where the first six aminoacid positions (M-R-E-I-L-H) are extremely conserved among all plant β-tubulins. RNA extracted from 5 days old seedlings was used as a template for cDNA synthesis with a T7-tagged primer. 5′ and 3’ cDNA regions were independently amplified with MREI/Rin2 and Fin2/T7 primer combination to generate two overlapping cDNA portions (Fig. 1a).
The degenerated MREI primer was first tested on flax genomic DNA. In combination with the reverse primer Rex1, it revealed a band pattern shifted in size by 300 bp with respect to a control amplification done with the forward primer Fex1, commonly used in the CE-TBP profiling protocol (Fig. 1b). As predicted, cDNA amplification of the 5’terminus of the ORF, spanning from the start codon to the middle of the third exon, originated a unique sharp band on agarose gel, since all plant β-tubulin cDNA have the same length in this region (Fig. 1c). Conversely, 3’RACE produced a smeared band due to the heterogeneous length of the third exon and the 3’UTR of β-tubulin genes, together with some shorter, non-specific amplicon (Fig. 1c). Amplified bands were eluted from gel and cloned into E.coli for sequencing. After nucleotide sequence assembly, tentative contigs were obtained for 14 different β-tubulin cDNA sequences, spanning from the beginning of the coding sequence to the poly-A tail. Contig accuracy was confirmed by cloning and sequencing of fourteen different cDNA fragments amplified with the MREI forward primer and Reverse Specific Primers (RSP), designed at the most downstream position of each 3’UTR (Additional file 1). cDNA sequences ranged in length from 1463 to 1665 nucleotides.
The flax β-tubulin gene family
List of Lusb-tub cDNA and corresponding genomic sequences
Sequence length (pb)
cDNA (GenBank ID)
ORF length (excluding stop)
TBP amplicon 1
TBP amplicon 2
Protein length (AA)
Putative genomic locus (GenBank ID)
Putative transcript (Phytozome ID)
Accordance of the putative β-tubulin gene family with the TBP profile
Expression profiles along fibre differentiation
β-tubulin gene expression was studied by RT-qPCR in stem sections at increasing distance from the apex and in hypocotyls at different days after sowing, taken as model systems for the progression of fibre differentiation. Stems of 30 days old plants (fast growth stage) were deprived of leaves and sectioned in order to distinguish three different phases of phloem fibre growth, according to what reported in published studies [2, 6]. 1) The top section (TS), easily distinguished by short internodes, characterized by fast elongating phloem cells with primary cell wall and transversally oriented cortical MT; 2) the medium section (MS), around the SP, characterized by the gradual cessation of fibre cell elongation and the onset of secondary cell wall thickening, with a prevalence of helical arrays of cortical MT, due to their gradual reorientation; 3) the basal, or bottom section (BS), close to the base of the stem, containing fully differentiated fibres, with thickened, gelatinous type cell walls and longitudinally oriented MT. The apical section (A), which includes the apical meristem with actively dividing and early differentiating cells, was analysed as a control. In elongating hypocotyls, cell division has almost ceased during the pre-germinative stage and fibre differentiation proceeds similarly to the stem, with cell elongation preceding in time secondary cell wall deposition and thickening . Hypocotyls were collected at four different stages, DAS 3, 5, 7 and 12, toward the end of elongation. This study provided an overview on the modulation of β-tubulin gene expression during plant growth, useful to identify some interesting candidate for further investigations on fibre differentiation.
Due to the very high degree of sequence similarity between ohnologs, gene-specific primers were designed in the 3’untranslated regions of each tubulin cDNA (Additional file 3). Expression of the β-Galactosidase 1 gene (LuBGal1), known to be up-regulated during cell wall galactan formation in fibre cells [4, 6, 45], was used as a marker for the fibre developmental stage in both the stem and hypocotyl.
Comparison of the same data, reported as fold changes of each target gene with respect to the apex in the stem and to DAS3 in the hypocotyl, highlighted similarity in the modulation of the β-tubulin gene expression pattern, despite different absolute levels (Fig. 6c and d). Genes that were more stably up-regulated during fibre growth in the stem, LusTub1, 2 and 4, were similarly up-regulated in the hypocotyl. LusTub5, 3 and 6, transiently up-regulated during the elongation stage, were progressively down-regulated as cell wall thickening progressed. LusTub7 transcript level decreased in both organs with respect to the apex or the 3-DAS stage, in parallel with fibre elongation and differentiation. This observation is in agreement with a strictly controlled expression of β-tubulin genes during stem and hypocotyls elongation, possibly related to fibre differentiation.
β-tubulin gene expression in reproductive organs
Expression of duplicated β-tubulin genes
Percent nucleotide identity between ohnologs
0.5 kb UP
1 kb UP
The most evident feature of the flax β-tubulin gene family, based on cDNA sequence alignment, is the presence of five pairs of duplicated genes (LusTub1, 2, 3, 4 and 6), plus one triplet, likely due to a further single-gene duplication (LusTub7), and one singleton (LusTub5). This configuration is in agreement with data from the flax genome assembly, inferring the occurrence of a WGD in the flax lineage, dated around 5–9 MYA . This evidence is also supported by karyotype studies that assigned larger chromosome sets (2n = 30) to both L. usitatissimum and L. bienne, the likely wild ancestor of cultivated flax, with respect to other close relatives such as L. grandiflorum (2n = 16) . This supports the occurrence of chromosome repatterning and diploydisation following tetraploydisation. The high level of retention of duplicated β-tubulin genes after the WGD event, with just one probable gene loss, the ohnolog of LusTub5, is in accordance with similar occurrence in other plants. Tubulin gene family expansion following genome duplication have been reported for the genome of cotton (G. hirsutum), a quite recent allotetraployd (1MYA) accounting for at least 19 β-tubulin genes , as well as for the Populus and Salix genera, both having 10 pairs of highly similar β-tubulin genes [23, 55], reminiscent of the ancient salicoid WGD event, estimated 60–65 MYA . In the flax genome, gene expansion was reported for other gene families, such as β-galactosidase, cellulose synthase and chitinase-like [6, 8, 57], but duplicate gene retention was lower and characterized by selective gene losses and differential expansion of specific subfamilies. Our observation is supported by recent models of duplicated gene evolution suggesting that rates of duplicate gene retention vary among functional groups. Duplicate genes retained after WGD have the tendency to persist and proliferate following subsequent duplication events, while singleton showed an opposite tendency to be repeatedly restored to the original status . The tendency to high retention rates was found to be more pronounced among genes encoding long complex proteins and proteins with high number of interactions. This is consistent with the gene balance hypothesis, which predicts that the fate of duplicate genes largely depends on maintaining a stoichiometric balance among members of macromolecular complexes [58, 59]. Such explanation is particularly suited for tubulins, whose unbalanced synthesis may induce lethal phenotype . Duplicated genes are not only retained but were found to accumulate mutations at a slower rate than singleton genes in rice subspecies and A. thaliana ecotypes . In the genome of Populus, relative rates of non-synonymous substitution were actually lower for genes with paralogs from the salicoid and eurosid whole-genome duplication events than for genes with no paralogs .
Tubulins in particular are very ancient, highly conserved proteins, since even single AA substitutions in a single isotype have high chance to confer lethal or dominant negative phenotypes which may severely alter growth morphology in plants  or lead to serious neurodevelopmental diseases, defined as tubulinopathies, in humans . Therefore, strong purifying selection is likely to act toward the fixation of tubulin duplicates. Passive forces rather than adaptive mechanism are likely to act to keep tubulin gene family large. Nevertheless, tubulin genes can slowly evolve and acquire slightly different properties. In flax, aminoacid substitutions between tubulin ohnologs were minimal, up to 4–5, and clustered in the CTT. Although protruding from the microtubule core, this highly polymorphic domain can provide distinctive features to the proteins. It is in fact the binding site for a large set of microtubule associated proteins (MAPs; ) and the target for multiple post-translational modifications such as tyrosination, phosphorylation or glutamylation, as largely documented by immunological data (reviewed in ), although mass spectrometry-based analysis has recently questioned the abundance of such modifications in some plants or plant tissues [64, 65]. Although not affecting dimer formation, CTT is also a key determinant for microtubule assembly, influencing the kinetics of polymerization/depolymerisation . Therefore, mutations in this domain can be well tolerated and may gradually accumulate providing slightly different properties to tubulin isotypes and to their corresponding microtubules.
Not only duplicated flax β-tubulin genes show minimal sequence divergence, but they also retain the same expression pattern, with some little bias toward one other or the other ohnolog observed in few cases, in the tissues we investigated. LusTub4a and LusTub4b, showing the higher, although modest in absolute values, difference at transcriptional level encode for identical proteins, but carry divergent regulatory regions. This suggests that differential expression may precede protein modification in tubulin gene evolution. In the aspen genome, in which the last whole genome duplication event dates back to 60–65 MYA, the ten pairs of duplicated β-tubulin genes still retain similar expression patterns, with few exceptions (i.ePtrTub 9 and PtrTub10), despite a greater divergence of their aminoacid sequences . A more detailed study of duplicated tubulin genes in different plants could reveal new interesting insights on the evolution mechanism of this fundamental gene family. The functional significance of the presence of numerous slightly different tubulin isoforms in plants is still poorly understood, despite the many attempts of assigning them functions.
Model systems, characterized by synchronously developing cells, may help to unravel this matter. Mesophyll cells in barley leaf sections have been used to study the differential and gradual modulation of steady state mRNA level of the five alpha tubulin genes that occur in parallel to the formation of differently oriented cortical microtubule arrays . The cotton fibre, a trichome elongating from a single seed epidermal cell, is another useful model used for the identification of tubulin genes related to cell elongation and cellulose deposition [24–26]. One of these studies also led to the isolation of a specific cotton β-tubulin gene able to induce abnormal elongation in fission yeast .
Flax stems and hypocotyls have been successfully used as a model to identify transcriptional changes that might correlate to the different developing stage of xylem and phloem fibres [6, 9, 45]. Since secondary cell wall deposition in xylem fibres begins early in hypocotyl development and in the top section in the stem, it was expected that a substantial portion of the differentially transcribed genes detected would have mainly reflected phloem-specific developmental events. Microarray analysis showed that actin, tubulin and other cytoskeletal proteins are enhanced in developing fibres, but no information was reported about differential expression.
In this work, we observed that the relative level of β-tubulin transcripts was modulated during stem and hypocotyls elongation. Using the apex as a control, we observed a more or less pronounced up-regulation of most genes in the TS, which is characterised by an extensive coordinated growth of fibre cells with other cell types. This process is associated with a transverse orientation of cortical MT that assist the correct positioning of CESA complexes in the membrane. Some of the tubulin genes, LusTub3, 5 and 6, were clearly down-regulated in the next stem sections, MS and BS, characterised by cortical MT re-orientation, concomitant with the arrest of elongation and gradual cell wall thickening in phloem fibres, while only cambium derivatives are still elongating. These two sections showed similar tubulin expression profiles, in accordance to a more general picture also envisaged by microarray analysis . Conversely, up-regulation of LusTub1, 2 and 4 in TS was stably maintained across the stem. LusTub7 paralogs showed an opposite behaviour, being mostly expressed in the apex and down-regulated in all the other sections. Interestingly, aspen PtrTub1 and PtrTub2, the most expressed β-tubulin genes in the apex, belong to the same class II as LusTub7. Comparison of large expression datasets in different plants could provide hints for some functional difference among tubulin classes. In the hypocotyl, a generalized increase in tubulin gene expression was observed at DAS5 with respect to DAS3, which paralleled the gradual transition from cell fibre elongation to cell wall thickening, as documented by the LuBGal1 expression. Similar to what observed in the stem, LusTub3, 5 and 6up-regulation was transient and rapidly declined with fibre development, while the increase of LusTub4 mRNA was more stable, and LusTub1 and 2 transcript level dropped at DAS12. LusTub7 paralogs were gradually downregulated with respect to early stages. Similarity in the relative expression pattern between stem and hypocotyl is in accordance with their correlation to fibre development. This could be further confirmed by expression studies on isolated fibres. How the different expression pattern observed influences microtubule array dynamics and cell function remains to be elucidated. Although post-transcriptional mechanisms are known to control tubulin levels, it is reasonable to infer that differential gene expression leads to the formation of different intracellular pools of tubulins, which may be further modified by post-translational modifications. It is then likely that microtubule arrays, randomly assembled from such different pools, have different kinetic properties and/or interact with different subsets of associated MAPs, making them better suited for specific functions (i.e, CESA positioning rather than vesicular transport for secondary cell wall deposition).
This can also apply to what observed in flax anthers. Complete pollen maturation at anthesis brings about changes in the steady-state transcript level of β-tubulin, with significant up-regulation of LusTub 3 and 4, paralleled by down-regulation of LusTub2 and 5, whose transcript become barely detectable. The newly synthesized tubulin pool might be more suited to build MTs more specialized for fine delivery of vesicles carrying proteins/lipids at the tip membrane to sustain polarized pollen tube growth.
We first described the flax β-tubulin gene family as composed by fourteen expressed members, reminiscent of the last WGD event during flax speciation. ML evolutionary tree including the complete set of β-tubulin genes of eight angiosperm species suggests that the multiple gene family originated by subsequent duplications of four ancestral genes. Gene expression analysis in stem, hypocotyl and flower suggest that different tubulin pools may originate MT with possible specialized functions.
Functional redundancy of tubulins has since now made difficult to dissect possible functions by simple knockout or overexpression experiments. Complete characterization of the gene family and an overview of the expression pattern during flax fibre elongation are the pre-requisite for further studies aimed to elucidate functional differences among similar tubulins. Further confirmation of the observed expression pattern can be obtained in isolated fibres. Genome editing studies with CRISPR-Cas9, may reveal new insights on such tubulin specificities. A future goal could be the modulation of the expression of tubulin genes specifically involved in cell elongation, in order to increase fibre length in flax. Tubulin genes that are up- or down-regulated during fibre elongation, such as LusTub3, 5 and 6, are good candidate for further studies.
We wish to thank Andrea Frascotti for collecting flax flowers pre and post-anthesis.
The project was financially supported by Regione Lombardia/CNR–Institutional Agreement n. 14,840/RCC 2011–2013. Research Project: VeLiCa - From ancient crops materials and products for the future.
Availability of data and materials
Flax β-tubulin cDNA sequences are available from GenBank under the accessions KM196480- KM196493.Other datasets used and/or analysed during the current study are available from the corresponding author on reasonable request.
LM Conceived the study and wrote the manuscript; GP, FG and SG performed the experiments and discussed results; LB performed phylogenetic analysis; FG and LM analyzed the data; DB supervised the work and critically revised the manuscript. All Authors read and accepted the final 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.
- Ageeva M, Petrovska B, Kieft H, Sal'nikov V, Snegireva A, van Dam J, van Veenendaal W, Emons A, Gorshkova T, van Lammeren A. Intrusive growth of flax phloem fibers is of intercalary type. Planta. 2005;222(4):565–74.View ArticlePubMedGoogle Scholar
- Gorshkova T, Sal'nikova V, Chemikosova S, Ageeva M, Pavlencheva N, van Dam J. The snap point: a transition point in Linum usitatissimum bast fiber development. Ind Crop Prod. 2003;18(3):213–21.View ArticleGoogle Scholar
- Gorshkova T, Gurjanov O, Mikshina P, Ibragimova N, Mokshina N, Salnikov V, Ageeva M, Amenitskii S, Chernova T, Chemikosova S. Specific type of secondary cell wall formed by plant fibers. Russ J Plant Physiol. 2010;57(3):328–41.View ArticleGoogle Scholar
- Roach MJ, Mokshina NY, Badhan A, Snegireva AV, Hobson N, Deyholos MK, Gorshkova TA. Development of cellulosic secondary walls in flax fibers requires beta-galactosidase. Plant Physiol. 2011;156(3):1351–63.View ArticlePubMedPubMed CentralGoogle Scholar
- Roach M, Deyholos M. Microarray analysis of developing flax hypocotyls identifies novel transcripts correlated with specific stages of phloem fibre differentiation. Ann Bot. 2008;102(3):317–30.View ArticlePubMedPubMed CentralGoogle Scholar
- Roach M, Deyholos M. Microarray analysis of flax (Linum usitatissimum L.) stems identifies transcripts enriched in fibre-bearing phloem tissues. Mol Gen Genomics. 2007;278(2):149–65.View ArticleGoogle Scholar
- Hobson N, Deyholos M. Genomic and expression analysis of the flax (Linum usitatissimum) family of glycosyl hydrolase 35 genes. BMC Genomics. 2013;14:344.Google Scholar
- Mokshina N, Gorshkova T, Deyholos M. Chitinase-like (CTL) and cellulose synthase (CESA) gene expression in gelatinous-type cellulosic walls of flax (Linum usitatissimum L.) bast fibers. PLoS One. 2014;9(6)Google Scholar
- Mokshina N, Gorshkov O, Ibragimova N, Chernova T, Gorshkova T. Cellulosic fibres of flax recruit both primary and secondary cell wall cellulose synthases during deposition of thick tertiary cell walls and in the course of graviresponse. Funct Plant Biol. doi:10.1071/FP17105
- Paredez A, Somerville C, Ehrhardt D. Visualization of cellulose synthase demonstrates functional association with microtubules. Science. 2006;312(5779):1491–5.View ArticlePubMedGoogle Scholar
- Bringmann M, Li E, Sampathkumar A, Kocabek T, Hauser MT, Persson S. POM-POM2/cellulose synthase interacting1 is essential for the functional association of cellulose synthase and microtubules in Arabidopsis. Plant Cell. 2012;24(1):163–77.View ArticlePubMedPubMed CentralGoogle Scholar
- Li S, Lei L, Somerville C, Gu Y. Cellulose synthase interactive protein 1 (CSI1) links microtubules and cellulose synthase complexes. Proc Natl Acad Sci U S A. 2012;109(1):185–90.View ArticlePubMedGoogle Scholar
- Chu Z, Chen H, Zhang Y, Zhang Z, Zheng N, Yin B, Yan H, Zhu L, Zhao X, Yuan M, et al. Knockout of the AtCESA2 gene affects microtubule orientation and causes abnormal cell expansion in Arabidopsis. Plant Physiol. 2007;143(1):213–24.View ArticlePubMedPubMed CentralGoogle Scholar
- Domozych D, Sorensen I, Popper Z, Ochs J, Andreas A, Fangel J, Pielach A, Sacks C, Brechka H, Ruisi-Besares P, et al. Pectin metabolism and assembly in the cell wall of the charophyte green algaPenium margaritaceum. Plant Physiol. 2014;165(1):105–18.View ArticlePubMedPubMed CentralGoogle Scholar
- Kong Z, Ioki M, Braybrook S, Li S, Ye Z, Lee Y, Hotta T, Chang A, Tian J, Wang G, et al. Kinesin-4 functions in vesicular transport on cortical microtubules and regulates cell wall mechanics during cell elongation in plants. Mol Plant. 2015;8(7):1011–23.View ArticlePubMedGoogle Scholar
- Yoneda A, Ito T, Higaki T, Kutsuna N, Saito T, Ishimizu T, Osada H, Hasezawa S, Matsui M, Demura T. Cobtorin target analysis reveals that pectin functions in the deposition of cellulose microfibrils in parallel with cortical microtubules. Plant J. 2010;64(4):657–67.View ArticlePubMedGoogle Scholar
- Zhu C, Ganguly A, Baskin T, McClosky D, Anderson C, Foster C, Meunier K, Okamoto R, Berg H, Dixit R. The fragile fiber1 kinesin contributes to cortical microtubule-mediated trafficking of cell wall components. Plant Physiol. 2015;167(3):780–92.View ArticlePubMedPubMed CentralGoogle Scholar
- Chen X, Grandont L, Li H, Hauschild R, Paque S, Abuzeineh A, Ralcusova H, Benkova E, Perrot-Rechenmann C, Friml J. Inhibition of cell expansion by rapid ABP1-mediated auxin effect on microtubules. Nature. 2014;516(7529):90–U206.PubMedPubMed CentralGoogle Scholar
- Bustos M, Guiltinan M, Cyr R, Ahdoot D, Fosket D. Light regulation of beta-tubulin gene-expression during internode development in soybean (Glycine max [l] Merr). Plant Physiol. 1989;91(3):1157–61.View ArticlePubMedPubMed CentralGoogle Scholar
- Qin X, Giani S, Breviario D. Molecular cloning of three rice alpha-tubulin isotypes: differential expression in tissues and during flower development. Biochim Et Biophys Acta-Gene Struct and Expr. 1997;1354(1):19–23.View ArticleGoogle Scholar
- Radchuk V. The transcriptome of the tubulin gene family in plants. In: The plant cytoskeleton: a key tool for agro-biotechnology. New York, NY: Springer Dordrecht; 2008.Google Scholar
- Spokevicius A, Southerton S, MacMillan C, Qiu D, Gan S, Tibbits J, Moran G, Bossinger G. Beta-tubulin affects cellulose microfibril orientation in plant secondary fibre cell walls. Plant J. 2007;51(4):717–26.View ArticlePubMedGoogle Scholar
- Oakley R, Wang Y, Ramakrishna W, Harding S, Tsai C. Differential expansion and expression of alpha- and beta-tubulin gene families in Populus. Plant Physiol. 2007;145(3):961–73.View ArticlePubMedPubMed CentralGoogle Scholar
- Li X, Cai L, Cheng N, Liu J. Molecular characterization of the cotton GhTUB1 gene that is preferentially expressed in fiber. Plant Physiol. 2002;130(2):666–74.View ArticlePubMedPubMed CentralGoogle Scholar
- He X, Qin Y, Xu Y, Hu C, Zhu Y. Molecular cloning, expression profiling, and yeast complementation of 19 beta-tubulin cDNAs from developing cotton ovules. J Exp Bot. 2008;59(10):2687–95.View ArticlePubMedPubMed CentralGoogle Scholar
- Ji S, Lu Y, Li J, Wei G, Liang X, Zhu Y. A beta-tubulin-like cDNA expressed specifically in elongating cotton fibers induces longitudinal growth of fission yeast. Biochem Biophys Res Commun. 2002;296(5):1245–50.View ArticlePubMedGoogle Scholar
- Onelli E, Idilli AI, Moscatelli A. Emerging roles for microtubules in angiosperm pollen tube growth highlight new research cues. Front Plant Sci. 2015;6:51.View ArticlePubMedPubMed CentralGoogle Scholar
- Yu Y, Li Y, Li L, Lin J, Zheng C, Zhang L. Overexpression of PwTUA1, a pollen-specific tubulin gene, increases pollen tube elongation by altering the distribution of alpha-tubulin and promoting vesicle transport. J Exp Bot. 2009;60(9):2737–49.View ArticlePubMedPubMed CentralGoogle Scholar
- Breviario D, Baird W, Sangoi S, Hilu K, Blumetti P, Giani S. High polymorphism and resolution in targeted fingerprinting with combined ss-tubulin introns. Mol Breed. 2007;20(3):249–59.View ArticleGoogle Scholar
- Gavazzi F, Braglia L, Mastromauro F, Giani S, Morello L, Breviario D. The tubulin-based-polymorphism method provides a simple and effective alternative to the genomic profiling of grape. PLoS One. 2016;11(9)Google Scholar
- Huis R, Hawkins S, Neutelings G. Selection of reference genes for quantitative gene expression normalization in flax (Linum usitatissimum L). BMC Plant Biol. 2010;10Google Scholar
- Hellemans J, Mortier G, De Paepe A, Speleman F, Vandesompele J. qBase relative quantification framework and software for management and automated analysis of real-time quantitative PCR data. Genome Biol. 2007;8(2)Google Scholar
- Livak K, Schmittgen T. Analysis of relative gene expression data using real-time quantitative PCR and the 2(-ΔΔC T ) method. Methods. 2001;25(4):402–8.View ArticlePubMedGoogle Scholar
- Pfaffl M. A new mathematical model for relative quantification in real-time RT-PCR. Nucleic Acids Res. 2001;29(9)Google Scholar
- Tamura K, Peterson D, Peterson N, Stecher G, Nei M, Kumar S. MEGA5: molecular evolutionary genetics analysis using maximum likelihood, evolutionary distance, and maximum parsimony methods. Mol Biol Evol. 2011;28(10):2731–9.View ArticlePubMedPubMed CentralGoogle Scholar
- Bardini M, Lee D, Donini P, Mariani A, Giani S, Toschi M, Lowe C, Breviario D. Tubulin-based polymorphism (TBP): a new tool, based on functionally relevant sequences, to assess genetic diversity in plant species. Genome. 2004;47(2):281–91.View ArticlePubMedGoogle Scholar
- Wang Z, Hobson N, Galindo L, Zhu S, Shi D, McDill J, Yang L, Hawkins S, Neutelings G, Datla R, et al. The genome of flax (Linum usitatissimum) assembled de novo from short shotgun sequence reads. Plant J. 2012;72(3):461–73.View ArticlePubMedGoogle Scholar
- Wolfe K. Robustness - it's not where you think it is. Nat Genet. 2000;25(1):3–4.View ArticlePubMedGoogle Scholar
- Breviario D, Giani S, Morello L. Multiple tubulins: evolutionary aspects and biological implications. Plant J. 2013;75(2):202–18.View ArticlePubMedGoogle Scholar
- Gavazzi F, Casazza A, Depedro C, Mastromauro F, Breviario D. Technical improvement of the TBP (tubulin-based polymorphism) method for plant species detection, based on capillary electrophoresis. Electrophoresis. 2012;33(18):2840–51.View ArticlePubMedGoogle Scholar
- Hashimoto T. Dissecting the cellular functions of plant microtubules using mutant tubulins. Cytoskeleton. 2013;70(4):191–200.View ArticlePubMedGoogle Scholar
- Lowe J, Li H, Downing K, Nogales E. Refined structure of alpha beta-tubulin at 3.5 a resolution. J Mol Biol. 2001;313(5):1045–57.View ArticlePubMedGoogle Scholar
- Selvaa Kumar C, Gadewal N, Mohammed SM. Seminal role of deletion of amino acid residues in H1-S2 and S-loop regions in eukaryotic β-tubulin investigated from docking and dynamics perspective. J Theor Biol. 2015;378:79–88.View ArticlePubMedGoogle Scholar
- Gorshkova T, Morvan C. Secondary cell-wall assembly in flax phloem fibers: role of galactans. Planta. 2006;223(2):149–58.View ArticlePubMedGoogle Scholar
- Hotte N, Deyholos M. A flax fibre proteome: identification of proteins enriched in bast fibres. BMC Plant Biol. 2008;8Google Scholar
- Carpenter J, Ploense S, Snustad D, Silflow C. Preferential expression of an alpha-tubulin gene of arabidopsis in pollen. Plant Cell. 1992;4(5):557–71.PubMedPubMed CentralGoogle Scholar
- Kim Y, An G. Pollen-specific expression of the Arabidopsis thaliana Alpha 1-tubulin promoter assayed by beta-glucuronidase, chloramphenicol acetyltransferase and diphtheria toxin reporter genes. Transgenic Res. 1992;1(4):188–94.View ArticlePubMedGoogle Scholar
- Rogers H, Greenland A, Hussey P. Fourmembers of the maize beta-tubulin gene family are expressed in the male gametophyte. Plant J. 1993;4(5):875–82.View ArticlePubMedGoogle Scholar
- Villemur R, Haas N, Joyce C, Snustad D, Silflow C. Characterization of four new beta-tubulin genes and their expression during male flower development in maize (Zea-Mays L). Plant Mol Biol. 1994;24(2):295–315.View ArticlePubMedGoogle Scholar
- Schewe L, Sawhney V, Davis A. Ontogeny of floral organs in flax (Linum usitatissimum; linaceae). Am J Bot. 2011;98(7):1077–85.View ArticlePubMedGoogle Scholar
- Lynch M, Conery J. The evolutionary fate and consequences of duplicate genes. Science. 2000;290(5494):1151–5.View ArticlePubMedGoogle Scholar
- Ohno S. Evolution by gene duplication. New York: Springer-Verlag; 1970.View ArticleGoogle Scholar
- Chapman B, Bowers J, Feltus F, Paterson A. Buffering of crucial functions by paleologous duplicated genes may contribute cyclicality to angiosperm genome duplication. Proc Natl Acad Sci U S A. 2006;103(8):2730–5.View ArticlePubMedPubMed CentralGoogle Scholar
- Muravenko OV, Bol'sheva NL, OIu I, Nosova IV, Rachinskaia OA, Samatadze TE, Zelenin AV. Karyogenomics of species of the genus Linum L. Genetika. 2010;46(10):1339–42.PubMedGoogle Scholar
- Rao G, Zeng Y, He C, Zhang J. Characterization and putative post-translational regulation of alpha- and beta-tubulin gene families in Salix arbutifolia. Sci Rep. 2016;6Google Scholar
- Tuskan G, DiFazio S, Jansson S, Bohlmann J, Grigoriev I, Hellsten U, Putnam N, Ralph S, Rombauts S, Salamov A, et al. The genome of black cottonwood, Populus Trichocarpa (Torr. & gray). Science. 2006;313(5793):1596–604.View ArticlePubMedGoogle Scholar
- Pydiura NA, GYa B, Galinouski DV, Yemets AI, PirkoYaV PTA, Anisimova NV, Khotyleva LV, Kilchevski AV, YaB B. Bioinformatic search for cellulose synthase genes in flax (Linum Usitatissimum) and their phylogenetic analysis. Cytol Genet. 2015;49(5):279–87. https://doi.org/10.3103/S0095452715050084.View ArticleGoogle Scholar
- Freeling M, Thomas B. Gene-balanced duplications, like tetraploidy, provide predictable drive to increase morphological complexity. Genome Res. 2006;16(7):805–14.View ArticlePubMedGoogle Scholar
- Rodgers-Melnick E, Mane S, Dharmawardhana P, Slavov G, Crasta O, Strauss S, Brunner A, DiFazio S. Contrasting patterns of evolution following whole genome versus tandem duplication events in Populus. Genome Res. 2012;22(1):95–105.View ArticlePubMedPubMed CentralGoogle Scholar
- Anthony R, Hussey P. Suppression of endogenous alpha and beta tubulin synthesis in transgenic maize calli overexpressing alpha and beta tubulins. Plant J. 1998;16(3):297–304.View ArticlePubMedGoogle Scholar
- Romaniello R, Arrigoni F, Bassi M, Borgatti R. Mutations in alpha- and beta-tubulin encoding genes: implications in brain malformations. Brain Dev. 2015;37(3):273–80.View ArticlePubMedGoogle Scholar
- Gardiner J. The evolution and diversification of plant microtubule-associated proteins. Plant J. 2013;75(2):219–29.View ArticlePubMedGoogle Scholar
- Parrotta L, Cresti M, Cai G. Accumulation and post-translational modifications of plant tubulins. Plant Biol. 2014;16(3):521–7.View ArticlePubMedGoogle Scholar
- Hotta T, Fujita S, Uchimura S, Noguchi M, Demura T, Muto E, Hashimoto T. Affinity purification and characterization of functional tubulin from cell suspension cultures of arabidopsis and tobacco. Plant Physiol. 2016;170(3):1189–205.PubMedPubMed CentralGoogle Scholar
- Hu H, Gu X, Xue L, Swamy P, Harding S, Tsai C. Tubulin C-Terminal post-translational modifications do not occur in wood forming tissue of Populus. Front Plant Sci. 2016;7Google Scholar
- Freedman H, Luchko T, Luduena R, Tuszynski J. Molecular dynamics modeling of tubulin C-terminal tail interactions with the microtubule surface. Proteins-Structure Function and Bioinformatics. 2011;79(10):2968–82.View ArticleGoogle Scholar
- Schroder J, Stenger H, Wernicke W. Alpha-tubulin genes are differentially expressed during leaf cell development in barley (Hordeum vulgare L). Plant Mol Biol. 2001;45(6):723–30.View ArticlePubMedGoogle Scholar