Establishment of a cloning-free CRISPR/Cas9 protocol to generate large deletions in the bovine MDBK cell line

The CRISPR/Cas9 technique applied to modify the cattle genome has value in increasing animal health and welfare. Here, we established a simple, fast, and efficient cloning-free CRISPR/Cas9 protocol for large deletions of genomic loci in the frequently used model bovine MDBK cell line. The main advantages of our protocol are as follows: (i) pre-screening of the sgRNA efficiency with a fast and simple cleavage assay, (ii) reliable detection of genomic edits primarily by PCR and confirmed by DNA sequencing, and (iii) single cell sorting with FACS providing specific genetic information from modified cells of interest. Therefore, our method could be successfully applied in different studies, including functional validation of any genetic or regulatory elements.

Characterizing structural variants based on graph-genotyping provides insights into pig domestication and local adaption

Structural variants (SVs), such as deletions (DELs) and insertions (INSs), contribute substantially to pig genetic diversity and phenotypic variation. Using a library of SVs discovered from long-read primary assemblies and short-read sequenced genomes, we map pig genomic SVs with a graph-based method for re-genotyping SVs in 402 genomes. Our results demonstrate that those SVs harboring specific trait-associated genes may greatly shape pig domestication and local adaptation. Further characterization of SVs reveals that some population-stratified SVs may alter the transcription of genes by affecting regulatory elements. We identify that the genotypes of two DELs (296-bp DEL, chr7: 52,172,101-52,172,397; 278-bp DEL, chr18: 23,840,143-23,840,421) located in muscle-specific enhancers are associated with the expression of target genes related to meat quality (FSD2) and muscle fiber hypertrophy (LMOD2 and WASL) in pigs. Our results highlight the role of SVs in domestic porcine evolution, and the identified candidate functional genes and SVs are valuable resources for future genomic research and breeding programs in pigs.

Co-culture platform for tuning of cancer receptor density allows for evaluation of bispecific immune cell engagers

Cancer immunotherapy, where a patient's immune system is harnessed to eradicate cancer cells selectively, is a leading strategy for cancer treatment. However, successes with immune checkpoint inhibitors (ICI) are hampered by reported systemic and organ-specific toxicities and by two-thirds of the patients being non-responders or subsequently acquiring resistance to approved ICIs. Hence substantial efforts are invested in discovering novel targeted immunotherapies aimed at reduced side-effects and improved potency. One way is utilizing the dual targeting feature of bispecific antibodies, which have made them increasingly popular for cancer immunotherapy. Easy and predictive screening methods for activation ranking of candidate drugs in tumor contra non-tumor environments are however lacking. Herein, we present a cell-based assay mimicking the tumor microenvironment by co-culturing B cells with engineered human embryonic kidney 293 T cells (HEK293T), presenting a controllable density of platelet-derived growth factor receptor β (PDGFRβ). A target density panel with three different surface protein levels on HEK293T cells was established by genetic constructs carrying regulatory elements limiting RNA translation of PDGFRβ. We employed a bispecific antibody-affibody construct called an AffiMab capable of binding PDGFRβ on cancer cells and CD40 expressed by B cells as a model. Specific activation of CD40-mediated signaling of immune cells was demonstrated with the two highest receptor-expressing cell lines, Level 2/3 and Level 4, while low-to-none in the low-expressing cell lines. The concept of receptor tuning and the presented co-culture protocol may be of general utility for assessing and developing novel bi-specific antibodies for immuno-oncology applications.

FORGEdb: a tool for identifying candidate functional variants and uncovering target genes and mechanisms for complex diseases

