High-resolution comprehensive radiation hybrid maps of the porcine chromosomes 2p and 9p compared with the human chromosome 11

High-resolution autosomal radiation hybrid maps of the pig genome and their contribution to the genome sequence assembly

BACKGROUND

The release of the porcine genome sequence offers great perspectives for Pig genetics and genomics, and more generally will contribute to the understanding of mammalian genome biology and evolution. The process of producing a complete genome sequence of high quality, while facilitated by high-throughput sequencing technologies, remains a difficult task. The porcine genome was sequenced using a combination of a hierarchical shotgun strategy and data generated with whole genome shotgun. In addition to the BAC contig map used for the clone-by-clone approach, genomic mapping resources for the pig include two radiation hybrid (RH) panels at two different resolutions. These two panels have been used extensively for the physical mapping of pig genes and markers prior to the availability of the pig genome sequence.

RESULTS

In order to contribute to the assembly of the pig genome, we genotyped the two radiation hybrid (RH) panels with a SNP array (the Illumina porcineSNP60 array) and produced high density physical RH maps for each pig autosome. We first present the methods developed to obtain high density RH maps with 38,379 SNPs from the SNP array genotyping. We then show how they were useful to identify problems in a draft of the pig genome assembly, and how the RH maps enabled the problems to be corrected in the porcine genome sequence. Finally, we used the RH maps to predict the position of 2,703 SNPs and 1,328 scaffolds currently unplaced on the porcine genome assembly.

CONCLUSIONS

A complete process, from genotyping of a high density SNP array on RH panels, to the construction of genome-wide high density RH maps, and finally their exploitation for validating and improving a genome assembly is presented here. The study includes the cross-validation of RH based findings with independent information from genetic data and comparative mapping with the Human genome. Several additional resources are also provided, in particular the predicted genomic location of currently unplaced SNPs and associated scaffolds summing up to a total of 72 megabases, that can be useful for the exploitation of the pig genome assembly.

A high-resolution comparative map of porcine chromosome 4 (SSC4)

We used the IMNpRH2(12,000-rad) RH and IMpRH(7,000-rad) panels to integrate 2019 transcriptome (RNA-seq)-generated contigs with markers from the porcine genetic and radiation hybrid (RH) maps and bacterial artificial chromosome finger-printed contigs, into 1) parallel framework maps (LOD ≥ 10) on both panels for swine chromosome (SSC) 4, and 2) a high-resolution comparative map of SSC4, thus and human chromosomes (HSA) 1 and 8. A total of 573 loci were anchored and ordered on SSC4 closing gaps identified in the porcine sequence assembly Sscrofa9. Alignment of the SSC4 RH with the genetic map identified five microsatellites incorrectly mapped around the centromeric region in the genetic map. Further alignment of the RH and comparative maps with the genome sequence identified four additional regions of discrepancy that are also suggestive of errors in assembly, three of which were resolved through conserved synteny with blocks on HSA1 and HSA8.

An integrated RH map of porcine chromosome 10

BACKGROUND

Whole genome radiation hybrid (WG-RH) maps serve as "scaffolds" to significantly improve the orientation of small bacterial artificial chromosome (BAC) contigs, order genes within the contigs and assist assembly of a sequence-ready map for virtually any species. Here, we report the construction of a porcine: human comparative map for pig (Sus scrofa) chromosome 10 (SSC10) using the IMNpRH2(12,000-rad) porcine WG-RH panel, integrated with the IMpRH(7000-rad) WG-RH, genetic and BAC fingerprinted contig (FPC) maps.

RESULTS

Map vectors from the IMNpRH2(12,000-rad) and IMpRH(7,000-rad) panels were merged to construct parallel framework (FW) maps, within which FW markers common to both panels have an identical order. This strategy reduced map discrepancies between the two panels and significantly improved map accuracy. A total of 216 markers, including 50 microsatellites (MSs), 97 genes and ESTs, and 69 BAC end sequences (BESs), were ordered within two linkage groups at two point (2 pt) LOD score of 8. One linkage group covers SSC10p with accumulated map distances of 738.2 cR(7,000) and 1814.5 cR(12,000), respectively. The second group covers SSC10q at map distances of 1336.9 cR(7,000) and 3353.6 cR(12,000), yielding an overall average map resolution of 16.4 kb/cR(12,000) or 393.5 kb per marker on SSC10. This represents an approximately 2.5-fold increase in map resolution over the IMpRH(7,000-rad) panel. Based on 127 porcine markers that have homologous sequences in the human genome, a detailed comparative map between SSC10 and human (Homo sapiens) chromosome (HSA) 1, 9 and 10 was built.

CONCLUSION

This initial comparative RH map of SSC10 refines the syntenic regions between SSC10 and HSA1, 9 and 10. It integrates the IMNpRH2(12,000-rad) and IMpRH(7,000-rad), genetic and BAC FPC maps and provides a scaffold to close potential gaps between contigs prior to genome sequencing and assembly. This map is also useful in fine mapping of QTLs on SSC10.