Long-term genome diploidization in allopolyploid Nicotiana section Repandae (Solanaceae)

Chromosome-level genome assemblies of Nicotiana tabacum, Nicotiana sylvestris, and Nicotiana tomentosiformis

The Solanaceae species Nicotiana tabacum, an economically important crop plant cultivated worldwide, is an allotetraploid species that appeared about 200,000 years ago as the result of the hybridization of diploid ancestors of Nicotiana sylvestris and Nicotiana tomentosiformis. The previously published genome assemblies for these three species relied primarily on short-reads, and the obtained pseudochromosomes only partially covered the genomes. In this study, we generated annotated de novo chromosome-level genomes of N. tabacum, N. sylvestris, and N. tomentosiformis, which contain 3.99 Gb, 2.32 Gb, and 1.74 Gb, respectively of sequence data, with 97.6%, 99.5%, and 95.9% aligned in chromosomes, and represent 99.2%, 98.3%, and 98.5% of the near-universal single-copy orthologs Solanaceae genes. The completion levels of these chromosome-level genomes for N. tabacum, N. sylvestris, and N. tomentosiformis are comparable to other reference Solanaceae genomes, enabling more efficient synteny-based cross-species research.

COL2-dependent photoperiodic floral induction in Nicotiana sylvestris seems to be lost in the N. sylvestris × N. tomentosiformis hybrid N. tabacum

Introduction

Plants are sessile organisms that maximize reproductive success by adapting to their environment. One of the key steps in the reproductive phase of angiosperms is flower development, requiring the perception of multiple endogenous and exogenous signals integrated via a complex regulatory network. Key floral regulators, including the main transcription factor of the photoperiodic pathway (CONSTANS, CO) and the central floral pathway integrator (FLOWERING LOCUS T, FT), are known in many species.

Methods and results

We identified several CO-like (COL) proteins in tobacco (Nicotiana tabacum). The NtCOL2a/b proteins in the day-neutral plant N. tabacum were most closely related to Arabidopsis CO. We characterized the diurnal expression profiles of corresponding genes in leaves under short-day (SD) and long-day (LD) conditions and confirmed their expression in phloem companion cells. Furthermore, we analyzed the orthologs of NtCOL2a/b in the maternal LD ancestor (N. sylvestris) and paternal, facultative SD ancestor (N. tomentosiformis) of N. tabacum and found that they were expressed in the same diurnal manner. NtCOL2a/b overexpression or knock-out using the CRISPR/Cas9 system did not support a substantial role for the CO homologs in the control of floral transition in N. tabacum. However, NsCOL2 overexpression induced flowering in N. sylvestris under typically non-inductive SD conditions, correlating with the upregulation of the endogenous NsFTd gene.

Discussion

Our results suggest that NsFTd is transcriptionally regulated by NsCOL2 and that this COL2-dependent photoperiodic floral induction seems to be lost in N. tabacum, providing insight into the diverse genetics of photoperiod-dependent flowering in different Nicotiana species.

Phylotranscriptomics supports numerous polyploidization events and phylogenetic relationships in Nicotiana

Introduction

Nicotiana L. (Solanaceae) is of great scientific and economic importance, and polyploidization has been pivotal in shaping this genus. Despite many previous studies on the Nicotiana phylogenetic relationship and hybridization, evidence from whole genome data is still lacking.

Methods

In this study, we obtained 995 low-copy genes and plastid transcript fragments from the transcriptome datasets of 26 Nicotiana species, including all sections. We reconstructed the phylogenetic relationship and phylogenetic network of diploid species.

Results

The incongruence among gene trees showed that the formation of N. sylvestris involved incomplete lineage sorting. The nuclear-plastid discordance and nuclear introgression absence indicated that organelle capture from section Trigonophyllae was involved in forming section Petunioides. Furthermore, we analyzed the evolutionary origin of polyploid species and dated the time of hybridization events based on the analysis of PhyloNet, sequence similarity search, and phylogeny of subgenome approaches. Our results highly evidenced the hybrid origins of five polyploid sections, including sections Nicotiana, Repandae, Rusticae, Polydicliae, and Suaveolentes. Notably, we provide novel insights into the hybridization event of section Polydicliae and Suaveolentes. The section Polydicliae formed from a single hybridization event between maternal progenitor N. attenuata and paternal progenitor N. undulata; the N. sylvestris (paternal progenitor) and the N. glauca (maternal progenitor) were involved in the formation of section Suaveolentes.

Discussion

This study represents the first exploration of Nicotiana polyploidization events and phylogenetic relationships using the high-throughput RNA-seq approach. It will provide guidance for further studies in molecular systematics, population genetics, and ecological adaption studies in Nicotiana and other related species.

Lack of ITS sequence homogenization in Erysimum species (Brassicaceae) with different ploidy levels

The internal transcribed spacers (ITS) exhibit concerted evolution by the fast homogenization of these sequences at the intragenomic level. However, the rate and extension of this process are unclear and might be conditioned by the number and divergence of the different ITS copies. In some cases, such as hybrid species and polyploids, ITS sequence homogenization appears incomplete, resulting in multiple haplotypes within the same organism. Here, we studied the dynamics of concerted evolution in 85 individuals of seven plant species of the genus Erysimum (Brassicaceae) with multiple ploidy levels. We estimated the rate of concerted evolution and the degree of sequence homogenization separately for ITS1 and ITS2 and whether these varied with ploidy. Our results showed incomplete sequence homogenization, especially for polyploid samples, indicating a lack of concerted evolution in these taxa. Homogenization was usually higher in ITS2 than in ITS1, suggesting that concerted evolution operates more efficiently on the former. Furthermore, the hybrid origin of several species appears to contribute to the maintenance of high haplotype diversity, regardless of the level of ploidy. These findings indicate that sequence homogenization of ITS is a dynamic and complex process that might result in varying intra- and inter-genomic diversity levels.

Comparative Chloroplast Genomes of Nicotiana Species (Solanaceae): Insights Into the Genetic Variation, Phylogenetic Relationship, and Polyploid Speciation

Nicotiana L. is a genus rich in polyploidy, which represents an ideal natural system for investigating speciation, biodiversity, and phytogeography. Despite a wealth of phylogenetic work on this genus, a robust evolutionary framework with a dated molecular phylogeny for the genus is still lacking. In this study, the 19 complete chloroplast genomes of Nicotiana species were assembled, and five published chloroplast genomes of Nicotiana were retrieved for comparative analyses. The results showed that the 24 chloroplast genomes of Nicotiana, ranging from 155,327 bp (N. paniculata) to 156,142 bp (N. heterantha) in size, exhibited typical quadripartite structure. The chloroplast genomes were rather conserved in genome structure, GC content, RNA editing sites, and gene content and order. The higher GC content observed in the IR regions could be a result of the presence of abundant rRNA and tRNA genes, which contained a relatively higher GC content. A total of seven hypervariable regions, as new molecular markers for phylogenetic analysis, were uncovered. Based on 78 protein-coding genes, we constructed a well-supported phylogenetic tree, which was largely in agreement with previous studies, except for a slight conflict in several sections. Chloroplast phylogenetic results indicated that the progenitors of diploid N. sylvestris, N. knightiana, and the common ancestor of N. sylvestris and N. glauca might have donated the maternal genomes of allopolyploid N. tabacum, N. rustica, and section Repandae, respectively. Meanwhile, the diploid section Noctiflorae lineages (N. glauca) acted as the most likely maternal progenitor of section Suaveolentes. Molecular dating results show that the polyploid events range considerably in ~0.12 million (section Nicotiana) to ~5.77 million (section Repandae) years ago. The younger polyploids (N. tabacum and N. rustica) were estimated to have arisen ~0.120 and ~0.186 Mya, respectively. The older polyploids (section Repandae and Suaveolentes) were considered to have originated from a single polyploid event at ~5.77 and ~4.49 Mya, respectively. In summary, the comparative analysis of chloroplast genomes of Nicotiana species has not only revealed a series of new insights into the genetic variation and phylogenetic relationships in Nicotiana but also provided rich genetic resources for speciation and biodiversity research in the future.

Molecular and Cytogenetic Analysis of rDNA Evolution in Crepis Sensu Lato

