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.

Phylogeny and evolution of Asparagaceae subfamily Nolinoideae: new insights from plastid phylogenomics.

BACKGROUND AND AIMS: Asparagaceae subfamily Nolinoideae is an economically important plant group, but the deep relationships and evolutionary history of the lineage remain poorly understood. Based on a large data set including 37 newly sequenced samples and publicly available plastomes, this study aims to better resolve the inter-tribal relationships of Nolinoideae, and to rigorously examine the tribe-level monophyly of Convallarieae, Ophiopogoneae and Polygonateae. METHODS: Maximum likelihood (ML) and Bayesian inference (BI) methods were used to infer phylogenetic relationships of Nolinoideae at the genus level and above. The diversification history of Nolinoideae was explored using molecular dating. KEY RESULTS: Both ML and BI analyses identically recovered five clades within Nolinoideae, respectively corresponding to Dracaeneae + Rusceae, Polygonateae + Theropogon, Ophiopogoneae, Nolineae, and Convallarieae excluding Theropogon, and most deep nodes were well supported. As Theropogon was embedded in Polygonateae, the plastome phylogeny failed to resolve Convallarieae and Polygonateae as reciprocally monophyletic. Divergence time estimation showed that the origins of most Nolinoideae genera were dated to the Miocene and Pliocene. The youthfulness of Nolinoideae genera is well represented in the three herbaceous tribes (Convallarieae, Ophiopogoneae and Polygonateae) chiefly distributed in temperate areas of the Northern Hemisphere, as the median stem ages of all 14 genera currently belonging to them were estimated at <12.37 Ma. CONCLUSIONS: This study recovered a robust backbone phylogeny, providing new insights for better understanding the evolution and classification of Nolinoideae. Compared with the deep relationships recovered by a previous study based on transcriptomic data, our data suggest that ancient hybridization or incomplete lineage sorting may have occurred in the early diversification of Nolinoideae. Our findings will provide important reference for further study of the evolutionary complexity of Nolinoideae using nuclear genomic data. The recent origin of these herbaceous genera currently belonging to Convallarieae, Ophiopogoneae and Polygonateae provides new evidence to support the hypothesis that the global expansion of temperate habitats caused by the climate cooling over the past 15 million years may have dramatically driven lineage diversification and speciation in the Northern Hemisphere temperate flora.

The Agavoideae: an emergent model clade for CAM evolutionary biology.

Crassulacean acid metabolism - or CAM photosynthesis - was described in the early to mid-20th century, and our understanding of this metabolic pathway was later expanded upon through detailed biochemical analyses of carbon balance. Soon after, scientists began to study the ecophysiological implications of CAM, and a large part of this early work was conducted in the genus Agave, in the subfamily Agavoideae of the family Asparagaceae. Today, the Agavoideae continues to be important for the study of CAM photosynthesis, from the ecophysiology of CAM species, to the evolution of the CAM phenotype and to the genomics underlying CAM traits. Here we review past and current work on CAM in the Agavoideae, in particular highlighting the work of Park Nobel in Agave, and focusing on the powerful comparative system the Agavoideae has become for studying the origins of CAM. We also highlight new genomics research and the potential for studying intraspecific variation within species of the Agavoideae, particularly species in the genus Yucca. The Agavoideae has served as an important model clade for CAM research for decades, and undoubtedly will continue to help push our understanding of CAM biology and evolution in the future.

Differential timing of gene expression and recruitment in independent origins of CAM in the Agavoideae (Asparagaceae).

Crassulacean acid metabolism (CAM) photosynthesis has evolved repeatedly across the plant tree of life, however our understanding of the genetic convergence across independent origins remains hampered by the lack of comparative studies. Here, we explore gene expression profiles in eight species from the Agavoideae (Asparagaceae) encompassing three independent origins of CAM. Using comparative physiology and transcriptomics, we examined the variable modes of CAM in this subfamily and the changes in gene expression across time of day and between well watered and drought-stressed treatments. We further assessed gene expression and the molecular evolution of genes encoding phosphoenolpyruvate carboxylase (PPC), an enzyme required for primary carbon fixation in CAM. Most time-of-day expression profiles are largely conserved across all eight species and suggest that large perturbations to the central clock are not required for CAM evolution. By contrast, transcriptional response to drought is highly lineage specific. Yucca and Beschorneria have CAM-like expression of PPC2, a copy of PPC that has never been shown to be recruited for CAM in angiosperms. Together the physiological and transcriptomic comparison of closely related C3 and CAM species reveals similar gene expression profiles, with the notable exception of differential recruitment of carboxylase enzymes for CAM function.

