The Plastome Sequences of Triticum sphaerococcum (ABD) and Triticum turgidum subsp. durum (AB) Exhibit Evolutionary Changes,Structural Characterization,Comparative Analysis,Phylogenomics and Time Divergence

The mechanism and course of Triticum plastome evolution is currently unknown; thus, it remains unclear how Triticum plastomes evolved during recent polyploidization. Here, we report the complete plastomes of two polyploid wheat species, Triticum sphaerococcum (AABBDD) and Triticum turgidum subsp. durum (AABB), and compare them with 19 available and complete Triticum plastomes to create the first map of genomic structural variation. Both T. sphaerococcum and T. turgidum subsp. durum plastomes were found to have a quadripartite structure, with plastome lengths of 134,531 bp and 134,015 bp, respectively. Furthermore, diploid (AA), tetraploid (AB, AG) and hexaploid (ABD, AGA^m) Triticum species plastomes displayed a conserved gene content and commonly harbored an identical set of annotated unique genes. Overall, there was a positive correlation between the number of repeats and plastome size. In all plastomes, the number of tandem repeats was higher than the number of palindromic and forward repeats. We constructed a Triticum phylogeny based on the complete plastomes and 42 shared genes from 71 plastomes. We estimated the divergence of Hordeum vulgare from wheat around 11.04–11.9 million years ago (mya) using a well-resolved plastome tree. Similarly, Sitopsis species diverged 2.8–2.9 mya before Triticum urartu (AA) and Triticum monococcum (AA). Aegilops speltoides was shown to be the maternal donor of polyploid wheat genomes and diverged ~0.2–0.9 mya. The phylogeny and divergence time estimates presented here can act as a reference framework for future studies of Triticum evolution.     

Plastid Genomes of the Early Vascular Plant Genus Selaginella Have Unusual Direct Repeat Structures and Drastically Reduced Gene Numbers

The early vascular plants in the genus Selaginella, which is the sole genus of the Selaginellaceae family, have an important place in evolutionary history, along with ferns, as such plants are valuable resources for deciphering plant evolution. In this study, we sequenced and assembled the plastid genome (plastome) sequences of two Selaginella tamariscina individuals, as well as Selaginella stauntoniana and Selaginella involvens. Unlike the inverted repeat (IR) structures typically found in plant plastomes, Selaginella species had direct repeat (DR) structures, which were confirmed by Oxford Nanopore long-read sequence assembly. Comparative analyses of 19 lycophytes, including two Huperzia and one Isoetes species, revealed unique phylogenetic relationships between Selaginella species and related lycophytes, reflected by structural rearrangements involving two rounds of large inversions that resulted in dynamic changes between IR and DR blocks in the plastome sequence. Furthermore, we present other uncommon characteristics, including a small genome size, drastic reductions in gene and intron numbers, a high GC content, and extensive RNA editing. Although the 16 Selaginella species examined may not fully represent the genus, our findings suggest that Selaginella plastomes have undergone unique evolutionary events yielding genomic features unparalleled in other lycophytes, ferns, or seed plants.     

Structural Mutations in the Organellar Genomes of Valeriana sambucifolia f. dageletiana (Nakai. ex Maekawa) Hara Show Dynamic Gene Transfer

Valeriana sambucifolia f. dageletiana (Nakai. ex Maekawa) Hara is a broad-leaved valerian endemic to Ulleung Island, a noted hot spot of endemism in Korea. However, despite its widespread pharmacological use, this plant remains comparatively understudied. Plant cells generally contain two types of organellar genomes (the plastome and the mitogenome) that have undergone independent evolution, which accordingly can provide valuable information for elucidating the phylogenetic relationships and evolutionary histories of terrestrial plants. Moreover, the extensive mega-data available for plant genomes, particularly those of plastomes, can enable researchers to gain an in-depth understanding of the transfer of genes between different types of genomes. In this study, we analyzed two organellar genomes (the 155,179 bp plastome and the 1,187,459 bp mitogenome) of V. sambucifolia f. dageletiana and detected extensive changes throughout the plastome sequence, including rapid structural mutations associated with inverted repeat (IR) contraction and genetic variation. We also described features characterizing the first reported mitogenome sequence obtained for a plant in the order Dipsacales and confirmed frequent gene transfer in this mitogenome. We identified eight non-plastome-originated regions (NPRs) distributed within the plastome of this endemic plant, for six of which there were no corresponding sequences in the current nucleotide sequence databases. Indeed, one of these unidentified NPRs unexpectedly showed certain similarities to sequences from bony fish. Although this is ostensibly difficult to explain, we suggest that this surprising association may conceivably reflect the occurrence of gene transfer from a bony fish to the plastome of an ancestor of V. sambucifolia f. dageletiana mediated by either fungi or bacteria.     

