Dcode.org anthology of comparative genomic tools

From Beginning to End: Expanding the SERINC3 Interactome Through an in silico Analysis

The serine incorporator (SERINC) family of proteins are a family of multipass transmembrane proteins associated with biosynthesis of serine-containing phospholipids and sphingolipids. Humans have 5 paralogs, SERINC1-5, which have been linked to disease including variable expression in tumor lines and possessing activity as restriction factors against HIV-1. Despite recent studies, the cellular function of SERINC proteins have yet to be fully elucidated. The goal of this study as to investigate the role of SERINC3 by expanding upon its interactome. We used a variety of bioinformatic tools to identify cellular factors that interact with SERINC3 and assessed how sequence variation might alter these interactions. Analysis of the promoter region indicates that SERINC3 is putatively regulated by transcription factors involved in tissue-specific development. Analysis of the unique 3′-untranslated region of one variant of HsSERINC3 revealed that this region serves as a conserved site of regulation by both RNA binding proteins and miRNA. In addition, SERINC3 is putatively regulated at the protein level by several posttranslational modifications. Our results show that extra-membrane portions of SERINC3 are subject to variation in the coding sequence as well as areas of relatively low conservation. Overall, our data suggest that regions of low homology as well as presence of variations in the nucleotide and protein sequences of HsSERINC3 suggest that these variations may lead to aberrant function and alternative regulatory mechanisms in homologs. The functional consequences of these sequence and structural variations need to be explored systematically to fully appreciate the role of SERINC3 in both health and disease.

New Findings on LMO7 Transcripts, Proteins and Regulatory Regions in Human and Vertebrate Model Organisms and the Intracellular Distribution in Skeletal Muscle Cells

LMO7 is a multifunctional PDZ–LIM protein that can interact with different molecular partners and is found in several intracellular locations. The aim of this work was to shed light on LMO7 evolution, alternative transcripts, protein structure and gene regulation through multiple in silico analyses. We also explored the intracellular distribution of the LMO7 protein in chicken and zebrafish embryonic skeletal muscle cells by means of confocal fluorescence microscopy. Our results revealed a single LMO7 gene in mammals, sauropsids, Xenopus and in the holostean fish spotted gar while two lmo7 genes (lmo7a and lmo7b) were identified in teleost fishes. In addition, several different transcripts were predicted for LMO7 in human and in major vertebrate model organisms (mouse, chicken, Xenopus and zebrafish). Bioinformatics tools revealed several structural features of the LMO7 protein including intrinsically disordered regions. We found the LMO7 protein in multiple intracellular compartments in chicken and zebrafish skeletal muscle cells, such as membrane adhesion sites and the perinuclear region. Curiously, the LMO7 protein was detected within the nuclei of muscle cells in chicken but not in zebrafish. Our data showed that a conserved regulatory element may be related to muscle-specific LMO7 expression. Our findings uncover new and important information about LMO7 and open new challenges to understanding how the diverse regulation, structure and distribution of this protein are integrated into highly complex vertebrate cellular milieux, such as skeletal muscle cells.

Effect of rehabilitation on biologic and transcriptomic responses after hospital-acquired deconditioning: a prospective longitudinal feasibility study

PURPOSE

The objective of this study is to explore the transcriptomic and biologic variables characterizing the longitudinal rehabilitation intervention of patients with hospital-acquired deconditioning (HAD).

METHODS

This prospective clinical trial recruited HAD patients (n = 10) who spent ≥3 weeks hospitalized and then received inpatient rehabilitation. Functional improvement was measured using the Functional Independence Measure (FIM). Transcriptomic and biological variables were recorded at rehabilitation admission and 1, 2, 4, and 6 weeks post-admission. RNA sequencing studied the temporal changes of gene expression in leukocytes. Between-subject transcriptome comparisons were performed using principle component analysis. Within-subject changes in gene expression were analyzed using a gene ontology hierarchical clustering to identify common biological terms. Heart rate, weight, albumin, creatinine, and complete blood counts were analyzed.

RESULTS

Patients average age was 50.6 ± 7.2, FIM increased during inpatient rehabilitation (p = 0.01), weight increased (p = 0.01), lymphocytes decreased (p = 0.05), neutrophil increased (0.03) and red cell distribution width decreased (p = 0.05). The temporal profiles of gene expression revealed within-patient homogeneity and between-patients heterogeneity. The biological terms "bone morphogenesis" and "muscle cell development" were the most significantly enriched differentially expressed genes.

CONCLUSION

Transcriptomic and biologic markers paralleled the functional improvements of HAD patients during inpatient rehabilitation. Transcriptomic analyses were consistent with the cohort heterogeneity. Enrichment of the biological pathways bone morphogenesis and muscle cell development constituted evidence at the gene expression level of the effect of rehabilitation. Larger studies of various rehabilitation patient groups may increase gene expression profile homogeneity. Objective transcriptomic and biologic markers have the potential to improve the rehabilitation of HAD patients.IMPLICATIONS FOR REHABILITATIONNovel gene expression methods are increasingly being integrated into clinical practice and may apply to rehabilitation.Patients with hospital-acquired deconditioning (HAD) enriched gene expression of pathways targeted by inpatient rehabilitation such as bone morphogenesis and muscle cell development.The gene expression paralleled functional improvement of HAD patients.These data demonstrated the feasibility of molecular methods to identify markers of rehabilitation success in HAD patients.