Plastid phylogenomics of the Sansevieria Clade of Dracaena (Asparagaceae) resolves a recent radiation.

Best known as low maintenance houseplants, sansevierias are a diverse group of flowering plants native to Africa, Madagascar, the Arabian Peninsula, and the Indian subcontinent. Traditionally recognised as a distinct genus, Sansevieria was recently merged with the larger genus Dracaena based on molecular phylogenetic data. Within the Sansevieria Clade of Dracaena, taxonomic uncertainties remain despite attempts to unravel the relationships between the species. To investigate the evolutionary relationships, morphological evolution and biogeographical history in the group, we aim to reconstruct a robust dated phylogenetic hypothesis. Using genome skimming, a chloroplast genome (cpDNA) dataset and a nuclear ribosomal (nrDNA) dataset were generated. The sampling included representatives of all sections and informal groups previously described in Sansevieria based on morphology. Analysis of the cpDNA dataset using a maximum likelihood approach resulted in a well-supported phylogeny. The time-calibrated phylogeny indicated a recent radiation with five main clades emerging in the Pliocene. Two strongly supported clades align with previously defined groups, i.e., Sansevieria section Dracomima, characterised by the Dracomima-type inflorescence, and the Zeylanica informal group, native to the Indian subcontinent. Other previously defined groups were shown to be polyphyletic; a result of convergent evolution of the identifying characters. Switches between flat and cylindrical leaves occurred multiple times in the evolution of the Sansevieria Clade. Similarly, the Cephalantha-type inflorescence has originated multiple times from an ancestor with a Sansevieria-type inflorescence. Analysis of the nrDNA dataset resulted in a phylogenetic hypothesis with low resolution, yet it supported the same two groups confirmed by the cpDNA dataset. This study furthers our understanding of the evolution of the Sansevieria Clade, which will benefit taxonomic and applied research, and aid conservation efforts.

Comparative transcriptome profiling of Polianthes tuberosa during a compatible interaction with root-knot nematode Meloidogyne incognita.

BACKGROUND: The root-knot nematode (RKN; Meloidogyne spp.) is the most destructive plant parasitic nematode known to date. RKN infections, especially those caused by Meloidogyne incognita, are one of the most serious diseases of tuberose. METHODS AND RESULTS: To investigate the molecular mechanism in the host-pathogen interactions, the Illumina sequencing platform was employed to generate comparative transcriptome profiles of uninfected and Meloidogyne incognita-infected tuberose plants, during early, mid, and late infection stage. A total of 7.5 GB (49 million reads) and 9.3 GB (61 million reads) of high-quality data was generated for the control and infected samples, respectively. These reads were combined and assembled using the Trinity assembly program which clustered them into 1,25,060 unigenes. A total of 85,360 validated CDS were obtained from the combined transcriptome whereas 6,795 CDS and 7,778 CDS were found in the data for the control and infected samples, respectively. Gene ontology terms were assigned to 958 and 1,310 CDSs from the control and infected data, respectively. The KAAS pathway analysis revealed that 1,248 CDS in the control sample and 1,482 CDS in the infected sample were enriched with KEGG pathways. The major proportions of CDS were annotated for carbohydrate metabolism, signal transduction and translation related pathways in control and infected samples. Of the 8,289 CDS commonly expressed between the control and infected plants, 256 were significantly upregulated and 129 were significantly downregulated in the infected plants. CONCLUSIONS: Collectively, our results provide a comprehensive gene expression changes in tuberose during its association with RKNs and point to candidate genes that are involved in nematode stress signaling for further investigation. This is the first report addressing genes associated with M. incognita-tuberose interaction and the results have important implications for further characterization of RKN resistance genes in tuberose.

Intracellular DNA transfer events restricted to the genus Convallaria within the Asparagaceae family: Possible mechanisms and potential as genetic markers for biographical studies.

