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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

A molecular phylogeny of the genus Drimia (Asparagaceae: Scilloideae: Urgineeae) in India inferred from non-coding chloroplast and nuclear ribosomal DNA sequences.

The evolutionary history of the medicinally important bulbous geophyte Drimia (subfamily: Scilloideae) has long been considered as a matter of debate in the monocot systematics. In India the genus is represented by a species complex, however, the taxonomic delimitation among them is ill-defined till date. In the present study, a comprehensive phylogenetic relationship among Indian species of this genus has been inferred for the first time based on chloroplast DNA trnL intron, rps16-trnK intergenic spacer, atpB-rbcL intergenic spacer and ribosomal DNA ITS1-5.8S-ITS2 sequences, leaf morphology, anatomy, stomatal characteristics and pollen exine ornamentations. The present findings revealed the monophyletic origin of the Indian members of Drimia and grouped them into two possible lineages (clade- I and II). The phylogenetic tree based on cpDNA concatenated sequences further resolved the clade-I into two distinct subclades (I and II) and clarified the intraspecies relationship among the studied members. The present study suggested a strong relationship between the molecular phylogeny and the morphological characteristics of the species studied. A possible trend of evolution of two important traits: 'type of palisade cells' in leaf and 'pollen exine patterns' among the members of Drimia in India was also suggested.

Exogenous 24-epibrassinolide enhanced benzene detoxification in Chlorophytum comosum via overexpression and conjugation by glutathione.

Benzene, a hydrophobic xenobiotic, induces cell damage in both humans and plants. Due to its volatilization, benzene is an airborne environmental problem. The potential of an exogenous bioactive brassinosteroid phytohormone to enhance benzene removal for phytoremediation was investigated. Chlorophytum comosum had higher brassinosteroids content under benzene stress. Plant treated with 24-epibrassinolide (EBR) removed significantly more gaseous benzene than untreated plants under both light and dark conditions at an initial benzene of 12.75 µmol in the systematic chambers (P < 0.05). Although benzene increased malondialdehyde in plant tissue, EBR-treated plants lowered this lipid peroxidation by enhancing their antioxidant content and increasing benzene detoxification-related genes expression, including ascorbic acid (AsA), homogentisate phytyltransferase (HPT), and glutathione synthethase (GS). This contributed to maintaining higher photosynthetic performances. Moreover, EBR-treated plants had higher gene expression of ferredoxin-NADP reductase (FNR) and glucose-6-phosphate 1-dehydrogenase (G6PDH), thus promoting NADPH biosynthesis to cope with benzene under light and dark conditions, respectively. Further, higher glutathione biosynthesis promoted more glutathione conjugate of benzene products including S-phenylcysteine (SPC) in EBR-treated plants. Hence, application of exogenous EBR as foliar spray provided for enhanced benzene detoxification via antioxidant content, benzene detoxification-related genes and benzene conjugation products with glutathione (GSH) and consequently greater gaseous benzene removal.

Genetic structuring of the coastal herb Arthropodium cirratum (Asparagaceae) is shaped by low gene flow, hybridization and prehistoric translocation.

We examined the genetic structuring of rengarenga (Arthropodium cirratum; Asparagaceae), an endemic New Zealand coastal herb, using nuclear microsatellite markers. This species was brought into cultivation by Maori within the last 700-800 years for its edible roots and was transplanted beyond its natural distribution as part of its cultivation. We found very high levels of genetic structuring in the natural populations (FST = 0.84), indicating low levels of gene flow. Reduced genetic diversity was found in the translocated populations, suggesting a large loss of genetic diversity early in the domestication process. The data indicates that rengarenga was brought into cultivation independently at least three times, with the sources of these introductions located within a narrow area encompassing about 250km of coastline. Hybridization was inferred between A. cirratum and the closely related A. bifurcatum, despite A. birfucatum not occurring in the vicinity.

Transcriptome analysis of Polianthes tuberosa during floral scent formation.

Polianthes tuberosa is a popular ornamental plant. Its floral scent volatiles mainly consist of terpenes and benzenoids that emit a charming fragrance. However, our understanding of the molecular mechanism responsible for the floral scent of P. tuberosa is limited. Using transcriptome sequencing and de novo assembly, a total of 228,706,703 high-quality reads were obtained, which resulted in the identification of 96,906 unigenes (SRA Accession Number SRP126470, TSA Acc. No. GGEA00000000). Approximately 41.85% of the unigenes were functionally annotated using public databases. A total of 4,694 differentially expressed genes (DEGs)were discovered during flowering. Gas chromatography-mass spectrometry analysis revealed that the majority of the volatiles comprised benzenoids and small amounts of terpenes. Homology analysis identified 13 and 17 candidate genes associated with terpene and benzenoid biosynthesis, respectively. Among these, PtTPS1, PtDAHPSs, PtPAL1, and PtBCMT2 might play important roles in regulating the formation of floral volatiles. The data generated by transcriptome sequencing provide a critical resource for exploring concrete characteristics as well as for supporting functional genomics studies. The results of the present study also lay the foundation for the elucidation of the molecular mechanism underlying the regulation of floral scents in monocots.

Long-term monitoring of B-chromosome invasion and neutralization in a population of Prospero autumnale (Asparagaceae).

B chromosomes have been reported in about 15% of eukaryotes, but long-term dynamics of B chromosomes in a single natural population has rarely been analyzed. Prospero autumnale plants collected in 1981 and 1983 at Cuesta de La Palma population had shown the presence of B chromosomes. We analyze here seven additional samples collected between 1987 and 2015, and show that B frequency increased significantly during the 1980s and showed minor fluctuations between 2005 and 2015. A mother-offspring analysis of B chromosome transmission, at population level, showed significant drive on the male side (kB  = 0.65) and significant drag on the female side (kB  = 0.33), with average B transmission rate being very close to the Mendelian rate (0.5). No significant effects of B chromosomes were observed on a number of vigor and fertility-related traits. Within a parasite/host framework, these results suggest that B chromosomes' drive on the male side is the main pathway for B chromosome invasion, whereas B chromosome drag on the female side might be the main manifestation of host genome resistance in this species. Prospero autumnale thus illuminates a novel evolutionary pathway for B chromosome neutralization by means of a decrease in B transmission through the nondriving sex.

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.