Putative Regulatory Sequences Research Articles

Functional and regulatory diversification of Period genes responsible for circadian rhythm in vertebrates.

The Period genes (Per) play essential roles in modulating the molecular circadian clock timing in a broad range of species, which regulates the physiological and cellular rhythms through the transcription-translation feedback loop. While the Period gene paralogs are widely observed among vertebrates, the evolutionary history and the functional diversification of Per genes across vertebrates are not well known. In this study, we comprehensively investigated the evolution of Per genes at the copy number and sequence levels, including de novo binding motif discovery by comparative genomics. We also determined the lineage-specific transcriptome landscape across tissues and developmental stages and phenotypic effects in public RNA-seq data sets of model species. We observed multiple lineage-specific gain and loss events of Per genes, though no simple association was observed between ecological factors and Per gene numbers in each species. Among salmonid fish species, the per3 gene has been lost in the majority, whereas those retaining the per3 gene exhibit not a signature of relaxed selective constraint but rather a signature of intensified selection. We also determined the signature of adaptive diversification of the CRY-binding region in Per1 and Per3, which modulates the circadian rhythm. We also discovered putative regulatory sequences, which are lineage-specific, suggesting that these cis-regulatory elements may have evolved rapidly and divergently across different lineages. Collectively, our findings revealed the evolution of Per genes and their fine-tuned contribution to the plastic and precise regulation of circadian rhythms in various vertebrate taxa.

Phenotypic and genetic insights into efflux pump mechanism in Mycoplasma anserisalpingitidis.

Mycoplasma anserisalpingitidis is one of the most important waterfowl-pathogenic mycoplasmas. Due to inadequate antibiotic treatment, many strains with high minimal inhibitory concentration (MIC) values for multiple drugs have been isolated lately. Decreased antibiotic susceptibility in several Mycoplasma species are known to be associated with mutations in topoisomerase and ribosomal genes, but other strategies such as active efflux pump mechanisms were also described. The scope of this study was the phenotypic and genetic characterization of the active efflux mechanism in M. anserisalpingitidis. We measured the MIC values in the presence and absence of different efflux pump inhibitors (EPIs), such as carbonyl cyanide m-chlorophenylhydrazine (CCCP), orthovanadate (OV), and reserpine (RSP). Moreover, bioinformatic tools were utilized to detect putative regulatory sequences of membrane transport proteins coding genes, while comparative genome analysis was performed to reveal potential markers of antibiotic resistance. Out of the three examined EPIs, CCCP decreased the MICs at least two-fold below the original MICs (in 23 cases out of 36 strains). In the presence of OV or RSP, MIC value differences could be seen only if modified dilution series (10% decrease steps were used instead of two-fold dilutions) were applied (in 24/36 cases with OV and 9/36 with RSP). During comparative genome analysis, non-synonymous single nucleotide polymorphisms (nsSNPs) were identified in genes encoding ABC membrane transport proteins, which were displayed in higher percentages in M. anserisalpingitidis strains with increased MICs. In terms of other genes, a nsSNP was identified in DNA gyrase subunit A (gyrA) gene which can be related to decreased susceptibility to enrofloxacin. The present study is the first to highlight the importance of efflux pump mechanisms in M. anserisalpingitidis. Considering the observed effects of the EPI CCCP against this bacterium, it can be assumed, that the use of EPIs would increase the efficiency of targeted antibiotic therapy in the future control of this pathogen. However, further research is required to obtain a more comprehensive understanding of efflux pump mechanism in this bacterium.

Functional and comparative genomics reveals conserved noncoding sequences in the nitrogen-fixing clade.

