Search

We present a simple algorithm for reconstruction of haplotypes from a sample of multilocus genotypes. The algorithm is aimed specifically for analysis of very large pedigrees for small chromosomal segments, where recombination frequency within the chromosomal segment can be assumed to be zero. The algorithm was tested both on simulated pedigrees of 155 individuals in a family structure of three generations and on real data of 1149 animals from the Israeli Holstein dairy cattle population, including 406 bulls with genotypes, but no females with genotypes. The rate of haplotype resolution for the simulated data was >91% with a standard deviation of 2%. With 20% missing data, the rate of haplotype resolution was 67.5% with a standard deviation of 1.3%. In both cases all recovered haplotypes were correct. In the real data, allele origin was resolved for 22% of the heterozygous genotypes, even though 70% of the genotypes were missing. Haplotypes were resolved for 36% of the males. Computing time was insignificant for both data sets. Despite the intricacy of large-scale real pedigree genotypes, the proposed algorithm provides a practical rule-based solution for resolving haplotypes for small chromosomal segments in commercial animal populations.

We previously localized a quantitative trait locus (QTL) on chromosome 6 affecting milk fat and protein concentration to a 4-cM confidence interval, centered on the microsatellite BM143. We characterized the genes and sequence variation in this region and identified common haplotypes spanning five polymorphic sites in the genes IBSP, SPP1, PKD2, and ABCG2 for two sires heterozygous for this QTL. Expression of SPP1 and ABCG2 in the bovine mammary gland increased from parturition through lactation. SPP1 and all the coding exons of ABCG2 and PKD2 were sequenced for these two sires. The single nucleotide change capable of encoding a substitution of tyrosine-581 to serine (Y581S) in the ABCG2 transporter was the only polymorphism corresponding to the segregation status of all 3 heterozygous and 15 homozygous sires for the QTL in the Israeli and U.S. Holstein populations. The allele substitution fixed effects on the genetic evaluations of 335 Israeli sires were -341 kg milk, +0.16% fat, and +0.13% protein (F-value = 200). No other polymorphism gave significant effect for fat and protein concentration in models that also included Y581S. The allele substitution effects on the genetic evaluations of 670 cows, daughters of two heterozygous sires, were -226 kg milk, 0.09% fat, and 0.08% protein (F-value = 394), with partial dominance towards the 581S homozygotes. We therefore propose that Y581S in ABCG2 is the causative site for this QTL.

A cluster of genes coding for proteins of the extracellular matrix (ECM) containing sequence motifs essential for integrin-receptor interactions is located on HSA4q21 and on BTA6, within the critical region of a quantitative trait locus (QTL) affecting milk protein production. Genes within this cluster are involved in the formation of bone and lobuloalveolar structures in mammary gland and in kidney function. We cloned a bovine gene neighboring this ECM cluster, termed FAM13A1, the first member of a novel gene family (FAM13). A short predominant 5.1-kb mRNA variant capable of encoding 697 amino acids is transcribed from 18-exon orthologous genes in human and mouse. All putative protein orthologs contained a bipartite nuclear-localization signal and two coil-coiled domains. We detected two other FAM13 paralogs, C5ORF5 and FAM13C1, on HSA5q31 and HSA10q21, respectively. All FAM13 paralogs produce transcripts that are complementary to adjacent genes, suggesting that antisense transcription may regulate their functions. The structure of a longer 5.9-kb variant, FAM13A1_v2, that has an extra 7-exon putative RhoGAP domain at the N-terminus provides further clues to FAM13 function. Analysis of 400 bulls revealed a population-wide linkage disequilibrium between FAM13A1 polymorphisms and the QTL.

Motivation: Sequencing of a bi-allelic PCR product, which contains an allele with a deletion/insertion mutation results in a superimposed tracefile following the site of this shift mutation. A trace file of this type hampers the use of current computer programs for base calling. ShiftDetector analyses a sequencing trace file in order to discover if it is a superimposed sequence of two molecules that differ in a shift mutation of 1 to 25 bases. The program calculates a probability score for the existence of such a shift and reconstructs the sequence of the original molecule., Availability: ShiftDetector is available from http://cowry.agri.huji.ac.il

