Sequence comparison of human and mouse genes reveals a homologous block structure in the promoter regions

Retrotrans-genomics identifies aberrant THE1B endogenous retrovirus fusion transcripts in the pathogenesis of sarcoidosis

Transposon-like human element 1B (THE1B) originates from ancient retroviral sequences integrated into the primate genome approximately 50 million years ago, now accounting for at least 27,233 copies in the human genome, suggesting their extensive influence on human genomic architecture. Here we report identification of 19 THE1B fusion transcripts through short- and long-read RNA-seq analysis, 15 of which are previously unmapped, showing elevated expression in 16 individuals with sarcoid myopathy (SM), as compared to 400 controls with various other muscle diseases. Analysis of publicly available RNA-seq data indicated a correlation between the reduced expression of eight THE1B fusion transcripts and clinical improvement in individuals with cutaneous sarcoidosis receiving tofacitinib treatment. Single-cell or single-nucleus RNA-seq analyses of sarcoidosis not only confirmed these transcripts but also revealed a novel read-through transcript, SIRPB1-SIRPD, and TREM2.1, predominantly in granuloma-associated macrophages. The expression profiles of THE1B fusion transcripts in tuberculosis (TB) significantly differed from SM in single-cell RNA-seq data, suggesting that the differences between TB's caseous granulomas and sarcoidosis's non-caseous granulomas might be linked to disparate expression patterns of THE1B fusion transcripts. Our retrotrans-genomics approach has not only identified the genomic landscape of sarcoidosis but also provided new insights into its etiology.

Deep learning-based assessment of missense variants in the COG4 gene presented with bilateral congenital cataract

OBJECTIVE

We compared the protein structure and pathogenicity of clinically relevant variants of the COG4 gene with AlphaFold2 (AF2), Alpha Missense (AM), and ThermoMPNN for the first time.

METHODS AND ANALYSIS

The sequences of clinically relevant Cog4 missense variants (one novel identified p.Y714F and three pre-existing p.G512R, p.R729W and p.L769R from Uniprot Q9H9E3) were imported into AF2 for protein structural prediction, and the pathogenicity was estimated using AM and ThermoMPNN. Different pathogenicity metrics were aggregated with principal component analysis (PCA) and further analysed at three levels (amino acid position, substitution and post-translation) based on all possible Cog4 missense variants (n=14 915).

RESULTS

Localised protein structural impact including change of conformation and amino acid polarity, breakage of hydrogen bond and salt-bridge, and formation of alpha-helix were identified among clinically relevant Cog4 variants. The global structural comparison with multidimensional scaling demonstrated variants with similar protein structures (AF2) tended to exhibit similar clinical and biological phenotypes. The Cog4 p.Y714F variant exhibited greater protein structural similarity to mutated Cog4 found in Saul‒Wilson syndrome (p.G512R) and shared similar clinical phenotype (congenital cataract and psychomotor retardation). PCA of included pathogenic metrics demonstrated p.Y714F occurred at a critical position in Cog4 amino acid sequence with disrupted post-translational phosphorylation.

CONCLUSION

Deep learning algorithms, including AF2, AM and ThermoMPNN, can be useful for evaluating variant of uncertain significance (VUS) by structural and pathogenicity prediction. Despite classified as VUS (American College of Medical Genetics and Genomics criteria: PM1, PP4), the pathogenicity in this Cog4 variant cannot be ruled out and warrants further investigation.

Therapeutic advances in the targeting of ROR1 in hematological cancers

Receptor tyrosine kinases (RTKs) are key cell surface receptors involved in cell communication and signal transduction, with great importance in cell growth, differentiation, survival, and metabolism. Dysregulation of RTKs, such as EGFR, VEGFR, HER2 or ROR, could lead to various diseases, particularly cancers. ROR1 has emerged as a promising target in hematological malignancies. The development of ROR1 targeted therapies is continuously growing leading to remarkable novel therapeutical approaches using mAbs, antibody-drug conjugates, several small molecules or CAR T cells which have shown encouraging preclinical results. In the hematological field, mAbs, small molecules, BiTEs or CAR T cell therapies displayed promising outcomes with the clinical trials data encouraging the use of anti-ROR1 therapies. This paper aims to offer a comprehensive analysis of the current landscape of ROR1-targeted therapies in hematological malignancies marking the innovative approaches with promising preclinical and clinical. Offering a better understanding of structural and functional aspects of ROR1 could lead to new perspectives in targeting a wide spectrum of malignancies.

Genomic binding of NF-Y in mouse and human cells

NF-Y is a Transcription Factor that regulates transcription through binding to the CCAAT-box. To understand its strategy, we analyzed 16 ChIP-seq datasets from human and mouse cells. Shared loci, mostly located in promoters of expressed genes of cell cycle, metabolism and gene expression pathways, are associated with histone marks of active chromatin and specific modules of TFs. Other peaks are in enhancers and Transposable Elements -TE- of retroviral origin in human and mouse. We evaluated the relationship with USF1, a common synergistic partner in promoters and MLT1 TEs, upon NF-YB inactivation: USF1 binding decreases in promoters, modestly in MLT1, suggesting a pioneering role of NF-Y in formers, not in the latters. These data define a common set of NF-Y functional targets across different mammalian cell types, suggesting a pioneering role in promoters with respect to TEs.

From zinc homeostasis to disease progression: Unveiling the neurodegenerative puzzle

Zinc is a crucial trace element in the human body, playing a role in various physiological processes such as oxidative stress, neurotransmission, protein synthesis, and DNA repair. The zinc transporters (ZnTs) family members are responsible for exporting intracellular zinc, while Zrt- and Irt-like proteins (ZIPs) are involved in importing extracellular zinc. These processes are essential for maintaining cellular zinc homeostasis. Imbalances in zinc metabolism have been linked to the development of neurodegenerative diseases. Disruptions in zinc levels can impact the survival and activity of neurons, thereby contributing to the progression of neurodegenerative diseases through mechanisms like cell apoptosis regulation, protein phase separation, ferroptosis, oxidative stress, and neuroinflammation. Therefore, conducting a systematic review of the regulatory network of zinc and investigating the relationship between zinc dysmetabolism and neurodegenerative diseases can enhance our understanding of the pathogenesis of these diseases. Additionally, it may offer new insights and approaches for the treatment of neurodegenerative diseases.

The Identification of Large Rearrangements Involving Intron 2 of the CDH1 Gene in BRCA1/2 Negative and Breast Cancer Susceptibility

E-cadherin, a CDH1 gene product, is a calcium-dependent cell-cell adhesion molecule playing a critical role in the establishment of epithelial architecture, maintenance of cell polarity, and differentiation. Germline pathogenic variants in the CDH1 gene are associated with hereditary diffuse gastric cancer (HDGC), and large rearrangements in the CDH1 gene are now being reported as well. Because CDH1 pathogenic variants could be associated with breast cancer (BC) susceptibility, CDH1 rearrangements could also impact it. The aim of our study is to identify rearrangements in the CDH1 gene in 148 BC cases with no BRCA1 and BRCA2 pathogenic variants. To do so, a zoom-in CGH array, covering the exonic, intronic, and flanking regions of the CDH1 gene, was used to screen our cohort. Intron 2 of the CDH1 gene was specifically targeted because it is largely reported to include several regulatory regions. As results, we detected one large rearrangement causing a premature stop in exon 3 of the CDH1 gene in a proband with a bilateral lobular breast carcinoma and a gastric carcinoma (GC). Two large rearrangements in the intron 2, a deletion and a duplication, were also reported only with BC cases without any familial history of GC. No germline rearrangements in the CDH1 coding region were detected in those families without GC and with a broad range of BC susceptibility. This study confirms the diversity of large rearrangements in the CDH1 gene. The rearrangements identified in intron 2 highlight the putative role of this intron in CDH1 regulation and alternative transcripts. Recurrent duplication copy number variations (CNV) are found in this region, and the deletion encompasses an alternative CDH1 transcript. Screening for large rearrangements in the CDH1 gene could be important for genetic testing of BC.

