Function and constraint in enhancers with multiple evolutionary origins

Motivation: Thousands of human gene regulatory enhancers are composed of sequences with multiple evolutionary origins. These evolutionarily "complex" enhancers consist of older "core" sequences and younger "derived" sequences. However, the functional relationship between the sequences of different evolutionary origins within complex enhancers is poorly understood. Results: We evaluated the function, selective pressures, and sequence variation across core and derived components of human complex enhancers. We find that both components are older than expected from the genomic background, and cores are enriched for derived sequences of similar evolutionary ages. Both components show strong evidence of biochemical activity in massively parallel report assays (MPRAs). However, core and derived sequences have distinct transcription factor (TF) binding preferences that are largely stable across evolutionary origins. Given these signatures of function, both core and derived sequences have substantial evidence of purifying selection. Nonetheless, derived sequences exhibit weaker purifying selection than adjacent cores. Derived sequences also tolerate more common genetic variation and are enriched compared to cores for eQTL associated with gene expression variability in human populations. Conclusions: Both core and derived sequences have strong evidence of gene regulatory function, but derived sequences have distinct constraint profiles, TF binding preferences, and tolerance to variation compared with cores. We propose that the step-wise integration of younger derived and older core sequences has generated regulatory substrates with robust activity and the potential for functional variation. Our analyses demonstrate that synthesizing study of enhancer evolution and function can aid interpretation of regulatory sequence activity and functional variation across human populations.

Using Machine Learning Approaches to Predict Target Gene Expression in Rice T-DNA Insertional Mutants

To change the expression of the flanking genes by inserting T-DNA into the genome is commonly used in rice functional gene research. However, whether the expression of a gene of interest is enhanced must be validated experimentally. Consequently, to improve the efficiency of screening activated genes, we established a model to predict gene expression in T-DNA mutants through machine learning methods. We gathered experimental datasets consisting of gene expression data in T-DNA mutants and captured the PROMOTER and MIDDLE sequences for encoding. In first-layer models, support vector machine (SVM) models were constructed with nine features consisting of information about biological function and local and global sequences. Feature encoding based on the PROMOTER sequence was weighted by logistic regression. The second-layer models integrated 16 first-layer models with minimum redundancy maximum relevance (mRMR) feature selection and the LADTree algorithm, which were selected from nine feature selection methods and 65 classified methods, respectively. The accuracy of the final two-layer machine learning model, referred to as TIMgo, was 99.3% based on fivefold cross-validation, and 85.6% based on independent testing. We discovered that the information within the local sequence had a greater contribution than the global sequence with respect to classification. TIMgo had a good predictive ability for target genes within 20 kb from the 35S enhancer. Based on the analysis of significant sequences, the G-box regulatory sequence may also play an important role in the activation mechanism of the 35S enhancer.

Wheat Breeding through Genetic and Physical Mapping 2

Following the success of the first topic, the special issue of “Wheat breeding through genetic and physical mapping 2” has been re-proposed in order to keep current the recent advancement in research on genetic and physical mapping of candidate genes for agronomically important traits, in studies of the regulatory sequence for biotic and abiotic stress resistance [...]

Conservation of A-to-I RNA editing in bowhead whale and pig

RNA editing is a post-transcriptional process in which nucleotide changes are introduced into an RNA sequence, many of which can contribute to proteomic sequence variation. The most common type of RNA editing, contributing to nearly 99% of all editing events in RNA, is A-to-I (adenosine-to-inosine) editing mediated by double-stranded RNA-specific adenosine deaminase (ADAR) enzymes. A-to-I editing at ‘recoding’ sites results in non-synonymous substitutions in protein-coding sequences. Here, we present studies of the conservation of A-to-I editing in selected mRNAs between pigs, bowhead whales, humans and two shark species. All examined mRNAs–NEIL1, COG3, GRIA2, FLNA, FLNB, IGFBP7, AZIN1, BLCAP, GLI1, SON, HTR2C and ADAR2 –showed conservation of A-to-I editing of recoding sites. In addition, novel editing sites were identified in NEIL1 and GLI1 in bowhead whales. The A-to-I editing site of human NEIL1 in position 242 was conserved in the bowhead and porcine homologues. A novel editing site was discovered in Tyr244. Differential editing was detected at the two adenosines in the NEIL1 242 codon in both pig and bowhead NEIL1 mRNAs in various tissues and organs. No conservation of editing of KCNB1 and EEF1A mRNAs was seen in bowhead whales. In silico analyses revealed conservation of five adenosines in ADAR2, some of which are subject to A-to-I editing in bowheads and pigs, and conservation of a regulatory sequence in GRIA2 mRNA that is responsible for recognition of the ADAR editing enzyme.

MaveDB v2: a curated community database with over three million variant effects from multiplexed functional assays

A central problem in genomics is understanding the effect of individual DNA variants. Multiplexed Assays of Variant Effect (MAVEs) can help address this challenge by measuring all possible single nucleotide variant effects in a gene or regulatory sequence simultaneously. Here we describe MaveDB v2, which has become the database of record for MAVEs. MaveDB now contains a large fraction of published studies, comprising over two hundred datasets and three million variant effect measurements. We created tools and APIs to streamline data submission and access, transforming MaveDB into a hub for the analysis and dissemination of these impactful datasets.

An evolutionarily conserved cis-regulatory element of Nkx3.2 contributes to early jaw joint morphology in zebrafish

