BioNano genome mapping of individual chromosomes supports physical mapping and sequence assembly in complex plant genomes

DNA Quantity and Quality Comparisons between Cryopreserved and FFPE Tumors from Matched Pan-Cancer Samples

Personalized cancer care requires molecular characterization of neoplasms. While the research community accepts frozen tissues as the gold standard analyte for molecular assays, the source of tissue for testing in clinical cancer care comes almost universally from formalin-fixed, paraffin-embedded tissue (FFPE). As newer technologies emerge for DNA characterization that requires higher molecular weight DNA, it was necessary to compare the quality of DNA in terms of DNA length between FFPE and cryopreserved samples. We hypothesized that cryopreserved samples would yield higher quantity and superior quality DNA compared to FFPE samples. We analyzed DNA metrics by performing a head-to-head comparison between FFPE and cryopreserved samples from 38 human tumors representing various cancer types. DNA quantity and purity were measured by UV spectrophotometry, and DNA from cryopreserved tissue demonstrated a 4.2-fold increase in DNA yield per mg of tissue (p-value < 0.001). DNA quality was measured on a fragment microelectrophoresis analyzer, and again, DNA from cryopreserved tissue demonstrated a 223% increase in the DNA quality number and a 9-fold increase in DNA fragments > 40,000 bp (p-value < 0.0001). DNA from the cryopreserved tissues was superior to the DNA from FFPE samples in terms of DNA yield and quality.

Structural Refinement by Direct Mapping Reveals Assembly Inconsistencies near Hi-C Junctions

High-throughput chromosome conformation capture (Hi-C) is widely used for scaffolding in de novo assembly because it produces highly contiguous genomes, but its indirect statistical approach can introduce connection errors. We employed optical mapping (Bionano Genomics) as an orthogonal scaffolding technology to assess the structural solidity of Hi-C reconstructed scaffolds. Optical maps were used to assess the correctness of five de novo genome assemblies based on long-read sequencing for contig generation and Hi-C for scaffolding. Hundreds of inconsistencies were found between the reconstructions generated using the Hi-C and optical mapping approaches. Manual inspection, exploiting raw long-read sequencing data and optical maps, confirmed that several of these conflicts were derived from Hi-C joining errors. Such misjoins were widespread, involved the connection of both small and large contigs, and even overlapped annotated genes. We conclude that the integration of optical mapping data after, not before, Hi-C-based scaffolding, improves the quality of the assembly and limits reconstruction errors by highlighting misjoins that can then be subjected to further investigation.

Flow Cytometric Analysis and Sorting of Plant Chromosomes

Flow cytometry offers a unique way of analyzing and manipulating plant chromosomes. During a rapid movement in a liquid stream, large populations can be classified in a short time according to their fluorescence and light scatter properties. Chromosomes whose optical properties differ from other chromosomes in a karyotype can be purified by flow sorting and used in a range of applications in cytogenetics, molecular biology, genomics, and proteomics. As the samples for flow cytometry must be liquid suspensions of single particles, intact chromosomes must be released from mitotic cells. This protocol describes a procedure for preparation of suspensions of mitotic metaphase chromosomes from meristem root tips and their flow cytometric analysis and sorting for various downstream applications.

Flow Sorting-Assisted Optical Mapping

Optical mapping-a technique that visualizes short sequence motives along DNA molecules of hundred kilobases to megabase in size-has found an important place in genome research. It is widely used to facilitate genome sequence assemblies and analyses of genome structural variations. Application of the technique is conditional on availability of highly pure ultra-long high-molecular-weight DNA (uHMW DNA), which is challenging to achieve in plants due to the presence of the cell wall, chloroplasts, and secondary metabolites, just as a high content of polysaccharides and DNA nucleases in some species. These obstacles can be overcome by employment of flow cytometry, enabling a fast and highly efficient purification of cell nuclei or metaphase chromosomes, which are afterward embedded in agarose plugs and used to isolate the uHMW DNA in situ. Here, we provide a detailed protocol for the flow sorting-assisted uHMW DNA preparation that has been successfully used to construct whole-genome as well as chromosomal optical maps for 20 plant species from several plant families.

Flow karyotyping of wheat-Aegilops additions facilitate dissecting the genomes of Ae. biuncialis and Ae. geniculata into individual chromosomes

Breeding of wheat adapted to new climatic conditions and resistant to diseases and pests is hindered by a limited gene pool due to domestication and thousands of years of human selection. Annual goatgrasses (Aegilops spp.) with M and U genomes are potential sources of the missing genes and alleles. Development of alien introgression lines of wheat may be facilitated by the knowledge of DNA sequences of Aegilops chromosomes. As the Aegilops genomes are complex, sequencing relevant Aegilops chromosomes purified by flow cytometric sorting offers an attractive route forward. The present study extends the potential of chromosome genomics to allotetraploid Ae. biuncialis and Ae. geniculata by dissecting their M and U genomes into individual chromosomes. Hybridization of FITC-conjugated GAA oligonucleotide probe to chromosomes suspensions of the two species allowed the application of bivariate flow karyotyping and sorting some individual chromosomes. Bivariate flow karyotype FITC vs. DAPI of Ae. biuncialis consisted of nine chromosome-populations, but their chromosome content determined by microscopic analysis of flow sorted chromosomes indicated that only 7M^b and 1U^b could be sorted at high purity. In the case of Ae. geniculata, fourteen chromosome-populations were discriminated, allowing the separation of nine individual chromosomes (1M^g, 3M^g, 5M^g, 6M^g, 7M^g, 1U^g, 3U^g, 6U^g, and 7U^g) out of the 14. To sort the remaining chromosomes, a partial set of wheat-Ae. biuncialis and a whole set of wheat-Ae. geniculata chromosome addition lines were also flow karyotyped, revealing clear separation of the GAA-rich Aegilops chromosomes from the GAA-poor A- and D-genome chromosomes of wheat. All of the alien chromosomes represented by individual addition lines could be isolated at purities ranging from 74.5% to 96.6% and from 87.8% to 97.7%, respectively. Differences in flow karyotypes between Ae. biuncialis and Ae. geniculata were analyzed and discussed. Chromosome-specific genomic resources will facilitate gene cloning and the development of molecular tools to support alien introgression breeding of wheat.

Reference quality genome sequence of Indian pomegranate cv. 'Bhagawa' (Punica granatum L.)

