Comparative analysis of full-length transcriptomes based on hybrid population reveals regulatory mechanisms of anthocyanin biosynthesis in sweet potato (Ipomoea batatas (L.) Lam)

BACKGROUND

Sweet potato (Ipomoea batatas (L.) Lam.) is a highly heterozygous autohexaploid crop with high yield and high anthocyanin content. Purple sweet potato is the main source of anthocyanins, and the mechanism of anthocyanin biosynthesis in storage roots has not been fully revealed.

RESULTS

In order to reveal the mechanism of anthocyanin biosynthesis and identify new homologous genes involved in anthocyanin biosynthesis in the storage roots of sweet potato, we used Ningzishu 1 and Jizishu 2 as parents to construct a F₁ hybrid population. Seven anthocyanin-containing lines and three anthocyanin-free lines were selected for full-length and second-generation transcriptome analyses. A total of 598,375 circular consensus sequencing reads were identified from full-length transcriptome sequencing. After analysis and correction of second-generation transcriptome data, 41,356 transcripts and 18,176 unigenes were obtained. Through a comparative analysis between anthocyanin-containing and anthocyanin-free groups 2329 unigenes were found to be significantly differentially expressed, of which 1235 were significantly up-regulated and 1094 were significantly down-regulated. GO enrichment analysis showed that the differentially expressed unigenes were significantly enriched in molecular function and biological process. KEGG enrichment analysis showed that the up-regulated unigenes were significantly enriched in the flavonoid biosynthesis and phenylpropanoid biosynthesis pathways, and the down-regulated unigenes were significantly enriched in the plant hormone signal transduction pathway. Weighted gene co-expression network analysis of differentially expressed unigenes revealed that anthocyanin biosynthesis genes were co-expressed with transcription factors such as MYB, bHLH and WRKY at the transcription level.

CONCLUSIONS

Our study will shed light on the regulatory mechanism of anthocyanin biosynthesis in sweet potato storage roots at the transcriptome level.

Asymmetrical canina meiosis is accompanied by the expansion of a pericentromeric satellite in non-recombining univalent chromosomes in the genus Rosa

BACKGROUND AND AIMS

Despite their abundant odd-ploidy (2n = 5x = 35), dogroses (Rosa sect. Caninae) are capable of sexual reproduction due to their unique meiosis. During canina meiosis, two sets of chromosomes form bivalents and are transmitted by male and female gametes, whereas the remaining chromosomes form univalents and are exclusively transmitted by the egg cells. Thus, the evolution of chromosomes is expected to be driven by their behaviour during meiosis.

METHODS

To gain insight into differential chromosome evolution, fluorescence in situ hybridization was conducted for mitotic and meiotic chromosomes in four dogroses (two subsections) using satellite and ribosomal DNA probes. By exploiting high-throughput sequencing data, we determined the abundance and diversity of the satellite repeats in the genus Rosa by analysing 20 pentaploid, tetraploid and diploid species in total.

KEY RESULTS

A pericentromeric satellite repeat, CANR4, was found in all members of the genus Rosa, including the basal subgenera Hulthemia and Hesperhodos. The satellite was distributed across multiple chromosomes (5-20 sites per mitotic cell), and its genomic abundance was higher in pentaploid dogroses (2.3 %) than in non-dogrose species (1.3 %). In dogrose meiosis, univalent chromosomes were markedly enriched in CANR4 repeats based on both the number and the intensity of the signals compared to bivalent-forming chromosomes. Single-nucleotide polymorphisms and cluster analysis revealed high intragenomic homogeneity of the satellite in dogrose genomes.

CONCLUSIONS

The CANR4 satellite arose early in the evolution of the genus Rosa. Its high content and extraordinary homogeneity in dogrose genomes is explained by its recent amplification in non-recombining chromosomes. We hypothesize that satellite DNA expansion may contribute to the divergence of univalent chromosomes in Rosa species with non-symmetrical meiosis.

SUPPRESSOR OF GAMMA RESPONSE 1 acts as a regulator coordinating crosstalk between DNA damage response and immune response in Arabidopsis thaliana

