Aneura pinguis is 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 of A. pinguis complex: A. pinguis A, 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 of A. pinguis acts positive selection. It suggests an advanced speciation of species. The phylogenetic analyses based on genomes show that A. pinguis is differentiated and forms three distinct clades. Moreover, on the cpDNA trees, Aneura mirabilis is nested among the cryptic species of A. pinguis. This indicates that the A. pinguis cryptic species do not derive directly from one common ancestor.

Background: Molecular research revealed that some of the European Calypogeia species described on the basis of morphological criteria are genetically heterogeneous and, in fact, are species complexes. DNA barcoding is already commonly used for correct identification of difficult to determine species, to disclose cryptic species, or detecting new taxa. Among liverworts, some DNA fragments, recommend as universal plant DNA barcodes, cause problems in amplification. Super-barcoding based on genomic data, makes new opportunities in a species identification. Results: On the basis of 22 individuals, representing 10 Calypogeia species, plastid genome was tested as a superbarcode. It is not effective in 100%, nonetheless its success of species discrimination (95.45%) is still conspicuous. It is not excluded that the above outcome may have been upset by cryptic speciation in C. suecica, as our results indicate. Having the sequences of entire plastomes of European Calypogeia species, we also discovered that the ndhB and ndhH genes and the trnT-trnL spacer identify species in 100%. Conclusions: This study shows that even if a super-barcoding is not effective in 100%, this method does not close the door to a traditional single-or multi-locus barcoding. Moreover, it avoids many complication resulting from the need to amplify selected DNA fragments. It seems that a good solution for species discrimination is a development of so-called "specific barcodes" for a given taxonomic group, based on plastome data.

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.

Liverwort mitogenomes are considered to be evolutionarily stable. A comparative analysis of four Calypogeia species revealed differences compared to previously sequenced liverwort mitogenomes. Such differences involve unexpected structural changes in the two genes, cox1 and atp1, which have lost three and two introns, respectively. The group I introns in the cox1 gene are proposed to have been lost by two-step localized retroprocessing, whereas one-step retroprocessing could be responsible for the disappearance of the group II introns in the atp1 gene. These cases represent the first identified losses of introns in mitogenomes of leafy liverworts (Jungermanniopsida) contrasting the stability of mitochondrial gene order with certain changes in the gene content and intron set in liverworts.

Oil bodies are the unique feature of most liverworts. Their shape, color and distribution pattern in leaf and underleaf cells are important taxonomic features of the genus Calypogeia. Most species of the genus Calypogeia have pellucid and colorless oil bodies, whereas colored, including gray to pale brown, purple-brown or blue oil bodies, are rare. To date, C. azurea was the only species with blue oil bodies to have been considered as a species of the Holarctic range. This species has been noted in various parts of the northern hemisphere–from North America, through Europe to the Far East. The aim of this study was to determine the genetic diversity of C. azurea from different parts of its distribution range and to ascertain whether blue oil bodies appeared once or several times in the evolution of the genus Calypogeia. The phylogenetic analyses based on four plastid regions (rbcL, trnG, trnL, trnH-psbA) and one nuclear region (ITS2) revealed that C. azurea is presently a paraphyletic taxon, with other Calypogeia species nested among C. azurea accessions that were clustered into four different clades. Based on the level of genetic divergence (1.03–2.17%) and the observed morphological, ecological and geographical differences, the evaluated clades could be regarded as previously unrecognized species. Four species were identified: C. azurea Stotler & Crotz (a European species corresponding to the holotype), two new species from Pacific Asia—C. orientalis Buczkowska & Bakalin and C. sinensis Bakalin & Buczkowska, and a North American species which, due to the lack of identifiable morphological features, must be regarded as the cryptic species of C. azurea with a provisional name of C. azurea species NA.