On the NF-Y regulome as in ENCODE (2019)

NF-Y is a trimeric Transcription Factor -TF- which binds with high selectivity to the conserved CCAAT element. Individual ChIP-seq analysis as well as ENCODE have progressively identified locations shared by other TFs. Here, we have analyzed data introduced by ENCODE over the last five years in K562, HeLa-S3 and GM12878, including several chromatin features, as well RNA-seq profiling of HeLa cells after NF-Y inactivation. We double the number of sequence-specific TFs and co-factors reported. We catalogue them in 4 classes based on co-association criteria, infer target genes categorizations, identify positional bias of binding sites and gene expression changes. Larger and novel co-associations emerge, specifically concerning subunits of repressive complexes as well as RNA-binding proteins. On the one hand, these data better define NF-Y association with single members of major classes of TFs, on the other, they suggest that it might have a wider role in the control of mRNA production.

The ongoing ENCODE consortium represents a useful compendium of locations of TFs, chromatin marks, gene expression data. In previous reports, we identified modules of CCAAT-binding NF-Y with individual TFs. Here, we analyzed all 363 factors currently present: 68 with enrichment of CCAAT in their locations, 38 with overlap of peaks. New sequence-specific TFs, co-activators and co-repressors are reported. Co-association patterns correspond to specific targeted genes categorizations and gene expression changes, as assessed by RNA-seq after NF-Y inactivation. These data widen and better define a coherent model of synergy of NF-Y with selected groups of TFs and co-factors.

Magnetic Resonance Imaging for Translational Research in Oncology

The translation of results from the preclinical to the clinical setting is often anything other than straightforward. Indeed, ideas and even very intriguing results obtained at all levels of preclinical research, i.e., in vitro, on animal models, or even in clinical trials, often require much effort to validate, and sometimes, even useful data are lost or are demonstrated to be inapplicable in the clinic. In vivo, small-animal, preclinical imaging uses almost the same technologies in terms of hardware and software settings as for human patients, and hence, might result in a more rapid translation. In this perspective, magnetic resonance imaging might be the most translatable technique, since only in rare cases does it require the use of contrast agents, and when not, sequences developed in the lab can be readily applied to patients, thanks to their non-invasiveness. The wide range of sequences can give much useful information on the anatomy and pathophysiology of oncologic lesions in different body districts. This review aims to underline the versatility of this imaging technique and its various approaches, reporting the latest preclinical studies on thyroid, breast, and prostate cancers, both on small laboratory animals and on human patients, according to our previous and ongoing research lines.

Characterization of gene promoters in pig: conservative elements, regulatory motifs and evolutionary trend

It is vital to understand the conservation and evolution of gene promoter sequences in order to understand environmental adaptation. The level of promoter conservation varies greatly between housekeeping (HK) and tissue-specific (TS) genes, denoting differences in the strength of the evolutionary constraints. Here, we analyzed promoter conservation and evolution to exploit differential regulation between HK and TS genes. The analysis of conserved elements showed CpG islands, short tandem repeats and G-quadruplex sequences are highly enriched in HK promoters relative to TS promoters. In addition, the type and density of regulatory motifs in TS promoters are much higher than HK promoters, indicating that TS genes show more complex regulatory patterns than HK genes. Moreover, the evolutionary dynamics of promoters showed similar evolutionary trend to coding sequences. HK promoters suffer more stringent selective pressure in the long-term evolutionary process. HK genes tend to show increased upstream sequence conservation due to stringent selection pressures acting on the promoter regions. The specificity of TS gene expression may be due to complex regulatory motifs acting in different tissues or conditions. The results from this study can be used to deepen our understanding of adaptive evolution.

Isoforms of Base Excision Repair Enzymes Produced by Alternative Splicing

Transcripts of many enzymes involved in base excision repair (BER) undergo extensive alternative splicing, but functions of the corresponding alternative splice variants remain largely unexplored. In this review, we cover the studies describing the common alternatively spliced isoforms and disease-associated variants of DNA glycosylases, AP-endonuclease 1, and DNA polymerase beta. We also discuss the roles of alternative splicing in the regulation of their expression, catalytic activities, and intracellular transport.

Porcine Models of Pancreatic Cancer

Pancreatic cancer is the fourth most common cause of cancer-related deaths in both men and women. The 5-year survival rate for metastatic pancreatic cancer is only 8%. There remains a need for improved early diagnosis and therapy for pancreatic cancer. Murine models are the current standard for preclinical study of pancreatic cancer. However, mice may not accurately reflect human biology because of a variety of differences between the two species. Remarkably, only 5-8% of anti-cancer drugs that have emerged from preclinical studies and entered clinical studies have ultimately been approved for clinical use. The cause of this poor approval rate is multi-factorial, but may in part be due to use of murine models that have limited accuracy with respect to human disease. Murine models also have limited utility in the development of diagnostic or interventional technology that require a human-sized model. So, at present, there remains a need for improved animal models of pancreatic cancer. The rationale for a porcine model of pancreatic cancer is (i) to enable development of diagnostic/therapeutic devices for which murine models have limited utility; and (ii) to have a highly predictive preclinical model in which anti-cancer therapies can be tested and optimized prior to a clinical trial. Recently, pancreatic tumors were induced in transgenic Oncopigs and porcine pancreatic ductal cells were transformed that contain oncogenic KRAS and p53-null mutations. Both techniques to induce pancreatic tumors in pigs are undergoing further refinement and expansion. The Oncopig currently is commercially available, and it is conceivable that other porcine models of pancreatic cancer may be available for general use in the near future.

Nonhuman Primate Optogenetics: Recent Advances and Future Directions

Optogenetics is the use of genetically coded, light-gated ion channels or pumps (opsins) for millisecond resolution control of neural activity. By targeting opsin expression to specific cell types and neuronal pathways, optogenetics can expand our understanding of the neural basis of normal and pathological behavior. To maximize the potential of optogenetics to study human cognition and behavior, optogenetics should be applied to the study of nonhuman primates (NHPs). The homology between NHPs and humans makes these animals the best experimental model for understanding human brain function and dysfunction. Moreover, for genetic tools to have translational promise, their use must be demonstrated effectively in large, wild-type animals such as Rhesus macaques. Here, we review recent advances in primate optogenetics. We highlight the technical hurdles that have been cleared, challenges that remain, and summarize how optogenetic experiments are expanding our understanding of primate brain function.

Normal and Cancerous Tissues Release Extrachromosomal Circular DNA (eccDNA) into the Circulation

