Evolutionary sequence comparisons using high-density oligonucleotide arrays.

We explored the utility of high-density oligonucleotide arrays (DNA chips) for obtaining sequence information from homologous genes in closely related species. Orthologues of the human BRCA1 exon 11, all approximately 3.4 kb in length and ranging from 98.2% to 83.5% nucleotide identity, were subjected to hybridization-based and conventional dideoxysequencing analysis. Retrospective guidelines for identifying high-fidelity hybridization-based sequence calls were formulated based upon dideoxysequencing results. Prospective application of these rules yielded base-calling with at least 98.8% accuracy over orthologous sequence tracts shown to have approximately 99% identity. For higher primate sequences with greater than 97% nucleotide identity, base-calling was made with at least 99.91% accuracy covering a minimum of 97% of the sequence. Using a second-tier confirmatory hybridization chip strategy, shown in several cases to confirm the identity of predicted sequence changes, the complete sequence of the chimpanzee, gorilla and orangutan orthologues should be deducible solely through hybridization-based methodologies. Analysis of less highly conserved orthologues can still identify conserved nucleotide tracts of at least 15 nucleotides and can provide useful information for designing primers. DNA-chip based assays can be a valuable new technology for obtaining high-throughput cost-effective sequence information from related genomes.

Determination of ancestral alleles for human single-nucleotide polymorphisms using high-density oligonucleotide arrays.

Here we report the application of high-density oligonucleotide array (DNA chip)-based analysis to determine the distant history of single nucleotide polymorphisms (SNPs) in current human populations. We analysed orthologues for 397 human SNP sites (identified in CEPH pedigrees from Amish, Venezuelan and Utah populations) from 23 common chimpanzee, 19 pygmy chimpanzee and 11 gorilla genomic DNA samples. From this data we determined 214 proposed ancestral alleles (the sequence found in the last common ancestor of humans and chimpanzees). In a diverse human population set, we found that SNP alleles with higher frequencies were more likely to be ancestral than less frequently occurring alleles. There were, however, exceptions. We also found three shared human/pygmy chimpanzee polymorphisms, all involving CpG dinucleotides, and two shared human/gorilla polymorphisms, one involving a CpG dinucleotide. We demonstrate that microarray-based assays allow rapid comparative sequence analysis of intra- and interspecies genetic variation.

HRPT2, encoding parafibromin, is mutated in hyperparathyroidism-jaw tumor syndrome.

We report here the identification of a gene associated with the hyperparathyroidism-jaw tumor (HPT-JT) syndrome. A single locus associated with HPT-JT (HRPT2) was previously mapped to chromosomal region 1q25-q32. We refined this region to a critical interval of 12 cM by genotyping in 26 affected kindreds. Using a positional candidate approach, we identified thirteen different heterozygous, germline, inactivating mutations in a single gene in fourteen families with HPT-JT. The proposed role of HRPT2 as a tumor suppressor was supported by mutation screening in 48 parathyroid adenomas with cystic features, which identified three somatic inactivating mutations, all located in exon 1. None of these mutations were detected in normal controls, and all were predicted to cause deficient or impaired protein function. HRPT2 is a ubiquitously expressed, evolutionarily conserved gene encoding a predicted protein of 531 amino acids, for which we propose the name parafibromin. Our findings suggest that HRPT2 is a tumor-suppressor gene, the inactivation of which is directly involved in predisposition to HPT-JT and in development of some sporadic parathyroid tumors.

A common nonsense mutation in EphB2 is associated with prostate cancer risk in African American men with a positive family history.

