Mapping of transcription start sites of human retina expressed genes

Polyglutamine-expanded ATXN7 alters a specific epigenetic signature underlying photoreceptor identity gene expression in SCA7 mouse retinopathy

BACKGROUND

Spinocerebellar ataxia type 7 (SCA7) is a neurodegenerative disorder that primarily affects the cerebellum and retina. SCA7 is caused by a polyglutamine expansion in the ATXN7 protein, a subunit of the transcriptional coactivator SAGA that acetylates histone H3 to deposit narrow H3K9ac mark at DNA regulatory elements of active genes. Defective histone acetylation has been presented as a possible cause for gene deregulation in SCA7 mouse models. However, the topography of acetylation defects at the whole genome level and its relationship to changes in gene expression remain to be determined.

METHODS

We performed deep RNA-sequencing and chromatin immunoprecipitation coupled to high-throughput sequencing to examine the genome-wide correlation between gene deregulation and alteration of the active transcription marks, e.g. SAGA-related H3K9ac, CBP-related H3K27ac and RNA polymerase II (RNAPII), in a SCA7 mouse retinopathy model.

RESULTS

Our analyses revealed that active transcription marks are reduced at most gene promoters in SCA7 retina, while a limited number of genes show changes in expression. We found that SCA7 retinopathy is caused by preferential downregulation of hundreds of highly expressed genes that define morphological and physiological identities of mature photoreceptors. We further uncovered that these photoreceptor genes harbor unusually broad H3K9ac profiles spanning the entire gene bodies and have a low RNAPII pausing. This broad H3K9ac signature co-occurs with other features that delineate superenhancers, including broad H3K27ac, binding sites for photoreceptor specific transcription factors and expression of enhancer-related non-coding RNAs (eRNAs). In SCA7 retina, downregulated photoreceptor genes show decreased H3K9 and H3K27 acetylation and eRNA expression as well as increased RNAPII pausing, suggesting that superenhancer-related features are altered.

CONCLUSIONS

Our study thus provides evidence that distinctive epigenetic configurations underlying high expression of cell-type specific genes are preferentially impaired in SCA7, resulting in a defect in the maintenance of identity features of mature photoreceptors. Our results also suggest that continuous SAGA-driven acetylation plays a role in preserving post-mitotic neuronal identity.

Investigation of promoter regions, motifs, and CpG islands in the regulation of gene expression in Trametes hirsuta strain 072

Background

In silico analysis of transcription start sites, promoter regions, transcription factors and their binding sites, and CpG islands for the Trametes hirsuta strain 072 genome were performed to understand the regulation mechanisms of gene expression and its genetic variations in the genomes. Therefore, a computational survey was carried out for the Trametes hirsuta strain 072 genome with the open reading frames from the National Center for Biotechnology Information database. Seventeen functional sequences were used to analyze promoter regions and their regulatory elements.

Result

The present study revealed that 94% of Trametes hirsuta strain 072 genes contained more than two TSSs. Among these identified TSSs, a TSS with the highest predictive score was considered to determine a promoter region of the genes. Moreover, a total of five common candidate motifs such as MotI, MotII, MotIII, MotIV, and MotV were identified. Among these motifs, motif IV was investigated as the common promoter motif for 41.17% of genes that serve as binding sites for transcription factors (TFs) involved in the expression regulation of Trametes hirsuta strain 072 genes. Motif IV was also compared to registered motifs in publically available databases to see if they are similar to known regulatory motifs for TF using TOMTOM web server. Hence, it was revealed that MotIV might serve as the binding site mainly for the leucine zipper TF gene family to regulate a gene expression of Trametes hirsuta strain 072. Regarding CpG island determination, it was concluded that there is no CpG island in both promoter and gene body regions of the Trametes hirsuta strain 072 genome.

Conclusions

This study provides a better insight into further molecular characterization which aimed to efficiently exploit a white rot fungus, Trametes hirsuta strain 072, for several biotechnological applications aimed to revitalize a severely contaminated environment.

In silico analysis of the promoter region of olfactory receptors in cattle (Bos indicus) to understand its gene regulation

