Abundance, arrangement, and function of sequence motifs in the chicken promoters

Chicken Erythrocyte: Epigenomic Regulation of Gene Activity

The chicken genome is one-third the size of the human genome and has a similarity of sixty percent when it comes to gene content. Harboring similar genome sequences, chickens' gene arrangement is closer to the human genomic organization than it is to rodents. Chickens have been used as model organisms to study evolution, epigenome, and diseases. The chicken nucleated erythrocyte's physiological function is to carry oxygen to the tissues and remove carbon dioxide. The erythrocyte also supports the innate immune response in protecting the chicken from pathogens. Among the highly studied aspects in the field of epigenetics are modifications of DNA, histones, and their variants. In understanding the organization of transcriptionally active chromatin, studies on the chicken nucleated erythrocyte have been important. Through the application of a variety of epigenomic approaches, we and others have determined the chromatin structure of expressed/poised genes involved in the physiological functions of the erythrocyte. As the chicken erythrocyte has a nucleus and is readily isolated from the animal, the chicken erythrocyte epigenome has been studied as a biomarker of an animal's long-term exposure to stress. In this review, epigenomic features that allow erythroid gene expression in a highly repressive chromatin background are presented.

The chicken pan-genome reveals gene content variation and a promoter region deletion in IGF2BP1 affecting body size

Domestication and breeding have reshaped the genomic architecture of chicken, but the retention and loss of genomic elements during these evolutionary processes remain unclear. We present the first chicken pan-genome constructed using 664 individuals, which identified an additional ∼66.5 Mb sequences that are absent from the reference genome (GRCg6a). The constructed pan-genome encoded 20,491 predicated protein-coding genes, of which higher expression level are observed in conserved genes relative to dispensable genes. Presence/absence variation (PAV) analyses demonstrated that gene PAV in chicken was shaped by selection, genetic drift, and hybridization. PAV-based GWAS identified numerous candidate mutations related to growth, carcass composition, meat quality, or physiological traits. Among them, a deletion in the promoter region of IGF2BP1 affecting chicken body size is reported, which is supported by functional studies and extra samples. This is the first time to report the causal variant of chicken body size QTL located at chromosome 27 which was repeatedly reported. Therefore, the chicken pan-genome is a useful resource for biological discovery and breeding. It improves our understanding of chicken genome diversity and provides materials to unveil the evolution history of chicken domestication.

Mathematical Algorithm for Identification of Eukaryotic Promoter Sequences

Identification of promoter sequences in the eukaryotic genome, by computer methods, is an important task of bioinformatics. However, this problem has not been solved since the best algorithms have a false positive probability of 10−3–10−4 per nucleotide. As a result of full genome analysis, there may be more false positives than annotated gene promoters. The probability of a false positive should be reduced to 10−6–10−8 to reduce the number of false positives and increase the reliability of the prediction. The method for multi alignment of the promoter sequences was developed. Then, mathematical methods were developed for calculation of the statistically important classes of the promoter sequences. Five promoter classes, from the rice genome, were created. We developed promoter classes to search for potential promoter sequences in the rice genome with a false positive number less than 10−8 per nucleotide. Five classes of promoter sequences contain 1740, 222, 199, 167 and 130 promoters, respectively. A total of 145,277 potential promoter sequences (PPSs) were identified. Of these, 18,563 are promoters of known genes, 87,233 PPSs intersect with transposable elements, and 37,390 PPSs were found in previously unannotated sequences. The number of false positives for a randomly mixed rice genome is less than 10−8 per nucleotide. The method developed for detecting PPSs was compared with some previously used approaches. The developed mathematical method can be used to search for genes, transposable elements, and transcript start sites in eukaryotic genomes.

The chicken model organism for epigenomic research