BACKGROUND: The EphB2 gene was recently implicated as a prostate cancer (PC) tumour suppressor gene, with somatic inactivating mutations occurring in approximately 10% of sporadic tumours. We evaluated the contribution of EphB2 to inherited PC susceptibility in African Americans (AA) by screening the gene for germline polymorphisms. METHODS: Direct sequencing of the coding region of EphB2 was performed on 72 probands from the African American Hereditary Prostate Cancer Study (AAHPC). A case-control association analysis was then carried out using the AAHPC probands and an additional 183 cases of sporadic PC compared with 329 healthy AA male controls. In addition, we performed an ancestry adjusted association study where we adjusted for individual ancestry among all subjects, in order to rule out a spurious association due to population stratification. RESULTS: Ten coding sequence variants were identified, including the K1019X (3055A-->T) nonsense mutation which was present in 15.3% of the AAHPC probands but only 1.7% of 231 European American (EA) control samples. We observed that the 3055A-->T mutation significantly increased risk for prostate cancer over twofold (Fisher's two sided test, p = 0.003). The T allele was significantly more common among AAHPC probands (15.3%) than among healthy AA male controls (5.2%) (odds ratio 3.31; 95% confidence interval 1.5 to 7.4; p = 0.008). The ancestry adjusted analyses confirmed the association. CONCLUSIONS: Our data show that the K1019X mutation in the EphB2 gene differs in frequency between AA and EA, is associated with increased risk for PC in AA men with a positive family history, and may be an important genetic risk factor for prostate cancer in AA.

The human RGL (RalGDS-like) gene: cloning, expression analysis and genomic organization.

Ral GDP dissociation stimulator (RalGDS) and its family members RGL, RLF and RGL2 are involved in Ras and Ral signaling pathways as downstream effector proteins. Here we report the precise localization and cloning of two forms of human RGL gene differing at the amino terminus. Transcript A, cloned from liver cDNA libraries has the same amino terminus as the mouse RGL, whereas transcript B cloned from brain has a substitution of 45 amino acids for the first nine amino acids. At the genomic level, exon 1 of transcript A is replaced by two alternative exons (1B1 and 1B2) in transcript B. Both forms share exons 2 through 18. The human RGL protein shares 94% amino acid identity with the mouse protein. Northern blot analysis shows that human RGL is expressed in a wide variety of tissues with strong expression being seen in the heart, brain, kidney, spleen and testis.

Cloning, genomic organization and expression of a putative human transmembrane protein related to the Caenorhabditis elegans M01F1.4 gene.

A novel human transcript CG-2 (C9ORF5), was isolated from the familial dysautonomia candidate region on 9q31 using a combination of cDNA selection and exon trapping. CG-2 was detected as a relatively abundant 8kb transcript in all adult and fetal tissues with the exception of adult thymus. Genomic analysis of CG-2 identified 18 exons that span more than 110kb. The gene encodes a 911-amino-acid protein with a predicted molecular weight of 101kDa and a hypothetical pI of 9.03. Sequence analysis of CG-2 indicates that it is likely to encode a transmembrane protein. Here, we assess CG-2 as a candidate for familial dysautonomia.

Cloning and characterization of an inversion breakpoint at 6q23.3 suggests a role for Map7 in sacral dysgenesis.

Here we report on a male patient with sacral dysgenesis (SD) and constitutional pericentric inversion of chromosome 6 (p11.2;q23.3). SD is a heterogeneous group of congenital anomalies with complex genetic etiology. Previously, a patient with sacral abnormalities and an interstitial deletion of 6q23-->q25 region has been described. We speculated that a susceptibility gene for SD lies in 6q23.3 region (disrupted in both patients), and therefore, cloning of the breakpoint in our patient would lead to the identification of the disrupted gene. We performed FISH analysis followed by Southern blot analysis and inverse PCR to clone the breakpoint. The 6p11.2 breakpoint mapped very close to the centromere, and the 6q23.3 breakpoint localized in the ninth intron of the MAP7 gene. We then evaluated the involvement of MAP7 in SD by further screening of the gene in several patients with a similar phenotype. Two nucleotide changes causing Ile257Asn and Glu571Ala substitutions in the protein, both affecting amino acid residues conserved in the mouse homolog, were identified in two patients. Both changes are either very rare polymorphisms or true mutations, since they were not detected in 167 normal individuals nor found in the SNP database. Therefore, our study suggests MAP7 as a candidate gene for SD. However, we were unable to detect any sacral defects in the MAP7 knockout mice.

Sequencing homopolymer tracts and repetitive elements.

We investigated the use of Taq dye primer and Taq terminator sequencing chemistry to optimize the quality of sequence data obtained from templates containing homopolymer tracts and repetitive elements. In direct side-by-side comparisons using the Applied Biosystems Model 373A Fluorescent Sequencer, the Taq terminator sequencing chemistry gave much cleaner and more consistent results on long homopolymer tracts and dinucleotide repeats. We also investigated various thermal cycling conditions and determined that higher annealing temperatures and longer denaturation times improved the ability to sequence through these problem templates.

Influence of galactosyl lactose on energy and protein digestibility, enteric microflora, and performance of weanling pigs.

In three experiments, weanling pigs were used to investigate the effect of feeding .5% galactosyl lactose in a simple corn-soy diet on nutrient digestibility, selected ileal microbial concentrations, fermentation acids, and performance. In Exp. 1, 20 21-d-old weaned pigs in two replicate trials were fed either a 19% CP corn and soybean meal diet or a similar diet containing .5% galactosyl lactose. Feces and urine were collected over a 3-wk period to determine apparent CP and energy digestibilities. Apparent digestibilities of CP and energy were not different between treatments. In Exp. 2, 24 14-d-old nursing pigs were cannulated at the terminal ileum and fed the diets described above after weaning at 21 d of age. Ileal samples were collected at 20, 23, 27, 30, 34, 37, and 41 d of age. Samples were assayed for VFA, lactate, and selected microflora. Ileal fermentation acid concentrations and selected microflora were affected by age of pig but were not different between feed treatments. In Exp. 3, 91 21- to 24-d-old weaned pigs were fed the control and test diets already described; intake, BW gain, and feed efficiency were determined for 9 wk after weaning. Intake was greater in pigs fed galactosyl lactose, but ADG and gain/feed did not differ between treatments. Inclusion of galactosyl lactose in weanling pig diets did not affect protein or energy digestibility, microflora, short-chain fatty acids, or performance of young pigs.

Effects of a direct-fed yeast culture on enteric microbial populations, fermentation acids, and performance of weanling pigs.

In three replicate trials, a total of 36 pigs that had been cannulated at the terminal ileum were used to determine the effects of a Saccharomyces cerevisiae culture in a phase feeding program (phase I was d 0 to 7 and phase II was d 8 to 21) on performance, ileal microflora, and short-chain fatty acids in weanling pigs. Pigs were cannulated at approximately 12 d of age, weaned at 17 d of age, and randomly assigned to one of three treatments: 1) a pelleted phase feeding program, 2) a similar program with the inclusion of a live S. cerevisiae culture (1 g/ kg), and 3) a nonpelleted feeding program otherwise similar to program 2. Ileal samples were collected at 17, 20, 24, 27, 31, 34, and 38 d of age, and samples were analyzed for total E. coli, streptococci, lactobacilli, yeast, short-chain fatty acids, pH, and dry matter. Performance data were also collected. At 41 d of age, pigs were killed and digesta were collected from various regions of the gastrointestinal tract. Total intake was less for pigs fed the control diet than for pigs fed the yeast diets, and overall gains tended to be greater for pigs fed diets including yeast. Treatment differences were not observed for ileal microflora or short-chain fatty acids in samples obtained from cannulas or from the various sites of the gastrointestinal tract. Inclusion of a live yeast culture in weanling pig diets affected intake and performance but did not alter tested intestinal microflora or net concentrations of fermentation products.