Identifications of transcription start sites (TSSs) and promoter regions are first step to understand the regulation mechanisms of gene expression and association with genetic variations in the regions. This analysis was conducted with the objectives to identify TSSs, determine the promoter regions, identify common candidate motifs and transcription factors (TFs), and search for CpG islands (CGIs) in cattle olfactory receptors (ORs) genes promoter regions. In the analysis, TSSs of cattle olfactory genes were first identified. The locations for 60% of the TSSs were below -500 bp relative to the start codon and five candidate motifs (MOR1, MOR2, MOR3, MOR4, and MOR5) were identified that are shared by at least 50% of the cattle ORs promoter input sequences from both strands. Among the five candidate motifs, MOR4 was revealed as the common promoter motif for 85.71% of cattle ORs genes that serves as binding sites for TFs involved in the expression regulation of these genes. MOR4 was also compared to registered motifs in publically available databases to see if they are similar to known regulatory motifs for TF by using the TOMTOM web application. Hence, it was revealed that MOR4 may serve as the binding site mainly for the Zinc finger (ZNF) TF gene family to regulate expression of cattle ORs genes. Further gene ontology analysis for MOR4 demonstrated ORs belong to the G-protein-coupled receptor superfamily and MOR4 tend to be located near the genes involved in the detection of chemical stimulus involved in sensory perception and in innate immune responses such as cytokine-mediated signaling. In silico digestion of cattle OR sequences was performed using restriction enzyme MspI. CGIs from OR10K1 and OR2L13 gene was found. In the present analysis, the poor CGIs observed might suggest their gene expression regulation pattern is in tissue specific manner.

Analysis of Pig Vomeronasal Receptor Type 1 (V1R) Promoter Region Reveals a Common Promoter Motif but Poor CpG Islands

Promoters are, generally, located immediately upstream of a transcription start site (TSS) and have a variety of regulatory motifs, such as transcription factors (TFs) and CpG islands (CGIs), that participate in the regulation of gene expression. Here analysis of the promoter region for pig vomeronasal receptor type 1 (V1R) was described. In the analysis, TSSs for pig V1R genes was first identified and five motifs (MV1, MV2, MV3, MV4, and MV5) were found that are shared by at least 50% of the pig V1R promoter input sequences from both strands. Among the five motifs, MV2 was identified as a common promoter motif shared by all (100%) pig V1R promoters. For further analysis, to better characterize and get deeper biological insight associated with MV2, TOMTOM web application was used. MV2 was compared to the known motif databases (such as JASPAR) to see if they are similar to a known regulatory motif (transcription factor). Hence, it was revealed that MV2 serves as the binding site mainly for the BetaBetaAlpha-zinc finger (BTB-ZF) transcription factor gene family to regulate expression of pig V1R genes. Moreover, it was shown that pig V1R promoters are CpG poor, suggesting that their gene expression regulation pattern is in tissue specific manner.

The c.-360 T>C mutation affects PGAM2 transcription activity and is linked with the water holding capacity of the longissimus lumborum muscle in pigs

The phosphoglycerate mutase 2 (PGAM2) gene encodes a key enzyme in the glycolytic process. This study examined a functional mutation in the PGAM2 gene and evaluated its relationship with water holding capacity (WHC). RT-qPCR analysis showed the PGAM2 mRNA level was significantly higher in the low-WHC group than in the high-WHC group (P<0.05). The c.-360 T>C mutation was identified through sequencing and found to have opposite allele distributions in the two groups. The allele was further genotyped in 170 Duroc×Large White×Yorkshire crossbred pigs using allele-specific PCR. The CC genotype was associated with lower WHC and higher PGAM2 mRNA levels, whereas the TT genotype corresponded to a higher WHC and lower PGAM2 mRNA levels (P<0.05). A luciferase activity assay also showed that the CC-genotype promoter had higher activity than the TT-genotype promoter (P<0.05). In conclusion, we discovered the c.-360 T>C mutation in the PGAM2 gene, which is a promising marker for improving pork WHC.

Where does transcription start? 5'-RACE adapted to next-generation sequencing