Sex chromosomes drive gene expression and regulatory dimorphisms in mouse embryonic stem cells

BACKGROUND

Pre-implantation embryos exhibit sexual dimorphisms in both primates and rodents. To determine whether these differences reflected sex-biased expression patterns, we generated transcriptome profiles for six 40,XX, six 40,XY, and two 39,X mouse embryonic stem (ES) cells by RNA sequencing.

RESULTS

We found hundreds of coding and non-coding RNAs that were differentially expressed between male and female cells. Surprisingly, the majority of these were autosomal and included RNA encoding transcription and epigenetic and chromatin remodeling factors. We showed differential Prdm14-responsive enhancer activity in male and female cells, correlating with the sex-specific levels of Prdm14 expression. This is the first time sex-specific enhancer activity in ES cells has been reported. Evaluation of X-linked gene expression patterns between our XX and XY lines revealed four distinct categories: (1) genes showing 2-fold greater expression in the female cells; (2) a set of genes with expression levels well above 2-fold in female cells; (3) genes with equivalent RNA levels in male and female cells; and strikingly, (4) a small number of genes with higher expression in the XY lines. Further evaluation of autosomal gene expression revealed differential expression of imprinted loci, despite appropriate parent-of-origin patterns. The 39,X lines aligned closely with the XY cells and provided insights into potential regulation of genes associated with Turner syndrome in humans. Moreover, inclusion of the 39,X lines permitted three-way comparisons, delineating X and Y chromosome-dependent patterns.

CONCLUSIONS

Overall, our results support the role of the sex chromosomes in establishing sex-specific networks early in embryonic development and provide insights into effects of sex chromosome aneuploidies originating at those stages.

The AIRE -230Y Polymorphism Affects AIRE Transcriptional Activity: Potential Influence on AIRE Function in the Thymus

BACKGROUND

The autoimmune regulator (AIRE) is expressed in the thymus, particularly in thymic medullary epithelial cells (mTECs), and is required for the ectopic expression of a diverse range of peripheral tissue antigens by mTECs, facilitating their ability to perform negative selection of auto-reactive immature T-cells. The expression profile of peripheral tissue antigens is affected not only by AIRE deficiency but also with variation of AIRE activity in the thymus.

METHOD AND RESULTS

Therefore we screened 591bp upstream of the AIRE transcription start site including AIRE minimal promoter for single nucleotide polymorphism (SNPs) and identified two SNPs -655R (rs117557896) and -230Y (rs751032) respectively. To study the effect of these variations on AIRE promoter activity we generated a Flp-In host cell line which was stably transfected with a single copy of the reporter vector. Relative promoter activity was estimated by comparing the luciferase specific activity for lysates of the different reporter AIRE promoter-reporter gene constructs including AIRE-655G AIRE-230C, AIRE-655G AIRE-230T and AIRE-655A AIRE-230C. The analysis showed that the commonest haplotype AIRE-655G AIRE-230C has the highest luciferase specific activity (p<0.001). Whereas AIRE-655G AIRE-230T has a luciferase specific activity value that approaches null. Both AIRE promoter polymorphic sites have one allele that forms a CpG methylation site which we determined can be methylated in methylation assays using the M.SssI CpG methyltransferase.

CONCLUSION

AIRE-230Y is in a conserved region of the promoter and is adjacent to a predicted WT1 transcription factor binding site, suggesting that AIRE-230Y affects AIRE expression by influencing the binding of biochemical factors to this region. Our findings show that AIRE-655GAIRE-230T haplotype could dramatically alter AIRE transcription and so have an effect on the process of negative selection and affect susceptibility to autoimmune conditions.

Identification and computational analysis of gene regulatory elements

Over the last two decades, advances in experimental and computational technologies have greatly facilitated genomic research. Next-generation sequencing technologies have made de novo sequencing of large genomes affordable, and powerful computational approaches have enabled accurate annotations of genomic DNA sequences. Charting functional regions in genomes must account for not only the coding sequences, but also noncoding RNAs, repetitive elements, chromatin states, epigenetic modifications, and gene regulatory elements. A mix of comparative genomics, high-throughput biological experiments, and machine learning approaches has played a major role in this truly global effort. Here we describe some of these approaches and provide an account of our current understanding of the complex landscape of the human genome. We also present overviews of different publicly available, large-scale experimental data sets and computational tools, which we hope will prove beneficial for researchers working with large and complex genomes.

GenomeVISTA--an integrated software package for whole-genome alignment and visualization