Pomegranate is an important fruit crop for ensuring livelihood and nutrition security in fragile semi-arid regions of the globe having limited irrigation resources. This is a high-value, nutritionally rich, and export-oriented agri-commodity that ensures high returns on investment to growers across the world. Although it is a valuable fruit crop, it has received only a limited genomics research outcome. To fast-track the pomegranate improvement program, de novo whole-genome sequencing of the main Indian cultivar 'Bhagawa' was initiated by the Indian Council of Agricultural Research-National Research Center on Pomegranate (ICAR-NRCP). We have demonstrated that a combination of commercially available technologies from Illumina, PacBio, 10X Genomics, and BioNano Genomics could be used efficiently for sequencing and reference-grade de novo assembly of the pomegranate genome. The research led to a final reference-quality genome assembly for 'Bhagawa' of 346.08 Mb in 342 scaffolds and an average N50 of 16.12 Mb and N90 of 1088.62 Kb. This assembly covered more than 98% of the estimated pomegranate genome size, 352.54 Mb. The LTR assembly index (LAI) value of 10 and 93.68% Benchmarking Universal Single-Copy Orthologs (BUSCO) completeness score over the 1,440 ortholog genes of the completed pomegranate genome indicates the quality of the assembled pomegranate genome. Furthermore, 29,435 gene models were discovered with a mean transcript length of 2,954 bp and a mean coding sequence length 1,090 bp. Four transcript data samples of pomegranate tissues were mapped over the assembled 'Bhagawa' genome up to 95% significant matches, indicating the high quality of the assembled genome. We have compared the 'Bhagawa' genome with the genomes of the pomegranate cultivars 'Dabenzi' and 'Taishanhong.' We have also performed whole-genome phylogenetic analysis using Computational Analysis of Gene Family Evolution (CAFE) and found that Eucalyptus grandis and pomegranate diverged 64 (60-70) million years ago. About 1,573 protein-coding resistance genes identified in the 'Bhagawa' genome were classified into 32 domains. In all, 314 copies of miRNA belonging to 26 different families were identified in the 'Bhagawa' genome. The reference-quality genome assembly of 'Bhagawa' is certainly a significant genomic resource for accelerated pomegranate improvement.

Meiotic Chromosomal Abnormality Detected in a Heterozygote of Elymus nutans

Elymus nutans is an allopolyploid with a genome constitution of StStYYHH (2n = 6x = 42). Highly frequent intergenomic translocations and chromosomal variations with repeat amplification and deletions in E. nutans have been identified in the previous studies. However, more complicated structural variations such as chromosomal inversions or intra-genomic translocations are still unknown in this species, so does the reason for the origin of the chromosomal variations. Heterozygotes with rearranged chromosomes always present irregular meiosis behaviors, which subsequently cause the secondary chromosome rearrangements. Investigation on the meiosis of heterozygotes, especially on the individual chromosome level, may provide the important clues to identify the more complicated chromosome structural variations in the populations, and clarify the origin of the chromosome variations. In this study, meiotic analysis was conducted on a heterozygote plant of Elymus nutans, which showed high intra- and inter-genome chromosomal variations, by sequential fluorescence in situ hybridization (FISH) and genomic in situ hybridization (GISH), with each chromosome clearly recognized. The results showed chromosomal abnormalities at every meiotic stage and abnormalities in frequency variations between different sub-genomes and different individual chromosomes. The abnormalities were revealed as univalent, fragment, rod, or Y shape bivalent in diakinesis; univalent and rod bivalent in metaphase I; lagged and segregated chromatid, bridge, fragment of the sister chromatid, fragment, bridge accompanied with fragment, and unequal segregated chromosome in anaphase I; bridge and lagged chromatid in ana-telophase II; and micronucleus at uninucleate stage. Generally, the St and H genomes harbor more abnormalities than the Y genome. Moreover, a paracentric inversion in 2St was exclusively determined, and another paracentric inversion in 6Y was tentatively identified. In addition, novel deletions were clearly detected in 3H, 4H, 1Y, and 3Y homologous chromosomes; in particular, de novo pericentric inversion in 3H was repeatedly identified in metaphase I. The study revealed the chromosomal inversions pre-existed in parents or populations, as well as de novo inversions and deletions originated in the meiosis of the heterozygote in E. nutans. Moreover, it indicated wide range of meiosis abnormalities on different stages and different chromosomes, and suggests that secondary rearrangements contribute much to the chromosome variations in E. nutans.

Fine structure and transcription dynamics of bread wheat ribosomal DNA loci deciphered by a multi-omics approach

Three out of four RNA components of ribosomes are encoded by 45S ribosomal DNA (rDNA) loci, which are organized as long head-to-tail tandem arrays of nearly identical units, spanning several megabases of sequence. Due to this structure, the rDNA loci are the major sources of gaps in genome assemblies, and gene copy number, sequence composition, and expression status of particular arrays remain elusive, especially in complex genomes harboring multiple loci. Here we conducted a multi-omics study to decipher the 45S rDNA loci in hexaploid bread wheat. Coupling chromosomal genomics with optical mapping, we reconstructed individual rDNA arrays, enabling locus-specific analyses of transcription activity and methylation status from RNA- and bisulfite-sequencing data. We estimated a total of 6,650 rDNA units in the bread wheat genome, with approximately 2,321, 3,910, 253, and 50 gene copies located in short arms of chromosomes 1B, 6B, 5D, and 1A, respectively. Only 1B and 6B loci contributed substantially to rRNA transcription at a roughly 2:1 ratio. The ratio varied among five tissues analyzed (embryo, coleoptile, root tip, primary leaf, mature leaf), being the highest (2.64:1) in mature leaf and lowest (1.72:1) in coleoptile. Cytosine methylation was considerably higher in CHG context in the silenced 5D locus as compared with the active 1B and 6B loci. In conclusion, a fine genomic organization and tissue-specific expression of rDNA loci were deciphered, for the first time, in a complex polyploid species. The results are discussed in the context of wheat evolution and transcription regulation.

Optical mapping of the Fusarium oxysporum f. sp. melongenae genome

Optical mapping approaches are widely preferred and applied in different branches of genomic studies because of their accuracy, low cost, and high efficiency. In the current study, a sequence orientation of the Fusarium oxysporum f. sp. melongenae (FOMG) genome that is deposited in GenBank NCBI under accession number MPIL00000000 was used as the reference genome, which we checked with Bionano Genomics (BNG) optical mapping approaches. Optical mapping produced 103 contigs, the longest of which was 3.05 Mb. The N50 value of optical map contigs is 0.85 Mb. The sequences of the FOMG reference genome and optical map mainly match each other. Results obtained in the current study indicate that optical mapping can be used to construct complete and gapless assemblies of the FOMG genome. It also can be applied to validate a previous genome assembly. This article is protected by copyright. All rights reserved.

Two gap-free reference genomes and a global view of the centromere architecture in rice

Rice (Oryza sativa), a major staple throughout the world and a model system for plant genomics and breeding, was the first crop genome sequenced almost two decades ago. However, reference genomes for all higher organisms to date contain gaps and missing sequences. Here, we report the assembly and analysis of gap-free reference genome sequences for two elite O. sativa xian/indica rice varieties, Zhenshan 97 and Minghui 63, which are being used as a model system for studying heterosis and yield. Gap-free reference genomes provide the opportunity for a global view of the structure and function of centromeres. We show that all rice centromeric regions share conserved centromere-specific satellite motifs with different copy numbers and structures. In addition, the similarity of CentO repeats in the same chromosome is higher than across chromosomes, supporting a model of local expansion and homogenization. Both genomes have over 395 non-TE genes located in centromere regions, of which ∼41% are actively transcribed. Two large structural variants at the end of chromosome 11 affect the copy number of resistance genes between the two genomes. The availability of the two gap-free genomes lays a solid foundation for further understanding genome structure and function in plants and breeding climate-resilient varieties.

Signal-based optical map alignment

