Genomic sequence analysis of Fugu rubripes CFTR and flanking genes in a 60 kb region conserving synteny with 800 kb of human chromosome 7

Exploiting species differences to understand the CFTR Cl- channel

The anion channel cystic fibrosis transmembrane conductance regulator (CFTR) is a unique ATP-binding cassette (ABC) transporter. CFTR plays a pivotal role in transepithelial ion transport as its dysfunction in the genetic disease cystic fibrosis (CF) dramatically demonstrates. Phylogenetic analysis suggests that CFTR first appeared in aquatic vertebrates fulfilling important roles in osmosensing and organ development. Here, we review selectively, knowledge of CFTR structure, function and pharmacology, gleaned from cross-species comparative studies of recombinant CFTR proteins, including CFTR chimeras. The data argue that subtle changes in CFTR structure can affect strongly channel function and the action of CF mutations.

Transcription-dependent spatial arrangements of CFTR and conserved adjacent loci are not conserved in human and murine nuclei

The human genes CFTR, ASZ1/GASZ, and CTTNBP2/CORTBP2 map to adjacent loci on chromosome 7q31 and display characteristic patterns of nuclear positioning, which strictly correlate with the state of activity. To address the evolutionary conservation of gene positioning, we investigated transcriptional activity and nuclear positioning of the highly conserved murine orthologs and of additional murine genes mapping to the region of conserved synteny on mouse chromosome 6. The results showed that all murine loci investigated constitutively localized in the nuclear interior irrespective of their functional state. Silenced loci did not display preferential association with the nuclear periphery or with chromocenters, respectively, and no differential positioning with respect to the chromosome 6 territory could be observed. This positional behavior of the murine loci was in striking contrast to the positioning of the human orthologs, and the results show that the transcription-dependent positioning of CFTR and adjacent loci has not been conserved. The findings reveal that the nuclear organization of conserved chromosomal regions can change rapidly during evolution and is not always as highly conserved as other features of chromosome organization. Furthermore, the results suggest that the way how nuclear positioning contributes to the regulation of conserved loci can be different in different vertebrate species.

Short fuzzy tandem repeats in genomic sequences, identification, and possible role in regulation of gene expression

MOTIVATION

Genomic sequences are highly redundant and contain many types of repetitive DNA. Fuzzy tandem repeats (FTRs) are of particular interest. They are found in regulatory regions of eukaryotic genes and are reported to interact with transcription factors. However, accurate assessment of FTR occurrences in different genome segments requires specific algorithm for efficient FTR identification and classification.

RESULTS

We have obtained formulas for P-values of FTR occurrence and developed an FTR identification algorithm implemented in TandemSWAN software. Using TandemSWAN we compared the structure and the occurrence of FTRs with short period length (up to 24 bp) in coding and non-coding regions including UTRs, heterochromatic, intergenic and enhancer sequences of Drosophila melanogaster and Drosophila pseudoobscura. Tandems with period three and its multiples were found in coding segments, whereas FTRs with periods multiple of six are overrepresented in all non-coding segment. Periods equal to 5-7 and 11-14 were characteristic of the enhancer regions and other non-coding regions close to genes.

AVAILABILITY

TandemSWAN web page, stand-alone version and documentation can be found at http://bioinform.genetika.ru/projects/swan/www/

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

The Multifunctional Fish Gill: Dominant Site of Gas Exchange, Osmoregulation, Acid-Base Regulation, and Excretion of Nitrogenous Waste

The fish gill is a multipurpose organ that, in addition to providing for aquatic gas exchange, plays dominant roles in osmotic and ionic regulation, acid-base regulation, and excretion of nitrogenous wastes. Thus, despite the fact that all fish groups have functional kidneys, the gill epithelium is the site of many processes that are mediated by renal epithelia in terrestrial vertebrates. Indeed, many of the pathways that mediate these processes in mammalian renal epithelial are expressed in the gill, and many of the extrinsic and intrinsic modulators of these processes are also found in fish endocrine tissues and the gill itself. The basic patterns of gill physiology were outlined over a half century ago, but modern immunological and molecular techniques are bringing new insights into this complicated system. Nevertheless, substantial questions about the evolution of these mechanisms and control remain.

Conserved and acquired features of neurogenin1 regulation