With the ubiquitous generation of complete genome assemblies for a variety of species, efficient tools for whole-genome alignment along with user-friendly visualization are critically important. Our VISTA family of tools for comparative genomics, based on algorithms for pairwise and multiple alignments of genomic sequences and whole-genome assemblies, has become one of the standard techniques for comparative analysis. Most of the VISTA programs have been implemented as Web-accessible servers and are extensively used by the biomedical community. In this manuscript, we introduce GenomeVISTA: a novel implementation that incorporates most features of the VISTA family--fast and accurate alignment, visualization capabilities, GUI and analytical tools within a stand-alone software package. GenomeVISTA thus provides flexibility and security for users who need to conduct whole-genome comparisons on their own computers.

AVAILABILITY AND IMPLEMENTATION

Implemented in Perl, C/C++ and Java, the source code is freely available for download at the VISTA Web site: http://genome.lbl.gov/vista/.

Identification of candidate transcription factor binding sites in the cattle genome

A resource that provides candidate transcription factor binding sites (TFBSs) does not currently exist for cattle. Such data is necessary, as predicted sites may serve as excellent starting locations for future omics studies to develop transcriptional regulation hypotheses. In order to generate this resource, we employed a phylogenetic footprinting approach—using sequence conservation across cattle, human and dog—and position-specific scoring matrices to identify 379,333 putative TFBSs upstream of nearly 8000 Mammalian Gene Collection (MGC) annotated genes within the cattle genome. Comparisons of our predictions to known binding site loci within the PCK1, ACTA1 and G6PC promoter regions revealed 75% sensitivity for our method of discovery. Additionally, we intersected our predictions with known cattle SNP variants in dbSNP and on the Illumina BovineHD 770k and Bos 1 SNP chips, finding 7534, 444 and 346 overlaps, respectively. Due to our stringent filtering criteria, these results represent high quality predictions of putative TFBSs within the cattle genome. All binding site predictions are freely available at http://bfgl.anri.barc.usda.gov/BovineTFBS/ or http://199.133.54.77/BovineTFBS.

A novel MSMB-related microprotein in the postovulatory egg coats of marsupials

BACKGROUND

Early marsupial conceptuses differ markedly from those of eutherian mammals, especially during cleavage and early blastocyst stages of development. Additionally, in marsupials the zona pellucida is surrounded by two acellular layers, the mucoid coat and shell, which are formed from secretions from the reproductive tract.

RESULTS

We report the identification of a novel postovulatory coat component in marsupials, which we call uterinesecreted microprotein (USM). USM belongs to a family of disulfide-rich microproteins of unconfirmed function that is found throughout deuterostomes and in some protostomes, and includes β-microseminoprotein (MSMB) and prostate-associated microseminoprotein (MSMP). We describe the evolution of this family in detail, including USM-related sequences in other vertebrates. The orthologue of USM in the tammar wallaby, USM1, is expressed by the endometrium with a dynamic temporal profile, possibly under the control of progesterone.

CONCLUSIONS

USM appears to have evolved in a mammalian ancestor specifically as a component of the postovulatory coats. By analogy with the known properties of MSMB, it may have roles in regulating sperm motility/survival or in the immune system. However, its C-terminal domain is greatly truncated compared with MSMB, suggesting a divergent function.

Comparative analysis of the ATRX promoter and 5' regulatory region reveals conserved regulatory elements which are linked to roles in neurodevelopment, alpha-globin regulation and testicular function

Background

ATRX is a tightly-regulated multifunctional protein with crucial roles in mammalian development. Mutations in the ATRX gene cause ATR-X syndrome, an X-linked recessive developmental disorder resulting in severe mental retardation and mild alpha-thalassemia with facial, skeletal and genital abnormalities. Although ubiquitously expressed the clinical features of the syndrome indicate that ATRX is not likely to be a global regulator of gene expression but involved in regulating specific target genes. The regulation of ATRX expression is not well understood and this is reflected by the current lack of identified upstream regulators. The availability of genomic data from a range of species and the very highly conserved 5' regulatory regions of the ATRX gene has allowed us to investigate putative transcription factor binding sites (TFBSs) in evolutionarily conserved regions of the mammalian ATRX promoter.

Results

We identified 12 highly conserved TFBSs of key gene regulators involved in biologically relevant processes such as neural and testis development and alpha-globin regulation.

Conclusions

Our results reveal potentially important regulatory elements in the ATRX gene which may lead to the identification of upstream regulators of ATRX and aid in the understanding of the molecular mechanisms that underlie ATR-X syndrome.

Conserved proximal promoter elements control repulsive guidance molecule c/hemojuvelin (Hfe2) gene transcription in skeletal muscle

