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.

Contiguous ∼16 Mb 1p36 deletion: Dominant features of classical distal 1p36 monosomy with haplo-lethality

Monosomy 1p36 results from heterozygous deletions of the terminal short chromosome 1 arm, the most common terminal deletion in humans. The microdeletion is split in two usually non-overlapping and clinically distinct classical distal and proximal 1p36 monosomy syndromes. Using comparative genome hybridization, MLPA and qPCR we identified the largest contiguous ∼16 Mb terminal 1p36 deletion reported to date. It covers both distal and proximal regions, causes a neonatally lethal variant with virtually exclusive features of distal 1p36 monosomy, highlighting the key importance of the gene-rich distal region for the "compound" 1p36 phenotype and a threshold deletion-size effect for haplo-lethality.

Molecular characterization and chromosomal assignment of the bovine glycinamide ribonucleotide formyltransferase (GART) gene on cattle chromosome 1q12.1-q12.2

The mammalian glycinamide ribonucleotide formyltransferase (GART) genes encode a trifunctional polypeptide involved in the de novo purine biosynthesis. We isolated a bacterial artificial chromosome (BAC) clone containing the bovine GART gene and determined the complete DNA sequence of the BAC clone. Cloning and characterization of the bovine GART gene revealed that the bovine gene consists of 23 exons spanning approximately 27 kb. RT-PCR amplification of bovine GART in different organs showed the expression of two GART transcripts in cattle similar to human and mouse. The GART transcripts encode two proteins of 1010 and 433 amino acids, respectively. Eleven single nucleotide polymorphisms (SNPs) were detected in a mutation scan of 24 unrelated animals of three different cattle breeds, including one SNP that affects the amino acid sequence of GART. The chromosomal localization of the gene was determined by fluorescence in situ hybridization. Comparative genome analysis between cattle, human and mouse indicates that the chromosomal location of the bovine GART gene is in agreement with a previously published mapping report.

A 4,103 marker integrated physical and comparative map of the horse genome

A comprehensive second-generation whole genome radiation hybrid (RH II), cytogenetic and comparative map of the horse genome (2n = 64) has been developed using the 5000rad horse x hamster radiation hybrid panel and fluorescence in situ hybridization (FISH). The map contains 4,103 markers (3,816 RH; 1,144 FISH) assigned to all 31 pairs of autosomes and the X chromosome. The RH maps of individual chromosomes are anchored and oriented using 857 cytogenetic markers. The overall resolution of the map is one marker per 775 kilobase pairs (kb), which represents a more than five-fold improvement over the first-generation map. The RH II incorporates 920 markers shared jointly with the two recently reported meiotic maps. Consequently the two maps were aligned with the RH II maps of individual autosomes and the X chromosome. Additionally, a comparative map of the horse genome was generated by connecting 1,904 loci on the horse map with genome sequences available for eight diverse vertebrates to highlight regions of evolutionarily conserved syntenies, linkages, and chromosomal breakpoints. The integrated map thus obtained presents the most comprehensive information on the physical and comparative organization of the equine genome and will assist future assemblies of whole genome BAC fingerprint maps and the genome sequence. It will also serve as a tool to identify genes governing health, disease and performance traits in horses and assist us in understanding the evolution of the equine genome in relation to other species.

The mutation causing the black-and-tan pigmentation phenotype of Mangalitza pigs maps to the porcine ASIP locus but does not affect its coding sequence

The gene for agouti signaling protein (ASIP) is centrally involved in the expression of coat color traits in animals. The Mangalitza pig breed is characterized by a black-and-tan phenotype with black dorsal pigmentation and yellow or white ventral pigmentation. We investigated a Mangalitza x Piétrain cross and observed a coat color segregation pattern in the F2 generation that can be explained by virtue of two alleles at the MC1R locus and two alleles at the ASIP locus. Complete linkage of the black-and-tan phenotype to microsatellite alleles at the ASIP locus on SSC 17q21 was observed. Corroborated by the knowledge of similar mouse coat color mutants, it seems therefore conceivable that the black-and-tan pigmentation of Mangalitza pigs is caused by an ASIP allele a(t), which is recessive to the wild-type allele A. Toward positional cloning of the a(t) mutation, a 200-kb genomic BAC/PAC contig of this chromosomal region has been constructed and subsequently sequenced. Full-length ASIP cDNAs obtained by RACE differed in their 5' untranslated regions, whereas they shared a common open reading frame. Comparative sequencing of all ASIP exons and ASIP cDNAs between Mangalitza and Piétrain pigs did not reveal any differences associated with the coat color phenotype. Relative qRT-PCR analyses showed different dorsoventral skin expression intensities of the five ASIP transcripts in black-and-tan Mangalitza. The a(t) mutation is therefore probably a regulatory ASIP mutation that alters its dorsoventral expression pattern.

Genome sequence, comparative analysis, and population genetics of the domestic horse