The variability and complexity of the transcription initiation process was examined by adapting RNA ligase-mediated rapid amplification of 5' cDNA ends (5'-RACE) to Next-Generation Sequencing (NGS). We oligo-labelled 5'-m(7)G-capped mRNA from two genes, the simple mono-exonic Beta-2-Adrenoceptor (ADRB2R)and the complex multi-exonic Glucocorticoid Receptor (GR, NR3C1), and detected a variability in TSS location that has received little attention up to now. Transcription was not initiated at a fixed TSS, but from loci of 4 to 10 adjacent nucleotides. Individual TSSs had frequencies from <0.001% to 38.5% of the total gene-specific 5' m(7)G-capped transcripts. ADRB2R used a single locus consisting of 4 adjacent TSSs. Unstimulated, the GR used a total of 358 TSSs distributed throughout 38 loci, that were principally in the 5' UTRs and were spliced using established donor and acceptor sites. Complete demethylation of the epigenetically sensitive GR promoter with 5-azacytidine induced one new locus and 127 TSSs, 12 of which were unique. We induced GR transcription with dexamethasone and Interferon-γ, adding one new locus and 185 additional TSSs distributed throughout the promoter region. In-vitro the TSS microvariability regulated mRNA translation efficiency and the relative abundance of the different GRN-terminal protein isoform levels.

Melanocortin 3 receptor has a 5' exon that directs translation of apically localized protein from the second in-frame ATG

Melanocortin-3 receptor (MC3R) is a canonical MSH receptor that plays an essential role in energy homeostasis. Variants in MC3R have been implicated in obesity in humans and mice. However, interpretation of the functional consequences of these variants is challenging because the translational start site of MC3R is unclear. Using 5' rapid amplification of cDNA ends, we discovered a novel upstream exon that extends the length of the 5' untranslated region (UTR) in MC3R without changing the open-reading frame. The full-length 5' UTR directs utilization of an evolutionarily conserved second in-frame ATG as the primary translation start site. MC3R synthesized from the second ATG is localized to apical membranes of polarized Madin-Darby canine kidney cells, consistent with its function as a cell surface mediator of melanocortin signaling. Expression of MC3R causes relocalization of melanocortin receptor accessory protein 2, an accessory factor for melanocortin-2 receptor, to the apical membrane, coincident with the location of MC3R. In contrast, protein synthesized from MC3R cDNAs lacking the 5' UTR displayed diffuse cytosolic distribution and has no effect on the distribution of melanocortin receptor accessory protein 2. Our findings demonstrate that a previously unannotated 5' exon directs translation of MC3R protein that localizes to apical membranes of polarized cells. Together, our work provides insight on the structure of human MC3R and reveals a new pathway for regulation of energy metabolism.

Lipolysis-stimulated lipoprotein receptor: a novel membrane protein of tricellular tight junctions

Tricellular tight junctions (tTJs) are specialized structural variants of tight junctions that restrict the free diffusion of solutes at the extracellular space of tricellular contacts. Their presence at cell corners, situated in the angles between three adjacent epithelial cells, was identified early by electron microscopy, but despite their potential importance, tTJs have been generally ignored in epithelial cell biology. Tricellulin was the first molecular component of tTJs shown to be involved in their formation and in epithelial barrier function. However, the precise molecular organization and function of tTJs are still largely unknown. Recently, we identified the lipolysis-stimulated lipoprotein receptor (LSR) as a tTJ-associated membrane protein. LSR recruits tricellulin to tTJs, suggesting that the LSR-tricellulin system plays a key role in tTJ formation. In this paper, we summarize the identification and characterization of LSR as a molecular component of tTJs.

Mitotic recombination of chromosome arm 17q as a cause of loss of heterozygosity of NF1 in neurofibromatosis type 1-associated glomus tumors