The majority of disease-associated variants identified through genome-wide association studies are located outside of protein-coding regions. Prioritizing candidate regulatory variants and gene targets to identify potential biological mechanisms for further functional experiments can be challenging. To address this challenge, we developed FORGEdb ( https://forgedb.cancer.gov/ ; https://forge2.altiusinstitute.org/files/forgedb.html ; and https://doi.org/10.5281/zenodo.10067458 ), a standalone and web-based tool that integrates multiple datasets, delivering information on associated regulatory elements, transcription factor binding sites, and target genes for over 37 million variants. FORGEdb scores provide researchers with a quantitative assessment of the relative importance of each variant for targeted functional experiments.

Prediction and Analysis of Transcription Factor Binding Sites: Practical Examples and Case Studies Using R Programming

Transcription factors (TFs) bind to specific regions of DNA known as transcription factor binding sites (TFBSs) and modulate gene expression by interacting with the transcriptional machinery. TFBSs are typically located upstream of target genes, within a few thousand base pairs of the transcription start site. The binding of TFs to TFBSs influences the recruitment of the transcriptional machinery, thereby regulating gene transcription in a precise and specific manner. This chapter provides practical examples and case studies demonstrating the extraction of upstream gene regions from the genome, identification of TFBSs using PWMEnrich R/Bioconductor package, interpretation of results, and preparation of publication-ready figures and tables. The EOMES promoter is used as a case study for single DNA sequence analysis, revealing potential regulation by the LHX9-FOXP1 complex during embryonic development. Additionally, an example is presented on how to investigate TFBSs in the upstream regions of a group of genes, using a case study of differentially expressed genes in response to human parainfluenza virus type 1 (HPIV1) infection and interferon-beta. Key regulators identified in this context include the STAT1:STAT2 heterodimer and interferon regulatory factor family proteins. The presented protocol is designed to be accessible to individuals with basic computer literacy. Understanding the interactions between TFs and TFBSs provides insights into the complex transcriptional regulatory networks that govern gene expression, with broad implications for several fields such as developmental biology, immunology, and disease research.

Epigenomic analysis of the myometrium during late implantation revealed regulatory elements in genes related to the cellular zinc homeostasis pathway in pigs

The myometrium, composed of the inner circular muscle (CM) and outer longitudinal muscle (LM), is crucial in establishing and maintaining early pregnancy. However, the molecular mechanisms involved are not well understood. In this study, we identified the transcriptomic features of the CM and LM collected from the mesometrial (M) and anti-mesometrial (AM) sides of the pig uterus on day 18 of pregnancy during the placentation initiation phase. Some genes in the cellular zinc ion level regulatory pathways (MT-1A, MT-1D, MT-2B, SLC30A2, and SLC39A2) were spatially and highly enriched in uterine CM at the mesometrial side. In addition, the histone modification profiles of H3K27ac and H3K4me3 in uterine CM and LM collected from the mesometrial side were characterized. Genomic regions associated with the expression of genes regulating the cellular zinc ion level were detected. Moreover, six highly linked variants in the H3K27ac-enriched region of the pig SLC30A2 gene were identified and found to be significantly associated with the total number born at the second parity (P < 0.05). In conclusion, the genes in the pathways of cellular zinc homeostasis and their regulatory elements identified have implications for pig reproduction trait improvement and warrant further investigations.

Differential usage of DNA modifications in neurons, astrocytes, and microglia

BACKGROUND

Cellular identity is determined partly by cell type-specific epigenomic profiles that regulate gene expression. In neuroscience, there is a pressing need to isolate and characterize the epigenomes of specific CNS cell types in health and disease. In this study, we developed an in vivo tagging mouse model (Camk2a-NuTRAP) for paired isolation of neuronal DNA and RNA without cell sorting and then used this model to assess epigenomic regulation, DNA modifications in particular, of gene expression between neurons and glia.

RESULTS

After validating the cell-specificity of the Camk2a-NuTRAP model, we performed TRAP-RNA-Seq and INTACT-whole genome oxidative bisulfite sequencing (WGoxBS) to assess the neuronal translatome and epigenome in the hippocampus of young mice (4 months old). WGoxBS findings were validated with enzymatic methyl-Seq (EM-Seq) and nanopore sequencing. Comparing neuronal data to microglial and astrocytic data from NuTRAP models, microglia had the highest global mCG levels followed by astrocytes and then neurons, with the opposite pattern observed for hmCG and mCH. Differentially modified regions between cell types were predominantly found within gene bodies and distal intergenic regions, rather than proximal promoters. Across cell types there was a negative correlation between DNA modifications (mCG, mCH, hmCG) and gene expression at proximal promoters. In contrast, a negative correlation of gene body mCG and a positive relationship between distal promoter and gene body hmCG with gene expression was observed. Furthermore, we identified a neuron-specific inverse relationship between mCH and gene expression across promoter and gene body regions.

CONCLUSIONS

Neurons, astrocytes, and microglia demonstrate different genome-wide levels of mCG, hmCG, and mCH that are reproducible across analytical methods. However, modification-gene expression relationships are conserved across cell types. Enrichment of differential modifications across cell types in gene bodies and distal regulatory elements, but not proximal promoters, highlights epigenomic patterning in these regions as potentially greater determinants of cell identity. These findings also demonstrate the importance of differentiating between mC and hmC in neuroepigenomic analyses, as up to 30% of what is conventionally interpreted as mCG can be hmCG, which often has a different relationship to gene expression than mCG.

A novel microdeletion of 517 kb downstream of the PAX6 gene in a Chinese family with congenital aniridia

BACKGROUND

To identify the disease-causing gene in a Chinese family affected with congenital aniridia.

METHODS

Patients underwent systematic ophthalmic examinations such as anterior segment photography, fundus photography, optical coherence tomography, and fundus fluorescein angiography. The proband was screened for pathogenic variants by whole exome sequencing (WES) and copy number variant (CNV) analysis. Real-time quantitative PCR (RT-qPCR) was applied to confirm the CNV results. Breakpoints were identified by long-range PCR followed by Sanger sequencing.

RESULTS

All seven members of this Chinese family, including four patients and three normal individuals, were recruited for this study. All patients showed bilateral congenital aniridia with nystagmus, except the son of the proband, who presented with bilateral partial coloboma of the iris. A novel heterozygous deletion (chr11:31,139,019-31,655,997) containing the 3' regulatory enhancers of the PAX6 gene was detected in this family. We also reviewed the reported microdeletions downstream of PAX6 in patients with aniridia.

CONCLUSIONS

We identified a novel microdeletion, 517 kb in size located about 133 kb downstream of the PAX6 gene, responsible for congenital aniridia in this Chinese family, which expands the spectrum of aniridia-associated mutations in PAX6.

Combined loss of CDH1 and downstream regulatory sequences drive early-onset diffuse gastric cancer and increase penetrance of hereditary diffuse gastric cancer

BACKGROUND

Germline CDH1 pathogenic or likely pathogenic variants cause hereditary diffuse gastric cancer (HDGC). Once a genetic cause is identified, stomachs' and breasts' surveillance and/or prophylactic surgery is offered to asymptomatic CDH1 carriers, which is life-saving. Herein, we characterized an inherited mechanism responsible for extremely early-onset gastric cancer and atypical HDGC high penetrance.

METHODS

Whole-exome sequencing (WES) re-analysis was performed in an unsolved HDGC family. Accessible chromatin and CDH1 promoter interactors were evaluated in normal stomach by ATAC-seq and 4C-seq, and functional analysis was performed using CRISPR-Cas9, RNA-seq and pathway analysis.

RESULTS

We identified a germline heterozygous 23 Kb CDH1-TANGO6 deletion in a family with eight diffuse gastric cancers, six before age 30. Atypical HDGC high penetrance and young cancer-onset argued towards a role for the deleted region downstream of CDH1, which we proved to present accessible chromatin, and CDH1 promoter interactors in normal stomach. CRISPR-Cas9 edited cells mimicking the CDH1-TANGO6 deletion display the strongest CDH1 mRNA downregulation, more impacted adhesion-associated, type-I interferon immune-associated and oncogenic signalling pathways, compared to wild-type or CDH1-deleted cells. This finding solved an 18-year family odyssey and engaged carrier family members in a cancer prevention pathway of care.

CONCLUSION

In this work, we demonstrated that regulatory elements lying down-stream of CDH1 are part of a chromatin network that control CDH1 expression and influence cell transcriptome and associated signalling pathways, likely explaining high disease penetrance and very young cancer-onset. This study highlights the importance of incorporating scientific-technological updates and clinical guidelines in routine diagnosis, given their impact in timely genetic diagnosis and disease prevention.

OGRE: calculate, visualize, and analyze overlap between genomic input regions and public annotations

BACKGROUND

Modern genome sequencing leads to an ever-growing collection of genomic annotations. Combining these elements with a set of input regions (e.g. genes) would yield new insights in genomic associations, such as those involved in gene regulation. The required data are scattered across different databases making a manual approach tiresome, unpractical, and prone to error. Semi-automatic approaches require programming skills in data parsing, processing, overlap calculation, and visualization, which most biomedical researchers lack. Our aim was to develop an automated tool providing all necessary algorithms, benefiting both bioinformaticians and researchers without bioinformatic training.

RESULTS

We developed overlapping annotated genomic regions (OGRE) as a comprehensive tool to associate and visualize input regions with genomic annotations. It does so by parsing regions of interest, mining publicly available annotations, and calculating possible overlaps between them. The user can thus identify location, type, and number of associated regulatory elements. Results are presented as easy to understand visualizations and result tables. We applied OGRE to recent studies and could show high reproducibility and potential new insights. To demonstrate OGRE's performance in terms of running time and output, we have conducted a benchmark and compared its features with similar tools.

CONCLUSIONS

OGRE's functions and built-in annotations can be applied as a downstream overlap association step, which is compatible with most genomic sequencing outputs, and can thus enrich pre-existing analyses pipelines. Compared to similar tools, OGRE shows competitive performance, offers additional features, and has been successfully applied to two recent studies. Overall, OGRE addresses the lack of tools for automatic analysis, local genomic overlap calculation, and visualization by providing an easy to use, end-to-end solution for both biologists and computational scientists.

Optimization of mammalian expression vector by cis-regulatory element combinations

The regulation of gene expression in mammalian cells by combining various cis-regulatory features has rarely been discussed. In this study, we constructed expression vectors containing various combinations of regulatory elements to examine the regulation of gene expression by different combinations of cis-regulatory elements. The effects of four promoters (CMV promoter, PGK promoter, Polr2a promoter, and EF-1α core promoter), two enhancers (CMV enhancer and SV40 enhancer), two introns (EF-1α intron A and hybrid intron), two terminators (CYC1 terminator and TEF terminator), and their different combinations on downstream gene expression were compared in various mammalian cells using fluorescence microscopy to observe fluorescence, quantitative real-time PCR (qRT-PCR), and western blot. The receptor binding domain (RBD) sequence from severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2) spike protein was used to replace the eGFP sequence in the expression vector and the RBD expression was detected by qRT-PCR and western blot. The results showed that protein expression can be regulated by optimizing the combination of cis-acting elements. The vector with the CMV enhancer, EF-1α core promoter, and TEF terminator was found to express approximately threefold higher eGFP than the unmodified vector in different animal cells, as well as 2.63-fold higher recombinant RBD protein than the original vector in HEK-293T cells. Moreover, we suggest that combinations of multiple regulatory elements capable of regulating gene expression do not necessarily exhibit synergistic effects to enhance expression further. Overall, our findings provide insights into biological applications that require the regulation of gene expression and will help to optimize expression vectors for biosynthesis and other fields. Additionally, we provide valuable insights into the production of RBD proteins, which may aid in developing reagents for diagnosis and treatment during the COVID-19 pandemic.

Relationships between some transcription factors and concordantly expressed drought stress-related genes in bread wheat

The challenge of climate change makes it mandatory to improve tolerance to drought stress in bread wheat (Triticum aestivum) via biotechnological approaches. Drought stress experiment was conducted followed by RNA-Seq analysis for leaves of two wheat cultivars namely Giza 168 and Gemmiza 10 with contrasting genotypes. Expression patterns of the regulated stress-related genes and concordantly expressed TFs were detected, then, validated via qPCR for two loss-of-function mutants in Arabidopsis background harboring mutated genes analogue to those in wheat. Drought-stress related genes were searched for concordantly expressed TFs and a total of eight TFs were shown to coexpress with 14 stress-related genes. Among these genes, one TF belongs to the zinc finger protein CONSTANS family and proved via qPCR to drive expression of a gene encoding a speculative TF namely zinc transporter 3-like and two other stress related genes encoding tryptophan synthase alpha chain and asparagine synthetase. Known functions of the two TFs under drought stress complement those of the two concordantly expressed stress-related genes, thus, it is likely that they are related. This study highlights the possibility to utilize metabolic engineering approaches to decipher and incorporate existing regulatory frameworks under drought stress in future breeding programs of bread wheat.

