Whole-genome sequencing of wild Siberian musk deer (Moschus moschiferus) provides insights into its genetic features.

BACKGROUND: Siberian musk deer, one of the seven species, is distributed in coniferous forests of Asia. Worldwide, the population size of Siberian musk deer is threatened by severe illegal poaching for commercially valuable musk and meat, habitat losses, and forest fire. At present, this species is categorized as Vulnerable on the IUCN Red List. However, the genetic information of Siberian musk deer is largely unexplored. RESULTS: Here, we produced 3.10 Gb draft assembly of wild Siberian musk deer with a contig N50 of 29,145 bp and a scaffold N50 of 7,955,248 bp. We annotated 19,363 protein-coding genes and estimated 44.44% of the genome to be repetitive. Our phylogenetic analysis reveals that wild Siberian musk deer is closer to Bovidae than to Cervidae. Comparative analyses showed that the genetic features of Siberian musk deer adapted in cold and high-altitude environments. We sequenced two additional genomes of Siberian musk deer constructed demographic history indicated that changes in effective population size corresponded with recent glacial epochs. Finally, we identified several candidate genes that may play a role in the musk secretion based on transcriptome analysis. CONCLUSIONS: Here, we present a high-quality draft genome of wild Siberian musk deer, which will provide a valuable genetic resource for further investigations of this economically important musk deer.

The red deer Cervus elaphus genome CerEla1.0: sequencing, annotating, genes, and chromosomes.

We present here the de novo genome assembly CerEla1.0 for the red deer, Cervus elaphus, an emblematic member of the natural megafauna of the Northern Hemisphere. Humans spread the species in the South. Today, the red deer is also a farm-bred animal and is becoming a model animal in biomedical and population studies. Stag DNA was sequenced at 74× coverage by Illumina technology. The ALLPATHS-LG assembly of the reads resulted in 34.7 × 103 scaffolds, 26.1 × 103 of which were utilized in Cer.Ela1.0. The assembly spans 3.4 Gbp. For building the red deer pseudochromosomes, a pre-established genetic map was used for main anchor points. A nearly complete co-linearity was found between the mapmarker sequences of the deer genetic map and the order and orientation of the orthologous sequences in the syntenic bovine regions. Syntenies were also conserved at the in-scaffold level. The cM distances corresponded to 1.34 Mbp uniformly along the deer genome. Chromosomal rearrangements between deer and cattle were demonstrated. 2.8 × 106 SNPs, 365 × 103 indels and 19368 protein-coding genes were identified in CerEla1.0, along with positions for centromerons. CerEla1.0 demonstrates the utilization of dual references, i.e., when a target genome (here C. elaphus) already has a pre-established genetic map, and is combined with the well-established whole genome sequence of a closely related species (here Bos taurus). Genome-wide association studies (GWAS) that CerEla1.0 (NCBI, MKHE00000000) could serve for are discussed.

De novo characterisation of the greenlip abalone transcriptome (Haliotis laevigata) with a focus on the heat shock protein 70 (HSP70) family.

Abalone (Haliotis) are economically important molluscs for fisheries and aquaculture industries worldwide. Despite this, genomic resources for abalone and molluscs are still limited. Here we present a description and functional annotation of the greenlip abalone (Haliotis laevigata) transcriptome. We present a focused analysis on the heat shock protein 70 (HSP70) family of genes with putative functions affecting temperature stress and immunity. A total of ~38 million paired end Illumina reads were obtained, resulting in a Trinity assembly of 222,172 contigs with minimum length of 200 base pairs and maximum length of 33 kilobases. The 20,702 contigs were annotated with gene descriptions by BLAST. We created a program to maximise the number of functionally annotated genes, and over 10,000 contigs were assigned Gene ontologies (GO terms). By using CateGOrizer, immunity related GO terms for stressors such as heat, hypoxia, oxidative stress and wounding received the highest counts. Twenty-six contigs with homology to the HSP70 family of genes were identified. Ninety-one putative single-nucleotide polymorphisms were observed in the abalone HSP70 contigs. Eleven of these were considered non-synonymous. The annotated transcriptome described in this study will be a useful basis for future work investigating the genetic response of abalone to stress.

