Using linkage maps to correct and scaffold de novo genome assemblies: methods, challenges, and computational tools

Telomere-to-telomere genome assembly of an allotetraploid pernicious weed, Echinochloa phyllopogon

Echinochloa phyllopogon is an allotetraploid pernicious weed species found in rice fields worldwide that often exhibit resistance to multiple herbicides. An accurate genome sequence is essential to comprehensively understand the genetic basis underlying the traits of this species. Here, the telomere-to-telomere genome sequence of E. phyllopogon was presented. Eighteen chromosome sequences spanning 1.0 Gb were constructed using the PacBio highly fidelity long technology. Of the 18 chromosomes, 12 sequences were entirely assembled into telomere-to-telomere and gap-free contigs, whereas the remaining six sequences were constructed at the chromosomal level with only eight gaps. The sequences were assigned to the A and B genome with total lengths of 453 and 520 Mb, respectively. Repetitive sequences occupied 42.93% of the A genome and 48.47% of the B genome, although 32,337, and 30,889 high-confidence genes were predicted in the A and B genomes, respectively. This suggested that genome extensions and gene disruptions caused by repeated sequence accumulation often occur in the B genome before polyploidization to establish a tetraploid genome. The highly accurate and comprehensive genome sequence could be a milestone in understanding the molecular mechanisms of the pernicious traits and in developing effective weed control strategies to avoid yield loss in rice production.

High-density linkage mapping and genetic dissection of resistance to broomrape (Orobanche crenata Forsk.) in pea (Pisum sativum L.)

Pea (Pisum sativum L.) is a widely cultivated legume of major importance for global food security and agricultural sustainability. Crenate broomrape (Orobanche crenata Forsk.) (Oc) is a parasitic weed severely affecting legumes, including pea, in the Mediterranean Basin and the Middle East. Previously, the identification of the pea line "ROR12", displaying resistance to Oc, was reported. Two-year field trials on a segregant population of 148 F₇ recombinant inbred lines (RILs), originating from a cross between "ROR12" and the susceptible cultivar "Sprinter", revealed high heritability (0.84) of the "ROR12" resistance source. Genotyping-by-sequencing (GBS) on the same RIL population allowed the construction of a high-density pea linkage map, which was compared with the pea reference genome and used for quantitative trait locus (QTL) mapping. Three QTLs associated with the response to Oc infection, named PsOcr-1, PsOcr-2, and PsOcr-3, were identified, with PsOcr-1 explaining 69.3% of the genotypic variance. Evaluation of the effects of different genotypic combinations indicated additivity between PsOcr-1 and PsOcr-2, and between PsOcr-1 and PsOcr-3, and epistasis between PsOcr-2 and PsOcr-3. Finally, three Kompetitive Allele Specific PCR (KASP) marker assays were designed on the single-nucleotide polymorphisms (SNPs) associated with the QTL significance peaks. Besides contributing to the development of pea genomic resources, this work lays the foundation for the obtainment of pea cultivars resistant to Oc and the identification of genes involved in resistance to parasitic Orobanchaceae.

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.

Smooth Descent: A ploidy-aware algorithm to improve linkage mapping in the presence of genotyping errors

Linkage mapping is an approach to order markers based on recombination events. Mapping algorithms cannot easily handle genotyping errors, which are common in high-throughput genotyping data. To solve this issue, strategies have been developed, aimed mostly at identifying and eliminating these errors. One such strategy is SMOOTH, an iterative algorithm to detect genotyping errors. Unlike other approaches, SMOOTH can also be used to impute the most probable alternative genotypes, but its application is limited to diploid species and to markers heterozygous in only one of the parents. In this study we adapted SMOOTH to expand its use to any marker type and to autopolyploids with the use of identity-by-descent probabilities, naming the updated algorithm Smooth Descent (SD). We applied SD to real and simulated data, showing that in the presence of genotyping errors this method produces better genetic maps in terms of marker order and map length. SD is particularly useful for error rates between 5% and 20% and when error rates are not homogeneous among markers or individuals. With a starting error rate of 10%, SD reduced it to ∼5% in diploids, ∼7% in tetraploids and ∼8.5% in hexaploids. Conversely, the correlation between true and estimated genetic maps increased by 0.03 in tetraploids and by 0.2 in hexaploids, while worsening slightly in diploids (∼0.0011). We also show that the combination of genotype curation and map re-estimation allowed us to obtain better genetic maps while correcting wrong genotypes. We have implemented this algorithm in the R package Smooth Descent.

High-Density Genetic Linkage Map of the Southern Blue-ringed Octopus (Octopodidae: Hapalochlaena maculosa)

Genetic linkage maps provide a useful resource for non-model genomes and can aid in genome reassembly to form more contiguous pseudo-chromosomes. We present the first linkage map of any cephalopod, H. maculosa, composed of 47 linkage groups (LG). A total of 2166 single nucleotide polymorphisms and 2455 presence–absence variant loci were utilised by Lep-Map3 in linkage map construction. The map length spans 2016.62 cM with an average marker distance of 0.85 cM. Integration of the recent H. maculosa genome allowed 1151 scaffolds comprising 34% of the total genomic sequence to be orientated and/or placed using 1278 markers across all 47 LG. The linkage map generated provides a new perspective on HOX gene distribution in octopods. In the H. maculosa linkage map three (SCR, LOX4 and POST1) of six identified HOX genes (HOX1/LAB, SCR, LOX2, LOX4, LOX5, POST1) were located within the same LG (LG 9). The generation of a linkage map for H. maculosa has provided a valuable resource for understanding the evolution of cephalopod genomes and will provide a base for future work.

The first high-density genetic map of common cockle (Cerastoderma edule) reveals a major QTL controlling shell color variation

Shell color shows broad variation within mollusc species and despite information on the genetic pathways involved in shell construction and color has recently increased, more studies are needed to understand its genetic architecture. The common cockle (Cerastoderma edule) is a valuable species from ecological and commercial perspectives which shows important variation in shell color across Northeast Atlantic. In this study, we constructed a high-density genetic map, as a tool for screening common cockle genome, which was applied to ascertain the genetic basis of color variation in the species. The consensus genetic map comprised 19 linkage groups (LGs) in accordance with the cockle karyotype (2n = 38) and spanned 1073 cM, including 730 markers per LG and an inter-marker distance of 0.13 cM. Five full-sib families showing segregation for several color-associated traits were used for a genome-wide association study and a major QTL on chromosome 13 associated to different color-traits was detected. Mining on this genomic region revealed several candidate genes related to shell construction and color. A genomic region previously reported associated with divergent selection in cockle distribution overlapped with this QTL suggesting its putative role on adaptation.

A Novel Strategy for Constructing an Integrated Linkage Map in an F1 Hybrid Population of Populus deltoides and Populus simonii

The genetic linkage maps of the traditional F₂ population in inbred lines were estimated from the frequency of recombination events in both parents, providing full genetic information for genetic and genomic studies. However, in outbred forest trees, it is almost impossible to generate the F₂ population because of their high heterozygosity and long generation times. We proposed a novel strategy to construct an integrated genetic linkage map that contained both parental recombination information, with restriction-site-associated DNA sequencing (RADSeq) data in an F1 hybrid population of Populus deltoides and Populus simonii. We selected a large number of specific RAD tags to construct the linkage map, each of which contained two SNPs, one heterozygous only in the female parent and the other heterozygous only in the male. Consequently, the integrated map contained a total of 1154 RAD tags and 19 linkage groups, with a total length of 5255.49 cM and an average genetic distance of 4.63 cM. Meanwhile, the two parent-specific linkage maps were also constructed with SNPs that were heterozygous in one parent and homozygous in the other. We found that the integrated linkage map was more consensus with the genomic sequences of P. simonii and P. deltoides. Additionally, the likelihood of the marker order in each linkage group of the integrated map was greater than that in both parental maps. The integrated linkage map was more accurate than the parent-specific linkage maps constructed in the same F₁ hybrid population, providing a powerful genetic resource for identifying the quantitative trait loci (QTLs) with dominant effects, assembling genomic sequences, and performing comparative genomics in related Populus species. More importantly, this novel strategy can be used in other outbred species to build an integrated linkage map.

