Large-scale comparative sequence analysis of the human and murine Bruton's tyrosine kinase loci reveals conserved regulatory domains

Comparison of inhibitory effects of irreversible and reversible Btk inhibitors on platelet function

All irreversible Bruton tyrosine kinase (Btk) inhibitors including ibrutinib and acalabrutinib induce platelet dysfunction and increased bleeding risk. New reversible Btk inhibitors were developed, like MK-1026. The mechanism underlying increased bleeding tendency with Btk inhibitors remains unclear. We investigated the effects of ibrutinib, acalabrutinib and MK-1026 on platelet function in healthy volunteers, patients and Btk-deficient mice, together with off-target effects on tyrosine kinase phosphorylation. All inhibitors suppressed GPVI- and CLEC-2-mediated platelet aggregation, activation and secretion in a dose-dependent manner. Only ibrutinib inhibited thrombus formation on vWF-co-coated surfaces, while on collagen this was not affected. In blood from Btk-deficient mice, collagen-induced thrombus formation under flow was reduced, but preincubation with either inhibitor was without additional effects. MK-1026 showed less off-target effects upon GPVI-induced TK phosphorylation as compared to ibrutinib and acalabrutinib. In ibrutinib-treated patients, GPVI-stimulated platelet activation, and adhesion on vWF-co-coated surfaces were inhibited, while CLEC-2 stimulation induced variable responses. The dual inhibition of GPVI and CLEC-2 signalling by Btk inhibitors might account for the increased bleeding tendency, with ibrutinib causing more high-grade bleedings due to additional inhibition of platelet-vWF interaction. As MK-1026 showed less off-target effects and only affected activation of isolated platelets, it might be promising for future treatment.

Effective, safe, and sustained correction of murine XLA using a UCOE-BTK promoter-based lentiviral vector

X-linked agammaglobulinemia (XLA) is an immune disorder caused by mutations in Bruton’s tyrosine kinase (BTK). BTK is expressed in B and myeloid cells, and its deficiency results in a lack of mature B cells and protective antibodies. We previously reported a lentivirus (LV) BTK replacement therapy that restored B cell development and function in Btk and Tec double knockout mice (a phenocopy of human XLA). In this study, with the goal of optimizing both the level and lineage specificity of BTK expression, we generated LV incorporating the proximal human BTK promoter. Hematopoietic stem cells from Btk^−/−Tec^−/− mice transduced with this vector rescued lineage-specific expression and restored B cell function in Btk^−/−Tec^−/− recipients. Next, we tested addition of candidate enhancers and/or ubiquitous chromatin opening elements (UCOEs), as well as codon optimization to improve BTK expression. An Eμ enhancer improved B cell rescue, but increased immunoglobulin G (IgG) autoantibodies. Addition of the UCOE avoided autoantibody generation while improving B cell development and function and reducing vector silencing. An optimized vector containing a truncated UCOE upstream of the BTK promoter and codon-optimized BTK cDNA resulted in stable, lineage-regulated BTK expression that mirrored endogenous BTK, making it a strong candidate for XLA therapy.

A genome-wide cis-regulatory element discovery method based on promoter sequences and gene co-expression networks

Background

Deciphering cis-regulatory networks has become an attractive yet challenging task. This paper presents a simple method for cis-regulatory network discovery which aims to avoid some of the common problems of previous approaches.

Results

Using promoter sequences and gene expression profiles as input, rather than clustering the genes by the expression data, our method utilizes co-expression neighborhood information for each individual gene, thereby overcoming the disadvantages of current clustering based models which may miss specific information for individual genes. In addition, rather than using a motif database as an input, it implements a simple motif count table for each enumerated k-mer for each gene promoter sequence. Thus, it can be used for species where previous knowledge of cis-regulatory motifs is unknown and has the potential to discover new transcription factor binding sites. Applications on Saccharomyces cerevisiae and Arabidopsis have shown that our method has a good prediction accuracy and outperforms a phylogenetic footprinting approach. Furthermore, the top ranked gene-motif regulatory clusters are evidently functionally co-regulated, and the regulatory relationships between the motifs and the enriched biological functions can often be confirmed by literature.

Conclusions

Since this method is simple and gene-specific, it can be readily utilized for insufficiently studied species or flexibly used as an additional step or data source for previous transcription regulatory networks discovery models.

Development of B-lineage predominant lentiviral vectors for use in genetic therapies for B cell disorders

Sustained, targeted, high-level transgene expression in primary B lymphocytes may be useful for gene therapy in B cell disorders. We developed several candidate B-lineage predominant self-inactivating lentiviral vectors (LV) containing alternative enhancer/promoter elements including: the immunoglobulin β (Igβ) (B29) promoter combined with the immunoglobulin µ enhancer (EµB29); and the endogenous BTK promoter with or without Eµ (EµBtkp or Btkp). LV-driven enhanced green fluorescent protein (eGFP) reporter expression was evaluated in cell lines and primary cells derived from human or murine hematopoietic stem cells (HSC). In murine primary cells, EµB29 and EµBtkp LV-mediated high-level expression in immature and mature B cells compared with all other lineages. Expression increased with B cell maturation and was maintained in peripheral subsets. Expression in T and myeloid cells was much lower in percentage and intensity. Similarly, both EµB29 and EµBtkp LV exhibited high-level activity in human primary B cells. In contrast to EµB29, Btkp and EµBtkp LV also exhibited modest activity in myeloid cells, consistent with the expression profile of endogenous Bruton's tyrosine kinase (Btk). Notably, EµB29 and EµBtkp activity was superior in all expression models to an alternative, B-lineage targeted vector containing the EµS.CD19 enhancer/promoter. In summary, EµB29 and EµBtkp LV comprise efficient delivery platforms for gene expression in B-lineage cells.

A systems approach to mapping transcriptional networks controlling surfactant homeostasis

Background

Pulmonary surfactant is required for lung function at birth and throughout life. Lung lipid and surfactant homeostasis requires regulation among multi-tiered processes, coordinating the synthesis of surfactant proteins and lipids, their assembly, trafficking, and storage in type II cells of the lung. The mechanisms regulating these interrelated processes are largely unknown.

Results

