Comparative analyses of multi-species sequences from targeted genomic regions.

The systematic comparison of genomic sequences from different organisms represents a central focus of contemporary genome analysis. Comparative analyses of vertebrate sequences can identify coding and conserved non-coding regions, including regulatory elements, and provide insight into the forces that have rendered modern-day genomes. As a complement to whole-genome sequencing efforts, we are sequencing and comparing targeted genomic regions in multiple, evolutionarily diverse vertebrates. Here we report the generation and analysis of over 12 megabases (Mb) of sequence from 12 species, all derived from the genomic region orthologous to a segment of about 1.8 Mb on human chromosome 7 containing ten genes, including the gene mutated in cystic fibrosis. These sequences show conservation reflecting both functional constraints and the neutral mutational events that shaped this genomic region. In particular, we identify substantial numbers of conserved non-coding segments beyond those previously identified experimentally, most of which are not detectable by pair-wise sequence comparisons alone. Analysis of transposable element insertions highlights the variation in genome dynamics among these species and confirms the placement of rodents as a sister group to the primates.

The human pregnane X receptor: genomic structure and identification and functional characterization of natural allelic variants.

The pregnane X receptor (PXR)/steroid and xenobiotic receptor (SXR) transcriptionally activates cytochrome P4503A4 (CYP3A4) when ligand activated by endobiotics and xenobiotics. We cloned the human PXR gene and analysed the sequence in DNAs of individuals whose CYP3A phenotype was known. The PXR gene spans 35 kb, contains nine exons, and mapped to chromosome 13q11-13. Thirty-eight single nucleotide polymorphisms (SNPs) were identified including six SNPs in the coding region. Three of the coding SNPs are non-synonymous creating new PXR alleles [PXR*2, P27S (79C to T); PXR*3, G36R (106G to A); and PXR*4, R122Q (4321G to A)]. The frequency of PXR*2 was 0.20 in African Americans and was never found in Caucasians. Hepatic expression of CYP3A4 protein was not significantly different between African Americans homozygous for PXR*1 compared to those with one PXR*2 allele. PXR*4 was a rare variant found in only one Caucasian person. Homology modelling suggested that R122Q, (PXR*4) is a direct DNA contact site variation in the third alpha-helix in the DNA binding domain. Compared with PXR*1, and variants PXR*2 and PXR*3, only the variant PXR*4 protein had significantly decreased affinity for the PXR binding sequence in electromobility shift assays and attenuated ligand activation of the CYP3A4 reporter plasmids in transient transfection assays. However, the person heterozygous for PXR*4 is normal for CYP3A4 metabolism phenotype. The relevance of each of the 38 PXR SNPs identified in DNA of individuals whose CYP3A basal and rifampin-inducible CYP3A4 expression was determined in vivo and/or in vitro was demonstrated by univariate statistical analysis. Because ligand activation of PXR and upregulation of a system of drug detoxification genes are major determinants of drug interactions, it will now be useful to extend this work to determine the association of these common PXR SNPs to human variation in induction of other drug detoxification gene targets.

Mutation of a new gene causes a unique form of Hermansky-Pudlak syndrome in a genetic isolate of central Puerto Rico.

Hermansky-Pudlak syndrome (HPS) is a rare autosomal recessive disorder characterized by oculocutaneous albinism and a storage pool deficiency due to an absence of platelet dense bodies. Lysosomal ceroid lipofuscinosis, pulmonary fibrosis and granulomatous colitis are occasional manifestations of the disease. HPS occurs with a frequency of one in 1800 in north-west Puerto Rico due to a founder effect. Several non-Puerto Rican patients also have mutations in HPS1, which produces a protein of unknown function. Another gene, ADTB3A, causes HPS in the pearl mouse and in two brothers with HPS-2 (refs. 11,12). ADTB3A encodes a coat protein involved in vesicle formation, implicating HPS as a disorder of membrane trafficking. We sought to identify other HPS-causing genes. Using homozygosity mapping on pooled DNA of 6 families from central Puerto Rico, we localized a new HPS susceptibility gene to a 1.6-cM interval on chromosome 3q24. The gene, HPS3, has 17 exons, and a putative 113.7-kD product expected to reveal how new vesicles form in specialized cells. The homozygous, disease-causing mutation is a large deletion and represents the second example of a founder mutation causing HPS on the small island of Puerto Rico. We also present an allele-specific assay for diagnosing individuals heterozygous or homozygous for this mutation.

The NIEHS Xenopus maternal EST project: interim analysis of the first 13,879 ESTs from unfertilized eggs.

