An integrated genetic and physical map of the bovine X chromosome

Dissecting the genomic regions of selection on the X chromosome in different cattle breeds

Mammalian X and Y chromosomes independently evolved from various autosomes approximately 300 million years ago (MYA). To fully understand the relationship between genomic composition and phenotypic diversity arising due to the course of evolution, we have scanned regions of selection signatures on the X chromosome in different cattle breeds. In this study, we have prepared the datasets of 184 individuals of different cattle breeds and explored the complete X chromosome by utilizing four within-population and two between-population methods. There were 23, 25, 30, 17, 17, and 12 outlier regions identified in Tajima's D, CLR, iHS, ROH, FST, and XP-EHH. Bioinformatics analysis showed that these regions harbor important candidate genes like AKAP4 for reproduction in Brown Swiss, MBTS2 for production traits in Brown Swiss and Guernsey, CXCR3 and CITED1 for health traits in Jersey and Nelore, and BMX and CD40LG for regulation of X chromosome inactivation in Nelore and Gir. We identified genes shared among multiple methods, such as TRNAC-GCA and IL1RAPL1, which appeared in Tajima's D, ROH, and iHS analyses. The gene TRNAW-CCA was found in ROH, CLR and iHS analyses. The X chromosome exhibits a distinctive interaction between demographic factors and genetic variations, and these findings may provide new insight into the X-linked selection in different cattle breeds.

Copy number variation of individual cattle genomes using next-generation sequencing

Copy number variations (CNVs) affect a wide range of phenotypic traits; however, CNVs in or near segmental duplication regions are often intractable. Using a read depth approach based on next-generation sequencing, we examined genome-wide copy number differences among five taurine (three Angus, one Holstein, and one Hereford) and one indicine (Nelore) cattle. Within mapped chromosomal sequence, we identified 1265 CNV regions comprising ~55.6-Mbp sequence--476 of which (~38%) have not previously been reported. We validated this sequence-based CNV call set with array comparative genomic hybridization (aCGH), quantitative PCR (qPCR), and fluorescent in situ hybridization (FISH), achieving a validation rate of 82% and a false positive rate of 8%. We further estimated absolute copy numbers for genomic segments and annotated genes in each individual. Surveys of the top 25 most variable genes revealed that the Nelore individual had the lowest copy numbers in 13 cases (~52%, χ(2) test; P-value <0.05). In contrast, genes related to pathogen- and parasite-resistance, such as CATHL4 and ULBP17, were highly duplicated in the Nelore individual relative to the taurine cattle, while genes involved in lipid transport and metabolism, including APOL3 and FABP2, were highly duplicated in the beef breeds. These CNV regions also harbor genes like BPIFA2A (BSP30A) and WC1, suggesting that some CNVs may be associated with breed-specific differences in adaptation, health, and production traits. By providing the first individualized cattle CNV and segmental duplication maps and genome-wide gene copy number estimates, we enable future CNV studies into highly duplicated regions in the cattle genome.

Analysis of recent segmental duplications in the bovine genome

BACKGROUND

Duplicated sequences are an important source of gene innovation and structural variation within mammalian genomes. We performed the first systematic and genome-wide analysis of segmental duplications in the modern domesticated cattle (Bos taurus). Using two distinct computational analyses, we estimated that 3.1% (94.4 Mb) of the bovine genome consists of recently duplicated sequences (>or= 1 kb in length, >or= 90% sequence identity). Similar to other mammalian draft assemblies, almost half (47% of 94.4 Mb) of these sequences have not been assigned to cattle chromosomes.

RESULTS

In this study, we provide the first experimental validation large duplications and briefly compared their distribution on two independent bovine genome assemblies using fluorescent in situ hybridization (FISH). Our analyses suggest that the (75-90%) of segmental duplications are organized into local tandem duplication clusters. Along with rodents and carnivores, these results now confidently establish tandem duplications as the most likely mammalian archetypical organization, in contrast to humans and great ape species which show a preponderance of interspersed duplications. A cross-species survey of duplicated genes and gene families indicated that duplication, positive selection and gene conversion have shaped primates, rodents, carnivores and ruminants to different degrees for their speciation and adaptation. We identified that bovine segmental duplications corresponding to genes are significantly enriched for specific biological functions such as immunity, digestion, lactation and reproduction.

CONCLUSION

Our results suggest that in most mammalian lineages segmental duplications are organized in a tandem configuration. Segmental duplications remain problematic for genome and assembly and we highlight genic regions that require higher quality sequence characterization. This study provides insights into mammalian genome evolution and generates a valuable resource for cattle genomics research.