Intracellular gene transfer among plant genomes is a common phenomenon. Due to their high conservation and high plastid membrane integrity, chloroplast (cp) genomes incorporate foreign genetic material very rarely. Convallaria is a small monocotyledonous genus consisting of C. keiskei, C. majalis and C. montana. Here, we characterized, analyzed and identified 3.3 and 3.7 kb of mitochondrial DNA sequences in the plastome (MCP) of C. majalis and C. montana, respectively. We identified 6 bp and 23 bp direct repeats and mitochondrial pseudogenes, with rps3, rps19 and rpl10 identified in the MCP region. Additionally, we developed novel plastid molecular genetic markers to differentiate Convallaria spp. based on 21 populations. BEAST and biogeographical analyses suggested that Convallaria separated into Eurasian and North American lineages during the middle Pliocene and originated in East Asia. Vicariance in the genus was followed by dispersal into Europe and southeastern North America. These analyses indicate that the MCP event was restricted to the genus Convallaria of Asparagaceae, in contrast to similar events that occurred in its common ancestors with other families of land plants. However, further mitochondrial and population studies are necessary to understand the integration of the MCP region and gene flow in the genus Convallaria.

Evolution of Lomandroideae: Multiple origins of polyploidy and biome occupancy in Australia.

Asparagaceae: Lomandroideae are a species-rich and economically important subfamily in the monocot order Asparagales, with a center of diversity in Australia. Lomandroideae are ecologically diverse, occupying mesic and arid biomes in Australia and possessing an array of key traits, including sexual dimorphism, storage organs and polyploidy that are potentially adaptive for survival in seasonally arid and fire-dependent habitats. The Lomandroideae phylogeny was reconstructed using maximum likelihood and Bayesian inference criteria, based on plastome data from genome-skimming to infer relationships. A fossil-calibrated chronogram provided a temporal framework for understanding trait transitions. Ancestral state reconstructions and phylogenetic comparative trait correlation analyses provided insights into the evolutionary and ecological drivers associated with Lomandroideae diversification. Lomandroideae diverged from the other Asparagaceae ca. 56.61 million years ago (95% highest posterior density values 70.31-45.34 million years) and the major lineages diversified since the Oligocene. The most recent common ancestor of the clade likely occupied the mesic biome, was hermaphroditic and geophytic. Biome occupancy transitions were correlated with polyploidy and the presence of storage roots. Polyploidy potentially serves as an "enabler" trait, generating novel phenotypes, which may confer tolerance to climatic ranges and soil conditions putatively required for expansion into and occupation of new arid biomes. Storage roots, as a key factor driving biome transitions, may have been associated with fire rather than with aridification events in the Australian flora. This study contributes significantly to our understanding of biome evolution by identifying polyploidy and storage organs as key factors associated with transitions in biome occupancy in this lineage.

Development of 19 novel microsatellite markers of lily-of-the-valley (Convallaria, Asparagaceae) from transcriptome sequencing.

The lily-of-the-valley Convallaria (Asparagaceae) consists of three herbaceous perennial species. The plants are commonly found in northern hemisphere, and are best-known for their ornamental and pharmaceutical value. In order to assess the genetic structure, diversity and demographic history of Convallaria species, 19 novel microsatellite markers were developed based on transcriptome data of C. keiskei. Polymorphism and cross-amplification of the markers were tested in three populations of C. keiskei and one population each of C. majalis and C. montana. The transferability rate in two species was both 89.5%. The average number of alleles detected per locus was 7.7, 3.3 and 2.7 in C. keiskei, C. majalis and C. montana, respectively, and the polymorphism information content correspondingly varied from 0.067 to 0.730, from 0.071 to 0.637 and from 0.195 to 0.680 at the population level. The observed and expected heterozygosity ranged from 0.000 to 1.000 and from 0.000 to 0.833, respectively. Seven of the 19 loci showed significant deviation from Hardy-Weinberg equilibrium. The availability of these markers will provide a useful molecular tool for further population genetics, phylogeographic and breeding studies of Convallaria species.

A novel R2R3-MYB from grape hyacinth, MaMybA, which is different from MaAN2, confers intense and magenta anthocyanin pigmentation in tobacco.

