Repeat turnover meets stable chromosomes: repetitive DNA sequences mark speciation and gene pool boundaries in sugar beet and wild beets.

Sugar beet and its wild relatives share a base chromosome number of nine and similar chromosome morphologies. Yet, interspecific breeding is impeded by chromosome and sequence divergence that is still not fully understood. Since repetitive DNAs are among the fastest evolving parts of the genome, we investigated, if repeatome innovations and losses are linked to chromosomal differentiation and speciation. We traced genome and chromosome-wide evolution across 13 beet species comprising all sections of the genera Beta and Patellifolia. For this, we combined short and long read sequencing, flow cytometry, and cytogenetics to build a comprehensive framework that spans the complete scale from DNA to chromosome to genome. Genome sizes and repeat profiles reflect the separation into three gene pools with contrasting evolutionary patterns. Among all repeats, satellite DNAs harbor most genomic variability, leading to fundamentally different centromere architectures, ranging from chromosomal uniformity in Beta and Patellifolia to the formation of patchwork chromosomes in Corollinae/Nanae. We show that repetitive DNAs are causal for the genome expansions and contractions across the beet genera, providing insights into the genomic underpinnings of beet speciation. Satellite DNAs in particular vary considerably between beet genomes, leading to the evolution of distinct chromosomal setups in the three gene pools, likely contributing to the barriers in beet breeding. Thus, with their isokaryotypic chromosome sets, beet genomes present an ideal system for studying the link between repeats, genomic variability, and chromosomal differentiation and provide a theoretical fundament for understanding barriers in any crop breeding effort.

Generation and characterisation of an Arabidopsis thaliana f3h/fls1/ans triple mutant that accumulates eriodictyol derivatives.