We integrated mRNA microarray data with array independent knowledge using Gene Ontology (GO) similarity analysis, promoter motif searching, protein interaction and literature mining to elucidate genetic networks regulating lipid related biological processes in lung. A Transcription factor (TF) - target gene (TG) similarity matrix was generated by integrating data from different analytic methods. A scoring function was built to rank the likely TF-TG pairs. Using this strategy, we identified and verified critical components of a transcriptional network directing lipogenesis, lipid trafficking and surfactant homeostasis in the mouse lung.

Conclusions

Within the transcriptional network, SREBP, CEBPA, FOXA2, ETSF, GATA6 and IRF1 were identified as regulatory hubs displaying high connectivity. SREBP, FOXA2 and CEBPA together form a common core regulatory module that controls surfactant lipid homeostasis. The core module cooperates with other factors to regulate lipid metabolism and transport, cell growth and development, cell death and cell mediated immune response. Coordinated interactions of the TFs influence surfactant homeostasis and regulate lung function at birth.

The correlations of the function and positional distribution of the cis-elements CArG around the TSS in the genes of Mus musculus

While many studies of cis-elements CArG bound by serum response factor (SRF) are in progress, little is known about the positional distribution of the functional CArG elements around the transcription start site (TSS) of genes that they influence. We use a validated CArG data set to calculate the distance distribution of functional CArG elements around the TSS. Distances between adjacent CArGs were also analyzed. We compare these distributions with those derived using a control set of randomly selected CArGs (that were not experimentally validated for function). Our results show that most functional CArG elements (108 of 152, 71%) exist upstream of the annotated TSS, with copy number increasing as one moves closer to the TSS. Moreover, the average number of the CArG elements in the CArG-containing genes is significantly more than that in the control genes. Our study extends earlier bioinformatic analyses of functional CArG elements and provides an application of comparative sequence data to the identification of transcription factor binding sites.

Genomic promoter analysis predicts functional transcription factor binding

Background. The computational identification of functional transcription factor binding sites (TFBSs) remains a major challenge of computational biology. Results. We have analyzed the conserved promoter sequences for the complete set of human RefSeq genes using our conserved transcription factor binding site (CONFAC) software. CONFAC identified 16296 human-mouse ortholog gene pairs, and of those pairs, 9107 genes contained conserved TFBS in the 3 kb proximal promoter and first intron. To attempt to predict in vivo occupancy of transcription factor binding sites, we developed a novel marginal effect isolator algorithm that builds upon Bayesian methods for multigroup TFBS filtering and predicted the in vivo occupancy of two transcription factors with an overall accuracy of 84%. Conclusion. Our analyses show that integration of chromatin immunoprecipitation data with conserved TFBS analysis can be used to generate accurate predictions of functional TFBS. They also show that TFBS cooccurrence can be used to predict transcription factor binding to promoters in vivo.

PipMaker: a World Wide Web server for genomic sequence alignments

PipMaker is a World-Wide Web site used to compare two long genomic sequences and identify conserved segments between them. This unit describes the use of the PipMaker server and explains the resulting output files. PipMaker provides an efficient method of aligning genomic sequences and returns a compact, but easy-to-interpret form of output, the percent identity plot (pip). For each aligning segment between two sequences the pip shows both the position relative to the first sequence and the degree of similarity. Optional annotations on the pip provide additional information to assist in the interpretation of the alignment. The default parameters of the underlying blastz alignment program are tuned for human-mouse alignments.

Proteasome-dependent autoregulation of Bruton tyrosine kinase (Btk) promoter via NF-kappaB

Bruton tyrosine kinase (Btk) is critical for B-cell development. Btk regulates a plethora of signaling proteins, among them nuclear factor-[kappa]B (NF-kappaB). Activation of NF-kappaB is a hallmark of B cells, and NF-kappaB signaling is severely compromised in Btk deficiency. We here present strong evidence indicating that NF-kappaB is required for efficient transcription of the Btk gene. First, we found that proteasome blockers and inhibitors of NF-kappaB signaling suppress Btk transcription and intracellular expression. Similar to Btk, proteasome inhibitors also reduced the expression of other members of this family of kinases, Itk, Bmx, and Tec. Second, 2 functional NF-kappaB-binding sites were found in the Btk promoter. Moreover, in live mice, by hydrodynamic transfection, we show that bortezomib (a blocker of proteasomes and NF-kappaB signaling), as well as NF-kappaB binding sequence-oligonucleotide decoys block Btk transcription. We also demonstrate that Btk induces NF-kappaB activity in mice. Collectively, we show that Btk uses a positive autoregulatory feedback mechanism to stimulate transcription from its own promoter via NF-kappaB.

Comparative gene mapping in cattle, Indian muntjac, and Chinese muntjac by fluorescence in situ hybridization

The Indian muntjac (Muntiacus muntjak vaginalis) has a karyotype of 2n = 6 in the female and 2n = 7 in the male. The karyotypic evolution of Indian muntjac via extensive tandem fusions and several centric fusions are well documented by molecular cytogenetic studies mainly utilizing chromosome paints. To achieve higher resolution mapping, a set of 42 different genomic clones coding for 37 genes and the nucleolar organizer region were used to examine homologies between the cattle (2n = 60), human (2n = 46), Indian muntjac (2n = 6/7) and Chinese muntjac (2n = 46) karyotypes. These genomic clones were mapped by fluorescence in situ hybridization (FISH). Localization of genes on all three pairs of M. m. vaginalis chromosomes and on the acrocentric chromosomes of M. reevesi allowed not only the analysis of the evolution of syntenic regions within the muntjac genus but also allowed a broader comparison of synteny with more distantly related species, such as cattle and human, to shed more light onto the evolving genome organization.

Organization of conserved elements near key developmental regulators in vertebrate genomes

Sequence conservation has traditionally been used as a means to target functional regions of complex genomes. In addition to its use in identifying coding regions of genes, the recent availability of whole genome data for a number of vertebrates has permitted high-resolution analyses of the noncoding "dark matter" of the genome. This has resulted in the identification of a large number of highly conserved sequence elements that appear to be preserved in all bony vertebrates. Further positional analysis of these conserved noncoding elements (CNEs) in the genome demonstrates that they cluster around genes involved in developmental regulation. This chapter describes the identification and characterization of these elements, with particular reference to their composition and organization.