A first-generation microsatellite-based genetic linkage map of the Siberian jay (Perisoreus infaustus): insights into avian genome evolution

BACKGROUND

Genomic resources for the majority of free-living vertebrates of ecological and evolutionary importance are scarce. Therefore, linkage maps with high-density genome coverage are needed for progress in genomics of wild species. The Siberian jay (Perisoreus infaustus; Corvidae) is a passerine bird which has been subject to lots of research in the areas of ecology and evolutionary biology. Knowledge of its genome structure and organization is required to advance our understanding of the genetic basis of ecologically important traits in this species, as well as to provide insights into avian genome evolution.

RESULTS

We describe the first genetic linkage map of Siberian jay constructed using 117 microsatellites and a mapping pedigree of 349 animals representing five families from a natural population breeding in western Finland from the years 1975 to 2006. Markers were resolved into nine autosomal and a Z-chromosome-specific linkage group, 10 markers remaining unlinked. The best-position map with the most likely positions of all significantly linked loci had a total sex-average size of 862.8 cM, with an average interval distance of 9.69 cM. The female map covered 988.4 cM, whereas the male map covered only 774 cM. The Z-chromosome linkage group comprised six markers, three pseudoautosomal and three sex-specific loci, and spanned 10.6 cM in females and 48.9 cM in males. Eighty-one of the mapped loci could be ordered on a framework map with odds of >1000:1 covering a total size of 809.6 cM in females and 694.2 cM in males. Significant sex specific distortions towards reduced male recombination rates were revealed in the entire best-position map as well as within two autosomal linkage groups. Comparative mapping between Siberian jay and chicken anchored 22 homologous loci on 6 different linkage groups corresponding to chicken chromosomes Gga1, 2, 3, 4, 5, and Z. Quite a few cases of intra-chromosomal rearrangements within the autosomes and three cases of inter-chromosomal rearrangement between the Siberian jay autosomal linkage groups (LG1, LG2 and LG3) and the chicken sex chromosome GgaZ were observed, suggesting a conserved synteny, but changes in marker order, within autosomes during about 100 million years of avian evolution.

CONCLUSION

The constructed linkage map represents a valuable resource for intraspecific genomics of Siberian jay, as well as for avian comparative genomic studies. Apart from providing novel insights into sex-specific recombination rates and patterns, the described maps - from a previously genomically uncharacterized superfamily (Corvidae) of passerine birds - provide new insights into avian genome evolution. In combination with high-resolution data on quantitative trait variability from the study population, they also provide a foundation for QTL-mapping studies.

Assessing the relative ages of admixture in the bovine hybrid zones of Africa and the Near East using X chromosome haplotype mosaicism

Historical hybridization events between the two subspecies of cattle, Bos taurus and B. indicus, have occurred in several regions of the world, while other populations have remained nonadmixed. We typed closely linked X chromosome microsatellites in cattle populations with differing histories of admixture from Africa, Europe, the Near East, and India. Haplotype breakdown will occur as admixed populations age, and longer ancestral haplotypes will remain intact in more recently admixed populations compared to older ones. We genotyped male animals from these populations, obtaining unambiguous haplotypes, and measured levels of linkage disequilibrium (LD) and ancestral mosaicism. Extensive LD, likely to be the result of ongoing admixture, was discovered in hybrid cattle populations from the perimeter of the tsetse zone in West Africa. A Bayesian method to assign microsatellite allele ancestry was used to designate the likely origin of each chromosomal segment and assess the relative ages of admixture in the populations. A gradient of the age of admixture in the African continent emerged, where older admixture has produced more fragmented haplotypes in the south, and longer intact haplotypes, indicating more recent hybridization, feature in the northwest.

Development and mapping of microsatellite markers derived from chicken chromosome-specific libraries

Chromosome-specific painting probes and libraries were developed for chicken Macrochromosomes 1, 2, 3, and 4 by chromosome microisolation and microcloning. Fluorescent in situ hybridization results using the painting probes on normal chicken metaphase chromosomes indicated the purity and specificity of each probe. Chromosome-specific libraries for chicken Macrochromosomes 1, 2, 3, and 4 were prepared in a phage vector. Fifty-two additional unique microsatellite markers of the (AC)n type were developed from these chromosome-specific libraries. These markers were mapped on the East Lansing reference population to increase the marker density on the four macrochromosomes. Results of the current study suggest that development of markers from chromosome-specific libraries is very useful for constructing high-density linkage maps for chicken macrochromosomes.

Genomic tools to improve parasite resistance

