BackgroundThere is a rapidly growing awareness that plant peptide signalling molecules are numerous and varied and they are known to play fundamental roles in angiosperm plant growth and development. Two closely related peptide signalling molecule families are the CLAVATA3-EMBRYO-SURROUNDING REGION (CLE) and CLE-LIKE (CLEL) genes, which encode precursors of secreted peptide ligands that have roles in meristem maintenance and root gravitropism. Progress in peptide signalling molecule research in gymnosperms has lagged behind that of angiosperms. We therefore sought to identify CLE and CLEL genes in gymnosperms and conduct a comparative analysis of these gene families with angiosperms.ResultsWe undertook a meta-analysis of the GenBank/EMBL/DDBJ gymnosperm EST database and the Picea abies and P. glauca genomes and identified 93 putative CLE genes and 11 CLEL genes among eight Pinophyta species, in the genera Cryptomeria, Pinus and Picea. The predicted conifer CLE and CLEL protein sequences had close phylogenetic relationships with their homologues in Arabidopsis. Notably, perfect conservation of the active CLE dodecapeptide in presumed orthologues of the Arabidopsis CLE41/44-TRACHEARY ELEMENT DIFFERENTIATION (TDIF) protein, an inhibitor of tracheary element (xylem) differentiation, was seen in all eight conifer species. We cloned the Pinus radiata CLE41/44-TDIF orthologues. These genes were preferentially expressed in phloem in planta as expected, but unexpectedly, also in differentiating tracheary element (TE) cultures. Surprisingly, transcript abundances of these TE differentiation-inhibitors sharply increased during early TE differentiation, suggesting that some cells differentiate into phloem cells in addition to TEs in these cultures. Applied CLE13 and CLE41/44 peptides inhibited root elongation in Pinus radiata seedlings. We show evidence that two CLEL genes are alternatively spliced via 3′-terminal acceptor exons encoding separate CLEL peptides.ConclusionsThe CLE and CLEL genes are found in conifers and they exhibit at least as much sequence diversity in these species as they do in other plant species. Only one CLE peptide sequence has been 100% conserved between gymnosperms and angiosperms over 300 million years of evolutionary history, the CLE41/44-TDIF peptide and its likely conifer orthologues. The preferential expression of these vascular development-regulating genes in phloem in conifers, as they are in dicot species, suggests close parallels in the regulation of secondary growth and wood formation in gymnosperm and dicot plants. Based on our bioinformatic analysis, we predict a novel mechanism of regulation of the expression of several conifer CLEL peptides, via alternative splicing resulting in the selection of alternative C-terminal exons encoding separate CLEL peptides.

There is a wide diversity of bioinformatic tools available for the assembly of next generation sequence and subsequence variant calling to identify genetic markers at scale. Integration of genomics tools such as genomic selection, association studies, pedigree analysis and analysis of genetic diversity, into operational breeding is a goal for New Zealand’s most widely planted exotic tree species, Pinus radiata. In the absence of full reference genomes for large megagenomes such as in conifers, RNA sequencing in a range of genotypes and tissue types, offers a rich source of genetic markers for downstream application. We compared nine different assembler and variant calling software combinations in a single transcriptomic library and found that Single Nucleotide Polymorphism (SNPs) discovery could vary by as much as an order of magnitude (8,061 SNPs up to 86,815 SNPs). The assembler with the best realignment of the packages trialled, Trinity, in combination with several variant callers was then applied to a much larger multi-genotype, multi-tissue transcriptome and identified 683,135 in silico SNPs across a predicted 449,951 exons when mapped to the Pinus taeda ver 1.01e reference.

Standard protocols for extracting genomic DNA from Pinus radiata D. Don needles, such as CTAB-based methods, can yield large quantities of DNA. However, final DNA purity can be an issue due to carry over of contaminants that can impede accurate high throughput genotyping. This study evaluated eight DNA extraction and purification protocols to determine which method provided the greatest improvement in call rates and accuracy when using the Sequenom iPLEX W Gold MassARRAY W genotyping technology. Of the methods tested, genomic DNA extracted using the Machery-Nagel NucleoSpin W -96 Plant II kit performed the best overall, and was more efficiently and accurately genotyped than genomic DNA extracted using the standard CTAB method. This study also demonstrated that the quality and assay performance of CTAB-extracted genomic DNA is greatly improved by further purification with the Qiagen W QIAquick 96 PCR Purification kit. Using these improvements, the Sequenom iPLEX W Gold MassARRAY W genotyping technology is now a viable option for genotyping plant genomes such as Pinus radiata.

Associations of single nucleotide polymorphisms (SNP) with the form traits, branch cluster frequency and stem straightness, were studied in three radiata pine (Pinus radiata D. Don) breeding populations. For branch cluster frequency, high genotype by environment (G×E) interactions were found between two sites of the POP1 trial series (Tarawera and Glenledi) and between two sites of the POP2 trial series (Tarawera and Woodhill). For stem straightness, high G×E interactions were found between two sites (Tarawera and Woodhill) in both the POP2 and POP3 trial series. Thirtytwo and 26 SNPs were associated with branch cluster frequency and stem straightness respectively, across a number of sites and trial series in a genetic model assuming heterogeneous additive genetic variances across sites. These SNPs had significant associations with branch cluster frequency or stem straightness either in one site only, across sites within a trial series or across different trial series. Seven SNPs had significant associations with branch cluster frequency and seven SNPs with stem straightness across two of three trial series. One SNP had significant association with stem straightness across all three trial series. Strong G×E interactions resulted in fewer significant SNPs in common across the sites among which there were large effect differences. The percentage of genetic variance explained by SNPs showing significant associations ranged from 0.23 to 8.76 % for branch cluster frequency and from 0.37 to 12.75 % for stem straightness. This study showed that SNP effects for a trait change across environments if G×E interactions exist for that trait.