Cis and trans effects differentially contribute to the evolution of promoters and enhancers

Contribution of the transcription factor SfGATAe to Bt Cry toxin resistance in Spodoptera frugiperda through reduction of ABCC2 expression

Insect resistance evolution poses a significant threat to the advantages of biopesticides and transgenic crops utilizing insecticidal Cry-toxins from Bacillus thuringiensis (Bt). However, there is limited research on the relationship between transcriptional regulation of specific toxin receptors in lepidopteran insects and their resistance to Bt toxins. Here, we report the positive regulatory role of the SfGATAe transcription factor on the expression of the ABCC2 gene in Spodoptera frugiperda. DNA regions in the SfABCC2 promoter that are vital for regulation by SfGATAe, utilizing DAP-seq technology and promoter deletion mapping. Through yeast one-hybrid assays, DNA pull-down experiments, and site-directed mutagenesis, we confirmed that the transcription factor SfGATAe regulates the core control site PBS2 in the ABCC2 target gene. Tissue-specific expression analysis has revealed that SfGATAe is involved in the regulation and expression of midgut cells in the fall armyworm. Silencing SfGATAe in fall armyworm larvae resulted in reduced expression of SfABCC2 and decreased sensitivity to Cry1Ac toxin. Overall, this study elucidated the regulatory mechanism of the transcription factor SfGATAe on the expression of the toxin receptor gene SfABCC2 and this transcriptional control mechanism impacts the resistance of the fall armyworm to Bt toxins.

Massively parallel screen uncovers many rare 3' UTR variants regulating mRNA abundance of cancer driver genes

Understanding the function of rare non-coding variants represents a significant challenge. Using MapUTR, a screening method, we studied the function of rare 3' UTR variants affecting mRNA abundance post-transcriptionally. Among 17,301 rare gnomAD variants, an average of 24.5% were functional, with 70% in cancer-related genes, many in critical cancer pathways. This observation motivated an interrogation of 11,929 somatic mutations, uncovering 3928 (33%) functional mutations in 155 cancer driver genes. Functional MapUTR variants were enriched in microRNA- or protein-binding sites and may underlie outlier gene expression in tumors. Further, we introduce untranslated tumor mutational burden (uTMB), a metric reflecting the amount of somatic functional MapUTR variants of a tumor and show its potential in predicting patient survival. Through prime editing, we characterized three variants in cancer-relevant genes (MFN2, FOSL2, and IRAK1), demonstrating their cancer-driving potential. Our study elucidates the function of tens of thousands of non-coding variants, nominates non-coding cancer driver mutations, and demonstrates their potential contributions to cancer.

Transforming our understanding of species-specific gene regulation

Mechanisms underlying phenotypic divergence across species remain unresolved. In this issue of Cell Genomics, Hansen, Fong, et al.1 systematically dissect human and rhesus macaque gene expression divergence by screening tens of thousands of orthologous elements for enhancer activity in lymphoblastoid cell lines, revealing a much greater role for trans divergence at levels equal to those of cis effects, counter to the prevailing consensus in the field.

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.

Super-enhancers: drivers of cells' identities and cells' debacles

Precise spatiotemporal regulations of gene expression are essential for determining cells' fates and functions. Enhancers are cis-acting DNA elements that act as periodic transcriptional thrusters and their activities are cell type specific. Clusters of enhancers, called super-enhancers, are more densely occupied by transcriptional activators than enhancers, driving stronger expression of their target genes, which have prominent roles in establishing and maintaining cellular identities. Here we review the current knowledge on the composition and structure of super-enhancers to understand how they robustly stimulate the expression of cellular identity genes. We also review their involvement in the development of various cell types and both noncancerous and cancerous disorders, implying the therapeutic interest of targeting them to fight against various diseases.

Dissecting the sequence and structural determinants guiding m6A deposition and evolution via inter- and intra-species hybrids

BACKGROUND

N6-methyladenosine (m6A) is the most abundant mRNA modification, and controls mRNA stability. m6A distribution varies considerably between and within species. Yet, it is unclear to what extent this variability is driven by changes in genetic sequences ('cis') or cellular environments ('trans') and via which mechanisms.

RESULTS

