Integration of fluorescence in situ hybridization and chromosome-length genome assemblies revealed synteny map for guinea pig, naked mole-rat, and human

Descriptions of karyotypes of many animal species are currently available. In addition, there has been a significant increase in the number of sequenced genomes and an ever-improving quality of genome assembly. To close the gap between genomic and cytogenetic data we applied fluorescent in situ hybridization (FISH) and Hi-C technology to make the first full chromosome-level genome comparison of the guinea pig (Cavia porcellus), naked mole-rat (Heterocephalus glaber), and human. Comparative chromosome maps obtained by FISH with chromosome-specific probes link genomic scaffolds to individual chromosomes and orient them relative to centromeres and heterochromatic blocks. Hi-C assembly made it possible to close all gaps on the comparative maps and to reveal additional rearrangements that distinguish the karyotypes of the three species. As a result, we integrated the bioinformatic and cytogenetic data and adjusted the previous comparative maps and genome assemblies of the guinea pig, naked mole-rat, and human. Syntenic associations in the two hystricomorphs indicate features of their putative ancestral karyotype. We postulate that the two approaches applied in this study complement one another and provide complete information about the organization of these genomes at the chromosome level.

Cytogenetic screening of a Canadian swine breeding nucleus using a newly developed karyotyping method named oligo-banding

BACKGROUND

The frequency of chromosomal rearrangements in Canadian breeding boars has been estimated at 0.91 to 1.64%. These abnormalities are widely recognized as a potential cause of subfertility in livestock production. Since artificial insemination is practiced in almost all intensive pig production systems, the use of elite boars carrying cytogenetic defects that have an impact on fertility can lead to major economic losses. To avoid keeping subfertile boars in artificial insemination centres and spreading chromosomal defects within populations, cytogenetic screening of boars is crucial. Different techniques are used for this purpose, but several issues are frequently encountered, i.e. environmental factors can influence the quality of results, the lack of genomic information outputted by these techniques, and the need for prior cytogenetic skills. The aim of this study was to develop a new pig karyotyping method based on fluorescent banding patterns.

RESULTS

The use of 207,847 specific oligonucleotides generated 96 fluorescent bands that are distributed across the 18 autosomes and the sex chromosomes. Tested alongside conventional G-banding, this oligo-banding method allowed us to identify four chromosomal translocations and a rare unbalanced chromosomal rearrangement that was not detected by conventional banding. In addition, this method allowed us to investigate chromosomal imbalance in spermatozoa.

CONCLUSIONS

The use of oligo-banding was found to be appropriate for detecting chromosomal aberrations in a Canadian pig nucleus and its convenient design and use make it an interesting tool for livestock karyotyping and cytogenetic studies.

Chromosome-level assembly of the Rangifer tarandus genome and validation of cervid and bovid evolution insights

BACKGROUND

Genome assembly into chromosomes facilitates several analyses including cytogenetics, genomics and phylogenetics. Despite rapid development in bioinformatics, however, assembly beyond scaffolds remains challenging, especially in species without closely related well-assembled and available reference genomes. So far, four draft genomes of Rangifer tarandus (caribou or reindeer, a circumpolar distributed cervid species) have been published, but none with chromosome-level assembly. This emblematic northern species is of high interest in ecological studies and conservation since most populations are declining.

RESULTS

We have designed specific probes based on Oligopaint FISH technology to upgrade the latest published reindeer and caribou chromosome-level genomes. Using this oligonucleotide-based method, we found six mis-assembled scaffolds and physically mapped 68 of the largest scaffolds representing 78% of the most recent R. tarandus genome assembly. Combining physical mapping and comparative genomics, it was possible to document chromosomal evolution among Cervidae and closely related bovids.

CONCLUSIONS

Our results provide validation for the current chromosome-level genome assembly as well as resources to use chromosome banding in studies of Rangifer tarandus.

Nitrous Oxide-Induced Metaphase Arrest: A Technique for Somatic Chromosome Analysis

Procedures to arrest metaphase chromosomes are used for determining chromosome numbers, chromosomal aberrations, and natural chromosome variation, as well as chromosome sorting. Here is described a technique of nitrous oxide gas treatment of freshly harvested root tips that is highly effective at producing an excellent mitotic index together with well-spread chromosomes. The details of the treatment and equipment used are provided. The metaphase spreads can be used directly for determining chromosome numbers or for in situ hybridization to reveal chromosomal features.

Evolution and co-evolution: insights into the divergence of plant heat shock factor genes

The Heat Shock Factor (Hsf) genes are widely distributed across the plant kingdom regulating the plant response to various abiotic stresses. In addition to natural selection, breeding and accelerated selection changed the structure and function of Hsf genes. 1076 Hsf genes from 30 genera from primitive algae to the most advanced plant species and major crop plants were used for phylogenetic analysis. The interspecific divergence was studied with 11 members of genus Oryza while intraspecific divergence was studied with sesame pan-genome adapted to diverse ecological niches. B2 genes in eudicots and monocots originated separately while A1 gave rise to the recently evolved Class-C genes and land colonization happened with evolution of A1 genes. An increase in the number of lineages in the Oryza clade with the evolution of AA genome indicated independent domestication and positive selection was observed in > 53% of loci whereas the highly conserved homologues were under purifying selection. The paralogous genes under positive selection exhibited more domain changes for diversified function and increased fitness. A significant co-evolving cluster involving amino acids Phenylalanine, Lysine and Valine played crucial role in maintaining hydrophobic core along with highly conserved Tryptophan residues. A mutation of Glutamic acid to Glutamine was observed in A8 genes of Lamiales affecting protein solvency. Breeding resulted in accumulation of mutations reducing the hydrophobicity of proteins and a further reduction in protein aggregation. This study identify genome duplications, non-neutral selection and co-evolving residues as causing drastic changes in the conserved domain of Hsf proteins.

Supplementary information

The online version contains supplementary material available at 10.1007/s12298-022-01183-7.

Assembled and annotated 26.5 Gbp coast redwood genome: a resource for estimating evolutionary adaptive potential and investigating hexaploid origin

Sequencing, assembly, and annotation of the 26.5 Gbp hexaploid genome of coast redwood (Sequoia sempervirens) was completed leading toward discovery of genes related to climate adaptation and investigation of the origin of the hexaploid genome. Deep-coverage short-read Illumina sequencing data from haploid tissue from a single seed were combined with long-read Oxford Nanopore Technologies sequencing data from diploid needle tissue to create an initial assembly, which was then scaffolded using proximity ligation data to produce a highly contiguous final assembly, SESE 2.1, with a scaffold N50 size of 44.9 Mbp. The assembly included several scaffolds that span entire chromosome arms, confirmed by the presence of telomere and centromere sequences on the ends of the scaffolds. The structural annotation produced 118,906 genes with 113 containing introns that exceed 500 Kbp in length and one reaching 2 Mb. Nearly 19 Gbp of the genome represented repetitive content with the vast majority characterized as long terminal repeats, with a 2.9:1 ratio of Copia to Gypsy elements that may aid in gene expression control. Comparison of coast redwood to other conifers revealed species-specific expansions for a plethora of abiotic and biotic stress response genes, including those involved in fungal disease resistance, detoxification, and physical injury/structural remodeling and others supporting flavonoid biosynthesis. Analysis of multiple genes that exist in triplicate in coast redwood but only once in its diploid relative, giant sequoia, supports a previous hypothesis that the hexaploidy is the result of autopolyploidy rather than any hybridizations with separate but closely related conifer species.