Summary Nitrogen is one of the most inaccessible plant nutrients, but certain species have overcome this limitation by establishing symbiotic interactions with nitrogen‐fixing bacteria in the root nodule. This root–nodule symbiosis (RNS) is restricted to species within a single clade of angiosperms, suggesting a critical, but undetermined, evolutionary event at the base of this clade.To identify putative regulatory sequences implicated in the evolution of RNS, we evaluated the genomes of 25 species capable of nodulation and identified 3091 conserved noncoding sequences (CNS) in the nitrogen‐fixing clade (NFC).We show that the chromatin accessibility of 452 CNS correlates significantly with the regulation of genes responding to lipochitooligosaccharides in Medicago truncatula. These included 38 CNS in proximity to 19 known genes involved in RNS. Five such regions are upstream of MtCRE1, Cytokinin Response Element 1, required to activate a suite of downstream transcription factors necessary for nodulation in M. truncatula. Genetic complementation of an Mtcre1 mutant showed a significant decrease of nodulation in the absence of the five CNS, when they are driving the expression of a functional copy of MtCRE1.CNS identified in the NFC may harbor elements required for the regulation of genes controlling RNS in M. truncatula.

Targeted Next-Generation Sequencing Reveals Mutations in Non-coding Regions and Potential Regulatory Sequences of Calpain-3 Gene in Polish Limb-Girdle Muscular Dystrophy Patients.

Limb–girdle muscular dystrophy type R1 (LGMDR1) is caused by mutations in CAPN3 and is the most common type of recessive LGMD. Even with the use of whole-exome sequencing (WES), only one mutant allele of CAPN3 is found in a significant number of LGMDR patients. This points to a role of non-coding, intronic or regulatory, sequence variants in the disease pathogenesis. Targeted sequencing of the whole CAPN3 gene including not only intronic, 3′ and 5′ UTRs but also potential regulatory regions was performed in 27 patients suspected with LGMDR1. This group included 13 patients with only one mutated CAPN3 allele detected previously with exome sequencing. A second rare variant in the non-coding part of CAPN3 was found in 11 of 13 patients with previously identified single mutation. Intronic mutations were found in 10 cases, with c.1746-20C>G variant present in seven patients. In addition, a large deletion of exons 2–8 was found in one patient. In the patients with no causative mutation previously found, we detected rare CAPN3 variants in 5 out of 10 patients and in two of them in a compound heterozygous state. Rare variants within putative regulatory sequences distant from the CAPN3 gene were found in 15 patients, although in 11 of these cases, other variants are deemed causative. The results indicate that intronic mutations are common in Polish LGMDR patients, and testing for non-coding mutations in CAPN3 should be performed in apparently single heterozygous patients.

Global properties of regulatory sequences are predicted by transcription factor recognition mechanisms

BackgroundMammalian genomes contain millions of putative regulatory sequences, which are delineated by binding of multiple transcription factors. The degree to which spacing and orientation constraints among transcription factor binding sites contribute to the recognition and identity of regulatory sequence is an unresolved but important question that impacts our understanding of genome function and evolution. Global mechanisms that underlie phenomena including the size of regulatory sequences, their uniqueness, and their evolutionary turnover remain poorly described.ResultsHere, we ask whether models incorporating different degrees of spacing and orientation constraints among transcription factor binding sites are broadly consistent with several global properties of regulatory sequence. These properties include length, sequence diversity, turnover rate, and dominance of specific TFs in regulatory site identity and cell type specification. Models with and without spacing and orientation constraints are generally consistent with all observed properties of regulatory sequence, and with regulatory sequences being fundamentally small (~ 1 nucleosome). Uniqueness of regulatory regions and their rapid evolutionary turnover are expected under all models examined. An intriguing issue we identify is that the complexity of eukaryotic regulatory sites must scale with the number of active transcription factors, in order to accomplish observed specificity.ConclusionsModels of transcription factor binding with or without spacing and orientation constraints predict that regulatory sequences should be fundamentally short, unique, and turn over rapidly. We posit that the existence of master regulators may be, in part, a consequence of evolutionary pressure to limit the complexity and increase evolvability of regulatory sites.

Molecular Characterization of Choroideremia-Associated Deletions Reveals an Unexpected Regulation of CHM Gene Transcription.

