Regulatory Region Sequences Research Articles

Repeated Shifts in Sociality Are Associated With Fine-tuning of Highly Conserved and Lineage-Specific Enhancers in a Socially Flexible Bee.

Comparative genomic studies of social insects suggest that changes in gene regulation are associated with evolutionary transitions in social behavior, but the activity of predicted regulatory regions has not been tested empirically. We used self-transcribing active regulatory region sequencing, a high-throughput enhancer discovery tool, to identify and measure the activity of enhancers in the socially variable sweat bee, Lasioglossum albipes. We identified over 36,000 enhancers in the L. albipes genome from 3 social and 3 solitary populations. Many enhancers were identified in only a subset of L. albipes populations, revealing rapid divergence in regulatory regions within this species. Population-specific enhancers were often proximal to the same genes across populations, suggesting compensatory gains and losses of regulatory regions may preserve gene activity. We also identified 1,182 enhancers with significant differences in activity between social and solitary populations, some of which are conserved regulatory regions across species of bees. These results indicate that social trait variation in L. albipes is associated with the fine-tuning of ancient enhancers as well as lineage-specific regulatory changes. Combining enhancer activity with population genetic data revealed variants associated with differences in enhancer activity and identified a subset of differential enhancers with signatures of selection associated with social behavior. Together, these results provide the first empirical map of enhancers in a socially flexible bee and highlight links between cis-regulatory variation and the evolution of social behavior.

Unraveling the Genetics of Shared Clinical and Serological Manifestations in Patients With Systemic Inflammatory Autoimmune Diseases.

Systemic inflammatory autoimmune diseases (SIADs) such as systemic lupus erythematosus (SLE), primary Sjögren disease (pSS), and idiopathic inflammatory myopathies (myositis) are complex conditions characterized by shared circulating autoantibodies and clinical manifestations, including skin rashes, among others. This study was aimed at elucidating the genetics underlying these common features. We performed targeted DNA sequencing of coding and regulatory regions from approximately 1,900 immune-related genes in a large cohort of 2,292 well-characterized Scandinavian patients with SIADs with SLE, pSS, and myositis as well as 1,252 controls. A gene-based functionally weighted genetic score for aggregate testing of all genetic variants, including rare variants, was complemented by in silico functional analyses and in vitro reporter experiments. Case-control association analysis detected known and potentially novel genetic loci in agreement with previous genetic and transcriptomics findings linked to the SIAD autoimmune background. Intriguingly, case-case comparisons between patient subgroups with and without specific autoantibodies revealed that the subgroups defined by antinuclear antibodies and anti-double-stranded DNA antibodies have unique genetic profiles reflecting their heterogeneity. When focusing on clinical features, we overall showed that dual-specificity phosphatase 1 (DUSP1) protective genetic variants lead to increased gene expression and potentially to anti-inflammatory effects on the SIAD-associated skin phenotype. This is consistent with recent genetic findings on eczema and with the previously reported down-regulation of the MAPK signaling-related gene DUSP1 in other skin disorders. Together, this suggests common molecular mechanisms potentially underlying overlapping clinical manifestations shared among different disorders and informs clinical heterogeneity, which could be translated to improve disease diagnostic and treatment, also in more generalized disease frameworks.

Principle and application of self-transcribing active regulatory region sequencing in enhancer discovery research.

Self-transcribing active regulatory region sequencing (STARR-seq) is a high-throughput sequencing method capable of simultaneously discovering and validating all enhancers within the genome. In this method, candidate sequences are inserted into plasmid vectors and electroporated into cells. Acting as both enhancers and target genes, the self-transcription of these sequences will also be enhanced by themselves. By sequencing the transcriptome and comparing the results with the non-inserted control, the locations and activity of enhancers can be determined. In traditional enhancer discovery strategies, the chromatin open regions and transcription active regions were sequenced and predicted as enhancers. However, the activity of these putative enhancers could only be validated one by one without a high-throughput method. STARR-seq solved this limitation, allowing simultaneous enhancers discovery and activity validation in a high-throughput manner. Since the introduction of STARR-seq, it has been widely used to discover enhancers and validate enhancer activity in a number of organisms and cells. In this review, we present the traditional enhancer prediction methods and the basic principles, development history, specific applications of STARR-seq, and its future prospects, aiming to provide a reference for researchers in related fields conducting enhancer studies.