Repulsive guidance molecule c (RGMc; gene symbol: Hfe2) plays a critical role in iron metabolism. Inactivating mutations cause juvenile hemochromatosis, a severe iron overload disorder. Understanding mechanisms controlling RGMc biosynthesis has been hampered by minimal information about the RGMc gene. Here we define the structure, examine the evolution, and establish mechanisms of regulation of the mouse RGMc gene. RGMc is a 4-exon gene that undergoes alternative RNA splicing to yield 3 mRNAs with 5' different untranslated regions. Gene transcription is induced during myoblast differentiation, producing all 3 mRNAs. We identify 3 critical promoter elements responsible for transcriptional activation in skeletal muscle, comprising paired E-boxes, a putative Stat and/or Ets element, and a MEF2 site, and muscle transcription factors myogenin and MEF2C stimulate RGMc promoter function in non-muscle cells. As these elements are conserved in RGMc genes from multiple species, our results suggest that RGMc has been a muscle-enriched gene throughout its evolutionary history.

Pax6 regulation of Math5 during mouse retinal neurogenesis

Activation of the bHLH factor Math5 (Atoh7) is an initiating event for mammalian retinal neurogenesis, as it is critically required for retinal ganglion cell formation. However, the cis-regulatory elements and trans-acting factors that control Math5 expression are largely unknown. Using a combination of transgenic mice and bioinformatics, we identified a phylogenetically conserved regulatory element that is required to activate Math5 transcription during early retinal neurogenesis. This element drives retinal expression in vivo, in a cross-species transgenic assay. Previously, Pax6 was shown to be necessary for the initiation of Math5 mRNA expression. We extend this finding by showing that the Math5 retinal enhancer also requires Pax6 for its activation, via Pax6 binding to a highly conserved binding site. In addition, our data reveal that other retinal factors are required for accurate regulation of Math5 by Pax6.

AhR-mediated gene expression in the developing mouse telencephalon

We hypothesize that TCDD-induced developmental neurotoxicity is modulated through an AhR-dependent interaction with key regulatory neuronal differentiation pathways during telencephalon development. To test this hypothesis we examined global gene expression in both dorsal and ventral telencephalon tissues in E13.5 AhR-/- and wildtype mice exposed to TCDD or vehicle. Consistent with previous biochemical, pathological and behavioral studies, our results suggest TCDD initiated changes in gene expression in the developing telencephalon are primarily AhR-dependent, as no statistically significant gene expression changes are evident after TCDD exposure in AhR-/- mice. Based on a gene regulatory network for neuronal specification in the developing telencephalon, the present analysis suggests differentiation of GABAergic neurons in the ventral telencephalon is compromised in TCDD exposed and AhR-/- mice. In addition, our analysis suggests Sox11 may be directly regulated by AhR based on gene expression and comparative genomics analyses. In conclusion, this analysis supports the hypothesis that AhR has a specific role in the normal development of the telencephalon and provides a mechanistic framework for neurodevelopmental toxicity of chemicals that perturb AhR signaling.

Evolution of MIR168 paralogs in Brassicaceae

Background

In plants, expression of ARGONAUTE1 (AGO1), the catalytic subunit of the RNA-Induced Silencing Complex responsible for post-transcriptional gene silencing, is controlled through a feedback loop involving the miR168 microRNA. This complex auto-regulatory loop, composed of miR168-guided AGO1-catalyzed cleavage of AGO1 mRNA and AGO1-mediated stabilization of miR168, was shown to ensure the maintenance of AGO1 homeostasis that is pivotal for the correct functioning of the miRNA pathway.

Results

We applied different approaches to studying the genomic organization and the structural and functional evolution of MIR168 homologs in Brassicaeae. A whole genome comparison of Arabidopsis and poplar, phylogenetic footprinting and phylogenetic reconstruction were used to date the duplication events originating MIR168 homologs in these genomes. While orthology was lacking between Arabidopsis and poplar MIR168 genes, we successfully isolated orthologs of both loci present in Arabidopsis (MIR168a and MIR168b) from all the Brassicaceae species analyzed, including the basal species Aethionema grandiflora, thus indicating that (1) independent duplication events took place in Arabidopsis and poplar lineages and (2) the origin of MIR168 paralogs predates both the Brassicaceae radiation and the Arabidopsis alpha polyploidization. Different phylogenetic footprints, corresponding to known functionally relevant regions (transcription starting site and double-stranded structures responsible for microRNA biogenesis and function) or for which functions could be proposed, were found to be highly conserved among MIR168 homologs. Comparative predictions of the identified microRNAs also indicate extreme conservation of secondary structure and thermodynamic stability.

Conclusion

We used a comparative phylogenetic footprinting approach to identify the structural and functional constraints that shaped MIR168 evolution in Brassicaceae. Although their duplication happened at least 40 million years ago, we found evidence that both MIR168 paralogs have been maintained throughout the evolution of Brassicaceae, most likely functionally as indicated by the extremely high conservation of functionally relevant regions, predicted secondary structure and thermodynamic profile. Interestingly, the expression patterns observed in Arabidopsis indicate that MIR168b underwent partial subfunctionalization as determined by the experimental characterization of its expression pattern provided in this study. We found further evolutionary evidence that pre-miR168 lower stem (the RNA-duplex structure adjacent to the miR-miR* stem) is significantly longer than animal lower stems and probably plays a relevant role in multi-step miR168 biogenesis.