The sequencing of expressed sequence tags (ESTs) from Xenopus laevis has lagged behind efforts on many other common experimental organisms and man, partly because of the pseudotetraploid nature of the Xenopus genome. Nonetheless, large collections of Xenopus ESTs would be useful in gene discovery, oligonucleotide-based knockout studies, gene chip analyses of normal and perturbed development, mapping studies in the related diploid frog X. tropicalis, and for other reasons. We have created a normalized library of cDNAs from unfertilized Xenopus eggs. These cells contain all of the information necessary for the first several cell divisions in the early embryo, as well as much of the information needed for embryonic pattern formation and cell fate determination. To date, we have successfully sequenced 13,879 ESTs out of 16,607 attempts (83.6% success rate), with an average sequence read length of 508 bp. Using a fragment assembly program, these ESTs were assembled into 8,985 'contigs' comprised of up to 11 ESTs each. When these contigs were used to search publicly available databases, 46.2% bore no relationship to protein or DNA sequences in the database at the significance level of 1e-6. Examination of a sample of 100 of the assembled contigs revealed that most ( approximately 87%) were comprised of two apparent allelic variants. Expression profiles of 16 of the most prominent contigs showed that 12 exhibited some degree of zygotic expression. These findings have implications for sequence-specific applications for Xenopus ESTs, particularly the use of allele-specific oligonucleotides for knockout studies, differential hybridization techniques such as gene chip analysis, and the establishment of accurate nomenclature and databases for this species.

Human and mouse alpha-synuclein genes: comparative genomic sequence analysis and identification of a novel gene regulatory element.

The human alpha-synuclein gene (SNCA) encodes a presynaptic nerve terminal protein that was originally identified as a precursor of the non-beta-amyloid component of Alzheimer's disease plaques. More recently, mutations in SNCA have been identified in some cases of familial Parkinson's disease, presenting numerous new areas of investigation for this important disease. Molecular studies would benefit from detailed information about the long-range sequence context of SNCA. To that end, we have established the complete genomic sequence of the chromosomal regions containing the human and mouse alpha-synuclein genes, with the objective of using the resulting sequence information to identify conserved regions of biological importance through comparative sequence analysis. These efforts have yielded approximately 146 and approximately 119 kb of high-accuracy human and mouse genomic sequence, respectively, revealing the precise genetic architecture of the alpha-synuclein gene in both species. A simple repeat element upstream of SNCA/Snca has been identified and shown to be necessary for normal expression in transient transfection assays using a luciferase reporter construct. Together, these studies provide valuable data that should facilitate more detailed analysis of this medically important gene.

The genomic region encompassing the nephropathic cystinosis gene (CTNS): complete sequencing of a 200-kb segment and discovery of a novel gene within the common cystinosis-causing deletion.

Nephropathic cystinosis is an autosomal recessive disorder caused by the defective transport of cystine out of lysosomes. Recently, the causative gene (CTNS) was identified and presumed to encode an integral membrane protein called cystinosin. Many of the disease-associated mutations in CTNS are deletions, including one >55 kb in size that represents the most common cystinosis allele encountered to date. In an effort to determine the precise genomic organization of CTNS and to gain sequence-based insight about the DNA within and flanking cystinosis-associated deletions, we mapped and sequenced the region of human chromosome 17p13 encompassing CTNS. Specifically, a bacterial artificial chromosome (BAC)-based physical map spanning CTNS was constructed by sequence-tagged site (STS)-content mapping. The resulting BAC contig provided the relative order of 43 STSs. Two overlapping BACs, which together contain all of the CTNS exons as well as extensive amounts of flanking DNA, were selected and subjected to shotgun sequencing. A total of 200,237 bp of contiguous, high-accuracy sequence was generated. Analysis of the resulting data revealed a number of interesting features about this genomic region, including the long-range organization of CTNS, insight about the breakpoints and intervening DNA associated with the common cystinosis-causing deletion, and structural information about five genes neighboring CTNS (human ortholog of rat vanilloid receptor subtype 1 gene, CARKL, TIP-1, P2X5, and HUMINAE). In particular, sequence analysis detected the presence of a novel gene (CARKL) residing within the most common cystinosis-causing deletion. This gene encodes a previously unknown protein that is predicted to function as a carbohydrate kinase. Interestingly, both CTNS and CARKL are absent in nearly half of all cystinosis patients (i.e., those homozygous for the common deletion). [The sequence data described in this paper have been submitted to the GenBank data library under accession nos. AF168787 and AF163573.]

