About the Authors:
Mehrshad Zeinalabedini
Roles Conceptualization, Data curation, Formal analysis, Methodology, Project administration, Software, Supervision, Visualization, Writing – original draft
* E-mail: [email protected] (MZ); [email protected] (NAKS)
Affiliation: Department of Systems and Synthetic biology, Agricultural Biotechnology Research Institute of Iran, Agricultural Research, Education and Extension Organization (AREEO), Karaj, Iran
ORCID logo https://orcid.org/0000-0002-3436-4334
Nayer Azam Khoshkholgh Sima
Roles Project administration, Resources, Supervision
* E-mail: [email protected] (MZ); [email protected] (NAKS)
Affiliation: Department of Molecular Plant Physiology, Agricultural Biotechnology Research Institute of Iran, Agricultural Research, Education and Extension Organization (AREEO), Karaj, Iran
Mohammad Reza Ghaffari
Roles Data curation, Formal analysis, Resources, Writing – review & editing
Affiliation: Department of Systems and Synthetic biology, Agricultural Biotechnology Research Institute of Iran, Agricultural Research, Education and Extension Organization (AREEO), Karaj, Iran
ORCID logo https://orcid.org/0000-0002-7139-8613
Ali Ebadi
Roles Investigation, Validation, Writing – original draft
Affiliation: Department of Molecular Plant Physiology, Agricultural Biotechnology Research Institute of Iran, Agricultural Research, Education and Extension Organization (AREEO), Karaj, Iran
ORCID logo https://orcid.org/0000-0002-5304-0706
Maryam Farsi
Roles Data curation, Investigation
Affiliation: Department of Systems and Synthetic biology, Agricultural Biotechnology Research Institute of Iran, Agricultural Research, Education and Extension Organization (AREEO), Karaj, Iran
Introduction
The world’s population is 7 billion in 2018 and is expected to rise by 35% and reach 9 billion in 2050. Food and Agriculture Organization of United Nations (FAO) estimates crop production must increase by at least 60% in the future to provide food security for the world’s growing population [1,2]. However, the use of genotypes with high-yield and more resistant to biotic and abiotic stress has gradually replaced indigenous varieties and local genotypes, endangering access and exploitation of resources and food production for the future [3]. Thus, low genetic diversity may reduce the opportunity to identify and use new sources responding to future challenges, including new pests and pathogens, as well as climate changes [4]. Today, soil salinity has increased due to unconventional cultivation and irrigation. Moreover, the level of arable land increased dramatically from 8 million to more than 220 million hectares, of which about 45 million hectares are salt-affected to varying degrees [5,6].
Halophytes, which include more than 600 taxa of the various genus and species [7], complete their life cycle under high salt concentration (at least 200 mM NaCl) [8,9]. Salicornia is a succulent halophyte, which grows naturally on mangrove swamps as well as seashores. It has also been evaluated as a vegetable, forage, and oilseed crop in the agronomic field trials [10]. There is a morphological variation across the genus distribution range of Salicornia sp. Some economically feasible applications have been suggested for Salicornia species. They also have suitable biomass for human and domestics consumption [11,12]. Furthermore, Salicornia seed contains 26%–33% oil, similar to safflower oil, along with 30%–33% protein [13].
Members of the Salicornioideae are promising candidates for saline agriculture [14]. For instance, Salicornia bigelovii Torr., a true halophyte, is a potential new crop due to the oil content, fresh vegetable, forage, and its seeds in coastal and saline lands. In recent years, Salicornia has been cultivated for various purposes around the world, from the production of Salicornia oil by a Saudi company to its cultivation by the UAE Institute of Science and Technology to produce aircraft biofuels [14]. The highest saline soils in Asia are China, India, Pakistan, and Iran [15]. Iran is located in arid and semi-arid geographical areas rich in the distribution of naturally occurring materials such as halite (NaCl) and gypsum (CaSO4 covering 34% of the regions of the country, especially in the central, south, and plain of Khuzestan in varying degrees of salinity [7]. Furthermore, about 20%-50% of the arable land is also affected by salinity [8,9].
Salicornia sp. is widespread across Iran, mainly in central, western, and northern parts. Thus, native Salicornia species selection might be a viable agricultural development strategy against Iran’s climate change. Pilot cultivation of Salicornia in different regions of Iran and study of its forage quality in collaboration with Animal Science Research Institute of IRAN (ASRI) showed that this plant could supply up to 21% of the volume of forage for feeding of light livestock such as sheep (Unpublished).
The number of Salicornia species is between 25 to 30 species worldwide [16]. The taxonomy classification of Salicornia is very complicated as there is no general published morphological descriptor for all accepted species. Previously, some morphological characteristics like growth form, inflorescence, branching of main and lateral stems, lateral flowers and their conditions relative to the main flowers, fruit formation, flowering, and fruiting characteristics as well as dry biomass were used to identify some species of Salicornia with high seed and biomass yield [17–19]. However, a low number of morphological characteristics, phenotypic flexibility, breeding, and hybridization systems are amongst the factors that make Salicornia challenging to distinguish precisely in an area [20]. Thus, the commercial importance and lack of description make Salicornia a priority for the circumscription of species, subtypes, ecotypes, and natural hybridization to develop it as a new crop [21]. Currently, few studies were conducted on taxonomic and phylogenetic analysis of the Salicornia genus. Teege et al. [21] studied the genetic relationships of intra and inter-variation of the two taxa of Salicornia prompumbens and Salicornia stricta species using AFLP markers [17]. Further, the phylogenetic relationships between eight taxa belonging to the Salicornia genus were examined using the ITS markers [21]. Likewise, the phylogenetic and geographic relationship in different germplasms of Salicornia was evaluated using six EST markers [21].
Moreover, the genetic diversity of some populations of Salicornia was investigated using RAPD markers [22]. All obtained results showed that molecular markers have some limitations/disadvantages inaccurate taxonomic classifications. DNA sequencing technology has recently facilitated identifying identical genotypes with different morphotypes or different genotypes with the same morphotypes [23]. DNA barcoding, which is a new and efficient tool based on the conserved DNA regions of organisms, has provided an efficient and precise method for species-level identifications [24]. DNA barcodes were developed by Consortium for the Barcode of Life (CBOL) as a global standard for the identification of biological species. To date, no comprehensive study has been performed on the utility of universal DNA barcodes for all proposed gene regions (either alone or in combination) in Salicornia. The overall aim of our study was to employ universality of three coding plastid regions (rbcL, matK and ycf), one non-coding plastid intergenic spacer regions (trnH-psbA), and the internal transcribed spacer of nuclear-encoded ribosomal DNA (ITS2) to identify the Salicornia species using molecular data and ecological niche. An integrated model was then developed to assess Salicornia distribution and best habitat based on climate factors and estimate critical parameters in land management and biodiversity conservation. We assume that our results will open up a new window for the future development of Salicornia farms.
Materials and methods
Sample collection
The previously collected records by the Agricultural Biotechnology Research Institute of Iran (ABRII) were used to identify distribution localities for sampling [25]. Field surveys were then done to collect samples from populations throughout entire natural distribution of Salicornia sp in Iran. Twenty-six different geographical locations were included in this research (Fig 1). A total of 300 samples were collected from different geographical areas ranging from 2 to 5 samples per location (Table 1). A minimum of 91 required samples was selected for sequencing. At least 3–5 samples were successfully sequenced and analyzed.
[Figure omitted. See PDF.]
Fig 1. Distribution map of Salicornia germplasms in Iran.
The gateway to astronaut photography of earth was used to create the map (https://eol.jsc.nasa.gov/).
https://doi.org/10.1371/journal.pone.0241162.g001
[Figure omitted. See PDF.]
Table 1. Geographic information of Salicornia populations that were collected in this study.
https://doi.org/10.1371/journal.pone.0241162.t001
Molecular analysis
DNA extraction and barcoding.
DNA extraction was carried out using the Core Bio DNA extraction kit (Core Bio, South Korea). The quality and quantity of DNA were determined using 0.8% agarose gel and NanoDrop ND-1000 spectrophotometer (Nanodrop Technologies, Wilmington, DE, USA), respectively. DNA barcoding was carried out using core-barcode consisting of two plastid coding regions, (rbcL+matK) and supplemented with one non-coding intergenic spacer (trnH-psbA) from the chloroplast genome, one internal transcribed spacer (ITS2) from the nuclear genome, and one plastid coding region (ycf) [24,26] (see S1 Table). For all markers, 50 μl reactions containing 1 U Taq DNA polymerase (Invitrogen or Promega), 2.5 μl Taq DNA polymerase buffer (10X), 1.5 μl MgCl2 (25 mM), 0.5 μl dNTPs (10 mM), 0.35 μl primers (20 pmol), 1.0 μl template DNA and 18.85 μl dH2O were carried out in a Veriti ABI thermal cycler (ABI Inc, USA). Then, the PCR products were purified using PureLINK PCR product purification kit (Thermo Fisher Scientific). The purified PCR product was then directly sequenced in both directions to minimize PCR artifacts, ambiguities, and base-calling errors on an automated ABI Prism 3730 XL Genetic Analyzer machine using ABI BigDye v3.1 Terminator Sequencing chemistry (Macrogen Inc, South Korea). In this study, the GC content threshold above 65% was considered high GC content. GC content is a good predictor of PCR success and can significantly impact the PCR success rate and sequencing error rate. It is also an indicator of binding ability and specificity (primer binding and specificity) to optimize the DNA barcoding process. Therefore, the sequenced samples with above 65% GC content were cloned and then sequenced.
The sequencing results were then sorted and trimmed using FinchTV software (https://digitalworldbiology.com/FinchTV) and analyzed by protein BLAST (BLASTP), nucleotide (BLAST) BLASTN [27], and LALIGN software (https://embnet.vital-it.ch/software/LALIGN_form.html). Finally, DNA barcode data compared to GenBank or otherwise publicly available in the BOLD database based on different integration strategies (coding, non-coding, plastid, and total combination) using a similarity-based method. The sequence similarity-based method was briefly utilized for species identification, using simple and optimized BLAST methods for different marker integration strategies (single, coding, non-coding, plastid, and total combination). First, sequences were queried using megablast using Geneious software (Biomatters) and NCBI (national center for biotechnology information) BLASTn against the nucleotide database. A simple method was then considered the top 10 hits, and the top 100 records were used to put extra weight on the identity value using "max score*(query cover/identity)" formula using the optimized BLAST method. Finally, taxonomic identifications were allocated based on the combination of the identity score (High identity: X≥95%; Medium identity: 90% ≤ X ≤ 95%; Low identity: X ≤ 90%) and the number of species within 1% deviation of the calculated similarity score (S1 Fig). Plant Sequences database (BOLD systems) and ITS2 database (V5) were also used as a validation tool for the core barcodes (matK, rbcL) and ITS2 [28,29]. Phylogenetic tree was obtained by MEGA X version 10.2.2 (https://www.megasoftware.net/) by analyzing all the sequence data.
Analysis of genetic diversity
Analysis of molecular variance (AMOVA) was carried out using Arlequin version 3.5 [30]. The significance level of Fst statistics was also performed using the software by a nonparametric permutation procedure with 1023 randomizations. Statistical calculations and graphics for Fst were conducted using Arlequin and R version 3.6.1 (cran.r-project.org) [31]. Further, the online version of Automatic Barcode Gap Discovery (ABGD) (https://bioinfo.mnhn.fr/abi/public/abgd/abgdweb.html) [32] was used to generate the genetic distance histograms and ranked distance based on K2P distance. The ABGD detects a gap. Then, when a gap exists in the distribution of divergence corresponded to differences between intra-specific and inter-specific distances, the method will work well for species delimitation. The following setting was used to perform the analysis: relative gap width X = 1.5, K2P distance, and intraspecific divergence (P) values range from 0.001 to 0.100. Further parameter values were employed based on the software defaults. The analyses were performed at the web server. In the ABGD system, specimens identified as one group could be considered as one species [32].
Species distribution modeling
The modeling program MAXENT v.3.3.3 was used in the default settings with a random test percentage of 25% of the input localities set aside for modeling Iranian Salicornia sp. potential distribution. The maximum entropy algorithm calculates the ecological niche modeling, determines the potential natural distribution of areas, and limits climatic factors using presence-only sampling data and environmental data layers [33].
Briefly, nineteen environmental variables (BIOCLIM) were provided by the web-based platform, WorldClim database (www.WorldClim.org) version 1.3, October 2004 [34], which includes climate databases from different sources. Environmental data were then further added from the FAO Map on Global Ecological Zones to build up the potential natural distribution model [35]. Then, the climate data information downloaded from www.WorldClim.org, was converted into a ESRI ASCII GIS (.asc) file by DIVA-GIS for the software. The data were first filtered, and outliers were then detected in DIVA-GIS (www.diva-gis.org) [36]. Then, all occurrence records and extreme values (minimum of 3 out of 19 bioclimatic variables examined) were rechecked for the inconsistency of their coordinates with administrative area level 1 and their climatic parameters based on the Reverse Jackknife method [37].
Modeling validation
The AUC of the receiver operating characteristic (ROC) was used to assess the constructed model [38]. The AUC value statistics as an independent threshold was ranged from 0.5 to 1, where a value below 0.5 was interpreted as a random prediction; 0.5–0.7 showed poor model performance; 0.7–0.9 showed moderate model performance, and a value above 0.9 was considered to have ’good’ discrimination abilities [39,40].
Results and discussion
Universality and evaluation of DNA barcodes for Salicornia sp
The success rates of PCR amplification and sequencing of the five barcode markers and their phylogenetic tree using single barcodes are shown in Fig 2 and S1 File. The three genomic and plastid DNA regions, including ITS2, trnH-psbA spacers and ycf genes and a combination of matK-rbcL were used as acceptable standard barcodes. Among the markers tested, rbcL had the highest amplification and recovery rates (98.90%), followed by trnH-psbA (82.42%), matK (80.21%), ycf (69.23%) and the rate for ITS2 was the lowest (65.93%). A low percentage amplification rate and recovery for ITS2 were due to the lack of suitable public primer, the incongruence of multiple copies, and or technical problems. Although the amplification success rate of ITS2 is low, it has been commonly used in studies of systematic and evolutionary botany. For example, ITS2 to be the most suitable region of more than 6,600 plant samples belonging to 4,800 species from 753 distinct genera tested, with 92.7 percent discrimination ability at the species level [41]. The sequence alignment method comprised 60 ITS2 sequences, 73 matK sequences, 90 sequences of rbcL, 75 sequences of trnH-psbA, and 63 sequences of ycf.
[Figure omitted. See PDF.]
Fig 2. The success rates of PCR amplification and sequencing of the five-barcode fragment and 26 different geographical locations.
https://doi.org/10.1371/journal.pone.0241162.g002
Furthermore, simple and optimized BLAST were carried for each of these sequences based on NCBI and ITS database. The efficiency of every single marker, as well as their combination for taxonomic identification, is presented in Fig 3. The highest level of identification using the single marker and simple BLAST at the species level, species group, genus, and family groups were obtained for ITS2 (74.55%), rbcL (98.78%), matK (62.19%) and ycf (83.08%) respectively while the discrimination efficiency using optimized BLAST was 70% (matK and trnH-psbA), 29.63% (ITS2), 94.44% (rbcL) and 53.96% (ycf) respectively (Fig 3). Likewise, data integration showed a higher identification at the level of the species group [42]. This group included S. persica, S. europea, S. patula, S. brachiate, S. herbacea and S. maritime suggesting species such as S. europea, S. patula and S. herbacea are similar or synonymous (The Plant List (2013). Version 1.1. Published on the Internet; http://www.theplantlist.org/ (accessed 1st January)).
[Figure omitted. See PDF.]
Fig 3.
Comparison of DNA barcoding efficiency using selected regions; a) simple and b) optimized BLAST based on NCBI and ITS database.
https://doi.org/10.1371/journal.pone.0241162.g003
Specific barcode for Salicornia sp
Automatic barcode gap discovery (ABGD) analysis, which assigns the sequences into potential species based on the barcode gap when the divergence within the same species is smaller than that among organisms from different species, detected a barcoding gap between the intraspecific and interspecific distance of barcode markers. Based on the distance-based approach as implemented in the software ABGD, different groups as candidate species were produced for rbcL, matk, ycf, psbA gene sequences. We analyzed K2P for distinct barcode gaps. No barcode gap is visible, including all samples and markers studied. For example, the study of matK identified the presence of two hypothetical species in all recursive partitions with prior intraspecific genetic divergence values between 0.46% and 10%, a result we considered four species with intraspecific divergence values below 0.28%. A barcoding gap between 3% and 16% of pairwise distances was also observed by the ABGD procedure. By the study of other markers, the two species hypothesis was also obtained. It should be noted that the ABGD algorithm yielded the same results when applied in the three implemented models (JC, K2P, and Simple Distance). Histograms of sequence divergence values and ranked distances among barcode sequences in Salicornia complex are shown in S2 Fig. The highest discrimination success in all combinations obtained for the S. europea species followed by S. brachiata (Fig 4). Further, the low success rate of markers might be related to the spreading and distribution of Salicornia germplasm species in Iran. This index is inverted with the amount of intra and inter-species gene flow. According to this theory, populations with a low distribution index are geographically or physiologically related. The low distribution index may cause the variants of neutral mutations to spread slowly throughout the entire population. Thus, the required time for lineage sorting for each gene locus increase by the natural genetic flow. In this case, the identification of specific barcodes is difficult for each species [43].
[Figure omitted. See PDF.]
Fig 4. Performances of matK, rbcL, trnH-psbA, ITS2 and ycf in resolving Salicornia species as different non-coding, coding, core and cpDNA combinations.
https://doi.org/10.1371/journal.pone.0241162.g004
Another consequence is that an increase in the penetration coefficient of a species in the adjacent species is caused by genetic flow, affecting the DNA barcoding method’s efficiency rate. In this study, the barriers and geographic distances, self-pollination, and cleistogamy reduced the gene flow rate among the plausible species of Salicornia populations. Later, the Iranian Salicornia germplasm has been studied with respect to the different nuclear and plastid sites with different heritability patterns in each group. As previously reported, the inheritance of plastid DNA is maternal in most angiosperms. Thus plastid haplotypes disperse through seeds [44]. In contrast, nuclear variants are distributed by seed and pollen and exhibited by biparental inheritance. The expected distribution range of pollen is higher relative to seed. Thus the genetic differentiation of the populations should be higher through plastid markers [45–47].
DNA barcodes variation survey on Salicornia population
Here, molecular variance analysis (AMOVA) and Fst index for core-barcode plastid (matK, rbcL) and nuclear marker (ITS2) showed that molecular variance within populations of the same species is greater than between populations indicating a higher genetic diversity (Fig 5). The largest difference between and within populations was observed in Dayer and Shahrbabbak genotypes based on ITS2 and matK markers, respectively. However, the rbcL fragment was not suitable for differentiation between populations. The highest Nei’s distance was obtained in Golmanxana and Qaregheshlaq genotypes based on ITS2 and rbcL markers. Moreover, the highest Fst index was observed in the Golmanxana genotype for the ITS2 marker, suggesting a low genetic flow rate compared to other genotypes. A higher Fst index conserves the allelic diversity needed for the conservation and exploitation of genetic resources. Further, the distribution of seed and pollen varies among populations, and in most cases, the genetic differentiation between populations depends on the geographical barriers and distance [48].
[Figure omitted. See PDF.]
Fig 5. Inter and intra molecular variance analysis and Fst index in Salicornia germplasm in Iran, based on plastid and nuclear markers.
https://doi.org/10.1371/journal.pone.0241162.g005
Mapping species diversity
The spatial analysis helps in better understanding and more accurate identification of biodiversity and developing strategies for identifying, managing, and exploiting genetic resources [49]. Outputs then provide vital information for prioritizing of protection and the Agro-Ecological Areas [50] finding. These findings may suggest that the new species can be grown in an environment of outside of its endemic region if additional water or soil nutrients are available. It should be noted that combining spatial results with other data is very useful in the management and exploitation of genetic reserves more efficiently [51]. The basis of spatial analysis in biodiversity is observed data from the sampled collection areas. In our study, the geographical data included identification codes, taxonomic names, geographic characteristics, and sampling locations (Fig 6).
[Figure omitted. See PDF.]
Fig 6. Distribution of the Salicornia populations in Iran based on the sampling areas.
The regions with the highest and lowest samples are depicted in red and dark green, respectively. The map was created using free software DIVA-GIS Version 7.5 (www.diva-gis.org).
https://doi.org/10.1371/journal.pone.0241162.g006
This information is commonly used for spatial diversity and distribution analysis. Biodiversity of plants is studied at three levels comprising of species-level and genetic level in a population (ecosystem). Here, Salicornia was studied at the species level or alpha diversity. In this case, the species level was observed in species diversity and evaluated in being absent or present conditions in each specific area. Further, the determination of diversity (including species, genotype, or ecotype) in different subunits was one of the most crucial aims in this study. The sub-units were areas that Salicornia species were found in a previous study. Consequently, diversity was studied at the alpha level, and the sampling areas were mapped using the observed number for each sample and their distributions (Fig 7). The regions were then classified into five groups based on their distributions.
[Figure omitted. See PDF.]
Fig 7.
Geographical distribution of Salicornia in Iran (a); the probability of presence increase from blue to red, (b) the effect of climates parameters on Salicornia distribution in Iran. The map was created using free open source MAXENT software version 3.4.1 (https://biodiversityinformatics.amnh.org/open_source/maxent/).
https://doi.org/10.1371/journal.pone.0241162.g007
The area with the highest and the lowest number of samples is depicted in red and dark green, respectively, in Fig 7. For example, if a plant protection program is considered, the areas with the highest alpha diversity (red) are in the top priority.
Ecological niche modelling of Salicornia sp
We used MAXENT to predict the potential distribution of the species of Salicornia populations based on the climatic, environmental variables. In the current study, 19 different climate variables affecting the distribution of species were evaluated [52] (Fig 8 and S2 File). As shown in Fig 8, isothermality (diurnal range mean/temperature annual range) *100, mean temperature of driest quarter, precipitation of warmest quarter, the maximum temperature of the warmest month, altitude, temperature seasonality, mean temperature of coldest quarter, minimum temperature of the coldest month, and annual mean temperature are the most critical variables to the Salicornia MAXENT model, based on jackknife test AUC. The average test AUC for each replicate runs were 0.994 with 0.002 standard deviation indicating MAXENT produces highly accurate predictions of AUC value, for all species.
[Figure omitted. See PDF.]
Fig 8.
Prediction model of Salicornia distribution (a) and limiting factors (b), based on geographical data and cluster analysis of sampled material. Different niches for protection and cultivation of Salicornia was shown with red to green color. The most competent areas for natural identification or cultivation of Salicornia when optimum water and soil parameter conditions are available shown in red colors while dark green colors show the least potential areas. The map was created using free software DIVA-GIS Version 7.5 (www.diva-gis.org).
https://doi.org/10.1371/journal.pone.0241162.g008
We also studied ecological niches on the conservation and utilization of genetic resources. This concept has been used to determine the priority of areas suitable for the conservation of wild species or the banking of natural genetic resources. Ecological niches are occupied area by a species in natural environmental conditions. Initial niches were also identified in this study. These areas have favorable climatic conditions for Salicornia growth and can be used to found these species.
Further, in the field trial condition, the coverage of the primary and ecological niches was confirmed. However, in some primary niches, Salicornia samples were not found due to the negative effect of climates or ecological parameters. Further, Geographic Information Systems (GIS) was used to model ecosystem nets based on available environmental data of each collected sample [15]. This model was drawn using DIVA and MAXENT software (Fig 8). The map showed different niches for the protection and cultivation of Salicornia. The most competent areas for natural identification or cultivation of Salicornia are shown in red colors where optimum water and soil parameters are available. However, the dark green colors show the least potential regions. Together, the model exhibited saline soil, and inefficient land might be used for industrial cultivation of Salicornia in Iran. Murray‐smith et al. (2009), previously applied MAXENT modeling and DIVA-GIS analysis to identify priority areas for the conservation of Myrtaceae. Their model showed observed species occurrences and predicted species occurrences and indicated complementarity analysis congruent in identifying areas with the most endemic species.
Conclusion
We validated genes in the field of DNA barcoding in Salicornia plants using matK, rbcL, trnH-psbA, ycf and ITS2, identifying species groups. Among the genetic markers tested, rbcL had the highest amplification and recoverability rates (98.90%), followed by trnH-psbA (82.42%), matK (80.21%), ycf (69.23%) and the rate for ITS2 was the lowest (65.93%). Data integration showed identification at the level of the species group was higher. This group included S. persica, S. europea, S. patula, S. brachiate, S. herbacea, and S. maritime. Molecular variance analysis and Fst index for plastid and nuclear markers showed genetic differences within the population of the same species is greater than between the populations indicating a high genetic diversity among the genotypes of each population. We also studied ecological niches on conservation and utilization of genetic resources to determine the priority of areas suitable for conservation of wild species or natural genetic resources banking. Our results provided valuable information on the diversity of specific geographical regions, conservation status of existing species, prioritization of conservation areas, and selection of regions for Agro-Ecological, which might be led to the development of industrial agriculture.
Supporting information
S1 Fig. Taxonomic identifications based on the combination of the identity score value.
(High identity: X≥95%; Medium identity: 90% ≤ X ≤ 95%; Low identity: X ≤ 90%).
https://doi.org/10.1371/journal.pone.0241162.s001
(TIF)
S2 Fig.
Automatic Barcode Gap Discovery (ABGD) analysis of barcoding markers, (a) Histogram of distances, (b) graph of ranked distances, (c) Automatic partition.
https://doi.org/10.1371/journal.pone.0241162.s002
(TIF)
S1 Table. DNA barcode primers sequences used in this study.
https://doi.org/10.1371/journal.pone.0241162.s003
(DOCX)
S1 File. The phylogenetic tree using a single fragment of ITS2, matK, rbcl, trn and ycf.
https://doi.org/10.1371/journal.pone.0241162.s004
(PDF)
S2 File. The primary data used in DIVA and MAXENT software.
https://doi.org/10.1371/journal.pone.0241162.s005
(PDF)
Citation: Zeinalabedini M, Khoshkholgh Sima NA, Ghaffari MR, Ebadi A, Farsi M (2021) Application of DNA barcodes and spatial analysis in conservation genetics and modeling of Iranian Salicornia genetic resources. PLoS ONE 16(4): e0241162. https://doi.org/10.1371/journal.pone.0241162
1. Bartels D, Sunkar R. Drought and Salt Tolerance in Plants. CRC Crit Rev Plant Sci. 2005;24: 23–58.
2. Witcombe JR, Gyawali S, Subedi M, Virk DS, Joshi KD. Plant breeding can be made more efficient by having fewer, better crosses. BMC Plant Biol. 2013;13: 1–12. pmid:23280064
3. Godfray HCJ, Beddington JR, Crute IR, Haddad L, Lawrence D, Muir JF, et al. Food security: the challenge of feeding 9 billion people. Science (80-). 2010;327: 812–8. pmid:20110467
4. Flowers TJ, Hajibagheri MA, Clipson NJW. Halophytes. Q Rev Biol. 1986;61: 313–337.
5. Flowers TJ. Improving crop salt tolerance. J Exp Bot. 2004;55: 307–319. pmid:14718494
6. Hadi MR, Khiyam-Nekoie SM, Khavarinejad R, Khoshkholgh Sima NA, Yavari P. Accumulation and role of ions (Ca2+, Mg2+, SO4-2) on salt tolerance in Triticum turgidum L. J Biol Sci. 2008;8: 143–148. Available: https://scholar.google.com/scholar?cluster=16416560739638406201&hl=en&oi=scholarr.
7. Öztürk M, Güvensen A, Sakçali S, Görk G. Halophyte plant diversity in the Irano-Turanian phytogeographical region of Turkey. Biosaline Agriculture and High Salinity Tolerance. Basel: Birkhäuser Basel; 2008. pp. 141–155.
8. Flowers TJ, Galal HK, Bromham L. Evolution of halophytes: multiple origins of salt tolerance in land plants. Funct Plant Biol. 2010;37: 604.
9. Hammer K, Khoshbakht K. Towards a ‘red list’ for crop plant species. Genet Resour Crop Evol. 2005;52: 249–265.
10. Ventura Y, Wuddineh WA, Shpigel M, Samocha TM, Klim BC, Cohen S, et al. Effects of day length on flowering and yield production of Salicornia and Sarcocornia species. Sci Hortic (Amsterdam). 2011;130: 510–516.
11. Katschnig D, Broekman R, Experimental JR-E and, 2013 undefined. Salt tolerance in the halophyte Salicornia dolichostachya Moss: Growth, morphology and physiology. Environ Exp Bot. 2013;92: 32–42. Available: http://www.sciencedirect.com/science/article/pii/S0098847212000895.
12. Glenn EP, Anday T, Chaturvedi R, Martinez-Garcia R, Pearlstein S, Soliz D, et al. Three halophytes for saline-water agriculture: An oilseed, a forage and a grain crop. Environ Exp Bot. 2013;92: 110–121.
13. Liu X, Xia Y, Wang F, Sun M, Jin Z, Science GW-F, et al. Analysis of Fatty Acid Compositions of Salicornia Europaea L. Seed Oil. Food Sci. 2005;2. Available: http://en.cnki.com.cn/Article_en/CJFDTOTAL-SPKX200502042.htm.
14. Gunning D. Cultivating Salicornia europaea (Marsh Samphire). Dublin, Ireland: Irish Sea Fisheries Board; 2016.
15. Kadereit G, Ball P, Beer S, Mucina L, Sokoloff D, Teege P, et al. A taxonomic nightmare comes true: phylogeny and biogeography of glassworts (Salicornia L., Chenopodiaceae). Taxon. 2007;56: 1143–1170. Available: https://www.ingentaconnect.com/contentone/iapt/tax/2007/00000056/00000004/art00013?crawler=true.
16. Muscolo A, Panuccio MR, Piernik A. Ecology, Distribution and Ecophysiology of Salicornia Europaea L. Springer, Dordrecht; 2014. pp. 233–240.
17. Petit RJ, Duminil J, Fineschi S, Hampe A, Salvini D, Vendramin GG. INVITED REVIEW: Comparative organization of chloroplast, mitochondrial and nuclear diversity in plant populations. Molecular Ecology. John Wiley & Sons, Ltd (10.1111); 2004. pp. 689–701. https://doi.org/10.1111/j.1365-294X.2004.02410.x pmid:15723661
18. Kadereit G, Piirainen M, Lambinon J, Vanderpoorten A. Cryptic taxa should have names: Reflections in the glasswort genus Salicornia (Amaranthaceae). Taxon. 2012;61: 1227–1239. Available: https://www.ingentaconnect.com/contentone/iapt/tax/2012/00000061/00000006/art00005?crawler=true.
19. Kaur G, Kumar S, Nayyar H, Upadhyaya HD. Cold Stress Injury during the Pod-Filling Phase in Chickpea (Cicer arietinum L.): Effects on Quantitative and Qualitative Components of Seeds. J Agron Crop Sci. 2008;194: 457–464.
20. Vorasoot N, Songsri P, Akkasaeng C, Jogloy S, Patanothai A. Effect of water stress on yield and agronomic characters of peanut (Arachis hypogaea L.). J Sci Technol. 2003;25: 283–288. Available: https://www.researchgate.net/profile/Chutipong_Akkasaeng/publication/26483170_Effect_of_water_stress_on_yield_and_agronomic_characters_of_peanut_Arachis_hypogaea_L/links/0c96051526ae30ecd8000000/Effect-of-water-stress-on-yield-and-agronomic-characters-of.
21. Teege P, Kadereit JW, Kadereit G. Tetraploid European Salicornia species are best interpreted as ecotypes of multiple origin. Flora—Morphol Distrib Funct Ecol Plants. 2011;206: 910–920.
22. Ingrouille MJ, Pearson J. The pattern of morphological variation in the Salicornia europaea L. aggregate (Chenopodiaceae). Watsonia. 1987;16: 269–281. Available: http://archive.bsbi.org.uk/Wats16p269.pdf.
23. Neves SS, Forrest LL. Plant DNA Sequencing for Phylogenetic Analyses: From Plants to Sequences. Humana Press; 2011. pp. 183–235. https://doi.org/10.1007/978-1-61779-276-2_10
24. Hollingsworth PM, Graham SW, Little DP. Choosing and Using a Plant DNA Barcode. Steinke D, editor. PLoS One. 2011;6: e19254. pmid:21637336
25. Reiahisamani N, Esmaeili M, Khoshkholgh Sima NA, Zaefarian F, Zeinalabedini M. Assessment of the oil content of the seed produced by Salicornia L., along with its ability to produce forage in saline soils. Genet Resour Crop Evol. 2018;65: 1879–1891.
26. Dong W, Liu J, Yu J, Wang L, Zhou S. Highly Variable Chloroplast Markers for Evaluating Plant Phylogeny at Low Taxonomic Levels and for DNA Barcoding. Moustafa A, editor. PLoS One. 2012;7: e35071. pmid:22511980
27. Altschul SF, Gish W, Miller W, Myers EW, Lipman DJ. Basic local alignment search tool. J Mol Biol. 1990;215: 403–410. pmid:2231712
28. Kearse M, Moir R, Wilson A, Stones-Havas S, Cheung M, Sturrock S, et al. Geneious Basic: An integrated and extendable desktop software platform for the organization and analysis of sequence data. Bioinformatics. 2012;28: 1647–1649. pmid:22543367
29. Ghorbani A, Saeedi Y, de Boer HJ. Unidentifiable by morphology: DNA barcoding of plant material in local markets in Iran. Bussmann R, editor. PLoS One. 2017;12: e0175722. pmid:28419161
30. Excoffier L, Lischer HE. Arlequin suite ver 3.5: a new series of programs to perform population genetics analyses under Linux and Windows. Mol Ecol Resour. 2010;10: 564–567. pmid:21565059
31. Weir BS, Cockerham CC. Estimating F -statistics for the analysis of population structure. Evolution (N Y). 1984;38: 1358–1370. pmid:28563791
32. Puillandre N, Lambert A, Brouillet S, Achaz G. ABGD, Automatic Barcode Gap Discovery for primary species delimitation. Mol Ecol. 2012. pmid:21883587
33. Warren DL, Seifert SN. Ecological niche modeling in Maxent: the importance of model complexity and the performance of model selection criteria. Ecol Appl. 2011;21: 335–342. pmid:21563566
34. Escudero A, Iriondo JM, Torres ME. Spatial analysis of genetic diversity as a tool for plant conservation. Biol Conserv. 2003;113: 351–365. (03)00122-8.
35. MacDicken K, Jonsson Ö, Piña L, Maulo S, Contessa V. Global forest resources assessment 2015: how are the world’s forests changing? 2016 [cited 30 Jul 2019]. Available: http://agris.fao.org/agris-search/search.do?recordID=XF2017001127.
36. Hijmans RJ, Guarino L, Cruz M, Rojas E. Computer tools for spatial analysis of plant genetic resources data: 1. DIVA-GIS. Plant Genetic Resources Newsletter. 2001. Available: https://books.google.com/books?hl=en&lr=&id=6QdBUMOtx0YC&oi=fnd&pg=PA15&dq=Computer+tools+for+spatial+analysis+of+plant+genetic+resources+data:+1.+DIVA-GIS.+Plant+Genetic+Resources+Newsletter&ots=pTiot0bO9v&sig=uK8ZvGqhL84JF125O2r0-meCbNE.
37. Chapman A. Principles and methods of data cleaning. GBIF; 2005. Available: https://books.google.com/books?hl=en&lr=&id=44gDJTEJoVIC&oi=fnd&pg=PA1&dq=Reverse+Jackknife+method+(Chapman+2005&ots=vSlT4ce3mp&sig=o8vwrA8B135r3bRtzzNwspnn05s.
38. Nyfløt LT, Sandven I, Stray-Pedersen B, Pettersen S, Al-Zirqi I, Rosenberg M, et al. Risk factors for severe postpartum hemorrhage: A case-control study. BMC Pregnancy Childbirth. 2017;17: 1–9. pmid:28049520
39. Wang F, Wang D, Guo G, Hu Y, Wei J, Liu J. Species delimitation of the Dermacentor ticks based on phylogenetic clustering and niche modeling. PeerJ. 2019;7: e6911. pmid:31123639
40. Peterson AT, Soberón J., Pearson RG, Anderson RP, Martínez-Meyer E, Nakamura M, et al. Ecological niches and geographic distributions (MPB-49). Princeton University Press; 2011. Available: https://books.google.com/books?hl=en&lr=&id=wyFnf9x9Vi0C&oi=fnd&pg=PR9&dq=Peterson+AT,+Soberón+J,+Pearson+RG,+Anderson+RP,+Martinez-Meyer+E,+Nakamura+M,+Araujo+MB.+2011.+Ecological+niches+and+geographic+distributions.+&ots=yntF45AZ2i&sig=29L9GkNxV9vfZytDCWdiyxk0lHA.
41. Chen S, Yao H, Han J, Liu C, Song J, Shi L, et al. Validation of the ITS2 Region as a Novel DNA Barcode for Identifying Medicinal Plant Species. Gilbert MTP, editor. PLoS One. 2010;5: e8613. pmid:20062805
42. Moritz C, Cicero C. DNA Barcoding: Promise and Pitfalls. Godfray Charles, editor. PLoS Biol. 2004;2: e354. pmid:15486587
43. Elith J, Phillips SJ, Hastie T, Dudík M, Chee YE, Yates CJ. A statistical explanation of MaxEnt for ecologists. Divers Distrib. 2011;17: 43–57. Available: https://onlinelibrary.wiley.com/doi/abs/10.1111/j.1472-4642.2010.00725.x%4010.1111/%28ISSN%291472-4642.species-distribution-models-in-conservation-biogeography.
44. Xu ZL, Ali Z, Yi JX, He XL, Zhang DY, Yu GH, et al. Expressed sequence tag-simple sequence repeat-based molecular variance in two Salicornia (Amaranthaceae) populations. Genet Mol Res. 2011;10: 1262–1276. Available: https://www.geneticsmr.com/sites/default/files/articles/year2011/vol10-2/pdf/gmr1321.pdf. pmid:21732290
45. Thakur P, Kumar S, Malik JA, Berger JD, Nayyar H. Cold stress effects on reproductive development in grain crops: An overview. Environ Exp Bot. 2010;67: 429–443.
46. Valentini A, Pompanon F, Taberlet P. DNA barcoding for ecologists. Trends Ecol Evol. 2009;24: 110–117. pmid:19100655
47. Ennos RA. Estimating the relative rates of pollen and seed migration among plant populations. Heredity (Edinb). 1994;72: 250–259.
48. Slatkin M. ene flow and the geographic structure of natural populations. Science (80-). 1987;236: 787–792.
49. Riccioli F, Fratini R, Boncinelli F, El Asmar T, El Asmar J-P, Casini L. Spatial analysis of selected biodiversity features in protected areas: a case study in Tuscany region. Land use policy. 2016;57: 540–554.
50. Vermeulen SJ, Campbell BM, Ingram JSI. Climate Change and Food Systems. Annu Rev Environ Resour. 2012;37: 195–222.
51. Kearney M, Porter W. Mechanistic niche modelling: combining physiological and spatial data to predict species’ ranges. Ecol Lett. 2009;12: 334–350. pmid:19292794
52. Esquinas-Alcázar J. Protecting crop genetic diversity for food security: Political, ethical and technical challenges. Nat Rev Genet. 2005;6: 946–953. pmid:16341075
You have requested "on-the-fly" machine translation of selected content from our databases. This functionality is provided solely for your convenience and is in no way intended to replace human translation. Show full disclaimer
Neither ProQuest nor its licensors make any representations or warranties with respect to the translations. The translations are automatically generated "AS IS" and "AS AVAILABLE" and are not retained in our systems. PROQUEST AND ITS LICENSORS SPECIFICALLY DISCLAIM ANY AND ALL EXPRESS OR IMPLIED WARRANTIES, INCLUDING WITHOUT LIMITATION, ANY WARRANTIES FOR AVAILABILITY, ACCURACY, TIMELINESS, COMPLETENESS, NON-INFRINGMENT, MERCHANTABILITY OR FITNESS FOR A PARTICULAR PURPOSE. Your use of the translations is subject to all use restrictions contained in your Electronic Products License Agreement and by using the translation functionality you agree to forgo any and all claims against ProQuest or its licensors for your use of the translation functionality and any output derived there from. Hide full disclaimer
© 2021 Zeinalabedini et al. This is an open access article distributed under the terms of the Creative Commons Attribution License: http://creativecommons.org/licenses/by/4.0/ (the “License”), which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited. Notwithstanding the ProQuest Terms and Conditions, you may use this content in accordance with the terms of the License.
Abstract
[...]the use of genotypes with high-yield and more resistant to biotic and abiotic stress has gradually replaced indigenous varieties and local genotypes, endangering access and exploitation of resources and food production for the future [3]. [...]low genetic diversity may reduce the opportunity to identify and use new sources responding to future challenges, including new pests and pathogens, as well as climate changes [4]. Today, soil salinity has increased due to unconventional cultivation and irrigation. [...]the level of arable land increased dramatically from 8 million to more than 220 million hectares, of which about 45 million hectares are salt-affected to varying degrees [5,6]. Salicornia sp. is widespread across Iran, mainly in central, western, and northern parts. [...]native Salicornia species selection might be a viable agricultural development strategy against Iran’s climate change. [...]a low number of morphological characteristics, phenotypic flexibility, breeding, and hybridization systems are amongst the factors that make Salicornia challenging to distinguish precisely in an area [20]. [...]the commercial importance and lack of description make Salicornia a priority for the circumscription of species, subtypes, ecotypes, and natural hybridization to develop it as a new crop [21].
You have requested "on-the-fly" machine translation of selected content from our databases. This functionality is provided solely for your convenience and is in no way intended to replace human translation. Show full disclaimer
Neither ProQuest nor its licensors make any representations or warranties with respect to the translations. The translations are automatically generated "AS IS" and "AS AVAILABLE" and are not retained in our systems. PROQUEST AND ITS LICENSORS SPECIFICALLY DISCLAIM ANY AND ALL EXPRESS OR IMPLIED WARRANTIES, INCLUDING WITHOUT LIMITATION, ANY WARRANTIES FOR AVAILABILITY, ACCURACY, TIMELINESS, COMPLETENESS, NON-INFRINGMENT, MERCHANTABILITY OR FITNESS FOR A PARTICULAR PURPOSE. Your use of the translations is subject to all use restrictions contained in your Electronic Products License Agreement and by using the translation functionality you agree to forgo any and all claims against ProQuest or its licensors for your use of the translation functionality and any output derived there from. Hide full disclaimer