BACKGROUND: The primary pigments in flowers are anthocyanins, the biosynthesis of which is mainly regulated by R2R3-MYBs. Muscari armeniacum is an ornamental garden plant with deep cobalt blue flowers containing delphinidin-based anthocyanins. An anthocyanin-related R2R3-MYB MaAN2 has previously been identified in M. armeniacum flowers; here, we also characterized a novel R2R3-MYB MaMybA, to determine its function and highlight similarities and differences between MaMybA and MaAN2. RESULTS: In this study, a novel anthocyanin-related R2R3-MYB gene was isolated from M. armeniacum flowers and functionally identified. A sequence alignment showed that MaMybA contained motifs typically conserved with MaAN2 and its orthologs. However, the shared identity of the entire amino acid sequence between MaMybA and MaAN2 was 43.5%. Phylogenetic analysis showed that they were both clustered into the AN2 subgroup of the R2R3-MYB family, but not in the same branch. We also identified a IIIf bHLH protein, MabHLH1, in M. armeniacum flowers. A bimolecular fluorescence complementation assay showed that MabHLH1 interacted with MaMybA or MaAN2 in vivo; a dual luciferase assay indicated that MaMybA alone or in interaction with MabHLH1 could regulate the expression of MaDFR and AtDFR, but MaAN2 required MabHLH1 to do so. When overexpressing MaMybA in Nicotiana tabacum 'NC89', the leaves, petals, anthers, and calyx of transgenic tobacco showed intense and magenta anthocyanin pigments, whereas those of OE-MaAN2 plants had lighter pigmentation. However, the ovary wall and seed skin of OE-MaMybA tobacco were barely pigmented, while those of OE-MaAN2 tobacco were reddish-purple. Moreover, overexpressing MaMybA in tobacco obviously improved anthocyanin pigmentation, compared to the OE-MaAN2 and control plants, by largely upregulating anthocyanin biosynthetic and endogenous bHLH genes. Notably, the increased transcription of NtF3'5'H in OE-MaMybA tobacco might lead to additional accumulation of delphinidin 3-rutinoside, which was barely detected in OE-MaAN2 and control plants. We concluded that the high concentration of anthocyanin and the newly produced Dp3R caused the darker color of OE-MaMybA compared to OE-MaAN2 tobacco. CONCLUSION: The newly identified R2R3-MYB transcription factor MaMybA functions in anthocyanin biosynthesis, but has some differences from MaAN2; MaMybA could also be useful in modifying flower color in ornamental plants.

Shifts in gene expression profiles are associated with weak and strong Crassulacean acid metabolism.

PREMISE OF THE STUDY: The relative ease of high throughput sequencing is facilitating comprehensive phylogenomic and gene expression studies, even for nonmodel groups. To date, however, these two approaches have not been merged; while phylogenomic methods might use transcriptome sequences to resolve relationships, assessment of gene expression patterns in a phylogenetic context is less common. Here we analyzed both carbon assimilation and gene expression patterns of closely related species within the Agavoideae (Asparagaceae) to elucidate changes in gene expression across weak and strong phenotypes for Crassulacean acid metabolism (CAM). METHODS: Gene expression patterns were compared across four genera: Agave (CAM), which is paraphyletic with Polianthes (weak CAM) and Manfreda (CAM), and Beschorneria (weak CAM). RNA-sequencing was paired with measures of gas exchange and titratable acidity. Climate niche space was compared across the four lineages to examine abiotic factors and their correlation to CAM. KEY RESULTS: Expression of homologous genes showed both shared and variable patterns in weak and strong CAM species. Network analysis highlights that despite shared expression patterns, highly connected genes differ between weak and strong CAM, implicating shifts in regulatory gene function as key for the evolution of CAM. Variation in carbohydrate metabolism between weak and strong CAM supports the importance of sugar turnovers for CAM physiology. CONCLUSIONS: Integration of phylogenetics and RNA-sequencing provides a powerful tool to study the evolution of CAM photosynthesis across closely related but photosynthetically variable species. Our findings regarding shared or shifted gene expression and regulation of CAM via carbohydrate metabolism have important implications for efforts to engineer the CAM pathway into C3 food and biofuel crops.

Small RNA profiling for identification of miRNAs involved in regulation of saponins biosynthesis in Chlorophytum borivilianum.

