Building human genome maps with radiation hybrids

Sequence analysis in Bos taurus reveals pervasiveness of X-Y arms races in mammalian lineages

Studies of Y Chromosome evolution have focused primarily on gene decay, a consequence of suppression of crossing-over with the X Chromosome. Here, we provide evidence that suppression of X-Y crossing-over unleashed a second dynamic: selfish X-Y arms races that reshaped the sex chromosomes in mammals as different as cattle, mice, and men. Using super-resolution sequencing, we explore the Y Chromosome of Bos taurus (bull) and find it to be dominated by massive, lineage-specific amplification of testis-expressed gene families, making it the most gene-dense Y Chromosome sequenced to date. As in mice, an X-linked homolog of a bull Y-amplified gene has become testis-specific and amplified. This evolutionary convergence implies that lineage-specific X-Y coevolution through gene amplification, and the selfish forces underlying this phenomenon, were dominatingly powerful among diverse mammalian lineages. Together with Y gene decay, X-Y arms races molded mammalian sex chromosomes and influenced the course of mammalian evolution.

Avian W and mammalian Y chromosomes convergently retained dosage-sensitive regulators

After birds diverged from mammals, different ancestral autosomes evolved into sex chromosomes in each lineage. In birds, females are ZW and males ZZ, but in mammals females are XX and males XY. We sequenced the chicken W chromosome, compared its gene content with our reconstruction of the ancestral autosomes, and followed the evolutionary trajectory of ancestral W-linked genes across birds. Avian W chromosomes evolved in parallel with mammalian Y chromosomes, preserving ancestral genes through selection to maintain the dosage of broadly-expressed regulators of key cellular processes. We propose that, like the human Y chromosome, the chicken W chromosome is essential for embryonic viability of the heterogametic sex. Unlike other sequenced sex chromosomes, the chicken W did not acquire and amplify genes specifically expressed in reproductive tissues. We speculate that the pressures that drive the acquisition of reproduction related genes on sex chromosomes may be specific to the male germ line.

Sequencing the mouse Y chromosome reveals convergent gene acquisition and amplification on both sex chromosomes

We sequenced the MSY (male-specific region of the Y chromosome) of the C57BL/6J strain of the laboratory mouse Mus musculus. In contrast to theories that Y chromosomes are heterochromatic and gene poor, the mouse MSY is 99.9% euchromatic and contains about 700 protein-coding genes. Only 2% of the MSY derives from the ancestral autosomes that gave rise to the mammalian sex chromosomes. Instead, all but 45 of the MSY's genes belong to three acquired, massively amplified gene families that have no homologs on primate MSYs but do have acquired, amplified homologs on the mouse X chromosome. The complete mouse MSY sequence brings to light dramatic forces in sex chromosome evolution: lineage-specific convergent acquisition and amplification of X-Y gene families, possibly fueled by antagonism between acquired X-Y homologs. The mouse MSY sequence presents opportunities for experimental studies of a sex-specific chromosome in its entirety, in a genetically tractable model organism.

Mammalian Y chromosomes retain widely expressed dosage-sensitive regulators

The human X and Y chromosomes evolved from an ordinary pair of autosomes, but millions of years ago genetic decay ravaged the Y chromosome, and only three percent of its ancestral genes survived. We reconstructed the evolution of the Y chromosome across eight mammals to identify biases in gene content and the selective pressures that preserved the surviving ancestral genes. Our findings indicate that survival was non-random, and in two cases, convergent across placental and marsupial mammals. We conclude that the Y chromosome's gene content became specialized through selection to maintain the ancestral dosage of homologous X-Y gene pairs that function as broadly expressed regulators of transcription, translation and protein stability. We propose that beyond its roles in testis determination and spermatogenesis, the Y chromosome is essential for male viability, and plays unappreciated roles in Turner syndrome and in phenotypic differences between the sexes in health and disease.

Strict evolutionary conservation followed rapid gene loss on human and rhesus Y chromosomes

The human X and Y chromosomes evolved from an ordinary pair of autosomes during the past 200–300 million years^1–3. Due to genetic decay, the human MSY (male-specific region of Y chromosome) retains only three percent of the ancestral autosomes’ genes^4,5. This evolutionary decay was driven by a series of five “stratification” events. Each event suppressed X-Y crossing over within a chromosome segment or “stratum”, incorporated that segment into the MSY, and subjected its genes to the erosive forces that attend the absence of crossing over^2,6. The last of these events occurred 30 million years ago (mya), or 5 million years before the human and Old World monkey (OWM) lineages diverged. Although speculation abounds regarding ongoing decay and looming extinction of the human Y chromosome^7–10, remarkably little is known about how many MSY genes were lost in the human lineage in the 25 million years that have followed its separation from the OWM lineage. To explore this question, we sequenced the MSY of the rhesus macaque, an OWM, and compared it to the human MSY. We discovered that, during the last 25 million years, MSY gene loss in the human lineage was limited to the youngest stratum (stratum 5), which comprises three percent of the human MSY. Within the older strata, which collectively comprise the bulk of the human MSY, gene loss evidently ceased more than 25 mya. Likewise, the rhesus MSY has not lost any older genes (from strata 1–4) during the past 25 million years, despite major structural differences from the human MSY. The rhesus MSY is simpler, with few amplified gene families or palindromes that might enable intrachromosomal recombination and repair. We present an empirical reconstruction of human MSY evolution in which each stratum transitioned from rapid, exponential loss of ancestral genes to strict conservation through purifying selection.