ALS is imprinted in the chromatin accessibility of blood cells

Amyotrophic Lateral Sclerosis (ALS) is a complex and incurable neurodegenerative disorder in which genetic and epigenetic factors contribute to the pathogenesis of all forms of ALS. The interplay of genetic predisposition and environmental footprints generates epigenetic signatures in the cells of affected tissues, which then alter transcriptional programs. Epigenetic modifications that arise from genetic predisposition and systemic environmental footprints should in theory be detectable not only in affected CNS tissue but also in the periphery. Here, we identify an ALS-associated epigenetic signature ('epiChromALS') by chromatin accessibility analysis of blood cells of ALS patients. In contrast to the blood transcriptome signature, epiChromALS includes also genes that are not expressed in blood cells; it is enriched in CNS neuronal pathways and it is present in the ALS motor cortex. By combining simultaneous ATAC-seq and RNA-seq with single-cell sequencing in PBMCs and motor cortex from ALS patients, we demonstrate that epigenetic changes associated with the neurodegenerative disease can be found in the periphery, thus strongly suggesting a mechanistic link between the epigenetic regulation and disease pathogenesis.

(Tn5-)FISH-based imaging in the era of 3D/spatial genomics

3D genomics mainly focuses on the 3D position of single genes at the cell level, while spatial genomics focuses more on the tissue level. In this exciting new era of 3D/spatial genomics, half-century old FISH and its derivative methods, including Tn5-FISH, play important roles. In this review, we introduce the Tn5-FISH we developed recently, and present six different applications published by our collaborators and us, based on (Tn5-)FISH, which can be either general BAC clone-based FISH or Tn5-FISH. In these interesting cases, (Tn5-)FISH demonstrated its vigorous ability of targeting sub-chromosomal structures across different diseases and cell lines (leukemia, mESCs (mouse embryonic stem cells), and differentiation cell lines). Serving as an effective tool to image genomic structures at the kilobase level, Tn5-FISH holds great potential to detect chromosomal structures in a high-throughput manner, thus bringing the dawn for new discoveries in the great era of 3D/spatial genomics.

Mini-review: Gene regulatory network benefits from three-dimensional chromatin conformation and structural biology

Gene regulatory networks are now at the forefront of precision biology, which can help researchers better understand how genes and regulatory elements interact to control cellular gene expression, offering a more promising molecular mechanism in biological research. Interactions between the genes and regulatory elements involve different promoters, enhancers, transcription factors, silencers, insulators, and long-range regulatory elements, which occur at a ∼10 µm nucleus in a spatiotemporal manner. In this way, three-dimensional chromatin conformation and structural biology are critical for interpreting the biological effects and the gene regulatory networks. In the review, we have briefly summarized the latest processes in three-dimensional chromatin conformation, microscopic imaging, and bioinformatics, and we have presented the outlook and future directions for these three aspects.

Ensemble learning based assessment of the role of transcription factors in gene expression

Cancer cells are formed when the associated, active genes fail to function the way they are meant to function. Multiple genes collectively control cell growth by activating a proper set of genes. Regulation of gene expression is controlled through the combined effort of multiple regulatory elements. Transcription of each gene is affected differently according to the combinatorial patterns of regulatory elements bound in the nearby regions. Identifying and analysing such patterns will give a better insight into the cell function. The main focus of this study is on developing a computational model to predict the functional role of transcriptional factors residing between divergent gene pairs. Acute Myeloid Leukaemia (AML) gene expression data from GEO and the two TFs EP300 and CTCF binding data calibrated in k562 cell line from ENCODE consortium are taken as a case study.

In Embryo Gene Reporter Assays for Evaluation of Cis-Regulatory Regions

Gene expression reporter assays measure the relevance of cis-regulatory elements and DNA-binding proteins in modulating transcriptional activity. Commonly, they are performed in cell lines. However, regulation of transcriptional activity during development is complex and dynamic, and not many cell lines reproduce the embryonic conditions. Thus, conclusions derived from cell line data provide limited information about embryonic development. On the other hand, one of the major hurdles for embryonic assays is delivering reporter plasmids in a tissue-specific manner. In this sense, the chick embryo is a good model system to perform these assays. Electroporation of chick embryos provides temporal and spatially controlled plasmid delivery. Further, it is a well-established, easy, and an economical procedure. Here, we describe in detail how to measure in the chick neural tube (1) enhancer activity with GFP, (2) enhancer activity with luciferase, and (3) 3'UTR activity with luciferase.

Characterizing chromatin interactions of regulatory elements and nucleosome positions, using Hi-C, Micro-C, and promoter capture Micro-C

BACKGROUND

Regulatory elements such as promoters, enhancers, and insulators interact each other to mediate molecular processes. To capture chromatin interactions of regulatory elements, 3C-derived methods such as Hi-C and Micro-C are developed. Here, we generated and analyzed Hi-C, Micro-C, and promoter capture Micro-C datasets with different sequencing depths to study chromatin interactions of regulatory elements and nucleosome positions in human prostate cancer cells.

RESULTS

Compared to Hi-C, Micro-C identifies more high-resolution loops, including ones around structural variants. By evaluating the effect of sequencing depth, we revealed that more than 2 billion reads of Micro-C are needed to detect chromatin interactions at 1 kb resolution. Moreover, we found that deep-sequencing identifies additional long-range loops that are longer than 1 Mb in distance. Furthermore, we found that more than 50% of the loops are involved in insulators while less than 10% of the loops are promoter-enhancer loops. To comprehensively capture chromatin interactions that promoters are involved in, we performed promoter capture Micro-C. Promoter capture Micro-C identifies loops near promoters with a lower amount of sequencing reads. Sequencing of 160 million reads of promoter capture Micro-C resulted in reaching a plateau of identifying loops. However, there was still a subset of promoters that are not involved in loops even after deep-sequencing. By integrating Micro-C with NOMe-seq and ChIP-seq, we found that active promoters involved in loops have a more accessible region with lower levels of DNA methylation and more highly phased nucleosomes, compared to active promoters that are not involved in loops.

CONCLUSION

We determined the required sequencing depth for Micro-C and promoter capture Micro-C to generate high-resolution chromatin interaction maps and loops. We also investigated the effect of sequencing coverage of Hi-C, Micro-C, and promoter capture Micro-C on detecting chromatin loops. Our analyses suggest the presence of distinct regulatory element groups, which are differently involved in nucleosome positions and chromatin interactions. This study does not only provide valuable insights on understanding chromatin interactions of regulatory elements, but also present guidelines for designing research projects on chromatin interactions among regulatory elements.

Comparative epigenomics reveals the impact of ruminant-specific regulatory elements on complex traits

BACKGROUND

Insights into the genetic basis of complex traits and disease in both human and livestock species have been achieved over the past decade through detection of genetic variants in genome-wide association studies (GWAS). A majority of such variants were found located in noncoding genomic regions, and though the involvement of numerous regulatory elements (REs) has been predicted across multiple tissues in domesticated animals, their evolutionary conservation and effects on complex traits have not been fully elucidated, particularly in ruminants. Here, we systematically analyzed 137 epigenomic and transcriptomic datasets of six mammals, including cattle, sheep, goats, pigs, mice, and humans, and then integrated them with large-scale GWAS of complex traits.

RESULTS

