1
|
Darnet E, Teixeira B, Schaller H, Rogez H, Darnet S. Elucidating the Mesocarp Drupe Transcriptome of Açai ( Euterpe oleracea Mart.): An Amazonian Tree Palm Producer of Bioactive Compounds. Int J Mol Sci 2023; 24:ijms24119315. [PMID: 37298279 DOI: 10.3390/ijms24119315] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/29/2023] [Revised: 05/13/2023] [Accepted: 05/16/2023] [Indexed: 06/12/2023] Open
Abstract
Euterpe oleracea palm, endemic to the Amazon region, is well known for açai, a fruit violet beverage with nutritional and medicinal properties. During E. oleracea fruit ripening, anthocyanin accumulation is not related to sugar production, contrarily to grape and blueberry. Ripened fruits have a high content of anthocyanins, isoprenoids, fibers, and proteins, and are poor in sugars. E. oleracea is proposed as a new genetic model for metabolism partitioning in the fruit. Approximately 255 million single-end-oriented reads were generated on an Ion Proton NGS platform combining fruit cDNA libraries at four ripening stages. The de novo transcriptome assembly was tested using six assemblers and 46 different combinations of parameters, a pre-processing and a post-processing step. The multiple k-mer approach with TransABySS as an assembler and Evidential Gene as a post-processer have shown the best results, with an N50 of 959 bp, a read coverage mean of 70x, a BUSCO complete sequence recovery of 36% and an RBMT of 61%. The fruit transcriptome dataset included 22,486 transcripts representing 18 Mbp, of which a proportion of 87% had significant homology with other plant sequences. Approximately 904 new EST-SSRs were described, and were common and transferable to Phoenix dactylifera and Elaeis guineensis, two other palm trees. The global GO classification of transcripts showed similar categories to that in P. dactylifera and E. guineensis fruit transcriptomes. For an accurate annotation and functional description of metabolism genes, a bioinformatic pipeline was developed to precisely identify orthologs, such as one-to-one orthologs between species, and to infer multigenic family evolution. The phylogenetic inference confirmed an occurrence of duplication events in the Arecaceae lineage and the presence of orphan genes in E. oleracea. Anthocyanin and tocopherol pathways were annotated entirely. Interestingly, the anthocyanin pathway showed a high number of paralogs, similar to in grape, whereas the tocopherol pathway exhibited a low and conserved gene number and the prediction of several splicing forms. The release of this exhaustively annotated molecular dataset of E. oleracea constitutes a valuable tool for further studies in metabolism partitioning and opens new great perspectives to study fruit physiology with açai as a model.
Collapse
Affiliation(s)
- Elaine Darnet
- Centre for Valorization of Amazonian Bioactive Compounds (CVACBA) & Institute of Biological Sciences, Federal University of Pará (UFPA), Belém 66075-750, PA, Brazil
- International Associated Laboratory PALMHEAT, Frech Scientific Research National Center (CNRS)/UFPA, 75016 Paris, France
| | - Bruno Teixeira
- Centre for Valorization of Amazonian Bioactive Compounds (CVACBA) & Institute of Biological Sciences, Federal University of Pará (UFPA), Belém 66075-750, PA, Brazil
| | - Hubert Schaller
- International Associated Laboratory PALMHEAT, Frech Scientific Research National Center (CNRS)/UFPA, 75016 Paris, France
- Plant Isoprenoid Biology, Institute of Molecular Biology of Plants of the Scientific Research National Center, Strasbourg University, 67081 Strasbourg, France
| | - Hervé Rogez
- Centre for Valorization of Amazonian Bioactive Compounds (CVACBA) & Institute of Biological Sciences, Federal University of Pará (UFPA), Belém 66075-750, PA, Brazil
| | - Sylvain Darnet
- Centre for Valorization of Amazonian Bioactive Compounds (CVACBA) & Institute of Biological Sciences, Federal University of Pará (UFPA), Belém 66075-750, PA, Brazil
- International Associated Laboratory PALMHEAT, Frech Scientific Research National Center (CNRS)/UFPA, 75016 Paris, France
- Plant Isoprenoid Biology, Institute of Molecular Biology of Plants of the Scientific Research National Center, Strasbourg University, 67081 Strasbourg, France
| |
Collapse
|
2
|
Ashwood LM, Elnahriry KA, Stewart ZK, Shafee T, Naseem MU, Szanto TG, van der Burg CA, Smith HL, Surm JM, Undheim EAB, Madio B, Hamilton BR, Guo S, Wai DCC, Coyne VL, Phillips MJ, Dudley KJ, Hurwood DA, Panyi G, King GF, Pavasovic A, Norton RS, Prentis PJ. Genomic, functional and structural analyses elucidate evolutionary innovation within the sea anemone 8 toxin family. BMC Biol 2023; 21:121. [PMID: 37226201 DOI: 10.1186/s12915-023-01617-y] [Citation(s) in RCA: 5] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/08/2023] [Accepted: 05/09/2023] [Indexed: 05/26/2023] Open
Abstract
BACKGROUND The ShK toxin from Stichodactyla helianthus has established the therapeutic potential of sea anemone venom peptides, but many lineage-specific toxin families in Actiniarians remain uncharacterised. One such peptide family, sea anemone 8 (SA8), is present in all five sea anemone superfamilies. We explored the genomic arrangement and evolution of the SA8 gene family in Actinia tenebrosa and Telmatactis stephensoni, characterised the expression patterns of SA8 sequences, and examined the structure and function of SA8 from the venom of T. stephensoni. RESULTS We identified ten SA8-family genes in two clusters and six SA8-family genes in five clusters for T. stephensoni and A. tenebrosa, respectively. Nine SA8 T. stephensoni genes were found in a single cluster, and an SA8 peptide encoded by an inverted SA8 gene from this cluster was recruited to venom. We show that SA8 genes in both species are expressed in a tissue-specific manner and the inverted SA8 gene has a unique tissue distribution. While the functional activity of the SA8 putative toxin encoded by the inverted gene was inconclusive, its tissue localisation is similar to toxins used for predator deterrence. We demonstrate that, although mature SA8 putative toxins have similar cysteine spacing to ShK, SA8 peptides are distinct from ShK peptides based on structure and disulfide connectivity. CONCLUSIONS Our results provide the first demonstration that SA8 is a unique gene family in Actiniarians, evolving through a variety of structural changes including tandem and proximal gene duplication and an inversion event that together allowed SA8 to be recruited into the venom of T. stephensoni.
Collapse
Affiliation(s)
- Lauren M Ashwood
- School of Biology and Environmental Science, Faculty of Science, Queensland University of Technology, Brisbane, QLD, 4000, Australia.
- Cancer Program, QIMR Berghofer Medical Research Institute, Brisbane, QLD, 4006, Australia.
| | - Khaled A Elnahriry
- Medicinal Chemistry, Monash Institute of Pharmaceutical Sciences, Monash University, Parkville, VIC, 3052, Australia
| | - Zachary K Stewart
- Centre for Agriculture and the Bioeconomy, Queensland University of Technology, Brisbane, QLD, 4000, Australia
| | - Thomas Shafee
- Department of Animal Plant & Soil Sciences, La Trobe University, Melbourne, Australia
- Swinburne University of Technology, Melbourne, VIC, Australia
| | - Muhammad Umair Naseem
- Department of Biophysics and Cell Biology, Faculty of Medicine, University of Debrecen, 4032, Debrecen, Hungary
| | - Tibor G Szanto
- Department of Biophysics and Cell Biology, Faculty of Medicine, University of Debrecen, 4032, Debrecen, Hungary
| | - Chloé A van der Burg
- School of Biology and Environmental Science, Faculty of Science, Queensland University of Technology, Brisbane, QLD, 4000, Australia
- Department of Anatomy, School of Biomedical Sciences, University of Otago, Dunedin, 9016, New Zealand
| | - Hayden L Smith
- School of Biology and Environmental Science, Faculty of Science, Queensland University of Technology, Brisbane, QLD, 4000, Australia
| | - Joachim M Surm
- Department of Ecology, Evolution and Behavior, Alexander Silberman Institute of Life Sciences, The Hebrew University of Jerusalem, 9190401, Jerusalem, Israel
| | - Eivind A B Undheim
- Department of Biosciences, Centre for Ecological and Evolutionary Synthesis, University of Oslo, Blindern, PO Box 1066, 0316, Oslo, Norway
- Centre for Advanced Imaging, The University of Queensland, St Lucia, QLD, 4072, Australia
| | - Bruno Madio
- Institute for Molecular Bioscience, The University of Queensland, St Lucia, QLD, 4072, Australia
| | - Brett R Hamilton
- Centre for Advanced Imaging, The University of Queensland, St Lucia, QLD, 4072, Australia
- Centre for Microscopy and Microanalysis, The University of Queensland, St Lucia, QLD, 4072, Australia
| | - Shaodong Guo
- Institute for Molecular Bioscience, The University of Queensland, St Lucia, QLD, 4072, Australia
| | - Dorothy C C Wai
- Medicinal Chemistry, Monash Institute of Pharmaceutical Sciences, Monash University, Parkville, VIC, 3052, Australia
| | - Victoria L Coyne
- Research Infrastructure, Central Analytical Research Facility, Queensland University of Technology, Brisbane, QLD, 4000, Australia
| | - Matthew J Phillips
- School of Biology and Environmental Science, Faculty of Science, Queensland University of Technology, Brisbane, QLD, 4000, Australia
| | - Kevin J Dudley
- School of Biology and Environmental Science, Faculty of Science, Queensland University of Technology, Brisbane, QLD, 4000, Australia
- Research Infrastructure, Central Analytical Research Facility, Queensland University of Technology, Brisbane, QLD, 4000, Australia
| | - David A Hurwood
- School of Biology and Environmental Science, Faculty of Science, Queensland University of Technology, Brisbane, QLD, 4000, Australia
- Centre for Agriculture and the Bioeconomy, Queensland University of Technology, Brisbane, QLD, 4000, Australia
| | - Gyorgy Panyi
- Department of Biophysics and Cell Biology, Faculty of Medicine, University of Debrecen, 4032, Debrecen, Hungary
| | - Glenn F King
- Institute for Molecular Bioscience, The University of Queensland, St Lucia, QLD, 4072, Australia
- ARC Centre for Innovations in Peptide and Protein Science, The University of Queensland, St Lucia, QLD, 4072, Australia
| | - Ana Pavasovic
- School of Biomedical Sciences, Faculty of Health, Queensland University of Technology, Brisbane, QLD, 4000, Australia
| | - Raymond S Norton
- Medicinal Chemistry, Monash Institute of Pharmaceutical Sciences, Monash University, Parkville, VIC, 3052, Australia
- ARC Centre for Fragment-Based Design, Monash University, Parkville, VIC, 3052, Australia
| | - Peter J Prentis
- School of Biology and Environmental Science, Faculty of Science, Queensland University of Technology, Brisbane, QLD, 4000, Australia
- Centre for Agriculture and the Bioeconomy, Queensland University of Technology, Brisbane, QLD, 4000, Australia
| |
Collapse
|
3
|
Visser EA, Kampmann TP, Wegrzyn JL, Naidoo S. Multispecies comparison of host responses to Fusarium circinatum challenge in tropical pines show consistency in resistance mechanisms. PLANT, CELL & ENVIRONMENT 2023; 46:1705-1725. [PMID: 36541367 DOI: 10.1111/pce.14522] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 01/27/2022] [Accepted: 12/18/2022] [Indexed: 06/17/2023]
Abstract
Fusarium circinatum poses a threat to both commercial and natural pine forests. Large variation in host resistance exists between species, with many economically important species being susceptible. Development of resistant genotypes could be expedited and optimised by investigating the molecular mechanisms underlying host resistance and susceptibility as well as increasing the available genetic resources. RNA-seq data, from F. circinatum inoculated and mock-inoculated ca. 6-month-old shoot tissue at 3- and 7-days postinoculation, was generated for three commercially important tropical pines, Pinus oocarpa, Pinus maximinoi and Pinus greggii. De novo transcriptomes were assembled and used to investigate the NLR and PR gene content within available pine references. Host responses to F. circinatum challenge were investigated in P. oocarpa (resistant) and P. greggii (susceptible), in comparison to previously generated expression profiles from Pinus tecunumanii (resistant) and Pinus patula (susceptible). Expression results indicated crosstalk between induced salicylate, jasmonate and ethylene signalling is involved in host resistance and compromised in susceptible hosts. Additionally, higher constitutive expression of sulfur metabolism and flavonoid biosynthesis in resistant hosts suggest involvement of these metabolites in resistance.
Collapse
Affiliation(s)
- Erik A Visser
- Department of Biochemistry, Genetics and Microbiology, Forestry and Agricultural Biotechnology Institute (FABI), University of Pretoria, Pretoria, South Africa
| | - Tamanique P Kampmann
- Department of Biochemistry, Genetics and Microbiology, Forestry and Agricultural Biotechnology Institute (FABI), University of Pretoria, Pretoria, South Africa
| | - Jill L Wegrzyn
- Department of Ecology and Evolutionary Biology, University of Connecticut, Storrs, Connecticut, USA
| | - Sanushka Naidoo
- Department of Biochemistry, Genetics and Microbiology, Forestry and Agricultural Biotechnology Institute (FABI), University of Pretoria, Pretoria, South Africa
| |
Collapse
|
4
|
Alkaloid production and response to natural adverse conditions in Peganum harmala: in silico transcriptome analyses. BIOTECHNOLOGIA 2022; 103:355-384. [PMID: 36685700 PMCID: PMC9837557 DOI: 10.5114/bta.2022.120706] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/22/2021] [Revised: 07/25/2022] [Accepted: 09/16/2022] [Indexed: 01/06/2023] Open
Abstract
Peganum harmala is a valuable wild plant that grows and survives under adverse conditions and produces pharmaceutical alkaloid metabolites. Using different assemblers to develop a transcriptome improves the quality of assembled transcriptome. In this study, a concrete and accurate method for detecting stress-responsive transcripts by comparing stress-related gene ontology (GO) terms and public domains was designed. An integrated transcriptome for P. harmala including 42 656 coding sequences was created by merging de novo assembled transcriptomes. Around 35 000 transcripts were annotated with more than 90% resemblance to three closely related species of Citrus, which confirmed the robustness of the assembled transcriptome; 4853 stress-responsive transcripts were identified. CYP82 involved in alkaloid biosynthesis showed a higher number of transcripts in P. harmala than in other plants, indicating its diverse alkaloid biosynthesis attributes. Transcription factors (TFs) and regulatory elements with 3887 transcripts comprised 9% of the transcriptome. Among the TFs of the integrated transcriptome, cystein2/histidine2 (C2H2) and WD40 repeat families were the most abundant. The Kyoto Encyclopedia of Genes and Genomes (KEGG) MAPK (mitogen-activated protein kinase) signaling map and the plant hormone signal transduction map showed the highest assigned genes to these pathways, suggesting their potential stress resistance. The P. harmala whole-transcriptome survey provides important resources and paves the way for functional and comparative genomic studies on this plant to discover stress-tolerance-related markers and response mechanisms in stress physiology, phytochemistry, ecology, biodiversity, and evolution. P. harmala can be a potential model for studying adverse environmental cues and metabolite biosynthesis and a major source for the production of various alkaloids.
Collapse
|
5
|
Sheikh-Assadi M, Naderi R, Salami SA, Kafi M, Fatahi R, Shariati V, Martinelli F, Cicatelli A, Triassi M, Guarino F, Improta G, Claros MG. Normalized Workflow to Optimize Hybrid De Novo Transcriptome Assembly for Non-Model Species: A Case Study in Lilium ledebourii (Baker) Boiss. PLANTS 2022; 11:plants11182365. [PMID: 36145766 PMCID: PMC9503428 DOI: 10.3390/plants11182365] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 07/23/2022] [Revised: 08/21/2022] [Accepted: 09/07/2022] [Indexed: 11/16/2022]
Abstract
A high-quality transcriptome is required to advance numerous bioinformatics workflows. Nevertheless, the effectuality of tools for de novo assembly and real precision assembled transcriptomes looks somewhat unexplored, particularly for non-model organisms with complicated (very long, heterozygous, polyploid) genomes. To disclose the performance of various transcriptome assembly programs, this study built 11 single assemblies and analyzed their performance on some significant reference-free and reference-based criteria. As well as to reconfirm the outputs of benchmarks, 55 BLAST were performed and compared using 11 constructed transcriptomes. Concisely, normalized benchmarking demonstrated that Velvet–Oases suffer from the worst results, while the EvidentialGene strategy can provide the most comprehensive and accurate transcriptome of Lilium ledebourii (Baker) Boiss. The BLAST results also confirmed the superiority of EvidentialGene, so it could capture even up to 59% more (than Velvet–Oases) unique gene hits. To promote assembly optimization, with the help of normalized benchmarking, PCA and AHC, it is emphasized that each metric can only provide part of the transcriptome status, and one should never settle for just a few evaluation criteria. This study supplies a framework for benchmarking and optimizing the efficiency of assembly approaches to analyze RNA-Seq data and reveals that selecting an inefficient assembly strategy might result in less identification of unique gene hits.
Collapse
Affiliation(s)
- Morteza Sheikh-Assadi
- Department of Horticultural Science, Faculty of Agricultural Science and Engineering, University of Tehran, Karaj 31587-77871, Iran
- Correspondence: (M.S.-A.); (R.N.)
| | - Roohangiz Naderi
- Department of Horticultural Science, Faculty of Agricultural Science and Engineering, University of Tehran, Karaj 31587-77871, Iran
- Correspondence: (M.S.-A.); (R.N.)
| | - Seyed Alireza Salami
- Department of Horticultural Science, Faculty of Agricultural Science and Engineering, University of Tehran, Karaj 31587-77871, Iran
| | - Mohsen Kafi
- Department of Horticultural Science, Faculty of Agricultural Science and Engineering, University of Tehran, Karaj 31587-77871, Iran
| | - Reza Fatahi
- Department of Horticultural Science, Faculty of Agricultural Science and Engineering, University of Tehran, Karaj 31587-77871, Iran
| | - Vahid Shariati
- NIGEB Genome Center, National Institute of Genetic Engineering and Biotechnology, Tehran 14965/161, Iran
| | - Federico Martinelli
- Department of Biology, University of Florence, 50019 Sesto Fiorentino, Italy
| | - Angela Cicatelli
- Department of Chemistry and Biology “A. Zambelli”, University of Salerno, 84084 Fisciano, Italy
| | - Maria Triassi
- Department of Public Health, University of Naples “Federico II”, 80131 Naples, Italy
| | - Francesco Guarino
- Department of Chemistry and Biology “A. Zambelli”, University of Salerno, 84084 Fisciano, Italy
| | - Giovanni Improta
- Department of Public Health, University of Naples “Federico II”, 80131 Naples, Italy
| | - Manuel Gonzalo Claros
- Molecular Biology and Biochemistry Department, University of Málaga, 29071 Málaga, Spain
- CIBER de Enfermedades Raras (CIBERER), 29071 Málaga, Spain
- Institute of Biomedical Research in Málaga (IBIMA), IBIMA-RARE, 29010 Málaga, Spain
- Instituto de Hortofruticultura Subtropical y Mediterránea (IHSM-UMA-CSIC), 29010 Málaga, Spain
| |
Collapse
|
6
|
Zamora-Ballesteros C, Martín-García J, Suárez-Vega A, Diez JJ. Genome-wide identification and characterization of Fusarium circinatum-responsive lncRNAs in Pinus radiata. BMC Genomics 2022; 23:194. [PMID: 35264109 PMCID: PMC8908662 DOI: 10.1186/s12864-022-08408-9] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/09/2021] [Accepted: 02/14/2022] [Indexed: 12/14/2022] Open
Abstract
Background One of the most promising strategies of Pine Pitch Canker (PPC) management is the use of reproductive plant material resistant to the disease. Understanding the complexity of plant transcriptome that underlies the defence to the causal agent Fusarium circinatum, would greatly facilitate the development of an accurate breeding program. Long non-coding RNAs (lncRNAs) are emerging as important transcriptional regulators under biotic stresses in plants. However, to date, characterization of lncRNAs in conifer trees has not been reported. In this study, transcriptomic identification of lncRNAs was carried out using strand-specific paired-end RNA sequencing, from Pinus radiata samples inoculated with F. circinatum at an early stage of infection. Results Overall, 13,312 lncRNAs were predicted through a bioinformatics approach, including long intergenic non-coding RNAs (92.3%), antisense lncRNAs (3.3%) and intronic lncRNAs (2.9%). Compared with protein-coding RNAs, pine lncRNAs are shorter, have lower expression, lower GC content and harbour fewer and shorter exons. A total of 164 differentially expressed (DE) lncRNAs were identified in response to F. circinatum infection in the inoculated versus mock-inoculated P. radiata seedlings. The predicted cis-regulated target genes of these pathogen-responsive lncRNAs were related to defence mechanisms such as kinase activity, phytohormone regulation, and cell wall reinforcement. Co-expression network analysis of DE lncRNAs, DE protein-coding RNAs and lncRNA target genes also indicated a potential network regulating pectinesterase activity and cell wall remodelling. Conclusions This study presents the first comprehensive genome-wide analysis of P. radiata lncRNAs and provides the basis for future functional characterizations of lncRNAs in relation to pine defence responses against F. circinatum. Supplementary Information The online version contains supplementary material available at 10.1186/s12864-022-08408-9.
Collapse
Affiliation(s)
- Cristina Zamora-Ballesteros
- Department of Vegetal Production and Forest Resources, University of Valladolid, Av Madrid 44, 34004 Palencia, Spain. .,Sustainable Forest Management Research Institute, University of Valladolid-INIA, 34004 Palencia, Spain.
| | - Jorge Martín-García
- Department of Vegetal Production and Forest Resources, University of Valladolid, Av Madrid 44, 34004 Palencia, Spain.,Sustainable Forest Management Research Institute, University of Valladolid-INIA, 34004 Palencia, Spain
| | - Aroa Suárez-Vega
- Departamento de Producción Animal, Facultad de Veterinaria, University of León, Campus de Vegazana s/n, 24071 León, Spain
| | - Julio Javier Diez
- Department of Vegetal Production and Forest Resources, University of Valladolid, Av Madrid 44, 34004 Palencia, Spain.,Sustainable Forest Management Research Institute, University of Valladolid-INIA, 34004 Palencia, Spain
| |
Collapse
|
7
|
Amaral J, Valledor L, Alves A, Martín-García J, Pinto G. Studying tree response to biotic stress using a multi-disciplinary approach: The pine pitch canker case study. FRONTIERS IN PLANT SCIENCE 2022; 13:916138. [PMID: 36160962 PMCID: PMC9501998 DOI: 10.3389/fpls.2022.916138] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 04/08/2022] [Accepted: 08/18/2022] [Indexed: 05/09/2023]
Abstract
In an era of climate change and global trade, forests sustainability is endangered by several biotic threats. Pine pitch canker (PPC), caused by Fusarium circinatum, is one of the most important disease affecting conifers worldwide. To date, no effective control measures have been found for this disease. Earlier studies on PPC were mainly focused on the pathogen itself or on determining the levels of susceptibility of different hosts to F. circinatum infection. However, over the last years, plenty of information on the mechanisms that may explain the susceptibility or resistance to PPC has been published. This data are useful to better understand tree response to biotic stress and, most importantly, to aid the development of innovative and scientific-based disease control measures. This review gathers and discusses the main advances on PPC knowledge, especially focusing on multi-disciplinary studies investigating the response of pines with different levels of susceptibility to PPC upon infection. After an overview of the general knowledge of the disease, the importance of integrating information from physiological and Omics studies to unveil the mechanisms behind PPC susceptibility/resistance and to develop control strategies is explored. An extensive review of the main host responses to PPC was performed, including changes in water relations, signalling (ROS and hormones), primary metabolism, and defence (resin, phenolics, and PR proteins). A general picture of pine response to PPC is suggested according to the host susceptibility level and the next steps and gaps on PPC research are pointed out.
Collapse
Affiliation(s)
- Joana Amaral
- Centre for Environmental and Marine Studies (CESAM), Department of Biology, University of Aveiro, Aveiro, Portugal
- *Correspondence: Joana Amaral,
| | - Luis Valledor
- Department of Organisms and Systems Biology, University of Oviedo, Oviedo, Spain
- University Institute of Biotechnology of Asturias, University of Oviedo, Oviedo, Spain
| | - Artur Alves
- Centre for Environmental and Marine Studies (CESAM), Department of Biology, University of Aveiro, Aveiro, Portugal
| | - Jorge Martín-García
- Department of Vegetal Production and Forest Resources, University of Valladolid, Palencia, Spain
- Sustainable Forest Management Research Institute, University of Valladolid-INIA, Palencia, Spain
| | - Glória Pinto
- Centre for Environmental and Marine Studies (CESAM), Department of Biology, University of Aveiro, Aveiro, Portugal
- Glória Pinto,
| |
Collapse
|
8
|
Zamora-Ballesteros C, Pinto G, Amaral J, Valledor L, Alves A, Diez JJ, Martín-García J. Dual RNA-Sequencing Analysis of Resistant ( Pinus pinea) and Susceptible ( Pinus radiata) Hosts during Fusarium circinatum Challenge. Int J Mol Sci 2021; 22:5231. [PMID: 34063405 PMCID: PMC8156185 DOI: 10.3390/ijms22105231] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/23/2021] [Revised: 05/11/2021] [Accepted: 05/13/2021] [Indexed: 12/13/2022] Open
Abstract
Fusarium circinatum causes one of the most important diseases of conifers worldwide, the pine pitch canker (PPC). However, no effective field intervention measures aiming to control or eradicate PPC are available. Due to the variation in host genetic resistance, the development of resistant varieties is postulated as a viable and promising strategy. By using an integrated approach, this study aimed to identify differences in the molecular responses and physiological traits of the highly susceptible Pinus radiata and the highly resistant Pinus pinea to F. circinatum at an early stage of infection. Dual RNA-Seq analysis also allowed to evaluate pathogen behavior when infecting each pine species. No significant changes in the physiological analysis were found upon pathogen infection, although transcriptional reprogramming was observed mainly in the resistant species. The transcriptome profiling of P. pinea revealed an early perception of the pathogen infection together with a strong and coordinated defense activation through the reinforcement and lignification of the cell wall, the antioxidant activity, the induction of PR genes, and the biosynthesis of defense hormones. On the contrary, P. radiata had a weaker response, possibly due to impaired perception of the fungal infection that led to a reduced downstream defense signaling. Fusarium circinatum showed a different transcriptomic profile depending on the pine species being infected. While in P. pinea, the pathogen focused on the degradation of plant cell walls, active uptake of the plant nutrients was showed in P. radiata. These findings present useful knowledge for the development of breeding programs to manage PPC.
Collapse
Affiliation(s)
- Cristina Zamora-Ballesteros
- Sustainable Forest Management Research Institute, University of Valladolid—INIA, 34004 Palencia, Spain; (J.J.D.); (J.M.-G.)
- Department of Vegetal Production and Forest Resources, University of Valladolid, 34004 Palencia, Spain
| | - Gloria Pinto
- Centre for Environmental and Marine Studies, CESAM, Department of Biology, University of Aveiro, 3810-193 Aveiro, Portugal; (G.P.); (J.A.); (A.A.)
| | - Joana Amaral
- Centre for Environmental and Marine Studies, CESAM, Department of Biology, University of Aveiro, 3810-193 Aveiro, Portugal; (G.P.); (J.A.); (A.A.)
| | - Luis Valledor
- Department of Organisms and Systems Biology, University of Oviedo, 33071 Oviedo, Spain;
| | - Artur Alves
- Centre for Environmental and Marine Studies, CESAM, Department of Biology, University of Aveiro, 3810-193 Aveiro, Portugal; (G.P.); (J.A.); (A.A.)
| | - Julio J. Diez
- Sustainable Forest Management Research Institute, University of Valladolid—INIA, 34004 Palencia, Spain; (J.J.D.); (J.M.-G.)
- Department of Vegetal Production and Forest Resources, University of Valladolid, 34004 Palencia, Spain
| | - Jorge Martín-García
- Sustainable Forest Management Research Institute, University of Valladolid—INIA, 34004 Palencia, Spain; (J.J.D.); (J.M.-G.)
- Department of Vegetal Production and Forest Resources, University of Valladolid, 34004 Palencia, Spain
| |
Collapse
|
9
|
Rodrigues AM, Miguel C, Chaves I, António C. Mass spectrometry-based forest tree metabolomics. MASS SPECTROMETRY REVIEWS 2021; 40:126-157. [PMID: 31498921 DOI: 10.1002/mas.21603] [Citation(s) in RCA: 7] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 03/07/2019] [Accepted: 08/05/2019] [Indexed: 05/24/2023]
Abstract
Research in forest tree species has advanced slowly when compared with other agricultural crops and model organisms, mainly due to the long-life cycles, large genome sizes, and lack of genomic tools. Additionally, trees are complex matrices, and the presence of interferents (e.g., oleoresins and cellulose) challenges the analysis of tree tissues with mass spectrometry (MS)-based analytical platforms. In this review, advances in MS-based forest tree metabolomics are discussed. Given their economic and ecological significance, particular focus is given to Pinus, Quercus, and Eucalyptus forest tree species to better understand their metabolite responses to abiotic and biotic stresses in the current climate change scenario. Furthermore, MS-based metabolomics technologies produce large and complex datasets that require expertize to adequately manage, process, analyze, and store the data in dedicated repositories. To ensure that the full potential of forest tree metabolomics data are translated into new knowledge, these data should comply with the FAIR principles (i.e., Findable, Accessible, Interoperable, and Re-usable). It is essential that adequate standards are implemented to annotate metadata from forest tree metabolomics studies as is already required by many science and governmental agencies and some major scientific publishers. © 2019 John Wiley & Sons Ltd. Mass Spec Rev 40:126-157, 2021.
Collapse
Affiliation(s)
- Ana Margarida Rodrigues
- Plant Metabolomics Laboratory, GreenIT-Bioresources for Sustainability, Instituto de Tecnologia Química e Biológica António Xavie, Universidade Nova de Lisboa (ITQB NOVA) Avenida da República, Oeiras, 2780-157, Portugal
| | - Célia Miguel
- Forest Genomics & Molecular Genetics Lab, BioISI-Biosystems & Integrative Sciences Institute, Faculty of Sciences, University of Lisboa, 1749-016, Lisboa, Portugal
- Instituto de Biologia Experimental e Tecnológica (iBET), 2780-157, Oeiras, Portugal
| | - Inês Chaves
- Forest Genomics & Molecular Genetics Lab, BioISI-Biosystems & Integrative Sciences Institute, Faculty of Sciences, University of Lisboa, 1749-016, Lisboa, Portugal
- Instituto de Biologia Experimental e Tecnológica (iBET), 2780-157, Oeiras, Portugal
| | - Carla António
- Plant Metabolomics Laboratory, GreenIT-Bioresources for Sustainability, Instituto de Tecnologia Química e Biológica António Xavie, Universidade Nova de Lisboa (ITQB NOVA) Avenida da República, Oeiras, 2780-157, Portugal
| |
Collapse
|
10
|
Šķipars V, Ruņģis D. Transcript Dynamics in Wounded and Inoculated Scots Pine. Int J Mol Sci 2021; 22:ijms22041505. [PMID: 33546141 PMCID: PMC7913219 DOI: 10.3390/ijms22041505] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/29/2020] [Revised: 01/28/2021] [Accepted: 01/30/2021] [Indexed: 11/16/2022] Open
Abstract
Comparative transcriptome analysis provides a useful tool for the exploration of plant-pathogen interaction by allowing in-depth comparison of gene expression between unaffected, inoculated and wounded organisms. Here we present the results of comparative transcriptome analysis in genetically identical one-year-old Scots pine ramets after wounding and inoculation with Heterobasidion annosum. We identified 230 genes that were more than 2-fold upregulated in inoculated samples (compared to controls) and 116 downregulated genes. Comparison of inoculated samp les with wounded samples identified 32 differentially expressed genes (30 were upregulated after inoculation). Several of the genes upregulated after inoculation are involved in protection from oxidative stress, while genes involved in photosynthesis, water transport and drought stress tolerance were downregulated. An NRT3 family protein was the most upregulated transcript in response to both inoculation and wounding, while a U-box domain-containing protein gene was the most upregulated gene comparing inoculation to wounding. The observed transcriptome dynamics suggest involvement of auxin, ethylene, jasmonate, gibberellin and reactive oxygen species pathways and cell wall modification regulation in response to H. annosum infection. The results are compared to methyl jasmonate induced transcriptome dynamics.
Collapse
|
11
|
Dong N, Bandura J, Zhang Z, Wang Y, Labadie K, Noel B, Davison A, Koene JM, Sun HS, Coutellec MA, Feng ZP. Ion channel profiling of the Lymnaea stagnalis ganglia via transcriptome analysis. BMC Genomics 2021; 22:18. [PMID: 33407100 PMCID: PMC7789530 DOI: 10.1186/s12864-020-07287-2] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/04/2020] [Accepted: 11/28/2020] [Indexed: 12/27/2022] Open
Abstract
BACKGROUND The pond snail Lymnaea stagnalis (L. stagnalis) has been widely used as a model organism in neurobiology, ecotoxicology, and parasitology due to the relative simplicity of its central nervous system (CNS). However, its usefulness is restricted by a limited availability of transcriptome data. While sequence information for the L. stagnalis CNS transcripts has been obtained from EST libraries and a de novo RNA-seq assembly, the quality of these assemblies is limited by a combination of low coverage of EST libraries, the fragmented nature of de novo assemblies, and lack of reference genome. RESULTS In this study, taking advantage of the recent availability of a preliminary L. stagnalis genome, we generated an RNA-seq library from the adult L. stagnalis CNS, using a combination of genome-guided and de novo assembly programs to identify 17,832 protein-coding L. stagnalis transcripts. We combined our library with existing resources to produce a transcript set with greater sequence length, completeness, and diversity than previously available ones. Using our assembly and functional domain analysis, we profiled L. stagnalis CNS transcripts encoding ion channels and ionotropic receptors, which are key proteins for CNS function, and compared their sequences to other vertebrate and invertebrate model organisms. Interestingly, L. stagnalis transcripts encoding numerous putative Ca2+ channels showed the most sequence similarity to those of Mus musculus, Danio rerio, Xenopus tropicalis, Drosophila melanogaster, and Caenorhabditis elegans, suggesting that many calcium channel-related signaling pathways may be evolutionarily conserved. CONCLUSIONS Our study provides the most thorough characterization to date of the L. stagnalis transcriptome and provides insights into differences between vertebrates and invertebrates in CNS transcript diversity, according to function and protein class. Furthermore, this study provides a complete characterization of the ion channels of Lymnaea stagnalis, opening new avenues for future research on fundamental neurobiological processes in this model system.
Collapse
Affiliation(s)
- Nancy Dong
- Department of Physiology, University of Toronto, 3308 MSB, 1 King's College Circle, Toronto, ON, M5S 1A8, Canada
| | - Julia Bandura
- Department of Physiology, University of Toronto, 3308 MSB, 1 King's College Circle, Toronto, ON, M5S 1A8, Canada
| | - Zhaolei Zhang
- Donnelly Centre for Cellular and Biomolecular Research and Department of Molecular Genetics, University of Toronto, Toronto, ON, M5S 3E1, Canada
| | - Yan Wang
- Department of Ecology and Evolutionary Biology, University of Toronto, Toronto, Ontario, M5S 3B2, Canada
- Department of Biological Sciences, University of Toronto Scarborough, Toronto, Ontario, M1C 1A4, Canada
| | - Karine Labadie
- Genoscope, Institut de biologie François Jacob, Commissariat à l'Energie Atomique (CEA), Université Paris-Saclay, BP5706, 91057, Evry, France
| | - Benjamin Noel
- Génomique Métabolique, Genoscope, Institut François Jacob, CEA, CNRS, University of Evry, Université Paris-Saclay, 91057, Evry, France
| | - Angus Davison
- School of Life Sciences, University of Nottingham, University Park, Nottingham, UK, NG7 2RD, UK
| | - Joris M Koene
- Department of Ecological Science, Faculty of Science, Vrije Universiteit, Amsterdam, The Netherlands
| | - Hong-Shuo Sun
- Department of Physiology, University of Toronto, 3308 MSB, 1 King's College Circle, Toronto, ON, M5S 1A8, Canada
- Department of Surgery, University of Toronto, Toronto, Ontario, M5S 1A8, Canada
| | | | - Zhong-Ping Feng
- Department of Physiology, University of Toronto, 3308 MSB, 1 King's College Circle, Toronto, ON, M5S 1A8, Canada.
| |
Collapse
|
12
|
Modern Approaches for Transcriptome Analyses in Plants. ADVANCES IN EXPERIMENTAL MEDICINE AND BIOLOGY 2021; 1346:11-50. [DOI: 10.1007/978-3-030-80352-0_2] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 10/19/2022]
|
13
|
Marasco V, Fusani L, Pola G, Smith S. Data on the de novo transcriptome assembly for the migratory bird, the Common quail ( Coturnix coturnix). Data Brief 2020; 32:106041. [PMID: 32775561 PMCID: PMC7394769 DOI: 10.1016/j.dib.2020.106041] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/11/2020] [Revised: 07/13/2020] [Accepted: 07/14/2020] [Indexed: 10/25/2022] Open
Abstract
The Common or European quail (Coturnix coturnix) is a Galliform bird of ecological importance for research in the field of animal migration. The Common quail is also a game bird, thus, of great interest for recreational activities and conservation management. Here, we generated a high quality de novo transcriptome for the Common quail for which no reference genome is to date publicly available. The transcriptome was obtained from a population of Common quail originated from captive founders raised under laboratory conditions. Paired-end RNA-Sequencing reads were obtained from extracted total RNA of brain tissue punches (preoptic-hypothalamic region) of 23 quails, which yielded to 5.5-11.2 million reads per individual bird for a total of 236 million reads. After assembly optimization, we used a stringent filtering analysis pipeline to remove redundant and low confidence transcripts. The final transcriptome consisted of 22,293 transcripts of which 21,551 (97%) were provided with annotation data. Our data offers a high quality pipeline for compiling transcriptomes of complex non-genomic species. Our data also provide a robust reference for gene expression studies in this species or other related Galliform species, including the Japanese quail.
Collapse
Affiliation(s)
- Valeria Marasco
- Konrad Lorenz Institute of Ethology, Department of Interdisciplinary Life Sciences, University of Veterinary Medicine Vienna, Savoyenstrasse 1a, Vienna 1160, Austria
| | - Leonida Fusani
- Konrad Lorenz Institute of Ethology, Department of Interdisciplinary Life Sciences, University of Veterinary Medicine Vienna, Savoyenstrasse 1a, Vienna 1160, Austria.,Department of Behavioural and Cognitive Biology, University of Vienna, Althanstrasse 14, Vienna 1090, Austria
| | - Gianni Pola
- Istituto Sperimentale Zootecnico per la Sicilia, via roccazzo 85, Palermo 90135, Italy
| | - Steve Smith
- Konrad Lorenz Institute of Ethology, Department of Interdisciplinary Life Sciences, University of Veterinary Medicine Vienna, Savoyenstrasse 1a, Vienna 1160, Austria
| |
Collapse
|
14
|
Maldonado JA, Firneno TJ, Roelke CE, Rains ND, Mwgiri J, Fujita MK. Transcriptome sequencing reveals signatures of positive selection in the Spot-Tailed Earless Lizard. PLoS One 2020; 15:e0234504. [PMID: 32542006 PMCID: PMC7295237 DOI: 10.1371/journal.pone.0234504] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/17/2020] [Accepted: 05/26/2020] [Indexed: 11/21/2022] Open
Abstract
The continual loss of threatened biodiversity is occurring at an accelerated pace. High-throughput sequencing technologies are now providing opportunities to address this issue by aiding in the generation of molecular data for many understudied species of high conservation interest. Our overall goal of this study was to begin building the genomic resources to continue investigations and conservation of the Spot-Tailed Earless lizard. Here we leverage the power of high-throughput sequencing to generate the liver transcriptome for the Northern Spot-Tailed Earless Lizard (Holbrookia lacerata) and Southern Spot-Tailed Earless Lizard (Holbrookia subcaudalis), which have declined in abundance in the past decades, and their sister species, the Common Lesser Earless Lizard (Holbrookia maculata). Our efforts produced high quality and robust transcriptome assemblies validated by 1) quantifying the number of processed reads represented in the transcriptome assembly and 2) quantifying the number of highly conserved single-copy orthologs that are present in our transcript set using the BUSCO pipeline. We found 1,361 1-to-1 orthologs among the three Holbrookia species, Anolis carolinensis, and Sceloporus undulatus. We carried out dN/dS selection tests using a branch-sites model and identified a dozen genes that experienced positive selection in the Holbrookia lineage with functions in development, immunity, and metabolism. Our single-copy orthologous sequences additionally revealed significant pairwise sequence divergence (~.73%) between the Northern H. lacerata and Southern H. subcaudalis that further supports the recent elevation of the Southern Spot-Tailed Earless Lizard to full species.
Collapse
Affiliation(s)
- Jose A. Maldonado
- Department of Biology, Amphibian and Reptile Diversity Research Center, The University of Texas at Arlington, Arlington, TX, United States of America
| | - Thomas J. Firneno
- Department of Biology, Amphibian and Reptile Diversity Research Center, The University of Texas at Arlington, Arlington, TX, United States of America
| | - Corey E. Roelke
- Department of Biology, Amphibian and Reptile Diversity Research Center, The University of Texas at Arlington, Arlington, TX, United States of America
| | - Nathan D. Rains
- Texas Parks and Wildlife Department, Austin, Texas, TX, United States of America
| | - Juliet Mwgiri
- Department of Biology, Amphibian and Reptile Diversity Research Center, The University of Texas at Arlington, Arlington, TX, United States of America
| | - Matthew K. Fujita
- Department of Biology, Amphibian and Reptile Diversity Research Center, The University of Texas at Arlington, Arlington, TX, United States of America
- * E-mail:
| |
Collapse
|
15
|
Sadat-Hosseini M, Bakhtiarizadeh MR, Boroomand N, Tohidfar M, Vahdati K. Combining independent de novo assemblies to optimize leaf transcriptome of Persian walnut. PLoS One 2020; 15:e0232005. [PMID: 32343733 PMCID: PMC7188282 DOI: 10.1371/journal.pone.0232005] [Citation(s) in RCA: 17] [Impact Index Per Article: 4.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/27/2019] [Accepted: 04/06/2020] [Indexed: 12/22/2022] Open
Abstract
Transcriptome resources can facilitate to increase yield and quality of walnuts. Finding the best transcriptome assembly has not been the subject of walnuts research as yet. This research generated 240,179,782 reads from 11 walnut leaves according to cDNA libraries. The reads provided a complete de novo transcriptome assembly. Fifteen different transcriptome assemblies were constructed from five different well-known assemblers used in scientific literature with different k-mer lengths (Bridger, BinPacker, SOAPdenovo-Trans, Trinity and SPAdes) as well as two merging approaches (EvidentialGene and Transfuse). Based on the four quality metrics of assembly, the results indicated an efficiency in the process of merging the assemblies after being generated by de novo assemblers. Finally, EvidentialGene was recognized as the best assembler for the de novo assembly of the leaf transcriptome in walnut. Among a total number of 183,191 transcripts which were generated by EvidentialGene, there were 109,413 transcripts capable of protein potential (59.72%) and 104,926 were recognized as ORFs (57.27%). In addition, 79,185 transcripts were predicted to exist with at least one hit to the Pfam database. A number of 3,931 transcription factors were identified by BLAST searching against PlnTFDB. Furthermore, 6,591 of the predicted peptide sequences contained signaling peptides, while 92,704 contained transmembrane domains. Comparison of the assembled transcripts with transcripts of the walnut and published genome assembly for the 'Chandler' cultivar using the BLAST algorithm led to identify a total number of 27,304 and 19,178 homologue transcripts, respectively. De novo transcriptomes in walnut leaves can be developed for the future studies in functional genomics and genetic studies of walnuts.
Collapse
Affiliation(s)
- Mohammad Sadat-Hosseini
- Department of Horticulture, College of Aburaihan, University of Tehran, Tehran, Iran
- Department of Horticulture, Faculty of Agriculture, University of Jiroft, Jiroft, Iran
| | | | - Naser Boroomand
- Department of Soil Science, Faculty of Agriculture, Shahid Bahonar University of Kerman, Kerman, Iran
| | - Masoud Tohidfar
- Department of Plant Biotechnology, Faculty of Life Science and Biotechnology, Shahid Beheshti University, Tehran, Iran
| | - Kourosh Vahdati
- Department of Horticulture, College of Aburaihan, University of Tehran, Tehran, Iran
| |
Collapse
|
16
|
Hernandez-Escribano L, Visser EA, Iturritxa E, Raposo R, Naidoo S. The transcriptome of Pinus pinaster under Fusarium circinatum challenge. BMC Genomics 2020; 21:28. [PMID: 31914917 PMCID: PMC6950806 DOI: 10.1186/s12864-019-6444-0] [Citation(s) in RCA: 16] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/26/2019] [Accepted: 12/30/2019] [Indexed: 01/19/2023] Open
Abstract
BACKGROUND Fusarium circinatum, the causal agent of pitch canker disease, poses a serious threat to several Pinus species affecting plantations and nurseries. Although Pinus pinaster has shown moderate resistance to F. circinatum, the molecular mechanisms of defense in this host are still unknown. Phytohormones produced by the plant and by the pathogen are known to play a crucial role in determining the outcome of plant-pathogen interactions. Therefore, the aim of this study was to determine the role of phytohormones in F. circinatum virulence, that compromise host resistance. RESULTS A high quality P. pinaster de novo transcriptome assembly was generated, represented by 24,375 sequences from which 17,593 were full length genes, and utilized to determine the expression profiles of both organisms during the infection process at 3, 5 and 10 days post-inoculation using a dual RNA-sequencing approach. The moderate resistance shown by Pinus pinaster at the early time points may be explained by the expression profiles pertaining to early recognition of the pathogen, the induction of pathogenesis-related proteins and the activation of complex phytohormone signaling pathways that involves crosstalk between salicylic acid, jasmonic acid, ethylene and possibly auxins. Moreover, the expression of F. circinatum genes related to hormone biosynthesis suggests manipulation of the host phytohormone balance to its own benefit. CONCLUSIONS We hypothesize three key steps of host manipulation: perturbing ethylene homeostasis by fungal expression of genes related to ethylene biosynthesis, blocking jasmonic acid signaling by coronatine insensitive 1 (COI1) suppression, and preventing salicylic acid biosynthesis from the chorismate pathway by the synthesis of isochorismatase family hydrolase (ICSH) genes. These results warrant further testing in F. circinatum mutants to confirm the mechanism behind perturbing host phytohormone homeostasis.
Collapse
Affiliation(s)
- Laura Hernandez-Escribano
- Instituto Nacional de Investigación y Tecnología Agraria y Alimentaria, Centro de Investigación Forestal (INIA-CIFOR), Madrid, Spain
- Departamento de Biotecnología-Biología Vegetal, Escuela Técnica Superior de Ingeniería Agronómica, Alimentaria y de Biosistemas, Universidad Politécnica de Madrid, Madrid, Spain
| | - Erik A Visser
- Department of Biochemistry, Genetics and Microbiology, Forestry and Agricultural Biotechnology Institute (FABI), Centre for Bioinformatics and Computational Biology, University of Pretoria, Pretoria, South Africa
| | - Eugenia Iturritxa
- NEIKER, Granja Modelo de Arkaute, Apdo 46, 01080, Vitoria-Gasteiz, Spain
| | - Rosa Raposo
- Instituto Nacional de Investigación y Tecnología Agraria y Alimentaria, Centro de Investigación Forestal (INIA-CIFOR), Madrid, Spain
- Instituto de Gestión Forestal Sostenible (iuFOR), Universidad de Valladolid/INIA, Valladolid, Spain
| | - Sanushka Naidoo
- Department of Biochemistry, Genetics and Microbiology, Forestry and Agricultural Biotechnology Institute (FABI), Centre for Bioinformatics and Computational Biology, University of Pretoria, Pretoria, South Africa.
| |
Collapse
|
17
|
Visser EA, Wegrzyn JL, Steenkamp ET, Myburg AA, Naidoo S. Dual RNA-Seq Analysis of the Pine- Fusarium circinatum Interaction in Resistant ( Pinus tecunumanii) and Susceptible ( Pinus patula) Hosts. Microorganisms 2019; 7:E315. [PMID: 31487786 PMCID: PMC6780516 DOI: 10.3390/microorganisms7090315] [Citation(s) in RCA: 15] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/03/2019] [Revised: 08/15/2019] [Accepted: 08/20/2019] [Indexed: 12/15/2022] Open
Abstract
Fusarium circinatum poses a serious threat to many pine species in both commercial and natural pine forests. Knowledge regarding the molecular basis of pine-F. circinatum host-pathogen interactions could assist efforts to produce more resistant planting stock. This study aimed to identify molecular responses underlying resistance against F. circinatum. A dual RNA-seq approach was used to investigate host and pathogen expression in F. circinatum challenged Pinus tecunumanii (resistant) and Pinus patula (susceptible), at three- and seven-days post inoculation. RNA-seq reads were mapped to combined host-pathogen references for both pine species to identify differentially expressed genes (DEGs). F. circinatum genes expressed during infection showed decreased ergosterol biosynthesis in P. tecunumanii relative to P. patula. For P. tecunumanii, enriched gene ontologies and DEGs indicated roles for auxin-, ethylene-, jasmonate- and salicylate-mediated phytohormone signalling. Correspondingly, key phytohormone signaling components were down-regulated in P. patula. Key F. circinatum ergosterol biosynthesis genes were expressed at lower levels during infection of the resistant relative to the susceptible host. This study further suggests that coordination of phytohormone signaling is required for F. circinatum resistance in P. tecunumanii, while a comparatively delayed response and impaired phytohormone signaling contributes to susceptibility in P. patula.
Collapse
Affiliation(s)
- Erik A Visser
- Department of Biochemistry, Genetics and Microbiology, Forestry and Agricultural Biotechnology Institute (FABI), Centre for Bioinformatics and Computational Biology, University of Pretoria, Private Bag X20, Pretoria 0028, South Africa
| | - Jill L Wegrzyn
- Department of Ecology and Evolutionary Biology, University of Connecticut, Storrs, CT 06269, USA
| | - Emma T Steenkamp
- Department of Biochemistry, Genetics and Microbiology, Forestry and Agricultural Biotechnology Institute (FABI), Centre for Bioinformatics and Computational Biology, University of Pretoria, Private Bag X20, Pretoria 0028, South Africa
| | - Alexander A Myburg
- Department of Biochemistry, Genetics and Microbiology, Forestry and Agricultural Biotechnology Institute (FABI), Centre for Bioinformatics and Computational Biology, University of Pretoria, Private Bag X20, Pretoria 0028, South Africa
| | - Sanushka Naidoo
- Department of Biochemistry, Genetics and Microbiology, Forestry and Agricultural Biotechnology Institute (FABI), Centre for Bioinformatics and Computational Biology, University of Pretoria, Private Bag X20, Pretoria 0028, South Africa.
| |
Collapse
|
18
|
Eguiluz M, Kulcheski FR, Margis R, Guzman F. De novo assembly of Vriesea carinata leaf transcriptome to identify candidate cysteine-proteases. Gene 2019; 691:96-105. [PMID: 30630096 DOI: 10.1016/j.gene.2018.12.053] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/10/2018] [Revised: 12/13/2018] [Accepted: 12/20/2018] [Indexed: 11/30/2022]
Abstract
Vriesea carinata is an endemic bromeliad from the Brazilian Atlantic Forest. It has trichome and tank system in their leaves which allows to absorb water and nutrients. It belongs to Bromeliaceae family, which includes several species highly enriched of cysteine-proteases (CysPs). These proteolytic enzymes regulate processes as senescence, cell differentiation, pathogen-linked programmed cell death and mobilization of proteins. Although, their biological importance, there are not genomic resources in V. carinata that can help to identify and understand their molecular mechanisms involved in different biological processes. Thus high-throughput transcriptome sequencing of V. carinata is necessary to generate sequences for the purpose of gene discovery and functional genomic studies. In the present study, we sequenced and assembled the V. carinata transcriptome to the identification of CysPs. A total of 43,232 contigs were assembled for the leaf tissue. BLAST analysis indicated that 23,803 contigs exhibited similarity to non-redundant Viridiplantae proteins. 28.24% of the contigs were classified into the COG database, and gene ontology categorized them into 61 functional groups. A metabolic pathway analysis with KEGG revealed 9679 contigs assigned to 31 metabolic pathways. Among 16 full-length CysPs identified, 11 were evaluated in respect to their expression patterns in the leaf apex, base and inflorescence tissues. The results showed differential expression levels of legumain, metacaspase, pyroglutamyl and papain-like CysPs depending of the leaf region. These results provide a global overview of V. carinata gene functions and expression activities of CysPs in those tissues.
Collapse
Affiliation(s)
- M Eguiluz
- PPGBM, Departamento de Genética, Universidade Federal do Rio Grande do Sul-UFRGS, Porto Alegre, Brazil; Facultad de Ciencias, Universidad Peruana Cayetano Heredia, Lima, Peru
| | - F R Kulcheski
- PPGBCD, Departamento de Biologia Celular, Genética e Embriologia, Universidade Federal de Santa Catarina-UFSC, Florianópolis, Brazil
| | - R Margis
- PPGBM, Departamento de Genética, Universidade Federal do Rio Grande do Sul-UFRGS, Porto Alegre, Brazil; PPGBCM, Centro de Biotecnologia, Universidade Federal do Rio Grande do Sul-UFRGS, Porto Alegre, Brazil; Departamento de Biofísica, Universidade Federal do Rio Grande do Sul-UFRGS, Porto Alegre, Brazil
| | - F Guzman
- PPGBCM, Centro de Biotecnologia, Universidade Federal do Rio Grande do Sul-UFRGS, Porto Alegre, Brazil.
| |
Collapse
|
19
|
Rodrigues AS, De Vega JJ, Miguel CM. Comprehensive assembly and analysis of the transcriptome of maritime pine developing embryos. BMC PLANT BIOLOGY 2018; 18:379. [PMID: 30594130 PMCID: PMC6310951 DOI: 10.1186/s12870-018-1564-2] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 05/06/2018] [Accepted: 11/22/2018] [Indexed: 05/05/2023]
Abstract
BACKGROUND There are clear differences in embryo development between angiosperm and gymnosperm species. Most of the current knowledge on gene expression and regulation during plant embryo development has derived from studies on angiosperms species, in particular from the model plant Arabidopsis thaliana. The few published studies on transcript profiling of conifer embryogenesis show the existence of many putative embryo-specific transcripts without an assigned function. In order to extend the knowledge on the transcriptomic expression during conifer embryogenesis, we sequenced the transcriptome of zygotic embryos for several developmental stages that cover most of Pinus pinaster (maritime pine) embryogenesis. RESULTS Total RNA samples collected from five zygotic embryo developmental stages were sequenced with Illumina technology. A de novo transcriptome was assembled as no genome sequence is yet published for Pinus pinaster. The transcriptome of reference for the period of zygotic embryogenesis in maritime pine contains 67,429 transcripts, which likely encode 58,527 proteins. The annotation shows a significant percentage, 31%, of predicted proteins exclusively present in pine embryogenesis. Functional categories and enrichment analysis of the differentially expressed transcripts evidenced carbohydrate transport and metabolism over-representation in early embryo stages, as highlighted by the identification of many putative glycoside hydrolases, possibly associated with cell wall modification, and carbohydrate transport transcripts. Moreover, the predominance of chromatin remodelling events was detected in early to middle embryogenesis, associated with an active synthesis of histones and their post-translational modifiers related to increased transcription, as well as silencing of transposons. CONCLUSIONS Our results extend the understanding of gene expression and regulation during zygotic embryogenesis in conifers and are a valuable resource to support further improvements in somatic embryogenesis for vegetative propagation of conifer species. Specific transcripts associated with carbohydrate metabolism, monosaccharide transport and epigenetic regulation seem to play an important role in pine early embryogenesis and may be a source of reliable molecular markers for early embryogenesis.
Collapse
Affiliation(s)
- Andreia S. Rodrigues
- Instituto de Biologia Experimental e Tecnológica (iBET), Apartado 12, 2780-901 Oeiras, Portugal
- Instituto de Tecnologia Química e Biológica António Xavier, Universidade Nova de Lisboa (ITQB NOVA), Av. da República, 2780-157 Oeiras, Portugal
| | - José J. De Vega
- Earlham Institute, Norwich Research Park, Norwich, NR4 7UZ UK
| | - Célia M. Miguel
- Instituto de Biologia Experimental e Tecnológica (iBET), Apartado 12, 2780-901 Oeiras, Portugal
- Instituto de Tecnologia Química e Biológica António Xavier, Universidade Nova de Lisboa (ITQB NOVA), Av. da República, 2780-157 Oeiras, Portugal
- Universidade de Lisboa, Faculdade de Ciências, BioISI - Biosystems & Integrative Sciences Institute, Campo Grande, 1749-016 Lisbon, Portugal
| |
Collapse
|
20
|
Fay JV, Watkins CJ, Shrestha RK, Litwiñiuk SL, Talavera Stefani LN, Rojas CA, Argüelles CF, Ferreras JA, Caccamo M, Miretti MM. Yerba mate (Ilex paraguariensis, A. St.-Hil.) de novo transcriptome assembly based on tissue specific genomic expression profiles. BMC Genomics 2018; 19:891. [PMID: 30526481 PMCID: PMC6286616 DOI: 10.1186/s12864-018-5240-6] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/02/2018] [Accepted: 11/12/2018] [Indexed: 12/27/2022] Open
Abstract
Background The most common infusion in southern Latin-American countries is prepared with dried leaves of Ilex paraguariensis A. St.-Hil., an aboriginal ancestral beverage known for its high polyphenols concentration currently consumed in > 90% of homes in Argentina, in Paraguay and Uruguay. The economy of entire provinces heavily relies on the production, collection and manufacture of Ilex paraguariensis, the fifth plant species with highest antioxidant activity. Polyphenols are associated to relevant health benefits including strong antioxidant properties. Despite its regional relevance and potential biotechnological applications, little is known about functional genomics and genetics underlying phenotypic variation of relevant traits. By generating tissue specific transcriptomic profiles, we aimed to comprehensively annotate genes in the Ilex paraguariensis phenylpropanoid pathway and to evaluate differential expression profiles. Results In this study we generated a reliable transcriptome assembly based on a collection of 15 RNA-Seq libraries from different tissues of Ilex paraguariensis. A total of 554 million RNA-Seq reads were assembled into 193,897 transcripts, where 24,612 annotated full-length transcripts had complete ORF. We assessed the transcriptome assembly quality, completeness and accuracy using BUSCO and TransRate; consistency was also evaluated by experimentally validating 11 predicted genes by PCR and sequencing. Functional annotation against KEGG Pathway database identified 1395 unigenes involved in biosynthesis of secondary metabolites, 531 annotated transcripts corresponded to the phenylpropanoid pathway. The top 30 differentially expressed genes among tissue revealed genes involved in photosynthesis and stress response. These significant differences were then validated by qRT-PCR. Conclusions Our study is the first to provide data from whole genome gene expression profiles in different Ilex paraguariensis tissues, experimentally validating in-silico predicted genes key to the phenylpropanoid (antioxidant) pathway. Our results provide essential genomic data of potential use in breeding programs for polyphenol content. Further studies are necessary to assess if the observed expression variation in the phenylpropanoid pathway annotated genes is related to variations in leaves’ polyphenol content at the population scale. These results set the current reference for Ilex paraguariensis genomic studies and provide a substantial contribution to research and biotechnological applications of phenylpropanoid secondary metabolites. Electronic supplementary material The online version of this article (10.1186/s12864-018-5240-6) contains supplementary material, which is available to authorized users.
Collapse
Affiliation(s)
- Jessica V Fay
- Grupo de Investigación en Genética Aplicada (GIGA), Facultad de Ciencias Exactas Químicas y Naturales, Instituto de Biología Subtropical (IBS UNaM-CONICET), Universidad Nacional de Misiones, Jujuy 1745, CP3300, Posadas, Misiones, Argentina
| | - Christopher J Watkins
- The Genome Analysis Centre, Norwich Research Park, Norwich, NR4 7UH, UK.,Present address: Earlham Institute, Norwich Research Park, Norwich, NR4 7UZ, UK
| | - Ram K Shrestha
- The Genome Analysis Centre, Norwich Research Park, Norwich, NR4 7UH, UK.,Present address: Earlham Institute, Norwich Research Park, Norwich, NR4 7UZ, UK
| | - Sergio L Litwiñiuk
- Grupo de Investigación en Genética Aplicada (GIGA), Facultad de Ciencias Exactas Químicas y Naturales, Instituto de Biología Subtropical (IBS UNaM-CONICET), Universidad Nacional de Misiones, Jujuy 1745, CP3300, Posadas, Misiones, Argentina
| | - Liliana N Talavera Stefani
- Grupo de Investigación en Genética Aplicada (GIGA), Facultad de Ciencias Exactas Químicas y Naturales, Instituto de Biología Subtropical (IBS UNaM-CONICET), Universidad Nacional de Misiones, Jujuy 1745, CP3300, Posadas, Misiones, Argentina
| | - Cristian A Rojas
- Universidad Federal de la Integración Latinoamericana, Foz de Iguazú, PR, Brazil
| | - Carina F Argüelles
- Grupo de Investigación en Genética Aplicada (GIGA), Facultad de Ciencias Exactas Químicas y Naturales, Instituto de Biología Subtropical (IBS UNaM-CONICET), Universidad Nacional de Misiones, Jujuy 1745, CP3300, Posadas, Misiones, Argentina
| | - Julian A Ferreras
- Grupo de Investigación en Genética Aplicada (GIGA), Facultad de Ciencias Exactas Químicas y Naturales, Instituto de Biología Subtropical (IBS UNaM-CONICET), Universidad Nacional de Misiones, Jujuy 1745, CP3300, Posadas, Misiones, Argentina
| | - Mario Caccamo
- The Genome Analysis Centre, Norwich Research Park, Norwich, NR4 7UH, UK.,Present address: NIAB, Huntingdon Road, Cambridge, CB3 0LE, UK
| | - Marcos M Miretti
- Grupo de Investigación en Genética Aplicada (GIGA), Facultad de Ciencias Exactas Químicas y Naturales, Instituto de Biología Subtropical (IBS UNaM-CONICET), Universidad Nacional de Misiones, Jujuy 1745, CP3300, Posadas, Misiones, Argentina.
| |
Collapse
|
21
|
Payá-Milans M, Olmstead JW, Nunez G, Rinehart TA, Staton M. Comprehensive evaluation of RNA-seq analysis pipelines in diploid and polyploid species. Gigascience 2018; 7:5168871. [PMID: 30418578 PMCID: PMC6275443 DOI: 10.1093/gigascience/giy132] [Citation(s) in RCA: 16] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/19/2018] [Accepted: 10/21/2018] [Indexed: 11/12/2022] Open
Abstract
Background The usual analysis of RNA sequencing (RNA-seq) reads is based on an existing reference genome and annotated gene models. However, when a reference for the sequenced species is not available, alternatives include using a reference genome from a related species or reconstructing transcript sequences with de novo assembly. In addition, researchers are faced with many options for RNA-seq data processing and limited information on how their decisions will impact the final outcome. Using both a diploid and polyploid species with a distant reference genome, we have tested the influence of different tools at various steps of a typical RNA-seq analysis workflow on the recovery of useful processed data available for downstream analysis. Findings At the preprocessing step, we found error correction has a strong influence on de novo assembly but not on mapping results. After trimming, a greater percentage of reads could be used in downstream analysis by selecting gentle quality trimming performed with Skewer instead of strict quality trimming with Trimmomatic. This availability of reads correlated with size, quality, and completeness of de novo assemblies and with number of mapped reads. When selecting a reference genome from a related species to map reads, outcome was significantly improved when using mapping software tolerant of greater sequence divergence, such as Stampy or GSNAP. Conclusions The selection of bioinformatic software tools for RNA-seq data analysis can maximize quality parameters on de novo assemblies and availability of reads in downstream analysis.
Collapse
Affiliation(s)
- Miriam Payá-Milans
- Department of Entomology and Plant Pathology, University of Tennessee, 370 PBB, 2505 EJ Chapman Blvd, Knoxville, TN, 37996, United States
| | - James W Olmstead
- Horticultural Sciences Department, University of Florida, 2550 Hull Rd, PO Box 110690, Gainesville, FL, 32611, United States
| | - Gerardo Nunez
- Horticultural Sciences Department, University of Florida, 2550 Hull Rd, PO Box 110690, Gainesville, FL, 32611, United States
| | - Timothy A Rinehart
- Thad Cochran Southern Horticultural Laboratory, USDA-Agricultural Research Service, PO Box 287, Poplarville, MS, 39470, United States.,Crop Production and Protection, USDA-Agricultural Research Service, 5601 Sunnyside Ave, Beltsville, MD, 20705, United States
| | - Margaret Staton
- Department of Entomology and Plant Pathology, University of Tennessee, 370 PBB, 2505 EJ Chapman Blvd, Knoxville, TN, 37996, United States
| |
Collapse
|
22
|
Visser EA, Wegrzyn JL, Myburg AA, Naidoo S. Defence transcriptome assembly and pathogenesis related gene family analysis in Pinus tecunumanii (low elevation). BMC Genomics 2018; 19:632. [PMID: 30139335 PMCID: PMC6108113 DOI: 10.1186/s12864-018-5015-0] [Citation(s) in RCA: 24] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/03/2018] [Accepted: 08/14/2018] [Indexed: 12/23/2022] Open
Abstract
BACKGROUND Fusarium circinatum is a pressing threat to the cultivation of many economically important pine tree species. Efforts to develop effective disease management strategies can be aided by investigating the molecular mechanisms involved in the host-pathogen interaction between F. circinatum and pine species. Pinus tecunumanii and Pinus patula are two closely related tropical pine species that differ widely in their resistance to F. circinatum challenge, being resistant and susceptible respectively, providing the potential for a useful pathosystem to investigate the molecular responses underlying resistance to F. circinatum. However, no genomic resources are available for P. tecunumanii. Pathogenesis-related proteins are classes of proteins that play important roles in plant-microbe interactions, e.g. chitinases; proteins that break down the major structural component of fungal cell walls. Generating a reference sequence for P. tecunumanii and characterizing pathogenesis related gene families in these two pine species is an important step towards unravelling the pine-F. circinatum interaction. RESULTS Eight reference based and 12 de novo assembled transcriptomes were produced, for juvenile shoot tissue from both species. EvidentialGene pipeline redundancy reduction, expression filtering, protein clustering and taxonomic filtering produced a 50 Mb shoot transcriptome consisting of 28,621 contigs for P. tecunumanii and a 72 Mb shoot transcriptome consisting of 52,735 contigs for P. patula. Predicted protein sequences encoded by the assembled transcriptomes were clustered with reference proteomes from 92 other species to identify pathogenesis related gene families in P. patula, P. tecunumanii and other pine species. CONCLUSIONS The P. tecunumanii transcriptome is the first gene catalogue for the species, representing an important resource for studying resistance to the pitch canker pathogen, F. circinatum. This study also constitutes, to our knowledge, the largest index of gymnosperm PR-genes to date.
Collapse
Affiliation(s)
- Erik A. Visser
- Department of Biochemistry, Genetics and Microbiology, Forestry and Agricultural Biotechnology Institute (FABI), Genomics Research Institute (GRI), University of Pretoria, Private bag X20, Pretoria, 0028 South Africa
| | - Jill L. Wegrzyn
- Department of Ecology and Evolutionary Biology, University of Connecticut, Storrs, CT 06269 USA
| | - Alexander A. Myburg
- Department of Biochemistry, Genetics and Microbiology, Forestry and Agricultural Biotechnology Institute (FABI), Genomics Research Institute (GRI), University of Pretoria, Private bag X20, Pretoria, 0028 South Africa
| | - Sanushka Naidoo
- Department of Biochemistry, Genetics and Microbiology, Forestry and Agricultural Biotechnology Institute (FABI), Genomics Research Institute (GRI), University of Pretoria, Private bag X20, Pretoria, 0028 South Africa
| |
Collapse
|
23
|
Naidoo S, Christie N, Acosta JJ, Mphahlele MM, Payn KG, Myburg AA, Külheim C. Terpenes associated with resistance against the gall wasp, Leptocybe invasa, in Eucalyptus grandis. PLANT, CELL & ENVIRONMENT 2018; 41:1840-1851. [PMID: 29710389 DOI: 10.1111/pce.13323] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 10/09/2017] [Revised: 04/13/2018] [Accepted: 04/17/2018] [Indexed: 05/24/2023]
Abstract
Leptocybe invasa is an insect pest causing gall formation on oviposited shoot tips and leaves of Eucalyptus trees leading to leaf deformation, stunting, and death in severe cases. We previously observed different constitutive and induced terpenes, plant specialized metabolites that may act as attractants or repellents to insects, in a resistant and susceptible clone of Eucalyptus challenged with L. invasa. We tested the hypothesis that specific terpenes are associated with pest resistance in a Eucalyptus grandis half-sib population. Insect damage was scored over 2 infestation cycles, and leaves were harvested for near-infrared reflectance (NIR) and terpene measurements. We used Bayesian model averaging for terpene selection and obtained partial least squares NIR models to predict terpene content and L. invasa infestation damage. In our optimal model, 29% of the phenotypic variation could be explained by 7 terpenes, and the monoterpene combination, limonene, α-terpineol, and 1,8-cineole, could be predicted with an NIR prediction ability of .67. Bayesian model averaging supported α-pinene, γ-terpinene, and iso-pinocarveol as important for predicting L. invasa infestation. Susceptibility was associated with increased γ-terpinene and α-pinene, which may act as a pest attractant, whereas reduced susceptibility was associated with iso-pinocarveol, which may act to recruit parasitoids or have direct toxic effects.
Collapse
Affiliation(s)
- Sanushka Naidoo
- Department of Biochemistry, Genetics and Microbiology, Forestry and Agricultural Biotechnology Institute (FABI), University of Pretoria, Private bag x20, Pretoria, 0028, South Africa
| | - Nanette Christie
- Department of Biochemistry, Genetics and Microbiology, Forestry and Agricultural Biotechnology Institute (FABI), University of Pretoria, Private bag x20, Pretoria, 0028, South Africa
| | - Juan J Acosta
- Camcore, Department of Forestry and Environmental Resources, North Carolina State University, Raleigh, NC, 27695-8008, USA
| | | | - Kitt G Payn
- Mondi Forests, Trahar Technology Centre, P.O. Box 12, Hilton, 3245, South Africa
| | - Alexander A Myburg
- Department of Biochemistry, Genetics and Microbiology, Forestry and Agricultural Biotechnology Institute (FABI), University of Pretoria, Private bag x20, Pretoria, 0028, South Africa
| | - Carsten Külheim
- Research School of Biology, Australian National University, 46 Sullivans Creek Road, Canberra, 2601, Australian Capital Territory, Australia
| |
Collapse
|
24
|
Evkaikina AI, Berke L, Romanova MA, Proux-Wéra E, Ivanova AN, Rydin C, Pawlowski K, Voitsekhovskaja OV. The Huperzia selago Shoot Tip Transcriptome Sheds New Light on the Evolution of Leaves. Genome Biol Evol 2018; 9:2444-2460. [PMID: 28957460 PMCID: PMC5622374 DOI: 10.1093/gbe/evx169] [Citation(s) in RCA: 17] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Accepted: 08/28/2017] [Indexed: 12/21/2022] Open
Abstract
Lycopodiophyta—consisting of three orders, Lycopodiales, Isoetales and Selaginellales, with different types of shoot apical meristems (SAMs)—form the earliest branch among the extant vascular plants. They represent a sister group to all other vascular plants, from which they differ in that their leaves are microphylls—that is, leaves with a single, unbranched vein, emerging from the protostele without a leaf gap—not megaphylls. All leaves represent determinate organs originating on the flanks of indeterminate SAMs. Thus, leaf formation requires the suppression of indeterminacy, that is, of KNOX transcription factors. In seed plants, this is mediated by different groups of transcription factors including ARP and YABBY. We generated a shoot tip transcriptome of Huperzia selago (Lycopodiales) to examine the genes involved in leaf formation. Our H. selago transcriptome does not contain any ARP homolog, although transcriptomes of Selaginella spp. do. Surprisingly, we discovered a YABBY homolog, although these transcription factors were assumed to have evolved only in seed plants. The existence of a YABBY homolog in H. selago suggests that YABBY evolved already in the common ancestor of the vascular plants, and subsequently was lost in some lineages like Selaginellales, whereas ARP may have been lost in Lycopodiales. The presence of YABBY in the common ancestor of vascular plants would also support the hypothesis that this common ancestor had a simplex SAM. Furthermore, a comparison of the expression patterns of ARP in shoot tips of Selaginella kraussiana (Harrison CJ, etal. 2005. Independent recruitment of a conserved developmental mechanism during leaf evolution. Nature 434(7032):509–514.) and YABBY in shoot tips of H. selago implies that the development of microphylls, unlike megaphylls, does not seem to depend on the combined activities of ARP and YABBY. Altogether, our data show that Lycopodiophyta are a diverse group; so, in order to understand the role of Lycopodiophyta in evolution, representatives of Lycopodiales, Selaginellales, as well as of Isoetales, have to be examined.
Collapse
Affiliation(s)
- Anastasiia I Evkaikina
- Laboratory of Molecular and Ecological Physiology, Komarov Botanical Institute, Russian Academy of Sciences, St. Petersburg, Russia
| | - Lidija Berke
- Department of Plant Sciences, Wageningen University, The Netherlands
| | - Marina A Romanova
- Department of Botany, St. Petersburg State University, St. Petersburg, Russia
| | - Estelle Proux-Wéra
- Department of Plant Protection Biology, Swedish University of Agricultural Sciences, Alnarp, Sweden.,Science for Life Laboratory, Department of Biochemistry and Biophysics, Stockholm University, Solna, Sweden
| | - Alexandra N Ivanova
- Laboratory of Anatomy and Morphology, Komarov Botanical Institute, Russian Academy of Sciences, St. Petersburg, Russia
| | - Catarina Rydin
- Department of Ecology, Environment and Plant Sciences, Stockholm University, Stockholm, Sweden
| | - Katharina Pawlowski
- Department of Ecology, Environment and Plant Sciences, Stockholm University, Stockholm, Sweden
| | - Olga V Voitsekhovskaja
- Laboratory of Molecular and Ecological Physiology, Komarov Botanical Institute, Russian Academy of Sciences, St. Petersburg, Russia
| |
Collapse
|
25
|
Fox H, Doron-Faigenboim A, Kelly G, Bourstein R, Attia Z, Zhou J, Moshe Y, Moshelion M, David-Schwartz R. Transcriptome analysis of Pinus halepensis under drought stress and during recovery. TREE PHYSIOLOGY 2018; 38:423-441. [PMID: 29177514 PMCID: PMC5982726 DOI: 10.1093/treephys/tpx137] [Citation(s) in RCA: 29] [Impact Index Per Article: 4.8] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 04/10/2017] [Revised: 08/24/2017] [Accepted: 10/12/2017] [Indexed: 05/09/2023]
Abstract
Forest trees use various strategies to cope with drought stress and these strategies involve complex molecular mechanisms. Pinus halepensis Miller (Aleppo pine) is found throughout the Mediterranean basin and is one of the most drought-tolerant pine species. In order to decipher the molecular mechanisms that P. halepensis uses to withstand drought, we performed large-scale physiological and transcriptome analyses. We selected a mature tree from a semi-arid area with suboptimal growth conditions for clonal propagation through cuttings. We then used a high-throughput experimental system to continuously monitor whole-plant transpiration rates, stomatal conductance and the vapor pressure deficit. The transcriptomes of plants were examined at six physiological stages: pre-stomatal response, partial stomatal closure, minimum transpiration, post-irrigation, partial recovery and full recovery. At each stage, data from plants exposed to the drought treatment were compared with data collected from well-irrigated control plants. A drought-stressed P. halepensis transcriptome was created using paired-end RNA-seq. In total, ~6000 differentially expressed, non-redundant transcripts were identified between drought-treated and control trees. Cluster analysis has revealed stress-induced down-regulation of transcripts related to photosynthesis, reactive oxygen species (ROS)-scavenging through the ascorbic acid (AsA)-glutathione cycle, fatty acid and cell wall biosynthesis, stomatal activity, and the biosynthesis of flavonoids and terpenoids. Up-regulated processes included chlorophyll degradation, ROS-scavenging through AsA-independent thiol-mediated pathways, abscisic acid response and accumulation of heat shock proteins, thaumatin and exordium. Recovery from drought induced strong transcription of retrotransposons, especially the retrovirus-related transposon Tnt1-94. The drought-related transcriptome illustrates this species' dynamic response to drought and recovery and unravels novel mechanisms.
Collapse
Affiliation(s)
- Hagar Fox
- Institute of Plant Sciences, Volcani Center, ARO, Bet Dagan 50250, Israel
- Institute of Plant Sciences and Genetics in Agriculture, The Robert H. Smith Faculty of Agriculture, Food and Environment, The Hebrew University of Jerusalem, Rehovot 76100, Israel
| | | | - Gilor Kelly
- Institute of Plant Sciences, Volcani Center, ARO, Bet Dagan 50250, Israel
| | - Ronny Bourstein
- Institute of Plant Sciences and Genetics in Agriculture, The Robert H. Smith Faculty of Agriculture, Food and Environment, The Hebrew University of Jerusalem, Rehovot 76100, Israel
| | - Ziv Attia
- Institute of Plant Sciences and Genetics in Agriculture, The Robert H. Smith Faculty of Agriculture, Food and Environment, The Hebrew University of Jerusalem, Rehovot 76100, Israel
| | - Jing Zhou
- Institute of Plant Sciences, Volcani Center, ARO, Bet Dagan 50250, Israel
| | - Yosef Moshe
- Institute of Plant Sciences, Volcani Center, ARO, Bet Dagan 50250, Israel
| | - Menachem Moshelion
- Institute of Plant Sciences and Genetics in Agriculture, The Robert H. Smith Faculty of Agriculture, Food and Environment, The Hebrew University of Jerusalem, Rehovot 76100, Israel
| | | |
Collapse
|
26
|
Orgeur M, Martens M, Börno ST, Timmermann B, Duprez D, Stricker S. A dual transcript-discovery approach to improve the delimitation of gene features from RNA-seq data in the chicken model. Biol Open 2018; 7:bio.028498. [PMID: 29183907 PMCID: PMC5827264 DOI: 10.1242/bio.028498] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/20/2023] Open
Abstract
The sequence of the chicken genome, like several other draft genome sequences, is presently not fully covered. Gaps, contigs assigned with low confidence and uncharacterized chromosomes result in gene fragmentation and imprecise gene annotation. Transcript abundance estimation from RNA sequencing (RNA-seq) data relies on read quality, library complexity and expression normalization. In addition, the quality of the genome sequence used to map sequencing reads, and the gene annotation that defines gene features, must also be taken into account. A partially covered genome sequence causes the loss of sequencing reads from the mapping step, while an inaccurate definition of gene features induces imprecise read counts from the assignment step. Both steps can significantly bias interpretation of RNA-seq data. Here, we describe a dual transcript-discovery approach combining a genome-guided gene prediction and a de novo transcriptome assembly. This dual approach enabled us to increase the assignment rate of RNA-seq data by nearly 20% as compared to when using only the chicken reference annotation, contributing therefore to a more accurate estimation of transcript abundance. More generally, this strategy could be applied to any organism with partial genome sequence and/or lacking a manually-curated reference annotation in order to improve the accuracy of gene expression studies.
Collapse
Affiliation(s)
- Mickael Orgeur
- Freie Universität Berlin, Institut für Chemie und Biochemie, Thielallee 63, 14195 Berlin, Germany.,Max Planck Institute for Molecular Genetics, Development and Disease Group, Ihnestrasse 63-73, 14195 Berlin, Germany.,Sorbonne Universités, UPMC Univ. Paris 06, CNRS UMR 7622, Inserm U1156, IBPS-Developmental Biology Laboratory, 9 Quai Saint-Bernard, 75252 Paris Cedex 05, France
| | - Marvin Martens
- Sorbonne Universités, UPMC Univ. Paris 06, CNRS UMR 7622, Inserm U1156, IBPS-Developmental Biology Laboratory, 9 Quai Saint-Bernard, 75252 Paris Cedex 05, France
| | - Stefan T Börno
- Max Planck Institute for Molecular Genetics, Development and Disease Group, Ihnestrasse 63-73, 14195 Berlin, Germany
| | - Bernd Timmermann
- Max Planck Institute for Molecular Genetics, Development and Disease Group, Ihnestrasse 63-73, 14195 Berlin, Germany
| | - Delphine Duprez
- Sorbonne Universités, UPMC Univ. Paris 06, CNRS UMR 7622, Inserm U1156, IBPS-Developmental Biology Laboratory, 9 Quai Saint-Bernard, 75252 Paris Cedex 05, France
| | - Sigmar Stricker
- Freie Universität Berlin, Institut für Chemie und Biochemie, Thielallee 63, 14195 Berlin, Germany .,Max Planck Institute for Molecular Genetics, Development and Disease Group, Ihnestrasse 63-73, 14195 Berlin, Germany
| |
Collapse
|
27
|
Carruthers M, Yurchenko AA, Augley JJ, Adams CE, Herzyk P, Elmer KR. De novo transcriptome assembly, annotation and comparison of four ecological and evolutionary model salmonid fish species. BMC Genomics 2018; 19:32. [PMID: 29310597 PMCID: PMC5759245 DOI: 10.1186/s12864-017-4379-x] [Citation(s) in RCA: 52] [Impact Index Per Article: 8.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/12/2017] [Accepted: 12/11/2017] [Indexed: 12/21/2022] Open
Abstract
Background Salmonid fishes exhibit high levels of phenotypic and ecological variation and are thus ideal model systems for studying evolutionary processes of adaptive divergence and speciation. Furthermore, salmonids are of major interest in fisheries, aquaculture, and conservation research. Improving understanding of the genetic mechanisms underlying traits in these species would significantly progress research in these fields. Here we generate high quality de novo transcriptomes for four salmonid species: Atlantic salmon (Salmo salar), brown trout (Salmo trutta), Arctic charr (Salvelinus alpinus), and European whitefish (Coregonus lavaretus). All species except Atlantic salmon have no reference genome publicly available and few if any genomic studies to date. Results We used paired-end RNA-seq on Illumina to generate high coverage sequencing of multiple individuals, yielding between 180 and 210 M reads per species. After initial assembly, strict filtering was used to remove duplicated, redundant, and low confidence transcripts. The final assemblies consisted of 36,505 protein-coding transcripts for Atlantic salmon, 35,736 for brown trout, 33,126 for Arctic charr, and 33,697 for European whitefish and are made publicly available. Assembly completeness was assessed using three approaches, all of which supported high quality of the assemblies: 1) ~78% of Actinopterygian single-copy orthologs were successfully captured in our assemblies, 2) orthogroup inference identified high overlap in the protein sequences present across all four species (40% shared across all four and 84% shared by at least two), and 3) comparison with the published Atlantic salmon genome suggests that our assemblies represent well covered (~98%) protein-coding transcriptomes. Thorough comparison of the generated assemblies found that 84-90% of transcripts in each assembly were orthologous with at least one of the other three species. We also identified 34-37% of transcripts in each assembly as paralogs. We further compare completeness and annotation statistics of our new assemblies to available related species. Conclusion New, high-confidence protein-coding transcriptomes were generated for four ecologically and economically important species of salmonids. This offers a high quality pipeline for such complex genomes, represents a valuable contribution to the existing genomic resources for these species and provides robust tools for future investigation of gene expression and sequence evolution in these and other salmonid species. Electronic supplementary material The online version of this article (10.1186/s12864-017-4379-x) contains supplementary material, which is available to authorized users.
Collapse
Affiliation(s)
- Madeleine Carruthers
- Institute of Biodiversity, Animal Health & Comparative Medicine, College of Medical, Veterinary & Life Sciences, University of Glasgow, G12 8QQ, Glasgow, UK
| | - Andrey A Yurchenko
- Institute of Biodiversity, Animal Health & Comparative Medicine, College of Medical, Veterinary & Life Sciences, University of Glasgow, G12 8QQ, Glasgow, UK
| | - Julian J Augley
- Glasgow Polyomics, Wolfson Wohl Cancer Research Centre, University of Glasgow, G61 1QH, Glasgow, UK.,Present Address: Fios Genomics Ltd., Nine Edinburgh Bioquarter, 9 Little France Road, Edinburgh, EH16 4UX, UK
| | - Colin E Adams
- Institute of Biodiversity, Animal Health & Comparative Medicine, College of Medical, Veterinary & Life Sciences, University of Glasgow, G12 8QQ, Glasgow, UK.,Scottish Centre for Ecology and the Natural Environment, University of Glasgow, Rowardennan, G63 0AW, UK
| | - Pawel Herzyk
- Glasgow Polyomics, Wolfson Wohl Cancer Research Centre, University of Glasgow, G61 1QH, Glasgow, UK.,Institute of Molecular, Cell & Systems Biology, College of Medical, Veterinary & Life Sciences, University of Glasgow, G12 8QQ, Glasgow, UK
| | - Kathryn R Elmer
- Institute of Biodiversity, Animal Health & Comparative Medicine, College of Medical, Veterinary & Life Sciences, University of Glasgow, G12 8QQ, Glasgow, UK.
| |
Collapse
|
28
|
Berger BA, Han J, Sessa EB, Gardner AG, Shepherd KA, Ricigliano VA, Jabaily RS, Howarth DG. The unexpected depths of genome-skimming data: A case study examining Goodeniaceae floral symmetry genes. APPLICATIONS IN PLANT SCIENCES 2017; 5:apps.1700042. [PMID: 29109919 PMCID: PMC5664964 DOI: 10.3732/apps.1700042] [Citation(s) in RCA: 18] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 04/21/2017] [Accepted: 09/07/2017] [Indexed: 05/20/2023]
Abstract
PREMISE OF THE STUDY The use of genome skimming allows systematists to quickly generate large data sets, particularly of sequences in high abundance (e.g., plastomes); however, researchers may be overlooking data in low abundance that could be used for phylogenetic or evo-devo studies. Here, we present a bioinformatics approach that explores the low-abundance portion of genome-skimming next-generation sequencing libraries in the fan-flowered Goodeniaceae. METHODS Twenty-four previously constructed Goodeniaceae genome-skimming Illumina libraries were examined for their utility in mining low-copy nuclear genes involved in floral symmetry, specifically the CYCLOIDEA (CYC)-like genes. De novo assemblies were generated using multiple assemblers, and BLAST searches were performed for CYC1, CYC2, and CYC3 genes. RESULTS Overall Trinity, SOAPdenovo-Trans, and SOAPdenovo implementing lower k-mer values uncovered the most data, although no assembler consistently outperformed the others. Using SOAPdenovo-Trans across all 24 data sets, we recovered four CYC-like gene groups (CYC1, CYC2, CYC3A, and CYC3B) from a majority of the species. Alignments of the fragments included the entire coding sequence as well as upstream and downstream regions. DISCUSSION Genome-skimming data sets can provide a significant source of low-copy nuclear gene sequence data that may be used for multiple downstream applications.
Collapse
Affiliation(s)
- Brent A. Berger
- Department of Biological Sciences, St. John’s University, 8000 Utopia Parkway, Queens, New York 11439 USA
| | - Jiahong Han
- Department of Biological Sciences, St. John’s University, 8000 Utopia Parkway, Queens, New York 11439 USA
| | - Emily B. Sessa
- Department of Biology, University of Florida, Box 118525, Gainesville, Florida 32611 USA
| | - Andrew G. Gardner
- Department of Biological Sciences, California State University, Stanislaus, One University Circle, Turlock, California 95382 USA
| | - Kelly A. Shepherd
- Western Australian Herbarium, Department of Biodiversity, Conservation and Attractions, 17 Dick Perry Avenue, Kensington 6151, Western Australia, Australia
| | - Vincent A. Ricigliano
- USDA-ARS Carl Hayden Bee Research Center, 2000 E. Allen Road, Tucson, Arizona 85719 USA
| | - Rachel S. Jabaily
- Department of Biology, Rhodes College, 2000 N. Parkway, Memphis, Tennessee 38112 USA
| | - Dianella G. Howarth
- Department of Biological Sciences, St. John’s University, 8000 Utopia Parkway, Queens, New York 11439 USA
| |
Collapse
|
29
|
Berger BA, Han J, Sessa EB, Gardner AG, Shepherd KA, Ricigliano VA, Jabaily RS, Howarth DG. The unexpected depths of genome-skimming data: A case study examining Goodeniaceae floral symmetry genes. APPLICATIONS IN PLANT SCIENCES 2017. [PMID: 29109919 DOI: 10.5061/dryad.0500c] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Subscribe] [Scholar Register] [Indexed: 05/08/2023]
Abstract
PREMISE OF THE STUDY The use of genome skimming allows systematists to quickly generate large data sets, particularly of sequences in high abundance (e.g., plastomes); however, researchers may be overlooking data in low abundance that could be used for phylogenetic or evo-devo studies. Here, we present a bioinformatics approach that explores the low-abundance portion of genome-skimming next-generation sequencing libraries in the fan-flowered Goodeniaceae. METHODS Twenty-four previously constructed Goodeniaceae genome-skimming Illumina libraries were examined for their utility in mining low-copy nuclear genes involved in floral symmetry, specifically the CYCLOIDEA (CYC)-like genes. De novo assemblies were generated using multiple assemblers, and BLAST searches were performed for CYC1, CYC2, and CYC3 genes. RESULTS Overall Trinity, SOAPdenovo-Trans, and SOAPdenovo implementing lower k-mer values uncovered the most data, although no assembler consistently outperformed the others. Using SOAPdenovo-Trans across all 24 data sets, we recovered four CYC-like gene groups (CYC1, CYC2, CYC3A, and CYC3B) from a majority of the species. Alignments of the fragments included the entire coding sequence as well as upstream and downstream regions. DISCUSSION Genome-skimming data sets can provide a significant source of low-copy nuclear gene sequence data that may be used for multiple downstream applications.
Collapse
Affiliation(s)
- Brent A Berger
- Department of Biological Sciences, St. John's University, 8000 Utopia Parkway, Queens, New York 11439 USA
| | - Jiahong Han
- Department of Biological Sciences, St. John's University, 8000 Utopia Parkway, Queens, New York 11439 USA
| | - Emily B Sessa
- Department of Biology, University of Florida, Box 118525, Gainesville, Florida 32611 USA
| | - Andrew G Gardner
- Department of Biological Sciences, California State University, Stanislaus, One University Circle, Turlock, California 95382 USA
| | - Kelly A Shepherd
- Western Australian Herbarium, Department of Biodiversity, Conservation and Attractions, 17 Dick Perry Avenue, Kensington 6151, Western Australia, Australia
| | - Vincent A Ricigliano
- USDA-ARS Carl Hayden Bee Research Center, 2000 E. Allen Road, Tucson, Arizona 85719 USA
| | - Rachel S Jabaily
- Department of Biology, Rhodes College, 2000 N. Parkway, Memphis, Tennessee 38112 USA
| | - Dianella G Howarth
- Department of Biological Sciences, St. John's University, 8000 Utopia Parkway, Queens, New York 11439 USA
| |
Collapse
|
30
|
Bauer E, Schmutzer T, Barilar I, Mascher M, Gundlach H, Martis MM, Twardziok SO, Hackauf B, Gordillo A, Wilde P, Schmidt M, Korzun V, Mayer KFX, Schmid K, Schön CC, Scholz U. Towards a whole-genome sequence for rye (Secale cereale L.). THE PLANT JOURNAL : FOR CELL AND MOLECULAR BIOLOGY 2017; 89:853-869. [PMID: 27888547 DOI: 10.1111/tpj.13436] [Citation(s) in RCA: 120] [Impact Index Per Article: 17.1] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 08/03/2016] [Revised: 11/08/2016] [Accepted: 11/21/2016] [Indexed: 05/18/2023]
Abstract
We report on a whole-genome draft sequence of rye (Secale cereale L.). Rye is a diploid Triticeae species closely related to wheat and barley, and an important crop for food and feed in Central and Eastern Europe. Through whole-genome shotgun sequencing of the 7.9-Gbp genome of the winter rye inbred line Lo7 we obtained a de novo assembly represented by 1.29 million scaffolds covering a total length of 2.8 Gbp. Our reference sequence represents nearly the entire low-copy portion of the rye genome. This genome assembly was used to predict 27 784 rye gene models based on homology to sequenced grass genomes. Through resequencing of 10 rye inbred lines and one accession of the wild relative S. vavilovii, we discovered more than 90 million single nucleotide variants and short insertions/deletions in the rye genome. From these variants, we developed the high-density Rye600k genotyping array with 600 843 markers, which enabled anchoring the sequence contigs along a high-density genetic map and establishing a synteny-based virtual gene order. Genotyping data were used to characterize the diversity of rye breeding pools and genetic resources, and to obtain a genome-wide map of selection signals differentiating the divergent gene pools. This rye whole-genome sequence closes a gap in Triticeae genome research, and will be highly valuable for comparative genomics, functional studies and genome-based breeding in rye.
Collapse
Affiliation(s)
- Eva Bauer
- Technical University of Munich, Plant Breeding, Liesel-Beckmann-Str. 2, 85354, Freising, Germany
| | - Thomas Schmutzer
- Leibniz Institute of Plant Genetics and Crop Plant Research (IPK) Gatersleben, Corrensstr. 3, 06466, Stadt Seeland, Germany
| | - Ivan Barilar
- Universität Hohenheim, Crop Biodiversity and Breeding Informatics, Fruwirthstr. 21, 70599, Stuttgart, Germany
| | - Martin Mascher
- Leibniz Institute of Plant Genetics and Crop Plant Research (IPK) Gatersleben, Corrensstr. 3, 06466, Stadt Seeland, Germany
| | - Heidrun Gundlach
- Helmholtz Zentrum München, Plant Genome and Systems Biology, Ingolstädter Landstraße 1, 85764, Neuherberg, Germany
| | - Mihaela M Martis
- Helmholtz Zentrum München, Plant Genome and Systems Biology, Ingolstädter Landstraße 1, 85764, Neuherberg, Germany
| | - Sven O Twardziok
- Helmholtz Zentrum München, Plant Genome and Systems Biology, Ingolstädter Landstraße 1, 85764, Neuherberg, Germany
| | - Bernd Hackauf
- Julius Kühn-Institute, Institute for Breeding Research on Agricultural Crops, Rudolf-Schick-Platz 3a, 18190, Sanitz, Germany
| | - Andres Gordillo
- KWS LOCHOW GMBH, Ferdinand-von-Lochow-Str. 5, 29303, Bergen, Germany
| | - Peer Wilde
- KWS LOCHOW GMBH, Ferdinand-von-Lochow-Str. 5, 29303, Bergen, Germany
| | - Malthe Schmidt
- KWS LOCHOW GMBH, Ferdinand-von-Lochow-Str. 5, 29303, Bergen, Germany
| | - Viktor Korzun
- KWS LOCHOW GMBH, Ferdinand-von-Lochow-Str. 5, 29303, Bergen, Germany
| | - Klaus F X Mayer
- Helmholtz Zentrum München, Plant Genome and Systems Biology, Ingolstädter Landstraße 1, 85764, Neuherberg, Germany
| | - Karl Schmid
- Universität Hohenheim, Crop Biodiversity and Breeding Informatics, Fruwirthstr. 21, 70599, Stuttgart, Germany
| | - Chris-Carolin Schön
- Technical University of Munich, Plant Breeding, Liesel-Beckmann-Str. 2, 85354, Freising, Germany
| | - Uwe Scholz
- Leibniz Institute of Plant Genetics and Crop Plant Research (IPK) Gatersleben, Corrensstr. 3, 06466, Stadt Seeland, Germany
| |
Collapse
|
31
|
Gonzalez-Ibeas D, Martinez-Garcia PJ, Famula RA, Delfino-Mix A, Stevens KA, Loopstra CA, Langley CH, Neale DB, Wegrzyn JL. Assessing the Gene Content of the Megagenome: Sugar Pine (Pinus lambertiana). G3 (BETHESDA, MD.) 2016; 6:3787-3802. [PMID: 27799338 PMCID: PMC5144951 DOI: 10.1534/g3.116.032805] [Citation(s) in RCA: 29] [Impact Index Per Article: 3.6] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 06/24/2016] [Accepted: 07/13/2016] [Indexed: 02/06/2023]
Abstract
Sugar pine (Pinus lambertiana Douglas) is within the subgenus Strobus with an estimated genome size of 31 Gbp. Transcriptomic resources are of particular interest in conifers due to the challenges presented in their megagenomes for gene identification. In this study, we present the first comprehensive survey of the P. lambertiana transcriptome through deep sequencing of a variety of tissue types to generate more than 2.5 billion short reads. Third generation, long reads generated through PacBio Iso-Seq have been included for the first time in conifers to combat the challenges associated with de novo transcriptome assembly. A technology comparison is provided here to contribute to the otherwise scarce comparisons of second and third generation transcriptome sequencing approaches in plant species. In addition, the transcriptome reference was essential for gene model identification and quality assessment in the parallel project responsible for sequencing and assembly of the entire genome. In this study, the transcriptomic data were also used to address questions surrounding lineage-specific Dicer-like proteins in conifers. These proteins play a role in the control of transposable element proliferation and the related genome expansion in conifers.
Collapse
Affiliation(s)
- Daniel Gonzalez-Ibeas
- Department of Ecology and Evolutionary Biology, University of Connecticut, Storrs, Connecticut 06269
| | | | - Randi A Famula
- Department of Plant Sciences, University of California, Davis, California 95616
| | - Annette Delfino-Mix
- United States Department of Agriculture Forest Service, Institute of Forest Genetics, Placerville, California 95667
| | - Kristian A Stevens
- Department of Evolution and Ecology, University of California, Davis, California 95616
| | - Carol A Loopstra
- Department of Ecosystem Science and Management, Texas A&M University, College Station, Texas 77843
| | - Charles H Langley
- Department of Evolution and Ecology, University of California, Davis, California 95616
| | - David B Neale
- Department of Plant Sciences, University of California, Davis, California 95616
| | - Jill L Wegrzyn
- Department of Ecology and Evolutionary Biology, University of Connecticut, Storrs, Connecticut 06269
| |
Collapse
|
32
|
Rai A, Nakamura M, Takahashi H, Suzuki H, Saito K, Yamazaki M. High-throughput sequencing and de novo transcriptome assembly of Swertia japonica to identify genes involved in the biosynthesis of therapeutic metabolites. PLANT CELL REPORTS 2016; 35:2091-111. [PMID: 27378356 DOI: 10.1007/s00299-016-2021-z] [Citation(s) in RCA: 24] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 04/26/2016] [Accepted: 06/29/2016] [Indexed: 05/08/2023]
Abstract
Here, we report potential transcripts involved in the biosynthesis of therapeutic metabolites in Swertia japonica , the first report of transcriptome assembly, and characterization of the medicinal plant from Swertia genus. Swertia genus, representing over 170 plant species including herbs such as S. chirata, S. hookeri, S. longifolia, S. japonica, among others, have been used as the traditional medicine in China, India, Korea, and Japan for thousands of years. Due to the lack of genomic and transcriptomic resources, little is known about the molecular basis involved in the biosynthesis of characteristic key bioactive metabolites. Here, we performed deep-transcriptome sequencing for the aerial tissues and the roots of S. japonica, generating over 2 billion raw reads with an average length of 101 bps. Using a combined approach of three popular assemblers, de novo transcriptome assembly for S. japonica was obtained, yielding 81,729 unigenes having an average length of 884 bps and N50 value of 1452 bps, of which 46,963 unigenes were annotated based on the sequence similarity against NCBI-nr protein database. Annotation of transcriptome assembly resulted in the identification of putative genes encoding all enzymes from the key therapeutic metabolite biosynthesis pathways. Transcript abundance analysis, gene ontology enrichment analysis, and KEGG pathway enrichment analysis revealed metabolic processes being up-regulated in the aerial tissues with respect to the roots of S. japonica. We also identified 37 unigenes as potential candidates involved in the glycosylation of bioactive metabolites. Being the first report of transcriptome assembly and annotation for any of the Swertia species, this study will be a valuable resource for future investigations on the biosynthetic pathways of therapeutic metabolites and their regulations.
Collapse
Affiliation(s)
- Amit Rai
- Graduate School of Pharmaceutical Sciences, Chiba University, Inohana 1-8-1, Chuo-ku, Chiba, 260-8675, Japan.
| | - Michimi Nakamura
- Graduate School of Pharmaceutical Sciences, Chiba University, Inohana 1-8-1, Chuo-ku, Chiba, 260-8675, Japan
| | - Hiroki Takahashi
- Medical Mycology Research Center, Chiba University, Chiba, Japan
| | - Hideyuki Suzuki
- Department of Research and Development, Kazusa DNA Research Institute, Chiba, Japan
| | - Kazuki Saito
- Graduate School of Pharmaceutical Sciences, Chiba University, Inohana 1-8-1, Chuo-ku, Chiba, 260-8675, Japan
| | - Mami Yamazaki
- Graduate School of Pharmaceutical Sciences, Chiba University, Inohana 1-8-1, Chuo-ku, Chiba, 260-8675, Japan.
| |
Collapse
|
33
|
Rai A, Kamochi H, Suzuki H, Nakamura M, Takahashi H, Hatada T, Saito K, Yamazaki M. De novo transcriptome assembly and characterization of nine tissues of Lonicera japonica to identify potential candidate genes involved in chlorogenic acid, luteolosides, and secoiridoid biosynthesis pathways. J Nat Med 2016; 71:1-15. [PMID: 27629269 PMCID: PMC5214891 DOI: 10.1007/s11418-016-1041-x] [Citation(s) in RCA: 34] [Impact Index Per Article: 4.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/07/2016] [Accepted: 09/01/2016] [Indexed: 11/18/2022]
Abstract
Lonicera japonica is one of the most important medicinal plants with applications in traditional Chinese and Japanese medicine for thousands of years. Extensive studies on the constituents of L. japonica extracts have revealed an accumulation of pharmaceutically active metabolite classes, such as chlorogenic acid, luteolin and other flavonoids, and secoiridoids, which impart characteristic medicinal properties. Despite being a rich source of pharmaceutically active metabolites, little is known about the biosynthetic enzymes involved, and their expression profile across different tissues of L. japonica. In this study, we performed de novo transcriptome assembly for L. japonica, representing transcripts from nine different tissues. A total of 22 Gbps clean RNA-seq reads from nine tissues of L. japonica were used, resulting in 243,185 unigenes, with 99,938 unigenes annotated based on a homology search using blastx against the NCBI-nr protein database. Unsupervised principal component analysis and correlation studies using transcript expression data from all nine tissues of L. japonica showed relationships between tissues, explaining their association at different developmental stages. Homologs for all genes associated with chlorogenic acid, luteolin, and secoiridoid biosynthesis pathways were identified in the L. japonica transcriptome assembly. Expression of unigenes associated with chlorogenic acid was enriched in stems and leaf-2, unigenes from luteolin were enriched in stems and flowers, while unigenes from secoiridoid metabolic pathways were enriched in leaf-1 and shoot apex. Our results showed that different tissues of L. japonica are enriched with sets of unigenes associated with specific pharmaceutically important metabolic pathways and, therefore, possess unique medicinal properties. The present study will serve as a resource for future attempts for functional characterization of enzyme coding genes within key metabolic processes.
Collapse
Affiliation(s)
- Amit Rai
- Department of Molecular Biology and Biotechnology, Graduate School of Pharmaceutical Sciences, Chiba University, Inohana 1-8-1, Chuo-ku, Chiba, 260-8675, Japan.
| | - Hidetaka Kamochi
- Department of Molecular Biology and Biotechnology, Graduate School of Pharmaceutical Sciences, Chiba University, Inohana 1-8-1, Chuo-ku, Chiba, 260-8675, Japan
| | | | - Michimi Nakamura
- Department of Molecular Biology and Biotechnology, Graduate School of Pharmaceutical Sciences, Chiba University, Inohana 1-8-1, Chuo-ku, Chiba, 260-8675, Japan
| | - Hiroki Takahashi
- Medical Mycology Research Center, Chiba University, Chiba, Japan
| | - Tomoki Hatada
- Department of Molecular Biology and Biotechnology, Graduate School of Pharmaceutical Sciences, Chiba University, Inohana 1-8-1, Chuo-ku, Chiba, 260-8675, Japan
| | - Kazuki Saito
- Department of Molecular Biology and Biotechnology, Graduate School of Pharmaceutical Sciences, Chiba University, Inohana 1-8-1, Chuo-ku, Chiba, 260-8675, Japan
| | - Mami Yamazaki
- Department of Molecular Biology and Biotechnology, Graduate School of Pharmaceutical Sciences, Chiba University, Inohana 1-8-1, Chuo-ku, Chiba, 260-8675, Japan.
| |
Collapse
|