Plants live in constantly changing and often unfavorable or stressful environments. Environmental changes induce biotic and abiotic stress, which, in turn, may cause genomic DNA damage. Hence, plants simultaneously suffer abiotic/biotic stress and DNA damage. However, little information is available on the signaling crosstalk that occurs between DNA damage and abiotic/biotic stresses. Arabidopsis thaliana SUPPRESSOR OF GAMMA RESPONSE1 (SOG1) is a pivotal transcription factor that regulates thousands of genes in response to DNA double-strand break (DSB), and we recently reported that SOG1 has a role in immune responses. In the present study, the effects of SOG1 overexpression on the DNA damage and immune responses were examined. Results found that SOG1 overexpression enhances the regulation of numerous downstream genes. Relative to the wild type plants, then, DNA damage responses were observed to be strongly induced. SOG1 overexpression also upregulates chitin (a major components of fungal cell walls) responsive genes in the presence of DSBs, implying that pathogen defense response is activated by DNA damage via SOG1. Further, SOG1 overexpression enhances fungal resistance. These results suggest that SOG1 regulates crosstalk between DNA damage response and the immune response and that plants have evolved a sophisticated defense network to contend with environmental stress.

Complete chloroplast genome sequence determination of Rheum species and comparative chloroplast genomics for the members of Rumiceae

Complete plastomes of Rheum species facilitated to clarify the phylogeny of Polygonaceae, and comparative chloroplast genomics contributed to develop genetic markers for the authentication of Rheum species. Rheum (Polygonaceae) is widely distributed throughout the temperate and subtropical areas of Asian interior. Rheum species are usually perennial herbs, and half of them are endemic to China with important medicinal properties. On account of similar morphological characteristics, species delimitation of Rheum still remains unclear. Chloroplast genomes of eight Rheum species, Rumex crispus and Oxyria digyna were characterized. Based on the comparison of genome structure of these species and the two published Rheum species, it is shown that plastome sequences of these species are relatively conserved with the same gene order, and three Sect. Palmata species remarkably showed high sequence similarities. Some hotspots could be used to discriminate the Rheum species, and 17 plastid genes were subject to positive selection. The phylogenetic analyses indicated that all the Polygonaceae species were clustered in the same group and showed that Rheum species, except for Rheum wittrockii, formed a monophyletic group with high maximum parsimony/maximum likelihood bootstrap support values and Bayesian posterior probabilities. The molecular dating based on plastomes indicated that the divergences within Polygonaceae species were dated to the Upper Cretaceous period [73.86-77.99 million years ago (Ma)]. The divergence of Sect. Palmata species was estimated to have occurred around 1.60 Ma, indicating that its diversification was affected by the repeated climatic fluctuation in the Quaternary.

Resequencing 93 accessions of coffee unveils independent and parallel selection during Coffea species divergence

Coffea arabica, C. canephora and C. excelsa, with differentiated morphological traits and distinct agro-climatic conditions, compose the majority of the global coffee plantation. To comprehensively understand their genetic diversity and divergence for future genetic improvement requires high-density markers. Here, we sequenced 93 accessions encompassing these three Coffea species, uncovering 15,367,960 single-nucleotide polymorphisms (SNPs). These SNPs are unequally distributed across different genomic regions and gene families, with two disease-resistant gene families showing the highest SNP density, suggesting strong balancing selection. Meanwhile, the allotetraploid C. arabica exhibits greater nucleotide diversity, followed by C. canephora and C. excelsa. Population divergence (F_ST), population stratification and phylogeny all support strong divergence among species, with C. arabica and its parental species C. canephora being closer genetically. Scanning of genomic islands with elevated F_ST and structure-disruptive SNPs contributing to species divergence revealed that most of the selected genes in each lineage are independent, with a few being selected in parallel for two or three species, such as genes in root hair cell development, flavonols accumulation and disease-resistant genes. Moreover, some of the SNPs associated with coffee lipids exhibit significantly biased allele frequency among species, being valuable for interspecific breeding. Overall, our study not only uncovers the key population genomic patterns among species but also contributes a substantial genomic resource for coffee breeding.

Genome size variation at constant chromosome number is not correlated with repetitive DNA dynamism in Anacyclus (Asteraceae)

BACKGROUND AND AIMS

Changes in the amount of repetitive DNA (dispersed and tandem repeats) are considered the main contributors to genome size variation across plant species in the absence of polyploidy. However, the study of repeatome dynamism in groups showing contrasting genomic features and complex evolutionary histories is needed to determine whether other processes underlying genome size variation may have been overlooked. The main aim here was to elucidate which mechanism best explains genome size evolution in Anacyclus (Asteraceae).

METHODS

