Use of bovine EST data and human genomic sequences to map 100 gene-specific bovine markers

A system to use bovine EST data in conjunction with human genomic sequence to improve the bovine linkage map over the entire genome or on specific chromosomes was evaluated. Bovine EST sequence was used to provide primer sequences corresponding to bovine genes, while human genomic sequence directed primer design to flank introns and produce amplicons of appropriate size for efficient direct sequencing. The sequence tagged sites (STS) produced in this way from the four sires of the MARC reference families were examined for single nucleotide polymorphisms (SNPs) that could be used to map the corresponding genes. With this approach, along with a primer/extension mass spectrometry SNP genotyping assay, 100 ESTs were placed on the bovine genetic linkage map. The first 70 were chosen at random from bovine EST–human genomic comparisons. An additional 30 ESTs were successfully mapped to bovine Chromosome 19 (BTA19), and comparison of the resulting BTA19 map to the position of the corresponding human orthologs on the HSA17 draft sequences revealed differences in the spacing and order of genes. Over 80% of successful amplicons contained SNPs, indicating that this is an efficient approach to generating EST-associated genetic markers. We have demonstrated the feasibility of constructing a linkage map based on SNPs associated with ESTs and the plausibility of utilizing EST, comparative mapping information, and human sequence data to target regions of the bovine genome for SNP marker development.     

Single nucleotide polymorphism (SNP) discovery and linkage mapping of bovine cytokine genes

Polymorphic markers at bovine gene loci facilitate the integration of cattle genetic maps with those of humans and mice. To this end, 31 single nucleotide polymorphism (SNP) markers were developed for seven bovine chemokine genes. Loci were amplified from bovine genomic DNA by the polymerase chain reaction, and candidate amplicons were sequenced to determine their identity. Amplified loci from 24 founding parents and select progeny from a beef cattle reference population were sequenced and analyzed for SNPs. SNP haplotype alleles were determined by examining segregation patterns and used to establish the locus position on the bovine linkage map. Loci for growth-related proteins (GRO3, GRO1, and GROX) were clustered with the related CXC chemokine genes, interleukin (IL) 8, and epithelial cell inflammatory protein 1, at 84 cM from the centromeric end of the bovine chromosome (BTA) 6 linkage group. Bovine loci for a cluster of IL8 receptors, a stromal cell-derived factor 1, interferon-γ, and tumor necrosis factor-α were mapped at 90, 55, 59, and 34 cM, respectively, from the centromeric ends of the BTA 2, 28, 5, and 23 linkage groups. The positions of these bovine loci were compared with those of orthologous loci on the human map to refine the boundaries of conserved synteny. These seven loci provide examples of SNP development in which the efficiency was largely dependent on the availability of bovine genomic or cDNA sequence. The polymorphic nature of these SNP haplotype markers suggests that they will be useful for mapping complex traits in cattle, such as resistance to infectious disease. Received: 30 April 1999 / Accepted: 12 July 1999     

A first-generation metric linkage disequilibrium map of bovine chromosome 6

We constructed a metric linkage disequilibrium (LD) map of bovine chromosome 6 (BTA6) on the basis of data from 220 SNPs genotyped on 433 Australian dairy bulls. This metric LD map has distances in LD units (LDUs) that are analogous to centimorgans in linkage maps. The LD map of BTA6 has a total length of 8.9 LDUs. Within the LD map, regions of high LD (represented as blocks) and regions of low LD (steps) are observed, when plotted against the integrated map in kilobases. At the most stringent block definition, namely a set of loci with zero LDU increase over the span of these markers, BTA6 comprises 40 blocks, accounting for 41% of the chromosome. At a slightly lower stringency of block definition (a set of loci covering a maximum of 0.2 LDUs on the LD map), up to 81% of BTA6 is spanned by 46 blocks and with 13 steps that are likely to reflect recombination hot spots. The mean swept radius (the distance over which LD is likely to be useful for mapping) is 13.3 Mb, confirming extensive LD in Holstein-Friesian dairy cattle, which makes such populations ideal for whole-genome association studies.     

A genome‐wide linkage map for the house sparrow (Passer domesticus) provides insights into the evolutionary history of the avian genome