Admixture and population stratification in African Caribbean populations.

Throughout biomedical research, there is growing interest in the use of ancestry informative markers (AIMs) to deconstruct racial categories into useful variables. Studies on recently admixed populations have shown significant population substructure due to differences in individual ancestry; however, few studies have examined Caribbean populations. Here we used a panel of 28 AIMs to examine the genetic ancestry of 298 individuals of African descent from the Caribbean islands of Jamaica, St. Thomas and Barbados. Differences in global admixture were observed, with Barbados having the highest level of West African ancestry (89.6%+/- 2.0) and the lowest levels of European (10.2%+/- 2.2) and Native American ancestry (0.2%+/- 2.0), while Jamaica possessed the highest levels of European (12.4%+/- 3.5) and Native American ancestry (3.2%+/- 3.1). St. Thomas, USVI had ancestry levels quite similar to African Americans in continental U.S. (86.8%+/- 2.2 West African, 10.6%+/- 2.3 European, and 2.6%+/- 2.1 Native American). Significant substructure was observed in the islands of Jamaica and St. Thomas but not Barbados (K=1), indicating that differences in population substructure exist across these three Caribbean islands. These differences likely stem from diverse colonial and historical experiences, and subsequent evolutionary processes. Most importantly, these differences may have significant ramifications for case-control studies of complex disease in Caribbean populations.