Using data from Illumina sequencing, we analysed the repetitive DNA in all species of Anacyclus, a genus with a reticulate evolutionary history, which displays significant genome size and karyotype diversity albeit presenting a stable chromosome number.

KEY RESULTS

By reconstructing ancestral genome size values, we inferred independent episodes of genome size expansions and contractions during the evolution of the genus. However, analysis of the repeatome revealed a similar DNA repeat composition across species, both qualitative and quantitative. Using comparative methods to study repeatome dynamics in the genus, we found no evidence for repeat activity causing genome size variation among species.

CONCLUSIONS

Our results, combined with previous cytogenetic data, suggest that genome size differences in Anacyclus are probably related to chromosome rearrangements involving losses or gains of chromosome fragments, possibly associated with homoploid hybridization. These could represent balanced rearrangements that do not disrupt gene dosage in merged genomes, for example via chromosome segment exchanges.

Molecular evolution of chloroplast genomes in Monsteroideae (Araceae)

This study provides broad insight into the chloroplast genomes of the subfamily Monsteroideae. The identified polymorphic regions may be suitable for designing unique and robust molecular markers for phylogenetic inference. Monsteroideae is the third largest subfamily (comprises 369 species) and one of the early diverging lineages of the monocot plant family Araceae. The phylogeny of this important subfamily is not well resolved at the species level due to scarcity of genomic resources and suitable molecular markers. Here, we report annotated chloroplast genome sequences of four Monsteroideae species: Spathiphyllum patulinervum, Stenospermation multiovulatum, Monstera adansonii, and Rhaphidophora amplissima. The quadripartite chloroplast genomes (size range 163,335-164,751 bp) consist of a pair of inverted repeats (25,270-25,931 bp), separating a small single copy region (21,448-22,346 bp) from a large single copy region (89,714-91,841 bp). The genomes contain 114 unique genes, including four rRNA genes, 80 protein-coding genes, and 30 tRNA genes. Gene features, amino acid frequencies, codon usage, GC contents, oligonucleotide repeats, and inverted repeats dynamics exhibit similarities among the four genomes. Higher rate of synonymous substitutions was observed as compared to non-synonymous substitutions in 76 protein-coding genes. Positive selection was observed in seven protein-coding genes, including psbK, ndhK, ndhD, rbcL, accD, rps8, and ycf2. Our included species of Araceae showed the monophyly in Monsteroideae and other subfamilies. We report 30 suitable polymorphic regions. The polymorphic regions identified here might be suitable for designing unique and robust markers for inferring the phylogeny and phylogeography among closely related species within the genus Spathiphyllum and among distantly related species within the subfamily Monsteroideae. The chloroplast genomes presented here are a valuable contribution towards understanding the molecular evolutionary dynamics in the family Araceae.

Non-destructive DNA extraction from herbarium specimens: a method particularly suitable for plants with small and fragile leaves

Protocols for DNA extraction from plants generally involve physical and chemical destruction of tissues. Use of these conventional methods precludes preservation of morphological information from herbarium specimens, especially for small plants with few leaves, and reduces the voucher value of specimens. Here, we developed a new, non-destructive DNA extraction protocol (Protocol 1) that only needs a small piece of leaf (< 25 mm²) to obtain DNA suitable for DNA sequencing from fragile herbarium specimens. The protocol was very simple and rapid; an extraction buffer was placed on the leaf surface of an intact specimen for 30 min at room temperature (20 °C). The quality of extracted DNA was checked by PCR amplification of two standard plant DNA barcode regions, the maturase K gene (matK, ca. 850 bp) and the ribulose-1,5-bisphosphatecarboxylase/oxygenase gene (rbcL, ca. 550 bp), for 14 vascular plant species encompassing various taxonomic groups. The protocol retrieved sequences from 80.0% of specimens for matK and 46.2% of specimens for rbcL. Placing of the extraction buffer onto specimens did not cause any tears or deformation, but caused discoloration in some plants. To improve DNA yield for specimens incompatible with Protocol 1, we developed an alternative protocol for DNA extraction with minimally invasive destruction of specimens (Protocol 2). In this protocol, a cut leaf was immersed in the extraction buffer for 30 min and stored subsequently in a fragment pocket on the specimen sheet. This alternative method retrieved matK sequences from 80.0% of specimens and rbcL sequences from 92.8% of specimens. The combination of Protocols 1 and 2 enabled us to obtain matK sequences from 90.0% of specimens and rbcL sequences form 92.8% of specimens. The new protocols facilitate the use of museum specimens for use of DNA of museum specimens while still preserving morphological information.