Enhancing radiation hybrid mapping through whole genome amplification

Radiation hybrid (RH) mapping is limited by the inherent genomic instability of RH clones entailing both, limited DNA sample amounts and genomic heterogeneity of the clones. Here the instability of RH clones is quantified and the suitability of the multiple strand displacement whole genome amplification method (WGA) for radiation hybrid mapping is assessed. To quantify the instability of RH clones, eleven clones of a 10,000(Rad) rhesus macaque radiation hybrid panel were passaged ten times and analyzed by interspersed repeat sequence specific quantitative PCR and by genotyping of 46 macaque chromosome 5 STS markers. The quantitative PCR data indicate an average loss of 55% of the donor DNA over 10 passages. Over the same period, a dropout of 46.2% of the STS markers was observed. These data indicate a genome wide half-life of the donor DNA of 8.7 passages and of 10.6 passages for the chromosome 5 markers. The genotyping data of the genomic RH DNA were compared to three sets of WGA experiments: 1) single wgaDNA amplifications, 2) six WGA replicates, and 3) re-amplification of wga DNA. The assays demonstrated concordance rates of 97.6%, 98% and 99.3%, respectively, and indicated the marker specificity of some repeated WGA dropouts. The study confirms that WGA is suitable for RH mapping studies should enable the accurate analysis of almost an infinite numbers of markers. WGA will allow the analysis of earliest RH clone passages, thus limiting their heterogeneity and RH mapping artifacts.

Cloning and expression of activation induced cytidine deaminase from Bos taurus

Activation induced cytidine deaminase is an enzyme crucial to somatic hypermutation and gene conversion, processes that are essential for the diversification of Ig V genes. The bovine Ig repertoire appears to be diversified by mechanisms that are significantly different to those that operate in humans and mice. This study set out to test the hypothesis that differences in the organization, coding sequence, expression or genomic location of the bovine AICDA gene enables the encoded enzyme to catalyse the unusual Ig diversification mechanism seen in cattle as well as conventional antigen-driven mutation. Characterization of bovine AICDA excluded the first two possibilities. AICDA expression was detected in lymphoid tissues from neonatal and older cattle, but AICDA cDNA could not be detected in muscle tissue. The pattern of gene expression did not therefore differ from that in other vertebrates. The AICDA cDNA was cloned and expressed successfully in Escherichia coli generating a phenotype consistent with the mutating action of this deaminase. Using a whole genome radiation hybrid panel, bovine AICDA was mapped to a region of bovine chromosome 5 syntenic with the location of human AICDA on chromosome 12. We conclude that the unusual nature of Ig diversification in cattle is unlikely to be attributable to the structure, sequence, activity or genomic location of bovine AICDA.

Assessment of AFLP marker behaviour in enriching STS radiation hybrid maps

Radiation hybrid (RH) mapping provides a powerful tool to build high-resolution maps of genomes. Here, we demonstrate the use of the AFLP technique for high-throughput typing of RH cell lines. Cattle were used as the model species because an RH panel was available to investigate the behaviour of AFLP markers within the microsatellite- and STS-based maps of this species. A total of 747 AFLP markers were typed on the TM112 RH radiation panel and 651 of these were assigned by two-point analysis to the 29 bovine autosomes and sex chromosomes. AFLP markers were added to the 1222 microsatellite and STS markers that were included in earlier RH maps. Multipoint maps were constructed for seven example chromosomes, which retained 248 microsatellite and STS markers, and added 123 AFLP markers at LOD 4. The addition of the AFLP markers increased the number of markers by 42.1% and the map length by 10.4%. The AFLP markers showed lower retention frequency (RF) values than the STS markers. The comparison of RF values in AFLP markers and their corresponding AFLP-derived STSs demonstrated that the lower RF values were due to the lower detection sensitivity of the AFLP technique. Despite these differences, AFLP and AFLP-derived STS markers mapped to identical or similar positions. These results demonstrate that it is possible to merge AFLP and microsatellite markers in the same map. The application of AFLP technology could permit the rapid construction of RH maps in species for which extensive genome information and large numbers of SNP and microsatellite markers are not available.

An original SERPINA3 gene cluster: elucidation of genomic organization and gene expression in the Bos taurus 21q24 region

Background

The superfamily of serine proteinase inhibitors (serpins) is involved in numerous fundamental biological processes as inflammation, blood coagulation and apoptosis. Our interest is focused on the SERPINA3 sub-family. The major human plasma protease inhibitor, α1-antichymotrypsin, encoded by the SERPINA3 gene, is homologous to genes organized in clusters in several mammalian species. However, although there is a similar genic organization with a high degree of sequence conservation, the reactive-centre-loop domains, which are responsible for the protease specificity, show significant divergences.

Results

