1
|
Khelghatibana F, Javan-Nikkhah M, Safaie N, Sobhani A, Shams S, Sari E. A reference transcriptome for walnut anthracnose pathogen, Ophiognomonia leptostyla, guides the discovery of candidate virulence genes. Fungal Genet Biol 2023; 169:103828. [PMID: 37657751 DOI: 10.1016/j.fgb.2023.103828] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/12/2023] [Revised: 08/13/2023] [Accepted: 08/28/2023] [Indexed: 09/03/2023]
Abstract
Despite the economic losses due to the walnut anthracnose, Ophiognomonia leptostyla is an orphan fungus with respect to genomic resources. In the present study, the transcriptome of O. leptostyla was assembled for the first time. RNA sequencing was conducted for the fungal mycelia grown in a liquid media, and the inoculated leaf samples of walnut with the fungal conidia sampled at 48, 96 and 144 h post inoculation (hpi). The completeness, correctness, and contiguity of the de novo transcriptome assemblies generated with Trinity, Oases, SOAPdenovo-Trans and Bridger were compared to identify a single superior reference assembly. In most of the assessment criteria including N50, Transrate score, number of ORFs with known description in gene bank, the percentage of reads mapped back to the transcript (RMBT), BUSCO score, Swiss-Prot coverage bin and RESM-EVAL score, the Bridger assembly was the superior and thus used as a reference for profiling the O. leptostyla transcriptome in liquid media vs. during walnut infection. The k-means clustering of transcripts resulted in four distinct transcription patterns across the three sampling time points. Most of the detected CAZy transcripts had elevated transcription at 96 hpi that is hypothetically concurrent with the start of intracellular growth. The in-silico analysis revealed 103 candidate effectors of which six were members of Necrosis and Ethylene Inducing Like Protein (NLP) gene family belonging to three distinct k-means clusters. This study provided a complex and temporal pattern of the CAZys and candidate effectors transcription during six days post O. leptostyla inoculation on walnut leaves, introducing a list of candidate virulence genes for validation in future studies.
Collapse
Affiliation(s)
- Fatemeh Khelghatibana
- Department of Plant Pathology, Iranian Research Institute of Plant Protection, Agricultural Research, Education and Extension Organization (AREEO), Tehran, Iran.
| | - Mohammad Javan-Nikkhah
- Department of Plant Protection, College of Agriculture and Natural Resources, University of Tehran, Karaj, Iran
| | - Naser Safaie
- Department of Plant Pathology, Faculty of Agriculture, Tarbiat Modares University, Tehran, Iran
| | - Ahmad Sobhani
- Agricultural Biotechnology Research Institute of Iran - Isfahan Branch, Agricultural Research, Education and Extension Organization (AREEO), Isfahan, Iran
| | - Somayeh Shams
- Department of Plant Production and Genetic Engineering, Faculty of Agriculture, University of Lorestan, Khorramabad, Iran
| | - Ehsan Sari
- Department of Microbiology and Plant Pathology, University of California, Riverside, CA, USA.
| |
Collapse
|
2
|
Fan S, Xu Y, Bai M, Luo F, Yu J, Yang G. Integrated Transcriptome and Metabolome Analysis Revealed the Causal Agent of Primary Bud Necrosis in 'Summer Black' Grape. Int J Mol Sci 2023; 24:10410. [PMID: 37373557 DOI: 10.3390/ijms241210410] [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/25/2023] [Revised: 06/13/2023] [Accepted: 06/16/2023] [Indexed: 06/29/2023] Open
Abstract
Primary bud necrosis of grape buds is a physiological disorder that leads to decreased berry yield and has a catastrophic impact on the double cropping system in sub-tropical areas. The pathogenic mechanisms and potential solutions remain unknown. In this study, the progression and irreversibility patterns of primary bud necrosis in 'Summer Black' were examined via staining and transmission electron microscopy observation. Primary bud necrosis was initiated at 60 days after bud break and was characterized by plasmolysis, mitochondrial swelling, and severe damage to other organelles. To reveal the underlying regulatory networks, winter buds were collected during primary bud necrosis progression for integrated transcriptome and metabolome analysis. The accumulation of reactive oxygen species and subsequent signaling cascades disrupted the regulation systems for cellular protein quality. ROS cascade reactions were related to mitochondrial stress that can lead to mitochondrial dysfunction, lipid peroxidation causing damage to membrane structure, and endoplasmic reticulum stress leading to misfolded protein aggregates. All these factors ultimately resulted in primary bud necrosis. Visible tissue browning was associated with the oxidation and decreased levels of flavonoids during primary bud necrosis, while the products of polyunsaturated fatty acids and stilbenes exhibited an increasing trend, leading to a shift in carbon flow from flavonoids to stilbene. Increased ethylene may be closely related to primary bud necrosis, while auxin accelerated cell growth and alleviated necrosis by co-chaperone VvP23-regulated redistribution of auxin in meristem cells. Altogether, this study provides important clues for further study on primary bud necrosis.
Collapse
Affiliation(s)
- Shaogang Fan
- College of Horticulture, Hunan Agricultural University, Changsha 410128, China
| | - Yanshuai Xu
- College of Horticulture, Hunan Agricultural University, Changsha 410128, China
| | - Miao Bai
- College of Horticulture, Hunan Agricultural University, Changsha 410128, China
| | - Feixiong Luo
- College of Horticulture, Hunan Agricultural University, Changsha 410128, China
| | - Jun Yu
- College of Horticulture, Hunan Agricultural University, Changsha 410128, China
| | - Guoshun Yang
- College of Horticulture, Hunan Agricultural University, Changsha 410128, China
| |
Collapse
|
3
|
Birt HWG, Pattison AB, Skarshewski A, Daniells J, Raghavendra A, Dennis PG. The core fungal microbiome of banana (Musa spp.). Front Microbiol 2023; 14:1127779. [PMID: 37065131 PMCID: PMC10098452 DOI: 10.3389/fmicb.2023.1127779] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/20/2022] [Accepted: 03/08/2023] [Indexed: 04/03/2023] Open
Abstract
Here, we report a metabarcoding (ITS2) study to define the common core fungal microbiome (mycobiome) of healthy Musa spp. (bananas and plantains). To identify a list of 21 core fungal taxa, we first characterised the effects of edaphic conditions and host genotype – two factors that are likely to differ between farms – on the diversity of fungal communities in bulk soil and seven plant compartments. This experiment facilitated shortlisting of core ‘candidates’, which were then elevated to full core status if also found to frequent a wide-range of field-grown Musa spp. and exhibit hub-like characteristics in network analyses. Subsequently, we conducted a meta-analysis of eleven publicly available datasets of Musa spp. associated fungi demonstrating that the core fungi identified in our study have close relatives in other countries. The diversity and composition of mycobiomes differed between plant compartments and soils, but not genotypes. The core mycobiome included Fusarium oxysporum and its relatives, which dominated all plant compartments, as well as members of the Sordariomycetes, Dothideomycetes, and Mortierellomycota. Our study provides a robust list of common core fungal taxa for Musa spp. Further studies may consider how changes in the frequencies and activities of these taxa influence host fitness and whether they can be managed to improve banana production.
Collapse
Affiliation(s)
- Henry W. G. Birt
- School of Earth and Environmental Sciences, The University of Queensland, Brisbane, QLD, Australia
| | - Anthony B. Pattison
- Department of Agriculture and Fisheries, Centre for Wet Tropics Agriculture, South Johnstone, QLD, Australia
| | - Adam Skarshewski
- School of Earth and Environmental Sciences, The University of Queensland, Brisbane, QLD, Australia
| | - Jeff Daniells
- Department of Agriculture and Fisheries, Centre for Wet Tropics Agriculture, South Johnstone, QLD, Australia
| | - Anil Raghavendra
- School of Earth and Environmental Sciences, The University of Queensland, Brisbane, QLD, Australia
| | - Paul G. Dennis
- School of Earth and Environmental Sciences, The University of Queensland, Brisbane, QLD, Australia
- *Correspondence: Paul G. Dennis,
| |
Collapse
|
4
|
Samal I, Bhoi TK, Majhi PK, Murmu S, Pradhan AK, Kumar D, Saini V, Paschapur AU, Raj MN, Ankur, Manik S, Behera PP, Mahanta DK, Komal J, Alam P, Balawi TA. Combatting insects mediated biotic stress through plant associated endophytic entomopathogenic fungi in horticultural crops. FRONTIERS IN PLANT SCIENCE 2023; 13:1098673. [PMID: 36743574 PMCID: PMC9894630 DOI: 10.3389/fpls.2022.1098673] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 11/15/2022] [Accepted: 12/16/2022] [Indexed: 09/12/2023]
Abstract
Horticultural production is a vital catalyst for economic growth, yet insect infestations reduce horticultural crop yield and quality. Pesticides and other pest control methods are used during planting to eliminate pests that cause direct and indirect losses. In such situations, endophytic entomo-pathogenic fungi (EEPF) can act as a potential tools for biological control. They protect plants by boosting growth, nutrition, morpho-physiology and salt or iron tolerance. Antixenosis, antibiosis and plant tolerance change insect performance and preferences. EEPF- plant colonisation slows herbivore development, food consumption, oviposition and larval survival. EEPF changes plant physio-chemical properties like volatile emission profile and secondary metabolite production to regulate insect pest defences. EEPF produces chitinases, laccases, amylases, and cellulases for plant defence. Recent studies focused on EEPF species' significance, isolation, identification and field application. Realizing their full potential is difficult due to insufficient mass production, storage stability and formulation. Genetic-molecular and bioinformatics can help to build EEPF-based biological control systems. Metagenomics helps study microbial EEPF taxonomy and function. Multi-omics and system biology can decode EEPF interactions with host plants and microorganisms. NGS (Next Generation Sequencing), comparative genomics, proteomics, transcriptomics, metabolomics, metatranscriptomics and microarrays are used to evaluate plant-EEPF relationships. IPM requires understanding the abiotic and biotic elements that influence plant-EEPF interaction and the physiological mechanisms of EEPF colonisation. Due to restricted research, there are hundreds of unexplored EEPFs, providing an urgent need to uncover and analyse them.
Collapse
Affiliation(s)
- Ipsita Samal
- Department of Entomology, Sri Sri University, Cuttack, Odisha, India
| | - Tanmaya Kumar Bhoi
- Forest Protection Division, Indian Council of Forestry Research and Education (ICFRE) - Arid Forest Research Institute (AFRI), Jodhpur, Rajasthan, India
| | - Prasanta Kumar Majhi
- Department of Plant Breeding and Genetics, Odisha University of Agriculture and Technology, Bhubaneswar, Odisha, India
| | - Sneha Murmu
- Division of Agricultural Bio-informatics, Indian Council of Agricultural Research (ICAR)- Indian Agricultural Statistics Research Institute, New Delhi, India
| | - Asit Kumar Pradhan
- Division, Social Science Division, Indian Council of Agricultural Research (ICAR)- National Rice Research Institute (NRRI), Cuttack, Odisha, India
| | - Dilip Kumar
- Division of Computer Application and IT, National Institute for Agricultural Economics and Policy Research (NIAP), New Delhi, National Capital Territory of Delhi, India
| | - Varun Saini
- Department of Entomology, Navsari Agricultural University, Navsari, Gujarat, India
| | - Amit Umesh Paschapur
- Crop Protection Division, Indian Council of Agricultural Research (ICAR) - Vivekananda Parvatiya Krishi Anusandhan Sansthan, Almora, Uttarakhand, India
| | - M Nikhil Raj
- Department of Entomology, Dr. Rajendra Prasad Central Agricultural University, Samastipur, Bihar, India
| | - Ankur
- Division of Entomology, Indian Council of Agricultural Research (ICAR-IARI)- Indian Agricultural Research Institute, New Delhi, India
| | - Suryakant Manik
- Department of Seed Science and Technology, Chaudhary Charan Singh Haryana Agricultural University, Hisar, Haryana, India
| | - Partha Pratim Behera
- Department of Plant Breeding and Genetics, Assam Agricultural University, Jorhat, Assam, India
| | - Deepak Kumar Mahanta
- Department of Entomology, Dr. Rajendra Prasad Central Agricultural University, Samastipur, Bihar, India
| | - J. Komal
- Department of Entomology, Navsari Agricultural University, Navsari, Gujarat, India
| | - Pravej Alam
- Department of Biology, College of Science and Humanities in Al-Kharj, Prince Sattam Bin Abdulaziz University, Al Kharj, Saudi Arabia
| | - Thamer Al Balawi
- Department of Biology, College of Science and Humanities, Prince Sattam Bin Abdulaziz University, Al Kharj, Saudi Arabia
| |
Collapse
|
5
|
High Variability of Fungal Communities Associated with the Functional Tissues and Rhizosphere Soil of Picea abies in the Southern Baltics. FORESTS 2022. [DOI: 10.3390/f13071103] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/16/2022]
Abstract
Climate change, which leads to higher temperatures, droughts, and storms, is expected to have a strong effect on both health of forest trees and associated biodiversity. The aim of this study was to investigate the diversity and composition of fungal communities associated with the functional tissues and rhizosphere soil of healthy-looking Picea abies to better understand these fungal communities and their potential effect on tree health in the process of climate change. The study sites included 30 P. abies stands, where needles, shoots, roots, and the rhizosphere soil was sampled. DNA was isolated from individual samples, amplified using ITS2 rRNA as a marker and subjected to high-throughput sequencing. The sequence analysis showed the presence of 232,547 high-quality reads, which following clustering were found to represent 2701 non-singleton fungal OTUs. The highest absolute richness of fungal OTUs was in the soil (1895), then in the needles (1049) and shoots (1002), and the lowest was in the roots (641). The overall fungal community was composed of Ascomycota (58.3%), Basidiomycota (37.2%), Zygomycota (2.5%), Chytridiomycota (1.6%), and Glomeromycota (0.4%). The most common fungi based on sequence read abundance were Aspergillus pseudoglaucus (7.9%), Archaeorhizomyces sp. (3.6%), and Rhinocladiella sp. (2.0%). Pathogens were relatively rare, among which the most common were Phacidium lacerum (1.7%), Cyphellophora sessilis (1.4%), and Rhizosphaera kalkhoffii (1.4%). The results showed that the detected diversity of fungal OTUs was generally high, but their relative abundance varied greatly among different study sites, thereby highlighting the complexity of interactions among the host trees, fungi, and local environmental conditions.
Collapse
|
6
|
Meshram S, Gogoi R, Bashyal BM, Kumar A, Mandal PK, Hossain F. Comparative Transcriptome Analysis of Fungal Pathogen Bipolaris maydis to Understand Pathogenicity Behavior on Resistant and Susceptible Non-CMS Maize Genotypes. Front Microbiol 2022; 13:837056. [PMID: 35572625 PMCID: PMC9100685 DOI: 10.3389/fmicb.2022.837056] [Citation(s) in RCA: 5] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/16/2021] [Accepted: 03/02/2022] [Indexed: 11/13/2022] Open
Abstract
Bipolaris maydis is pathogen of maize which causes maydis leaf blight disease. In India major losses occur due to the B. maydis race “O” pathogen, whereas in other parts of the world, major losses are due to the race “T” pathogen. In the present study, we conducted an in planta transcriptomics study of the B. maydis race “O” pathogen after infection on non-CMS maize resistant and susceptible genotypes by mRNA sequencing to understand the molecular basis of pathogenicity for better management of the pathogen. Approximately 23.4 GB of mRNA-seq data of B. maydis were obtained from both resistant and susceptible maize backgrounds for fungus. Differentially expressed genes (DEGs) analysis of B. maydis in two different genetic backgrounds suggested that the majority of highly DEGs were associated with mitochondrial, cell wall and chitin synthesis, sugar metabolism, peroxidase activity, mitogen-activated protein kinase (MAPK) activity, and shikimate dehydrogenase. KEGG analysis showed that the biosynthetic pathways for secondary metabolism, antibiotics, and carbon metabolism of fungus were highly enriched, respectively, in susceptible backgrounds during infection. Previous studies in other host pathogen systems suggest that these genes play a vital role in causing disease in their host plants. Our study is probably the first transcriptome study of the B. maydis race “O” pathogen and provides in-depth insight of pathogenicity on the host.
Collapse
Affiliation(s)
- Shweta Meshram
- Division of Plant Pathology, Indian Council of Agricultural Research (ICAR)-Indian Agricultural Research Institute, New Delhi, India
| | - Robin Gogoi
- Division of Plant Pathology, Indian Council of Agricultural Research (ICAR)-Indian Agricultural Research Institute, New Delhi, India
- *Correspondence: Robin Gogoi,
| | - Bishnu Maya Bashyal
- Division of Plant Pathology, Indian Council of Agricultural Research (ICAR)-Indian Agricultural Research Institute, New Delhi, India
| | - Aundy Kumar
- Division of Plant Pathology, Indian Council of Agricultural Research (ICAR)-Indian Agricultural Research Institute, New Delhi, India
| | - Pranab Kumar Mandal
- Indian Council of Agricultural Research (ICAR)-National Institute for Plant Biotechnology, New Delhi, India
| | - Firoz Hossain
- Division of Genetics, Indian Council of Agricultural Research (ICAR)-Indian Agricultural Research Institute, New Delhi, India
| |
Collapse
|
7
|
Aluko OO, Li C, Wang Q, Liu H. Sucrose Utilization for Improved Crop Yields: A Review Article. Int J Mol Sci 2021; 22:4704. [PMID: 33946791 PMCID: PMC8124652 DOI: 10.3390/ijms22094704] [Citation(s) in RCA: 39] [Impact Index Per Article: 13.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/02/2021] [Revised: 04/14/2021] [Accepted: 04/17/2021] [Indexed: 12/13/2022] Open
Abstract
Photosynthetic carbon converted to sucrose is vital for plant growth. Sucrose acts as a signaling molecule and a primary energy source that coordinates the source and sink development. Alteration in source-sink balance halts the physiological and developmental processes of plants, since plant growth is mostly triggered when the primary assimilates in the source leaf balance with the metabolic needs of the heterotrophic sinks. To measure up with the sink organ's metabolic needs, the improvement of photosynthetic carbon to synthesis sucrose, its remobilization, and utilization at the sink level becomes imperative. However, environmental cues that influence sucrose balance within these plant organs, limiting positive yield prospects, have also been a rising issue over the past few decades. Thus, this review discusses strategies to improve photosynthetic carbon assimilation, the pathways actively involved in the transport of sucrose from source to sink organs, and their utilization at the sink organ. We further emphasize the impact of various environmental cues on sucrose transport and utilization, and the strategic yield improvement approaches under such conditions.
Collapse
Affiliation(s)
- Oluwaseun Olayemi Aluko
- Tobacco Research Institute of Chinese Academy of Agricultural Sciences, Qingdao 266101, China; (O.O.A.); (C.L.)
- Graduate School of Chinese Academy of Agricultural Sciences, Beijing 100081, China
| | - Chuanzong Li
- Tobacco Research Institute of Chinese Academy of Agricultural Sciences, Qingdao 266101, China; (O.O.A.); (C.L.)
- Graduate School of Chinese Academy of Agricultural Sciences, Beijing 100081, China
| | - Qian Wang
- Tobacco Research Institute of Chinese Academy of Agricultural Sciences, Qingdao 266101, China; (O.O.A.); (C.L.)
| | - Haobao Liu
- Tobacco Research Institute of Chinese Academy of Agricultural Sciences, Qingdao 266101, China; (O.O.A.); (C.L.)
| |
Collapse
|
8
|
Hill R, Llewellyn T, Downes E, Oddy J, MacIntosh C, Kallow S, Panis B, Dickie JB, Gaya E. Seed Banks as Incidental Fungi Banks: Fungal Endophyte Diversity in Stored Seeds of Banana Wild Relatives. Front Microbiol 2021; 12:643731. [PMID: 33841366 PMCID: PMC8024981 DOI: 10.3389/fmicb.2021.643731] [Citation(s) in RCA: 6] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/18/2020] [Accepted: 02/19/2021] [Indexed: 01/19/2023] Open
Abstract
Seed banks were first established to conserve crop genetic diversity, but seed banking has more recently been extended to wild plants, particularly crop wild relatives (CWRs) (e.g., by the Millennium Seed Bank (MSB), Royal Botanic Gardens Kew). CWRs have been recognised as potential reservoirs of beneficial traits for our domesticated crops, and with mounting evidence of the importance of the microbiome to organismal health, it follows that the microbial communities of wild relatives could also be a valuable resource for crop resilience to environmental and pathogenic threats. Endophytic fungi reside asymptomatically inside all plant tissues and have been found to confer advantages to their plant host. Preserving the natural microbial diversity of plants could therefore represent an important secondary conservation role of seed banks. At the same time, species that are reported as endophytes may also be latent pathogens. We explored the potential of the MSB as an incidental fungal endophyte bank by assessing diversity of fungi inside stored seeds. Using banana CWRs in the genus Musa as a case-study, we sequenced an extended ITS-LSU fragment in order to delimit operational taxonomic units (OTUs) and used a similarity and phylogenetics approach for classification. Fungi were successfully detected inside just under one third of the seeds, with a few genera accounting for most of the OTUs-primarily Lasiodiplodia, Fusarium, and Aspergillus-while a large variety of rare OTUs from across the Ascomycota were isolated only once. Fusarium species were notably abundant-of significance in light of Fusarium wilt, a disease threatening global banana crops-and so were targeted for additional sequencing with the marker EF1α in order to delimit species and place them in a phylogeny of the genus. Endophyte community composition, diversity and abundance was significantly different across habitats, and we explored the relationship between community differences and seed germination/viability. Our results show that there is a previously neglected invisible fungal dimension to seed banking that could well have implications for the seed collection and storage procedures, and that collections such as the MSB are indeed a novel source of potentially useful fungal strains.
Collapse
Affiliation(s)
- Rowena Hill
- Department of Comparative Plant and Fungal Biology, Royal Botanic Gardens, Kew, Richmond, United Kingdom
- School of Biological and Chemical Sciences, Faculty of Science and Engineering, Queen Mary University of London, London, United Kingdom
| | - Theo Llewellyn
- Department of Comparative Plant and Fungal Biology, Royal Botanic Gardens, Kew, Richmond, United Kingdom
- Department of Life Sciences, Faculty of Natural Sciences, Imperial College London, London, United Kingdom
| | - Elizabeth Downes
- Department for Environment, Food and Rural Affairs, London, United Kingdom
| | - Joseph Oddy
- Department of Plant Science, Rothamsted Research, Harpenden, United Kingdom
| | - Catriona MacIntosh
- Department of Comparative Plant and Fungal Biology, Royal Botanic Gardens, Kew, Richmond, United Kingdom
- School of Life Sciences, University of Glasgow, Glasgow, United Kingdom
| | - Simon Kallow
- Collections Department, Royal Botanic Gardens, Kew, Millennium Seed Bank, Ardingly, United Kingdom
- Division of Crop Biotechnics, Department of Biosystems, Faculty of Bioscience Engineering, University of Leuven, Leuven, Belgium
| | - Bart Panis
- Bioversity International, Montpellier, France
| | - John B. Dickie
- Collections Department, Royal Botanic Gardens, Kew, Millennium Seed Bank, Ardingly, United Kingdom
| | - Ester Gaya
- Department of Comparative Plant and Fungal Biology, Royal Botanic Gardens, Kew, Richmond, United Kingdom
| |
Collapse
|