Search

The coconut rhinoceros beetle (Oryctes rhinoceros, CRB) is a serious pest of coconut and oil palms. It is native to South and Southeast Asia and was inadvertently introduced to Samoa in 1909. It has invaded many other Pacific countries throughout the last century. Oryctes rhinoceros nudivirus (OrNV), a natural pathogen of CRB in its native range, was successfully introduced as a classical biocontrol agent and has effectively suppressed invasive CRB populations for decades. However, resurgence of CRB has been recorded, with new invasions detected in several Pacific Island Countries and Territories. Additionally, new populations of CRB are emerging in some invaded areas that have a degree of resistance to the virus isolates commonly released for CRB biocontrol. Here, we designed a fast and reliable tool for distinguishing between different OrNV isolates that can help with the selection process to identify effective isolates for management of new CRB invasions. A comparison of 13 gene/gene region sequences within the OrNV genome of 16 OrNV isolates from native and invaded ranges allowed us to identify unique Single Nucleotide Polymorphisms (SNPs). With these SNPs, we developed an assay using multiplex PCR-amplicon-based nanopore sequencing to distinguish between OrNV isolates. We found that as few as four gene fragments were sufficient to identify 15 out of 20 OrNV isolates. This method can be used as a tool to monitor the establishment and distribution of OrNV isolates selected for release as biocontrol agents in CRB-infected areas., (© 2024. The Author(s).)

Oryctes rhinoceros nudivirus, a double-stranded DNA virus of the family Nudiviridae , is an important biocontrol agent of the coconut rhinoceros beetle (Coleoptera: Scarabaeidae). We present the genome sequences of six isolates of Oryctes rhinoceros nudivirus collected from the Philippines, Papua New Guinea, and Tanzania between the years 1977 and 2016., Competing Interests: The authors declare no conflict of interest.

This study compared the protein composition of M. longissimus thoracis of lambs from six commercial forage production systems in New Zealand. A total of 286 proteins were identified based on liquid chromatography-tandem mass spectrometry. First, a binomial model showed that different production groups could be distinguished based on abundances of 16 proteins. Second, pair-wise comparisons were performed to search for protein abundance differences in meat due to animal sex (ewe vs. wether), diet (perennial ryegrass vs. chicory), and age (4 vs. 6-8 months old). Greater abundance of some myofibrillar and sarcoplasmic proteins were observed in lamb loins from ewes compared to wethers. Chicory diet and older age at slaughter were associated with meat with lower abundance of some myofibrillar proteins, possibly due to a greater proportion of muscle glycolytic fibres. The proteins that showed significant differences in their abundances due to production factors could be further investigated to understand their influence on meat quality.

Curvature in mammalian fibers, such as wool and human hair, is an important feature of the functional trait of coat structure-it affects mechanical resilience and thermo-insulation. However, to examine the relationship between fiber curvature, ultrastructure and protein composition fiber diameter variability has to be minimal. To achieve this we utilised the progeny of straight-wool domestic sheep mutant rams (crimp mutants) and wild-type ewes. Proteomic and structural results of the resulting mutant/wild-type twin pairs confirmed that straight crimp mutant wool had a normal cuticle and the same cortical protein and ultrastructural building blocks as wild-type (crimpy) fibers but differed in the layout of its cortical cells and in the relative proportions of keratin (K) and keratin-associated proteins (KAPs). In the case of the crimp mutants (straight fibers), the orthocortex was distributed in a fragmented, annular ring, with some orthocortical cells near the central medulla, a pattern similar to that of straight hairs from humans and other mammals. Crimp mutant fibers were noted for the reduced abundance of some proteins in the high glycine-tyrosine class normally associated with the orthocortex, specifically the KAP6, KAP7, and KAP8 families, while proteins from the KAP16 and KAP19 were found in increased abundance. In addition to this, the type I keratin, K38, which is also associated with the orthocortex, was also found at lower abundance in the mutant fibers. Conversely, proteins from the ultra-high sulfur class normally associated with the paracortex, specifically the KAP4 and KAP9 families, were found in higher abundance., (© 2021 Wiley Periodicals LLC.)