BACKGROUND: MicroRNAs act as molecular regulator of cell signaling, plant growth and development, and regulate various primary and secondary plant metabolic processes. In the present study, deep sequencing of small RNAs was carried out to identify known and novel miRNAs from pharmaceutically important plant, Chlorophytum borivilianum. RESULTS: Total 442 known miRNAs and 5 novel miRNAs were identified from young leaf small RNA library. Experimental validation with stem loop RT-PCR confirmed the in silico identification. Based on transcriptome data of root and leaf of C. borivilianum, Oryza sativa, and Arabidopsis thaliana target gene prediction was done using psRNAtarget and mirRanda. BLAST2GO helped in localization of predicted targets and KEGG (Kyoto Encyclopedia for Genes and Genomes) pathway analysis concluded that miR9662, miR894, miR172, and miR166 might be involved in regulating saponin biosynthetic pathway. The correlation between miRNA and its target gene was further validated by RT-qPCR analysis. CONCLUSION: This study provides first elaborated glimpse of miRNA pool of C. borivilianum, which can help to understand the miRNA dependent regulation of saponin biosynthesis and to design further metabolic engineering experiment to enhance their contents in the plant.

Structural polymorphisms and distinct genomic composition suggest recurrent origin and ongoing evolution of B chromosomes in the Prospero autumnale complex (Hyacinthaceae).

Supernumerary B chromosomes (Bs) are genomic parasitic components, originating from the A complement via chromosomal rearrangements, which follow their own evolutionary trajectories. They often contain repetitive DNAs, some shared with regular chromosomes and some newly evolved. Genomic composition, origin and evolution of Bs have been analysed in the chromosomally variable Prospero autumnale complex. Two rDNAs and a satellite DNA (PaB6) from regular chromosomes were mapped to Bs of 26 plants from three diploid cytotypes, their hybrids and polyploid derivatives. In homoploid diploid hybrids, genomic in situ hybridization (GISH) allowed B painting with the parental DNAs. Bs were structurally variable and highly enriched in 5S rDNA and satDNA PaB6, and rarely in 35S rDNA. Eleven combinations of rDNA and PaB6 localization were observed. The quantities of PaB6 in Bs and regular chromosomes were not correlated, suggesting amplification mechanisms other than recombination. PaB6 and 5S rDNA amounts increased with increasing ploidy level. GISH revealed two independent origins of Bs. The structural variation, repeat content, repeat-type fluctuations and differing genomic affinities of Bs in different cytotypes suggest that they represent young proto-B chromosomes. Bs in P. autumnale probably form recurrently as by-products of the extensive genome restructuring within this chromosomally variable species complex.

Timing of rapid diversification and convergent origins of active pollination within Agavoideae (Asparagaceae).

PREMISE OF THE STUDY: Yucca species are ideal candidates for the study of coevolution due to the obligate mutualism they form with yucca moth pollinators (genera Tegeticula and Parategeticula). Yuccas are not the only species to exhibit a mutualism with yucca moths; the genus Hesperoyucca is pollinated by the California yucca moth (Tegeticula maculata). Relationships among yuccas, Hesperoyucca, and other members of subfamily Agavoideae are necessary to understand the evolution of this unique pollination syndrome. Here, we investigate evolutionary relationships of yuccas and closely related genera looking at the timing and origin of yucca moth pollination. METHODS: In this study, we sequenced the chloroplast genomes of 20 species in the subfamily Agavoideae (Asparagaceae) and three confamilial outgroup taxa to resolve intergeneric phylogenetic relationships of Agavoideae. We estimated divergence times using protein-coding genes from 67 chloroplast genomes sampled across monocots to determine the timing of the yucca moth pollination origin. KEY RESULTS: We confidently resolved intergeneric relationships in Agavoideae, demonstrating the origin of the yucca-yucca moth mutualism on two distinct lineages that diverged 27 million years ago. Comparisons of Yucca and Hesperoyucca divergence time to those of yucca moths (Tegeticula and Parategeticula, Prodoxidae) indicate overlapping ages for the origin of pollinating behavior in the moths and pollination by yucca moths in the two plant lineages. CONCLUSION: Whereas pollinating yucca moths have been shown to have a single origin within the Prodoxidae, there were independent acquisitions of active pollination on lineages leading to Yucca and Hesperoyucca within the Agavoideae.

Population Genetics of the Endemic Hawaiian Species Chrysodracon hawaiiensis and Chrysodracon auwahiensis (Asparagaceae): Insights from RAPD and ISSR Variation.

