Functional characterization of SIM1-associated enhancers

Rare variation in non-coding regions with evolutionary signatures contributes to autism spectrum disorder risk

Little is known about the role of non-coding regions in the etiology of autism spectrum disorder (ASD). We examined three classes of non-coding regions: human accelerated regions (HARs), which show signatures of positive selection in humans; experimentally validated neural VISTA enhancers (VEs); and conserved regions predicted to act as neural enhancers (CNEs). Targeted and whole-genome analysis of >16,600 samples and >4,900 ASD probands revealed that likely recessive, rare, inherited variants in HARs, VEs, and CNEs substantially contribute to ASD risk in probands whose parents share ancestry, which enriches for recessive contributions, but modestly contribute, if at all, in simplex family structures. We identified multiple patient variants in HARs near IL1RAPL1 and in VEs near OTX1 and SIM1 and showed that they change enhancer activity. Our results implicate both human-evolved and evolutionarily conserved non-coding regions in ASD risk and suggest potential mechanisms of how regulatory changes can modulate social behavior.

Identification of ancestral gnathostome Gli3 enhancers with activity in mammals

Abnormal expression of the transcriptional regulator and hedgehog (Hh) signaling pathway effector Gli3 is known to trigger congenital disease, most frequently affecting the central nervous system (CNS) and the limbs. Accurate delineation of the genomic cis-regulatory landscape controlling Gli3 transcription during embryonic development is critical for the interpretation of noncoding variants associated with congenital defects. Here, we employed a comparative genomic analysis on fish species with a slow rate of molecular evolution to identify seven previously unknown conserved noncoding elements (CNEs) in Gli3 intronic intervals (CNE15-21). Transgenic assays in zebrafish revealed that most of these elements drive activities in Gli3 expressing tissues, predominantly the fins, CNS, and the heart. Intersection of these CNEs with human disease associated SNPs identified CNE15 as a putative mammalian craniofacial enhancer, with conserved activity in vertebrates and potentially affected by mutation associated with human craniofacial morphology. Finally, comparative functional dissection of an appendage-specific CNE conserved in slowly evolving fish (elephant shark), but not in teleost (CNE14/hs1586) indicates co-option of limb specificity from other tissues prior to the divergence of amniotes and lobe-finned fish. These results uncover a novel subset of intronic Gli3 enhancers that arose in the common ancestor of gnathostomes and whose sequence components were likely gradually modified in other species during the process of evolutionary diversification.

Rare variation in noncoding regions with evolutionary signatures contributes to autism spectrum disorder risk

Little is known about the role of noncoding regions in the etiology of autism spectrum disorder (ASD). We examined three classes of noncoding regions: Human Accelerated Regions (HARs), which show signatures of positive selection in humans; experimentally validated neural Vista Enhancers (VEs); and conserved regions predicted to act as neural enhancers (CNEs). Targeted and whole genome analysis of >16,600 samples and >4900 ASD probands revealed that likely recessive, rare, inherited variants in HARs, VEs, and CNEs substantially contribute to ASD risk in probands whose parents share ancestry, which enriches for recessive contributions, but modestly, if at all, in simplex family structures. We identified multiple patient variants in HARs near IL1RAPL1 and in a VE near SIM1 and showed that they change enhancer activity. Our results implicate both human-evolved and evolutionarily conserved noncoding regions in ASD risk and suggest potential mechanisms of how changes in regulatory regions can modulate social behavior.

Regeneration and developmental enhancers are differentially compatible with minimal promoters

Enhancers and promoters are cis-regulatory elements that control gene expression. Enhancers are activated in a cell type-, tissue-, and condition-specific manner to stimulate promoter function and transcription. Zebrafish have emerged as a powerful animal model for examining the activities of enhancers derived from various species through transgenic enhancer assays, in which an enhancer is coupled with a minimal promoter. However, the efficiency of minimal promoters and their compatibility with multiple developmental and regeneration enhancers have not been systematically tested in zebrafish. Thus, we assessed the efficiency of six minimal promoters and comprehensively interrogated the compatibility of the promoters with developmental and regeneration enhancers. We found that the fos minimal promoter and Drosophila synthetic core promoter (DSCP) yielded high rates of leaky expression that may complicate the interpretation of enhancer assays. Notably, the adenovirus E1b promoter, the zebrafish lepb 0.8-kb (P0.8) and lepb 2-kb (P2) promoters, and a new zebrafish synthetic promoter (ZSP) that combines elements of the E1b and P0.8 promoters drove little or no ectopic expression, making them suitable for transgenic assays. We also found significant differences in compatibility among specific combinations of promoters and enhancers, indicating the importance of promoters as key regulatory elements determining the specificity of gene expression. Our study provides guidelines for transgenic enhancer assays in zebrafish to aid in the discovery of functional enhancers regulating development and regeneration.

Clinical Significance of ADAMTS19, BMP7, SIM1, and SFRP1 Promoter Methylation in Renal Clear Cell Carcinoma