We provide additional information by analyzing the situation of SERPINA3 in the bovine genome. A cluster of eight genes and one pseudogene sharing a high degree of identity and the same structural organization was characterized. Bovine SERPINA3 genes were localized by radiation hybrid mapping on 21q24 and only spanned over 235 Kilobases. For all these genes, we propose a new nomenclature from SERPINA3-1 to SERPINA3-8. They share approximately 70% of identity with the human SERPINA3 homologue. In the cluster, we described an original sub-group of six members with an unexpected high degree of conservation for the reactive-centre-loop domain, suggesting a similar peptidase inhibitory pattern. Preliminary expression analyses of these bovSERPINA3s showed different tissue-specific patterns and diverse states of glycosylation and phosphorylation. Finally, in the context of phylogenetic analyses, we improved our knowledge on mammalian SERPINAs evolution.

Conclusion

Our experimental results update data of the bovine genome sequencing, substantially increase the bovSERPINA3 sub-family and enrich the phylogenetic tree of serpins. We provide new opportunities for future investigations to approach the biological functions of this unusual subset of serine proteinase inhibitors.

Mapping, expression, and association study of the bovine PSMC1 gene

The 26S ATP-dependent protease is composed of a 20S catalytic proteasome and two PA700 regulatory modules; it plays a central role in many regulatory pathways, such as cell cycle regulation, differentiation, and apoptosis. The PA700 complex is composed of multiple subunits, including at least six related ATPases and approximately 15 non-ATPase polypeptides. PSMC1 (proteasome 26S subunit, ATPase, 1) is one of these ATPases. In this study, we amplified a fragment of 507 bp from intron 9 of the bovine PSMC1 gene and found a SNP (G/A) at position 216 in the PCR fragment. Genotyping of 138 animals from four beef breeds revealed that the average frequency for allele A (G-base) was 0.4271 (0.3269-0.5517); for allele B (A-base) it was 0.5729 (0.4483-0.6731). This SNP is significantly associated with average daily feed intake (P < 0.01), average daily gain, finishing average daily gain, body length, ratio of feed to meat, backfat thickness, and loin-muscle area (P < 0.05). Our experimental data showed that animals with an AA genotype have a significantly lower food intake, grow faster, are longer in the body, and have less backfat and bigger loin muscle; hence, their ratio of feed to meat is significantly lower. We believe that the PSMC1 SNP is a potential candidate marker for marker-assisted selection in these traits. We also found that the bovine PSMC1 gene was expressed mainly in lung, testis, and spleen. In addition, we mapped the bovine PSMC1 gene on BTA10 by an RH mapping method.

Evaluation of Compass as a comparative mapping tool for ESTs using horse radiation hybrid maps

Loci for 9322 equine expressed sequence tags (ESTs) were predicted using the Comparative Mapping by Annotation and Sequence Similarity (Compass) strategy in order to evaluate the programme's ability to make accurate locus predictions in species with comparative gene maps. Using human genome sequence information from Build 35 (May 2004) and published marker information from the radiation hybrid (RH) maps for equine chromosomes (ECA) 17 and X, 162 ESTs were predicted to locations on ECA17 and 328 ESTs to locations on ECAX by selection of the 'top blast hit'. The locations of 30 ESTs were assessed experimentally by RH mapping analysis to evaluate the accuracy of the Compass predictions. The data revealed that 53% (16 of 30) of the ESTs predicted on ECA17 and ECAX mapped to those chromosomes. Analysis of the results suggested the need to identify expressed orthologous sequences in order to generate more accurate predictions for ESTs. Locus predictions were reassessed with three modifications to the Compass strategy's orthologue selection parameters. Selection of the 'top gene hit' improved accuracy to 72% (21 of 29), while selection of the 'top expressed gene hit' improved accuracy to 86% (24 of 28). Using the default Compass parameters with the UniGene database improved prediction accuracy to 96% (22 of 23); however, this level of accuracy came with a substantial decrease in the total number of predictions. When used with optimized prediction parameters, the Compass strategy can be a practical in silico map location prediction tool for large EST sample sets from unsequenced animal genomes.

A new 4016-marker radiation hybrid map for porcine-human genome analysis

We constructed a 5000-rad comprehensive radiation hybrid (RH) map of the porcine (Sus scrofa) genome and compared the results with the human genome. Of 4475 typed markers, 4016 (89.7%) had LOD >5 compared with the markers used in our previous RH map by means of two-point analysis and were grouped onto the 19 porcine chromosomes (SSCs). All mapped markers had LOD >3 as determined by RHMAPPER analysis. The current map comprised 430 microsatellite (MS) framework markers, 914 other MS markers, and 2672 expressed sequence tags (ESTs). The whole-genome map was 8822.1 cR in length, giving an average marker density of 0.342 Mb/cR. The average retention frequency was 35.8%. Using BLAST searches of porcine ESTs against the RefSeq human nucleotide and amino acid sequences (release 22), we constructed high-resolution comparative maps of each SSC and each human chromosome (HSA). The average distance between ESTs in the human genome was 1.38 Mb. SSC contained 50 human chromosomal syntenic groups, and SSC11, SSC12, and SSC16 were only derived from the HSA13q, HSA17, and HSA5 regions, respectively. Among 38 porcine terminal regions, we found that at least 20 regions have been conserved between the porcine and human genomes; we also found four paralogous regions for the major histocompatibility complex (MHC) on SSC7, SSC2, SSC4, and SSC1.