Choroideremia (CHM) is a X-linked recessive chorioretinal dystrophy due to deficiency of the CHM gene product, i.e., Rab escort protein isoform 1 (REP1). To date, gene therapy for CHM has shown variable effectiveness, likely because the underlying pathogenic mechanisms as well as genotype-phenotype correlation are not yet fully known. Small nucleotide variants leading to premature termination codons (PTCs) are a major cause of CHM, but about 20% of patients has CHM gene deletions. To improve understanding of the disease mechanisms, we analyzed molecular features of seven deletions involving the CHM gene sequence. We mapped the deletion breakpoints by using polymerase chain reaction, sequencing and array comparative genomic hybridization; to identify rearrangement-promoting DNA sequences, we analyzed genomic architecture surrounding the breakpoint regions. Moreover, in some CHM patients with different mutation types, we measured transcript level of CHM and of CHML, encoding the REP2 isoform. Scattered along the whole CHM gene and in close proximity to the deletion breakpoints we found numerous repeat elements that generate a locus-specific rearrangement hot spot. Unexpectedly, patients with non-PTC variants had increased expression of the aberrant CHM mRNA; CHML expression was higher than normal in a patient lacking CHM and its putative regulatory sequences. This latest evidence suggests that mechanisms regulating CHM and CHML gene expression are worthy of further study, because their full knowledge could be also useful for developing effective therapies for this hitherto untreatable inherited retinal degeneration.

Enhancer Pleiotropy, Gene Expression, and the Architecture of Human Enhancer-Gene Interactions.

Enhancers are often studied as noncoding regulatory elements that modulate the precise spatiotemporal expression of genes in a highly tissue-specific manner. This paradigm has been challenged by recent evidence of individual enhancers acting in multiple tissues or developmental contexts. However, the frequency of these enhancers with high degrees of “pleiotropy” out of all putative enhancers is not well understood. Consequently, it is unclear how the variation of enhancer pleiotropy corresponds to the variation in expression breadth of target genes. Here, we use multi-tissue chromatin maps from diverse human tissues to investigate the enhancer–gene interaction architecture while accounting for 1) the distribution of enhancer pleiotropy, 2) the variations of regulatory links from enhancers to target genes, and 3) the expression breadth of target genes. We show that most enhancers are tissue-specific and that highly pleiotropy enhancers account for <1% of all putative regulatory sequences in the human genome. Notably, several genomic features are indicative of increasing enhancer pleiotropy, including longer sequence length, greater number of links to genes, increasing abundance and diversity of encoded transcription factor motifs, and stronger evolutionary conservation. Intriguingly, the number of enhancers per gene remains remarkably consistent for all genes (∼14). However, enhancer pleiotropy does not directly translate to the expression breadth of target genes. We further present a series of Gaussian Mixture Models to represent this organization architecture. Consequently, we demonstrate that a modest trend of more pleiotropic enhancers targeting more broadly expressed genes can generate the observed diversity of expression breadths in the human genome.

Comparative transcriptome of compatible and incompatible interaction of Erysiphe pisi and garden pea reveals putative defense and pathogenicity factors.

Comparative transcriptome analysis of Erysiphe pisi-infected pea (Pisum sativum) genotypes JI-2480 (resistant) and Arkel (susceptible) at 72 hours post-inoculation (hpi) was carried to detect molecular components involved in compatible and incompatible interactions. Differential gene expression was observed in Arkel and JI-2480 genotype at 72 hpi with E. pisi isolate (Ep01) using EdgeR software. Out of 32 217 transcripts, 2755 transcripts showed significantly altered gene expression in case of plants while 530 were related to E. pisi (P <0.05). The higher transcript number of differentially expressed genes demonstrated peak activity of pathogenicity genes in plants at 72 hpi. Glycolysis was observed to be the major pathway for energy source during fungal growth. Differential gene expression of plant transcripts revealed significant expression of putative receptor and regulatory sequences involved in defense in the resistant, JI-2480 compared to susceptible, Arkel genotype. Expression of genes involved in defense and hormonal signaling, genes related to hypersensitive response, reactive oxygen species and phenylpropanoid pathway in JI-2480 indicated their crucial role in disease resistance against E. pisi. Down-regulation of transcription factors like-WRKY-28 and up-regulation of several putative pattern recognition receptors in JI-2480 compared to Arkel also suggested activation of host-mediated defense responses against E. pisi in pea.