The house sparrow is an important model species for studying physiological, ecological and evolutionary processes in wild populations. Here, we present a medium density, genome wide linkage map for house sparrow (Passer domesticus) that has aided the assembly of the house sparrow reference genome, and that will provide an important resource for ongoing mapping of genes controlling important traits in the ecology and evolution of this species. Using a custom house sparrow 10 K iSelect Illumina SNP chip we have assigned 6,498 SNPs to 29 autosomal linkage groups, based on a mean of 430 informative meioses per SNP. The map was constructed by combining the information from linkage with that of the physical position of SNPs within scaffold sequences in an iterative process. Averaged between the sexes; the linkage map had a total length of 2,004 cM, with a longer map for females (2,240 cM) than males (1,801 cM). Additionally, recombination rates also varied along the chromosomes. Comparison of the linkage map to the reference genomes of zebra finch, collared flycatcher and chicken, showed a chromosome fusion of the two avian chromosomes 8 and 4A in house sparrow. Lastly, information from the linkage map was utilized to conduct analysis of linkage disequilibrium (LD) in eight populations with different effective population sizes (N_e) in order to quantify the background level LD. Together, these results aid the design of future association studies, facilitate the development of new genomic tools and support the body of research that describes the evolution of the avian genome.     

SNP Identification from RNA Sequencing and Linkage Map Construction of Rubber Tree for Anchoring the Draft Genome

Hevea brasiliensis, or rubber tree, is an important crop species that accounts for the majority of natural latex production. The rubber tree nuclear genome consists of 18 chromosomes and is roughly 2.15 Gb. The current rubber tree reference genome assembly consists of 1,150,326 scaffolds ranging from 200 to 531,465 bp and totalling 1.1 Gb. Only 143 scaffolds, totalling 7.6 Mb, have been placed into linkage groups. We have performed RNA-seq on 6 varieties of rubber tree to identify SNPs and InDels and used this information to perform target sequence enrichment and high throughput sequencing to genotype a set of SNPs in 149 rubber tree offspring from a cross between RRIM 600 and RRII 105 rubber tree varieties. We used this information to generate a linkage map allowing for the anchoring of 24,424 contigs from 3,009 scaffolds, totalling 115 Mb or 10.4% of the published sequence, into 18 linkage groups. Each linkage group contains between 319 and 1367 SNPs, or 60 to 194 non-redundant marker positions, and ranges from 156 to 336 cM in length. This linkage map includes 20,143 of the 69,300 predicted genes from rubber tree and will be useful for mapping studies and improving the reference genome assembly.     

A first-generation genetic linkage map of the European flat oyster Ostrea edulis (L.) based on AFLP and microsatellite markers

This study presents the first genetic linkage map for the European flat oyster Ostrea edulis . Two hundred and forty-six AFLP and 20 microsatellite markers were genotyped in a three-generation pedigree comprising two grandparents, two parents and 92 progeny. Chi-square goodness-of-fit tests revealed high segregation distortion, which was significant for 32.8% of markers. Sixteen microsatellites and 235 AFLPs (170 type 1:1 AFLPs and 65 type 3:1 AFLPs) were used to build sex-specific linkage maps using crimap software. The first parental map (P₁) consisted of 104 markers grouped in nine linkage groups, and spanned 471.2 cM with an average spacing of 4.86 cM. The second parental map (P₂) consisted of 117 markers grouped in 10 linkage groups (which equals the haploid chromosome number), and covered 450.0 cM with an average spacing of 4.21 cM. The estimated coverage of the genome was 82.4% for the P₁ map and 84.2% for the P₂ map. Eight linkage groups that were probably homologous between the two parents contained the same microsatellites and 3:1 AFLPs (segregating through both parents). Distorted markers were not randomly distributed across the genome and tended to cluster in a few linkage groups. Sex-specific differences in recombination rates were evident. This first-generation genetic linkage map for O. edulis represents a major step towards the mapping of QTL such as resistance to bonamiasis, a parasitosis that has drastically decreased populations of flat oysters since the 1960s.     

FISH mapping in cattle (Bos taurus L.) is not yet out of fashion

