Comparative analysis of the river buffalo somatostatin gene.

The somatostatin protein plays a crucial role in the regulation of multiple biological functions, such as growth, fat deposition, and nutrient absorption in vertebrates. Polymorphisms in the somatostatin gene have been associated with growth traits in livestock species, including cattle and goat. In this study, we conducted complete molecular characterization of the somatostatin gene in Bubalus bubalis (Murrah breed) by sequencing a Murrah BAC clone spanning 72,489 base pairs (bp) in length. The buffalo somatostatin gene contains 1481 bp organized into a 5'-untranslated region (135 bp), exon 1 (139 bp), intron 1 (839 bp), exon 2 (212 bp), and 3'UTR (156 bp). Comparative analysis between the buffalo somatostatin DNA coding sequence and the amino acid sequence with other bovids (cattle, goat, and sheep), horse, pig, human, rodents (mouse and rat), and chicken. Identity varied from 83-99% on the DNA sequence level and 88-100% on the protein level. In addition, a comparison of gene sequences between Murrah and Mediterranean breeds revealed 6 potential single-nucleotide polymorphisms (1 in exon 1 and 5 in intron 1), which were validated in different buffalo populations. This comparative analysis provides basic information for future studies of different buffalo herds using the position candidate gene approach, quantitative trait loci analysis, and polymorphisms associated with growth traits.

A high-resolution radiation hybrid map of the river buffalo major histocompatibility complex and comparison with BoLA.

The major histocompatibility complex (MHC) in mammals codes for antigen-presenting proteins. For this reason, the MHC is of great importance for immune function and animal health. Previous studies revealed this gene-dense and polymorphic region in river buffalo to be on the short arm of chromosome 2, which is homologous to cattle chromosome 23. Using cattle-derived STS markers and a river buffalo radiation hybrid (RH) panel (BBURH5000 ), we generated a high-resolution RH map of the river buffalo MHC region. The buffalo MHC RH map (cR5000 ) was aligned with the cattle MHC RH map (cR12000 ) to compare gene order. The buffalo MHC had similar organization to the cattle MHC, with class II genes distributed in two segments, class IIa and class IIb. Class IIa was closely associated with the class I and class III regions, and class IIb was a separate cluster. A total of 53 markers were distributed into two linkage groups based on a two-point LOD score threshold of ≥8. The first linkage group included 32 markers from class IIa, class I and class III. The second linkage group included 21 markers from class IIb. Bacterial artificial chromosome clones for seven loci were mapped by fluorescence in situ hybridization on metaphase chromosomes using single- and double-color hybridizations. The order of cytogenetically mapped markers in the region corroborated the physical order of markers obtained from the RH map and served as anchor points to align and orient the linkage groups.

Construction and preliminary characterization of a river buffalo bacterial artificial chromosome library.

River buffalo genome analyses have advanced significantly in the last decade, and the genome sequence of Bubalus bubalis will be available shortly. Nonetheless, large-insert DNA library resources such as bacterial artificial chromosomes (BAC) are still required for validation and accurate assembly of the genome sequence. We constructed a river buffalo BAC library containing 52,224 clones with an average insert size of 97 kb, representing 1.7 × coverage of the genome. This genomic resource for river buffalo will facilitate further studies in this economically important species allowing for instance, whole genome physical mapping and isolation of genes and gene clusters, contributing to the elucidation of gene organization and identification of regulatory elements.

Comparative RH maps of the river buffalo and bovine Y chromosomes.