Enhancers associated with unstable RNAs are rare in plants

Unstable transcripts have emerged as markers of active enhancers in vertebrates and shown to be involved in many cellular processes and medical disorders. However, their prevalence and role in plants is largely unexplored. Here, we comprehensively captured all actively initiating (nascent) transcripts across diverse crops and other plants using capped small (cs)RNA sequencing. We discovered that unstable transcripts are rare in plants, unlike in vertebrates, and when present, often originate from promoters. In addition, many ‘distal’ elements in plants initiate tissue-specific stable transcripts and are likely bona fide promoters of as-yet-unannotated genes or non-coding RNAs, cautioning against using reference genome annotations to infer putative enhancer sites. To investigate enhancer function, we integrated data from self-transcribing active regulatory region (STARR) sequencing. We found that annotated promoters and other regions that initiate stable transcripts, but not those marked by unstable or bidirectional unstable transcripts, showed stronger enhancer activity in this assay. Our findings underscore the blurred line between promoters and enhancers and suggest that cis-regulatory elements can encompass diverse structures and mechanisms in eukaryotes, including humans.

Constructing eRNA-mediated gene regulatory networks to explore the genetic basis of muscle and fat-relevant traits in pigs

BackgroundEnhancer RNAs (eRNAs) play a crucial role in transcriptional regulation. While significant progress has been made in understanding epigenetic regulation mediated by eRNAs, research on the construction of eRNA-mediated gene regulatory networks (eGRN) and the identification of critical network components that influence complex traits is lacking.ResultsHere, employing the pig as a model, we conducted a comprehensive study using H3K27ac histone ChIP-seq and RNA-seq data to construct eRNA expression profiles from multiple tissues of two distinct pig breeds, namely Enshi Black (ES) and Duroc. In addition to revealing the regulatory landscape of eRNAs at the tissue level, we developed an innovative network construction and refinement method by integrating RNA-seq, ChIP-seq, genome-wide association study (GWAS) signals and enhancer-modulating effects of single nucleotide polymorphisms (SNPs) measured by self-transcribing active regulatory region sequencing (STARR-seq) experiments. Using this approach, we unraveled eGRN that significantly influence the growth and development of muscle and fat tissues, and identified several novel genes that affect adipocyte differentiation in a cell line model.ConclusionsOur work not only provides novel insights into the genetic basis of economic pig traits, but also offers a generalizable approach to elucidate the eRNA-mediated transcriptional regulation underlying a wide spectrum of complex traits for diverse organisms.

Human gene regulatory evolution is driven by the divergence of regulatory element function in both cis and trans

Gene regulatory divergence between species can result from cis-acting local changes to regulatory element DNA sequences or global trans-acting changes to the regulatory environment. Understanding how these mechanisms drive regulatory evolution has been limited by challenges in identifying trans-acting changes. We present a comprehensive approach to directly identify cis- and trans-divergent regulatory elements between human and rhesus macaque lymphoblastoid cells using assay for transposase-accessible chromatin coupled to self-transcribing active regulatory region (ATAC-STARR) sequencing. In addition to thousands of cis changes, we discover an unexpected number (∼10,000) of trans changes and show that cis and trans elements exhibit distinct patterns of sequence divergence and function. We further identify differentially expressed transcription factors that underlie ∼37% of trans differences and trace how cis changes can produce cascades of trans changes. Overall, we find that most divergent elements (67%) experienced changes in both cis and trans, revealing a substantial role for trans divergence-alone and together with cis changes-in regulatory differences between species.

Identification of Highly Repetitive Enhancers with Long-range Regulation Potential in Barley via STARR-seq.