Neurofibromatosis type 1 (NF1) is a common, autosomal dominant, tumor-predisposition syndrome that arises secondary to mutations in NF1. Glomus tumors are painful benign tumors that originate from the glomus body in the fingers and toes due to biallelic inactivation of NF1. We karyotyped cultures from four previously reported and one new glomus tumor and hybridized tumor (and matching germline) DNA on Illumina HumanOmni1-Quad SNP arrays (≈ 1 × 10(6) SNPs). Two tumors displayed evidence of copy-neutral loss of heterozygosity of chromosome arm 17q not observed in the germline sample, consistent with a mitotic recombination event. One of these two tumors, NF1-G12, featured extreme polyploidy (near-tetraploidy, near-hexaploidy, or near-septaploidy) across all chromosomes. In the remaining four tumors, there were few cytogenetic abnormalities observed, and copy-number analysis was consistent with diploidy in all chromosomes. This is the first study of glomus tumors cytogenetics, to our knowledge, and the first to report biallelic inactivation of NF1 secondary to mitotic recombination of chromosome arm 17q in multiple NF1-associated glomus tumors. We have observed mitotic recombination in 22% of molecularly characterized NF1-associated glomus tumors, suggesting that it is a not uncommon mechanism in the reduction to homozygosity of the NF1 germline mutation in these tumors. In tumor NF1-G12, we hypothesize that mitotic recombination also "unmasked" (reduced to homozygosity) a hypomorphic germline allele in a gene on chromosome arm 17q associated with chromosomal instability, resulting in the extreme polyploidy.

Isoform diversity and its importance for axon regeneration

Axon regeneration is a fundamental problem facing neuroscientists and clinicians. Failure of axon regeneration is caused by both extrinsic and intrinsic mechanisms. New techniques to examine gene expression such as Next Generation Sequencing of the Transcriptome (RNA-Seq) drastically increase our knowledge of both gene expression complexity (RNA isoforms) and gene expression regulation. By utilizing RNA-Seq, gene expression can now be defined at the level of isoforms, an essential step for understanding the mechanisms governing cell identity, growth and ultimately cellular responses to injury and disease.

Variation in the IGF2 gene promoter region is associated with intramuscular fat content in porcine skeletal muscle

Intramuscular fat (IMF) and subcutaneous fat (back fat-BF) are two of the major fat depots in livestock. A QTN located in the insulin-like growth factor 2 gene (IGF2) has been associated with a desirable reduction in BF depth in pigs. Given that the lipid metabolism of intramuscular adipocytes differs from that of subcutaneous fat adipocytes, this study aimed to search for genetic variation in the IGF2 gene that may be associated with IMF, as well as BF, in diverse pig breeds. Four proximal promoter regions of the IGF2 gene were characterised and the association of IGF2 genetic variation with IMF and BF was assessed. Six promoter SNPs were identified in four promoter regions (P1-P4; sequence coverage 945, 866, 784 and 864 bp, respectively) in phenotypically diverse F1 cross populations. Three promoter SNPs were subsequently genotyped in three pure breeds (Pietrain = 98, Duroc = 99 and Large White = 98). All three SNPs were >95% monomorphic in the Pietrain and Duroc breeds but minor alleles were at moderate frequencies in the Large White breed. These SNPs were linked and one was located in a putative transcription factor binding site. Five haplotypes were inferred and three combined diplotypes tested for association with IMF and BF in the Large White. As expected haplotype 1 (likely in LD with the beneficial QTN allele) was superior for BF level. In contrast, the heterozygote diplotype of the most common haplotypes (1 and 2) was associated with higher IMF and marbling scores compared to either homozygote. Gene expression analysis of divergent animals showed that IGF2 was 1.89 fold up-regulated in muscle with higher compared to lower IMF content. These findings suggest that genetic variation in the promoter region of the IGF2 gene is associated with IMF content in porcine skeletal muscle and that greater expression of the IGF2 gene is associated with higher IMF content.

Regulatory polymorphisms in the bovine Ankyrin 1 gene promoter are associated with tenderness and intramuscular fat content

Background

Recent QTL and gene expression studies have highlighted ankyrins as positional and functional candidate genes for meat quality. Our objective was to characterise the promoter region of the bovine ankyrin 1 gene and to test polymorphisms for association with sensory and technological meat quality measures.

Results

Seven novel promoter SNPs were identified in a 1.11 kb region of the ankyrin 1 promoter in Angus, Charolais and Limousin bulls (n = 15 per breed) as well as 141 crossbred beef animals for which meat quality data was available. Eighteen haplotypes were inferred with significant breed variation in haplotype frequencies. The five most frequent SNPs and the four most frequent haplotypes were subsequently tested for association with sensory and technological measures of meat quality in the crossbred population. SNP1, SNP3 and SNP4 (which were subsequently designated regulatory SNPs) and SNP5 were associated with traits that contribute to sensorial and technological measurements of tenderness and texture; Haplotype 1 and haplotype 4 were oppositely correlated with traits contributing to tenderness (P < 0.05). While no single SNP was associated with intramuscular fat (IMF), a clear association with increased IMF and juiciness was observed for haplotype 2.