The impact of different negative training data on regulatory sequence predictions.

Regulatory regions, like promoters and enhancers, cover an estimated 5-15% of the human genome. Changes to these sequences are thought to underlie much of human phenotypic variation and a substantial proportion of genetic causes of disease. However, our understanding of their functional encoding in DNA is still very limited. Applying machine or deep learning methods can shed light on this encoding and gapped k-mer support vector machines (gkm-SVMs) or convolutional neural networks (CNNs) are commonly trained on putative regulatory sequences. Here, we investigate the impact of negative sequence selection on model performance. By training gkm-SVM and CNN models on open chromatin data and corresponding negative training dataset, both learners and two approaches for negative training data are compared. Negative sets use either genomic background sequences or sequence shuffles of the positive sequences. Model performance was evaluated on three different tasks: predicting elements active in a cell-type, predicting cell-type specific elements, and predicting elements' relative activity as measured from independent experimental data. Our results indicate strong effects of the negative training data, with genomic backgrounds showing overall best results. Specifically, models trained on highly shuffled sequences perform worse on the complex tasks of tissue-specific activity and quantitative activity prediction, and seem to learn features of artificial sequences rather than regulatory activity. Further, we observe that insufficient matching of genomic background sequences results in model biases. While CNNs achieved and exceeded the performance of gkm-SVMs for larger training datasets, gkm-SVMs gave robust and best results for typical training dataset sizes without the need of hyperparameter optimization.

Cationic porphyrins as destabilizer of a G-quadruplex located at the promoter of human MYH7 β gene

G-quadruplex (GQ) architecture is adopted by guanine rich sequences, present throughout the eukaryotic genome including promoter locations and telomeric ends. The in vivo presence indicates their involvement and role in various biological processes. Various small ligands have been developed to interact and stabilize/destabilize G-quadruplex structures. Cationic porphyrins are among the most studied ligands, reported to bind and stabilize G-quadruplexes. Herein, we report the recognition and destabilization of a parallel G-quadruplex by porphyrins (TMPyP3 and TMPyP4). This G-quadruplex forming 23-nt G-rich sequence is in the promoter region of Human Myosin Heavy Chain β gene (MYH7β). Presence of various putative regulatory sequence elements (TATA Box, CCAAT, SP-1) located in the vicinity of this quadruplex motif, highlight its regulatory implications. Biophysical methods as Circular Dichroism Spectroscopy, UV-Absorption Spectroscopy, UV-Thermal Denaturation and Fluorescence Spectroscopy (steady as well as Time Resolved) have been used for studying the interaction and binding parameters. It is proposed that porphyrins have a destabilizing effect on the G-quadruplexes with parallel topology and a stronger binding specifically via intercalation mode is needed to cause destabilization. The study deals with better understanding and insights of DNA-Drug interactions in biological systems.Communicated by Ramaswamy H. Sarma

Novel cis-acting regulatory elements in wild Oryza species impart improved rice bran quality by lowering the expression of phospholipase D alpha1 enzyme (OsPLDα1).