Gramene: A Resource for Comparative Analysis of Plants Genomes and Pathways

Gramene is an integrated bioinformatics resource for accessing, visualizing, and comparing plant genomes and biological pathways. Originally targeting grasses, Gramene has grown to host annotations for over 90 plant genomes including agronomically important cereals (e.g., maize, sorghum, wheat, teff), fruits and vegetables (e.g., apple, watermelon, clementine, tomato, cassava), specialty crops (e.g., coffee, olive tree, pistachio, almond), and plants of special or emerging interest (e.g., cotton, tobacco, cannabis, or hemp). For some species, the resource includes multiple varieties of the same species, which has paved the road for the creation of species-specific pan-genome browsers. The resource also features plant research models, including Arabidopsis and C4 warm-season grasses and brassicas, as well as other species that fill phylogenetic gaps for plant evolution studies. Its strength derives from the application of a phylogenetic framework for genome comparison and the use of ontologies to integrate structural and functional annotation data. This chapter outlines system requirements for end-users and database hosting, data types and basic navigation within Gramene, and provides examples of how to (1) explore Gramene's search results, (2) explore gene-centric comparative genomics data visualizations in Gramene, and (3) explore genetic variation associated with a gene locus. This is the first publication describing in detail Gramene's integrated search interface-intended to provide a simplified entry portal for the resource's main data categories (genomic location, phylogeny, gene expression, pathways, and external references) to the most complete and up-to-date set of plant genome and pathway annotations.

A near-complete genome sequence of mungbean (Vigna radiata L.) provides key insights into the modern breeding program