Here we dissect the determinants governing RNA methylation via interspecies and intraspecies hybrids in yeast and mammalian systems, coupled with massively parallel reporter assays and m6A-QTL reanalysis. We find that m6A evolution and variability is driven primarily in 'cis', via two mechanisms: (1) variations altering m6A consensus motifs, and (2) variation impacting mRNA secondary structure. We establish that mutations impacting RNA structure - even when distant from an m6A consensus motif - causally dictate methylation propensity. Finally, we demonstrate that allele-specific differences in m6A levels lead to allele-specific changes in gene expression.

CONCLUSIONS

Our findings define the determinants governing m6A evolution and diversity and characterize the consequences thereof on gene expression regulation.

Causal effects of lipid-lowering therapies on aging-related outcomes and risk of cancers: a drug-target Mendelian randomization study

BACKGROUND

Despite the widespread use of statins, newer lipid-lowering drugs have been emerging. It remains unclear how the long-term use of novel lipid-lowering drugs affects the occurrence of cancers and age-related diseases.

METHODS

A drug-target Mendelian randomization study was performed. Genetic variants of nine lipid-lowering drug-target genes (HMGCR, PCKS9, NPC1L1, LDLR, APOB, CETP, LPL, APOC3, and ANGPTL3) were extracted as exposures from the summary data of Global Lipids Genetics Consortium Genome-Wide Association Studies (GWAS). GWAS summary data of cancers and noncancerous diseases were used as outcomes. The inverse-variance weighted method was applied as the main statistical approach. Sensitivity tests were conducted to evaluate the robustness, pleiotropy, and heterogeneity of the results.

RESULTS

In addition to marked effects on decreased risks of atherosclerotic cardiovascular diseases, genetically proxied lipid-lowering variants of PCKS9, CETP, LPL, LDLR, and APOC3 were associated with longer human lifespans (q<0.05). Lipid-lowering variants of ANGPTL3 and LDLR were associated with reduced risks of colorectal cancer, and ANGPTL3 was also associated with lower risks of gastric cancer (q<0.05). Lipid-lowering LPL variants were associated with decreased risks of hypertension, type 2 diabetes, nonalcoholic fatty liver disease, and bladder cancer (q<0.05). Lipid-lowering variants of PCKS9 and HMGCR were associated with decreased risks of osteoporosis (q<0.05). Lipid-lowering APOB variants were associated with a decreased risk of thyroid cancer (q<0.05).

CONCLUSIONS

Our study provides genetic evidence that newer nonstatin lipid-lowering agents have causal effects on decreased risks of several common cancers and cardiometabolic diseases. These data provide genetic insights into the potential benefits of newer nonstatin therapies.

High-throughput functional analysis of regulatory variants using a massively parallel reporter assay

Association studies describe genetic associations between noncoding variants and disease susceptibility; however, they do not provide functional insight into the underlying molecular mechanisms of these variants. We present a protocol to assay the regulatory potential of thousands of noncoding variants using massively parallel reporter assays. We describe steps for oligo design, generating a plasmid pool, and extracting tag-seq libraries from cells to quantify the tested sequences. For complete details on the use and execution of this protocol, please refer to Oliveros and Delfosse et al.1.

Mild dehydration effects on the murine kidney single-nucleus transcriptome and chromatin accessibility

Daily, we may experience mild dehydration with a rise in plasma osmolality that triggers the release of vasopressin. Although the effect of dehydration is well characterized in collecting duct principal cells (CDPCs), we hypothesized that mild dehydration (<12 h) results in many kidney cell-specific changes in transcriptomes and chromatin accessibility. Single-nucleus (sn) multiome (RNA-assay for transposase-accessible chromatin) sequencing and bulk RNA sequencing of kidneys from male and female mice that were mildly water deprived or not were compared. Water-deprived mice had a significant increase in plasma osmolality. sn-multiome-seq resulted in 19,837 nuclei that were annotated into 33 clusters. In CDPCs, aquaporin 2 (Aqp2) and aquaporin 3 (Apq3) were greater in dehydrated mice, but there were novel genes like gremlin 2 (Grem2; a cytokine) that were increased compared with ad libitum mice. The transcription factor cAMP-responsive element modulator (Crem) was greater in CDPCs of dehydrated mice, and the Crem DNA motif was more accessible. There were hundreds of sex- and dehydration-specific differentially expressed genes (DEGs) throughout the kidney, especially in the proximal tubules and thin limbs. In male mice, DEGs were enriched in pathways related to lipid metabolism, whereas female DEGs were enriched in organic acid metabolism. Many highly expressed genes had a positive correlation with increased chromatin accessibility, and mild dehydration exerted many transcriptional changes that we detected at the chromatin level. Even with a rise in plasma osmolality, male and female kidneys have distinct transcriptomes suggesting that there may be diverse mechanisms used to remain in fluid balance.NEW & NOTEWORTHY The kidney consists of >30 cell types that work collectively to maintain fluid-electrolyte balance. Kidney single-nucleus transcriptomes and chromatin accessibility profiles from male and female control (ad libitum water and food) or mildly dehydrated mice (ad libitum food, water deprivation) were determined. Mild dehydration caused hundreds of cell- and sex-specific transcriptomic changes, even though the kidney function to conserve water was the same.