A metagenomic approach was used to study the gut microbiome of Haemonchus contortus field strains and that of its predilection site, the abomasum of Dohne Merino sheep. The abomasum contents and H. contortus were collected from 10 naturally infected Dohne Merino sheep. The H. contortus specimens were classified and sexually differentiated using morphometric characters and was further confirmed through molecular identification. We investigated differences and similarities between the bacterial composition of the adult male and female H. contortus gut microbiomes, which were both dominated by bacteria from the Escherichia, Shigella, Vibrio and Halomonas genera. Major abundance variations were identified between the shared adult male and female H. contortus microbiomes. The results also revealed that Succiniclasticum, Rikenellaceae RC9 gut group and Candidatus Saccharimonas were the predominant genera in the Dohne Merino abomasum. This study provides insight into the highly diverse bacterial composition of the H. contortus gut microbiome and the Dohne Merino abomasum which needs to be studied further to explore the complex interactions of different gastrointestinal nematode microbiomes with the host.

Bacillus anthracis, the aetiological agent of anthrax, is regarded as a highly monomorphic pathogen that presents a low genetic diversity using standard molecular techniques. Whole genome sequencing and single nucleotide polymorphisms (SNPs) are definitive signatures for subtyping of B. anthracis. Here we employed whole genome single nucleotide polymorphism (wgSNP) analysis to investigate the genetic diversity of B. anthracis in the historically endemic region of Northern Cape Province (NCP), South Africa. Twenty-six isolates from anthrax outbreaks that occurred between 1998 and 2008/9 in NCP as well as from Namibia-South Africa Transfontier Conservation area and Botswana were compared to global B. anthracis genomes. Most NCP B. anthracis strains (n = 22) clustered in the A.Br.003/004 (A.Br.101) branch and are closely related to the Zimbabwe and Mozambique strains (A.Br.102 branch). A total of 4923 parsimony informative-SNPs accurately established the A.Br.003/004 phylogenetic relationships of the NCP isolates into two distinct sub-clades and SNP markers designated as A.Br.172 and A.Br.173 were developed. Other NCP strains (n = 2) grouped in the A.Br.001/002 (Sterne) branch while strains (n = 2) from the Namibia-South Africa Transfontier Conservation area and Botswana clustered in A.Br.005/006 (Ancient A) branch. The sequenced B. anthracis strains (A0094, A0096 and A0097) that clustered in the A.Br.064 (V770) clade were isolated from Vaalbos National Park and similar strains have not been isolated. The B. anthracis A0088 strain cluster with the NCP strains in the A.Br.003/004 (A.Br.172) SNP branch which has been isolated in NCP, South Africa. This study highlights the phylogenetic structure of NCP B. anthracis strains with distinctive SNP branches important for forensic tracing and novel SNP discovery purposes. The sequenced strains will serve as a means to further trace the dissemination of B. anthracis outbreaks in NCP, South Africa, and on the continent, as well as for forensic tracking on a global scale., Competing Interests: Declaration of Competing Interest The authors declare no competing interests., (Copyright © 2020 Elsevier B.V. All rights reserved.)

Background: Rhodococci are industrially important soil-dwelling Gram-positive bacteria that are well known for both nitrile hydrolysis and oxidative metabolism of aromatics. Rhodococcus rhodochrous ATCC BAA-870 is capable of metabolising a wide range of aliphatic and aromatic nitriles and amides. The genome of the organism was sequenced and analysed in order to better understand this whole cell biocatalyst., Results: The genome of R. rhodochrous ATCC BAA-870 is the first Rhodococcus genome fully sequenced using Nanopore sequencing. The circular genome contains 5.9 megabase pairs (Mbp) and includes a 0.53 Mbp linear plasmid, that together encode 7548 predicted protein sequences according to BASys annotation, and 5535 predicted protein sequences according to RAST annotation. The genome contains numerous oxidoreductases, 15 identified antibiotic and secondary metabolite gene clusters, several terpene and nonribosomal peptide synthetase clusters, as well as 6 putative clusters of unknown type. The 0.53 Mbp plasmid encodes 677 predicted genes and contains the nitrile converting gene cluster, including a nitrilase, a low molecular weight nitrile hydratase, and an enantioselective amidase. Although there are fewer biotechnologically relevant enzymes compared to those found in rhodococci with larger genomes, such as the well-known Rhodococcus jostii RHA1, the abundance of transporters in combination with the myriad of enzymes found in strain BAA-870 might make it more suitable for use in industrially relevant processes than other rhodococci., Conclusions: The sequence and comprehensive description of the R. rhodochrous ATCC BAA-870 genome will facilitate the additional exploitation of rhodococci for biotechnological applications, as well as enable further characterisation of this model organism. The genome encodes a wide range of enzymes, many with unknown substrate specificities supporting potential applications in biotechnology, including nitrilases, nitrile hydratase, monooxygenases, cytochrome P450s, reductases, proteases, lipases, and transaminases.