Identification of differentially expressed genes in distinct skeletal muscles in cattle using cDNA microarray

The 788-gene microarray was manufactured using selected elements from three different cDNA libraries in order to identify molecular processes that determine phenotypic characteristics between loin (M. longissimus thoracis) and round (M. semimembranosus) muscles. Microarray analyses identified 24 differentially expressed genes between the two muscles investigated. Five of the genes were verified by quantitative RT-PCR and three of them were mapped on bovine chromosomes using 5,000 rad bovine radiation hybrid (RH) panel. The map locations indicated that they were mapped in the same chromosomal regions where IMF and growth QTLs were located, suggesting that they are most possible positional candidate genes for the traits.

Radiation hybrid mapping and comparative sequence analysis of bovine RIG-I and MAVS genes

Retinoic acid inducible gene I (RIG-I) and mitochondrial antiviral signaling (MAVS) proteins have recently been found to operate in a pathway for the detection and subsequent elimination of replicating viral genomes. Because of this innate immunity role, RIG-I and MAVS are candidates for studies of disease resistance. The objectives of this work were to (1) radiation hybrid (RH) map bovine RIG-I and MAVS and (2) perform comparative sequence analysis of partial genomic sequence from each gene. Using a bovine 5000(rad) RH panel, RIG-I was localized to BTA08 (LOD > 12) and MAVS was localized to BTA13 (LOD > 12). RIG-I exon 14 and partial MAVS exon five were sequenced in nine breeds and compared with available sequence from the Bovine Genome Project. RIG-I exon 14 and partial MAYS exon five were conserved in all samples examined. One T-A transversion SNP was found in intronic sequence downstream of RIG-I exon 14.

rh_tsp_map 3.0: end-to-end radiation hybrid mapping with improved speed and quality control

rh_tsp_map is a software package for computing radiation hybrid (RH) maps and for integrating physical and genetic maps. It solves the central mapping instances by reducing them to the traveling salesman problem (TSP) and using a modification of the CONCORDE package to solve the TSP instances. We present some of the features added between the initial rh_tsp_map version 1.0 and the current version 3.0, emphasizing the automation of many steps and addition of various checks designed to find problems with the input data. Iterations of improved input data followed by fast re-computation of the maps improves the quality of the final maps.

AVAILABILITY

rh_tsp_map source code and documentation including a tutorial is available at ftp://ftp.ncbi.nih.gov/pub/agarwala/rhmapping/rh_tsp_map.tar.gz. CONCORDE modified for RH mapping is available in the directory http://www.isye.gatech.edu/~wcook/rh/. The QSopt library needed for CONCORDE is available at http://www2.isye.gatech.edu/~wcook/qsopt/downloads/downloads.htm

A high-resolution comparative radiation hybrid map of equine chromosome 4q12-q22

In this study, we present a comprehensive 5000-rad radiation hybrid map of a 40-cM region on equine chromosome 4 (ECA4) that contains quantitative trait loci for equine osteochondrosis. We mapped 29 gene-associated sequence tagged site markers using primers designed from equine expressed sequence tags or BAC clones in the ECA4q12-q22 region. Three blocks of conserved synteny, showing two chromosomal breakpoints, were identified in the segment of ECA4q12-q22. Markers from other segments of HSA7q mapped to ECA13p and ECA4p, and a region of HSA7p was homologous to ECA13p. Therefore, we have improved the resolution of the human-equine comparative map, which allows the identification of candidate genes underlying traits of interest.

Molecular characterization of the equine collagen, type IX, alpha 2 (COL9A2) gene on horse chromosome 2p16-->p15

The mammalian collagen, type IX, alpha 2 gene (COL9A2) encodes the alpha-2 chain of type IX collagen and is located on horse chromosome 2p16-->p14 harbouring a quantitative trait locus for osteochondrosis. We isolated a bacterial artificial chromosome (BAC) clone containing the equine COL9A2 gene and determined the complete genomic sequence of this gene. Cloning and characterization of equine COL9A2 revealed that the equine gene consists of 32 exons spanning approximately 15 kb. The COL9A2 transcript encodes a single protein of 688 amino acids. Thirty two single nucleotide polymorphisms (SNPs) equally distributed in the gene were detected in a mutation scan of eight unrelated Hanoverian warmblood stallions, including one SNP that affects the amino acid sequence of COL9A2. Comparative analyses between horse, human, mouse and rat indicate that the chromosomal location of equine COL9A2 is in agreement with known chromosomal synteny relationships. The comparison of the gene structure and transcript revealed a high degree of conservation towards the other mammalian COL9A2 genes. We chose three informative SNPs for association and linkage disequilibrium tests in three to five paternal half-sib families of Hanoverian warmblood horses consisting of 44 to 75 genotyped animals. The test statistics did not reach the significance threshold of 5% and so we could not show an association of COL9A2 with equine osteochondrosis.