In genomics, optical mapping technology provides long-range contiguity information to improve genome sequence assemblies and detect structural variation. Originally a laborious manual process, Bionano Genomics platforms now offer high-throughput, automated optical mapping based on chips packed with nanochannels through which unwound DNA is guided and the fluorescent DNA backbone and specific restriction sites are recorded. Although the raw image data obtained is of high quality, the processing and assembly software accompanying the platforms is closed source and does not seem to make full use of data, labeling approximately half of the measured signals as unusable. Here we introduce two new software tools, independent of Bionano Genomics software, to extract and process molecules from raw images (OptiScan) and to perform molecule-to-molecule and molecule-to-reference alignments using a novel signal-based approach (OptiMap). We demonstrate that the molecules detected by OptiScan can yield better assemblies, and that the approach taken by OptiMap results in higher use of molecules from the raw data. These tools lay the foundation for a suite of open-source methods to process and analyze high-throughput optical mapping data. The Python implementations of the OptiTools are publicly available through http://www.bif.wur.nl/.

Genome reconstruction and haplotype phasing using chromosome conformation capture methodologies

Genomic analysis of individuals or organisms is predicated on the availability of high-quality reference and genotype information. With the rapidly dropping costs of high-throughput DNA sequencing, this is becoming readily available for diverse organisms and for increasingly large populations of individuals. Despite these advances, there are still aspects of genome sequencing that remain challenging for existing sequencing methods. This includes the generation of long-range contiguity during genome assembly, identification of structural variants in both germline and somatic tissues, the phasing of haplotypes in diploid organisms and the resolution of genome sequence for organisms derived from complex samples. These types of information are valuable for understanding the role of genome sequence and genetic variation on genome function, and numerous approaches have been developed to address them. Recently, chromosome conformation capture (3C) experiments, such as the Hi-C assay, have emerged as powerful tools to aid in these challenges for genome reconstruction. We will review the current use of Hi-C as a tool for aiding in genome sequencing, addressing the applications, strengths, limitations and potential future directions for the use of 3C data in genome analysis. We argue that unique features of Hi-C experiments make this data type a powerful tool to address challenges in genome sequencing, and that future integration of Hi-C data with alternative sequencing assays will facilitate the continuing revolution in genomic analysis and genome sequencing.

Comparison of structural variants detected by optical mapping with long-read next-generation sequencing

MOTIVATION

Recent studies have shown the potential of using long-read whole-genome sequencing (WGS) approaches and optical mapping (OM) for the detection of clinically relevant structural variants (SVs) in cancer research. Three main long-read WGS platforms are currently in use: Pacific Biosciences (PacBio), Oxford Nanopore Technologies (ONT) and 10x Genomics. Recently, whole-genome OM technology (Bionano Genomics) has been introduced into human diagnostics. Questions remain about the accuracy of these long-read sequencing platforms, how comparable/interchangeable they are when searching for SVs and to what extent they can be replaced or supplemented by OM. Moreover, no tool can effectively compare SVs obtained by OM and WGS.

RESULTS

This study compared optical maps of the breast cancer cell line SKBR3 with AnnotSV outputs from WGS platforms. For this purpose, a software tool with comparative and filtering features was developed. The majority of SVs up to a 50 kbp distance variance threshold found by OM were confirmed by all WGS platforms, and 99% of translocations and 80% of deletions found by OM were confirmed by both PacBio and ONT, with ∼70% being confirmed by 10x Genomics in combination with PacBio and/or ONT. Interestingly, long deletions (>100 kbp) were detected only by 10x Genomics. Regarding insertions, ∼72% was confirmed by PacBio and ONT, but none by 10x Genomics. Inversions and duplications detected by OM were not detected by WGS. Moreover, the tool enabled the confirmation of SVs that overlapped in the same gene(s) and was applied to the filtering of disease-associated SVs.

AVAILABILITY

https://github.com/novosadt/om-annotsv-svc.

Genome Sequencing and Assembly Strategies and a Comparative Analysis of the Genomic Characteristics in Penaeid Shrimp Species

Penaeid shrimp (family Penaeidae) represents one of the most economically and ecologically important groups of crustaceans. However, their genome sequencing and assembly have encountered extreme difficulties during the last 20 years. In this study, based on our previous genomic data, we investigated the genomic characteristics of four penaeid shrimp species and identified potential factors that result in their poor genome assembly, including heterozygosity, polyploidization, and repeats. Genome sequencing and comparison of somatic cells (diploid) of the four shrimp species and a single sperm cell (haploid) of Litopenaeus vannamei identified a common bimodal distribution of K-mer depths, suggesting either high heterozygosity or abundant homo-duplicated sequences present in their genomes. However, penaeids have not undergone whole-genome duplication as indicated by a series of approaches. Besides, the remarkable expansion of simple sequence repeats was another outstanding character of penaeid genomes, which also made the genome assembly highly fragmented. Due to this situation, we tried to assemble the genome of penaeid shrimp using various genome sequencing and assembly strategies and compared the quality. Therefore, this study provides new insights about the genomic characteristics of penaeid shrimps while improving their genome assemblies.

Chromosome genomics uncovers plant genome organization and function

The identification of causal genomic loci and their interactions underlying various traits in plants has been greatly aided by progress in understanding the organization of the nuclear genome. This provides clues to the responses of plants to environmental stimuli at the molecular level. Apart from other uses, these insights are needed to fully explore the potential of new breeding techniques that rely on genome editing. However, genome analysis and sequencing is not straightforward in the many agricultural crops and their wild relatives that possess large and complex genomes. Chromosome genomics streamlines this task by dissecting the genome to single chromosomes whose DNA is then used instead of nuclear DNA. This results in a massive and lossless reduction in DNA sample complexity, reduces the time and cost of the experiment, and simplifies data interpretation. Flow cytometric sorting of condensed mitotic chromosomes makes it possible to purify single chromosomes in large quantities, and as the DNA remains intact this process can be coupled successfully with many techniques in molecular biology and genomics. Since the first experiments with flow cytometric sorting in the late 1980s, numerous applications have been developed, and chromosome genomics has been having a significant impact in many areas of research, including the sequencing of complex genomes of important crops and gene cloning. This review discusses these applications, describes their contribution to advancements in plant genome analysis and gene cloning, and outlines future directions.

Complete representation of a tapeworm genome reveals chromosomes capped by centromeres, necessitating a dual role in segregation and protection

BACKGROUND

Chromosome-level assemblies are indispensable for accurate gene prediction, synteny assessment, and understanding higher-order genome architecture. Reference and draft genomes of key helminth species have been published, but little is yet known about the biology of their chromosomes. Here, we present the complete genome of the tapeworm Hymenolepis microstoma, providing a reference quality, end-to-end assembly that represents the first fully assembled genome of a spiralian/lophotrochozoan, revealing new insights into chromosome evolution.

RESULTS

Long-read sequencing and optical mapping data were added to previous short-read data enabling complete re-assembly into six chromosomes, consistent with karyology. Small genome size (169 Mb) and lack of haploid variation (1 SNP/3.2 Mb) contributed to exceptionally high contiguity with only 85 gaps remaining in regions of low complexity sequence. Resolution of repeat regions reveals novel gene expansions, micro-exon genes, and spliced leader trans-splicing, and illuminates the landscape of transposable elements, explaining observed length differences in sister chromatids. Syntenic comparison with other parasitic flatworms shows conserved ancestral linkage groups indicating that the H. microstoma karyotype evolved through fusion events. Strikingly, the assembly reveals that the chromosomes terminate in centromeric arrays, indicating that these motifs play a role not only in segregation, but also in protecting the linear integrity and full lengths of chromosomes.

CONCLUSIONS

