Comparative genomics reveals conservative evolution of the xylem transcriptome in vascular plants

<p>Abstract</p> <p>Background</p> <p>Wood is a valuable natural resource and a major carbon sink. Wood formation is an important developmental process in vascular plants which played a crucial role in plant evolution. Although genes involved in xylem formation have been investigated, the molecular mechanisms of xylem evolution are not well understood. We use comparative genomics to examine evolution of the xylem transcriptome to gain insights into xylem evolution.</p> <p>Results</p> <p>The xylem transcriptome is highly conserved in conifers, but considerably divergent in angiosperms. The functional domains of genes in the xylem transcriptome are moderately to highly conserved in vascular plants, suggesting the existence of a common ancestral xylem transcriptome. Compared to the total transcriptome derived from a range of tissues, the xylem transcriptome is relatively conserved in vascular plants. Of the xylem transcriptome, cell wall genes, ancestral xylem genes, known proteins and transcription factors are relatively more conserved in vascular plants. A total of 527 putative xylem orthologs were identified, which are unevenly distributed across the <it>Arabidopsis </it>chromosomes with eight hot spots observed. Phylogenetic analysis revealed that evolution of the xylem transcriptome has paralleled plant evolution. We also identified 274 conifer-specific xylem unigenes, all of which are of unknown function. These xylem orthologs and conifer-specific unigenes are likely to have played a crucial role in xylem evolution.</p> <p>Conclusions</p> <p>Conifers have highly conserved xylem transcriptomes, while angiosperm xylem transcriptomes are relatively diversified. Vascular plants share a common ancestral xylem transcriptome. The xylem transcriptomes of vascular plants are more conserved than the total transcriptomes. Evolution of the xylem transcriptome has largely followed the trend of plant evolution.</p

Generation and analysis of expressed sequence tags from six developing xylem libraries in Pinus radiata D. Don

<p>Abstract</p> <p>Background</p> <p>Wood is a major renewable natural resource for the timber, fibre and bioenergy industry. <it>Pinus radiata </it>D. Don is the most important commercial plantation tree species in Australia and several other countries; however, genomic resources for this species are very limited in public databases. Our primary objective was to sequence a large number of expressed sequence tags (ESTs) from genes involved in wood formation in radiata pine.</p> <p>Results</p> <p>Six developing xylem cDNA libraries were constructed from earlywood and latewood tissues sampled at juvenile (7 yrs), transition (11 yrs) and mature (30 yrs) ages, respectively. These xylem tissues represent six typical development stages in a rotation period of radiata pine. A total of 6,389 high quality ESTs were collected from 5,952 cDNA clones. Assembly of 5,952 ESTs from 5' end sequences generated 3,304 unigenes including 952 contigs and 2,352 singletons. About 97.0% of the 5,952 ESTs and 96.1% of the unigenes have matches in the UniProt and TIGR databases. Of the 3,174 unigenes with matches, 42.9% were not assigned GO (Gene Ontology) terms and their functions are unknown or unclassified. More than half (52.1%) of the 5,952 ESTs have matches in the Pfam database and represent 772 known protein families. About 18.0% of the 5,952 ESTs matched cell wall related genes in the MAIZEWALL database, representing all 18 categories, 91 of all 174 families and possibly 557 genes. Fifteen cell wall-related genes are ranked in the 30 most abundant genes, including <it>CesA</it>, <it>tubulin</it>, <it>AGP</it>, <it>SAMS</it>, <it>actin</it>, <it>laccase, CCoAMT, MetE</it>, <it>phytocyanin, pectate lyase</it>, <it>cellulase, SuSy</it>, <it>expansin</it>, <it>chitinase </it>and <it>UDP-glucose dehydrogenase</it>. Based on the PlantTFDB database 41 of the 64 transcription factor families in the poplar genome were identified as being involved in radiata pine wood formation. Comparative analysis of GO term abundance revealed a distinct transcriptome in juvenile earlywood formation compared to other stages of wood development.</p> <p>Conclusion</p> <p>The first large scale genomic resource in radiata pine was generated from six developing xylem cDNA libraries. Cell wall-related genes and transcription factors were identified. Juvenile earlywood has a distinct transcriptome, which is likely to contribute to the undesirable properties of juvenile wood in radiata pine. The publicly available resource of radiata pine will also be valuable for gene function studies and comparative genomics in forest trees.</p