Evolution of small RNA expression following hybridization and allopolyploidization: insights from Spartina species (Poaceae, Chloridoideae)

Differential expression of mi-RNAs targeting developmental processes and progressive downregulation of repeat-associated siRNAs following genome merger and genome duplication in the context of allopolyploid speciation in Spartina. The role of small RNAs on gene expression regulation and genome stability is arousing increased interest and is being explored in various plant systems. In spite of prominence of reticulate evolution and polyploidy that affects the evolutionary history of all plant lineages, very few studies analysed RNAi mechanisms with this respect. Here, we explored small RNAs diversity and expression in the context of recent allopolyploid speciation, using the Spartina system, which offers a unique opportunity to explore the immediate changes following hybridization and genome duplication. Small RNA-Seq analyses were conducted on hexaploid parental species (S. alterniflora and S. maritima), their F1 hybrid S. x townsendii, and the neoallododecaploid S. anglica. We identified 594 miRNAs, 2197 miRNA-target genes, and 3730 repeat-associated siRNAs (mostly targeting Class I/Copia-Ivana- Copia-SIRE and LINEs elements). For both mi- and ra-siRNAs, we detected differential expression patterns following genome merger and genome duplication. These misregulations include non-additive expression of miRNAs in the F1 hybrid and additional changes in the allopolyploid targeting developmental processes. Expression of repeat-associated siRNAs indicates a strengthen of transposable element repression during the allopolyploidization process. Altogether, these results confirm the central role small RNAs play in shaping regulatory changes in naturally formed recent allopolyploids.

Repetitive DNA landscape in essential A and supernumerary B chromosomes of Festuca pratensis Huds

Here, we characterized the basic properties of repetitive sequences in essential A and supernumerary B chromosomes of Festuca pratensis Huds. This was performed by comparative analysis of low-pass Illumina sequence reads of B chromosome lacking (-B) and B chromosome containing (+B) individuals of F. pratensis. 61% of the nuclear genome is composed of repetitive sequences. 43.1% of the genome are transposons of which DNA transposons and retrotransposons made up 2.3% and 40.8%, respectively. LTR retrotransposons are the most abundant mobile elements and contribute to 40.7% of the genome and divided into Ty3-gypsy and Ty1-copia super families with 32.97% and 7.78% of the genome, respectively. Eighteen different satellite repeats were identified making up 3.9% of the genome. Five satellite repeats were used as cytological markers for chromosome identification and genome analysis in the genus Festuca. Four satellite repeats were identified on B chromosomes among which Fp-Sat48 and Fp-Sat253 were specific to the B chromosome of F. pratensis.

DNA content equivalence in haploid and diploid maize leaves

The qPCR assay developed to differentiate haploid and diploid maize leaf samples was unsuccessful due to DNA content difference. Haploid cells are packed more closely together with less cellular expansion. Increased ploidy content (> 2 N) directly correlates with increased cell size in plants, but few studies have examined cell morphology in plants with reduced ploidy (i.e., haploids). To pioneer a scalable new ploidy test, we compared DNA content and cellular morphology of haploid and diploid maize leaves. The amount of genomic DNA recovered from standardized leaf-punch samples was equivalent between these two ploidy types, while both epidermal and mesophyll cell types were smaller in haploid plants. Pavement cells had a substantially smaller size than mesophyll cells, and this effect was more pronounced in the abaxial epidermis. Interveinal distance and guard cell size were significantly reduced in haploids, but the cell percentage comprising stomata did not change. These results confirm the direct correlation between ploidy content and cell size in plants, and suggest that reduced cell expansion predominantly explains DNA content equivalence between haploid and diploid samples, confounding efforts to develop a haploid detection method using DNA content.

Identification and Monitoring of Amomi Fructus and its Adulterants Based on DNA Barcoding Analysis and Designed DNA Markers

