1
|
Wu ZY, Chapman MA, Liu J, Milne RI, Zhao Y, Luo YH, Zhu GF, Cadotte MW, Luan MB, Fan PZ, Monro AK, Li ZP, Corlett RT, Li DZ. Genomic variation, environmental adaptation, and feralization in ramie, an ancient fiber crop. PLANT COMMUNICATIONS 2024:100942. [PMID: 38720463 DOI: 10.1016/j.xplc.2024.100942] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 11/19/2023] [Revised: 12/20/2023] [Accepted: 05/06/2024] [Indexed: 06/29/2024]
Abstract
Feralization is an important evolutionary process, but the mechanisms behind it remain poorly understood. Here, we use the ancient fiber crop ramie (Boehmeria nivea (L.) Gaudich.) as a model to investigate genomic changes associated with both domestication and feralization. We first produced a chromosome-scale de novo genome assembly of feral ramie and investigated structural variations between feral and domesticated ramie genomes. Next, we gathered 915 accessions from 23 countries, comprising cultivars, major landraces, feral populations, and the wild progenitor. Based on whole-genome resequencing of these accessions, we constructed the most comprehensive ramie genomic variation map to date. Phylogenetic, demographic, and admixture signal detection analyses indicated that feral ramie is of exoferal or exo-endo origin, i.e., descended from hybridization between domesticated ramie and the wild progenitor or ancient landraces. Feral ramie has higher genetic diversity than wild or domesticated ramie, and genomic regions affected by natural selection during feralization differ from those under selection during domestication. Ecological analyses showed that feral and domesticated ramie have similar ecological niches that differ substantially from the niche of the wild progenitor, and three environmental variables are associated with habitat-specific adaptation in feral ramie. These findings advance our understanding of feralization, providing a scientific basis for the excavation of new crop germplasm resources and offering novel insights into the evolution of feralization in nature.
Collapse
Affiliation(s)
- Zeng-Yuan Wu
- Germplasm Bank of Wild Species & Yunnan Key Laboratory of Crop Wild Relatives Omics, Kunming Institute of Botany, Chinese Academy of Sciences, Kunming, Yunnan 650201, China
| | - Mark A Chapman
- School of Biological Sciences, University of Southampton, Southampton SO17 1BJ, UK
| | - Jie Liu
- Germplasm Bank of Wild Species & Yunnan Key Laboratory of Crop Wild Relatives Omics, Kunming Institute of Botany, Chinese Academy of Sciences, Kunming, Yunnan 650201, China; CAS Key Laboratory for Plant Diversity and Biogeography of East Asia, Kunming Institute of Botany, Chinese Academy of Sciences, Kunming, Yunnan 650201, China.
| | - Richard I Milne
- Institute of Molecular Plant Sciences, School of Biological Sciences, University of Edinburgh, Edinburgh EH9 3JH, UK
| | - Ying Zhao
- Germplasm Bank of Wild Species & Yunnan Key Laboratory of Crop Wild Relatives Omics, Kunming Institute of Botany, Chinese Academy of Sciences, Kunming, Yunnan 650201, China
| | - Ya-Huang Luo
- CAS Key Laboratory for Plant Diversity and Biogeography of East Asia, Kunming Institute of Botany, Chinese Academy of Sciences, Kunming, Yunnan 650201, China
| | - Guang-Fu Zhu
- CAS Key Laboratory for Plant Diversity and Biogeography of East Asia, Kunming Institute of Botany, Chinese Academy of Sciences, Kunming, Yunnan 650201, China
| | - Marc W Cadotte
- Department of Biological Sciences, University of Toronto-Scarborough, Toronto, Ontario, Canada
| | - Ming-Bao Luan
- Institute of Bast Fiber Crops, Chinese Academy of Agricultural Sciences, Changsha, Hunan 410205, China.
| | - Peng-Zhen Fan
- Germplasm Bank of Wild Species & Yunnan Key Laboratory of Crop Wild Relatives Omics, Kunming Institute of Botany, Chinese Academy of Sciences, Kunming, Yunnan 650201, China
| | - Alex K Monro
- Royal Botanic Gardens Kew, Richmond, Surrey TW9 3AE, UK
| | - Zhi-Peng Li
- CAS Key Laboratory for Plant Diversity and Biogeography of East Asia, Kunming Institute of Botany, Chinese Academy of Sciences, Kunming, Yunnan 650201, China
| | - Richard T Corlett
- Royal Botanic Gardens Kew, Richmond, Surrey TW9 3AE, UK; Center for Integrative Conservation and Yunnan Key Laboratory for the Conservation of Tropical Rainforests and Asian Elephants, Xishuangbanna Tropical Botanical Garden, Chinese Academy of Sciences, Menglun, Yunnan 666303, China
| | - De-Zhu Li
- Germplasm Bank of Wild Species & Yunnan Key Laboratory of Crop Wild Relatives Omics, Kunming Institute of Botany, Chinese Academy of Sciences, Kunming, Yunnan 650201, China; CAS Key Laboratory for Plant Diversity and Biogeography of East Asia, Kunming Institute of Botany, Chinese Academy of Sciences, Kunming, Yunnan 650201, China.
| |
Collapse
|
2
|
Habibi N, Al Salameen F, Vyas N, Rahman M, Kumar V, Shajan A, Zakir F, Razzack NA, Al Doaij B. Genome survey and genetic characterization of Acacia pachyceras O. Schwartz. FRONTIERS IN PLANT SCIENCE 2023; 14:1062401. [PMID: 36875582 PMCID: PMC9979705 DOI: 10.3389/fpls.2023.1062401] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 10/05/2022] [Accepted: 01/19/2023] [Indexed: 06/18/2023]
Abstract
Acacia pachyceras O. Schwartz (Leguminoseae), a woody tree growing in Kuwait is critically endangered. High throughput genomic research is immediately needed to formulate effective conservation strategies for its rehabilitation. We therefore, performed a genome survey analysis of the species. Whole genome sequencing generated ~97 Gb of raw reads (92x coverage) with a per base quality score above Q30. The k-mer analysis (17 mer) revealed its genome to be 720Mb in size with an average guanine-cytosine (GC) ratio of 35%. The assembled genome was analyzed for repeat regions (45.4%-interspersed repeats; 9%-retroelements; 2%-DNA transposons). BUSCO assessment of completeness of genome identified 93% of assembly to be complete. Gene alignments in BRAKER2 yielded 34,374 transcripts corresponding to 33,650 genes. Average length of coding sequences and protein sequences were recorded as 1,027nts and 342aa, respectively. GMATA software filtered a total of 901,755 simple sequence repeats (SSRs) regions against which 11,181 unique primers were designed. A subset of 110 SSR primers were PCR validated and demonstrated for its application in genetic diversity analysis of Acacia. The SSR primers successfully amplified A. gerrardii seedlings DNA depicting cross transferability among species. The principal coordinate analysis and the split decomposition tree (bootstrapping runs of 1000 replicates) distributed the Acacia genotypes into two clusters. The flow cytometry analysis revealed the A. pachyceras genome to be polyploid (6x). The DNA content was predicted as 2.46 pg, 1.23 pg, and 0.41 pg corresponding to 2C DNA, 1C DNA and 1Cx DNA, respectively. The results provide a base for further high throughput genomic studies and molecular breeding for its conservation.
Collapse
Affiliation(s)
- Nazima Habibi
- Environment and Life Science Research Centre, Kuwait Institute for Scientific Research, Kuwait, Kuwait
| | - Fadila Al Salameen
- Environment and Life Science Research Centre, Kuwait Institute for Scientific Research, Kuwait, Kuwait
| | - Nishant Vyas
- Department of Immunology, Logical Life Sciences, Pune, India
| | - Muhammad Rahman
- Environment and Life Science Research Centre, Kuwait Institute for Scientific Research, Kuwait, Kuwait
| | - Vinod Kumar
- Environment and Life Science Research Centre, Kuwait Institute for Scientific Research, Kuwait, Kuwait
| | - Anisha Shajan
- Environment and Life Science Research Centre, Kuwait Institute for Scientific Research, Kuwait, Kuwait
| | - Farhana Zakir
- Environment and Life Science Research Centre, Kuwait Institute for Scientific Research, Kuwait, Kuwait
| | - Nasreem Abdul Razzack
- Environment and Life Science Research Centre, Kuwait Institute for Scientific Research, Kuwait, Kuwait
| | - Bashayer Al Doaij
- Environment and Life Science Research Centre, Kuwait Institute for Scientific Research, Kuwait, Kuwait
| |
Collapse
|
3
|
Establishing MinION Sequencing and Genome Assembly Procedures for the Analysis of the Rooibos (Aspalathus linearis) Genome. PLANTS 2022; 11:plants11162156. [PMID: 36015459 PMCID: PMC9416007 DOI: 10.3390/plants11162156] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 06/25/2022] [Revised: 08/08/2022] [Accepted: 08/14/2022] [Indexed: 11/17/2022]
Abstract
While plant genome analysis is gaining speed worldwide, few plant genomes have been sequenced and analyzed on the African continent. Yet, this information holds the potential to transform diverse industries as it unlocks medicinally and industrially relevant biosynthesis pathways for bioprospecting. Considering that South Africa is home to the highly diverse Cape Floristic Region, local establishment of methods for plant genome analysis is essential. Long-read sequencing is becoming standard procedure for plant genome research, as these reads can span repetitive regions of the DNA, substantially facilitating reassembly of a contiguous genome. With the MinION, Oxford Nanopore offers a cost-efficient sequencing method to generate long reads; however, DNA purification protocols must be adapted for each plant species to generate ultra-pure DNA, essential for these analyses. Here, we describe a cost-effective procedure for the extraction and purification of plant DNA and evaluate diverse genome assembly approaches for the reconstruction of the genome of rooibos (Aspalathus linearis), an endemic South African medicinal plant widely used for tea production. We discuss the pros and cons of nine tested assembly programs, specifically Redbean and NextDenovo, which generated the most contiguous assemblies, and Flye, which produced an assembly closest to the predicted genome size.
Collapse
|
4
|
Becher H, Sampson J, Twyford AD. Measuring the Invisible: The Sequences Causal of Genome Size Differences in Eyebrights ( Euphrasia) Revealed by k-mers. FRONTIERS IN PLANT SCIENCE 2022; 13:818410. [PMID: 35968114 PMCID: PMC9372453 DOI: 10.3389/fpls.2022.818410] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 11/19/2021] [Accepted: 06/20/2022] [Indexed: 06/15/2023]
Abstract
Genome size variation within plant taxa is due to presence/absence variation, which may affect low-copy sequences or genomic repeats of various frequency classes. However, identifying the sequences underpinning genome size variation is challenging because genome assemblies commonly contain collapsed representations of repetitive sequences and because genome skimming studies by design miss low-copy number sequences. Here, we take a novel approach based on k-mers, short sub-sequences of equal length k, generated from whole-genome sequencing data of diploid eyebrights (Euphrasia), a group of plants that have considerable genome size variation within a ploidy level. We compare k-mer inventories within and between closely related species, and quantify the contribution of different copy number classes to genome size differences. We further match high-copy number k-mers to specific repeat types as retrieved from the RepeatExplorer2 pipeline. We find genome size differences of up to 230Mbp, equivalent to more than 20% genome size variation. The largest contributions to these differences come from rDNA sequences, a 145-nt genomic satellite and a repeat associated with an Angela transposable element. We also find size differences in the low-copy number class (copy number ≤ 10×) of up to 27 Mbp, possibly indicating differences in gene space between our samples. We demonstrate that it is possible to pinpoint the sequences causing genome size variation within species without the use of a reference genome. Such sequences can serve as targets for future cytogenetic studies. We also show that studies of genome size variation should go beyond repeats if they aim to characterise the full range of genomic variants. To allow future work with other taxonomic groups, we share our k-mer analysis pipeline, which is straightforward to run, relying largely on standard GNU command line tools.
Collapse
Affiliation(s)
- Hannes Becher
- School of Biological Sciences, Institute of Evolutionary Biology, University of Edinburgh, Edinburgh, United Kingdom
| | - Jacob Sampson
- School of Biological Sciences, Institute of Evolutionary Biology, University of Edinburgh, Edinburgh, United Kingdom
| | - Alex D. Twyford
- School of Biological Sciences, Institute of Evolutionary Biology, University of Edinburgh, Edinburgh, United Kingdom
- Royal Botanic Garden Edinburgh, Edinburgh, United Kingdom
| |
Collapse
|
5
|
Dai SF, Zhu XG, Hutang GR, Li JY, Tian JQ, Jiang XH, Zhang D, Gao LZ. Genome Size Variation and Evolution Driven by Transposable Elements in the Genus Oryza. FRONTIERS IN PLANT SCIENCE 2022; 13:921937. [PMID: 35874017 PMCID: PMC9301470 DOI: 10.3389/fpls.2022.921937] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 04/16/2022] [Accepted: 05/16/2022] [Indexed: 05/08/2023]
Abstract
Genome size variation and evolutionary forces behind have been long pursued in flowering plants. The genus Oryza, consisting of approximately 25 wild species and two cultivated rice, harbors eleven extant genome types, six of which are diploid (AA, BB, CC, EE, FF, and GG) and five of which are tetraploid (BBCC, CCDD, HHJJ, HHKK, and KKLL). To obtain the most comprehensive knowledge of genome size variation in the genus Oryza, we performed flow cytometry experiments and estimated genome sizes of 166 accessions belonging to 16 non-AA genome Oryza species. k-mer analyses were followed to verify the experimental results of the two accessions for each species. Our results showed that genome sizes largely varied fourfold in the genus Oryza, ranging from 279 Mb in Oryza brachyantha (FF) to 1,203 Mb in Oryza ridleyi (HHJJ). There was a 2-fold variation (ranging from 570 to 1,203 Mb) in genome size among the tetraploid species, while the diploid species had 3-fold variation, ranging from 279 Mb in Oryza brachyantha (FF) to 905 Mb in Oryza australiensis (EE). The genome sizes of the tetraploid species were not always two times larger than those of the diploid species, and some diploid species even had larger genome sizes than those of tetraploids. Nevertheless, we found that genome sizes of newly formed allotetraploids (BBCC-) were almost equal to totaling genome sizes of their parental progenitors. Our results showed that the species belonging to the same genome types had similar genome sizes, while genome sizes exhibited a gradually decreased trend during the evolutionary process in the clade with AA, BB, CC, and EE genome types. Comparative genomic analyses further showed that the species with different rice genome types may had experienced dissimilar amplification histories of retrotransposons, resulting in remarkably different genome sizes. On the other hand, the closely related rice species may have experienced similar amplification history. We observed that the contents of transposable elements, long terminal repeats (LTR) retrotransposons, and particularly LTR/Gypsy retrotransposons varied largely but were significantly correlated with genome sizes. Therefore, this study demonstrated that LTR retrotransposons act as an active driver of genome size variation in the genus Oryza.
Collapse
Affiliation(s)
- Shuang-feng Dai
- Institution of Genomics and Bioinformatics, South China Agricultural University, Guangzhou, China
| | - Xun-ge Zhu
- Plant Germplasm and Genomics Center, Kunming Institute of Botany, Chinese Academy of Sciences, Kunming, China
| | - Ge-rang Hutang
- Plant Germplasm and Genomics Center, Kunming Institute of Botany, Chinese Academy of Sciences, Kunming, China
| | - Jia-yue Li
- Institution of Genomics and Bioinformatics, South China Agricultural University, Guangzhou, China
| | - Jia-qi Tian
- Institution of Genomics and Bioinformatics, South China Agricultural University, Guangzhou, China
| | - Xian-hui Jiang
- Institution of Genomics and Bioinformatics, South China Agricultural University, Guangzhou, China
| | - Dan Zhang
- College of Tropical Crops, Hainan University, Haikou, China
| | - Li-zhi Gao
- Institution of Genomics and Bioinformatics, South China Agricultural University, Guangzhou, China
- Plant Germplasm and Genomics Center, Kunming Institute of Botany, Chinese Academy of Sciences, Kunming, China
- College of Tropical Crops, Hainan University, Haikou, China
- *Correspondence: Li-zhi Gao,
| |
Collapse
|
6
|
Tomaszewska P, Pellny TK, Hernández LM, Mitchell RAC, Castiblanco V, de Vega JJ, Schwarzacher T, Heslop-Harrison P(J. Flow Cytometry-Based Determination of Ploidy from Dried Leaf Specimens in Genomically Complex Collections of the Tropical Forage Grass Urochloa s. l. Genes (Basel) 2021; 12:957. [PMID: 34201593 PMCID: PMC8306847 DOI: 10.3390/genes12070957] [Citation(s) in RCA: 9] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/26/2021] [Revised: 06/09/2021] [Accepted: 06/16/2021] [Indexed: 01/31/2023] Open
Abstract
Urochloa (including Brachiaria, Megathyrus and some Panicum) tropical grasses are native to Africa and are now, after selection and breeding, planted worldwide, particularly in South America, as important forages with huge potential for further sustainable improvement and conservation of grasslands. We aimed to develop an optimized approach to determine ploidy of germplasm collection of this tropical forage grass group using dried leaf material, including approaches to collect, dry and preserve plant samples for flow cytometry analysis. Our methods enable robust identification of ploidy levels (coefficient of variation of G0/G1 peaks, CV, typically <5%). Ploidy of some 348 forage grass accessions (ploidy range from 2x to 9x), from international genetic resource collections, showing variation in basic chromosome numbers and reproduction modes (apomixis and sexual), were determined using our defined standard protocol. Two major Urochloa agamic complexes are used in the current breeding programs at CIAT and EMBRAPA: the 'brizantha' and 'humidicola' agamic complexes are variable, with multiple ploidy levels. Some U. brizantha accessions have odd level of ploidy (5x), and the relative differences in fluorescence values of the peak positions between adjacent cytotypes is reduced, thus more precise examination of this species is required. Ploidy measurement of U. humidicola revealed aneuploidy.
Collapse
Affiliation(s)
- Paulina Tomaszewska
- Department of Genetics and Genome Biology, University of Leicester, University Road, Leicester LE1 7RH, UK; (T.S.); (P.H.-H.)
| | - Till K. Pellny
- Rothamsted Research, Harpenden, Hertfordshire AL5 2JQ, UK; (T.K.P.); (R.A.C.M.)
| | - Luis M. Hernández
- International Center for Tropical Agriculture (CIAT), Cali 6713, Colombia; (L.M.H.); (V.C.)
| | | | - Valheria Castiblanco
- International Center for Tropical Agriculture (CIAT), Cali 6713, Colombia; (L.M.H.); (V.C.)
| | - José J. de Vega
- Earlham Institute, Norwich Research Park, Norwich NR4 7UZ, UK;
| | - Trude Schwarzacher
- Department of Genetics and Genome Biology, University of Leicester, University Road, Leicester LE1 7RH, UK; (T.S.); (P.H.-H.)
| | - Pat (J.S.) Heslop-Harrison
- Department of Genetics and Genome Biology, University of Leicester, University Road, Leicester LE1 7RH, UK; (T.S.); (P.H.-H.)
| |
Collapse
|
7
|
Phan NT, Orjuela J, Danchin EGJ, Klopp C, Perfus‐Barbeoch L, Kozlowski DK, Koutsovoulos GD, Lopez‐Roques C, Bouchez O, Zahm M, Besnard G, Bellafiore S. Genome structure and content of the rice root-knot nematode ( Meloidogyne graminicola). Ecol Evol 2020; 10:11006-11021. [PMID: 33144944 PMCID: PMC7593179 DOI: 10.1002/ece3.6680] [Citation(s) in RCA: 21] [Impact Index Per Article: 5.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/22/2020] [Revised: 07/04/2020] [Accepted: 07/17/2020] [Indexed: 12/15/2022] Open
Abstract
Discovered in the 1960s, Meloidogyne graminicola is a root-knot nematode species considered as a major threat to rice production. Yet, its origin, genomic structure, and intraspecific diversity are poorly understood. So far, such studies have been limited by the unavailability of a sufficiently complete and well-assembled genome. In this study, using a combination of Oxford Nanopore Technologies and Illumina sequencing data, we generated a highly contiguous reference genome (283 scaffolds with an N50 length of 294 kb, totaling 41.5 Mb). The completeness scores of our assembly are among the highest currently published for Meloidogyne genomes. We predicted 10,284 protein-coding genes spanning 75.5% of the genome. Among them, 67 are identified as possibly originating from horizontal gene transfers (mostly from bacteria), which supposedly contribute to nematode infection, nutrient processing, and plant defense manipulation. Besides, we detected 575 canonical transposable elements (TEs) belonging to seven orders and spanning 2.61% of the genome. These TEs might promote genomic plasticity putatively related to the evolution of M. graminicola parasitism. This high-quality genome assembly constitutes a major improvement regarding previously available versions and represents a valuable molecular resource for future phylogenomic studies of Meloidogyne species. In particular, this will foster comparative genomic studies to trace back the evolutionary history of M. graminicola and its closest relatives.
Collapse
Affiliation(s)
- Ngan Thi Phan
- IRD‐CIRAD‐University of MontpellierUMR Interactions Plantes Microorganismes Environnement (IPME)MontpellierFrance
| | - Julie Orjuela
- IRD‐CIRAD‐University of MontpellierUMR Interactions Plantes Microorganismes Environnement (IPME)MontpellierFrance
| | | | - Christophe Klopp
- Plateforme BioInfo GenotoulUR875INRAECastanet‐Tolosan cedexFrance
| | | | - Djampa K. Kozlowski
- Institut Sophia AgrobiotechINRAECNRSUniversité Côte d’AzurSophia AntipolisFrance
| | | | | | | | - Margot Zahm
- Plateforme BioInfo GenotoulUR875INRAECastanet‐Tolosan cedexFrance
| | | | - Stéphane Bellafiore
- IRD‐CIRAD‐University of MontpellierUMR Interactions Plantes Microorganismes Environnement (IPME)MontpellierFrance
| |
Collapse
|