Cell-free circulating linear DNA is being explored for noninvasive diagnosis and management of tumors and fetuses, the so-called liquid biopsy. Previously, we observed the presence of small extrachromosomal circular DNA (eccDNA), called microDNA, in the nuclei of mammalian tissues and cell lines. Now, we demonstrate that cell-free microDNA derived from uniquely mapping regions of the genome is detectable in plasma and serum from both mice and humans and that they are significantly longer (30%-60% >250 bases) than cell-free circulating linear DNA (∼150 bases). Tumor-derived human microDNA is detected in the mouse circulation in a mouse xenograft model of human ovarian cancer. Comparing the microDNA from paired tumor and normal lung tissue specimens reveals that the tumors contain longer microDNA. Consistent with human cancers releasing microDNA into the circulation, serum and plasma samples (12 lung and 11 ovarian cancer) collected prior to surgery are enriched for longer cell-free microDNA compared with samples from the same patient obtained several weeks after surgical resection of the tumor. Thus, circular DNA in the circulation is a previously unexplored pool of nucleic acids that could complement miRNAs and linear DNA for diagnosis and for intercellular communication.Implications: eccDNA derived from chromosomal genomic sequence, first discovered in the nuclei of cells, are detected in the circulation, are longer than linear cell-free DNA, and are released from normal tissue and tumors into the circulation. Mol Cancer Res; 15(9); 1197-205. ©2017 AACR.

Generation and characterization of two immortalized human osteoblastic cell lines useful for epigenetic studies

Different model systems using osteoblastic cell lines have been developed to help understand the process of bone formation. Here, we report the establishment of two human osteoblastic cell lines obtained from primary cultures upon transduction of immortalizing genes. The resulting cell lines had no major differences to their parental lines in their gene expression profiles. Similar to primary osteoblastic cells, osteocalcin transcription increased following 1,25-dihydroxyvitamin D3 treatment and the immortalized cells formed a mineralized matrix, as detected by Alizarin Red staining. Moreover, these human cell lines responded by upregulating ALPL gene expression after treatment with the demethylating agent 5-aza-2'-deoxycytidine (AzadC), as shown before for primary osteoblasts. We further demonstrate that these cell lines can differentiate in vivo, using a hydroxyapatite/tricalcium phosphate composite as a scaffold, to produce bone matrix. More importantly, we show that these cells respond to demethylating treatment, as shown by the increase in SOST mRNA levels, the gene encoding sclerostin, upon treatment of the recipient mice with AzadC. This also confirms, in vivo, the role of DNA methylation in the regulation of SOST expression previously shown in vitro. Altogether our results show that these immortalized cell lines constitute a particularly useful model system to obtain further insight into bone homeostasis, and particularly into the epigenetic mechanisms regulating sclerostin production.

A Five-Gene Expression Signature Predicts Clinical Outcome of Ovarian Serous Cystadenocarcinoma

Ovarian serous cystadenocarcinoma is a common malignant tumor of female genital organs. Treatment is generally less effective as patients are usually diagnosed in the late stage. Therefore, a well-designed prognostic marker provides valuable data for optimizing therapy. In this study, we analyzed 303 samples of ovarian serous cystadenocarcinoma and the corresponding RNA-seq data. We observed the correlation between gene expression and patients' survival and eventually established a risk assessment model of five factors using Cox proportional hazards regression analysis. We found that the survival time in high-risk patients was significantly shorter than in low-risk patients in both training and testing sets after Kaplan-Meier analysis. The AUROC value was 0.67 when predicting the survival time in testing set, which indicates a relatively high specificity and sensitivity. The results suggest diagnostic and therapeutic applications of our five-gene model for ovarian serous cystadenocarcinoma.

A high definition look at the NF-Y regulome reveals genome-wide associations with selected transcription factors

NF-Y is a trimeric transcription factor (TF), binding the CCAAT box element, for which several results suggest a pioneering role in activation of transcription. In this work, we integrated 380 ENCODE ChIP-Seq experiments for 154 TFs and cofactors with sequence analysis, protein–protein interactions and RNA profiling data, in order to identify genome-wide regulatory modules resulting from the co-association of NF-Y with other TFs. We identified three main degrees of co-association with NF-Y for sequence-specific TFs. In the most relevant one, we found TFs having a significant overlap with NF-Y in their DNA binding loci, some with a precise spacing of binding sites with respect to the CCAAT box, others (FOS, Sp1/2, RFX5, IRF3, PBX3) mostly lacking their canonical binding site and bound to arrays of well spaced CCAAT boxes. As expected, NF-Y binding also correlates with RNA Pol II General TFs and with subunits of complexes involved in the control of H3K4 methylations. Co-association patterns are confirmed by protein–protein interactions, and correspond to specific functional categorizations and expression level changes of target genes following NF-Y inactivation. These data define genome-wide rules for the organization of NF-Y-centered regulatory modules, supporting a model of distinct categorization and synergy with well defined sets of TFs.

Cancer specificity of promoters of the genes controlling cell proliferation

Violation of proliferation control is a common feature of cancer cells. We put forward the hypothesis that promoters of genes involved in the control of cell proliferation should possess intrinsic cancer specific activity. We cloned promoter regions of CDC6, POLD1, CKS1B, MCM2, and PLK1 genes into pGL3 reporter vector and studied their ability to drive heterologous gene expression in transfected cancer cells of different origin and in normal human fibroblasts. Each promoter was cloned in short (335-800 bp) and long (up to 2.3 kb) variants to cover probable location of core and whole promoter regulatory elements. Cloned promoters were significantly more active in cancer cells than in normal fibroblasts that may indicate their cancer specificity. Both versions of CDC6 promoters were shown to be most active while the activities of others were close to that of BIRC5 gene (survivin) gene promoter. Long and short variants of each cloned promoter demonstrated very similar cancer specificity with the exception of PLK1-long promoter that was substantially more specific than its short variant and other promoters under study. The data indicate that most of the important cis-regulatory transcription elements responsible for intrinsic cancer specificity are located in short variants of the promoters under study. CDC6 short promoter may serve as a promising candidate for transcription targeted cancer gene therapy.

Secreted histidyl-tRNA synthetase splice variants elaborate major epitopes for autoantibodies in inflammatory myositis

Inflammatory and debilitating myositis and interstitial lung disease are commonly associated with autoantibodies (anti-Jo-1 antibodies) to cytoplasmic histidyl-tRNA synthetase (HisRS). Anti-Jo-1 antibodies from different disease-afflicted patients react mostly with spatially separated epitopes in the three-dimensional structure of human HisRS. We noted that two HisRS splice variants (SVs) include these spatially separated regions, but each SV lacks the HisRS catalytic domain. Despite the large deletions, the two SVs cross-react with a substantial population of anti-Jo-l antibodies from myositis patients. Moreover, expression of at least one of the SVs is up-regulated in dermatomyositis patients, and cell-based experiments show that both SVs and HisRS can be secreted. We suggest that, in patients with inflammatory myositis, anti-Jo-1 antibodies may have extracellular activity.

The novel long noncoding RNA linc00467 promotes cell survival but is down-regulated by N-Myc