The Barley Chloroplast Mutator (cpm) Mutant: All Roads Lead to the Msh1 Gene

The barley chloroplast mutator (cpm) is a nuclear gene mutant that induces a wide spectrum of cytoplasmically inherited chlorophyll deficiencies. Plastome instability of cpm seedlings was determined by identification of a particular landscape of polymorphisms that suggests failures in a plastome mismatch repair (MMR) protein. In Arabidopsis, MSH genes encode proteins that are in charge of mismatch repair and have anti-recombination activity. In this work, barley homologs of these genes were identified, and their sequences were analyzed in control and cpm mutant seedlings. A substitution, leading to a premature stop codon and a truncated MSH1 protein, was identified in the Msh1 gene of cpm plants. The relationship between this mutation and the presence of chlorophyll deficiencies was established in progenies from crosses and backcrosses. These results strongly suggest that the mutation identified in the Msh1 gene of the cpm mutant is responsible for the observed plastome instabilities. Interestingly, comparison of mutant phenotypes and molecular changes induced by the barley cpm mutant with those of Arabidopsis MSH1 mutants revealed marked differences.     

Genomic Resources of Three Pulsatilla Species Reveal Evolutionary Hotspots,Species-Specific Sites and Variable Plastid Structure in the Family Ranunculaceae

Background: The European continent is presently colonized by nine species of the genus Pulsatilla, five of which are encountered only in mountainous regions of southwest and south-central Europe. The remaining four species inhabit lowlands in the north-central and eastern parts of the continent. Most plants of the genus Pulsatilla are rare and endangered, which is why most research efforts focused on their biology, ecology and hybridization. The objective of this study was to develop genomic resources, including complete plastid genomes and nuclear rRNA clusters, for three sympatric Pulsatilla species that are most commonly found in Central Europe. The results will supply valuable information about genetic variation, which can be used in the process of designing primers for population studies and conservation genetics research. The complete plastid genomes together with the nuclear rRNA cluster can serve as a useful tool in hybridization studies. Methodology/principal findings: Six complete plastid genomes and nuclear rRNA clusters were sequenced from three species of Pulsatilla using the Illumina sequencing technology. Four junctions between single copy regions and inverted repeats and junctions between the identified locally-collinear blocks (LCB) were confirmed by Sanger sequencing. Pulsatilla genomes of 120 unique genes had a total length of approximately 161–162 kb, and 21 were duplicated in the inverted repeats (IR) region. Comparative plastid genomes of newly-sequenced Pulsatilla and the previously-identified plastomes of Aconitum and Ranunculus species belonging to the family Ranunculaceae revealed several variations in the structure of the genome, but the gene content remained constant. The nuclear rRNA cluster (18S-ITS1-5.8S-ITS2-26S) of studied Pulsatilla species is 5795 bp long. Among five analyzed regions of the rRNA cluster, only Internal Transcribed Spacer 2 (ITS2) enabled the molecular delimitation of closely-related Pulsatilla patens and Pulsatilla vernalis. Conclusions/significance: The determination of complete plastid genome and nuclear rRNA cluster sequences in three species of the genus Pulsatilla is an important contribution to our knowledge of the evolution and phylogeography of those endangered taxa. The resulting data can be used to identify regions that are particularly useful for barcoding, phylogenetic and phylogeographic studies. The investigated taxa can be identified at each stage of development based on their species-specific SNPs. The nuclear and plastid genomic resources enable advanced studies on hybridization, including identification of parent species, including their roles in that process. The identified nonsynonymous mutations could play an important role in adaptations to changing environments. The results of the study will also provide valuable information about the evolution of the plastome structure in the family Ranunculaceae.     