Introgression can introduce novel genetic variation at a faster rate than mutation alone and result in adaptive introgression when adaptive alleles are maintained in the recipient genome over time by natural selection. A previous study from our group demonstrated adaptive introgression from Populus balsamifera into P. trichocarpa in a target genomic region. Here we expand our local ancestry analysis to the whole genome of both parents to provide a comprehensive view of introgression patterns and to identify additional candidate regions for adaptive introgression genomewide. Populus trichocarpa is a large, fast-growing tree of mild coastal regions of the Pacific Northwest, whereas P. balsamifera is a smaller stature tree of continental and boreal regions with intense winter cold. The species hybridize where they are parapatric. We detected asymmetric patterns of introgression across the whole genome of these two poplar species adapted to contrasting environments, with stronger introgression from P. balsamifera to P. trichocarpa than vice versa. Admixed P. trichocarpa individuals contained more genomic regions with unusually high levels of introgression (19 regions) and also the largest introgressed genome fragment (1.02 Mb) compared with admixed P. balsamifera (nine regions). Our analysis also revealed numerous candidate regions for adaptive introgression with strong signals of selection, notably related to disease resistance, and enriched for genes that may play crucial roles in survival and adaptation. Furthermore, we detected a potential overrepresentation of subtelomeric regions in P. balsamifera introgressed into P. trichocarpa and possible protection of sex-determining regions from interspecific gene flow., (© 2018 John Wiley & Sons Ltd.)

Introgression can be an important source of adaptive phenotypes, although conversely it can have deleterious effects. Evidence for adaptive introgression is accumulating but information on the genetic architecture of introgressed traits lags behind. Here we determine trait architecture in Populus trichocarpa under introgression from P. balsamifera using admixture mapping and phenotypic analyses. Our results reveal that admixture is a key driver of clinal adaptation and suggest that the northern range extension of P. trichocarpa depends, at least in part, on introgression from P. balsamifera. However, admixture with P. balsamifera can lead to potentially maladaptive early phenology, and a reduction in growth and disease resistance in P. trichocarpa. Strikingly, an introgressed chromosome 9 haplotype block from P. balsamifera restores the late phenology and high growth parental phenotype in admixed P. trichocarpa. This epistatic restorer block may be strongly advantageous in maximizing carbon assimilation and disease resistance in the southernmost populations where admixture has been detected. We also confirm a previously demonstrated case of adaptive introgression in chromosome 15 and show that introgression generates a transgressive chlorophyll-content phenotype. We provide strong support that introgression provides a reservoir of genetic variation associated with adaptive characters that allows improved survival in new environments., (© 2017 The Authors. New Phytologist © 2017 New Phytologist Trust.)

Alkenes are linear hydrocarbons with one or more double bonds. Despite their potential as biofuels and precursors for specialty chemicals, the underlying biochemistry and genetics of alkene biosynthesis in plants remain elusive. Here, we report on a screen of natural accessions of poplar ( Populus trichocarpa ), revealing that the leaf cuticular waxes are predominantly composed of alkanes and alkenes. Interestingly, the accumulation of alkenes increases with leaf development, is limited to the abaxial side of the leaf, and is impaired in a few accessions. Among other genes, a β-ketoacyl CoA synthase gene ( PotriKCS1 ) was downregulated in leaves from non-alkene-producing accessions. We demonstrated biochemically that PotriKCS1 elongates monounsaturated fatty acids and is responsible for the recruitment of unsaturated substrates to the cuticular wax. Moreover, we found significant associations between the presence of alkenes and tree growth and resistance to leaf spot. These findings highlight the crucial role of cuticular waxes as the first point of contact with the environment, and they provide a foundation for engineering long-chain monounsaturated oils in other species., (© 2017 American Society of Plant Biologists. All rights reserved.)