Hybrid Epigenomes Reveal Extensive Local Genetic Changes to Chromatin Accessibility Contribute to Divergence in Embryonic Gene Expression Between Species

Chromatin accessibility plays an important role in shaping gene expression, yet little is known about the genetic and molecular mechanisms that influence the evolution of chromatin configuration. Both local (cis) and distant (trans) genetic influences can in principle influence chromatin accessibility and are based on distinct molecular mechanisms. We, therefore, sought to characterize the role that each of these plays in altering chromatin accessibility in 2 closely related sea urchin species. Using hybrids of Heliocidaris erythrogramma and Heliocidaris tuberculata, and adapting a statistical framework previously developed for the analysis of cis and trans influences on the transcriptome, we examined how these mechanisms shape the regulatory landscape at 3 important developmental stages, and compared our results to similar analyses of the transcriptome. We found extensive cis- and trans-based influences on evolutionary changes in chromatin, with cis effects generally larger in effect. Evolutionary changes in accessibility and gene expression are correlated, especially when expression has a local genetic basis. Maternal influences appear to have more of an effect on chromatin accessibility than on gene expression, persisting well past the maternal-to-zygotic transition. Chromatin accessibility near gene regulatory network genes appears to be distinctly regulated, with trans factors appearing to play an outsized role in the configuration of chromatin near these genes. Together, our results represent the first attempt to quantify cis and trans influences on evolutionary divergence in chromatin configuration in an outbred natural study system and suggest that chromatin regulation is more genetically complex than was previously appreciated.

Deep flanking sequence engineering for efficient promoter design using DeepSEED

Designing promoters with desirable properties is essential in synthetic biology. Human experts are skilled at identifying strong explicit patterns in small samples, while deep learning models excel at detecting implicit weak patterns in large datasets. Biologists have described the sequence patterns of promoters via transcription factor binding sites (TFBSs). However, the flanking sequences of cis-regulatory elements, have long been overlooked and often arbitrarily decided in promoter design. To address this limitation, we introduce DeepSEED, an AI-aided framework that efficiently designs synthetic promoters by combining expert knowledge with deep learning techniques. DeepSEED has demonstrated success in improving the properties of Escherichia coli constitutive, IPTG-inducible, and mammalian cell doxycycline (Dox)-inducible promoters. Furthermore, our results show that DeepSEED captures the implicit features in flanking sequences, such as k-mer frequencies and DNA shape features, which are crucial for determining promoter properties.

Environmentally robust cis-regulatory changes underlie rapid climatic adaptation

Changes in gene expression are thought to play a major role in adaptive evolution. While it is known that gene expression is highly sensitive to the environment, very few studies have determined the influence of genetic and environmental effects on adaptive gene expression differences in natural populations. Here, we utilize allele-specific expression to characterize cis and trans gene regulatory divergence in temperate and tropical house mice in two metabolic tissues under two thermal conditions. First, we show that gene expression divergence is pervasive between populations and across thermal conditions, with roughly 5 to 10% of genes exhibiting genotype-by-environment interactions. Second, we found that most expression divergence was due to cis-regulatory changes that were stable across temperatures. In contrast, patterns of expression plasticity were largely attributable to trans-effects, which showed greater sensitivity to temperature. Nonetheless, we found a small subset of temperature-dependent cis-regulatory changes, thereby identifying loci underlying expression plasticity. Finally, we performed scans for selection in wild house mice to identify genomic signatures of rapid adaptation. Genomic outliers were enriched in genes with evidence for cis-regulatory divergence. Notably, these genes were associated with phenotypes that affected body weight and metabolism, suggesting that cis-regulatory changes are a possible mechanism for adaptive body size evolution between populations. Our results show that gene expression plasticity, largely controlled in trans, may facilitate the colonization of new environments, but that evolved changes in gene expression are largely controlled in cis, illustrating the genetic and nongenetic mechanisms underlying the establishment of populations in new environments.