The telencephalon shows vast morphological variations among different vertebrate groups. The transcription factor neurogenin1(ngn1) controls neurogenesis in the mouse pallium and is also expressed in the dorsal telencephalon of the evolutionary distant zebrafish. The upstream regions of the zebrafish and mammalian ngn1 loci harbour several stretches of conserved sequences. Here, we show that the upstream region of zebrafish ngn1 is capable of faithfully recapitulating endogenous expression in the zebrafish and mouse telencephalon. A single conserved regulatory region is essential for dorsal telencephalic expression in the zebrafish, and for expression in the dorsal pallium of the mouse. However, a second conserved region that is inactive in the fish telencephalon is necessary for expression in the lateral pallium of mouse embryos. This regulatory region, which drives expression in the zebrafish diencephalon and hindbrain, is dependent on Pax6 activity and binds recombinant Pax6 in vitro. Thus, the regulatory elements of ngn1 appear to be conserved among vertebrates, with certain differences being incorporated in the utilisation of these enhancers, for the acquisition of more advanced features in amniotes. Our data provide evidence for the co-option of regulatory regions as a mechanism of evolutionary diversification of expression patterns, and suggest that an alteration in Pax6expression was crucial in neocortex evolution.

Theatre: A software tool for detailed comparative analysis and visualization of genomic sequence

Theatre is a web-based computing system designed for the comparative analysis of genomic sequences, especially with respect to motifs likely to be involved in the regulation of gene expression. Theatre is an interface to commonly used sequence analysis tools and biological sequence databases to determine or predict the positions of coding regions, repetitive sequences and transcription factor binding sites in families of DNA sequences. The information is displayed in a manner that can be easily understood and can reveal patterns that might not otherwise have been noticed. In addition to web-based output, Theatre can produce publication quality colour hardcopies showing predicted features in aligned genomic sequences. A case study using the p53 promoter region of four mammalian species and two fish species is described. Unlike the mammalian sequences the promoter regions in fish have not been previously predicted or characterized and we report the differences in the p53 promoter region of four mammals and that predicted for two fish species. Theatre can be accessed at http://www.hgmp.mrc.ac.uk/Registered/Webapp/theatre/.

A comparative genomic analysis of the cow, pig, and human CFTR genes identifies potential intronic regulatory elements

The identification of sequences within noncoding regions of genes that are conserved between several species may indicate potential regulatory elements. This is important for genes with complex control mechanisms such as the cystic fibrosis transmembrane conductance regulator (CFTR). CFTR demonstrates similar patterns of temporal and spatial expression in human and sheep, but these differ significantly in mouse cftr. The complete sheep CFTR sequence is unavailable so we annotated BAC clones encompassing the CFTR gene from two other artiodactyl species (cow and pig) for comparative sequence analysis. Regions of introns 2, 3, 10, 17a, 18, and 21 and 3' flanking sequence corresponding to human CFTR DNase I hypersensitive sites (DHS) showed high homology in the cow and pig. Cross-species sequence conservation also enabled finer mapping of other human DHS, including those in introns 1, 16, and 20. Additional potential regulatory elements not associated with human DHS were also identified.

Evolutionary constraints on the Disrupted in Schizophrenia locus

The Disrupted in Schizophrenia (DISC) locus on human chromosome 1q42 has been strongly implicated by genetic studies as a susceptibility locus for major mental illnesses. In humans the locus is transcriptionally complex, with multiple alternate splicing events, antisense transcription, and intergenic splicing all evident. We have compared the genomic sequence and transcription maps of this locus between human, mouse, pufferfish (Fugu rubripes), and, in part, zebrafish (Danio rerio). The order and orientation of EGLN1, TSNAX, and DISC1 genes are conserved between mammals and F. rubripes. Intergenic splicing and short intergenic transcripts are not found to be conserved features. DISC2, a putative noncoding transcript partially antisense to DISC1, is not conserved in mouse or F. rubripes. Alternate splice forms of the protein-coding DISC1 gene are conserved even though the genomic structure is not. The amino acid sequence of DISC1 is diverging rapidly, although a putative nuclear localization signal and discrete blocks of coiled coil are specifically conserved features.

Cystic fibrosis transmembrane conductance regulator in teleost fish