Rice bran oil is good quality edible oil, rich in antioxidants and comprised typically of oleic-linoleic type fatty acids. However, presence of a highly lipolytic enzyme Phospholipase D alpha1 (OsPLDα1) increases free fatty acid content in the oil which further leads to stale flavor and rancidity of the oil, making it unfit for human consumption. In this study, we compared the upstream regions of OsPLDα1 orthologs across 34 accessions representing 5 wild Oryza species and 8 cultivars, to uncover sequence variations and identify cis-elements involved in differential transcription of orthologs. Alignment of the upstream sequences to the Nipponbare OsPLDα1 reference sequence revealed the presence of 39 SNPs. Phylogenetic analysis showed that all the selected cultivars and wild species accessions are closely related to the reference except for three accessions of O. rufipogon (IRGC89224, IRGC104425, and IRGC105902). Furthermore, using exon-specific qRT-PCR, OsPLDα1 expression patterns in immature grains indicated significant differences in transcript abundance between the wild species accessions. In comparison to the control, lowest gene expression was observed in IRGC89224 accession (0.20-fold) followed by IRGC105902 (0.26-fold) and IRGC104425 (0.41-fold) accessions. In-silico analysis of the OsPLDα1 promoter revealed that the copy number variations of CGCGBOXAT, GT1CONSENSUS, IBOXCORE, NODCON2GM, OSE2ROOTNODULE, SURECOREATSULTR11, and SORLIP1AT cis-elements play an important role in the transcriptional activities of orthologous genes. Owing to the presence of ARFAT and SEBF elements only in the IRGC89224 accession, which had the lowest gene expression as well, these putative upstream regulatory sequences have been identified as novel cis-elements which may act as repressors in regulating the OsPLDα1 gene expression. The accessions identified with low OsPLDα1 expressions could be further deployed as potential donors of ideal OsPLDα1 allele for transfer of the desired trait into elite rice cultivars.

Mitochondrial genomic variation drives differential nuclear gene expression in discrete regions of Drosophila gene and protein interaction networks

BackgroundMitochondria perform many key roles in their eukaryotic hosts, from integrating signaling pathways through to modulating whole organism phenotypes. The > 1 billion years of nuclear and mitochondrial gene co-evolution has necessitated coordinated expression of gene products from both genomes that maintain mitochondrial, and more generally, eukaryotic cellular function. How mitochondrial DNA (mtDNA) variation modifies host fitness has proved a challenging question but has profound implications for evolutionary and medical genetics. In Drosophila, we have previously shown that recently diverged mtDNA haplotypes within-species can have more impact on organismal phenotypes than older, deeply diverged haplotypes from different species. Here, we tested the effects of mtDNA haplotype variation on gene expression in Drosophila under standardized conditions. Using the Drosophila Genetic Reference Panel (DGRP), we constructed a panel of mitonuclear genotypes that consists of factorial variation in nuclear and mtDNA genomes, with mtDNAs originating in D. melanogaster (2x haplotypes) and D. simulans (2x haplotypes).ResultsWe show that mtDNA haplotype variation unequivocally alters nuclear gene expression in both females and males, and mitonuclear interactions are pervasive modifying factors for gene expression. There was appreciable overlap between the sexes for mtDNA-sensitive genes, and considerable transcriptional variation attributed to particular mtDNA contrasts. These genes are generally found in low-connectivity gene co-expression networks, occur in gene clusters along chromosomes, are often flanked by non-coding RNA, and are under-represented among housekeeping genes. Finally, we identify the giant (gt) transcription factor motif as a putative regulatory sequence associated with mtDNA-sensitive genes.ConclusionsThere are predictive conditions for nuclear genes that are influenced by mtDNA variation.

Identification of Putative UL54 (ICP27) Transcription Regulatory Sequences Binding to Oct-1, v-Myb, Pax-6 and Hairy in Herpes Simplex Viruses.

An oncolytic herpes simplex virus (oHSV) has proven amenable in oncolytic virotherapy and was approved to treat melanoma. The immediate-early (IE) protein ICP27 encoded by gene UL54 is essential for HSV infection. Post-transcriptional modification of UL54 would increase tumor targeting of oHSVs. However, UL54 gene transcription regulatory sequences and factors were not reported yet. Here we isolated a new strain LXMW of type 1 HSV (HSV-1-LXMW) in China and found it's closely related to HSV-1 strains Patton and H129 in the US by the first and next generation DNA sequencing viral DNA phylogenetic analysis. Using a weight matrix-based program Match, we found the UL54 transcription regulatory sequences binding to the transcription factors Oct-1, v-Myb and Pax-6 in HSV-1-LXMW, while the sequences binding to Oct-1 and Hairy in a HSV-2 strain. Further validation showed that HSV-1 and HSV-2 shared the common sequence binding to Oct-1, but had unique sequences to bind v-Myb and Pax-6, or Hairy, respectively, by DNA sequence alignment of total 11 HSV strains. The published results howed that the expression of transcription factors is consistent with the tissue tropism of HSV-1 and HSV-2. In the current article a new HSV-1 strain LXMW was isolated and its putative HSV UL54 transcription regulatory sequences and factors were identified for the first time. Our findings highlight the new understanding of the principles of transcriptional regulation in HSV biology and oncolytic virotherapy.