Enhancers are DNA sequences that can strengthen transcription initiation. However, the global identification of plant enhancers is complicated due to uncertainty in the distance and orientation of enhancers, especially in species with large genomes. In this study, we performed self-transcribing active regulatory region sequencing (STARR-seq) for the first time to identify enhancers across the barley genome. A total of 7323 enhancers were successfully identified, and among 45 randomly selected enhancers, over 75% were effective as validated by a dual-luciferase reporter assay system in the lower epidermis of tobacco leaves. Interestingly, up to 53.5% of the barley enhancers were repetitive sequences, especially transposable elements (TEs), thus reinforcing the vital role of repetitive enhancers in gene expression. Both the common active mark H3K4me3 and repressive mark H3K27me3 were abundant among the barley STARR-seq enhancers. In addition, the functional range of barley STARR-seq enhancers seemed much broader than that of rice or maize and extended to ±100 kb of the gene body, and this finding was consistent with the high expression levels of genes in the genome. This study specifically depicts the unique features of barley enhancers and provides available barley enhancers for further utilization.

Molecular Cloning and Functional Characterization of the 5′ Regulatory Region of the SLC11A1 Gene from Yaks

The solute transport protein family 11 A1 (SLC11A1), also recognized as natural resistance-associated macrophage protein 1 (NRAMP1), represents a transmembrane protein encoded by the SLC11A1 gene. A variety of prior investigations have illuminated its involvement in conferring resistance or susceptibility to bacterial agents, positioning it as a promising candidate gene for breeding disease-resistant animals. Yaks (Bos grunniens), renowned inhabitants of the Qinghai-Tibet Plateau in China, stand as robust ruminants distinguished by their adaptability and formidable disease resistance. Notwithstanding these unique traits, there is scant literature on the SLC11A1 gene in the yak population. Our inquiry commences with the cloning of the 5′ regulatory region sequence of the Zhongdian yak SLC11A1 gene. We employ bioinformatics tools to identify transcription factor binding sites, delineating pivotal elements like enhancers and cis-acting elements. To ascertain the promoter activity of this region, we amplify four distinct promoter fragments within the 5′ regulatory region of the yak SLC11A1 gene. Subsequently, we design a luciferase reporter gene vector containing four site-specific deletion mutations and perform transient transfection experiments. Through these experiments, we measure and compare the activity of disparate gene fragments located within the 5′ regulatory region, revealing regions bearing promoter functionality and discerning key regulatory elements. Our findings validate the promoter functionality of the 5′ regulatory region, offering preliminary insights into the core and principal regulatory segments of this promoter. Notably, we identified single nucleotide polymorphisms (SNPs) that may be associated with important regulatory elements such as NF-1 and NF-1/L. This study provides a theoretical framework for in-depth research on the function and expression regulation mechanism of the yak SLC11A1 gene.

Identification of human genetic variants modulating the course of COVID-19 infection with importance in other viral infections.

Introduction: COVID-19 has been a major focus of scientific research since early 2020. Due to its societal, economic, and clinical impact worldwide, research efforts aimed, among other questions, to address the effect of host genetics in susceptibility and severity of COVID-19. Methods: We, therefore, performed next-generation sequencing of coding and regulatory regions of 16 human genes, involved in maintenance of the immune system or encoding receptors for viral entry into the host cells, in a subset of 60 COVID-19 patients from the General Hospital Tešanj, Bosnia and Herzegovina, classified into three groups of clinical conditions of different severity ("mild," "moderate," and "severe"). Results: We confirmed that the male sex and older age are risk factors for severe clinical picture and identified 13 variants on seven genes (CD55, IL1B, IL4, IRF7, DDX58, TMPRSS2, and ACE2) with potential functional significance, either as genetic markers of modulated susceptibility to SARS-CoV-2 infection or modifiers of the infection severity. Our results include variants reported for the first time as potentially associated with COVID-19, but further research and larger patient cohorts are required to confirm their effect. Discussion: Such studies, focused on candidate genes and/or variants, have a potential to answer the questions regarding the effect of human genetic makeup on the expected infection outcome. In addition, loci we identified here were previously reported to have clinical significance in other diseases and viral infections, thus confirming a general, broader significance of COVID-19-related research results following the end of the pandemic period.