Analysis of expressed immunoglobulin heavy chain genes in familial B-CLL.

In this study, we wished to determine whether familial chronic lymphocytic leukemia of B-cell phenotype (CLL) shares with sporadic B-CLL the same immunoglobulin (Ig) heavy chain variable region (VH) gene usage and occurrence of somatic mutation, to gain insight into the pathogenetic relatedness of these epidemiologically distinct forms of CLL. We therefore analyzed the expressed Ig heavy chain genes in 23 cases (11 families) of familial CLL, and compared these results with data previously reported for sporadic CLL. In addition, we assessed the relationship of the occurrence of somatic mutation to several clinical and phenotypic features. The distribution of V genes among these cases was similar to that observed in sporadic CLL: VH3 > VH1 > VH4. Thirteen of the 23 cases (57%) showed germ line VH gene sequences, whereas somatic mutations were detected in 10 cases (43%). The average mutation frequency of these latter 10 cases of was 6.7% (ranging from 1.7% to 8.8%), and evidence of antigen selection was noted in 6. Intraclonal variation, followed by clonal evolution and the appearance of a second clone over a 20-year period was observed in 1 case, suggesting that mutations can continue to accumulate after neoplastic transformation. The presence of somatic mutations correlated with age at presentation, low white blood cell (WBC) count, and low fluorescence intensity of surface CD5, and the potential significance of these relationships is discussed. Our data indicate that familial and sporadic B-CLL display a similar pattern of immunoglobulin gene usage and frequency of somatic mutation, and are consistent with a common ontogeny and immunogenetic origin for these 2 epidemiologically distinct forms of CLL. (Blood. 2000;95:1413-1419)

Comparative genomic sequence analysis of the human and mouse cystic fibrosis transmembrane conductance regulator genes.

The identification of the cystic fibrosis transmembrane conductance regulator gene (CFTR) in 1989 represents a landmark accomplishment in human genetics. Since that time, there have been numerous advances in elucidating the function of the encoded protein and the physiological basis of cystic fibrosis. However, numerous areas of cystic fibrosis biology require additional investigation, some of which would be facilitated by information about the long-range sequence context of the CFTR gene. For example, the latter might provide clues about the sequence elements responsible for the temporal and spatial regulation of CFTR expression. We thus sought to establish the sequence of the chromosomal segments encompassing the human CFTR and mouse Cftr genes, with the hope of identifying conserved regions of biologic interest by sequence comparison. Bacterial clone-based physical maps of the relevant human and mouse genomic regions were constructed, and minimally overlapping sets of clones were selected and sequenced, eventually yielding approximately 1.6 Mb and approximately 358 kb of contiguous human and mouse sequence, respectively. These efforts have produced the complete sequence of the approximately 189-kb and approximately 152-kb segments containing the human CFTR and mouse Cftr genes, respectively, as well as significant amounts of flanking DNA. Analyses of the resulting data provide insights about the organization of the CFTR/Cftr genes and potential sequence elements regulating their expression. Furthermore, the generated sequence reveals the precise architecture of genes residing near CFTR/Cftr, including one known gene (WNT2/Wnt2) and two previously unknown genes that immediately flank CFTR/Cftr.

Mutations in the gene encoding KRIT1, a Krev-1/rap1a binding protein, cause cerebral cavernous malformations (CCM1).

Cerebral cavernous malformations (CCM) are congenital vascular anomalies of the brain that can cause significant neurological disabilities, including intractable seizures and hemorrhagic stroke. One locus for autosomal dominant CCM ( CCM1 ) maps to chromosome 7q21-q22. Recombination events in linked family members define a critical region of approximately 2 Mb and a shared disease haplotype associated with a presumed founder effect in families of Mexican-American descent points to a potentially smaller region of interest. Using a genomic sequence-based positional cloning strategy, we have identified KRIT1, encoding a protein that interacts with the Krev-1/rap1a tumor suppressor, as the CCM1 gene. Seven different KRIT1 mutations have been identified in 23 distinct CCM1 families. The identical mutation is present in 16 of 21 Mexican-American families analyzed, substantiating a founder effect in this population. Other Mexican-American and non-Hispanic Caucasian CCM1 kindreds harbor other KRIT1 mutations. Identification of a common Mexican-American mutation has potential clinical significance for presymptomatic diagnosis of CCM in this population. In addition, these data point to a key role for the Krev-1/rap1a signaling pathway in angiogenesis and cerebrovascular disease.