Despite strong conservation of canonical telomeres, our results show that they can be substituted by more complex, species-specific sequences, as represented by centromeres. The assembly provides a robust platform for investigations that require complete genome representation.

3D Printing-Enabled DNA Extraction for Long-Read Genomics

Long-read genomics technologies such as nanopore sequencing and genome mapping in nanochannels extract genomic information in the kilobase to megabase pair range from single DNA molecules, thereby overcoming read-length limitations in next-generation DNA sequencing. Long-read technologies start with long DNA molecules as the input and thus benefit from universal sample preparation methods that are fast and shear-free and present a scope of automation and direct upstream integration. We describe a 3D printing-assisted poly(dimethylysiloxane)-based DNA sample preparation device, where diffusive chemical lysis followed by electrophoresis produces circa 100 ng of long DNA directly from cells with less than 5 min of labor. Assessment of the product DNA by confinement in nanochannels reveals that the DNA sizes are commensurate with the requirements for long-read single-molecule technologies. Microfluidics not only expedites sample preparation, but also offers the opportunity for integration with genomics technologies to eliminate DNA fragmentation and loss during transfer to the genomic device.

Karyotypic Evolution of Sauropsid Vertebrates Illuminated by Optical and Physical Mapping of the Painted Turtle and Slider Turtle Genomes

Recent sequencing and software enhancements have advanced our understanding of the evolution of genomic structure and function, especially addressing novel evolutionary biology questions. Yet fragmentary turtle genome assemblies remain a challenge to fully decipher the genetic architecture of adaptive evolution. Here, we use optical mapping to improve the contiguity of the painted turtle (Chrysemys picta) genome assembly and use de novo fluorescent in situ hybridization (FISH) of bacterial artificial chromosome (BAC) clones, BAC-FISH, to physically map the genomes of the painted and slider turtles (Trachemys scripta elegans). Optical mapping increased C. picta's N50 by ~242% compared to the previous assembly. Physical mapping permitted anchoring ~45% of the genome assembly, spanning 5544 genes (including 20 genes related to the sex determination network of turtles and vertebrates). BAC-FISH data revealed assembly errors in C. picta and T. s. elegans assemblies, highlighting the importance of molecular cytogenetic data to complement bioinformatic approaches. We also compared C. picta's anchored scaffolds to the genomes of other chelonians, chicken, lizards, and snake. Results revealed a mostly one-to-one correspondence between chromosomes of painted and slider turtles, and high homology among large syntenic blocks shared with other turtles and sauropsids. Yet, numerous chromosomal rearrangements were also evident across chelonians, between turtles and squamates, and between avian and non-avian reptiles.

Advances in optical mapping for genomic research

Recent advances in optical mapping have allowed the construction of improved genome assemblies with greater contiguity. Optical mapping also enables genome comparison and identification of large-scale structural variations. Association of these large-scale genomic features with biological functions is an important goal in plant and animal breeding and in medical research. Optical mapping has also been used in microbiology and still plays an important role in strain typing and epidemiological studies. Here, we review the development of optical mapping in recent decades to illustrate its importance in genomic research. We detail its applications and algorithms to show its specific advantages. Finally, we discuss the challenges required to facilitate the optimization of optical mapping and improve its future development and application.

Genome assembly provides insights into the genome evolution and flowering regulation of orchardgrass

Orchardgrass (Dactylis glomerata L.) is an important forage grass for cultivating livestock worldwide. Here, we report an ~1.84-Gb chromosome-scale diploid genome assembly of orchardgrass, with a contig N50 of 0.93 Mb, a scaffold N50 of 6.08 Mb and a super-scaffold N50 of 252.52 Mb, which is the first chromosome-scale assembled genome of a cool-season forage grass. The genome includes 40 088 protein-coding genes, and 69% of the assembled sequences are transposable elements, with long terminal repeats (LTRs) being the most abundant. The LTRretrotransposons may have been activated and expanded in the grass genome in response to environmental changes during the Pleistocene between 0 and 1 million years ago. Phylogenetic analysis reveals that orchardgrass diverged after rice but before three Triticeae species, and evolutionarily conserved chromosomes were detected by analysing ancient chromosome rearrangements in these grass species. We also resequenced the whole genome of 76 orchardgrass accessions and found that germplasm from Northern Europe and East Asia clustered together, likely due to the exchange of plants along the 'Silk Road' or other ancient trade routes connecting the East and West. Last, a combined transcriptome, quantitative genetic and bulk segregant analysis provided insights into the genetic network regulating flowering time in orchardgrass and revealed four main candidate genes controlling this trait. This chromosome-scale genome and the online database of orchardgrass developed here will facilitate the discovery of genes controlling agronomically important traits, stimulate genetic improvement of and functional genetic research on orchardgrass and provide comparative genetic resources for other forage grasses.

Chlorella vulgaris genome assembly and annotation reveals the molecular basis for metabolic acclimation to high light conditions

Chlorella vulgaris is a fast-growing fresh-water microalga cultivated on the industrial scale for applications ranging from food to biofuel production. To advance our understanding of its biology and to establish genetics tools for biotechnological manipulation, we sequenced the nuclear and organelle genomes of Chlorella vulgaris 211/11P by combining next generation sequencing and optical mapping of isolated DNA molecules. This hybrid approach allowed us to assemble the nuclear genome in 14 pseudo-molecules with an N50 of 2.8 Mb and 98.9% of scaffolded genome. The integration of RNA-seq data obtained at two different irradiances of growth (high light, HL versus low light, LL) enabled us to identify 10 724 nuclear genes, coding for 11 082 transcripts. Moreover, 121 and 48 genes, respectively, were found in the chloroplast and mitochondrial genome. Functional annotation and expression analysis of nuclear, chloroplast and mitochondrial genome sequences revealed particular features of Chlorella vulgaris. Evidence of horizontal gene transfers from chloroplast to mitochondrial genome was observed. Furthermore, comparative transcriptomic analyses of LL versus HL provided insights into the molecular basis for metabolic rearrangement under HL versus LL conditions leading to enhanced de novo fatty acid biosynthesis and triacylglycerol accumulation. The occurrence of a cytosolic fatty acid biosynthetic pathway could be predicted and its upregulation upon HL exposure was observed, consistent with the increased lipid amount under HL conditions. These data provide a rich genetic resource for future genome editing studies, and potential targets for biotechnological manipulation of Chlorella vulgaris or other microalgae species to improve biomass and lipid productivity.

Haplotype-resolved genomes of geminivirus-resistant and geminivirus-susceptible African cassava cultivars

Background

Cassava is an important food crop in tropical and sub-tropical regions worldwide. In Africa, cassava production is widely affected by cassava mosaic disease (CMD), which is caused by the African cassava mosaic geminivirus that is transmitted by whiteflies. Cassava breeders often use a single locus, CMD2, for introducing CMD resistance into susceptible cultivars. The CMD2 locus has been genetically mapped to a 10-Mbp region, but its organization and genes as well as their functions are unknown.

Results

We report haplotype-resolved de novo assemblies and annotations of the genomes for the African cassava cultivar TME (tropical Manihot esculenta), which is the origin of CMD2, and the CMD-susceptible cultivar 60444. The assemblies provide phased haplotype information for over 80% of the genomes. Haplotype comparison identified novel features previously hidden in collapsed and fragmented cassava genomes, including thousands of allelic variants, inter-haplotype diversity in coding regions, and patterns of diversification through allele-specific expression. Reconstruction of the CMD2 locus revealed a highly complex region with nearly identical gene sets but limited microsynteny between the two cultivars.