The chicken model organism has advanced the areas of developmental biology, virology, immunology, oncology, epigenetic regulation of gene expression, conservation biology, and genomics of domestication. Further, the chicken model organism has aided in our understanding of human disease. Through the recent advances in high-throughput sequencing and bioinformatic tools, researchers have successfully identified sequences in the chicken genome that have human orthologs, improving mammalian genome annotation. In this review, we highlight the importance of chicken as an animal model in basic and pre-clinical research. We will present the importance of chicken in poultry epigenetics and in genomic studies that trace back to their ancestor, the last link between human and chicken in the tree of life. There are still many genes of unknown function in the chicken genome yet to be characterized. By taking advantage of recent sequencing technologies, it is possible to gain further insight into the chicken epigenome.

Atypical chromatin structure of immune-related genes expressed in chicken erythrocytes

The major biological role of red blood cells is to carry oxygen to the tissues in the body. However, another role of the erythroid cell is to participate in the immune response. Mature erythrocytes from chickens express Toll-like receptors and several cytokines in response to stimulation of the immune system. We previously reported the application of a biochemical fractionation protocol to isolate highly enriched transcribed DNA from polychromatic erythrocytes from chickens. In conjunction with next-generation DNA, RNA sequencing, chromatin immunoprecipitation-DNA sequencing, and formaldehyde-assisted isolation of regulatory elements (FAIRE) sequencing, we identified the active chromosomal compartments and determined their structural signatures in relation to expression levels. Here, we present the detailed chromatin characteristics of erythroid genes participating in the innate immune response. Our studies revealed an atypical chromatin structure for several genes coding for Toll-like receptors, interleukins, and interferon regulatory factors. The body of these genes had nucleosome-free regions intermingled with nucleosomes modified with H3K4me3 and H3K27ac, suggesting a dynamic unstable chromatin structure. We further show that human genes involved in cell identity have gene bodies with the same chromatin-instability features as the chicken polychromatic erythrocyte genes participating in the innate immune response.

The evolution of CpG density and lifespan in conserved primate and mammalian promoters

Gene promoters are evolutionarily conserved across holozoans and enriched in CpG sites, the target for DNA methylation. As animals age, the epigenetic pattern of DNA methylation degrades, with highly methylated CpG sites gradually becoming demethylated while CpG islands increase in methylation. Across vertebrates, aging is a trait that varies among species. We used this variation to determine whether promoter CpG density correlates with species’ maximum lifespan. Human promoter sequences were used to identify conserved regions in 131 mammals and a subset of 28 primate genomes. We identified approximately 1000 gene promoters (5% of the total), that significantly correlated CpG density with lifespan. The correlations were performed via the phylogenetic least squares method to account for trait similarity by common descent using phylogenetic branch lengths. Gene set enrichment analysis revealed no significantly enriched pathways or processes, consistent with the hypothesis that aging is not under positive selection. However, within both mammals and primates, 95% of the promoters showed a positive correlation between increasing CpG density and species lifespan, and two thirds were shared between the primate subset and mammalian datasets. Thus, these genes may require greater buffering capacity against age-related dysregulation of DNA methylation in longer-lived species.

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.

Systematic analysis of transcription start sites in avian development

Cap Analysis of Gene Expression (CAGE) in combination with single-molecule sequencing technology allows precision mapping of transcription start sites (TSSs) and genome-wide capture of promoter activities in differentiated and steady state cell populations. Much less is known about whether TSS profiling can characterize diverse and non-steady state cell populations, such as the approximately 400 transitory and heterogeneous cell types that arise during ontogeny of vertebrate animals. To gain such insight, we used the chick model and performed CAGE-based TSS analysis on embryonic samples covering the full 3-week developmental period. In total, 31,863 robust TSS peaks (>1 tag per million [TPM]) were mapped to the latest chicken genome assembly, of which 34% to 46% were active in any given developmental stage. ZENBU, a web-based, open-source platform, was used for interactive data exploration. TSSs of genes critical for lineage differentiation could be precisely mapped and their activities tracked throughout development, suggesting that non-steady state and heterogeneous cell populations are amenable to CAGE-based transcriptional analysis. Our study also uncovered a large set of extremely stable housekeeping TSSs and many novel stage-specific ones. We furthermore demonstrated that TSS mapping could expedite motif-based promoter analysis for regulatory modules associated with stage-specific and housekeeping genes. Finally, using Brachyury as an example, we provide evidence that precise TSS mapping in combination with Clustered Regularly Interspaced Short Palindromic Repeat (CRISPR)-on technology enables us, for the first time, to efficiently target endogenous avian genes for transcriptional activation. Taken together, our results represent the first report of genome-wide TSS mapping in birds and the first systematic developmental TSS analysis in any amniote species (birds and mammals). By facilitating promoter-based molecular analysis and genetic manipulation, our work also underscores the value of avian models in unravelling the complex regulatory mechanism of cell lineage specification during amniote development.