A second generation radiation hybrid map to aid the assembly of the bovine genome sequence

BACKGROUND

Several approaches can be used to determine the order of loci on chromosomes and hence develop maps of the genome. However, all mapping approaches are prone to errors either arising from technical deficiencies or lack of statistical support to distinguish between alternative orders of loci. The accuracy of the genome maps could be improved, in principle, if information from different sources was combined to produce integrated maps. The publicly available bovine genomic sequence assembly with 6x coverage (Btau_2.0) is based on whole genome shotgun sequence data and limited mapping data however, it is recognised that this assembly is a draft that contains errors. Correcting the sequence assembly requires extensive additional mapping information to improve the reliability of the ordering of sequence scaffolds on chromosomes. The radiation hybrid (RH) map described here has been contributed to the international sequencing project to aid this process.

RESULTS

An RH map for the 30 bovine chromosomes is presented. The map was built using the Roslin 3000-rad RH panel (BovGen RH map) and contains 3966 markers including 2473 new loci in addition to 262 amplified fragment-length polymorphisms (AFLP) and 1231 markers previously published with the first generation RH map. Sequences of the mapped loci were aligned with published bovine genome maps to identify inconsistencies. In addition to differences in the order of loci, several cases were observed where the chromosomal assignment of loci differed between maps. All the chromosome maps were aligned with the current 6x bovine assembly (Btau_2.0) and 2898 loci were unambiguously located in the bovine sequence. The order of loci on the RH map for BTA 5, 7, 16, 22, 25 and 29 differed substantially from the assembled bovine sequence. From the 2898 loci unambiguously identified in the bovine sequence assembly, 131 mapped to different chromosomes in the BovGen RH map.

CONCLUSION

Alignment of the BovGen RH map with other published RH and genetic maps showed higher consistency in marker order and chromosome assignment than with the current 6x sequence assembly. This suggests that the bovine sequence assembly could be significantly improved by incorporating additional independent mapping information.

An integrated radiation hybrid map of bovine chromosome 18 that refines a critical region associated with multiple ocular defects in cattle

Congenital multiple ocular defects (MOD) of Japanese black cattle is a hereditary ocular disorder with an autosomal recessive mode of inheritance showing developmental defects of the lens, retina and iris, persistent embryonic eye vascularization and microphthalmia. The MOD locus has been mapped by linkage analysis to a 6.6-cM interval on the proximal end of bovine chromosome 18, which corresponds to human chromosome 16q and mouse chromosome 8. To refine the MOD region in cattle, we constructed an integrated radiation hybrid (RH) map of the proximal region of bovine chromosome 18, which consisted of 17 genes and 10 microsatellite markers, using the SUNbRH7000 panel. Strong conservation of gene order was found among the corresponding chromosomal regions in cattle, human and mouse. The MOD-critical region was fine mapped to a 59.5-cR region that corresponds to a 6.3-Mb segment of human chromosome 16 and a 4.8-Mb segment of mouse chromosome 8. Several positional candidate genes, including FOXC2 and USP10, were identified in this region.

Identification of genetic variation and putative regulatory regions in bovine CARD15

Mutations in caspase recruitment domain 15 (CARD15) are associated with susceptibility to Crohn's disease and Blau Syndrome. We performed comparative analyses of the bovine, murine, and human CARD15 transcripts to elucidate functionality of bovine CARD15 and examine its potential role in bovine disease resistance. Comparative analyses of intronic sequence across seven divergent species were performed to identify putative regulatory element binding motifs. High levels of interspecies conservation in sequence, genomic structure, and protein domains were detected indicating common functionality for CARD15 in cattle, human, and mouse. We identified species-specific regulatory elements in the 5' and 3' untranslated regions, suggesting that modes of regulation may have diverged across species. Thirty-one conserved putative regulatory element binding motifs were identified in the CARD15 intronic sequence of seven species. To assess the extent of genetic diversity within bovine CARD15, 41 individuals from two subspecies were sequenced and screened for polymorphisms. Thirty-six single nucleotide polymorphisms (SNPs) were identified. Finally, 20 subspecies-specific haplotypes were predicted with 7 and 13 unique haplotypes explaining the diversity within B. taurus taurus and B. taurus indicus animals, respectively. Strong evidence for a simple causal relationship between these SNP loci and their haplotypes with Johne's disease was not detected.

Identification and gene expression of the bovine C-type lectin Dectin-1

C-type lectin receptors (CTLR) are cell-surface signalling molecules that recognize a range of highly conserved pathogen molecules and instigate the appropriate immune response. Here, we report the cloning, sequencing, mapping and expression pattern of the bovine C-type lectin domain family 7, member A (CLEC7A; synonyms CLCSF12, Dectin-1). We identified two isoforms, similar to the human system, with a long and short neck. Overall, the organization of the two bovine CLEC7A genes is similar to that of humans and mice. The CLEC7A gene maps on Bos taurus chromosome 5 (BTA5). mRNA transcripts for CLEC7A were detected in bone-marrow cells, monocytes, macrophages and dendritic cells and NK cells, but not in CD4(+) T-cells or CD21(+) B-cells. The increased knowledge of the genomic organization of the bovine CTLR genes may promote our understanding of their evolution and help in the identification of bovine genes underlying disease-resistance traits.