Conclusions

The genome maps of the CMD2 locus in both 60444 and TME3, together with the newly annotated genes, will help the identification of the causal genetic basis of CMD2 resistance to geminiviruses. Our de novo cassava genome assemblies will also facilitate genetic mapping approaches to narrow the large CMD2 region to a few candidate genes for better informed strategies to develop robust geminivirus resistance in susceptible cassava cultivars.

Electronic supplementary material

The online version of this article (10.1186/s12915-019-0697-6) contains supplementary material, which is available to authorized users.

Accessing a Russian Wheat Aphid Resistance Gene in Bread Wheat by Long-Read Technologies

Russian wheat aphid (RWA) ( Kurdjumov) is a serious invasive pest of small-grain cereals and many grass species. An efficient strategy to defy aphid attacks is to identify sources of natural resistance and transfer resistance genes into susceptible crop cultivars. Revealing the genes helps understand plant defense mechanisms and engineer plants with durable resistance to the pest. To date, more than 15 RWA resistance genes have been identified in wheat ( L.) but none of them has been cloned. Previously, we genetically mapped the RWA resistance gene into an interval of 0.83 cM on the short arm of chromosome 7D and spanned it with five bacterial artificial chromosome (BAC) clones. Here, we used a targeted strategy combining traditional approaches toward gene cloning (genetic mapping and sequencing of BAC clones) with novel technologies, including optical mapping and long-read nanopore sequencing. The latter, with reads spanning the entire length of a BAC insert, enabled us to assemble the whole region, a task that was not achievable with short reads. Long-read optical mapping validated the DNA sequence in the interval and revealed a difference in the locus organization between resistant and susceptible genotypes. The complete and accurate sequence of the region facilitated the identification of new markers and precise annotation of the interval, revealing six high-confidence genes. Identification of as the most likely candidate opens an avenue for its validation through functional genomics approaches.

The Dark Matter of Large Cereal Genomes: Long Tandem Repeats

Reference genomes of important cereals, including barley, emmer wheat and bread wheat, were released recently. Their comparison with genome size estimates obtained by flow cytometry indicated that the assemblies represent not more than 88-98% of the complete genome. This work is aimed at identifying the missing parts in two cereal genomes and proposing techniques to make the assemblies more complete. We focused on tandemly organised repetitive sequences, known to be underrepresented in genome assemblies generated from short-read sequence data. Our study found arrays of three tandem repeats with unit sizes of 1242 to 2726 bp present in the bread wheat reference genome generated from short reads. However, this and another wheat genome assembly employing long PacBio reads failed in integrating correctly the 2726-bp repeat in the pseudomolecule context. This suggests that tandem repeats of this size, frequently incorporated in unassigned scaffolds, may contribute to shrinking of pseudomolecules without reducing size of the entire assembly. We demonstrate how this missing information may be added to the pseudomolecules with the aid of nanopore sequencing of individual BAC clones and optical mapping. Using the latter technique, we identified and localised a 470-kb long array of 45S ribosomal DNA absent from the reference genome of barley.

GOOGA: A platform to synthesize mapping experiments and identify genomic structural diversity

Understanding genomic structural variation such as inversions and translocations is a key challenge in evolutionary genetics. We develop a novel statistical approach to comparative genetic mapping to detect large-scale structural mutations from low-level sequencing data. The procedure, called Genome Order Optimization by Genetic Algorithm (GOOGA), couples a Hidden Markov Model with a Genetic Algorithm to analyze data from genetic mapping populations. We demonstrate the method using both simulated data (calibrated from experiments on Drosophila melanogaster) and real data from five distinct crosses within the flowering plant genus Mimulus. Application of GOOGA to the Mimulus data corrects numerous errors (misplaced sequences) in the M. guttatus reference genome and confirms or detects eight large inversions polymorphic within the species complex. Finally, we show how this method can be applied in genomic scans to improve the accuracy and resolution of Quantitative Trait Locus (QTL) mapping.

Characterization and evolutionary dynamics of complex regions in eukaryotic genomes

Complex regions in eukaryotic genomes are typically characterized by duplications of chromosomal stretches that often include one or more genes repeated in a tandem array or in relatively close proximity. Nevertheless, the repetitive nature of these regions, together with the often high sequence identity among repeats, have made complex regions particularly recalcitrant to proper molecular characterization, often being misassembled or completely absent in genome assemblies. This limitation has prevented accurate functional and evolutionary analyses of these regions. This is becoming increasingly relevant as evidence continues to support a central role for complex genomic regions in explaining human disease, developmental innovations, and ecological adaptations across phyla. With the advent of long-read sequencing technologies and suitable assemblers, the development of algorithms that can accommodate sample heterozygosity, and the adoption of a pangenomic-like view of these regions, accurate reconstructions of complex regions are now within reach. These reconstructions will finally allow for accurate functional and evolutionary studies of complex genomic regions, underlying the generation of genotype-phenotype maps of unprecedented resolution.

Microfluidic long DNA sample preparation from cells

A number of outstanding problems in genomics, such as identifying structural variations and sequencing through centromeres and telomeres, stand poised to benefit tremendously from emerging long-read genomics technologies such as nanopore sequencing and genome mapping in nanochannels. However, optimal application of these new genomics technologies requires facile methods for extracting long DNA from cells. These sample preparation tools should be amenable to automation and minimize fragmentation of the long DNA molecules by shear. We present one such approach in a poly(dimethylsiloxane) device, where gel-based high molecular weight DNA extraction and continuous flow purification in a 3D cell culture-inspired geometry is followed by electrophoretic extraction of the long DNA from the miniaturized gel. Molecular combing reveals that the device produces molecules that are typically in excess of 100 kilobase pairs in size, with the longest molecule extending up to 4 megabase pairs. The microfluidic format reduces the standard day-long and labor-intensive DNA extraction process to 4 hours, making it a promising prototype platform for routine long DNA sample preparation.

Development of a Sequence-Based Reference Physical Map of Pea (Pisum sativum L.)

Whole genome profiling (WGP) is a sequence-based physical mapping technology and uses sequence tags generated by next generation sequencing for construction of bacterial artificial chromosome (BAC) contigs of complex genomes. The physical map provides a framework for assembly of genome sequence and information for localization of genes that are difficult to find through positional cloning. To address the challenges of accurate assembly of the pea genome (∼4.2 GB of which approximately 85% is repetitive sequences), we have adopted the WGP technology for assembly of a pea BAC library. Multi-dimensional pooling of 295,680 BAC clones and sequencing the ends of restriction fragments of pooled DNA generated 1,814 million high quality reads, of which 825 million were deconvolutable to 1.11 million unique WGP sequence tags. These WGP tags were used to assemble 220,013 BACs into contigs. Assembly of the BAC clones using the modified Fingerprinted Contigs (FPC) program has resulted in 13,040 contigs, consisting of 213,719 BACs, and 6,294 singleton BACs. The average contig size is 0.33 Mbp and the N₅₀ contig size is 0.62 Mbp. WGP^TM technology has proved to provide a robust physical map of the pea genome, which would have been difficult to assemble using traditional restriction digestion based methods. This sequence-based physical map will be useful to assemble the genome sequence of pea. Additionally, the 1.1 million WGP tags will support efficient assignment of sequence scaffolds to the BAC clones, and thus an efficient sequencing of BAC pools with targeted genome regions of interest.