Experimental validation of novel genes predicted in the un-annotated regions of the Arabidopsis genome

Background

Several lines of evidence support the existence of novel genes and other transcribed units which have not yet been annotated in the Arabidopsis genome. Two gene prediction programs which make use of comparative genomic analysis, Twinscan and EuGene, have recently been deployed on the Arabidopsis genome. The ability of these programs to make use of sequence data from other species has allowed both Twinscan and EuGene to predict over 1000 genes that are intergenic with respect to the most recent annotation release. A high throughput RACE pipeline was utilized in an attempt to verify the structure and expression of these novel genes.

Results

1,071 un-annotated loci were targeted by RACE, and full length sequence coverage was obtained for 35% of the targeted genes. We have verified the structure and expression of 378 genes that were not present within the most recent release of the Arabidopsis genome annotation. These 378 genes represent a structurally diverse set of transcripts and encode a functionally diverse set of proteins.

Conclusion

We have investigated the accuracy of the Twinscan and EuGene gene prediction programs and found them to be reliable predictors of gene structure in Arabidopsis. Several hundred previously un-annotated genes were validated by this work. Based upon this information derived from these efforts it is likely that the Arabidopsis genome annotation continues to overlook several hundred protein coding genes.

Comparative genomics: a tool to functionally annotate human DNA

The availability of an increasing number of vertebrate genomes has enabled comparative methods to infer functional sequences based on evolutionary constraint. Although this has proved powerful for gene identification, significant progress has also been made in uncovering gene regulatory sequences such as distant acting transcriptional enhancers. These pursuits have led to the development of a variety of valuable databases and resources that should serve as a routine toolbox for biological discovery.

Sequence analyses to study the evolutionary history and cis-regulatory elements of Hedgehog genes

Sequence analysis and comparative genomics are powerful tools to gain knowledge on multiple aspects of gene and protein regulation and function. These have been widely used to understand the evolutionary history and the biochemistry of Hedgehog (Hh) proteins, and the molecular control of Hedgehog gene expression. Here, we report on some of the methods available to retrieve protein and genomic sequences. We describe how protein sequence comparison can produce information on the evolutionary history of Hh proteins. Moreover, we describe the use of genomic sequence analysis including phylogenetic footprinting and transcription factor-binding site search tools, techniques that allow for the characterization of cis-regulatory elements of developmental genes such as the Hedgehog genes.

Btk expression is controlled by Oct and BOB.1/OBF.1

BOB.1/OBF.1 is a lymphocyte-restricted transcriptional coactivator. It binds together with the Oct1 and Oct2 transcription factors to DNA and enhances their transactivation potential. Mice deficient for the transcriptional coactivator BOB.1/OBF.1 show several defects in differentiation, function and signaling of B cells. In search of BOB.1/OBF.1 regulated genes we identified Btk—a cytoplasmic tyrosine kinase—as a direct target of BOB.1/OBF.1. Analyses of the human as well as murine Btk promoters revealed a non-consensus octamer site close to the start site of transcription. Here we show that Oct proteins together with BOB.1/OBF.1 are able to form ternary complexes on these sites in vitro and in vivo. This in turn leads to the induction of Btk promoter activity in synergism with the transcription factor PU.1. Btk, like BOB.1/OBF.1, plays a critical role in B cell development and B cell receptor signalling. Therefore the down-regulation of Btk expression in BOB.1/OBF.1-deficient B cells could be related to the functional and developmental defects observed in BOB.1/OBF.1-deficient mice.

The human desmin locus: gene organization and LCR-mediated transcriptional control

Locus control regions (LCRs) are defined by their ability to confer reproducible physiological levels of transgene expression in mice and therefore thought to possess the ability to generate dominantly a transcriptionally active chromatin structure. We report the first characterization of a muscle-cell-specific LCR, which is linked to the human desmin gene (DES). The DES LCR consists of five regions of muscle-specific DNase I hypersensitivity (HS) localized between -9 and -18 kb 5' of DES and reproducibly drives full physiological levels of expression in all muscle cell types. The DES LCR DNase I HS regions are highly conserved between humans and other mammals and can potentially bind a broad range of muscle-specific and ubiquitous transcription factors. Bioinformatics and direct molecular analysis show that the DES locus consists of three muscle-specific (DES) or muscle preferentially expressed genes (APEG1 and SPEG, the human orthologue of murine striated-muscle-specific serine/threonine protein kinase, Speg). The DES LCR may therefore regulate expression of SPEG and APEG1 as well as DES.

Condition-specific coregulation with cis-regulatory motifs and modules in the mouse genome

Deciphering genetic regulatory codes remains a challenge. Here, we present an effective approach to identifying in vivo condition-specific coregulation with cis-regulatory motifs and modules in the mouse genome. A resampling-based algorithm was adopted to cluster our microarray data of a stress response, which generated 35 tight clusters with unique expression patterns containing 811 genes of 5652 genes significantly altered. Database searches identified many known motifs within the 3-kb regulatory regions of 40 genes from 3 clusters and modules with six to nine motifs that were commonly shared by 60-100% of these genes. The upstream regulatory region contained the highest frequency of these common motifs. CisModule program predictions were comparable with the results from database searches and found four potentially novel motifs. This result indicates that these motifs and modules could be responsible for gene coregulation of the stress response in the lacrimal gland.

Trade-offs in detecting evolutionarily constrained sequence by comparative genomics

As whole-genome sequencing efforts extend beyond more traditional model organisms to include a deep diversity of species, comparative genomic analyses will be further empowered to reveal insights into the human genome and its evolution. The discovery and annotation of functional genomic elements is a necessary step toward a detailed understanding of our biology, and sequence comparisons have proven to be an integral tool for that task. This review is structured to broadly reflect the statistical challenges in discriminating these functional elements from the bulk of the genome that has evolved neutrally. Specifically, we review the comparative genomics literature in terms of specificity, sensitivity, and phylogenetic scope, as well as the trade-offs that relate these factors in standard analyses. We consider the impact of an expanding diversity of orthologous sequences on our ability to resolve functional elements. This impact is assessed through both recent comparative analyses of deep alignments and mathematical modeling.