In vitro analysis of promoter activity in Müller cells

PURPOSE

Rational modification of promoter architecture is necessary for manipulation of transgene activity and requires accurate deciphering of regulatory control elements. Identification of minimally sized promoters is critical to the design of viral vectors for gene therapy. To this end, we evaluated computational methods for predicting short DNA sequences capable of driving gene expression in Müller cells.

METHODS

We measured enhanced green fluorescent protein (eGFP) expression levels driven by "full-length" promoters, and compared these data with computationally identified shorter promoter elements from the same genes. We cloned and screened over 90 sequences from nine Müller cell-associated genes: CAR2, CD44, GFAP, GLUL, PDGFRA, RLBP1, S100B, SLC1A3, and vimentin (VIM). We PCR-amplified the "full-length" promoter (~1500 bp), the proximal promoter (~500 bp), and the most proximal evolutionarily conserved region (ECR; 95-871 bp) for each gene, both with and without their respective 5' untranslated regions (UTRs), from C57BL/6J mouse genomic DNA. We selected and cloned additional ECRs from more distal genomic regions (both 5' and 3') of the VIM and CD44 genes, using both mouse and rat (Sprague-Dawley) genomic DNA as templates. PCR products were cloned into the pFTMGW or pFTM3GW lentiviral transfer vectors. Plasmid constructs were transfected into rat (wMC) or human (MIO-M1) Müller cells, and eGFP expression levels were evaluated by fluorescence microscopy and flow cytometry. Selected constructs were also examined in NIH/3T3 and Neuro-2a cells.

RESULTS

Several ECRs from the nine Müller cell-associated genes were able to drive reporter gene expression as well as their longer counterparts. Preliminary comparisons of ECRs from the VIM and CD44 genes suggested that inclusion of UTRs in promoter constructs resulted in increased transgene expression levels. Systematic comparison of promoter activity from nine Müller cell-expressed genes supported this finding, and characteristic regulation profiles were evident among the different genes tested. Importantly, individual cloned promoter sequences were capable of driving distinct levels of transgene expression, resulting in up to eightfold more cells expressing eGFP with up to 3.8-fold higher mean fluorescence intensity (MFI). Furthermore, combining constructs into single regulatory "units" modulated transgene expression, suggesting that secondary gene sequences provided in cis may be used to fine-tune gene expression levels.

CONCLUSIONS

In this study, we demonstrate that computational and empirical methods, when used in combination, can efficiently identify short promoters that are active in cultured Müller cells. In addition, the pFTM3GW vector can be used to study the effects of combined promoter elements. We anticipate that these methods will expedite the design and testing of synthetic/chimeric promoter constructs that should be useful for both in vitro and in vivo applications.

Intergenic bidirectional promoter and cooperative regulation of the IIx and IIb MHC genes in fast skeletal muscle

This study investigated the dynamic regulation of IIx-IIb MHC genes in the fast white medial gastrocnemius (WMG) muscle in response to intermittent resistance exercise training (RE), a model associated with a rapid shift from IIb to IIx expression (11). We investigated the effect of 4 days of RE on the transcriptional activity across the skeletal MHC gene locus in the WMG in female Sprague-Dawley rats. Our results show that RE resulted in significant shifts from IIb to IIx observed at both the pre-mRNA and mRNA levels. An antisense RNA (xII NAT) was detected in the intergenic (IG) region between IIx and IIb, extending across the entire IIx gene and into its promoter. The expression of the xII NAT was positively correlated with IIb pre-mRNA (R = +0.8), and negatively correlated with IIx pre-mRNA (R = -0.8). Transcription mapping of the IIx-IIb IG region revealed the generation of sense IIb and xII NATs from a single promoter region. This bidirectional promoter is highly conserved among species and contains several regulatory elements that may be implicated in its regulation. These results suggest that the IIx and the IIb genes are physically and functionally linked via the bidirectional promoter. In order for the IIx MHC gene to be regulated, a feedback mechanism from the IG xII NAT is needed. In conclusion, the IG bidirectional promoter generating antisense RNA appears to be essential for the coordinated regulation of the skeletal muscle MHC genes during dynamic phenotype shifts.

Characterization of the proneural gene regulatory network during mouse telencephalon development

Background

The proneural proteins Mash1 and Ngn2 are key cell autonomous regulators of neurogenesis in the mammalian central nervous system, yet little is known about the molecular pathways regulated by these transcription factors.

Results

Here we identify the downstream effectors of proneural genes in the telencephalon using a genomic approach to analyze the transcriptome of mice that are either lacking or overexpressing proneural genes. Novel targets of Ngn2 and/or Mash1 were identified, such as members of the Notch and Wnt pathways, and proteins involved in adhesion and signal transduction. Next, we searched the non-coding sequence surrounding the predicted proneural downstream effector genes for evolutionarily conserved transcription factor binding sites associated with newly defined consensus binding sites for Ngn2 and Mash1. This allowed us to identify potential novel co-factors and co-regulators for proneural proteins, including Creb, Tcf/Lef, Pou-domain containing transcription factors, Sox9, and Mef2a. Finally, a gene regulatory network was delineated using a novel Bayesian-based algorithm that can incorporate information from diverse datasets.