Physical mapping of genes by fluorescence in situ hybridization (FISH) seems to be out of fashion in species whose assembled genome sequences are available. However, in this work we evidence the existence of errors in gene location in the Btau_4.0 assembly. We show thatDFNA5 andCHCHD6 genes are located on BTA4 and BTA22, respectively, instead of BTA10 and BTA3, as displayed by Btau_4.0. This report emphasizes the need to verify the data on physical localization of genes in the cattle genome (at least by taking into account comparative data reported in available papers) and the need to improve the cattle genome assembly. Our results indicate that FISH mapping in cattle is still useful.     

Characterization of a 320-kb region containing the HEXA gene on bovine chromosome 10 and analysis of its association with BSE susceptibility

Bovine spongiform encephalopathy (BSE) belongs to a group of neurodegenerative diseases known as transmissible prion diseases. Recently, variants in the promoter region of the prion protein ( PRNP ) gene have been shown to have a considerable effect on the susceptibility to BSE. However, a previous genome scan revealed other putative BSE-susceptibility loci. Here, we analysed such a region on BTA10, which contains the functional candidate gene HEXA . Three hundred and twenty kilobases that, besides HEXA , also contain ARIH1 , BRUNOL6 and PARP6 were characterized and screened for polymorphisms. Genotyping of 38 SNPs in Holstein–Friesian animals from the UK (350 diseased and 270 controls) revealed two intronic SNPs that were associated with BSE incidence, with experiment-wise P -values of 3.5 × 10⁻³ and 7.7 × 10⁻³ respectively. Both SNPs were in strong linkage disequilibrium and the rare alleles had a protective effect. These alleles were contained in a haplotype dubbed 'UK-protective' that was significantly overrepresented in the controls with a permuted P -value of 2 × 10⁻³. An association study in German Holstein animals (73 diseased and 627 controls) revealed an opposite effect of the 'UK-protective' haplotype in this population, i.e. it was overrepresented in the diseased animals, although not significant after correction for multiple testing. These findings indicate a causal variant for BSE susceptibility on BTA10 in linkage disequilibrium with the markers studied. Candidate gene analyses of the surrounding region and additional association studies will help to clarify the origin of the protective effects and to identify causal variants for BSE susceptibility on BTA10.     

Consensus and comprehensive linkage maps of bovine chromosome 24

This study describes development of a consensus genetic linkage map of bovine chromosome 24 (BTA24). Eight participating laboratories contributed data for 58 unique markers including a total of 25 409 meioses. Eighteen markers, which were typed in more than one reference population, were used as potential anchors to generate a consensus framework map. The framework map contained 16 loci ordered with odds greater than 1000:1 and spanned 79.3 cM. Remaining markers were included in a comprehensive map relative to these anchors. The resulting BTA24 comprehensive map was 98.3 cM in length. Average marker intervals were 6.1 and 2.5 cM for framework and comprehensive maps, respectively. Marker order was generally consistent with previously reported BTA24 linkage maps. Only one discrepancy was found when comparing the comprehensive map with the published USDA-MARC linkage map. Integration of genetic information from different maps provides a high-resolution BTA24 linkage map.     

A Multi-Population Consensus Genetic Map Reveals Inconsistent Marker Order among Maps Likely Attributed to Structural Variations in the Apple Genome

Genetic maps serve as frameworks for determining the genetic architecture of quantitative traits, assessing structure of a genome, as well as aid in pursuing association mapping and comparative genetic studies. In this study, a dense genetic map was constructed using a high-throughput 1,536 EST-derived SNP GoldenGate genotyping platform and a global consensus map established by combining the new genetic map with four existing reliable genetic maps of apple. The consensus map identified markers with both major and minor conflicts in positioning across all five maps. These major inconsistencies among marker positions were attributed either to structural variations within the apple genome, or among mapping populations, or genotyping technical errors. These also highlighted problems in assembly and anchorage of the reference draft apple genome sequence in regions with known segmental duplications. Markers common across all five apple genetic maps resulted in successful positioning of 2875 markers, consisting of 2033 SNPs and 843 SSRs as well as other specific markers, on the global consensus map. These markers were distributed across all 17 linkage groups, with an average of 169±33 marker per linkage group and with an average distance of 0.70±0.14 cM between markers. The total length of the consensus map was 1991.38 cM with an average length of 117.14±24.43 cM per linkage group. A total of 569 SNPs were mapped onto the genetic map, consisting of 140 recombinant individuals, from our recently developed apple Oligonucleotide pool assays (OPA). The new functional SNPs, along with the dense consensus genetic map, will be useful for high resolution QTL mapping of important traits in apple and for pursuing comparative genetic studies in Rosaceae.     