BACKGROUND: Flavonoids are plant specialised metabolites, which derive from phenylalanine and acetate metabolism. They possess a variety of beneficial characteristics for plants and humans. Several modification steps in the synthesis of tricyclic flavonoids cause for the amazing diversity of flavonoids in plants. The 2-oxoglutarate-dependent dioxygenases (2-ODDs) flavanone 3-hydroxylase (F3H, synonym FHT), flavonol synthase (FLS) and anthocyanidin synthase (ANS, synonym leucoanthocyanidin dioxygenase (LDOX)), catalyse oxidative modifications to the central C ring. They are highly similar and have been shown to catalyse, at least in part, each other's reactions. FLS and ANS have been identified as bifunctional enzymes in many species, including Arabidopsis thaliana, stressing the capability of plants to bypass missing or mutated reaction steps on the way to flavonoid production. However, little is known about such bypass reactions and the flavonoid composition of plants lacking all three central flavonoid 2-ODDs. RESULTS: To address this issue, we generated a f3h/fls1/ans mutant, as well as the corresponding double mutants and investigated the flavonoid composition of this mutant collection. The f3h/fls1/ans mutant was further characterised at the genomic level by analysis of a nanopore DNA sequencing generated genome sequence assembly and at the transcriptomic level by RNA-Seq analysis. The mutant collection established, including the novel double mutants f3h/fls1 and f3h/ans, was used to validate and analyse the multifunctionalities of F3H, FLS1, and ANS in planta. Metabolite analyses revealed the accumulation of eriodictyol and additional glycosylated derivatives in mutants carrying the f3h mutant allele, resulting from the conversion of naringenin to eriodictyol by flavonoid 3'-hydroxylase (F3'H) activity. CONCLUSIONS: We describe the in planta multifunctionality of the three central flavonoid 2-ODDs from A. thaliana and identify a bypass in the f3h/fls1/ans triple mutant that leads to the formation of eriodictyol derivatives. As (homo-)eriodictyols are known as bitter taste maskers, the annotated eriodictyol (derivatives) and in particular the observations made on their in planta production, could provide valuable insights for the creation of novel food supplements.

Genomic characterization of a nematode tolerance locus in sugar beet.

BACKGROUND: Infection by beet cyst nematodes (BCN, Heterodera schachtii) causes a serious disease of sugar beet, and climatic change is expected to improve the conditions for BCN infection. Yield and yield stability under adverse conditions are among the main breeding objectives. Breeding of BCN tolerant sugar beet cultivars offering high yield in the presence of the pathogen is therefore of high relevance. RESULTS: To identify causal genes providing tolerance against BCN infection, we combined several experimental and bioinformatic approaches. Relevant genomic regions were detected through mapping-by-sequencing using a segregating F2 population. DNA sequencing of contrasting F2 pools and analyses of allele frequencies for variant positions identified a single genomic region which confers nematode tolerance. The genomic interval was confirmed and narrowed down by genotyping with newly developed molecular markers. To pinpoint the causal genes within the potential nematode tolerance locus, we generated long read-based genome sequence assemblies of the tolerant parental breeding line Strube U2Bv and the susceptible reference line 2320Bv. We analyzed continuous sequences of the potential locus with regard to functional gene annotation and differential gene expression upon BCN infection. A cluster of genes with similarity to the Arabidopsis thaliana gene encoding nodule inception protein-like protein 7 (NLP7) was identified. Gene expression analyses confirmed transcriptional activity and revealed clear differences between susceptible and tolerant genotypes. CONCLUSIONS: Our findings provide new insights into the genomic basis of plant-nematode interactions that can be used to design and accelerate novel management strategies against BCN.

Complete pan-plastome sequences enable high resolution phylogenetic classification of sugar beet and closely related crop wild relatives.

BACKGROUND: As the major source of sugar in moderate climates, sugar-producing beets (Beta vulgaris subsp. vulgaris) have a high economic value. However, the low genetic diversity within cultivated beets requires introduction of new traits, for example to increase their tolerance and resistance attributes - traits that often reside in the crop wild relatives. For this, genetic information of wild beet relatives and their phylogenetic placements to each other are crucial. To answer this need, we sequenced and assembled the complete plastome sequences from a broad species spectrum across the beet genera Beta and Patellifolia, both embedded in the Betoideae (order Caryophyllales). This pan-plastome dataset was then used to determine the wild beet phylogeny in high-resolution. RESULTS: We sequenced the plastomes of 18 closely related accessions representing 11 species of the Betoideae subfamily and provided high-quality plastome assemblies which represent an important resource for further studies of beet wild relatives and the diverse plant order Caryophyllales. Their assembly sizes range from 149,723 bp (Beta vulgaris subsp. vulgaris) to 152,816 bp (Beta nana), with most variability in the intergenic sequences. Combining plastome-derived phylogenies with read-based treatments based on mitochondrial information, we were able to suggest a unified and highly confident phylogenetic placement of the investigated Betoideae species. Our results show that the genus Beta can be divided into the two clearly separated sections Beta and Corollinae. Our analysis confirms the affiliation of B. nana with the other Corollinae species, and we argue against a separate placement in the Nanae section. Within the Patellifolia genus, the two diploid species Patellifolia procumbens and Patellifolia webbiana are, regarding the plastome sequences, genetically more similar to each other than to the tetraploid Patellifolia patellaris. Nevertheless, all three Patellifolia species are clearly separated. CONCLUSION: In conclusion, our wild beet plastome assemblies represent a new resource to understand the molecular base of the beet germplasm. Despite large differences on the phenotypic level, our pan-plastome dataset is highly conserved. For the first time in beets, our whole plastome sequences overcome the low sequence variation in individual genes and provide the molecular backbone for highly resolved beet phylogenomics. Hence, our plastome sequencing strategy can also guide genomic approaches to unravel other closely related taxa.

Functional and evolutionary analysis of the Arabidopsis 4R-MYB protein SNAPc4 as part of the SNAP complex.

Transcription initiation of the genes coding for small nuclear RNA (snRNA) has been extensively analyzed in humans and fruit fly, but only a single ortholog of a snRNA-activating protein complex (SNAPc) subunit has so far been characterized in plants. The genome of the model plant Arabidopsis thaliana encodes orthologs of all three core SNAPc subunits, including A. thaliana SNAP complex 4 (AtSNAPc4)-a 4R-MYB-type protein with four-and-a-half adjacent MYB repeat units. We report the conserved role of AtSNAPc4 as subunit of a protein complex involved in snRNA gene transcription and present genetic evidence that AtSNAPc4 is an essential gene in gametophyte and zygote development. We present experimental evidence that the three A. thaliana SNAPc subunits assemble into a SNAP complex and demonstrate the binding of AtSNAPc4 to snRNA promoters. In addition, co-localization studies show a link between AtSNAPc4 accumulation and Cajal bodies, known to aggregate at snRNA gene loci in humans. Moreover, we show the strong evolutionary conservation of single-copy 4R-MYB/SNAPc4 genes in a broad range of eukaryotes and present additional shared protein features besides the MYB domain, suggesting a conservation of the snRNA transcription initiation machinery along the course of the eukaryotic evolution.

Large scale genomic rearrangements in selected Arabidopsis thaliana T-DNA lines are caused by T-DNA insertion mutagenesis.

BACKGROUND: Experimental proof of gene function assignments in plants is based on mutant analyses. T-DNA insertion lines provided an invaluable resource of mutants and enabled systematic reverse genetics-based investigation of the functions of Arabidopsis thaliana genes during the last decades. RESULTS: We sequenced the genomes of 14 A. thaliana GABI-Kat T-DNA insertion lines, which eluded flanking sequence tag-based attempts to characterize their insertion loci, with Oxford Nanopore Technologies (ONT) long reads. Complex T-DNA insertions were resolved and 11 previously unknown T-DNA loci identified, resulting in about 2 T-DNA insertions per line and suggesting that this number was previously underestimated. T-DNA mutagenesis caused fusions of chromosomes along with compensating translocations to keep the gene set complete throughout meiosis. Also, an inverted duplication of 800 kbp was detected. About 10 % of GABI-Kat lines might be affected by chromosomal rearrangements, some of which do not involve T-DNA. Local assembly of selected reads was shown to be a computationally effective method to resolve the structure of T-DNA insertion loci. We developed an automated workflow to support investigation of long read data from T-DNA insertion lines. All steps from DNA extraction to assembly of T-DNA loci can be completed within days. CONCLUSIONS: Long read sequencing was demonstrated to be an effective way to resolve complex T-DNA insertions and chromosome fusions. Many T-DNA insertions comprise not just a single T-DNA, but complex arrays of multiple T-DNAs. It is becoming obvious that T-DNA insertion alleles must be characterized by exact identification of both T-DNA::genome junctions to generate clear genotype-to-phenotype relations.

Transcriptomic analysis of temporal shifts in berry development between two grapevine cultivars of the Pinot family reveals potential genes controlling ripening time.

BACKGROUND: Grapevine cultivars of the Pinot family represent clonally propagated mutants with major phenotypic and physiological differences, such as different colour or shifted ripening time, as well as changes in important viticultural traits. Specifically, the cultivars 'Pinot Noir' (PN) and 'Pinot Noir Precoce' (PNP, early ripening) flower at the same time, but vary in the beginning of berry ripening (veraison) and, consequently, harvest time. In addition to genotype, seasonal climatic conditions (i.e. high temperatures) also affect ripening times. To reveal possible regulatory genes that affect the timing of veraison onset, we investigated differences in gene expression profiles between PN and PNP throughout berry development with a closely meshed time series and over two separate years. RESULTS: The difference in the duration of berry formation between PN and PNP was quantified to be approximately two weeks under the growth conditions applied, using plant material with a proven PN and PNP clonal relationship. Clusters of co-expressed genes and differentially expressed genes (DEGs) were detected which reflect the shift in the timing of veraison onset. Functional annotation of these DEGs fit to observed phenotypic and physiological changes during berry development. In total, we observed 3,342 DEGs in 2014 and 2,745 DEGs in 2017 between PN and PNP, with 1,923 DEGs across both years. Among these, 388 DEGs were identified as veraison-specific and 12 were considered as berry ripening time regulatory candidates. The expression profiles revealed two candidate genes for ripening time control which we designated VviRTIC1 and VviRTIC2 (VIT_210s0071g01145 and VIT_200s0366g00020, respectively). These genes likely contribute the phenotypic differences observed between PN and PNP. CONCLUSIONS: Many of the 1,923 DEGs show highly similar expression profiles in both cultivars if the patterns are aligned according to developmental stage. In our work, putative genes differentially expressed between PNP and PN which could control ripening time as well as veraison-specific genes were identified. We point out connections of these genes to molecular events during berry development and discuss potential candidate genes which may control ripening time. Two of these candidates were observed to be differentially expressed in the early berry development phase. Several down-regulated genes during berry ripening are annotated as auxin response factors / ARFs. Conceivably, general changes in auxin signaling may cause the earlier ripening phenotype of PNP.

The transition to flowering in winter rapeseed during vernalization.

Flowering time is a major determinant of adaptation, fitness and yield in the allopolyploid species rapeseed (Brassica napus). Despite being a close relative to Arabidopsis thaliana, little is known about the timing of floral transition and the genes that govern this process. Winter, semi-winter and spring type plants have important life history characteristics that differ in vernalization requirements for flowering and are important for growing rapeseed in different regions of the world. In this study, we investigated the timing of vernalization-driven floral transition in winter rapeseed and the effect of photoperiod and developmental age on flowering time and vernalization responsiveness. Microscopy and whole transcriptome analyses at the shoot apical meristems of plants grown under controlled conditions showed that floral transition is initiated within few weeks of vernalization. Certain Bna.SOC1 and Bna.SPL5 homeologs were among the induced genes, suggesting that they are regulating the timing of cold-induced floral transition. Moreover, the flowering response of plants with shorter pre-vernalization period correlated with a delayed expression of Bna.SOC1 and Bna.SPL5 genes. In essence, this study presents a detailed analysis of vernalization-driven floral transition and the aspects of juvenility and dormancy and their effect on flowering time in rapeseed.

The genome of the recently domesticated crop plant sugar beet (Beta vulgaris).

Sugar beet (Beta vulgaris ssp. vulgaris) is an important crop of temperate climates which provides nearly 30% of the world's annual sugar production and is a source for bioethanol and animal feed. The species belongs to the order of Caryophylalles, is diploid with 2n = 18 chromosomes, has an estimated genome size of 714-758 megabases and shares an ancient genome triplication with other eudicot plants. Leafy beets have been cultivated since Roman times, but sugar beet is one of the most recently domesticated crops. It arose in the late eighteenth century when lines accumulating sugar in the storage root were selected from crosses made with chard and fodder beet. Here we present a reference genome sequence for sugar beet as the first non-rosid, non-asterid eudicot genome, advancing comparative genomics and phylogenetic reconstructions. The genome sequence comprises 567 megabases, of which 85% could be assigned to chromosomes. The assembly covers a large proportion of the repetitive sequence content that was estimated to be 63%. We predicted 27,421 protein-coding genes supported by transcript data and annotated them on the basis of sequence homology. Phylogenetic analyses provided evidence for the separation of Caryophyllales before the split of asterids and rosids, and revealed lineage-specific gene family expansions and losses. We sequenced spinach (Spinacia oleracea), another Caryophyllales species, and validated features that separate this clade from rosids and asterids. Intraspecific genomic variation was analysed based on the genome sequences of sea beet (Beta vulgaris ssp. maritima; progenitor of all beet crops) and four additional sugar beet accessions. We identified seven million variant positions in the reference genome, and also large regions of low variability, indicating artificial selection. The sugar beet genome sequence enables the identification of genes affecting agronomically relevant traits, supports molecular breeding and maximizes the plant's potential in energy biotechnology.

Diversification, evolution and methylation of short interspersed nuclear element families in sugar beet and related Amaranthaceae species.

Short interspersed nuclear elements (SINEs) are non-autonomous non-long terminal repeat retrotransposons which are widely distributed in eukaryotic organisms. While SINEs have been intensively studied in animals, only limited information is available about plant SINEs. We analysed 22 SINE families from seven genomes of the Amaranthaceae family and identified 34 806 SINEs, including 19 549 full-length copies. With the focus on sugar beet (Beta vulgaris), we performed a comparative analysis of the diversity, genomic and chromosomal organization and the methylation of SINEs to provide a detailed insight into the evolution and age of Amaranthaceae SINEs. The lengths of consensus sequences of SINEs range from 113 nucleotides (nt) up to 224 nt. The SINEs show dispersed distribution on all chromosomes but were found with higher incidence in subterminal euchromatic chromosome regions. The methylation of SINEs is increased compared with their flanking regions, and the strongest effect is visible for cytosines in the CHH context, indicating an involvement of asymmetric methylation in the silencing of SINEs.

Enhancing the GABI-Kat Arabidopsis thaliana T-DNA Insertion Mutant Database by Incorporating Araport11 Annotation.

SimpleSearch provides access to a database containing information about T-DNA insertion lines of the GABI-Kat collection of Arabidopsis thaliana mutants. These mutants are an important tool for reverse genetics, and GABI-Kat is the second largest collection of such T-DNA insertion mutants. Insertion sites were deduced from flanking sequence tags (FSTs), and the database contains information about mutant plant lines as well as insertion alleles. Here, we describe improvements within the interface (available at http://www.gabi-kat.de/db/genehits.php) and with regard to the database content that have been realized in the last five years. These improvements include the integration of the Araport11 genome sequence annotation data containing the recently updated A. thaliana structural gene descriptions, an updated visualization component that displays groups of insertions with very similar insertion positions, mapped confirmation sequences, and primers. The visualization component provides a quick way to identify insertions of interest, and access to improved data about the exact structure of confirmed insertion alleles. In addition, the database content has been extended by incorporating additional insertion alleles that were detected during the confirmation process, as well as by adding new FSTs that have been produced during continued efforts to complement gaps in FST availability. Finally, the current database content regarding predicted and confirmed insertion alleles as well as primer sequences has been made available as downloadable flat files.

ReadXplorer 2-detailed read mapping analysis and visualization from one single source.

MOTIVATION: The vast amount of already available and currently generated read mapping data requires comprehensive visualization, and should benefit from bioinformatics tools offering a wide spectrum of analysis functionality from just one source. Appropriate handling of multiple mapped reads during mapping analyses remains an issue that demands improvement. RESULTS: The capabilities of the read mapping analysis and visualization tool ReadXplorer were vastly enhanced. Here, we present an even finer granulated read mapping classification, improving the level of detail for analyses and visualizations. The spectrum of automatic analysis functions has been broadened to include genome rearrangement detection as well as correlation analysis between two mapping data sets. Existing functions were refined and enhanced, namely the computation of differentially expressed genes, the read count and normalization analysis and the transcription start site detection. Additionally, ReadXplorer 2 features a highly improved support for large eukaryotic data sets and a command line version, enabling its integration into workflows. Finally, the new version is now able to display any kind of tabular results from other bioinformatics tools. AVAILABILITY AND IMPLEMENTATION: http://www.readxplorer.org CONTACT: readxplorer@computational.bio.uni-giessen.deSupplementary information: Supplementary data are available at Bioinformatics online.

The AtMYB12 activation domain maps to a short C-terminal region of the transcription factor.

The Arabidopsis thaliana R2R3-MYB transcription factor MYB12 is a light-inducible, flavonol-specific activator of flavonoid biosynthesis. The transactivation activity of the AtMYB12 protein was analyzed using a C-terminal deletion series in a transient A. thaliana protoplast assay with the goal of mapping the activation domain (AD). Although the deletion of the last 46 C-terminal amino acids did not affect the activation capacity, the deletion of the last 98 amino acids almost totally abolished transactivation of two different target promoters. A domain swap experiment using the yeast GAL4 DNA-binding domain revealed that the region from positions 282 to 328 of AtMYB12 was sufficient for transactivation. In contrast to the R2R3-MYB ADs known thus far, that of AtMYB12 is not located at the rearmost C-terminal end of the protein. The AtMYB12 AD is conserved in other experimentally proven R2R3-MYB flavonol regulators from different species.

TRANSPARENT TESTA 13 is a tonoplast P3A -ATPase required for vacuolar deposition of proanthocyanidins in Arabidopsis thaliana seeds.

Intracellular pH homeostasis is essential for all living cells. In plants, pH is usually maintained by three structurally distinct and differentially localized types of proton pump: P-type H(+) -ATPases in the plasma membrane, and multimeric vacuolar-type H(+) -ATPases (V-ATPases) and vacuolar H(+) -pyrophosphatases (H(+) -PPases) in endomembranes. Here, we show that reduced accumulation of proanthocyanidins (PAs) and hence the diminished brown seed coloration found in the Arabidopsis thaliana mutant transparent testa 13 (tt13) is caused by disruption of the gene encoding the P3A -ATPase AHA10. Identification of the gene encoded by the tt13 locus completes the molecular characterization of the classical set of transparent testa mutants. Cells of the tt13 seed coat endothelium do not contain PA-filled central vacuoles as observed in the wild-type. tt13 phenocopies tt12, a mutant that is defective in vacuolar import of the PA precursor epicatechin. Our data show that vacuolar loading with PA precursors depends on TT13. Consistent with the tt13 phenotype, but in contrast to other isoforms of P-type H(+) -ATPases, TT13 localizes to the tonoplast. PA accumulation in tt13 is partially restored by expression of the tonoplast localized H(+) -PPase VHP1. Our findings indicate that the P3A -ATPase TT13 functions as a proton pump in the tonoplast of seed coat endothelium cells, and generates the driving force for TT12-mediated transport of PA precursors to the vacuole.

Rapid gene identification in sugar beet using deep sequencing of DNA from phenotypic pools selected from breeding panels.

BACKGROUND: The combination of bulk segregant analysis (BSA) and next generation sequencing (NGS), also known as mapping by sequencing (MBS), has been shown to significantly accelerate the identification of causal mutations for species with a reference genome sequence. The usual approach is to cross homozygous parents that differ for the monogenic trait to address, to perform deep sequencing of DNA from F2 plants pooled according to their phenotype, and subsequently to analyze the allele frequency distribution based on a marker table for the parents studied. The method has been successfully applied for EMS induced mutations as well as natural variation. Here, we show that pooling genetically diverse breeding lines according to a contrasting phenotype also allows high resolution mapping of the causal gene in a crop species. The test case was the monogenic locus causing red vs. green hypocotyl color in Beta vulgaris (R locus). RESULTS: We determined the allele frequencies of polymorphic sequences using sequence data from two diverging phenotypic pools of 180 B. vulgaris accessions each. A single interval of about 31 kbp among the nine chromosomes was identified which indeed contained the causative mutation. CONCLUSIONS: By applying a variation of the mapping by sequencing approach, we demonstrated that phenotype-based pooling of diverse accessions from breeding panels and subsequent direct determination of the allele frequency distribution can be successfully applied for gene identification in a crop species. Our approach made it possible to identify a small interval around the causative gene. Sequencing of parents or individual lines was not necessary. Whenever the appropriate plant material is available, the approach described saves time compared to the generation of an F2 population. In addition, we provide clues for planning similar experiments with regard to pool size and the sequencing depth required.

Repeat Composition of CenH3-chromatin and H3K9me2-marked heterochromatin in Sugar Beet (Beta vulgaris).

BACKGROUND: Sugar beet (Beta vulgaris) is an important crop of temperate climate zones, which provides nearly 30 % of the world's annual sugar needs. From the total genome size of 758 Mb, only 567 Mb were incorporated in the recently published genome sequence, due to the fact that regions with high repetitive DNA contents (e.g. satellite DNAs) are only partially included. Therefore, to fill these gaps and to gain information about the repeat composition of centromeres and heterochromatic regions, we performed chromatin immunoprecipitation followed by sequencing (ChIP-Seq) using antibodies against the centromere-specific histone H3 variant of sugar beet (CenH3) and the heterochromatic mark of dimethylated lysine 9 of histone H3 (H3K9me2). RESULTS: ChIP-Seq analysis revealed that active centromeres containing CenH3 consist of the satellite pBV and the Ty3-gypsy retrotransposon Beetle7, while heterochromatin marked by H3K9me2 exhibits heterogeneity in repeat composition. H3K9me2 was mainly associated with the satellite family pEV, the Ty1-copia retrotransposon family Cotzilla and the DNA transposon superfamily of the En/Spm type. In members of the section Beta within the genus Beta, immunostaining using the CenH3 antibody was successful, indicating that orthologous CenH3 proteins are present in closely related species within this section. CONCLUSIONS: The identification of repetitive genome portions by ChIP-Seq experiments complemented the sugar beet reference sequence by providing insights into the repeat composition of poorly characterized CenH3-chromatin and H3K9me2-heterochromatin. Therefore, our work provides the basis for future research and application concerning the sugar beet centromere and repeat-rich heterochromatic regions characterized by the presence of H3K9me2.

Evolution of tonoplast P-ATPase transporters involved in vacuolar acidification.

Petunia mutants (Petunia hybrida) with blue flowers defined a novel vacuolar proton pump consisting of two interacting P-ATPases, PH1 and PH5, that hyper-acidify the vacuoles of petal cells. PH5 is similar to plasma membrane H(+) P3A -ATPase, whereas PH1 is the only known eukaryoticP3B -ATPase. As there were no indications that this tonoplast pump is widespread in plants, we investigated the distribution and evolution of PH1 and PH5. We combined database mining and phylogenetic and synteny analyses of PH1- and PH5-like proteins from all kingdoms with functional analyses (mutant complementation and intracellular localization) of homologs from diverse angiosperms. We identified functional PH1 and PH5 homologs in divergent angiosperms. PH5 homologs evolved from plasma membrane P3A -ATPases, acquiring an N-terminal tonoplast-sorting sequence and new cellular function before angiosperms appeared. PH1 is widespread among seed plants and related proteins are found in some groups of bacteria and fungi and in one moss, but is absent in most algae, suggesting that its evolution involved several cases of gene loss and possibly horizontal transfer events. The distribution of PH1 and PH5 in the plant kingdom suggests that vacuolar acidification by P-ATPases appeared in gymnosperms before flowers. This implies that, next to flower color determination, vacuolar hyper-acidification is required for yet unknown processes.

Natural variation in flavonol accumulation in Arabidopsis is determined by the flavonol glucosyltransferase BGLU6.

Flavonols are colourless secondary metabolites, primarily regarded as UV-protection pigments that are deposited in plants in their glycosylated forms. The glycosylation of flavonols is mainly catalysed by UDP-sugar-dependent glycosyltransferases (UGTs). Although the structures of flavonol glycosides accumulating in Arabidopsis thaliana are known, many genes involved in the flavonol glycosylation pathway are yet to be discovered. The flavonol glycoside profiles of seedlings from 81 naturally occurring A. thaliana accessions were screened using high performance thin layer chromatography. A qualitative variation in flavonol 3-O-gentiobioside 7-O-rhamnoside (F3GG7R) content was identified. Ler × Col-0 recombinant inbred line mapping and whole genome association mapping led to the identification of a glycoside hydrolase family 1-type gene, At1g60270/BGLU6, that encodes a homolog of acyl-glucose-dependent glucosyltransferases involved in the glycosylation of anthocyanins, possibly localized in the cytoplasm, and that is co-expressed with genes linked to phenylpropanoid biosynthesis. A causal single nucleotide polymorphism introducing a premature stop codon in non-producer accessions was found to be absent in the producers. Several other naturally occurring loss-of-function alleles were also identified. Two independent bglu6 T-DNA insertion mutants from the producer accessions showed loss of F3GG7R. Furthermore, bglu6 mutant lines complemented with the genomic Ler BGLU6 gene confirmed that BGLU6 is essential for production of F3GGR7. We have thus identified an accession-specific gene that causes a qualitative difference in flavonol glycoside accumulation in A. thaliana strains. This gene encodes a flavonol 3-O-glucoside: 6″-O-glucosyltransferase that does not belong to the large canonical family of flavonol glycosyltransferases that use UDP-conjugates as the activated sugar donor substrate.

SMRT sequencing only de novo assembly of the sugar beet (Beta vulgaris) chloroplast genome.

BACKGROUND: Third generation sequencing methods, like SMRT (Single Molecule, Real-Time) sequencing developed by Pacific Biosciences, offer much longer read length in comparison to Next Generation Sequencing (NGS) methods. Hence, they are well suited for de novo- or re-sequencing projects. Sequences generated for these purposes will not only contain reads originating from the nuclear genome, but also a significant amount of reads originating from the organelles of the target organism. These reads are usually discarded but they can also be used for an assembly of organellar replicons. The long read length supports resolution of repetitive regions and repeats within the organelles genome which might be problematic when just using short read data. Additionally, SMRT sequencing is less influenced by GC rich areas and by long stretches of the same base. RESULTS: We describe a workflow for a de novo assembly of the sugar beet (Beta vulgaris ssp. vulgaris) chloroplast genome sequence only based on data originating from a SMRT sequencing dataset targeted on its nuclear genome. We show that the data obtained from such an experiment are sufficient to create a high quality assembly with a higher reliability than assemblies derived from e.g. Illumina reads only. The chloroplast genome is especially challenging for de novo assembling as it contains two large inverted repeat (IR) regions. We also describe some limitations that still apply even though long reads are used for the assembly. CONCLUSIONS: SMRT sequencing reads extracted from a dataset created for nuclear genome (re)sequencing can be used to obtain a high quality de novo assembly of the chloroplast of the sequenced organism. Even with a relatively small overall coverage for the nuclear genome it is possible to collect more than enough reads to generate a high quality assembly that outperforms short read based assemblies. However, even with long reads it is not always possible to clarify the order of elements of a chloroplast genome sequence reliantly which we could demonstrate with Fosmid End Sequences (FES) generated with Sanger technology. Nevertheless, this limitation also applies to short read sequencing data but is reached in this case at a much earlier stage during finishing.

The CHH motif in sugar beet satellite DNA: a modulator for cytosine methylation.

Methylation of DNA is important for the epigenetic silencing of repetitive DNA in plant genomes. Knowledge about the cytosine methylation status of satellite DNAs, a major class of repetitive DNA, is scarce. One reason for this is that arrays of tandemly arranged sequences are usually collapsed in next-generation sequencing assemblies. We applied strategies to overcome this limitation and quantified the level of cytosine methylation and its pattern in three satellite families of sugar beet (Beta vulgaris) which differ in their abundance, chromosomal localization and monomer size. We visualized methylation levels along pachytene chromosomes with respect to small satellite loci at maximum resolution using chromosome-wide fluorescent in situ hybridization complemented with immunostaining and super-resolution microscopy. Only reduced methylation of many satellite arrays was obtained. To investigate methylation at the nucleotide level we performed bisulfite sequencing of 1569 satellite sequences. We found that the level of methylation of cytosine strongly depends on the sequence context: cytosines in the CHH motif show lower methylation (44-52%), while CG and CHG motifs are more strongly methylated. This affects the overall methylation of satellite sequences because CHH occurs frequently while CG and CHG are rare or even absent in the satellite arrays investigated. Evidently, CHH is the major target for modulation of the cytosine methylation level of adjacent monomers within individual arrays and contributes to their epigenetic function. This strongly indicates that asymmetric cytosine methylation plays a role in the epigenetic modification of satellite repeats in plant genomes.