Conclusion

Together, these data shed light on the molecular pathways regulated by proneural genes and demonstrate that the integration of experimentation with bioinformatics can guide both hypothesis testing and hypothesis generation.

Genetic analysis of chromosome 20-related posterior polymorphous corneal dystrophy: genetic heterogeneity and exclusion of three candidate genes

PURPOSE

Posterior polymorphous corneal dystrophy (PPCD) is a genetically heterogeneous autosomal dominant condition which maps to the pericentromeric region of chromosome 20. Mutations in the VSX1 transcription factor have been reported in patients affected with PPCD, keratoconus, or a combination of both phenotypes. However, no mutation was identified in the coding region of VSX1 in the family used for the original mapping. To clarify the genetic basis of PPCD1, a thorough analysis was performed on the original PPCD1 family and two other PPCD1-linked families. As part of the analysis, the expression profile, transcript variants, and evolutionary conserved regions of VSX1, a key candidate gene within the linkage interval, were characterized.

METHODS

Haplotype analysis was performed using highly informative markers on the pericentromeric region of chromosome 20. VSX1 transcript variants were identified using RT-PCR and characterized by 3'RACE assay. Temporal expression profile of VSX1 was evaluated using semi-quantitative real-time RT-PCR on human tissues. Evolutionary conserved regions (ECRs) were identified in the vicinity of VSX1 using publicly available sequence alignments (UCSC and rVista) and sequenced for mutation analysis.

RESULTS

Recombination events were identified that narrow the PPCD1-disease interval from 20 to 16.44 cM. This smaller interval includes the CHED1 locus and a recently described PPCD locus in Czech families. The three strongest candidate genes of the PPCD1-CHED1 overlap region (RBBP9, ZNF133, SLC24A3) did not show any mutations in our PPCD1-linked families. Semi-quantitative real-time RT-PCR detected VSX1 expression in neonatal human cornea. Six transcript variants of VSX1 were characterized. Four of the transcript variants spliced to two novel exons downstream of the gene. Mutation analysis of the PPCD1-linked families did not reveal any mutations in the full genomic sequence of VSX1 (considering all splice variants) or in the six cis- regulatory modules predicted in the vicinity of VSX1 (100 kb).

CONCLUSIONS

This is the first documentation of VSX1 expression in human neonatal cornea. We provide evidence for genetic heterogeneity of chromosome 20-related PPCD and refinement of the original PPCD1 interval. The full genomic sequence of VSX1 and coding exons of three other candidate genes were excluded from being pathogenic in the original PPCD1 family.

Comparative and evolutionary genomics of globin genes in fish

Sequencing genomes of model organisms is a great challenge for biological sciences. In the past decade, scientists have developed a large number of methods to align and compare sequenced genomes. The analysis of a given sequence provides much information on the genome structure but to a lesser extent on the function. Comparative genomics are a useful tool for functional and evolutionary annotation of genomes. In principle, comparison of genomic sequences may allow for identification of the evolutionary selection (negative or positive) that the functional sequences have been subjected to over time. Positively selected genome regions are the most important ones for evolution, because most changes are adaptive and often induce biological differences in organisms. The draft genomes of five fish species have recently become available. We herewith review and discuss some new insights into comparative genomics in fish globin genes. Special attention will be given to a complementary methodological approach to comparative genomics, fluorescence in situ hybridization (FISH). Internet resources for analyzing sequence alignments and annotations and new bioinformatic tools to address critical problems are thoroughly discussed.

A highly conserved Wnt-dependent TCF4 binding site within the proximal enhancer of the anti-myogenic Msx1 gene supports expression within Pax3-expressing limb bud muscle precursor cells

The product of the Msx1 gene is a potent inhibitor of muscle differentiation. Msx1 is expressed in muscle precursor cells of the limb bud that also express Pax3. It is thought that Msx1 may facilitate distal migration by delaying myogenesis in these cells. Despite the role played by Msx1 in inhibiting muscle differentiation, nothing is known of the mechanisms that support the expression of the Msx1 gene within limb bud muscle precursor cells. In the present study we have used a combination of comparative genomics, mouse transgenic analysis, in situ hybridisation and immunohistochemistry to identify a highly conserved and tissue-specific regulatory sub-domain within the previously characterised Msx1 gene proximal enhancer element that supports the expression of the Msx1 gene in Pax3-expressing mouse limb pre-muscle masses. Furthermore, using a combination of in situ hybridisation, in vivo ChIP assay and transgenic explant culture analysis we provide evidence that Msx1 expression in limb bud muscle precursor cells is dependent on the canonical Wnt/TCF signalling pathway that is important in muscle shape formation. The results of these studies provide evidence of a mechanistic link between the Wnt/TCF and the Msx1/Pax3/MyoD pathways within limb bud muscle precursor cells.