Uterine tumours are a phenotypic manifestation of the hyperparathyroidism-jaw tumour syndrome.

The hyperparathyroidism-jaw tumour (HPT-JT) syndrome is an autosomal dominant disorder characterized by parathyroid tumours, which are frequently carcinomas, and ossifying jaw fibromas. In addition, some patients may develop renal tumours and cysts. The gene causing HPT-JT, which is referred to as HRPT2 and is located on chromosome 1q31.2, encodes a 531 amino acid protein called PARAFIBROMIN. To date 42 mutations, of which 22 are germline, have been reported and 97% of these are inactivating and consistent with a tumour suppressor role for HRPT2. We have investigated another four HPT-JT families for germline mutations, searched for additional clinical phenotypes, and examined for a genotype-phenotype correlation. Mutations were found in two families. One family had a novel deletional-insertion at codon 669, and the other had a 2 bp insertion at codon 679, which has been reported in four other unrelated patients. These five unrelated patients and their families with the same mutation were not found to develop the same tumours, thereby indicating an absence of a genotype-phenotype correlation. An analysis of 33 HPT-JT kindreds revealed that affected women in 13 HPT-JT families suffered from menorrhagia in their second to fourth decades. This often required hysterectomy, which revealed the presence of uterine tumours. This resulted in a significantly reduced maternal transmission of the disease. Thus, the results of our analysis expand the spectrum of HPT-JT-associated tumours to include uterine tumours, and these may account for the decreased reproductive fitness in females from HPT-JT families.

An IL-2 response element in the human IL-2 receptor alpha chain promoter is a composite element that binds Stat5, Elf-1, HMG-I(Y) and a GATA family protein.

Expression of the human interleukin-2 (IL-2) receptor alpha chain gene is potently upregulated by its own ligand, IL-2. In this study, we characterize an essential upstream IL-2 response element that contains both consensus and non-consensus GAS motifs, two putative Ets binding sites (EBS), one of which overlaps the consensus GAS motif, and a GATA motif, which overlaps the non-consensus GAS motif. We demonstrate that although the individual components of this element do not respond to IL-2, together they form a composite element capable of conferring IL-2 responsiveness to a heterologous promoter. Multiple factors including Stat5, Elf-1, HMG-I(Y) and GATA family proteins bind to the IL-2 response element and mutation of any one of these binding sites diminishes the activity of this element. An unidentified Ets family protein binds to the EBS overlapping the consensus GAS motif and appears to negatively regulate the human IL-2R alpha promoter. Thus, IL-2-induced IL-2R alpha promoter activity requires a complex upstream element, which appears to contain binding sites for both positive and negative regulatory factors.

Occurrence of embryotoxicity in mouse embryos following in utero exposure to 2'-deoxycoformycin (pentostatin).

Previous investigations had shown that i.p. injection of 2'-deoxycoformycin (dCF; pentostatin; 5 mg/kg) on either E7 or E8 into pregnant mice results in a 61-81% resorption rate at E17. The incidence of visible gross malformations among the surviving conceptuses was exceptionally low (3%) at the time of necropsy on E17 and was unrelated to dCF dose (Knudsen et al., Teratology, 40:5-626, '89; Teratology, 45:91-103, '92). These findings demonstrated the embryotoxicity of dCF but provided no clues as to the site(s) of dCF action. To define the lesion site(s), we have now examined embryos at 72 h (E10), 96 h (E11), and 120 h (E12) following administration of a highly embryotoxic dose of 5 mg dCF/kg to dams on E7. Deoxycoformycin caused multiple abnormalities and growth retardation, and the temporal sequence between maximal abnormal embryo incidence and resorption frequency was established. The quantitative data show that the maximal occurrence of abnormal embryos on E10 (71%) was followed by a maximal resorption rate on E12 (78%). There was a strong correlation (r = -0.82; P < 0.05) between the rapid decline of percent abnormal embryos over E10-E12 and the simultaneous increase in resorption rate, with linear regression analysis showing nearly equal but opposite slopes (-31.2% vs. +35.8% per gestational day, respectively). This suggests that one or more of the abnormalities seen at E10 is associated with the death and resorption of the embryo at E12. The dCF treatment perturbed a wide spectrum of developmental events, including neural tube closure, craniofacial and limb development, turning of the embryo, and growth retardation. None of the individual abnormalities, however, can quantitatively account for the high percentage of dead and resorbed embryos. Therefore, the specific cause of dCF-induced embryolethality is not clear. There is evidence both for direct dCF toxicity at specific embryonic sites as well as for a generalized retardation in the rate of development.