Characterization of the human and mouse unconventional myosin XV genes responsible for hereditary deafness DFNB3 and shaker 2.

Mutations in myosin XV are responsible for congenital profound deafness DFNB3 in humans and deafness and vestibular defects in shaker 2 mice. By combining direct cDNA analyses with a comparison of 95.2 kb of genomic DNA sequence from human chromosome 17p11.2 and 88.4 kb from the homologous region on mouse chromosome 11, we have determined the genomic and mRNA structures of the human (MYO15) and mouse (Myo15) myosin XV genes. Our results indicate that full-length myosin XV transcripts contain 66 exons, are >12 kb in length, and encode 365-kDa proteins that are unique among myosins in possessing very long approximately 1200-aa N-terminal extensions preceding their conserved motor domains. The tail regions of the myosin XV proteins contain two MyTH4 domains, two regions with similarity to the membrane attachment FERM domain, and a putative SH3 domain. Northern and dot blot analyses revealed that myosin XV is expressed in the pituitary gland in both humans and mice. Myosin XV transcripts were also observed by in situ hybridization within areas corresponding to the sensory epithelia of the cochlea and vestibular systems in the developing mouse inner ear. Immunostaining of adult mouse organ of Corti revealed that myosin XV protein is concentrated within the cuticular plate and stereocilia of cochlear sensory hair cells. These results indicate a likely role for myosin XV in the formation or maintenance of the unique actin-rich structures of inner ear sensory hair cells.

Identification and detection of the common 65-kb deletion breakpoint in the nephropathic cystinosis gene (CTNS).

The most common mutation in the cystinosis gene, CTNS, is a 65-kb deletion thought to have originated in Germany. Although homozygotes for this deletion are detectable by the absence of the D17S829 polymorphic marker, no method exists to identify heterozygotes. We identified the 65-kb deletion breakpoints and used flanking PCR primers to amplify a 423-bp fragment present only in the deletion alleles. Using this method, we determined that 121 of 216 (56%) cystinosis alleles examined bore the 65-kb deletion. We found no non-Europeans with the deletion, and the deletion size and breakpoints appeared identical in all patients studied, supporting the concept of a founder effect. The addition of D17S829 primers (266 bp apart) to the PCR created a multiplex PCR system useful for diagnosing cystinosis patients homozygous and heterozygous for the 65-kb deletion.

An effective approach for analyzing "prefinished" genomic sequence data.

Ongoing efforts to sequence the human genome are already generating large amounts of data, with substantial increases anticipated over the next few years. In most cases, a shotgun sequencing strategy is being used, which rapidly yields most of the primary sequence in incompletely assembled sequence contigs ("prefinished" sequence) and more slowly produces the final, completely assembled sequence ("finished" sequence). Thus, in general, prefinished sequence is produced in excess of finished sequence, and this trend is certain to continue and even accelerate over the next few years. Even at a prefinished stage, genomic sequence represents a rich source of important biological information that is of great interest to many investigators. However, analyzing such data is a challenging and daunting task, both because of its sheer volume and because it can change on a day-by-day basis. To facilitate the discovery and characterization of genes and other important elements within prefinished sequence, we have developed an analytical strategy and system that uses readily available software tools in new combinations. Implementation of this strategy for the analysis of prefinished sequence data from human chromosome 7 has demonstrated that this is a convenient, inexpensive, and extensible solution to the problem of analyzing the large amounts of preliminary data being produced by large-scale sequencing efforts. Our approach is accessible to any investigator who wishes to assimilate additional information about particular sequence data en route to developing richer annotations of a finished sequence.

Correction of deafness in shaker-2 mice by an unconventional myosin in a BAC transgene.

The shaker-2 mouse mutation, the homolog of human DFNB3, causes deafness and circling behavior. A bacterial artificial chromosome (BAC) transgene from the shaker-2 critical region corrected the vestibular defects, deafness, and inner ear morphology of shaker-2 mice. An unconventional myosin gene, Myo15, was discovered by DNA sequencing of this BAC. Shaker-2 mice were found to have an amino acid substitution at a highly conserved position within the motor domain of this myosin. Auditory hair cells of shaker-2 mice have very short stereocilia and a long actin-containing protrusion extending from their basal end. This histopathology suggests that Myo15 is necessary for actin organization in the hair cells of the cochlea.

Association of unconventional myosin MYO15 mutations with human nonsyndromic deafness DFNB3.