Systematic characterization of regulatory variants of blood pressure genes

High blood pressure (BP) is the major risk factor for cardiovascular disease. Genome-wide association studies have identified genetic variants for BP, but functional insights into causality and related molecular mechanisms lag behind. We functionally characterize 4,608 genetic variants in linkage with 135 BP loci in vascular smooth muscle cells and cardiomyocytes by massively parallel reporter assays. High densities of regulatory variants at BP loci (i.e., ULK4, MAP4, CFDP1, PDE5A) indicate that multiple variants drive genetic association. Regulatory variants are enriched in repeats, alter cardiovascular-related transcription factor motifs, and spatially converge with genes controlling specific cardiovascular pathways. Using heuristic scoring, we define likely causal variants, and CRISPR prime editing finally determines causal variants for KCNK9, SFXN2, and PCGF6, which are candidates for developing high BP. Our systems-level approach provides a catalog of functionally relevant variants and their genomic architecture in two trait-relevant cell lines for a better understanding of BP gene regulation.

Applications of advanced technologies for detecting genomic structural variation

Chromosomal structural variation (SV) encompasses a heterogenous class of genetic variants that exerts strong influences on human health and disease. Despite their importance, many structural variants (SVs) have remained poorly characterized at even a basic level, a discrepancy predicated upon the technical limitations of prior genomic assays. However, recent advances in genomic technology can identify and localize SVs accurately, opening new questions regarding SV risk factors and their impacts in humans. Here, we first define and classify human SVs and their generative mechanisms, highlighting characteristics leveraged by various SV assays. We next examine the first-ever gapless assembly of the human genome and the technical process of assembling it, which required third-generation sequencing technologies to resolve structurally complex loci. The new portions of that "telomere-to-telomere" and subsequent pangenome assemblies highlight aspects of SV biology likely to develop in the near-term. We consider the strengths and limitations of the most promising new SV technologies and when they or longstanding approaches are best suited to meeting salient goals in the study of human SV in population-scale genomics research, clinical, and public health contexts. It is a watershed time in our understanding of human SV when new approaches are expected to fundamentally change genomic applications.

Chromosomal Instability in Genome Evolution: From Cancer to Macroevolution

The integrity of the genome is crucial for the survival of all living organisms. However, genomes need to adapt to survive certain pressures, and for this purpose use several mechanisms to diversify. Chromosomal instability (CIN) is one of the main mechanisms leading to the creation of genomic heterogeneity by altering the number of chromosomes and changing their structures. In this review, we will discuss the different chromosomal patterns and changes observed in speciation, in evolutional biology as well as during tumor progression. By nature, the human genome shows an induction of diversity during gametogenesis but as well during tumorigenesis that can conclude in drastic changes such as the whole genome doubling to more discrete changes as the complex chromosomal rearrangement chromothripsis. More importantly, changes observed during speciation are strikingly similar to the genomic evolution observed during tumor progression and resistance to therapy. The different origins of CIN will be treated as the importance of double-strand breaks (DSBs) or the consequences of micronuclei. We will also explain the mechanisms behind the controlled DSBs, and recombination of homologous chromosomes observed during meiosis, to explain how errors lead to similar patterns observed during tumorigenesis. Then, we will also list several diseases associated with CIN, resulting in fertility issues, miscarriage, rare genetic diseases, and cancer. Understanding better chromosomal instability as a whole is primordial for the understanding of mechanisms leading to tumor progression.

Machine learning dissection of human accelerated regions in primate neurodevelopment