The acquisition of movable jaws was a major event during vertebrate evolution. The role of NK3 homeobox 2 (Nkx3.2) transcription factor in patterning the primary jaw joint of gnathostomes (jawed vertebrates) is well known, however knowledge about its regulatory mechanism is lacking. In this study, we report a proximal enhancer element of Nkx3.2 that is deeply conserved in gnathostomes but undetectable in the jawless hagfish. This enhancer is active in the developing jaw joint region of the zebrafish Danio rerio, and was thus designated as jaw joint regulatory sequence 1 (JRS1). We further show that JRS1 enhancer sequences from a range of gnathostome species, including a chondrichthyan and mammals, have the same activity in the jaw joint as the native zebrafish enhancer, indicating a high degree of functional conservation despite the divergence of cartilaginous and bony fish lineages or the transition of the primary jaw joint into the middle ear of mammals. Finally, we show that deletion of JRS1 from the zebrafish genome using CRISPR/Cas9 leads to a transient jaw joint deformation and partial fusion. Emergence of this Nkx3.2 enhancer in early gnathostomes may have contributed to the origin and shaping of the articulating surfaces of vertebrate jaws.

Alu insertion variants alter gene transcript levels

Alu are high copy number interspersed repeats that have accumulated near genes during primate and human evolution. They are a pervasive source of structural variation in modern humans. Impacts that Alu insertions may have on gene expression are not well understood, although some have been associated with expression quantitative trait loci (eQTLs). Here, we directly test regulatory effects of polymorphic Alu insertions in isolation of other variants on the same haplotype. To screen insertion variants for those with such effects, we used ectopic luciferase reporter assays and evaluated 110 Alu insertion variants, including more than 40 with a potential role in disease risk. We observed a continuum of effects with significant outliers that up- or down-regulate luciferase activity. Using a series of reporter constructs, which included genomic context surrounding the Alu, we can distinguish between instances in which the Alu disrupts another regulator and those in which the Alu introduces new regulatory sequence. We next focused on three polymorphic Alu loci associated with breast cancer that display significant effects in the reporter assay. We used CRISPR to modify the endogenous sequences, establishing cell lines varying in the Alu genotype. Our findings indicate that Alu genotype can alter expression of genes implicated in cancer risk, including PTHLH, RANBP9, and MYC. These data show that commonly occurring polymorphic Alu elements can alter transcript levels and potentially contribute to disease risk.

TCF7L1 regulates cytokine response and neuroendocrine differentiation of prostate cancer

Neuroendocrine differentiation (NED) is associated with WNT signaling activation and can be significantly observed after failure of androgen-deprivation therapy (ADT) for prostatic adenocarcinomas. Cytokine signaling is stimulated in NED prostate cancer; however, how ADT-upregulated WNT signaling promotes activation of cytokine signaling and contributes to NED of prostate cancer is poorly understood. In this study, we identified ADT-mediated upregulation of transcription factor 7 like 1 (TCF7L1), which increases the cytokine response and enhances NED of prostate cancer through interleukin (IL)-8/C-X-C motif chemokine receptor type 2 (CXCR2) signaling activation. ADT induced the secretion of WNT4 which upon engagement of TCF7L1 in prostate cancer cells, enhanced IL-8 and CXCR2 expressions. TCF7L1 directly binds to the regulatory sequence region of IL-8 and CXCR2 through WNT4 activation, thus upregulating IL-8/CXCR2 signaling-driven NED and cell motility. Analysis of prostate tissue samples collected from small-cell neuroendocrine prostate cancer (SCPC) and castration-resistant prostate cancer (CRPC) tumors showed an increased intensity of nuclear TCF7L1 associated with CXCR2. Our results suggest that induction of WNT4/TCF7L1 results in increased NED and malignancy in prostate cancer that is linked to dysregulation of androgen receptor signaling and activation of the IL-8/CXCR2 pathway.

Long-range regulatory effects of Neandertal DNA in modern humans

The admixture between modern humans and Neandertals has resulted in ~2% of the genomes of present-day non-Africans still being composed of Neandertal DNA. Association studies have shown that introgressed DNA still significantly influences skin and hair traits, immunity and behavioral phenotypes in people today. Several of the phenotype-associated archaic variants had links to regulatory effects as well. In general, analyses of allele-specific expression, regulatory sequence composition and cis-eQTL have demonstrated a significant contribution of this introgressed DNA to the transcriptomic landscape in people today. However, little is known about the impact of Neandertal DNA on trans-eQTLs - long-range regulatory effects that have been shown to explain ~20% of expression variation. Here, we used blood eQTL results from >30,000 individuals from eQTLGen Consortium. The cohort size allowed for a robust identification of trans-eQTLs and in addition enabled quantifying the role of transcription factors (TF) in mediating long-range regulatory effects. In our study we used this information to (i) annotate trans-eQTLs that are linked to Neandertal variants and (ii) predict long-range regulatory effects that are induced by Neandertal DNA by screening for the predicted target genes of TFs that are cis-eQTLs linked to Neandertal variants. We show that both trans-eQTL-associated Neandertal variants and those predicted to have long-range regulatory effects affect genes in genomic regions devoid of Neandertal DNA. In addition, both types of variants included candidates for local adaptation and show associations with autoimmune disorders, a severe Covid-19 phenotype, blood cell type composition and anthropometric measures. Our results suggest that the regulatory reach of Neandertal DNA goes beyond the 40% of genomic sequence that it still covers in present-day non-Africans and that via this mechanism Neandertal DNA additionally influences the phenotypic variation in people today.

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

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