The worst subtype of neuroblastoma is caused by MYCN oncogene amplification and N-Myc oncoprotein over-expression. Long noncoding RNAs (lncRNAs) are emerging as critical regulators of gene expression and tumourigenesis. While Myc oncoproteins are well-known to exert tumourigenic effects by regulating the expression of protein-coding genes and microRNAs, little is known about which lncRNAs are Myc targets and whether the Myc target lncRNAs play a role in Myc-induced oncogenesis. Here we performed differential gene expression studies using lncRNA microarray in neuroblastoma cells after transfection with control or N-Myc-specific small interfering RNA (siRNA), and identified N-Myc target lncRNAs including the novel lncRNA linc00467, the expression and function of which were completely unknown. RT-PCR, chromatin immunoprecipitation and luciferase assays showed that N-Myc suppressed linc00467 gene expression through direct binding to the linc00467 gene promoter and reducing linc00467 promoter activity. While N-Myc suppressed the expression of RD3, the protein-coding gene immediately down-stream of linc00467 gene, through direct binding to the RD3 gene promoter and reducing RD3 promoter activity, linc00467 reduced RD3 mRNA expression. Moreover, Affymetrix microarray analysis revealed that one of genes significantly up-regulated by linc00467 siRNA was the tumour suppressor gene DKK1. Importantly, knocking-down linc00467 expression with siRNA in neuroblastoma cells reduced the number of viable cells and increased the percentage of apoptotic cells, and co-transfection with DKK1 siRNA blocked the effects. These findings therefore demonstrate that N-Myc-mediated suppression of linc00467 gene transcription counterintuitively blocks N-Myc-mediated reduction in RD3 mRNA expression, and reduces neuroblastoma cell survival by inducing DKK1 expression.

Comparison of genome-wide DNA methylation in urothelial carcinomas of patients with and without arsenic exposure

BACKGROUND

Arsenic is a well-documented carcinogen of human urothelial carcinoma (UC) with incompletely understood mechanisms.

OBJECTIVES

This study aimed to compare the genome-wide DNA methylation profiles of arsenic-induced UC (AsUC) and non-arsenic-induced UC (Non-AsUC), and to assess associations between site-specific methylation levels and cumulative arsenic exposure.

METHODS

Genome-wide DNA methylation profiles in 14 AsUC and 14 non-AsUC were analyzed by Illumina Infinium methylation27 BeadChip and validated by bisulfite pyrosequencing. Mean methylation levels (β¯) in AsUC and non-AsUC were compared by their ratio (β¯ ratio) and difference (Δβ¯). Associations between site-specific methylation levels in UC and cumulative arsenic exposure were examined.

RESULTS

Among 27,578 methylation sites analyzed, 231 sites had β¯ ratio >2 or <0.5 and 45 sites had Δβ¯ >0.2 or <-0.2. There were 13 sites showing statistically significant (q<0.05) differences in β¯ between AsUC and non-AsUC including 12 hypermethylation sites in AsUC and only one hypermethylation site in non-AsUC. Significant associations between cumulative arsenic exposure and DNA methylation levels of 28 patients were observed in nine CpG sites of nine gens including PDGFD (Spearman rank correlation, 0.54), CTNNA2 (0.48), KCNK17 (0.52), PCDHB2 (0.57), ZNF132 (0.48), DCDC2 (0.48), KLK7 (0.48), FBXO39 (0.49), and NPY2R (0.45). These associations remained statistically significant for CpG sites in CTNNA2, KLK7, NPY2R, ZNF132 and KCNK17 in 20 non-smoking women after adjustment for tumor stage and age.

CONCLUSIONS

Significant associations between cumulative arsenic exposure and methylation level of CTNNA2, KLK7, NPY2R, ZNF132 and KCNK17 were found in smoking-unrelated urothelial carcinoma. Arsenic exposure may cause urothelial carcinomas through the hypermethylation of genes involved in cell adhesion, proteolysis, transcriptional regulation, neuronal pathway, and ion transport. The findings of this study, which are limited by its small sample size and moderate dose-response relation, remain to be validated by further studies with large sample sizes.

Targeting the Y/CCAAT box in cancer: YB-1 (YBX1) or NF-Y?

The Y box is an important sequence motif found in promoters and enhancers containing a CCAAT box - one of the few elements enriched in promoters of large sets of genes overexpressed in cancer. The search for the transcription factor(s) acting on it led to the biochemical purification of the nuclear factor Y (NF-Y) heterotrimer, and to the cloning - through the screening of expression libraries - of Y box-binding protein 1 (YB-1), an oncogene, overexpressed in aggressive tumors and associated with drug resistance. These two factors have been associated with Y/CCAAT-dependent activation of numerous growth-related genes, notably multidrug resistance protein 1. We review two decades of data indicating that NF-Y ultimately acts on Y/CCAAT in cancer cells, a notion recently confirmed by genome-wide data. Other features of YB-1, such as post-transcriptional control of mRNA biology, render it important in cancer biology.

NF-Y and the transcriptional activation of CCAAT promoters

The CCAAT box promoter element and NF-Y, the transcription factor (TF) that binds to it, were among the first cis-elements and trans-acting factors identified; their interplay is required for transcriptional activation of a sizeable number of eukaryotic genes. NF-Y consists of three evolutionarily conserved subunits: a dimer of NF-YB and NF-YC which closely resembles a histone, and the "innovative" NF-YA. In this review, we will provide an update on the functional and biological features that make NF-Y a fundamental link between chromatin and transcription. The last 25 years have witnessed a spectacular increase in our knowledge of how genes are regulated: from the identification of cis-acting sequences in promoters and enhancers, and the biochemical characterization of the corresponding TFs, to the merging of chromatin studies with the investigation of enzymatic machines that regulate epigenetic states. Originally identified and studied in yeast and mammals, NF-Y - also termed CBF and CP1 - is composed of three subunits, NF-YA, NF-YB and NF-YC. The complex recognizes the CCAAT pentanucleotide and specific flanking nucleotides with high specificity (Dorn et al., 1997; Hatamochi et al., 1988; Hooft van Huijsduijnen et al, 1987; Kim & Sheffery, 1990). A compelling set of bioinformatics studies clarified that the NF-Y preferred binding site is one of the most frequent promoter elements (Suzuki et al., 2001, 2004; Elkon et al., 2003; Mariño-Ramírez et al., 2004; FitzGerald et al., 2004; Linhart et al., 2005; Zhu et al., 2005; Lee et al., 2007; Abnizova et al., 2007; Grskovic et al., 2007; Halperin et al., 2009; Häkkinen et al., 2011). The same consensus, as determined by mutagenesis and SELEX studies (Bi et al., 1997), was also retrieved in ChIP-on-chip analysis (Testa et al., 2005; Ceribelli et al., 2006; Ceribelli et al., 2008; Reed et al., 2008). Additional structural features of the CCAAT box - position, orientation, presence of multiple Transcriptional Start Sites - were previously reviewed (Dolfini et al., 2009) and will not be considered in detail here.

Large-scale motif discovery using DNA Gray code and equiprobable oligomers

Motivation: How to find motifs from genome-scale functional sequences, such as all the promoters in a genome, is a challenging problem. Word-based methods count the occurrences of oligomers to detect excessively represented ones. This approach is known to be fast and accurate compared with other methods. However, two problems have hampered the application of such methods to large-scale data. One is the computational cost necessary for clustering similar oligomers, and the other is the bias in the frequency of fixed-length oligomers, which complicates the detection of significant words.

