1
|
Silva SR, Miranda VFO, Michael TP, Płachno BJ, Matos RG, Adamec L, Pond SLK, Lucaci AG, Pinheiro DG, Varani AM. The phylogenomics and evolutionary dynamics of the organellar genomes in carnivorous Utricularia and Genlisea species (Lentibulariaceae). Mol Phylogenet Evol 2023; 181:107711. [PMID: 36693533 DOI: 10.1016/j.ympev.2023.107711] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/03/2022] [Revised: 01/13/2023] [Accepted: 01/18/2023] [Indexed: 01/22/2023]
Abstract
Utricularia and Genlisea are highly specialized carnivorous plants whose phylogenetic history has been poorly explored using phylogenomic methods. Additional sampling and genomic data are needed to advance our phylogenetic and taxonomic knowledge of this group of plants. Within a comparative framework, we present a characterization of plastome (PT) and mitochondrial (MT) genes of 26 Utricularia and six Genlisea species, with representatives of all subgenera and growth habits. All PT genomes maintain similar gene content, showing minor variation across the genes located between the PT junctions. One exception is a major variation related to different patterns in the presence and absence of ndh genes in the small single copy region, which appears to follow the phylogenetic history of the species rather than their lifestyle. All MT genomes exhibit similar gene content, with most differences related to a lineage-specific pseudogenes. We find evidence for episodic positive diversifying selection in PT and for most of the Utricularia MT genes that may be related to the current hypothesis that bladderworts' nuclear DNA is under constant ROS oxidative DNA damage and unusual DNA repair mechanisms, or even low fidelity polymerase that bypass lesions which could also be affecting the organellar genomes. Finally, both PT and MT phylogenetic trees were well resolved and highly supported, providing a congruent phylogenomic hypothesis for Utricularia and Genlisea clade given the study sampling.
Collapse
Affiliation(s)
- Saura R Silva
- UNESP - São Paulo State University, School of Agricultural and Veterinarian Sciences, Department of Agricultural and Environmental Biotechnology, Campus Jaboticabal, CEP 14884-900 SP, Brazil.
| | - Vitor F O Miranda
- UNESP - São Paulo State University, School of Agricultural and Veterinarian Sciences, Department of Biology, Laboratory of Plant Systematics, Campus Jaboticabal, CEP 14884-900 SP, Brazil.
| | - Todd P Michael
- Plant Molecular and Cellular Biology Laboratory, The Salk Institute for Biological Studies, La Jolla, CA 92037, USA.
| | - Bartosz J Płachno
- Department of Plant Cytology and Embryology, Institute of Botany, Faculty of Biology, Jagiellonian University in Kraków, Gronostajowa 9 St., 30-387 Cracow, Poland.
| | - Ramon G Matos
- UNESP - São Paulo State University, School of Agricultural and Veterinarian Sciences, Department of Biology, Laboratory of Plant Systematics, Campus Jaboticabal, CEP 14884-900 SP, Brazil.
| | - Lubomir Adamec
- Department of Experimental and Functional Morphology, Institute of Botany CAS, Dukelská 135, CZ-379 01 Třeboň, Czech Republic.
| | - Sergei L K Pond
- Institute for Genomics and Evolutionary Medicine, Temple University, Philadelphia, PA 19122, USA.
| | - Alexander G Lucaci
- Institute for Genomics and Evolutionary Medicine, Temple University, Philadelphia, PA 19122, USA.
| | - Daniel G Pinheiro
- UNESP - São Paulo State University, School of Agricultural and Veterinarian Sciences, Department of Agricultural and Environmental Biotechnology, Campus Jaboticabal, CEP 14884-900 SP, Brazil.
| | - Alessandro M Varani
- UNESP - São Paulo State University, School of Agricultural and Veterinarian Sciences, Department of Agricultural and Environmental Biotechnology, Campus Jaboticabal, CEP 14884-900 SP, Brazil.
| |
Collapse
|
2
|
De Novo Hybrid Assembly of the Salvia miltiorrhiza Mitochondrial Genome Provides the First Evidence of the Multi-Chromosomal Mitochondrial DNA Structure of Salvia Species. Int J Mol Sci 2022; 23:ijms232214267. [PMID: 36430747 PMCID: PMC9694629 DOI: 10.3390/ijms232214267] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/29/2022] [Revised: 11/09/2022] [Accepted: 11/15/2022] [Indexed: 11/19/2022] Open
Abstract
Salvia miltiorrhiza has been an economically important medicinal plant. Previously, an S. miltiorrhiza mitochondrial genome (mitogenome) assembled from Illumina short reads, appearing to be a single circular molecule, has been published. Based on the recent reports on the plant mitogenome structure, we suspected that this conformation does not accurately represent the complexity of the S. miltiorrhiza mitogenome. In the current study, we assembled the mitogenome of S. miltiorrhiza using the PacBio and Illumina sequencing technologies. The primary structure of the mitogenome contained two mitochondrial chromosomes (MC1 and MC2), which corresponded to two major conformations, namely, Mac1 and Mac2, respectively. Using two approaches, including (1) long reads mapping and (2) polymerase chain reaction amplification followed by Sanger sequencing, we observed nine repeats that can mediate recombination. We predicted 55 genes, including 33 mitochondrial protein-coding genes (PCGs), 3 rRNA genes, and 19 tRNA genes. Repeat analysis identified 112 microsatellite repeats and 3 long-tandem repeats. Phylogenetic analysis using the 26 shared PCGs resulted in a tree that was congruent with the phylogeny of Lamiales species in the APG IV system. The analysis of mitochondrial plastid DNA (MTPT) identified 16 MTPTs in the mitogenome. Moreover, the analysis of nucleotide substitution rates in Lamiales showed that the genes atp4, ccmB, ccmFc, and mttB might have been positively selected. The results lay the foundation for future studies on the evolution of the Salvia mitogenome and the molecular breeding of S. miltiorrhiza.
Collapse
|
3
|
The complete mitochondrial genome of carnivorous Genlisea tuberosa (Lentibulariaceae): Structure and evolutionary aspects. Gene 2022; 824:146391. [PMID: 35259463 DOI: 10.1016/j.gene.2022.146391] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/28/2021] [Revised: 02/18/2022] [Accepted: 02/28/2022] [Indexed: 11/21/2022]
Abstract
Sequenced genomic data for carnivorous plants are scarce, especially regarding the mitogenomes (MTs) and further studies are crucial to obtain a better understanding of the topic. In this study, we sequenced and characterized the mitochondrial genome of the tuberous carnivorous plant Genlisea tuberosa, being the first of its genus to be sequenced. The genome comprises 729,765 bp, encoding 80 identified genes of which 36 are protein-coding, 40 tRNA, four rRNA genes, and three pseudogenes. An intronic region from the cox1 gene was identified that encodes an endonuclease enzyme that is present in the other sequenced species of Lentibulariaceae. Chloroplast genes (pseudogene and complete) inserted in the MT genome were identified, showing possible horizontal transfer between organelles. In addition, 50 pairs of long repeats from 94 to 274 bp are present, possibly playing an important role in the maintenance of the MT genome. Phylogenetic analysis carried out with 34 coding mitochondrial genes corroborated the positioning of the species listed here within the family. The molecular dynamism in the mitogenome (e.g. the loss or pseudogenization of genes, insertion of foreign genes, the long repeats as well as accumulated mutations) may be reflections of the carnivorous lifestyle where a significant part of cellular energy was shifted for the adaptation of leaves into traps molding the mitochondrial DNA. The sequence and annotation of G. tuberosa's MT will be useful for further studies and serve as a model for evolutionary and taxonomic clarifications of the group as well as improving our comprehension of MT evolution.
Collapse
|
4
|
Transcriptional Landscape and Splicing Efficiency in Arabidopsis Mitochondria. Cells 2021; 10:cells10082054. [PMID: 34440822 PMCID: PMC8392254 DOI: 10.3390/cells10082054] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/30/2021] [Revised: 07/23/2021] [Accepted: 07/30/2021] [Indexed: 12/18/2022] Open
Abstract
Plant mitochondrial transcription is initiated from multiple promoters without an apparent motif, which precludes their identification in other species based on sequence comparisons. Even though coding regions take up only a small fraction of plant mitochondrial genomes, deep RNAseq studies uncovered that these genomes are fully or nearly fully transcribed with significantly different RNA read depth across the genome. Transcriptomic analysis can be a powerful tool to understand the transcription process in diverse angiosperms, including the identification of potential promoters and co-transcribed genes or to study the efficiency of intron splicing. In this work, we analyzed the transcriptional landscape of the Arabidopsis mitochondrial genome (mtDNA) based on large-scale RNA sequencing data to evaluate the use of RNAseq to study those aspects of the transcription process. We found that about 98% of the Arabidopsis mtDNA is transcribed with highly different RNA read depth, which was elevated in known genes. The location of a sharp increase in RNA read depth upstream of genes matched the experimentally identified promoters. The continuously high RNA read depth across two adjacent genes agreed with the known co-transcribed units in Arabidopsis mitochondria. Most intron-containing genes showed a high splicing efficiency with no differences between cis and trans-spliced introns or between genes with distinct splicing mechanisms. Deep RNAseq analyses of diverse plant species will be valuable to recognize general and lineage-specific characteristics related to the mitochondrial transcription process.
Collapse
|
5
|
Garcia LE, Edera AA, Palmer JD, Sato H, Sanchez-Puerta MV. Horizontal gene transfers dominate the functional mitochondrial gene space of a holoparasitic plant. THE NEW PHYTOLOGIST 2021; 229:1701-1714. [PMID: 32929737 DOI: 10.1111/nph.16926] [Citation(s) in RCA: 15] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 06/26/2020] [Accepted: 09/02/2020] [Indexed: 06/11/2023]
Abstract
Although horizontal gene transfer (HGT) is common in angiosperm mitochondrial DNAs (mtDNAs), few cases of functional foreign genes have been identified. The one outstanding candidate for large-scale functional HGT is the holoparasite Lophophytum mirabile, whose mtDNA has lost most native genes but contains intact foreign homologs acquired from legume host plants. To investigate the extent to which this situation results from functional replacement of native by foreign genes, functional mitochondrial gene transfer to the nucleus, and/or loss of mitochondrial biochemical function in the context of extreme parasitism, we examined the Lophophytum mitochondrial and nuclear transcriptomes by deep paired-end RNA sequencing. Most foreign mitochondrial genes in Lophophytum are highly transcribed, accurately spliced, and efficiently RNA edited. By contrast, we found no evidence for functional gene transfer to the nucleus or loss of mitochondrial functions in Lophophytum. Many functional replacements occurred via the physical replacement of native genes by foreign genes. Some of these events probably occurred as the final act of HGT itself. Lophophytum mtDNA has experienced an unprecedented level of functional replacement of native genes by foreign copies. This raises important questions concerning population-genetic and molecular regimes that underlie such a high level of foreign gene takeover.
Collapse
Affiliation(s)
- Laura E Garcia
- IBAM, Universidad Nacional de Cuyo, CONICET, Facultad de Ciencias Agrarias, Almirante Brown 500, Chacras de Coria, M5528AHB, Argentina
- Facultad de Ciencias Exactas y Naturales, Universidad Nacional de Cuyo, Padre Jorge Contreras 1300, Mendoza, M5502JMA, Argentina
| | - Alejandro A Edera
- IBAM, Universidad Nacional de Cuyo, CONICET, Facultad de Ciencias Agrarias, Almirante Brown 500, Chacras de Coria, M5528AHB, Argentina
| | - Jeffrey D Palmer
- Department of Biology, Indiana University, Bloomington, IN, 47405, USA
| | - Hector Sato
- Facultad de Ciencias Agrarias (UNJu), Cátedra de Botánica General-Herbario JUA, Alberdi 47, Jujuy, CP 4600, Argentina
| | - M Virginia Sanchez-Puerta
- IBAM, Universidad Nacional de Cuyo, CONICET, Facultad de Ciencias Agrarias, Almirante Brown 500, Chacras de Coria, M5528AHB, Argentina
- Facultad de Ciencias Exactas y Naturales, Universidad Nacional de Cuyo, Padre Jorge Contreras 1300, Mendoza, M5502JMA, Argentina
| |
Collapse
|
6
|
Feng Y, Xiang X, Akhter D, Pan R, Fu Z, Jin X. Mitochondrial Phylogenomics of Fagales Provides Insights Into Plant Mitogenome Mosaic Evolution. FRONTIERS IN PLANT SCIENCE 2021; 12:762195. [PMID: 34733309 PMCID: PMC8558628 DOI: 10.3389/fpls.2021.762195] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 08/21/2021] [Accepted: 09/17/2021] [Indexed: 05/13/2023]
Abstract
Fagales are an order of woody plants and comprise more than 1,100 species, most of which produce economically important timbers, nuts, and fruits. Their nuclear and plastid genomes are well-sequenced and provided valuable resources to study their phylogeny, breeding, resistance, etc. However, little is known about the mitochondrial genomes (mitogenomes), which hinder a full understanding of their genome evolution. In this study, we assembled complete mitogenomes of 23 species, covering five of the seven families of Fagales. These mitogenomes had similar gene sets but varied 2.4 times in size. The mitochondrial genes were highly conserved, and their capacity in phylogeny was challenging. The mitogenomic structure was extremely dynamic, and synteny among species was poor. Further analyses of the Fagales mitogenomes revealed extremely mosaic characteristics, with horizontal transfer (HGT)-like sequences from almost all seed plant taxa and even mitoviruses. The largest mitogenome, Carpinus cordata, did not have large amounts of specific sequences but instead contained a high proportion of sequences homologous to other Fagales. Independent and unequal transfers of third-party DNA, including nuclear genome and other resources, may partially account for the HGT-like fragments and unbalanced size expansions observed in Fagales mitogenomes. Supporting this, a mitochondrial plasmid-like of nuclear origin was found in Carpinus. Overall, we deciphered the last genetic materials of Fagales, and our large-scale analyses provide new insights into plant mitogenome evolution and size variation.
Collapse
Affiliation(s)
- Yanlei Feng
- Institute of Biology, Westlake Institute for Advanced Study, Hangzhou, China
- College of Agriculture and Biotechnology, Institute of Crop Science, Zhejiang University, Hangzhou, China
- School of Life Sciences, Westlake University, Hangzhou, China
- Yanlei Feng
| | - Xiaoguo Xiang
- School of Life Sciences, Nanchang University, Nanchang, China
| | - Delara Akhter
- College of Agriculture and Biotechnology, Institute of Crop Science, Zhejiang University, Hangzhou, China
- Department of Genetics and Plant Breeding, Sylhet Agricultural University, Sylhet, Bangladesh
| | - Ronghui Pan
- College of Agriculture and Biotechnology, Institute of Crop Science, Zhejiang University, Hangzhou, China
| | - Zhixi Fu
- College of Life Science, Sichuan Normal University, Chengdu, China
- Sustainable Development Research Center of Resources and Environment of Western Sichuan, Sichuan Normal University, Chengdu, China
- *Correspondence: Zhixi Fu
| | - Xiaohua Jin
- Institute of Botany, The Chinese Academy of Sciences, Beijing, China
- Xiaohua Jin
| |
Collapse
|
7
|
Zhou J, Cao Y, Liu X, Jiang H, Li W. Bladder entrance of microplastic likely induces toxic effects in carnivorous macrophyte Utricularia aurea Lour. ENVIRONMENTAL SCIENCE AND POLLUTION RESEARCH INTERNATIONAL 2020; 27:32124-32131. [PMID: 32564318 DOI: 10.1007/s11356-020-09529-y] [Citation(s) in RCA: 11] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 03/02/2020] [Accepted: 05/29/2020] [Indexed: 06/11/2023]
Abstract
The global distribution of microplastic (particle size < 5 mm) is of growing concern, especially in aquatic environments where it may cause adverse effects on resident organisms. To date, however, few studies have focused on the impacts of microplastic on aquatic plants. Here, we conducted a microcosm study to investigate the toxic effects of microplastic on the carnivorous aquatic macrophyte Utricularia aurea Lour. Based on microscopic images and Raman spectrum analysis, we found that most polyvinyl chloride (PVC) particles were smaller than the valve of U. aurea bladders, thus allowing entrance into the plant, but this was not so for polyethylene (PE) particles. Furthermore, PVC (50 mg L-1) had significantly negative effects on growth and physiological parameters such as macrophyte length, chlorophyll content, and fluorescence, whereas, at the same concentration, PE had no such effects. Further analysis revealed that after bladder removal, the macrophytes did not respond to PVC particle toxicity. Thus, intake of microplastics (i.e., PVC) through bladders is likely responsible for inducing toxic effects to the growth and physiological parameters of U. aurea.
Collapse
Affiliation(s)
- Jingzhe Zhou
- Key Laboratory of Aquatic Botany and Watershed Ecology, Wuhan Botanical Garden, Chinese Academy of Sciences, Wuhan, 430074, China
- University of Chinese Academy of Sciences, Beijing, 100049, China
| | - Yu Cao
- Key Laboratory of Aquatic Botany and Watershed Ecology, Wuhan Botanical Garden, Chinese Academy of Sciences, Wuhan, 430074, China.
- University of Chinese Academy of Sciences, Beijing, 100049, China.
| | - Xiaoning Liu
- Institute of Hydroecology, State Key Laboratory of Water Resources and Hydropower Engineering Science, Wuhan University, Wuhan, 430072, China
| | - Hongsheng Jiang
- Key Laboratory of Aquatic Botany and Watershed Ecology, Wuhan Botanical Garden, Chinese Academy of Sciences, Wuhan, 430074, China
- University of Chinese Academy of Sciences, Beijing, 100049, China
| | - Wei Li
- Key Laboratory of Aquatic Botany and Watershed Ecology, Wuhan Botanical Garden, Chinese Academy of Sciences, Wuhan, 430074, China
- University of Chinese Academy of Sciences, Beijing, 100049, China
| |
Collapse
|
8
|
Abstract
Ever since its discovery, the double-stranded DNA contained in the mitochondria of eukaryotes has fascinated researchers because of its bacterial endosymbiotic origin, crucial role in encoding subunits of the respiratory complexes, compact nature, and specific inheritance mechanisms. In the last few years, high-throughput sequencing techniques have accelerated the sequencing of mitochondrial genomes (mitogenomes) and uncovered the great diversity of organizations, gene contents, and modes of replication and transcription found in living eukaryotes. Some early divergent lineages of unicellular eukaryotes retain certain synteny and gene content resembling those observed in the genomes of alphaproteobacteria (the inferred closest living group of mitochondria), whereas others adapted to anaerobic environments have drastically reduced or even lost the mitogenome. In the three main multicellular lineages of eukaryotes, mitogenomes have pursued diverse evolutionary trajectories in which different types of molecules (circular versus linear and single versus multipartite), gene structures (with or without self-splicing introns), gene contents, gene orders, genetic codes, and transfer RNA editing mechanisms have been selected. Whereas animals have evolved a rather compact mitochondrial genome between 11 and 50 Kb in length with a highly conserved gene content in bilaterians, plants exhibit large mitochondrial genomes of 66 Kb to 11.3 Mb with large intergenic repetitions prone to recombination, and fungal mitogenomes have intermediate sizes of 12 to 236 Kb.
Collapse
Affiliation(s)
- Rafael Zardoya
- Departamento de Biodiversidad y Biología Evolutiva, Museo Nacional de Ciencias Naturales (MNCN-CSIC), Madrid, Spain
| |
Collapse
|
9
|
Methods and Tools for Plant Organelle Genome Sequencing, Assembly, and Downstream Analysis. Methods Mol Biol 2020; 2107:49-98. [PMID: 31893443 DOI: 10.1007/978-1-0716-0235-5_4] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/03/2023]
Abstract
Organelles play an important role in a eukaryotic cell. Among them, the two organelles, chloroplast and mitochondria, are responsible for the critical function of photosynthesis and aerobic respiration. Organellar genomes are also very important for plant systematic studies. Here we have described the methods for isolation of the mitochondrial and plastid DNA and its subsequent sequencing with the help of NGS technology. We have also discussed in detail the various tools available for assembly, annotation, and visualization of the organelle genome sequence.
Collapse
|
10
|
Silva SR, Moraes AP, Penha HA, Julião MHM, Domingues DS, Michael TP, Miranda VFO, Varani AM. The Terrestrial Carnivorous Plant Utricularia reniformis Sheds Light on Environmental and Life-Form Genome Plasticity. Int J Mol Sci 2019; 21:E3. [PMID: 31861318 PMCID: PMC6982007 DOI: 10.3390/ijms21010003] [Citation(s) in RCA: 10] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/23/2019] [Revised: 12/13/2019] [Accepted: 12/15/2019] [Indexed: 12/22/2022] Open
Abstract
Utricularia belongs to Lentibulariaceae, a widespread family of carnivorous plants that possess ultra-small and highly dynamic nuclear genomes. It has been shown that the Lentibulariaceae genomes have been shaped by transposable elements expansion and loss, and multiple rounds of whole-genome duplications (WGD), making the family a platform for evolutionary and comparative genomics studies. To explore the evolution of Utricularia, we estimated the chromosome number and genome size, as well as sequenced the terrestrial bladderwort Utricularia reniformis (2n = 40, 1C = 317.1-Mpb). Here, we report a high quality 304 Mb draft genome, with a scaffold NG50 of 466-Kb, a BUSCO completeness of 87.8%, and 42,582 predicted genes. Compared to the smaller and aquatic U. gibba genome (101 Mb) that has a 32% repetitive sequence, the U. reniformis genome is highly repetitive (56%). The structural differences between the two genomes are the result of distinct fractionation and rearrangements after WGD, and massive proliferation of LTR-retrotransposons. Moreover, GO enrichment analyses suggest an ongoing gene birth-death-innovation process occurring among the tandem duplicated genes, shaping the evolution of carnivory-associated functions. We also identified unique patterns of developmentally related genes that support the terrestrial life-form and body plan of U. reniformis. Collectively, our results provided additional insights into the evolution of the plastic and specialized Lentibulariaceae genomes.
Collapse
Affiliation(s)
- Saura R. Silva
- Departamento de Tecnologia, Faculdade de Ciências Agrárias e Veterinárias, UNESP—Universidade Estadual Paulista, Jaboticabal 14884-900, Brazil; (S.R.S.); (H.A.P.); (M.H.M.J.)
| | - Ana Paula Moraes
- Centro de Ciências Naturais e Humanas, Universidade Federal do ABC, São Bernardo do Campo 09606-070, Brazil;
| | - Helen A. Penha
- Departamento de Tecnologia, Faculdade de Ciências Agrárias e Veterinárias, UNESP—Universidade Estadual Paulista, Jaboticabal 14884-900, Brazil; (S.R.S.); (H.A.P.); (M.H.M.J.)
| | - Maria H. M. Julião
- Departamento de Tecnologia, Faculdade de Ciências Agrárias e Veterinárias, UNESP—Universidade Estadual Paulista, Jaboticabal 14884-900, Brazil; (S.R.S.); (H.A.P.); (M.H.M.J.)
| | - Douglas S. Domingues
- Departamento de Botânica, Instituto de Biociências, UNESP—Universidade Estadual Paulista, Rio Claro 13506-900, Brazil;
| | | | - Vitor F. O. Miranda
- Departamento de Biologia Aplicada à Agropecuária, Faculdade de Ciências Agrárias e Veterinárias, UNESP—Universidade Estadual Paulista, Jaboticabal 14884-900, Brazil
| | - Alessandro M. Varani
- Departamento de Tecnologia, Faculdade de Ciências Agrárias e Veterinárias, UNESP—Universidade Estadual Paulista, Jaboticabal 14884-900, Brazil; (S.R.S.); (H.A.P.); (M.H.M.J.)
| |
Collapse
|
11
|
Starr EP, Nuccio EE, Pett-Ridge J, Banfield JF, Firestone MK. Metatranscriptomic reconstruction reveals RNA viruses with the potential to shape carbon cycling in soil. Proc Natl Acad Sci U S A 2019; 116:25900-25908. [PMID: 31772013 PMCID: PMC6926006 DOI: 10.1073/pnas.1908291116] [Citation(s) in RCA: 105] [Impact Index Per Article: 21.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/18/2022] Open
Abstract
Viruses impact nearly all organisms on Earth, with ripples of influence in agriculture, health, and biogeochemical processes. However, very little is known about RNA viruses in an environmental context, and even less is known about their diversity and ecology in soil, 1 of the most complex microbial systems. Here, we assembled 48 individual metatranscriptomes from 4 habitats within a planted soil sampled over a 22-d time series: Rhizosphere alone, detritosphere alone, rhizosphere with added root detritus, and unamended soil (4 time points and 3 biological replicates). We resolved the RNA viral community, uncovering a high diversity of viral sequences. We also investigated possible host organisms by analyzing metatranscriptome marker genes. Based on viral phylogeny, much of the diversity was Narnaviridae that may parasitize fungi or Leviviridae, which may infect Proteobacteria. Both host and viral communities appear to be highly dynamic, and rapidly diverged depending on experimental conditions. The viral and host communities were structured based on the presence of root litter. Clear temporal dynamics by Leviviridae and their hosts indicated that viruses were replicating. With this time-resolved analysis, we show that RNA viruses are diverse, abundant, and active in soil. When viral infection causes host cell death, it may mobilize cell carbon in a process that may represent an overlooked component of soil carbon cycling.
Collapse
Affiliation(s)
- Evan P Starr
- Department of Plant and Microbial Biology, University of California, Berkeley, CA 94720
| | - Erin E Nuccio
- Physical and Life Sciences Directorate, Lawrence Livermore National Laboratory, Livermore, CA 94550
| | - Jennifer Pett-Ridge
- Physical and Life Sciences Directorate, Lawrence Livermore National Laboratory, Livermore, CA 94550
| | - Jillian F Banfield
- Department of Earth and Planetary Science, University of California, Berkeley, CA 94720;
- Earth Sciences Division, Lawrence Berkeley National Laboratory, Berkeley, CA 94720
- Department of Environmental Science, Policy, and Management, University of California, Berkeley, CA 94720
- Chan Zuckerberg Biohub, San Francisco, CA 94158
- Innovative Genomics Institute, Berkeley, CA 94720
| | - Mary K Firestone
- Earth Sciences Division, Lawrence Berkeley National Laboratory, Berkeley, CA 94720;
- Department of Environmental Science, Policy, and Management, University of California, Berkeley, CA 94720
| |
Collapse
|
12
|
Silva SR, Pinheiro DG, Penha HA, Płachno BJ, Michael TP, Meer EJ, Miranda VFO, Varani AM. Intraspecific Variation within the Utricularia amethystina Species Morphotypes Based on Chloroplast Genomes. Int J Mol Sci 2019; 20:E6130. [PMID: 31817365 PMCID: PMC6940893 DOI: 10.3390/ijms20246130] [Citation(s) in RCA: 17] [Impact Index Per Article: 3.4] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/18/2019] [Revised: 11/09/2019] [Accepted: 11/27/2019] [Indexed: 01/09/2023] Open
Abstract
Utricularia amethystina Salzm. ex A.St.-Hil. & Girard (Lentibulariaceae) is a highly polymorphic carnivorous plant taxonomically rearranged many times throughout history. Herein, the complete chloroplast genomes (cpDNA) of three U. amethystina morphotypes: purple-, white-, and yellow-flowered, were sequenced, compared, and putative markers for systematic, populations, and evolutionary studies were uncovered. In addition, RNA-Seq and RNA-editing analysis were employed for functional cpDNA evaluation. The cpDNA of three U. amethystina morphotypes exhibits typical quadripartite structure. Fine-grained sequence comparison revealed a high degree of intraspecific genetic variability in all morphotypes, including an exclusive inversion in the psbM and petN genes in U. amethystina yellow. Phylogenetic analyses indicate that U. amethystina morphotypes are monophyletic. Furthermore, in contrast to the terrestrial Utricularia reniformis cpDNA, the U. amethystina morphotypes retain all the plastid NAD(P)H-dehydrogenase (ndh) complex genes. This observation supports the hypothesis that the ndhs in terrestrial Utricularia were independently lost and regained, also suggesting that different habitats (aquatic and terrestrial) are not related to the absence of Utricularia ndhs gene repertoire as previously assumed. Moreover, RNA-Seq analyses recovered similar patterns, including nonsynonymous RNA-editing sites (e.g., rps14 and petB). Collectively, our results bring new insights into the chloroplast genome architecture and evolution of the photosynthesis machinery in the Lentibulariaceae.
Collapse
Affiliation(s)
- Saura R. Silva
- Departamento de Tecnologia, Faculdade de Ciências Agrárias e Veterinárias, Jaboticabal, Universidade Estadual Paulista (Unesp), Sao Paulo 14884-900, Brazil; (D.G.P.); (H.A.P.)
| | - Daniel G. Pinheiro
- Departamento de Tecnologia, Faculdade de Ciências Agrárias e Veterinárias, Jaboticabal, Universidade Estadual Paulista (Unesp), Sao Paulo 14884-900, Brazil; (D.G.P.); (H.A.P.)
| | - Helen A. Penha
- Departamento de Tecnologia, Faculdade de Ciências Agrárias e Veterinárias, Jaboticabal, Universidade Estadual Paulista (Unesp), Sao Paulo 14884-900, Brazil; (D.G.P.); (H.A.P.)
| | - Bartosz J. Płachno
- Department of Plant Cytology and Embryology, Institute of Botany, Faculty of Biology, Jagiellonian University in Kraków, 30-387 Krakow, Poland;
| | | | | | - Vitor F. O. Miranda
- Departamento de Biologia Aplicada à Agropecuária, Faculdade de Ciências Agrárias e Veterinárias, Jaboticabal, Universidade Estadual Paulista (Unesp), Sao Paulo 14884-900, Brazil
| | - Alessandro M. Varani
- Departamento de Tecnologia, Faculdade de Ciências Agrárias e Veterinárias, Jaboticabal, Universidade Estadual Paulista (Unesp), Sao Paulo 14884-900, Brazil; (D.G.P.); (H.A.P.)
| |
Collapse
|
13
|
Kozik A, Rowan BA, Lavelle D, Berke L, Schranz ME, Michelmore RW, Christensen AC. The alternative reality of plant mitochondrial DNA: One ring does not rule them all. PLoS Genet 2019; 15:e1008373. [PMID: 31469821 PMCID: PMC6742443 DOI: 10.1371/journal.pgen.1008373] [Citation(s) in RCA: 144] [Impact Index Per Article: 28.8] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/24/2019] [Revised: 09/12/2019] [Accepted: 08/16/2019] [Indexed: 01/27/2023] Open
Abstract
Plant mitochondrial genomes are usually assembled and displayed as circular maps based on the widely-held view across the broad community of life scientists that circular genome-sized molecules are the primary form of plant mitochondrial DNA, despite the understanding by plant mitochondrial researchers that this is an inaccurate and outdated concept. Many plant mitochondrial genomes have one or more pairs of large repeats that can act as sites for inter- or intramolecular recombination, leading to multiple alternative arrangements (isoforms). Most mitochondrial genomes have been assembled using methods unable to capture the complete spectrum of isoforms within a species, leading to an incomplete inference of their structure and recombinational activity. To document and investigate underlying reasons for structural diversity in plant mitochondrial DNA, we used long-read (PacBio) and short-read (Illumina) sequencing data to assemble and compare mitochondrial genomes of domesticated (Lactuca sativa) and wild (L. saligna and L. serriola) lettuce species. We characterized a comprehensive, complex set of isoforms within each species and compared genome structures between species. Physical analysis of L. sativa mtDNA molecules by fluorescence microscopy revealed a variety of linear, branched, and circular structures. The mitochondrial genomes for L. sativa and L. serriola were identical in sequence and arrangement and differed substantially from L. saligna, indicating that the mitochondrial genome structure did not change during domestication. From the isoforms in our data, we infer that recombination occurs at repeats of all sizes at variable frequencies. The differences in genome structure between L. saligna and the two other Lactuca species can be largely explained by rare recombination events that rearranged the structure. Our data demonstrate that representations of plant mitochondrial genomes as simple, circular molecules are not accurate descriptions of their true nature and that in reality plant mitochondrial DNA is a complex, dynamic mixture of forms.
Collapse
Affiliation(s)
- Alexander Kozik
- Genome Center and Department of Plant Sciences, University of California, Davis, California, United States of America
| | - Beth A. Rowan
- Genome Center and Department of Plant Sciences, University of California, Davis, California, United States of America
| | - Dean Lavelle
- Genome Center and Department of Plant Sciences, University of California, Davis, California, United States of America
| | - Lidija Berke
- Wageningen University & Research, PB Wageningen, Gelderland, The Netherlands
| | - M. Eric Schranz
- Wageningen University & Research, PB Wageningen, Gelderland, The Netherlands
| | - Richard W. Michelmore
- Genome Center and Department of Plant Sciences, University of California, Davis, California, United States of America
| | - Alan C. Christensen
- School of Biological Sciences, University of Nebraska - Lincoln, Lincoln, Nebraska, United States of America
| |
Collapse
|
14
|
Roossinck MJ. Evolutionary and ecological links between plant and fungal viruses. THE NEW PHYTOLOGIST 2019; 221:86-92. [PMID: 30084143 DOI: 10.1111/nph.15364] [Citation(s) in RCA: 43] [Impact Index Per Article: 8.6] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 01/30/2018] [Accepted: 06/21/2018] [Indexed: 06/08/2023]
Abstract
Contents Summary 86 I. Introduction 86 II. Lineages shared by plant and fungal viruses 87 III. Virus transmission between plants and fungi 90 IV. Additional plant virus families identified in fungi by metagenomics 91 Acknowledgements 91 References 91 SUMMARY: Plants and microorganisms have been interacting in both positive and negative ways for millions of years. They are also frequently infected with viruses that can have positive or negative impacts. A majority of virus families with members that infect fungi have counterparts that infect plants, and in some cases the phylogenetic analyses of these virus families indicate transmission between the plant and fungal kingdoms. These similarities reflect the host relationships; fungi are evolutionarily more closely related to animals than to plants but share very few viral signatures with animal viruses. The details of several of these interactions are described, and the evolutionary implications of viral cross-kingdom interactions and horizontal gene transfer are proposed.
Collapse
Affiliation(s)
- Marilyn J Roossinck
- Department of Plant Pathology and Environmental Microbiology, Center for Infectious Disease Dynamics, Pennsylvania State University, University Park, PA, 16802, USA
| |
Collapse
|
15
|
Abstract
RNA editing is a fundamental biochemical process relating to the modification of nucleotides in messenger RNAs of functional genes in cells. RNA editing leads to re-establishment of conserved amino acid residues for functional proteins in nuclei, chloroplasts, and mitochondria. Identification of RNA editing factors that contributes to target site recognition increases our understanding of RNA editing mechanisms. Significant progress has been made in recent years in RNA editing studies for both animal and plant cells. RNA editing in nuclei and mitochondria of animal cells and in chloroplast of plant cells has been extensively documented and reviewed. RNA editing has been also extensively documented on plant mitochondria. However, functional diversity of RNA editing factors in plant mitochondria is not overviewed. Here, we review the biological significance of RNA editing, recent progress on the molecular mechanisms of RNA editing process, and function diversity of editing factors in plant mitochondrial research. We will focus on: (1) pentatricopeptide repeat proteins in Arabidopsis and in crop plants; (2) the progress of RNA editing process in plant mitochondria; (3) RNA editing-related RNA splicing; (4) RNA editing associated flower development; (5) RNA editing modulated male sterile; (6) RNA editing-regulated cell signaling; and (7) RNA editing involving abiotic stress. Advances described in this review will be valuable in expanding our understanding in RNA editing. The diverse functions of RNA editing in plant mitochondria will shed light on the investigation of molecular mechanisms that underlies plant development and abiotic stress tolerance.
Collapse
|
16
|
Sirová D, Bárta J, Šimek K, Posch T, Pech J, Stone J, Borovec J, Adamec L, Vrba J. Hunters or farmers? Microbiome characteristics help elucidate the diet composition in an aquatic carnivorous plant. MICROBIOME 2018; 6:225. [PMID: 30558682 PMCID: PMC6297986 DOI: 10.1186/s40168-018-0600-7] [Citation(s) in RCA: 18] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 07/24/2018] [Accepted: 11/18/2018] [Indexed: 05/09/2023]
Abstract
BACKGROUND Utricularia are rootless aquatic carnivorous plants which have recently attracted the attention of researchers due to the peculiarities of their miniaturized genomes. Here, we focus on a novel aspect of Utricularia ecophysiology-the interactions with and within the complex communities of microorganisms colonizing their traps and external surfaces. RESULTS Bacteria, fungi, algae, and protozoa inhabit the miniature ecosystem of the Utricularia trap lumen and are involved in the regeneration of nutrients from complex organic matter. By combining molecular methods, microscopy, and other approaches to assess the trap-associated microbial community structure, diversity, function, as well as the nutrient turn-over potential of bacterivory, we gained insight into the nutrient acquisition strategies of the Utricularia hosts. CONCLUSIONS We conclude that Utricularia traps can, in terms of their ecophysiological function, be compared to microbial cultivators or farms, which center around complex microbial consortia acting synergistically to convert complex organic matter, often of algal origin, into a source of utilizable nutrients for the plants.
Collapse
Affiliation(s)
- Dagmara Sirová
- Biology Centre CAS, Institute of Hydrobiology, Na Sádkách 7, CZ-37005, České Budějovice, Czech Republic.
- Faculty of Science, University of South Bohemia, Branišovská 1760, CZ-37005, České Budějovice, Czech Republic.
| | - Jiří Bárta
- Faculty of Science, University of South Bohemia, Branišovská 1760, CZ-37005, České Budějovice, Czech Republic
| | - Karel Šimek
- Biology Centre CAS, Institute of Hydrobiology, Na Sádkách 7, CZ-37005, České Budějovice, Czech Republic
- Faculty of Science, University of South Bohemia, Branišovská 1760, CZ-37005, České Budějovice, Czech Republic
| | - Thomas Posch
- Limnological Station, Department of Plant and Microbial Biology, University of Zurich, CH-8802, Kilchberg, Switzerland
| | - Jiří Pech
- Faculty of Science, University of South Bohemia, Branišovská 1760, CZ-37005, České Budějovice, Czech Republic
| | - James Stone
- Department of Biology and Wildlife, University of Alaska Fairbanks, Fairbanks, AK-99775, USA
- Institute of Experimental Botany CAS, Rozvojová 263, CZ-16502, Praha 6-Lysolaje, Czech Republic
| | - Jakub Borovec
- Biology Centre CAS, Institute of Hydrobiology, Na Sádkách 7, CZ-37005, České Budějovice, Czech Republic
| | - Lubomír Adamec
- Institute of Botany CAS, Dukelská 135, CZ-37982, Třeboň, Czech Republic
| | - Jaroslav Vrba
- Biology Centre CAS, Institute of Hydrobiology, Na Sádkách 7, CZ-37005, České Budějovice, Czech Republic
- Faculty of Science, University of South Bohemia, Branišovská 1760, CZ-37005, České Budějovice, Czech Republic
| |
Collapse
|
17
|
Gruzdev EV, Mardanov AV, Beletsky AV, Ravin NV, Skryabin KG. The complete mitochondrial genome of the carnivorous flowering plant Nepenthes X Ventrata. MITOCHONDRIAL DNA PART B-RESOURCES 2018; 3:1259-1260. [PMID: 33474484 PMCID: PMC7800874 DOI: 10.1080/23802359.2018.1532353] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Indexed: 11/07/2022]
Abstract
Carnivorous plants have the ability to capture and digest small animals as a source of additional nutrients, which allows them to grow in nutrient-poor habitats. This study reports the complete mitochondrial genome sequence of a pitcher plant Nepenthes x ventrata. It was 520,764 bp in size with a GC content of 44.17% and contained 37 protein-coding genes, 2 pseudogenes, 18 tRNA genes and 3 rRNA genes. Four tRNA genes and the rps11 gene were probably transferred to mitochondrion form the chloroplast genome. Phylogenetic analysis confirmed that N. x ventrata belongs to the order Caryophyllales.
Collapse
Affiliation(s)
- Eugene V Gruzdev
- Institute of Bioengineering, Research Center of Biotechnology of the Russian Academy of Sciences, Moscow, Russia
| | - Andrey V Mardanov
- Institute of Bioengineering, Research Center of Biotechnology of the Russian Academy of Sciences, Moscow, Russia
| | - Alexey V Beletsky
- Institute of Bioengineering, Research Center of Biotechnology of the Russian Academy of Sciences, Moscow, Russia
| | - Nikolai V Ravin
- Institute of Bioengineering, Research Center of Biotechnology of the Russian Academy of Sciences, Moscow, Russia
| | - Konstantin G Skryabin
- Institute of Bioengineering, Research Center of Biotechnology of the Russian Academy of Sciences, Moscow, Russia
| |
Collapse
|
18
|
Lee HO, Choi JW, Baek JH, Oh JH, Lee SC, Kim CK. Assembly of the Mitochondrial Genome in the Campanulaceae Family Using Illumina Low-Coverage Sequencing. Genes (Basel) 2018; 9:E383. [PMID: 30061537 PMCID: PMC6116063 DOI: 10.3390/genes9080383] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/28/2018] [Revised: 07/24/2018] [Accepted: 07/25/2018] [Indexed: 11/16/2022] Open
Abstract
Platycodongrandiflorus (balloon flower) and Codonopsislanceolata (bonnet bellflower) are important herbs used in Asian traditional medicine, and both belong to the botanical family Campanulaceae. In this study, we designed and implemented a de novo DNA sequencing and assembly strategy to map the complete mitochondrial genomes of the first two members of the Campanulaceae using low-coverage Illumina DNA sequencing data. We produced a total of 28.9 Gb of paired-end sequencing data from the genomic DNA of P.grandiflorus (20.9 Gb) and C.lanceolata (8.0 Gb). The assembled mitochondrial genome of P.grandiflorus was found to consist of two circular chromosomes; the master circle contains 56 genes, and the minor circle contains 42 genes. The C.lanceolata mitochondrial genome consists of a single circle harboring 54 genes. Using a comparative genome structure and a pattern of repeated sequences, we show that the P.grandiflorus minor circle resulted from a recombination event involving the direct repeats of the master circle. Our dataset will be useful for comparative genomics and for evolutionary studies, and will facilitate further biological and phylogenetic characterization of species in the Campanulaceae.
Collapse
Affiliation(s)
- Hyun-Oh Lee
- Phyzen Genomics Institute, Seongnam 13558, Korea.
- Department of Plant Science, Seoul National University, Seoul 08826, Korea.
| | - Ji-Weon Choi
- Postharvest Technology Division, National Institute of Horticultural and Herbal Science, Wanju 55365, Korea.
| | - Jeong-Ho Baek
- Gene Engineering Division, National Institute of Agricultural Sciences, RDA, Jeonju 54874, Korea.
| | - Jae-Hyeon Oh
- Genomics Division, National Institute of Agricultural Sciences, RDA, Jeonju 54874, Korea.
| | | | - Chang-Kug Kim
- Genomics Division, National Institute of Agricultural Sciences, RDA, Jeonju 54874, Korea.
| |
Collapse
|
19
|
Silva SR, Michael TP, Meer EJ, Pinheiro DG, Varani AM, Miranda VFO. Comparative genomic analysis of Genlisea (corkscrew plants-Lentibulariaceae) chloroplast genomes reveals an increasing loss of the ndh genes. PLoS One 2018; 13:e0190321. [PMID: 29293597 PMCID: PMC5749785 DOI: 10.1371/journal.pone.0190321] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/18/2017] [Accepted: 12/12/2017] [Indexed: 11/18/2022] Open
Abstract
In the carnivorous plant family Lentibulariaceae, all three genome compartments (nuclear, chloroplast, and mitochondria) have some of the highest rates of nucleotide substitutions across angiosperms. While the genera Genlisea and Utricularia have the smallest known flowering plant nuclear genomes, the chloroplast genomes (cpDNA) are mostly structurally conserved except for deletion and/or pseudogenization of the NAD(P)H-dehydrogenase complex (ndh) genes known to be involved in stress conditions of low light or CO2 concentrations. In order to determine how the cpDNA are changing, and to better understand the evolutionary history within the Genlisea genus, we sequenced, assembled and analyzed complete cpDNA from six species (G. aurea, G. filiformis, G. pygmaea, G. repens, G. tuberosa and G. violacea) together with the publicly available G. margaretae cpDNA. In general, the cpDNA structure among the analyzed Genlisea species is highly similar. However, we found that the plastidial ndh genes underwent a progressive process of degradation similar to the other terrestrial Lentibulariaceae cpDNA analyzed to date, but in contrast to the aquatic species. Contrary to current thinking that the terrestrial environment is a more stressful environment and thus requiring the ndh genes, we provide evidence that in the Lentibulariaceae the terrestrial forms have progressive loss while the aquatic forms have the eleven plastidial ndh genes intact. Therefore, the Lentibulariaceae system provides an important opportunity to understand the evolutionary forces that govern the transition to an aquatic environment and may provide insight into how plants manage water stress at a genome scale.
Collapse
Affiliation(s)
- Saura R. Silva
- Universidade Estadual Paulista (Unesp), Botucatu, Instituto de Biociências, São Paulo, Brazil
| | - Todd P. Michael
- J. Craig Venter Institute, La Jolla, CA, United States of America
| | - Elliott J. Meer
- 10X Genomics, Pleasanton, California, United States of America
| | - Daniel G. Pinheiro
- Universidade Estadual Paulista (Unesp), Faculdade de Ciências Agrárias e Veterinárias, Jaboticabal, Departamento de Tecnologia, São Paulo, Brazil
| | - Alessandro M. Varani
- Universidade Estadual Paulista (Unesp), Faculdade de Ciências Agrárias e Veterinárias, Jaboticabal, Departamento de Tecnologia, São Paulo, Brazil
| | - Vitor F. O. Miranda
- Universidade Estadual Paulista (Unesp), Faculdade de Ciências Agrárias e Veterinárias, Jaboticabal, Departamento de Biologia Aplicada à Agropecuária, São Paulo, Brazil
| |
Collapse
|