Construction of a high-density SNP genetic map in flue-cured tobacco based on SLAF-seq

Tobacco (Nicotiana tabacum L., 2n = 48) is an important agronomic crop and model plant. Flue-cured tobacco is the most important type and accounts for approximately 80 % of tobacco production worldwide. The low genetic diversity of flue-cured tobacco impedes the construction of a high-density genetic linkage map using simple sequence repeat (SSR) markers and warrants the exploitation of single nucleotide polymorphic (SNP) markers from genomic regions. In this article, initially using specific locus-amplified fragment sequencing, we discovered 10,891 SNPs that were subsequently used as molecular markers for genetic map construction. Combined with SSR markers, a final high-density genetic map was generated containing 4215 SNPs and 194 SSRs distributed on 24 linkage groups (LGs). The genetic map was 2662.43 cM in length, with an average distance of 0.60 cM between adjacent markers. Furthermore, by mapping the SNP markers to the ancestral genomes of Nicotiana tomentosiformis and Nicotiana sylvestris, a large number of genome rearrangements were identified as occurring after the polyploidization event. Finally, using this novel integrated map and mapping population, two major quantitative trait loci (QTLs) were identified for flue-curing and mapped to the LG6 of tobacco. This is the first report of SNP markers and a SNP-based linkage map being developed in tobacco. The high-density genetic map and QTLs related to tobacco curing will support gene/QTL fine mapping, genome sequence assembly and molecular breeding in tobacco.     

Recombination rates across porcine autosomes inferred from high-density linkage maps

Studies of the variation in recombination rate across the genome provide a better understanding of evolutionary genomics and are also an important step towards mapping and dissecting complex traits in domestic animals. With the recent completion of the porcine genome sequence and the availability of a high‐density porcine single nucleotide polymorphism (SNP) array, it is now possible to construct a high‐density porcine linkage map and estimate recombination rate across the genome. A total of 416 animals were genotyped with the Porcine SNP60BeadChip, and high‐density chromosome linkage maps were constructed using CRI‐MAP, assuming the physical order of the Sscrofa10 assembly. The total linkage map length was 2018.79 cM, using 658 meioses and 14 503 SNPs. The estimated average recombination rate across the porcine autosomes was 0.86 cM/Mb. However, a large variation in recombination rate was observed among chromosomes. The estimated average recombination rates (cM/Mb) per chromosome ranged from 0.48 in SSC1 to 1.48 in SSC10, displaying a significant negative correlation with the chromosome sizes. In addition, the analysis of the variation in the recombination rates taking 1‐Mb sliding windows has allowed us to demonstrate the variation in recombination rates within chromosomes. In general, a larger recombination rate was observed in the extremes than in the centre of the chromosome. Finally, the ratio between female and male recombination rates was also inferred, obtaining a value of 1.38, with the heterogametic sex having the least recombination.     

Comparison of the power between microsatellite and single-nucleotide polymorphism markers for linkage and linkage disequilibrium mapping of an electrophysiological phenotype

We performed linkage and linkage disequilibrium (LD) mapping analyses to compare the power between microsatellite and single nucleotide polymorphism (SNP) markers. Chromosome-wide analyses were performed for a quantitative electrophysiological phenotype, ttth1, on chromosome 7. Multipoint analysis of microsatellite markers using the variance component (VC) method showed the highest LOD score of 4.20 at 162 cM, near D7S509 (163.7 cM). Two-point analysis of SNPs using the VC method yielded the highest LOD score of 3.98 in the Illumina SNP data and 3.45 in the Affymetrix SNP data around 152-153 cM. In family-based single SNP and SNP haplotype LD analysis, we identified seven SNPs associated with ttth1. We searched for any potential candidate genes in the location of the seven SNPs. The SNPs rs1476640 and rs768055 are located in the FLJ40852 gene (a hypothetical protein), and SNP rs1859646 is located in the TAS2R5 gene (a taste receptor). The other four SNPs are not located in any known or annotated genes. We found the high density SNP scan to be superior to microsatellites because it is effective in downstream fine mapping due to a better defined linkage region. Our study proves the utility of high density SNP in genome-wide mapping studies.     