DFNB3, a locus for nonsyndromic sensorineural recessive deafness, maps to a 3-centimorgan interval on human chromosome 17p11.2, a region that shows conserved synteny with mouse shaker-2. A human unconventional myosin gene, MYO15, was identified by combining functional and positional cloning approaches in searching for shaker-2 and DFNB3. MYO15 has at least 50 exons spanning 36 kilobases. Sequence analyses of these exons in affected individuals from three unrelated DFNB3 families revealed two missense mutations and one nonsense mutation that cosegregated with congenital recessive deafness.

A physical map of human chromosome 7: an integrated YAC contig map with average STS spacing of 79 kb.

The construction of highly integrated and annotated physical maps of human chromosomes represents a critical goal of the ongoing Human Genome Project. Our laboratory has focused on developing a physical map of human chromosome 7, a approximately 170-Mb segment of DNA that corresponds to an estimated 5% of the human genome. Using a yeast artificial chromosome (YAC)-based sequence-tagged site (STS)-content mapping strategy, 2150 chromosome 7-specific STSs have been established and mapped to a collection of YACs highly enriched for chromosome 7 DNA. The STSs correspond to sequences generated from a variety of DNA sources, with particular emphasis placed on YAC insert ends, genetic markers, and genes. The YACs include a set of relatively nonchimeric clones from a human-hamster hybrid cell line as well as clones isolated from total genomic libraries. For map integration, we have localized 260 STSs corresponding to Genethon genetic markers and 259 STSs corresponding to markers orders by radiation hybrid (RH) mapping on our YAC contigs. Analysis of the data with the program SEGMAP results in the assembly of 22 contigs that are "anchored" on the Genethon genetic map, the RH map, and/or the cytogenetic map. These 22 contigs are ordered relative to one another, are (in all but 3 cases) oriented relative to the centromere and telomeres, and contain > 98% of the mapped STSs. The largest anchored YAC contig, accounting for most of 7p, contains 634 STSs and 1260 YACs. An additional 14 contigs, accounting for approximately 1.5% of the mapped STSs, are assembled but remain unanchored on either the genetic or RH map. Therefore, these 14 "orphan" contigs are not ordered relative to other contigs. In our contig maps, adjacent STSs are connected by two or more YACs in > 95% of cases. With 2150 mapped STSs, our map provides an average STS spacing of approximately 79 kb. The physical map we report here exceeds the goal of 100-kb average STS spacing and should provide an excellent framework for systematic sequencing of the chromosome.

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.

Branchpoint and polypyrimidine tract mutations mediating the loss and partial recovery of the Moloney murine sarcoma virus MuSVts110 thermosensitive splicing phenotype.

Balanced splicing of retroviral RNAs is mediated by weak signals at the 3' splice site (ss) acting in concert with other cis elements. Moloney murine sarcoma virus MuSVts110 shows a similar balance between unspliced and spliced RNAs, differing only in that the splicing of its RNA is, in addition, growth temperature sensitive. We have generated N-nitroso-N-methylurea (NMU)-treated MuSVts110 revertants in which splicing was virtually complete at all temperatures and have investigated the molecular basis of this reversion on the assumption that the findings would reveal cis-acting elements controlling MuSVts110 splicing thermosensitivity. In a representative revertant (NMU-20), we found that complete splicing was conferred by a G-to-A substitution generating a consensus branchpoint (BP) signal (-CCCUGGC- to -CCCUGAC- [termed G(-25)A]) at -25 relative to the 3' ss. Weakening this BP to -CCCGAC- [G(-25)A,U(-27)C] moderately reduced splicing at the permissive temperature and sharply inhibited splicing at the originally nonpermissive temperature, arguing that MuSVts110 splicing thermosensitivity depends on a suboptimal BP-U2 small nuclear RNA interaction. This conclusion was supported by results indicating that lengthening the short MuSVts110 polypyrimidine tract and altering its uridine content doubled splicing efficiency at permissive temperatures and nearly abrogated splicing thermosensitivity. In vitro splicing experiments showed that MuSVts110 G(-25)A RNA intermediates were far more efficiently ligated than RNAs carrying the wild-type BP, the G(-25)A,U (-27)C BP, or the extended polypyrimidine tract. The efficiency of ligation in vitro roughly paralleled splicing efficiency in vivo [G(-25)A BP > extended polypyrimidine tract > G(-25)A,U(-27)C BP > wild-type BP]. These results suggest that MuSVts110 RNA splicing is balanced by cis elements similar to those operating in other retroviruses and, in addition, that its splicing thermosensitivity is a response to the presence of multiple suboptimal splicing signals.