Comparative genomics analysis and genome assembly integration with the recombination landscape contribute to Takifugu bimaculatus assembly refinement

Takifugu genus has been brought to the fore in scientific and practical research due to its compact genome, explosive speciation progress and economic value. Here we updated the chromosome-level genome of Takifugu bimaculatus by an ultra-high-density linkage map, a classic and accurate way of chromosome assembly. The map constituted a robust assembly frame, with 92.2% (372.77 Mb) of the draft genome cumulatively placed. With intraspecies and interspecies comparative genomic analysis, we developed a criterion to quantify the differences between assemblies and established a novel way to integrate information from multiple assemblies. The integrated assembly rectified potential mis-assemblies, greatly improving the genome contiguity and correctness. Our results rendered profound information on the genetic recombination of T. bimaculatus and provided new insights into effective genome assembly. The consolidated assembly will be a contributory tool of T. bimaculatus and broadly across the Takifugu by providing a convincing reference for genomic research.

Integrating Genetic and Chromosome Maps of Allium cepa: From Markers Visualization to Genome Assembly Verification

The ability to directly look into genome sequences has opened great opportunities in plant breeding. Yet, the assembly of full-length chromosomes remains one of the most difficult problems in modern genomics. Genetic maps are commonly used in de novo genome assembly and are constructed on the basis of a statistical analysis of the number of recombinations. This may affect the accuracy of the ordering and orientation of scaffolds within the chromosome, especially in the region of recombination suppression. Moreover, it is impossible to assign contigs lacking DNA markers. Here, we report the use of Tyr-FISH to determine the position of the short DNA sequence of markers and non-mapped unique copy sequence on the physical chromosomes of a large-genome onion (Allium cepa L.). In order to minimize potential background masking of the target signal, we improved our earlier developed pipeline for probe design. A total of 23 markers were located on physical chromosomes 2 and 6. The order of markers was corrected by the integration of genetic, pseudochromosome maps and cytogenetic maps. Additionally, the position of the mlh1 gene, which was not on the genetic map, was defined on physical chromosome 2. Tyr-FISH mapping showed that the order of 23.1% (chromosome 2) and 27.3% (chromosome 6) of the tested genes differed between physical chromosomes and pseudochromosomes. The results can be used for the improvement of pseudochromosome 2 and 6 assembly. The present study aims to demonstrate the value of the in situ visualization of DNA sequences in chromosome-scaffold genome assembly.

A Novel Strategy to Reveal the Landscape of Crossovers in an F1 Hybrid Population of Populus deltoides and Populus simonii

Although the crossover (CO) patterns of different species have been extensively investigated, little is known about the landscape of CO patterns in Populus because of its high heterozygosity and long-time generation. A novel strategy was proposed to reveal the difference of CO rate and interference between Populus deltoides and Populus simonii using their F1 hybrid population. We chose restriction site-associated DNA (RAD) tags that contained two SNPs, one only receiving the CO information from the female P. deltoides and the other from the male P. simonii. These RAD tags allowed us to investigate the CO patterns between the two outbred species, instead of using the traditional backcross populations in inbred lines. We found that the CO rate in P. deltoides was generally greater than that in P. simonii, and that the CO interference was a common phenomenon across the two genomes. The COs landscape of the different Populus species facilitates not only to understand the evolutionary mechanism for adaptability but also to rebuild the statistical model for precisely constructing genetic linkage maps that are critical in genome assembly in Populus. Additionally, the novel strategy could be applied in other outbred species for investigating the CO patterns.

De novo chromosome level assembly of a plant genome from long read sequence data

Recent advances in the sequencing and assembly of plant genomes have allowed the generation of genomes with increasing contiguity and sequence accuracy. Chromosome level genome assemblies using sequence contigs generated from long read sequencing have involved the use of proximity analysis (Hi-C) or traditional genetic maps to guide the placement of sequence contigs within chromosomes. The development of highly accurate long reads by repeated sequencing of circularized DNA (HiFi; PacBio) has greatly increased the size of contigs. We now report the use of HiFiasm to assemble the genome of Macadamia jansenii, a genome that has been used as a model to test sequencing and assembly. This achieved almost complete chromosome level assembly from the sequence data alone without the need for higher level chromosome map information. Eight of the 14 chromosomes were represented by a single large contig (six with telomere repeats at both ends) and the other six assembled from two to four main contigs. The small number of chromosome breaks appears to be the result of highly repetitive regions including ribosomal genes that cannot be assembled by these approaches. De novo assembly of near complete chromosome level plant genomes now appears possible using these sequencing and assembly tools. Further targeted strategies might allow these remaining gaps to be closed.

Chromosome-level genome assembly reveals genomic architecture of northern range expansion in the mountain pine beetle, Dendroctonus ponderosae Hopkins (Coleoptera: Curculionidae)

Genome sequencing methods and assembly tools have improved dramatically since the 2013 publication of draft genome assemblies for the mountain pine beetle, Dendroctonus ponderosae Hopkins (Coleoptera: Curculionidae). We conducted proximity ligation library sequencing and scaffolding to improve contiguity, and then used linkage mapping and recent bioinformatic tools for correction and further improvement. The new assemblies have dramatically improved contiguity and gaps compared to the originals: N50 values increased 26- to 36-fold, and the number of gaps were reduced by half. Ninety per cent of the content of the assemblies is now contained in 12 and 11 scaffolds for the female and male assemblies, respectively. Based on linkage mapping information, the 12 largest scaffolds in both assemblies represent all 11 autosomal chromosomes and the neo-X chromosome. These assemblies now have nearly chromosome-sized scaffolds and will be instrumental for studying genomic architecture, chromosome evolution, population genomics, functional genomics, and adaptation in this and other pest insects. We also identified regions in two chromosomes, including the ancestral-X portion of the neo-X chromosome, with elevated differentiation between northern and southern Canadian populations.

Efficient construction of a linkage map and haplotypes for Mentha suaveolens using sequence capture

The sustainability of many crops is hindered by the lack of genomic resources and a poor understanding of natural genetic diversity. Particularly, application of modern breeding requires high-density linkage maps that are integrated into a highly contiguous reference genome. Here, we present a rapid method for deriving haplotypes and developing linkage maps, and its application to Mentha suaveolens, one of the diploid progenitors of cultivated mints. Using sequence-capture via DNA hybridization to target single nucleotide polymorphisms (SNPs), we successfully genotyped ∼5000 SNPs within the genome of >400 individuals derived from a self cross. After stringent quality control, and identification of nonredundant SNPs, 1919 informative SNPs were retained for linkage map construction. The resulting linkage map defined a total genetic space of 942.17 cM divided among 12 linkage groups, ranging from 56.32 to 122.61 cM in length. The linkage map is in good agreement with pseudomolecules from our preliminary genome assembly, proving this resource effective for the correction and validation of the reference genome. We discuss the advantages of this method for the rapid creation of linkage maps.