Using machine learning (ML), we interrogated the function of all human-chimpanzee variants in 2,645 human accelerated regions (HARs), finding 43% of HARs have variants with large opposing effects on chromatin state and 14% on neurodevelopmental enhancer activity. This pattern, consistent with compensatory evolution, was confirmed using massively parallel reporter assays in chimpanzee and human neural progenitor cells. The species-specific enhancer activity of HARs was accurately predicted from the presence and absence of transcription factor footprints in each species. Despite these striking cis effects, activity of a given HAR sequence was nearly identical in human and chimpanzee cells. This suggests that HARs did not evolve to compensate for changes in the trans environment but instead altered their ability to bind factors present in both species. Thus, ML prioritized variants with functional effects on human neurodevelopment and revealed an unexpected reason why HARs may have evolved so rapidly.

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 ATAC-STARR-seq. 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 >50% 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-to regulatory differences between species.

Hybrid epigenomes reveal extensive local genetic changes to chromatin accessibility contribute to divergence in embryonic gene expression between species

Chromatin accessibility plays an important role in shaping gene expression patterns across development and evolution; however, little is known about the genetic and molecular mechanisms that influence chromatin configuration itself. Because cis and trans influences can both theoretically influence the accessibility of the epigenome, we sought to better characterize the role that both of these mechanisms play in altering chromatin accessibility in two closely related sea urchin species. Using hybrids of the two species, and adapting a statistical framework previously developed for the analysis of cis and trans influences on the transcriptome, we examined how these mechanisms shape the regulatory landscape at three important developmental stages, and compared our results to similar patterns in the transcriptome. We found extensive cis- and trans-based influences on evolutionary changes in chromatin, with cis effects slightly more numerous and larger in effect. Genetic mechanisms influencing gene expression and chromatin configuration are correlated, but differ in several important ways. Maternal influences also appear to have more of an effect on chromatin accessibility than on gene expression, persisting well past the maternal-to-zygotic transition. Furthermore, chromatin accessibility near GRN genes appears to be regulated differently than the rest of the epigenome, and indicates that trans factors may play an outsized role in the configuration of chromatin near these genes. Together, our results represent the first attempt to quantify cis and trans influences on evolutionary divergence in chromatin configuration in an outbred natural study system, and suggest that the regulation of chromatin is more genetically complex than was previously appreciated.

Current advances in primate genomics: novel approaches for understanding evolution and disease

Primate genomics holds the key to understanding fundamental aspects of human evolution and disease. However, genetic diversity and functional genomics data sets are currently available for only a few of the more than 500 extant primate species. Concerted efforts are under way to characterize primate genomes, genetic polymorphism and divergence, and functional landscapes across the primate phylogeny. The resulting data sets will enable the connection of genotypes to phenotypes and provide new insight into aspects of the genetics of primate traits, including human diseases. In this Review, we describe the existing genome assemblies as well as genetic variation and functional genomic data sets. We highlight some of the challenges with sample acquisition. Finally, we explore how technological advances in single-cell functional genomics and induced pluripotent stem cell-derived organoids will facilitate our understanding of the molecular foundations of primate biology.

A trans-omics assessment of gene-gene interaction in early-stage NSCLC

Epigenome-wide gene-gene (G × G) interactions associated with non-small-cell lung cancer (NSCLC) survival may provide insights into molecular mechanisms and therapeutic targets. Hence, we proposed a three-step analytic strategy to identify significant and robust G × G interactions that are relevant to NSCLC survival. In the first step, among 49 billion pairs of DNA methylation probes, we identified 175 775 G × G interactions with P_Bonferroni  ≤ 0.05 in the discovery phase of epigenomic analysis; among them, 15 534 were confirmed with P ≤ 0.05 in the validation phase. In the second step, we further performed a functional validation for these G × G interactions at the gene expression level by way of a two-phase (discovery and validation) transcriptomic analysis, and confirmed 25 significant G × G interactions enriched in the 6p21.33 and 6p22.1 regions. In the third step, we identified two G × G interactions using the trans-omics analysis, which had significant (P ≤ 0.05) epigenetic cis-regulation of transcription and robust G × G interactions at both the epigenetic and transcriptional levels. These interactions were cg14391855 × cg23937960 (β_interaction  = 0.018, P = 1.87 × 10^-12 ), which mapped to RELA × HLA-G (β_interaction  = 0.218, P = 8.82 × 10^-11 ) and cg08872738 × cg27077312 (β_interaction  = -0.010, P = 1.16 × 10^-11 ), which mapped to TUBA1B × TOMM40 (β_interaction =-0.250, P = 3.83 × 10^-10 ). A trans-omics mediation analysis revealed that 20.3% of epigenetic effects on NSCLC survival were significantly (P = 0.034) mediated through transcriptional expression. These statistically significant trans-omics G × G interactions can also discriminate patients with high risk of mortality. In summary, we identified two G × G interactions at both the epigenetic and transcriptional levels, and our findings may provide potential clues for precision treatment of NSCLC.

