Comparative Analyses of Complete Chloroplast Genomes and Karyotypes of Allotetraploid Iris koreana and Its Putative Diploid Parental Species (Iris Series Chinenses, Iridaceae)

Complete plastid genome of Iris orchioides and comparative analysis with 19 Iris plastomes

Iris is a cosmopolitan genus comprising approximately 280 species distributed throughout the Northern Hemisphere. Although Iris is the most diverse group in the Iridaceae, the number of taxa is debatable owing to various taxonomic issues. Plastid genomes have been widely used for phylogenetic research in plants; however, only limited number of plastid DNA markers are available for phylogenetic study of the Iris. To understand the genomic features of plastids within the genus, including its structural and genetic variation, we newly sequenced and analyzed the complete plastid genome of I. orchioides and compared it with those of 19 other Iris taxa. Potential plastid markers for phylogenetic research were identified by computing the sequence divergence and phylogenetic informativeness. We then tested the utility of the markers with the phylogenies inferred from the markers and whole-plastome data. The average size of the plastid genome was 152,926 bp, and the overall genomic content and organization were nearly identical among the 20 Iris taxa, except for minor variations in the inverted repeats. We identified 10 highly informative regions (matK, ndhF, rpoC2, ycf1, ycf2, rps15-ycf, rpoB-trnC, petA-psbJ, ndhG-ndhI and psbK-trnQ) and inferred a phylogeny from each region individually, as well as from their concatenated data. Remarkably, the phylogeny reconstructed from the concatenated data comprising three selected regions (rpoC2, ycf1 and ycf2) exhibited the highest congruence with the phylogeny derived from the entire plastome dataset. The result suggests that this subset of data could serve as a viable alternative to the complete plastome data, especially for molecular diagnoses among closely related Iris taxa, and at a lower cost.

Comparative microscopic investigations of leaf epidermis in four Ajuga species from Korea

The genus Ajuga is widely distributed in temperate to subtropical regions, and four species are currently recognized in Korea (A. decumbens, A. multiflora, A. nipponensis, and A. spectabilis), but epidermal anatomical differences across these species have never been described. A comparative study of the leaf micromorphological characteristics of Korean Ajuga species was performed using light microscopy (LM) and scanning electron microscopy (SEM) to elucidate their taxonomic usefulness and to assess leaf micromorphological diversity. Considerable diversity in epidermal and stomatal anatomy was observed across Korean Ajuga species. Species had both hypostomatic or amphistomatic leaves, with anomocytic, anisocytic, diactyic, or actinocytic stomatal complexes. Guard cell length across species ranged from 17.66 ± 0.57 μm to 32.50 ± 2.38 μm and correlated with genome size. Abnormal stomata were frequently observed in three species (A. decumbens, A. multiflora, and A. nipponensis) but not in A. spectabilis. Three types of glandular trichomes were found: peltate in all species, short-stalked in all species, and long-stalked glandular trichomes in A. multiflora. Among the investigated leaf micromophological characters, trichome type, epidermal cell shape, and stomatal morphology were all taxonomically informative traits at a species level. RESEARCH HIGHLIGHTS: A comprehensive micromorphological description of the leaf surface is provided for Korean Ajuga species using scanning electron microscopic (SEM) and light microscopic (LM) analyses. The diverse range of stomatal development and the occurrence of polymorphic stomatal types are documented for the first time in Korean Ajuga species. The great diversity in stomatal and trichome morphology in Korean Ajuga species are taxonomically useful traits for species identification.

Comparative Analysis of Chloroplast Genomes of "Tiantai Wu-Yao" (Lindera aggregata) and Taxa of the Same Genus and Different Genera

Lindera aggregata is a species of the Lauraceae family, which has important medicinal, economic and ornamental values. In this study, we sequenced, assembled and annotated the chloroplast genome of L. aggregata and reannotated and corrected eight unverified annotations in the same genus. The chloroplast genomes taxa from Lindera and from different genera of Lauraceae were compared and analyzed, and their phylogenetic relationship and divergence time were speculated. All the 36 chloroplast genomes had typical quadripartite structures that ranged from 150,749 to 154,736 bp in total length. These genomes encoded 111-112 unique genes, including 78-79 protein-coding genes, 29-30 tRNA and 4 rRNA. Furthermore, there were 78-97 SSRs loci in these genomes, in which mononucleotide repeats were the most abundant; there were 24-49 interspersed repeats, and forward repeat types were the most frequent. The codon bias patterns of all species tended to use codons ending with A or U. Five and six highly variable regions were identified within genus and between genera, respectively, and three common regions (ycf1, ndhF-rpl32 and rpl32-trnL) were identified, which can be used as important DNA markers for phylogeny and species identification. According to the evaluation of the Ka/Ks ratio, most of the genes were under purifying selection, and only 10 genes were under positive selection. Finally, through the construction of the evolutionary tree of 39 chloroplast genomes, the phylogenetic relationship of Lauraceae was clarified and the evolutionary relationship of Lindera was revealed. The species of genus Lindera experienced rapid adaptive radiation from Miocene to Pleistocene. The results provided valuable insights for the study of chloroplast genomes in the Lauraceae family, especially in the genus Lindera.