Sequencing and Chromosome-Scale Assembly of Plant Genomes, Brassica rapa as a Use Case

With the rise of long-read sequencers and long-range technologies, delivering high-quality plant genome assemblies is no longer reserved to large consortia. Not only sequencing techniques, but also computer algorithms have reached a point where the reconstruction of assemblies at the chromosome scale is now feasible at the laboratory scale. Current technologies, in particular long-range technologies, are numerous, and selecting the most promising one for the genome of interest is crucial to obtain optimal results. In this study, we resequenced the genome of the yellow sarson, Brassica rapa cv. Z1, using the Oxford Nanopore PromethION sequencer and assembled the sequenced data using current assemblers. To reconstruct complete chromosomes, we used and compared three long-range scaffolding techniques, optical mapping, Omni-C, and Pore-C sequencing libraries, commercialized by Bionano Genomics, Dovetail Genomics, and Oxford Nanopore Technologies, respectively, or a combination of the three, in order to evaluate the capability of each technology.

Optimizing experimental design for genome sequencing and assembly with Oxford Nanopore Technologies

High quality reference genome sequences are the core of modern genomics. Oxford Nanopore Technologies (ONT) produces inexpensive DNA sequences, but has high error rates, which make sequence assembly and analysis difficult as genome size and complexity increases. Robust experimental design is necessary for ONT genome sequencing and assembly, but few studies have addressed eukaryotic organisms. Here, we present novel results using simulated and empirical ONT and DNA libraries to identify best practices for sequencing and assembly for several model species. We find that the unique error structure of ONT libraries causes errors to accumulate and assembly statistics plateau as sequence depth increases. High-quality assembled eukaryotic sequences require high-molecular-weight DNA extractions that increase sequence read length, and computational protocols that reduce error through pre-assembly correction and read selection. Our quantitative results will be helpful for researchers seeking guidance for de novo assembly projects.

The Assembled and Annotated Genome of the Fairy-Ring Fungus Marasmius oreades

Marasmius oreades is a basidiomycete fungus that grows in so called "fairy rings," which are circular, underground mycelia common in lawns across temperate areas of the world. Fairy rings can be thought of as natural, long-term evolutionary experiments. As each ring has a common origin and expands radially outwards over many years, different sectors will independently accumulate mutations during growth. The genotype can be followed to the next generation, as mushrooms producing the sexual spores are formed seasonally at the edge of the ring. Here, we present new genomic data from 95 single-spore isolates of the species, which we used to construct a genetic linkage map and an updated version of the genome assembly. The 44-Mb assembly was anchored to 11 linkage groups, producing chromosome-length scaffolds. Gene annotation revealed 13,891 genes, 55% of which contained a pfam domain. The repetitive fraction of the genome was 22%, and dominated by retrotransposons and DNA elements of the KDZ and Plavaka groups. The level of assembly contiguity we present is so far rare in mushroom-forming fungi, and we expect studies of genomics, transposons, phylogenetics, and evolution to be facilitated by the data we present here of the iconic fairy-ring mushroom.

The Diversity of Plant Sex Chromosomes Highlighted through Advances in Genome Sequencing

For centuries, scientists have been intrigued by the origin of dioecy in plants, characterizing sex-specific development, uncovering cytological differences between the sexes, and developing theoretical models. Through the invention and continued improvements in genomic technologies, we have truly begun to unlock the genetic basis of dioecy in many species. Here we broadly review the advances in research on dioecy and sex chromosomes. We start by first discussing the early works that built the foundation for current studies and the advances in genome sequencing that have facilitated more-recent findings. We next discuss the analyses of sex chromosomes and sex-determination genes uncovered by genome sequencing. We synthesize these results to find some patterns are emerging, such as the role of duplications, the involvement of hormones in sex-determination, and support for the two-locus model for the origin of dioecy. Though across systems, there are also many novel insights into how sex chromosomes evolve, including different sex-determining genes and routes to suppressed recombination. We propose the future of research in plant sex chromosomes should involve interdisciplinary approaches, combining cutting-edge technologies with the classics to unravel the patterns that can be found across the hundreds of independent origins.

Male-specific Y-chromosomal regions in waterhemp (Amaranthus tuberculatus) and Palmer amaranth (Amaranthus palmeri)

Amaranthus tuberculatus and Amaranthus palmeri are agronomically important weed species, both with stable dioecious reproductive systems. An understanding of the genetic basis of sex determination may lead to new methods of managing these troublesome weeds. Previous research identified genomic sequences associated with maleness in each species. Male-specific sequences were used to identify genomic regions in both species that are believed to contain sex-determining genes, i.e. the male-specific Y (MSY) region. These regions were compared to understand if sex determination is controlled via the same physiological pathway and if dioecy evolved independently. A contiguously assembled candidate MSY region identified in Amaranthus palmeri is approximately 1.3 Mb with 121 predicted gene models. In Amaranthus tuberculatus, several contigs, with combined length of 4.6 Mb and with 147 gene models, were identified as belonging to the MSY region. Synteny was not detected between the two species' candidate MSY regions but they shared two predicted genes. With lists of candidate genes for sex determination containing fewer than 200 in each species, future research can address whether sex determination is controlled via similar physiological pathways and whether dioecy has indeed evolved independently in these species.

Construction of a High-Density Genetic Map of Acca sellowiana (Berg.) Burret, an Outcrossing Species, Based on Two Connected Mapping Populations

Acca sellowiana, known as feijoa or pineapple guava, is a diploid, (2n = 2x = 22) outcrossing fruit tree species native to Uruguay and Brazil. The species stands out for its highly aromatic fruits, with nutraceutical and therapeutic value. Despite its promising agronomical value, genetic studies on this species are limited. Linkage genetic maps are valuable tools for genetic and genomic studies, and constitute essential tools in breeding programs to support the development of molecular breeding strategies. A high-density composite genetic linkage map of A. sellowiana was constructed using two genetically connected populations: H5 (TCO × BR, N = 160) and H6 (TCO × DP, N = 184). Genotyping by sequencing (GBS) approach was successfully applied for developing single nucleotide polymorphism (SNP) markers. A total of 4,921 SNP markers were identified using the reference genome of the closely related species Eucalyptus grandis, whereas other 4,656 SNPs were discovered using a de novo pipeline. The individual H5 and H6 maps comprised 1,236 and 1,302 markers distributed over the expected 11 linkage groups, respectively. These two maps spanned a map length of 1,593 and 1,572 cM, with an average inter-marker distance of 1.29 and 1.21 cM, respectively. A large proportion of markers were common to both maps and showed a high degree of collinearity. The composite map consisted of 1,897 SNPs markers with a total map length of 1,314 cM and an average inter-marker distance of 0.69. A novel approach for the construction of composite maps where the meiosis information of individuals of two connected populations is captured in a single estimator is described. A high-density, accurate composite map based on a consensus ordering of markers provides a valuable contribution for future genetic research and breeding efforts in A. sellowiana. A novel mapping approach based on an estimation of multipopulation recombination fraction described here may be applied in the construction of dense composite genetic maps for any other outcrossing diploid species.

An ultra-high density SNP-based linkage map for enhancing the pikeperch (Sander lucioperca) genome assembly to chromosome-scale