The genus Chrysodracon has six endemic species in the Hawaii Islands. Chrysodracon hawaiiensis is endemic to Hawaii Island and was described as a distinct species in 1980. It was listed as an endangered species on the International Union for the Conservation of Nature and Natural Resources (IUCN) Red List in 1997. This woody plant species was, at one time, common in exposed dry forests, but it became very rare due to grazing pressure and human development. The tree species Chrysodracon auwahiensis (C. auwahiensis), endemic to Maui and Molokai, still has large adult populations in dry lands of the islands, but unfortunately no regeneration from seed has been reported in those areas for many years. The two endemic species were examined using the molecular technique of random amplified polymorphic DNA (RAPD) and inter simple sequence repeats (ISSR) to determine the genetic structure of the populations and the amount of variation. Both species possess similar genetic structure. Larger and smaller populations of both species contain similar levels of genetic diversity as determined by the number of polymorphic loci, estimated heterozygosity, and Shannon's index of genetic diversity. Although population diversity of Chrysodracon hawaiiensis (C. hawaiiensis) is thought to have remained near pre-disturbance levels, population size continues to decline as recruitment is either absent or does not keep pace with senescence of mature plants. Conservation recommendations for both species are suggested.

Comparative and phylogenetic analysis of the complete chloroplast genomes of six Polygonatum species (Asparagaceae).

Polygonatum Miller belongs to the tribe Polygonateae of Asparagaceae. The horizontal creeping fleshy roots of several species in this genus serve as traditional Chinese medicine. Previous studies have mainly reported the size and gene contents of the plastomes, with little information on the comparative analysis of the plastid genomes of this genus. Additionally, there are still some species whose chloroplast genome information has not been reported. In this study, the complete plastomes of six Polygonatum were sequenced and assembled, among them, the chloroplast genome of P. campanulatum was reported for the first time. Comparative and phylogenetic analyses were then conducted with the published plastomes of three related species. Results indicated that the whole plastome length of the Polygonatum species ranged from 154,564 bp (P. multiflorum) to 156,028 bp (P. stenophyllum) having a quadripartite structure of LSC and SSC separated by two IR regions. A total of 113 unique genes were detected in each of the species. Comparative analysis revealed that gene content and total GC content in these species were highly identical. No significant contraction or expansion was observed in the IR boundaries among all the species except P. sibiricum1, in which the rps19 gene was pseudogenized owing to incomplete duplication. Abundant long dispersed repeats and SSRs were detected in each genome. There were five remarkably variable regions and 14 positively selected genes were identified among Polygonatum and Heteropolygonatum. Phylogenetic results based on chloroplast genome strongly supported the placement of P. campanulatum with alternate leaves in sect. Verticillata, a group characterized by whorled leaves. Moreover, P. verticillatum and P. cyrtonema were displayed as paraphyletic. This study revealed that the characters of plastomes in Polygonatum and Heteropolygonatum maintained a high degree of similarity. Five highly variable regions were found to be potential specific DNA barcodes in Polygonatum. Phylogenetic results suggested that leaf arrangement was not suitable as a basis for delimitation of subgeneric groups in Polygonatum and the definitions of P. cyrtonema and P. verticillatum require further study.

AFLP analysis and zebra disease resistance identification of 40 sisal genotypes in China.

Sisal is the most important fiber crop in tropical and subtropical areas in China and the world. Zebra disease is a serious threat to the main cultivar Agave hybrid No.11648 (H.11648) worldwide. To select germplasm materials with zebra disease resistance for breeding, the fluorescent amplified fragment length polymorphism (AFLP) technique was used to make a cluster analysis of the genetic relationships of 40 sisal genotypes grown in China, and Phytophthora nicotianae was used to inoculate the 40 genotypes to identify their resistance to zebra disease. As a result, the similarity coefficient among 40 sisal genotypes was found to be 0.44-0.83 and the 40 genotypes show different levels of disease resistance. According to the AFLP analysis, the disease resistance and chromosomal ploidy, it can be reasoned that, A. attenuata var. marginata, Dong 109, Nan ya 1 and A. attenuata are suitable for hybridization with H.11648 to breed a new disease-resistant variety.

Evolutionary Comparison of the Complete Chloroplast Genomes in Convallaria Species and Phylogenetic Study of Asparagaceae.