Complete Chloroplast Genomes and the Phylogenetic Analysis of Three Native Species of Paeoniaceae from the Sino-Himalayan Flora Subkingdom

Paeonia delavayi var. lutea, Paeonia delavayi var. angustiloba, and Paeonia ludlowii are Chinese endemics that belong to the Paeoniaceae family and have vital medicinal and ornamental value. It is often difficult to classify Paeoniaceae plants based on their morphological characteristics, and the limited genomic information has strongly hindered molecular evolution and phylogenetic studies of Paeoniaceae. In this study, we sequenced, assembled, and annotated the chloroplast genomes of P. delavayi var. lutea, P. delavayi var. angustiloba, and P. ludlowii. The chloroplast genomes of these strains were comparatively analyzed, and their phylogenetic relationships and divergence times were inferred. These three chloroplast genomes exhibited a typical quadripartite structure and were 152,687-152,759 bp in length. Each genome contains 126-132 genes, including 81-87 protein-coding genes, 37 transfer RNAs, and 8 ribosomal RNAs. In addition, the genomes had 61-64 SSRs, with mononucleotide repeats being the most abundant. The codon bias patterns of the three species tend to use codons ending in A/U. Six regions of high variability were identified (psbK-psbL, trnG-UCC, petN-psbM, psbC, rps8-rpl14, and ycf1) that can be used as DNA molecular markers for phylogenetic and taxonomic analysis. The Ka/Ks ratio indicates positive selection for the rps18 gene associated with self-replication. The phylogenetic analysis of 99 chloroplast genomes from Saxifragales clarified the phylogenetic relationships of Paeoniaceae and revealed that P. delavayi var. lutea, P. delavayi var. angustiloba, and P. ludlowii are monophyletic groups and sisters to P. delavayi. Divergence time estimation revealed two evolutionary divergences of Paeoniaceae species in the early Oligocene and Miocene. Afterward, they underwent rapid adaptive radiation from the Pliocene to the early Pleistocene when P. delavayi var. lutea, P. delavayi var. angustiloba, and P. ludlowii formed. The results of this study enrich the chloroplast genomic information of Paeoniaceae and reveal new insights into the phylogeny of Paeoniaceae.

Does the evolution of micromorphology accompany chromosomal changes on dysploid and polyploid levels in the Barnardia japonica complex (Hyacinthaceae)?

BACKGROUND

Chromosome number and genome size changes via dysploidy and polyploidy accompany plant diversification and speciation. Such changes often impact also morphological characters. An excellent system to address the questions of how extensive and structured chromosomal changes within one species complex affect the phenotype is the monocot species complex of Barnardia japonica. This taxon contains two well established and distinct diploid cytotypes differing in base chromosome numbers (AA: x = 8, BB: x = 9) and their allopolyploid derivatives on several ploidy levels (from 3x to 6x). This extensive and structured genomic variation, however, is not mirrored by gross morphological differentiation.

RESULTS

The current study aims to analyze the correlations between the changes of chromosome numbers and genome sizes with palynological and leaf micromorphological characters in diploids and selected allopolyploids of the B. japonica complex. The chromosome numbers varied from 2n = 16 and 18 (2n = 25 with the presence of supernumerary B chromosomes), and from 2n = 26 to 51 in polyploids on four different ploidy levels (3x, 4x, 5x, and 6x). Despite additive chromosome numbers compared to diploid parental cytotypes, all polyploid cytotypes have experienced genome downsizing. Analyses of leaf micromorphological characters did not reveal any diagnostic traits that could be specifically assigned to individual cytotypes. The variation of pollen grain sizes correlated positively with ploidy levels.

CONCLUSIONS

This study clearly demonstrates that karyotype and genome size differentiation does not have to be correlated with morphological differentiation of cytotypes.

Comparative Genomic Analysis Uncovers the Chloroplast Genome Variation and Phylogenetic Relationships of Camellia Species

Camellia is the largest genus in the family Theaceae. Due to phenotypic diversity, frequent hybridization, and polyploidization, an understanding of the phylogenetic relationships between Camellia species remains challenging. Comparative chloroplast (cp) genomics provides an informative resource for phylogenetic analyses of Camellia. In this study, 12 chloroplast genome sequences from nine Camellia species were determined using Illumina sequencing technology via de novo assembly. The cp genome sizes ranged from 156,545 to 157,021 bp and were organized into quadripartite regions with the typical angiosperm cp genomes. Each genome harbored 87 protein-coding, 37 transfer RNA, and 8 ribosomal RNA genes in the same order and orientation. Differences in long and short sequence repeats, SNPs, and InDels were detected across the 12 cp genomes. Combining with the complete cp sequences of seven other species in the genus Camellia, a total of nine intergenic sequence divergent hotspots and 14 protein-coding genes with high sequence polymorphism were identified. These hotspots, especially the InDel (~400 bp) located in atpH-atpI region, had sufficient potential to be used as barcode markers for further phylogenetic analysis and species identification. Principal component and phylogenetic analysis suggested that regional constraints, rather than functional constraints, strongly affected the sequence evolution of the cp genomes in this study. These cp genomes could facilitate the development of new molecular markers, accurate species identification, and investigations of the phylogenomic relationships of the genus Camellia.