Pikeperch (Sander lucioperca) is a fish species with growing economic significance in the aquaculture industry. However, successful positioning of pikeperch in large-scale aquaculture requires advances in our understanding of its genome organization. In this study, an ultra-high density linkage map for pikeperch comprising 24 linkage groups and 1,023,625 single nucleotide polymorphisms markers was constructed after genotyping whole-genome sequencing data from 11 broodstock and 363 progeny, belonging to 6 full-sib families. The sex-specific linkage maps spanned a total of 2985.16 cM in females and 2540.47 cM in males with an average inter-marker distance of 0.0030 and 0.0026 cM, respectively. The sex-averaged map spanned a total of 2725.53 cM with an average inter-marker distance of 0.0028 cM. Furthermore, the sex-averaged map was used for improving the contiguity and accuracy of the current pikeperch genome assembly. Based on 723,360 markers, 706 contigs were anchored and oriented into 24 pseudomolecules, covering a total of 896.48 Mb and accounting for 99.47% of the assembled genome size. The overall contiguity of the assembly improved with a scaffold N50 length of 41.06 Mb. Finally, an updated annotation of protein-coding genes and repetitive elements of the enhanced genome assembly is provided at NCBI.

High-Quality SNP Linkage Maps Improved QTL Mapping and Genome Assembly in Populus

With the advances in high-throughput sequencing technologies and the development of new software for extracting single nucleotide polymorphisms (SNPs) across a mapping population, it is possible to construct high-quality genetic maps with thousands of SNPs in outbred forest trees. Two parent-specific linkage maps were constructed with restriction site-associated DNA sequencing data from an F₁ hybrid population derived from Populus deltoides and Populus simonii, and applied in QTL mapping and genome assembly. The female P. deltoides map contained 4018 SNPs, which were divided into 19 linkage groups under a wide range of LOD thresholds from 7 to 55. The male P. simonii map showed similar characteristics, consisting of 2097 SNPs, which also belonged to 19 linkage groups under LOD thresholds of 7 to 29. The SNP order of each linkage group was optimal among different ordering results from several available software. Moreover, the linkage maps allowed the detection of 39 QTLs underlying tree height and 47 for diameter at breast height. In addition, the linkage maps improved the anchoring of 689 contigs of P. simonii to chromosomes. The 2 parental genetic maps of Populus are of high quality, especially in terms of SNP data quality, the SNP order within linkage groups, and the perfect match between the number of linkage groups and the karyotype of Populus, as well as the excellent performances in QTL mapping and genome assembly. Both approaches for extracting and ordering SNPs could be applied to other species for constructing high-quality genetic maps.

High-Density Linkage Maps Based on Genotyping-by-Sequencing (GBS) Confirm a Chromosome-Level Genome Assembly and Reveal Variation in Recombination Rate for the Pacific Oyster Crassostrea gigas

Studies of linkage and linkage mapping have advanced genetic and biological knowledge for over 100 years. In addition to their growing role, today, in mapping phenotypes to genotypes, dense linkage maps can help to validate genome assemblies. Previously, we showed that 40% of scaffolds in the first genome assembly for the Pacific oyster Crassostrea gigas were chimeric, containing single nucleotide polymorphisms (SNPs) mapping to different linkage groups. Here, we merge 14 linkage maps constructed of SNPs generated from genotyping-by-sequencing (GBS) methods with five, previously constructed linkage maps, to create a compendium of nearly 69 thousand SNPs mapped with high confidence. We use this compendium to assess a recently available, chromosome-level assembly of the C. gigas genome, mapping SNPs in 275 of 301 contigs and comparing the ordering of these contigs, by linkage, to their assembly by Hi-C sequencing methods. We find that, while 26% of contigs contain chimeric blocks of SNPs, i.e., adjacent SNPs mapping to different linkage groups than the majority of SNPs in their contig, these apparent misassemblies amount to only 0.08% of the genome sequence. Furthermore, nearly 90% of 275 contigs mapped by linkage and sequencing are assembled identically; inconsistencies between the two assemblies for the remaining 10% of contigs appear to result from insufficient linkage information. Thus, our compilation of linkage maps strongly supports this chromosome-level assembly of the oyster genome. Finally, we use this assembly to estimate, for the first time in a Lophotrochozoan, genome-wide recombination rates and causes of variation in this fundamental process.

Assembly of whole-chromosome pseudomolecules for polyploid plant genomes using outbred mapping populations

Despite advances in sequencing technologies, assembly of complex plant genomes remains elusive due to polyploidy and high repeat content. Here we report PolyGembler for grouping and ordering contigs into pseudomolecules by genetic linkage analysis. Our approach also provides an accurate method with which to detect and fix assembly errors. Using simulated data, we demonstrate that our approach is of high accuracy and outperforms three existing state-of-the-art genetic mapping tools. Particularly, our approach is more robust to the presence of missing genotype data and genotyping errors. We used our method to construct pseudomolecules for allotetraploid lawn grass utilizing PacBio long reads in combination with restriction site-associated DNA sequencing, and for diploid Ipomoea trifida and autotetraploid potato utilizing contigs assembled from Illumina reads in combination with genotype data generated by single-nucleotide polymorphism arrays and genotyping by sequencing, respectively. We resolved 13 assembly errors for a published I. trifida genome assembly and anchored eight unplaced scaffolds in the published potato genome.

Lep-Anchor: automated construction of linkage map anchored haploid genomes

MOTIVATION

Linkage mapping provides a practical way to anchor de novo genome assemblies into chromosomes and to detect chimeric or otherwise erroneous contigs. Such anchoring improves with higher number of markers and individuals, as long as the mapping software can handle all the information. Recent software Lep-MAP3 can robustly construct linkage maps for millions of genotyped markers and on thousands of individuals, providing optimal maps for genome anchoring. For such large datasets, automated and robust genome anchoring tool is especially valuable and can significantly reduce intensive computational and manual work involved.

RESULTS

Here, we present a software Lep-Anchor (LA) to anchor genome assemblies automatically using dense linkage maps. As the main novelty, it takes into account the uncertainty of the linkage map positions caused by low recombination regions, cross type or poor mapping data quality. Furthermore, it can automatically detect and cut chimeric contigs, and use contig-contig, single read or alternative genome assembly alignments as additional information on contig order and orientations and to collapse haplotype contigs. We demonstrate the performance of LA using real data and show that it outperforms ALLMAPS on anchoring completeness and speed. Accuracy-wise LA and ALLMAPS are about equal, but at the expense of lower completeness of ALLMAPS. The software Chromonomer was faster than the other two methods but has major limitations and is lower in accuracy. We also show that with additional information, such as contig-contig and read alignments, the anchoring completeness can be improved by up to 70% without significant loss in accuracy. Based on simulated data, we conclude that the anchoring accuracy can be improved by utilizing information about map position uncertainty. Accuracy is the rate of contigs in correct orientation and completeness is the number contigs with inferred orientation.

AVAILABILITY AND IMPLEMENTATION

Lep-Anchor is available with the source code under GNU general public license from http://sourceforge.net/projects/lep-anchor. All the scripts and code used to produce the reported results are included with Lep-Anchor.

A high-density SNP-based linkage map using genotyping-by-sequencing and its utilization for improved genome assembly of chickpea (Cicer arietinum L.)