The genus Convallaria (Asparagaceae) comprises three herbaceous perennial species that are widely distributed in the understory of temperate deciduous forests in the Northern Hemisphere. Although Convallaria species have high medicinal and horticultural values, studies related to the phylogenetic analysis of this genus are few. In the present study, we assembled and reported five complete chloroplast (cp) sequences of three Convallaria species (two of C. keiskei Miq., two of C. majalis L., and one of C. montana Raf.) using Illumina paired-end sequencing data. The cp genomes were highly similar in overall size (161,365-162,972 bp), and all consisted of a pair of inverted repeats (IR) regions (29,140-29,486 bp) separated by a large single-copy (LSC) (85,183-85,521 bp) and a small single-copy (SSC) region (17,877-18,502 bp). Each cp genome contained the same 113 unique genes, including 78 protein-coding genes, 30 transfer RNA genes, and 4 ribosomal RNA genes. Gene content, gene order, AT content and IR/SC boundary structure were nearly identical among all of the Convallaria cp genomes. However, their lengths varied due to contraction/expansion at the IR/LSC borders. Simple sequence repeat (SSR) analyses indicated that the richest SSRs are A/T mononucleotides. Three highly variable regions (petA-psbJ, psbI-trnS and ccsA-ndhD) were identified as valuable molecular markers. Phylogenetic analysis of the family Asparagaceae using 48 cp genome sequences supported the monophyly of Convallaria, which formed a sister clade to the genus Rohdea. Our study provides a robust phylogeny of the Asparagaceae family. The complete cp genome sequences will contribute to further studies in the molecular identification, genetic diversity, and phylogeny of Convallaria.

A novel R3 MYB transcriptional repressor, MaMYBx, finely regulates anthocyanin biosynthesis in grape hyacinth.

R3-MYBs negatively regulate anthocyanin pigmentation in plants. However, how R3-MYB repressors finely modulate anthocyanin biosynthesis in cooperation with R2R3-MYB activators remains unclear in monocots. We previously identified two anthocyanin-related R2R3-MYB activators (MaMybA and MaAN2) in grape hyacinth (Muscari spp.). Here, we isolated a R3-MYB repressor, MaMYBx, and characterized its role in anthocyanin biosynthesis using genetic and biochemical markers. The temporal expression pattern of MaMYBx was similar to that of MaMybA and MaAN2, and it was correlated with anthocyanin accumulation during flower development. MaMYBx could be activated either by MaMybA alone or by MaMybA/MaAN2 and cofactor MabHLH1, and it suppressed its own activation and that of MaMybA promoters mediated by MaMybA/MaAN2 and MabHLH1. Like MaMybA, MaMYBx interacted with MabHLH1. MaDFR and MaANS transcription and anthocyanin accumulation mediated by MaMybA/MaAN2 and MabHLH1 were inhibited by MaMYBx. Overexpression of MaMYBx in tobacco greatly reduced flower pigmentation and repressed the expression of late structural and regulatory anthocyanin pathway genes. Thus, MaMYBx finely regulates anthocyanin biosynthesis by binding to MabHLH1 and disrupting the R2R3 MYB-bHLH complex in grape hyacinth. The regulatory network of transcriptional activators and repressors modulating anthocyanin biosynthesis is conserved within monocots. MaMYBx seems a potentially valuable target for flower color modification in ornamental plants.

Phylogeography of Camassia quamash in western North America: postglacial colonization and transport by indigenous peoples.

Recent human activities have spread numerous plant species across the globe, yet it is unclear to what degree historical human activities influenced plant dispersal. In western North America, Camassia quamash was one of the most important food plants for indigenous peoples, who transported its propagules either intentionally or accidentally. We investigated how human and natural dispersal might have contributed to the current pattern of spatial genetic structure in C. quamash by performing phylogeographical surveys at two geographical scales. We sequenced two noncoding regions of chloroplast deoxyribonucleic acid (DNA) in 226 individuals from 53 populations of C. quamash as well as 126 individuals from 21 populations of the non-food plant Zigadenus venenosus. Contrary to the expectation of anthropogenic transport, C. quamash populations did not exhibit weaker genetic structure than Z. venenosus populations. We also failed to find convincing evidence for signatures of transport. Instead, our data showed strong effects of past glaciation and geographical barriers of the mountains in the Cascade Range, Olympic Peninsula and Vancouver Island. West of the Cascades, the species appears to have largely migrated northward from a southern refugium after deglaciation, whereas few populations having a highly divergent haplotype might have survived in southwestern Washington. Our data suggest that despite substantial ethnobotanical evidence for anthropogenic transport, the current pattern of genetic structure of C. quamash does not show any detectable signatures of transport by indigenous peoples and is better understood as the result of natural dispersal processes.