Comparative mapping of bovine chromosome 27 with human chromosome 8 near a dairy form QTL in cattle

In the absence of a complete and annotated bovine genome sequence, detailed human-bovine comparative maps are one of the most effective tools for identification of positional candidate genes contributing to quantitative trait loci (QTL) in cattle. In the present study, eight genes from human chromosome 8 were selected for mapping in cattle to improve breakpoint resolution and confirm gene order on the comparative map near the 40 cM region of the BTA27 linkage map where a QTL affecting dairy form had previously been identified. The resulting map identified ADRB3 as a positional candidate gene for the QTL contributing to the dairy form trait based on its estimated position between 40 and 45 cM on the linkage map. It is also a functional candidate gene due to its role in fat metabolism, and polymorphisms in the ADRB3 gene associated with obesity and metabolic disease in humans, as well as, carcass fat in sheep. Further studies are underway to investigate the existence of polymorphisms in the bovine ADRB3 gene and their association with traits related to fat deposition in cattle.     

A genetic linkage map and comparative genome analysis of common carp (Cyprinus carpio L.) using microsatellites and SNPs

A genetic linkage map is a powerful research tool for mapping traits of interest and is essential to understanding genome evolution. The aim of this study is to provide an expanded genetic linkage map of common carp to effectively carry out quantitative trait loci analysis and conduct comparative mapping analysis between lineages. Here, we constructed a genetic linkage map of common carp (Cyprinus carpio L.) using microsatellite and single-nucleotide polymorphism (SNP) markers in a 159 sibling family. A total of 246 microsatellites and 306 SNP polymorphic markers were genotyped in this family. Linkage analysis using JoinMap 4.0 organized 427 markers (186 microsatellites and 241 SNPs) to 50 linkage groups, ranging in size from 1.4 to 130.1 cM. Each group contained 2-30 markers. The linkage map covered a genetic distance of 2,039.2 cM and the average interval for markers within the linkage groups was approximately 6.4 cM. In addition, comparative genome analysis within five model teleost fish revealed a high percentage (74.7%) of conserved loci corresponding to zebrafish chromosomes. In most cases, each zebrafish chromosome comprised two common carp linkage groups. The comparative analysis also revealed independent chromosome rearrangements in common carp and zebrafish. The linkage map will be of great assistance in mapping genes of interest and serve as a reference to approach comparative mapping and enable further insights into the comprehensive investigations of genome evolution of common carp.     

Confirmation and refinement of a QTL on BTA5 affecting milk production traits in the Fleckvieh dual purpose cattle breed

We analysed a QTL affecting milk yield (MY), milk protein yield (PY) and milk fat yield (FY) in the dual purpose cattle breed Fleckvieh on BTA5. Twenty-six microsatellite markers covering 135 cM were selected to analyse nine half-sib families containing 605 sons in a granddaughter design. We thereby assigned two new markers to the public linkage map using the CRI-MAP program. Phenotypic records were daughter yield deviations (DYD) originating from the routinely performed genetic evaluations of breeding animals. To determine the position of the QTL, three different approaches were applied: interval mapping (IM), linkage analysis by variance component analysis (LAVC), and combined linkage disequilibrium (LD) and linkage (LDL) analysis. All three methods mapped the QTL in the same marker interval ( BM2830-ETH152 ) with the greatest test-statistic value at 118, 119.33 and 119.33 cM respectively. The positive QTL allele simultaneously increases DYD in the first lactation by 272 kg milk, 7.1 kg milk protein and 7.0 kg milk fat. Although the mapping accuracy and the significance of a QTL effect increased from IM over LAVC to LDL, the confidence interval was large (13, 20 and 24 cM for FY, MY and PY respectively) for the positional cloning of the causal gene. The estimated averages of pair wise marker LD with a distance <5 cM were low (0.107) and reflect the large effective population size of the Fleckvieh subpopulation analysed. This low level of LD suggests a need for increase in marker density in following fine mapping steps.     

A genetic and physical map of bovine chromosome 3