Many plant genes are known to be involved in the development of cambium and wood, but how the expression and functional interaction of these genes determine the unique biology of wood remains largely unknown. We used the soc1ful loss of function mutant - the woodiest genotype known in the otherwise herbaceous model plant Arabidopsis - to investigate the expression and interactions of genes involved in secondary growth (wood formation). Detailed anatomical observations of the stem in combination with mRNA sequencing were used to assess transcriptome remodeling during xylogenesis in wild-type and woody soc1ful plants. To interpret the transcriptome changes, we constructed functional gene association networks of differentially expressed genes using the STRING database. This analysis revealed functionally enriched gene association hubs that are differentially expressed in herbaceous and woody tissues. In particular, we observed the differential expression of genes related to mechanical stress and jasmonate biosynthesis/signaling during wood formation in soc1ful plants that may be an effect of greater tension within woody tissues. Our results suggest that habit shifts from herbaceous to woody life forms observed in many angiosperm lineages could have evolved convergently by genetic changes that modulate the gene expression and interaction network, and thereby redeploy the conserved wood developmental program., (© 2016 The Authors. The Plant Journal published by Society for Experimental Biology and John Wiley & Sons Ltd.)

Natural hybrid zones in forest trees provide systems to study the transfer of adaptive genetic variation by introgression. Previous landscape genomic studies in Populus trichocarpa, a keystone tree species, indicated genomic footprints of admixture with its sister species Populus balsamifera and identified candidate genes for local adaptation. Here, we explored the patterns of introgression and signals of local adaptation in P. trichocarpa and P. balsamifera, employing genome resequencing data from three chromosomes in pure species and admixed individuals from wild populations. Local ancestry analysis in admixed P. trichocarpa revealed a telomeric region in chromosome 15 with P. balsamifera ancestry, containing several candidate genes for local adaptation. Genomic analyses revealed signals of selection in certain genes in this region (e.g. PRR5, COMT1), and functional analyses based on gene expression variation and correlations with adaptive phenotypes suggest distinct functions of the introgressed alleles. In contrast, a block of genes in chromosome 12 paralogous to the introgressed region showed no signs of introgression or signatures of selection. We hypothesize that the introgressed region in chromosome 15 has introduced modular or cassette-like variation into P. trichocarpa. These linked adaptive mutations are associated with a block of genes in chromosome 15 that appear to have undergone neo- or subfunctionalization relative to paralogs in a duplicated region on chromosome 12 that show no signatures of adaptive variation. The association between P. balsamifera introgressed alleles with the expression of adaptive traits in P. trichocarpa supports the hypothesis that this is a case of adaptive introgression in an ecologically important foundation species., (© 2016 John Wiley & Sons Ltd.)

Background: Recent studies show that galling Hymenoptera and Diptera are able to synthesize the plant hormone indole-3-acetic acid (auxin) from tryptophan and that plant response to insect-produced auxin is implicated in gall formation. We examined the leaf transcriptome of galled and ungalled leaves of individuals of the Hawaiian endemic plant Metrosideros polymorpha (Myrtaceae) subject to infestation by psyllid (Hemiptera) gall-makers in the genus Trioza (Triozidae)., Results: Transcript libraries were sequenced using Illumina technology and the reads assembled de novo into contigs. Functional identification of contigs followed a two-step procedure, first identifying contigs by comparison to the completely sequenced genome of the related Eucalyptus, followed by identifying the equivalent Arabidopsis gene using a pre-computed mapping between Eucalyptus and Arabidopsis genes. This allowed us to use the rich functional annotation of the Arabidopsis genome to assess the transcriptional landscape of galling in Metrosideros. Comparing galled and ungalled leaves, we find a highly significant enrichment of expressed genes with a gene ontology (GO) annotation to auxin response in the former. One gene consistently expressed in all galled trees examined but not detected in any libraries from ungalled leaves was the Metrosideros version of SMALL AUXIN UPREGULATED (SAUR) 67 which appears to be a marker for leaf-galling in Metrosideros., Conclusions: We conclude that an auxin response is involved in galling by Metrosideros psyllids. The possibility should therefore be considered that psyllids (like other insects examined) are able to synthesize auxin.