Closing a gap in the physical map of the ovine major histocompatibility complex.

A 184 kb gap in an ovine MHC physical map was successfully closed by identification of two overlapping clones (304C7 and 222G18) from a Chinese fine wool merino sheep BAC library. The location and tiling path of the two clones were confirmed by BAC-end sequencing and PCR amplification of loci in overlapping regions. Full-length sequencing of the clones identified 13 novel ovine genes in the gap between loci Notch4 and Btnl2, and eight of them belonging to the Butyrophilin-like (Btn-like or Btnl) gene family. The scattered distribution of the Btnl gene cluster at the gap provided a clue to explain the difficulties previously experienced in closing the gap. Completed BAC contigs of the ovine MHC will facilitate sequencing of the entire ovine leukocyte antigen (OLA) region, providing detailed information for comparative studies of MHC evolution.

A BAC contig map over the proximal approximately 3.3 Mb region of horse chromosome 21.

A total of 207 BAC clones containing 155 loci were isolated and arranged into a map of linearly ordered overlapping clones over the proximal part of horse chromosome 21 (ECA21), which corresponds to the proximal half of the short arm of human chromosome 19 (HSA19p) and part of HSA5. The clones form two contigs - each corresponding to the respective human chromosomes - that are estimated to be separated by a gap of approximately 200 kb. Of the 155 markers present in the two contigs, 141 (33 genes and 108 STS) were generated and mapped in this study. The BACs provide a 4-5x coverage of the region and span an estimated length of approximately 3.3 Mb. The region presently contains one mapped marker per 22 kb on average, which represents a major improvement over the previous resolution of one marker per 380 kb obtained through the generation of a dense RH map for this segment. Dual color fluorescence in situ hybridization on metaphase and interphase chromosomes verified the relative order of some of the BACs and helped to orient them accurately in the contigs. Despite having similar gene order and content, the equine region covered by the contigs appears to be distinctly smaller than the corresponding region in human (3.3 Mb vs. 5.5-6 Mb) because the latter harbors a host of repetitive elements and gene families unique to humans/primates. Considering limited representation of the region in the latest version of the horse whole genome sequence EquCab2, the dense map developed in this study will prove useful for the assembly and annotation of the sequence data on ECA21 and will be instrumental in rapid search and isolation of candidate genes for traits mapped to this region.

A 2.5-Mb contig constructed from Angus, Longhorn and horned Hereford DNA spanning the polled interval on bovine chromosome 1.

The polled locus has been mapped by genetic linkage analysis to the proximal region of bovine chromosome 1. As an intermediate step in our efforts to identify the polled locus and the underlying causative mutation for the polled phenotype, we have constructed a BAC-based physical map of the interval containing the polled locus. Clones containing genes and markers in the critical interval were isolated from the TAMBT (constructed from Angus and Longhorn genomic DNA) and CHORI-240 (constructed from horned Hereford genomic DNA) BAC libraries and ordered based on fingerprinting and the presence or absence of 80 STS markers. A single contig spanning 2.5 Mb was assembled. Comparison of the physical order of STSs to the corresponding region of human chromosome 21 revealed the same order of genes within the polled critical interval. This contig of overlapping BAC clones from horned and polled breeds is a useful resource for SNP discovery and characterization of positional candidate genes.

An expressed sequence tag analysis of the chicken reproductive tract transcriptome.