Using 40 ChIP-seq datasets of H3K4me3 and H3K27ac, we detected 68,479, 58,562, 63,273, 97,244, 111,881, and 87,049 REs in the liver of cattle, sheep, goats, pigs, humans and mice, respectively. We then systematically characterized the dynamic functional landscapes of these REs by integrating multi-omics datasets, including gene expression, chromatin accessibility, and DNA methylation. We identified a core set (n = 6359) of ruminant-specific REs that are involved in liver development, metabolism, and immune processes. Genes with more complex cis-REs exhibited higher gene expression levels and stronger conservation across species. Furthermore, we integrated expression quantitative trait loci (eQTLs) and GWAS from 44 and 52 complex traits/diseases in cattle and humans, respectively. These results demonstrated that REs with different degrees of evolutionary conservation across species exhibited distinct enrichments for GWAS signals of complex traits.

CONCLUSIONS

We systematically annotated genome-wide functional REs in liver across six mammals and demonstrated the evolution of REs and their associations with transcriptional output and conservation. Detecting lineage-specific REs allows us to decipher the evolutionary and genetic basis of complex phenotypes in livestock and humans, which may benefit the discovery of potential biomedical models for functional variants and genes of specific human diseases.

The Neuroepigenome: Implications of Chemical and Physical Modifications of Genomic DNA in Schizophrenia

Schizophrenia is a chronic mental illness with a substantial genetic component. To unfold the complex etiology of schizophrenia, it is important to understand the interplay between genetic and nongenetic factors. Genetic factors involve variation in the DNA sequences of protein-coding genes, which directly contribute to phenotypic traits, and variation in noncoding sequences, which comprise 98% of the genome and contain DNA elements known to play a role in regulating gene expression. The epigenome refers to the chemical modifications on both DNA and the structural proteins that package DNA into the nucleus, which together regulate gene expression in specific cell types, conditions, and developmental stages. The dynamic nature of the epigenome makes it an ideal tool to investigate the relationship between inherited genetic mutations associated with schizophrenia and altered gene regulation throughout the course of brain development. In this review, we focus on the current understanding of the role of epigenetic marks and their three-dimensional nuclear organization in the developmental trajectory of distinct brain cell types to decipher the complex gene regulatory mechanisms that are disrupted in schizophrenia.

Development of dual reporter vector system for estimating translational activity of regulatory elements

BACKGROUND

For the needs of modern biotechnology, a quantitative approach to the control of regulatory elements at all stages of gene expression has long become indispensable. Such a control regime is impossible without a quantitative analysis of the role of each regulatory element or pattern used. Therefore, it seems important to modify and develop the accuracy, reproducibility, and availability of methods for quantifying the contribution of each regulatory code to the implementation of genetic information.

RESULTS

A new vector system for transient expression in plants is described; this system is intended for quantitative analysis of the contribution of regulatory elements to transcription and translation efficiencies. The proposed vector comprises two expression cassettes carrying reporter genes (of the Clostridium thermocellum thermostable lichenase and E. coli β-glucuronidase) under the control of different promoters. Herewith we also propose a new method for quantification of the effect of tested regulatory elements on expression, which relies on assessment of the enzyme activities of reporter proteins taking into account the transcription of their genes.

CONCLUSIONS

In our view, this approach makes it possible to precisely determine the amounts of reporter proteins and their transcripts at all stages of expression. The efficiency of the proposed system has been validated by the analysis of the roles of known translation enhancers at the stages of transcription and translation.

Genetic variants of the hypoxia-inducible factor 3 alpha subunit (Hif3a) gene in the Fat and Lean mouse selection lines

BACKGROUND

Adipose tissue hypoxia and members of the hypoxia-inducible factor alpha (HIFA) are involved in development of obesity. However, the mechanism and functions of HIF3A, one of three HIFA paralogs, in fat deposition have not been sufficiently studied.

METHODS AND RESULTS

In the present study, we investigated whether Hif3a sequence variants are associated with divergent fat deposition in mouse selection lines for fatness and leanness. Sequencing and RFLP were used to analyse sequence variants within Hif3a. To identify candidate regulatory variants, we performed literature screening and used databases and bioinformatics tools like Ensembl, MethPrimer, TargetScanMouse, miRDB, PolyAsite, RISE, LncRRIsearch, RNAfold, PredictProtein, CAIcal, and switches.ELM Resource. There are 90 sequence variants in Hif3a between the two mouse lines. While most Fat line variants locate within intronic regions, Lean line variants are mainly in 3' UTR. We constructed a map of Hif3a potential regulatory regions and identified 39 regulatory variants by integrating data on constrained and regulatory elements, CpGs, and miRNAs and lncRNAs binding sites. Moreover, 3' UTR and two exonic variants may influence mRNA stability, translation rate and protein functionality. We propose as priority candidates for further functional studies a missense (rs37398126) and synonymous (rs37739792) variants, and intronic (rs47471302) variant that overlap conserved element in promoter region and predicted lncRNAs binding site.

CONCLUSION

The results indicate a potential involvement of Hif3a in fat deposition. Additionally, approach used in the present study may serve as a general guideline for constructing an integrative gene map for prioritizing candidate gene variants with phenotypic effects.

Mapping Chromatin Accessibility in Human Naïve Pluripotent Stem Cells Using ATAC-Seq

Regulatory elements, such as promoters and enhancers, typically show reduced nucleosome occupancy, which is a feature that is commonly referred to as "open chromatin". The distribution of open chromatin sites can provide important clues about the transcription factors and regulatory networks that drive gene expression in a given cell. Assay for Transposase-Accessible Chromatin using sequencing (ATAC-seq) is a rapid and robust method for mapping open chromatin sites. ATAC-seq data can also discern the binding sites of nucleosomes and transcription factors. In this chapter, we describe how to produce and assess the quality of ATAC-seq libraries that are generated from naïve human pluripotent stem cells.

A cis-regulatory-directed pipeline for the identification of genes involved in cardiac development and disease

BACKGROUND

Congenital heart diseases are the major cause of death in newborns, but the genetic etiology of this developmental disorder is not fully known. The conventional approach to identify the disease-causing genes focuses on screening genes that display heart-specific expression during development. However, this approach would have discounted genes that are expressed widely in other tissues but may play critical roles in heart development.

RESULTS

We report an efficient pipeline of genome-wide gene discovery based on the identification of a cardiac-specific cis-regulatory element signature that points to candidate genes involved in heart development and congenital heart disease. With this pipeline, we retrieve 76% of the known cardiac developmental genes and predict 35 novel genes that previously had no known connectivity to heart development. Functional validation of these novel cardiac genes by RNAi-mediated knockdown of the conserved orthologs in Drosophila cardiac tissue reveals that disrupting the activity of 71% of these genes leads to adult mortality. Among these genes, RpL14, RpS24, and Rpn8 are associated with heart phenotypes.

CONCLUSIONS

Our pipeline has enabled the discovery of novel genes with roles in heart development. This workflow, which relies on screening for non-coding cis-regulatory signatures, is amenable for identifying developmental and disease genes for an organ without constraining to genes that are expressed exclusively in the organ of interest.

Gene expression profiling and functional analysis of ventricular tissues from murine transverse aortic constriction

BACKGROUND

Transverse aortic constriction (TAC) model is widely used to study pressure overload-induced cardiac remodeling. However, the conserved transcriptional features of TAC model and the underlying regulatory mechanisms remain unclear.

METHODS

In this study, we screened out the high-quality microarray data for ventricular tissue from murine TAC model. The transcriptional changes in ventricular tissue were analyzed by identifying the common differently expressed genes (DEGs) and enriched gene sets. We also analyzed the protein-protein interaction and mRNA-mRNA association of DEGs. Furthermore, the potential regulatory elements of the DEGs were explored through comparative analysis between mouse and human.

RESULTS

265 common DEGs and 45 enriched canonical pathways were identified in murine TAC model. 201 DEGs had the protein-protein interaction, whereas 96 DEGs had mRNA-mRNA association. 99 transcription factor (TF)-mRNA and 2997 microRNA (miRNA)-mRNA regulatory relationships were retrieved.

CONCLUSIONS

In pressure overload-induced cardiac remodeling, inflammation, fibrosis, metabolic remodeling and hypoxia were significant features. Approaches to intervene these phenomena may have therapeutic values. TFs and miRNAs are important regulator elements of DEGs in both mouse and human. Examination of miRNAs is a promising tool to detect the occurrence of pressure overload-induced cardiac remodeling in patients.