Amomi Fructus is one of the traditional medicines derived from the ripe fruits of the Zingiberaceae family of plants, which include Amomum villosum, A. villosum var. xanthioides, and A. longiligulare. Owing to their highly similar morphological traits, several kinds of adulterants of Amomi Fructus have been reported. Therefore, accurate and reliable methods of identification are necessary in order to ensure drug safety and quality. We performed DNA barcoding using five regions (ITS, matK, rbcL, rpoB, and trnL-F intergenic spacer) of 23 Amomi Fructus samples and 22 adulterants. We designed specific DNA markers for Amomi Fructus based on the single nucleotide polymorphisms (SNPs) in the ITS. Amomi Fructus was well separated from the adulterants and was classified with the species of origin based on the detected SNPs from the DNA barcoding results. The AVF1/ISR DNA marker for A. villosum produced a 270 bases amplified product, while the ALF1/ISF DNA marker produced a 350 bases product specific for A. longiligulare. Using these DNA markers, the monitoring of commercially distributed Amomi Fructus was performed, and the monitoring results were confirmed by ITS analysis. This method identified samples that were from incorrect origins, and a new species of adulterant was also identified. These results confirmed the accuracy and efficiency of the designed DNA markers; this method may be used as an efficient tool for the identification and verification of Amomi Fructus.

Diversity of repetitive sequences within compact genomes of Phaseolus L. beans and allied genera Cajanus L. and Vigna Savi

Repetitive sequences are ubiquitous and fast-evolving elements responsible for size variation and large-scale organization of plant genomes. Within tribe Phaseoleae (Fabaceae), some genera, such as Phaseolus, Vigna, and Cajanus, show small genome and mostly stable chromosome number. Here, we applied a combined computational and cytological approach to study the organization and diversification of repetitive elements in some species of these genera. Sequences were classified in terms of type and repetitiveness and the most abundant were mapped to chromosomes. We identified long terminal repeat (LTR) retrotransposons, especially Ogre and Chromovirus elements, making up most of genomes, other than P. acutifolius and Vigna species. Satellite DNAs (SatDNAs) were less representative, but highly diverse among species, showing a clear phylogenetic relationship. In situ localization revealed preferential location at pericentromeres and centromeres for both types of sequences, suggesting a heterogeneous composition, especially for centromeres. Few elements showed subterminal accumulation. Copy number variation among chromosomes within and among species was observed for all nine identified SatDNAs. Altogether, our data pointed two main elements (Ty3/Gypsy retrotransponsons and SatDNAs) to the diversification on the repetitive landscape in Phaseoleae, with a typical set of repeats in each species. The high turnover of these sequences, however, did not affect total genome size.

Genome-wide identification and expression analysis of SWEET gene family in Litchi chinensis reveal the involvement of LcSWEET2a/3b in early seed development

BACKGROUND

SWEETs (Sugar Will Eventually be Exported transporters) function as sugar efflux transporters that perform diverse physiological functions, including phloem loading, nectar secretion, seed filling, and pathogen nutrition. The SWEET gene family has been identified and characterized in a number of plant species, but little is known about in Litchi chinensis, which is an important evergreen fruit crop.

RESULTS

In this study, 16 LcSWEET genes were identified and nominated according to its homologous genes in Arabidopsis and grapevine. Multiple sequence alignment showed that the 7 alpha-helical transmembrane domains (7-TMs) were basically conserved in LcSWEETs. The LcSWEETs were divided into four clades (Clade I to Clade IV) by phylogenetic tree analysis. A total of 8 predicted motifs were detected in the litchi LcSWEET genes. The 16 LcSWEET genes were unevenly distributed in 9 chromosomes and there was one pairs of segmental duplicated events by synteny analysis. The expression patterns of the 16 LcSWEET genes showed higher expression levels in reproductive organs. The temporal and spatial expression patterns of LcSWEET2a and LcSWEET3b indicated they play central roles during early seed development.

CONCLUSIONS

The litchi genome contained 16 SWEET genes, and most of the genes were expressed in different tissues. Gene expression suggested that LcSWEETs played important roles in the growth and development of litchi fruits. Genes that regulate early seed development were preliminarily identified. This work provides a comprehensive understanding of the SWEET gene family in litchi, laying a strong foundation for further functional studies of LcSWEET genes and improvement of litchi fruits.

Intergenomic gene transfer in diploid and allopolyploid Gossypium

BACKGROUND

Intergenomic gene transfer (IGT) between nuclear and organellar genomes is a common phenomenon during plant evolution. Gossypium is a useful model to evaluate the genomic consequences of IGT for both diploid and polyploid species. Here, we explore IGT among nuclear, mitochondrial, and plastid genomes of four cotton species, including two allopolyploids and their model diploid progenitors (genome donors, G. arboreum: A₂ and G. raimondii: D₅).