Developing allantois is a primary site of 2'-deoxycoformycin toxicity.

In this study, an assessment of normal mouse allantoic development and its sensitivity to 2'-(R)-deoxycoformycin (dCF; Pentostatin) exposure were examined. Both dissecting microscopy and scanning electron microscopy were used to describe the normal growth and morphogenesis of the mouse allantois over gestational days 7-10 as a preliminary step in evaluating potential abnormal allantoic ontogeny and its effect on umbilical cord and placental development. Two abnormal allantoic/umbilical cord phenotypes were observed subsequent to injecting pregnant mice with 5 mg dCF/kg, i.p., on gestational day 7 (GD 7) and evaluating litters on GD 10, 11, and 12. Abnormal phenotypes included: (1) an allantois which extended approximately halfway across the exocoelom but failed to establish a functional contact with the chorion; and (2) a phenotype characterized by reduced expansion of the allantois across the chorionic surface, a very thin umbilical cord, and aberrant vascularization throughout the structure. Both abnormal phenotypes exhibited either an agenesis or hypogenesis of the umbilical cord and chorioallantoic plate, respectively. Neither abnormal phenotype, however, exhibited errors in the directionality of allantoic growth toward the chorion nor in the formation of aberrant contacts between allantois and adjacent yolk sac or amnionic mesenchyme. Statistical interpretation of the experimental data strongly suggested that abnormalities in allantoic/umbilical cord development were directly associated with embryolethality as evidenced by a decline in the frequency of abnormal allantoic/umbilical cord phenotypes over GD 10-12 (73, 36, and 4%; respectively) and a concomitant increase in the frequency of implantation site resorptions over the same time period (7, 47, and 78%). These results strongly suggest that the developing allantois is very sensitive to the effects of dCF exposure, and that interference with its development leads to embryolethality by GD 12.

Isolation and characterization of the human homeobox gene HOX D1.

Homeobox genes, first identified in Drosophila, encode transcription factors that regulate embryonic development along the anteroposterior axis of an organism. Vertebrate homeobox genes are described on the basis of their homology to the genes found within the Drosophila Antennapedia and Bithorax homeotic gene complexes. Mammals possess four paralogous homeobox (HOX) gene clusters, HOX A, HOX B, HOX C and HOX D, each located on different chromosomes, consisting of 9 to 11 genes arranged in tandem. We report the characterization of the human HOX D1 gene. This gene consists of two exons, encoding a 328 amino acid protein, separated by an intron of 354 bp. The human HOX D1 protein is one amino acid longer (328 amino acids) than the mouse protein (327 amino acids) and is 82% identical to the mouse HOX D1 homolog. The DNA binding homeodomain region of the human protein exhibits a 97% and 80% identity between mouse Hoxd1 and Drosophila labial homeodomains, respectively. The exon/intron and intron/exon splice junctions are conserved in position between human and mouse genes. Determination of the human HOX D1 gene structure permits the use of PCR based analysis of this gene for the assessment of mutations, for diseases that link to the HOXD cluster (such as Duanes Retraction Syndrome (DRS)), or polymorphisms associated with human variation. Molecular characterization of the HOXD1 gene may also permit analysis of the functional role of this gene in human neurogenisis.

Genomic organization of the murine Miller-Dieker/lissencephaly region: conservation of linkage with the human region.