Genotyping-by-sequencing (GBS) allows rapid identification of markers for use in development of linkage maps, which expedite efficient breeding programs. In the present study, we have utilized GBS approach to identify and genotype single-nucleotide polymorphism (SNP) markers in an inter-specific RIL population of Cicer arietinum L. X C. reticulatum. A total of 141,639 raw SNPs were identified using the TASSEL-GBS pipeline. After stringent filtering, 8208 candidate SNPs were identified of which ~ 37% were localized in the intragenic regions followed by genic regions (~ 30%) and intergenic regions (~ 27%). We then utilized 6920 stringent selected SNPs from present study and 6714 SNPs and microsatellite markers available from previous studies for construction of linkage map. The resulting high-density linkage map comprising of eight linkage groups contained 13,590 markers which spanned 1299.14 cM of map length with an average marker density of 0.095 cM. Further, the derived linkage map was used to improve the available assembly of desi chickpea genome by anchoring 443 previously unplaced scaffolds onto eight linkage groups. The present efforts have refined anchoring of the desi chickpea genome assembly to 55.57% of the ~ 520 Mb of assembled desi genome. To the best of our knowledge, the linkage map generated in the present study represents one of the most dense linkage map developed for the crop till date. It will serve as a valuable resource for fine mapping and positional cloning of important quantitative trait loci (QTLs) associated with agronomical traits and also for anchoring and ordering of future genome sequence assemblies.

Chromosome-Level Assembly of the Common Lizard (Zootoca vivipara) Genome

Squamate reptiles exhibit high variation in their phenotypic traits and geographical distributions and are therefore fascinating taxa for evolutionary and ecological research. However, genomic resources are very limited for this group of species, consequently inhibiting research efforts. To address this gap, we assembled a high-quality genome of the common lizard, Zootoca vivipara (Lacertidae), using a combination of high coverage Illumina (shotgun and mate-pair) and PacBio sequencing data, coupled with RNAseq data and genetic linkage map generation. The 1.46-Gb genome assembly has a scaffold N50 of 11.52 Mb with N50 contig size of 220.4 kb and only 2.96% gaps. A BUSCO analysis indicates that 97.7% of the single-copy Tetrapoda orthologs were recovered in the assembly. In total, 19,829 gene models were annotated to the genome using a combination of ab initio and homology-based methods. To improve the chromosome-level assembly, we generated a high-density linkage map from wild-caught families and developed a novel analytical pipeline to accommodate multiple paternity and unknown father genotypes. We successfully anchored and oriented almost 90% of the genome on 19 linkage groups. This annotated and oriented chromosome-level reference genome represents a valuable resource to facilitate evolutionary studies in squamate reptiles.

Genomic, Transcriptomic and Epigenomic Tools to Study the Domestication of Plants and Animals: A Field Guide for Beginners

In the last decade, genomics and the related fields of transcriptomics and epigenomics have revolutionized the study of the domestication process in plants and animals, leading to new discoveries and new unresolved questions. Given that some domesticated taxa have been more studied than others, the extent of genomic data can range from vast to nonexistent, depending on the domesticated taxon of interest. This review is meant as a rough guide for students and academics that want to start a domestication research project using modern genomic tools, as well as for researchers already conducting domestication studies that are interested in following a genomic approach and looking for alternate strategies (cheaper or more efficient) and future directions. We summarize the theoretical and technical background needed to carry out domestication genomics, starting from the acquisition of a reference genome and genome assembly, to the sampling design for population genomics, paleogenomics, transcriptomics, epigenomics and experimental validation of domestication-related genes. We also describe some examples of the aforementioned approaches and the relevant discoveries they made to understand the domestication of the studied taxa.

CSA: A high-throughput chromosome-scale assembly pipeline for vertebrate genomes

BACKGROUND

Easy-to-use and fast bioinformatics pipelines for long-read assembly that go beyond the contig level to generate highly continuous chromosome-scale genomes from raw data remain scarce.

RESULT

Chromosome-Scale Assembler (CSA) is a novel computationally highly efficient bioinformatics pipeline that fills this gap. CSA integrates information from scaffolded assemblies (e.g., Hi-C or 10X Genomics) or even from diverged reference genomes into the assembly process. As CSA performs automated assembly of chromosome-sized scaffolds, we benchmark its performance against state-of-the-art reference genomes, i.e., conventionally built in a laborious fashion using multiple separate assembly tools and manual curation. CSA increases the contig lengths using scaffolding, local re-assembly, and gap closing. On certain datasets, initial contig N50 may be increased up to 4.5-fold. For smaller vertebrate genomes, chromosome-scale assemblies can be achieved within 12 h using low-cost, high-end desktop computers. Mammalian genomes can be processed within 16 h on compute-servers. Using diverged reference genomes for fish, birds, and mammals, we demonstrate that CSA calculates chromosome-scale assemblies from long-read data and genome comparisons alone. Even contig-level draft assemblies of diverged genomes are helpful for reconstructing chromosome-scale sequences. CSA is also capable of assembling ultra-long reads.

CONCLUSIONS

CSA can speed up and simplify chromosome-level assembly and significantly lower costs of large-scale family-level vertebrate genome projects.

Development and characterization of genomic resources for a non-model marine teleost, the red snapper (Lutjanus campechanus, Lutjanidae): Construction of a high-density linkage map, anchoring of genome contigs and comparative genomic analysis

The red snapper Lutjanus campechanus is an exploited reef fish of major economic importance in the Gulf of Mexico region. Studies of genome wide genetic variation are needed to understand the structure of wild populations and develop breeding programs for aquaculture but interpretation of these genome scans is limited by the absence of reference genome. In this work, the first draft of a reference genome was developed and characterized for the red snapper. P-454 and Illumina sequencing were conducted to produce paired-end reads that were assembled into reference contigs and scaffolds. The current assembly spans over 770 Mb, representing an estimated 69% of the red snapper genome in 67,254 scaffolds (N50 = 16,803 bp). The genome contigs were applied to map double digest Restriction-Site Associated DNA Tags and characterize Single Nucleotide Polymorphisms (SNPs) in five outbred full-sib families. The identified SNPs and 97 microsatellite loci were used to generate a high-density linkage map that includes 7,420 markers distributed across 24 linkage groups and spans 1,346.64 cM with an average inter–marker distance of 0.18 cM. Sex-specific maps revealed a 1.10:1 female to male map length ratio. A total of 4,422 genome contigs (10.5% of the assembly) were anchored to the map and used in a comparative genomic analysis of the red snapper and two model teleosts. Red snapper showed a high degree of chromosome level syntenic conservation with both medaka and spotted green puffer and a near one to one correspondence between the 24 red snapper linkage groups and corresponding medaka chromosomes was observed. This work established the first draft of a reference genome for a lutjanid fish. The obtained genomic resources will serve as a framework for the interpretation of genome scans during studies of wild populations and captive breeding programs.

Maximising recombination across macadamia populations to generate linkage maps for genome anchoring

The Proteaceae genus Macadamia has a recent history of domestication as a commercial nut crop. We aimed to establish the first sequence-based haploid-correlated reference genetic linkage maps for this primarily outcrossing perennial tree crop, with marker density suitable for genome anchoring. Four first generation populations were used to maximise the segregation patterns available within full-sib, biparental and self-pollinated progeny. This allowed us to combine segregation data from overlapping subsets of >4,000 informative sequence-tagged markers to increase the effective coverage of the karyotype represented by the recombinant crossover events detected. All maps had 14 linkage groups, corresponding to the Macadamia haploid chromosome number, and enabled the anchoring and orientation of sequence scaffolds to construct a pseudo-chromosomal genome assembly for macadamia. Comparison of individual maps indicated a high level of congruence, with minor discrepancies satisfactorily resolved within the integrated maps. The combined set of maps significantly improved marker density and the proportion (70%) of the genome sequence assembly anchored. Overall, increasing our understanding of the genetic landscape and genome for this nut crop represents a substantial advance in macadamia genetics and genomics. The set of maps, large number of sequence-based markers and the reconstructed genome provide a toolkit to underpin future breeding that should help to extend the macadamia industry as well as provide resources for the long term conservation of natural populations in eastern Australia of this unique genus.