Homogeneous Nature of Malaysian Marine Fish Epinephelus fuscoguttatus (Perciformes; Serranidae): Evidence Based on Molecular Markers,Morphology and Fourier Transform Infrared Analysis

Taxonomic confusion exists within the genus Epinephelus due to the lack of morphological specializations and the overwhelming number of species reported in several studies. The homogenous nature of the morphology has created confusion in the Malaysian Marine fish species Epinephelus fuscoguttatus and Epinephelus hexagonatus. In this study, the partial DNA sequence of the 16S gene and mitochondrial nucleotide sequences of two gene regions, Cytochrome Oxidase Subunit I and III were used to investigate the phylogenetic relationship between them. In the phylogenetic trees, E. fuscoguttatus was monophyletic with E. hexagonatus species and morphology examination shows that no significant differences were found in the morphometric features between these two taxa. This suggests that E. fuscoguttatus is not distinguishable from E. hexagonatus species, and that E. fuscoguttatus have been identified to be E. hexagonatus species is likely attributed to differences in environment and ability to camouflage themselves under certain conditions. Interestingly, this finding was also supported by Principal Component Analysis on Attenuated Total Reflectance–Fourier-transform Infrared (ATR-FTIR) data analysis. Molecular, morphological and meristic characteristics were combined with ATR-FTIR analysis used in this study offer new perspectives in fish species identification. To our knowledge, this is the first report of an extensive genetic population study of E. fuscoguttatus in Malaysia and this understanding will play an important role in informing genetic stock-specific strategies for the management and conservation of this highly valued fish.     

Ovule Development and in Planta Transformation of Paphiopedilum Maudiae by Agrobacterium-Mediated Ovary-Injection

In this paper, the development of the Paphiopedilum Maudiae embryo sac at different developmental stages after pollination was assessed by confocal laser scanning microscopy. The mature seeds of P. Maudiae consisted of an exopleura and a spherical embryo, but without an endosperm, while the inner integument cells were absorbed by the developing embryo. The P. Maudiae embryo sac exhibited an Allium type of development. The time taken for the embryo to develop to a mature sac was 45-50 days after pollination (DAP) and most mature embryo sacs had completed fertilization and formed zygotes by about 50–54 DAP. In planta transformation was achieved by injection of the ovaries by Agrobacterium, resulting in 38 protocorms or seedlings after several rounds of hygromycin selection, corresponding to 2, 7, 5, 1, 3, 4, 9, and 7 plantlets from Agrobacterium-mediated ovary-injection at 30, 35, 42, 43, 45, 48, 50, and 53 DAP, respectively. Transformation efficiency was highest at 50 DAP (2.54%), followed by 2.48% at 53 DAP and 2.45% at 48 DAP. Four randomly selected hygromycin-resistant plants were GUS-positive after PCR analysis. Semi-quantitative PCR and quantitative real-time PCR analysis revealed the expression of the hpt gene in the leaves of eight hygromycin-resistant seedlings following Agrobacterium-mediated ovary-injection at 30, 35, 42, 43, 45, 48, 50, and 53 DAP, while hpt expression was not detected in the control. The best time to inject P. Maudiae ovaries in planta with Agrobacterium is 48-53 DAP, which corresponds to the period of fertilization. This protocol represents the first genetic transformation protocol for any Paphiopedilum species and will allow for expanded molecular breeding programs to introduce useful and interesting genes that can expand its ornamental and horticulturally important characteristics.     

Morphological and Genetic Diversities of Habenaria radiata (Orchidaceae) in the Kinki Area,Japan