Results: We introduce a method that uses a DNA Gray code and equiprobable oligomers, which solve the clustering problem and the oligomer bias, respectively. Our method can analyze 18 000 sequences of ~1 kbp long in 30 s. We also show that the accuracy of our method is superior to that of a leading method, especially for large-scale data and small fractions of motif-containing sequences.

Availability: The online and stand-alone versions of the application, named Hegma, are available at our website: http://www.genome.ist.i.kyoto-u.ac.jp/~ichinose/hegma/

Contact:ichinose@i.kyoto-u.ac.jp; o.gotoh@i.kyoto-u.ac.jp

Benefits of random-priming: exhaustive survey of a cDNA library from lung tissue of a SARS patient

The severe acute respiratory syndrome (SARS) leads to severe injury in the lungs with multiple factors, though the pathogenesis is still largely unclear. This paper describes the particular analyses of the transcriptome of human lung tissue that was infected by SARS‐associated coronavirus (SARS‐CoV). Random primers were used to produce ESTs from total RNA samples of the lung tissue. The result showed a high diversity of the transcripts, covering much of the human genome, including loci which do not contain protein coding sequences. 10,801 ESTs were generated and assembled into 267 contigs plus 7,659 singletons. Sequences matching to SARS‐CoV RNAs and other pneumonia‐related microbes were found. The transcripts were well classified by functional annotation. Among the 7,872 assembled sequences that were identified as from human genome, 578 non‐coding genes were revealed by BLAST search. The transcripts were mapped to the human genome with the restriction of identity = 100%, which found a candidate pool of 448 novel transcriptional loci where EST transcriptional signal was never found before. Among these, 13 loci were never reported to be transcriptional by other detection methods such as gene chips, tiling arrays, and paired‐end ditags (PETs). The result showed that random‐priming cDNA library is valid for the investigation of transcript diversity in the virus‐infected tissue. The EST data could be a useful supplemental source for SARS pathology researches. J. Med. Virol. 83:574–586, 2011. © 2011 Wiley‐Liss, Inc.

Genome-wide DNA methylation analysis reveals phytoestrogen modification of promoter methylation patterns during embryonic stem cell differentiation

BACKGROUND

Environmental challenges during development affect the fetal epigenome, but the period(s) vulnerable to epigenetic dysregulation is(are) not clear. By employing a soy phytoestrogen, genistein, that is known to alter the epigenetic states of the A(vy) allele during embryogenesis, we have explored the sensitive period for epigenetic regulation. The post-implantation period, when de novo DNA methylation actively proceeds, is amenable to in vitro analysis using a mouse embryonic stem (ES) cell differentiation system.

METHODS AND FINDINGS

Mouse ES cells were differentiated in the presence or absence of genistein, and DNA methylation patterns on day 10 were compared by microarray-based promoter methylation analysis coupled with a methylation-sensitive endonuclease (HpaII/McrBC)-dependent enrichment procedure. Moderate changes in methylation levels were observed in a subset of promoters following genistein treatment. Detailed investigation of the Ucp1 and Sytl1 promoters further revealed that genistein does not affect de novo methylation occurring between day 0 and day 4, but interferes with subsequent regulatory processes and leads to a decrease in methylation level for both promoters.

CONCLUSION

Genistein perturbed the methylation pattern of differentiated ES cells after de novo methylation. Our observations suggest that, for a subset of genes, regulation after de novo DNA methylation in the early embryo may be sensitive to genistein.

A B56gamma mutation in lung cancer disrupts the p53-dependent tumor-suppressor function of protein phosphatase 2A

Earlier studies have shown both p53-dependent and -independent tumor-suppressive functions of B56gamma-specific protein phosphatase 2A (B56gamma-PP2A). In the absence of p53, B56gamma-PP2A can inhibit cell proliferation and cell transformation by an unknown mechanism. In the presence of p53, on DNA damage, a complex including B56gamma-PP2A and p53 is formed, which leads to Thr55 dephosphorylation of p53, induction of the p53 transcriptional target p21 and inhibition of cell proliferation. In spite of its significance in inhibition of cell proliferation, no B56gamma mutations have been linked to human cancer to date. In this study, we first differentiate between the p53-dependent and -independent functions of B56gamma-PP2A by identifying a domain of the B56gamma protein required for interaction with p53. Within this region, we identify a B56gamma mutation, F395C, in lung cancer that disrupts the B56gamma-p53 interaction. More importantly, we show that F395C is unable to promote p53 Thr55 dephosphorylation, transcriptional activation of p21 and the p53-dependent tumor-suppressive function of PP2A. This finding provides a mechanistic basis for the p53-dependent and -independent functions of B56gamma-PP2A and establishes a critical link between B56gamma-PP2A p53-dependent tumor-suppressive function and tumorigenesis.

Functional analysis of a novel cis-acting regulatory region within the human ankyrin gene (ANK-1) promoter

The characterization of atypical mutations in loci associated with diseases is a powerful tool to discover novel regulatory elements. We previously identified a dinucleotide deletion in the human ankyrin-1 gene (ANK-1) promoter that underlies ankyrin-deficient hereditary spherocytosis. The presence of the deletion was associated with a decrease in promoter function both in vitro and in vivo establishing it as a causative hereditary spherocytosis mutation. The dinucleotide deletion is located in the 5' untranslated region of the ANK-1 gene and disrupts the binding of TATA binding protein and TFIID, components of the preinitiation complex. We hypothesized that the nucleotides surrounding the mutation define an uncharacterized regulatory sequence. To test this hypothesis, we generated a library of more than 16,000 ANK-1 promoters with degenerate sequence around the mutation and cloned the functional promoter sequences after cell-free transcription. We identified the wild type and three additional sequences, from which we derived a consensus. The sequences were shown to be functional in cell-free transcription, transient-transfection, and transgenic mouse assays. One sequence increased ANK-1 promoter function 5-fold, while randomly chosen sequences decreased ANK-1 promoter function. Our results demonstrate a novel functional motif in the ANK-1 promoter.

Novel neuronal proteolipid protein isoforms encoded by the human myelin proteolipid protein 1 gene

The human myelin proteolipid protein 1 gene (hPLP1), which encodes the major structural myelin proteins of the central nervous system (CNS), is classically described as expressed in the oligodendrocytes, the CNS myelinating cells. We identified two new exons in the intron 1 of the hPLP1 gene that lead to the expression of additional mRNA and protein isoforms mainly expressed in neurons instead of oligodendrocytes. Those novel neuronal PLP isoforms are detected as soon as human fetal development and their concomitant expression is specific of the human species. As classical PLP proteins, the novel protein isoforms seem to be addressed to the plasma membrane. These results suggest for the first time that PLP may have functions in humans not only in oligodendrocytes but also in neurons and could be implicated in axono-glial communication. Moreover, this neuronal expression of the hPLP1 gene might explain the neuronal dysfunctions in patients carrying hPLP1 gene mutations.

Heterogeneous patterns of gene-expression diversification in mammalian gene duplicates

