Danilo Trabuco Amaral1,2, Juliana Rodrigues Bombonato1,2, Sónia Cristina da Silva Andrade3, Evandro Marsola Moraes1, Fernando Faria Franco4. 1. Department of Biology, Center for Human and Biological Sciences, Universidade Federal de São Carlos (UFSCar), Rodovia João Leme dos Santos, Km 110, SP264, Sorocaba, 18052-780, Brazil. 2. Graduate Program in Comparative Biology, Faculty of Philosophy, Sciences and Languages of Ribeirão Preto, Universidade de São Paulo (USP), Ribeirão Preto, Brazil. 3. Department of Genetics and Evolutionary Biology, Instituto de Biociências, Universidade de São Paulo (USP), São Paulo, Brazil. 4. Department of Biology, Center for Human and Biological Sciences, Universidade Federal de São Carlos (UFSCar), Rodovia João Leme dos Santos, Km 110, SP264, Sorocaba, 18052-780, Brazil. franco@ufscar.br.
Abstract
MAIN CONCLUSION: The first South American cactus nuclear genome assembly associated with comparative genomic analyses provides insights into nuclear and plastid genomic features, such as size, transposable elements, and metabolic processes related to cactus development. Here, we assembled the partial genome, plastome, and transcriptome of Cereus fernambucensis (Cereeae, Cactaceae), a representative species of the South American core Cactoideae. We accessed other genomes and transcriptomes available for cactus species to compare the heterozygosity level, genome size, transposable elements, orthologous genes, and plastome structure. These estimates were obtained from the literature or using the same pipeline adopted for C. fermabucensis. In addition to the C. fernambucensis plastome, we also performed de novo plastome assembly of Pachycereus pringlei, Stenocereus thurberi, and Pereskia humboldtii based on the sequences available in public databases. We estimated a genome size of ~ 1.58 Gb for C. fernambucensis, the largest genome among the compared species. The genome heterozygosity was 0.88% in C. fernambucensis but ranged from 0.36 (Carnegiea gigantea) to 17.4% (Lophocereus schottii) in the other taxa. The genome lengths of the studied cacti are constituted by a high amount of transposable elements, ranging from ~ 57 to ~ 67%. Putative satellite DNAs are present in all species, excepting C. gigantea. The plastome of C. fernambucensis was ~ 104 kb, showing events of translocation, inversion, and gene loss. We observed a low number of shared unique orthologs, which may suggest gene duplication events and the simultaneous expression of paralogous genes. We recovered 37 genes that have undergone positive selection along the Cereus branch that are associated with different metabolic processes, such as improving photosynthesis during drought stress and nutrient absorption, which may be related to the adaptation to xeric areas of the Neotropics.
MAIN CONCLUSION: The first South American cactus nuclear genome assembly associated with comparative genomic analyses provides insights into nuclear and plastid genomic features, such as size, transposable elements, and metabolic processes related to cactus development. Here, we assembled the partial genome, plastome, and transcriptome of Cereus fernambucensis (Cereeae, Cactaceae), a representative species of the South American core Cactoideae. We accessed other genomes and transcriptomes available for cactus species to compare the heterozygosity level, genome size, transposable elements, orthologous genes, and plastome structure. These estimates were obtained from the literature or using the same pipeline adopted for C. fermabucensis. In addition to the C. fernambucensis plastome, we also performed de novo plastome assembly of Pachycereus pringlei, Stenocereus thurberi, and Pereskia humboldtii based on the sequences available in public databases. We estimated a genome size of ~ 1.58 Gb for C. fernambucensis, the largest genome among the compared species. The genome heterozygosity was 0.88% in C. fernambucensis but ranged from 0.36 (Carnegiea gigantea) to 17.4% (Lophocereus schottii) in the other taxa. The genome lengths of the studied cacti are constituted by a high amount of transposable elements, ranging from ~ 57 to ~ 67%. Putative satellite DNAs are present in all species, excepting C. gigantea. The plastome of C. fernambucensis was ~ 104 kb, showing events of translocation, inversion, and gene loss. We observed a low number of shared unique orthologs, which may suggest gene duplication events and the simultaneous expression of paralogous genes. We recovered 37 genes that have undergone positive selection along the Cereus branch that are associated with different metabolic processes, such as improving photosynthesis during drought stress and nutrient absorption, which may be related to the adaptation to xeric areas of the Neotropics.