Hardwood Tree Genomics: Unlocking Woody Plant Biology

Woody perennial angiosperms (i.e., hardwood trees) are polyphyletic in origin and occur in most angiosperm orders. Despite their independent origins, hardwoods have shared physiological, anatomical, and life history traits distinct from their herbaceous relatives. New high-throughput DNA sequencing platforms have provided access to numerous woody plant genomes beyond the early reference genomes of Populus and Eucalyptus, references that now include willow and oak, with pecan and chestnut soon to follow. Genomic studies within these diverse and undomesticated species have successfully linked genes to ecological, physiological, and developmental traits directly. Moreover, comparative genomic approaches are providing insights into speciation events while large-scale DNA resequencing of native collections is identifying population-level genetic diversity responsible for variation in key woody plant biology across and within species. Current research is focused on developing genomic prediction models for breeding, defining speciation and local adaptation, detecting and characterizing somatic mutations, revealing the mechanisms of gender determination and flowering, and application of systems biology approaches to model complex regulatory networks underlying quantitative traits. Emerging technologies such as single-molecule, long-read sequencing is being employed as additional woody plant species, and genotypes within species, are sequenced, thus enabling a comparative ("evo-devo") approach to understanding the unique biology of large woody plants. Resource availability, current genomic and genetic applications, new discoveries and predicted future developments are illustrated and discussed for poplar, eucalyptus, willow, oak, chestnut, and pecan.

Discovery of multi-megabase polymorphic inversions by chromosome conformation capture sequencing in large-genome plant species

Chromosomal inversions occur in natural populations of many species, and may underlie reproductive isolation and local adaptation. Traditional methods of inversion discovery are labor-intensive and lack sensitivity. Here, we report the use of three-dimensional contact probabilities between genomic loci as assayed by chromosome-conformation capture sequencing (Hi-C) to detect multi-megabase polymorphic inversions in four barley genotypes. Inversions are validated by fluorescence in situ hybridization and Bionano optical mapping. We propose Hi-C as a generally applicable method for inversion discovery in natural populations.

Genomic approaches for studying crop evolution

Understanding how crop plants evolved from their wild relatives and spread around the world can inform about the origins of agriculture. Here, we review how the rapid development of genomic resources and tools has made it possible to conduct genetic mapping and population genetic studies to unravel the molecular underpinnings of domestication and crop evolution in diverse crop species. We propose three future avenues for the study of crop evolution: establishment of high-quality reference genomes for crops and their wild relatives; genomic characterization of germplasm collections; and the adoption of novel methodologies such as archaeogenetics, epigenomics, and genome editing.

Optical and physical mapping with local finishing enables megabase-scale resolution of agronomically important regions in the wheat genome

BACKGROUND

Numerous scaffold-level sequences for wheat are now being released and, in this context, we report on a strategy for improving the overall assembly to a level comparable to that of the human genome.

RESULTS

Using chromosome 7A of wheat as a model, sequence-finished megabase-scale sections of this chromosome were established by combining a new independent assembly using a bacterial artificial chromosome (BAC)-based physical map, BAC pool paired-end sequencing, chromosome-arm-specific mate-pair sequencing and Bionano optical mapping with the International Wheat Genome Sequencing Consortium RefSeq v1.0 sequence and its underlying raw data. The combined assembly results in 18 super-scaffolds across the chromosome. The value of finished genome regions is demonstrated for two approximately 2.5 Mb regions associated with yield and the grain quality phenotype of fructan carbohydrate grain levels. In addition, the 50 Mb centromere region analysis incorporates cytological data highlighting the importance of non-sequence data in the assembly of this complex genome region.

CONCLUSIONS

Sufficient genome sequence information is shown to now be available for the wheat community to produce sequence-finished releases of each chromosome of the reference genome. The high-level completion identified that an array of seven fructosyl transferase genes underpins grain quality and that yield attributes are affected by five F-box-only-protein-ubiquitin ligase domain and four root-specific lipid transfer domain genes. The completed sequence also includes the centromere.

Shifting the limits in wheat research and breeding using a fully annotated reference genome

An annotated reference sequence representing the hexaploid bread wheat genome in 21 pseudomolecules has been analyzed to identify the distribution and genomic context of coding and noncoding elements across the A, B, and D subgenomes. With an estimated coverage of 94% of the genome and containing 107,891 high-confidence gene models, this assembly enabled the discovery of tissue- and developmental stage-related coexpression networks by providing a transcriptome atlas representing major stages of wheat development. Dynamics of complex gene families involved in environmental adaptation and end-use quality were revealed at subgenome resolution and contextualized to known agronomic single-gene or quantitative trait loci. This community resource establishes the foundation for accelerating wheat research and application through improved understanding of wheat biology and genomics-assisted breeding.

Intact DNA purified from flow-sorted nuclei unlocks the potential of next-generation genome mapping and assembly in Solanum species

Next-generation genome mapping through nanochannels (Bionano optical mapping) of plant genomes brings genome assemblies to the ‘nearly-finished’ level for reliable and detailed gene annotations and assessment of structural variations. Despite the recent progress in its development, researchers face the technical challenges of obtaining sufficient high molecular weight (HMW) nuclear DNA due to cell walls which are difficult to disrupt and to the presence of cytoplasmic polyphenols and polysaccharides that co-precipitate or are covalently bound to DNA and might cause oxidation and/or affect the access of nicking enzymes to DNA, preventing downstream applications. Here we describe important improvements for obtaining HMW DNA that we tested on Solanum crops and wild relatives. The methods that we further elaborated and refined focus on

•

Improving flexibility of using different tissues as source materials, like fast-growing root tips and young leaves from seedlings or in vitro plantlets.

•

Obtaining nuclei suspensions through either lab homogenizers or by chopping.

•

Increasing flow sorting efficiency using DAPI (4′,6-diamidino-2-phenylindole) and PI (propidium iodide) DNA stains, with different lasers (UV or 488 nm) and sorting platforms such as the FACSAria and FACSVantage flow sorters, thus making it appropriate for more laboratories working on plant genomics.

The obtained nuclei are embedded into agarose plugs for processing and isolating uncontaminated HMW DNA, which is a prerequisite for nanochannel-based next-generation optical mapping strategies.

Chromosome-scale pseudomolecules refined by optical, physical and genetic maps in flax

Genomes of varying sizes have been sequenced with next-generation sequencing platforms. However, most reference sequences include draft unordered scaffolds containing chimeras caused by mis-scaffolding. A BioNano genome (BNG) optical map was constructed to improve the previously sequenced flax genome (Linum usitatissimum L., 2n = 30, about 373 Mb), which consisted of 3852 scaffolds larger than 1 kb and totalling 300.6 Mb. The high-resolution BNG map of cv. CDC Bethune totalled 317 Mb and consisted of 251 BNG contigs with an N50 of 2.15 Mb. A total of 622 scaffolds (286.6 Mb, 94.9%) aligned to 211 BNG contigs (298.6 Mb, 94.2%). Of those, 99 scaffolds, diagnosed to contain assembly errors, were refined into 225 new scaffolds. Using the newly refined scaffold sequences and the validated bacterial artificial chromosome-based physical map of CDC Bethune, the 211 BNG contigs were scaffolded into 94 super-BNG contigs (N50 of 6.64 Mb) that were further assigned to the 15 flax chromosomes using the genetic map. The pseudomolecules total about 316 Mb, with individual chromosomes of 15.6 to 29.4 Mb, and cover 97% of the annotated genes. Evidence from the chromosome-scale pseudomolecules suggests that flax has undergone palaeopolyploidization and mesopolyploidization events, followed by rearrangements and deletions or fusion of chromosome arms from an ancient progenitor with a haploid chromosome number of eight.