Evolution in agriculture: the application of evolutionary approaches to the management of biotic interactions in agro-ecosystems

Anthropogenic impacts increasingly drive ecological and evolutionary processes at many spatio-temporal scales, demanding greater capacity to predict and manage their consequences. This is particularly true for agro-ecosystems, which not only comprise a significant proportion of land use, but which also involve conflicting imperatives to expand or intensify production while simultaneously reducing environmental impacts. These imperatives reinforce the likelihood of further major changes in agriculture over the next 30–40 years. Key transformations include genetic technologies as well as changes in land use. The use of evolutionary principles is not new in agriculture (e.g. crop breeding, domestication of animals, management of selection for pest resistance), but given land-use trends and other transformative processes in production landscapes, ecological and evolutionary research in agro-ecosystems must consider such issues in a broader systems context. Here, we focus on biotic interactions involving pests and pathogens as exemplars of situations where integration of agronomic, ecological and evolutionary perspectives has practical value. Although their presence in agro-ecosystems may be new, many traits involved in these associations evolved in natural settings. We advocate the use of predictive frameworks based on evolutionary models as pre-emptive management tools and identify some specific research opportunities to facilitate this. We conclude with a brief discussion of multidisciplinary approaches in applied evolutionary problems

Transcriptome profiling of Pinus radiata juvenile wood with contrasting stiffness identifies putative candidate genes involved in microfibril orientation and cell wall mechanics

<p>Abstract</p> <p>Background</p> <p>The mechanical properties of wood are largely determined by the orientation of cellulose microfibrils in secondary cell walls. Several genes and their allelic variants have previously been found to affect microfibril angle (MFA) and wood stiffness; however, the molecular mechanisms controlling microfibril orientation and mechanical strength are largely uncharacterised. In the present study, cDNA microarrays were used to compare gene expression in developing xylem with contrasting stiffness and MFA in juvenile <it>Pinus radiata </it>trees in order to gain further insights into the molecular mechanisms underlying microfibril orientation and cell wall mechanics.</p> <p>Results</p> <p>Juvenile radiata pine trees with higher stiffness (HS) had lower MFA in the earlywood and latewood of each ring compared to low stiffness (LS) trees. Approximately 3.4 to 14.5% out of 3, 320 xylem unigenes on cDNA microarrays were differentially regulated in juvenile wood with contrasting stiffness and MFA. Greater variation in MFA and stiffness was observed in earlywood compared to latewood, suggesting earlywood contributes most to differences in stiffness; however, 3-4 times more genes were differentially regulated in latewood than in earlywood. A total of 108 xylem unigenes were differentially regulated in juvenile wood with HS and LS in at least two seasons, including 43 unigenes with unknown functions. Many genes involved in cytoskeleton development and secondary wall formation (cellulose and lignin biosynthesis) were preferentially transcribed in wood with HS and low MFA. In contrast, several genes involved in cell division and primary wall synthesis were more abundantly transcribed in LS wood with high MFA.</p> <p>Conclusions</p> <p>Microarray expression profiles in <it>Pinus radiata </it>juvenile wood with contrasting stiffness has shed more light on the transcriptional control of microfibril orientation and the mechanical properties of wood. The identified candidate genes provide an invaluable resource for further gene function and association genetics studies aimed at deepening our understanding of cell wall biomechanics with a view to improving the mechanical properties of wood.</p