Functional promoter testing using a modified lentiviral transfer vector

PURPOSE

The importance of retinal glial cells in the maintenance of retinal health and in retinal degenerations has not been fully explored. Several groups have suggested that secretion of neurotrophic proteins from the retina's primary glial cell type, the Müller cell, holds promise for treating retinal degenerations. Tight regulation of transgene expression in Müller cells is likely to be critical to the efficacy of long-term neuroprotective therapies, due to the genetic heterogeneity and progressive nature of retinal disease. To this end, we developed a modified lentiviral (LV) transfer vector (pFTMGW) to accelerate the testing and evaluation of novel transcriptional regulatory elements. This vector facilitates identification and characterization of regulatory elements in terms of size, cell specificity and ability to control transgene expression levels.

METHODS

A synthetic multiple cloning site (MCS) which can accept up to five directionally cloned DNA regulatory elements was inserted immediately upstream of an enhanced green fluorescent protein (eGFP) reporter. A cytomegalovirus (CMV) promoter, required for tat-independent viral packaging, is located around 2 kb upstream of the eGFP reporter and is capable of directing transgene expression. A synthetic transcription blocker (TB) was inserted to insulate the MCS/eGFP from the CMV promoter. We evaluated eGFP expression from pFTMGW and control constructs using flow cytometry and quantitative reverse transcriptase polymerase chain reaction (RT-PCR). We also tested and compared the activity and cell specificity of a computationally identified promoter fragment from the rat vimentin gene (Vim409) in transfection and lentiviral infection experiments using fluorescence microscopy.

RESULTS

Transfection data, quantitative RT-PCR, and flow cytometry show that around 85% of expression from the CMV promoter was blocked by the TB element, allowing direct evaluation of expression from the Vim409 candidate promoter cloned into the MCS. Lentiviruses generated from this construct containing the Vim409 promoter (without the TB element) drove robust eGFP expression in Müller cells in vitro and in vivo.

CONCLUSIONS

The TB element efficiently prevented eGFP expression by the upstream CMV promoter and the novel MCS facilitated testing of an evolutionarily conserved regulatory element. Additional sites allow for combinatorial testing of additional promoter, enhancer, and/or repressor elements in various configurations. This modified LV transfer vector is an effective tool for expediting functional analysis of gene regulatory elements in Müller glia, and should prove useful for promoter analyses in other cell types and tissues.

Enhancer identification through comparative genomics

With the availability of genomic sequence from numerous vertebrates, a paradigm shift has occurred in the identification of distant-acting gene regulatory elements. In contrast to traditional gene-centric studies in which investigators randomly scanned genomic fragments that flank genes of interest in functional assays, the modern approach begins electronically with publicly available comparative sequence datasets that provide investigators with prioritized lists of putative functional sequences based on their evolutionary conservation. However, although a large number of tools and resources are now available, application of comparative genomic approaches remains far from trivial. In particular, it requires users to dynamically consider the species and methods for comparison depending on the specific biological question under investigation. While there is currently no single general rule to this end, it is clear that when applied appropriately, comparative genomic approaches exponentially increase our power in generating biological hypotheses for subsequent experimental testing. It is anticipated that cardiac-related genes and the identification of their distant-acting transcriptional enhancers are particularly poised to benefit from these modern capabilities.

Phylogenetic analysis of the promoter region of the CD40L gene in primates and other mammals

CD40L is a type II membrane protein comprised of 261 amino acids. CD40L plays a crucial role in the immune system where it is primarily expressed on activated T cells and triggers immunoglobulin class switching. The genetic disease X-linked hypergammaglobulinemia (HIGM1, XHIGM or XHIM) is caused by mutations in the CD40L gene. Individuals with HIGM1 are susceptible to recurrent infections to pathogens and a relationship has been shown to exist with malaria [Sabeti, P., Usen, S., Farhadian, S., Jallow, M., Doherty, T., Newport, M., Pinder, M., Ward, R., Kwiatkowski, D., 2002a. CD40L association with protection from severe malaria. Genes Immun. 3, 286-291]. In this paper, we phylogenetically examine the promoter region of CD40L in primates and other mammals via phylogenetic shadowing. This analysis revealed several regions of the CD40L promoter that were highly constrained and thereby inferred to be functional. These constrained regions confirmed many known regulatory sites. In addition, a novel, highly constrained upstream region was also identified which had an NF-AT recognition motif. These analyses also showed that the different mammal groups do not share an exactly similar set of promoter binding sites and taxon-specific promoters have evolved.

VISTA Enhancer Browser--a database of tissue-specific human enhancers