Variability of Inverted Repeats in All Available Genomes of Bacteria.

Noncanonical secondary structures in nucleic acids have been studied intensively in recent years. Important biological roles of cruciform structures formed by inverted repeats (IRs) have been demonstrated in diverse organisms, including humans. Using Palindrome analyser, we analyzed IRs in all accessible bacterial genome sequences to determine their frequencies, lengths, and localizations. IR sequences were identified in all species, but their frequencies differed significantly across various evolutionary groups. We detected 242,373,717 IRs in all 1,565 bacterial genomes. The highest mean IR frequency was detected in the Tenericutes (61.89 IRs/kbp) and the lowest mean frequency was found in the Alphaproteobacteria (27.08 IRs/kbp). IRs were abundant near genes and around regulatory, tRNA, transfer-messenger RNA (tmRNA), and rRNA regions, pointing to the importance of IRs in such basic cellular processes as genome maintenance, DNA replication, and transcription. Moreover, we found that organisms with high IR frequencies were more likely to be endosymbiotic, antibiotic producing, or pathogenic. On the other hand, those with low IR frequencies were far more likely to be thermophilic. This first comprehensive analysis of IRs in all available bacterial genomes demonstrates their genomic ubiquity, nonrandom distribution, and enrichment in genomic regulatory regions. IMPORTANCE Our manuscript reports for the first time a complete analysis of inverted repeats in all fully sequenced bacterial genomes. Thanks to the availability of unique computational resources, we were able to statistically evaluate the presence and localization of these important regulatory sequences in bacterial genomes. This work revealed a strong abundance of these sequences in regulatory regions and provides researchers with a valuable tool for their manipulation.

Cis-regulatory atlas of primary human CD4+ T cells

Cis-regulatory elements (CRE) are critical for coordinating gene expression programs that dictate cell-specific differentiation and homeostasis. Recently developed self-transcribing active regulatory region sequencing (STARR-Seq) has allowed for genome-wide annotation of functional CREs. Despite this, STARR-Seq assays are only employed in cell lines, in part, due to difficulties in delivering reporter constructs. Herein, we implemented and validated a STARR-Seq–based screen in human CD4+ T cells using a non-integrating lentiviral transduction system. Lenti-STARR-Seq is the first example of a genome-wide assay of CRE function in human primary cells, identifying thousands of functional enhancers and negative regulatory elements (NREs) in human CD4+ T cells. We find an unexpected difference in nucleosome organization between enhancers and NRE: enhancers are located between nucleosomes, whereas NRE are occupied by nucleosomes in their endogenous locations. We also describe chromatin modification, eRNA production, and transcription factor binding at both enhancers and NREs. Our findings support the idea of silencer repurposing as enhancers in alternate cell types. Collectively, these data suggest that Lenti-STARR-Seq is a successful approach for CRE screening in primary human cell types, and provides an atlas of functional CREs in human CD4+ T cells.

An unbiased AAV-STARR-seq screen revealing the enhancer activity map of genomic regions in the mouse brain in vivo

Enhancers are important cis-regulatory elements controlling cell-type specific expression patterns of genes. Furthermore, combinations of enhancers and minimal promoters are utilized to construct small, artificial promoters for gene delivery vectors. Large-scale functional screening methodology to construct genomic maps of enhancer activities has been successfully established in cultured cell lines, however, not yet applied to terminally differentiated cells and tissues in a living animal. Here, we transposed the Self-Transcribing Active Regulatory Region Sequencing (STARR-seq) technique to the mouse brain using adeno-associated-viruses (AAV) for the delivery of a highly complex screening library tiling entire genomic regions and covering in total 3 Mb of the mouse genome. We identified 483 sequences with enhancer activity, including sequences that were not predicted by DNA accessibility or histone marks. Characterizing the expression patterns of fluorescent reporters controlled by nine candidate sequences, we observed differential expression patterns also in sparse cell types. Together, our study provides an entry point for the unbiased study of enhancer activities in organisms during health and disease.