Cis-Regulatory Control of Mammalian Sex Determination

Sex determination is the process by which an initial bipotential gonad adopts either a testicular or ovarian cell fate. The inability to properly complete this process leads to a group of developmental disorders classified as disorders of sex development (DSD). To date, dozens of genes were shown to play roles in mammalian sex determination, and mutations in these genes can cause DSD in humans or gonadal sex reversal/dysfunction in mice. However, exome sequencing currently provides genetic diagnosis for only less than half of DSD patients. This points towards a major role for the non-coding genome during sex determination. In this review, we highlight recent advances in our understanding of non-coding, cis-acting gene regulatory elements and discuss how they may control transcriptional programmes that underpin sex determination in the context of the 3-dimensional folding of chromatin. As a paradigm, we focus on the Sox9 gene, a prominent pro-male factor and one of the most extensively studied genes in gonadal cell fate determination.

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

Background

Mammalian 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.

Results

Here, 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.

Conclusions

Models 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.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13059-021-02503-y.

Regulatory regions in natural transposable element insertions drive interindividual differences in response to immune challenges in Drosophila

Background

Variation in gene expression underlies interindividual variability in relevant traits including immune response. However, the genetic variation responsible for these gene expression changes remains largely unknown. Among the non-coding variants that could be relevant, transposable element insertions are promising candidates as they have been shown to be a rich and diverse source of cis-regulatory elements.

Results

In this work, we use a population genetics approach to identify transposable element insertions likely to increase the tolerance of Drosophila melanogaster to bacterial infection by affecting the expression of immune-related genes. We identify 12 insertions associated with allele-specific expression changes in immune-related genes. We experimentally validate three of these insertions including one likely to be acting as a silencer, one as an enhancer, and one with a dual role as enhancer and promoter. The direction in the change of gene expression associated with the presence of several of these insertions is consistent with an increased survival to infection. Indeed, for one of the insertions, we show that this is the case by analyzing both natural populations and CRISPR/Cas9 mutants in which the insertion is deleted from its native genomic context.

Conclusions

We show that transposable elements contribute to gene expression variation in response to infection in D. melanogaster and that this variation is likely to affect their survival capacity. Because the role of transposable elements as regulatory elements is not restricted to Drosophila, transposable elements are likely to play a role in immune response in other organisms as well.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13059-021-02471-3.

Molecular and computational approaches to map regulatory elements in 3D chromatin structure

Epigenetic marks do not change the sequence of DNA but affect gene expression in a cell-type specific manner by altering the activities of regulatory elements. Development of new molecular biology assays, sequencing technologies, and computational approaches enables us to profile the human epigenome in three-dimensional structure genome-wide. Here we describe various molecular biology techniques and bioinformatic tools that have been developed to measure the activities of regulatory elements and their chromatin interactions. Moreover, we list currently available three-dimensional epigenomic data sets that are generated in various human cell types and tissues to assist in the design and analysis of research projects.

Differential transcriptional regulation of the NANOG gene in chicken primordial germ cells and embryonic stem cells

Background

NANOG is a core transcription factor (TF) in embryonic stem cells (ESCs) and primordial germ cells (PGCs). Regulation of the NANOG gene by TFs, epigenetic factors, and autoregulatory factors is well characterized in ESCs, and transcriptional regulation of NANOG is well established in these cells. Although NANOG plays a key role in germ cells, the molecular mechanism underlying its transcriptional regulation in PGCs has not been studied. Therefore, we investigated the mechanism that regulates transcription of the chicken NANOG (cNANOG) gene in PGCs and ESCs.

Results

We first identified the transcription start site of cNANOG by 5′-rapid amplification of cDNA ends PCR analysis. Then, we measured the promoter activity of various 5′ flanking regions of cNANOG in chicken PGCs and ESCs using the luciferase reporter assay. cNANOG expression required transcriptional regulatory elements, which were positively regulated by POU5F3 (OCT4) and SOX2 and negatively regulated by TP53 in PGCs. The proximal region of the cNANOG promoter contains a positive transcriptional regulatory element (CCAAT/enhancer-binding protein (CEBP)-binding site) in ESCs. Furthermore, small interfering RNA-mediated knockdown demonstrated that POU5F3, SOX2, and CEBP played a role in cell type-specific transcription of cNANOG.

Conclusions

We show for the first time that different trans-regulatory elements control transcription of cNANOG in a cell type-specific manner. This finding might help to elucidate the mechanism that regulates cNANOG expression in PGCs and ESCs.

Self-organizing maps with variable neighborhoods facilitate learning of chromatin accessibility signal shapes associated with regulatory elements

BACKGROUND

Assigning chromatin states genome-wide (e.g. promoters, enhancers, etc.) is commonly performed to improve functional interpretation of these states. However, computational methods to assign chromatin state suffer from the following drawbacks: they typically require data from multiple assays, which may not be practically feasible to obtain, and they depend on peak calling algorithms, which require careful parameterization and often exclude the majority of the genome. To address these drawbacks, we propose a novel learning technique built upon the Self-Organizing Map (SOM), Self-Organizing Map with Variable Neighborhoods (SOM-VN), to learn a set of representative shapes from a single, genome-wide, chromatin accessibility dataset to associate with a chromatin state assignment in which a particular RE is prevalent. These shapes can then be used to assign chromatin state using our workflow.

RESULTS

We validate the performance of the SOM-VN workflow on 14 different samples of varying quality, namely one assay each of A549 and GM12878 cell lines and two each of H1 and HeLa cell lines, primary B-cells, and brain, heart, and stomach tissue. We show that SOM-VN learns shapes that are (1) non-random, (2) associated with known chromatin states, (3) generalizable across sets of chromosomes, and (4) associated with magnitude and multimodality. We compare the accuracy of SOM-VN chromatin states against the Clustering Aggregation Tool (CAGT), an unsupervised method that learns chromatin accessibility signal shapes but does not associate these shapes with REs, and we show that overall precision and recall is increased when learning shapes using SOM-VN as compared to CAGT. We further compare enhancer state assignments from SOM-VN in signals above a set threshold to enhancer state assignments from Predicting Enhancers from ATAC-seq Data (PEAS), a deep learning method that assigns enhancer chromatin states to peaks. We show that the precision-recall area under the curve for the assignment of enhancer states is comparable to PEAS.

CONCLUSIONS

Our work shows that the SOM-VN workflow can learn relationships between REs and chromatin accessibility signal shape, which is an important step toward the goal of assigning and comparing enhancer state across multiple experiments and phenotypic states.

Interrogating the Accessible Chromatin Landscape of Eukaryote Genomes Using ATAC-seq

The ATAC-seq assay has emerged as the most useful, versatile, and widely adaptable method for profiling accessible chromatin regions and tracking the activity of cis-regulatory elements (cREs) in eukaryotes. Thanks to its great utility, it is now being applied to map active chromatin in the context of a very wide diversity of biological systems and questions. In the course of these studies, considerable experience working with ATAC-seq data has accumulated and a standard set of computational tasks that need to be carried for most ATAC-seq analyses has emerged. Here, we review and provide examples of common such analytical procedures (including data processing, quality control, peak calling, identifying differentially accessible open chromatin regions, and variable transcription factor (TF) motif accessibility) and discuss recommended optimal practices.

Insights into gene regulation: From regulatory genomic elements to DNA-protein and protein-protein interactions

Transcription is orchestrated by non-coding regulatory elements embedded in chromatin, which exist within the larger context of chromosome topology. Here, we review recent insights into the functions of non-coding regulatory elements and their protein interactors during transcription control. A picture emerges in which the topological environment constraints enhancer-promoter interactions and specific enhancer-bound proteins with distinct promoter-compatibilities refine target promoter choice. Such compatibilities are encoded within the sequences of enhancers and promoters and realized by diverse transcription factors and cofactors with distinct biochemical activities. An emerging property of transcription factors and cofactors is the formation of nuclear microenvironments or membraneless compartments that can have properties of phase-separated liquids. These environments are able to selectively enrich certain proteins and small molecules over others. Further investigation into the interaction of transcriptional regulators with themselves and regulatory DNA elements will help reveal the complexities of gene regulation within the context of the nucleus.