Although Crepis was the first model plant group in which chromosomal changes were considered to play an important role in speciation, their chromosome structure and evolution have been barely investigated using molecular cytogenetic methods. The aim of the study was to provide a better understanding of the patterns and directions of Crepis chromosome evolution, using comparative analyses of rDNA loci number and localisation. The chromosome base number and chromosomal organisation of 5S and 35S rDNA loci were analysed in the phylogenetic background for 39 species of Crepis, which represent the evolutionary lineages of Crepis sensu stricto and Lagoseris, including Lapsana communis. The phylogenetic relationships among all the species were inferred from nrITS and newly obtained 5S rDNA NTS sequences. Despite high variations in rDNA loci chromosomal organisation, most species had a chromosome with both rDNA loci within the same (usually short) chromosomal arm. The comparative analyses revealed several independent rDNA loci number gains and loci repositioning that accompanied diversification and speciation in Crepis. Some of the changes in rDNA loci patterns were reconstructed for the same evolutionary lineages as descending dysploidy.

Tracing the Evolution of the Angiosperm Genome from the Cytogenetic Point of View

Cytogenetics constitutes a branch of genetics that is focused on the cellular components, especially chromosomes, in relation to heredity and genome structure, function and evolution. The use of modern cytogenetic approaches and the latest microscopes with image acquisition and processing systems enables the simultaneous two- or three-dimensional, multicolour visualisation of both single-copy and highly-repetitive sequences in the plant genome. The data that is gathered using the cytogenetic methods in the phylogenetic background enable tracing the evolution of the plant genome that involve changes in: (i) genome sizes; (ii) chromosome numbers and morphology; (iii) the content of repetitive sequences and (iv) ploidy level. Modern cytogenetic approaches such as FISH using chromosome- and genome-specific probes have been widely used in studies of the evolution of diploids and the consequences of polyploidy. Nowadays, modern cytogenetics complements analyses in other fields of cell biology and constitutes the linkage between genetics, molecular biology and genomics.

The dynamics of genome size and GC contents evolution in genus Nicotiana

Hybridization and Polyploidization are most common of the phenomenon observed in plants, especially in the genus Nicotiana leading to the duplication of genome. Although genomic changes associated with these events has been studied at various levels but the genome size and GC content variation is less understood because of absence of sufficient genomic data. In this study the flow cytometry technique was used to uncover the genome size and GC contents of 46 Nicotiana species and we compared the genomic changes associated with the hybridization events along evolutionary time scale. The genome size among Nicotiana species varied between 3.28 pg and 11.88 pg whereas GC contents varied between 37.22% and 51.25%. The tetraploid species in genus Nicotiana including section Polydiclae, Repandae, Nicotiana, Rustica and Sauveolentes revealed both up and downsizing in their genome sizes when compared to the sum of genomes of their ancestral species. The genome sizes of three homoploid hybrids were found near their ancestral species. Loss of large genome sequence was observed in the evolutionary more aged species (>10 Myr) as compared to the recently evolved one's (<0.2 Myr). The GC contents were found homogenous with a mean difference of 2.46% among the Nicotiana species. It is concluded that genome size change appeared in either direction whereas the GC contents were found more homogenous in genus Nicotiana.

Origin and Fates of TERT Gene Copies in Polyploid Plants

The gene coding for the telomerase reverse transcriptase (TERT) is essential for the maintenance of telomeres. Previously we described the presence of three TERT paralogs in the allotetraploid plant Nicotiana tabacum, while a single TERT copy was identified in the paleopolyploid model plant Arabidopsis thaliana. Here we examine the presence, origin and functional status of TERT variants in allotetraploid Nicotiana species of diverse evolutionary ages and their parental genome donors, as well as in other diploid and polyploid plant species. A combination of experimental and in silico bottom-up analyses of TERT gene copies in Nicotiana polyploids revealed various patterns of retention or loss of parental TERT variants and divergence in their functions. RT-qPCR results confirmed the expression of all the identified TERT variants. In representative plant and green algal genomes, our synteny analyses show that their TERT genes were located in a conserved locus that became advantageous after the divergence of eudicots, and the gene was later translocated in several plant groups. In various diploid and polyploid species, translocation of TERT became fixed in target loci that show ancient synapomorphy.

Molecular characterization of the 5S rDNA non-transcribed spacer and reconstruction of phylogenetic relationships in Capsicum

Abstract Capsicum includes ca. 41 species of chili peppers. In this original report we PCR amplified, cloned, sequenced and characterized the 5S rDNA non-transcribed spacer -NTS- in 23 taxa of nine clades of Capsicum, divergent at geographical origin and fruit and chromosome traits, and compared the NTS features throughout Solanaceae. According to GC content, inner variability and regulatory elements, the NTS organizes into three distinct structural regions; genetic variability at the NTS in Capsicum and related genus clusters into defined taxa hierarchies. Based on the reconstruction of a maximum-likelihood phylogenetic tree and phylogenetic networks, NTS sequences of Capsicum and related taxa grouped into well recognized categories -genus, section, clade, species, variety-. An evolutionary scenario arose from combined genetic and phylogenetic NTS data, in which monophyly and lineage diversification over time of Capsicum are addressed. Our analysis is original to include all domesticated species of Capsicum prevailing in germplasm collections and breeding programs, together with a large group of wild taxa that demanded further genetic characterization. The NTS set up as a double purpose marker in Capsicum, to directly evaluate genetic variability and reconstruct phylogenetic relationships to a broad extent, and constitutes a valuable tool for germplasm characterization and evolutionary studies within Solanaceae.

Identification and characterization of a tobacco metal tolerance protein, NtMTP2

Metal tolerance proteins (MTPs) from the CDF (Cation Diffusion Facilitator) family are efflux transporters that play a crucial role in metal homeostasis by maintaining optimal metal concentrations in the cytoplasm. Here, a novel tobacco NtMTP2 transporter was cloned and characterized. It encodes a 512 aa protein containing all specific CDF family domains. A GFP-NtMTP2 fusion protein localizes to the tonoplast in tobacco cells. NtMTP2 expression in yeast conferred tolerance to Co and Ni, indicating that the protein mediates transport of both metals, but not Zn, Mn, Cu, Fe, or Cd. Nonetheless, the expression level was not affected by Co or Ni, except for an increase in leaves at high Co concentrations. Its expression in plant parts remained stable during development, but increased in the leaves of older plants. Analysis of tobacco expressing a promoter-GUS construct indicates that the main sites of promoter activity are the conductive tissue throughout the plant and the palisade parenchyma in leaves. Our results suggest that NtMTP2 is a tonoplast transporter mediating sequestration of Co and Ni into vacuoles and an important housekeeping protein that controls the basal availability of micronutrients and plays a role in the sequestration of metal excess, specifically in leaves.

The fate of 35S rRNA genes in the allotetraploid grass Brachypodium hybridum

Nucleolar dominance (ND) consists of the reversible silencing of 35S/45S rDNA loci inherited from one of the ancestors of an allopolyploid. The molecular mechanisms by which one ancestral rDNA set is selected for silencing remain unclear. We applied a combination of molecular (Southern blot hybridization and reverse-transcription cleaved amplified polymorphic sequence analysis), genomic (analysis of variants) and cytogenetic (fluorescence in situ hybridization) approaches to study the structure, expression and epigenetic landscape of 35S rDNA in an allotetraploid grass that exhibits ND, Brachypodium hybridum (genome composition DDSS), and its putative progenitors, Brachypodium distachyon (DD) and Brachypodium stacei (SS). In progenitor genomes, B. stacei showed a higher intragenomic heterogeneity of rDNA compared with B. distachyon. In all studied accessions of B. hybridum, there was a reduction in the copy number of S homoeologues, which was accompanied by their inactive transcriptional status. The involvement of DNA methylation in CG and CHG contexts in the silencing of the S-genome rDNA loci was revealed. In the B. hybridum allotetraploid, ND is stabilized towards the D-genome units, irrespective of the polyphyletic origin of the species, and does not seem to be influenced by homoeologous 35S rDNA ratios and developmental stage.

Plastid genomics of Nicotiana (Solanaceae): insights into molecular evolution, positive selection and the origin of the maternal genome of Aztec tobacco (Nicotiana rustica)