Despite the known existence of distant-acting cis-regulatory elements in the human genome, only a small fraction of these elements has been identified and experimentally characterized in vivo. This paucity of enhancer collections with defined activities has thus hindered computational approaches for the genome-wide prediction of enhancers and their functions. To fill this void, we utilize comparative genome analysis to identify candidate enhancer elements in the human genome coupled with the experimental determination of their in vivo enhancer activity in transgenic mice [L. A. Pennacchio et al. (2006) Nature, in press]. These data are available through the VISTA Enhancer Browser (). This growing database currently contains over 250 experimentally tested DNA fragments, of which more than 100 have been validated as tissue-specific enhancers. For each positive enhancer, we provide digital images of whole-mount embryo staining at embryonic day 11.5 and an anatomical description of the reporter gene expression pattern. Users can retrieve elements near single genes of interest, search for enhancers that target reporter gene expression to a particular tissue, or download entire collections of enhancers with a defined tissue specificity or conservation depth. These experimentally validated training sets are expected to provide a basis for a wide range of downstream computational and functional studies of enhancer function.

Regulation of mouse inducible costimulator (ICOS) expression by Fyn-NFATc2 and ERK signaling in T cells

The inducible costimulator (ICOS), a member of the CD28 family of costimulatory molecules, is rapidly induced upon T cell activation. Although the critical role of ICOS in costimulating T cell responses is well documented, little is known of the intracellular signaling pathways and mechanisms that regulate ICOS expression. Here, we report that Fyn, NFAT, and ERK signaling influence ICOS expression as various chemical inhibitors, such as PP2 that targets Src kinases, U0126 that targets MEK1/2, and cyclosporin A or FK506 that targets calcineurin and thereby affects NFAT, attenuate T cell receptor-mediated ICOS induction. Moreover, ectopic expression of NFATc2 or a constitutively active MEK2 amplifies ICOS transcription and transactivates a 288-bp core region of the icos promoter in luciferase reporter assays. We also identify a site on the icos promoter that is sensitive to ERK signaling and further show that NFATc2 can bind the icos promoter in vivo and that this binding is diminished when Fyn signaling is ablated. The normal activation of ERK but reduced nuclear translocation of NFATc2 in Fyn(-/-) CD4(+) T cells further suggest that Fyn and NFATc2 act in a common axis, separate from that involving ERK, to drive ICOS transcription. Taken together, our findings indicate that Fyn-calcineurin-NFATc2 and MEK2-ERK1/2 are two independent signaling pathways that cooperate to control T cell receptor-mediated ICOS induction.

Approaches to lipid metabolism gene identification and characterization in the postgenomic era

The availability of genomic resources has already had a tremendous impact on biomedical research. In this review, we describe how whole genome sequence and high-throughput functional genomics projects have facilitated the identification and characterization of important genes in lipid metabolism and disease. We review key approaches and lipid genes identified in the first years of this century and discuss how genomic resources are likely to streamline gene identification and functional characterization in the future.

New technologies, new findings, and new concepts in the study of vertebrate cis-regulatory sequences

All vertebrates share a similar early embryonic body plan and use the same regulatory genes for their development. The availability of numerous sequenced vertebrate genomes and significant advances in bioinformatics have resulted in the finding that the genomic regions of many of these developmental regulatory genes also contain highly conserved noncoding sequence. In silico discovery of conserved noncoding regions and of transcription factor binding sites as well as the development of methods for high throughput transgenesis in Xenopus and zebrafish are dramatically increasing the speed with which regulatory elements can be discovered, characterized, and tested in the context of whole live embryos. We review here some of the recent technological developments that will likely lead to a surge in research on how vertebrate genomes encode regulation of transcriptional activity, how regulatory sequences constrain genomic architecture, and ultimately how vertebrate form has evolved.

Array2BIO: from microarray expression data to functional annotation of co-regulated genes

BACKGROUND

There are several isolated tools for partial analysis of microarray expression data. To provide an integrative, easy-to-use and automated toolkit for the analysis of Affymetrix microarray expression data we have developed Array2BIO, an application that couples several analytical methods into a single web based utility.

RESULTS

Array2BIO converts raw intensities into probe expression values, automatically maps those to genes, and subsequently identifies groups of co-expressed genes using two complementary approaches: (1) comparative analysis of signal versus control and (2) clustering analysis of gene expression across different conditions. The identified genes are assigned to functional categories based on Gene Ontology classification and KEGG protein interaction pathways. Array2BIO reliably handles low-expressor genes and provides a set of statistical methods for quantifying expression levels, including Benjamini-Hochberg and Bonferroni multiple testing corrections. An automated interface with the ECR Browser provides evolutionary conservation analysis for the identified gene loci while the interconnection with Crème allows prediction of gene regulatory elements that underlie observed expression patterns.

CONCLUSION

We have developed Array2BIO - a web based tool for rapid comprehensive analysis of Affymetrix microarray expression data, which also allows users to link expression data to Dcode.org comparative genomics tools and integrates a system for translating co-expression data into mechanisms of gene co-regulation. Array2BIO is publicly available at http://array2bio.dcode.org.