RESULTS

Extensive IGT events exist for both diploid and allotetraploid cotton (Gossypium) species, with the nuclear genome being the predominant recipient of transferred DNA followed by the mitochondrial genome. The nuclear genome has integrated 100 times more foreign sequences than the mitochondrial genome has in total length. In the nucleus, the integrated length of chloroplast DNA (cpDNA) was between 1.87 times (in diploids) to nearly four times (in allopolyploids) greater than that of mitochondrial DNA (mtDNA). In the mitochondrion, the length of nuclear DNA (nuDNA) was typically three times than that of cpDNA. Gossypium mitochondrial genomes integrated three nuclear retrotransposons and eight chloroplast tRNA genes, and incorporated chloroplast DNA prior to divergence between the diploids and allopolyploid formation. For mitochondrial chloroplast-tRNA genes, there were 2-6 bp conserved microhomologies flanking their insertion sites across distantly related genera, which increased to 10 bp microhomologies for the four cotton species studied. For organellar DNA sequences, there are source hotspots, e.g., the atp6-trnW intergenic region in the mitochondrion and the inverted repeat region in the chloroplast. Organellar DNAs in the nucleus were rarely expressed, and at low levels. Surprisingly, there was asymmetry in the survivorship of ancestral insertions following allopolyploidy, with most numts (nuclear mitochondrial insertions) decaying or being lost whereas most nupts (nuclear plastidial insertions) were retained.

CONCLUSIONS

This study characterized and compared intracellular transfer among nuclear and organellar genomes within two cultivated allopolyploids and their ancestral diploid cotton species. A striking asymmetry in the fate of IGTs in allopolyploid cotton was discovered, with numts being preferentially lost relative to nupts. Our results connect intergenomic gene transfer with allotetraploidy and provide new insight into intracellular genome evolution.

A novel miniature transposon-like element discovered in the coding sequence of a gene that encodes for 5-formyltetrahydrofolate in wheat

BACKGROUND

Transposable elements (TEs) comprise over 80% of the wheat genome and usually possess unique features for specific super-families and families. However, the role of TEs in wheat evolution and reshaping the wheat genome remains largely unclear.

RESULTS

In this study, we discovered a miniature (307 bp in length) TE-like sequence in exon 6 of a gene that encodes for 5-formyltetrahydrofolate, in two accessions of wild emmer wheat (T. turgidum ssp. dicoccoides) and has interfered with the gene translation by creating a shorter reading frame as a result of a stop codon. The sequence that was termed Mariam, does not show any structural similarity to known TEs. It does not possess terminal inverted repeats (TIRs) that would allow us to assign this element to one of the TIR DNA super-families, and it does not possess characteristic features of SINE, such as a Pol-III promotor or a poly-A tail. In-silico analysis of five publicly available genome drafts of Triticum and Aegilops species revealed that Mariam element appears in a very low copy number (1-3 insertions) in diploid wheat species and ~ 12 insertions in tetraploid and hexaploidy wheat species. In addition, Mariam element was found to be unique to wheat, as it was not found in other plant genomes. The dynamic nature of Mariam in the wheat genome was assessed by site-specific PCR analysis and revealed that it retained activity in wild emmer populations in a population-specific manner.

CONCLUSIONS

This study provides additional insight into the evolutionary impact of TEs in wheat.

A taxonomic monograph of Ipomoea integrated across phylogenetic scales

Taxonomic monographs have the potential to make a unique contribution to the understanding of global biodiversity. However, such studies, now rare, are often considered too daunting to undertake within a realistic time frame, especially as the world's collections have doubled in size in recent times. Here, we report a global-scale monographic study of morning glories (Ipomoea) that integrated DNA barcodes and high-throughput sequencing with the morphological study of herbarium specimens. Our approach overhauled the taxonomy of this megadiverse group, described 63 new species and uncovered significant increases in net diversification rates comparable to the most iconic evolutionary radiations in the plant kingdom. Finally, we show that more than 60 species of Ipomoea, including sweet potato, independently evolved storage roots in pre-human times, indicating that the storage root is not solely a product of human domestication but a trait that predisposed the species for cultivation. This study demonstrates how the world's natural history collections can contribute to global challenges in the Anthropocene.

A 3,000-year-old Egyptian emmer wheat genome reveals dispersal and domestication history