Floral organs have evolved from leaves for reproduction, and the morphological analyses help to understand the plant diversity and evolution. Habenaria radiata (syn. Pecteilis radiata) is a terrestrial orchid living in wetlands in Japan, Russia, South Korea, and China. The habitats of this plant in Japan have been reduced because of environmental destruction and overexploitation, and thus it is on the Red List of Japan as a Near Threatened species. One of the three petals of the H. radiata flower is called a lip or labellum, which resembles a flying white bird, egret, or white heron, with its proposed function being to attract pollinators. To understand the diversity of H. radiata plants in different areas, we examined the lip morphology and phylogeny of populations from eight habitats in the Kinki area, Japan. The complex shapes of the lips were quantified and presented as a radar chart, enabling characterization of the morphological difference among populations. Phylogenetic analysis with microsatellite markers that we generated showed the variation of genetic diversity among populations, suggesting the different degrees of inbreeding, outbreeding, and vegetative propagation. Our approach offers a basic method to characterize the morphological and genetic diversity in natural populations.     

Lineage-Specific Variation in IR Boundary Shift Events,Inversions, and Substitution Rates among Caprifoliaceae s.l. (Dipsacales) Plastomes

Caprifoliaceae s.l. plastid genomes (plastomes) show that one inversion and two inverted repeat boundary shifts occurred in the common ancestor of this family, after which the plastomes are generally conserved. This study reports plastome sequences of five additional species, Fedia cornucopiae, Valeriana fauriei, and Valerianella locusta from the subfamily Valerianoideae, as well as Dipsacus japonicus and Scabiosa comosa from the subfamily Dipsacoideae. Combined with the published plastomes, these plastomes provide new insights into the structural evolution of plastomes within the family. Moreover, the three plastomes from the subfamily Valerianoideae exhibited accelerated nucleotide substitution rates, particularly at synonymous sites, across the family. The patterns of accD sequence divergence in the family are dynamic with structural changes, including interruption of the conserved domain and increases in nonsynonymous substitution rates. In particular, the Valeriana accD gene harbors a large insertion of amino acid repeat (AAR) motifs, and intraspecific polymorphism with a variable number of AARs in the Valeriana accD gene was detected. We found a correlation between intron losses and increased ratios of nonsynonymous to synonymous substitution rates in the clpP gene with intensified positive selection. In addition, two Dipsacoideae plastomes revealed the loss of the plastid-encoded rps15, and a potential functional gene transfer to the nucleus was confirmed.     

Genome-Wide SNP Markers for Genotypic and Phenotypic Differentiation of Melon (Cucumis melo L.) Varieties Using Genotyping-by-Sequencing

Melon (Cucumis melo L.) is an economically important horticultural crop with abundant morphological and genetic variability. Complex genetic variations exist even among melon varieties and remain unclear to date. Therefore, unraveling the genetic variability among the three different melon varieties, muskmelon (C. melo subsp. melo), makuwa (C. melo L. var. makuwa), and cantaloupes (C. melo subsp. melo var. cantalupensis), could provide a basis for evolutionary research. In this study, we attempted a systematic approach with genotyping-by-sequencing (GBS)-derived single nucleotide polymorphisms (SNPs) to reveal the genetic structure and diversity, haplotype differences, and marker-based varieties differentiation. A total of 6406 GBS-derived SNPs were selected for the diversity analysis, in which the muskmelon varieties showed higher heterozygote SNPs. Linkage disequilibrium (LD) decay varied significantly among the three melon varieties, in which more rapid LD decay was observed in muskmelon (r² = 0.25) varieties. The Bayesian phylogenetic tree provided the intraspecific relationships among the three melon varieties that formed, as expected, individual clusters exhibiting the greatest genetic distance based on the posterior probability. The haplotype analysis also supported the phylogeny result by generating three major networks for 48 haplotypes. Further investigation for varieties discrimination allowed us to detect a total of 52 SNP markers that discriminated muskmelon from makuwa varieties, of which two SNPs were converted into cleaved amplified polymorphic sequence markers for practical use. In addition to these markers, the genome-wide association study identified two SNPs located in the genes on chromosome 6, which were significantly associated with the phenotypic traits of melon seed. This study demonstrated that a systematic approach using GBS-derived SNPs could serve to efficiently classify and manage the melon varieties in the genebank.     

Phylogenetic Relationships and Adaptation in Deep-Sea Mussels: Insights from Mitochondrial Genomes

