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 with each other. However, the resolving power of these markers is still insufficient to reliably discriminate some closely related species. To achieve this aim more reliably there is a need for application of modern genomic approaches in combination with traditional ones.