A human-horse comparative map based on equine BAC end sequences

In an effort to increase the density of sequence-based markers for the horse genome we generated 9473 BAC end sequences (BESs) from the CHORI-241 BAC library with an average read length of 677 bp. BLASTN searches with the BESs revealed 4036 meaningful hits (E <or= 10(-5)) in the human genome that provide useful markers for the human-horse comparative map. The 4036 BLASTN hits allowed the anchoring of 3079 BAC clones to the human genome, on average one corresponding equine BAC clone per megabase of human DNA. We used the BLASTN anchored BESs for an in silico prediction of the gene content and chromosome assignment of comparatively mapped equine BAC clones. As a first verification of our in silico mapping strategy we placed 19 equine BESs with matches to HSA6 onto the RH map. All markers were assigned to the predicted localizations on ECA10, ECA20, and ECA31, respectively.

Sequence analysis of a 212 kb defensin gene cluster on ECA 27q17

Defensins are a family of evolutionary ancient antimicrobial peptides consisting of three sub-families: alpha-, beta- and theta-defensins. This investigation was focused on the genomic characterization of equine beta-defensins and the investigation of the potential clustering of beta-defensin genes in the equine genome. Six genomic BAC clones were isolated from the CHORI-241 library and one of these was mapped by FISH to ECA 27q17. This location was confirmed by RH-mapping. The contiguous 212 kb sequence of this clone was determined. Sequence analysis revealed the identification of ten pseudogenes and nine genes, six of which were highly homologous to human beta-defensin DEFB4. Clustering of the beta-defensin genes was confirmed and the order of the genes on the analyzed BAC was related to the corresponding defensin cluster on HSA 8. The knowledge about the sequence and the genomic structure of the equine beta-defensin genes will improve the classification of different paralogous defensin genes and is a prerequisite for subsequent functional studies. Additionally, the first alpha-defensin-like sequence outside the groups of primates, lagomorphs and rodents (glires) was identified.

A gene-based high-resolution comparative radiation hybrid map as a framework for genome sequence assembly of a bovine chromosome 6 region associated with QTL for growth, body composition, and milk performance traits

Background

A number of different quantitative trait loci (QTL) for various phenotypic traits, including milk production, functional, and conformation traits in dairy cattle as well as growth and body composition traits in meat cattle, have been mapped consistently in the middle region of bovine chromosome 6 (BTA6). Dense genetic and physical maps and, ultimately, a fully annotated genome sequence as well as their mutual connections are required to efficiently identify genes and gene variants responsible for genetic variation of phenotypic traits. A comprehensive high-resolution gene-rich map linking densely spaced bovine markers and genes to the annotated human genome sequence is required as a framework to facilitate this approach for the region on BTA6 carrying the QTL.

Results

Therefore, we constructed a high-resolution radiation hybrid (RH) map for the QTL containing chromosomal region of BTA6. This new RH map with a total of 234 loci including 115 genes and ESTs displays a substantial increase in loci density compared to existing physical BTA6 maps. Screening the available bovine genome sequence resources, a total of 73 loci could be assigned to sequence contigs, which were already identified as specific for BTA6. For 43 loci, corresponding sequence contigs, which were not yet placed on the bovine genome assembly, were identified. In addition, the improved potential of this high-resolution RH map for BTA6 with respect to comparative mapping was demonstrated. Mapping a large number of genes on BTA6 and cross-referencing them with map locations in corresponding syntenic multi-species chromosome segments (human, mouse, rat, dog, chicken) achieved a refined accurate alignment of conserved segments and evolutionary breakpoints across the species included.

Conclusion

The gene-anchored high-resolution RH map (1 locus/300 kb) for the targeted region of BTA6 presented here will provide a valuable platform to guide high-quality assembling and annotation of the currently existing bovine genome sequence draft to establish the final architecture of BTA6. Hence, a sequence-based map will provide a key resource to facilitate prospective continued efforts for the selection and validation of relevant positional and functional candidates underlying QTL for milk production and growth-related traits mapped on BTA6 and on similar chromosomal regions from evolutionary closely related species like sheep and goat. Furthermore, the high-resolution sequence-referenced BTA6 map will enable precise identification of multi-species conserved chromosome segments and evolutionary breakpoints in mammalian phylogenetic studies.

Bovine NALP5, NALP8, and NALP9 Genes: Assignment to a QTL Region and the Expression in Adult Tissues, Oocytes, and Preimplantation Embryos