Improving Nelumbo nucifera genome assemblies using high-resolution genetic maps and BioNano genome mapping reveals ancient chromosome rearrangements

Genetic and physical maps are powerful tools to anchor fragmented draft genome assemblies generated from next-generation sequencing. Currently, two draft assemblies of Nelumbo nucifera, the genomes of 'China Antique' and 'Chinese Tai-zi', have been released. However, there is presently no information on how the sequences are assembled into chromosomes in N. nucifera. The lack of physical maps and inadequate resolution of available genetic maps hindered the assembly of N. nucifera chromosomes. Here, a linkage map of N. nucifera containing 2371 bin markers [217 577 single nucleotide polymorphisms (SNPs)] was constructed using restriction-site associated DNA sequencing data of 181 F₂ individuals and validated by adding 197 simple sequence repeat (SSR) markers. Additionally, a BioNano optical map covering 86.20% of the 'Chinese Tai-zi' genome was constructed. The draft assembly of 'Chinese Tai-zi' was improved based on the BioNano optical map, showing an increase of the scaffold N50 from 0.989 to 1.48 Mb. Using a combination of multiple maps, 97.9% of the scaffolds in the 'Chinese Tai-zi' draft assembly and 97.6% of the scaffolds in the 'China Antique' draft assembly were anchored into pseudo-chromosomes, and the centromere regions along the pseudo-chromosomes were identified. An evolutionary scenario was proposed to reach the modern N. nucifera karyotype from the seven ancestral eudicot chromosomes. The present study provides the highest-resolution linkage map, the optical map and chromosome level genome assemblies for N. nucifera, which are valuable for the breeding and cultivation of N. nucifera and future studies of comparative and evolutionary genomics in angiosperms.

A guide to sequence your favorite plant genomes

With the rapid development of sequencing technology and the plummeting cost, assembling whole genomes from non-model plants will soon become routine for plant systematists and evolutionary biologists. Here we summarize and compare several of the latest genome sequencing and assembly approaches, offering a practical guide on how to approach a genome project. We also highlight certain precautions that need to be taken before investing time and money into a genome project.

Post genomics era for orchid research

Among 300,000 species in angiosperms, Orchidaceae containing 30,000 species is one of the largest families. Almost every habitats on earth have orchid plants successfully colonized, and it indicates that orchids are among the plants with significant ecological and evolutionary importance. So far, four orchid genomes have been sequenced, including Phalaenopsis equestris, Dendrobium catenatum, Dendrobium officinale, and Apostaceae shengen. Here, we review the current progress and the direction of orchid research in the post genomics era. These include the orchid genome evolution, genome mapping (genome-wide association analysis, genetic map, physical map), comparative genomics (especially receptor-like kinase and terpene synthase), secondary metabolomics, and genome editing.

Fine organization of genomic regions tagged to the 5S rDNA locus of the bread wheat 5B chromosome

BACKGROUND

The multigene family encoding the 5S rRNA, one of the most important structurally-functional part of the large ribosomal subunit, is an obligate component of all eukaryotic genomes. 5S rDNA has long been a favored target for cytological and phylogenetic studies due to the inherent peculiarities of its structural organization, such as the tandem arrays of repetitive units and their high interspecific divergence. The complex polyploid nature of the genome of bread wheat, Triticum aestivum, and the technically difficult task of sequencing clusters of tandem repeats mean that the detailed organization of extended genomic regions containing 5S rRNA genes remains unclear. This is despite the recent progress made in wheat genomic sequencing. Using pyrosequencing of BAC clones, in this work we studied the organization of two distinct 5S rDNA-tagged regions of the 5BS chromosome of bread wheat.

RESULTS

Three BAC-clones containing 5S rDNA were identified in the 5BS chromosome-specific BAC-library of Triticum aestivum. Using the results of pyrosequencing and assembling, we obtained six 5S rDNA- containing contigs with a total length of 140,417 bp, and two sets (pools) of individual 5S rDNA sequences belonging to separate, but closely located genomic regions on the 5BS chromosome. Both regions are characterized by the presence of approximately 70-80 copies of 5S rDNA, however, they are completely different in their structural organization. The first region contained highly diverged short-type 5S rDNA units that were disrupted by multiple insertions of transposable elements. The second region contained the more conserved long-type 5S rDNA, organized as a single tandem array. FISH using probes specific to both 5S rDNA unit types showed differences in the distribution and intensity of signals on the chromosomes of polyploid wheat species and their diploid progenitors.

CONCLUSION

A detailed structural organization of two closely located 5S rDNA-tagged genomic regions on the 5BS chromosome of bread wheat has been established. These two regions differ in the organization of both 5S rDNA and the neighboring sequences comprised of transposable elements, implying different modes of evolution for these regions.

Multiple hybrid de novo genome assembly of finger millet, an orphan allotetraploid crop

Finger millet (Eleusine coracana (L.) Gaertn) is an important crop for food security because of its tolerance to drought, which is expected to be exacerbated by global climate changes. Nevertheless, it is often classified as an orphan/underutilized crop because of the paucity of scientific attention. Among several small millets, finger millet is considered as an excellent source of essential nutrient elements, such as iron and zinc; hence, it has potential as an alternate coarse cereal. However, high-quality genome sequence data of finger millet are currently not available. One of the major problems encountered in the genome assembly of this species was its polyploidy, which hampers genome assembly compared with a diploid genome. To overcome this problem, we sequenced its genome using diverse technologies with sufficient coverage and assembled it via a novel multiple hybrid assembly workflow that combines next-generation with single-molecule sequencing, followed by whole-genome optical mapping using the Bionano Irys® system. The total number of scaffolds was 1,897 with an N50 length >2.6 Mb and detection of 96% of the universal single-copy orthologs. The majority of the homeologs were assembled separately. This indicates that the proposed workflow is applicable to the assembly of other allotetraploid genomes.

Strategies for optimizing BioNano and Dovetail explored through a second reference quality assembly for the legume model, Medicago truncatula

Background

Third generation sequencing technologies, with sequencing reads in the tens- of kilo-bases, facilitate genome assembly by spanning ambiguous regions and improving continuity. This has been critical for plant genomes, which are difficult to assemble due to high repeat content, gene family expansions, segmental and tandem duplications, and polyploidy. Recently, high-throughput mapping and scaffolding strategies have further improved continuity. Together, these long-range technologies enable quality draft assemblies of complex genomes in a cost-effective and timely manner.

Results