Species of the genus Nicotiana (Solanaceae), commonly referred to as tobacco plants, are often cultivated as non-food crops and garden ornamentals. In addition to the worldwide production of tobacco leaves, they are also used as evolutionary model systems due to their complex development history tangled by polyploidy and hybridization. Here, we assembled the plastid genomes of five tobacco species: N. knightiana, N. rustica, N. paniculata, N. obtusifolia and N. glauca. De novo assembled tobacco plastid genomes had the typical quadripartite structure, consisting of a pair of inverted repeat (IR) regions (25,323-25,369 bp each) separated by a large single-copy (LSC) region (86,510-86,716 bp) and a small single-copy (SSC) region (18,441-18,555 bp). Comparative analyses of Nicotiana plastid genomes with currently available Solanaceae genome sequences showed similar GC and gene content, codon usage, simple sequence and oligonucleotide repeats, RNA editing sites, and substitutions. We identified 20 highly polymorphic regions, mostly belonging to intergenic spacer regions (IGS), which could be suitable for the development of robust and cost-effective markers for inferring the phylogeny of the genus Nicotiana and family Solanaceae. Our comparative plastid genome analysis revealed that the maternal parent of the tetraploid N. rustica was the common ancestor of N. paniculata and N. knightiana, and the later species is more closely related to N. rustica. Relaxed molecular clock analyses estimated the speciation event between N. rustica and N. knightiana appeared 0.56 Ma (HPD 0.65-0.46). Biogeographical analysis supported a south-to-north range expansion and diversification for N. rustica and related species, where N. undulata and N. paniculata evolved in North/Central Peru, while N. rustica developed in Southern Peru and separated from N. knightiana, which adapted to the Southern coastal climatic regimes. We further inspected selective pressure on protein-coding genes among tobacco species to determine if this adaptation process affected the evolution of plastid genes. These analyses indicate that four genes involved in different plastid functions, including DNA replication (rpoA) and photosynthesis (atpB, ndhD and ndhF), came under positive selective pressure as a result of specific environmental conditions. Genetic mutations in these genes might have contributed to better survival and superior adaptations during the evolutionary history of tobacco species.

Expression patterns of NbrgsCaM family genes in Nicotiana benthamiana and their potential roles in development and stress responses

rgsCaM has been reported as a calmodulin-like (CML) factor induced by viral infection in Nicotiana. There are three CMLs that belong to the rgsCaM family in Arabidopsis thaliana. In this study, we found a total of 5 NbrgsCaM coding sequences in N. benthamiana genome. We analyzed transcription patterns of NbrgsCaMs in transgenic plants expressing a β-glucuronidase (GUS) under the promoter of NbrgsCaMs by histochemistry staining and RT-qPCR. Similar to their Arabidopsis homologs, most NbrgsCaMs have an overlapping but distinct expression pattern in response to developmental and environmental changes. Specifically, the NbrgsCaM4 promoter exhibited robust activity and showed distinct regulatory response to viral infection, developmental stages and other abiotic stimuli. Overall, these findings provide clues for further understanding of the NbrgsCaM family genes in regulating plant growth and development under biotic stress and environmental stimulation.

Genome-wide association study of leaf chemistry traits in tobacco

Leaf chemistry traits are some of the key factors influencing tobacco quality, which can be significantly reduced by lower chemical components in cured leaf. To improve tobacco quality through breeding, genetic diversity analysis, population structure analysis, and genome-wide association studies were performed in a panel of 347 tobacco germplasms and the markers associated with five leaf chemistry traits, including total sugar (TS), reducing sugar (RS), total nitrogen (TN), nicotine (NIC), and total potassium (TP) contents were identified. Four groups were classified at a genetic distance of 0.316 by genetic diversity analysis based on coefficient parameter NEI72 using a program NTSYS-pc2.10e, whereas four well-differentiated subpopulations were postulated in the 347 tobacco accessions. A total of 47 target trait-associated SNPs was detected in at least three environments as well as the best linear unbiased predictions (BLUPs) across all environments, among which two, two, four, six, and one highly suggestive associated SNPs were repeatedly detected in all environments and BLUPs for TS, RS, TN, NIC, and TP, respectively. On the basis of the phenotypic effects of the alleles corresponding to suggestive associated SNPs, five tobacco accessions harboring favorable alleles with elite phenotypic performance in leaf chemistry traits were identified. The results could facilitate quality tobacco breeding for higher leaf chemistry trait contents through molecular marker-assisted approaches.

Reconstructing phylogenetic relationships based on repeat sequence similarities

A recent phylogenetic method based on genome-wide abundance of different repeat types proved to be useful in reconstructing the evolutionary history of several plant and animal groups. Here, we demonstrate that an alternative information source from the repeatome can also be employed to infer phylogenetic relationships among taxa. Specifically, this novel approach makes use of the repeat sequence similarity matrices obtained from the comparative clustering analyses of RepeatExplorer 2, which are subsequently transformed to between-taxa distance matrices. These pairwise matrices are used to construct neighbour-joining trees for each of the top most-abundant clusters and they are finally summarized in a consensus network. This methodology was tested on three groups of angiosperms and one group of insects, resulting in congruent evolutionary hypotheses compared to more standard systematic analyses based on commonly used DNA markers. We propose that the combined application of these phylogenetic approaches based on repeat abundances and repeat sequence similarities could be helpful to understand mechanisms governing genome and repeatome evolution.

Parental origin of the allotetraploid tobacco Nicotiana benthamiana

Nicotiana section Suaveolentes is an almost all-Australian clade of allopolyploid tobacco species including the important plant model Nicotiana benthamiana. The homology relationships of this clade and its formation are not completely understood. To address this gap, we assessed phylogenies of all individual genes of N. benthamiana and the well studied N. tabacum (section Nicotiana) and their homologues in six diploid Nicotiana species. We generated sets of 44 424 and 65 457 phylogenetic trees of N. benthamiana and N. tabacum genes, respectively, each collectively called a phylome. Members of Nicotiana sections Noctiflorae and Sylvestres were represented as the species closest to N. benthamiana in most of the gene trees. Analyzing the gene trees of the phylome we: (i) dated the hybridization event giving rise to N. benthamiana to 4-5 MyA, and (ii) separated the subgenomes. We assigned 1.42 Gbp of the genome sequence to section Noctiflorae and 0.97 Gbp to section Sylvestres based on phylome analysis. In contrast, read mapping of the donor species did not succeed in separating the subgenomes of N. benthamiana. We show that the maternal progenitor of N. benthamiana was a member of section Noctiflorae, and confirm a member of section Sylvestres as paternal subgenome donor. We also demonstrate that the advanced stage of long-term genome diploidization in N. benthamiana is reflected in its subgenome organization. Taken together, our results underscore the usefulness of phylome analysis for subgenome characterization in hybrid species.

Construction of a high-density genetic map with whole genome sequencing in Nicotiana tabacum L

Tobacco (Nicotiana tabacum L.) is an essential commercial crop and an ideal model plant for biological mechanism studies. As an allopolyploid species, tobacco harbors a massive and complex genome, which makes the application of molecular markers complicated and challenging. In our study, we performed whole-genome sequencing of an intraspecific recombinant inbred line (RIL) population, a F1 generation and their parents. With the Nicotiana tabacum (K326 cultivar) genome as reference, a total of 45,081 markers were characterized to construct the genetic map, which spanned a genetic distance of 3486.78 cM. Evaluation of a two-dimensional heat map proved the high quality of the genetic map. We utilized these markers to anchor scaffolds and analyzed the ancestral genome origin of linkage groups (LGs). Furthermore, such a high-density genetic map will be applied for quantitative trait locus (QTL) detection, gene localization, genome-wide association studies (GWAS), and marker-assisted breeding in tobacco.

Repetitive DNA Restructuring Across Multiple Nicotiana Allopolyploidisation Events Shows a Lack of Strong Cytoplasmic Bias in Influencing Repeat Turnover

Allopolyploidy is acknowledged as an important force in plant evolution. Frequent allopolyploidy in Nicotiana across different timescales permits the evaluation of genome restructuring and repeat dynamics through time. Here we use a clustering approach on high-throughput sequence reads to identify the main classes of repetitive elements following three allotetraploid events, and how these are inherited from the closest extant relatives of the maternal and paternal subgenome donors. In all three cases, there was a lack of clear maternal, cytoplasmic bias in repeat evolution, i.e., lack of a predicted bias towards maternal subgenome-derived repeats, with roughly equal contributions from both parental subgenomes. Different overall repeat dynamics were found across timescales of <0.5 (N. rustica L.), 4 (N. repanda Willd.) and 6 (N. benthamiana Domin) Ma, with nearly additive, genome upsizing, and genome downsizing, respectively. Lower copy repeats were inherited in similar abundance to the parental subgenomes, whereas higher copy repeats contributed the most to genome size change in N. repanda and N. benthamiana. Genome downsizing post-polyploidisation may be a general long-term trend across angiosperms, but at more recent timescales there is species-specific variance as found in Nicotiana.

Divergence of duplicated genes by repeated partitioning of splice forms and subcellular localization