Construction of a High-Density Genetic Linkage Map and QTL Mapping for Growth-Related Traits in Takifugu bimaculatus

Takifugu bimaculatus is a euryhaline species, distributed ranging from the southern Yellow Sea to the South China Sea. Their tolerance to a wide range of salinity and temperature, coupled with a desirable firm texture, makes T. bimaculatus a strong candidate for Takifugu aquaculture in subtropics areas. Due to the increasing demand in markets and emerging of the Takifugu aquaculture industry, close attention has been paid to improvement on the T. bimaculatus production. In aquaculture, the great effort has been put into marker-assisted selective breeding, and efficient improvement was realized. However, few genetic resources on T. bimaculatus are provided so far. Aiming at understanding the genetic basis underlying important economic growth traits, facilitating genetic improvement and enriching the genetic resource in T. bimaculatus, we constructed the first genetic linkage map for T. bimaculatus via double digestion restriction-site association DNA sequencing and conducted quantitative traits locus (QTL) mapping for growth-related traits. The map comprised 1976 single nucleotide polymorphism markers distributed on 22 linkage groups (LG), with a total genetic distance of 2039.74 cM. Based on the linkage map, a chromosome-level assembly was constructed whereby we carried out comparative genomics analysis, verifying the high accuracy on contigs ordering of the linkage map. On the other hand, 18 QTLs associated with growth traits were detected on LG6, LG7, LG8, LG10, LG20, and LG21 with phenotypical variance ranging from 15.1 to 56.4%. Candidate genes participating in cartilage development, fat accumulation, and other growth-related regulation activities were identified from these QTLs, including col11a1, foxa2, and thrap3. The linkage map provided a solid foundation for chromosomes assembly and refinement. QTLs reported here unraveled the genomic architecture of some growth traits, which will advance the investigation of aquaculture breeding efforts in T. bimaculatus.

Evolutionary superscaffolding and chromosome anchoring to improve Anopheles genome assemblies

Background

New sequencing technologies have lowered financial barriers to whole genome sequencing, but resulting assemblies are often fragmented and far from ‘finished’. Updating multi-scaffold drafts to chromosome-level status can be achieved through experimental mapping or re-sequencing efforts. Avoiding the costs associated with such approaches, comparative genomic analysis of gene order conservation (synteny) to predict scaffold neighbours (adjacencies) offers a potentially useful complementary method for improving draft assemblies.

Results

We evaluated and employed 3 gene synteny-based methods applied to 21 Anopheles mosquito assemblies to produce consensus sets of scaffold adjacencies. For subsets of the assemblies, we integrated these with additional supporting data to confirm and complement the synteny-based adjacencies: 6 with physical mapping data that anchor scaffolds to chromosome locations, 13 with paired-end RNA sequencing (RNAseq) data, and 3 with new assemblies based on re-scaffolding or long-read data. Our combined analyses produced 20 new superscaffolded assemblies with improved contiguities: 7 for which assignments of non-anchored scaffolds to chromosome arms span more than 75% of the assemblies, and a further 7 with chromosome anchoring including an 88% anchored Anopheles arabiensis assembly and, respectively, 73% and 84% anchored assemblies with comprehensively updated cytogenetic photomaps for Anopheles funestus and Anopheles stephensi.

Conclusions

Experimental data from probe mapping, RNAseq, or long-read technologies, where available, all contribute to successful upgrading of draft assemblies. Our evaluations show that gene synteny-based computational methods represent a valuable alternative or complementary approach. Our improved Anopheles reference assemblies highlight the utility of applying comparative genomics approaches to improve community genomic resources.

Genome of an iconic Australian bird: High-quality assembly and linkage map of the superb fairy-wren (Malurus cyaneus)

The superb fairy-wren, Malurus cyaneus, is one of the most iconic Australian passerine species. This species belongs to an endemic Australasian clade, Meliphagides, which diversified early in the evolution of the oscine passerines. Today, the oscine passerines comprise almost half of all avian species diversity. Despite the rapid increase of available bird genome assemblies, this part of the avian tree has not yet been represented by a high-quality reference. To rectify that, we present the first high-quality genome assembly of a Meliphagides representative: the superb fairy-wren. We combined Illumina shotgun and mate-pair sequences, PacBio long-reads, and a genetic linkage map from an intensively sampled pedigree of a wild population to generate this genome assembly. Of the final assembled 1.07-Gb genome, 975 Mb (90.4%) was anchored onto 25 pseudochromosomes resulting in a final superscaffold N50 of 68.11 Mb. This high-quality bird genome assembly is one of only a handful which is also accompanied by a genetic map and recombination landscape. In comparison to other pedigree-based bird genetic maps, we find that the fairy-wren genetic map more closely resembles those of Taeniopygia guttata and Parus major maps, unlike the Ficedula albicollis map which more closely resembles that of Gallus gallus. Lastly, we also provide a predictive gene and repeat annotation of the genome assembly. This new high-quality, annotated genome assembly will be an invaluable resource not only regarding the superb fairy-wren species and relatives but also broadly across the avian tree by providing a novel reference point for comparative genomic analyses.

Construction of Genetic Linkage Maps in Multiparental Populations

Construction of genetic linkage maps has become a routine step for mapping quantitative trait loci (QTL), particularly in animal and plant breeding populations. Many multiparental populations have recently been produced to increase genetic diversity and QTL mapping resolution. However, few software packages are available for map construction in these populations. In this paper, we build a general framework for the construction of genetic linkage maps from genotypic data in diploid populations, including bi- and multiparental populations, cross-pollinated (CP) populations, and breeding pedigrees. The framework is implemented as an automatic pipeline called magicMap, where the maximum multilocus likelihood approach utilizes genotypic information efficiently. We evaluate magicMap by extensive simulations and eight real datasets: one biparental, one CP, four multiparent advanced generation intercross (MAGIC), and two nested association mapping (NAM) populations, the number of markers ranging from a few hundred to tens of thousands. Not only is magicMap the only software capable of accommodating all of these designs, it is more accurate and robust to missing genotypes and genotyping errors than commonly used packages.

A novel high-density grapevine (Vitis vinifera L.) integrated linkage map using GBS in a half-diallel population

A half-diallel population involving five elite grapevine cultivars was generated and genotyped by GBS, and highly-informative segregation data was used to construct a high-density genetic map for Vitis vinifera L. Grapevine is one of the most relevant fruit crops in the world. Deeper genetic knowledge could assist modern grapevine breeding programs to develop new wine grape varieties able to face climate change effects. To assist in the rapid identification of markers for crop yield components, grape quality traits and adaptation potential, we generated a large Vitis vinifera L. population (N = 624) by crossing five red wine cultivars in a half-diallel scheme, which was subsequently sequenced by an efficient GBS procedure. A high number of fully informative genetic variants was detected using a novel mapping approach capable of reconstructing local haplotypes from adjacent biallelic SNPs, which were subsequently used to construct the densest consensus genetic map available for the cultivated grapevine to date. This 1378.3-cM map integrates 10 bi-parental consensus maps and orders 4437 markers in 3353 unique positions on 19 chromosomes. Markers are well distributed all along the grapevine reference genome, covering up to 98.8% of its genomic sequence. Additionally, a good agreement was observed between genetic and physical orders, adding confidence in the quality of this map. Collectively, our results pave the way for future genetic studies (such as fine QTL mapping) aimed to understand the complex relationship between genotypic and phenotypic variation in the cultivated grapevine. In addition, the method used (which efficiently delivers a high number of fully informative markers) could be of interest to other outbred organisms, notably perennial fruit crops.