Background

Clear cell renal cell carcinoma (ccRCC) is the most common subtype of kidney tumors, accounting for the majority of deaths from genitourinary cancers. The currently used nomograms for predicting patient outcomes are based on clinical-pathological characteristics only; however, a significant number of ccRCC survivors with similar radiological and histological features still demonstrate a different clinical course of the disease. This study aimed at the identification of novel DNA methylation biomarkers for the monitoring of patients with ccRCC.

Methods

Gene expression profiling by SurePrint G3 Human GE 8×60K Microarrays was performed in 4 ccRCC tissues and adjacent non-cancerous renal tissue (NRT) samples. Four down-regulated genes were selected for further DNA methylation status analysis in 123 ccRCC and 45 NRT samples using methylation-specific PCR (MSP).

Results

DNA methylation changes of ADAMTS19, BMP7, SIM1, and SFRP1 were cancer-specific with significantly (P<0.050) higher methylation frequency (37%, 20%, 18%, and 42%, respectively) in tumor tissues. The methylated status of at least one gene was significantly related to various clinical-pathological parameters, including tumor size, Fuhrman and WHO/ISUP grades, intravascular invasion, and necrosis. Moreover, the methylated status of multimarker panel ADAMTS19, BMP7 & SFRP1 was predictive for poorer overall survival (HR, 4.11; 95% CI, 1.22–13.86).

Conclusion

In conclusion, DNA methylation of the three-gene panel consisting of ADAMTS19, BMP7 & SFRP1 supposedly predicts the outcome of patients diagnosed with ccRCC and possibly might be used to enrich the current prognostic tools.

Comparative genomic analysis of human GLI2 locus using slowly evolving fish revealed the ancestral gnathostome set of early developmental enhancers

Background

The zinc finger‐containing transcription factor Gli2, is a key mediator of Hedgehog (Hh) signaling and participates in embryonic patterning of various organs including the central nervous system (CNS) and limbs. Abnormal expression of Gli2 can impede the transcription of Hh target genes through disruption of proper balance between Gli2 and Gli3 functions. Therefore, delineation of enhancers that are required for complementary roles of Glis would allow the interrogation of those pathogenic variants that cause gene dysregulation, and a corresponding abnormal phenotype. Previously, we reported tissue‐specific enhancers for Gli family including Gli2 through direct tetrapod‐teleost comparisons.

Results

Here, we employed the sequence alignments of slowly evolving spotted gar and elephant shark and have identified six novel conserved noncoding elements in human GLI2 containing locus. Zebrafish‐based transgenic assays revealed that combined action of these autonomous CNEs reflects many aspects of Gli2 specific endogenous transcriptional activity, including CNS and pectoral fins.

Conclusion

Taken together with our previous findings, this study suggests that Hh‐signaling controlled deployment of Gli2 activity in embryonic patterning arose in the common ancestor of gnathostomes. These GLI2 specific cis‐regulatory modules will help to identify DNA variants that probably reside outside of coding intervals and are associated with congenital anomalies.

We performed a phylogenetic footprint analyses of human GLI2 containing locus by incorporating relatively slowly evolving gar and elephant shark genomes and have identified multiple novel conserved non‐coding elements (CNEs) that were not predicted by direct human‐teleostcomparisons.

Comparative analyses suggest that majority of the GLI2 associated CNEs identified in the present data and reported previously arose in the common ancestor of gnathostomes but lost in teleosts, presumably because of fast teleost sequence evolution.

Functional testing of GLI2 associated CNEs by employing zebrafish based transgenic reporter assays revealed their tissue specific cis‐regulatory potential that corresponds with the results based on whole‐mount in situ hybridization analysis of gli2 mRNA in zebrafish.

The delineated set of GLI2 associated enhancers can be further interrogated to determine their role in canonical Hh signaling, gene dysregulation, and a corresponding congenital anomaly.

CRISPR-mediated activation of a promoter or enhancer rescues obesity caused by haploinsufficiency

A wide range of human diseases result from haploinsufficiency, where the function of one of the two gene copies is lost. Here, we targeted the remaining functional copy of a haploinsufficient gene using CRISPR-mediated activation (CRISPRa) in Sim1 and Mc4r heterozygous mouse models to rescue their obesity phenotype. Transgenic-based CRISPRa targeting of the Sim1 promoter or its distant hypothalamic enhancer up-regulated its expression from the endogenous functional allele in a tissue-specific manner, rescuing the obesity phenotype in Sim1 heterozygous mice. To evaluate the therapeutic potential of CRISPRa, we injected CRISPRa-recombinant adeno-associated virus into the hypothalamus, which led to reversal of the obesity phenotype in Sim1 and Mc4r haploinsufficient mice. Our results suggest that endogenous gene up-regulation could be a potential strategy to treat altered gene dosage diseases.

Genetic variation in the SIM1 locus is associated with erectile dysfunction