Background: Copy number variations (CNVs) are modifications in DNA structure comprising of deletions, duplications, insertions and complex multi-site variants. Although CNVs are proven to be involved in a variety of phenotypic discrepancies, the full extent and consequence of CNVs is yet to be understood. To date, no such genomic characterization has been performed in indigenous South African Nguni cattle. Nguni cattle are recognized for their ability to sustain harsh environmental conditions while exhibiting enhanced resistance to disease and parasites and are thought to comprise of up to nine different ecotypes., Methods: Illumina BovineSNP50 Beadchip data was utilized to investigate genomic population structure and the prevalence of CNVs in 492 South African Nguni cattle. PLINK, ADMIXTURE, R, gPLINK and Haploview software was utilized for quality control, population structure and haplotype block determination. PennCNV hidden Markov model identified CNVs and genes contained within and 10 Mb downstream from reported CNVs. PANTHER and Ensembl databases were subsequently utilized for gene annotation analyses., Results: Population structure analyses on Nguni cattle revealed 5 sub-populations with a possible sub-structure evident at K equal to 8. Four hundred and thirty three CNVs that formed 334 CNVRs ranging from 30 kb to 1 Mb in size are reported. Only 231 of the 492 animals demonstrated CNVRs. Two hundred and eighty nine genes were observed within CNVRs identified. Of these 149, 28, 44, 2 and 14 genes were unique to sub-populations A, B, C, D and E respectively. Gene ontology analyses demonstrated a number of pathways to be represented by respective genes, including immune response, response to abiotic stress and biological regulation processess., Conclusions: CNVs may explain part of the phenotypic diversity and the enhanced adaptation evident in Nguni cattle. Genes involved in a number of cellular components, biological processes and molecular functions are reported within CNVRs identified. The significance of such CNVRs and the possible effect thereof needs to be ascertained and may hold interesting insight into the functional and adaptive consequence of CNVs in cattle.

All species of the genus Populus (poplar, aspen) are dioecious, suggesting an ancient origin of this trait. Despite some empirical counter examples, theory suggests that nonrecombining sex-linked regions should quickly spread, eventually becoming heteromorphic chromosomes. In contrast, we show using whole-genome scans that the sex-associated region in Populus trichocarpa is small and much younger than the age of the genus. This indicates that sex determination is highly labile in poplar, consistent with recent evidence of 'turnover' of sex-determination regions in animals. We performed whole-genome resequencing of 52 P. trichocarpa (black cottonwood) and 34 Populus balsamifera (balsam poplar) individuals of known sex. Genomewide association studies in these unstructured populations identified 650 SNPs significantly associated with sex. We estimate the size of the sex-linked region to be ~100 kbp. All SNPs significantly associated with sex were in strong linkage disequilibrium despite the fact that they were mapped to six different chromosomes (plus 3 unmapped scaffolds) in version 2.2 of the reference genome. We show that this is likely due to genome misassembly. The segregation pattern of sex-associated SNPs revealed this to be an XY sex-determining system. Estimated divergence times of X and Y haplotype sequences (6-7 Ma) are much more recent than the divergence of P. trichocarpa (poplar) and Populus tremuloides (aspen). Consistent with this, in P. tremuloides, we found no XY haplotype divergence within the P. trichocarpa sex-determining region. These two species therefore have a different genomic architecture of sex, suggestive of at least one turnover event in the recent past., (© 2015 John Wiley & Sons Ltd.)

Tension wood has distinct physical and chemical properties, including altered fibre properties, cell wall composition and ultrastructure. It serves as a good system for investigating the genetic regulation of secondary cell wall biosynthesis and wood formation. The reference genome sequence for Eucalyptus grandis allows investigation of the global transcriptional reprogramming that accompanies tension wood formation in this global wood fibre crop. We report the first comprehensive analysis of physicochemical wood property changes in tension wood of Eucalyptus measured in a hybrid (E. grandis × Eucalyptus urophylla) clone, as well as genome-wide gene expression changes in xylem tissues 3 wk post-induction using RNA sequencing. We found that Eucalyptus tension wood in field-grown trees is characterized by an increase in cellulose, a reduction in lignin, xylose and mannose, and a marked increase in galactose. Gene expression profiling in tension wood-forming tissue showed corresponding down-regulation of monolignol biosynthetic genes, and differential expression of several carbohydrate active enzymes. We conclude that alterations of cell wall traits induced by tension wood formation in Eucalyptus are a consequence of a combination of down-regulation of lignin biosynthesis and hemicellulose remodelling, rather than the often proposed up-regulation of the cellulose biosynthetic pathway., (© 2014 University of Pretoria New Phytologist © 2014 New Phytologist Trust.)