LTR-retrotransposon transcriptome modulation in response to endotoxin-induced stress in PBMCs

BackgroundHuman Endogenous Retroviruses (HERVs) and Mammalian apparent LTR-retrotransposons (MaLRs) represent the 8% of our genome and are distributed among our 46 chromosomes. These LTR-retrotransposons are thought to be essentially silent except in cancer, autoimmunity and placental development. Their Long Terminal Repeats (LTRs) constitute putative promoter or polyA regulatory sequences. In this study, we used a recently described high-density microarray which can be used to study HERV/MaLR transcriptome including 353,994 HERV/MaLR loci and 1559 immunity-related genes.ResultsWe described, for the first time, the HERV transcriptome in peripheral blood mononuclear cells (PBMCs) using a cellular model mimicking inflammatory response and monocyte anergy observed after septic shock. About 5.6% of the HERV/MaLR repertoire is transcribed in PBMCs. Roughly one-tenth [5.7–13.1%] of LTRs exhibit a putative constitutive promoter or polyA function while one-quarter [19.5–27.6%] may shift from silent to active. Evidence was given that some HERVs/MaLRs and genes may share similar regulation control under lipopolysaccharide (LPS) stimulation conditions. Stimulus-dependent response confirms that HERV expression is tightly regulated in PBMCs. Altogether, these observations make it possible to integrate 62 HERVs/MaLRs and 26 genes in 11 canonical pathways and suggest a link between HERV expression and immune response. The transcriptional modulation of HERVs located close to genes such as OAS2/3 and IFI44/IFI44L or at a great distance from genes was discussed.ConclusionThis microarray-based approach revealed the expression of about 47,466 distinct HERV loci and identified 951 putative promoter LTRs and 744 putative polyA LTRs in PBMCs. HERV/MaLR expression was shown to be tightly modulated under several stimuli including high-dose and low-dose LPS and Interferon-γ (IFN-γ). HERV incorporation at the crossroads of immune response pathways paves the way for further functional studies and analyses of the HERV transcriptome in altered immune responses in vivo such as in sepsis.

Organophosphate pesticide chlorpyrifos impairs STAT1 signaling to induce dopaminergic neurotoxicity: Implications for mitochondria mediated oxidative stress signaling events