A 3204-bp full-length cDNA of bovine NALP9 was cloned and its genomic organization was analyzed. The 2988-bp open reading frame covers 9 exons and encodes a deduced protein of 996 amino acids containing Pyrin, Nacht and leucine-rich repeat domains like the human NALP gene family members. Mapping with the WGRH5000 panel and fluorescence in situ hybridization assigned NALP9 in close vicinity to BM2078 (LOD score 25.71; distance 0.0 cR5000) on bovine chromosome 18, BTA18q25-q26, within a previously identified QTL region for reproductive traits flanked by the bovine marker BM2078 and TGLA227. BAC contig analysis revealed that NALP9, NALP8, and NALP5 map in this QTL region. Temporospatial expression of these members of the NALP gene family was monitored. Among the adult tissues examined, transcripts of NALP8 and NALP9 were detected exclusively in testis and ovary, whereas transcripts of the NALP5 gene are limited to the ovary. The transcripts of NALP9, NALP8, and NALP5 were detected in oocytes before and after in vitro maturation and with a gradual decline from 2-cell to 8-cell stage, suggesting no reactivation at the time of bovine maternal to embryonic transition. Assignment to a QTL region for reproductive traits and preferential expression of NALP9, NALP8, and NALP5 in oocyte, germinal lineage, and gonad cells may suggest their functional relevance to reproduction and possible contribution to phenotypic variation.

Constructing the parental linkage phase and the genetic map over distances <1 cM using pooled haploid DNA

A new statistical approach for construction of the genetic linkage map and estimation of the parental linkage phase based on allele frequency data from pooled gametic (sperm or egg) samples is introduced. This method can be applied for estimation of recombination fractions (over distances <1 cM) and ordering of large numbers (even hundreds) of closely linked markers. This method should be extremely useful in species with a long generation interval and a large genome size such as in dairy cattle or in forest trees; the conifer species have haploid tissues available in megagametophytes. According to Mendelian expectation, two parental alleles should occur in gametes in 1:1 proportions, if segregation distortion does not occur. However, due to mere sampling variation, the observed proportions may deviate from their expected value in practice. These deviations and their dependence along the chromosome can provide information on the parental linkage phase and on the genetic linkage map. Usefulness of the method is illustrated with simulations. The role of segregation distortion as a source of these deviations is also discussed. The software implementing this method is freely available for research purposes from the authors.

A comprehensive radiation hybrid map of the bovine genome comprising 5593 loci

A bovine whole genome 7000-rad radiation hybrid (RH) panel, SUNbRH(7000-rad), was constructed to build a high-resolution RH map. The Shirakawa-USDA linkage map served as a scaffold to construct a framework map of 3216 microsatellites on which 2377 ESTs were ordered. The resulting RH map provided essentially complete coverage across the genome, with 1 cR7000 corresponding to 114 kb, and a cattle-human comparative map of 1716 bovine genes and sequences annotated in the human genome, which covered 79 and 72% of the bovine and human genomes, respectively. We then integrated the bovine RH and comparative maps with BAC fingerprint information in to construct a detailed, BAC-based physical map covering a reported 40-cM quantitative trait locus region for intramuscular fat or "marbling" on BTA 4. In summary, the new, high-resolution SUNbRH7000-rad, comparative, Shirakawa-USDA linkage, and BAC fingerprint maps provide a set of genomic tools for fine mapping regions of interest in cattle.

Characterization and RH mapping of bovine microsatellites generated from a microdissected BTA20-specific DNA library

Bovine chromosome 20 (BTA20) is associated with several quantitative trait loci (QTL) for meat tenderness, birth weight, milk yield and composition. Fine mapping of these QTL requires the development of additional informative markers to increase the resolution of the BTA20 genetic and physical maps. A BTA20-specific library was constructed by means of microdissection and microcloning, and screened for dinucleotide repeats with (CA)16 and (GT)16 oligos. A total of 60 new microsatellites (MS) were developed and characterized for polymorphism using the U.S. Department of Agriculture (USDA)/Meat Animal Research Center (MARC) bovine reference family, of which 53 markers were informative in this family. The number of alleles for these loci varied from 1 to 14, with an average of 6.5. Thirty-three of these MSs, together with 105 markers previously mapped to BTA20, were scored on a 7000-rad cattle-hamster whole-genome radiation hybrid panel (SUNbRH), resulting in a high-resolution RH7000 rad map for BTA20.

High-resolution comparative mapping between human chromosomes 4 and 8 and bovine chromosome 27 provides genes and segments serving as positional candidates for udder health in cattle

To get more information about the order of genes located in Bos taurus (BTA) chromosome 27 segments, supposed to harbor loci influencing clinical mastitis and somatic cell count, and to identify genes that serve as positional candidates for the mentioned traits, we constructed a high-resolution, comparative, and comprehensive gene map for BTA27. The map includes 57 loci in a 5000-rad cattle-hamster whole genome radiation hybrid panel supported by 50 syntenic assignments in a cattle-murine somatic hybrid cell panel. Thirty-eight new loci (36 genes, 2 microsatellites) together with repeated mappings of 5 genes and 7 microsatellites and integration of existing data from 7 microsatellites were used to generate a comprehensive RH5000 map. The RH map, constructed at lod score criterion 8 using the software RHMAP v.3.0, consisted of three linkage groups 23, 22, and 590 cR5000 in length. Gene assignments on BTA27 and the localization of 8 more genes on BTA8 and BTA14 previously predicted on BTA8/BTA27 and BTA14/BTA27 narrowed down significantly the chromosome break points between the three cattle chromosomes and segments on Homo sapiens chromosomes HSA4 and HSA8. Defined evolutionary break points increase the accuracy of comparative in silico mapping of further human genes in conserved chromosome segments of BTA27.