Mungbean (Vigna radiata L.), a fast-growing legume species, is an important source of carbohydrates and proteins in developing countries of Asia. Here, we constructed a near-complete genome sequence of mungbean with a scaffold N50 value of 5.2 Mb and only a 0.4% gap, with a total scaffold size of 475 Mb. We identified several misassembled pseudomolecules (Chr03, Chr04, Chr05, and Chr08) in the previous draft assembly; Chr03, Chr04, and Chr08 were assembled into one chromosome, and Chr05 was broken into two chromosomes in the improved reference genome assembly, thus providing more accurate linkage information to breeders. Additionally, using an ultra-high-resolution linkage map constructed based on resequencing data, we identified several quantitative trait loci (QTLs) and the underlying candidate genes affecting synchronous pod maturity (SPM). Mungbean homologs of two soybean ([Glycine max (L.) Merr.] flowering genes, E3 (phytochrome A) and J (early flowering 3), were identified as candidate genes for the QTLs, and the candidate genes for plant height, node number, and SPM showed critical nucleotide substitutions between the reference cultivar and other genotypes (landraces and wild accessions). Based on the analysis of genetic diversity among 276 accessions collected from 23 countries, we identified 36 selective sweep regions and observed that the overall genetic diversity of cultivars decreased to 30% of that in wild accessions postdomestication. The near-complete genome sequence of mungbean represents an important resource for genome-assisted improvement in the mungbean breeding program.

Conditional stomatal closure in a fern shares molecular features with flowering plant active stomatal responses

Stomata evolved as plants transitioned from water to land, enabling carbon dioxide uptake and water loss to be controlled. In flowering plants, the most recently divergent land plant lineage, stomatal pores actively close in response to drought. In this response, the phytohormone abscisic acid (ABA) triggers signaling cascades that lead to ion and water loss in the guard cells of the stomatal complex, causing a reduction in turgor and pore closure. Whether this stimulus-response coupling pathway acts in other major land plant lineages is unclear, with some investigations reporting that stomatal closure involves ABA but others concluding that closure is passive. Here, we show that in the model fern Ceratopteris richardii active stomatal closure is conditional on sensitization by pre-exposure to either low humidity or exogenous ABA and is promoted by ABA. RNA-seq analysis and de novo transcriptome assembly reconstructed the protein-coding complement of the C. richardii genome, with coverage comparable to other plant models, enabling transcriptional signatures of stomatal sensitization and closure to be inferred. In both cases, changes in abundance of homologs of ABA, Ca²⁺, and ROS-related signaling components were observed, suggesting that the closure-response pathway is conserved in ferns and flowering plants. These signatures further suggested that sensitization is achieved by lowering the threshold required for a subsequent closure-inducing signal to trigger a response. We conclude that the canonical signaling network for active stomatal closure functioned in at least a rudimentary form in the stomata of the last common ancestor of ferns and flowering plants.

Fast and efficient Rmap assembly using the Bi-labelled de Bruijn graph

Genome wide optical maps are high resolution restriction maps that give a unique numeric representation to a genome. They are produced by assembling hundreds of thousands of single molecule optical maps, which are called Rmaps. Unfortunately, there are very few choices for assembling Rmap data. There exists only one publicly-available non-proprietary method for assembly and one proprietary software that is available via an executable. Furthermore, the publicly-available method, by Valouev et al. (Proc Natl Acad Sci USA 103(43):15770-15775, 2006), follows the overlap-layout-consensus (OLC) paradigm, and therefore, is unable to scale for relatively large genomes. The algorithm behind the proprietary method, Bionano Genomics' Solve, is largely unknown. In this paper, we extend the definition of bi-labels in the paired de Bruijn graph to the context of optical mapping data, and present the first de Bruijn graph based method for Rmap assembly. We implement our approach, which we refer to as RMAPPER, and compare its performance against the assembler of Valouev et al. (Proc Natl Acad Sci USA 103(43):15770-15775, 2006) and Solve by Bionano Genomics on data from three genomes: E. coli, human, and climbing perch fish (Anabas Testudineus). Our method was able to successfully run on all three genomes. The method of Valouev et al. (Proc Natl Acad Sci USA 103(43):15770-15775, 2006) only successfully ran on E. coli. Moreover, on the human genome RMAPPER was at least 130 times faster than Bionano Solve, used five times less memory and produced the highest genome fraction with zero mis-assemblies. Our software, RMAPPER is written in C++ and is publicly available under GNU General Public License at https://github.com/kingufl/Rmapper .

The Histone Acetyltransferase GCN5 and the Associated Coactivators ADA2: From Evolution of the SAGA Complex to the Biological Roles in Plants

Transcription of protein-encoding genes starts with forming a pre-initiation complex comprised of RNA polymerase II and several general transcription factors. To activate gene expression, transcription factors must overcome repressive chromatin structure, which is accomplished with multiprotein complexes. One such complex, SAGA, modifies the nucleosomal histones through acetylation and other histone modifications. A prototypical histone acetyltransferase (HAT) known as general control non-repressed protein 5 (GCN5), was defined biochemically as the first transcription-linked HAT with specificity for histone H3 lysine 14. In this review, we analyze the components of the putative plant SAGA complex during plant evolution, and current knowledge on the biological role of the key components of the HAT module, GCN5 and ADA2b in plants, will be summarized.

Chromosome-scale genome assembly for the duckweed Spirodela intermedia, integrating cytogenetic maps, PacBio and Oxford Nanopore libraries

Duckweeds are small, free-floating, morphologically highly reduced organisms belonging to the monocot order Alismatales. They display the most rapid growth among flowering plants, vary ~ 14-fold in genome size and comprise five genera. Spirodela is the phylogenetically oldest genus with only two mainly asexually propagating species: S. polyrhiza (2n = 40; 160 Mbp/1C) and S. intermedia (2n = 36; 160 Mbp/1C). This study combined comparative cytogenetics and de novo genome assembly based on PacBio, Illumina and Oxford Nanopore (ON) reads to obtain the first genome reference for S. intermedia and to compare its genomic features with those of the sister species S. polyrhiza. Both species' genomes revealed little more than 20,000 putative protein-coding genes, very low rDNA copy numbers and a low amount of repetitive sequences, mainly Ty3/gypsy retroelements. The detection of a few new small chromosome rearrangements between both Spirodela species refined the karyotype and the chromosomal sequence assignment for S. intermedia.

The Evolution of Chromosome Numbers: Mechanistic Models and Experimental Approaches

Chromosome numbers have been widely used to describe the most fundamental genomic attribute of an organism or a lineage. Although providing strong phylogenetic signal, chromosome numbers vary remarkably among eukaryotes at all levels of taxonomic resolution. Changes in chromosome numbers regularly serve as indication of major genomic events, most notably polyploidy and dysploidy. Here, we review recent advancements in our ability to make inferences regarding historical events that led to alterations in the number of chromosomes of a lineage. We first describe the mechanistic processes underlying changes in chromosome numbers, focusing on structural chromosomal rearrangements. Then, we focus on experimental procedures, encompassing comparative cytogenomics and genomics approaches, and on computational methodologies that are based on explicit models of chromosome-number evolution. Together, these tools offer valuable predictions regarding historical events that have changed chromosome numbers and genome structures, as well as their phylogenetic and temporal placements.

Kohdista: an efficient method to index and query possible Rmap alignments

Background

Genome-wide optical maps are ordered high-resolution restriction maps that give the position of occurrence of restriction cut sites corresponding to one or more restriction enzymes. These genome-wide optical maps are assembled using an overlap-layout-consensus approach using raw optical map data, which are referred to as Rmaps. Due to the high error-rate of Rmap data, finding the overlap between Rmaps remains challenging.

Results

We present Kohdista, which is an index-based algorithm for finding pairwise alignments between single molecule maps (Rmaps). The novelty of our approach is the formulation of the alignment problem as automaton path matching, and the application of modern index-based data structures. In particular, we combine the use of the Generalized Compressed Suffix Array (GCSA) index with the wavelet tree in order to build Kohdista. We validate Kohdista on simulated E. coli data, showing the approach successfully finds alignments between Rmaps simulated from overlapping genomic regions.

Conclusion

we demonstrate Kohdista is the only method that is capable of finding a significant number of high quality pairwise Rmap alignments for large eukaryote organisms in reasonable time.

The C-Fern (Ceratopteris richardii) genome: insights into plant genome evolution with the first partial homosporous fern genome assembly

Ferns are notorious for possessing large genomes and numerous chromosomes. Despite decades of speculation, the processes underlying the expansive genomes of ferns are unclear, largely due to the absence of a sequenced homosporous fern genome. The lack of this crucial resource has not only hindered investigations of evolutionary processes responsible for the unusual genome characteristics of homosporous ferns, but also impeded synthesis of genome evolution across land plants. Here, we used the model fern species Ceratopteris richardii to address the processes (e.g., polyploidy, spread of repeat elements) by which the large genomes and high chromosome numbers typical of homosporous ferns may have evolved and have been maintained. We directly compared repeat compositions in species spanning the green plant tree of life and a diversity of genome sizes, as well as both short- and long-read-based assemblies of Ceratopteris. We found evidence consistent with a single ancient polyploidy event in the evolutionary history of Ceratopteris based on both genomic and cytogenetic data, and on repeat proportions similar to those found in large flowering plant genomes. This study provides a major stepping-stone in the understanding of land plant evolutionary genomics by providing the first homosporous fern reference genome, as well as insights into the processes underlying the formation of these massive genomes.

Re-sequencing and optical mapping reveals misassemblies and real inversions on Corynebacterium pseudotuberculosis genomes

The number of draft genomes deposited in Genbank from the National Center for Biotechnology Information (NCBI) is higher than the complete ones. Draft genomes are assemblies that contain fragments of misassembled regions (gaps). Such draft genomes present a hindrance to the complete understanding of the biology and evolution of the organism since they lack genomic information. To overcome this problem, strategies to improve the assembly process are developed continuously. Also, the greatest challenge to the assembly progress is the presence of repetitive DNA regions. This article highlights the use of optical mapping, to detect and correct assembly errors in Corynebacterium pseudotuberculosis. We also demonstrate that choosing a reference genome should be done with caution to avoid assembly errors and loss of genetic information.

Mechanistic insights into the evolution of DUF26-containing proteins in land plants

Large protein families are a prominent feature of plant genomes and their size variation is a key element for adaptation. However, gene and genome duplications pose difficulties for functional characterization and translational research. Here we infer the evolutionary history of the DOMAIN OF UNKNOWN FUNCTION (DUF) 26-containing proteins. The DUF26 emerged in secreted proteins. Domain duplications and rearrangements led to the appearance of CYSTEINE-RICH RECEPTOR-LIKE PROTEIN KINASES (CRKs) and PLASMODESMATA-LOCALIZED PROTEINS (PDLPs). The DUF26 is land plant-specific but structural analyses of PDLP ectodomains revealed strong similarity to fungal lectins and thus may constitute a group of plant carbohydrate-binding proteins. CRKs expanded through tandem duplications and preferential retention of duplicates following whole genome duplications, whereas PDLPs evolved according to the dosage balance hypothesis. We propose that new gene families mainly expand through small-scale duplications, while fractionation and genetic drift after whole genome multiplications drive families towards dosage balance.

Shining a Spotlight on DNA: Single-Molecule Methods to Visualise DNA

The ability to watch single molecules of DNA has revolutionised how we study biological transactions concerning nucleic acids. Many strategies have been developed to manipulate DNA molecules to investigate mechanical properties, dynamics and protein–DNA interactions. Imaging methods using small molecules and protein-based probes to visualise DNA have propelled our understanding of complex biochemical reactions involving DNA. This review focuses on summarising some of the methodological developments made to visualise individual DNA molecules and discusses how these probes have been used in single-molecule biophysical assays.

Error correcting optical mapping data

Optical mapping is a unique system that is capable of producing high-resolution, high-throughput genomic map data that gives information about the structure of a genome . Recently it has been used for scaffolding contigs and for assembly validation for large-scale sequencing projects, including the maize, goat, and Amborella genomes. However, a major impediment in the use of this data is the variety and quantity of errors in the raw optical mapping data, which are called Rmaps. The challenges associated with using Rmap data are analogous to dealing with insertions and deletions in the alignment of long reads. Moreover, they are arguably harder to tackle since the data are numerical and susceptible to inaccuracy. We develop cOMet to error correct Rmap data, which to the best of our knowledge is the only optical mapping error correction method. Our experimental results demonstrate that cOMet has high prevision and corrects 82.49% of insertion errors and 77.38% of deletion errors in Rmap data generated from the Escherichia coli K-12 reference genome. Out of the deletion errors corrected, 98.26% are true errors. Similarly, out of the insertion errors corrected, 82.19% are true errors. It also successfully scales to large genomes, improving the quality of 78% and 99% of the Rmaps in the plum and goat genomes, respectively. Last, we show the utility of error correction by demonstrating how it improves the assembly of Rmap data. Error corrected Rmap data results in an assembly that is more contiguous and covers a larger fraction of the genome.

Generating a high-confidence reference genome map of the Greater Duckweed by integration of cytogenomic, optical mapping, and Oxford Nanopore technologies

Duckweeds are the fastest growing angiosperms and have the potential to become a new generation of sustainable crops. Although a seed plant, Spirodela polyrhiza clones rarely flower and multiply mainly through vegetative propagation. Whole-genome sequencing using different approaches and clones yielded two reference maps. One for clone 9509, supported in its assembly by optical mapping of single DNA molecules, and one for clone 7498, supported by cytogenetic assignment of 96 fingerprinted bacterial artificial chromosomes (BACs) to its 20 chromosomes. However, these maps differ in the composition of several individual chromosome models. We validated both maps further to resolve these differences and addressed whether they could be due to chromosome rearrangements in different clones. For this purpose, we applied sequential multicolor fluorescence in situ hybridization (mcFISH) to seven S. polyrhiza clones, using 106 BACs that were mapped onto the 39 pseudomolecules for clone 7498. Furthermore we integrated high-depth Oxford Nanopore (ON) sequence data for clone 9509 to validate and revise the previously assembled chromosome models. We found no major structural rearrangements between these seven clones, identified seven chimeric pseudomolecules and Illumina assembly errors in the previous maps, respectively. A new S. polyrhiza genome map with high contiguity was produced with the ON sequence data and genome-wide synteny analysis supported the occurrence of two Whole Genome Duplication events during its evolution. This work generated a high confidence genome map for S. polyrhiza at the chromosome scale, and illustrates the complementarity of independent approaches to produce whole-genome assemblies in the absence of a genetic map.

Members of the Arabidopsis FORKED1-LIKE gene family act to localize PIN1 in developing veins

The reticulate leaf vein pattern typical of angiosperms is proposed to have been a driving force for their evolutionary success. Vein pattern is established through auxin canalization via the auxin efflux protein PINFORMED1 (PIN1). During formation of vein loops, PIN1 cellular localization is increasingly restricted to either the basal side of cells in the lower domain or to the apical side in the upper domain. We previously identified the gene FORKED1 (FKD1) to be required for PIN1 asymmetric localization and for the formation of closed vein loops. FKD1 encodes a plant-specific protein with a domain of unknown function (DUF828) and a Pleckstrin-like homology domain. The Arabidopsis genome encodes eight similar proteins, which we term the FORKED1-LIKE (FL) gene family. Five FL family members localize primarily to the trans-Golgi network or the Golgi, and several co-localize with FKD1-green flourescent protein (GFP) and RABA1c, suggesting action in the secretory pathway. While single FL gene family mutations do not result in vein pattern defects, triple mutants with mutations in FKD1, FL2, and FL3 result in a more symmetric PIN1 localization and a highly disconnected vein pattern. Our data suggest that FL genes act redundantly with FKD1 in the secretory pathway to establish appropriate PIN1 localization in provascular tissue.

Rapid and Efficient FISH using Pre-Labeled Oligomer Probes

Fluorescence in situ hybridization (FISH) is used to visualize the distribution of DNA elements within a genome. Conventional methods for FISH take 1-2 days. Here, we developed a simplified, rapid FISH technique using pre-labeled oligonucleotide probes (PLOPs) and tested the procedure using 18 PLOPs from 45S and 5S rDNA, Arabidopsis-type telomere, and newly-identified Panax ginseng-specific tandem repeats. The 16 developed rDNA PLOPs can be universally applied to plants and animals. The telomere PLOPs can be utilized in most plants with Arabidopsis-type telomeres. The ginseng-specific PLOP can be used to distinguish P. ginseng from related Panax species. Differential labeling of PLOPs allowed us to simultaneously visualize different target loci while reducing the FISH hybridization time from ~16 h to 5 min. PLOP-FISH is efficient, reliable, and rapid, making it ideal for routine analysis, especially of newly sequenced genomes using either universal or specific targets, such as novel tandem repeats identified from whole-genome sequencing data.

Functional diversification of Arabidopsis SEC1-related SM proteins in cytokinetic and secretory membrane fusion

Sec1/Munc18 (SM) proteins contribute to membrane fusion by interacting with Qa-SNAREs or nascent trans-SNARE complexes. Gymnosperms and the basal angiosperm Amborella have only a single SEC1 gene related to the KEULE gene in Arabidopsis However, the genomes of most angiosperms including Arabidopsis encode three SEC1-related SM proteins of which only KEULE has been functionally characterized as interacting with the cytokinesis-specific Qa-SNARE KNOLLE during cell-plate formation. Here we analyze the closest paralog of KEULE named SEC1B. In contrast to the cytokinesis defects of keule mutants, sec1b mutants are homozygous viable. However, the keule sec1b double mutant was nearly gametophytically lethal, displaying collapsed pollen grains, which suggests substantial overlap between SEC1B and KEULE functions in secretion-dependent growth. SEC1B had a strong preference for interaction with the evolutionarily ancient Qa-SNARE SYP132 involved in secretion and cytokinesis, whereas KEULE interacted with both KNOLLE and SYP132. This differential interaction with Qa-SNAREs is likely conferred by domains 1 and 2a of the two SM proteins. Comparative analysis of all four possible combinations of the relevant SEC1 Qa-SNARE double mutants revealed that in cytokinesis, the interaction of SEC1B with KNOLLE plays no role, whereas the interaction of KEULE with KNOLLE is prevalent and functionally as important as the interactions of both SEC1B and KEU with SYP132 together. Our results suggest that functional diversification of the two SEC1-related SM proteins during angiosperm evolution resulted in enhanced interaction of SEC1B with Qa-SNARE SYP132, and thus a predominant role of SEC1B in secretion.

Convergent adaptive evolution in marginal environments: unloading transposable elements as a common strategy among mangrove genomes

Several clades of mangrove trees independently invade the interface between land and sea at the margin of woody plant distribution. As phenotypic convergence among mangroves is common, the possibility of convergent adaptation in their genomes is quite intriguing. To study this molecular convergence, we sequenced multiple mangrove genomes. In this study, we focused on the evolution of transposable elements (TEs) in relation to the genome size evolution. TEs, generally considered genomic parasites, are the most common components of woody plant genomes. Analyzing the long terminal repeat-retrotransposon (LTR-RT) type of TE, we estimated their death rates by counting solo-LTRs and truncated elements. We found that all lineages of mangroves massively and convergently reduce TE loads in comparison to their nonmangrove relatives; as a consequence, genome size reduction happens independently in all six mangrove lineages; TE load reduction in mangroves can be attributed to the paucity of young elements; the rarity of young LTR-RTs is a consequence of fewer births rather than access death. In conclusion, mangrove genomes employ a convergent strategy of TE load reduction by suppressing element origination in their independent adaptation to a new environment.

Visible laser-assisted reduction of plasmonic Ag nanoparticles with narrow-band optical absorption for colored holographic reconstruction

Noble metal plasmonic resonance has been utilized in optical data storage widely for its excellent photo-transformation efficiency. TiO₂ nanoporous films deposited with Ag nanoparticles present outstanding polarization-response and color-modulation ability. However, the low exposure-sensitivity at single wavelength inhibits their application in optical information processing, which is urgent to be improved by innovative methods. Here, we report that Ag nanoparticles were deposited efficiently via continuous laser irradiation in the TiO₂ nanoporous film treated by tannic acid, presenting high-efficient monochromic absorption property. As a result, two sets of holograms were recorded sequentially at the same point of the film with orthogonal circular polarization configurations. The colored reconstruction of the mixed holograms was achieved by utilizing laser polarization state as chrominance segmentation channel. Our method provides a distinctive route for enhancing the photo-energy conversion efficiency of plasmonic nanoparticles, and paves a way to develop advanced display device.

Mapping Genomic Scaffolds to Chromosomes Using Laser Capture Microdissection in Application to Hawaiian Picture-Winged Drosophila

Next-generation sequencing technologies have led to a decreased cost and an increased throughput in genome sequencing. Yet, many genome assemblies based on short sequencing reads have been assembled only to the scaffold level due to the lack of sufficient chromosome mapping information. Traditional ways of mapping scaffolds to chromosomes require a large amount of laboratory work and time to generate genetic and/or physical maps. To address this problem, we conducted a rapid technique which uses laser capture microdissection and enables mapping scaffolds of de novo genome assemblies directly to chromosomes in Hawaiian picture-winged Drosophila. We isolated and sequenced intact chromosome arms from larvae of D. differens. By mapping the reads of each chromosome to the recently assembled scaffolds from 3 Hawaiian picture-winged Drosophila species, at least 67% of the scaffolds were successfully assigned to chromosome arms. Even though the scaffolds are not ordered within a chromosome, the fast-generated chromosome information allows for chromosome-related analyses after genome assembling. We utilize this new information to test the faster-X evolution effect for the first time in these Hawaiian picture-winged Drosophila species.

A comparative transcriptomic analysis reveals the core genetic components of salt and osmotic stress responses in Braya humilis

Braya humilis is a member of the Euclidieae tribe within the family Brassicaceae. This species exhibits a broad range of adaptations to different climatic zones and latitudes as it has a distribution that ranges from northern Asia to the arctic-alpine regions of northern North America. In China, B. humilis is mainly found on the Qinghai-Tibetan Plateau (QTP) and in adjacent arid regions. In this study, we sequenced a sample from an arid region adjacent to the QTP using the Illumina platform generating a total of 46,485 highly accurate unigenes, of which 78.41% were annotated by BLASTing versus public protein databases. The B. humilis transcriptome is characterized by a high level of sequence conservation compared with its close relative, Arabidopsis thaliana. We also used reciprocal blast to identify shared orthologous genes between B. humilis and four other sequenced Brassicaceae species (i.e. A. thaliana, A. lyrata, Capsella rubella, and Thellungiella parvula). To enable precise characterization of orthologous genes, the early-diverging basal angiosperm Amborella trichopoda was also included. A total of 6,689 orthologous genes were identified before stricter criteria for the determination of e-values, amino acid hit lengths, and identity values was applied to further reduce this list. This led to a final list of 381 core orthologous genes for B. humilis; 39 out of these genes are involved in salt and osmotic stress responses and estimations of nonsynonymous/synonymous substitution ratios for this species and A. thaliana orthologs show that these genes are under purifying selection in B. humilis. Expression of six genes was detected in B. humilis seedlings under salt and osmotic stress treatments. Comparable expression patterns to their counterparts in Arabidopsis suggest that these orthologous genes are both sequence and functional conservation. The results of this study demonstrate that the environmental adaptations of B. humilis are mainly the results of preexisting genetic components. Future work will be required to characterize the expression patterns of these orthologous genes in natural populations and will provide further insights into the adaptive mechanisms underlying the wide range of B. humilis adaptations.

Reconstruction of chromosome rearrangements between the two most ancestral duckweed species Spirodela polyrhiza and S. intermedia

The monophyletic duckweeds comprising five genera within the monocot order Alismatales are neotenic, free-floating, aquatic organisms with fast vegetative propagation. Some species are considered for efficient biomass production, for life stock feeding, and for (simultaneous) wastewater phytoremediation. The ancestral genus Spirodela consists of only two species, Spirodela polyrhiza and Spirodela intermedia, both with a similar small genome (~160 Mbp/1C). Reference genome drafts and a physical map of 96 BACs on the 20 chromosome pairs of S. polyrhiza strain 7498 are available and provide useful tools for further evolutionary studies within and between duckweed genera. Here we applied sequential comparative multicolor fluorescence in situ hybridization (mcFISH) to address homeologous chromosomes in S. intermedia (2n = 36), to detect chromosome rearrangements between both species and to elucidate the mechanisms which may have led to the chromosome number alteration after their evolutionary separation. Ten chromosome pairs proved to be conserved between S. polyrhiza and S. intermedia, the remaining ones experienced, depending on the assumed direction of evolution, translocations, inversion, and fissions, respectively. These results represent a first step to unravel karyotype evolution among duckweeds and are anchor points for future genome assembly of S. intermedia.

WUSCHEL-RELATED HOMEOBOX4 (WOX4)-like genes regulate cambial cell division activity and secondary growth in Populus trees

Plant secondary growth derives from the meristematic activity of the vascular cambium. In Arabidopsis thaliana, cell divisions in the cambium are regulated by the transcription factor WOX4, a key target of the CLAVATA3 (CLV3)/EMBRYO SURROUNDING REGION (ESR)-RELATED 41 (CLE41) signaling pathway. However, function of the WOX4-like genes in plants that are dependent on a much more prolific secondary growth, such as trees, remains unclear. Here, we investigate the role of WOX4 and CLE41 homologs for stem secondary growth in Populus trees. In Populus, PttWOX4 genes are specifically expressed in the cambial region during vegetative growth, but not after growth cessation and during dormancy, possibly involving a regulation by auxin. In PttWOX4a/b RNAi trees, primary growth was not affected whereas the width of the vascular cambium was severely reduced and secondary growth was greatly diminished. Our data show that in Populus trees, PttWOX4 genes control cell division activity in the vascular cambium, and hence growth in stem girth. This activity involves the positive regulation of PttWOX4a/b through PttCLE41-related genes. Finally, expression profiling suggests that the CLE41 signaling pathway is an evolutionarily conserved program for the regulation of vascular cambium activity between angiosperm and gymnosperm tree species.

Noise reduction in single time frame optical DNA maps

In optical DNA mapping technologies sequence-specific intensity variations (DNA barcodes) along stretched and stained DNA molecules are produced. These “fingerprints” of the underlying DNA sequence have a resolution of the order one kilobasepairs and the stretching of the DNA molecules are performed by surface adsorption or nano-channel setups. A post-processing challenge for nano-channel based methods, due to local and global random movement of the DNA molecule during imaging, is how to align different time frames in order to produce reproducible time-averaged DNA barcodes. The current solutions to this challenge are computationally rather slow. With high-throughput applications in mind, we here introduce a parameter-free method for filtering a single time frame noisy barcode (snap-shot optical map), measured in a fraction of a second. By using only a single time frame barcode we circumvent the need for post-processing alignment. We demonstrate that our method is successful at providing filtered barcodes which are less noisy and more similar to time averaged barcodes. The method is based on the application of a low-pass filter on a single noisy barcode using the width of the Point Spread Function of the system as a unique, and known, filtering parameter. We find that after applying our method, the Pearson correlation coefficient (a real number in the range from -1 to 1) between the single time-frame barcode and the time average of the aligned kymograph increases significantly, roughly by 0.2 on average. By comparing to a database of more than 3000 theoretical plasmid barcodes we show that the capabilities to identify plasmids is improved by filtering single time-frame barcodes compared to the unfiltered analogues. Since snap-shot experiments and computational time using our method both are less than a second, this study opens up for high throughput optical DNA mapping with improved reproducibility.

The impact of third generation genomic technologies on plant genome assembly

Since the introduction of next generation sequencing, plant genome assembly projects do not need to rely on dedicated research facilities or community-wide consortia anymore, even individual research groups can sequence and assemble the genomes they are interested in. However, such assemblies are typically not based on the entire breadth of genomic technologies including genetic and physical maps and their contiguities tend to be low compared to the full-length gold standard reference sequences. Recently emerging third generation genomic technologies like long-read sequencing or optical mapping promise to bridge this quality gap and enable simple and cost-effective solutions for chromosomal-level assemblies.

Dioecy in Amborella trichopoda: evidence for genetically based sex determination and its consequences for inferences of the breeding system in early angiosperms

BACKGROUND AND AIMS

This work aimed to gain insight into the breeding system at the base of living angiosperms through both character state reconstructions and the study of sex ratios and phenotypes in the likely sister to all other living angiosperms, Amborella trichopoda .

METHODS

Sex phenotypes were mapped onto a phylogeny of basally diverging angiosperms using maximum parsimony. In parallel, sex ratios and phenotypes were studied over two consecutive flowering seasons in an ex situ population of A. trichopoda , while the sex ratio of an in situ population was also assessed.

KEY RESULTS

Parsimony analyses failed to resolve the breeding system present at the base of living angiosperms, but indicated the importance of A. trichopoda for the future elucidation of this question. The ex situ A. trichopoda population studied showed a primary sex ratio close to 1:1, though sex ratio bias was found in the in situ population studied. Instances of sexual instability were quantified in both populations.

CONCLUSIONS

Sex ratio data support the presence of genetic sex determination in A. trichopoda , whose further elucidation may guide inferences on the breeding system at the base of living angiosperms. Sexual instability in A. trichopoda suggests the operation of epigenetic mechanisms, and the evolution of dioecy via a gynodioecious intermediate.

Ensembl Plants: Integrating Tools for Visualizing, Mining, and Analyzing Plant Genomic Data

Ensembl Plants ( http://plants.ensembl.org ) is an integrative resource presenting genome-scale information for 39 sequenced plant species. Available data includes genome sequence, gene models, functional annotation, and polymorphic loci; for the latter, additional information including population structure, individual genotypes, linkage, and phenotype data is available for some species. Comparative data is also available, including genomic alignments and "gene trees," which show the inferred evolutionary history of each gene family represented in the resource. Access to the data is provided through a genome browser, which incorporates many specialist interfaces for different data types, through a variety of programmatic interfaces, and via a specialist data mining tool supporting rapid filtering and retrieval of bulk data. Genomic data from many non-plant species, including those of plant pathogens, pests, and pollinators, is also available via the same interfaces through other divisions of Ensembl.Ensembl Plants is updated 4-6 times a year and is developed in collaboration with our international partners in the Gramene ( http://www.gramene.org ) and transPLANT projects ( http://www.transplantdb.eu ).

Comparative phylogenomics of symbiotic associations

89 I. 89 II. 90 III. 90 IV. 91 V. 92 VI. 93 References 93 SUMMARY: Understanding the genetic bases of complex traits has been a main challenge in biology for decades. Comparative phylogenomics offers an opportunity to identify candidate genes associated with these complex traits. This approach initially developed in prokaryotes consists in looking at shared coevolution between genes and traits. It thus requires a precise reconstruction of the trait evolution, a large genomic sampling in the clades of interest and an accurate definition of orthogroups. Recently, with the growing body of sequenced plant genomes, comparative genomics has been successfully applied to plants to study the widespread arbuscular mycorrhizal symbiosis. Here I will use these findings to illustrate the main principles of comparative phylogenomic approaches and propose directions to improve our understanding of symbiotic associations.

Species-Specific Genome Sequence Databases: A Practical Review

This chapter presents a use case illustrating the search for homologues of a known protein in species-specific genome sequence databases. The results from different species-specific resources are compared to each other and to results obtained from a more general genome sequence database (Phytozome). Various options and settings relevant when searching these databases are discussed. For example, it is shown how the choice of reference sequence set in a given database influences the results one obtains. The provided examples illustrate some problems and pitfalls related to interpreting results obtained from species-specific genome sequence databases.

Genome Mapping in Plant Comparative Genomics

Genome mapping produces fingerprints of DNA sequences to construct a physical map of the whole genome. It provides contiguous, long-range information that complements and, in some cases, replaces sequencing data. Recent advances in genome-mapping technology will better allow researchers to detect large (>1kbp) structural variations between plant genomes. Some molecular and informatics complications need to be overcome for this novel technology to achieve its full utility. This technology will be useful for understanding phenotype responses due to DNA rearrangements and will yield insights into genome evolution, particularly in polyploids. In this review, we outline recent advances in genome-mapping technology, including the processes required for data collection and analysis, and applications in plant comparative genomics.

Identification of peanut (Arachis hypogaea) chromosomes using a fluorescence in situ hybridization system reveals multiple hybridization events during tetraploid peanut formation

The cultivated peanut Arachis hypogaea (AABB) is thought to have originated from the hybridization of Arachis duranensis (AA) and Arachis ipaënsis (BB) followed by spontaneous chromosome doubling. In this study, we cloned and analyzed chromosome markers from cultivated peanut and its wild relatives. A fluorescence in situ hybridization (FISH)-based karyotyping cocktail was developed with which to study the karyotypes and chromosome evolution of peanut and its wild relatives. Karyotypes were constructed in cultivated peanut and its two putative progenitors using our FISH-based karyotyping system. Comparative karyotyping analysis revealed that chromosome organization was highly conserved in cultivated peanut and its two putative progenitors, especially in the B genome chromosomes. However, variations existed between A. duranensis and the A genome chromosomes in cultivated peanut, especially for the distribution of the interstitial telomere repeats (ITRs). A search of additional A. duranensis varieties from different geographic regions revealed both numeric and positional variations of ITRs, which were similar to the variations in tetraploid peanut varieties. The results provide evidence for the origin of cultivated peanut from the two diploid ancestors, and also suggest that multiple hybridization events of A. ipaënsis with different varieties of A. duranensis may have occurred during the origination of peanut.

The phylogeny of C/S1 bZIP transcription factors reveals a shared algal ancestry and the pre-angiosperm translational regulation of S1 transcripts

Basic leucine zippers (bZIPs) form a large plant transcription factor family. C and S1 bZIP groups can heterodimerize, fulfilling crucial roles in seed development and stress response. S1 sequences also harbor a unique regulatory mechanism, termed Sucrose-Induced Repression of Translation (SIRT). The conservation of both C/S1 bZIP interactions and SIRT remains poorly characterized in non-model species, leaving their evolutionary origin uncertain and limiting crop research. In this work, we explored recently published plant sequencing data to establish a detailed phylogeny of C and S1 bZIPs, investigating their intertwined role in plant evolution, and the origin of SIRT. Our analyses clarified C and S1 bZIP orthology relationships in angiosperms, and identified S1 sequences in gymnosperms. We experimentally showed that the gymnosperm orthologs are regulated by SIRT, tracing back the origin of this unique regulatory mechanism to the ancestor of seed plants. Additionally, we discovered an earlier S ortholog in the charophyte algae Klebsormidium flaccidum, together with a C ortholog. This suggests that C and S groups originated by duplication from a single algal proto-C/S ancestor. Based on our observations, we propose a model wherein the C/S1 bZIP dimer network evolved in seed plants from pre-existing C/S bZIP interactions.

Application of a modified drop method for high-resolution pachytene chromosome spreads in two Phalaenopsis species

BACKGROUND

Preparation of good chromosome spreads without cytoplasmic contamination is the crucial step in cytogenetic mapping. To date, cytogenetic research in the Orchidaceae family has been carried out solely on mitotic metaphase chromosomes. Well-spread meiotic pachytene chromosomes can provide higher resolution and fine detail for analysis of chromosomal structure and are also beneficial for chromosomal FISH (fluorescence in situ hybridization) mapping. However, an adequate method for the preparation of meiotic pachytene chromosomes in orchid species has not yet been reported.

RESULTS

Two Taiwanese native Phalaenopsis species were selected to test the modified drop method for preparation of meiotic pachytene chromosomes from pollinia. In this modified method, pollinia were ground and treated with an enzyme mixture to completely remove cell walls. Protoplasts were resuspended in ethanol/glacial acetic acid and dropped onto a wet inclined slide of 30° from a height of 0.5 m. The sample was then flowed down the inclined plane to spread the chromosomes. Hundreds of pachytene chromosomes with little to no cytoplasmic contamination were well spread on each slide. We also showed that the resolution of 45S rDNA-containing chromosomes at the pachytene stage was up to 20 times higher than that at metaphase. Slides prepared following this modified drop method were amenable to FISH mapping of both 45S and 5S rDNA on pachytene chromosomes and, after FISH, the chromosomal structure remained intact for further analysis.

CONCLUSION

This modified drop method is suitable for pachytene spreads from pollinia of Phalaenopsis orchids. The large number and high-resolution pachytene spreads, with little or no cytoplasmic contamination, prepared by the modified drop method could be used for FISH mapping of DNA fragments to accelerate the integration of cytogenetic and molecular research in Phalaenopsis orchids.

Advancing Eucalyptus Genomics: Cytogenomics Reveals Conservation of Eucalyptus Genomes

The genus Eucalyptus encloses several species with high ecological and economic value, being the subgenus Symphyomyrtus one of the most important. Species such as E. grandis and E. globulus are well characterized at the molecular level but knowledge regarding genome and chromosome organization is very scarce. Here we characterized and compared the karyotypes of three economically important species, E. grandis, E. globulus, and E. calmadulensis, and three with ecological relevance, E. pulverulenta, E. cornuta, and E. occidentalis, through an integrative approach including genome size estimation, fluorochrome banding, rDNA FISH, and BAC landing comprising genes involved in lignin biosynthesis. All karyotypes show a high degree of conservation with pericentromeric 35S and 5S rDNA loci in the first and third pairs, respectively. GC-rich heterochromatin was restricted to the 35S rDNA locus while the AT-rich heterochromatin pattern was species-specific. The slight differences in karyotype formulas and distribution of AT-rich heterochromatin, along with genome sizes estimations, support the idea of Eucalyptus genome evolution by local expansions of heterochromatin clusters. The unusual co-localization of both rDNA with AT-rich heterochromatin was attributed mainly to the presence of silent transposable elements in those loci. The cinnamoyl CoA reductase gene (CCR1) previously assessed to linkage group 10 (LG10) was clearly localized distally at the long arm of chromosome 9 establishing an unexpected correlation between the cytogenetic chromosome 9 and the LG10. Our work is novel and contributes to the understanding of Eucalyptus genome organization which is essential to develop successful advanced breeding strategies for this genus.

MaGuS: a tool for quality assessment and scaffolding of genome assemblies with Whole Genome Profiling™ Data

Background

Scaffolding is an essential step in the genome assembly process. Current methods based on large fragment paired-end reads or long reads allow an increase in contiguity but often lack consistency in repetitive regions, resulting in fragmented assemblies. Here, we describe a novel tool to link assemblies to a genome map to aid complex genome reconstruction by detecting assembly errors and allowing scaffold ordering and anchoring.

Results

We present MaGuS (map-guided scaffolding), a modular tool that uses a draft genome assembly, a Whole Genome Profiling™ (WGP) map, and high-throughput paired-end sequencing data to estimate the quality and to enhance the contiguity of an assembly. We generated several assemblies of the Arabidopsis genome using different scaffolding programs and applied MaGuS to select the best assembly using quality metrics. Then, we used MaGuS to perform map-guided scaffolding to increase contiguity by creating new scaffold links in low-covered and highly repetitive regions where other commonly used scaffolding methods lack consistency.

Conclusions

MaGuS is a powerful reference-free evaluator of assembly quality and a WGP map-guided scaffolder that is freely available at https://github.com/institut-de-genomique/MaGuS. Its use can be extended to other high-throughput sequencing data (e.g., long-read data) and also to other map data (e.g., genetic maps) to improve the quality and the contiguity of large and complex genome assemblies.

Electronic supplementary material

The online version of this article (doi:10.1186/s12859-016-0969-x) contains supplementary material, which is available to authorized users.

The role of ABC genes in shaping perianth phenotype in the basal angiosperm Magnolia

It is generally accepted that the genus Magnolia is characterised by an undifferentiated perianth, typically organised into three whorls of nearly identical tepals. In some species, however, we encountered interesting and significant perianth modifications. In Magnolia acuminata, M. liliiflora and M. stellata the perianth elements of the first whorl are visually different from the others. In M. stellata the additional, spirally arranged perianth elements are present above the first three whorls, which suggests that they have been formed within the domain of stamen primordia. In these three species, we analysed expression patterns of the key flower genes (AP1, AGL6, AP3, PI, AG) responsible for the identity of flower elements and correlated them with results of morphological and anatomical investigations. In all studied species the elements of the first whorl lacked the identity of petals (lack of AP3 and PI expression) but also that of leaves (presence of AGL6 expression), and this seems to prove their sepal character. The analysis of additional perianth elements of M. stellata, spirally arranged on the elongated floral axis, revealed overlapping and reduced activity of genes involved in specification of the identity of the perianth (AGL6) but also of generative parts (AG), even though no clear gradient of morphological changes could be observed. In conclusion, Magnolia genus is capable of forming, in some species, a perianth differentiated into a calyx (sepals) and corolla (petals). Spirally arranged, additional perianth elements of M. stellata, despite activity of AG falling basipetally, resemble petals.

The genome sequence of the outbreeding globe artichoke constructed de novo incorporating a phase-aware low-pass sequencing strategy of F1 progeny

Globe artichoke (Cynara cardunculus var. scolymus) is an out-crossing, perennial, multi-use crop species that is grown worldwide and belongs to the Compositae, one of the most successful Angiosperm families. We describe the first genome sequence of globe artichoke. The assembly, comprising of 13,588 scaffolds covering 725 of the 1,084 Mb genome, was generated using ~133-fold Illumina sequencing data and encodes 26,889 predicted genes. Re-sequencing (30×) of globe artichoke and cultivated cardoon (C. cardunculus var. altilis) parental genotypes and low-coverage (0.5 to 1×) genotyping-by-sequencing of 163 F₁ individuals resulted in 73% of the assembled genome being anchored in 2,178 genetic bins ordered along 17 chromosomal pseudomolecules. This was achieved using a novel pipeline, SOILoCo (Scaffold Ordering by Imputation with Low Coverage), to detect heterozygous regions and assign parental haplotypes with low sequencing read depth and of unknown phase. SOILoCo provides a powerful tool for de novo genome analysis of outcrossing species. Our data will enable genome-scale analyses of evolutionary processes among crops, weeds, and wild species within and beyond the Compositae, and will facilitate the identification of economically important genes from related species.

Gramene: A Resource for Comparative Analysis of Plants Genomes and Pathways

Gramene is an integrated informatics resource for accessing, visualizing, and comparing plant genomes and biological pathways. Originally targeting grasses, Gramene has grown to host annotations for economically important and research model crops, including wheat, potato, tomato, banana, grape, poplar, and Chlamydomonas. Its strength derives from the application of a phylogenetic framework for genome comparison and the use of ontologies to integrate structural and functional annotation data. This chapter outlines system requirements for end users and database hosting, data types and basic navigation within Gramene, and provides examples of how to (1) view a phylogenetic tree for a family of transcription factors, (2) explore genetic variation in the orthologues of a gene with a known trait association, and (3) upload, visualize, and privately share end user data into a new genome browser track.Moreover, this is the first publication describing Gramene's new web interface-intended to provide a simplified portal to the most complete and up-to-date set of plant genome and pathway annotations.

Ensembl Plants: Integrating Tools for Visualizing, Mining, and Analyzing Plant Genomics Data

Ensembl Plants ( http://plants.ensembl.org ) is an integrative resource presenting genome-scale information for a growing number of sequenced plant species (currently 33). Data provided includes genome sequence, gene models, functional annotation, and polymorphic loci. Various additional information are provided for variation data, including population structure, individual genotypes, linkage, and phenotype data. In each release, comparative analyses are performed on whole genome and protein sequences, and genome alignments and gene trees are made available that show the implied evolutionary history of each gene family. Access to the data is provided through a genome browser incorporating many specialist interfaces for different data types, and through a variety of additional methods for programmatic access and data mining. These access routes are consistent with those offered through the Ensembl interface for the genomes of non-plant species, including those of plant pathogens, pests, and pollinators.Ensembl Plants is updated 4-5 times a year and is developed in collaboration with our international partners in the Gramene ( http://www.gramene.org ) and transPLANT projects ( http://www.transplantdb.org ).

The map-based genome sequence of Spirodela polyrhiza aligned with its chromosomes, a reference for karyotype evolution

Duckweeds are aquatic monocotyledonous plants of potential economic interest with fast vegetative propagation, comprising 37 species with variable genome sizes (0.158-1.88 Gbp). The genomic sequence of Spirodela polyrhiza, the smallest and the most ancient duckweed genome, needs to be aligned to its chromosomes as a reference and prerequisite to study the genome and karyotype evolution of other duckweed species. We selected physically mapped bacterial artificial chromosomes (BACs) containing Spirodela DNA inserts with little or no repetitive elements as probes for multicolor fluorescence in situ hybridization (mcFISH), using an optimized BAC pooling strategy, to validate its physical map and correlate it with its chromosome complement. By consecutive mcFISH analyses, we assigned the originally assembled 32 pseudomolecules (supercontigs) of the genomic sequences to the 20 chromosomes of S. polyrhiza. A Spirodela cytogenetic map containing 96 BAC markers with an average distance of 0.89 Mbp was constructed. Using a cocktail of 41 BACs in three colors, all chromosome pairs could be individualized simultaneously. Seven ancestral blocks emerged from duplicated chromosome segments of 19 Spirodela chromosomes. The chromosomally integrated genome of S. polyrhiza and the established prerequisites for comparative chromosome painting enable future studies on the chromosome homoeology and karyotype evolution of duckweed species.

Misassembly detection using paired-end sequence reads and optical mapping data

Motivation: A crucial problem in genome assembly is the discovery and correction of misassembly errors in draft genomes. We develop a method called misSEQuel that enhances the quality of draft genomes by identifying misassembly errors and their breakpoints using paired-end sequence reads and optical mapping data. Our method also fulfills the critical need for open source computational methods for analyzing optical mapping data. We apply our method to various assemblies of the loblolly pine, Francisella tularensis, rice and budgerigar genomes. We generated and used stimulated optical mapping data for loblolly pine and F.tularensis and used real optical mapping data for rice and budgerigar.

Results: Our results demonstrate that we detect more than 54% of extensively misassembled contigs and more than 60% of locally misassembled contigs in assemblies of F.tularensis and between 31% and 100% of extensively misassembled contigs and between 57% and 73% of locally misassembled contigs in assemblies of loblolly pine. Using the real optical mapping data, we correctly identified 75% of extensively misassembled contigs and 100% of locally misassembled contigs in rice, and 77% of extensively misassembled contigs and 80% of locally misassembled contigs in budgerigar.

Availability and implementation:misSEQuel can be used as a post-processing step in combination with any genome assembler and is freely available at http://www.cs.colostate.edu/seq/.

Contact:muggli@cs.colostate.edu

Supplementary information:Supplementary data are available at Bioinformatics online.

Gramene 2016: comparative plant genomics and pathway resources

Gramene (http://www.gramene.org) is an online resource for comparative functional genomics in crops and model plant species. Its two main frameworks are genomes (collaboration with Ensembl Plants) and pathways (The Plant Reactome and archival BioCyc databases). Since our last NAR update, the database website adopted a new Drupal management platform. The genomes section features 39 fully assembled reference genomes that are integrated using ontology-based annotation and comparative analyses, and accessed through both visual and programmatic interfaces. Additional community data, such as genetic variation, expression and methylation, are also mapped for a subset of genomes. The Plant Reactome pathway portal (http://plantreactome.gramene.org) provides a reference resource for analyzing plant metabolic and regulatory pathways. In addition to ∼ 200 curated rice reference pathways, the portal hosts gene homology-based pathway projections for 33 plant species. Both the genome and pathway browsers interface with the EMBL-EBI's Expression Atlas to enable the projection of baseline and differential expression data from curated expression studies in plants. Gramene's archive website (http://archive.gramene.org) continues to provide previously reported resources on comparative maps, markers and QTL. To further aid our users, we have also introduced a live monthly educational webinar series and a Gramene YouTube channel carrying video tutorials.