Mitochondrial genomes (mitogenomes) are an excellent source of information for phylogenetic and evolutionary studies, but their application in marine invertebrates is limited. In the present study, we utilized mitogenomes to elucidate the phylogeny and environmental adaptation in deep-sea mussels (Mytilidae: Bathymodiolinae). We sequenced and assembled seven bathymodioline mitogenomes. A phylogenetic analysis integrating the seven newly assembled and six previously reported bathymodioline mitogenomes revealed that these bathymodiolines are divided into three well-supported clades represented by five Gigantidas species, six Bathymodiolus species, and two “Bathymodiolus” species, respectively. A Common interval Rearrangement Explorer (CREx) analysis revealed a gene order rearrangement in bathymodiolines that is distinct from that in other shallow-water mytilids. The CREx analysis also suggested that reversal, transposition, and tandem duplications with subsequent random gene loss (TDRL) may have been responsible for the evolution of mitochondrial gene orders in bathymodiolines. Moreover, a comparison of the mitogenomes of shallow-water and deep-sea mussels revealed that the latter lineage has experienced relaxed purifying selection, but 16 residues of the atp6, nad4, nad2, cob, nad5, and cox2 genes have underwent positive selection. Overall, this study provides new insights into the phylogenetic relationships and mitogenomic adaptations of deep-sea mussels     

Nuclear Microsatellite Primers for the Endangered Relict Fir,Abies pinsapo (Pinaceae) and Cross-Amplification in Related Mediterranean Species

Twelve nuclear microsatellite primers (nSSR) were developed for the endangered species Abies pinsapo Boiss. to enable the study of gene flow and genetic structure in the remaining distribution areas. Microsatellite primers were developed using next-generation sequencing (454) data from a single Abies pinsapo individual. Primers were applied to thirty individuals from the three extant localities. The number of alleles per locus ranged from one to four. Cross-amplification was tested for other Abies species from the Mediterranean Basin, and most of the loci showed higher polymorphisms in the Mediterranean species than in A. pinsapo. These microsatellite markers provide tools for conservation genetic studies in Abies pinsapo as well other Abies species from the Mediterranean Basin.     

Comparative Studies of 5S rDNA Profiles and Cyt b Sequences in two Onychostoma Species (Cyprinidae)

Onychostoma barbatulum and O. alticorpus, two primarily freshwater cyprinid fish, have similar morphological characters and partially overlapping ecological habitats. In order to explore the genetic differences between these two species, chromosomal characteristics and genetic variations were examined by fluorescence in situ hybridization (FISH) of 5S rDNA and cytochrome (Cyt) b gene analysis. Ten specimens of O. barbatulum and O. alticorpus were collected from the Nanzihsian Stream in southern Taiwan. FISH revealed that the 5S rDNA loci of O. barbatulum and O. alticorpus were found at a pericentromeric and subtelomeric position, respectively, in a pair of submetacentric chromosomes. Cyt b genes were amplified and sequenced from five individuals of each species. Intraspecific genetic distances ranged from 0.001–0.004 in O. barbatulum and from 0.001–0.006 in O. alticorpus. Genetic distances between these two species ranged from 0.132–0.142. The phylogenetic tree showed these two species are not sister species. In conclusion, FISH cytogenetic information and Cyt b gene analyses indicated that these two species have significantly different genetic characteristics; nevertheless, their morphological similarities may be due to environmental adaptation.     

Genetic Diversity and Population Differentiation of Calanthe tsoongiana,a Rare and Endemic Orchid in China

Calanthe tsoongiana is a rare terrestrial orchid endemic to China, and this species has experienced severe habitat loss and fragmentation. Inter-simple sequence repeat (ISSR) markers were employed to assess the genetic diversity and differentiation of six populations of C. tsoongiana. Based on 124 discernible fragments yielded by eleven selected primers, high genetic diversity was revealed at the species level; however, genetic diversity at the population level was relatively low. High-level genetic differentiation among populations was detected based on analysis of molecular variance (AMOVA), indicating potential limited gene flow. No significant relationship was observed between genetic and geographic distances among the sampled populations. These results suggested that restricted gene flow might be due to habitat fragmentation and reduced population size as a result of human activities. Based on the findings, several conservation strategies were proposed for the preservation of this threatened species.     