Radiation hybrid maps were constructed for river buffalo and cattle Y chromosomes. A total of 41 cattle-derived Y-chromosome molecular markers were selected and tested with 2 previously described 5,000-rad whole-genome radiation hybrid (RH) panels (river buffalo - BBURH(5000) and cattle - BTARH(5000)) for generation of maps. Among the initial 41 selected markers, a subset of 26 markers generated PCR products suitable for scoring with the BBURH(5000) panel. Of these, 19 markers (73%) were distributed in 1 linkage group spanning 341.3 cR. Retention frequencies (RF) for individual markers ranged from 17.8% for SMCY to 56.7% for BTY1, with an average RF of 37.6%. From the selected markers, 37 generated reliable scores using the BTARH(5000) panel. The newly constructed BTAY RH map contains 28 markers distributed within 1 linkage group. Twenty-four of these markers had been previously mapped on BTAY using a 7,000-rad cattle-hamster WG-RH panel and 4 markers were mapped for the first time (ZFY, SeqRep, RepSeqS4 and BTY1). The length of the BTAY RH map was estimated to be 602.4 cR. Retention frequencies for individual mapped markers ranged from 10% (INRA126) to 63.3% (SeqRep), with an average RF of 35.3%. RH marker positions along the Y chromosome were compared between BBUY and BTAY, which revealed differences in the order of some of the markers. The BBUY pseudoautosomal region (PAR) is delineated by 3 BTAY PAR markers (MAF45, TGLA325 and UMN2008). These markers are telomeric in both species but are not found in the same order. Here we have demonstrated the effective use of bovine Y chromosome markers for the development of the first BBUY RH map. Likewise, these set of markers can be used for comparative assessment of Y chromosomes in other members of the Bovidae family.

Mapping MHC genes in river buffalo.

The major histocompatibility complex (MHC) contains a set of genes necessary for antigen presentation in the immune system. This gene dense and polymorphic region of the mammalian genome is of considerable interest due to the role of MHC genes in immune function and animal health. Previous cytogenetic studies have indicated that the MHC in river buffalo resides on the short arm of chromosome 2 (BBU2). A 5000-rad radiation hybrid mapping panel was recently generated to enable construction of a whole genome map of river buffalo. To this end, the aims of this project were to elucidate the general organization of the MHC on BBU2, and to compare gene order within this region to the MHC in cattle. PCR primers were selected from the bovine gene map and used with the BBURH5000 panel to map a set of ten MHC class II genes in river buffalo. Analysis indicates that these genes fall into two linkage groups, consistent with organization of the MHC in cattle. This comparison of buffalo and bovine MHC gene order provides the first insight into the organization of the MHC on river buffalo chromosome 2.

Mapping the BDA20 gene to BTAX.

The lipocalin family is a large group of proteins that exhibits great structural and functional variation both within and among species, including a significant number of animal-derived aeroallergens, such as the bovine BDA20 (major cow dander allergen). This protein is classified as an occupational allergen causing asthma and other work-related allergic disorders among dairy farmers. Using a somatic cell panel the BDA20 gene was assigned to the bovine X chromosome (BTAX) with a significant concordant value of 97% to the previously mapped reference marker MAF45. A radiation hybrid (RH) mapping approach confirmed the assignment of BDA20 to BTAX. Two-point LOD scores showed that BDA20 is linked to XBM451 with a LOD score of 22.1 for a theta value of 0.03.

First radiation hybrid map of the river buffalo X chromosome (BBUX) and comparison with BTAX.

We report the first radiation hybrid map of the river buffalo X chromosome generated from a recently constructed river buffalo (Bubalus bubalis) whole-genome radiation hybrid panel (BBURH(5000)). This map contains a total of 33 cattle-derived markers, including 10 genes, four ESTs and 19 microsatellites. The markers are distributed in two linkage groups: LG1 contains eight markers spanning 125.6 cR, and LG2 contains 25 markers spanning 366.3 cR. LG1 contains six markers in common with bovine sequence assembly build 3.1. With the exception of BMS2152, the order of these markers on our BBUX map is shuffled when compared to the cow X chromosome (Bos taurus; BTAX). From LG2, two markers (AMELX and BL22) map to a more distal portion of BTAX compared to BBUX. In addition, two pairs of LG2 markers exhibit inversions compared to BTAX (ILSTS017 and ATRX; XBM38 and PPEF1). Alternatively, when compared to the most recent bovine RH map (Bov-Gen 3000rads), BL1098 and BMS2227 from LG1 as well as PLS3 and BMS1820 from LG2 showed inverted positions on the BBUX map. These discrepancies in buffalo and cattle maps may reflect evolutionary divergence of the chromosomes or mapping errors in one of the two species. Although the set of mapped markers does not cover the entire X chromosome, this map is a starting point for the construction of a high-resolution map, which is necessary for characterization of small rearrangements that might have occurred between the Bubalus bubalis and Bos taurus X chromosomes.