An Ultra-Dense Haploid Genetic Map for Evaluating the Highly Fragmented Genome Assembly of Norway Spruce (Picea abies)

Norway spruce (Picea abies (L.) Karst.) is a conifer species of substanital economic and ecological importance. In common with most conifers, the P. abies genome is very large (∼20 Gbp) and contains a high fraction of repetitive DNA. The current P. abies genome assembly (v1.0) covers approximately 60% of the total genome size but is highly fragmented, consisting of >10 million scaffolds. The genome annotation contains 66,632 gene models that are at least partially validated (www.congenie.org), however, the fragmented nature of the assembly means that there is currently little information available on how these genes are physically distributed over the 12 P. abies chromosomes. By creating an ultra-dense genetic linkage map, we anchored and ordered scaffolds into linkage groups, which complements the fine-scale information available in assembly contigs. Our ultra-dense haploid consensus genetic map consists of 21,056 markers derived from 14,336 scaffolds that contain 17,079 gene models (25.6% of the validated gene models) that we have anchored to the 12 linkage groups. We used data from three independent component maps, as well as comparisons with previously published Picea maps to evaluate the accuracy and marker ordering of the linkage groups. We demonstrate that approximately 3.8% of the anchored scaffolds and 1.6% of the gene models covered by the consensus map have likely assembly errors as they contain genetic markers that map to different regions within or between linkage groups. We further evaluate the utility of the genetic map for the conifer research community by using an independent data set of unrelated individuals to assess genome-wide variation in genetic diversity using the genomic regions anchored to linkage groups. The results show that our map is sufficiently dense to enable detailed evolutionary analyses across the P. abies genome.

High-Density Linkage Map and QTLs for Growth in Snapper (Chrysophrys auratus)

Characterizing the genetic variation underlying phenotypic traits is a central objective in biological research. This research has been hampered in the past by the limited genomic resources available for most non-model species. However, recent advances in sequencing technologies and related genotyping methods are rapidly changing this. Here we report the use of genome-wide SNP data from the ecologically and commercially important marine fish species Chrysophrys auratus (snapper) to 1) construct the first linkage map for this species, 2) scan for growth QTL, and 3) search for putative candidate genes in the surrounding QTL regions. The newly constructed linkage map contained ∼11K SNP markers and is one of the densest maps to date in the fish family Sparidae. Comparisons with genome scaffolds of the recently assembled snapper genome indicated that marker placement was mostly consistent between the scaffolds and linkage map (R = 0.7), but that at fine scales (< 5 cM) some precision limitations occurred. Of the 24 linkage groups, which likely reflect the 24 chromosomes of this species, three were found to contain QTL with genome-wide significance for growth-related traits. A scan of 13 candidate growth genes located the growth hormone, myogenin, and parvalbumin genes within 5.3, 9.6, and 25.0 cM of these QTL, respectively. The linkage map and QTL found in this study will advance the investigation of genome structure and aquaculture breeding efforts in this and related species.

Kermit: linkage map guided long read assembly

BACKGROUND

With long reads getting even longer and cheaper, large scale sequencing projects can be accomplished without short reads at an affordable cost. Due to the high error rates and less mature tools, de novo assembly of long reads is still challenging and often results in a large collection of contigs. Dense linkage maps are collections of markers whose location on the genome is approximately known. Therefore they provide long range information that has the potential to greatly aid in de novo assembly. Previously linkage maps have been used to detect misassemblies and to manually order contigs. However, no fully automated tools exist to incorporate linkage maps in assembly but instead large amounts of manual labour is needed to order the contigs into chromosomes.

RESULTS

We formulate the genome assembly problem in the presence of linkage maps and present the first method for guided genome assembly using linkage maps. Our method is based on an additional cleaning step added to the assembly. We show that it can simplify the underlying assembly graph, resulting in more contiguous assemblies and reducing the amount of misassemblies when compared to de novo assembly.

CONCLUSIONS

We present the first method to integrate linkage maps directly into genome assembly. With a modest increase in runtime, our method improves contiguity and correctness of genome assembly.

New Approaches for Genome Assembly and Scaffolding

Affordable, high-throughput DNA sequencing has accelerated the pace of genome assembly over the past decade. Genome assemblies from high-throughput, short-read sequencing, however, are often not as contiguous as the first generation of genome assemblies. Whereas early genome assembly projects were often aided by clone maps or other mapping data, many current assembly projects forego these scaffolding data and only assemble genomes into smaller segments. Recently, new technologies have been invented that allow chromosome-scale assembly at a lower cost and faster speed than traditional methods. Here, we give an overview of the problem of chromosome-scale assembly and traditional methods for tackling this problem. We then review new technologies for chromosome-scale assembly and recent genome projects that used these technologies to create highly contiguous genome assemblies at low cost.

A European Whitefish Linkage Map and Its Implications for Understanding Genome-Wide Synteny Between Salmonids Following Whole Genome Duplication

Genomic datasets continue to increase in number due to the ease of production for a wider selection of species including non-model organisms. For many of these species, especially those with large or polyploid genomes, highly contiguous and well-annotated genomes are still rare due to the complexity and cost involved in their assembly. As a result, a common starting point for genomic work in non-model species is the production of a linkage map. Dense linkage maps facilitate the analysis of genomic data in a variety of ways, from broad scale observations regarding genome structure e.g., chromosome number and type or sex-related structural differences, to fine scale patterns e.g., recombination rate variation and co-localization of differentiated regions. Here we present both sex-averaged and sex-specific linkage maps for Coregonus sp. "Albock", a member of the European whitefish lineage (C. lavaretus spp. complex), containing 5395 single nucleotide polymorphism (SNP) loci across 40 linkage groups to facilitate future investigation into the genomic basis of whitefish adaptation and speciation. The map was produced using restriction-site associated digestion (RAD) sequencing data from two wild-caught parents and 156 F1 offspring. We discuss the differences between our sex-averaged and sex-specific maps and identify genome-wide synteny between C. sp. "Albock" and Atlantic Salmon (Salmo salar), which have diverged following the salmonid-specific whole genome duplication. Our analysis confirms that many patterns of synteny observed between Atlantic Salmon and Oncorhynchus and Salvelinus species are also shared by members of the Coregoninae subfamily. We also show that regions known for their species-specific rediploidization history can pose challenges for synteny identification since these regions have diverged independently in each salmonid species following the salmonid-specific whole genome duplication. The European whitefish map provided here will enable future studies to understand the distribution of loci of interest, e.g., FST outliers, along the whitefish genome as well as assisting with the de novo assembly of a whitefish reference genome.

gmRAD: an integrated SNP calling pipeline for genetic mapping with RADseq across a hybrid population