Gene duplication is a major mechanism for molecular evolutionary innovation. Young gene duplicates typically exhibit elevated rates of protein evolution and, according to a number of recent studies, increased expression divergence. However, the nature of these changes is still poorly understood. To gain novel insights into the functional consequences of gene duplication, we have undertaken an in-depth analysis of a large data set of gene families containing primate- and/or rodent-specific gene duplicates. We have found a clear tendency toward an increase in protein, promoter, and expression divergence with increasing number of duplication events undergone by each gene since the human-mouse split. In addition, gene duplication is significantly associated with a reduction in expression breadth and intensity. Interestingly, it is possible to identify three main groups regarding the evolution of gene expression following gene duplication. The first group, which comprises around 25% of the families, shows patterns compatible with tissue-expression partitioning. The second and largest group, comprising 33-53% of the families, shows broad expression of one of the gene copies and reduced, overlapping, expression of the other copy or copies. This can be attributed, in most cases, to loss of expression in several tissues of one or more gene copies. Finally, a substantial number of families, 19-35%, maintain a very high level of tissue-expression overlap (>0.8) after tens of millions of years of evolution. These families may have been subject to selection for increased gene dosage.

Coding region structural heterogeneity and turnover of transcription start sites contribute to divergence in expression between duplicate genes

BACKGROUND

Gene expression divergence is one manifestation of functional differences between duplicate genes. Although rapid accumulation of expression divergence between duplicate gene copies has been observed, the driving mechanisms behind this phenomenon have not been explored in detail.

RESULTS

We examine which factors influence expression divergence between human duplicate genes, utilizing the latest genome-wide data sets. We conclude that the turnover of transcription start sites between duplicate genes occurs rapidly after gene duplication and that gene pairs with shared transcription start sites have significantly higher expression similarity than those without shared transcription start sites. Moreover, we find that most (55%) duplicate gene pairs do not retain the same coding sequence structure between the two duplicate copies and this also contributes to divergence in their expression. Furthermore, the proportion of aligned sequences in cis-regulatory regions between the two copies is positively correlated with expression similarity. Surprisingly, we find no effect of copy-specific transposable element insertions on the divergence of duplicate gene expression.

CONCLUSIONS

Our results suggest that turnover of transcription start sites, structural heterogeneity of coding sequences, and divergence of cis-regulatory regions between copies play a pivotal role in determining the expression divergence of duplicate genes.

NF-Y substitutes H2A-H2B on active cell-cycle promoters: recruitment of CoREST-KDM1 and fine-tuning of H3 methylations

The CCAAT box is a frequent promoter element, as illustrated by bioinformatic analysis, and it is bound by NF-Y, a trimer with H2A-H2B-like subunits. We developed a MNase I-based ChIP protocol on homogeneous cell populations to study cell-cycle promoters at the single nucleosome level. We analyzed histone methylations and the association of enzymatic activities. Two novel results emerged: (i) H3-H4 are present on core promoters under active conditions, with the expected cohort of ‘positive’ modifications; H2A-H2B are removed and substituted by NF-Y. Through the use of a dominant negative mutant we show that NF-Y is important for H3K36me3 deposition and for elongation, not recruitment of Pol II; (ii) H3K4 methylations are highly dynamic and H3K4me1 is a crucial positive mark. Functional siRNA inactivation and treatment with Tranylcypromine determined that KDM1 (LSD1) plays a positive role in transcription, specifically of G2/M genes. It requires CoREST, which is recruited on active promoters through direct interactions with NF-Y. These data are the first in vivo indication of a crucial interplay between core histones and ‘deviant’ histone-fold such as NF-Y, leading to fine-tuning of histone methylations.

An NF-Y-dependent switch of positive and negative histone methyl marks on CCAAT promoters

BACKGROUND

Histone tails have a plethora of different post-translational modifications, which are located differently in "open" and "closed" parts of genomes. H3K4me3/H3K79me2 and H4K20me3 are among the histone marks associated with the early establishment of active and inactive chromatin, respectively. One of the most widespread promoter elements is the CCAAT box, bound by the NF-Y trimer. Two of NF-Y subunits have an H2A-H2B-like structure.

PRINCIPAL FINDINGS

We established the causal relationship between NF-Y binding and positioning of methyl marks, by ChIP analysis of mouse and human cells infected with a dominant negative NF-YA: a parallel decrease in NF-Y binding, H3K4me3, H3K79me2 and transcription was observed in promoters that are dependent upon NF-Y. On the contrary, changes in the levels of H3K9-14ac were more subtle. Components of the H3K4 methylating MLL complex are not recruited in the absence of NF-Y. As for repressed promoters, NF-Y removal leads to a decrease in the H4K20me3 mark and deposition of H3K4me3.

CONCLUSIONS

Two relevant findings are reported: (i) NF-Y gains access to its genomic locations independently from the presence of methyl histone marks, either positive or negative; (ii) NF-Y binding has profound positive or negative consequences on the deposition of histone methyl marks. Therefore NF-Y is a fundamental switch at the heart of decision between gene activation and repression in CCAAT regulated genes.

Weak correlation between sequence conservation in promoter regions and in protein-coding regions of human-mouse orthologous gene pairs

Background

Interspecies sequence comparison is a powerful tool to extract functional or evolutionary information from the genomes of organisms. A number of studies have compared protein sequences or promoter sequences between mammals, which provided many insights into genomics. However, the correlation between protein conservation and promoter conservation remains controversial.

Results

We examined promoter conservation as well as protein conservation for 6,901 human and mouse orthologous genes, and observed a very weak correlation between them. We further investigated their relationship by decomposing it based on functional categories, and identified categories with significant tendencies. Remarkably, the 'ribosome' category showed significantly low promoter conservation, despite its high protein conservation, and the 'extracellular matrix' category showed significantly high promoter conservation, in spite of its low protein conservation.

Conclusion

Our results show the relation of gene function to protein conservation and promoter conservation, and revealed that there seem to be nonparallel components between protein and promoter sequence evolution.

Housekeeping genes tend to show reduced upstream sequence conservation

Mammalian housekeeping genes show significantly lower promoter sequence conservation, especially upstream of position -500 with respect to the transcription start site, than genes expressed in a subset of tissues.

Background

Understanding the constraints that operate in mammalian gene promoter sequences is of key importance to understand the evolution of gene regulatory networks. The level of promoter conservation varies greatly across orthologous genes, denoting differences in the strength of the evolutionary constraints. Here we test the hypothesis that the number of tissues in which a gene is expressed is related in a significant manner to the extent of promoter sequence conservation.

Results

We show that mammalian housekeeping genes, expressed in all or nearly all tissues, show significantly lower promoter sequence conservation, especially upstream of position -500 with respect to the transcription start site, than genes expressed in a subset of tissues. In addition, we evaluate the effect of gene function, CpG island content and protein evolutionary rate on promoter sequence conservation. Finally, we identify a subset of transcription factors that bind to motifs that are specifically over-represented in housekeeping gene promoters.

Conclusion

This is the first report that shows that the promoters of housekeeping genes show reduced sequence conservation with respect to genes expressed in a more tissue-restricted manner. This is likely to be related to simpler gene expression, requiring a smaller number of functional cis-regulatory motifs.

Co-regulation of Gremlin and Notch signalling in diabetic nephropathy