Evaluation of the RET regulatory landscape reveals the biological relevance of a HSCR-implicated enhancer

Evolutionary sequence conservation is now a relatively common approach for the prediction of functional DNA sequences. However, the fraction of conserved non-coding sequences with regulatory potential is still unknown. In this study, we focus on elucidating the regulatory landscape of RET, a crucial developmental gene within which we have recently identified a regulatory Hirschsprung disease (HSCR) susceptibility variant. We report a systematic examination of conserved non-coding sequences (n=45) identified in a 220 kb interval encompassing RET. We demonstrate that most of these conserved elements are capable of enhancer or suppressor activity in vitro, and the majority of the elements exert cell type-dependent control. We show that discrete sequences within regulatory elements can bind nuclear protein in a cell type-dependent manner that is consistent with their identified in vitro regulatory control. Finally, we focused our attention on the enhancer implicated in HSCR to demonstrate that this element drives reporter expression in cell populations of the excretory system and central nervous system (CNS) and peripheral nervous system (PNS), consistent with expression of the endogenous RET protein. Importantly, this sequence also modulates expression in the enteric nervous system consistent with its proposed role in HSCR.

Highly conserved non-coding sequences are associated with vertebrate development

In addition to protein coding sequence, the human genome contains a significant amount of regulatory DNA, the identification of which is proving somewhat recalcitrant to both in silico and functional methods. An approach that has been used with some success is comparative sequence analysis, whereby equivalent genomic regions from different organisms are compared in order to identify both similarities and differences. In general, similarities in sequence between highly divergent organisms imply functional constraint. We have used a whole-genome comparison between humans and the pufferfish, Fugu rubripes, to identify nearly 1,400 highly conserved non-coding sequences. Given the evolutionary divergence between these species, it is likely that these sequences are found in, and furthermore are essential to, all vertebrates. Most, and possibly all, of these sequences are located in and around genes that act as developmental regulators. Some of these sequences are over 90% identical across more than 500 bases, being more highly conserved than coding sequence between these two species. Despite this, we cannot find any similar sequences in invertebrate genomes. In order to begin to functionally test this set of sequences, we have used a rapid in vivo assay system using zebrafish embryos that allows tissue-specific enhancer activity to be identified. Functional data is presented for highly conserved non-coding sequences associated with four unrelated developmental regulators (SOX21, PAX6, HLXB9, and SHH), in order to demonstrate the suitability of this screen to a wide range of genes and expression patterns. Of 25 sequence elements tested around these four genes, 23 show significant enhancer activity in one or more tissues. We have identified a set of non-coding sequences that are highly conserved throughout vertebrates. They are found in clusters across the human genome, principally around genes that are implicated in the regulation of development, including many transcription factors. These highly conserved non-coding sequences are likely to form part of the genomic circuitry that uniquely defines vertebrate development.

R11: a cis-regulatory node of the sea urchin embryo gene network that controls early expression of SpDelta in micromeres

A gene regulatory network (GRN) controls the process by which the endomesoderm of the sea urchin embryo is specified. In this GRN, the program of gene expression unique to the skeletogenic micromere lineage is set in train by activation of the pmar1 gene. Through a double repression system, this gene is responsible for localization of expression of downstream regulatory and signaling genes to cells of this lineage. One of these genes, delta, encodes a Notch ligand, and its expression in the right place and time is crucial to the specification of the endomesoderm. Here we report a cis-regulatory element R11 that is responsible for localizing the expression of delta by means of its response to the pmar1 repression system. R11 was identified as an evolutionarily conserved genomic sequence located about 13 kb downstream of the last exon of the delta gene. We demonstrate here that this cis-regulatory element is able to drive the expression of a reporter gene in the same cells and at the same time that the endogenous delta gene is expressed, and that temporally, spatially, and quantitatively it responds to the pmar1 repression system just as predicted for the delta gene in the endomesoderm GRN. This work illustrates the application of cis-regulatory analysis to the validation of predictions of the GRN model. In addition, we introduce new methodological tools for quantitative measurement of the output of expression constructs that promise to be of general value for cis-regulatory analysis in sea urchin embryos.

Microarray and comparative genomics-based identification of genes and gene regulatory regions of the mouse immune system

Background

In this study we have built and mined a gene expression database composed of 65 diverse mouse tissues for genes preferentially expressed in immune tissues and cell types. Using expression pattern criteria, we identified 360 genes with preferential expression in thymus, spleen, peripheral blood mononuclear cells, lymph nodes (unstimulated or stimulated), or in vitro activated T-cells.

Results

Gene clusters, formed based on similarity of expression-pattern across either all tissues or the immune tissues only, had highly significant associations both with immunological processes such as chemokine-mediated response, antigen processing, receptor-related signal transduction, and transcriptional regulation, and also with more general processes such as replication and cell cycle control. Within-cluster gene correlations implicated known associations of known genes, as well as immune process-related roles for poorly described genes. To characterize regulatory mechanisms and cis-elements of genes with similar patterns of expression, we used a new version of a comparative genomics-based cis-element analysis tool to identify clusters of cis-elements with compositional similarity among multiple genes. Several clusters contained genes that shared 5–6 cis-elements that included ETS and zinc-finger binding sites. cis-Elements AP2 EGRF ETSF MAZF SP1F ZF5F and AREB ETSF MZF1 PAX5 STAT were shared in a thymus-expressed set; AP4R E2FF EBOX ETSF MAZF SP1F ZF5F and CREB E2FF MAZF PCAT SP1F STAT cis-clusters occurred in activated T-cells; CEBP CREB NFKB SORY and GATA NKXH OCT1 RBIT occurred in stimulated lymph nodes.

Conclusion

This study demonstrates a series of analytic approaches that have allowed the implication of genes and regulatory elements that participate in the differentiation, maintenance, and function of the immune system. Polymorphism or mutation of these could adversely impact immune system functions.

Comparative genomics