Analysis of the chicken reproductive tract transcriptome is important in comparative biology for analysis of reproductive tract development and evolution. In addition, molecular analysis of the reproductive tract is important for identification of genes affecting fertility in the poultry industry. We sampled the chicken reproductive tract (ovary, oviduct, and testis) transcriptome, generating 5,328 expressed sequence tags that assembled into 4,518 contigs. We identified 475 contigs with no match in the current expressed sequence tag databases or in GenBank. The novel contigs included 31 with no match to the current assembly of the chicken genome, 119 representing spliced transcripts, and 309 that were unspliced. More detailed molecular characterization of the 428 novel contigs present in the assembly will be important to gene discovery and annotation of the chicken and other vertebrate genomes.

A comparative map of bovine chromosome 19 based on a combination of mapping on a bacterial artificial chromosome scaffold map, a whole genome radiation hybrid panel and the human draft sequence.

We have constructed a medium density physical map of bovine chromosome 19 using a combination of mapping loci on both a bovine bacterial artificial chromosome (BAC) scaffold map and a whole genome radiation hybrid (WGRH) panel. The resulting map contains 70 loci spanning the length of bovine chromosome 19. Three contiguous groups of BACs were identified on the basis of multiple loci mapping to individual BAC clones. Bovine chromosome 19 was found in this study to be comprised almost entirely from regions of human chromosome 17, with a small region putatively assigned to human chromosome 10. Fourteen breakpoints between the bovine and human chromosomes were detected, with a possibility of five more based on ordering of the WGRH map.

Generation of a 5.5-Mb BAC/PAC contig of pig chromosome 6q1.2 and its integration with existing RH, genetic and comparative maps.

We generated a sequence-ready BAC/PAC contig spanning approximately 5.5 Mb on porcine chromosome 6q1.2, which represents a very gene-rich genome region. STS content mapping was used as the main strategy for the assembly of the contig and a total of 6 microsatellite markers, 53 gene-related STS and 116 STS corresponding to BAC and PAC end sequences were analyzed. The contig comprises 316 BAC and PAC clones covering the region between the genes GPI and LIPE. The correct contig assembly was verified by RH-mapping of STS markers and comparative mapping of BAC/PAC end sequences using BLAST searches. The use of microsatellite primer pairs allowed the integration of the physical maps with the genetic map of this region. Comparative mapping of the porcine BAC/PAC contig with respect to the gene-rich region on the human chromosome 19q13.1 map revealed a completely conserved gene order of this segment, however, physical distances differ somewhat between HSA19q13.1 and SSC6q1.2. Three major differences in DNA content between human and pig are found in two large intergenic regions and in one region of a clustered gene family, respectively. While there is a complete conservation of gene order between pig and human, the comparative analysis with respect to the rodent species mouse and rat shows one breakpoint where a genome segment is inverted.

An ordered BAC contig map of the equine major histocompatibility complex.

A physical map of ordered bacterial artificial chromosome (BAC) clones was constructed to determine the genetic organization of the horse major histocompatibility complex. Human, cattle, pig, mouse, and rat MHC gene sequences were compared to identify highly conserved regions which served as source templates for the design of overgo primers. Thirty-five overgo probes were designed from 24 genes and used for hybridization screening of the equine USDA CHORI 241 BAC library. Two hundred thirty-eight BAC clones were assembled into two contigs spanning the horse MHC region. The first contig contains the MHC class II region and was reduced to a minimum tiling path of nine BAC clones that span approximately 800 kb and contain at least 20 genes. A minimum tiling path of a second contig containing the class III/I region is comprised of 14 BAC clones that span approximately 1.6 Mb and contain at least 34 genes. Fluorescence in situ hybridization (FISH) using representative clones from each of the three regions of the MHC localized the contigs onto ECA20q21 and oriented the regions relative to one another and the centromere. Dual-colored FISH revealed that the class I region is proximal to the centromere, the class II region is distal, and the class III region is located between class I and II. These data indicate that the equine MHC is a single gene-dense region similar in structure and organization to the human MHC and is not disrupted as in ruminants and pigs.

Integration of animal linkage and BAC contig maps using overgo hybridization.