Prediction of putative regulatory elements in the subgenomic promoters of cucumber green mottle mosaic virus and their interactions with the RNA dependent RNA polymerase domain

Characterization of the subgenomic RNA (sgRNA) promoter of many plant viruses is important to understand the expression of downstream genes and also to configure their genome into a suitable virus gene-vector system. Cucumber green mottle mosaic virus (CGMMV, genus Tobamovirus) is one of the RNA viruses, which is extensively being exploited as the suitable gene silencing and protein expression vector. Even though, characters of the sgRNA promoters (SGPs) of CGMMV are yet to be addressed. In the present study, we predicted the SGP for the movement protein (MP) and coat protein (CP) of CGMMV. Further, we identified the key regulatory elements in the SGP regions of MP and CP, and their interactions with the core RNA dependent RNA polymerase (RdRp) domain of CGMMV was deciphered. The modeled structure of core RdRp contains two palm (1-41 aa, and 63-109 aa), one finger (42-62 aa) subdomains with three conserved RdRp motifs that played important role in binding to the SGP nucleic acids. RdRp strongly preferred the double helix form of the stem region in the stem and loop (SL) structures, and the internal bulge elements. In MP-SGP, a total of six elements was identified; of them, the affinity of binding to - 26 nt to - 17 nt site (CGCGGAAAAG) was higher through the formation of strong hydrogen bonds with LYS16, TYR17, LYS19, SER20, etc. of the motif A in the palm subdomain of RdRp. Similar strong interactions were noticed in the internal bulge (CAACUUU) located at + 33 to + 39 nt adjacent to the translation start site (TLSS) (+ 1). These could be proposed as the putative core promoter elements in MP-SGP. Likewise, total five elements were predicted within - 114 nt to + 144 nt region of CP-SGP with respect to CP-TLSS. Of them, RdRp preferred to bind at the small hairpin located at - 60 nt to - 43 nt (UUGGAGGUUUAGCCUCCA) in the upstream region, and at the complex duplex structure spanning between + 99 and + 114 nt in the downstream region, thus indicating the distribution of core promoter within - 60 nt to + 114 nt region of CP-SGP with respect to TLSS (+ 1) of the CP; whereas, the - 114 nt to + 144 nt region of CP-SGP might be necessary for the full activity of the CP-SGP. Our in silico prediction certifies the gravity of these nucleotide stretches as the RNA regulatory elements and identifies their potentiality for binding with of palm and finger sub-domain of RdRp. Identification of such elements will be helpful to anticipate the critical length of the SGPs. Our finding will not only be helpful to delineate the SGPs of CGMMV but also their subsequent application in the efficient construction of virus gene-vector for the expression of foreign protein in plant.

Identification of atheroprone shear stress responsive regulatory elements in endothelial cells

AIMS

Oscillatory shear stress (OSS) is an atheroprone haemodynamic force that occurs in areas of vessel irregularities and is implicated in the pathogenesis of atherosclerosis. Changes in signalling and transcriptional programme in response to OSS have been vigorously studied; however, the underlying changes in the chromatin landscape controlling transcription remain to be elucidated. Here, we investigated the changes in the regulatory element (RE) landscape of endothelial cells under atheroprone OSS conditions in an in vitro model.

METHODS AND RESULTS

Analyses of H3K27ac chromatin immunoprecipitation-Seq enrichment and RNA-Seq in primary human umbilical vein endothelial cells 6 h after onset of OSS identified 2806 differential responsive REs and 33 differentially expressed genes compared with control cells kept under static conditions. Furthermore, gene ontology analyses of putative RE-associated genes uncovered enrichment of WNT/HIPPO pathway and cytoskeleton reorganization signatures. Transcription factor (TF) binding motif analysis within RE sequences identified over-representation of ETS, Zinc finger, and activator protein 1 TF families that regulate cell cycle, proliferation, and apoptosis, implicating them in the development of atherosclerosis. Importantly, we confirmed the activation of EGR1 as well as the YAP/TAZ complex early (6 h) after onset of OSS in both cultured human vein and artery endothelial cells and, by undertaking luciferase assays, functionally verified their role in RE activation in response to OSS.

CONCLUSIONS

Based on the identification and verification of specific responsive REs early upon OSS exposure, we propose an expanded mechanism of how OSS might contribute to the development of atherosclerosis.

Identification of integrin β6 gene promoter and analysis of its transcription regulation in colon cancer cells

BACKGROUND

The integrin β6 gene, which is expressed in epithelial cancer, plays a pivotal role in various aspects of cancer progression. The present research for integrin β6 regulation mainly focuses on the post-transcription and translation related regulation mechanism and its role in tumorigenesis. The mechanisms of how the integrin β6 gene is regulated transcriptionally, and the promoter and transcription factors responsible for basic transcription of integrin β6 gene remain unknown.

AIM

To clone and characterize the integrin β6 promoter.

METHODS

Software analysis was used to predict the region of integrin β6 promoter. Luciferase reporter plasmids, which contained the integrin β6 promoter, were constructed. Element deletion analysis was performed to identify the location of core promoter and binding sites for transcription factors.

RESULTS

The regulatory elements for the transcription of the integrin β6 gene were located between -286 and -85 and contained binding sites for transcription factors such as STAT3 and Ets-1.

CONCLUSION

For the first time, we found the region of β6 core promoter and demonstrated the binding sites for transcription factors such as Ets-1 and STAT3, which are important for integrin β6 promoter transcription activity. These findings are important for investigating the mechanism of integrin β6 activation in cancer progression.

Cell fate decisions during the development of the peripheral nervous system in the vertebrate head

Sensory placodes and neural crest cells are among the key cell populations that facilitated the emergence and diversification of vertebrates throughout evolution. Together, they generate the sensory nervous system in the head: both form the cranial sensory ganglia, while placodal cells make major contributions to the sense organs-the eye, ear and olfactory epithelium. Both are instrumental for integrating craniofacial organs and have been key to drive the concentration of sensory structures in the vertebrate head allowing the emergence of active and predatory life forms. Whereas the gene regulatory networks that control neural crest cell development have been studied extensively, the signals and downstream transcriptional events that regulate placode formation and diversity are only beginning to be uncovered. Both cell populations are derived from the embryonic ectoderm, which also generates the central nervous system and the epidermis, and recent evidence suggests that their initial specification involves a common molecular mechanism before definitive neural, neural crest and placodal lineages are established. In this review, we will first discuss the transcriptional networks that pattern the embryonic ectoderm and establish these three cell fates with emphasis on sensory placodes. Second, we will focus on how sensory placode precursors diversify using the specification of otic-epibranchial progenitors and their segregation as an example.

In silico identification of novel transcription factors associated with CYP27B1 transcriptional regulation in LPS-challenged mononuclear phagocytes

Renal and extrarenal production of the active form of vitamin D, 1,25-dihydroxyvitamin D (1,25(OH)2D), is catalyzed by CYP27B1, an enzyme also called 1-α-hydroxylase. The overproduction of 1,25(OH)2D has been described in granulomatous diseases. High circulating concentrations of 1,25(OH)2D can lead to hypercalcemia. The aim of this work was to characterize the transcriptional regulation of CYP27B1 in human mononuclear phagocytes exposed to LPS due to its relevance to understanding the hypercalcemia and ectopic calcifications associated with chronic inflammatory diseases such as tuberculosis and other granulomatous diseases. The human CYP27B1 promoter analysis identified binding sites for published TF, SNPs, novel putative TFBS and conserved sites compared to mice. Then, using microarray data, a meta-analysis was performed to obtain a global view of the gene expression in LPS-challenged dendritic cells, monocytes and macrophages. Finally, two experiments, GSE40885 and time series GSE19765, were analyzed in-depth using differential expression analysis which permitted the identification of TF co-expressed with CYP27B1. This work allowed us to formulate a CYP27B1 transcriptional regulation model for LPS-challenged monocytes/macrophages. The importance of two TF families, NFKB and CEBPB, was confirmed. Data also suggests that PLAGL2 and STAT4 which are novel TF could participate in the CYP27B1 transcriptional regulation in cells exposed to LPS. These TF, in turn, would be interacting with regions that present polymorphisms in the general population which might explain the pathological phenotypes associated with altered vitamin D metabolism.