A second-generation capture panel for cost-effective sequencing of genome regulatory regions in wheat and relatives.

As genome resources for wheat (Triticum L.) expand at a rapid pace, it is important to update targeted sequencing tools to incorporate improved sequence assemblies and regions of previously unknown significance. Here, we developed an updated regulatory region enrichment capture for wheat and other Triticeae species. The core target space includes sequences from 2-Kbp upstream of each gene predicted in the Chinese Spring wheat genome (IWGSC RefSeq Annotation v1.0) and regions of open chromatin identified with an assay for transposase-accessible chromatin using sequencing from wheat leaf and root samples. To improve specificity, we aggressively filtered candidate repetitive sequences using a combination of nucleotide basic local alignment search tool (BLASTN) searches to the Triticeae Repetitive Sequence Database (TREP), identification of regions with read over-coverage from previous target enrichment experiments, and k-mer frequency analyses. The final design comprises 216.5 Mbp of predicted hybridization space in hexaploid wheat and showed increased specificity and coverage of targeted sequences relative to previous protocols. Test captures on hexaploid and tetraploid wheat and other diploid cereals show that the assay has broad potential utility for cost-effective promoter and open chromatin resequencing and general-purpose genotyping of various Triticeae species.

STARR-seq for high-throughput identification of plant enhancers.

Self-transcribing active regulatory region sequencing (STARR-seq) is a convenient method for genome-wide identification of plant enhancers. The basic principle is that an enhancer is able to activate the reporter gene expression (e.g., GFP) in a well-designed transient expression vector. The constructed vectors, which harbor the enhancer candidates, are transiently transformed into tobacco leaves or various protoplasts. The resulting reporter transcripts are enriched for sequencing. In the original design, the locations and activities of the enhancers are predicted by mapping to the reference genome and enrichment analysis. Recent research in plants showed the robustness of adding barcode sequences to link transcripts to the candidates placed upstream. The enhancer activities can, thus, be evaluated as the relative reporter expression level driven by the candidate versus the negative control.

De-erosion of X chromosome dosage compensation by the editing of XIST regulatory regions restores the differentiation potential in hPSCs

Human pluripotent stem cells (hPSCs) regularly and irreversibly show the erosion of X chromosome inactivation (XCI) by long non-coding RNA (lncRNA) XIST silencing, causing challenges in various applications of female hPSCs. Here, we report reliable methods to reactivate XIST with monoallelic expression in female hPSCs. Surprisingly, we find that the editing of XIST regulatory regions by Cas9-mediated non-homologous end joining is sufficient for the reactivation of XIST by endogenous systems. Proliferated hPSCs with XIST reactivation show XCI from an eroded X chromosome, suggesting that hPSCs with normal dosage compensation might lead to a growth advantage. Furthermore, the use of targeting vectors, including the XIST regulatory region sequences and selection cassette, enables XIST reactivation in hPSCs with high efficiency. XIST-reactivated hPSCs can show the restoration of differentiation potential. Thus, our findings demonstrate that XIST re-expression is a beneficial method to maximize the use of female hPSCs in various applications, such as proper disease modeling.

A second generation capture panel for cost-effective sequencing of genome regulatory regions in wheat and relatives

A second generation capture panel for cost-effective sequencing of genome regulatory regions in wheat and relatives

A second generation capture panel for cost-effective sequencing of genome regulatory regions in wheat and relatives

A second generation capture panel for cost-effective sequencing of genome regulatory regions in wheat and relatives

A second generation capture panel for cost-effective sequencing of genome regulatory regions in wheat and relatives

A second generation capture panel for cost-effective sequencing of genome regulatory regions in wheat and relatives

A second generation capture panel for cost-effective sequencing of genome regulatory regions in wheat and relatives

A second generation capture panel for cost-effective sequencing of genome regulatory regions in wheat and relatives

A second generation capture panel for cost-effective sequencing of genome regulatory regions in wheat and relatives

A second generation capture panel for cost-effective sequencing of genome regulatory regions in wheat and relatives