Molecular evidence and ecological niche modeling reveal an extensive hybrid zone among three Bursera species (section Bullockia)

Complete chloroplast genomes of three copal trees (Bursera: Bullockia): comparative analysis and phylogenetic relationships

BACKGROUND

Bursera trees are conspicuous elements of the tropical dry forests in the Neotropics that have significant cultural value due to their fragrant resins (incense), wood sources (handcrafts), and ecological benefits. Despite their relevance, genetic resources developed for the genus are scarce.

METHODS AND RESULTS

We obtained the complete chloroplast (Cp) genome sequence, analyzed the genome structure, and performed functional annotation of three Bursera species of the Bullockia section: Bursera cuneata, B. palmeri, and B. bipinnata. The Cp genome sizes ranged from 159,824 to 159,872 bp in length, including a large single-copy (LSC) region from 87,668 to 87,656 bp, a small single-copy (SSC) from 18,581 to 18,571 bp, and two inverted repeats regions (IRa and IRb) of 26,814 bp each. The three Cp genomes consisted of 135 genes, of which 90 were protein-coding, 37 tRNAs, and 8 rRNAs. The Cp genomes were relatively conserved, with the LSC region exhibiting the greatest nucleotide divergence (psbJ, trnQ-UCC, trnG-UCC, and petL genes), whereas few changes were observed in the IR border regions. Between 589 and 591 simple sequence repeats were identified. Analysis of phylogenetic relationships using our data for each Cp region (LSC, SSC, IRa, and IRb) and of seven species within Burseraceae confirmed that Commiphora is the sister genus of Bursera. Only the phylogenetic trees based on the SSC and LSC regions resolved the close relationship between B. bipinnata and B. palmeri.

CONCLUSION

Our work contributes to the development of Bursera's genomic resources for taxonomic, evolutionary, and ecological-genetic studies.

Development of Genomic Resources in Mexican Bursera (Section: Bullockia: Burseraceae): Genome Assembly, Annotation, and Marker Discovery for Three Copal Species

Bursera comprises ~100 tropical shrub and tree species, with the center of the species diversification in Mexico. The genomic resources developed for the genus are scarce, and this has limited the study of the gene flow, local adaptation, and hybridization dynamics. In this study, based on ~155 million Illumina paired-end reads per species, we performed a de novo genome assembly and annotation of three Bursera species of the Bullockia section: Bursera bipinnata, Bursera cuneata, and Bursera palmeri. The total lengths of the genome assemblies were 253, 237, and 229 Mb for B. cuneata, B. palmeri, and B. bipinnata, respectively. The assembly of B. palmeri retrieved the most complete and single-copy BUSCOs (87.3%) relative to B. cuneata (86.5%) and B. bipinnata (76.6%). The ab initio gene prediction recognized between 21,000 and 32,000 protein-coding genes. Other genomic features, such as simple sequence repeats (SSRs), were also detected. Using the de novo genome assemblies as a reference, we identified single-nucleotide polymorphisms (SNPs) for a set of 43 Bursera individuals. Moreover, we mapped the filtered reads of each Bursera species against the chloroplast genomes of five Burseraceae species, obtaining consensus sequences ranging from 156 to 160 kb in length. Our work contributes to the generation of genomic resources for an important but understudied genus of tropical-dry-forest species.