G-quadruplexes in human promoters: A challenge for therapeutic applications

BACKGROUND

G-rich sequences undergo unique structural equilibria to form G-quadruplexes (G4) both in vitro and in cell systems. Several pathologies emerged to be directly related to G4 occurrence at defined genomic portions. Additionally, G-rich sequences are significantly represented around transcription start sites (TSS) thus leading to the hypothesis of a gene regulatory function for G4. Thus, the tuning of G4 formation has been proposed as a new powerful tool to regulate gene expression to treat related pathologies. However, up-to date this approach did not provide any new really efficient treatment.

SCOPE OF REVIEW

Here, we summarize the most recent advances on the correlation between the structural features of G4 in human promoters and the role these systems physiologically exert. In particular we focus on the effect of G4 localization among cell compartments and along the promoters in correlation with protein interaction networks and epigenetic state. Finally the intrinsic structural features of G4 at promoters are discussed to unveil the contribution of different G4 structural modules in this complex architecture.

MAJOR CONCLUSIONS

It emerges that G4s play several roles in the intriguing and complex mechanism of gene expression, being able to produce opposite effects on the same target. This reflects the occurrence of a highly variegate network of several components working simultaneously.

GENERAL SIGNIFICANCE

The resulting picture is still fuzzy but some points of strength are definitely emerging, which prompts all of us to strengthen our efforts in view of a selective control of gene expression through G4 modulation. This article is part of a Special Issue entitled "G-quadruplex" Guest Editor: Dr. Concetta Giancola and Dr. Daniela Montesarchio.

FoxP2 directly regulates the reelin receptor VLDLR developmentally and by singing

Mutations of the transcription factor FOXP2 cause a severe speech and language disorder. In songbirds, FoxP2 is expressed in the medium spiny neurons (MSNs) of the avian basal ganglia song nucleus, Area X, which is crucial for song learning and adult song performance. Experimental downregulation of FoxP2 in Area X affects spine formation, prevents neuronal plasticity induced by social context and impairs song learning. Direct target genes of FoxP2 relevant for song learning and song production are unknown. Here we show that a lentivirally mediated FoxP2 knockdown in Area X of zebra finches downregulates the expression of VLDLR, one of the two reelin receptors. Zebra finch FoxP2 binds to the promoter of VLDLR and activates it, establishing VLDLR as a direct FoxP2 target. Consistent with these findings, VLDLR expression is co-regulated with FoxP2 as a consequence of adult singing and during song learning. We also demonstrate that knockdown of FoxP2 affects glutamatergic transmission at the corticostriatal MSN synapse. These data raise the possibility that the regulatory relationship between FoxP2 and VLDLR guides structural plasticity towards the subset of FoxP2-positive MSNs in an activity dependent manner via the reelin pathway.

Evolutionary Footprints of Short Tandem Repeats in Avian Promoters

Short tandem repeats (STRs) or microsatellites are well-known sequence elements that may change the spacing between transcription factor binding sites (TFBSs) in promoter regions by expansion or contraction of repetitive units. Some of these mutations have the potential to contribute to phenotypic diversity by altering patterns of gene expression. To explore how repetitive sequence motifs within promoters have evolved in avian lineages under mutation-selection balance, more than 400 evolutionary conserved STRs (ecSTRs) were identified in this study by comparing the 2 kb upstream promoter sequences of chicken against those of other birds (turkey, duck, zebra finch, and flycatcher). The rate of conservation was significantly higher in AG dinucleotide repeats than in AC or AT repeats, with the expansion of AG motifs being noticeably constrained in passerines. Analysis of the relative distance between ecSTRs and TFBSs revealed a significantly higher rate of conserved TFBSs in the vicinity of ecSTRs in both chicken-duck and chicken-passerine comparisons. Our comparative study provides a novel insight into which intrinsic factors have influenced the degree of constraint on repeat expansion/contraction during avian promoter evolution.