The gills and intestinal epithelia of teleost fish express cystic fibrosis transmembrane conductance regulator (CFTR), and utilize this low conductance anion channel in the apical membrane for ion secretion in seawater gill and in the basolateral membrane for ion absorption in freshwater gill. Similarly, in the intestine CFTR is present in the basolateral membrane for intestinal absorption and also in the apical membrane of secreting intestine. The expression of CFTR and the directed trafficking of the protein to the apical or basolateral membrane is salinity-dependent. The CFTR gene has been cloned and sequenced from several teleost species and although all the major elements in the human gene are present, including two nucleotide binding domains that are common to all ATP binding cassette (ABC) transporters, the sequences are divergent compared to shark or human. In euryhaline fish adapting to seawater, CFTR, localized immunocytochemically, redistributes slowly from a basolateral location to the apical membrane while ion secretory capacity increases. The facility with which teleosts regulate CFTR expression and activation during salinity adaptation make this system an appealing model for the expression and trafficking operation of this labile gene product.

Na(+), Cl(-), Ca(2+) and Zn(2+) transport by fish gills: retrospective review and prospective synthesis

The secondary active Cl(-) secretion in seawater (SW) teleost fish gills and elasmobranch rectal gland involves basolateral Na(+),K(+)-ATPase and NKCC, apical membrane CFTR anion channels, and a paracellular Na(+)-selective conductance. In freshwater (FW) teleost gill, the mechanism of NaCl uptake is more controversial and involves apical V-type H(+)-ATPase linked to an apical Na(+) channel, apical Cl(-)-HCO-3 exchange and basolateral Na(+),K(+)-ATPase. Ca(2+) uptake (in FW and SW) is via Ca(2+) channels in the apical membrane and Ca(2+)-ATPase in the basolateral membrane. Mainly this transport occurs in mitochondria rich (MR) chloride cells, but there is a role for the pavement cells also. Future research will likely expand in two major directions, molded by methodology: first in physiological genomics of all the transporters, including their expression, trafficking, operation, and regulation at the molecular level, and second in biotelemetry to examine multivariable components in behavioral physiological ecology, thus widening the integration of physiology from the molecular to the environmental levels while deepening understanding at all levels.

Genomic organization, transcript variants and comparative analysis of the human nucleoporin 155 (NUP155) gene

Nucleoporin 155 (Nup155) is a major component of the nuclear pore complex (NPC) involved in cellular nucleo-cytoplasmic transport. We have acquired the complete sequence and interpreted the genomic organization of the Nup155 orthologos from human (Homo sapiens) and pufferfish (Fugu rubripes), which are approximately 80 and 8 kb in length, respectively. The human gene is ubiquitously expressed in many tissues analyzed and has two major transcript variants, resulted from an alternative usage of the 5' cryptic or consensus splice donor in intron 1 and two polyadenylation signals. We have also cloned DNA complementary to RNAs of the Nup155 orthologs from Fugu and mouse. Comparative analysis of the Nup155 orthologs in many species, including H. sapiens, Mus musculus, Rattus norvegicus, F. rubripes, Arabidopsis thaliana, Drosophila melanogaster, and Saccharomyces cerevisiae, has revealed two paralogs in S. cerevisiae but only a single gene with increasing number of introns in more complex organisms. The amino acid sequences of the Nup155 orthologos are highly conserved in the evolution of eukaryotes. Different gene orders in the human and Fugu genomic regions harboring the Nup155 orthologs advocate cautious interpretation of synteny in comparative genomic analysis even within the vertebrate lineage.

Vertebrate genome sequencing: building a backbone for comparative genomics

The human genome sequence provides a reference point from which we can compare ourselves with other organisms. Interspecies comparison is a powerful tool for inferring function from genomic sequence and could ultimately lead to the discovery of what makes humans unique. To date, most comparative sequencing has focused on pair-wise comparisons between human and a limited number of other vertebrates, such as mouse. Targeted approaches now exist for mapping and sequencing vertebrate bacterial artificial chromosomes (BACs) from numerous species, allowing rapid and detailed molecular and phylogenetic investigation of multi-megabase loci. Such targeted sequencing is complementary to current whole-genome sequencing projects, and would benefit greatly from the creation of BAC libraries from a diverse range of vertebrates.

Comparative genomics and regulatory evolution: conservation and function of the Chs and Apetala3 promoters