As a step toward functional annotation of genes required for floral initiation and development within the Eucalyptus genome, we used short read sequencing to analyze transcriptomes of floral buds from early and late developmental stages, and compared these with transcriptomes of diverse vegetative tissues, including leaves, roots, and stems. A subset of 4807 genes (13% of protein-coding genes) were differentially expressed between floral buds of either stage and vegetative tissues. A similar proportion of genes were differentially expressed among all tissues. A total of 479 genes were differentially expressed between early and late stages of floral development. Gene function enrichment identified 158 gene ontology classes that were overrepresented in floral tissues, including 'pollen development' and 'aromatic compound biosynthetic process'. At least 40 floral-dominant genes lacked functional annotations and thus may be novel floral transcripts. We analyzed several genes and gene families in depth, including 49 putative biomarkers of floral development, the MADS-box transcription factors, 'S-domain'-receptor-like kinases, and selected gene family members with phosphatidylethanolamine-binding protein domains. Expanded MADS-box gene subfamilies in Eucalyptus grandis included SUPPRESSOR OF OVEREXPRESSION OF CO 1 (SOC1), SEPALLATA (SEP) and SHORT VEGETATIVE PHASE (SVP) Arabidopsis thaliana homologs. These data provide a rich resource for functional and evolutionary analysis of genes controlling eucalypt floral development, and new tools for breeding and biotechnology., (© 2014 The Authors. New Phytologist © 2014 New Phytologist Trust.)

NAC domain transcription factors regulate many developmental processes and stress responses in plants and vary widely in number and family structure. We analysed the characteristics and evolution of the NAC gene family of Eucalyptus grandis, a fast-growing forest tree in the rosid order Myrtales. NAC domain genes identified in the E. grandis genome were subjected to amino acid sequence, phylogenetic and motif analyses. Transcript abundance in developing tissues and abiotic stress conditions in E. grandis and E. globulus was quantified using RNA-seq and reverse transcription quantitative PCR (RT-qPCR). One hundred and eighty-nine E. grandis NAC (EgrNAC) proteins, arranged into 22 subfamilies, are extensively duplicated in subfamilies associated with stress response. Most EgrNAC genes form tandem duplicate arrays that frequently carry signatures of purifying selection. Sixteen amino acid motifs were identified in EgrNAC proteins, eight of which are enriched in, or unique to, Eucalyptus. New candidates for the regulation of normal and tension wood development and cold responses were identified. This first description of a Myrtales NAC domain family reveals an unique history of tandem duplication in stress-related subfamilies that has likely contributed to the adaptation of eucalypts to the challenging Australian environment. Several new candidates for the regulation of stress, wood formation and tree-specific development are reported., (© 2014 The Authors. New Phytologist © 2014 New Phytologist Trust.)

The R2R3-MYB family, one of the largest transcription factor families in higher plants, controls a wide variety of plant-specific processes including, notably, phenylpropanoid metabolism and secondary cell wall formation. We performed a genome-wide analysis of this superfamily in Eucalyptus, one of the most planted hardwood trees world-wide. A total of 141 predicted R2R3-MYB sequences identified in the Eucalyptus grandis genome sequence were subjected to comparative phylogenetic analyses with Arabidopsis thaliana, Oryza sativa, Populus trichocarpa and Vitis vinifera. We analysed features such as gene structure, conserved motifs and genome location. Transcript abundance patterns were assessed by RNAseq and validated by high-throughput quantitative PCR. We found some R2R3-MYB subgroups with expanded membership in E. grandis, V. vinifera and P. trichocarpa, and others preferentially found in woody species, suggesting diversification of specific functions in woody plants. By contrast, subgroups containing key genes regulating lignin biosynthesis and secondary cell wall formation are more conserved across all of the species analysed. In Eucalyptus, R2R3-MYB tandem gene duplications seem to disproportionately affect woody-preferential and woody-expanded subgroups. Interestingly, some of the genes belonging to woody-preferential subgroups show higher expression in the cambial region, suggesting a putative role in the regulation of secondary growth., (© 2014 The Authors New Phytologist © 2014 New Phytologist Trust.)