Nine Israeli Holstein sire families with 2978 daughters were analyzed for quantitative trait loci effects on chromosome 6 for five milk production traits by a daughter design. All animals were genotyped for 2 markers. The three families with significant effects were genotyped for up to 10 additional markers spanning positions 0-122 cM of BTA6. Two sires were segregating for a locus affecting protein and fat percentage near position 55 cM with an estimated substitution effect of 0.18% protein, which is equivalent to one phenotypic standard deviation. This locus was localized to a confidence interval of 4 cM. One of these sires was also heterozygous for a locus affecting milk, fat, and protein production near the centromere. The hypothesis of two segregating loci was verified by multiple regression analysis. A third sire was heterozygous for a locus affecting milk and protein percentage near the telomeric end of the chromosome. Possible candidates for the major quantitative gene near position 55 cM were determined by comparative mapping. IBSP and SSP1 were used as anchors for the orthologous region on human chromosome 4. Twelve genes were detected within a 2-Mbp sequence. None of these genes have been previously associated with lactogenesis.

Partial cDNAs of different isoforms of protein phosphatase 2Cbeta (PP2Cbeta or PPM1B) have been characterized in mammals. We disclose here the full cDNAs of two major PP2Cbeta isoforms from human, rat and mouse. These cDNAs (2.6 and 3.3 kb) are able to encode 53 kDa (PP2Cbetal) and 43 kDa (PP2Cbetas) polypeptides, respectively. The isoforms are co-expressed ubiquitously with the highest level in skeletal muscle, as assessed by Northern-blot analysis. Western and in situ analyses using monoclonal antibodies against PP2Cbeta confirmed the existence of two isoforms in the cytoplasm. Comparative sequence analysis revealed that both cDNAs consist of six exons with an alternate usage of the 3' exons that underlies the differences between them. The genomic structure of PP2Cbeta is similar to that of other PP2C paralogs and includes a non-coding first exon followed by a large intron and a large second exon that encoded most of the catalytic domain. Both variants of the ending exon include large non-coding regions. All non-translated regions (NTRs) are highly conserved between the orthologous genes, indicating their regulatory function. The 5'-NTR is long (379 bp), includes upstream start codons and is predicted to contain stable secondary structures. Such features inhibit translation initiation by the scanning mechanism. Introduction of this NTR element into a bi-luciferase expression-cassette enabled expression of the second cistron, suggesting that it might serve as an internal ribosome entry site, or it contains a cryptic promoter. Overexpression of PP2Cbeta under CMV-promoter in 293 cells led to cell-growth arrest or cell death., (Copyright 2001 Academic Press.)

We have recently cloned the mouse activity-dependent neuroprotective protein (ADNP). Here, we disclose the cloning of human ADNP (hADNP) from a fetal brain cDNA library. Comparative sequence analysis of these two ADNP orthologs indicated 90% identity at the mRNA level. Several single nucleotide polymorphic sites were noticed. The deduced protein structure contained nine zinc fingers, a proline-rich region, a nuclear bipartite localization signal, and a homeobox domain profile, suggesting a transcription factor function. Further comparative analysis identified an ADNP paralog (33% identity and 46% similarity), indicating that these genes belong to a novel protein family with a nine-zinc finger motif followed by a homeobox domain. The hADNP gene structure spans approximately 40 kilobases and includes five exons and four introns with alternative splicing of an untranslated second exon. The hADNP gene was mapped to chromosome 20q12-13.2, a region associated with aggressive tumor growth, frequently amplified in many neoplasias, including breast, bladder, ovarian, pancreatic, and colon cancers. hADNP mRNA is abundantly expressed in distinct normal tissues, and high expression levels were encountered in malignant cells. Down-regulation of ADNP by antisense oligodeoxynucleotides up-regulated the tumor suppressor p53 and reduced the viability of intestinal cancer cells by 90%. Thus, ADNP is implicated in maintaining cell survival, perhaps through modulation of p53.

We have developed an elimination test to identify chromosomal regions that contain tumor inhibitory genes. Monochromosomal human/mouse microcell hybrids are generated and passaged through SCID mice. Derived tumors are then analyzed for deletions on the transgenomic chromosome. Using this strategy, we have previously identified a 1.6-cM common eliminated region 1 (CER1) on human 3p21. 3. We now report that CER1 contains 14 markers that are deleted in 19 SCID-derived tumors. A 1-Mb PAC contig that spans CER1 was assembled. Five chemokine receptor genes (CCR1, CCR3, CCR2, CCR5, and CCR6) were localized in CER1 in a 225-kb cluster. The lactotransferrin gene (LTF, or lactoferrin, LF), which reportedly has tumor inhibitory activity, also maps to CER1. Our results create a basis for characterization and further functional testing of genes within CER1., (Copyright 1999 Academic Press.)

