Aneura pinguis is a thalloid liverwort species with broad geographical distribution. It is composed of cryptic species, however, the number of cryptic species within A. pinguis is not known. Five cpDNA regions (matK, rbcL, rpoC1, trnH-psbA and trnL-trnF) and the entire nuclear ITS region were studied in 130 samples of A. pinguis from different geographical regions. The relationships between the cryptic species of A. pinguis, A. maxima and A. mirabilis were analyzed. All of the examined samples were clustered into 10 clades corresponding to 10 cryptic species of A. pinguis (marked A to J). Aneura mirabilis and A. maxima were nested among different cryptic species of A. pinguis, which indicates that A. pinguis is a paraphyletic taxon. Subgroups were found in cryptic species A, B, C and E. As single barcodes, all tested DNA regions had 100% discriminant power and fulfilled DNA barcode criteria for species identification; however, the only combination detected in all subgroups was trnL-trnF with trnH-psbA or ITS2. The distances between cryptic species were 11- to 35-fold higher than intraspecific distances. In all analyzed DNA regions, the distances between most pairs of cryptic A. pinguis species were higher than between A. maxima and A. mirabilis. All cryptic species of A. pinguis clearly differed in their habitat preferences, which suggests that habitat adaptation could be the main driving force behind cryptic speciation within this taxon.

Aneura pinguisis known as a species complex with several morphologically indiscernible species, which are often reproductively isolated from each other and show distinguishable genetic differences. Genetic dissimilarity of cryptic species may be detected by genomes comparison. This study presents the first complete sequences of chloroplast and mitochondrial genomes of six cryptic species ofA. pinguiscomplex:A. pinguisA, B, C, E, F, J. These genomes have been compared to each other in order to reconstruct phylogenetic relationships and to gain better understanding of the evolutionary process of cryptic speciation in this complex. The chloroplast genome with the nucleotide diversity 0.05111 and 1537 indels is by far more variable than mitogenome with π value 0.00233 and number of indels 1526. Tests of selection evidenced that on about 36% of chloroplast genes and on 10% of mitochondrial genes ofA. pinguisacts positive selection. It suggests an advanced speciation of species. The phylogenetic analyses based on genomes show thatA. pinguisis differentiated and forms three distinct clades. Moreover, on the cpDNA trees,Aneura mirabilisis nested among the cryptic species ofA. pinguis. This indicates that theA. pinguiscryptic species do not derive directly from one common ancestor.

DNA barcoding is a standard and efficient method, frequently used for identification, discrimination and discovery of new species. Although this approach is very useful for classifying the world's biodiversity, little is known about its usefulness in barcoding at lower taxonomic level and its discrimination rate for closely related species, like conifers. In this study, we compared the genetic variation of eight chloroplast DNA barcode regions ( matK, rbcL, trnH-psbA, trnL-trnF, rpl20-rps18, trnV, ycf1, ycf2 ) in 17 conifers - three closely related pines from Pinus mugo complex and 14 more distant conifers representing two genera and four sections of the Pinaceae family. The discrimination rate for a single and for multiple DNA barcode regions analyzed in this study was estimated using the Tree-Building and PWG-Distance methods. The usefulness of the DNA barcoding approach for analyzing and resolving taxonomic inconsistency among closely related and more phylogenetically distant conifers was evaluated and discussed.