Gene duplication is a prominent and recurrent process in plant genomes. Among the possible fates of duplicated genes, subfunctionalization refers to duplicates taking on different parts of the function or expression pattern of the ancestral gene. This partitioning could be accompanied by changes in subcellular localization of the protein products. When alternative splicing of gene products leads to protein products with different subcellular localizations, we propose that after gene duplication there will be partitioning of the alternatively spliced forms such that the products of each duplicate are localized to only one of the original locations, which we refer to as sublocalization. We identified the plastid ascorbate peroxidase (cpAPX) genes across angiosperms and analyzed their duplication history, alternative splicing, and subcellular targeting patterns to identify cases of sublocalization. We found angiosperms typically have one cpAPX gene that generates both thylakoidal APX (tAPX) and stromal APX (sAPX) through alternative splicing. We identified several independent lineage-specific sublocalization cases with specialized paralogues of tAPX and sAPX. We determined that the sublocalization happened through two types of sequence evolution patterns. Our findings suggest that the divergence through sublocalization is key to the retention of paralogous cpAPX genes in angiosperms.

Evolutionary history of the Pasque-flowers (Pulsatilla, Ranunculaceae): Molecular phylogenetics, systematics and rDNA evolution

Pulsatilla (Anemoneae, Ranunculaceae) is sister to Anemone s.s. and contains ca 40 perennial species of considerable horticultural and medical importance. We sequenced 31 of those species, plus nine subspecies, two cultivars and six outgroups, for two nuclear regions (high-copy nrITS and low-copy MLH1) and three plastid regions (rbcL, accD-psaI, trnL intron) in order to generate the first comprehensive species-level phylogeny of the genus. Phylogenetic trees were constructed using both concatenation-based (maximum likelihood and Bayesian inference) and coalescence methods. The better supported among the internal nodes were subjected to molecular clock dating and ancestral area reconstruction, and karyotypic characters identified by us using Fluorescence In Situ Hybridization were mapped across the tree. The preferred species tree from the coalescence analysis formed the basis of a new infrageneric classification based on monophyly plus degree of divergence. The earliest divergent of the three subgenera, Kostyczewianae, is represented by only a single species that is morphologically intermediate between Anemone s.s. and 'core' Pulsatilla. Subgenus Pulsatilla is considerably richer in species than its sister subgenus Preonanthus and contains three monophyletic sections. Species possessing nodding flowers and pectinately dissected leaves are phylogenetically derived compared with groups possessing erect flowers and palmately lobed leaves. Pulsatilla separated from Anemone s.s. at ca 25 Ma. Our results indicate a central Asian mountain origin of the genus and an initial diversification correlated with late Tertiary global cooling plus regional mountain uplift, aridification and consequent expansion of grasslands. The more rapid and extensive diversification within subgenus Pulsatilla began at ca 3 Ma and continued throughout the Quaternary, driven not only by major perturbations in global climate but also by well-documented polyploidy.

Change in nuclear DNA content and pollen size with polyploidisation in the sweet potato (Ipomoea batatas, Convolvulaceae) complex

Genome size evolution and its relationship with pollen grain size has been investigated in sweet potato (Ipomoea batatas), an economically important crop which is closely related to diploid and tetraploid species, assessing the nuclear DNA content of 22 accessions from five Ipomoea species, ten sweet potato varieties and two outgroup taxa. Nuclear DNA amounts were determined using flow cytometry. Pollen grains were studied using scanning and transmission electron microscopy. 2C DNA content of hexaploid I. batatas ranged between 3.12-3.29 pg; the mean monoploid genome size being 0.539 pg (527 Mbp), similar to the related diploid accessions. In tetraploid species I. trifida and I. tabascana, 2C DNA content was, respectively, 2.07 and 2.03 pg. In the diploid species closely related to sweet potato e.g. I. ×leucantha, I. tiliacea, I. trifida and I. triloba, 2C DNA content was 1.01-1.12 pg. However, two diploid outgroup species, I. setosa and I. purpurea, were clearly different from the other diploid species, with 2C of 1.47-1.49 pg; they also have larger chromosomes. The I. batatas genome presents 60.0% AT bases. DNA content and ploidy level were positively correlated within this complex. In I. batatas and the more closely related species I. trifida, the genome size and ploidy levels were correlated with pollen size. Our results allow us to propose alternative or complementary hypotheses to that currently proposed for the formation of hexaploid Ipomoea batatas.

Ribosomal DNA loci derived from Brachypodium stacei are switched off for major parts of the life cycle of Brachypodium hybridum

Nucleolar dominance is an epigenetic phenomenon that occurs in some plant and animal allopolyploids and hybrids, whereby only one ancestral set of 35S rRNA genes retains the ability to form the nucleolus while the rDNA loci derived from the other progenitor are transcriptionally silenced. There is substantial evidence that nucleolar dominance is regulated developmentally. This study focuses upon the establishment and/or maintenance of nucleolar dominance during different stages of development in the model grass allotetraploid Brachypodium hybridum. Fluorescence in situ hybridization with a 25S rDNA probe to cells in three-dimensional cytogenetic preparations showed that nucleolar dominance is present not only in root meristematic and differentiated cells of this species, but also in male meiocytes at prophase I, tetrads of microspores, and different embryonic tissues. The inactive state of Brachypodium stacei-originated rDNA loci was confirmed by silver staining. Only B. distachyon-derived 35S rDNA loci formed nucleoli in the aforementioned tissues, whereas B. stacei-like loci remained highly condensed and thus transcriptionally suppressed. The establishment of nucleolar dominance during earlier stages of B. hybridum embryo development cannot be ruled out. However, we propose that gradual pseudogenization of B. stacei-like loci in the evolution of the allotetraploid seems to be more likely.

The impact of genome evolution on the allotetraploid Nicotiana rustica - an intriguing story of enhanced alkaloid production

BACKGROUND

Nicotiana rustica (Aztec tobacco), like common tobacco (Nicotiana tabacum), is an allotetraploid formed through a recent hybridization event; however, it originated from completely different progenitor species. Here, we report the comparative genome analysis of wild type N. rustica (5 Gb; 2n = 4x = 48) with its three putative diploid progenitors (2.3-3 Gb; 2n = 2x =24), Nicotiana undulata, Nicotiana paniculata and Nicotiana knightiana.

RESULTS

In total, 41% of N. rustica genome originated from the paternal donor (N. undulata), while 59% originated from the maternal donor (N. paniculata/N. knightiana). Chloroplast genome and gene analyses indicated that N. knightiana is more closely related to N. rustica than N. paniculata. Gene clustering revealed 14,623 ortholog groups common to other Nicotiana species and 207 unique to N. rustica. Genome sequence analysis indicated that N. knightiana is more closely related to N. rustica than N. paniculata, and that the higher nicotine content of N. rustica leaves is the result of the progenitor genomes combination and of a more active transport of nicotine to the shoot.

CONCLUSIONS

The availability of four new Nicotiana genome sequences provide insights into how speciation impacts plant metabolism, and in particular alkaloid transport and accumulation, and will contribute to better understanding the evolution of Nicotiana species.

Satellite DNA in Paphiopedilum subgenus Parvisepalum as revealed by high-throughput sequencing and fluorescent in situ hybridization

BACKGROUND

Satellite DNA is a rapidly diverging, largely repetitive DNA component of many eukaryotic genomes. Here we analyse the evolutionary dynamics of a satellite DNA repeat in the genomes of a group of Asian subtropical lady slipper orchids (Paphiopedilum subgenus Parvisepalum and representative species in the other subgenera/sections across the genus). A new satellite repeat in Paphiopedilum subgenus Parvisepalum, SatA, was identified and characterized using the RepeatExplorer pipeline in HiSeq Illumina reads from P. armeniacum (2n = 26). Reconstructed monomers were used to design a satellite-specific fluorescent in situ hybridization (FISH) probe. The data were also analysed within a phylogenetic framework built using the internal transcribed spacer (ITS) sequences of 45S nuclear ribosomal DNA.

RESULTS

SatA comprises c. 14.5% of the P. armeniacum genome and is specific to subgenus Parvisepalum. It is composed of four primary monomers that range from 230 to 359 bp and contains multiple inverted repeat regions with hairpin-loop motifs. A new karyotype of P. vietnamense (2n = 28) is presented and shows that the chromosome number in subgenus Parvisepalum is not conserved at 2n = 26, as previously reported. The physical locations of SatA sequences were visualised on the chromosomes of all seven Paphiopedilum species of subgenus Parvisepalum (2n = 26-28), together with the 5S and 45S rDNA loci using FISH. The SatA repeats were predominantly localisedin the centromeric, peri-centromeric and sub-telocentric chromosome regions, but the exact distribution pattern was species-specific.

CONCLUSIONS

We conclude that the newly discovered, highly abundant and rapidly evolving satellite sequence SatA is specific to Paphiopedilum subgenus Parvisepalum. SatA and rDNA chromosomal distributions are characteristic of species, and comparisons between species reveal that the distribution patterns generate a strong phylogenetic signal. We also conclude that the ancestral chromosome number of subgenus Parvisepalum and indeed of all Paphiopedilum could be either 2n = 26 or 28, if P. vietnamense is sister to all species in the subgenus as suggested by the ITS data.