The organophosphate (OP) pesticide chlorpyrifos (CPF), used in agricultural settings, induces developmental and neurological impairments. Recent studies using in vitro cell culture models have reported CPF exposure to have a positive association with mitochondria-mediated oxidative stress response and dopaminergic cell death; however, the mechanism by which mitochondrial reactive oxygen species (ROS) contribute to dopaminergic cell death remains unclear. Therefore, we hypothesized that STAT1, a transcription factor, causes apoptotic dopaminergic cell death via mitochondria-mediated oxidative stress mechanisms. Here we show that exposure of dopaminergic neuronal cells such as N27 cells (immortalized murine mesencephalic dopaminergic cells) to CPF resulted in a dose-dependent increase in apoptotic cell death as measured by MTS assay and DNA fragmentation. Similar effects were observed in CPF-treated human dopaminergic neuronal cells (LUHMES cells), with an associated increase in mitochondrial dysfunction. Moreover, CPF (10 μM) induced time-dependent increase in STAT1 activation coincided with the collapse of mitochondrial transmembrane potential, increase in ROS generation, proteolytic cleavage of protein kinase C delta (PKCδ), inhibition of the mitochondrial basal oxygen consumption rate (OCR), with a concomitant reduction in ATP-linked OCR and reserve capacity, increase in Bax/Bcl-2 ratio and enhancement of autophagy. Additionally, by chromatin immunoprecipitation (ChIP), we demonstrated that STAT1 bound to a putative regulatory sequence in the NOX1 and Bax promoter regions in response to CPF in N27 cells. Interestingly, overexpression of non-phosphorylatable STAT1 mutants (STAT1Y701F and STAT1S727A) but not STAT1 WT construct attenuated the cleavage of PKCδ and ultimately cell death in CPF-treated cells. Furthermore, small interfering RNA knockdown demonstrated STAT1 to be a critical regulator of autophagy and mitochondria-mediated proapoptotic cell signaling events after CPF treatment in N27 cells. Finally, oral administration of CPF (5 mg/kg) in postnatal rats (PNDs 27–61) induced motor deficits, and nigrostriatal dopaminergic neurodegeneration with a concomitant induction of STAT1-dependent proapoptotic cell signaling events. Conversely, co-treatment with mitoapocynin (a mitochondrially-targeted antioxidant) and CPF rescued motor deficits, and restored dopaminergic neuronal survival via abrogation of STAT1-dependent proapoptotic cell signaling events. Taken together, our study identifies a novel mechanism by which STAT1 regulates mitochondria-mediated oxidative stress response, PKCδ activation and autophagy. In this context, the phosphorylation of Tyrosine 701 and Serine 727 in STAT1 was found to be essential for PKCδ cleavage. By attenuating mitochondrial-derived ROS, mitoapocynin may have therapeutic applications for reversing CPF-induced dopaminergic neurotoxicity and associated neurobehavioral deficits as well as neurodegenerative diseases.

Identification of rare noncoding sequence variants in gamma-aminobutyric acid A receptor, alpha 4 subunit in autism spectrum disorder.

Alterations of the gamma-aminobutyric acid (GABA) signaling system has been strongly linked to the pathophysiology of autism spectrum disorder (ASD). Genetic associations of common variants in GABA receptor subunits, in particular GABRA4 on chromosome 4p12, with ASD have been replicated by several studies. Moreover, molecular investigations have identified altered transcriptional and translational levels of this gene and protein in brains of ASD individuals. Since the genotyped common variants are likely not the functional variants contributing to the molecular consequences or underlying ASD phenotype, this study aims to examine rare sequence variants in GABRA4, including those outside the protein coding regions of the gene. We comprehensively re-sequenced the entire protein coding and noncoding portions of the gene and putative regulatory sequences in 82 ASD individuals and 55 developmentally typical pediatric controls, all homozygous for the most significant previously associated ASD risk allele (G/G at rs1912960). We identified only a single common, coding variant, and no association of any single marker or set of variants with ASD. Functional annotation of noncoding variants identified several rare variants in putative regulatory sites. Finally, a rare variant unique to ASD cases, in an evolutionary conserved site of the 3'UTR, shows a trend toward decreasing gene expression. Hence, GABRA4 rare variants in noncoding DNA may be variants of modest physiological effects in ASD etiology.

Cis-regulatory landscapes of four cell types of the retina.

The retina is composed of ∼50 cell-types with specific functions for the process of vision. Identification of the cis-regulatory elements active in retinal cell-types is key to elucidate the networks controlling this diversity. Here, we combined transcriptome and epigenome profiling to map the regulatory landscape of four cell-types isolated from mouse retinas including rod and cone photoreceptors as well as rare inter-neuron populations such as horizontal and starburst amacrine cells. Integration of this information reveals sequence determinants and candidate transcription factors for controlling cellular specialization. Additionally, we refined parallel reporter assays to enable studying the transcriptional activity of large collection of sequences in individual cell-types isolated from a tissue. We provide proof of concept for this approach and its scalability by characterizing the transcriptional capacity of several hundred putative regulatory sequences within individual retinal cell-types. This generates a catalogue of cis-regulatory regions active in retinal cell types and we further demonstrate their utility as potential resource for cellular tagging and manipulation.