The genomes from three mammals (human, mouse, and rat), two worms, and several yeasts have been sequenced, and more genomes will be completed in the near future for comparison with those of the major model organisms. Scientists have used various methods to align and compare the sequenced genomes to address critical issues in genome function and evolution. This review covers some of the major new insights about gene content, gene regulation, and the fraction of mammalian genomes that are under purifying selection and presumed functional. We review the evolutionary processes that shape genomes, with particular attention to variation in rates within genomes and along different lineages. Internet resources for accessing and analyzing the treasure trove of sequence alignments and annotations are reviewed, and we discuss critical problems to address in new bioinformatic developments in comparative genomics.

Methods in comparative genomics: genome correspondence, gene identification and regulatory motif discovery

In Kellis et al. (2003), we reported the genome sequences of S. paradoxus, S. mikatae, and S. bayanus and compared these three yeast species to their close relative, S. cerevisiae. Genomewide comparative analysis allowed the identification of functionally important sequences, both coding and noncoding. In this companion paper we describe the mathematical and algorithmic results underpinning the analysis of these genomes. (1) We present methods for the automatic determination of genome correspondence. The algorithms enabled the automatic identification of orthologs for more than 90% of genes and intergenic regions across the four species despite the large number of duplicated genes in the yeast genome. The remaining ambiguities in the gene correspondence revealed recent gene family expansions in regions of rapid genomic change. (2) We present methods for the identification of protein-coding genes based on their patterns of nucleotide conservation across related species. We observed the pressure to conserve the reading frame of functional proteins and developed a test for gene identification with high sensitivity and specificity. We used this test to revisit the genome of S. cerevisiae, reducing the overall gene count by 500 genes (10% of previously annotated genes) and refining the gene structure of hundreds of genes. (3) We present novel methods for the systematic de novo identification of regulatory motifs. The methods do not rely on previous knowledge of gene function and in that way differ from the current literature on computational motif discovery. Based on genomewide conservation patterns of known motifs, we developed three conservation criteria that we used to discover novel motifs. We used an enumeration approach to select strongly conserved motif cores, which we extended and collapsed into a small number of candidate regulatory motifs. These include most previously known regulatory motifs as well as several noteworthy novel motifs. The majority of discovered motifs are enriched in functionally related genes, allowing us to infer a candidate function for novel motifs. Our results demonstrate the power of comparative genomics to further our understanding of any species. Our methods are validated by the extensive experimental knowledge in yeast and will be invaluable in the study of complex genomes like that of the human.

CONFAC: automated application of comparative genomic promoter analysis to DNA microarray datasets