Restriction site-associated DNA sequencing (RADseq) is a powerful technology that has been extensively applied in population genetics, phylogenetics and genetic mapping. Although many software packages are available for ecological and evolutionary studies, a few effective tools are available for extracting genotype data with RADseq for genetic mapping, a prerequisite for quantitative trait locus mapping, comparative genomics and genome scaffold assembly. Here, we present an integrated pipeline called gmRAD for generating single nucleotide polymorphism (SNP) genotypes from RADseq data, de novo, across a genetic mapping population derived by crossing two parents. As an analytical strategy, the software takes five steps to implement the whole algorithms, including clustering the first (forward) reads of each parent, building two parental references, generating parental SNP catalogs, calling SNP genotypes across all individuals and filtering the genotype data for genetic linkage mapping. All the steps can be completed with a simple command line, but they can be also performed optionally if prerequisite files are available. To validate its application, we also performed a real data analysis with RADseq data from an F1 hybrid population derived by crossing Populus deltoides and Populus simonii. The software gmRAD is freely available at https://github.com/tongchf/gmRAD.

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.

Genetic parameters and quantitative trait loci analysis associated with body size and timing at metamorphosis into glass eels in captive-bred Japanese eels (Anguilla japonica)

The Japanese eel (Anguilla japonica) is among the most important aquaculture fish species in Eastern Asia. The present study aimed to identify the genetic parameters underlying body size and the timing at metamorphosis from leptocephali to glass eels in captive-bred Japanese eels, with the intent to foster sustainable development. Larvae from a partly factorial cross (14 sires × 11 dams) were reared until the point of metamorphosis into glass eels. In these organisms, we observed moderate heritability and mild genetic correlations among traits related to body size (h² = 0.16–0.33) and timing at metamorphosis (h² = 0.36–0.41). In an F₁ full-sib family, quantitative trait loci (QTL) mapping for these traits identified one significant (genome-wide P < 0.05) and five suggestive QTLs (chromosome-wide P < 0.05). These results suggest that in the Japanese eel, metamorphic traits exhibit a polygenic genetic structure comprising many QTLs with small effects. In addition, we updated the genetic linkage map for the Japanese eel and integrated it with our newly constructed de novo genome assembly. The information and tools generated from this study will contribute to the development of freshwater eel genetics and genomics.

High-Resolution Genetic Map and QTL Analysis of Growth-Related Traits of Hevea brasiliensis Cultivated Under Suboptimal Temperature and Humidity Conditions

Rubber tree (Hevea brasiliensis) cultivation is the main source of natural rubber worldwide and has been extended to areas with suboptimal climates and lengthy drought periods; this transition affects growth and latex production. High-density genetic maps with reliable markers support precise mapping of quantitative trait loci (QTL), which can help reveal the complex genome of the species, provide tools to enhance molecular breeding, and shorten the breeding cycle. In this study, QTL mapping of the stem diameter, tree height, and number of whorls was performed for a full-sibling population derived from a GT1 and RRIM701 cross. A total of 225 simple sequence repeats (SSRs) and 186 single-nucleotide polymorphism (SNP) markers were used to construct a base map with 18 linkage groups and to anchor 671 SNPs from genotyping by sequencing (GBS) to produce a very dense linkage map with small intervals between loci. The final map was composed of 1,079 markers, spanned 3,779.7 cM with an average marker density of 3.5 cM, and showed collinearity between markers from previous studies. Significant variation in phenotypic characteristics was found over a 59-month evaluation period with a total of 38 QTLs being identified through a composite interval mapping method. Linkage group 4 showed the greatest number of QTLs (7), with phenotypic explained values varying from 7.67 to 14.07%. Additionally, we estimated segregation patterns, dominance, and additive effects for each QTL. A total of 53 significant effects for stem diameter were observed, and these effects were mostly related to additivity in the GT1 clone. Associating accurate genome assemblies and genetic maps represents a promising strategy for identifying the genetic basis of phenotypic traits in rubber trees. Then, further research can benefit from the QTLs identified herein, providing a better understanding of the key determinant genes associated with growth of Hevea brasiliensis under limiting water conditions.

The State of "Omics" Research for Farmed Penaeids: Advances in Research and Impediments to Industry Utilization

Elucidating the underlying genetic drivers of production traits in agricultural and aquaculture species is critical to efforts to maximize farming efficiency. "Omics" based methods (i.e., transcriptomics, genomics, proteomics, and metabolomics) are increasingly being applied to gain unprecedented insight into the biology of many aquaculture species. While the culture of penaeid shrimp has increased markedly, the industry continues to be impeded in many regards by disease, reproductive dysfunction, and a poor understanding of production traits. Extensive effort has been, and continues to be, applied to develop critical genomic resources for many commercially important penaeids. However, the industry application of these genomic resources, and the translation of the knowledge derived from "omics" studies has not yet been completely realized. Integration between the multiple "omics" resources now available (i.e., genome assemblies, transcriptomes, linkage maps, optical maps, and proteomes) will prove critical to unlocking the full utility of these otherwise independently developed and isolated resources. Furthermore, emerging "omics" based techniques are now available to address longstanding issues with completing keystone genome assemblies (e.g., through long-read sequencing), and can provide cost-effective industrial scale genotyping tools (e.g., through low density SNP chips and genotype-by-sequencing) to undertake advanced selective breeding programs (i.e., genomic selection) and powerful genome-wide association studies. In particular, this review highlights the status, utility and suggested path forward for continued development, and improved use of "omics" resources in penaeid aquaculture.

Highly dense linkage maps from 31 full-sibling families of turbot (Scophthalmus maximus) provide insights into recombination patterns and chromosome rearrangements throughout a newly refined genome assembly

Highly dense linkage maps enable positioning thousands of landmarks useful for anchoring the whole genome and for analysing genome properties. Turbot is the most important cultured flatfish worldwide and breeding programs in the fifth generation of selection are targeted to improve growth rate, obtain disease resistant broodstock and understand sex determination to control sex ratio. Using a Restriction-site Associated DNA approach, we genotyped 18,214 single nucleotide polymorphism in 1,268 turbot individuals from 31 full-sibling families. Individual linkage maps were combined to obtain a male, female and species consensus maps. The turbot consensus map contained 11,845 markers distributed across 22 linkage groups representing a total normalised length of 3,753.9 cM. The turbot genome was anchored to this map, and scaffolds representing 96% of the assembly were ordered and oriented to obtain the expected 22 megascaffolds according to its karyotype. Recombination rate was lower in males, especially around centromeres, and pairwise comparison of 44 individual maps suggested chromosome polymorphism at specific genomic regions. Genome comparison across flatfish provided new evidence on karyotype reorganisations occurring across the evolution of this fish group.

Applications of genotyping by sequencing in aquaculture breeding and genetics

Selective breeding is increasingly recognized as a key component of sustainable production of aquaculture species. The uptake of genomic technology in aquaculture breeding has traditionally lagged behind terrestrial farmed animals. However, the rapid development and application of sequencing technologies has allowed aquaculture to narrow the gap, leading to substantial genomic resources for all major aquaculture species. While high-density single-nucleotide polymorphism (SNP) arrays for some species have been developed recently, direct genotyping by sequencing (GBS) techniques have underpinned many of the advances in aquaculture genetics and breeding to date. In particular, restriction-site associated DNA sequencing (RAD-Seq) and subsequent variations have been extensively applied to generate population-level SNP genotype data. These GBS techniques are not dependent on prior genomic information such as a reference genome assembly for the species of interest. As such, they have been widely utilized by researchers and companies focussing on nonmodel aquaculture species with relatively small research communities. Applications of RAD-Seq techniques have included generation of genetic linkage maps, performing genome-wide association studies, improvements of reference genome assemblies and, more recently, genomic selection for traits of interest to aquaculture like growth, sex determination or disease resistance. In this review, we briefly discuss the history of GBS, the nuances of the various GBS techniques, bioinformatics approaches and application of these techniques to various aquaculture species.