1
|
Simon SJ, Furches A, Chhetri H, Evans L, Abeyratne CR, Jones P, Wimp G, Macaya-Sanz D, Jacobson D, Tschaplinski TJ, Tuskan GA, DiFazio SP. Genetic underpinnings of arthropod community distributions in Populus trichocarpa. New Phytol 2024; 242:1307-1323. [PMID: 38488269 DOI: 10.1111/nph.19660] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 10/18/2023] [Accepted: 02/21/2024] [Indexed: 04/12/2024]
Abstract
Community genetics seeks to understand the mechanisms by which natural genetic variation in heritable host phenotypes can encompass assemblages of organisms such as bacteria, fungi, and many animals including arthropods. Prior studies that focused on plant genotypes have been unable to identify genes controlling community composition, a necessary step to predict ecosystem structure and function as underlying genes shift within plant populations. We surveyed arthropods within an association population of Populus trichocarpa in three common gardens to discover plant genes that contributed to arthropod community composition. We analyzed our surveys with traditional single-trait genome-wide association analysis (GWAS), multitrait GWAS, and functional networks built from a diverse set of plant phenotypes. Plant genotype was influential in structuring arthropod community composition among several garden sites. Candidate genes important for higher level organization of arthropod communities had broadly applicable functions, such as terpenoid biosynthesis and production of dsRNA binding proteins and protein kinases, which may be capable of targeting multiple arthropod species. We have demonstrated the ability to detect, in an uncontrolled environment, individual genes that are associated with the community assemblage of arthropods on a host plant, further enhancing our understanding of genetic mechanisms that impact ecosystem structure.
Collapse
Affiliation(s)
- Sandra J Simon
- Department of Biology, West Virginia University, Morgantown, WV, 26506, USA
| | - Anna Furches
- Biosciences Division and Center for Bioenergy Innovation, Oak Ridge National Laboratory, Oak Ridge, TN, 37831, USA
- The Bredesen Center for Interdisciplinary Research and Graduate Education, University of Tennessee, Knoxville, TN, 37996, USA
| | - Hari Chhetri
- Department of Biology, West Virginia University, Morgantown, WV, 26506, USA
- Computational Systems Biology Group, Oak Ridge National Laboratory, Oak Ridge, TN, 37831, USA
| | - Luke Evans
- Institute for Behavioral Genetics, University of Colorado Boulder, Boulder, CO, 80309, USA
| | | | - Piet Jones
- Biosciences Division and Center for Bioenergy Innovation, Oak Ridge National Laboratory, Oak Ridge, TN, 37831, USA
- The Bredesen Center for Interdisciplinary Research and Graduate Education, University of Tennessee, Knoxville, TN, 37996, USA
| | - Gina Wimp
- Department of Biology, Georgetown University, Washington, DC, 20057, USA
| | - David Macaya-Sanz
- Department of Biology, West Virginia University, Morgantown, WV, 26506, USA
| | - Daniel Jacobson
- Biosciences Division and Center for Bioenergy Innovation, Oak Ridge National Laboratory, Oak Ridge, TN, 37831, USA
- The Bredesen Center for Interdisciplinary Research and Graduate Education, University of Tennessee, Knoxville, TN, 37996, USA
| | - Timothy J Tschaplinski
- Biosciences Division and Center for Bioenergy Innovation, Oak Ridge National Laboratory, Oak Ridge, TN, 37831, USA
| | - Gerald A Tuskan
- Biosciences Division and Center for Bioenergy Innovation, Oak Ridge National Laboratory, Oak Ridge, TN, 37831, USA
| | - Stephen P DiFazio
- Department of Biology, West Virginia University, Morgantown, WV, 26506, USA
| |
Collapse
|
2
|
Macaya-Sanz D, Witzell J, Collada C, Gil L, Martín JA. Core endophytic mycobiome in Ulmus minor and its relation to Dutch elm disease resistance. Front Plant Sci 2023; 14:1125942. [PMID: 36925756 PMCID: PMC10011445 DOI: 10.3389/fpls.2023.1125942] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Grants] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 12/16/2022] [Accepted: 02/13/2023] [Indexed: 06/18/2023]
Abstract
The core microbiota of plants exerts key effects on plant performance and resilience to stress. The aim of this study was to identify the core endophytic mycobiome in U. minor stems and disentangle associations between its composition and the resistance to Dutch elm disease (DED). We also defined its spatial variation within the tree and among distant tree populations. Stem samples were taken i) from different heights of the crown of a 168-year-old elm tree, ii) from adult elm trees growing in a common garden and representing a gradient of resistance to DED, and iii) from trees growing in two distant natural populations, one of them with varying degrees of vitality. Endophyte composition was profiled by high throughput sequencing of the first internal transcribed spacer region (ITS1) of the ribosomal DNA. Three families of yeasts (Buckleyzymaceae, Trichomeriaceae and Bulleraceae) were associated to DED-resistant hosts. A small proportion (10%) of endophytic OTUs was almost ubiquitous throughout the crown while tree colonization by most fungal taxa followed stochastic patterns. A clear distinction in endophyte composition was found between geographical locations. By combining all surveys, we found evidence of a U. minor core mycobiome, pervasive within the tree and ubiquitous across locations, genotypes and health status.
Collapse
Affiliation(s)
- David Macaya-Sanz
- Departamento de Ecología y Genética Forestal, Instituto de Ciencias Forestales (ICIFOR-INIA), CSIC, Madrid, Spain
| | - Johanna Witzell
- Department of Forestry and Wood Technology, Linnaeus University, Växjö, Sweden
| | - Carmen Collada
- Departamento de Sistemas y Recursos Naturales, Escuela Técnica Superior de Ingeniería (ETSI) Montes, Forestal y del Medio Natural, Universidad Politécnica de Madrid, Madrid, Spain
| | - Luis Gil
- Departamento de Sistemas y Recursos Naturales, Escuela Técnica Superior de Ingeniería (ETSI) Montes, Forestal y del Medio Natural, Universidad Politécnica de Madrid, Madrid, Spain
| | - Juan A. Martín
- Departamento de Sistemas y Recursos Naturales, Escuela Técnica Superior de Ingeniería (ETSI) Montes, Forestal y del Medio Natural, Universidad Politécnica de Madrid, Madrid, Spain
| |
Collapse
|
3
|
Abeyratne CR, Macaya-Sanz D, Zhou R, Barry KW, Daum C, Haiby K, Lipzen A, Stanton B, Yoshinaga Y, Zane M, Tuskan GA, DiFazio SP. High-resolution mapping reveals hotspots and sex-biased recombination in Populus trichocarpa. G3 (Bethesda) 2022; 13:6762080. [PMID: 36250890 PMCID: PMC9836356 DOI: 10.1093/g3journal/jkac269] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 05/10/2022] [Accepted: 09/28/2022] [Indexed: 12/14/2022]
Abstract
Fine-scale meiotic recombination is fundamental to the outcome of natural and artificial selection. Here, dense genetic mapping and haplotype reconstruction were used to estimate recombination for a full factorial Populus trichocarpa cross of 7 males and 7 females. Genomes of the resulting 49 full-sib families (N = 829 offspring) were resequenced, and high-fidelity biallelic SNP/INDELs and pedigree information were used to ascertain allelic phase and impute progeny genotypes to recover gametic haplotypes. The 14 parental genetic maps contained 1,820 SNP/INDELs on average that covered 376.7 Mb of physical length across 19 chromosomes. Comparison of parental and progeny haplotypes allowed fine-scale demarcation of cross-over regions, where 38,846 cross-over events in 1,658 gametes were observed. Cross-over events were positively associated with gene density and negatively associated with GC content and long-terminal repeats. One of the most striking findings was higher rates of cross-overs in males in 8 out of 19 chromosomes. Regions with elevated male cross-over rates had lower gene density and GC content than windows showing no sex bias. High-resolution analysis identified 67 candidate cross-over hotspots spread throughout the genome. DNA sequence motifs enriched in these regions showed striking similarity to those of maize, Arabidopsis, and wheat. These findings, and recombination estimates, will be useful for ongoing efforts to accelerate domestication of this and other biomass feedstocks, as well as future studies investigating broader questions related to evolutionary history, perennial development, phenology, wood formation, vegetative propagation, and dioecy that cannot be studied using annual plant model systems.
Collapse
Affiliation(s)
| | - David Macaya-Sanz
- Department of Forest Ecology & Genetics, CIFOR-INIA, CSIC, Madrid 28040, Spain
| | - Ran Zhou
- Warnell School of Forestry and Natural Resources, Department of Genetics, and Department of Plant Biology, University of Georgia, Athens, GA 30602, USA
| | - Kerrie W Barry
- Department of Energy Joint Genome Institute, Berkeley, CA 94720, USA
| | - Christopher Daum
- Department of Energy Joint Genome Institute, Berkeley, CA 94720, USA
| | | | - Anna Lipzen
- Department of Energy Joint Genome Institute, Berkeley, CA 94720, USA
| | | | - Yuko Yoshinaga
- Department of Energy Joint Genome Institute, Berkeley, CA 94720, USA
| | - Matthew Zane
- Department of Energy Joint Genome Institute, Berkeley, CA 94720, USA
| | - Gerald A Tuskan
- Biosciences Division, Center for Bioenergy Innovation, Oak Ridge National Laboratory, Oak Ridge, TN 37830, USA
| | - Stephen P DiFazio
- Corresponding author: Department of Biology, West Virginia University, Morgantown, WV 26506, USA.
| |
Collapse
|
4
|
Harman-Ware AE, Happs RM, Macaya-Sanz D, Doeppke C, Muchero W, DiFazio SP. Abundance of Major Cell Wall Components in Natural Variants and Pedigrees of Populus trichocarpa. Front Plant Sci 2022; 13:757810. [PMID: 35185975 PMCID: PMC8850957 DOI: 10.3389/fpls.2022.757810] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Grants] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 08/12/2021] [Accepted: 01/04/2022] [Indexed: 06/14/2023]
Abstract
The rapid analysis of biopolymers including lignin and sugars in lignocellulosic biomass cell walls is essential for the analysis of the large sample populations needed for identifying heritable genetic variation in biomass feedstocks for biofuels and bioproducts. In this study, we reported the analysis of cell wall lignin content, syringyl/guaiacyl (S/G) ratio, as well as glucose and xylose content by high-throughput pyrolysis-molecular beam mass spectrometry (py-MBMS) for >3,600 samples derived from hundreds of accessions of Populus trichocarpa from natural populations, as well as pedigrees constructed from 14 parents (7 × 7). Partial Least Squares (PLS) regression models were built from the samples of known sugar composition previously determined by hydrolysis followed by nuclear magnetic resonance (NMR) analysis. Key spectral features positively correlated with glucose content consisted of m/z 126, 98, and 69, among others, deriving from pyrolyzates such as hydroxymethylfurfural, maltol, and other sugar-derived species. Xylose content positively correlated primarily with many lignin-derived ions and to a lesser degree with m/z 114, deriving from a lactone produced from xylose pyrolysis. Models were capable of predicting glucose and xylose contents with an average error of less than 4%, and accuracy was significantly improved over previously used methods. The differences in the models constructed from the two sample sets varied in training sample number, but the genetic and compositional uniformity of the pedigree set could be a potential driver in the slightly better performance of that model in comparison with the natural variants. Broad-sense heritability of glucose and xylose composition using these data was 0.32 and 0.34, respectively. In summary, we have demonstrated the use of a single high-throughput method to predict sugar and lignin composition in thousands of poplar samples to estimate the heritability and phenotypic plasticity of traits necessary to develop optimized feedstocks for bioenergy applications.
Collapse
Affiliation(s)
- Anne E. Harman-Ware
- Renewable Resources and Enabling Sciences Center, Center for Bioenergy Innovation, National Renewable Energy Laboratory, Golden, CO, United States
| | - Renee M. Happs
- Renewable Resources and Enabling Sciences Center, Center for Bioenergy Innovation, National Renewable Energy Laboratory, Golden, CO, United States
| | - David Macaya-Sanz
- Department of Forest Ecology and Genetics, INIA-CIFOR, Madrid, Spain
- Department of Biology, West Virginia University, Morgantown, WV, United States
| | - Crissa Doeppke
- Renewable Resources and Enabling Sciences Center, Center for Bioenergy Innovation, National Renewable Energy Laboratory, Golden, CO, United States
| | - Wellington Muchero
- Biosciences Division, Center for Bioenergy Innovation, Oak Ridge National Laboratory, Oak Ridge, TN, United States
| | - Stephen P. DiFazio
- Department of Biology, West Virginia University, Morgantown, WV, United States
| |
Collapse
|
5
|
Harman-Ware AE, Macaya-Sanz D, Abeyratne CR, Doeppke C, Haiby K, Tuskan GA, Stanton B, DiFazio SP, Davis MF. Accurate determination of genotypic variance of cell wall characteristics of a Populus trichocarpa pedigree using high-throughput pyrolysis-molecular beam mass spectrometry. Biotechnol Biofuels 2021; 14:59. [PMID: 33676543 PMCID: PMC7937246 DOI: 10.1186/s13068-021-01908-y] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Grants] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 04/16/2020] [Accepted: 02/17/2021] [Indexed: 06/12/2023]
Abstract
BACKGROUND Pyrolysis-molecular beam mass spectrometry (py-MBMS) analysis of a pedigree of Populus trichocarpa was performed to study the phenotypic plasticity and heritability of lignin content and lignin monomer composition. Instrumental and microspatial environmental variability were observed in the spectral features and corrected to reveal underlying genetic variance of biomass composition. RESULTS Lignin-derived ions (including m/z 124, 154, 168, 194, 210 and others) were highly impacted by microspatial environmental variation which demonstrates phenotypic plasticity of lignin composition in Populus trichocarpa biomass. Broad-sense heritability of lignin composition after correcting for microspatial and instrumental variation was determined to be H2 = 0.56 based on py-MBMS ions known to derive from lignin. Heritability of lignin monomeric syringyl/guaiacyl ratio (S/G) was H2 = 0.81. Broad-sense heritability was also high (up to H2 = 0.79) for ions derived from other components of the biomass including phenolics (e.g., salicylates) and C5 sugars (e.g., xylose). Lignin and phenolic ion abundances were primarily driven by maternal effects, and paternal effects were either similar or stronger for the most heritable carbohydrate-derived ions. CONCLUSIONS We have shown that many biopolymer-derived ions from py-MBMS show substantial phenotypic plasticity in response to microenvironmental variation in plantations. Nevertheless, broad-sense heritability for biomass composition can be quite high after correcting for spatial environmental variation. This work outlines the importance in accounting for instrumental and microspatial environmental variation in biomass composition data for applications in heritability measurements and genomic selection for breeding poplar for renewable fuels and materials.
Collapse
Affiliation(s)
- Anne E Harman-Ware
- Renewable Resources and Enabling Sciences Center, National Renewable Energy Laboratory, Golden, CO, 80401, USA.
| | - David Macaya-Sanz
- Department of Biology, West Virginia University, Morgantown, WV, 26506, USA
| | | | - Crissa Doeppke
- Renewable Resources and Enabling Sciences Center, National Renewable Energy Laboratory, Golden, CO, 80401, USA
| | | | - Gerald A Tuskan
- Biosciences Division, Oak Ridge National Laboratory, Oak Ridge, TN, 37830, USA
| | | | - Stephen P DiFazio
- Department of Biology, West Virginia University, Morgantown, WV, 26506, USA
| | - Mark F Davis
- Biosciences Center, National Renewable Energy Laboratory, Golden, CO, 80401, USA
| |
Collapse
|
6
|
Zhou R, Macaya-Sanz D, Schmutz J, Jenkins JW, Tuskan GA, DiFazio SP. Sequencing and Analysis of the Sex Determination Region of Populus trichocarpa. Genes (Basel) 2020; 11:E843. [PMID: 32722098 PMCID: PMC7465354 DOI: 10.3390/genes11080843] [Citation(s) in RCA: 15] [Impact Index Per Article: 3.8] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/01/2020] [Revised: 07/17/2020] [Accepted: 07/22/2020] [Indexed: 12/19/2022] Open
Abstract
The ages and sizes of a sex-determination region (SDR) are difficult to determine in non-model species. Due to the lack of recombination and enrichment of repetitive elements in SDRs, the quality of assembly with short sequencing reads is universally low. Unique features present in the SDRs help provide clues about how SDRs are established and how they evolve in the absence of recombination. Several Populus species have been reported with a male heterogametic configuration of sex (XX/XY system) mapped on chromosome 19, but the exact location of the SDR has been inconsistent among species, and thus far, none of these SDRs has been fully assembled in a genomic context. Here we identify the Y-SDR from a Y-linked contig directly from a long-read PacBio assembly of a Populus trichocarpa male individual. We also identified homologous gene sequences in the SDR of P. trichocarpa and the SDR of the W chromosome in Salix purpurea. We show that inverted repeats (IRs) found in the Y-SDR and the W-SDR are lineage-specific. We hypothesize that, although the two IRs are derived from the same orthologous gene within each species, they likely have independent evolutionary histories. Furthermore, the truncated inverted repeats in P. trichocarpa may code for small RNAs that target the homologous gene for RNA-directed DNA methylation. These findings support the hypothesis that diverse sex-determining systems may be achieved through similar evolutionary pathways, thereby providing a possible mechanism to explain the lability of sex-determination systems in plants in general.
Collapse
Affiliation(s)
- Ran Zhou
- Department of Biology, West Virginia University, Morgantown, WV 26506-6057, USA; (R.Z.); (D.M.-S.)
| | - David Macaya-Sanz
- Department of Biology, West Virginia University, Morgantown, WV 26506-6057, USA; (R.Z.); (D.M.-S.)
| | - Jeremy Schmutz
- HudsonAlpha Institute of Biotechnology, 601 Genome Way Northwest, Huntsville, AL 35806, USA; (J.S.); (J.W.J.)
- Department of Energy Joint Genome Institute, 2800 Mitchell Drive, Walnut Creek, CA 94598, USA;
| | - Jerry W. Jenkins
- HudsonAlpha Institute of Biotechnology, 601 Genome Way Northwest, Huntsville, AL 35806, USA; (J.S.); (J.W.J.)
| | - Gerald A. Tuskan
- Department of Energy Joint Genome Institute, 2800 Mitchell Drive, Walnut Creek, CA 94598, USA;
- Center for Bioenergy Innovation, Biosciences Division, Oak Ridge National Lab, Oak Ridge, TN 37830, USA
| | - Stephen P. DiFazio
- Department of Biology, West Virginia University, Morgantown, WV 26506-6057, USA; (R.Z.); (D.M.-S.)
| |
Collapse
|
7
|
Zhou R, Macaya-Sanz D, Carlson CH, Schmutz J, Jenkins JW, Kudrna D, Sharma A, Sandor L, Shu S, Barry K, Tuskan GA, Ma T, Liu J, Olson M, Smart LB, DiFazio SP. A willow sex chromosome reveals convergent evolution of complex palindromic repeats. Genome Biol 2020; 21:38. [PMID: 32059685 PMCID: PMC7023750 DOI: 10.1186/s13059-020-1952-4] [Citation(s) in RCA: 49] [Impact Index Per Article: 12.3] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/23/2019] [Accepted: 02/03/2020] [Indexed: 12/19/2022] Open
Abstract
Background Sex chromosomes have arisen independently in a wide variety of species, yet they share common characteristics, including the presence of suppressed recombination surrounding sex determination loci. Mammalian sex chromosomes contain multiple palindromic repeats across the non-recombining region that show sequence conservation through gene conversion and contain genes that are crucial for sexual reproduction. In plants, it is not clear if palindromic repeats play a role in maintaining sequence conservation in the absence of homologous recombination. Results Here we present the first evidence of large palindromic structures in a plant sex chromosome, based on a highly contiguous assembly of the W chromosome of the dioecious shrub Salix purpurea. The W chromosome has an expanded number of genes due to transpositions from autosomes. It also contains two consecutive palindromes that span a region of 200 kb, with conspicuous 20-kb stretches of highly conserved sequences among the four arms that show evidence of gene conversion. Four genes in the palindrome are homologous to genes in the sex determination regions of the closely related genus Populus, which is located on a different chromosome. These genes show distinct, floral-biased expression patterns compared to paralogous copies on autosomes. Conclusion The presence of palindromes in sex chromosomes of mammals and plants highlights the intrinsic importance of these features in adaptive evolution in the absence of recombination. Convergent evolution is driving both the independent establishment of sex chromosomes as well as their fine-scale sequence structure.
Collapse
Affiliation(s)
- Ran Zhou
- Department of Biology, West Virginia University, Morgantown, WV, 26506-6057, USA
| | - David Macaya-Sanz
- Department of Biology, West Virginia University, Morgantown, WV, 26506-6057, USA
| | - Craig H Carlson
- Horticulture Section, School of Integrative Plant Science, Cornell University, New York State Agricultural Experiment Station, Geneva, NY, 14456, USA
| | - Jeremy Schmutz
- HudsonAlpha Institute of Biotechnology, Huntsville, AL, USA.,Department of Energy Joint Genome Institute, Walnut Creek, California, USA
| | | | - David Kudrna
- Arizona Genomics Institute, School of Plant Sciences, University of Arizona, Tucson, AZ, USA
| | - Aditi Sharma
- Department of Energy Joint Genome Institute, Walnut Creek, California, USA
| | - Laura Sandor
- Department of Energy Joint Genome Institute, Walnut Creek, California, USA
| | - Shengqiang Shu
- Department of Energy Joint Genome Institute, Walnut Creek, California, USA
| | - Kerrie Barry
- Department of Energy Joint Genome Institute, Walnut Creek, California, USA
| | - Gerald A Tuskan
- Biosciences Division, Oak Ridge National Laboratory, Oak Ridge, TN, 37831, USA.,DOE-Center for Bioenergy Innovation (CBI), Oak Ridge National Laboratory, Oak Ridge, TN, 37831, USA
| | - Tao Ma
- Key Laboratory of Bio-Resource and Eco-Environment of Ministry of Education, College of Life Sciences, Sichuan University, Chengdu, 610065, China
| | - Jianquan Liu
- Key Laboratory of Bio-Resource and Eco-Environment of Ministry of Education, College of Life Sciences, Sichuan University, Chengdu, 610065, China.,State Key Laboratory of Grassland Agro-Ecosystem, Institute of Innovation Ecology & College of Life Sciences, Lanzhou University, Lanzhou, 730000, China
| | - Matthew Olson
- Department of Biological Sciences, Texas Tech University, Box 43131, Lubbock, TX, 79409-3131, USA
| | - Lawrence B Smart
- Horticulture Section, School of Integrative Plant Science, Cornell University, New York State Agricultural Experiment Station, Geneva, NY, 14456, USA
| | - Stephen P DiFazio
- Department of Biology, West Virginia University, Morgantown, WV, 26506-6057, USA.
| |
Collapse
|
8
|
Chhetri HB, Furches A, Macaya-Sanz D, Walker AR, Kainer D, Jones P, Harman-Ware AE, Tschaplinski TJ, Jacobson D, Tuskan GA, DiFazio SP. Genome-Wide Association Study of Wood Anatomical and Morphological Traits in Populus trichocarpa. Front Plant Sci 2020; 11:545748. [PMID: 33013968 PMCID: PMC7509168 DOI: 10.3389/fpls.2020.545748] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 03/26/2020] [Accepted: 08/21/2020] [Indexed: 05/04/2023]
Abstract
To understand the genetic mechanisms underlying wood anatomical and morphological traits in Populus trichocarpa, we used 869 unrelated genotypes from a common garden in Clatskanie, Oregon that were previously collected from across the distribution range in western North America. Using GEMMA mixed model analysis, we tested for the association of 25 phenotypic traits and nine multitrait combinations with 6.741 million SNPs covering the entire genome. Broad-sense trait heritabilities ranged from 0.117 to 0.477. Most traits were significantly correlated with geoclimatic variables suggesting a role of climate and geography in shaping the variation of this species. Fifty-seven SNPs from single trait GWAS and 11 SNPs from multitrait GWAS passed an FDR threshold of 0.05, leading to the identification of eight and seven nearby candidate genes, respectively. The percentage of phenotypic variance explained (PVE) by the significant SNPs for both single and multitrait GWAS ranged from 0.01% to 6.18%. To further evaluate the potential roles of candidate genes, we used a multi-omic network containing five additional data sets, including leaf and wood metabolite GWAS layers and coexpression and comethylation networks. We also performed a functional enrichment analysis on coexpression nearest neighbors for each gene model identified by the wood anatomical and morphological trait GWAS analyses. Genes affecting cell wall composition and transport related genes were enriched in wood anatomy and stomatal density trait networks. Signaling and metabolism related genes were also common in networks for stomatal density. For leaf morphology traits (leaf dry and wet weight) the networks were significantly enriched for GO terms related to photosynthetic processes as well as cellular homeostasis. The identified genes provide further insights into the genetic control of these traits, which are important determinants of the suitability and sustainability of improved genotypes for lignocellulosic biofuel production.
Collapse
Affiliation(s)
- Hari B. Chhetri
- Department of Biology, West Virginia University, Morgantown, WV, United States
| | - Anna Furches
- Biosciences Division, and The Center for Bioenergy Innovation, Oak Ridge National Laboratory, Oak Ridge, TN, United States
- The Bredesen Center for Interdisciplinary Research and Graduate Education, University of Tennessee, Knoxville, TN, United States
| | - David Macaya-Sanz
- Department of Biology, West Virginia University, Morgantown, WV, United States
| | - Alejandro R. Walker
- Department of Oral Biology, College of Dentistry, University of Florida, Gainesville, FL, United States
| | - David Kainer
- Biosciences Division, and The Center for Bioenergy Innovation, Oak Ridge National Laboratory, Oak Ridge, TN, United States
| | - Piet Jones
- Biosciences Division, and The Center for Bioenergy Innovation, Oak Ridge National Laboratory, Oak Ridge, TN, United States
- The Bredesen Center for Interdisciplinary Research and Graduate Education, University of Tennessee, Knoxville, TN, United States
| | - Anne E. Harman-Ware
- Biosciences Center, and National Bioenergy Center, National Renewable Energy Laboratory, Golden, CO, United States
| | - Timothy J. Tschaplinski
- Biosciences Division, and The Center for Bioenergy Innovation, Oak Ridge National Laboratory, Oak Ridge, TN, United States
| | - Daniel Jacobson
- Biosciences Division, and The Center for Bioenergy Innovation, Oak Ridge National Laboratory, Oak Ridge, TN, United States
- The Bredesen Center for Interdisciplinary Research and Graduate Education, University of Tennessee, Knoxville, TN, United States
| | - Gerald A. Tuskan
- Biosciences Division, and The Center for Bioenergy Innovation, Oak Ridge National Laboratory, Oak Ridge, TN, United States
| | - Stephen P. DiFazio
- Department of Biology, West Virginia University, Morgantown, WV, United States
- *Correspondence: Stephen P. DiFazio,
| |
Collapse
|
9
|
Chhetri HB, Macaya-Sanz D, Kainer D, Biswal AK, Evans LM, Chen JG, Collins C, Hunt K, Mohanty SS, Rosenstiel T, Ryno D, Winkeler K, Yang X, Jacobson D, Mohnen D, Muchero W, Strauss SH, Tschaplinski TJ, Tuskan GA, DiFazio SP. Multitrait genome-wide association analysis of Populus trichocarpa identifies key polymorphisms controlling morphological and physiological traits. New Phytol 2019; 223:293-309. [PMID: 30843213 DOI: 10.1111/nph.15777] [Citation(s) in RCA: 55] [Impact Index Per Article: 11.0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 12/07/2018] [Accepted: 02/22/2019] [Indexed: 05/08/2023]
Abstract
Genome-wide association studies (GWAS) have great promise for identifying the loci that contribute to adaptive variation, but the complex genetic architecture of many quantitative traits presents a substantial challenge. We measured 14 morphological and physiological traits and identified single nucleotide polymorphism (SNP)-phenotype associations in a Populus trichocarpa population distributed from California, USA to British Columbia, Canada. We used whole-genome resequencing data of 882 trees with more than 6.78 million SNPs, coupled with multitrait association to detect polymorphisms with potentially pleiotropic effects. Candidate genes were validated with functional data. Broad-sense heritability (H2 ) ranged from 0.30 to 0.56 for morphological traits and 0.08 to 0.36 for physiological traits. In total, 4 and 20 gene models were detected using the single-trait and multitrait association methods, respectively. Several of these associations were corroborated by additional lines of evidence, including co-expression networks, metabolite analyses, and direct confirmation of gene function through RNAi. Multitrait association identified many more significant associations than single-trait association, potentially revealing pleiotropic effects of individual genes. This approach can be particularly useful for challenging physiological traits such as water-use efficiency or complex traits such as leaf morphology, for which we were able to identify credible candidate genes by combining multitrait association with gene co-expression and co-methylation data.
Collapse
Affiliation(s)
- Hari B Chhetri
- Department of Biology, West Virginia University, Morgantown, WV, 26506, USA
| | - David Macaya-Sanz
- Department of Biology, West Virginia University, Morgantown, WV, 26506, USA
| | - David Kainer
- Biosciences Division, Oak Ridge National Laboratory, Oak Ridge, TN, 37831, USA
| | - Ajaya K Biswal
- Department of Biochemistry and Molecular Biology, University of Georgia, Athens, GA, 30602, USA
- Complex Carbohydrate Research Center, University of Georgia, Athens, GA, 30602, USA
| | - Luke M Evans
- Department of Biology, West Virginia University, Morgantown, WV, 26506, USA
| | - Jin-Gui Chen
- Biosciences Division, Oak Ridge National Laboratory, Oak Ridge, TN, 37831, USA
| | | | - Kimberly Hunt
- ArborGen, Inc., 2011 Broadbank Ct., Ridgeville, SC, 29472, USA
| | - Sushree S Mohanty
- Department of Biochemistry and Molecular Biology, University of Georgia, Athens, GA, 30602, USA
| | - Todd Rosenstiel
- Department of Biology, Portland State University, Portland, OR, 97207, USA
| | - David Ryno
- Complex Carbohydrate Research Center, University of Georgia, Athens, GA, 30602, USA
| | - Kim Winkeler
- ArborGen, Inc., 2011 Broadbank Ct., Ridgeville, SC, 29472, USA
| | - Xiaohan Yang
- Biosciences Division, Oak Ridge National Laboratory, Oak Ridge, TN, 37831, USA
| | - Daniel Jacobson
- Biosciences Division, Oak Ridge National Laboratory, Oak Ridge, TN, 37831, USA
| | - Debra Mohnen
- Department of Biochemistry and Molecular Biology, University of Georgia, Athens, GA, 30602, USA
- Complex Carbohydrate Research Center, University of Georgia, Athens, GA, 30602, USA
| | - Wellington Muchero
- Biosciences Division, Oak Ridge National Laboratory, Oak Ridge, TN, 37831, USA
| | - Steven H Strauss
- Department of Forest Ecosystems & Society, Oregon State University, Corvallis, OR, 97331, USA
| | | | - Gerald A Tuskan
- Biosciences Division, Oak Ridge National Laboratory, Oak Ridge, TN, 37831, USA
| | - Stephen P DiFazio
- Department of Biology, West Virginia University, Morgantown, WV, 26506, USA
| |
Collapse
|
10
|
Weighill D, Macaya-Sanz D, DiFazio SP, Joubert W, Shah M, Schmutz J, Sreedasyam A, Tuskan G, Jacobson D. Wavelet-Based Genomic Signal Processing for Centromere Identification and Hypothesis Generation. Front Genet 2019; 10:487. [PMID: 31214244 PMCID: PMC6554479 DOI: 10.3389/fgene.2019.00487] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/02/2018] [Accepted: 05/06/2019] [Indexed: 12/14/2022] Open
Abstract
Various ‘omics data types have been generated for Populus trichocarpa, each providing a layer of information which can be represented as a density signal across a chromosome. We make use of genome sequence data, variants data across a population as well as methylation data across 10 different tissues, combined with wavelet-based signal processing to perform a comprehensive analysis of the signature of the centromere in these different data signals, and successfully identify putative centromeric regions in P. trichocarpa from these signals. Furthermore, using SNP (single nucleotide polymorphism) correlations across a natural population of P. trichocarpa, we find evidence for the co-evolution of the centromeric histone CENH3 with the sequence of the newly identified centromeric regions, and identify a new CENH3 candidate in P. trichocarpa.
Collapse
Affiliation(s)
- Deborah Weighill
- The Bredesen Center for Interdisciplinary Research and Graduate Education, University of Tennessee, Knoxville, TN, United States.,Biosciences Division, Oak Ridge National Laboratory, Oak Ridge, TN, United States
| | - David Macaya-Sanz
- Department of Biology, West Virginia University, Morgantown, WV, United States
| | | | - Wayne Joubert
- Oak Ridge Leadership Computing Facility, Oak Ridge National Laboratory, Oak Ridge, TN, United States
| | - Manesh Shah
- Biosciences Division, Oak Ridge National Laboratory, Oak Ridge, TN, United States
| | - Jeremy Schmutz
- Department of Energy Joint Genome Institute, Walnut Creek, CA, United States.,HudsonAlpha Institute for Biotechnology, Huntsville, AL, United States
| | | | - Gerald Tuskan
- Biosciences Division, Oak Ridge National Laboratory, Oak Ridge, TN, United States
| | - Daniel Jacobson
- The Bredesen Center for Interdisciplinary Research and Graduate Education, University of Tennessee, Knoxville, TN, United States.,Biosciences Division, Oak Ridge National Laboratory, Oak Ridge, TN, United States
| |
Collapse
|
11
|
Weighill D, Jones P, Bleker C, Ranjan P, Shah M, Zhao N, Martin M, DiFazio S, Macaya-Sanz D, Schmutz J, Sreedasyam A, Tschaplinski T, Tuskan G, Jacobson D. Multi-Phenotype Association Decomposition: Unraveling Complex Gene-Phenotype Relationships. Front Genet 2019; 10:417. [PMID: 31134130 PMCID: PMC6522845 DOI: 10.3389/fgene.2019.00417] [Citation(s) in RCA: 14] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/19/2018] [Accepted: 04/16/2019] [Indexed: 01/18/2023] Open
Abstract
Various patterns of multi-phenotype associations (MPAs) exist in the results of Genome Wide Association Studies (GWAS) involving different topologies of single nucleotide polymorphism (SNP)-phenotype associations. These can provide interesting information about the different impacts of a gene on closely related phenotypes or disparate phenotypes (pleiotropy). In this work we present MPA Decomposition, a new network-based approach which decomposes the results of a multi-phenotype GWAS study into three bipartite networks, which, when used together, unravel the multi-phenotype signatures of genes on a genome-wide scale. The decomposition involves the construction of a phenotype powerset space, and subsequent mapping of genes into this new space. Clustering of genes in this powerset space groups genes based on their detailed MPA signatures. We show that this method allows us to find multiple different MPA and pleiotropic signatures within individual genes and to classify and cluster genes based on these SNP-phenotype association topologies. We demonstrate the use of this approach on a GWAS analysis of a large population of 882 Populus trichocarpa genotypes using untargeted metabolomics phenotypes. This method should prove invaluable in the interpretation of large GWAS datasets and aid in future synthetic biology efforts designed to optimize phenotypes of interest.
Collapse
Affiliation(s)
- Deborah Weighill
- The Bredesen Center for Interdisciplinary Research and Graduate Education, University of Tennessee, Knoxville, TN, United States.,Biosciences Division, Oak Ridge National Laboratory, Oak Ridge, TN, United States
| | - Piet Jones
- The Bredesen Center for Interdisciplinary Research and Graduate Education, University of Tennessee, Knoxville, TN, United States.,Biosciences Division, Oak Ridge National Laboratory, Oak Ridge, TN, United States
| | - Carissa Bleker
- The Bredesen Center for Interdisciplinary Research and Graduate Education, University of Tennessee, Knoxville, TN, United States.,Biosciences Division, Oak Ridge National Laboratory, Oak Ridge, TN, United States
| | - Priya Ranjan
- Biosciences Division, Oak Ridge National Laboratory, Oak Ridge, TN, United States.,Department of Plant Sciences, The University of Tennessee Institute of Agriculture, University of Tennessee, Knoxville, TN, United States
| | - Manesh Shah
- Biosciences Division, Oak Ridge National Laboratory, Oak Ridge, TN, United States
| | - Nan Zhao
- Department of Plant Sciences, The University of Tennessee Institute of Agriculture, University of Tennessee, Knoxville, TN, United States
| | - Madhavi Martin
- Biosciences Division, Oak Ridge National Laboratory, Oak Ridge, TN, United States
| | - Stephen DiFazio
- Department of Biology, West Virginia University, Morgantown, WV, United States
| | - David Macaya-Sanz
- Department of Biology, West Virginia University, Morgantown, WV, United States
| | - Jeremy Schmutz
- Department of Energy Joint Genome Institute, Walnut Creek, CA, United States.,HudsonAlpha Institute for Biotechnology, Huntsville, AL, United States
| | | | - Timothy Tschaplinski
- Biosciences Division, Oak Ridge National Laboratory, Oak Ridge, TN, United States
| | - Gerald Tuskan
- Biosciences Division, Oak Ridge National Laboratory, Oak Ridge, TN, United States
| | - Daniel Jacobson
- The Bredesen Center for Interdisciplinary Research and Graduate Education, University of Tennessee, Knoxville, TN, United States.,Biosciences Division, Oak Ridge National Laboratory, Oak Ridge, TN, United States
| |
Collapse
|
12
|
Furches A, Kainer D, Weighill D, Large A, Jones P, Walker AM, Romero J, Gazolla JGFM, Joubert W, Shah M, Streich J, Ranjan P, Schmutz J, Sreedasyam A, Macaya-Sanz D, Zhao N, Martin MZ, Rao X, Dixon RA, DiFazio S, Tschaplinski TJ, Chen JG, Tuskan GA, Jacobson D. Finding New Cell Wall Regulatory Genes in Populus trichocarpa Using Multiple Lines of Evidence. Front Plant Sci 2019; 10:1249. [PMID: 31649710 PMCID: PMC6791931 DOI: 10.3389/fpls.2019.01249] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 02/15/2019] [Accepted: 09/09/2019] [Indexed: 05/05/2023]
Abstract
Understanding the regulatory network controlling cell wall biosynthesis is of great interest in Populus trichocarpa, both because of its status as a model woody perennial and its importance for lignocellulosic products. We searched for genes with putatively unknown roles in regulating cell wall biosynthesis using an extended network-based Lines of Evidence (LOE) pipeline to combine multiple omics data sets in P. trichocarpa, including gene coexpression, gene comethylation, population level pairwise SNP correlations, and two distinct SNP-metabolite Genome Wide Association Study (GWAS) layers. By incorporating validation, ranking, and filtering approaches we produced a list of nine high priority gene candidates for involvement in the regulation of cell wall biosynthesis. We subsequently performed a detailed investigation of candidate gene GROWTH-REGULATING FACTOR 9 (PtGRF9). To investigate the role of PtGRF9 in regulating cell wall biosynthesis, we assessed the genome-wide connections of PtGRF9 and a paralog across data layers with functional enrichment analyses, predictive transcription factor binding site analysis, and an independent comparison to eQTN data. Our findings indicate that PtGRF9 likely affects the cell wall by directly repressing genes involved in cell wall biosynthesis, such as PtCCoAOMT and PtMYB.41, and indirectly by regulating homeobox genes. Furthermore, evidence suggests that PtGRF9 paralogs may act as transcriptional co-regulators that direct the global energy usage of the plant. Using our extended pipeline, we show multiple lines of evidence implicating the involvement of these genes in cell wall regulatory functions and demonstrate the value of this method for prioritizing candidate genes for experimental validation.
Collapse
Affiliation(s)
- Anna Furches
- Biosciences Division, and The Center for Bioenergy Innovation, Oak Ridge National Laboratory, Oak Ridge, TN, United States
- The Bredesen Center for Interdisciplinary Research and Graduate Education, University of Tennessee, Knoxville, TN, United States
| | - David Kainer
- Biosciences Division, and The Center for Bioenergy Innovation, Oak Ridge National Laboratory, Oak Ridge, TN, United States
| | - Deborah Weighill
- Biosciences Division, and The Center for Bioenergy Innovation, Oak Ridge National Laboratory, Oak Ridge, TN, United States
- The Bredesen Center for Interdisciplinary Research and Graduate Education, University of Tennessee, Knoxville, TN, United States
| | - Annabel Large
- Biosciences Division, and The Center for Bioenergy Innovation, Oak Ridge National Laboratory, Oak Ridge, TN, United States
- Oak Ridge Associated Universities (ORAU), Oak Ridge, TN, United States
- Department of Chemical and Biomolecular Engineering, University of Tennessee, Knoxville, TN, United States
| | - Piet Jones
- Biosciences Division, and The Center for Bioenergy Innovation, Oak Ridge National Laboratory, Oak Ridge, TN, United States
- The Bredesen Center for Interdisciplinary Research and Graduate Education, University of Tennessee, Knoxville, TN, United States
| | - Angelica M. Walker
- Biosciences Division, and The Center for Bioenergy Innovation, Oak Ridge National Laboratory, Oak Ridge, TN, United States
- Oak Ridge Associated Universities (ORAU), Oak Ridge, TN, United States
- Department of Computer Science, Johns Hopkins University, Baltimore, MD, United States
- Department of Biology, Johns Hopkins University, Baltimore, MD, United States
| | - Jonathon Romero
- Biosciences Division, and The Center for Bioenergy Innovation, Oak Ridge National Laboratory, Oak Ridge, TN, United States
- The Bredesen Center for Interdisciplinary Research and Graduate Education, University of Tennessee, Knoxville, TN, United States
| | | | - Wayne Joubert
- Oak Ridge Leadership Computing Facility, Oak Ridge National Laboratory, Oak Ridge, TN, United States
| | - Manesh Shah
- Biosciences Division, and The Center for Bioenergy Innovation, Oak Ridge National Laboratory, Oak Ridge, TN, United States
| | - Jared Streich
- Biosciences Division, and The Center for Bioenergy Innovation, Oak Ridge National Laboratory, Oak Ridge, TN, United States
| | - Priya Ranjan
- Biosciences Division, and The Center for Bioenergy Innovation, Oak Ridge National Laboratory, Oak Ridge, TN, United States
- Department of Plant Sciences, The University of Tennessee Institute of Agriculture, University of Tennessee, Knoxville, TN, United States
| | - Jeremy Schmutz
- Joint Genome Institute, Walnut Creek, CA, United States
- HudsonAlpha Institute for Biotechnology, Huntsville, AL, United States
| | | | - David Macaya-Sanz
- Department of Biology, West Virginia University, Morgantown, WV, United States
| | - Nan Zhao
- Department of Plant Sciences, The University of Tennessee Institute of Agriculture, University of Tennessee, Knoxville, TN, United States
| | - Madhavi Z. Martin
- Biosciences Division, and The Center for Bioenergy Innovation, Oak Ridge National Laboratory, Oak Ridge, TN, United States
| | - Xiaolan Rao
- BioDiscovery Institute and Department of Biological Sciences, University of North Texas, Denton, TX, United States
| | - Richard A. Dixon
- BioDiscovery Institute and Department of Biological Sciences, University of North Texas, Denton, TX, United States
| | - Stephen DiFazio
- Department of Biology, West Virginia University, Morgantown, WV, United States
| | - Timothy J. Tschaplinski
- Biosciences Division, and The Center for Bioenergy Innovation, Oak Ridge National Laboratory, Oak Ridge, TN, United States
| | - Jin-Gui Chen
- Biosciences Division, and The Center for Bioenergy Innovation, Oak Ridge National Laboratory, Oak Ridge, TN, United States
| | - Gerald A. Tuskan
- Biosciences Division, and The Center for Bioenergy Innovation, Oak Ridge National Laboratory, Oak Ridge, TN, United States
| | - Daniel Jacobson
- Biosciences Division, and The Center for Bioenergy Innovation, Oak Ridge National Laboratory, Oak Ridge, TN, United States
- The Bredesen Center for Interdisciplinary Research and Graduate Education, University of Tennessee, Knoxville, TN, United States
- *Correspondence: Daniel Jacobson,
| |
Collapse
|
13
|
Zhou R, Macaya-Sanz D, Rodgers-Melnick E, Carlson CH, Gouker FE, Evans LM, Schmutz J, Jenkins JW, Yan J, Tuskan GA, Smart LB, DiFazio SP. Characterization of a large sex determination region in Salix purpurea L. (Salicaceae). Mol Genet Genomics 2018; 293:1437-1452. [PMID: 30022352 DOI: 10.1007/s00438-018-1473-y] [Citation(s) in RCA: 47] [Impact Index Per Article: 7.8] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/16/2018] [Accepted: 07/09/2018] [Indexed: 12/30/2022]
Abstract
Dioecy has evolved numerous times in plants, but heteromorphic sex chromosomes are apparently rare. Sex determination has been studied in multiple Salix and Populus (Salicaceae) species, and P. trichocarpa has an XY sex determination system on chromosome 19, while S. suchowensis and S. viminalis have a ZW system on chromosome 15. Here we use whole genome sequencing coupled with quantitative trait locus mapping and a genome-wide association study to characterize the genomic composition of the non-recombining portion of the sex determination region. We demonstrate that Salix purpurea also has a ZW system on chromosome 15. The sex determination region has reduced recombination, high structural polymorphism, an abundance of transposable elements, and contains genes that are involved in sex expression in other plants. We also show that chromosome 19 contains sex-associated markers in this S. purpurea assembly, along with other autosomes. This raises the intriguing possibility of a translocation of the sex determination region within the Salicaceae lineage, suggesting a common evolutionary origin of the Populus and Salix sex determination loci.
Collapse
Affiliation(s)
- Ran Zhou
- Department of Biology, West Virginia University, 53 Campus Drive, Morgantown, WV, 26506-6057, USA
| | - David Macaya-Sanz
- Department of Biology, West Virginia University, 53 Campus Drive, Morgantown, WV, 26506-6057, USA
| | - Eli Rodgers-Melnick
- Department of Biology, West Virginia University, 53 Campus Drive, Morgantown, WV, 26506-6057, USA
| | - Craig H Carlson
- Horticulture Section, School of Integrative Plant Science, Cornell University, New York State Agricultural Experiment Station, Geneva, NY, 14456, USA
| | - Fred E Gouker
- Horticulture Section, School of Integrative Plant Science, Cornell University, New York State Agricultural Experiment Station, Geneva, NY, 14456, USA
| | - Luke M Evans
- Department of Biology, West Virginia University, 53 Campus Drive, Morgantown, WV, 26506-6057, USA
| | - Jeremy Schmutz
- HudsonAlpha Institute of Biotechnology, 601 Genome Way Northwest, Huntsville, AL, 35806, USA.,Department of Energy Joint Genome Institute, 2800 Mitchell Drive, Walnut Creek, CA, 94598, USA
| | - Jerry W Jenkins
- HudsonAlpha Institute of Biotechnology, 601 Genome Way Northwest, Huntsville, AL, 35806, USA
| | - Juying Yan
- Department of Energy Joint Genome Institute, 2800 Mitchell Drive, Walnut Creek, CA, 94598, USA
| | - Gerald A Tuskan
- Department of Energy Joint Genome Institute, 2800 Mitchell Drive, Walnut Creek, CA, 94598, USA.,Biosciences Division, Oak Ridge National Lab, Oak Ridge, USA
| | - Lawrence B Smart
- Horticulture Section, School of Integrative Plant Science, Cornell University, New York State Agricultural Experiment Station, Geneva, NY, 14456, USA
| | - Stephen P DiFazio
- Department of Biology, West Virginia University, 53 Campus Drive, Morgantown, WV, 26506-6057, USA.
| |
Collapse
|
14
|
Macaya-Sanz D, Chen J, Kalluri UC, Muchero W, Tschaplinski TJ, Gunter LE, Simon SJ, Biswal AK, Bryan AC, Payyavula R, Xie M, Yang Y, Zhang J, Mohnen D, Tuskan GA, DiFazio SP. Agronomic performance of Populus deltoides trees engineered for biofuel production. Biotechnol Biofuels 2017; 10:253. [PMID: 29213313 PMCID: PMC5707814 DOI: 10.1186/s13068-017-0934-6] [Citation(s) in RCA: 14] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 06/26/2017] [Accepted: 10/19/2017] [Indexed: 05/02/2023]
Abstract
BACKGROUND One of the major barriers to the development of lignocellulosic feedstocks is the recalcitrance of plant cell walls to deconstruction and saccharification. Recalcitrance can be reduced by targeting genes involved in cell wall biosynthesis, but this can have unintended consequences that compromise the agronomic performance of the trees under field conditions. Here we report the results of a field trial of fourteen distinct transgenic Populus deltoides lines that had previously demonstrated reduced recalcitrance without yield penalties under greenhouse conditions. RESULTS Survival and productivity of the trial were excellent in the first year, and there was little evidence for reduced performance of the transgenic lines with modified target gene expression. Surprisingly, the most striking phenotypic effects in this trial were for two empty-vector control lines that had modified bud set and bud flush. This is most likely due to somaclonal variation or insertional mutagenesis. Traits related to yield, crown architecture, herbivory, pathogen response, and frost damage showed few significant differences between target gene transgenics and empty vector controls. However, there were a few interesting exceptions. Lines overexpressing the DUF231 gene, a putative O-acetyltransferase, showed early bud flush and marginally increased height growth. Lines overexpressing the DUF266 gene, a putative glycosyltransferase, had significantly decreased stem internode length and slightly higher volume index. Finally, lines overexpressing the PFD2 gene, a putative member of the prefoldin complex, had a slightly reduced volume index. CONCLUSIONS This field trial demonstrates that these cell wall modifications, which decreased cell wall recalcitrance under laboratory conditions, did not seriously compromise first-year performance in the field, despite substantial challenges, including an outbreak of a stem boring insect (Gypsonoma haimbachiana), attack by a leaf rust pathogen (Melampsora spp.), and a late frost event. This bodes well for the potential utility of these lines as advanced biofuels feedstocks.
Collapse
Affiliation(s)
- David Macaya-Sanz
- Department of Biology, West Virginia University, Morgantown, WV 26506 USA
| | - Jin‐Gui Chen
- BioEnergy Science Center and Biosciences Division, Oak Ridge National Laboratory, Oak Ridge, TN 37831 USA
| | - Udaya C. Kalluri
- BioEnergy Science Center and Biosciences Division, Oak Ridge National Laboratory, Oak Ridge, TN 37831 USA
| | - Wellington Muchero
- BioEnergy Science Center and Biosciences Division, Oak Ridge National Laboratory, Oak Ridge, TN 37831 USA
| | - Timothy J. Tschaplinski
- BioEnergy Science Center and Biosciences Division, Oak Ridge National Laboratory, Oak Ridge, TN 37831 USA
| | - Lee E. Gunter
- BioEnergy Science Center and Biosciences Division, Oak Ridge National Laboratory, Oak Ridge, TN 37831 USA
| | - Sandra J. Simon
- Department of Biology, West Virginia University, Morgantown, WV 26506 USA
| | - Ajaya K. Biswal
- Department of Biochemistry and Molecular Biology, University of Georgia, Athens, GA 30602 USA
- Complex Carbohydrate Research Center, University of Georgia, Athens, GA 30602 USA
| | - Anthony C. Bryan
- BioEnergy Science Center and Biosciences Division, Oak Ridge National Laboratory, Oak Ridge, TN 37831 USA
| | - Raja Payyavula
- BioEnergy Science Center and Biosciences Division, Oak Ridge National Laboratory, Oak Ridge, TN 37831 USA
| | - Meng Xie
- BioEnergy Science Center and Biosciences Division, Oak Ridge National Laboratory, Oak Ridge, TN 37831 USA
| | - Yongil Yang
- BioEnergy Science Center and Biosciences Division, Oak Ridge National Laboratory, Oak Ridge, TN 37831 USA
| | - Jin Zhang
- BioEnergy Science Center and Biosciences Division, Oak Ridge National Laboratory, Oak Ridge, TN 37831 USA
| | - Debra Mohnen
- Department of Biochemistry and Molecular Biology, University of Georgia, Athens, GA 30602 USA
- Complex Carbohydrate Research Center, University of Georgia, Athens, GA 30602 USA
| | - Gerald A. Tuskan
- BioEnergy Science Center and Biosciences Division, Oak Ridge National Laboratory, Oak Ridge, TN 37831 USA
| | - Stephen P. DiFazio
- Department of Biology, West Virginia University, Morgantown, WV 26506 USA
| |
Collapse
|
15
|
Macaya-Sanz D, Heuertz M, Lindtke D, Vendramin GG, Lexer C, González-Martínez SC. Causes and consequences of large clonal assemblies in a poplar hybrid zone. Mol Ecol 2016; 25:5330-5344. [PMID: 27661461 DOI: 10.1111/mec.13850] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/13/2015] [Revised: 09/02/2016] [Accepted: 09/06/2016] [Indexed: 11/27/2022]
Abstract
Asexual reproduction is a common and fundamental mode of reproduction in plants. Although persistence in adverse conditions underlies most known cases of clonal dominance, proximal genetic drivers remain unclear, in particular for populations dominated by a few large clones. In this study, we studied a clonal population of the riparian tree Populus alba in the Douro river basin (northwestern Iberian Peninsula) where it hybridizes with Populus tremula, a species that grows in highly contrasted ecological conditions. We used 73 nuclear microsatellites to test whether genomic background (species ancestry) is a relevant cause of clonal success, and to assess the evolutionary consequences of clonal dominance by a few genets. Additional genotyping-by-sequencing data were produced to estimate the age of the largest clones. We found that a few ancient (over a few thousand years old) and widespread genets dominate the population, both in terms of clone size and number of sexual offspring produced. Interestingly, large clones possessed two genomic regions introgressed from P. tremula, which may have favoured their spread under stressful environmental conditions. At the population level, the spread of large genets was accompanied by an overall ancient (>0.1 Myr) but soft decline of effective population size. Despite this decrease, and the high clonality and dominance of sexual reproduction by large clones, the Douro hybrid zone still displays considerable genetic diversity and low inbreeding. This suggests that even in extreme cases as in the Douro, asexual and sexual dominance of a few large, geographically extended individuals does not threaten population survival.
Collapse
Affiliation(s)
- David Macaya-Sanz
- Department of Forest Ecology and Genetics, INIA-Forest Research Centre, Madrid, 28040, Spain.,Department of Biology, West Virginia University, Morgantown, WV, 26505, USA
| | | | - Dorothea Lindtke
- Unit of Ecology and Evolution, Department of Biology, University of Fribourg, Fribourg, 1700, Switzerland.,Department of Animal and Plant Sciences, University of Sheffield, Sheffield, S10 2TN, UK
| | - Giovanni G Vendramin
- Institute of Biosciences and Bioresources, National Research Council, Sesto Fiorentino (Florence), 50019, Italy
| | - Christian Lexer
- Unit of Ecology and Evolution, Department of Biology, University of Fribourg, Fribourg, 1700, Switzerland.,Department of Botany and Biodiversity Research, Faculty of Life Sciences, University of Vienna, Vienna, A-1030, Austria
| | - Santiago C González-Martínez
- Department of Forest Ecology and Genetics, INIA-Forest Research Centre, Madrid, 28040, Spain. .,BIOGECO, INRA, Univ. Bordeaux, Cestas, 33610, France.
| |
Collapse
|
16
|
Fillat Ú, Martín-Sampedro R, Macaya-Sanz D, Martín JA, Ibarra D, Martínez MJ, Eugenio ME. Screening of eucalyptus wood endophytes for laccase activity. Process Biochem 2016. [DOI: 10.1016/j.procbio.2016.02.006] [Citation(s) in RCA: 17] [Impact Index Per Article: 2.1] [Reference Citation Analysis] [What about the content of this article? (0)] [Track Full Text] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Indexed: 10/22/2022]
|
17
|
Blumenstein K, Macaya-Sanz D, Martín JA, Albrectsen BR, Witzell J. Phenotype MicroArrays as a complementary tool to next generation sequencing for characterization of tree endophytes. Front Microbiol 2015; 6:1033. [PMID: 26441951 PMCID: PMC4585013 DOI: 10.3389/fmicb.2015.01033] [Citation(s) in RCA: 22] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/17/2015] [Accepted: 09/11/2015] [Indexed: 12/31/2022] Open
Abstract
There is an increasing need to calibrate microbial community profiles obtained through next generation sequencing (NGS) with relevant taxonomic identities of the microbes, and to further associate these identities with phenotypic attributes. Phenotype MicroArray (PM) techniques provide a semi-high throughput assay for characterization and monitoring the microbial cellular phenotypes. Here, we present detailed descriptions of two different PM protocols used in our recent studies on fungal endophytes of forest trees, and highlight the benefits and limitations of this technique. We found that the PM approach enables effective screening of substrate utilization by endophytes. However, the technical limitations are multifaceted and the interpretation of the PM data challenging. For the best result, we recommend that the growth conditions for the fungi are carefully standardized. In addition, rigorous replication and control strategies should be employed whether using pre-configured, commercial microwell-plates or in-house designed PM plates for targeted substrate analyses. With these precautions, the PM technique is a valuable tool to characterize the metabolic capabilities of individual endophyte isolates, or successional endophyte communities identified by NGS, allowing a functional interpretation of the taxonomic data. Thus, PM approaches can provide valuable complementary information for NGS studies of fungal endophytes in forest trees.
Collapse
Affiliation(s)
- Kathrin Blumenstein
- Southern Swedish Forest Research Centre, Swedish University of Agricultural Sciences, AlnarpSweden
| | - David Macaya-Sanz
- Department of Natural Systems and Resources, School of Forest Engineers, Technical University of MadridMadrid, Spain
| | - Juan A. Martín
- Department of Natural Systems and Resources, School of Forest Engineers, Technical University of MadridMadrid, Spain
| | - Benedicte R. Albrectsen
- Department of Plant Physiology, Umeå Plant Science Centre, Umeå UniversityUmeå, Sweden
- Department of Plant and Environmental Sciences, University of CopenhagenCopenhagen, Denmark
| | - Johanna Witzell
- Southern Swedish Forest Research Centre, Swedish University of Agricultural Sciences, AlnarpSweden
- School of Forest Sciences, Faculty of Science and Forestry, University of Eastern Finland, JoensuuFinland
| |
Collapse
|
18
|
Lindtke D, González-Martínez SC, Macaya-Sanz D, Lexer C. Admixture mapping of quantitative traits in Populus hybrid zones: power and limitations. Heredity (Edinb) 2013; 111:474-85. [PMID: 23860234 PMCID: PMC3833683 DOI: 10.1038/hdy.2013.69] [Citation(s) in RCA: 32] [Impact Index Per Article: 2.9] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/13/2013] [Revised: 06/06/2013] [Accepted: 06/12/2013] [Indexed: 01/30/2023] Open
Abstract
Uncovering the genetic architecture of species differences is of central importance for understanding the origin and maintenance of biological diversity. Admixture mapping can be used to identify the number and effect sizes of genes that contribute to the divergence of ecologically important traits, even in taxa that are not amenable to laboratory crosses because of their long generation time or other limitations. Here, we apply admixture mapping to naturally occurring hybrids between two ecologically divergent Populus species. We map quantitative trait loci for eight leaf morphological traits using 77 mapped microsatellite markers from all 19 chromosomes of Populus. We apply multivariate linear regression analysis allowing the modeling of additive and non-additive gene action and identify several candidate genomic regions associated with leaf morphology using an information-theoretic approach. We perform simulation studies to assess the power and limitations of admixture mapping of quantitative traits in natural hybrid populations for a variety of genetic architectures and modes of gene action. Our results indicate that (1) admixture mapping has considerable power to identify the genetic architecture of species differences if sample sizes and marker densities are sufficiently high, (2) modeling of non-additive gene action can help to elucidate the discrepancy between genotype and phenotype sometimes seen in interspecific hybrids, and (3) the genetic architecture of leaf morphological traits in the studied Populus species involves complementary and overdominant gene action, providing the basis for rapid adaptation of these ecologically important forest trees.
Collapse
Affiliation(s)
- D Lindtke
- Unit of Ecology and Evolution, Department of
Biology, University of Fribourg, Fribourg,
Switzerland
| | - S C González-Martínez
- INIA, Forest Research Centre, Department of
Forest Ecology and Genetics, Madrid, Spain
| | - D Macaya-Sanz
- INIA, Forest Research Centre, Department of
Forest Ecology and Genetics, Madrid, Spain
- Technical University of Madrid, ETS Forestry
Engineering, Department of Silviculture, Madrid, Spain
| | - C Lexer
- Unit of Ecology and Evolution, Department of
Biology, University of Fribourg, Fribourg,
Switzerland
- Jodrell Laboratory, Royal Botanic
Gardens, Kew, Richmond, Surrey, UK
| |
Collapse
|
19
|
Macaya-Sanz D, Suter L, Joseph J, Barbará T, Alba N, González-Martínez SC, Widmer A, Lexer C. Genetic analysis of post-mating reproductive barriers in hybridizing European Populus species. Heredity (Edinb) 2011; 107:478-86. [PMID: 21587301 DOI: 10.1038/hdy.2011.35] [Citation(s) in RCA: 25] [Impact Index Per Article: 1.9] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/09/2022] Open
Abstract
Molecular genetic analyses of experimental crosses provide important information on the strength and nature of post-mating barriers to gene exchange between divergent populations, which are topics of great interest to evolutionary geneticists and breeders. Although not a trivial task in long-lived organisms such as trees, experimental interspecific recombinants can sometimes be created through controlled crosses involving natural F(1)'s. Here, we used this approach to understand the genetics of post-mating isolation and barriers to introgression in Populus alba and Populus tremula, two ecologically divergent, hybridizing forest trees. We studied 86 interspecific backcross (BC(1)) progeny and >350 individuals from natural populations of these species for up to 98 nuclear genetic markers, including microsatellites, indels and single nucleotide polymorphisms, and inferred the origin of the cytoplasm of the cross with plastid DNA. Genetic analysis of the BC(1) revealed extensive segregation distortions on six chromosomes, and >90% of these (12 out of 13) favored P. tremula donor alleles in the heterospecific genomic background. Since selection was documented during early diploid stages of the progeny, this surprising result was attributed to epistasis, cyto-nuclear coadaptation, heterozygote advantage at nuclear loci experiencing introgression or a combination of these. Our results indicate that gene flow across 'porous' species barriers affects these poplars and aspens beyond neutral, Mendelian expectations and suggests the mechanisms responsible. Contrary to expectations, the Populus sex determination region is not protected from introgression. Understanding the population dynamics of the Populus sex determination region will require tests based on natural interspecific hybrid zones.
Collapse
Affiliation(s)
- D Macaya-Sanz
- Department of Forest Ecology and Genetics, Center of Forest Research, CIFOR-INIA, Carretera de A Coruña, Madrid, Spain
| | | | | | | | | | | | | | | |
Collapse
|