A revision of the Old World Black Nightshades (Morelloid clade of Solanum L., Solanaceae)

The Morelloid clade, also known as the black nightshades or "Maurella" (Morella), is one of the 10 major clades within Solanum L. The pantropical clade consists of 75 currently recognised non-spiny herbaceous and suffrutescent species with simple or branched hairs with or without glandular tips, with a centre of distribution in the tropical Andes. A secondary centre of diversity is found in Africa, where a set of mainly polyploid taxa occur. A yet smaller set of species is found in Australasia and Europe, including Solanum nigrum L., the type of the genus Solanum. Due to the large number of published synonyms, combined with complex morphological variation, our understanding of species limits and diversity in the Morelloid clade has remained poor despite detailed morphological studies carried out in conjunction with breeding experiments. Here we provide the first taxonomic overview since the 19th century of the entire group in the Old World, including Africa, Asia, Australia, Europe and islands of the Pacific. Complete synonymy, morphological descriptions, distribution maps and common names and uses are provided for all 19 species occurring outside the Americas (i.e. Africa, Asia, Australia, Europe and islands of the Pacific). We treat 12 species native to the Old World, as well as 7 taxa that are putatively introduced and/or invasive in the region. The current knowledge of the origin of the polyploid species is summarised. A key to all of the species occurring in the Old World is provided, together with line drawings and colour figures to aid identification both in herbaria and in the field. Preliminary conservation assessments are provided for all species.

NtTPN1: a RPP8-like R gene required for Potato virus Y-induced veinal necrosis in tobacco

Potato virus Y (PVY) is one of the most damaging viruses of tobacco. In particular, aggressive necrotic strains (PVYN ) lead to considerable losses in yield. The main source of resistance against PVY is linked to the va locus. However, va-overcoming PVY isolates inducing necrotic symptoms were observed in several countries. In this context, it is important to find va-independent protection strategies. In a previous study, the phenotyping of 162 tobacco varieties revealed 10 accessions that do not carry the va allele and do not exhibit typical PVYN -induced veinal necrosis. Despite the absence of necrotic symptoms, normal viral accumulation in these plants suggests a va-independent mechanism of tolerance to PVYN -induced systemic veinal necrosis. Fine mapping of the genetic determinant(s) was performed in a segregating F2 population. The tolerance trait is inherited as a single recessive gene, and allelism tests demonstrated that eight of the 10 tolerant varieties carry the same determinant. Anchoring the linkage map to the tobacco genome physical map allowed the identification of a RPP8-like R gene, called NtTPN1 (for Nicotiana tabacum Tolerance to PVY-induced Necrosis1), with the same single-nucleotide polymorphism in the eight tolerant accessions. Functional assays using homozygous NtTPN1 EMS mutants confirmed the role of NtTPN1 in the tolerance phenotype. PVYN -induced systemic veinal necrosis in tobacco likely represents an inefficient defense response with hypersensitive response-like characteristics. The identification of NtTPN1 opens breeding options to minimize the impact of emerging and so far uncontrolled va-breaking necrotic PVY isolates.

Agrobacterium-Mediated Transformation in the Evolution of Plants

In most cases, the genetic engineering of plants uses Agrobacterium-mediated transformation to introduce novel genes. In nature, insertion of T-DNA into the plant genome and its subsequent transfer via sexual reproduction have been shown for several species in the genera Nicotiana, Ipomoea , and Linaria . A sequence homologous to T-DNA of the Ri plasmid of Agrobacterium rhizogenes was found in the genome of wild-type Nicotiana glauca (section Noctiflorae) more than 30 years ago and was named "cellular T-DNA" (cT-DNA). It comprises an imperfect inverted repeat and contains homologs of several T-DNA oncogenes (NgrolB, NgrolC, Ngorf13, Ngorf14) and an opine synthesis gene (Ngmis). Multiple cT-DNAs have also been found in species of the sections Tomentosae and Nicotiana of the genus Nicotiana. These ancient cT-DNA genes are still expressed, indicating that they may play a role in the evolution of these plants. In 2012-2013, cT-DNA was detected and characterized in Linaria vulgaris and L. genistifolia ssp. dalmatica. Their cT-DNA is present in two copies and organized as an imperfect direct tandem repeat, containing LvORF2, LvORF3, LvORF8, LvrolA, LvrolB, LvrolC, LvORF13, LvORF14, and the Lvmis genes. In 2015, cT-DNA was found in Ipomoea. Two types of T-DNA-like sequences were described within this genera, and their distribution varied among cultured hexaploid, tetraploid, and wild diploid forms. Thus, several independent T-DNA integration events occurred in the genomes of these three plant genera. We propose that the events of T-DNA insertion in the plant genome might have affected their evolution, resulting in the creation of new plant species. In this chapter, we focus on the structure and functions of cT-DNA in Linaria, Nicotiana, and Ipomoea and discuss their possible evolutionary role.

Functional Analysis of NtZIP4B and Zn Status-Dependent Expression Pattern of Tobacco ZIP Genes

Tobacco is frequently considered as a plant useful for phytoremediation of metal-contaminated soil, despite the mechanisms for regulation of uptake and accumulation being largely unknown. Here we cloned and characterized a new tobacco Zn and Cd transporter NtZIP4B from the ZIP family (ZRT-IRT-Like proteins). It complemented the Zn-uptake defective yeast mutant zrt1zrt2, and rendered the wild type DY1457 yeast more sensitive to Cd. Bioinformatic analysis and transient expression of the NtZIP4B-GFP fusion protein in tobacco leaves indicated its localization to the plasma membrane. Real-time q-PCR based analysis showed that it is expressed in all vegetative organs with the highest level in leaves. The Zn status determined transcript abundance; NtZIP4B was upregulated by Zn-deficiency and downregulated by Zn excess. At the tissue level, in roots NtZIP4B is expressed in the vasculature of the middle part of the roots and in surrounding tissues including the root epidermis; in leaves primarily in the vasculature. Bioinformatic analysis identified two copies of ZIP4 in tobacco, NtZIP4A and NtZIP4B with 97.57% homology at the amino acid level, with the same expression pattern for both, indicating a high degree of functional redundancy. Moreover, the present study provides new insights into the coordinated function of NtZIP1, NtZIP2, NtZIP4, NtZIP5, NtZIP8, NtIRT1, and NtIRT1-like in response to low-to-high Zn status. Leaves were the major site of NtZIP4, NtZIP5, and NtZIP8 expression, and roots for NtZIP1, NtZIP2, NtIRT1, and NtIRT1-like. Contrasting expression level in the apical and basal root parts indicates distinct roles in root-specific processes likely contributing to the regulation of Zn root-to-shoot translocation. In summary, new insight into the role of ZIP genes in Zn homeostasis pointing to their overlapping and complementary functions, offers opportunities for strategies to modify Zn and Cd root/shoot partition in tobacco.

Ecological Roles and Biological Activities of Specialized Metabolites from the Genus Nicotiana

Species of Nicotiana grow naturally in different parts of the world and have long been used both medicinally and recreationally by human societies. More recently in our history, Nicotiana tabacum has attracted interest as one of the most economically important industrial crops. Nicotiana species are frequently investigated for their bioactive natural products, and the ecological role of their specialized metabolites in responses to abiotic stress or biotic stress factors like pathogens and herbivores. The interest of tobacco companies in genetic information as well as the success of a few wild tobacco species as experimental model organisms have resulted in growing knowledge about the molecular biology and ecology of these plants and functional studies of the plant's natural products. Although a large number of reviews and books on biologically active natural products already exists, mostly from N. tabacum, we focus our attention on the ecological roles and biological activity of natural products, versus products from cured and processed material, in this Review. The studied compounds include alkaloids, aromatic compounds, flavonoids, volatiles, sesquiterpenoids, diterpenes alcohols, and sugar esters from trichomes of the plants, and recently characterized acyclic hydroxygeranyllinalool diterpene glycosides (HGL-DTGs). In this Review (1800s-2017), we describe the above-mentioned classes of natural products, emphasizing their biological activities and functions as they have been determined either in bioassay-guided purification approaches or in bioassays with plants in which the expression of specific biosynthetic genes has been genetically manipulated. Additionally, a review on the history, taxonomy, ecology, and medicinal application of different Nicotiana species growing around the globe presented in this Review may be of interest for pharmacognosists, natural products, and ecological chemists.

Evolutional dynamics of 45S and 5S ribosomal DNA in ancient allohexaploid Atropa belladonna