Somatic synonymous mutations in regulatory elements contribute to the genetic aetiology of melanoma

Background

Non-synonymous mutations altering tumor suppressor genes and oncogenes are widely studied. However, synonymous mutations, which do not alter the protein sequence, are rarely investigated in melanoma genome studies.

Methods

We explored the role of somatic synonymous mutations in melanoma samples from TCGA (The Cancer Genome Atlas). The pathogenic synonymous mutation and neutral synonymous mutation data were used to assess the significance of pathogenic synonymous mutations in melanoma likely to affect genetic regulatory elements using Fisher’s exact test. Poisson distribution probabilities of each gene were used to mine the genes with multiple potential functional synonymous mutations affecting regulatory elements.

Results

Concentrating on five types of genetic regulatory functions, we found that the mutational patterns of pathogenic synonymous mutations are mostly involved in exonic splicing regulators in near-splicing sites or inside DNase I hypersensitivity sites or non-optimal codon. Moreover, the sites of miRNA binding alteration exhibit a significantly lower rate of evolution than other sites. Finally, 12 genes were hit by recurrent potentially functional synonymous mutations, which showed statistical significance in the pathogenic mutations. Among them, nine genes (DNAH5, ADCY8, GRIN2A, KSR2, TECTA, RIMS2, XKR6, MYH1, SCN10A) have been reported to be mutated in melanoma, and other three genes (SLC9A2, CASR, SLC8A3) have a great potential to impact melanoma.

Conclusion

These findings confirm the functional consequences of somatic synonymous mutations in melanoma, emphasizing the significance of research in future studies.

Engineering salinity tolerance in plants: progress and prospects

There is a need to integrate conceptual framework based on the current understanding of salt stress responses with different approaches for manipulating and improving salt tolerance in crop plants. Soil salinity exerts significant constraints on global crop production, posing a serious challenge for plant breeders and biotechnologists. The classical transgenic approach for enhancing salinity tolerance in plants revolves by boosting endogenous defence mechanisms, often via a single-gene approach, and usually involves the enhanced synthesis of compatible osmolytes, antioxidants, polyamines, maintenance of hormone homeostasis, modification of transporters and/or regulatory proteins, including transcription factors and alternative splicing events. Occasionally, genetic manipulation of regulatory proteins or phytohormone levels confers salinity tolerance, but all these may cause undesired reduction in plant growth and/or yields. In this review, we present and evaluate novel and cutting-edge approaches for engineering salt tolerance in crop plants. First, we cover recent findings regarding the importance of regulatory proteins and transporters, and how they can be used to enhance salt tolerance in crop plants. We also evaluate the importance of halobiomes as a reservoir of genes that can be used for engineering salt tolerance in glycophytic crops. Additionally, the role of microRNAs as critical post-transcriptional regulators in plant adaptive responses to salt stress is reviewed and their use for engineering salt-tolerant crop plants is critically assessed. The potentials of alternative splicing mechanisms and targeted gene-editing technologies in understanding plant salt stress responses and developing salt-tolerant crop plants are also discussed.

Capturing functional epigenomes for insight into metabolic diseases

Background

Metabolic diseases such as obesity are known to be driven by both environmental and genetic factors. Although genome-wide association studies of common variants and their impact on complex traits have provided some biological insight into disease etiology, identified genetic variants have been found to contribute only a small proportion to disease heritability, and to map mainly to non-coding regions of the genome. To link variants to function, association studies of cellular traits, such as epigenetic marks, in disease-relevant tissues are commonly applied.

Scope of the review

We review large-scale efforts to generate genome-wide maps of coordinated epigenetic marks and their utility in complex disease dissection with a focus on DNA methylation. We contrast DNA methylation profiling methods and discuss the advantages of using targeted methods for single-base resolution assessments of methylation levels across tissue-specific regulatory regions to deepen our understanding of contributing factors leading to complex diseases.

Major conclusions

Large-scale assessments of DNA methylation patterns in metabolic disease-linked study cohorts have provided insight into the impact of variable epigenetic variants in disease etiology. In-depth profiling of epigenetic marks at regulatory regions, particularly at tissue-specific elements, will be key to dissect the genetic and environmental components contributing to metabolic disease onset and progression.

•

Changes in epigenetic marks have been linked to metabolic disease phenotypes.

•

Disease-linked sites of variable DNA methylation status are enriched in distal regulatory regions of disease-linked tissues.

•

Distal regulatory elements remain underrepresented in popular array-based methylation profiling technologies.

•

Novel next-generation capture methods provide cost-effective solutions to assess the impact of DNA methylation in metabolic diseases specifically at regulatory elements.

•

Improvements in methodologies to account for tissue heterogeneity and causality will be crucial in future epigenome-wide association studies.

Synthetic promoter for efficient and muscle-specific expression of exogenous genes

Synthetic promoters (SPs) have many advantages over their natural counterparts, especially with regard to transcriptional activity and tissue specificity. Here, we report a new strategy to construct SPs for efficient and muscle-specific gene expression. First, 19 nucleic acid motifs classified to 3 kinds of transcriptional regulatory elements were rationally selected. A recombinant promoter library was constructed by randomly assembling these motifs. Second, the transcriptional activities of ~1200 SPs were screened by intramuscular expression of several reporter genes in different cell lines for activity higher than that of the cytomegalovirus (CMV) promoter, with SP-301 finally identified as the strongest. A single intramuscular injection of mice with an SP-301 plasmid expressing mouse growth hormone releasing hormone accelerated mouse growth significantly over 24 days. Third, the muscle specificity of SP-301 was confirmed in transgenic mice. Finally, in comparison with the CMV promoter, SP-301 accelerated translocation and increased the level of plasmid in the nuclei of myoblast cells to a greater extent than in non-muscle cells. Altogether, the study has provided a more rational strategy to construct efficient and tissue-specific promoters, with the promoter SP-301 exhibiting promising potential for establishing an intramuscular gene expression system for therapeutic applications.

Identification of cis-regulatory regions responsible for developmental and hormonal regulation of HbHMGS1 in transgenic Arabidopsis thaliana

OBJECTIVES

3-Hydroxy-3-methylglutaryl coenzyme A (HMG-CoA) synthase (HMGS) is an important enzyme in mevalonate (MVA) pathway of isoprenoid biosynthesis, which regulates the rubber biosynthetic pathway in rubber tree (Hevea brasiliensis) in coordination with HMG-CoA reductase (HMGR). However, little information is available about the regulation of HMGS gene expression. To understand the mechanism controlling the HbHMGS1 gene expression, we characterized the HbHMGS1 promoter sequence in transgenic plants with the β-glucuronidase (GUS) reporter gene.

RESULTS

GUS activity analysis of the transgenic plants showed that the HbHMGS1 promoter is active in all organs of the transgenic Arabidopsis plants during various developmental stages (from 6 to 45-day-old). Deletion of different portions of the upstream HbHMGS1 promoter identified sequences responsible for either positive or negative regulation of the GUS expression. Particularly, the - 454 bp HbHMGS1 promoter resulted in a 2.19-fold increase in promoter activity compared with the CaMV 35S promoter, suggesting that the - 454 bp HbHMGS1 promoter is a super-strong near-constitutive promoter. In addition, a number of promoter regions important for the responsiveness to ethylene, methyl jasmonate (MeJA) and gibberellic acid (GA) were identified.

CONCLUSION

The - 454 bp HbHMGS1 promoter has great application potential in plant transformation studies as an alternative to the CaMV 35S promoter. The HbHMGS1 promoter may play important roles in regulating ethylene-, MeJA- and GA-mediated gene expression. The functional complexity of cis-elements revealed by this study remains to be elucidated.