Several human syndromes are associated with haploinsufficiency of chromosomal regions secondary to microdeletions. Isolated lissencephaly sequence (ILS), a human developmental disease characterized by a smooth cerebral surface (classical lissencephaly) and microscopic evidence of incomplete neuronal migration, is often associated with small deletions or translocations at chromosome 17p13.3. Miller-Dieker syndrome (MDS) is associated with larger deletions of 17p13.3 and consists of classical lissencephaly with additional phenotypes including facial abnormalities. We have isolated the murine homologs of three genes located inside and outside the MDS region: Lis1, Mnt/Rox, and 14-3-3 epsilon. These genes are all located on mouse chromosome 11B2, as determined by metaphase FISH, and the relative order and approximate gene distance was determined by interphase FISH analysis. The transcriptional orientation and intergenic distance of Lis1 and Mnt/Rox were ascertained by fragmentation analysis of a mouse yeast artificial chromosome containing both genes. To determine the distance and orientation of 14-3-3 epsilon with respect to Lis1 and Mnt/Rox, we introduced a super-rare cutter site (VDE) that is unique in the mouse genome into 14-3-3 epsilon by gene targeting. Using the introduced VDE site, the orientation of this gene was determined by pulsed field gel electrophoresis and Southern blot analysis. Our results demonstrate that the MDS region is conserved between human and mouse. This conservation of linkage suggests that the mouse can be used to model microdeletions that occur in ILS and MDS.

2006 expressed-sequence tags derived from human chromosome 7-enriched cDNA libraries.

The establishment and mapping of gene-specific DNA sequences greatly complement the ongoing efforts to map and sequence all human chromosomes. To facilitate our studies of human chromosome 7, we have generated and analyzed 2006 expressed-sequence tags (ESTs) derived from a collection of direct selection cDNA libraries that are highly enriched for human chromosome 7 gene sequences. Similarity searches indicate that approximately two-thirds of the ESTs are not represented by sequences in the public databases, including those in dbEST. In addition, a large fraction (68%) of the ESTs do not have redundant or overlapping sequences within our collection. Human DNA-specific sequence-tagged sites (STSs) have been developed from 190 of the ESTs. Remarkably, 180 (96%) of these STSs map to chromosome 7, demonstrating the robustness of chromosome enrichment in constructing the direct selection cDNA libraries. Thus far, 140 of these EST-specific STSs have been assigned unequivocally to YAC contigs that are distributed across the chromosome. Together, these studies provide > 2000 ESTs highly enriched for chromosome 7 gene sequences, 180 new chromosome 7 STSs corresponding to ESTs, and a definitive demonstration of the ability to enrich for chromosome-specific cDNAs by direct selection. Furthermore, the libraries, sequence data, and mapping information will contribute to the construction of a chromosome 7 transcript map.

Positional cloning utilizing genomic DNA microarrays: the Niemann-Pick type C gene as a model system.

A major obstacle in positional cloning is identifying the specific mutated gene from within a large physical contig. Here we describe the application of DNA microarray technology to a defined genomic region (physical map) to identify: (i) exons without a priori sequence data and (ii) the disease gene based on differential gene expression in a recessive disorder. The feasibility was tested using resources from the positional cloning of the Neimann-Pick Type C (NP-C) disease gene, NPC1. To identify NPC1 exons and optimize the technology, an array was generated from genomic fragments of the 110-kb bacterial artificial chromosome, 108N2, which encodes NPC1. First, as a test case for blindly identifying exons, fluorescently labeled NPC1 cDNA identified 108N2 fragments that contained NPC1 exons, many of which also contained intronic sequences and could be used to determine part of the NPC1 genomic structure. Second, to demonstrate that the NPC1 disease gene could be identified based upon differential gene expression, subarrays of 108N2 fragments were hybridized with fluorescently labeled cDNA probes generated from total RNA from hamster cell lines differentially expressing NPC1. A probe derived from the NP-C cell line CT60 did not detect NPC1 exons or other genomic fragments from 108N2. In contrast, several NPC1 exons were detected by a probe generated from the non-NP-C cell line 911D5A13, which was derived from CT60, and expressed NPC1 as a consequence of stable transduction with a YAC that contains NPC1 and encompasses 108N2. Thus, the array technology identified NPC1 as a candidate gene based on a physical contig and differential NPC1 expression between NP-C and non-NP-C cells. This technique should facilitate gene identification when a physical contig exists for a region of interest and mutations result in changes in the mRNA level of the disease gene or portions thereof.