A 1463 gene cattle-human comparative map with anchor points defined by human genome sequence coordinates

A second-generation 5000 rad radiation hybrid (RH) map of the cattle genome was constructed primarily using cattle ESTs that were targeted to gaps in the existing cattle-human comparative map, as well as to sparsely populated map intervals. A total of 870 targeted markers were added, bringing the number of markers mapped on the RH(5000) panel to 1913. Of these, 1463 have significant BLASTN hits (E < e(-5)) against the human genome sequence. A cattle-human comparative map was created using human genome sequence coordinates of the paired orthologs. One-hundred and ninety-five conserved segments (defined by two or more genes) were identified between the cattle and human genomes, of which 31 are newly discovered and 34 were extended singletons on the first-generation map. The new map represents an improvement of 20% genome-wide comparative coverage compared with the first-generation map. Analysis of gene content within human genome regions where there are gaps in the comparative map revealed gaps with both significantly greater and significantly lower gene content. The new, more detailed cattle-human comparative map provides an improved resource for the analysis of mammalian chromosome evolution, the identification of candidate genes for economically important traits, and for proper alignment of sequence contigs on cattle chromosomes.

A graph-theoretic approach to comparing and integrating genetic, physical and sequence-based maps

For many species, multiple maps are available, often constructed independently by different research groups using different sets of markers and different source material. Integration of these maps provides a higher density of markers and greater genome coverage than is possible using a single study. In this article, we describe a novel approach to comparing and integrating maps by using abstract graphs. A map is modeled as a directed graph in which nodes represent mapped markers and edges define the order of adjacent markers. Independently constructed graphs representing corresponding maps from different studies are merged on the basis of their common loci. Absence of a path between two nodes indicates that their order is undetermined. A cycle indicates inconsistency among the mapping studies with regard to the order of the loci involved. The integrated graph thus produced represents a complete picture of all of the mapping studies that comprise it, including all of the ambiguities and inconsistencies among them. The objective of this representation is to guide additional research aimed at interpreting these ambiguities and inconsistencies in locus order rather than presenting a "consensus order" that ignores these problems.

The bovine 5' AMPK gene family: mapping and single nucleotide polymorphism detection

The 5'-AMP-activated protein kinase (AMPK) family is an ancient stress response system whose primary function is regulation of cellular ATP. Activation of AMPK, which is instigated by environmental and nutritional stresses, initiates energy-conserving measures that protect the cell by inhibition and phosphorylation of key enzymes in energy-consuming biochemical pathways. The seven genes that comprise the bovine AMPK family were mapped in cattle by using a radiation hybrid panel. The seven genes mapped to six different cattle chromosomes, each with a LOD score greater than 10.0. PRKAA1 mapped to BTA 20, PRKAA2 and PRKAB2 to BTA 3, PRKAB1 to BTA 17, PRKAG1 to BTA 5, PRKAG2 to BTA 4, and PRKAG3 to BTA 2. Five of the seven genes mapped to regions expected from human/cattle comparative maps. PRKAB2 and PRKAG3, however, have not been mapped in humans. We predict these genes to be located on HSA 1 and 2, respectively. Additionally, one synonymous and one non-synonymous single nucleotide polymorphism (SNP) were detected in PRKAG3 in Bos taurus cattle. In an effort to determine ancestral origins, various herds of mixed breed cattle as well as other ruminant species were characterized for sequence variation in this region of PRKAG3. Owing to the physiological importance of this gene family, we believe that its individual genes are candidate genes for conferring resistance to diseases in cattle.

Mastitis increases mammary mRNA abundance of beta-defensin 5, toll-like-receptor 2 (TLR2), and TLR4 but not TLR9 in cattle

Coordination of the primary defense mechanisms against pathogens relies on the appropriate expression of pathogen recognition receptors (PRRs) triggering the early release of effector molecules of the innate immune system. To analyze the impact of this system on the counteraction of infections of the mammary gland (mastitis), we characterized the bovine gene encoding the key PRR Toll-like receptor 9 (TLR9) and mapped its precise position on chromosome BTA22. The sequence information was used to establish real-time PCR quantification assays to measure the mRNA abundances of TLR9, TLR2, and TLR4 together with those of beta-defensin 5 (BNBD5), an early bactericidal effector molecule of the innate system, in healthy and infected mammary glands. Mastitis strongly increased (4- to 13-fold) the mRNA abundances of all of these genes except TLR9. Slight subclinical infections already caused a substantial increase in the copy numbers, though they did so the least for TLR9. Induction was not systemic, since mRNA abundance was low in uninfected control quarters of the udder but high in the severely infected quarters of the same animal. The number of TLR2 copies correlated well with those of TLR4, indicating coordinated regulation of these two PRRs during infection of the udder. Their coordinated regulation explains our unexpected observation that pure Staphylococcus aureus infections caused a strong increase also in TLR4 mRNA abundance. In situ hybridizations revealed that BNBD5 is expressed predominantly in the mammary epithelial cells (MEC) of the infected gland. Our data therefore suggest a significant contribution of the innate immune system to counteract mastitis and attribute a prominent effector function to the MEC.