Diabetic nephropathy is currently the leading cause of end-stage renal disease worldwide, and occurs in approximately one third of all diabetic patients. The molecular pathogenesis of diabetic nephropathy has not been fully characterized and novel mediators and drivers of the disease are still being described. Previous data from our laboratory has identified the developmentally regulated gene Gremlin as a novel target implicated in diabetic nephropathy in vitro and in vivo. We used bioinformatic analysis to examine whether Gremlin gene sequence and structure could be used to identify other genes implicated in diabetic nephropathy. The Notch ligand Jagged1 and its downstream effector, hairy enhancer of split-1 (Hes1), were identified as genes with significant similarity to Gremlin in terms of promoter structure and predicted microRNA binding elements. This led us to discover that transforming growth factor-beta (TGFbeta1), a primary driver of cellular changes in the kidney during nephropathy, increased Gremlin, Jagged1 and Hes1 expression in human kidney epithelial cells. Elevated levels of Gremlin, Jagged1 and Hes1 were also detected in extracts from renal biopsies from diabetic nephropathy patients, but not in control living donors. In situ hybridization identified specific upregulation and co-expression of Gremlin, Jagged1 and Hes1 in the same tubuli of kidneys from diabetic nephropathy patients, but not controls. Finally, Notch pathway gene clustering showed that samples from diabetic nephropathy patients grouped together, distinct from both control living donors and patients with minimal change disease. Together, these data suggest that Notch pathway gene expression is elevated in diabetic nephropathy, co-incident with Gremlin, and may contribute to the pathogenesis of this disease.

Characterization of the equine glycogen debranching enzyme gene (AGL): Genomic and cDNA structure, localization, polymorphism and expression

Glycogen debranching enzyme (AGL) is a multifunctional enzyme acting in the glycogen degradation pathway. In humans, the AGL activity deficiency causes a type III glycogen storage disease (Cori-Forbes disease). One particularity of AGL gene expression lies in the multiple alternative splicing in its 5' region. The AGL gene was localized on ECA5q14-q15. The sequence of the equine cDNA was determined to be 7.5 kb in length with an open reading frame of 4602 bp. The gene is 69 kb long and contains 35 exons. The equine AGL gene has an ubiquitous expression and presents five tissue-dependent cDNA variants arising from alternative splicing of the first exons. The equine skeletal muscle and heart contain four out of six variants previously described in humans and the equine liver express three of these four human variants. We identified a new alternative splicing variant expressed in equine skeletal and heart muscles. All these mRNA variants most probably encode only two different protein isoforms of 1533 and 1377 amino-acids. Four SNPs were detected in the mRNA. The equine in silico promoter sequence reveals a structure similar to those of other mammalian species. The disposition of the transcription factor biding sites does not correlate to the transcription start sites of tissue-specific variants.

Distinct class of putative "non-conserved" promoters in humans: comparative studies of alternative promoters of human and mouse genes

Although recent studies have revealed that the majority of human genes are subject to regulation of alternative promoters, the biological relevance of this phenomenon remains unclear. We have also demonstrated that roughly half of the human RefSeq genes examined contain putative alternative promoters (PAPs). Here we report large-scale comparative studies of PAPs between human and mouse counterpart genes. Detailed sequence comparison of the 17,245 putative promoter regions (PPRs) in 5463 PAP-containing human genes revealed that PPRs in only a minor fraction of genes (807 genes) showed clear evolutionary conservation as one or more pairs. Also, we found that there were substantial qualitative differences between conserved and non-conserved PPRs, with the latter class being AT-rich PPRs of relative minor usage, enriched in repetitive elements and sometimes producing transcripts that encode small or no proteins. Systematic luciferase assays of these PPRs revealed that both classes of PPRs did have promoter activity, but that their strength ranges were significantly different. Furthermore, we demonstrate that these characteristic features of the non-conserved PPRs are shared with the PPRs of previously discovered putative non-protein coding transcripts. Taken together, our data suggest that there are two distinct classes of promoters in humans, with the latter class of promoters emerging frequently during evolution.

Anecdotes, data and regulatory modules

Beginning in the late 1980s, Eric Davidson's group at Cal Tech developed a modularity hypothesis of developmental gene regulation, showing that in an expanding number of cases, particular aspects of development were governed by compact 'modules' of transcription factor binding sites (TFBSs), and that these modules were separable, complex and interconnected. Davidson made no attempt to further generalize the hypothesis, but others took up the idea, transported it out of development and extended it to a general rule of clustering. Despite such misbegotten origins, the 'extended' modularity hypothesis--that TFBSs in general tend to come in compact clusters--has been highly productive, yet it has never been challenged with a large, diverse and unbiased dataset to see how universal it actually is. The aim of the present paper is to do so. Applying human-mouse-rat phylogenetic footprinting to neighbourhoods of a diverse set of TFBSs, including both developmental and non-developmental signals, we find that the extended hypothesis holds in at least 93.5% of cases. Based on this particular sample, we found a mean module length of 609 nucleotides containing, on an average, 24.5 presumptive regulatory signals of length greater than 5 and averaging 8.5 nucleotides each.

A scaffold for the Chinese hamster genome

Chinese hamster ovary (CHO) cells are a prevalent tool in biological research and are among the most widely used host cell lines for production of recombinant therapeutic proteins. While research in other organisms has been revolutionized through the development of DNA sequence-based tools, the lack of comparable genomic resources for the Chinese hamster has impeded similar work in CHO cell lines. A comparative genomics approach, based upon the completely sequenced mouse genome, can facilitate genomic work in this important organism. Using chromosome synteny to define regions of conserved linkage between Chinese hamster and mouse chromosomes, a working scaffold for the Chinese hamster genome has been developed. Mapping CHO and Chinese hamster sequences to the mouse genome creates direct access to relevant information in public databases. Additionally, mapping gene expression data onto a chromosome scaffold affords the ability to interpret information in a genomic context, potentially revealing important structural and regulatory features in the Chinese hamster genome. Further development of this genomic scaffold will provide opportunities to use biomolecular tools for research in CHO cell lines today and will be an asset to future efforts to sequence the Chinese hamster genome.

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.

Bioinformatics tools for modeling transcription factor target genes and epigenetic changes

The combinatorial control of gene regulatory switches involves both transcription factor (TF) complexes and associated epigenetic modifications to the chromatin template. The novel high-throughput technologies, such as Chromatin ImmunoPrecipitation ChIP-chip, have enabled genome-wide in vivo identification of TF target regulatory regions and related epigenetic modifications, which led to the view of highly dynamic TF-DNA interactions in activated or repressed promoters. Consequently, modeling and elucidating the combinatorial interaction of TFs and corresponding cis-regulatory modules in target promoters is of paramount interest. An estimated 5% of the genes in mammalian genomes code for TF proteins, and computational modeling of cis-regulatory logic would rapidly increase the pace of experimental confirmation of TF target promoters at the bench. The purpose of this chapter is to discuss the use of different bioinformatics tools for predicting the target genes of TFs of interest in mammalian genomes, and the application of these methods in the analysis of ChIP-chip experimental data. The author describes most commonly used databases and prediction programs that are available on the World Wide Web and demonstrate the use of some of these programs by an example. A list of these programs is provided along with their web Uniform Resource Locator (URLs) and guidelines for successful application are suggested.

A computational genomics approach to identify cis-regulatory modules from chromatin immunoprecipitation microarray data--a case study using E2F1