The avian tom1 (target of myb 1) gene has been previously characterized from v-myb-transformed cells. We report here cloning of the human and mouse tom1 orthologs. Both genes are expressed ubiquitously, with the highest levels in skeletal muscle, brain, and intestines, as assessed by Northern blot and mRNA in situ hybridization. The N-terminal domain of the TOM1 protein shares similarity with HGS (hepatocyte growth factor-regulated tyrosine kinase substrate) and STAM (signal-transducing adaptor molecule), which are associated with vesicular trafficking at the endosome. A putative coiled-coil domain was also detected in the central part of the TOM1 protein. This domain structure suggests that TOM1 is another member of a family of genes implicated in the trafficking regulation of growth-factor-receptor complexes that are destined for degradation in the lysosome. We also show that a human paralog of TOM1 (TOM1-like gene 1) exists. Furthermore, we provide a transcription map over a 190-kb contig of the TOM1 region. This map includes its distal neighbors HMOX1 and MCM5 and two proximal novel genes, one of which is a HMG-box-containing gene (HMG2L1), and the other of unknown function. Using a genomic PAC clone, we demonstrate that the mouse Tom1 and Hmox1 genes are part of an as yet undescribed syntenic group between mouse chromosome 8C1 and human chromosome 22q13.1., (Copyright 1999 Academic Press.)

Meningioma, a tumor of the meninges covering the central nervous system, shows frequent loss of material from human chromosome 22. Homozygous and heterozygous deletions in meningiomas defined a candidate region of >1 Mbp in 22q12.3-q13.1 and directed us to gene cloning in this segment. We characterized a new member of the N-acetylglucosaminyltransferase gene family, the LARGE gene. It occupies >664 kilobases and is one of the largest human genes. The predicted 756-aa N-acetylglucosaminyltransferase encoded by LARGE displays features that are absent in other glycosyltransferases. The human like-acetylglucosaminyltransferase polypeptide is much longer and contains putative coiled-coil domains. We characterized the mouse LARGE ortholog, which encodes a protein 97.75% identical with the human counterpart. Both genes reveal ubiquitous expression as assessed by Northern blot analysis and in situ histochemistry. Chromosomal mapping of the mouse gene reveals that mouse chromosome 8C1 corresponds to human 22q12.3-q13.1. Abnormal glycosylation of proteins and glycosphingolipids has been shown as a mechanism behind an increased potential of tumor formation and/or progression. Human tumors overexpress ganglioside GD3 (NeuAcalpha2,8NeuAcalpha2, 3Galbeta1,4Glc-Cer), which in meningiomas correlates with deletions on chromosome 22. It is the first time that a glycosyltransferase gene is involved in tumor-specific genomic rearrangements. An abnormal function of the human like-acetylglucosaminyltransferase protein may be linked to the development/progression of meningioma by altering the composition of gangliosides and/or by effect(s) on other glycosylated molecules in tumor cells.

Dermatofibrosarcoma protuberans (DFSP) and its juvenile form, giant-cell fibroblastoma (GCF), are uncommon infiltrative tumors of the dermis, which present unique cytogenetic features, such as the reciprocal translocation t(17;22) or, more commonly, supernumerary ring chromosomes containing sequences from chromosomes 17 and 22. We have recently shown that these aberrations are cytogenetic manifestations of gene fusions between the platelet-derived growth factor B-chain gene (PDGFB), the cellular equivalent of the v-sis oncogene, and the collagen type 1 alpha 1 gene (COL1A1), the major protein constituent of the extracellular matrix in connective tissue of skin. We now report characterization of COL1A1/PDGFB chimeric genes at the RNA and DNA sequence levels in a series of DFSPs and GCFs. All 16 tumors studied contained the COL1A1/PDGFB gene. The location of breakpoints within COL1A1 varied greatly, but was always limited to the region encoding the alpha-helical domain. The PDGFB segment of the chimeric transcript always starts with exon 2, placing PDGFB under the control of the COL1A1 promoter and removing all known elements negatively controlling PDGFB transcription and translation. Production of these aberrant transcripts in fibroblasts, the suspected cell of origin of DFSP/GCF, likely causes autocrine stimulation and cell proliferation. No specific function has yet been assigned to exon 2 of PDGFB, and this exon does not encode for the mature growth factor. Its retention in all chimeric COL1A1/PDGFB genes suggests that it is important for the normal processing of the PDGFB polypeptide.