The alignment of genome linkage maps, defined primarily by segregation of sequence-tagged site (STS) markers, with BAC contig physical maps and full genome sequences requires high throughput mechanisms to identify BAC clones that contain specific STS. A powerful technique for this purpose is multi-dimensional hybridization of "overgo" probes. The probes are chosen from available STS sequence data by selecting unique probe sequences that have a common melting temperature. We have hybridized sets of 216 overgo probes in subset pools of 36 overgos at a time to filter-spotted chicken BAC clone arrays. A four-dimensional pooling strategy, including one degree of redundancy, has been employed. This requires 24 hybridizations to completely assign BACs for all 216 probes. Results to date are consistent with about a 10% failure rate in overgo probe design and a 15-20% false negative detection rate within a group of 216 markers. Three complete rounds of overgo hybridization, each to sets of about 39,000 BACs (either BAMHI or ECORI partial digest inserts) generated a total of 1853 BAC alignments for 517 mapped chicken genome STS markers. These data are publicly available, and they have been used in the assembly of a first generation BAC contig map of the chicken genome.

Integration of chicken genomic resources to enable whole-genome sequencing.

Different genomic resources in chicken were integrated through the Wageningen chicken BAC library. First, a BAC anchor map was created by screening this library with two sets of markers: microsatellite markers from the consensus linkage map and markers created from BAC end sequencing in chromosome walking experiments. Second, HINdIII digestion fingerprints were created for all BACs of the Wageningen chicken BAC library. Third, cytogenetic positions of BACs were assigned by FISH. These integrated resources will facilitate further chromosome-walking experiments and whole-genome sequencing.

Positional candidate cloning of a QTL in dairy cattle: identification of a missense mutation in the bovine DGAT1 gene with major effect on milk yield and composition.

We recently mapped a quantitative trait locus (QTL) with a major effect on milk composition--particularly fat content--to the centromeric end of bovine chromosome 14. We subsequently exploited linkage disequilibrium to refine the map position of this QTL to a 3-cM chromosome interval bounded by microsatellite markers BULGE13 and BULGE09. We herein report the positional candidate cloning of this QTL, involving (1) the construction of a BAC contig spanning the corresponding marker interval, (2) the demonstration that a very strong candidate gene, acylCoA:diacylglycerol acyltransferase (DGAT1), maps to that contig, and (3) the identification of a nonconservative K232A substitution in the DGAT1 gene with a major effect on milk fat content and other milk characteristics.

Sequence-ready BAC contig, physical, and transcriptional map of a 2-Mb region overlapping the mouse chromosome 6 host-resistance locus Cmv1.

The host-resistance locus Cmv1 controls viral replication of mouse cytomegalovirus (MCMV) in the spleen of infected mice. Cmv1 maps on distal chromosome 6, very tightly linked to the Ly49 gene family within a 0.35-cM interval defined proximally by Cd94/Nkg2d and distally by D6Mit13/D6Mit111/D6Mit219/Prp/Kap. To facilitate the cloning of the gene, we have created a high-resolution physical map of the Cmv1 genetic interval that is based on long-range restriction mapping by pulsed-field gel electrophoresis, fluorescence in situ hybridization analysis of interphase nuclei, and the assembly of a cloned contig. A contig of BAC and YAC clones was assembled using probes derived from the minimal genetic interval. Individual clones from the region were validated by (1) restriction digest fingerprinting, (2) STS content mapping, (3) Southern hybridizations, and (4) sequencing and mapping of clone ends. This contig contains 25 YACs anchored by 71 STSs and 73 BACs anchored by 40 STSs. We also report the cloning of 31 new STSs and 18 new polymorphic markers. A minimum tiling path was defined that consists of either 4 YACs or 13 BACs covering 1.82 Mb between D6Ott8, the closest proximal marker, and D6Ott115, the closest distal marker. Gene distribution in the region includes 14 Ly49 genes as well as 3 new additional transcripts. This high-resolution, sequence-ready BAC contig provides a backbone for the identification of Cmv1 and its relationship with genes involved in innate immunity.