The advent of DNA microarray technology and the sequencing of multiple vertebrate genomes has provided a unique opportunity for the integration of comparative genomics with high-throughput gene expression analysis. Here we describe the conserved transcription factor binding site (CONFAC) software that enables the high-throughput identification of conserved transcription factor binding sites (TFBSs) in the regulatory regions of hundreds of genes at a time (http://morenolab.whitehead.emory.edu/cgi-bin/confac/login.pl). The CONFAC software compares non-coding regulatory sequences between human and mouse genomes to enable identification of conserved TFBSs that are significantly enriched in promoters of gene clusters from microarray analyses compared to sets of unchanging control genes using a Mann-Whitney U-test. Analysis of random gene sets demonstrated that using our approach, over 98% of TFBSs had false positive rates below 5%. As a proof-of-principle, we have validated the CONFAC software using gene sets from four separate microarray studies and identified TFBSs known to be functionally important for regulation of each of the four gene sets.

Comparative sequence and x-inactivation analyses of a domain of escape in human xp11.2 and the conserved segment in mouse

We have performed X-inactivation and sequence analyses on 350 kb of sequence from human Xp11.2, a region shown previously to contain a cluster of genes that escape X inactivation, and we compared this region with the region of conserved synteny in mouse. We identified several new transcripts from this region in human and in mouse, which defined the full extent of the domain escaping X inactivation in both species. In human, escape from X inactivation involves an uninterrupted 235-kb domain of multiple genes. Despite highly conserved gene content and order between the two species, Smcx is the only mouse gene from the conserved segment that escapes inactivation. As repetitive sequences are believed to facilitate spreading of X inactivation along the chromosome, we compared the repetitive sequence composition of this region between the two species. We found that long terminal repeats (LTRs) were decreased in the human domain of escape, but not in the majority of the conserved mouse region adjacent to Smcx in which genes were subject to X inactivation, suggesting that these repeats might be excluded from escape domains to prevent spreading of silencing. Our findings indicate that genomic context, as well as gene-specific regulatory elements, interact to determine expression of a gene from the inactive X-chromosome.

Comparative analysis of the alpha-like globin clusters in mouse, rat, and human chromosomes indicates a mechanism underlying breaks in conserved synteny

We have sequenced and fully annotated a 65,871-bp region of mouse Chromosome 17 including the Hba-ps4 alpha-globin pseudogene. Comparative sequence analysis with the functional alpha-globin loci at human Chromosome 16p13.3 and mouse Chromosome 11 shows that this segment of mouse Chromosome 17 contains a group of three alpha-like pseudogenes (Hba-psm-Hba-ps4-Hba-q3), similar to the duplicated sets found at the functional mouse cluster on Chromosome 11. In addition, exons 7 to 12 of the mLuc7L gene are present just downstream from the pseudogene cluster, indicating that this clone contains the region in which human 16p13.3 switches in synteny between mouse Chromosomes 11 and 17. Comparison of the sequences around the alpha-like clusters on the two mouse chromosomes reveals the presence of conserved tandem repeats. We propose that these repetitive elements have played a role in the fragmentation of the mouse alpha cluster during evolution.

zPicture: dynamic alignment and visualization tool for analyzing conservation profiles

Comparative sequence analysis has evolved as an essential technique for identifying functional coding and noncoding elements conserved throughout evolution. Here, we introduce zPicture, an interactive Web-based sequence alignment and visualization tool for dynamically generating conservation profiles and identifying evolutionarily conserved regions (ECRs). zPicture is highly flexible, because critical parameters can be modified interactively, allowing users to differentially predict ECRs in comparisons of sequences of different phylogenetic distances and evolutionary rates. We demonstrate the application of this module to identify a known regulatory element in the HOXD locus, in which functional ECRs are difficult to discern against the highly conserved genomic background. zPicture also facilitates transcription factor binding-site analysis via the rVista tool portal. We present an example of the HBB complex when zPicture/rVista combination specifically pinpoints to two ECRs containing GATA-1, NF-E2, and TAL1/E47 binding sites that were identified previously as transcriptional enhancers. In addition, zPicture is linked to the UCSC Genome Browser, allowing users to automatically extract sequences and gene annotations for any recorded locus. Finally, we describe how this tool can be efficiently applied to the analysis of nonvertebrate genomes, including those of microbial organisms.

A Distal Enhancer in the Interferon-γ (IFN-γ) Locus Revealed by Genome Sequence Comparison

Large-scale cross-species DNA sequence comparison has become a powerful tool to identify conserved cis-regulatory modules of genes. However, bioinformatic analysis alone cannot reveal how an evolutionarily conserved region regulates gene expression: whether it functions as an enhancer, silencer, or insulator; whether its function is cell-type restricted; and whether biologically relevant transcription factors bind to the element. Here we combine bioinformatics with wet-lab techniques to illustrate a general and systematic method of identifying functional conserved regulatory regions of genes. We applied this approach to the interferon-gamma (IFN-gamma) gene. Comparison of human and mouse IFN-gamma reveals a highly conserved non-coding sequence located approximately 5 kb 5' of the transcription start site. This region coincides with constitutive and inducible DNase I hypersensitivity sites present in IFN-gamma-producing Th1 cells but not in Th2 cells that do not produce IFN-gamma. Histone methylation at the 5' conserved non-coding sequences indicates a more accessible chromatin structure in Th1 cells compared with Th2 cells. This element binds two transcription factors known to be essential for IFN-gamma expression: nuclear factor of activated T cells, an inducible transcription factor, and T-box protein expressed in T cells, a cell lineage-restricted transcription factor. Together, these findings identify a highly conserved distal enhancer in the IFN-gamma cytokine locus and validate our approach as a successful method to detect cis-regulatory elements.

Analysis of the gene-dense major histocompatibility complex class III region and its comparison to mouse

In mammals, the Major Histocompatibility Complex class I and II gene clusters are separated by an approximately 700-kb stretch of sequence called the MHC class III region, which has been associated with susceptibility to numerous diseases. To facilitate understanding of this medically important and architecturally interesting portion of the genome, we have sequenced and analyzed both the human and mouse class III regions. The cross-species comparison has facilitated the identification of 60 genes in human and 61 in mouse, including a potential RNA gene for which the introns are more conserved across species than the exons. Delineation of global organization, gene structure, alternative splice forms, protein similarities, and potential cis-regulatory elements leads to several conclusions: (1) The human MHC class III region is the most gene-dense region of the human genome: >14% of the sequence is coding, approximately 72% of the region is transcribed, and there is an average of 8.5 genes per 100 kb. (2) Gene sizes, number of exons, and intergenic distances are for the most part similar in both species, implying that interspersed repeats have had little impact in disrupting the tight organization of this densely packed set of genes. (3) The region contains a heterogeneous mixture of genes, only a few of which have a clearly defined and proven function. Although many of the genes are of ancient origin, some appear to exist only in mammals and fish, implying they might be specific to vertebrates. (4) Conserved noncoding sequences are found primarily in or near the 5'-UTR or the first intron of genes, and seldom in the intergenic regions. Many of these conserved blocks are likely to be cis-regulatory elements.

A combinatorial network of evolutionarily conserved myelin basic protein regulatory sequences confers distinct glial-specific phenotypes

Myelin basic protein (MBP) is required for normal myelin compaction and is implicated in both experimental and human demyelinating diseases. In this study, as an initial step in defining the regulatory network controlling MBP transcription, we located and characterized the function of evolutionarily conserved regulatory sequences. Long-range human-mouse sequence comparison revealed over 1 kb of conserved noncoding MBP 5' flanking sequence distributed into four widely spaced modules ranging from 0.1 to 0.4 kb. We demonstrate first that a controlled strategy of transgenesis provides an effective means to assign and compare qualitative and quantitative in vivo regulatory programs. Using this strategy, single-copy reporter constructs, designed to evaluate the regulatory significance of modular and intermodular sequences, were introduced by homologous recombination into the mouse hprt (hypoxanthine-guanine phosphoribosyltransferase) locus. The proximal modules M1 and M2 confer comparatively low-level oligodendrocyte expression primarily limited to early postnatal development, whereas the upstream M3 confers high-level oligodendrocyte expression extending throughout maturity. Furthermore, constructs devoid of M3 fail to target expression to newly myelinating oligodendrocytes in the mature CNS. Mutation of putative Nkx6.2/Gtx sites within M3, although not eliminating oligodendrocyte targeting, significantly decreases transgene expression levels. High-level and continuous expression is conferred to myelinating or remyelinating Schwann cells by M4. In addition, when isolated from surrounding MBP sequences, M3 confers transient expression to Schwann cells elaborating myelin. These observations define the in vivo regulatory roles played by conserved noncoding MBP sequences and lead to a combinatorial model in which different regulatory modules are engaged during primary myelination, myelin maintenance, and remyelination.

Fragile site orthologs FHIT/FRA3B and Fhit/Fra14A2: evolutionarily conserved but highly recombinogenic

Common fragile sites are regions that show elevated susceptibility to DNA damage, leading to alterations that can contribute to cancer development. FRA3B, located at chromosome region 3p14.2, is the most frequently expressed human common fragile site, and allelic losses at FRA3B have been observed in many types of cancer. The FHIT gene, encompassing the FRA3B region, is a tumor-suppressor gene. To identify the features of FHIT/FRA3B that might contribute to fragility, sequences of the human FHIT and the flanking PTPRG gene were compared with those of murine Fhit and Ptprg. Human and mouse orthologous genes, FHIT and Fhit, are more highly conserved through evolution than PTPRG/Ptprg and yet contain more sequence elements that are exquisitely sensitive to genomic rearrangements, such as high-flexibility regions and long interspersed nuclear element 1s, suggesting that common fragile sites serve a function. The conserved AT-rich high-flexibility regions are the most characteristic of common fragile sites.

Recent advances in computational promoter analysis in understanding the transcriptional regulatory network

The computational approach to the study of transcriptional regulation networks has become more attractive and feasible with the rapid accumulation of complete genome sequences and the advance of high-throughput expression profiling technology. In this review, current computational approaches for understanding the transcriptional regulatory network, including promoter prediction, transcription factor binding site identification, combinatorial regulatory elements prediction, and transcription factor target gene identification, are discussed. The role of comparative genomics in transcription regulatory region analysis is also reviewed.

Sequencing and comparison of yeast species to identify genes and regulatory elements

Identifying the functional elements encoded in a genome is one of the principal challenges in modern biology. Comparative genomics should offer a powerful, general approach. Here, we present a comparative analysis of the yeast Saccharomyces cerevisiae based on high-quality draft sequences of three related species (S. paradoxus, S. mikatae and S. bayanus). We first aligned the genomes and characterized their evolution, defining the regions and mechanisms of change. We then developed methods for direct identification of genes and regulatory motifs. The gene analysis yielded a major revision to the yeast gene catalogue, affecting approximately 15% of all genes and reducing the total count by about 500 genes. The motif analysis automatically identified 72 genome-wide elements, including most known regulatory motifs and numerous new motifs. We inferred a putative function for most of these motifs, and provided insights into their combinatorial interactions. The results have implications for genome analysis of diverse organisms, including the human.

Comparative promoter analysis and its application in analysis of PTH-regulated gene expression

Taking advantage of the "working draft" of the human genome and the MIT shotgun assembly of the mouse genome, we performed a comparative promoter analysis of human RefSeq mRNA (sequences from GenBank's RefSeq database). By combining this analysis with a transcription factor (TF) binding site analysis using a TRANSFAC position weight matrix (PWM) search, 86% of non-specific TF sites were removed. Using a set of genes that are regulated by parathyroid hormone (PTH), a statistical analysis was performed on the conserved TF binding sites among a set of eight human and mouse genes. From among the eight genes tested, we obtained a set of 31 TFs, suggesting possible roles for associated genes in PTH-mediated pathways. All three known PTH-responsive TFs (AP1, RUNX2, CREB) were correctly predicted by this analysis as well as two other potential TFs (VDR and CEBP Delta). Additionally, a model was made to describe the TF site characteristic module of PTH-regulated genes. This model was then used to search all human RefSeq gene promoters with established human-mouse ortholog relationships to identify other PTH-regulated genes. This comparative approach combined with statistical analysis proved to be sufficiently specific to decipher critical TFs involved in PTH-regulated pathways.

Of mice and men: phylogenetic footprinting aids the discovery of regulatory elements

Phylogenetic footprinting is an approach to finding functionally important sequences in the genome that relies on detecting their high degrees of conservation across different species. A new study shows how much it improves the prediction of gene-regulatory elements in the human genome.

Covariation in frequencies of substitution, deletion, transposition, and recombination during eutherian evolution

Six measures of evolutionary change in the human genome were studied, three derived from the aligned human and mouse genomes in conjunction with the Mouse Genome Sequencing Consortium, consisting of (1) nucleotide substitution per fourfold degenerate site in coding regions, (2) nucleotide substitution per site in relics of transposable elements active only before the human-mouse speciation, and (3) the nonaligning fraction of human DNA that is nonrepetitive or in ancestral repeats; and three derived from human genome data alone, consisting of (4) SNP density, (5) frequency of insertion of transposable elements, and (6) rate of recombination. Features 1 and 2 are measures of nucleotide substitutions at two classes of "neutral" sites, whereas 4 is a measure of recent mutations. Feature 3 is a measure dominated by deletions in mouse, whereas 5 represents insertions in human. It was found that all six vary significantly in megabase-sized regions genome-wide, and many vary together. This indicates that some regions of a genome change slowly by all processes that alter DNA, and others change faster. Regional variation in all processes is correlated with, but not completely accounted for, by GC content in human and the difference between GC content in human and mouse.

Strategies and tools for whole-genome alignments

The availability of the assembled mouse genome makes possible, for the first time, an alignment and comparison of two large vertebrate genomes. We investigated different strategies of alignment for the subsequent analysis of conservation of genomes that are effective for assemblies of different quality. These strategies were applied to the comparison of the working draft of the human genome with the Mouse Genome Sequencing Consortium assembly, as well as other intermediate mouse assemblies. Our methods are fast and the resulting alignments exhibit a high degree of sensitivity, covering more than 90% of known coding exons in the human genome. We obtained such coverage while preserving specificity. With a view towards the end user, we developed a suite of tools and Web sites for automatically aligning and subsequently browsing and working with whole-genome comparisons. We describe the use of these tools to identify conserved non-coding regions between the human and mouse genomes, some of which have not been identified by other methods.

Leveraging the mouse genome for gene prediction in human: from whole-genome shotgun reads to a global synteny map

The availability of draft sequences for both the mouse and human genomes makes it possible, for the first time, to annotate whole mammalian genomes using comparative methods. TWINSCAN is a gene-prediction system that combines the methods of single-genome predictors like GENSCAN with information derived from genome comparison, thereby improving accuracy. Because TWINSCAN uses genomic sequence only, it is less biased toward highly and/or ubiquitously expressed genes than GENEWISE, GENOMESCAN, and other methods based on evidence derived from transcripts. We show that TWINSCAN improves gene prediction in human using intermediate products from various stages of the sequencing and analysis of the mouse genome, from low-redundancy, whole-genome shotgun reads to the draft assembly and the synteny map. TWINSCAN improves on the prior state of the art even when alignments from only 1X coverage of the mouse genome are available. Gene prediction accuracy improves steadily from 1X through 3X, more slowly from 3X to 4X, and relatively little thereafter. The assembly and the synteny map greatly speed the computations, however. Our human annotation using the mouse assembly is conservative, predicting only 25,622 genes, and appears to be one of the best de novo annotations of the human genome to date.

Initial sequencing and comparative analysis of the mouse genome

The sequence of the mouse genome is a key informational tool for understanding the contents of the human genome and a key experimental tool for biomedical research. Here, we report the results of an international collaboration to produce a high-quality draft sequence of the mouse genome. We also present an initial comparative analysis of the mouse and human genomes, describing some of the insights that can be gleaned from the two sequences. We discuss topics including the analysis of the evolutionary forces shaping the size, structure and sequence of the genomes; the conservation of large-scale synteny across most of the genomes; the much lower extent of sequence orthology covering less than half of the genomes; the proportions of the genomes under selection; the number of protein-coding genes; the expansion of gene families related to reproduction and immunity; the evolution of proteins; and the identification of intraspecies polymorphism.

Isolation and developmental expression analysis of Tbx22, the mouse homolog of the human X-linked cleft palate gene

Mutations in the TBX22 gene have been identified recently in patients with the X-linked cleft palate and ankyloglossia syndrome, suggesting that the TBX22 transcription factor plays an important role in palate development. However, because ankyloglossia has been reported in the majority of patients with TBX22 mutations, it has been speculated that the cleft palate phenotype is secondary to defective fetal tongue movement. To understand the role of TBX22 in disease pathogenesis and in normal development, it is necessary to carry out a detailed temporal and spatial gene expression analysis. We report here the isolation and developmental expression analysis of the mouse homolog Tbx22. The mouse Tbx22 gene encodes a putative protein of 517 amino acid residues, which shares 72% overall amino acid sequence identity with the human TBX22 protein. By using interspecific backcross analysis, we have localized the Tbx22 gene to mouse chromosome X, in a region syntenic to human chromosome Xq21, where the TBX22 gene resides, indicating that Tbx22 is the ortholog of human TBX22. Our in situ hybridization analysis shows that Tbx22 is expressed in a temporally and spatially highly restricted pattern during mouse palate and tongue development. Together with the mutant phenotypes in human patients, our data indicate a primary role for Tbx22 in both palate and tongue development.

Detection and visualization of compositionally similar cis-regulatory element clusters in orthologous and coordinately controlled genes

Evolutionarily conserved noncoding genomic sequences represent a potentially rich source for the discovery of gene regulatory regions. However, detecting and visualizing compositionally similar cis-element clusters in the context of conserved sequences is challenging. We have explored potential solutions and developed an algorithm and visualization method that combines the results of conserved sequence analyses (BLASTZ) with those of transcription factor binding site analyses (MatInspector) (http://trafac.chmcc.org). We define hits as the density of co-occurring cis-element transcription factor (TF)-binding sites measured within a 200-bp moving average window through phylogenetically conserved regions. The results are depicted as a Regulogram, in which the hit count is plotted as a function of position within each of the two genomic regions of the aligned orthologs. Within a high-scoring region, the relative arrangement of shared cis-elements within compositionally similar TF-binding site clusters is depicted in a Trafacgram. On the basis of analyses of several training data sets, the approach also allows for the detection of similarities in composition and relative arrangement of cis-element clusters within nonorthologous genes, promoters, and enhancers that exhibit coordinate regulatory properties. Known functional regulatory regions of nonorthologous and less-conserved orthologous genes frequently showed cis-element shuffling, demonstrating that compositional similarity can be more sensitive than sequence similarity. These results show that combining sequence similarity with cis-element compositional similarity provides a powerful aid for the identification of potential control regions.

A physical map of the mouse genome

A physical map of a genome is an essential guide for navigation, allowing the location of any gene or other landmark in the chromosomal DNA. We have constructed a physical map of the mouse genome that contains 296 contigs of overlapping bacterial clones and 16,992 unique markers. The mouse contigs were aligned to the human genome sequence on the basis of 51,486 homology matches, thus enabling use of the conserved synteny (correspondence between chromosome blocks) of the two genomes to accelerate construction of the mouse map. The map provides a framework for assembly of whole-genome shotgun sequence data, and a tile path of clones for generation of the reference sequence. Definition of the human-mouse alignment at this level of resolution enables identification of a mouse clone that corresponds to almost any position in the human genome. The human sequence may be used to facilitate construction of other mammalian genome maps using the same strategy.

The common fragile site FRA16D and its associated gene WWOX are highly conserved in the mouse at Fra8E1

Recently, several common fragile sites (CFSs) have been cloned and characterized, including the two most frequently observed in the human population, FRA3B and FRA16D. In addition to their high frequency of breakage, FRA3B and FRA16D colocalize with genes crossing large regions of breakage. At FRA3B, the fragile histidine triad (FHIT) gene spans more than 1 Mb, and at FRA16D, the WWOX gene spans more than 750 kb. It has also been shown that in Mus musculus, a CFS Fra14A2 and the mouse Fhit gene are conserved in the orthologous region of the genome. In this study, we positioned the ortholog to WWOX (Wox1) at chromosome band 8E1 in the mouse genome. To determine whether, like Fra14A2 and Fhit, Fra8E1 and Wox1 colocalized in the mouse, we prepared bacterial and yeast artificial chromosome probes, and we hybridized them to aphidicolin-treated mouse metaphase chromosomes. Our data demonstrate that Wox1 colocalizes with Fra8E1. Furthermore, the sequence from this region, including introns, is highly conserved over at least a 100-kb region. This evolutionary conservation suggests that the two most active CFSs share many features, and that CFSs and their associated genes may be necessary for cell survival.

Search for enhancers: teleost models in comparative genomic and transgenic analysis of cis regulatory elements

Homology searches between DNA sequences of evolutionary distant species (phylogenetic footprinting) offer a fast detection method for regulatory sequences. Because of the small size of their genomes, tetraodontid species such as the Japanese pufferfish and green spotted pufferfish have become attractive models for comparative genomics. A disadvantage of the tetraodontid species is, however, that they cannot be bred and manipulated routinely under laboratory conditions, so these species are less attractive for developmental and genetic analysis. In contrast, an increasing arsenal of transgene techniques with the developmental model species zebrafish and medaka are being used for functional analysis of cis regulatory sequences. The main disadvantage is the much larger genome. While comparison between many loci proved the suitability of phylogenetic footprinting using fish and mammalian sequences, fast rate of change in enhancer structure and gene duplication within teleosts may obscure detection of homologies. Here we discuss the contribution and potentials provided by different teleost models for the detection and functional analysis of conserved cis-regulatory elements.