We report a high-quality draft sequence of the genome of the horse (Equus caballus). The genome is relatively repetitive but has little segmental duplication. Chromosomes appear to have undergone few historical rearrangements: 53% of equine chromosomes show conserved synteny to a single human chromosome. Equine chromosome 11 is shown to have an evolutionary new centromere devoid of centromeric satellite DNA, suggesting that centromeric function may arise before satellite repeat accumulation. Linkage disequilibrium, showing the influences of early domestication of large herds of female horses, is intermediate in length between dog and human, and there is long-range haplotype sharing among breeds.

A comparative radiation hybrid map of sheep chromosome 10

Comparative radiation hybrid (RH) maps of individual ovine chromosomes are essential to identify genes governing traits of economic importance in sheep, a livestock species for which whole genome sequence data are not yet available. The USUoRH5000 radiation hybrid panel was used to generate a RH map of sheep chromosome 10 (OAR10) with 59 markers that span 1,422 cR over an estimated 92 Mb of the chromosome, thus providing markers every 2 Mb (equivalent to every 24 cR). The markers were derived from 46 BAC end sequences (BESs), a single EST, and 12 microsatellites. Comparative analysis showed that OAR10 shares remarkable conservation of gene order along the entire length of cattle chromosome 12 and that OAR10 contains four major homologous synteny blocks, each related to segments of the homologous human chromosome 13. Extending the comparison to the horse, dog, mouse, and chicken genome showed that these blocks share conserved synteny across species.

A 4 Mb high resolution BAC contig on bovine chromosome 1q12 and comparative analysis with human chromosome 21q22

The bovine RPCI-42 BAC library was screened to construct a sequence-ready ~4 Mb single contig of 92 BAC clones on BTA 1q12. The contig covers the region between the genes KRTAP8P1 and CLIC6. This genomic segment in cattle is of special interest as it contains the dominant gene responsible for the hornless or polled phenotype in cattle. The construction of the BAC contig was initiated by screening the bovine BAC library with heterologous cDNA probes derived from 12 human genes of the syntenic region on HSA 21q22. Contig building was facilitated by BAC end sequencing and chromosome walking. During the construction of the contig, 165 BAC end sequences and 109 single-copy STS markers were generated. For comparative mapping of 25 HSA 21q22 genes, genomic PCR primers were designed from bovine EST sequences and the gene-associated STSs mapped on the contig. Furthermore, bovine BAC end sequence comparisons against the human genome sequence revealed significant matches to HSA 21q22 and allowed the in silico mapping of two new genes in cattle. In total, 31 orthologues of human genes located on HSA 21q22 were directly mapped within the bovine BAC contig, of which 16 genes have been cloned and mapped for the first time in cattle. In contrast to the existing comparative bovine-human RH maps of this region, these results provide a better alignment and reveal a completely conserved gene order in this 4 Mb segment between cattle, human and mouse. The mapping of known polled linked BTA 1q12 microsatellite markers allowed the integration of the physical contig map with existing linkage maps of this region and also determined the exact order of these markers for the first time. Our physical map and transcript map may be useful for positional cloning of the putative polled gene in cattle.

Comparative expression profiling of E. coli and S. aureus inoculated primary mammary gland cells sampled from cows with different genetic predispositions for somatic cell score

BACKGROUND: During the past ten years many quantitative trait loci (QTL) affecting mastitis incidence and mastitis related traits like somatic cell score (SCS) were identified in cattle. However, little is known about the molecular architecture of QTL affecting mastitis susceptibility and the underlying physiological mechanisms and genes causing mastitis susceptibility. Here, a genome-wide expression analysis was conducted to analyze molecular mechanisms of mastitis susceptibility that are affected by a specific QTL for SCS on Bos taurus autosome 18 (BTA18). Thereby, some first insights were sought into the genetically determined mechanisms of mammary gland epithelial cells influencing the course of infection. METHODS: Primary bovine mammary gland epithelial cells (pbMEC) were sampled from the udder parenchyma of cows selected for high and low mastitis susceptibility by applying a marker-assisted selection strategy considering QTL and molecular marker information of a confirmed QTL for SCS in the telomeric region of BTA18. The cells were cultured and subsequently inoculated with heat-inactivated mastitis pathogens Escherichia coli and Staphylococcus aureus, respectively. After 1, 6 and 24 h, the cells were harvested and analyzed using the microarray expression chip technology to identify differences in mRNA expression profiles attributed to genetic predisposition, inoculation and cell culture. RESULTS: Comparative analysis of co-expression profiles clearly showed a faster and stronger response after pathogen challenge in pbMEC from less susceptible animals that inherited the favorable QTL allele 'Q' than in pbMEC from more susceptible animals that inherited the unfavorable QTL allele 'q'. Furthermore, the results highlighted RELB as a functional and positional candidate gene and related non-canonical Nf-kappaB signaling as a functional mechanism affected by the QTL. However, in both groups, inoculation resulted in up-regulation of genes associated with the Ingenuity pathways 'dendritic cell maturation' and 'acute phase response signaling', whereas cell culture affected biological processes involved in 'cellular development'. CONCLUSIONS: The results indicate that the complex expression profiling of pathogen challenged pbMEC sampled from cows inheriting alternative QTL alleles is suitable to study genetically determined molecular mechanisms of mastitis susceptibility in mammary epithelial cells in vitro and to highlight the most likely functional pathways and candidate genes underlying the QTL effect.