Tetraploid emmer wheat (Triticum turgidum ssp. dicoccon) is a progenitor of the world's most widely grown crop, hexaploid bread wheat (Triticum aestivum), as well as the direct ancestor of tetraploid durum wheat (T. turgidum subsp. turgidum). Emmer was one of the first cereals to be domesticated in the old world; it was cultivated from around 9700 BC in the Levant^1,2 and subsequently in south-western Asia, northern Africa and Europe with the spread of Neolithic agriculture^3,4. Here, we report a whole-genome sequence from a museum specimen of Egyptian emmer wheat chaff, ¹⁴C dated to the New Kingdom, 1130-1000 BC. Its genome shares haplotypes with modern domesticated emmer at loci that are associated with shattering, seed size and germination, as well as within other putative domestication loci, suggesting that these traits share a common origin before the introduction of emmer to Egypt. Its genome is otherwise unusual, carrying haplotypes that are absent from modern emmer. Genetic similarity with modern Arabian and Indian emmer landraces connects ancient Egyptian emmer with early south-eastern dispersals, whereas inferred gene flow with wild emmer from the Southern Levant signals a later connection. Our results show the importance of museum collections as sources of genetic data to uncover the history and diversity of ancient cereals.

An update and reassessment of fern and lycophyte diversity data in the Japanese Archipelago

The fern and lycophyte flora of Japan comprising 721 native taxa (including subspecies and varieties) plus 371 interspecific hybrids was reassessed using a nearly comprehensively sampled distribution map at 10 km resolution vouchered by 216,687 specimens, up-to-date cytotaxonomic information covering 74% of the taxa, and an rbcL sequence dataset covering 97.9% of the taxa. Spatial distribution of species richness and phylogenetic diversity was visualized. Apomixis was observed in 11.0% of the native taxa whose reproductive modes are known. The number of sexually reproducing polyploid taxa (n = 199) is less than sexual diploids (n = 241), and 30 of them are evidently allopolyploid, in contrast with the low number of possible autopolyploids (n = 4). Apomictic taxa were found to have smaller latitudinal ranges than sexual taxa or taxa with multiple reproductive modes. A morphological character dataset in Lucid format is provided for taxonomic identification of the native taxa.

An update and reassessment of fern and lycophyte diversity data in the Japanese Archipelago

The fern and lycophyte flora of Japan comprising 721 native taxa (including subspecies and varieties) plus 371 interspecific hybrids was reassessed using a nearly comprehensively sampled distribution map at 10 km resolution vouchered by 216,687 specimens, up-to-date cytotaxonomic information covering 74% of the taxa, and an rbcL sequence dataset covering 97.9% of the taxa. Spatial distribution of species richness and phylogenetic diversity was visualized. Apomixis was observed in 11.0% of the native taxa whose reproductive modes are known. The number of sexually reproducing polyploid taxa (n = 199) is less than sexual diploids (n = 241), and 30 of them are evidently allopolyploid, in contrast with the low number of possible autopolyploids (n = 4). Apomictic taxa were found to have smaller latitudinal ranges than sexual taxa or taxa with multiple reproductive modes. A morphological character dataset in Lucid format is provided for taxonomic identification of the native taxa.

UGT709G1: a novel uridine diphosphate glycosyltransferase involved in the biosynthesis of picrocrocin, the precursor of safranal in saffron (Crocus sativus)

Saffron, a spice derived from the dried red stigmas of Crocus sativus, is one of the oldest natural food additives. The flowers have long red stigmas, which store significant quantities of the glycosylated apocarotenoids crocins and picrocrocin. The apocarotenoid biosynthetic pathway in saffron starts with the oxidative cleavage of zeaxanthin, from which crocins and picrocrocin are derived. In the processed stigmas, picrocrocin is converted to safranal, giving saffron its typical aroma. By a targeted search for differentially expressed uridine diphosphate glycosyltransferases (UGTs) in Crocus transcriptomes, a novel apocarotenoid glucosyltransferase (UGT709G1) from saffron was identified. Biochemical analyses revealed that UGT709G1 showed a high catalytic efficiency toward 2,6,6-trimethyl-4-hydroxy-1-carboxaldehyde-1-cyclohexene (HTCC), making it suited for the biosynthesis of picrocrocin, the precursor of safranal. The role of UGT709G1 in picrocrocin/safranal biosynthesis was supported by the absence or presence of gene expression in a screening for HTCC and picrocrocin production in different Crocus species and by a combined transient expression assay with CsCCD2L in Nicotiana benthamiana leaves. The identification of UGT709G1 completes one of the most highly valued specialized metabolic biosynthetic pathways in plants and provides novel perspectives on the industrial production of picrocrocin to be used as a flavor additive or as a pharmacological constituent.

