Shuaibin Wang
Junping Gao
Zhaowu Li2
Kai Chen
Chen Feng- 1Tobacco Research Institute of Technology Centre, China Tobacco Hunan Industrial Corporation, Changsha, China
- 2Puai Medical College, Shaoyang University, Shaoyang, China
- 3Jiangxi Provincial Key Laboratory of ex-situ Plant Conservation and Utilization, Lushan Botanical Garden, Chinese Academy of Sciences, Jiujiang, China
Introduction: Nicotiana L. (Solanaceae) is of great scientific and economic importance, and polyploidization has been pivotal in shaping this genus. Despite many previous studies on the Nicotiana phylogenetic relationship and hybridization, evidence from whole genome data is still lacking.
Methods: In this study, we obtained 995 low-copy genes and plastid transcript fragments from the transcriptome datasets of 26 Nicotiana species, including all sections. We reconstructed the phylogenetic relationship and phylogenetic network of diploid species.
Results: The incongruence among gene trees showed that the formation of N. sylvestris involved incomplete lineage sorting. The nuclear–plastid discordance and nuclear introgression absence indicated that organelle capture from section Trigonophyllae was involved in forming section Petunioides. Furthermore, we analyzed the evolutionary origin of polyploid species and dated the time of hybridization events based on the analysis of PhyloNet, sequence similarity search, and phylogeny of subgenome approaches. Our results highly evidenced the hybrid origins of five polyploid sections, including sections Nicotiana, Repandae, Rusticae, Polydicliae, and Suaveolentes. Notably, we provide novel insights into the hybridization event of section Polydicliae and Suaveolentes. The section Polydicliae formed from a single hybridization event between maternal progenitor N. attenuata and paternal progenitor N. undulata; the N. sylvestris (paternal progenitor) and the N. glauca (maternal progenitor) were involved in the formation of section Suaveolentes.
Discussion: This study represents the first exploration of Nicotiana polyploidization events and phylogenetic relationships using the high-throughput RNA-seq approach. It will provide guidance for further studies in molecular systematics, population genetics, and ecological adaption studies in Nicotiana and other related species.
1 Introduction
Nicotiana L. (Solanaceae) is of great scientific and economic importance, containing the cultivated tobaccos (N. tabacum and N. rustica), the model plant (N. benthamiana), as well as some essential ornamentals (e.g., N. alata and N. sylvestris) (Wang and Bennetzen, 2015). The genus Nicotiana comprises about 87 species, including one recently reported new Australian species in the section Suaveolentes (N. paulineana) (Bally et al., 2021), nearly half of which are allotetraploids (Knapp et al., 2004). The classification of Nicotiana mainly relied on geographical distribution, morphological characters, and cytological investigations, which were first reported by Goodspeed (1956). It was subsequently updated by Knapp et al. (2004) based on the phylogenetic analysis and morphological description. The current classification of Nicotiana comprises three subgenera (Rustica, Tabacum, and Petunioides) and 13 sections, five of which contain polyploids formed by interspecific hybridization (seven diploid sections: Alatae, Noctiflorae, Petunioides, Undulatae, Paniculatae, Trigonophyllae, Tomentosae, and five polyploidy sections: Suaveolentes, Repandae, Nicotiana, Polydicliae, Rusticae) (Knapp et al., 2004; Leitch et al., 2008). The Nicotiana species are distributed across tropical and temperate regions and are primarily endemic to South America, North America, and Australia (Knapp et al., 2004), of which N. tabacum is one of the most widely cultivated non-food crops, having been spread worldwide by humans. About 75% of Nicotiana species occur naturally in America and 25% in Australia (Aoki and Ito, 2000; Clarkson et al., 2004). Interestingly, all native Australian Nicotiana species belong to section Suaveolentes (Bally et al., 2021). In addition, Nicotiana species exhibit a spectacular range of floral morphology and color, genome size, and karyotypic diversity (Leitch et al., 2008; Marks et al., 2011; Renny-Byfield et al., 2013; McCarthy et al., 2015; McCarthy et al., 2016) and the polyploids of Nicotiana formed at different stages of evolutionary divergence (Leitch et al., 2008). The genus Nicotiana is, therefore, an excellent system in which to take advantage of recent advances in the research of speciation, biodiversity, and phytogeography (Aoki and Ito, 2000; McCarthy et al., 2016).
Phylogenetic relationships have been the subject of study in this genus for around two decades based on the plastid markers (coding and noncoding) (Aoki and Ito, 2000; Clarkson et al., 2004), low-copy nuclear genes (Kitamura et al., 2001; Clarkson et al., 2017), nuclear ribosomal internal transcribed spacer (ITS) (Chase et al., 2003) and random amplified polymorphic DNA (RAPD) analysis (Khan and Narayan, 2007). The previous phylogenetic studies have provided new insights into interspecific relationships (Clarkson et al., 2004) and also led to a modification in the traditional classification of the genus Nicotiana (Knapp et al., 2004). However, these phylogenetic analyses are based on only several molecular markers or short sequences (glutamine synthetase gene, leafy/floricaula gene, ITS, trnL-F, trnS-G, ndhF, and matK) and provide limited resolution of relationships among Nicotiana species. And the deep relationships of Nicotiana have usually neither been resolved nor well-supported. Recent analyses using the complete chloroplast genomes of Nicotiana in 11 sections have recovered a nearly fully resolved phylogenetic relationship and deduced a potential maternal progenitor of polyploid species (Wang et al., 2022). However, evidence of phylogenetic relationships and diploid progenitors of polyploid species from the nuclear are still lacking. Therefore, sufficient evidence from whole genome data is needed to deduce the deep phylogenetic relationships and demonstrate the phylogenetic discordance among this genus.
Polyploidization has been pivotal in shaping this genus (Leitch et al., 2008; McCarthy et al., 2016). Approximately half of the Nicotiana species were considered natural tetraploid species of different ages (Goodspeed, 1956; Knapp et al., 2004; Leitch et al., 2008). The majority of the Nicotiana species possess 12 or 24 chromosome pairs, except for several diploid species in section Alatae with 9 or 10 pairs and several polyploidy species in section Suaveolentes with 15, 16, 18, 19, 20, 21, or 22 pairs (Khan and Narayan, 2007; Marks et al., 2011). So far, the morphological, distributional, cytogenetic, and molecular evidence has been used to discover the diploid progenitors for each tetraploid species (Chase et al., 2003; Lim et al., 2005; Kelly et al., 2013; Schiavinato et al., 2020). One of the most studied polyploids in this genus was the tetraploid N. tabacum (2n = 4x = 48), with different lines of evidence suggesting N. sylvestris (2n = 2x = 24, maternal donor) and N. tomentosiformis (2n = 2x = 24, paternal donor) as candidate parents based on morphological observation (Goodspeed, 1956), plastid genome comparison (Yukawa et al., 2006), RAPD analysis (Khan and Narayan, 2007), the genome size (Leitch et al., 2008), and whole genome sequencing (Sierro et al., 2014; Sierro et al., 2018). In addition, the species in section Repandae were proposed to have a hybrid origin between N. sylvestris and N. obtusifolia based on a phylogenetic context (Leitch et al., 2008). The species in the section Rusticae were suggested as hybrids between the ancestral species of N. paniculata of section Paniculatae and N. undulata of section Alatae based on the comparison of karyotype and genome size (Khan and Narayan, 2007; Leitch et al., 2008). The ancestral species of section Suaveolentes possibly were related to sections Acuminatae (which should be called Petunioides, following Knapp et al., 2004), Noctiflorae, and Alatae based on external morphology (Khan and Narayan, 2007). Still, this hypothesis has never been formally tested (Leitch et al., 2008), and the ancestry of allopolyploid species in section Polydicliae is unresolved (Khan and Narayan, 2007; Leitch et al., 2008). Thus, a genome-wide perspective on the origin and evolution of allopolyploid species, including estimation of divergence dates, has been lacking.
The rise of high-throughput sequencing techniques has produced massive amounts of genomic or transcriptomic data, providing an unprecedented opportunity for systematic and evolutionary studies in great depth (Lemmon and Lemmon, 2013). Notably, because of the relatively low cost of transcriptome sequencing compared with genome sequencing and the fact that phylotranscriptomics is almost as reliable as phylogenomics (Cheon et al., 2020), phylotranscriptomic analysis has emerged as the preferred method for studying evolutionary biology (Cheon et al., 2021). Recent studies based on nuclear genes, especially phylotranscriptomics, have been successful in resolving relationships of various scales from the genus (Yang et al., 2018; Chen et al., 2021) to angiosperm-wide (Yang et al., 2020), even gymnosperms (Liu et al., 2022) and ferns (Qi et al., 2018) plants. During the last decade, large amounts of transcriptomic data have been generated for this genus, which has provided new opportunities for studying the phylogenetic relationships and evolution of polyploids at the scale of whole genomes. The newer approaches based on transcriptomic studies applying to Nicotiana species will provide a more accurate evaluation of speciation and polyploid events of Nicotiana (Clarkson et al., 2017).
Nicotiana species have significant economic importance. Most notably, tobacco (N. tabacum) is a major cash crop widely used in the production of tobacco products (Wang and Bennetzen, 2015). Additionally, Nicotiana species have pharmaceutical and research value. They are extensively used in the preparation of medicines, including both traditional herbal remedies and modern pharmaceuticals (Wang and Bennetzen, 2015). Moreover, due to the rich genomic and genetic diversity of Nicotiana species, they serve as ideal model organisms for studying genetic engineering, molecular biology, and genetics (Wang and Bennetzen, 2015).
In this study, we obtained the transcriptome datasets from 26 Nicotiana species, including 17 diploid species in eight sections and nine allopolyploid species in five sections. Two outgroup species (Petunia axillaris and Petunia inflata) from the genus Petunia were used as the sister taxa in the family of Solanaceae. We performed the analysis of phylogeny and PhyloNet based on the low-copy nuclear genes and transcript fragments of plastid genomes, respectively. Our study aimed to address the following topics: (1) to re-examine the classification and phylogenetic relationships reported in previous studies of Nicotiana and provide a relationship strongly supported among the lineages of Nicotiana; (2) to assess the conflict between nuclear and plastid phylogenetic topology and the inconsistency of gene trees; (3) to investigate the potential parental origin for the species of five polyploid sections; (4) to estimate the divergence time of diploid species and the time of the interspecific hybridization events that gave rise to polyploid species. This study represents the first exploration of Nicotiana phylogeny and timing of diversification in allopolyploids utilizing the high-throughput RNA-seq approach. It will provide valuable insights for future research in molecular systematic, population genetic, and ecological adaptation studies in Nicotiana and other related species.
2 Materials and methods
2.1 Data source, reads trimming, and transcriptome assembly
Here the transcriptome datasets were obtained from 26 species of genus Nicotiana including 17 diploid Nicotiana species and nine Nicotiana allotetraploid species representing each of the 13 Nicotiana sections: Alatae (n = 12), Nicotiana (n = 24), Noctiflorae (n = 12), Paniculatae (n = 12), Petunioides (Noctiflorae) (n = 12), Polydicliae (n = 24), Repandae (n = 24), Rusticae (n = 24), Suaveolentes (n = 16-24), Sylvestres (n = 12), Tomentosae (n = 12), Trigonophyllae (n = 12), and Undulatae (n = 12) (Clarkson et al., 2004). These transcriptome sequencing reads were retrieved from the NCBI Sequence Read Archive (SRA) database using the fastq-dump software of the SRA Toolkit package (available from https://github.com/ncbi/sra-tools). Additionally, two species of the genus Petunia (P. axillaris and P. inflata) from the non-Nicotiana Solanaceae family were included as outgroups. The protein-coding sequences of P. axillaris and P. inflata (Bombarely et al., 2016) were retrieved from Solanaceae Genomics Network (available from https://solgenomics.net/).
For transcriptome analyses, the quality control of raw data was performed using the fastp v0.23.2 (Chen et al., 2018), and the adapter, short reads (min. read length 50 bp), reads containing N, and reads with low-quality score (min. quality 20) were removed. All subsequent analyses were based on the filtered clean data. The clean data was first aligned into the genome sequences of allotetraploid N. tabacum (version 4.5, available from https://solgenomics.net/ftp/genomes/Nicotiana_tabacum/edwards_et_al_2017/) (Edwards et al., 2017) using STAR v2.7.10a (Dobin et al., 2012) with the default parameters. The aligned BAM files from the same species were merged, sorted, and indexed using samtools v1.15.1 (Danecek et al., 2021). Then, the quality of BAM files were assessed using samtools with the ‘flagstat’ and ‘depth’ algorithms. The Trinity v2.14.0 (Grabherr et al., 2011) was used to perform genome-guided de novo transcriptome assembly based on the BAM file of each species with default parameters. The BUSCO v5.3.2 (Manni et al., 2021) was applied to evaluate the completeness of the assembly results based on the solanales homologous gene database (available from https://busco-data.ezlab.org/v5/data/lineages/). The TransDecoder v5.5.0 (available from http://transdecoder.github.io) was employed to predict the open reading frame (ORF) within the assemblies with a minimum protein length of 100 bp. A BLAST search against the Uniref90 proteins database (Suzek et al., 2007) using BLAST (Camacho et al., 2009) and an HMMER search against the Pfam protein domain database v35.0 to identify common protein domains using PfamScan (Mistry et al., 2020). TransDecoder leveraged the outputs generated above to ensure that those peptides with blast hits or domain hits were retained in the set of reported likely coding regions. Then, CD-HIT v4.8.1 (Huang et al., 2010) was applied to remove redundant sequences with a threshold value of 0.95 identities. The statistic of assembly result was evaluated using ‘TrinityStats.pl’ script in the Trinity utility (Grabherr et al., 2011).
2.2 Identification and filtering of orthogroups
The OrthoFinder v2.5.4 (Emms and Kelly, 2015) was employed to infer orthologous genes from the predicted non-redundant CDS of the Nicotiana species and outgroups with default settings. Its utility has been tested in previous phylogenomic studies of plants (Peng et al., 2020) and animals (Fukushima and Pollock, 2020). To further increase the robustness of phylogenetic analyses, we filtered the orthologous groups according to the following criteria: (1) the max gene number of each sample in any orthologous group < 5; (2) the max gene number of each diploid sample in any orthologous group < 3; (3) the sequences of coding regions with the length >300 bp; (4) coverage of at least 50% of the 26 samples in each orthologous group. Thus, a total of 995 low-copy orthologous were generated to investigate Nicotiana’s phylogenetic relationships. The gene number and the presence/absence of orthologous groups in each species were shown as a heatmap using Python script.
A significant difficulty in reconstructing phylogeny was that polyploidization events involved many species across the genus Nicotiana, resulting in hidden paralogs (remaining single-copy genes after the loss of distinct paralogs in different taxa), which should be avoided in phylogenetic analysis. Thus, low-copy orthologous from only the diploid species were selected to build the phylogenetic backbone structure. The low-copy orthologous groups, including polyploid species, were used to construct the phylogenetic network and infer the hybrid process.
2.3 Diploid gene tree and species tree reconstruction
The polyploid samples were excluded from the filtered low-copy orthologous groups. Thus, 17 Nicotiana ploidy species and two Petunia species as outgroups with one-to-one single-copy orthologous were extracted for phylogenetic analyses. The concatenation- and coalescent-based methods were applied to reconstruct Nicotiana’s gene and species tree, respectively. For the concatenation-based method, individual single-copy orthologous were aligned using MAFFT v7.490 (Katoh and Standley, 2013) with the ‘L-INS- I’ algorithm, and poorly aligned regions were removed using Gblock v0.91b (Talavera and Castresana, 2007) with parameters ‘-b4 = 5 -b5=h -t=d’. All trimmed alignments were concatenated as a supermatrix and then performed summary statistics using AMAS (Borowiec, 2016). Two conventional approaches were used to construct a phylogenetic tree: (1) A maximum-likelihood (ML) analysis was performed using RAxML-NG v1.1.0 (Kozlov et al., 2019), with 1,000 bootstraps to find the best-scoring ML tree and the best substitution model GTR+I+G4 under the Akaike information criterion (AIC) calculated by ModelTest-NG v0.1.7 (Darriba et al., 2019). (2) A Bayesian-inference (BI) analysis was performed using MrBayes v3.2.7a (Ronquist et al., 2012), with the parameters: nst = 6 and rates = gamma. Four independent Markov chain Monte Carlo (MCMC) chains were run for 1,000,000 generations with random initial trees. Trees were sampled per 100 generations. The first 25% of trees were discarded as burn-in, with the remaining trees being used for generating the consensus tree. Tracer v1.7.1 (available from http://beast.community/tracer) was used to assess the quality of the MCMC simulations and the stability of runs. For the coalescent-based method, a massively parallel tool ParGenes v1.2.0 (Morel et al., 2018; Zhang et al., 2018; Kozlov et al., 2019) was used for model selection. And species tree inference based on the single-copy orthologous with the parameters applied as: ‘pargenes.py -a coalescence -o pargenes -c 32 -d nt -m –use-astral -b 1000’. Finally, the concatenated and coalescent phylogenetic results were visualized using FigTree v1.4.4 (available from http://tree.bio.ed.ac.uk/software/figtree/).
2.4 Identification of conflict and concordance among gene trees
Phyparts v0.0.1 (Smith et al., 2015) was employed to summarize the conflict and concordance information by comparing the gene tree of each single-copy orthologous against the species tree described above. Specifically, the gene tree of each single-copy orthologous with bootstrap support (BS) values was exported from the outputs of ParGenes and optimally rooted with the outgroups (P. axillaris and P. inflata) using phyx v1.3 (Brown et al., 2017). Still, in cases where the outgroups were missing, the gene trees were rooted by Minimal Ancestor Deviation (MAD) using MADroot (Tria et al., 2017). The branches in each gene tree with less than 33% BS were considered uninformative and collapsed using Newick utilities (Junier and Zdobnov, 2010). Two hundred four remaining gene trees and the species tree were used as inputs for phyparts to summarize the conflict and concordance information. The results of support and conflict statistics between gene trees and species trees were visualized with phypartspiecharts.py (available from https://github.com/mossmatters/phyloscripts/).
Furthermore, to quantify branch support values for the species tree, Quartet Sampling v1.3.1 (QS) (Pease et al., 2018) analysis was conducted with 1,000 replicates, and the log‐likelihood cutoff was 2. QS was a method to analyze molecular phylogeny by calculating branch support using repeated sampling of quartets. It provided four values in the outputs: QC (the Quartet Concordance score), QD (the Quartet Differential score), QI (the Quartet Informativeness score), and QF (the Quartet Fidelity score). The QS outputs were visualized as a figure with plot_QC_ggtree.R (available from https://github.com/ShuiyinLIU/QS_visualization). The results from phyparts and QS could provide alternative evidence for evaluating the discordance of gene trees.
2.5 Phylogenetic network analysis of diploid species
Nuclear-plastid discordance indicated a hybrid origin for section Petunioides, comprising N. acuminata, N. miersii, and N. attenuata (see Results). Moreover, the discordance of species and gene trees also suggested a hybridization or incomplete lineage sorting (ILS) for the section Sylvestres (see Results). We then used the maximum pseudolikelihood (MPL) approach implemented in PhyloNet v3.8.2 (Wen et al., 2018) to assess corroborative evidence supporting these conclusions. A total of 596 gene trees covering all 17 Nicotiana ploidy species and two outgroups were used to infer the phylogenetic network with the command ‘InferNetwork_MPL’ (Yu and Nakhleh, 2015). The section Sylvestres and section Petunioides were set as a potential hybrid, respectively. Uncertain nodes were bypassed in the gene trees by applying a bootstrap support threshold of 30 using the -b flag. The networks were visualized using Dendroscpoe v3.82 (Huson and Scornavacca, 2012).
In addition, the NeighborNet method implemented in SplitsTree v4.18.2 (Huson and Bryant, 2005) were used to reconstruct phylogenetic networks for the concatenated dataset. The K2P model was used for distance analysis, and support values at each node were estimated by running 1,000 bootstrap replicates.
2.6 Inference on the origin of allopolyploid species
Two strategies were used to infer the hybridization process within polyploid species of Nicotiana: (1) Inference based on the phylogenetic network. The phylogenetic network of each allopolyploid species was carried out using PhyloNet v.3.8.2 (Wen et al., 2018) with the command ‘InferNetwork_MP_Allopp’ under the MDC criterion (Yan et al., 2021). Network searches were performed using only nodes in the gene trees with BS support of at least 30%, allowing for one hybridization event. The networks were visualized using Dendroscpoe v3.82 (Huson and Scornavacca, 2012). (2) Inference based on the phylogenetic relationship of subgenomes of the allopolyploid. The first step was to split the CDS sequences of each allopolyploid species into two subsets based on the sequence similarity between polyploid and all diploid species. Specifically, each low-copy orthologous group was used to perform a sequence similarity search using BLAST+ (Camacho et al., 2009) against itself with the default parameters. The diploid species with the best hit of the CDS sequences in each allopolyploid species were recorded using the custom Python script from the blast tabular output. The total number of diploid species with the best hit in the CDS sequences of each polyploid species was counted and visualized as a heatmap. The diploid species most closely related to the progenitors of each allopolyploid would yield the highest number of best hits. Therefore, we sought to identify the diploid progenitors’ representatives for each allopolyploid. Then, each CDS sequence of allopolyploid species was assigned to the maternal or paternal group based on the putative results from the heatmap and previous PhyloNet analysis. The CDS sequences with ambiguous classification were discarded. The second step was to reconstruct the phylogenetic relationship based on the subgenomes of the polyploid species. Specifically, the phylogenetic relationship was inferred based on the concatenated datasets of the CDS sequences of diploid species and the classified CDS sequences of polyploid species using the ML analysis method described above. In addition, the inferred hybridization among the section Sylvestres and section Petunioides was also evaluated based on the sequence similarity search in the second strategy described above.
2.7 Plastid transcript assembly and phylogenetic analysis
Plastid genomes inherit maternally in the genus Nicotiana (Svab and Maliga, 2007). Thus the plastid phylogenetic tree could be used to determine the maternal progenitor of tetraploid species. Concretely, the clean reads were mapped into a previously published plastid genome (N. sylvestris: NC_007500.1) (Yukawa et al., 2006) with only one copy of inverted repeat regions using STAR v2.7.10a (Dobin et al., 2012), following the unmapped reads were filtered using samtools v1.15.1 (Danecek et al., 2021) with the parameter: -F 12. Genome-guided de novo transcriptome assembly was performed based on the filtered BAM file of each species. CD-HIT v4.8.1 (Huang et al., 2010) was used to exclude similar sequences with the parameter: -c 0.95. The local collinear regions among the filtered assemblies were identified using Mugsy v1.2.3 (Angiuoli and Salzberg, 2011). Only the conserved collinear regions with a length > 100 bp and coverage of at least 50% of the samples were extracted from the output of Mugsy. The filtered collinear regions were aligned and filtered according to the above-mentioned method. Then the combined supermatrix was used to perform phylogeny analysis using RAxML-NG v1.1.0 (Kozlov et al., 2019) with 1,000 bootstraps.
2.8 Dating of the divergence time and interspecific hybridization event
The divergence times within the species tree were inferred using BEAST v.2.6.7 (Bouckaert et al., 2019) optimized for OpenGL graphics. The dating analysis of concatenated single-copy genes of diploid species and the classified genes of allopolyploid species was performed with a strict clock, HKY substitution model, gamma site heterogeneity model, estimated base frequencies, and an ML starting tree. A Calibrated Yule model was specified as the tree prior. As no reliable fossils were available to calibrate the internal nodes of the Nicotiana, one secondary calibration from a recently published dated phylogeny of the Solanaceae (Särkinen et al., 2013) was used to calibrate the crown age of Solanaceae at 30.4 (95% HPD 26.3-34.0) million years ago (Ma). Two independent MCMC analyses of 10 million generations with 10% burn-in and sampling every 1,000 generations were conducted to evaluate the credibility of posterior distributions of parameters. The log files from BEAST were analyzed with Tracer v.1.7.0 (available from http://beast.community/tracer) to evaluate and ensure convergence, effective sample size (ESS) values, density plots, and trace plots. A maximum clade credibility tree with median heights was generated with TreeAnnotator v.1.8.4 (Bouckaert et al., 2019). The final tree was visualized using FigTree v1.4.4.
3 Results
3.1 RNAseq data, transcriptomes assembly, and ortholog identification
Here, a total of ~300 Gb transcriptome data from 26 Nicotiana species were collected, covering all 13 Nicotiana sections (Table 1). After read trimming, the number of clean reads per species ranged from 29.4 million to 454.6 million, with an average of 139.2 million.
Table 1 Summary information of the 26 Nicotiana species and the genome-guided transcriptome assemblies.
For the genome-guide transcript assembly, the clean reads of each species were mapped to the N. tabacum genome, which resulted in an average of 81.6% primary mapped reads, with the fewest in N. otophora (66.9%) and the most in N. tabacum (93.8%) (Table 1). These assemblies produced between 77,747 and 323,807 contigs (≥ 300 bp) for each species, with an N50 length ranging from 757 bp to 2,756 bp and a total length ranging from 45.7 Mb to 337.8 Mb (Table 1). A total of 26,388 to 71,557 unigenes were detected in each species, with N50 ranging from 897 to 1,515 bp and GC content from 42.4% to 43.4%. The assembly completeness evaluation results showed that all assembled unigenes had relatively high BUSCO scores, ranging from 50.0% (N. obtusifolia) to 96.6% (N. attenuata) among the 26 species (Table 1). These values, below a fully satisfactory BUSCO score, could be explained by the absence of tissue diversity and (or) enough data. Nevertheless, we considered that these genome-guided assembled gene sets would provide a reasonably good representation of the transcriptomes of Nicotiana species. For the plastid transcript assembly, the clean reads of each species were mapped to the plastid genome of N. sylvestris. The unmapped reads were filtered from the mapping results subsequently, which resulted in an average of 0.22 million mapped reads and an average of 85.6% reads coverage rate (>5x), with the fewest in N. petunioides (68.32%) and the most in N. velutina (99.7%) (Supplementary Table S1). After the transcript assemblies, ORF prediction, and de-redundancy, there were between 12 and 144 unigenes detected in each species, with N50 ranging from 873 to 19,600 bp and GC content ranging from 36.50% to 38.97% (Supplementary Table S1).
Based on the combined datasets from the unigenes of Nicotiana species and the outgroups, a total of 73,634 orthologous groups were identified, of which 995 low-copy genes passed our filtering criteria as mentioned above. Among these low-copy orthologous groups, the number of genes in each species ranges from 752 to 1,185. A total of 776 orthologous groups contained at least one of the outgroups, 558 shared in all 26 Nicotiana species, and 510 shared in all 28 species (26 Nicotiana specie and two outgroups) (Supplementary Figure S1). All the filtered low-copy genes were used to perform the phylogeny analysis, detect the parental progenitor, and date the hybridization event.
3.2 Diploid phylogenetic inference
We obtained the sequences of 995 genes with at least 273 bp in length for the nuclear phylogenetic dataset of diploid species. The alignment length for these single-copy genes ranged from 249 to 5,192 bp, with a mean length of 813 bp. After concatenation, the aligned 995-gene super matrix reached 808,952 bp in length, with 138,331 (17.1%) variable sites, 85,084 (10.5%) parsimony informative sites, and 154,142 (19.1%) missing sites (gaps and undetermined characters) (Supplementary Table S2).
For the nuclear phylogenetic dataset, the phylogenetic relationships reconstructed by ML and BI methods had identical topologies that separated the genus into two major clades with strongly supported (MLBS=100, PP=1.0) (Figure 1A). The first of them was integrated by a clade of section Tomentosae, sister of section Trigonophyllae, and a clade where section Paniculatae was sister to the clade of section Undulatae. The second was integrated by a clade, where section Alatae and Sylvestres were recovered as successive sister species of the clade of section Petunioides plus section Noctiflorae (Figure 1A). The coalescent-based species tree inferred from the diploid nuclear dataset yielded a concordant phylogenetic relationship (ASTRAL LPP=1.00) with the concatenation analysis, except for the N. sylvestris clade with a slightly less confident (ASTRAL LPP=0.94) (Figure 1A).
Figure 1 Tanglegram illustrating the nuclear-plastid discordance and the conflict signal among the gene trees in Nicotiana diploid species. (A) The diploid nuclear phylogeny recovered from concatenation- or coalescent-based methods based on the 995 single-copy genes (left). (B) The diploid plastid phylogeny recovered from the maximum-likelihood method based on the concatenated local collinear regions of transcript fragments (right). Bootstrap percentages were indicated beside the branches, and only values less than 100 were shown. (C) Patterns of gene-tree concordance and conflict of Nicotiana based on the phyparts analysis (left). The tree topology used was inferred by ASTRAL. The pie charts at each node show the proportion of genes in concordance (blue), conflict (green: a single dominant alternative; red: all other conflicting trees), and without enough information (gray). The numbers above and below each branch were the numbers of concordant and conflicting genes at each bipartition, respectively. (D) Information on Nicotiana’s gene-tree concordance and conflict based on the result from quartet sampling (right). Branch labels show quartet concordance (QC), quartet differential (QD), and quartet informativeness (QI), respectively, for each relationship. Corresponding clades in the nuclear and plastid phylogeny were colored. Asterisks on the branches of nuclear phylogeny (left) indicated local posterior probabilities of 0.94 in the coalescent-based species tree. Sections were classified according to Knapp et al. (2004) and labeled to the left.
For the plastid phylogenetic dataset, the concatenated transcript fragment from plastid genomes reached 45,728 bp in length, with 2,103 (4.6%) variable sites, 1,125 (2.5%) parsimony informative sites, and 4,448 (9.72%) missing sites (Supplementary Table S2). The diploid plastid genomic phylogeny showed that the genus Nicotiana was monophyletic. Unlike nuclear phylogeny, within this plastid phylogenetic topology, three main clades with strong support were recovered (Figure 1B). The first was section Tomentosae as the basal-most clade. A clade of section Petunioides integrated the second, sister of section Trigonophyllae, and a clade of section Paniculatae, sister of section Undulatae. The third was integrated by a clade, where section Alatae and Sylvestres were recovered as successive sister species of section Noctiflorae.
Phylogenetic discordance was observed between topologies inferred from the nuclear and plastid concatenated datasets (Figures 1A, B). In nuclear phylogeny, section Petunioides (N. acuminata, N. miersii, and N. attenuata) were strongly supported (MLBS = 100) to be sister to section Noctiflorae (N. noctiflora, N. petunioides, and N. glauca). Section Petunioides were placed as sister to section Trigonophyllae (N. obtusifolia) with strongly supported (MLBS = 92) in plastid phylogeny. This result showed the hybridization of section Petunioides between section Noctiflorae and Trigonophyllae, which might explain the discordance of nuclear and plastid phylogenetic topologies. In addition, the section Tomentosae (N. otophora and N. tomentosiformis) was recovered as a basal taxon in plastid phylogeny. In construction, two main groups were recovered within the nuclear phylogeny.
3.3 Gene tree concordance and conflict
Our species tree inferred from coalescent nuclear data suggested that two major clades were identified in the phylogeny of the genus Nicotiana with full support (node A). Likewise, the result of Quartet Sampling (QS) demonstrated that node A was confirmed with full support (1/-/1; i.e., all informative quartets support that lineage) and phyparts result with almost all the informative gene trees (750 out of 751) support this topology (Figures 1C, D). Although the nuclear phylogenetic tree confirmed N. sylvestris sister to species of the clade of section Petunioides plus section Noctiflorae with full support, this clade (node B) was supported by only 28% of informative quartets with a skewed frequency for alternative discordant topologies (QS score = 0.28/0/1). Likewise, the result of phyparts supported this clade with only 148 out of the 716 informative gene trees (20.7%). This result revealed that ILS or hybridization might explain this phylogenomic discordance.
All the conflict nodes between nuclear and plastid phylogenetic topology showed relatively low informative gene trees supported in the phyparts result (Figure 1C). For example, nodes C and D have alternative discordant topologies. The result of phyparts supported this clade with only 275 of the 739 informative gene trees (37.2%) and 253 of the 655 informative gene trees (38.6%), respectively. In contrast, the QS result demonstrated that all these two nodes related to three sections were confirmed with full support (1/-/1) (Figure 1D).
3.4 Diploid phylogenetic networks
The phylogenetic discordances of nuclear-plastid and gene trees showed a complex evolutionary history among N. sylvestris and the species of section Petunioides. Using the N. sylvestris as potential hybrid species, our PhyloNet analysis indicated that N. sylvestris was a potential hybrid species between the section Alatae and the common ancestor of section Noctiflorae and Petunioides, or at least introgressed with section Alatae (Figure 2A). When using the species of section Petunioides as potential hybrid species, our PhyloNet analysis did not support the hybrid origin of N. sylvestris and the species of section Petunioides (data not shown). The analysis of sequence similarity search in N. sylvestris showed that a total of 264 genes in section Petunioides, 119 in section Noctiflorae, and 249 in section Alatae were similar to the gene sequences in N. sylvestris with the best hit (Figure 2B). In the sequence similarity search of section Petunioides, the best hit mainly occurred in N. sylvestris, section Noctiflorae, and Alatae, not in the section Trigonophyllae (Figure 2B).
Figure 2 The phylogenetic network analysis of diploid species. (A) Phylogenetic networks were inferred by setting N. sylvestris as the hybrid species using the InferNetwork_ML method in PhyloNet. Blue branches indicated lineages involved in reticulated histories, and numerical values were the inheritance probabilities for each reticulation. (B) The heatmap showed the best hit number of the CDS sequences in N. sylvestris and section Petunioides against the other diploid species. (C) Phylogenetic networks generated by NeighborNet method in SplitTree4.
Phylogenetic networks reconstructed using NeighborNet revealed apparent clustering among the Nicotiana diploid sections. When rooted in the genus Petunia, the topology of phylogenetic networks (Figure 2C) was very similar to the ML and Bayes trees (Figure 1A). However, N. sylvestris was alternative splits connecting to the section Alatae or the common ancestor of section Petunioides and Noctiflorae (Figure 2C), which showed uncertainty regarding the phylogenetic placement of this clade.
3.5 Putative diploid progenitors of polyploid species
Based on the InferNetwork_MP_Allopp approach implemented in PhyloNet, we inferred the allopolyploid network from gene trees under the MDC criterion (Figures 3A–E). For the N. tabacum hybridization test, the PhyloNet result showed a hybridization event between the N. sylvestris and N. tomentosiformis (Figure 3A). For the N. stocktonii, the PhyloNet result showed a hybridization event between the N. obtusifolia and N. sylvestris (Figure 3B). For the N. rustica, the PhyloNet result showed a hybridization event between the N. knightiana and N. undulata (Figure 3C). For the N. clevelandii, the PhyloNet detected a hybridization between the N. attenuata and N. undulata into the clade N. clevelandii (Figure 3D). For the hybridization tests of the N. velutina, N. cavicola, N. rosulata, N. benthamiana, and N. amplexicaulis, the PhyloNet detected a hybridization between the N. glauca and N. sylvestris into the clade of section Suaveolentes (Figure 3E).
Figure 3 The diploid origin analysis of polyploid species. (A–E) Phylogenetic networks were inferred by setting each polyploid species as the hybrid species using the InferNetwork_MP_Allopp method in PhyloNet. Blue branches indicated lineages involved in reticulated histories. (F) The heatmap showed the best hit number of the CDS sequences in the polyploid species against the other diploid species. (G) The nuclear phylogeny recovered from the maximum-likelihood method based on the classified low-copy genes. Bootstrap percentages were indicated beside the branches, and only values less than 100 were shown. Sections were classified according to Knapp et al. (2004) and labeled to the right. ♀: maternal origin; ♂: paternal origin.
For each polyploid species, we counted the number of diploid species with the best hit of the CDS sequences (Figure 3F). In the CDS sequences of tetraploid N. tabacum, 341 and 456 CDS sequences had the best hit against the N. sylvestris and N. tomentosiformis, respectively, which suggested that N. sylvestris and N. tomentosiformis were the putative diploid progenitors of N. tabacum. Similarly, in tetraploid N. Stocktonii, a total of 217 and 241 CDS sequences against the N. sylvestris and N. obtusifolia with the best hit, respectively, suggested that N. sylvestris and N. obtusifolia were the putative diploid progenitors of N. Stocktonii. In the tetraploid N. rustica, 370 CDS sequences against the N. sylvestris with the best hit and 212 and 185 CDS sequences against the N. paniculata and N. knightiana, respectively, supported the hybrid simulation of N. sylvestris and N. paniculata and (or) N. knightiana into N. rustica. In the tetraploid N. clevelandii, 191 CDS sequences against the N. attenuata with the best hit. And 235 and 123 CDS sequences against the N. undulata and N. paniculata, respectively. This result supported the hybrid simulation between N. attenuata and N. undulata or N. paniculata into N. clevelandii. In the tetraploid N. Velutina, N. cavicola, N. rosulata, N. benthamiana, and N. amplexicaulis, 233–258 CDS sequences against the N. sylvestris with the best hit, and another parent could be origin from section Petunioides, Noctiflorae or Alatae.
Based on the classified sequences of polyploid species and transcript fragments of plastid, we construct the phylogenetic relationship of nuclear and plastid, respectively, including all diploid and polyploid species. In the plastid phylogenetic relationship (Supplementary Figure S2), the N. tabacum was sister to the N. sylvestris; the N. rustica was sister to the N. knightiana; the N. clevelandii was sister to the N. obtusifolia; both N. stocktonii and section Suaveolentes was sister to the section Noctiflorae. In the nuclear phylogenetic relationship (Figure 3G), the N. tabacum was sister to the N. tomentosiformis and N. sylvestris, respectively; the N. clevelandii was sister to the N. undulata and N. attenuata, respectively; the N. stocktonii was sister to the N. obtusifolia and N. sylvestris, respectively; and the N. rustica was sister to N. undulata and the common ancestor of N. knightiana and N. paniculata, respectively. The results of phylogenetic relationships provided clues about the parent origin of polyploid species. We concluded the hybridization event by combining the phylogenetic relationships of nuclear and plastid (Figures 3G, S2). Namely, the N. tabacum formed from a single hybridization event between extant relatives of maternal progenitor N. sylvestris and paternal progenitor N. tomentosiformis; the N. stocktonii formed from a single hybridization event between maternal progenitor N. sylvestris and paternal progenitor N. obtusifolia; the N. rustica formed from a single hybridization event between paternal progenitor N. undulata and maternal progenitor common ancestor of N. knightiana and N. paniculata; the N. clevelandii formed from a single hybridization event between maternal progenitor N. attenuata and paternal progenitor N. undulata; the section Suaveolentes could be formed from a single hybridization event between maternal progenitor N. sylvestris and paternal progenitor of N. glauca.
3.6 Dating the time of hybridization event among Nicotiana species
To date the time of hybridization events that had led to the formation of polyploid species, we then used the classified CDS based on the sequence similarity search to obtain the time distance between the hybrid and the parental species (a hybridization date) (Figure 4; Supplementary Table S3). The divergence between the genus Nicotiana and Petunia was dated to c. 30.2 Ma (95% Highest Posterior Density (HPD) = 25.6–35.2 Ma). The diversification of Nicotiana was inferred to begin at c. 9.24 Ma (95% HPD = 7.84–10.77 Ma). The N. tabacum, as the only species in section Nicotiana, was the youngest allotetraploid section with age estimates of c. 0.42 Ma from its maternal progenitor N. sylvestris (95% HPD = 0.35-0.50 Ma) and c. 0.58 Ma from its paternal progenitor N. tomentosiformis (95% HPD = 0.49–0.68 Ma). Section Rusticae was a monotypic section containing only N. rustica and yielded age estimates of c.1.52 Ma from its maternal progenitor, the common ancestor of N. paniculata and N. knightiana (95% HPD = 1.28-1.76 Ma) and c. 1.30 Ma from its paternal progenitor N. undulata (95% HPD = 1.09-1.52 Ma). Section Polydicliae consists of two species, of which N. clevelandii yielded age estimates of c. 3.71 Ma from its maternal progenitor, the common ancestor of section Petunioides (95% HPD = 3.11-4.30 Ma) and c. 3.73 Ma from its paternal progenitor N. undulata (95% HPD = 3.12-4.30 Ma). Section Repandae contained four species, of which N. Stocktonii yielded age estimates of c. 5.02 Ma from its maternal progenitor N. sylvestris (95% HPD = 4.25-5.84 Ma) and c. 3.39 Ma from its paternal progenitor N. obtusifolia (95% HPD = 2.86-3.94 Ma). Section Suaveolentes was the oldest and most species-rich allotetraploid section with age estimates of c. 6.81 Ma from its maternal progenitor, the common ancestor of N. sylvestris (95% HPD = 5.76-7.92 Ma).
Figure 4 The phylogenetic tree showed the topology and divergence time for 26 Nicotiana species. Divergence times were indicated by light blue bars at the internodes; the range of these bars indicates 95% of the highest posterior density (HPD) interval of the divergence time. Numbers at the internodes indicate the mean divergence time. The geological timescale was illustrated at the bottom. ♀: maternal origin; ♂: paternal origin.
4 Discussion
4.1 Strongly supported diploid phylogeny and nuclear-plastid phylogenetic discordance
The Molecular Phylogeny of the Nicotiana genus has been researched for more than two decades (Aoki and Ito, 2000). The previous phylogenetic analyses using a combination of the internal transcribed spacer (ITS) and several plastid markers suggested that section Tomentosae was the base taxa and the section Petunioides was sister to the MCRA of section Noctiflorae, Alatae, and Sylvestres (Leitch et al., 2008). Recently, the nuclear phylogenetic tree based on GS and LFY genes was used to investigate the timing of diversification (Clarkson et al., 2017). However, these phylogeny relationships in the genus Nicotiana were inferred from several nuclear or plastid makers, which lack enough support and reliable results. Our recent analysis of whole plastid genomes provided a well-supported phylogenetic relationship of 11 sections in Nicotiana (Wang et al., 2022), which supported the section Tomentosae as the base clade of all others, and the section Petunioides was sister to the section Trigonophyllae.
Over the past ten years, high‐throughput transcriptome sequencing has provided an unprecedented volume of available genetic data. The transcriptome data have been widely used for reconstructing the phylogenetic relationship of plants (Leebens-Mack et al., 2019) regardless of the tissue origin of the transcriptomes (Cheon et al., 2020). We used gene and species tree approaches to construct a diploid phylogeny of Nicotiana based on RNAseq data that includes representatives from all 13 sections recognized in the Nicotiana genus. Our phylogenetic analyses of nuclear and plastid datasets produced mostly harmonious and well‐supported relationships (i.e., ≥92 BP) among major lineages in Nicotiana, including those not well resolved in previous studies. According to the nuclear phylogeny, the earliest divergence in Nicotiana involves two major clades (Figure 1A). The first major clade was formed by the most recent common ancestor (MRCA) of sections Trigonophyllae, Tomentosae, Undulatae, and Paniculatae. Section Trigonophyllae was sister to section Tomentosae, while section Undulatae was sister to section Paniculatae. The second major clade was formed by the MRCA of sections Petunioides, Noctiflorae, Sylvestres, and Alatae. All of these relationships are strongly supported, except for the placement of section Sylvestres (ASTRAL LPP = 0.94/ASTRAL BS = 90/concatenated BS > 99; 148/568 informative gene trees). According to the plastid phylogeny, section Tomentosae was sister to all remaining sections, with a grade formed by the MRCA of section Trigonophyllae, Petunioides, Undulatae, Paniculatae and the MRCA of section Alatae, Sylvestres, and Noctiflorae.
Our nuclear analyses revealed two major clades, which were not indicated by previous molecular or morphological analyses (McCarthy et al., 2016; Clarkson et al., 2017). The plastid phylogenetic tree in this study showed a consistent placement with our previous result based on the whole plastid genomes (Wang et al., 2022). Furthermore, the availability of genome‐scale data allowed us to examine the consistency of phylogenetic signals in the nuclear and plastid genomes for the first time, and several incongruent have been identified between nuclear and plastid phylogeny, especially for section Petunioides and Sylvestres (Figure 1). Incomplete lineage sorting and hybridization are the two main evolutionary processes that could lead to incongruent topologies between nuclear and organelle genomes (Willyard et al., 2009; Kao et al., 2022). It is often difficult to disentangle these processes. For section Sylvestres, incomplete lineage sorting should count for the incongruent topologies between gene trees. In addition, our analyses show an inconsistent placement of section Petunioides between nuclear and plastid phylogeny (Figures 1A, B), which is likely the result of an interspecific hybrid origin of this section between section Noctiflorae and Trigonophyllae in more early time. In contrast, the hybrid signal and nuclear introgression were lacking in section Petunioides based on the sequence similarity search and PhyloNet analysis, respectively (Figure 2B), which means that organelle capture can explain this observation. Using multiple gene trees enabled us to detect evidence for hybridization events between diploid species and resolved the phylogeny more robustly than in the previous studies.
4.2 Inferring polyploid parentage of Nicotiana species
As the genus Nicotiana contains several groups of tetraploids that formed at different times from different diploid progenitors (Kelly et al., 2013), it provides an ideal system for understanding polyploidization, a pervasive and powerful evolutionary force in plants. We identified the most likely diploid progenitors of five allopolyploid sections using the combined approaches of the phylogenetic network, sequence similarity search, and phylogenetic tree of subgenomes (Figure 3). Our results provide novel insights into parental species for the Nicotiana polyploids. Below we discussed the putative parental species of Nicotiana polyploids and proposed the results obtained in this study (Figure 5).
Figure 5 Cladogram summary of the polyploidization events and phylogenetic relationships in Nicotiana. Sections of allotetraploid origin were indicated by dashed black lines and solid black lines from their maternal lineages and paternal lineages, respectively. The section Petunioides involved organelle capture event was indicated by red dashed lines. The N. sylvestris involved incomplete lineage sorting (ILS) event was indicated by solid green lines. The time of hybridization events was noted under the allotetraploid sections. Sections were classified according to Knapp et al. (2004) and labeled to the left.
4.2.1 Section Nicotiana
N. tabacum (common tobacco), the only species in section Nicotiana, was most widely grown commercially for tobacco production. Its diploid ancestors (N. sylvestris as maternal progenitor and N. tomentosiformis as paternal progenitor) and the details of the hybridization have been well characterized (Sierro et al., 2013; Sierro et al., 2014). This study also provided three pieces of evidence to support that the diploid N. sylvestris and diploid N. tomentosiformis were the maternal and paternal progenitors of N. tabacum, respectively (Figure 5), which validated the reliability and accuracy of these strategies (Figures 3A, G).
4.2.2 Section Repandae
Section Repandae consists of four allopolyploid species (N. Nudicaulis, N. repanda, N. Stocktonii, and N. Nesophila; 2n = 4x = 48) (Knapp et al., 2004; Clarkson et al., 2005). It has been thought that section Repandae formed from a single hybridization event between extant relatives of maternal progenitor N. sylvestris and paternal progenitor N. obtusifolia. Subsequently, four allopolyploid species were produced following speciation (Chase et al., 2003; Clarkson et al., 2010; Kelly et al., 2013; Dodsworth et al., 2017). Based on our strategies, the diploid N. sylvestris and N. obtusifolia were recognized as the maternal and paternal progenitors of section Repandae, respectively (Figure 5), consistent with the previous results (Dodsworth et al., 2017).
4.2.3 Section rustica
Like common tobacco, N. rustica (Aztec tobacco) was the only species in section Rusticae and was an allotetraploid native to South America formed through a recent hybridization event (Knapp et al., 2004). Based on morphology, karyotype analyses, and breeding experiments, Goodspeed (1956) proposed that the diploid parents of N. rustica were ancestors of N. paniculata and N. undulata. Based on the comparative nuclear genome analysis, Sierro et al. (2018) found that the tetraploid species N. rustica inherited about 41% of its genome from its paternal progenitor, N. undulata, and 59% from its maternal progenitor, the common ancestor of N. paniculata and N. knightiana. Whereas, the pieces of evidence from comparative plastid genes and genome analysis revealed that the maternal parent of the tetraploid N. rustica was from the species of section Paniculatae, and the diploid N. knightiana was genetically closer than N. paniculata to N. rustica (Clarkson et al., 2004; Mehmood et al., 2020; Wang et al., 2022). The result from sequence similarity analysis support that both N. knightiana and N. paniculata might have donated the maternal genome of N. rustica (Figure 3B). Although the plastid genome of N. knightiana appears to be closer than that of N. paniculata to the N.rustica chloroplast genome, our analysis of PhyloNet and the phylogeny of subgenomes still suggested that a common ancestor of both N. knightiana and N. paniculata served as the maternal donor to N. rustica (Figure 5).
4.2.4 Section Polydicliae
Section Polydicliae consists of two species, N. quadrivalvis and N. clevelandii, the only allopolyploid section found in western North America (Goodspeed, 1956). In early research, plastid-based analyses indicated that a diploid species from section Trigonophyllae was the maternal genome donor of section Polydicliae (Clarkson et al., 2004). Based on the analysis of genome size, Leitch et al. (2008) found that the genome size of N. attenuata was most closely related to the paternal genome donor of section Polydicliae. Clarkson et al. (2010) suggested that section Polydiclieae was the product of a cross between the ancestors of section Trigonophyllae (maternal) and N. attenuata (paternal). Based on the strict consensus trees from the ADH and LFY/FLO datasets, Kelly et al. (2013) proposed that section Polydicliae formed from a single hybridization event between extant relatives of maternal progenitor N. obtusifolia and paternal progenitor N. attenuata. Subsequently, the analysis of floral evolution in section Polydicliae was performed based on this hybridization model (McCarthy et al., 2015; Bombarely et al., 2016). Based on the phylogenetic tree of five plastid loci and two nuclear genes, Clarkson et al. (2017) suggested a similar result. However, we obtained a different conclusion. In the plastid phylogenetic relationship (Supplementary Figure S2), the N. clevelandii was sister to the N. obtusifolia, the plastid of which was transferred to section Petunioides. Thus, we now believe that section Polydiclieae is the product of a cross between the ancestors of N. attenuata (maternal) and N. undulata (paternal) based on our strategies (Figure 5).
4.2.5 Section Suaveolentes
Section Suaveolentes was an almost all-Australian clade (the exception being N. africana of Namibia) of allopolyploid species, including the vital plant model N. benthamiana (Schiavinato et al., 2020). Karyotypic variation within this section was very enrichment from n = 15 to n = 32 chromosomes (Marks et al., 2011). It seems likely that section Suaveolentes has explosive radiation of taxa occurred, primarily accompanied by diploid reductions probably due to fusions of chromosomes (Clarkson et al., 2004). As the oldest Nicotiana polyploids (Clarkson et al., 2017), the diploid progenitors of section Suaveolentes were poorly understood (Marks et al., 2011). Goodspeed (1956) proposed that several diploid sections of Nicotiana, namely sections Alatae, Sylvestres, Noctiflorae, and Petunioides, were involved in the formation of the allotetraploid section Suaveolentes. Kelly et al. (2013) reconstructed the evolutionary origin of sect Suaveolentes using four regions from the nuclear and plastid genome. They identified a likely scenario where a member of the N. Sylvestres acted as the paternal progenitor, and a member of either section Petunioides or Noctiflorae was the maternal progenitor. Recently, Schiavinato et al. (2020) showed that the maternal progenitor of N. benthamiana was a member of section Noctiflorae and confirmed a member of section Sylvestres as a paternal subgenome donor. Our analysis based on the PhyloNet approaches supported a scenario where the N. sylvestris acted as the paternal progenitor, and the N. glauca of section Noctiflorae acted as the maternal progenitor in the formation of section Suaveolentes (Figure 3E), in line with previous findings (Schiavinato et al., 2020). But the analysis from the sequence similarity search and the phylogeny of the subgenome only support that N. sylvestris could have been involved in its formation (Figure 3G). Lastly, we summarized that the N. sylvestris (paternal progenitor) and the N. glauca (maternal progenitor) were involved in the formation of section Suaveolentes (Figure 5).
5 Conclusion
In conclusion, this study sheds light on the genetic diversity, phylogenetic relationships, and evolutionary history of Nicotiana species. The findings provide valuable insights into the classification and phylogenetic relationships within the genus. The identification of parental origins and the estimation of divergence times of polyploid species contribute to our understanding of speciation and hybridization events. Furthermore, the application of high-throughput RNA-seq technology in this study demonstrates its efficacy in phylogenetic studies and paves the way for future molecular systematic and population genetic investigations. The comprehensive dataset and analytical approaches used in this study serve as a valuable resource for further research in Nicotiana and related species.
Data availability statement
The original contributions presented in the study are included in the article/Supplementary Material. Further inquiries can be directed to the corresponding author.
Author contributions
SW and JG conceived and designed the study. SW and ZL conducted the bioinformatics analysis. ZL, WP, and KC assisted in data collection. SW, JG, and CF wrote and revised the manuscript. All authors read and approved the final manuscript.
Funding
This work was supported by the Special Project for Lushan Plants [2021ZWZX18] and the Key Funding of CNTC [No. 110201801030 (JY-07)] and CNTC [No. 110202101003 (JY-03)].
Conflict of interest
SW, JG, KC, and WP were employed by China Tobacco Hunan Industrial Corporation.
The remaining authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
This study received funding from CNTC No. 110201801030 JY-07 and CNTC No. 110202101003 JY-03. The funder was not involved in the study design, collection, analysis, interpretation of data, the writing of this article, or the decision to submit it for publication. All authors declare no other competing interests.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fpls.2023.1205683/full#supplementary-material
References
Angiuoli, S. V., Salzberg, S. L. (2011). Mugsy: fast multiple alignment of closely related whole genomes. Bioinformatics 27, 334–342. doi: 10.1093/bioinformatics/btq665
Aoki, S., Ito, M. (2000). Molecular phylogeny of Nicotiana (Solanaceae) based on the nucleotide sequence of the matK gene. Plant Biol. 2, 316–324. doi: 10.1055/s-2000-3710
Bally, J., Marks, C. E., Jung, H., Jia, F., Roden, S., Cooper, T., et al. (2021). Nicotiana paulineana, a new Australian species in Nicotiana section Suaveolentes. Aust. Syst. Bot. 34, 477–484. doi: 10.1071/SB20025
Bombarely, A., Moser, M., Amrad, A., Bao, M., Bapaume, L., Barry, C. S., et al. (2016). Insight into the evolution of the Solanaceae from the parental genomes of Petunia hybrida. Nat. Plants 2, 16074. doi: 10.1038/nplants.2016.74
Borowiec, M. L. (2016). AMAS: a fast tool for alignment manipulation and computing of summary statistics. PeerJ 4, e1660. doi: 10.7717/peerj.1660
Bouckaert, R., Vaughan, T. G., Barido-Sottani, J., Duchêne, S., Fourment, M., Gavryushkina, A., et al. (2019). BEAST 2.5: An advanced software platform for Bayesian evolutionary analysis. PloS Comput. Biol. 15, e1006650. doi: 10.1371/journal.pcbi.1006650
Brown, J. W., Walker, J. F., Smith, S. A. (2017). Phyx: phylogenetic tools for unix. Bioinformatics 33, 1886–1888. doi: 10.1093/bioinformatics/btx063
Camacho, C., Coulouris, G., Avagyan, V., Ma, N., Papadopoulos, J., Bealer, K., et al. (2009). BLAST+: architecture and applications. BMC Bioinf. 10, 421. doi: 10.1186/1471-2105-10-421
Chase, M. W., Knapp, S., Cox, A. V., Clarkson, J. J., Butsko, Y., Joseph, J., et al. (2003). Molecular systematics, GISH and the origin of hybrid taxa in Nicotiana (Solanaceae). Ann. Bot. 92, 107–127. doi: 10.1093/aob/mcg087
Chen, C.-L., Zhang, L., Li, J.-L., Mao, X.-X., Zhang, L.-S., Hu, Q.-J., et al. (2021). Phylotranscriptomics reveals extensive gene duplication in the subtribe Gentianinae (Gentianaceae). J. Syst. Evol. 59, 1198–1208. doi: 10.1111/jse.12651
Chen, S., Zhou, Y., Chen, Y., Bioinformatics, J. G. J. (2018). fastp: an ultra-fast all-in-one FASTQ preprocessor. Bioinformatics 34, i884–i890. doi: 10.1093/bioinformatics/bty560
Cheon, S., Lee, S.-G., Hong, H.-H., Lee, H.-G., Kim, K. Y., Park, C. (2021). A guide to phylotranscriptomic analysis for phycologists. Algae 36, 333–340. doi: 10.4490/algae.2021.36.12.7
Cheon, S., Zhang, J., Park, C. (2020). Is phylotranscriptomics as reliable as phylogenomics? Mol. Biol. Evol. 37, 3672–3683. doi: 10.1093/molbev/msaa181
Clarkson, J. J., Dodsworth, S., Chase, M. W. (2017). Time-calibrated phylogenetic trees establish a lag between polyploidisation and diversification in Nicotiana (Solanaceae). Plant Syst. Evol. 303, 1001–1012. doi: 10.1007/s00606-017-1416-9
Clarkson, J. J., Kelly, L. J., Leitch, A. R., Knapp, S., Chase, M. W. (2010). Nuclear glutamine synthetase evolution in Nicotiana: Phylogenetics and the origins of allotetraploid and homoploid (diploid) hybrids. Mol. Phylogenet. Evol. 55, 99–112. doi: 10.1016/j.ympev.2009.10.003
Clarkson, J., Knapp, S., Garcia, V. F., Olmstead, R. G., Leitch, A., Chase, M. (2004). Phylogenetic relationships in Nicotiana (Solanaceae) inferred from multiple plastid DNA regions. Mol. Phylogenet. Evol. 33, 75–90. doi: 10.1016/j.ympev.2004.05.002
Clarkson, J. J., Lim, K. Y., Kovarik, A., Chase, M. W., Knapp, S., Leitch, A. R. (2005). Long-term genome diploidization in allopolyploid Nicotiana section Repandae (Solanaceae). New Phytol. 168, 241–252. doi: 10.1111/j.1469-8137.2005.01480.x
Danecek, P., Bonfield, J. K., Liddle, J., Marshall, J., Ohan, V., Pollard, M. O., et al. (2021). Twelve years of SAMtools and BCFtools. GigaScience 10, giab008. doi: 10.1093/gigascience/giab008
Darriba, D., Posada, D., Kozlov, A. M., Stamatakis, A., Morel, B., Flouri, T. (2019). ModelTest-NG: A new and scalable tool for the selection of DNA and protein evolutionary models. Mol. Biol. Evol. 37, 291–294. doi: 10.1093/molbev/msz189
Dobin, A., Davis, C. A., Schlesinger, F., Drenkow, J., Zaleski, C., Jha, S., et al. (2012). STAR: ultrafast universal RNA-seq aligner. Bioinformatics 29, 15–21. doi: 10.1093/bioinformatics/bts635
Dodsworth, S., Jang, T.-S., Struebig, M., Chase, M. W., Weiss-Schneeweiss, H., Leitch, A. R. (2017). Genome-wide repeat dynamics reflect phylogenetic distance in closely related allotetraploid Nicotiana (Solanaceae). Plant Syst. Evol. 303, 1013–1020. doi: 10.1007/s00606-016-1356-9
Edwards, K. D., Fernandez-Pozo, N., Drake-Stowe, K., Humphry, M., Evans, A. D., Bombarely, A., et al. (2017). A reference genome for Nicotiana tabacum enables map-based cloning of homeologous loci implicated in nitrogen utilization efficiency. BMC Genomics 18, 448. doi: 10.1186/s12864-017-3791-6
Emms, D. M., Kelly, S. (2015). OrthoFinder: solving fundamental biases in whole genome comparisons dramatically improves orthogroup inference accuracy. Genome Biol. 16, 157. doi: 10.1186/s13059-015-0721-2
Fukushima, K., Pollock, D. D. (2020). Amalgamated cross-species transcriptomes reveal organ-specific propensity in gene expression evolution. Nat. Commun. 11, 4459. doi: 10.1038/s41467-020-18090-8
Goodspeed, T. H. (1956). The genus Nicotiana. J. Am. Pharm. Assoc. 45, 193–193. doi: 10.1002/jps.3030450326
Grabherr, M. G., Haas, B. J., Yassour, M., Levin, J. Z., Thompson, D. A., Amit, I., et al. (2011). Full-length transcriptome assembly from RNA-Seq data without a reference genome. Nat. Biotechnol. 29, 644–U130. doi: 10.1038/nbt.1883
Huang, Y., Niu, B., Gao, Y., Fu, L., Li, W. (2010). CD-HIT Suite: a web server for clustering and comparing biological sequences. Bioinformatics 26, 680–682. doi: 10.1093/bioinformatics/btq003
Huson, D. H., Bryant, D. (2005). Application of phylogenetic networks in evolutionary studies. Mol. Biol. Evol. 23, 254–267. doi: 10.1093/molbev/msj030
Huson, D. H., Scornavacca, C. (2012). Dendroscope 3: An interactive tool for rooted phylogenetic trees and networks. Syst. Biol. 61, 1061–1067. doi: 10.1093/sysbio/sys062
Junier, T., Zdobnov, E. M. (2010). The Newick utilities: high-throughput phylogenetic tree processing in the Unix shell. Bioinformatics 26, 1669–1670. doi: 10.1093/bioinformatics/btq243
Kao, T.-T., Wang, T.-H., Ku, C. (2022). Rampant nuclear–mitochondrial–plastid phylogenomic discordance in globally distributed calcifying microalgae. New Phytol. 235, 1394–1408. doi: 10.1111/nph.18219
Katoh, K., Standley, D. M. (2013). MAFFT multiple sequence alignment software version 7: Improvements in performance and usability. Mol. Biol. Evol. 30, 772–780. doi: 10.1093/molbev/mst010
Kelly, L. J., Leitch, A. R., Clarkson, J. J., Knapp, S., Chase, M. W. (2013). Reconstructing the complex evolutionary origin of wild allopolyploid tobaccos (Nicotiana section suaveolentes). Evolution 67, 80–94. doi: 10.1111/j.1558-5646.2012.01748.x
Khan, M. Q., Narayan, R. K. J. (2007). Phylogenetic diversity and relationships among species of genus Nicotiana using RAPDs analysis. Afr. J. Biotechnol. 6, 148–162. doi: 10.5897/AJB06.442
Kitamura, S., Inoue, M., Shikazono, N., Tanaka, A. (2001). Relationships among Nicotiana species revealed by the 5S rDNA spacer sequence and fluorescence in situ hybridization. Theor. Appl. Genet. 103, 678–686. doi: 10.1007/s001220100643
Knapp, S., Chase, M. W., Clarkson, J. J. (2004). Nomenclatural changes and a new sectional classification in Nicotiana (Solanaceae). Taxon 53, 73–82. doi: 10.2307/4135490
Kozlov, A. M., Darriba, D., Flouri, T., Morel, B., Stamatakis, A. (2019). RAxML-NG: a fast, scalable and user-friendly tool for maximum likelihood phylogenetic inference. Bioinformatics 35, 4453–4455. doi: 10.1093/bioinformatics/btz305
Leebens-Mack, J. H., Barker, M. S., Carpenter, E. J., Deyholos, M. K., Gitzendanner, M. A., Graham, S. W., et al. (2019). One thousand plant transcriptomes and the phylogenomics of green plants. Nature 574, 679–685. doi: 10.1038/s41586-019-1693-2
Leitch, I. J., Hanson, L., Lim, K. Y., Kovarik, A., Chase, M. W., Clarkson, J. J., et al. (2008). The ups and downs of genome size evolution in polyploid species of Nicotiana (Solanaceae). Ann. Bot. 101, 805–814. doi: 10.1093/aob/mcm326
Lemmon, E. M., Lemmon, A. R. (2013). High-throughput genomic data in systematics and phylogenetics. Annu. Rev. Ecol. Evol. Syst. 44, 99–121. doi: 10.1146/annurev-ecolsys-110512-135822
Lim, K. Y., Matyasek, R., Kovarik, A., Fulnecek, J., Leitch, A. R. (2005). Molecular cytogenetics and tandem repeat sequence evolution in the allopolyploid Nicotiana rustica compared with diploid progenitors N. paniculata and N. undulata. Cytogenetic Genome Res. 109, 298–309. doi: 10.1159/000082413
Liu, Y.-Y., Jin, W.-T., Wei, X.-X., Wang, X.-Q. (2022). Phylotranscriptomics reveals the evolutionary history of subtropical East Asian white pines: further insights into gymnosperm diversification. Mol. Phylogenet. Evol. 168, 107403. doi: 10.1016/j.ympev.2022.107403
Manni, M., Berkeley, M. R., Seppey, M., Simão, F. A., Zdobnov, E. M. (2021). BUSCO update: Novel and streamlined workflows along with broader and deeper phylogenetic coverage for scoring of eukaryotic, prokaryotic, and viral genomes. Mol. Biol. Evol. 38, 4647–4654. doi: 10.1093/molbev/msab199
Marks, C. E., Ladiges, P. Y., Newbigin, E. (2011). Karyotypic variation in Nicotiana section Suaveolentes. Genet. Resour. Crop Evol. 58, 797–803. doi: 10.1007/s10722-011-9724-3
McCarthy, E. W., Arnold, S. E. J., Chittka, L., Le Comber, S. C., Verity, R., Dodsworth, S., et al. (2015). The effect of polyploidy and hybridization on the evolution of floral colour in Nicotiana (Solanaceae). Ann. Bot. 115, 1117–1131. doi: 10.1093/aob/mcv048
McCarthy, E. W., Chase, M. W., Knapp, S., Litt, A., Leitch, A. R., Le Comber, S. C. (2016). Transgressive phenotypes and generalist pollination in the floral evolution of Nicotiana polyploids. Nat. Plants 2, 16119. doi: 10.1038/nplants.2016.119
Mehmood, F., Abdullah, Ubaid, Z., Shahzadi, I., Ahmed, I., Waheed, M. T., et al. (2020). Plastid genomics of Nicotiana (Solanaceae): insights into molecular evolution, positive selection and the origin of the maternal genome of Aztec tobacco (Nicotiana rustica). PeerJ 8, e9552. doi: 10.1101/2020.01.13.905158
Mistry, J., Chuguransky, S., Williams, L., Qureshi, M., Salazar, G. A., Sonnhammer, E. L. L., et al. (2020). Pfam: The protein families database in 2021. Nucleic Acids Res. 49, D412–D419. doi: 10.1093/nar/gkaa913
Morel, B., Kozlov, A. M., Stamatakis, A. (2018). ParGenes: a tool for massively parallel model selection and phylogenetic tree inference on thousands of genes. Bioinformatics 35, 1771–1773. doi: 10.1093/bioinformatics/bty839
Pease, J. B., Brown, J. W., Walker, J. F., Hinchliff, C. E., Smith, S. A. (2018). Quartet Sampling distinguishes lack of support from conflicting support in the green plant tree of life. Am. J. Bot. 105, 385–403. doi: 10.1002/ajb2.1016
Peng, J.-C., Ma, X.-G., Wang, Y.-H., Sun, H. (2020). New insights into the evolutionary history of Megacodon: Evidence from a newly discovered species. Plant Diversity 42, 198–208. doi: 10.1016/j.pld.2020.05.003
Qi, X., Kuo, L.-Y., Guo, C., Li, H., Li, Z., Qi, J., et al. (2018). A well-resolved fern nuclear phylogeny reveals the evolution history of numerous transcription factor families. Mol. Phylogenet. Evol. 127, 961–977. doi: 10.1016/j.ympev.2018.06.043
Renny-Byfield, S., Kovarik, A., Kelly, L. J., Macas, J., Novak, P., Chase, M. W., et al. (2013). Diploidization and genome size change in allopolyploids is associated with differential dynamics of low- and high-copy sequences. Plant J. 74, 829–839. doi: 10.1111/tpj.12168
Ronquist, F., Teslenko, M., van der Mark, P., Ayres, D. L., Darling, A., Höhna, S., et al. (2012). MrBayes 3.2: efficient Bayesian phylogenetic inference and model choice across a large model space. Syst. Biol. 61, 539–542. doi: 10.1093/sysbio/sys029
Särkinen, T., Bohs, L., Olmstead, R. G., Knapp, S. (2013). A phylogenetic framework for evolutionary study of the nightshades (Solanaceae): a dated 1000-tip tree. BMC Evol. Biol. 13, 214. doi: 10.1186/1471-2148-13-214
Schiavinato, M., Marcet-Houben, M., Dohm, J. C., Gabaldón, T., Himmelbauer, H. (2020). Parental origin of the allotetraploid tobacco Nicotiana benthamiana. Plant J. 102, 541–554. doi: 10.1111/tpj.14648
Sierro, N., Battey, J. N. D., Bovet, L., Liedschulte, V., Ouadi, S., Thomas, J., et al. (2018). The impact of genome evolution on the allotetraploid Nicotiana rustica – an intriguing story of enhanced alkaloid production. BMC Genomics 19, 855. doi: 10.1186/s12864-018-5241-5
Sierro, N., Battey, J. N. D., Ouadi, S., Bakaher, N., Bovet, L., Willig, A., et al. (2014). The tobacco genome sequence and its comparison with those of tomato and potato. Nat. Commun. 5, 3833. doi: 10.1038/ncomms4833
Sierro, N., Battey, J. N., Ouadi, S., Bovet, L., Goepfert, S., Bakaher, N., et al. (2013). Reference genomes and transcriptomes of Nicotiana sylvestris and Nicotiana tomentosiformis. Genome Biol. 14, R60. doi: 10.1186/gb-2013-14-6-r60
Smith, S. A., Moore, M. J., Brown, J. W., Yang, Y. (2015). Analysis of phylogenomic datasets reveals conflict, concordance, and gene duplications with examples from animals and plants. BMC Evol. Biol. 15, 150. doi: 10.1186/s12862-015-0423-0
Suzek, B. E., Huang, H., McGarvey, P., Mazumder, R., Wu, C. H. (2007). UniRef: comprehensive and non-redundant UniProt reference clusters. Bioinformatics 23, 1282–1288. doi: 10.1093/bioinformatics/btm098
Svab, Z., Maliga, P. (2007). Exceptional transmission of plastids and mitochondria from the transplastomic pollen parent and its impact on transgene containment. Proc. Natl. Acad. Sci. 104, 7003–7008. doi: 10.1073/pnas.0700063104
Talavera, G., Castresana, J. (2007). Improvement of phylogenies after removing divergent and ambiguously aligned blocks from protein sequence alignments. System. Biol. 56, 564–577. doi: 10.1080/10635150701472164
Tria, F. D. K., Landan, G., Dagan, T. (2017). Phylogenetic rooting using minimal ancestor deviation. Nat. Ecol. Evol. 1, 193. doi: 10.1038/s41559-017-0193
Wang, X., Bennetzen, J. (2015). Current status and prospects for the study of Nicotiana genomics, genetics, and nicotine biosynthesis genes. Mol. Genet. Genomics 290, 1–11. doi: 10.1007/s00438-015-0989-7
Wang, S., Gao, J., Chao, H., Li, Z., Pu, W., Wang, Y., et al. (2022). Comparative chloroplast genomes of Nicotiana species (Solanaceae): Insights into the genetic variation, phylogenetic relationship, and polyploid speciation. Front. Plant Sci. 13. doi: 10.3389/fpls.2022.899252
Wen, D., Yu, Y., Zhu, J., Nakhleh, L. (2018). Inferring phylogenetic networks using PhyloNet. Syst. Biol. 67, 735–740. doi: 10.1093/sysbio/syy015
Willyard, A., Cronn, R., Liston, A. (2009). Reticulate evolution and incomplete lineage sorting among the ponderosa pines. Mol. Phylogenet. Evol. 52, 498–511. doi: 10.1016/j.ympev.2009.02.011
Yan, Z., Cao, Z., Liu, Y., Ogilvie, H. A., Nakhleh, L. (2021). Maximum parsimony inference of phylogenetic networks in the presence of polyploid complexes. Syst. Biol. 71, 706–720. doi: 10.1093/sysbio/syab081
Yang, Y., Moore, M. J., Brockington, S. F., Mikenas, J., Olivieri, J., Walker, J. F., et al. (2018). Improved transcriptome sampling pinpoints 26 ancient and more recent polyploidy events in Caryophyllales, including two allopolyploidy events. New Phytol. 217, 855–870. doi: 10.1111/nph.14812
Yang, L., Su, D., Chang, X., Foster, C. S. P., Sun, L., Huang, C.-H., et al. (2020). Phylogenomic insights into deep phylogeny of Angiosperms based on broad nuclear gene sampling. Plant Commun. 1, 100027. doi: 10.1016/j.xplc.2020.100027
Yu, Y., Nakhleh, L. (2015). A maximum pseudo-likelihood approach for phylogenetic networks. BMC Genomics 16, S10. doi: 10.1186/1471-2164-16-S10-S10
Yukawa, M., Tsudzuki, T., Sugiura, M. (2006). The chloroplast genome of Nicotiana sylvestris and Nicotiana tomentosiformis: complete sequencing confirms that the Nicotiana sylvestris progenitor is the maternal genome donor of Nicotiana tabacum. Mol. Genet. Genomics 275, 367–373. doi: 10.1007/s00438-005-0092-6
Keywords: Nicotiana L., transcriptomic analysis, polyploid formation, evolutionary relationship, divergence time
Citation: Wang S, Gao J, Li Z, Chen K, Pu W and Feng C (2023) Phylotranscriptomics supports numerous polyploidization events and phylogenetic relationships in Nicotiana. Front. Plant Sci. 14:1205683. doi: 10.3389/fpls.2023.1205683
Received: 14 April 2023; Accepted: 04 July 2023;
Published: 28 July 2023.
Edited by:
Muhammad Aamir Manzoor, Shanghai Jiao Tong University, ChinaReviewed by:
Jingjian Li, South China Agricultural University, ChinaMuhammad Anwar, Hainan University, China
Copyright © 2023 Wang, Gao, Li, Chen, Pu and Feng. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Chen Feng, ZmVuZ2NAbHNiZy5jbg==
†These authors have contributed equally to this work