Transcriptome sequencing of <it>Eucalyptus camaldulensis</it> seedlings subjected to water stress reveals functional single nucleotide polymorphisms and genes under selection

Abstract Background Water stress limits plant survival and production in many parts of the world. Identification of genes and alleles responding to water stress conditions is important in breeding plants better adapted to drought. Currently there are no studies examining the transcriptome wide gene and allelic expression patterns under water stress conditions. We used RNA sequencing (RNA-seq) to identify the candidate genes and alleles and to explore the evolutionary signatures of selection. Results We studied the effect of water stress on gene expression in Eucalyptus camaldulensis seedlings derived from three natural populations. We used reference-guided transcriptome mapping to study gene expression. Several genes showed differential expression between control and stress conditions. Gene ontology (GO) enrichment tests revealed up-regulation of 140 stress-related gene categories and down-regulation of 35 metabolic and cell wall organisation gene categories. More than 190,000 single nucleotide polymorphisms (SNPs) were detected and 2737 of these showed differential allelic expression. Allelic expression of 52% of these variants was correlated with differential gene expression. Signatures of selection patterns were studied by estimating the proportion of nonsynonymous to synonymous substitution rates (Ka/Ks). The average Ka/Ks ratio among the 13,719 genes was 0.39 indicating that most of the genes are under purifying selection. Among the positively selected genes (Ka/Ks > 1.5) apoptosis and cell death categories were enriched. Of the 287 positively selected genes, ninety genes showed differential expression and 27 SNPs from 17 positively selected genes showed differential allelic expression between treatments. Conclusions Correlation of allelic expression of several SNPs with total gene expression indicates that these variants may be the cis-acting variants or in linkage disequilibrium with such variants. Enrichment of apoptosis and cell death gene categories among the positively selected genes reveals the past selection pressures experienced by the populations used in this study.</p

Transcriptome sequencing of Eucalyptus camaldulensis seedlings subjected to water stress reveals functional single nucleotide polymorphisms and genes under selection

BACKGROUND: Water stress limits plant survival and production in many parts of the world. Identification of genes and alleles responding to water stress conditions is important in breeding plants better adapted to drought. Currently there are no studies examining the transcriptome wide gene and allelic expression patterns under water stress conditions. We used RNA sequencing (RNA-seq) to identify the candidate genes and alleles and to explore the evolutionary signatures of selection. RESULTS: We studied the effect of water stress on gene expression in Eucalyptus camaldulensis seedlings derived from three natural populations. We used reference-guided transcriptome mapping to study gene expression. Several genes showed differential expression between control and stress conditions. Gene ontology (GO) enrichment tests revealed up-regulation of 140 stress-related gene categories and down-regulation of 35 metabolic and cell wall organisation gene categories. More than 190,000 single nucleotide polymorphisms (SNPs) were detected and 2737 of these showed differential allelic expression. Allelic expression of 52% of these variants was correlated with differential gene expression. Signatures of selection patterns were studied by estimating the proportion of nonsynonymous to synonymous substitution rates (Ka/Ks). The average Ka/Ks ratio among the 13,719 genes was 0.39 indicating that most of the genes are under purifying selection. Among the positively selected genes (Ka/Ks > 1.5) apoptosis and cell death categories were enriched. Of the 287 positively selected genes, ninety genes showed differential expression and 27 SNPs from 17 positively selected genes showed differential allelic expression between treatments. CONCLUSIONS: Correlation of allelic expression of several SNPs with total gene expression indicates that these variants may be the cis-acting variants or in linkage disequilibrium with such variants. Enrichment of apoptosis and cell death gene categories among the positively selected genes reveals the past selection pressures experienced by the populations used in this study

Eucalypt MADS-Box Genes Expressed in Developing Flowers