Genomic sequencing was combined with searches of databases for identification of active genes on human chromosome 22. A cosmid from 22q13, located in the telomeric vicinity of the PDGFB (platelet-derived growth factor B-chain) gene, was fully sequenced. Using an expressed sequence tag-based approach we characterized human (SYNGR1) and mouse (Syngr1) orthologs of the previously cloned rat synaptogyrin gene (RATSYNGR1). The human SYNGR1 gene reveals three (SYNGR1a, SYNGR1b, SYNGR1c) alternative transcript forms of 4.5, 1.3 and 0.9 kb, respectively. The transcription of SYNGR1 starts from two different promoters, and leads to predicted proteins with different N- and C-terminal ends. The most abundant SYNGR1 a transcript, the 4.5-kb form, which corresponds to RATSYNGR1, is highly expressed in neurons of the central nervous system and at much lower levels in other tissues, as determined by in situ hybridization histochemistry. The levels of SYNGR1b and SYNGR1c transcripts are low and limited to heart, skeletal muscle, ovary and fetal liver. We also characterized two additional members of this novel synaptogyrin gene family in human (SYNGR2 and SYNGR3), and one in mouse (Syngr2). The human SYNGR2 gene transcript of 1.6 kb is expressed at high levels in all tissues, except brain. The 2.2-kb SYNGR3 transcript was detected in brain and placenta only. The human SYNGR2 and SYNGR3 genes were mapped by fluorescence in situ hybridization to 17qtel and 16ptel, respectively. The human SYNGR2 gene has a processed pseudogene localized in 15q11. All predicted synaptogyrin proteins contain four strongly conserved transmembrane domains, which is consistent with the M-shaped topology. The C-terminal polypeptide ends are variable in length, display a low degree of sequence similarity between family members, and are therefore likely to convey the functional specificity of each protein.

Multiple endocrine neoplasia type 1 (MEN1) is tightly linked to the muscle-type glycogen phosphorylase (PYGM) gene in 11q13. This region of the human genome contains additional disease-related loci implicated in the development of insulin-dependent diabetes mellitus, familial paraganglioma type 2, spinocerebellar ataxia type 5, Bardet-Biedl syndrome and translocation t(11;17) described in B-cell non-Hodgkin's lymphoma. We approached cloning of candidate disease genes from 11q13 by large-scale genomic sequencing. We obtained > 106 kb of sequence around the PYGM gene and established a transcriptional map that includes: (i) two genes previously localized to 11q13, PYGM and a zinc-finger protein (ZFM1) gene; (ii) the germinal center kinase (GCK, human B-lymphocyte serine/threonine protein kinase) gene; (iii) a novel human CDC25-like (HCDC25L) gene; (iv) a dystrophia myotonica protein kinase-like (DMPKL) gene; and (v) a novel ubiquitously expressed gene of unknown function (germinal center kinase- neighboring gene, GCKNG).

Usher syndrome is a group of autosomal recessive disorders associated with congenital sensorineural deafness and progressive visual loss due to retinitis pigmentosa. Sixteen members of the small inbred Samaritan isolate with autosomal recessive deafness were studied in 10 related sibships. DNA samples from 59 individuals including parents and affected and nonaffected sibs were typed for markers on chromosomes 1q and 11q for which linkage has recently been established for Usher syndrome types II and I. Statistically significant linkage was observed with four markers on 11q (D11S533, D11S527, OMP, and INT2) with a maximum six-point location score of 11.61 at the D11S533 locus. Analysis of haplotypes supports the notion that the mutation arose only once in an ancestral chromosome carrying a specific haplotype. The availability of markers closely linked to the disease locus allows indirect genotype analysis and identifies all carriers of the gene within the community. Furthermore, the detection of complete linkage disequilibrium between the D11S533 marker and the Usher gene suggests that these loci are either identical or adjacent and narrows the critical region to which physical mapping efforts are currently directed.

A novel method for efficient and rapid isolation of dsRNA molecules was developed. The dsRNA content of Ustilago maydis was reexamined; two distinct dsRNA classes were identified. Class I includes the dsRNA segments reported earlier for U. maydis virus systems and class II includes unencapsidated dsRNA molecules that were barely detected by the conventional extraction methods despite their high titer. Segments of the class II, some of which are reported for the first time, were further characterized; all the segments are independent of the killer system and other encapsidated dsRNA molecules. These segments are cytoplasmically transmitted and, in sharp contrast with class I-encapsidated dsRNA segments, their relative copy number decreases rapidly while entering the stationary phase.