Genetic mapping and development of molecular markers for a candidate gene locus controlling rind color in watermelon

ClCG08G017810 (ClCGMenG) encoding a 2-phytyl-1,4-beta-naphthoquinone methyltransferase protein is associated with formation of dark green versus light green rind color in watermelon. Rind color is an important agronomic trait in watermelon [Citrullus lanatus (Thunb.) Matsum. and Nakai], but the underlying molecular mechanism for this trait is not fully known. In the present study, we identified a single locus on chromosome 8 accounting for watermelon rind color (dark green vs. light green). Genetic analysis of F₁, F₂, and BC₁ populations derived from two parental lines (9904 with dark green rind and Handel with light green rind) revealed that the watermelon rind color (dark green vs. light green) is controlled by a single locus, and dark green is dominant to light green rind. Initial mapping revealed a region of interest spanning 2.07 Mb on chromosome 8. Genetic mapping with CAPS and SNP markers narrowed down the candidate region to 31.4 kb. Gene annotation of the corresponding region in the reference genome revealed the ClCG08G017810 gene sequence encoding the 2-phytyl-1,4-beta-naphthoquinone methyltransferase protein. The sequence alignment of the candidate gene with the two parental lines suggested a nonsynonymous SNP mutation in the coding region of ClCG08G017810, converting an arginine (R) to glycine (G). The SNP might be associated with rind color of 103 watermelon germplasm lines investigated in this study. The qRT-PCR analysis revealed higher expression of ClCG08G017810 in dark green rind than in light green rind. Therefore, ClCG08G017810 is a candidate gene associated with watermelon rind color. The present study facilitates marker-assisted selection useful for the development of cultivars with desirable rind color.

Increase in lysophosphatidate acyltransferase activity in oilseed rape (Brassica napus) increases seed triacylglycerol content despite its low intrinsic flux control coefficient

Lysophosphatidate acyltransferase (LPAAT) catalyses the second step of the Kennedy pathway for triacylglycerol (TAG) synthesis. In this study we expressed Trapaeolum majus LPAAT in Brassica napus (B. napus) cv 12075 to evaluate the effects on lipid synthesis and estimate the flux control coefficient for LPAAT. We estimated the flux control coefficient of LPAAT in a whole plant context by deriving a relationship between it and overall lipid accumulation, given that this process is a exponential. Increasing LPAAT activity resulted in greater TAG accumulation in seeds of between 25% and 29%; altered fatty acid distributions in seed lipids (particularly those of the Kennedy pathway); and a redistribution of label from 14 C-glycerol between phosphoglycerides. Greater LPAAT activity in seeds led to an increase in TAG content despite its low intrinsic flux control coefficient on account of the exponential nature of lipid accumulation that amplifies the effect of the small flux increment achieved by increasing its activity. We have also developed a novel application of metabolic control analysis likely to have broad application as it determines the in planta flux control that a single component has upon accumulation of storage products.

Meiotic DNA Repair in the Nucleolus Employs a Nonhomologous End-Joining Mechanism

Ribosomal RNA genes are arranged in large arrays with hundreds of rDNA units in tandem. These highly repetitive DNA elements pose a risk to genome stability since they can undergo nonallelic exchanges. During meiosis, DNA double-strand breaks (DSBs) are induced as part of the regular program to generate gametes. Meiotic DSBs initiate homologous recombination (HR), which subsequently ensures genetic exchange and chromosome disjunction. In Arabidopsis (Arabidopsis thaliana), we demonstrate that all 45S rDNA arrays become transcriptionally active and are recruited into the nucleolus early in meiosis. This shields the rDNA from acquiring canonical meiotic chromatin modifications and meiotic cohesin and allows only very limited meiosis-specific DSB formation. DNA lesions within the rDNA arrays are repaired in an RAD51-independent but LIG4-dependent manner, establishing that nonhomologous end-joining maintains rDNA integrity during meiosis. Utilizing ectopically integrated rDNA repeats, we validate our findings and demonstrate that the rDNA constitutes an HR-refractory genome environment.