ITS secondary structure reconstruction to resolve taxonomy and phylogeny of the Betula L. genus

Andrii S. Tarieiev1, Oliver Gailing1,2 and Konstantin V. Krutovsky1,2,3,4,5

1 Department of Forest Genetics and Forest Tree Breeding, Georg-August University of Göttingen, Göttingen, Germany 2 Center for Integrated Breeding Research, Georg-August University of Göttingen, Göttingen, Germany 3 Laboratory of Forest Genomics, Genome Research and Education Center, Institute of Fundamental Biology and Biotechnology, Siberian Federal University, Krasnoyarsk, Russia 4 Laboratory of Population Genetics, N.I. Vavilov Institute of General Genetics, Russian Academy of Sciences, Moscow, Russia 5 Department of Ecosystem Science and Management, Texas A&M University, College Station, TX, United States of America

ABSTRACT The taxonomy and phylogeny of the Betula L. genus remain unresolved and are very difficult to assess due to several factors, especially because of frequent hybridization among different species. In the current study, we used nucleotide sequences of two internal transcribed spacer regions (ITS1 and ITS2), which are commonly used as phylogenetic markers. In addition to their nucleotide variation we reconstructed their secondary structure and used it to resolve phylogenetic relationships of some birch species. We explored whether consideration of secondary structure in phylogenetic analyses based on neighbor-joining, maximum parsimony, maximum likelihood, and Bayesian inference methods would help us obtain more solid support of the reconstructed phylogenetic trees. The results were not unambiguous. There were only a few clades with higher support when secondary structure was included into analysis. The phylogenetic trees generated using different methods were mostly in agreement Submitted 16 October 2020 with each other. However, the resolving power of these markers is still insufficient to Accepted 12 January 2021 reliably discriminate some closely related species. To achieve this aim more reliably Published 23 March 2021 there is a need for application of modern genomic approaches in combination with Corresponding author traditional ones. Konstantin V. Krutovsky, [email protected] goettingen.de Subjects Computational Biology, Evolutionary Studies, Genetics, Plant Science, Taxonomy Academic editor Keywords Betula, Taxonomy, Phylogeny, ITS, Secondary structure, Phylogenetics Gabriele Casazza Additional Information and Declarations can be found on INTRODUCTION page 16 Birches (Betula L. genus) are common pioneer trees and shrubs of temperate and boreal DOI 10.7717/peerj.10889 zones in the Northern Hemisphere. Trees classified within this genus are highly diverse, Copyright especially in the Old World. Birches occupy habitats predominantly in cool and moist 2021 Tarieiev et al. regions, including peatlands, stream banks, lakeshores, cool, damp woods, and moist Distributed under slopes in cool coves. Most of them represent organisms adapted to low temperatures and Creative Commons CC-BY 4.0 able to grow under different extreme conditions (bogs, high altitude mountains, sands, OPEN ACCESS

How to cite this article Tarieiev AS, Gailing O, Krutovsky KV. 2021. ITS secondary structure reconstruction to resolve taxonomy and phylogeny of the Betula L. genus. PeerJ 9:e10889 http://doi.org/10.7717/peerj.10889 subarctic zone, etc.) and on very poor soils, frequently as the first colonizing pioneer woody plants (Ashburner & McAllister, 2013). The taxonomy and exact phylogenetic relations in the Betula L. genus are still unresolved and very difficult to assess because of several reasons. Number of taxa and their range vary greatly in different studies or sources from ∼30 (De Jong, 1993) to ∼120 (World Checklist of Selected Plant Families; Ashburner & McAllister, 2013) species. Birches comprise taxonomically a notoriously complex group because of their high morphological variation, phenotypic plasticity, and frequent hybridization (especially introgressive one) (Johnsson, 1945; Johnsson, 1949; Brown, Kennedy & Williams, 1982; Anamthawat-Jonsson & Tomasson, 1999; Anamthawat-Jonsson & Thorsson, 2003; Thorsson, Salmela & Anamthawat-Jonsson, 2001; Thorsson et al., 2010; Palmé et al., 2004; Bona, Petrova & Jadwiszczak, 2018). Earlier classifications were mostly based on morphological traits (Regel, 1865; Schneider, 1915; Schneider, 1916; De Jong, 1993; Walters, 1993; Chen, Manchester & Sun, 1999; Skvortsov, 2002; Tzvelev, 2002; Ashburner & McAllister, 2013). The internal classification of the genus (subgenera, sections, and subsections) was changed several times, but even till now it is still mostly based on morphological traits. However, they are often insufficient to discriminate different birch species because of their high phenotypic plasticity and frequent hybridization, which blurs boundaries between different taxa. Only a few birch taxa were studied using different molecular genetic markers, such as rbcL and ITS sequences (Chen, Manchester & Sun, 1999), nuclear ADH and chloroplast matK sequences (Jarvinen et al., 2004), nuclear ribosomal DNA sequences (Li, Shoupa & Chen, 2005), nuclear nitrate reductase DNA sequences (Li, Shoupa & Chen, 2007), other multiple chloroplast and nuclear regions (Singevar et al., 2020) and AFLP markers (Schenk et al., 2008), but sampling sizes were quite low in most studies. In addition, there was a recent study of the genetic structure and hybridization of some birch species across Eurasia that provided new data regarding the hybridization patterns and the origin of B. pubescens (Tsuda et al., 2017). Some birch taxa were recently reassessed with propositions to lower taxonomic ranks using also internal transcribed spacers (ITS) as molecular genetic markers (e.g., Kuneš et al., 2019; Tarieiev et al., 2019; Jadwiszczak et al., 2020; Linda et al., 2020). Genome-wide markers inferred by restriction site associated DNA sequencing (RAD- seq) were used to study the phylogeny of 27 diploid and 31 polyploid birch taxa (Wang et al., 2021) and specifically Costatae section with description of a new species, Betula buggsii N. Wang (Wang et al., 2020). Complete reference genome sequence assemblies are also available for three birch species, B. nana L. (Wang et al., 2013), B. pendula Roth (Salojärvi et al., 2017) and B. platyphylla (Chen et al., 2021) respectively. ITS are among the most common phylogenetic markers in plants including also several studies on birch genus (Chen, Manchester & Sun, 1999; Li, Shoupa & Chen, 2005; Wang et al., 2016; Kuneš et al., 2019; Tarieiev et al., 2019). There is a large number of ITS sequences deposited and freely available in the NCBI GenBank, ITS2, and EMBL databases. The most complex phylogenetic reconstruction of ∼60 birch taxa (233 sequences in total) was done relatively recently using ITS sequences as phylogenetic markers (Wang et al., 2016). However, this phylogenetic study was unable to discriminate some closely related and hybrid taxa. The objective of our study was to test whether inclusion of secondary structure

Tarieiev et al. (2021), PeerJ, DOI 10.7717/peerj.10889 2/21 of ITS would improve phylogenetic analysis. In our study we inferred secondary structure of both ITS1 and ITS2.

MATERIALS & METHODS ITS markers We analyzed 223 complete nucleotide sequences of ITS1 that represented 73 birch taxa* (57 species, including two species with hybrid origin, 6 subspecies, 9 varieties, and one form) and 263 sequences of ITS2 that represented 78 birch taxa* (60 species, including two with hybrid origin, one interspecific hybrid, 6 subspecies, 10 varieties, and one form) retrieved from the NCBI GenBank and compared for secondary structure reconstruction with the ITS2 Database (for ITS2 only; Ankenbrand et al., 2015) (see Table S1). In addition, we sequenced several ITS sequences representing the two endemic birch species from Ukraine B. borysthenica Klokov (Klokov, 1946) and B. klokovii Zaverucha (Zaverucha, 1964), dark-barked birches B. pubescens var. sibakademica Baranov (Kuzeneva) (≡B. pubescens f. sibakademica (Baranov) Tarieiev), (Baranov, 1924; Kuzeneva, 1936; Tarieiev et al., 2019) B. kotulae Zaverucha (=B. pendula f. obscura (Kotula) Tarieiev), (Zaverucha, 1964; Tarieiev et al., 2019) B. atrata Domin (Domin, 1927) and carpathian birch B. pubescens ssp. carpatica from Poland. We used specimens originating from natural habitats of Ukraine (especially endemic species from type localities (locus classicus) and type specimens listed in Olshansky et al., 2016), the Botanical Garden of the Georg-August University of Göttingen, and other herbaria–GOET, KW, LE and the former herbarium of Taras Shevchenko National University of Kyiv (see Table S2). Birch taxa are listed in Tables S1 and S2 according to the NCBI taxonomy with more detailed explanation provided in Data S1. However, the NCBI taxonomic ranks do not always reflect the current changes in taxonomy and, therefore, could be wrong or outdated in some cases. ITS1 and ITS2 from ribosomal gene clusters were amplified and sequenced in nine birch species either as a single fragment that included partial 18S, ITS1, 5.8S, ITS2, and partial 28S using universal PCR primers ITS5 and ITS4 (White et al., 1990) or in parts using ITS5 and ITS2 primer pairs modified for birch species to amplify the ITS1 region and ITS3 and ITS4 primer pairs to amplify the ITS2 region, respectively. PCR products were sequenced on an ABI 3130x Genetic Analyzer (Applied Biosystems, Thermo Fisher Scientific Inc., Waltham, MA, USA). Using BioEdit software (Hall, 1999) these sequences were checked visually for quality and aligned together with 223 and 263 complete sequences of ITS1 and ITS2 retrieved from the NCBI GenBank and ITS2 databases and presented in Tables S1 and S2, respectively. Secondary structure (Mai & Coleman, 1997; Coleman, 2000; Coleman, 2003; Coleman, 2007; Coleman, 2009; Müller et al., 2007; Wolf et al., 2013) was inferred for all analyzed sequences using mfold web server (Zuker, 2003). Obtained structures were checked for correctness (e.g., correct position of conservative motives, U-U mismatch, and 5.8S-28S rRNA hybrid stem formation were manually inspected for ITS2) and Gibbs energy value. Structures were compared with models proposed in the ITS2-Database (Koetschan et al., 2010; Koetschan et al., 2012) and aligned automatically

Tarieiev et al. (2021), PeerJ, DOI 10.7717/peerj.10889 3/21 with manual adjustment afterwards using the 4SALE program (Seibel et al., 2006; Seibel et al., 2008). Alignments are presented in Data S2 and S3. Obtained structures were visualized using VARNA v3-93 (Darty, Denise & Ponty, 2009), and then their features were colored using Inkscape v. 0.92.3 (Free Software Foundation, Inc.; https://inkscape.org). Phylogenetic relationships were reconstructed using Bayesian inference implemented in MrBayes 3.2.6. (Huelsenbeck & Ronquist, 2001; Ronquist & Huelsenbeck, 2003; Ronquist et al., 2012) on the same sequence set with exclusion of identical sequences. Secondary structure was inferred by applying two different models: doublet model for paired regions and 4 by 4 with GTR +I for unpaired ones, respectively. The analysis was performed in four chains with sampling frequency 100 and different number of iterations –20,000,000, 50,000,000, and 100,000,000. Burnin was performed first separately in each run by monitoring the standard error value, convergence was considered to be reached when the standard error decreased below 0.01. The obtained final phylogenetic tree was visualised using FigTree (http://tree.bio.ed.ac.uk/software/figtree) and TreeGraph2 (Stöver & Müller, 2010). To analyse the impact of rRNA secondary structure on the topology of the phylogenetic tree and support of the clades the following four methods of phylogenetic analysis were performed: (1) Bayesian analysis using an identical number of iterations and 4 by 4 GTR+I+G model suggested by the ModelTest in PAUP v.4.0a (http://paup.phylosolutions. com; Swofford, 2003) for all the regions, (2) neighbor-joining (NJ) clustering in ProfDist (Friedrich et al., 2005; Wolf et al., 2008) using GTR model, 1000 bootstrapping iterations and rate matrix Q for ITS2 in case of considering secondary structure or rate matrix Q for Chlorophyta in case of not considering secondary structure, (3) maximum parsimony (MP) analysis using PAUP v.4.0a with GTR+I+G model, 1000 bootstrap iterations, and 12-letter coding of secondary structure performed using 4SALE, 4) maximum likelihood (ML) analysis using MEGA X without considering secondary structure and R script provided with 4SALE package (http://4sale.bioapps.biozentrum.uni-wuerzburg.de/mlseqstr.html) together with 12-letter coding to consider the structure. It should be noted that 12-letter coding does not include IUPAC codes for ambiguous nucleotides, and, therefore, some information is lost, especially in case of considering hCBCs. Pairwise topology of phylogenetic trees was adjusted manually for better comparison by merging and rearranging some clusters with statistically insignificant support, and respective drawings were generated using Inkscape v. 0.92.3 (Free Software Foundation, Inc.; https://inkscape.org).

RESULTS Secondary structure reconstruction Structural variation was observed for different groups of birch taxa. The consensus structural models comprising all variants were reconstructed for the ITS1 and ITS2 markers (Fig. 1). The ITS1 structure consists of five helices: four main (H1, H2, H3, and H4) plus one additional (Ha). Additional helix Ha was present in four structural variants (Fig. 2): the A

Tarieiev et al. (2021), PeerJ, DOI 10.7717/peerj.10889 4/21 C A A W A C G U U C C Y U H3 w U G U S G U G g U C G U G C H1 C R G C C G C U C G C a S G A C C U R G C C G K G K C U C R U Y U H2 A R G Y G G A C U C Y G G Y G C G A Ha A G G C C A G G U C C G Y A C G U Y G C C Y G H G G A G Y A Y C R N U A G A G C C C G C Y A W Y K U U C A G G G A G Y A G U A G C Y A A G C S Y A R R Y C K G U

C A C Y

Y A G G C G c G Y A U U Y C A C C Y C M C U G K A G A ITS1 A Y A A G K A C Y G U G C G K U C Y H2 R C C W G a Y S C R A U U C Y U G A U 1 G C U C G U G G U A G C G G Y A W U G U A A C G A G G C C C Y C 218 C C C G G G C C C U C A A G Y C H3 G C A A G G U Y M G U Y G R U A U S U Y U ITS2 C U A U U A A G G A Y A U G C C U U G H4 G U G Y A A C A Y U G A U A U Y G C Y A G K G U Y c A G U C A A M A G C G Y W S G G C Y G A V hemicompensatory base change (hCBC) C G H C U G R G Y C A C C C G R C R C C R structural base change H 5.8S U A 28S Y W R U 1 233 C A G C R g C rRNA hybrid stem Y C N non-structural base change G C G G G C U C R c G C wobble base pair R A C G A U

G Y G C U C Y U G U canonical or wobble base pair A C U C G G C G G G G U C G A conditional bond U A A H4 H1 A K A H A

Figure 1 Consensus models of the secondary structure comprising all variants of ITS1 and ITS2 mark- ers for the entire Betula L. genus. Indels are depicted by small letters. Full-size DOI: 10.7717/peerj.10889/fig-1

variant was common in most samples (Fig. 2A), B was specific for B. globispica (AB243883.1, KT308904.1, KT308905.1, AY761111.1), B. chinensis (KT308917.1, KT308918.1), B. fargesii (KT308906) from section Asperae, subsection Chinenses, B. glandulosa (AY761110.1) from section Apterocaryon, B. lenta (AY761115.1, AY352330.1), B. medwediewii (AY761120.1), and B. uber (AY352334) from section Lentae (Fig. 2B), C was specific for B. alnoides (AY763114.1) and B. luminifera (KT308944.1) from the section Acuminatae (Fig. 2C), and D occurred only in two haplotypes of B. apoiensis (AB243908.1 and AB243912.1)(Fig. 2D). The H1 helix was the most variable in ITS1 and represented by seven structural variants (Fig. 3). However, the last two ones (F and G, respectively) could be reduced to the first and second one (A and B) depending on the nucleotide in ambiguous sites N and R. The H2 helix was represented by two structural variants A and B (Fig. 4). B was specific only for B. insignis (KT308927.1, KT308928.1, KT308929.1). The H3 helix was the same for all analyzed sequences and without any structural variation. Four structural variants were

Tarieiev et al. (2021), PeerJ, DOI 10.7717/peerj.10889 5/21 ITS1 additional helix Ha C C C C C C C C A C C G A C C A G C G G G G U G A G G Y A U G C A C A G A U G A A A A A C A A C A C A A C A G Y G C G A A C A U A U A U G C

Y G C G C G A U

G U G U G U C G

R U A U A U G U

C C G A C C G A C C G A C A U S 1 A 28 1 B 28 1 C 28 1 D 26

Figure 2 Structural variants of additional helix Ha. Full-size DOI: 10.7717/peerj.10889/fig-2 ITS1 Helix 1 U C U C U C U C A U U C W U C U A U A U U C A U G Y U C G C G C G C A U R C C G C G C A C G G C C G C G G C G C G C G C C G G C G C G C G C G C G C G C G C G C G C G C G C G C G C G C G C G C G C G C G C G C G C G C G U G U K K G U G G G C G U G U R U G U G U G U U U G U Y G Y G Y G U G G U C G N G

G Y G Y G C G C U G G U G C

G C G C G C G C G C G C G C

G C G Y G C G C G C G C G C

G Y G Y G C G C G C G C G C

G C G C G C G C G C G C G C

C U G A C U G A Y U G A C U G A G G C G C U G A C U G A 1 A 36 1 B 36 1 C 34 1 D 36 1 E 34 1 F/A 36 1G/B36

Figure 3 Structural variants of helix 1 (H1) in ITS1. Full-size DOI: 10.7717/peerj.10889/fig-3

observed for the H4 helix (Fig. 5). The B variant was species-specific for B. nigra (see detailed information about distribution of other variants in Data S4). The ITS2 marker had no structural variation for H2 and H4 helices. Five structural variants were observed for the H1 helix (Fig. 6) (see also details in Data S3). The helix 3 (H3) was the most variable with 10 structural variants in total (Fig. 7). Further details

Tarieiev et al. (2021), PeerJ, DOI 10.7717/peerj.10889 6/21 ITS1 Helix 2 U Y C G U C A U C G

C G A U A C C G G A A C A U G A G C A U

G C G C

G C G C

A U G C

Y G A U

G C C G

G Y G C

A Y G A C G C U 1 A 32 1 B 30

Figure 4 Structural variants of helix 2 (H2) in ITS1. Full-size DOI: 10.7717/peerj.10889/fig-4

regarding distribution of different variants among analyzed sequences are presented in Data S5. There were no compensatory base changes (CBCs) detected in both structural models of markers across all samples. However, there were several hemicompensatory base changes (hCBCs): in total, 16 hCBCs sites were detected in ITS1 and 15 in ITS2, respectively (more detailed information is in Data S6 and S7). The exact distribution of those hCBCs is depicted in Figs. 1–7 in red color.

Neighbor-joining (NJ) phylogenetic analysis NJ phylogenetic tree with consideration of the secondary ITS structure showed higher statistical support of several clusters and more detailed clustering compared to the NJ trees without such consideration (42 vs. 35 supported clusters, respectively, considering 50% as a cut-off for bootstrap support; Fig. 8, Figs. S1–S2). In addition, 14 clusters that were supported in both trees had higher bootstrap values when the secondary structure was considered. However, the resolution power remained insufficient to delineate some of closely related species and hybrids. In general, clustering was mainly in agreement with genus taxonomy, but not always.

Tarieiev et al. (2021), PeerJ, DOI 10.7717/peerj.10889 7/21 ITS1 Helix 4 A A A A A A A A A C A A G A C C G C G G G G G G G G G G G C U G G C U C U G U G C U U G U G C U U G C U C U Y G C U C G C G C G G C G C G C G C C G C G C G C U C G C U C U G G C U G G G G N G G C G G C C G Y C G G C G G G C C A C C G G G U G C U G U U G C R U C A C A C A C G C C G C G A G G G G G G G U U U Y U A U U V G C A G G G G G A A A A A A A A G G A G G U A U U A U A A A A A C C C C A A A A A A A A G G G G U U U U C C C C U U U U A A A A G G G G A A A A U U U U U U C U U C U U C U U C

U A U A U A U A U U U U M G C G A G C G A A A A A A A A

A A A A A A A A G C G C G C G C

A G C A A G C A A G C A A G C A 1 A 77 1 B 77 1 C 77 1 D 77

Figure 5 Structural variants of helix 4 (H4) in ITS1. Full-size DOI: 10.7717/peerj.10889/fig-5

ITS2 Helix 1

A A A A C C A A G G G G U U G G A U U G G A A A A C C A G A A A H G A A G C G Y G C C C G C G G U A Y G A U G U G G U K U C A C G G U C U C C C C G A A U A C U A C G A G G G A C A G R U A U G U G C C U A A G C G G C G U G G A G C A A U U C C U U G G C G C C U G A G G U C A G G G G C C C A A C G G G G C G C G C C C C C C G C C C G C C G C Y C C U C C C C C C C C C C U H U A U A A A A A A A A A A G G G G G W A A A C Y C U C U C U C U C U C G C G C G C G C G C G

C G C G C G C G C G C G

C G C G C G C G C G C G

C G C G C G C G C G C G

G C G Y G C G C G U G C

U A U A U A U A U A U A

G U G A G U G A G U G A G U G A G U G A G U G A 1 A 53 1 B 53 1 C 52 1 D 53 1 E 53 1 F 52

Figure 6 Structural variants of helix 1 (H1) in ITS2. Full-size DOI: 10.7717/peerj.10889/fig-6

Tarieiev et al. (2021), PeerJ, DOI 10.7717/peerj.10889 8/21 ITS2 Helix 3 C U A C A C C A C A C A C A C A A W A G A U W A C A A A A A A C A A A A U A A A G A A G G A A G G A G U G A A A A A U C C U U U G A U U U C C G A C C C U C U G C C U U U U U C U U Y C C C C C C C C G C C C U C C G C G G U G C G G G U C C G C G U G U G U G U G G U G U G U U G U C G U U C U C G U U C C G C G U U C U U U G G G U G G C G G G G U G U G C G U G U G U G G G G G U G G G G G U U G G G U G U G G G G U U G U G G U U C U G C G C U C C C G U C C C U C C U C C C U C C U C U U C U C U C C U C U C A C C C A C A C A A A U A A C A C C A C A C A C C A C C G A C A C G C G A A C C G C G A C C G C G A C C C C G R U C G A U A U A U A U A U A U A U A U G C A U G C G C G C G C G C G C G C G C C G G C C G C G C G C G C G C G C G C G R U C G A U A U A U A U A U A U A U A U C G A U C G C G C G C G C G C G C G Y C C G C C C C C C C G C C C C C Y C C G G C G G G G C C G G G G C G G G G G G G U G U U U U U G U Y U C G C C C C C U U G A U G G G G G G G C G G G G G C G G G C N C C C G C C C A C M G C C C C C G C C C C G C G G C C G C G C G C C G G G G C G C G G C G C C C G C C G C G C Y C C G C G C C A C A G A C A A C A Y A G A G C G C G G G G G C G A G G G G G A U G U A U G U U U C U C G C G C C C U C Y C C G C U C C C G C G U G C G G G C G G G G C G C G G G C G G C G C C G G C G C G C G U C G U G U U U U U Y G C G C G C C U U C C U C C C C G C C G C C C C G Y C U U U C U U U Y U U C C U C C C C C C U C U C C U C C C C U C U C U U U U U U C G M G A G A C G A G A G A G A A U G A U G A U A U A U A A U A U A U A U G A A G G Y G C G C A U G C G C G C G C G C G C Y G C G C G G C C G C G C G C G C G C W

A R U G A G U G A G U G A G U G A G C G G G A G U G A G U G A G U G A G U G A G U G 1 A 85 1 B 85 1 C 84 1 D 85 1 E 85 1 F 85 1 G 84 1 H 84 1 IJ85 1 85

Figure 7 Structural variants of helix 3 (H3) in ITS2. Full-size DOI: 10.7717/peerj.10889/fig-7

Maximum parsimony (MP) phylogenetic analysis MP phylogenetic trees demonstrated the same pattern as the NJ trees: there were more supported clusters when the secondary structure was considered compared to the MP tree without such consideration: 27 vs. 18, respectively (Fig. 9, Figs. S3–S4). In addition, twelve clusters that were supported in both trees demonstrated higher bootstrap values when the secondary structure was considered. However, the resolution power of the MP analysis was lower than in the NJ analysis. Maximum likelihood (ML) phylogenetic analysis MP phylogenetic trees demonstrated 34 supported clusters when secondary structure was considered, but only 28 ones when traditional analysis without consideration of secondary structure was used (Fig. 10, Figs. S5–S6). In addition, twelve clusters that were supported in both trees demonstrated higher bootstrap values when the secondary structure was considered. Bayesian phylogenetic analysis Bayesian phylogeny reconstructions showed some changes in phylogenetic trees; both changes in the topology and support values (posterior probability) were observed when the secondary structure of ITS was considered. Clusters supported with posterior probability values >95% are highlighted in red in Fig. 11, Figs. S7–S11. Several Bayesian trees were constructed using different number of iterations (20,000,000, 50,000,000, and 100,000,000). At 20,000,000 iterations the tree with consideration of the secondary structure had greater statistical support values for several clades (Figs. S10, S11) compared to the tree without consideration of the secondary structure (Fig. S9). However, with increasing iteration number this difference vanished, and for 100,000,000 iterations the situation is opposite –the tree without structure had greater statistical support (Fig. 11, Figs. S7, S8). At the same time, the stability of the tree in terms of topology, number

Tarieiev et al. (2021), PeerJ, DOI 10.7717/peerj.10889 9/21 1 ·············································································································································AY352329 Alnus viridis································································································································· 0.99 ············································································································································AY352306 Alnus cordata······························································································································· 0.64 ······················································································································KT308934, KT308935, AY761112, AY761113 B. grossa·············································································································· 0.64 ···············································································································································AB243892 B. grossa·································································································································· 0.59 ··································································································AY761115, AY352330, KT308936, B. lenta; KT308937, KT308938 B. lenta f. uber······························································································· 0.85 ·······························································································AY761120, KT308930, KT308931 B. medwediewii; KT308932, KT308933 B. megrelica·························································································· ···········································································································································KT308927 B. insignis····································································································································· 0.99 0.67 0.52 ···············································································································KT308928 B. insignis; KT308929 B. insignis ssp. fansipanensis··········································································································

0.86 ·······························································································································KT308925, AY761100 B. alleghaniensis························································································································· 0.86 0.64 ·········································································································································KT308926 B. murrayana··································································································································· 0.64

··································································································································KT308911, KT308912 B. bomiensis···························································································································

0.74 0.99 0.52 0.57 1 0.51

··········································································································································KT308943 B. luminifera··································································································································· 0.51 ··········································································································································AY761116 B. luminifera··································································································································· 0.57 0.61 0.73 0.74 0.88 ··········································································································································AY761117 B. luminifera··································································································································· 0.61 0.53 ··········································································································································AY763113 B. luminifera··································································································································· 0.72 0.61 0.84 ·········································································································································KT308942 B. hainanensis·································································································································· ·······································································································································KT308941 B. cylindrostachya······························································································································ 0.75 ···········································································································································KT308944 B. luminifera·································································································································· 0.91 ············································································································································AY763114 B. alnoides··································································································································· ·····································································································································AY761101, KT308940 B. alnoides··························································································································· ·············································································································································AJ783641 B. alnoides·································································································································· ·······································································································AB243890, KT308945, KT308946, AY761118, AY761119 B. maximowicziana································································································

0.63 0.94 0.56 0.58 0.72 0.5 0.66 ···········································································································································AY761110 B. glandulosa································································································································· 0.68

0.58 ······························································································································AY761106, KT308907, KT308908 B. corylifolia················································································································ 0.91 0.9 ············································································································································AB243885 B. corylifolia································································································································· ············································································································································AB243884 B. corylifolia································································································································· 0.54 0.86 0.89

0.96 0.91

0.99 0.86 0.99 0.77 0.58 0.67 0.66 0.61 0.77

0.97 0.99

0.57 ··································································································KT308951 B. utilis var. jaquemontii; KT308952, KT308953, KT308961 B. ashburneri·························································································· 0.55 ·····························································································································KT308949 B. utilis; AY761099 B. albosinensis··················································································································· 0.5 ················································································································································KT308948 B. utilis···································································································································· ·································································································································FJ011780 B. utilis; FJ011773 B. davurica····················································································································· ················································································································································AY761134 B. utilis···································································································································· ·····························································································································KT308923, KT308950 B. utilis var. occidentalis················································································································· ········································································································································KT308955 B. utilis var. prattii······························································································································ ···································································································································KT308924, KT308954 B. albosinensis······················································································································· ······························································································································KT308947 B. albosinensis var. septentrionalis·················································································································· ··············································································································AB243886 - AB243888, KT308956, KT308957, AY761108 B. ermanii···································································································· ··············································································································································AB243889 B. ermanii·································································································································· 0.52

0.55

····································································································································AY761121, KT308978 B. michauxii·························································································································· 0.72 ···············································································································································KT309025 B. humilis·································································································································· ·······································································································································KT309024, KT309026 B. humilis·························································································································· ········································································································································AJ783643, AY761114 B. humilis·························································································································· ···············································································································································AY761109 B. fruticosa································································································································ ················································································································································AJ783642 B. pumila································································································································· 0.69 ················································································································································AY761131 B. pumila································································································································· 0.77 ··························································································································KT309021 B. pumila; KT309019 B. nana ssp. exilis·············································································································· ·································································································································AY352336, AY761122, KT309020 B. nana···················································································································· ······················································································································AB243897, AB243898, KT309022, KT309023 B. ovalifolia·········································································································· ············································································································································AB243895 B. middendorffii····························································································································· ·············································································································································KT309017 B. glandulosa······························································································································ ·······································································································································KT308966, AY761132 B. raddeana······················································································································· 0.64 ···············································································································································FJ011770 B. davurica································································································································

0.6 ··············································································································································AB243894 B. davurica································································································································ 0.84 ·····························································································································KT308963 B. dahurica; AY761105 B. chinensis················································································································· ·······················································································································AB243893, FJ011771, FJ011772, KT308962 B. dahurica·········································································································· 0.63 0.51 ·············································································································································KT308995 B. glandulosa······························································································································ 0.6

0.54 0.67 0.69

0.65 0.65 0.56

Subgenus Betula Aspera Acuminata Nipponobetula

Section Costatae Apterocaryon Betula Dahuricae Asperae Lentae Acuminatae Nipponobetula

Subsection Asperae Chinenses B. × minor: B. glandulosa (Apterocaryon) × B. cordifolia (Betula)

Figure 8 Two neighbor-joining (NJ) phylogenetic consensus trees based on 1,000 iterations without (left) and with (right) consideration of ITS secondary structure generated using ProfDistS with GTR model and rate matrix Q. They are presented also separately in Figs. S1–S2, respectively. Clusters with bootstrap value less than 0.5 (50%) were collapsed. Subgenera, sections, and subsections are marked by different colors according to Ashburner & McAllister (2013). Identical taxa are connected by lines. Full-size DOI: 10.7717/peerj.10889/fig-8

Tarieiev et al. (2021), PeerJ, DOI 10.7717/peerj.10889 10/21 1 ·····································································································································AY352306 Alnus cordata······································································································································ 1 ······································································································································AY352329 Alnus viridis······································································································································· 0.60 ·······································································································································AB243892 B. grossa········································································································································· 0.63 ··················································································································KT308934, KT308935, AY761112, AY761113 B. grossa················································································································· 0.56 ·································································································AY761115, AY352330, KT308936, B. lenta; KT308937, KT308938 B. lenta f. uber································································································ 0.79 ·····························································································AY761120, KT308930, KT308931 B. medwediewii; KT308932, KT308933 B. megrelica····························································································

0.99 ·······································································································································KT308927 B. insignis········································································································································· 1 0.55 ············································································································KT308928 B. insignis; KT308929 B. insignis ssp. fansipanensis············································································································ ···························································································································KT308925, AY761100 B. alleghaniensis····························································································································· 0.86 0.88 ·····································································································································KT308926 B. murrayana······································································································································

····························································································································KT308911, KT308912 B. bomiensis·································································································································

0.72 1 1

0.59 0.62 0.81

0.60

0.74 0.87

·····························································································································AB243883, KT308905 B. globispica································································································································ ······································································································································KT308918 B. chinensis······································································································································· 0.60 0.94 ································································································KT308917 B. chinensis; KT308906 B. fargesii; KT308904, AY761111 B. globispica································································································ 0.80 0.60 ·····································································································································AY761110 B. glandulosa······································································································································· 0.73

0.85 ······························································································································AB243884, AB243885 B. corylifolia································································································································ 0.87 ·······················································································································AY761106, KT308907, KT308908 B. corylifolia························································································································ 0.89 ····························································································································AB243881, AB243882 B. chichibuensis···························································································································· 0.87

0.93 0.95

0.99 0.70 1

0.63 0.71 0.62 0.61

0.97 0.99

······························································································ KT308951 B. utilis var. jaquemonti; KT308952, KT308953, KT308961 B. ashburneri······························································································· ·························································································································KT308949 B. utilis; AY761099 B. albosinensis······················································································································ ·································································································································AY761134, KT308948 B. utilis···································································································································· ····························································································································FJ011780 B. utilis; FJ011773 B. davurica·························································································································· ·························································································································KT308923, KT308950 B. utilis var. occidentalis····················································································································· ···································································································································KT308955 B. utilis var. prattii··································································································································· ······························································································································KT308924, KT308954 B. albosinensis···························································································································· ·························································································································KT308947 B. albosinensis var. septentrionalis······················································································································ ········································································································································KT308956 B. ermanii········································································································································ 0.56 ········································································································································KT308957 B. ermanii········································································································································ ·························································································································AB243886 - AB243890, AY761108 B. ermanii······················································································································· ······························································································································AY761121, KT308978 B. michauxii································································································································ ·································································································································KT309024, KT309026 B. humilis································································································································ 0.64 ········································································································································KT309025 B. humilis········································································································································ ·································································································································AJ783643, AY761114 B. humilis································································································································· ·······································································································································AY761109 B. fruticosa········································································································································ ·································································································································AJ783642, AY761131 B. pumila································································································································· ·····················································································································KT309021 B. pumila; KT309019 B. nana ssp. exilis···················································································································· ··························································································································AY352336, AY761122, KT309020 B. nana·························································································································· ················································································································AB243897, AB243898, KT309022, KT309023 B. ovalifolia··············································································································· ····································································································································AB243895 B. middendorffii···································································································································· ·····································································································································KT309017 B. glandulosa······································································································································ ·······························································································································KT308966, AY761132 B. raddeana······························································································································· ·····································································································································KT308995 B. glandulosa······································································································································ 0.63 0.74 0.63 ······················································································································AY761126, KT309011, KT309012 B. papyrifera························································································································ ·····································································································································AF432067 B. papyrifera······································································································································· ····································································································································KT308988 B. × caerulea········································································································································ ·····································································································································KT309013 B. papyrifera······································································································································· ··························································································································KT309014 B. papyrifera var. commutata·····························································································································

····································································································KT309010 B. populifolia; KT308987 B. × caerulea; KT308985 B. × minor········································································································ ······················································································································AY761129, KT308994, KT309009 B. populifolia························································································································ ·······················································································································KT309027 B. occidentalis; KT309018 B. nana························································································································ ·····················································································································AY761125, DQ397523, KT309028 B. occidentalis······················································································································ ·····································································································································AY761128 B. platyphylla······································································································································· ·····································································································································KT308984 B. microphylla······································································································································ ·····································································································································KT308979 B. × utahensis····································································································································· ·······························································································································KT309015, KT309016 B. cordifolia······························································································································· ······································································································································KT308967 B. halophila········································································································································ ····································································································································KT308989 B. tianschanica····································································································································· ·······························································································································AM503889, AY761127 B. pendula································································································································ ···························································································································KT308991 B. pendula ssp. mandshurica···························································································································· ·················································AY352332 B. pendula; KT308992, KT308997 B. turkestanica; KT308993 B. obscura; KT308998, KT309001, KT309002, KT309007 B. pendula ssp. pendula············································ 0.51 ········································································································································FJ011777 B. pendula········································································································································ ········································································································AJ006445 B. pendula; MH178109 B. pendula f. obscura (B. kotulae)·········································································································· KT308968 B. browicziana; KT308990 B. pendula; KT308996, KT308999, KT309005, KT309008 B. pendula ssp. mandshurica; KT309003, KT309004 B. pendula ssp. szechuanica; KT309000, KT309006 B. pendula ssp. pendula ············································AB243891, FJ011778 B. platyphylla ; JN247411 B. pendula ; MH014809 B. borysthenica; MH178103 - MH178108, MH231207 B. pendula f. obscura (B. kotulae)·········································· ··························································································································KT308972 B. pubescens ssp. celtiberica···························································································································· ······································································KT308969, KT308970 B. pubescens var. pubescens; KT308977 B. pubescens ssp. celtiberica; KT308986 B. middendorffii···································································· ························································································································MH014808 B. borysthenica KW0064215 type························································································································ 0.61 ···································································································AJ783644 B. populifolia; MH178101, MH231212 B. pubescens f. sibakademica································································································ 0.64 ·····································································································································AY761130 B. pubescens······································································································································ ······KT308971, KT308983 B. pubescens ssp. litwinowii; KT308973, KT308976, KT308980 B. pubescens var. pumila; KT308974 KT308975 B. pubescens var. fragrans; KT308981, KT308982 B. pubescens ssp. pubescens······

Subgenus Betula Aspera Acuminata Nipponobetula

Section Costatae Apterocaryon Betula Dahuricae Asperae Lentae Acuminatae Nipponobetula

Subsection Asperae Chinenses B. apoiensis: B. × minor: B. ermannii (Costatae) × B. fruticosa (Apterocaryon ) B. glandulosa (Apterocaryon) × B. cordifolia (Betula)

Figure 9 Two maximum parsimony (MP) phylogenetic consensus trees based on 1000 iterations and generated using PAUP without (left) and with (right) consideration of ITS secondary structure recorded as 12-letter code. They are presented also separately in Figs. S3–S4, respectively. Clusters with bootstrap value less than 0.5 (50%) were collapsed. Subgenera, sections, and subsections are marked by different colors according to Ashburner & McAllister (2013). Identical taxa are connected by lines. Full-size DOI: 10.7717/peerj.10889/fig-9

Tarieiev et al. (2021), PeerJ, DOI 10.7717/peerj.10889 11/21 1 ·····································································································································AY352306 Alnus cordata········································································································································ 1 ······································································································································AY352329 Alnus viridis········································································································································· 0.633 ·······································································································································AB243892 B. grossa··········································································································································· 0.614 ·····················································································································KT308934, KT308935, AY761112, AY761113 B. grossa················································································································ 0.672 ··································································································AY761115, AY352330, KT308936, B. lenta; KT308937, KT308938 B. lenta f. uber································································································ 0.873 ······························································································AY761120, KT308930, KT308931 B. medwediewii; KT308932, KT308933 B. megrelica····························································································· 0.698 0.999 ·······································································································································KT308927 B. insignis·········································································································································· 1 0.542 ·············································································································KT308928 B. insignis; KT308929 B. insignis ssp. fansipanensis·············································································································

0.92 ·····························································································································KT308925, AY761100 B. alleghaniensis···························································································································· 0.972 0.762 ·····································································································································KT308926 B. murrayana········································································································································ 0.64

································································································································KT308911, KT308912 B. bomiensis······························································································································· 0.561 0.514 0.705 0.999 1

0.654 0.657 0.641 0.552 0.783

0.687 0.608 0.822 0.728 0.62 0.877

0.801 0.873

·······························································································································AB243883, KT308905 B. globispica······························································································································· ·······································································································································KT308918 B. chinensis········································································································································ 0.635 0.981 ·······························································································KT308917 B. chinensis; KT308906 B. fargesii; KT308904, AY761111 B. globispica·································································································· 0.802 ······································································································································AY761110 B. glandulosa······································································································································· 0.658 0.725

0.947 0.905 0.723

0.972 ·····························································································································AB243881, AB243882 B. chichibuensis····························································································································· 0.936 ·····················································································································AY761104, KT308915, KT308916 B. chichibuensis······················································································································ 0.901 0.938 ········································································································································KT308913 B. delavayi········································································································································ 0.992 ········································································································································KT308921 B. delavayi········································································································································ 0.791 0.99 0.711 0.501 ········································································································································AY761103 B. calcicola········································································································································ 0.527 0.702 0.735 ········································································································································KT308909 B. potaninii········································································································································ ·······································································································KT308914 B. calcicola; KT308910 B. potaninii; AY761107 B. delavayi··········································································································· 0.998 ·····························································································································AB243880, FJ011779 B. schmidtii···································································································································· 0.994 ······················································································································KT308919, KT308920, AY761133 B. schmidtii··························································································································· AY761107 B. delavayi·····················································································································

·

0.545

0.529 ·····································································································································AB243899 B. apoiensis········································································································································· ······················································································································AB243900 · AB243914, AY761102 B. apoiensis························································································································· ·····································································································································AB243915 B. apoiensis········································································································································· 0.58 ······························································································································KT309022, KT309023 B. ovalifolia·································································································································· ······························································································································AB243897, AB243898 B. ovalifolia·································································································································· ······························································································································AY761121, KT308978 B. michauxii·································································································································· ································································································································KT309024, KT309026 B. humilis·································································································································· 0.691 ·······································································································································KT309025 B. humilis··········································································································································· ·································································································································AJ783643, AY761114 B. humilis·································································································································· ·······································································································································AY761109 B. fruticosa········································································································································· ········································································································································AJ783642 B. pumila··········································································································································· ········································································································································AY761131 B. pumila··········································································································································· 0.59 ·················································································································KT309021 B. pumila; KT309019 B. nana ssp. exilis························································································································· ·······················································································································AY352336, AY761122, KT309020 B. nana······························································································································· ·································································································································AB243895 B. middendorffi········································································································································· ···························································································································KT308995, KT309017 B. glandulosa··································································································································

···························································································································· AB243894, FJ011770 B. davurica···································································································································· 0.623

0.521 ·····················································································································AB243893, FJ011771, FJ011772 B. dahurica····························································································································· 0.72 ······································································································································KT308962 B. dahurica·········································································································································· ···················································································································KT308963 B. dahurica; AY761105 B. chinensis·····························································································································

0.636 0.628

Subgenus Betula Aspera Acuminata Nipponobetula

Section Costatae Apterocaryon Betula Dahuricae Asperae Lentae Acuminatae Nipponobetula

Subsection Asperae Chinenses B. apoiensis: B. × minor: B. ermannii (Costatae) × B. fruticosa (Apterocaryon ) B. glandulosa (Apterocaryon) × B. cordifolia (Betula)

Figure 10 Two maximum likelihood (ML) phylogenetic consensus trees based on 1,000 iterations gen- erated without consideration of ITS secondary structure using MEGA X (left) and with consideration of ITS secondary structure recorded as 12-letter code using 4SALE and phangorn script in R package (right) . They are presented also separately in Figs. S5–S6, respectively. Clusters with bootstrap value less than 0.5 (50%) were collapsed. Subgenera, sections, and subsections are marked according to Ashburner & McAllister (2013). Identical taxa are connected by lines. Full-size DOI: 10.7717/peerj.10889/fig-10

Tarieiev et al. (2021), PeerJ, DOI 10.7717/peerj.10889 12/21 ················································································AY352306 Alnus cordata··················································································· ·······························································································································AY352329 Alnus viridis····································································································································· 0.987 ···································································AY761115, AY352330, KT308936 B. lenta; KT308937, KT308938 B. lenta f. uber······································································· 0.983 ····································································AY761120, KT308930, KT308931 B. medwediewii; KT308932, KT308933 B. megrelica········································································

0.882 1 ······························································KT308928 B. insignis; KT308929 B. insignis ssp. fansipanensis·································································· 1 0.847 ··························································································KT308927 B. insignis······························································································ ····································································································KT308926 B. murrayana······································································································· 0.999 1 ······························································································KT308925, AY761100 B. alleghaniensis·································································································

0.632 ·······························AY352331 B. nigra················································· 1 1 ····················KT308965 AY761124, AJ783646 B. nigra························· 0.851 ····································KT308964 B. nigra·······························

··································································································KT308911, KT308912 B. bomiensis·······································································································

······················································································AB243883, KT308905 B. globispica·························································································· 1 1 ·························································································KT308918 B. chinensis····························································································· ·······················································KT308917 B. chinensis; KT308906 B. fargesii; KT308904, AY761111 B. globispica···························································· 0.898 0.976 ···················································································AY761110 B. glandulosa·······················································································

0.873 0.997 ···································································AB243881 B. chichibuensis······································································ 0.999 0.96 ················································AB243882, KT308915, KT308916, AY761104 B. chichibuensis····················································· 1 1 ··························································KT308913 B. delavayi······························································ 1 ··············································KT308921 B. delavayi··················································· 1

0.993 0.521 ···················KT308914 B. calcicola; KT308910 B. potaninii; AY761107 B. delavayi·························· 0.501 0.999 ···················································KT308909 B. potaninii························································ ··············································AY761103 B. calcicola·················································· 1 ···········································································AB243880, FJ011779 B. schmidtii·················································································· 1 ···································································AY761133, KT308919, KT308920 B. schmidtii·········································································· ············································································································AY761121, KT308978 B. michauxii···················································································································

·····································································KT308948 B. utilis; KT308951 B. utilis var. jacquemontii; KT308961, KT308952, KT308953 B. ashburneri············································································ 0.828 ·············································································································KT308950 B. utilis var. occidentalis····················································································································· 0.851 ·······································································································AY761099 B. albosinensis; KT308949 B. utilis············································································································· 0.746 ················································································································FJ011780 B. utilis; ··································································································································AY761134 B. utilis········································································································································ ······································································································································KT308958 B. costata··········································································································································· ··········································································································································KT308955 B. utilis var. prattii················································································································································ ············································································································································KT308954 B. albosinensis················································································································································· ·······························································································································KT308947 B. albosinensis var. septentrionalis···································································································································· ···········································································································································AB243888 B. ermanii··············································································································································· ················································································································AB243886, AB243887, KT308956, KT308957, AY761108 B. ermanii······················································································································ 0.811 ·······································································································································AB243889 B. ermanii··········································································································································· 0.829

1 ~

···········································································································AB243890, AY761118, AY761119, KT308945, KT308946 B. maximowicziana················································································································ ···············································································································································AJ783641 B. alnoides·····································································································································

································································································································AB243894 B. davurica······················································································································· 0.882 ···············································································································AB243893, FJ011771, FJ011772, KT308962 B. dahurica······································································································ 0.976 ·····················································································································KT308963 B. dahurica; AY761105 B. chinensis············································································································· ······················································································································FJ011770 B. davurica············································································································· ·······································································································································AY761132, KT308966 B. raddeana····························································································································

~1 0.575 ····················································································································KT309021 B. pumila; KT309019 B. nana ssp. exilis·················································································································· ····································································································································AY761131 B. pumila································································································································· ·················································································································································AJ783642 B. pumila······································································································································ ··································································································································AY352336, AY761122, KT309020 B. nana························································································································ ·····························································································································KT309024, KT309026 B. humilis···················································································································· ································································································································AJ783643, AY761114, KT309025 B. humilis························································································································ ················································································································································AY761109 B. fruticosa···································································································································· ·············································································································································KT308995 B. glandulosa···································································································································· ········································································································································KT309017 B. glandulosa······························································································································ ············································································································································AB243895 B. middendorffii·································································································································· ···············································································································································AB243898 B. ovalifolia····································································································································

·············································································································································AB243897 B. ovalifolia·························································································································· 0.682 0.929 ········································································································AB243914 B. apoiensis···················································································· 0.98 ··············································································································AB243912 B. apoiensis·························································································· 0.837 ··············································································································KT309022, KT309023 B. ovailfolia··································································································· 0.965 ················································································································AB243915 B. apoiensis·····································································································

0.706 ······························································································································AB243909 B. apoiensis·················································································································· 0.728 0.774 ·······································································································AB243910 B. apoiensis········································································································ ················································································································AB243913 B. apoiensis··················································································································

····························································································································AB243911 B. apoiensis··················································································································· 0.51 ····································································································································AB243908 B. apoiensis·························································································································· ··········································································································································AY761102 B. apoiensis······························································································································

y y

Subgenus Betula Aspera Acuminata Nipponobetula

Section Costatae Apterocaryon Betula Dahuricae Asperae Lentae Acuminatae Nipponobetula

Subsection Asperae Chinenses B. apoiensis: B. × minor: B. ermannii (Costatae) × B. fruticosa (Apterocaryon ) B. glandulosa (Apterocaryon) × B. cordifolia (Betula)

Figure 11 Two Bayesian phylogenetic trees based on 100000000 iterations including 10% burnin and 4 by 4 model for all regions without consideration of ITS secondary structure (left) and doublet model for paired regions, but 4 by 4 with GTR +I for unpaired ones with consideration of ITS secondary struc- ture (right). They are presented also separately in Figs. S7–S8, respectively. Supported clades (> 0.95) are highlighted in red. Subgenera, sections, and subsections are marked according to Ashburner & McAllister (2013). Identical taxa are connected by lines. Full-size DOI: 10.7717/peerj.10889/fig-11

Tarieiev et al. (2021), PeerJ, DOI 10.7717/peerj.10889 13/21 of clusters and support values was greater for the trees with structure. However, despite the fact that consideration of the secondary structure increased the stability of clusters, many closely related taxa inside the genus were still not resolved. Even subgenera, sections, and subsections (marked in Fig. 11, Figs. S7–S11 according to Ashburner & McAllister (2013) using the same symbols as in Wang et al. (2016) could not be clearly distinguished since there are some cases where taxa from different sections and subsections have totally identical ITS sequences.

DISCUSSION Reconstructed secondary structure of ITS1 and ITS1 reveals a common pattern for the whole genus Betula L. The secondary structure for both of these spacers is more conservative than their sequences, and some common patterns are present across very large groups (Coleman, 2007). At the same time, secondary structure reconstruction reveals several types of nucleotide base changes, such as structural, non-structural, hemicompensatory (hCBC) and compensatory (CBC) changes, and, therefore, they can considerably enhance phylogenetic studies (Keller et al., 2010). The approach is described for a wide range of land plants, algae and animals (Coleman, 2009). In case of birch, structural, non-structural and hemicompensatory (hCBC) changes were detected, but compensatory base changes (CBCs) were not found. This was expected because CBCs are usually observed for fully reproductively isolated species (Coleman, 2000; Coleman, 2009; Müller et al., 2007), and, therefore, they could be used in delineation of biological species according to Mayr’s biological species concept (Mayr, 1942). Meanwhile, birches have very weak reproductive barriers with common hybridization events, and a number of taxa have hybrid origin (Ashburner & McAllister, 2013; Anamthawat-Jonsson & Thorsson, 2003; Anamthawat- Jonsson & Tomasson, 1999; Brown, Kennedy & Williams, 1982; Johnsson, 1945; Johnsson, 1949; Thorsson, Salmela & Anamthawat-Jonsson, 2001; Thorsson et al., 2010; Palmé et al., 2004; Thomson et al., 2015; Wang et al., 2014; Zohren et al., 2016). Therefore, the biological species concept is hardly applicable to many birch species. There are several other species concepts that could better accommodate this genus, but there is currently no strict consensus on this. Genetic clustering is a common approach to delineate birch species, but it could be spurious in case of species with high gene flow that can easily hybridize (Tsuda et al., 2017). Some birch taxa, even from different subgenera, sections, and subsections, have identical ITS sequences and, therefore, cannot be resolved using ITS markers. Since most of the sequences were retrieved from the GenBank, it is not always possible to verify the exact origin of the samples and correct species identification. The errors could occur due to morphological misidentification, and it is often impossible to check the accuracy of morphological identification of those specimens. However, identity of ITS sequences in species that have very distinctive morphology could be real, especially if it was confirmed by sequence data obtained by different independent research groups for the same species. It is likely that birch is still in the process of speciation, and reproductive barriers are not established yet or weak isolation is beneficial for survival in extreme environments.

Tarieiev et al. (2021), PeerJ, DOI 10.7717/peerj.10889 14/21 This reproductive feature is likely under positive selection, which could also explain the relatively low variation of ITS sequences within diverse groups. There are circa 15 known hybrids between different subgenera, sections, and subsections (Ashburner & McAllister, 2013). However, most of them were not included in this study. The exceptions are B. × minor (Tuck.) Fernald, which is a hybrid between B. cordifolia Regel (2n = 2x = 28, sect. Betula) and B. glandulosa Michx (2n = 2x = 28, sect. Apterocaryon), and B. apoiensis which is likely a hybrid between B. ermannii Cham. (2n = 4x = 56, sect. Costatae) and B. fruticosa Pallas (2n = 4x = 56, sect. Apterocaryon)(Ashburner & McAllister, 2013). This could explain a variety of different haplotypes of B. apoiensis with different clustering in the dendrograms. It should also be noted that those taxonomic categories such as subgenera, sections and subsections within birch genus were based only on morphological traits and do not necessarily reflect genetic differences. Multiple copies of ITS regions of different origin could be present in allopolyploids and interspecific hybrids, but only one of them could be amplified by PCR, which can also result in spurious phylogenetic trees. Recent genome-wide RAD-seq studies (Wang et al., 2020; Wang et al., 2021) solved some specific problems in birch taxonomy and phylogeny, especially clarifying parentage of some tetraploid species, enhancing resolution for diploid birch species with proposing changes on the level of sections and subsections and describing the new species B. buggsii that is morphologically similar to B. utilis ssp. albosinensis. However, some problems especially with polyploids with high ploidy levels remain. Therefore, we believe that further genome wide analyses can help to resolve the phylogenetic and systematic uncertainties in genus Betula even better.

CONCLUSIONS The general models of secondary structure of the ITS1 and ITS2 sequences across the whole genus Betula L. were presented for the first time in this study. Species-specific structural variants were detected for B. insignis and B. nigra. In combination with sequence variation it provides higher confidence in clustering, but only when the number of iterations is relatively low or moderate. With increasing the number of iterations the number of supported clusters and support values for the tree considering secondary structure remain more or less the same, while for the tree without considering secondary structure the number of supported clusters and support values tend to increase. In case of very closely related and hybrid taxa, ITS markers did not provide sufficient discrimination, as well as even some subsections and sections could not be clearly separated. In general, our data were in agreement with the relatively recently published phylogenetic study of the entire genus Betula L. based on the nucleotide variation of the ITS markers (Wang et al., 2016), but structural variation can additionally be used to distinguish or support some taxa.

ACKNOWLEDGEMENTS We thank Prof. Dr. Sergei L. Mosyakin and Dr. Nataliia Shiyan for herbarium type material of B. borysthenica, B. klokovii, and B. kotulae from KW herbarium, Institute of Botany, National Academy of Sciences of Ukraine; Prof. Dr. Dmitry V. Geltman and Dr. Maryna

Tarieiev et al. (2021), PeerJ, DOI 10.7717/peerj.10889 15/21 Legchenko for authentic material of B. pubescens var. sibakademica from LE herbarium, Botanical Institute of Russian Academy of Sciences (BIN RAS); Dr. Dmytro Iakushenko for the plant material of B. pubescens ssp. carpatica collected in Poland; Dr. Oksana Abduloieva for herbarium material from the (former) herbarium of Taras Shevchenko University of Kyiv; Dr. Thomas Pröschold for valuable advice and help with Bayesian inference analysis; Dr. Matthias Wolf for useful advice concerning ProfDist; Oleksandra Dolynska for technical assistance during experiments.

ADDITIONAL INFORMATION AND DECLARATIONS

Funding The research was funded by an EMA2ALRAKIS II internship, a Rufford Small Grant Foundation grant, and a DAAD (German Academic Exchange Service) PhD scholarship for Andrii Tarieiev. Support was also provided by the German Research Foundation (DFG) and the Open Access Publication Funds of the University of Göttingen. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript. Grant Disclosures The following grant information was disclosed by the authors: EMA2 ALRAKIS II internship, Rufford Small Grant Foundation grant. DAAD (German Academic Exchange Service). Competing Interests The authors declare there are no competing interests. Author Contributions • Andrii S. Tarieiev conceived and designed the experiments, performed the experiments, analyzed the data, prepared figures and/or tables, authored or reviewed drafts of the paper, and approved the final draft. • Oliver Gailing analyzed the data, authored or reviewed drafts of the paper, and approved the final draft. • Konstantin V. Krutovsky conceived and designed the experiments, analyzed the data, authored or reviewed drafts of the paper, and approved the final draft. Data Availability The following information was supplied regarding data availability: Raw data is available in the Supplemental Files. Supplemental Information Supplemental information for this article can be found online at http://dx.doi.org/10.7717/ peerj.10889#supplemental-information.

Tarieiev et al. (2021), PeerJ, DOI 10.7717/peerj.10889 16/21 REFERENCES Anamthawat-Jónsson K, Thórsson AT. 2003. Natural hybridization in birch: triploid hybrids between Betula nana and B. pubescens. Plant Cell Tissue and Organ Culture 75:99–107 DOI 10.1023/A:1025063123552. Anamthawat-Jónsson K, Tomasson T. 1999. High frequency of triploid birch hybrid by Betula nana seed parent. Hereditas 130:191–193. Ankenbrand MJ, Keller A, Wolf M, Schultz J, Förster F. 2015. ITS2 database V: Twice as much. Molecular Biology and Evolution 32(11):3030–3032 DOI 10.1093/molbev/msv174. Ashburner K, McAllister HA. 2013. The Genus Betula: a taxonomic revision of birches. Kew: Kew publishing, Royal Botanical Gardens. Baranov VI. 1924. Notes about birches from forest steppe of Western Siberia. Reports of Western-Siberian Department of Russian Geographical Society 4:47–54. Bona A, Petrova G, Jadwiszczak KA. 2018. Unfavourable habitat conditions can facilitate hybridisation between the endangered Betula humilis and its widespread relatives B. pendula and B. pubescens. Plant Ecology & Diversity 11(3):295–306. Brown IR, Kennedy D, Williams DA. 1982. The occurrence of natural hybrids between Betula pendula Roth and B. pubescens. Ehrh. Watsonia 14:133–145. Chen S, Wang Y, Yu L, Zheng T, Wang S, Yue Z, Jiang J, Kumari S, Zheng C, Tang H, Li J, Li Y, Chen J, Zhang W, Kuang H, Robertson JS, Zhao PX, Li H, Shu S, Yordanov YS, Huang H, Goodstein DM, Gai Y, Qi Q, Min JM, Xu CY, Wang SB, Qu G-Z, Paterson AH, Sankoff D, Wei H, Liu G, Yang C. 2021. Genome sequence and evolution of Betula platyphylla. Horticulture Research 8:37 DOI 10.1038/s41438-021-00481-7. Chen ZD, Manchester SR, Sun HY. 1999. Phylogeny and evolution of the Betulaceae as inferred from DNA sequence, morphology, and paleobotany. American Journal of Botany 86:1168–1181 DOI 10.2307/2656981. Coleman AW. 2000. The significance of a coincidence between evolutionary landmarks found in mating affinity and a DNA sequence. Protist 151:1–9 DOI 10.1078/1434-4610-00002. Coleman AW. 2003. ITS2 is a double-edged tool for eukaryote evolutionary compar- isons. Trends in Genetics 19:370–375 DOI 10.1016/S0168-9525(03)00118-5. Coleman AW. 2007. Pan-eukaryote ITS2 homologies revealed by RNA secondary structure. Nucleic Acids Research 35:3322–3329 DOI 10.1093/nar/gkm233. Coleman AW. 2009. Is there a molecular key to the level of biological species in eukaryotes? A DNA guide. Molecular Phylogenetics and Evolution 50:197–203 DOI 10.1016/j.ympev.2008.10.008. Darty K, Denise A, Ponty Y. 2009. VARNA: interactive drawing and editing of the RNA secondary structure. Bioinformatics 25(15):1974–1975 DOI 10.1093/bioinformatics/btp250. De Jong PC. 1993. An introduction to Betula: its morphology, evolution, classification and distribution, with a survey of recent work. In: Hunt D, ed. Proceedings of the

Tarieiev et al. (2021), PeerJ, DOI 10.7717/peerj.10889 17/21 IDS Betula symposium, 199. Richmond, UK: International Dendrology Society, Richmond, UK, 2–4. Domin K. 1927. Betula atrata Domin, nová bříza z Českomoravské vysočiny. Rozpravy České akademie věd a umění 36(52):1–10. Friedrich J, Dandekar T, Wolf M, Müller T. 2005. ProfDist: a tool for the construction of large phylogenetic trees based on profile distances. Bioinformatics 21(9):2108–2109 DOI 10.1093/bioinformatics/bti289. Hall TA. 1999. BioEdit: a user-friendly biological sequence alignment editor and analysis program for Windows 95/98/NT. Nucleic Acids Symposium Series 41:95–98. Huelsenbeck JP, Ronquist F. 2001. MRBAYES: Bayesian inference of phylogeny. Bioinformatics 17:754–755 DOI 10.1093/bioinformatics/17.8.754. Jadwiszczak KA, Vetchinnikova LV, Bona A, Tyburski Ł, Kuznetsova TJ, Isidorov VA. 2020. Analyses of molecular markers and leaf morphology of two rare birches, Betula obscura and B. pendula var carelica. Annals of Forest Research 63(2):121–113 DOI 10.15287/afr.2020.1973. Järvinen P, Palmé A, Orlando Morales L, Lännenpää M, Keinänen M, Sopanen T, Lascoux M. 2004. Phylogenetic relationships of Betula species (Betulaceae) based on nuclear ADH and chloroplast matK sequences. American Journal of Botany 91(11):1834–1845 DOI 10.3732/ajb.91.11.1834. Johnsson H. 1945. Interspecific hybridization within the genus Betula. Hereditas 31:163–176. Johnsson H. 1949. Studies on birch hybrids. Hereditas 35(2):115–135. Keller A, Förster F, Müller T, Dandekar T, Schultz J, Wolf M. 2010. Including RNA secondary structures improves accuracy and robustness in reconstruction of phylogenetic trees. Biology Direct 5:4 DOI 10.1186/1745-6150-5-4. Klokov MV. 1946. New materials to knowledge Ukrainian flora. Botanical Journal AS Ukrainian SSR 3(1–2):17–18. Koetschan C, Förster F, Keller A, Schleicher T, Ruderisch B, Schwarz R, Müller T, Wolf M, Schultz J. 2010. The ITS2 Database III–sequences and structures for phylogeny. Nucleic Acids Research 38:D275–D279 DOI 10.1093/nar/gkp966. Koetschan C, Hackl T, Müller T, Wolf M, Schultz JFörsterF. 2012. ITS2 database IV: interactive taxon sampling for internal transcribed spacer 2 based phylogenies. Molecular Phylogenetics and Evolution 63:585–588 DOI 10.1016/j.ympev.2012.01.026. Kuneš I, Linda R, Fér T, Karlík P, Baláš M, Ešnerova J, Vítámvás J, Bílý J, Urfus T. 2019. Is Betula carpatica genetically distinctive? A morphometric, cytometric and molecular study of birches in the Bohemian Massif with a focus on Carpathian birch. PLOS ONE 14(10):e0224387 DOI 10.1371/journal.pone.0224387. Kuzeneva IO. 1936. Betuleae Döll. In: Komarov VL, ed. Flora of the USSR. vol. 5. Mosqua-Leningrad: Editio Academiae Scientarum URSS 296. Li J, Shoupa S, Chen Z. 2005. Phylogenetics of Betula (Betulaceae) inferred from sequences of nuclear ribosomal DNA. Rhodora 107(929):69–86.

Tarieiev et al. (2021), PeerJ, DOI 10.7717/peerj.10889 18/21 Li J, Shoupa S, Chen Z. 2007. Phylogenetic relationships of diploid species of Betula (Betulaceae) inferred from DNA sequences of nuclear nitrate reductase. Systematic Botany 32(2):357–365 DOI 10.1600/036364407781179699. Linda R, Kunes˘ I, Balás˘ M. 2020. Genetic and morphometric variability between populations of Betula ÃŮoycoviensis from Poland and Czechia: A revised view of the taxonomic treatment of the Ojców birch. PLoS ONE. e0243310 DOI 10.1371/journal.pone.0243310. Mai JC, Coleman AW. 1997. The internal transcribed spacer 2 exhibits a common secondary structure in green algae and flowering plants. Journal of Molecular Evolution 44:258–271 DOI 10.1007/PL00006143. Mayr E. 1942. Systematics and the origin of species. New York: Columbia Univ. Press. Müller T, Philippi N, Dandekar T, Schultz J, Wolf M. 2007. Distinguishing species. RNA 13:1469–1472 DOI 10.1261/rna.617107. Olshanskyi IG, Shyian NM, Tarieiev AS. 2016. Types of Betula (Betulaceae) names de- scribed from Ukraine, kept in the National Herbarium of Ukraine (KW). Ukrainian Botanical Journal 73(4):401–403 DOI 10.15407/ukrbotj73.04.401. Palmé AE, Su Q, Palsson S, Lascoux M. 2004. Extensive sharing of chloroplast haplo- types among European birches indicates hybridization among Betula pendula, B. pubescens and B. nana. Molecular Ecology 13:167–178 DOI 10.1046/j.1365-294X.2003.02034.x. Regel E. 1865. Bemerkungen über die Gattungen Betula und Alnus nebst Beschreibung einiger neuer Arten. Bulletin de la Société impériale des naturalistes de Moscou 38(2):388–434. Ronquist F, Huelsenbeck JP. 2003. MRBAYES 3: Bayesian phylogenetic inference under mixed models. Bioinformatics 19:1572–1574 DOI 10.1093/bioinformatics/btg180. Ronquist F, Teslenko M, van der Mark P, Ayres DL, Darling A, Höhna S, Larget B, Liu L, Suchard MA, Huelsenbeck JP. 2012. MRBAYES 3.2: Efficient Bayesian phylogenetic inference and model selection across a large model space. Systematic Biology 61:539–542 DOI 10.1093/sysbio/sys029. Salojärvi J, Smolander OP, Nieminen K, Rajaraman S, Safronov O, Safdari P, Lammin- mäki A, Immanen J, Lan T, Tanskanen J, Rastas P, Amiryousefi A, Jayaprakash B, Kammonen JI, Hagqvist R, Eswaran G, Ahonen VH, Serra JA, Asiegbu FO, de Dios Barajas-Lopez J, Blande D, Blokhina O, Blomster T, Broholm S, Brosché M, Cui F, Dardick C, Ehonen SE, Elomaa P, Escamez S, Fagerstedt KV, Fujii H, Gauthier A, Gollan PJ, Halimaa P, Heino PI, Himanen K, Hollender C, Kangasjärvi S, Kauppinen L, Kelleher CT, Kontunen-Soppela S, Koskinen JP, Kovalchuk A, Kärenlampi SO, Kärkönen AK, Lim K-J, Leppälä J, Macpherson L, Mikola J, Mouhu K, Mähönen AP, Niinemets Ü, Oksanen E, Overmyer K, Palva ET, Pazouki L, Pennanen V, Puhakainen T, Poczai P, Possen BJHM, Punkkinen M, Rahikainen MM, Rousi M, Ruonala R, Van der Schoot C, Shapiguzov A, Sierla M, Sipilä TP, Sutela S, Teeri TH, Tervahauta AI, Vaattovaara A, Vahala J, Vetchinnikova L, Welling A, Wrzaczek M, Xu E, Paulin LG, Schulman AH, Lascoux M, Albert VA, Auvinen P, Helariutta Y, Kangasjärvi J. 2017. Genome sequencing and population

Tarieiev et al. (2021), PeerJ, DOI 10.7717/peerj.10889 19/21 genomic analyses provide insights into the adaptive landscape of silver birch. Nature Genetics 49:904–912 DOI 10.1038/ng.3862. Schenk MF, Thienpont C-N, Koopman WJM, Gilissen LJWJ, Smulders MJM. 2008. Phylogenetic relationships in Betula (Betulaceae) based on AFLP markers. Tree Genetics and Genomes 4(4):911–924 DOI 10.1007/s11295-008-0162-0. Schneider CK. 1915. Bemerkungen zur Systematik der Gattung Betula. Österreichische botanische Zeitschrift 65(10–12):305–312 DOI 10.1007/BF01660999. Schneider CK. 1916. Betulaceae. In: Sargent CS, ed. Plantae Wilsonianae. vol. 2. Cambridge: The University Press, 423–508. Seibel PN, Müller T, Dandekar T, Schultz J, Wolf M. 2006. 4SALE - A tool for syn- chronous RNA sequence and secondary structure alignment and editing. BMC Bioinformatics 7:498 DOI 10.1186/1471-2105-7-498. Seibel PN, Müller T, Dandekar T, Wolf M. 2008. Synchronous visual analysis and editing of RNA sequence and secondary structure alignments using 4SALE. BMC Research Notes 1:91 DOI 10.1186/1756-0500-1-91. Singevar K, Moschner CR, Hartung E, Fladung M. 2020. Species determination and phylogenetic relationships of the genus Betula inferred from multiple chloroplast and nuclear regions reveal the high methyl salicylate-producing ability of the ancestor. Trees 34:1131–1146 DOI 10.1007/s00468-020-01984-x. Skvortsov AK. 2002. A new system of the genus Betula L. –the birch. Bulletin of Moscow Society of Naturalists 107:73–76. Stöver BC, Müller KF. 2010. TreeGraph 2: combining and visualizing evidence from dif- ferent phylogenetic analyses. BMC Bioinformatics 11:7 DOI 10.1186/1471-2105-11-7. Swofford DL. 2003. PAUP*. Phylogenetic analysis using parsimony (*and other methods). Version 4. Sunderland: Sinauer Associates. Tarieiev AS, Olshanskyi I, Gailing O, Krutovsky KV. 2019. Taxonomy of dark- and white-barked birches related to Betula pendula Roth and B. pubescens Ehrh. (Betu- laceae) in Ukraine based on both morphological traits and DNA markers. Botanical Journal of the Linnean Society 191(1):142–154 DOI 10.1093/botlinnean/boz031. Thomson AM, Dick CW, Pascoini AL, Dayanandan S. 2015. Despite introgres- sive hybridization, North American birches (Betula spp.) maintain strong dif- ferentiation at nuclear microsatellite loci. Tree Genetics and Genomes 11:101 DOI 10.1007/s11295-015-0922-6. Thórsson ÆTh, Pálsson S, Lascoux M, Anamthawat-Jónsson K. 2010. Introgression and phylogeography of Betula nana (diploid), B. pubescens (tetraploid) and their triploid hybrids in Iceland inferred from cpDNA haplotype variation. Journal of Biogeography 37:2098–2110 DOI 10.1111/j.1365-2699.2010.02353.x. Thórsson ÆTh, Salmela E, Anamthawat-Jónsson K. 2001. Morphological, cytological, and molecular evidence for introgressive hybridisation in birch. Journal of Heredity 92:404–408 DOI 10.1093/jhered/92.5.404. Tsuda Y, Semerikov V, Sebastiani F, Vendramin GG, Lascoux M. 2017. Multispecies genetic structure and hybridization in the Betula genus across Eurasia. Molecular Ecology 26:589–605 DOI 10.1111/mec.13885.

Tarieiev et al. (2021), PeerJ, DOI 10.7717/peerj.10889 20/21 Tzvelev NN. 2002. De generibus Betula L. et Alnus Mill (Betulaceae) in Europa Orientali. Novitas Systematicae Plantarum Vascularum 34:47–73. Walters SM. 1993. Betula L. In: Tutin TG, Heywood VH, Burges NA, Moore DM, Valentine DH, Walters SM, Webb DA, eds. Flora Europaea 2nd edition Volume 1 Psilotaceae to Platanaceae. Cambridge University Press, 57–58. Wang N, Borrell JS, Bodles WJ, Kuttapitiya A, Nichols RA, Buggs RJ. 2014. Molecular footprints of the Holocene retreat of dwarf birch in Britain. Molecular Ecology 23:2771–2782 DOI 10.1111/mec.12768. Wang L, Ding J, Borrell JS, McAllister HA, Wang F, Liu L, Wang N. 2020. Molecular and morphological analyses clarify species delimitation and reveal a new Betula species in section Costatae. bioRxiv. DOI 10.1101/2020.10.29.361519. Wang N, Kelly LJ, McAllister HA, Zohren J, Buggs RJA. 2021. Resolving phylogeny and polyploid parentage using genus-wide genome-wide sequence data from birch trees. Molecular Phylogenetics and Evolution. 107126 DOI 10.1016/j.ympev.2021.107126. Wang N, McAllister HA, Bartlett PR, Buggs RG. 2016. Molecular phylogeny and genome size evolution of the genus Betula (Betulaceae). Annals of Botany 117(6):1023–1035 DOI 10.1093/aob/mcw048. Wang N, Thomson M, Bodles WJ, Crawford RM, Hunt HV, Featherstone AW, Pellicer J, Buggs RJ. 2013. Genome sequence of dwarf birch (Betula nana) and cross-species RAD markers. Molecular Ecology 22:3098–3111 DOI 10.1111/mec.12131. White TJ, Bruns T, Lee S, Taylor J. 1990. Amplification and direct sequencing of fungal ribosomal RNA genes for phylogenetics. In: Innis MA, Gelfand DH, Sninsky JJ, White TJ, eds. PCR Protocols: a guide to methods and applications. New York: Academic Press, 315–322. Wolf M, Chen S, Song J, Ankenbrand M, Müller T. 2013. Compensatory base changes in ITS2 secondary structures correlate with the biological species concept despite intragenomic variability in ITS2 sequences - a proof of concept. PLOS ONE 8(6):e66726 DOI 10.1371/journal.pone.0066726. Wolf M, Ruderisch B, Dandekar T, Müller T. 2008. ProfdistS: (Profile-) Distance based phylogeny on sequence-structure alignments. Bioinformatics 24:2401–2402 DOI 10.1093/bioinformatics/btn453. Zaverucha BV. 1964. New and rare birch species of Ukrainian flora. Ukrainian Botanical Journal 21(5):78–86. Zohren J, Wang N, Kardailsky I, Borrel JS, Joecker A, Nichols RA, Buggs RJA. 2016. Unidirectional diploidtetraploid introgression among British birch trees with shifting ranges shown by RAD markers. Molecular Ecology 25(11):2413–2426 DOI 10.1111/mec.13644. Zuker M. 2003. Mfold web server for nucleic acid folding and hybridization prediction. Nucleic Acids Research 31(13):3406–3415 DOI 10.1093/nar/gkg595.

Tarieiev et al. (2021), PeerJ, DOI 10.7717/peerj.10889 21/21