BACKGROUND

Polyploid hybrids represent a rich natural resource to study molecular evolution of plant genes and genomes. Here, we applied a combination of karyological and molecular methods to investigate chromosomal structure, molecular organization and evolution of ribosomal DNA (rDNA) in nightshade, Atropa belladonna (fam. Solanaceae), one of the oldest known allohexaploids among flowering plants. Because of their abundance and specific molecular organization (evolutionarily conserved coding regions linked to variable intergenic spacers, IGS), 45S and 5S rDNA are widely used in plant taxonomic and evolutionary studies.

RESULTS

Molecular cloning and nucleotide sequencing of A. belladonna 45S rDNA repeats revealed a general structure characteristic of other Solanaceae species, and a very high sequence similarity of two length variants, with the only difference in number of short IGS subrepeats. These results combined with the detection of three pairs of 45S rDNA loci on separate chromosomes, presumably inherited from both tetraploid and diploid ancestor species, example intensive sequence homogenization that led to substitution/elimination of rDNA repeats of one parent. Chromosome silver-staining revealed that only four out of six 45S rDNA sites are frequently transcriptionally active, demonstrating nucleolar dominance. For 5S rDNA, three size variants of repeats were detected, with the major class represented by repeats containing all functional IGS elements required for transcription, the intermediate size repeats containing partially deleted IGS sequences, and the short 5S repeats containing severe defects both in the IGS and coding sequences. While shorter variants demonstrate increased rate of based substitution, probably in their transition into pseudogenes, the functional 5S rDNA variants are nearly identical at the sequence level, pointing to their origin from a single parental species. Localization of the 5S rDNA genes on two chromosome pairs further supports uniparental inheritance from the tetraploid progenitor.

CONCLUSIONS

The obtained molecular, cytogenetic and phylogenetic data demonstrate complex evolutionary dynamics of rDNA loci in allohexaploid species of Atropa belladonna. The high level of sequence unification revealed in 45S and 5S rDNA loci of this ancient hybrid species have been seemingly achieved by different molecular mechanisms.

Natural Agrobacterium Transformants: Recent Results and Some Theoretical Considerations

Agrobacterium rhizogenes causes hairy root growth on a large number of plant species. It does so by transferring specific DNA fragments (T-DNA) from its root-inducing plasmid (pRi) into plant cells. Expression of T-DNA genes leads to abnormal root growth and production of specific metabolites (opines) which are taken up by the bacterium and used for its growth. Recent work has shown that several Nicotiana, Linaria, and Ipomoea species contain T-DNA genes from A. rhizogenes in their genomes. Plants carrying such T-DNAs (called cellular T-DNA or cT-DNA) can be considered as natural transformants. In the Nicotiana genus, seven different T-DNAs are found originating from different Agrobacterium strains, and in the Tomentosae section no <4 successive insertion events took place. In several cases cT-DNA genes were found to be expressed. In some Nicotiana tabacum cultivars the opine synthesis gene TB-mas2' is expressed in the roots. These cultivars were found to produce opines. Here we review what is known about natural Agrobacterium transformants, develop a theoretical framework to analyze this unusual phenomenon, and provide some outlines for further research.

Related allopolyploids display distinct floral pigment profiles and transgressive pigments

PREMISE OF THE STUDY

Both polyploidy and shifts in floral color have marked angiosperm evolution. Here, we investigate the biochemical basis of the novel and diverse floral phenotypes seen in allopolyploids in Nicotiana (Solanaceae) and examine the extent to which the merging of distinct genomes alters flavonoid pigment production.

METHODS

We analyzed flavonol and anthocyanin pigments from Nicotiana allopolyploids of different ages (N. tabacum, 0.2 million years old; several species from Nicotiana section Repandae, 4.5 million years old; and five lines of first-generation synthetic N. tabacum) as well as their diploid progenitors.

KEY RESULTS

Allopolyploid floral pigment profiles tend not to overlap with their progenitors or related allopolyploids, and allopolyploids produce transgressive pigments that are not present in either progenitor. Differences in floral color among N. tabacum accessions seems mainly to be due to variation in cyanidin concentration, but changes in flavonol concentrations among accessions are also present.

CONCLUSIONS

Competition for substrates within the flavonoid biosynthetic pathway to make either flavonols or anthocyanins may drive the differences seen among related allopolyploids. Some of the pigment differences observed in allopolyploids may be associated with making flowers more visible to nocturnal pollinators.

Genome-wide repeat dynamics reflect phylogenetic distance in closely related allotetraploid Nicotiana (Solanaceae)

Nicotiana sect. Repandae is a group of four allotetraploid species originating from a single allopolyploidisation event approximately 5 million years ago. Previous phylogenetic analyses support the hypothesis of N. nudicaulis as sister to the other three species. This is concordant with changes in genome size, separating those with genome downsizing (N. nudicaulis) from those with genome upsizing (N. repanda, N. nesophila, N. stocktonii). However, a recent analysis reflecting genome dynamics of different transposable element families reconstructed greater similarity between N. nudicaulis and the Revillagigedo Island taxa (N. nesophila and N. stocktonii), thereby placing N. repanda as sister to the rest of the group. This could reflect a different phylogenetic hypothesis or the unique evolutionary history of these particular elements. Here we re-examine relationships in this group and investigate genome-wide patterns in repetitive DNA, utilising high-throughput sequencing and a genome skimming approach. Repetitive DNA clusters provide support for N. nudicaulis as sister to the rest of the section, with N. repanda sister to the two Revillagigedo Island species. Clade-specific patterns in the occurrence and abundance of particular repeats confirm the original (N. nudicaulis (N. repanda (N. nesophila + N. stocktonii))) hypothesis. Furthermore, overall repeat dynamics in the island species N. nesophila and N. stocktonii confirm their similarity to N. repanda and the distinctive patterns between these three species and N. nudicaulis. Together these results suggest that broad-scale repeat dynamics do in fact reflect evolutionary history and could be predicted based on phylogenetic distance.

rDNA Loci Evolution in the Genus Glechoma (Lamiaceae)

Glechoma L. (Lamiaceae) is distributed in eastern Asia and Europe. Understanding chromosome evolution in Glechoma has been strongly hampered by its small chromosomes, constant karyotype and polyploidy. Here phylogenetic patterns and chromosomal variation in Glechoma species are considered, using genome sizes, chromosome mapping of 5S and 35S rDNAs by fluorescence in situ hybridisation (FISH), and phylogenetic analyses of internal transcribed spacers (nrITS) of 35S rDNA and 5S rDNA NTS sequences. Species and populations of Glechoma are tetraploid (2n = 36) with base chromosome number of x = 9. Four chromosomes carry pericentric 5S rDNA sites in their short arms in all the species. Two to four of these chromosomes also carry 35S rDNA in subterminal regions of the same arms. Two to four other chromosomes have 35S rDNA sites, all located subterminally within short arms; one individual possessed additional weak pericentric 35S rDNA signals on three other chromosomes. Five types of rDNA locus distribution have been defined on the basis of 35S rDNA variation, but none is species-specific, and most species have more than one type. Glechoma hederacea has four types. Genome size in Glechoma ranges from 0.80 to 0.94 pg (1C), with low levels of intrapopulational variation in all species. Phylogenetic analyses of ITS and NTS sequences distinguish three main clades coinciding with geographical distribution: European (G. hederacea–G. hirsuta), Chinese and Korean (G. longituba), and Japanese (G. grandis). The paper presents the first comparative cytogenetic analyses of Glechoma species including karyotype structure, rDNA location and number, and genome size interpreted in a phylogenetic context. The observed variation suggests that the genus is still in genomic flux. Genome size, but not rDNA loci number and distribution, provides a character for species delimitation which allows better inferences of interspecific relationships to be made, in the absence of well-defined morphological differentiation.

Unlocking the Karyological and Cytogenetic Diversity of Iris from Lebanon: Oncocyclus Section Shows a Distinctive Profile and Relative Stasis during Its Continental Radiation