Phylogenetic Relationships between Four Salix L. Species Based on DArT Markers

The objectives of this study were to evaluate the usefulness of DArT markers in genotypic identification of willow species and describe genetic relationships between four willow species: Salix viminalis, S. purpurea, S. alba and S. triandra. The experimental plant material comprised 53 willow genotypes of these four species, which are popularly grown in Poland. DArT markers seem to identify Salix species with a high degree of accuracy. As a result, the examined species were divided into four distinct groups which corresponded to the four analyzed species. In our study, we observed that S. triandra was very different genetically from the other species, including S. alba which is generally classified into the same subgenus of Salix. The above corroborates the findings of other authors who relied on molecular methods to reveal that the classification of S. triandra to the subgenus Salix was erroneous. The Principal Coordinate Analysis (PCoA) and the neighbor-joining dendrogram also confirmed the clear division of the studied willow genotypes into four clusters corresponding to individual species. This confirmed the usefulness of DArT markers in taxonomic analyses and identification of willow species.     

First Report of Crown Gall of Kiwifruit (Actinidia deliciosa) Caused by Agrobacterium fabacearum in China and the Establishment of Loop-Mediated Isothermal Amplification Technique

Kiwifruit is moderately sweet and sour and quite popular among consumers; it has been widely planted in some areas of the world. In 2019, the crown gall disease of kiwifruit was discovered in the main kiwifruit-producing area of Guizhou Province, China. This disease can weaken and eventually cause the death of the tree. The phylogeny, morphological and biological characteristics of the bacteria were described, and were related to diseases. The pathogenicity of this species follows the Koch hypothesis, confirming that A. fabacearum is the pathogen of crown gall disease of kiwifruit in China. In this study, Loop-mediated isothermal amplification (LAMP) analysis for genome-specific gene sequences was developed for the specific detection of A. fabacearum. The detection limit of the LAMP method is 5 × 10⁻⁷ ng/μL, which has high sensitivity. At the same time, the amplified product is stained with SYBR Green I after the reaction is completed, so that the amplification can be detected with the naked eye. LAMP analysis detected the presence of A. fabacearum in the roots and soil samples of the infected kiwifruit plant. The proposed LAMP detection technology in this study offers the advantages of ease of operation, visibility of results, rapidity, accuracy and high sensitivity, making it suitable for the early diagnosis of crown gall disease of kiwifruit.     

Structure and Principal Components Analyses Reveal an Intervarietal Fusion in Malaysian Mistletoe Fig (Ficus deltoidea Jack) Populations

Genetic structure and biodiversity of the medicinal plant Ficus deltoidea have rarely been scrutinized. To fill these lacunae, five varieties, consisting of 30 F. deltoidea accessions were collected across the country and studied on the basis of molecular and morphological data. Molecular analysis of the accessions was performed using nine Inter Simple Sequence Repeat (ISSR) markers, seven of which were detected as polymorphic markers. ISSR-based clustering generated four clusters supporting the geographical distribution of the accessions to some extent. The Jaccard’s similarity coefficient implied the existence of low diversity (0.50–0.75) in the studied population. STRUCTURE analysis showed a low differentiation among the sampling sites, while a moderate varietal differentiation was unveiled with two main populations of F. deltoidea. Our observations confirmed the occurrence of gene flow among the accessions; however, the highest degree of this genetic interference was related to the three accessions of FDDJ10, FDTT16 and FDKT25. These three accessions may be the genetic intervarietal fusion points of the plant’s population. Principal Components Analysis (PCA) relying on quantitative morphological characteristics resulted in two principal components with Eigenvalue >1 which made up 89.96% of the total variation. The cluster analysis performed by the eight quantitative characteristics led to grouping the accessions into four clusters with a Euclidean distance ranged between 0.06 and 1.10. Similarly, a four-cluster dendrogram was generated using qualitative traits. The qualitative characteristics were found to be more discriminating in the cluster and PCA analyses, while ISSRs were more informative on the evolution and genetic structure of the population.