Wood formation is a complex developmental process governed by genetic and environmental stimuli. Populus and Eucalyptus are fast-growing, high-yielding tree genera that represent ecologically and economically important species suitable for generating significant lignocellulosic biomass. Comparative analysis of the developing xylem and leaf transcriptomes of Populus trichocarpa and Eucalyptus grandis together with phylogenetic analyses identified clusters of homologous genes preferentially expressed during xylem formation in both species. A conserved set of 336 single gene pairs showed highly similar xylem preferential expression patterns, as well as evidence of high functional constraint. Individual members of multi-gene orthologous clusters known to be involved in secondary cell wall biosynthesis also showed conserved xylem expression profiles. However, species-specific expression as well as opposite (xylem versus leaf) expression patterns observed for a subset of genes suggest subtle differences in the transcriptional regulation important for xylem development in each species. Using sequence similarity and gene expression status, we identified functional homologs likely to be involved in xylem developmental and biosynthetic processes in Populus and Eucalyptus. Our study suggests that, while genes involved in secondary cell wall biosynthesis show high levels of gene expression conservation, differential regulation of some xylem development genes may give rise to unique xylem properties., (© 2015 The Authors. New Phytologist © 2015 New Phytologist Trust.)

Background: Carbohydrate metabolism is a key feature of vascular plant architecture, and is of particular importance in large woody species, where lignocellulosic biomass is responsible for bearing the bulk of the stem and crown. Since Carbohydrate Active enZymes (CAZymes) in plants are responsible for the synthesis, modification and degradation of carbohydrate biopolymers, the differences in gene copy number and regulation between woody and herbaceous species have been highlighted previously. There are still many unanswered questions about the role of CAZymes in land plant evolution and the formation of wood, a strong carbohydrate sink., Results: Here, twenty-two publically available plant genomes were used to characterize the frequency, diversity and complexity of CAZymes in plants. We find that a conserved suite of CAZymes is a feature of land plant evolution, with similar diversity and complexity regardless of growth habit and form. In addition, we compared the diversity and levels of CAZyme gene expression during wood formation in trees using mRNA-seq data from two distantly related angiosperm tree species Eucalyptus grandis and Populus trichocarpa, highlighting the major CAZyme classes involved in xylogenesis and lignocellulosic biomass production., Conclusions: CAZyme domain ratio across embryophytes is maintained, and the diversity of CAZyme domains is similar in all land plants, regardless of woody habit. The stoichiometric conservation of gene expression in woody and non-woody tissues of Eucalyptus and Populus are indicative of gene balance preservation.

As molecular phylogenetic analyses incorporate ever-greater numbers of loci, cases of cytonuclear discordance - the phenomenon in which nuclear gene trees deviate significantly from organellar gene trees - are being reported more frequently. Plant examples of topological discordance, caused by recent hybridization between extant species, are well known. However, examples of branch-length discordance are less reported in plants relative to animals. We use a combination of de novo assembly and reference-based mapping using short-read shotgun sequences to construct a robust phylogeny of the plastome for multiple individuals of all the common Populus species in North America. We demonstrate a case of strikingly high plastome divergence, in contrast to little nuclear genome divergence, in two closely related balsam poplars, Populus balsamifera and Populus trichocarpa (Populus balsamifera ssp. trichocarpa). Previous studies with nuclear loci indicate that the two species (or subspecies) diverged since the late Pleistocene, whereas their plastomes indicate deep divergence, dating to at least the Pliocene (6-7 Myr ago). Our finding is in marked contrast to the estimated Pleistocene divergence of the nuclear genomes, previously calculated at 75 000 yr ago, suggesting plastid capture from a 'ghost lineage' of a now-extinct North American poplar., (© 2014 The Authors. New Phytologist © 2014 New Phytologist Trust.)