Conclusion

The conclusion from this study is that alleles defining haplotypes 2 and 4 could usefully contribute to marker SNP panels used to select individuals with improved IMF/juiciness or tenderness in a genome-assisted selection framework.

Mapping of the transcription start site (TSS) and identification of SNPs in the bovine neuropeptide Y (NPY) gene

BACKGROUND

Neuropeptide Y is a key neurotransmitter of the central nervous system which plays a vital role in the feed energy homeostasis in mammals. Mutations in the regulatory and coding regions of the bovine NPY gene can potentially affect the neuronal regulation of appetite and feeding behaviour in cattle. The objectives of this experiment were to: a) fully characterize the bovine NPY gene transcript and b) identify the SNP diversity in both coding and non-coding regions of the NPY gene in a panel of Bos taurus and B. indicus cattle.

RESULTS

Bovine NPY gene consists of four exons (99, 188, 81 and 195 nucleotides) and three introns. The promoter region of the NPY gene consists of TATA and GC boxes which are separated from the transcription start site (TSS) by 29 and ~100 nt, respectively. Analyses of the tissue specific expression of the bovine NPY gene revealed the presence of highly abundant NPY gene transcripts in the arcuate nucleus, cerebral and cerebellar regions of the bovine brain. We identified a total of 59 SNPs in the 8.4 kb region of the bovine NPY gene. Seven out of nine total SNPs in the promoter region affect binding of putative transcription factors. A high level of nucleotide diversity was evident in the promoter regions (2.84 x 10(-3)) compared to the exonic (1.44 x 10(-3)), intronic (1.30 x 10(-3)) and 3' untranslated (1.26 x 10(-3)) regions.

CONCLUSION

The SNPs identified in different regions of bovine NPY gene may serve as a basis for understanding the regulation of the expression of the bovine NPY gene under a variety of physiological conditions and identification of genotypes with high feed energy conversion efficiency.

Bioinformatic identification of novel putative photoreceptor specific cis-elements

Background

Cell specific gene expression is largely regulated by different combinations of transcription factors that bind cis-elements in the upstream promoter sequence. However, experimental detection of cis-elements is difficult, expensive, and time-consuming. This provides a motivation for developing bioinformatic methods to identify cis-elements that could prioritize future experimental studies. Here, we use motif discovery algorithms to predict transcription factor binding sites involved in regulating the differences between murine rod and cone photoreceptor populations.

Results

To identify highly conserved motifs enriched in promoters that drive expression in either rod or cone photoreceptors, we assembled a set of murine rod-specific, cone-specific, and non-photoreceptor background promoter sequences. These sets were used as input to a newly devised motif discovery algorithm called Iterative Alignment/Modular Motif Selection (IAMMS). Using IAMMS, we predicted 34 motifs that may contribute to rod-specific (19 motifs) or cone-specific (15 motifs) expression patterns. Of these, 16 rod- and 12 cone-specific motifs were found in clusters near the transcription start site. New findings include the observation that cone promoters tend to contain TATA boxes, while rod promoters tend to be TATA-less (exempting Rho and Cnga1). Additionally, we identify putative sites for IL-6 effectors (in rods) and RXR family members (in cones) that can explain experimental data showing changes to cell-fate by activating these signaling pathways during rod/cone development. Two of the predicted motifs (NRE and ROP2) have been confirmed experimentally to be involved in cell-specific expression patterns. We provide a full database of predictions as additional data that may contain further valuable information. IAMMS predictions are compared with existing motif discovery algorithms, DME and BioProspector. We find that over 60% of IAMMS predictions are confirmed by at least one other motif discovery algorithm.

Conclusion

We predict novel, putative cis-elements enriched in the promoter of rod-specific or cone-specific genes. These are candidate binding sites for transcription factors involved in maintaining functional differences between rod and cone photoreceptor populations.