Despite being an important target of conservation concern and horticultural interest, Lebanese irises yet have a confusing taxonomic history and species’ delimitation is often considered problematic, more especially among royal irises (Iris section Oncocyclus). Indeed, these irises of exceptionally large and spectacular flowers have radiated across Caucasus and eastern Mediterranean giving rise to a number of strict endemic taxa, many of them being considered under threat. Whilst efforts have mostly focused on clarifying the evolutionary relationships in the group based on morphological and molecular data, karyological and cytogenetic characters have been comparatively overlooked. In this study, we established for the first time the physical mapping of 35S rDNA loci and heterochromatin, and obtained karyo-morphological data for ten Lebanese Iris species belonging to four sections (Iris, Limniris, Oncocyclus and Scorpiris). Our results evidenced distinctive genomic profiles for each one of the sections, where Oncocyclus irises, while having the lowest chromosome numbers, exhibit both the highest number of 35S loci and CMA3+ sites. The continental radiation of royal irises has been accompanied by a relative karyological and cytogenetic stasis, even though some changes were observed regarding karyotype formula and asymmetry indexes. In addition to that, our results enabled taxonomic differentiation between I. germanica and I. mesopotamica–two taxa currently considered as synonyms–and highlighted the need for further studies on populations of I. persica and I. wallasiae in the Eastern Mediterranean Region.

Transgressive phenotypes and generalist pollination in the floral evolution of Nicotiana polyploids

Polyploidy is an important driving force in angiosperm evolution, and much research has focused on genetic, epigenetic and transcriptomic responses to allopolyploidy. Nicotiana is an excellent system in which to study allopolyploidy because half of the species are allotetraploids of different ages, allowing us to examine the trajectory of floral evolution over time. Here, we study the effects of allopolyploidy on floral morphology in Nicotiana, using corolla tube measurements and geometric morphometrics to quantify petal shape. We show that polyploid morphological divergence from the intermediate phenotype expected (based on progenitor morphology) increases with time for floral limb shape and tube length, and that most polyploids are distinct or transgressive in at least one trait. In addition, we show that polyploids tend to evolve shorter and wider corolla tubes, suggesting that allopolyploidy could provide an escape from specialist pollination via reversion to more generalist pollination strategies.

Cytomolecular Analysis of Ribosomal DNA Evolution in a Natural Allotetraploid Brachypodium hybridum and Its Putative Ancestors-Dissecting Complex Repetitive Structure of Intergenic Spacers

Nucleolar dominance is an epigenetic phenomenon associated with nuclear 35S rRNA genes and consists in selective suppression of gene loci inherited from one of the progenitors in the allopolyploid. Our understanding of the exact mechanisms that determine this process is still fragmentary, especially in case of the grass species. This study aimed to shed some light on the molecular basis of this genome-specific inactivation of 35S rDNA loci in an allotetraploid Brachypodium hybridum (2n = 30), which arose from the interspecific hybridization between two diploid ancestors that were very similar to modern B. distachyon (2n = 10) and B. stacei (2n = 20). Using fluorescence in situ hybridization with 25S rDNA and chromosome-specific BAC clones as probes we revealed that the nucleolar dominance is present not only in meristematic root-tip cells but also in differentiated cell fraction of B. hybridum. Additionally, the intergenic spacers (IGSs) from both of the putative ancestors and the allotetraploid were sequenced and analyzed. The presumptive transcription initiation sites, spacer promoters and repeated elements were identified within the IGSs. Two different length variants, 2.3 and 3.5 kb, of IGSs were identified in B. distachyon and B. stacei, respectively, however only the IGS that had originated from B. distachyon-like ancestor was present in the allotetraploid. The amplification pattern of B. hybridum IGSs suggests that some genetic changes occurred in inactive B. stacei-like rDNA loci during the evolution of the allotetraploid. We hypothesize that their preferential silencing is an effect of structural changes in the sequence rather than just the result of the sole inactivation at the epigenetic level.

Cytotype diversity and genome size variation in Knautia (Caprifoliaceae, Dipsacoideae)

BACKGROUND

Polyploidisation is one of the most important mechanisms in the evolution of angiosperms. As in many other genera, formation of polyploids has significantly contributed to diversification and radiation of Knautia (Caprifoliaceae, Dipsacoideae). Comprehensive studies of fine- and broad-scale patterns of ploidy and genome size (GS) variation are, however, still limited to relatively few genera and little is known about the geographic distribution of ploidy levels within these genera. Here, we explore ploidy and GS variation in Knautia based on a near-complete taxonomic and comprehensive geographic sampling.

RESULTS

Genome size is a reliable indicator of ploidy level in Knautia, even if monoploid genome downsizing is observed in the polyploid cytotypes. Twenty-four species studied are diploid, 16 tetraploid and two hexaploid, whereas ten species possess two, and two species possess three ploidy levels. Di- and tetraploids are distributed across most of the distribution area of Knautia, while hexaploids were sampled in the Balkan and Iberian Peninsulas and the Alps.

CONCLUSIONS

We show that the frequency of polyploidisation is unevenly distributed in Knautia both in a geographic and phylogenetic context. Monoploid GS varies considerably among three evolutionary lineages (sections) of Knautia, but also within sections Trichera and Tricheroides, as well as within some of the species. Although the exact causes of this variation remain elusive, we demonstrate that monoploid GS increases significantly towards the limits of the genus' distribution.

The effect of polyploidy and hybridization on the evolution of floral colour in Nicotiana (Solanaceae)

BACKGROUND AND AIMS

Speciation in angiosperms can be accompanied by changes in floral colour that may influence pollinator preference and reproductive isolation. This study investigates whether changes in floral colour can accompany polyploid and homoploid hybridization, important processes in angiosperm evolution.

METHODS

Spectral reflectance of corolla tissue was examined for 60 Nicotiana (Solanaceae) accessions (41 taxa) based on spectral shape (corresponding to pigmentation) as well as bee and hummingbird colour perception in order to assess patterns of floral colour evolution. Polyploid and homoploid hybrid spectra were compared with those of their progenitors to evaluate whether hybridization has resulted in floral colour shifts.

KEY RESULTS

Floral colour categories in Nicotiana seem to have arisen multiple times independently during the evolution of the genus. Most younger polyploids displayed an unexpected floral colour, considering those of their progenitors, in the colour perception of at least one pollinator type, whereas older polyploids tended to resemble one or both of their progenitors.

CONCLUSIONS

Floral colour evolution in Nicotiana is weakly constrained by phylogeny, and colour shifts do occur in association with both polyploid and homoploid hybrid divergence. Transgressive floral colour in N. tabacum has arisen by inheritance of anthocyanin pigmentation from its paternal progenitor while having a plastid phenotype like its maternal progenitor. Potentially, floral colour evolution has been driven by, or resulted in, pollinator shifts. However, those polyploids that are not sympatric (on a regional scale) with their progenitor lineages are typically not divergent in floral colour from them, perhaps because of a lack of competition for pollinators.

Do phosphoinositides regulate membrane water permeability of tobacco protoplasts by enhancing the aquaporin pathway?

MAIN CONCLUSION

Enhancing the membrane content of PtdInsP 2 , the already-recognized protein-regulating lipid, increased the osmotic water permeability of tobacco protoplasts, apparently by increasing the abundance of active aquaporins in their membranes. While phosphoinositides are implicated in cell volume changes and are known to regulate some ion channels, their modulation of aquaporins activity has not yet been reported for any organism. To examine this, we compared the osmotic water permeability (P f) of protoplasts isolated from tobacco (Nicotiana tabacum) cultured cells (NT1) with different (genetically lowered or elevated relative to controls) levels of inositol trisphosphate (InsP3) and phosphatidyl inositol [4,5] bisphosphate (PtdInsP2). To achieve this, the cells were transformed with, respectively, the human InsP3 5-phosphatase ('Ptase cells') or human phosphatidylinositol (4) phosphate 5-kinase ('PIPK cells'). The mean P f of the PIPK cells was several-fold higher relative to that of controls and Ptase cells. Three results favor aquaporins over the membrane matrix as underlying this excessive P f: (1) transient expression of the maize aquaporin ZmPIP2;4 in the PIPK cells increased P f by 12-30 μm s(-1), while in the controls only by 3-4 μm s(-1). (2) Cytosol acidification-known to inhibit aquaporins-lowered the P f in the PIPK cells down to control levels. (3) The transcript of at least one aquaporin was elevated in the PIPK cells. Together, the three results demonstrate the differences between the PIPK cells and their controls, and suggest a hitherto unobserved regulation of aquaporins by phosphoinositides, which could occur through direct interaction or indirect phosphoinositides-dependent cellular effects.

Persistence of subgenomes in paleopolyploid cotton after 60 my of evolution