Three MADS-box genes were identified from a cDNA library derived from young flowers of Eucalyptus grandis W. Hill ex Maiden. The three egm genes are single-copy genes and are expressed almost exclusively in flowers. The egm1 and egm3 genes shared strongest homology with other plant MADS-box genes, which mediate between the floral meristem and the organ-identity genes. The egm3 gene was also expressed strongly in the receptacle or floral tube, which surrounds the carpels in the eucalypt flower and bears the sepals, petals, and numerous stamens. There appeared to be a group of genes in eucalypts with strong homology with the 3′ region of the egm1 gene. The egm2 gene was expressed in eucalypt petals and stamens and was most homologous to MADS-box genes, which belong to the globosa group of genes, which regulate organogenesis of the second and third floral whorls. The possible role of these three genes in eucalypt floral development is discussed

Fasciclin-like arabinogalactan proteins: specialization for stem biomechanics and cell wall architecture in Arabidopsis and Eucalyptus

The ancient cell adhesion fasciclin (FAS) domain is found in bacteria, fungi, algae, insects and animals, and occurs in a large family of fasciclin-like arabinogalactan proteins (FLAs) in higher plants. Functional roles for FAS-containing proteins have been determined for insects, algae and vertebrates; however, the biological functions of the various higher-plant FLAs are not clear. Expression of some FLAs has been correlated with the onset of secondary-wall cellulose synthesis in Arabidopsis stems, and also with wood formation in the stems and branches of trees, suggesting a biological role in plant stems. We examined whether FLAs contribute to plant stem biomechanics. Using phylogenetic, transcript abundance and promoter-GUS fusion analyses, we identified a conserved subset of single FAS domain FLAs (group A FLAs) in Eucalyptus and Arabidopsis that have specific and high transcript abundance in stems, particularly in stem cells undergoing secondary-wall deposition, and that the phylogenetic conservation appears to extend to other dicots and monocots. Gene-function analyses revealed that Arabidopsis T-DNA knockout double mutant stems had altered stem biomechanics with reduced tensile strength and a reduced tensile modulus of elasticity, as well as altered cell-wall architecture and composition, with increased cellulose microfibril angle and reduced arabinose, galactose and cellulose content. Using materials engineering concepts, we relate the effects of these FLAs on cell-wall composition with stem biomechanics. Our results suggest that a subset of single FAS domain FLAs contributes to plant stem strength by affecting cellulose deposition, and to the stem modulus of elasticity by affecting the integrity of the cell-wall matrix

Gene expression in Eucalyptus branch wood with marked variation in cellulose microfibril orientation and lacking G-layers

• In response to gravitational stresses, angiosperm trees form tension wood in the upper sides of branches and leaning stems in which cellulose content is higher, microfibrils are typically aligned closely with the fibre axis and the fibres often have

Quantitative trait locus (QTL) analysis of wood quality traits in Eucalyptus nitens

To identify the chromosomal regions affecting wood quality traits, we conducted a genome-wide quantitative trait locus (QTL) analysis of wood quality traits in Eucalyptus nitens. This information is important to exploit the full potential of the impending Eucalyptus genome sequence. A three generational mapping population consisting of 296 progeny trees was used to identify QTL associated with several wood quality traits in E. nitens. Thirty-six QTL positions for cellulose content, pulp yield, lignin content, density, and microfibril angle (MFA) were identified across different linkage groups. On linkage groups (LG)2 and 8, cellulose QTL cluster with pulp yield and extractives QTL while on LG4 and 10 cellulose and pulp yield QTLs cluster together. Similarly, on LG4, 5, and 6 QTL for lignin traits were clustered together. At two positions, QTL for MFA, a physical trait related to wood stiffness, were clustered with QTL for lignin traits. Several cell wall candidate genes were co-located to QTL positions affecting different traits. Comparative QTL analysis with Eucalyptus globulus revealed two common QTL regions for cellulose and pulp yield. The QTL positions identified in this study provide a resource for identifying wood quality genes using the impending Eucalyptus genome sequence. Candidate genes identified in this study through co-location to QTL regions may be useful in association studies