The natural genetic variability of the ruminant immune system provides a feasible means to control gastrointestinal (GI) parasite infection without anthelmintics. However, the paradigm of traditional selection has not been effectively applied to the moderately heritable traits of parasite resistance (h approximately equal to 0.3) due to the difficulty and expense of gathering accurate phenotypes in a commercial production setting. These characteristics make host traits related to GI nematode infection ideal candidates for genomics-based research. To initiate explanation of important allelic differences, economic trait loci (ETL) are being identified and mapped using a resource population of Angus cattle segregating for GI nematode resistance and susceptibility to the two most common nematode parasites of US cattle, Ostertagia ostertagi and Cooperia oncophora. The population is composed of five generations of half-sib progeny with complete phenotypic records produced from controlled infections. To detect the genomic locations of the three distinct phenotypic traits being expressed (innately immune, acquired immune, and immunologically non-responsive), genotypes have been generated for DNA markers (N=199) spaced at regular intervals (approximately 20cm intervals) throughout the entire genome (3000cm). Although initial ETL detection may be limited by half-sib family size, the unique structure of this population provides additional statistical power for refining map position of potential ETL. After allele frequency and contribution to phenotype are determined in this population, marker tests associated with ETL most beneficial for controlling parasite infection can be accurately used for selection. Comparative map and functional genomic information from humans and other species of biomedical importance will be utilized in further investigations to elucidate the genes underlying ETL.

Partial deletion of the bovine ED1 gene causes anhidrotic ectodermal dysplasia in cattle

Anhidrotic ectodermal dysplasia (ED1) is characterized by hypotrichosis, reduced number of sweat glands, and incisior anodontia in human, mouse, and cattle. In affected humans and mice, mutations in the ED1 gene coding for ectodysplasin 1 are found. Ectodysplasin 1 is a novel trimeric transmembrane protein with an extracellular TNF-like signaling domain that is believed to be involved in the formation of hair follicles and tooth buds during fetal development. We report the construction of a 480-kb BAC contig harboring the complete bovine ED1 gene on BTA Xq22-Xq24. Physical mapping and sequence analysis of the coding parts of the ED1 gene revealed that a large genomic region including exon 3 of the ED1 gene is deleted in cattle with anhidrotic ectodermal dysplasia in a family of German Holstein cattle with three affected maternal half sibs.

An Enhanced Linkage Map of the Sheep Genome Comprising More Than 1000 Loci

A medium-density linkage map of the ovine genome has been developed. Marker data for 550 new loci were generated and merged with the previous sheep linkage map. The new map comprises 1093 markers representing 1062 unique loci (941 anonymous loci, 121 genes) and spans 3500 cM (sex-averaged) for the autosomes and 132 cM (female) on the X chromosome. There is an average spacing of 3.4 cM between autosomal loci and 8.3 cM between highly polymorphic [polymorphic information content (PIC) ≥ 0.7] autosomal loci. The largest gap between markers is 32.5 cM, and the number of gaps of >20 cM between loci, or regions where loci are missing from chromosome ends, has been reduced from 40 in the previous map to 6. Five hundred and seventy-three of the loci can be ordered on a framework map with odds of >1000 : 1. The sheep linkage map contains strong links to both the cattle and goat maps. Five hundred and seventy-two of the loci positioned on the sheep linkage map have also been mapped by linkage analysis in cattle, and 209 of the loci mapped on the sheep linkage map have also been placed on the goat linkage map. Inspection of ruminant linkage maps indicates that the genomic coverage by the current sheep linkage map is comparable to that of the available cattle maps. The sheep map provides a valuable resource to the international sheep, cattle, and goat gene mapping community.

Development and mapping of microsatellites from a microdissected BTA 11-specific DNA library

A chromosome-specific library was developed for Bos taurus autosome 11 by chromosome microdissection and microcloning using a bovine primary fibroblast culture, obtained from a t(X;23) heifer, that spontaneously developed a translocation chromosome involving bovine chromosome 11. The library was screened using (AC)12 oligos, positive clones selected, sequenced and primers developed to generate bovine chromosome 11-specific microsatellite markers. This study suggests that chromosome-specific libraries have great potential for development of microsatellite markers for the construction of marker-saturated linkage maps for each chromosome.

Evaluation of sire predicted transmitting abilities for evidence of X-chromosomal inheritance in north american sire families