Enhancer Function and Evolutionary Roles of Human Accelerated Regions

Human accelerated regions (HARs) are the fastest-evolving sequences in the human genome. When HARs were discovered in 2006, their function was mysterious due to scant annotation of the noncoding genome. Diverse technologies, from transgenic animals to machine learning, have consistently shown that HARs function as gene regulatory enhancers with significant enrichment in neurodevelopment. It is now possible to quantitatively measure the enhancer activity of thousands of HARs in parallel and model how each nucleotide contributes to gene expression. These strategies have revealed that many human HAR sequences function differently than their chimpanzee orthologs, though individual nucleotide changes in the same HAR may have opposite effects, consistent with compensatory substitutions. To fully evaluate the role of HARs in human evolution, it will be necessary to experimentally and computationally dissect them across more cell types and developmental stages.

Genome-wide identification of functional enhancers and their potential roles in pig breeding

Background

The pig is an economically important livestock species and is a widely applied large animal model in medical research. Enhancers are critical regulatory elements that have fundamental functions in evolution, development and disease. Genome-wide quantification of functional enhancers in the pig is needed.

Results

We performed self-transcribing active regulatory region sequencing (STARR-seq) in the porcine kidney epithelial PK15 and testicular ST cell lines, and reliably identified 2576 functional enhancers. Most of these enhancers were located in repetitive sequences and were enriched within silent and lowly expressed genes. Enhancers poorly overlapped with chromatin accessibility regions and were highly enriched in chromatin with the repressive histone modification H3K9me3, which is different from predicted pig enhancers detected using ChIP-seq for H3K27ac or/and H3K4me1 modified histones. This suggests that most pig enhancers identified with STARR-seq are endogenously repressed at the chromatin level and may function during cell type-specific development or at specific developmental stages. Additionally, the PPP3CA gene is associated with the loin muscle area trait and the QKI gene is associated with alkaline phosphatase activity that may be regulated by distal functional enhancers.

Conclusions

In summary, we generated the first functional enhancer map in PK15 and ST cells for the pig genome and highlight its potential roles in pig breeding.

Supplementary Information

The online version contains supplementary material available at 10.1186/s40104-022-00726-y.

Prediction and Experimental Validation of a New Salinity-Responsive Cis-Regulatory Element (CRE) in a Tilapia Cell Line

Transcriptional regulation is a major mechanism by which organisms integrate gene x environment interactions. It can be achieved by coordinated interplay between cis-regulatory elements (CREs) and transcription factors (TFs). Euryhaline tilapia (Oreochromis mossambicus) tolerate a wide range of salinity and thus are an appropriate model to examine transcriptional regulatory mechanisms during salinity stress in fish. Quantitative proteomics in combination with the transcription inhibitor actinomycin D revealed 19 proteins that are transcriptionally upregulated by hyperosmolality in tilapia brain (OmB) cells. We searched the extended proximal promoter up to intron1 of each corresponding gene for common motifs using motif discovery tools. The top-ranked motif identified (STREME1) represents a binding site for the Forkhead box TF L1 (FoxL1). STREME1 function during hyperosmolality was experimentally validated by choosing two of the 19 genes, chloride intracellular channel 2 (clic2) and uridine phosphorylase 1 (upp1), that are enriched in STREME1 in their extended promoters. Transcriptional induction of these genes during hyperosmolality requires STREME1, as evidenced by motif mutagenesis. We conclude that STREME1 represents a new functional CRE that contributes to gene x environment interactions during salinity stress in tilapia. Moreover, our results indicate that FoxL1 family TFs are contribute to hyperosmotic induction of genes in euryhaline fish.