Advances in high-throughput technologies, such as ChIP-chip, and the completion of human and mouse genomic sequences now allow analysis of the mechanisms of gene regulation on a systems level. In this study, we have developed a computational genomics approach (termed ChIPModules), which begins with experimentally determined binding sites and integrates positional weight matrices constructed from transcription factor binding sites, a comparative genomics approach, and statistical learning methods to identify transcriptional regulatory modules. We began with E2F1 binding site information obtained from ChIP-chip analyses of ENCODE regions, from both HeLa and MCF7 cells. Our approach not only distinguished targets from nontargets with a high specificity, but it also identified five regulatory modules for E2F1. One of the identified modules predicted a colocalization of E2F1 and AP-2alpha on a set of target promoters with an intersite distance of <270 bp. We tested this prediction using ChIP-chip assays with arrays containing approximately 14,000 human promoters. We found that both E2F1 and AP-2alpha bind within the predicted distance to a large number of human promoters, demonstrating the strength of our sequence-based, unbiased, and universal protocol. Finally, we have used our ChIPModules approach to develop a database that includes thousands of computationally identified and/or experimentally verified E2F1 target promoters.

Genomic structure, alternative maturation and tissue expression of the human BBOX1 gene

Gamma-butyrobetaine hydroxylase (BBOX1) is the enzyme responsible for the biosynthesis of l-carnitine, a key molecule of fatty acid metabolism. This cytosolic dimeric protein belongs to the dioxygenase family. In human, enzyme activity has been detected in kidney, liver and brain. The human gene encoding gamma-butyrobetaine hydroxylase is located on chromosome 11. Although the protein structure and activity have been extensively described, little information is available concerning BBOX1 structure and expression. In this study, the organization of the human gene was determined. The structure and functions of the 5'- and 3'-untranslated regions of the human BBOX1 mRNA were characterized in kidney, liver and brain. Our experiments revealed that the transcription initiation of the human BBOX1 gene might occur at 3 different exons, and that the expression level of each type of transcript is organ-specific. We showed that the use of 3 different promoters is responsible for the 5'-end heterogeneity. Investigations on BBOX1 mRNA maturation highlighted an alternative polyadenylation mechanism that generates two 3'-untranslated regions differing by their length. This alternative polyadenylation exhibited a tissue specificity.

Comparative genomic analysis of transcription regulation elements involved in human map kinase G-protein coupling pathway

The identification of cis-elements (motifs) in the regulatory regions of higher eukaryotes is an important and challenging problem in computational biology. Eukaryotic transcriptional regulatory mechanisms pose several difficulties for promoter analysis: including a high variance in the motif locations, frequently large divergence from motif consensus patterns, and a large amount of repetitive elements (confusing to many motif finding procedures). One promising approach to this difficult problem involves cross-species comparison. In this work we analyzed the full-length regulatory regions of genes involved in the G-protein coupling MAP kinase pathway and compared the results with ribosomal genes using human, mouse and rat genomic data. We found 19 high likely transcription factors (TFs) candidates for MAPK and 12 TFs for the ribosomal dataset. In the case of the MAPK dataset, regulatory regions of genes functionally grouped as receptors and MAP-core genes were found mostly highly conserved across the three species.

Analysis of promoter regions of co-expressed genes identified by microarray analysis

Background

The use of global gene expression profiling to identify sets of genes with similar expression patterns is rapidly becoming a widespread approach for understanding biological processes. A logical and systematic approach to study co-expressed genes is to analyze their promoter sequences to identify transcription factors that may be involved in establishing specific profiles and that may be experimentally investigated.

Results

We introduce promoter clustering i.e. grouping of promoters with respect to their high scoring motif content, and show that this approach greatly enhances the identification of common and significant transcription factor binding sites (TFBS) in co-expressed genes. We apply this method to two different dataset, one consisting of micro array data from 108 leukemias (AMLs) and a second from a time series experiment, and show that biologically relevant promoter patterns may be obtained using phylogenetic foot-printing methodology. In addition, we also found that 15% of the analyzed promoter regions contained transcription factors start sites for additional genes transcribed in the opposite direction.

Conclusion

Promoter clustering based on global promoter features greatly improve the identification of shared TFBS in co-expressed genes. We believe that the outlined approach may be a useful first step to identify transcription factors that contribute to specific features of gene expression profiles.

dbQSNP: a database of SNPs in human promoter regions with allele frequency information determined by single-strand conformation polymorphism-based methods

We present a database, dbQSNP (http://qsnp.gen.kyushu-u.ac.jp/), that provides sequence and allele frequency information for single-nucleotide polymorphisms (SNPs) located in the promoter regions of human genes, which were defined by the 5' ends of full-length cDNA clones. We searched for the SNPs in these regions by sequencing or single-strand conformation polymorphism (SSCP) analysis. The allele frequencies of the identified SNPs in two ethnic groups were quantified by SSCP analyses of pooled DNA samples. The accuracy of our estimation is supported by strong correlations between the frequencies in our data and those in other databases for the same ethnic groups. The frequencies vary considerably between the two ethnic groups studied, suggesting the need for population-based collections and allele frequency determination of SNPs, in, e.g., association studies of diseases. We show profiles of SNP densities that are characteristic of transcription start site regions. A fraction of the SNPs revealed a significantly different allele frequency between the groups, suggesting differential selection of the genes involved.

Genome-wide analysis of core promoter elements from conserved human and mouse orthologous pairs

Background

The canonical core promoter elements consist of the TATA box, initiator (Inr), downstream core promoter element (DPE), TFIIB recognition element (BRE) and the newly-discovered motif 10 element (MTE). The motifs for these core promoter elements are highly degenerate, which tends to lead to a high false discovery rate when attempting to detect them in promoter sequences.

Results

In this study, we have performed the first analysis of these core promoter elements in orthologous mouse and human promoters with experimentally-supported transcription start sites. We have identified these various elements using a combination of positional weight matrices (PWMs) and the degree of conservation of orthologous mouse and human sequences – a procedure that significantly reduces the false positive rate of motif discovery. Our analysis of 9,010 orthologous mouse-human promoter pairs revealed two combinations of three-way synergistic effects, TATA-Inr-MTE and BRE-Inr-MTE. The former has previously been putatively identified in human, but the latter represents a novel synergistic relationship.

Conclusion

Our results demonstrate that DNA sequence conservation can greatly improve the identification of functional core promoter elements in the human genome. The data also underscores the importance of synergistic occurrence of two or more core promoter elements. Furthermore, the sequence data and results presented here can help build better computational models for predicting the transcription start sites in the promoter regions, which remains one of the most challenging problems.

Genome-wide prediction of mammalian enhancers based on analysis of transcription-factor binding affinity

Understanding the regulation of human gene expression requires knowledge of the "second genetic code," which consists of the binding specificities of transcription factors (TFs) and the combinatorial code by which TF binding sites are assembled to form tissue-specific enhancer elements. Using a novel high-throughput method, we determined the DNA binding specificities of GLIs 1-3, Tcf4, and c-Ets1, which mediate transcriptional responses to the Hedgehog (Hh), Wnt, and Ras/MAPK signaling pathways. To identify mammalian enhancer elements regulated by these pathways on a genomic scale, we developed a computational tool, enhancer element locator (EEL). We show that EEL can be used to identify Hh and Wnt target genes and to predict activated TFs based on changes in gene expression. Predictions validated in transgenic mouse embryos revealed the presence of multiple tissue-specific enhancers in mouse c-Myc and N-Myc genes, which has implications for organ-specific growth control and tumor-type specificity of oncogenes.