Report of a chimeric origin of transposable elements in a bovine-coding gene.

Despite the wide distribution of transposable elements (TEs) in mammalian genomes, part of their evolutionary significance remains to be discovered. Today there is a substantial amount of evidence showing that TEs are involved in the generation of new exons in different species. In the present study, we searched 22,805 genes and reported the occurrence of TE-cassettes in coding sequences of 542 cow genes using the RepeatMasker program. Despite the significant number (542) of genes with TE insertions in exons only 14 (2.6%) of them were translated into protein, which we characterized as chimeric genes. From these chimeric genes, only the FAST kinase domains 3 (FASTKD3) gene, present on chromosome BTA 20, is a functional gene and showed evidence of the exaptation event. The genome sequence analysis showed that the last exon coding sequence of bovine FASTKD3 is approximately 85% similar to the ART2A retrotransposon sequence. In addition, comparison among FASTKD3 proteins shows that the last exon is very divergent from those of Homo sapiens, Pan troglodytes and Canis familiares. We suggest that the gene structure of bovine FASTKD3 gene could have originated by several ectopic recombinations between TE copies. Additionally, the absence of TE sequences in all other species analyzed suggests that the TE insertion is clade-specific, mainly in the ruminant lineage.

Report of a chimeric origin of transposable elements in a bovine-coding gene

Despite the wide distribution of transposable elements (TEs) in mammalian genomes, part of their evolutionary significance remains to be discovered. Today there is a substantial amount of evidence showing that TEs are involved in the generation of new exons in different species. In the present study, we searched 22,805 genes and reported the occurrence of TE-cassettes in coding sequences of 542 cow genes using the RepeatMasker program. Despite the significant number (542) of genes with TE insertions in exons only 14 (2.6%) of them were translated into protein, which we characterized as chimeric genes. From these chimeric genes, only the FAST kinase domains 3 (FASTKD3) gene, present on chromosome BTA 20, is a functional gene and showed evidence of the exaptation event. The genome sequence analysis showed that the last exon coding sequence of bovine FASTKD3 is ~85% similar to the ART2A retrotransposon sequence. In addition, comparison among FASTKD3 proteins shows that the last exon is very divergent from those of Homo sapiens, Pan troglodytes and Canis familiares. We suggest that the gene structure of bovine FASTKD3 gene could have originated by several ectopic recombinations between TE copies. Additionally, the absence of TE sequences in all other species analyzed suggests that the TE insertion is clade-specific, mainly in the ruminant lineage.

A radiation hybrid map of river buffalo (Bubalus bubalis) chromosome 1 (BBU1).

The largest chromosome in the river buffalo karyotype, BBU1, is a submetacentric chromosome with reported homology between BBU1q and bovine chromosome 1 and between BBU1p and BTA27. We present the first radiation hybrid map of this chromosome containing 69 cattle derived markers including 48 coding genes, 17 microsatellites and four ESTs distributed in two linkage groups spanning a total length of 1330.1 cR(5000). The RH map was constructed based on analysis of a recently developed river buffalo-hamster whole genome radiation hybrid (BBURH(5000)) panel. The retention frequency of individual markers across the panel ranged from 17.8 to 52.2%. With few exceptions, the order of markers within linkage groups is identical to the order established for corresponding cattle RH maps. The BBU1 map provides a starting point for comparison of gene order rearrangements between river buffalo chromosome 1 and its bovine homologs.

Preliminary radiation hybrid map for river buffalo chromosome 6 and comparison to bovine chromosome 3.