LUPA: a European initiative taking advantage of the canine genome architecture for unravelling complex disorders in both human and dogs

The domestic dog offers a unique opportunity to explore the genetic basis of disease, morphology and behaviour. Humans share many diseases with our canine companions, making dogs an ideal model organism for comparative disease genetics. Using newly developed resources, genome-wide association studies in dog breeds are proving to be exceptionally powerful. Towards this aim, veterinarians and geneticists from 12 European countries are collaborating to collect and analyse the DNA from large cohorts of dogs suffering from a range of carefully defined diseases of relevance to human health. This project, named LUPA, has already delivered considerable results. The consortium has collaborated to develop a new high density single nucleotide polymorphism (SNP) array. Mutations for four monogenic diseases have been identified and the information has been utilised to find mutations in human patients. Several complex diseases have been mapped and fine mapping is underway. These findings should ultimately lead to a better understanding of the molecular mechanisms underlying complex diseases in both humans and their best friend.

Comparative human-horse sequence analysis of the CYP3A subfamily gene cluster

Cytochrome P450 enzymes (CYP450s) represent a superfamily of haem-thiolate proteins. CYP450s are most abundant in the liver, a major site of drug metabolism, and play key roles in the metabolism of a variety of substrates, including drugs and environmental contaminants. Interaction of two or more different drugs with the same enzyme can account for adverse effects and failure of therapy. Human CYP3A4 metabolizes about 50% of all known drugs, but little is known about the orthologous CYP450s in horses. We report here the genomic organization of the equine CYP3A gene cluster as well as a comparative analysis with the human CYP3A gene cluster. The equine CYP450 genes of the 3A family are located on ECA 13 between 6.97-7.53 Mb, in a region syntenic to HSA 7 99.05-99.35 Mb. Seven potential, closely linked equine CYP3A genes were found, in contrast to only four genes in the human genome. RNA was isolated from an equine liver sample, and the approximately 1.5-kb coding sequence of six CYP3A genes could be amplified by RT-PCR. Sequencing of the RT-PCR products revealed numerous hitherto unknown single nucleotide polymorphisms (SNPs) in these six CYP3A genes, and one 6-bp deletion compared to the reference sequence (EquCab2.0). The presence of the variants was confirmed in a sample of genomic DNA from the same horse. In conclusion, orthologous genes for the CYP3A family exist in horses, but their number differs from those of the human CYP3A gene family. CYP450 genes of the same family show high homology within and between mammalian species, but can be highly polymorphic.

Genome-Wide Survey of Pseudogenes in 80 Fully Re-sequenced Arabidopsis thaliana Accessions

Pseudogenes (Ψs), including processed and non-processed Ψs, are ubiquitous genetic elements derived from originally functional genes in all studied genomes within the three kingdoms of life. However, systematic surveys of non-processed Ψs utilizing genomic information from multiple samples within a species are still rare. Here a systematic comparative analysis was conducted of Ψs within 80 fully re-sequenced Arabidopsis thaliana accessions, and 7546 genes, representing ~28% of the genomic annotated open reading frames (ORFs), were found with disruptive mutations in at least one accession. The distribution of these Ψs on chromosomes showed a significantly negative correlation between Ψs/ORFs and their local gene densities, suggesting a higher proportion of Ψs in gene desert regions, e.g. near centromeres. On the other hand, compared with the non-Ψ loci, even the intact coding sequences (CDSs) in the Ψ loci were found to have shorter CDS length, fewer exon number and lower GC content. In addition, a significant functional bias against the null hypothesis was detected in the Ψs mainly involved in responses to environmental stimuli and biotic stress as reported, suggesting that they are likely important for adaptive evolution to rapidly changing environments by pseudogenization to accumulate successive mutations.

Phylogenetic analysis of the gag region encoding the matrix protein of small ruminant lentiviruses: comparative analysis and molecular epidemiological applications

Little sequence information exists on the matrix-protein (MA) encoding region of small ruminant lentiviruses (SRLV). Fifty-two novel sequences were established and permitted a first phylogenetic analysis of this region of the SRLV genome. The variability of the MA encoding region is higher compared to the gag region encoding the capsid protein and surprisingly close to that reported for the env gene. In contrast to primate lentiviruses, the deduced amino acid sequences of the N- and C-terminal domains of MA are variable. This permitted to pinpoint a basic domain in the N-terminal domain that is conserved in all lentiviruses and likely to play an important functional role. Additionally, a seven amino acid insertion was detected in all MVV strains, which may be used to differentiate CAEV and MVV isolates. A molecular epidemiology analysis based on these sequences indicates that the Italian lentivirus strains are closely related to each other and to the CAEV-CO strain, a prototypic strain isolated three decades ago in the US. This suggests a common origin of the SRLV circulating in the monitored flocks, possibly related to the introduction of infected goats in a negative population. Finally, this study shows that the MA region is suitable for phylogenetic studies and may be applied to monitor SRLV eradication programs.

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 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.