This paper reports a map of nine polymorphic microsatellite markers previously assigned to bovine chromosome 3 (BTA3) by somatic cell genetics. The linkage group covers 101 cM on the chromosome with an average intermarker distance of 13-9 cM. One marker (INRA200) was isolated from a peak of flow sorted chromosomes 2 and 3. Another marker (INRA197) was derived from a cosmid. The localization of the cosmid by in situ hybridization enabled the orientation of the linkage group on BTA3. Markers were relatively evenly spaced and consequently can be used to complement other mapping data about this chromosome. This establishes a framework of polymorphic markers that can be used to search for quantitative trait loci (QTL).     

Levels of linkage disequilibrium in a wild bird population

Population-based mapping approaches are attractive for tracing the genetic background to phenotypic traits in wild species, given that it is often difficult to gather extensive and well-defined pedigrees needed for quantitative trait locus analysis. However, the feasibility of association or hitch-hiking mapping is dependent on the degree of linkage disequilibrium (LD) in the population, on which there is yet limited information for wild species. Here we use single nucleotide polymorphism (SNP) markers from 23 genes in a recently established linkage map of the Z chromosome of the collared flycatcher, to study the extent of LD in a natural bird population. In most but not all cases we find SNPs within the same intron (less than 500 bp) to be in perfect LD. However, LD then decays to background level at a distance 1cM or 400-500 kb. Although LD seems more extensive than in other species, if the observed pattern is representative for other regions of the genome and turns out to be a general feature of natural bird populations, dense marker maps might be needed for genome scans aimed at identifying association between marker and trait loci.     

Construction of bovine whole-genome radiation hybrid and linkage maps using high-throughput genotyping

High-density whole-genome maps are essential for ordering genes or markers and aid in the assembly of genome sequence. To increase the density of markers on the bovine radiation hybrid map, and hence contribute to the assembly of the bovine genome sequence, an Illumina BeadStation was used to simultaneously type large numbers of markers on the Roslin-Cambridge 3000 rad bovine-hamster whole-genome radiation hybrid panel (WGRH3000). In five multiplex reactions, 6738 sequence tagged site (STS) markers were successfully typed on the WGRH3000 panel DNA. These STSs harboured SNPs that were developed as a result of the bovine genome sequencing initiative. Typically, the most time consuming and expensive part of creating high-density radiation hybrid (RH) maps is genotyping the markers on the RH panel with conventional approaches. Using the method described in this article, we have developed a high-density whole-genome RH map with 4690 loci and a linkage map with 2701 loci, with direct comparison to the bovine whole-genome sequence assembly (Btau_2.0) in a fraction of the time it would have taken with conventional typing and genotyping methods.     

Whole genome scan to detect quantitative trait loci for bovine milk protein composition

The objective of this study was to perform a whole genome scan to detect quantitative trait loci (QTL) for milk protein composition in 849 Holstein–Friesian cows originating from seven sires. One morning milk sample was analysed for the major milk proteins using capillary zone electrophoresis. A genetic map was constructed with 1341 single nucleotide polymorphisms, covering 2829 centimorgans (cM) and 95% of the cattle genome. The chromosomal regions most significantly related to milk protein composition ( P _genome < 0.05) were found on Bos taurus autosomes (BTA) 6, 11 and 14. The QTL on BTA6 was found at about 80 cM, and affected α_S1-casein, α_S2-casein, β-casein and κ-casein. The QTL on BTA11 was found at 124 cM, and affected β-lactoglobulin, and the QTL on BTA14 was found at 0 cM, and affected protein percentage. The proportion of phenotypic variance explained by the QTL was 3.6% for β-casein and 7.9% for κ-casein on BTA6, 28.3% for β-lactoglobulin on BTA11, and 8.6% for protein percentage on BTA14. The QTL affecting α_S2-casein on BTA6 and 17 showed a significant interaction. We investigated the extent to which the detected QTL affecting milk protein composition could be explained by known polymorphisms in β-casein , κ -casein , β-lactoglobulin and DGAT1 genes. Correction for these polymorphisms decreased the proportion of phenotypic variance explained by the QTL previously found on BTA6, 11 and 14. Thus, several significant QTL affecting milk protein composition were found, of which some QTL could partially be explained by polymorphisms in milk protein genes.