DNA sequence variations of chalcone synthase (Chs) and Apetala3 gene promoters from 22 cruciferous plant species were analyzed to identify putative conserved regulatory elements. Our comparative approach confirmed the existence of numerous conserved sequences which may act as regulatory elements in both investigated promoters. To confirm the correct identification of a well-conserved UV-light-responsive promoter region, a subset of Chs promoter fragments were tested in Arabidopsis thaliana protoplasts. All promoters displayed similar light responsivenesses, indicating the general functional relevance of the conserved regulatory element. In addition to known regulatory elements, other highly conserved regions were detected which are likely to be of functional importance. Phylogenetic trees based on DNA sequences from both promoters (gene trees) were compared with the hypothesized phylogenetic relationships (species trees) of these taxa. The data derived from both promoter sequences were congruent with the phylogenies obtained from coding regions of other nuclear genes and from chloroplast DNA sequences. This indicates that promoter sequence evolution generally is reflective of species phylogeny. Our study also demonstrates the great value of comparative genomics and phylogenetics as a basis for functional analysis of promoter action and gene regulation.

Genomic organization and expression of the doublesex-related gene cluster in vertebrates and detection of putative regulatory regions for DMRT1

Genes related to the Drosophila melanogaster doublesex and Caenorhabditis elegans mab-3 genes are conserved in human. They are identified by a DNA-binding homology motif, the DM domain, and constitute a gene family (DMRTs). Unlike the invertebrate genes, whose role in the sex-determination process is essentially understood, the function of the different vertebrate DMRT genes is not as clear. Evidence has accumulated for the involvement of DMRT1 in male sex determination and differentiation. DMRT2 (known as terra in zebrafish) seems to be a critical factor for somitogenesis. To contribute to a better understanding of the function of this important gene family, we have analyzed DMRT1, DMRT2, and DMRT3 from the genome model organism Fugu rubripes and the medakafish, a complementary model organism for genetics and functional studies. We found conservation of synteny of human chromosome 9 in F. rubripes and an identical gene cluster organization of the DMRTs in both fish. Although expression analysis and gene linkage mapping in medaka exclude a function for any of the three genes in the primary step of male sex determination, comparison of F. rubripes and human sequences uncovered three putative regulatory regions that might have a role in more downstream events of sex determination and human XY sex reversal.

A combined analysis of the cystic fibrosis transmembrane conductance regulator: implications for structure and disease models

Over the past decade, nearly 1,000 variants have been identified in the cystic fibrosis transmembrane conductance regulator (CFTR) gene in classic and atypical cystic fibrosis (CF) patients worldwide, and an enormous wealth of information concerning the structure and function of the protein has also been accumulated. These data, if evaluated together in a sequence comparison of all currently available CFTR homologs, are likely to refine the global structure-function relationship of the protein, which will, in turn, facilitate interpretation of the identified mutations in the gene. Based on such a combined analysis, we had recently defined a "functional R domain" of the CFTR protein. First, presenting two full-length cDNA sequences (termed sCFTR-I and sCFTR-II) from the Atlantic salmon (Salmo salar) and an additional partial coding sequence from the eastern gray kangaroo (Macropus giganteus), this study went further to refine the boundaries of the two nucleotide-binding domains (NBDs) and the COOH-terminal tail (C-tail), wherein NBD1 was defined as going from P439 to G646, NBD2 as going from A1225 to E1417, and the C-tail as going from E1418 to L1480. This approach also provided further insights into the differential roles of the two halves of CFTR and highlighted several well-conserved motifs that may be involved in inter- or intramolecular interactions. Moreover, a serious concern that a certain fraction of missense mutations identified in the CFTR gene may not have functional consequences was raised. Finally, phylogenetic analysis of all the full-length CFTR amino acid sequences and an extended set of exon 13--coding nucleotide sequences reinforced the idea that the rabbit may represent a better CF model than the mouse and strengthened the assertion that a long-branch attraction artifact separates the murine rodents from the rabbit and the guinea pig, the other Glires.

The sheep genome contributes to localization of control elements in a human gene with complex regulatory mechanisms

Genes that show complex tissue-specific and temporal control by regulatory elements located outside their promoters present a considerable challenge to identify the sequences involved. The rapid accumulation of genomic sequence information for a number of species has enabled a comparative phylogenetic approach to find important regulatory elements. For some genes, which show a similar pattern of expression in humans and rodents, genomic sequence information for these two species may be sufficient. Others, such as the cystic fibrosis transmembrane conductance regulator (CFTR) gene, show significant divergence in expression patterns between mouse and human, necessitating phylogenetic approaches involving additional species. The ovine CFTR gene has a temporal and spatial expression pattern that is very similar to that of human CFTR. Comparative genomic sequence analysis of ovine and human CFTR identified high levels of homology between the core elements in several potential regulatory elements defined as DNase I hypersensitive sites in human CFTR. These data provide a case for the power of an artiodactyl genome to contribute to the understanding of human genetic disease.

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.