Divergence of paralogous growth-hormone-encoding genes and their promoters in Salmonidae

In many fish species, including salmonids, the growth-hormone is encoded by two duplicated paralogous genes, gh1 and gh2. Both genes were already in place at the time of divergence of species in this group. A comparison of the entire sequence of these genes of salmonids has shown that their conserved regions are associated with exons, while their most variable regions correspond to introns. Introns C and D include putative regulatory elements (sites Pit-1, CRE, and ERE), that are also conserved. In chars, the degree of polymorphism of gh2 gene is 2-3 times as large as that in gh1 gene. However, a comparison across all Salmonidae species would not extent this observation to other species. In both these chars' genes, the promoters are conserved mainly because they correspond to putative regulatory sequences (TATA box, binding sites for the pituitary transcription factor Pit-1 (F1-F4), CRE, GRE and RAR/RXR elements). The promoter of gh2 gene has a greater degree of polymorphism compared with gh1 gene promoter in all investigated species of salmonids. The observed differences in the rates of accumulation of changes in growth hormone encoding paralogs could be explained by differences in the intensity of selection.

Regulation of ULK1 Expression and Autophagy by STAT1

Autophagy involves the lysosomal degradation of cytoplasmic contents for regeneration of anabolic substrates during nutritional or inflammatory stress. Its initiation occurs rapidly after inactivation of the protein kinase mammalian target of rapamycin (mTOR) (or mechanistic target of rapamycin), leading to dephosphorylation of Unc-51-like kinase 1 (ULK1) and autophagosome formation. Recent studies indicate that mTOR can, in parallel, regulate the activity of stress transcription factors, including signal transducer and activator of transcription-1 (STAT1). The current study addresses the role of STAT1 as a transcriptional suppressor of autophagy genes and autophagic activity. We show that STAT1-deficient human fibrosarcoma cells exhibited enhanced autophagic flux as well as its induction by pharmacological inhibition of mTOR. Consistent with enhanced autophagy initiation, ULK1 mRNA and protein levels were increased in STAT1-deficient cells. By chromatin immunoprecipitation, STAT1 bound a putative regulatory sequence in the ULK1 5'-flanking region, the mutation of which increased ULK1 promoter activity, and rendered it unresponsive to mTOR inhibition. Consistent with an anti-apoptotic effect of autophagy, rapamycin-induced apoptosis and cytotoxicity were blocked in STAT1-deficient cells but restored in cells simultaneously exposed to the autophagy inhibitor ammonium chloride. In vivo, skeletal muscle ULK1 mRNA and protein levels as well as autophagic flux were significantly enhanced in STAT1-deficient mice. These results demonstrate a novel mechanism by which STAT1 negatively regulates ULK1 expression and autophagy.

CRISPR-DO for genome-wide CRISPR design and optimization.

Despite the growing popularity in using CRISPR/Cas9 technology for genome editing and gene knockout, its performance still relies on well-designed single guide RNAs (sgRNA). In this study, we propose a web application for the Design and Optimization (CRISPR-DO) of guide sequences that target both coding and non-coding regions in spCas9 CRISPR system across human, mouse, zebrafish, fly and worm genomes. CRISPR-DO uses a computational sequence model to predict sgRNA efficiency, and employs a specificity scoring function to evaluate the potential of off-target effect. It also provides information on functional conservation of target sequences, as well as the overlaps with exons, putative regulatory sequences and single-nucleotide polymorphisms (SNPs). The web application has a user-friendly genome-browser interface to facilitate the selection of the best target DNA sequences for experimental design. CRISPR-DO is available at http://cistrome.org/crispr/ CONTACT: qiliu@tongji.edu.cn or hanxu@jimmy.harvard.edu or xsliu@jimmy.harvard.eduSupplementary information: Supplementary data are available at Bioinformatics online.