Eucalypts are the world's most widely planted hardwood trees. Their outstanding diversity, adaptability and growth have made them a global renewable resource of fibre and energy. We sequenced and assembled >94% of the 640-megabase genome of Eucalyptus grandis. Of 36,376 predicted protein-coding genes, 34% occur in tandem duplications, the largest proportion thus far in plant genomes. Eucalyptus also shows the highest diversity of genes for specialized metabolites such as terpenes that act as chemical defence and provide unique pharmaceutical oils. Genome sequencing of the E. grandis sister species E. globulus and a set of inbred E. grandis tree genomes reveals dynamic genome evolution and hotspots of inbreeding depression. The E. grandis genome is the first reference for the eudicot order Myrtales and is placed here sister to the eurosids. This resource expands our understanding of the unique biology of large woody perennials and provides a powerful tool to accelerate comparative biology, breeding and biotechnology.

Background: F1 hybrid clones of Eucalyptus grandis and E. urophylla are widely grown for pulp and paper production in tropical and subtropical regions. Volume growth and wood quality are priority objectives in Eucalyptus tree improvement. The molecular basis of quantitative variation and trait expression in eucalypt hybrids, however, remains largely unknown. The recent availability of a draft genome sequence (http://www.phytozome.net) and genome-wide genotyping platforms, combined with high levels of genetic variation and high linkage disequilibrium in hybrid crosses, greatly facilitate the detection of quantitative trait loci (QTLs) as well as underlying candidate genes for growth and wood property traits. In this study, we used Diversity Arrays Technology markers to assess the genetic architecture of volume growth (diameter at breast height, DBH) and wood basic density in four-year-old progeny of an interspecific backcross pedigree of E. grandis and E. urophylla. In addition, we used Illumina RNA-Seq expression profiling in the E. urophylla backcross family to identify cis- and trans-acting polymorphisms (eQTLs) affecting transcript abundance of genes underlying QTLs for wood basic density., Results: A total of five QTLs for DBH and 12 for wood basic density were identified in the two backcross families. Individual QTLs for DBH and wood basic density explained 3.1 to 12.2% of phenotypic variation. Candidate genes underlying QTLs for wood basic density on linkage groups 8 and 9 were found to share trans-acting eQTLs located on linkage groups 4 and 10, which in turn coincided with QTLs for wood basic density suggesting that these QTLs represent segregating components of an underlying transcriptional network., Conclusion: This is the first demonstration of the use of next-generation expression profiling to quantify transcript abundance in a segregating tree population and identify candidate genes potentially affecting wood property variation. The QTLs identified in this study provide a resource for identifying candidate genes and developing molecular markers for marker-assisted breeding of volume growth and wood basic density. Our results suggest that integrated analysis of transcript and trait variation in eucalypt hybrids can be used to dissect the molecular basis of quantitative variation in wood property traits.

Background: De novo assembly of transcript sequences produced by short-read DNA sequencing technologies offers a rapid approach to obtain expressed gene catalogs for non-model organisms. A draft genome sequence will be produced in 2010 for a Eucalyptus tree species (E. grandis) representing the most important hardwood fibre crop in the world. Genome annotation of this valuable woody plant and genetic dissection of its superior growth and productivity will be greatly facilitated by the availability of a comprehensive collection of expressed gene sequences from multiple tissues and organs., Results: We present an extensive expressed gene catalog for a commercially grown E. grandis × E. urophylla hybrid clone constructed using only Illumina mRNA-Seq technology and de novo assembly. A total of 18,894 transcript-derived contigs, a large proportion of which represent full-length protein coding genes were assembled and annotated. Analysis of assembly quality, length and diversity show that this dataset represent the most comprehensive expressed gene catalog for any Eucalyptus tree. mRNA-Seq analysis furthermore allowed digital expression profiling of all of the assembled transcripts across diverse xylogenic and non-xylogenic tissues, which is invaluable for ascribing putative gene functions., Conclusions: De novo assembly of Illumina mRNA-Seq reads is an efficient approach for transcriptome sequencing and profiling in Eucalyptus and other non-model organisms. The transcriptome resource (Eucspresso, http://eucspresso.bi.up.ac.za/) generated by this study will be of value for genomic analysis of woody biomass production in Eucalyptus and for comparative genomic analysis of growth and development in woody and herbaceous plants.