MAP4K4 controlled transcription factor POUM1 regulates PxABCG1 expression influencing Cry1Ac resistance in Plutella xylostella (L.)

Deciphering the molecular mechanisms of insect resistance to Bacillus thuringiensis (Bt) based biotechnology products including Bt sprays and Bt crops is critical for the long-term application of Bt technology. Previously, we established that down-regulation of the ABC transporter gene PxABCG1, trans-regulated by the MAPK signaling pathway, contributed to high-level resistance to Bt Cry1Ac toxin in diamondback moth, Plutella xylostella (L.). However, the underlying transcriptional regulatory mechanism was unknown. Herein, we identified putative binding sites (PBSs) of the transcription factor (TF) POUM1 in the PxABCG1 promoter and used a dual-luciferase reporter assay (DLRA) and yeast one-hybrid (Y1H) assay to reveal that POUM1 activates PxABCG1 via interaction with one of these sites. The expression of POUM1 was significantly decreased in the midgut tissue of Cry1Ac-resistant P. xylostella strains compared to a Cry1Ac-susceptible P. xylostella strain. Silencing of POUM1 expression resulted in reduced expression of the PxABCG1 gene and an increase in larval tolerance to Bt Cry1Ac toxin in the Cry1Ac-susceptible P. xylostella strain. Furthermore, silencing of PxMAP4K4 expression increased the expression of both POUM1 and PxABCG1 genes in the Cry1Ac-resistant P. xylostella strain. These results indicate that the POUM1 induces PxABCG1 expression, while the activated MAPK cascade represses PxABCG1 expression thus reducing Cry1Ac susceptibility in P. xylostella. This result deepens our understanding of the transcriptional regulatory mechanism of midgut Cry receptor genes and the molecular basis of the evolution of Bt resistance in insects.

Transcriptomic Variations and Network Hubs Controlling Seed Size and Weight During Maize Seed Development

To elucidate the mechanisms underlying seed development in maize, comprehensive RNA-seq analyses were conducted on Zhengdan1002 (ZD1002), Zhengdan958 (ZD958), and their parental lines during seven seed developmental stages. We found that gene expression levels were largely nonadditive in hybrids and that cis-only or trans × cis pattern played a large role in hybrid gene regulation during seed developmental stage. Weighted gene co-expression network (WGCNA) analysis showed that 36 modules were highly correlated (r = -0.90-0.92, p < 0.05) with kernel weight, length, and width during seed development. Forty-five transcription factors and 38 ribosomal protein genes were identified as major hub genes determining seed size/weight. We also described a network hub, Auxin Response Factor 12 of maize (ZmARF12), a member of a family of transcription factor that mediate gene expression in response to auxin, potentially links auxin signal pathways, cell division, and the size of the seeds. The ZmARF12 mutant exhibited larger seed size and higher grain weight. ZmARF12 transcription was negatively associated with cell division during seed development, which was confirmed by evaluating the yield of protoplasts that isolated from the kernels of the mutant and other inbred lines. Transient knock-down of ZmARF12 in maize plants facilitated cell expansion and division, whereas transient silencing of its potential interactor ZmIAA8 impaired cell division. ZmIAA8 expression was repressed in the ZmARF12 over-expressed protoplasts. The mutant phenotype and the genetics studies presented here illustrated evidence that ZmARF12 is a cell division repressor, and potentially determines the final seed size.

Evolution of a New Testis-Specific Functional Promoter Within the Highly Conserved Map2k7 Gene of the Mouse

Map2k7 (synonym Mkk7) is a conserved regulatory kinase gene and a central component of the JNK signaling cascade with key functions during cellular differentiation. It shows complex transcription patterns, and different transcript isoforms are known in the mouse (Mus musculus). We have previously identified a newly evolved testis-specific transcript for the Map2k7 gene in the subspecies M. m. domesticus. Here, we identify the new promoter that drives this transcript and find that it codes for an open reading frame (ORF) of 50 amino acids. The new promoter was gained in the stem lineage of closely related mouse species but was secondarily lost in the subspecies M. m. musculus and M. m. castaneus. A single mutation can be correlated with its transcriptional activity in M. m. domesticus, and cell culture assays demonstrate the capability of this mutation to drive expression. A mouse knockout line in which the promoter region of the new transcript is deleted reveals a functional contribution of the newly evolved promoter to sperm motility and the spermatid transcriptome. Our data show that a new functional transcript (and possibly protein) can evolve within an otherwise highly conserved gene, supporting the notion of regulatory changes contributing to the emergence of evolutionary novelties.