This study tested for differences between genetic merits of sons and daughters of sires and for evidence of segregating quantitative trait loci on the X chromosomes of North American Holsteins. Son PTA adjusted for sire PTA was used as the dependent variable to test for biases and for genes that were passed from sire to daughter but not to son. The test of variability across sires of sons merely indicated an unaccounted source of variation, for which genes on X chromosomes might be responsible. Critical values for this test and power were determined by simulation for a variety of populations and traits differing in heritability, size of the X chromosome effect, and allelic frequency. Simulated genes on the X chromosome were detected with high power at intermediate frequencies of the favorable allele. The power of the test increased as the size of the effect increased and as genetic variance attributed to autosomes decreased. The test was then applied to recently evaluated data from US and Canadian Holstein populations. Genetic evaluations for >17,000 bulls from the US and >9000 from Canada were included. Results suggested that little extra variation was present for some traits formally evaluated in North America, but that genes on the X chromosome were unlikely to be the cause.

Comparative cytogenetic mapping reveals chromosome rearrangements between the X chromosomes of two closely related mammalian species (cattle and goats)

Cytogenetic localization of 24 BACs containing type I (genes and ESTs) and type II (microsatellites) markers were used to construct cytogenetic maps of caprine (CHI) and bovine (BTA) X chromosomes. Comparison of these two maps revealed that the distal region of the goat X long arm (CHI Xq38-->q42) was located inside the bovine X chromosome, between PGK1 (BTA Xq25) and DVEPC137 (BTA Xq12). The marker order was globally conserved without any pericentric inversion, as previously postulated in the literature. The caprine centromere was found between DVEPC053 and DVEPC102 (belonging to the same band in the bovine X: BTA Xq41), whereas the bovine centromere was between DVEPC076 and DVEPC132, belonging to the same region of the caprine X chromosome (CHI Xq31-->q33). The pseudoautosomal region was situated at the tip of the bovine X long arm and on the tiny short arm of the caprine X chromosome. In the non-pseudoautosomal (NPA) region, the synteny of coding sequences was well conserved between the human species and the two ruminant species, but the gene order was dramatically divergent. It is suggested that the 24 BACs of this study could constitute a new tool to measure phylogenetic distances between different mammalian species by comparing chromosome rearrangements inside the NPA region of the X.

A second-generation linkage map of the bovine genome

We report a bovine linkage map constructed with 1236 polymorphic DNA markers and 14 erythrocyte antigens and serum proteins. The 2990-cM map consists of a sex-specific, X chromosome linkage group and 29 sex-averaged, autosomal linkage groups with an average interval size of 2.5 cM. The map contains 627 new markers and 623 previously linked markers, providing a basis for integrating the four published bovine maps. Orientation and chromosomal assignment of all the linkage groups, except BTA20 and BTA22, was provided by 88 markers that were assigned previously to chromosomes. This map provides sufficient marker density for genomic scans of populations segregating quantitative trait loci (QTL) and subsequent implementation of marker-assisted selection (MAS) mating schemes.

Targeted development of microsatellite markers from the defined region of bovine chromosome 6q21-31

A methodical strategy for the isolation of microsatellite markers specific for targeted regions of bovine chromosomes is presented. The procedure involves directed microdissection of one defined subchromosomal area, its DOP-PCR-amplification and cloning. With this approach, a library specific to the BTA 6q21-31 chromosomal region was constructed. Eleven unique microsatellite-containing sequences were isolated, converted into sequence-tagged microsatellite sites, and characterized concerning their species-specific origin. Seven primer pairs generated bovine-specific PCR products and provided a set of microsatellite markers that generally revealed high informativity in the HF breed. Linkage analysis assigned six of them to their predefined subchromosomal origin on BTA 6 corresponding to the specific rehybridization signal of the DOP-PCR product generated from the microdissected chromosome area 6q21-31. The results underline the usefulness of the BTA 6q21-31 library for targeted isolation of unique sequences that are specific for the dissected chromosomal region as demonstrated here by the isolation of microsatellite markers.

A chromosome-specific microdissected library increases marker density on bovine chromosome 1

Genetic resolution of bovine chromosome 1 (BTA1) linkage group was significantly increased by screening for microsatellite clones a microdissected library constructed from a bovine cell line carrying a t(1;29) translocation. Eighty-five percent of the microsatellites (ms) (46/54) identified were informative in the USDA/MARC mapping population, and 96% of these ms (44/46) linked to BTA1 (LOD > 3.0). When merged with 40 existing BTA1 markers the genetic map spanned 153.8 cM (sex-averaged interval, 1.9 cM). The fourfold improvement in marker density of BTA1 provides a genetic map that enhances mapping of quantitative trait loci and implementation of marker assisted selection.