Navigating the non-coding genome in heart development and Congenital Heart Disease

Congenital Heart Disease (CHD) is characterised by a wide range of cardiac defects, from mild to life-threatening, which occur in babies worldwide. To date, there is no cure to CHD, however, progress in surgery has reduced its mortality allowing children affected by CHD to reach adulthood. In an effort to understand its genetic basis, several studies involving whole-genome sequencing (WGS) of patients with CHD have been undertaken and generated a great wealth of information. The majority of putative causative mutations identified in WGS studies fall into the non-coding part of the genome. Unfortunately, due to the lack of understanding of the function of these non-coding mutations, it is challenging to establish a causal link between the non-coding mutation and the disease. Thus, here we review the state-of-the-art approaches to interpret non-coding mutations in the context of CHD and address the following questions: What are the non-coding sequences important for cardiac function? Which technologies are used to identify them? Which resources are available to analyse them? What mutations are expected in these non-coding sequences? Learning from developmental process, what is their expected role in CHD?

Blepharophimosis, Ptosis, Epicanthus Inversus Syndrome: New Report with a 197-kb Deletion Upstream of FOXL2 and Review of the Literature

Blepharophimosis, ptosis, and epicanthus inversus syndrome (BPES) is due to heterozygous FOXL2 intragenic mutations in about 70% of the patients, whereas total or partial gene deletions account for a minority of cases. Alteration of FOXL2 regulatory elements has been rarely described in patients with BPES. In this study, a prepubertal girl with BPES due to a 197-kb de novo deletion of the regulatory elements upstream of FOXL2 is reported. This girl presented with additional clinical features such as a soft cleft palate and microcephaly; thus, this copy number variant might have other somatic effects. The present deletion encompasses 2 coding genes (MRPS22 and COPB2), whose homozygous mutations have been associated with microcephaly. In our case, the sequences of the non-deleted allele were normal, ruling out a compound genetic defect. Normal levels of new biomarkers of ovarian reserve (anti-müllerian hormone, inhibin B) likely indicate an early diagnosis of type 2 BPES, but an evolutive gonadal damage will be excluded only by long-term follow-up. Additional reports of microdeletions upstream of FOXL2 are needed to better define the underlying genetic mechanism and the related phenotypic spectrum; the ability of the new hormonal markers to predict ovarian function in adolescence and adulthood should be confirmed.

The regulation of Sox9 expression in the gonad

The bipotential nature of cell types in the early developing gonad and the process of sex determination leading to either testis or ovary differentiation makes this an interesting system in which to study transcriptional regulation of gene expression and cell fate decisions. SOX9 is a transcription factor with multiple roles during development, including being a key player in mediating testis differentiation and therefore subsequent male development. Loss of Sox9 expression in both humans and mice results in XY female development, whereas its inappropriate activation in XX embryonic gonads can give male development. Multiple cases of Disorders of Sex Development in human patients or sex reversal in mice and other vertebrates can be explained by mutations affecting upstream regulators of Sox9 expression, such as the product of the Y chromosome gene Sry that triggers testis differentiation. Other cases are due to mutations in the Sox9 gene itself, including its own regulatory region. Indeed, rearrangements in and around the Sox9 genomic locus indicate the presence of multiple critical enhancers and the complex nature of its regulation. Here we summarize what is known about the role of Sox9 and its regulation during gonad development, including recently discovered critical enhancers. We also discuss higher order chromatin organization and how this might be involved. We end with some interesting future directions that have the potential to further enrich our understanding on the complex, multi-layered regulation controlling Sox9 expression in the gonads.

Gonadal supporting cells acquire sex-specific chromatin landscapes during mammalian sex determination

Cis-regulatory elements are critical for the precise spatiotemporal regulation of genes during development. However, identifying functional regulatory sites that drive cell differentiation in vivo has been complicated by the high numbers of cells required for whole-genome epigenetic assays. Here, we identified putative regulatory elements during sex determination by performing ATAC-seq and ChIP-seq for H3K27ac in purified XX and XY gonadal supporting cells before and after sex determination in mice. We show that XX and XY supporting cells initiate sex determination with similar chromatin landscapes and acquire sex-specific regulatory elements as they commit to the male or female fate. To validate our approach, we identified a functional gonad-specific enhancer downstream of Bmp2, an ovary-promoting gene. This work increases our understanding of the complex regulatory network underlying mammalian sex determination and provides a powerful resource for identifying non-coding regulatory elements that could harbor mutations that lead to Disorders of Sexual Development.

Defining the transcription landscape of the Gram-negative marine bacterium Vibrio harveyi

Vibrio harveyi is a Gram-negative pathogenic bacterium ubiquitously present in natural aquatic systems. Although environmental adaptability in V. harveyi may be enabled by profound reprogramming of gene expression previously observed during responses to starvation, suboptimal temperatures and other stress factors, the key characteristics of V. harveyi transcripts and operons, such as their boundaries and size as well as location of small RNA genes, remain largely unknown. To reveal the main features of the V. harveyi transcriptome, total RNA of this organism was analyzed by differential RNA sequencing (dRNA-seq). Analysis of the dRNA-seq data made it possible to define the primary transcriptome of V. harveyi along with cis-acting regulatory elements (riboswitches and leader sequences) and new genes. The latter encode a number of putative polypeptides and new phylogenetically conserved antisense RNAs potentially involved in the post-transcriptional control of gene expression.

Reconsideration of the immunotherapeutic pediatric safe dose levels of aluminum

FDA regulations require safety testing of constituent ingredients in drugs (21 CFR 610.15). With the exception of extraneous proteins, no component safety testing is required for vaccines or vaccine schedules. The dosing of aluminum in vaccines is based on the production of antibody titers, not safety science. Here we estimate a Pediatric Dose Limit that considers body weight. We identify several serious historical missteps in past analyses of provisional safe levels of aluminum in vaccines, and provide updates relevant to infant aluminum exposure in the pediatric schedule considering pediatric body weight. When aluminum doses are estimated from Federal Regulatory Code given body weight, exposure from the current vaccine schedule are found to exceed our estimate of a weight-corrected Pediatric Dose Limit. Our calculations show that the levels of aluminum suggested by the currently used limits place infants at risk of acute, repeated, and possibly chronic exposures of toxic levels of aluminum in modern vaccine schedules. Individual adult exposures are on par with Provisional Tolerable Weekly Intake "limits", but some individuals may be aluminum intolerant due to genetics or previous exposures. Vaccination in neonates and low birth-weight infants must be re-assessed; other implications for the use of aluminum-containing vaccines, and additional limitations in our understanding of neurotoxicity and safety levels of aluminum in biologics are discussed.

Characterization of regulatory elements in OsRGLP2 gene promoter from different rice accessions through sequencing and in silico evaluation

Germins and germin-like proteins from cupin superfamily contribute resistance to heat denaturation, chemical degradation and against plant pathogens, further functions in plant growth and development. In this study, from three different Oryza sativa accessions KS-282 and Pak 7178 and Pak 7865, OsRGLP2 gene promoter region was amplified, sequenced and analyzed. Sequencing data was evaluated via different computational tools. The regulatory elements were predicted by Consite tool and mapping was done. Many transcription factors binding sites were discovered in OsRGLP2 gene promoter; among these factors, HFH-1 having a significant role in germination was picked for further investigation. To study the interaction between HFH-1 and corresponding regulatory factors, HADDOCK Webserver was used. Graphical models for the interactions of HFH-1 and related regulatory elements were studied by graphic molecular system PyMOL. Mapping of cis-acting regulatory elements in OsRGLP2 gene promoter from three rice accessions showed differences in their position and copy number. Important regulatory elements found in OsRGLP2 promoter region were TATA, CAAT Box, ARR1, GATA, AGAAA, CAAT and DNA-binding One Zinc Finger (Dof) factors, few of them contribute to the regulation of plant defensive system, light responses, developmental and growth activities. Furthermore, during DNA interaction studies, it was found that HFH-1 transcription factor participates in hydrogen bonds formation with thymine and adenine bases.