Erectile dysfunction is a common condition of men in middle and older ages. Twin studies suggest that about one-third of the risk is due to genetic factors, independent of other known erectile dysfunction risk factors. However, studies that have searched for specific genetic contributors have been limited due to small sample sizes, candidate gene approaches, and weak phenotyping. As a result, there are no confirmed genetic risk factors for erectile dysfunction. This study finds a specific genetic cause for erectile dysfunction.

Erectile dysfunction affects millions of men worldwide. Twin studies support the role of genetic risk factors underlying erectile dysfunction, but no specific genetic variants have been identified. We conducted a large-scale genome-wide association study of erectile dysfunction in 36,649 men in the multiethnic Kaiser Permanente Northern California Genetic Epidemiology Research in Adult Health and Aging cohort. We also undertook replication analyses in 222,358 men from the UK Biobank. In the discovery cohort, we identified a single locus (rs17185536-T) on chromosome 6 near the single-minded family basic helix-loop-helix transcription factor 1 (SIM1) gene that was significantly associated with the risk of erectile dysfunction (odds ratio = 1.26, P = 3.4 × 10⁻²⁵). The association replicated in the UK Biobank sample (odds ratio = 1.25, P = 6.8 × 10⁻¹⁴), and the effect is independent of known erectile dysfunction risk factors, including body mass index (BMI). The risk locus resides on the same topologically associating domain as SIM1 and interacts with the SIM1 promoter, and the rs17185536-T risk allele showed differential enhancer activity. SIM1 is part of the leptin–melanocortin system, which has an established role in body weight homeostasis and sexual function. Because the variants associated with erectile dysfunction are not associated with differences in BMI, our findings suggest a mechanism that is specific to sexual function.

Aberrant DNA methylation of DLX4 and SIM1 is a predictive marker for disease progression of uterine cervical low-grade squamous intraepithelial lesion

BACKGROUND

A few uterine cervical low-grade squamous intraepithelial lesions (LSILs) are known to progress with high-risk human papillomavirus (hrHPV).

METHODS

One hundred and thirteen patients were classified into four groups according to their cervical cytology, hrHPV infection, and follow up. Cytology samples were examined for aberrant DNA methylation of DLX4 and SIM1 genes and protein expressions. CaSki cells were treated with 5-Aza-2'-deoxycytidine (5-aza-dC).

RESULTS

Group 1 was negative for intraepithelial lesions or malignancies. LSIL in group 2 showed a continuance of LSIL for longer than 365 days and LSIL in group 3 showed an upgrading to high-grade (H) SIL or higher (HSIL+) within 365 days of LSIL diagnosis. Group 4 was squamous cell carcinoma. All but group 1 were infected with hrHPV. Significant differences existed in the frequency of DNA methylation between groups 2 and 3 (p = 0.044), between groups 3 and 4 (p = 0.020) for DLX4, and between groups 1 and 3 (p = 0.0003), and groups 2 and 3 (p = 0.005) for the SIM1 gene. DLX4 protein expression was significantly reduced in the DLX4 methylation positive tissues (p = 0.008). The 5-aza-dC treatment restored DLX4 mRNA expressions of CaSki cells (p < 0.005). The LSIL cases with DNA methylation of the SIM1 gene, or both genes, progressed faster to HSIL+ than did the others (p = 0.033 and p = 0.045, respectively).

CONCLUSION

Aberrant DNA methylation of the DLX4 and SIM1 genes should be a novel progression marker for uterine cervical LSIL with hrHPV infection.

Incomplete penetrance and phenotypic variability of 6q16 deletions including SIM1

6q16 deletions have been described in patients with a Prader-Willi-like (PWS-like) phenotype. Recent studies have shown that certain rare single-minded 1 (SIM1) loss-of-function variants were associated with a high intra-familial risk for obesity with or without features of PWS-like syndrome. Although SIM1 seems to have a key role in the phenotype of patients carrying 6q16 deletions, some data support a contribution of other genes, such as GRIK2, to explain associated behavioural problems. We describe 15 new patients in whom de novo 6q16 deletions were characterised by comparative genomic hybridisation or single-nucleotide polymorphism (SNP) array analysis, including the first patient with fetopathological data. This fetus showed dysmorphic facial features, cerebellar and cerebral migration defects with neuronal heterotopias, and fusion of brain nuclei. The size of the deletion in the 14 living patients ranged from 1.73 to 7.84 Mb, and the fetus had the largest deletion (14 Mb). Genotype-phenotype correlations confirmed the major role for SIM1 haploinsufficiency in obesity and the PWS-like phenotype. Nevertheless, only 8 of 13 patients with SIM1 deletion exhibited obesity, in agreement with incomplete penetrance of SIM1 haploinsufficiency. This study in the largest series reported to date confirms that the PWS-like phenotype is strongly linked to 6q16.2q16.3 deletions and varies considerably in its clinical expression. The possible involvement of other genes in the 6q16.2q16.3-deletion phenotype is discussed.