Here, we present high quality genome assemblies of the model legume plant, Medicago truncatula (R108) using PacBio, Dovetail Chicago (hereafter, Dovetail) and BioNano technologies. To test these technologies for plant genome assembly, we generated five assemblies using all possible combinations and ordering of these three technologies in the R108 assembly. While the BioNano and Dovetail joins overlapped, they also showed complementary gains in continuity and join numbers. Both technologies spanned repetitive regions that PacBio alone was unable to bridge. Combining technologies, particularly Dovetail followed by BioNano, resulted in notable improvements compared to Dovetail or BioNano alone. A combination of PacBio, Dovetail, and BioNano was used to generate a high quality draft assembly of R108, a M. truncatula accession widely used in studies of functional genomics. As a test for the usefulness of the resulting genome sequence, the new R108 assembly was used to pinpoint breakpoints and characterize flanking sequence of a previously identified translocation between chromosomes 4 and 8, identifying more than 22.7 Mb of novel sequence not present in the earlier A17 reference assembly.

Conclusions

Adding Dovetail followed by BioNano data yielded complementary improvements in continuity over the original PacBio assembly. This strategy proved efficient and cost-effective for developing a quality draft assembly compared to traditional reference assemblies.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-017-3971-4) contains supplementary material, which is available to authorized users.

An advanced reference genome of Trifolium subterraneum L. reveals genes related to agronomic performance

Subterranean clover is an important annual forage legume, whose diploidy and inbreeding nature make it an ideal model for genomic analysis in Trifolium. We reported a draft genome assembly of the subterranean clover TSUd_r1.1. Here we evaluate genome mapping on nanochannel arrays and generation of a transcriptome atlas across tissues to advance the assembly and gene annotation. Using a BioNano-based assembly spanning 512 Mb (93% genome coverage), we validated the draft assembly, anchored unplaced contigs and resolved misassemblies. Multiple contigs (264) from the draft assembly coalesced into 97 super-scaffolds (43% of genome). Sequences longer than >1 Mb increased from 40 to 189 Mb giving 1.4-fold increase in N50 with total genome in pseudomolecules improved from 73 to 80%. The advanced assembly was re-annotated using transcriptome atlas data to contain 31 272 protein-coding genes capturing >96% of the gene content. Functional characterization and GO enrichment confirmed gene expression for response to water deprivation, flavonoid biosynthesis and embryo development ending in seed dormancy, reflecting adaptation to the harsh Mediterranean environment. Comparative analyses across Papilionoideae identified 24 893 Trifolium-specific and 6325 subterranean-clover-specific genes that could be mined further for traits such as geocarpy and grazing tolerance. Eight key traits, including persistence, improved livestock health by isoflavonoid production in addition to important agro-morphological traits, were fine-mapped on the high-density SNP linkage map anchored to the assembly. This new genomic information is crucial to identify loci governing traits allowing marker-assisted breeding, comparative mapping and identification of tissue-specific gene promoters for biotechnological improvement of forage legumes.

Physical Map of the Short Arm of Bread Wheat Chromosome 3D

Bread wheat ( L.) is one of the most important crops worldwide. Although a reference genome sequence would represent a valuable resource for wheat improvement through genomics-assisted breeding and gene cloning, its generation has long been hampered by its allohexaploidy, high repeat content, and large size. As a part of a project coordinated by the International Wheat Genome Sequencing Consortium (IWGSC), a physical map of the short arm of wheat chromosome 3D (3DS) was prepared to facilitate reference genome assembly and positional gene cloning. It comprises 869 contigs with a cumulative length of 274.5 Mbp and represents 85.5% of the estimated chromosome arm size. Eighty-six Mbp of survey sequences from chromosome arm 3DS were assigned in silico to physical map contigs via next-generation sequencing of bacterial artificial chromosome pools, thus providing a high-density framework for physical map ordering along the chromosome arm. About 60% of the physical map was anchored in this single experiment. Finally, 1393 high-confidence genes were anchored to the physical map. Comparisons of gene space of the chromosome arm 3DS with genomes of closely related species [ (L.) P.Beauv., rice ( L.), and sorghum [ (L.) Moench] and homeologous wheat chromosomes provided information about gene movement on the chromosome arm.

BioNanoAnalyst: a visualisation tool to assess genome assembly quality using BioNano data

Background

Reference genome assemblies are valuable, as they provide insights into gene content, genetic evolution and domestication. The higher the quality of a reference genome assembly the more accurate the downstream analysis will be. During the last few years, major efforts have been made towards improving the quality of genome assemblies. However, erroneous and incomplete assemblies are still common. Complementary to DNA sequencing technologies, optical mapping has advanced genomic studies by facilitating the production of genome scaffolds and assessing structural variation. However, there are few tools available to comprehensively examine misassemblies in reference genome sequences using optical map data.

Results

We present BioNanoAnalyst, a software package to examine genome assemblies based on restriction endonuclease cut sites and optical map data. A graphical user interface (GUI) allows users to assess reference genome sequences on different computer platforms without the requirement of programming knowledge. The zoom function makes visualisation convenient, while a GFF3 format output file gives an option to directly visualise questionable assembly regions by location and nucleotides following import into a local genome browser.

Conclusions

BioNanoAnalyst is a tool to identify misassemblies in a reference genome sequence using optical map data. With the reported information, users can rapidly identify assembly errors and correct them using other software tools, which could facilitate an accurate downstream analysis.

Electronic supplementary material

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

Higher-order organisation of extremely amplified, potentially functional and massively methylated 5S rDNA in European pikes (Esox sp.)

Background

Pikes represent an important genus (Esox) harbouring a pre-duplication karyotype (2n = 2x = 50) of economically important salmonid pseudopolyploids. Here, we have characterized the 5S ribosomal RNA genes (rDNA) in Esox lucius and its closely related E. cisalpinus using cytogenetic, molecular and genomic approaches. Intragenomic homogeneity and copy number estimation was carried out using Illumina reads. The higher-order structure of rDNA arrays was investigated by the analysis of long PacBio reads. Position of loci on chromosomes was determined by FISH. DNA methylation was analysed by methylation-sensitive restriction enzymes.

Results

The 5S rDNA loci occupy exclusively (peri)centromeric regions on 30–38 acrocentric chromosomes in both E. lucius and E. cisalpinus. The large number of loci is accompanied by extreme amplification of genes (>20,000 copies), which is to the best of our knowledge one of the highest copy number of rRNA genes in animals ever reported. Conserved secondary structures of predicted 5S rRNAs indicate that most of the amplified genes are potentially functional. Only few SNPs were found in genic regions indicating their high homogeneity while intergenic spacers were more heterogeneous and several families were identified. Analysis of 10–30 kb-long molecules sequenced by the PacBio technology (containing about 40% of total 5S rDNA) revealed that the vast majority (96%) of genes are organised in large several kilobase-long blocks. Dispersed genes or short tandems were less common (4%). The adjacent 5S blocks were directly linked, separated by intervening DNA and even inverted. The 5S units differing in the intergenic spacers formed both homogeneous and heterogeneous (mixed) blocks indicating variable degree of homogenisation between the loci. Both E. lucius and E. cisalpinus 5S rDNA was heavily methylated at CG dinucleotides.

Conclusions

Extreme amplification of 5S rRNA genes in the Esox genome occurred in the absence of significant pseudogenisation suggesting its recent origin and/or intensive homogenisation processes. The dense methylation of units indicates that powerful epigenetic mechanisms have evolved in this group of fish to silence amplified genes. We discuss how the higher-order repeat structures impact on homogenisation of 5S rDNA in the genome.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-017-3774-7) contains supplementary material, which is available to authorized users.