The importance of whole-genome multiplication (WGM) in plant evolution has long been recognized. In flowering plants, WGM is both ubiquitous and in many lineages cyclical, each round followed by substantial gene loss (fractionation). This process may be biased with respect to duplicated chromosomes, often with overexpression of genes in less fractionated relative to more fractionated regions. This bias is hypothesized to arise through downregulation of gene expression through silencing of local transposable elements (TEs). We assess differences in gene expression between duplicated regions of the paleopolyploid cotton genome and demonstrate that the rate of fractionation is negatively correlated with gene expression. We examine recent hypotheses regarding the source of fractionation bias and show that TE-mediated, positional downregulation is absent in the modern cotton genome, seemingly excluding this phenomenon as the primary driver of biased gene loss. Nevertheless, the paleo subgenomes of diploid cotton are still distinguishable with respect to TE content, targeting of 24-nt-small interfering RNAs and GC content, despite approximately 60 My of evolution. We propose that repeat content per se and differential recombination rates may drive biased fractionation following WGM. These data highlight the likely importance of ancient genomic fractionation biases in shaping modern crop genomes.

Genomic repeat abundances contain phylogenetic signal

A large proportion of genomic information, particularly repetitive elements, is usually ignored when researchers are using next-generation sequencing. Here we demonstrate the usefulness of this repetitive fraction in phylogenetic analyses, utilizing comparative graph-based clustering of next-generation sequence reads, which results in abundance estimates of different classes of genomic repeats. Phylogenetic trees are then inferred based on the genome-wide abundance of different repeat types treated as continuously varying characters; such repeats are scattered across chromosomes and in angiosperms can constitute a majority of nuclear genomic DNA. In six diverse examples, five angiosperms and one insect, this method provides generally well-supported relationships at interspecific and intergeneric levels that agree with results from more standard phylogenetic analyses of commonly used markers. We propose that this methodology may prove especially useful in groups where there is little genetic differentiation in standard phylogenetic markers. At the same time as providing data for phylogenetic inference, this method additionally yields a wealth of data for comparative studies of genome evolution.

Differential genome evolution and speciation of Coix lacryma-jobi L. and Coix aquatica Roxb. hybrid guangxi revealed by repetitive sequence analysis and fine karyotyping

Abstract

Background

Coix, Sorghum and Zea are closely related plant genera in the subtribe Maydeae. Coix comprises 9–11 species with different ploidy levels (2n = 10, 20, 30, and 40). The exclusively cultivated C. lacryma-jobi L. (2n = 20) is widely used in East and Southeast Asia for food and medicinal applications. Three fertile cytotypes (2n = 10, 20, and 40) have been reported for C. aquatica Roxb. One sterile cytotype (2n = 30) closely related to C. aquatica has been recently found in Guangxi of China. This putative hybrid has been named C. aquatica HG (Hybrid Guangxi). The genome composition and the evolutionary history of C. lacryma-jobi and C. aquatica HG are largely unclear.

Results

About 76% of the genome of C. lacryma-jobi and 73% of the genome of C. aquatica HG are repetitive DNA sequences as shown by low coverage genome sequencing followed by similarity-based cluster analysis. In addition, long terminal repeat (LTR) retrotransposable elements are dominant repetitive sequences in these two genomes, and the proportions of many repetitive sequences in whole genome varied greatly between the two species, indicating evolutionary divergence of them. We also found that a novel 102 bp variant of centromeric satellite repeat CentX and two other satellites only appeared in C. aquatica HG. The results from FISH analysis with repeat probe cocktails and the data from chromosomes pairing in meiosis metaphase showed that C. lacryma-jobi is likely a diploidized paleotetraploid species and C. aquatica HG is possibly a recently formed hybrid. Furthermore, C. lacryma-jobi and C. aquatica HG shared more co-existing repeat families and higher sequence similarity with Sorghum than with Zea.

Conclusions

The composition and abundance of repetitive sequences are divergent between the genomes of C. lacryma-jobi and C. aquatica HG. The results from fine karyotyping analysis and chromosome pairing suggested diploidization of C. lacryma-jobi during evolution and C. aquatica HG is a recently formed hybrid. The genome-wide comparison of repetitive sequences indicated that the repeats in Coix were more similar to those in Sorghum than to those in Zea, which is consistent with the phylogenetic relationship reported by previous work.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-1025) contains supplementary material, which is available to authorized users.

The evolutionary fate of the horizontally transferred agrobacterial mikimopine synthase gene in the genera Nicotiana and Linaria

Few cases of spontaneously horizontally transferred bacterial genes into plant genomes have been described to date. The occurrence of horizontally transferred genes from the T-DNA of Agrobacterium rhizogenes into the plant genome has been reported in the genus Nicotiana and in the species Linaria vulgaris. Here we compare patterns of evolution in one of these genes (a gene encoding mikimopine synthase, mis) following three different events of horizontal gene transfer (HGT). As this gene plays an important role in Agrobacterium, and there are known cases showing that genes from pathogens can acquire plant protection function, we hypothesised that in at least some of the studied species we will find signs of selective pressures influencing mis sequence. The mikimopine synthase (mis) gene evolved in a different manner in the branch leading to Nicotiana tabacum and N. tomentosiformis, in the branch leading to N. glauca and in the genus Linaria. Our analyses of the genus Linaria suggest that the mis gene began to degenerate soon after the HGT. In contrast, in the case of N. glauca, the mis gene evolved under significant selective pressures. This suggests a possible role of mikimopine synthase in current N. glauca and its ancestor(s). In N. tabacum and N. tomentosiformis, the mis gene has a common frameshift mutation that disrupted its open reading frame. Interestingly, our results suggest that in spite of the frameshift, the mis gene could evolve under selective pressures. This sequence may still have some regulatory role at the RNA level as suggested by coverage of this sequence by small RNAs in N. tabacum.

Inactivation of two newly identified tobacco heavy metal ATPases leads to reduced Zn and Cd accumulation in shoots and reduced pollen germination

Cadmium (Cd) is a non-essential heavy metal, which is classified as a "known human carcinogen" by the International Agency for Research on Cancer (IARC). Understanding the mechanisms controlling Cd distribution in planta is essential to develop phytoremediation approaches as well as for food safety. Unlike most other plants, tobacco (Nicotiana tabacum) plants translocate most of the Cd taken up from the soil, out of the roots and into the shoots, leading to high Cd accumulation in tobacco shoots. Two orthologs of the Arabidopsis thaliana HMA2 and HMA4 Zn and Cd ATPases that are responsible for zinc (Zn) and Cd translocation from roots to shoots were identified in tobacco and sequenced. These genes, named NtHMAα and NtHMAβ, were more highly expressed in roots than in shoots. NtHMAα was expressed in the vascular tissues of both roots and leaves as well as in anthers. No visual difference was observed between wild-type plants and plants in which the NtHMAα and NtHMAβ genes were either mutated or silenced. These mutants showed reduced Zn and Cd accumulation in shoots as well as increased Cd tolerance. When both NtHMA genes were silenced, plant development was altered and pollen germination was severely impaired due to Zn deficiency. Interestingly, seeds from these lines also showed decreased Zn concentration but increased iron (Fe) concentration.

Species origin of genomic factors in Nicotiana nudicaulis Watson controlling hybrid lethality in interspecific hybrids between N. nudicaulis Watson and N. tabacum L

Hybrid lethality is expressed at 28°C in the cross Nicotiana nudicaulis × N. tabacum. The S subgenome of N. tabacum has been identified as controlling this hybrid lethality. To clarify the responsible genomic factor(s) of N. nudicaulis, we crossed N. trigonophylla (paternal progenitor of N. nudicaulis) with N. tabacum, because hybrids between N. sylvestris (maternal progenitor of N. nudicaulis) and N. tabacum are viable when grown in a greenhouse. In the cross N. trigonophylla×N. tabacum, approximately 50% of hybrids were vitrified, 20% were viable, and 20% were nonviable at 28°C. To reveal which subgenome of N. tabacum was responsible for these phenotypes, we crossed N. trigonophylla with two progenitors of N. tabacum, N. sylvestris (SS) and N. tomentosiformis (TT). In the cross N. sylvestris × N. trigonophylla, we confirmed that over half of hybrids of N. sylvestris × N. trigonophylla were vitrified, and none of the hybrids of N. trigonophylla × N. tomentosiformis were. The results imply that the S subgenome, encoding a gene or genes inducing hybrid lethality in the cross between N. nudicaulis and N. tabacum, has one or more genomic factors that induce vitrification. Furthermore, in vitrified hybrids of N. trigonophylla × N. tabacum and N. sylvestris × N. trigonophylla, we found that nuclear fragmentation, which progresses during expression of hybrid lethality, was accompanied by vitrification. This observation suggests that vitrification has a relationship to hybrid lethality. Based on these results, we speculate that when N. nudicaulis was formed approximately 5 million years ago, several causative genomic factors determining phenotypes of hybrid seedlings were inherited from N. trigonophylla. Subsequently, genome downsizing and various recombination-based processes took place. Some of the causative genomic factors were lost and some became genomic factor(s) controlling hybrid lethality in extant N. nudicaulis.