Genetic Correlates of Individual Differences in Sleep Behavior of Free-Living Great Tits (Parus major)

Within populations, free-living birds display considerable variation in observable sleep behaviors, reflecting dynamic interactions between individuals and their environment. Genes are expected to contribute to repeatable between-individual differences in sleep behaviors, which may be associated with individual fitness. We identified and genotyped polymorphisms in nine candidate genes for sleep, and measured five repeatable sleep behaviors in free-living great tits (Parus major), partly replicating a previous study in blue tits (Cyanistes caeruleus). Microsatellites in the CLOCK and NPAS2 clock genes exhibited an association with sleep duration relative to night length, and morning latency to exit the nest box, respectively. Furthermore, microsatellites in the NPSR1 and PCSK2 genes associated with relative sleep duration and proportion of time spent awake at night, respectively. Given the detection rate of associations in the same models run with random markers instead of candidate genes, we expected two associations to arise by chance. The detection of four associations between candidate genes and sleep, however, suggests that clock genes, a clock-related gene, or a gene involved in the melanocortin system, could play key roles in maintaining phenotypic variation in sleep behavior in avian populations. Knowledge of the genetic architecture underlying sleep behavior in the wild is important because it will enable ecologists to assess the evolution of sleep in response to selection.

Characteristics of Antisense Transcript Promoters and the Regulation of Their Activity

Recently, an increasing number of studies on natural antisense transcripts have been reported, especially regarding their classification, temporal and spatial expression patterns, regulatory functions and mechanisms. It is well established that natural antisense transcripts are produced from the strand opposite to the strand encoding a protein. Despite the pivotal roles of natural antisense transcripts in regulating the expression of target genes, the transcriptional mechanisms initiated by antisense promoters (ASPs) remain unknown. To date, nearly all of the studies conducted on this topic have focused on the ASP of a single gene of interest, whereas no study has systematically analyzed the locations of ASPs in the genome, ASP activity, or factors influencing this activity. This review focuses on elaborating on and summarizing the characteristics of ASPs to extend our knowledge about the mechanisms of antisense transcript initiation.

Human endogenous retroviruses sustain complex and cooperative regulation of gene-containing loci and unannotated megabase-sized regions

Background

Evidence suggests that some human endogenous retroviruses and endogenous retrovirus-like repeats (here collectively ERVs) regulate the expression of neighboring genes in normal and disease states; e.g. the human globin locus is regulated by an ERV9 that coordinates long-range gene switching during hematopoiesis and activates also intergenic transcripts. While complex transcription regulation is associated with integration of certain exogenous retroviruses, comparable regulation sustained by ERVs is less understood.

Findings

We analyzed ERV transcription using ERV9 consensus sequences and publically available RNA-sequencing, chromatin immunoprecipitation with sequencing (ChIP-seq) and cap analysis gene expression (CAGE) data from ENCODE. We discovered previously undescribed and advanced transcription regulation mechanisms in several human reference cell lines. We show that regulation by ERVs involves long-ranging activations including complex RNA splicing patterns, and transcription of large unannotated regions ranging in size from several hundred kb to around 1 Mb. Moreover, regulation was found to be cooperatively sustained in some loci by multiple ERVs and also non-LTR repeats.

Conclusion

Our analyses show that endogenous retroviruses sustain advanced transcription regulation in human cell lines, which shows similarities to complex insertional mutagenesis effects exerted by exogenous retroviruses. By exposing previously undescribed regulation effects, this study should prove useful for understanding fundamental transcription mechanisms resulting from evolutionary acquisition of retroviral sequence in the human genome.

Electronic supplementary material

The online version of this article (doi:10.1186/s12977-015-0161-9) contains supplementary material, which is available to authorized users.