We present the first radiation hybrid (RH) map of river buffalo (Bubalus bubalis) chromosome 6 (BBU6) developed with a recently constructed river buffalo whole-genome RH panel (BBURH(5000)). The preliminary map contains 33 cattle-derived markers, including 12 microsatellites, 19 coding genes and two ESTs, distributed across two linkage groups. Retention frequencies for markers ranged from 14.4% to 40.0%. Most of the marker orders within the linkage groups on BBU6 were consistent with the cattle genome sequence and RH maps. This preliminary RH map is the starting point for comparing gene order between river buffalo and cattle, presenting an opportunity for the examination of micro-rearrangements of these chromosomes. Also, resources for positional candidate cloning in river buffalo are enhanced.

Construction of a river buffalo (Bubalus bubalis) whole-genome radiation hybrid panel and preliminary RH mapping of chromosomes 3 and 10.

The buffalo (Bubalus bubalis) is a source of milk and meat, and also serves as a draft animal. In this study, a 5000-rad whole-genome radiation hybrid (RH) panel for river buffalo was constructed and used to build preliminary RH maps for BBU3 and BBU10 chromosomes. The preliminary maps contain 66 markers, including coding genes, cattle expressed sequence tags (ESTs) and microsatellite loci. The RH maps presented here are the starting point for mapping additional loci that will allow detailed comparative maps between buffalo, cattle and other species whose genomes may be mapped in the future. A large quantity of DNA has been prepared from the cell lines forming the river buffalo RH panel and will be made publicly available to the international community both for the study of chromosome evolution and for the improvement of traits important to the role of buffalo in animal agriculture.

A radiation hybrid map of river buffalo (Bubalus bubalis) chromosome 7 and comparative mapping to the cattle and human genomes.

A preliminary radiation hybrid (RH) map containing 50 loci on chromosome 7 of the domestic river buffalo Bubalus bubalis (BBU; 2n = 50) was constructed based on a comparative mapping approach. The RH map of BBU7 includes thirty-seven gene markers and thirteen microsatellites. All loci have been previously assigned to Bos taurus (BTA) chromosome BTA6, which is known for its association with several economically important milk production traits in cattle. The map consists of two linkage groups spanning a total length of 627.9 cR(5,000). Comparative analysis of the BBU7 RH(5,000) map with BTA6 in cattle gave new evidence for strong similarity between the two chromosomes over their entire length and exposed minor differences in locus order. Comparison of the BBU7 RH(5,000) map with the Homo sapiens (HSA) genome revealed similarity with a large chromosome segment of HSA4. Comparative analysis of loci in both species revealed more variability than previously known in gene order and several chromosome rearrangements including centromere relocation. The data obtained in our study define the evolutionarily conserved segment on BBU7 and HSA4 to be between 3.5 megabases (Mb) and 115.8 Mb in the HSA4 (genome build 36) DNA sequence.

Clarifications on breakpoints in HSAX and BTAX by comparative mapping of F9, HPRT, and XIST in cattle.

The coagulation factor IX gene (F9), the hypoxanthine phosphoribosyl transferase 1 gene (HPRT1), and the X-inactive specific transcript gene (XIST) were physically assigned in cattle to analyze chromosomal breakpoints on BTAX recently identified by radiation hybrid (RH) mapping experiments. Whereas the FISH assignment of XIST indicates a similar location on the q-arm of the human and cattle X chromosomes, the locus of HPRT1 supported the assumption of a chromosome rearrangement between the distal half of the q-arm of HSAX and the p-arm of BTAX identified by RH mapping. F9 previously located on the q-arm of BTAX was assigned to the p-arm of BTAX using RH mapping and FISH. The suggested new position of F9 close to HPRT1 supports the homology between HSAXq and BTAXp. The F9 locus corresponds with the gene order found in the homologous human chromosome segment. XIST was assigned on BTAXq23, HPRT1 and F9 were mapped to BTAXp22, and the verification of the location of F9 in a 5000 rad cattle-hamster whole genome radiation hybrid panel linked the gene to markers URB10 and HPRT1.