Inheritance of gene expression throughout fruit development in chili pepper

Gene expression is the primary molecular phenotype and can be estimated in specific organs or tissues at particular times. Here we analyzed genome-wide inheritance of gene expression in fruits of chili pepper (Capsicum annuum L.) in reciprocal crosses between a domesticated and a wild accession, estimating this parameter during fruit development. We defined a general hierarchical schema to classify gene expression inheritance which can be employed for any quantitative trait. We found that inheritance of gene expression is affected by both, the time of fruit development as well as the direction of the cross, and propose that such variations could be common in many developmental processes. We conclude that classification of inheritance patterns is important to have a better understanding of the mechanisms underlying gene expression regulation, and demonstrate that sets of genes with specific inheritance pattern at particular times of fruit development are enriched in different biological processes, molecular functions and cell components. All curated data and functions for analysis and visualization are publicly available as an R package.

The cis-regulatory effects of modern human-specific variants

The Neanderthal and Denisovan genomes enabled the discovery of sequences that differ between modern and archaic humans, the majority of which are noncoding. However, our understanding of the regulatory consequences of these differences remains limited, in part due to the decay of regulatory marks in ancient samples. Here, we used a massively parallel reporter assay in embryonic stem cells, neural progenitor cells, and bone osteoblasts to investigate the regulatory effects of the 14,042 single-nucleotide modern human-specific variants. Overall, 1791 (13%) of sequences containing these variants showed active regulatory activity, and 407 (23%) of these drove differential expression between human groups. Differentially active sequences were associated with divergent transcription factor binding motifs, and with genes enriched for vocal tract and brain anatomy and function. This work provides insight into the regulatory function of variants that emerged along the modern human lineage and the recent evolution of human gene expression.

Prospects Of Non-Coding Elements In Genomic Dna Based Gene Therapy

Gene therapy has made significant development since the commencement of the first clinical trials a few decades ago and has remained a dynamic area of research regardless of obstacles such as immune response and insertional mutagenesis. Progression in various technologies like next-generation sequencing (NGS) and nanotechnology has established the importance of non-coding segments of a genome, thereby taking gene therapy to the next level. In this review, we have summarized the importance of non-coding elements, highlighting the advantages of using full-length genomic DNA loci (gDNA) compared to complementary DNA (cDNA) or minigene, currently used in gene therapy. The focus of this review is to provide an overview of the advances and the future of potential use of gDNA loci in gene therapy, expanding the therapeutic repertoire in molecular medicine.

Cis and trans effects differentially contribute to the evolution of promoters and enhancers

BACKGROUND

Gene expression differences between species are driven by both cis and trans effects. Whereas cis effects are caused by genetic variants located on the same DNA molecule as the target gene, trans effects are due to genetic variants that affect diffusible elements. Previous studies have mostly assessed the impact of cis and trans effects at the gene level. However, how cis and trans effects differentially impact regulatory elements such as enhancers and promoters remains poorly understood. Here, we use massively parallel reporter assays to directly measure the transcriptional outputs of thousands of individual regulatory elements in embryonic stem cells and measure cis and trans effects between human and mouse.

RESULTS

Our approach reveals that cis effects are widespread across transcribed regulatory elements, and the strongest cis effects are associated with the disruption of motifs recognized by strong transcriptional activators. Conversely, we find that trans effects are rare but stronger in enhancers than promoters and are associated with a subset of transcription factors that are differentially expressed between human and mouse. While we find that cis-trans compensation is common within promoters, we do not see evidence of widespread cis-trans compensation at enhancers. Cis-trans compensation is inversely correlated with enhancer redundancy, suggesting that such compensation may often occur across multiple enhancers.

CONCLUSIONS

Our results highlight differences in the mode of evolution between promoters and enhancers in complex mammalian genomes and indicate that studying the evolution of individual regulatory elements is pivotal to understand the tempo and mode of gene expression evolution.