Functional annotation and investigation of the 10q24.33 melanoma risk locus identifies a common variant that influences transcriptional regulation of OBFC1.

The 10q24.33 locus is known to be associated with susceptibility to cutaneous malignant melanoma (CMM), but the mechanisms underlying this association have been not extensively investigated. We carried out an integrative genomic analysis of 10q24.33 using epigenomic annotations and in vitro reporter gene assays to identify regulatory variants. We found two putative functional single nucleotide polymorphisms (SNPs) in an enhancer and in the promoter of OBFC1, respectively, in neural crest and CMM cells, one, rs2995264, altering enhancer activity. The minor allele G of rs2995264 correlated with lower OBFC1 expression in 470 CMM tumors and was confirmed to increase the CMM risk in a cohort of 484 CMM cases and 1801 controls of Italian origin. Hi-C and chromosome conformation capture (3C) experiments showed the interaction between the enhancer-SNP region and the promoter of OBFC1 and an isogenic model characterized by CRISPR-Cas9 deletion of the enhancer-SNP region confirmed the potential regulatory effect of rs2995264 on OBFC1 transcription. Moreover, the presence of G-rs2995264 risk allele reduced the binding affinity of the transcription factor MEOX2. Biologic investigations showed significant cell viability upon depletion of OBFC1, specifically in CMM cells that were homozygous for the protective allele. Clinically, high levels of OBFC1 expression associated with histologically favorable CMM tumors. Finally, preliminary results suggested the potential effect of decreased OBFC1 expression on telomerase activity in tumorigenic conditions. Our results support the hypothesis that reduced expression of OBFC1 gene through functional heritable DNA variation can contribute to malignant transformation of normal melanocytes.

Evidence for a Pathogenic Role of CSMD1 in Childhood Apraxia of Speech.

Childhood apraxia of speech (CAS) is a pediatric motor speech disorder. The genetic etiology of this complex neurological condition is not yet well understood, although some genes have been linked to it. We describe the case of a boy with a severe and persistent motor speech disorder, consistent with CAS, and a coexisting language impairment.Whole exome sequencing in our case revealed a de novo and splicing mutation in the CSMD1 gene.

Neural crest-derived tumor neuroblastoma and melanoma share 1p13.2 as susceptibility locus that shows a long-range interaction with the SLC16A1 gene.

Neuroblastoma (NB) and malignant cutaneous melanoma (CMM) are neural crest cells (NCC)-derived tumors and may have a shared genetic basis, but this has not been investigated systematically by genome-wide association studies (GWAS). We took a three-staged approach to conduct cross-disease meta-analysis of GWAS for NB and CMM (2101 NB cases and 4202 controls; 12 874 CMM cases and 23 203 controls) to identify shared loci. Findings were replicated in 1403 NB cases and 1403 controls of European ancestry and in 636 NB, 508 CMM cases and 2066 controls of Italian origin. We found a cross-association at locus 1p13.2 (rs2153977, odds ratio = 0.91, P = 5.36 × 10-8). We also detected a suggestive (P < 10-7) NB-CMM cross-association at 2q37.1 with opposite effect on cancer risk. Pathway analysis of 110 NB-CMM risk loci with P < 10-4 demonstrated enrichment of biological processes such as cell migration, cell cycle, metabolism and immune response, which are essential of human NCC development, underlying both tumors. In vitro and in silico analyses indicated that the rs2153977-T protective allele, located in an NB and CMM enhancer, decreased expression of SLC16A1 via long-range loop formation and altered a T-box protein binding site. Upon depletion of SLC16A1, we observed a decrease of cellular proliferation and invasion in both NB and CMM cell lines, suggesting its role as oncogene. This is the largest study to date examining pleiotropy across two NC cell-derived tumors identifying 1p13.2 as common susceptibility locus for NB and CMM risk. We demonstrate that combining genome-wide association studies results across cancers with same origins can identify new loci common to neuroblastoma and melanoma arising from tissues which originate from neural crest cells. Our results also show 1p13.2 confer risk to neuroblastoma and melanoma by regulating SLC16A1.

Kinome multigenic panel identified novel druggable EPHB4-V871I somatic variant in high-risk neuroblastoma.

Neuroblastoma (NB) is the most common extracranial neoplasm in children. The overall outcome for high-risk NB patients is still unacceptable, therefore, it is critical to deeply understand molecular mechanisms associated with NB, which in turn can be utilized for developing drugs towards the treatment of NB. Protein kinases (TKs) play an essential role in the regulation of cell survival and proliferation. Different kinases, such as anaplastic lymphoma kinase (ALK), Aurora kinase, RET receptor tyrosine kinase, are potential therapeutic targets in various cancers, including NB. We analysed a cohort of 45 high-risk NB patients and 9 NB cell lines by a targeted-(t)NGS custom gene panel (genes codifying for the kinase domains of 90 TKs). We identified somatic variants in four TK genes (ALK, EPHB4, LMTK3 and EPHB6) in NB patients and we functionally characterized an interesting somatic variant, V871I, in EPHB4 gene. EPHB4 plays a crucial role in cardiovascular development and regulates vascularization in cancer-promoting angiogenesis, tumour growth and metastasis. Several EPHB4 mutations have previously been identified in solid and haematological tumour specimens but EPHB4 mutations were not described until now in NB. Interestingly, a re-analysis of public CGH-array showed that the EPHB4 gain is associated with advanced diseases in NB. We further demonstrated that higher EPHB4 expression is correlated to stage 4 of NB and with poor overall survival. Additionally, we also revealed that the EPHB4-V871I accounts for increased proliferation, migration and invasion properties in two NB cell lines by acting on VEGF, c-RAF and CDK4 target genes and by increasing the phosphorylation of ERK1-2 pathway. The use of two EPHB4 inhibitors, JI-101 and NVP-BHG712, was able to rescue the phenotype driven by the variant. Our study suggested that EPHB4 is a promising therapeutic target in high-risk NB.

ZNF687 Mutations in Severe Paget Disease of Bone Associated with Giant Cell Tumor.

Paget disease of bone (PDB) is a skeletal disorder characterized by focal abnormalities of bone remodeling, which result in enlarged and deformed bones in one or more regions of the skeleton. In some cases, the pagetic tissue undergoes neoplastic transformation, resulting in osteosarcoma and, less frequently, in giant cell tumor of bone (GCT). We performed whole-exome sequencing in a large family with 14 PDB-affected members, four of whom developed GCT at multiple pagetic skeletal sites, and we identified the c.2810C>G (p.Pro937Arg) missense mutation in the zinc finger protein 687 gene (ZNF687). The mutation precisely co-segregated with the clinical phenotype in all affected family members. The sequencing of seven unrelated individuals with GCT associated with PDB (GCT/PDB) identified the same mutation in all individuals, unravelling a founder effect. ZNF687 is highly expressed during osteoclastogenesis and osteoblastogenesis and is dramatically upregulated in the tumor tissue of individuals with GCT/PDB. Interestingly, our preliminary findings showed that ZNF687, indicated as a target gene of the NFkB transcription factor by ChIP-seq analysis, is also upregulated in the peripheral blood of PDB-affected individuals with (n = 5) or without (n = 6) mutations in SQSTM1, encouraging additional studies to investigate its potential role as a biomarker of PDB risk.

Fine mapping of 2q35 high-risk neuroblastoma locus reveals independent functional risk variants and suggests full-length BARD1 as tumor-suppressor.

A previous genome-wide association study (GWAS) identified common variation at the BARD1 locus as being highly associated with susceptibility to high-risk neuroblastoma, but the mechanisms underlying this association have been not extensively investigated. Here, we performed a fine mapping analysis of BARD1 locus (2q35) using GWAS data from 556 high-risk neuroblastoma patients and 2,575 controls of European-American ancestry, and identified two independent genome-wide neuroblastoma-associated loci. Functional single-nucleotide polymorphism (SNP) prioritization identified two causative variants that independently contributed to neuroblastoma risk, and each replicated robustly in multiple independent cohorts comprising 445 high-risk cases and 3,170 controls (rs17489363: combined p = 1.07 × 10-31 , OR:1.79, 95% CI:1.62-1.98 and rs1048108: combined p = 7.27 × 10-14 , OR:0.65, 95% CI:0.58-0.73). Particularly, the T risk allele of rs17489363 in the canonical promoter region of full-length BARD1 altered binding site of the transcription factor HSF1 and correlated with low expression of full-length BARD1 mRNA and protein. Low-level expression of full-length BARD1 associated with advanced neuroblastoma. In human neuroblastoma cells, attenuating full-length BARD1 increased proliferation and invasion capacity. In conclusion, we have identified two potentially causative SNPs at the BARD1 locus associated with predisposition to high-risk neuroblastoma, and have shown that full-length BARD1 may act as tumor suppressor.

The functional variant rs34330 of CDKN1B is associated with risk of neuroblastoma.

The genetic aetiology of sporadic neuroblastoma is still largely unknown. We have identified diverse neuroblastoma susceptibility loci by genomewide association studies (GWASs); however, additional SNPs that likely contribute to neuroblastoma susceptibility prompted this investigation for identification of additional variants that are likely hidden among signals discarded by the multiple testing corrections used in the analysis of genomewide data. There is evidence suggesting the CDKN1B, coding for the cycle inhibitor p27Kip1, is involved in neuroblastoma. We thus assess whether genetic variants of CDKN1B are associated with neuroblastoma. We imputed all possible genotypes across CDKN1B locus on a discovery case series of 2101 neuroblastoma patients and 4202 genetically matched controls of European ancestry. The most significantly associated rs34330 was analysed in an independent Italian cohort of 311 cases and 709 controls. In vitro functional analysis was carried out in HEK293T and in neuroblastoma cell line SHEP-2, both transfected with pGL3-CDKN1B-CC or pGL3-CDKN1B-TT constructs. We identified an association of the rs34330 T allele (-79C/T) with the neuroblastoma risk (Pcombined = 0.002; OR = 1.17). The risk allele (T) of this single nucleotide polymorphism led to a lower transcription rate in cells transfected with a luciferase reporter driven by the polymorphic p27Kip1 promoter (P < 0.05). Three independent sets of neuroblastoma tumours carrying -79TT genotype showed a tendency towards lower CDKN1B mRNA levels. Our study shows that a functional variant, associated with a reduced CDKN1B gene transcription, influences neuroblastoma susceptibility.

An 18 gene expression-based score classifier predicts the clinical outcome in stage 4 neuroblastoma.

BACKGROUND: The prognosis of children with metastatic stage 4 neuroblastoma (NB) has remained poor in the past decade. PATIENTS AND METHODS: Using microarray analyses of 342 primary tumors, we here developed and validated an easy to use gene expression-based risk score including 18 genes, which can robustly predict the outcome of stage 4 patients. RESULTS: This classifier was a significant predictor of overall survival in two independent validation cohorts [cohort 1 (n = 214): P = 6.3 × 10(-5); cohort 2 (n = 27): P = 3.1 × 10(-2)]. The prognostic value of the risk score was validated by multivariate analysis including the established markers age and MYCN status (P = 0.027). In the pooled validation cohorts (n = 241), integration of the risk score with the age and/or MYCN status identified subgroups with significantly differing overall survival (ranging from 35 to 100 %). CONCLUSION: Together, the 18-gene risk score classifier can identify patients with stage 4 NB with favorable outcome and may therefore improve risk assessment and treatment stratification of NB patients with disseminated disease.

CFDP1 is a neuroblastoma susceptibility gene that regulates transcription factors of the noradrenergic cell identity.

Pleiotropic genetic factors (e.g., DNA polymorphisms) may be involved in the initiation of neuroblastoma (NB) and coronary artery disease (CAD) given their common origin from defects in neural crest development. To discover novel NB susceptibility genes, we conducted a three-stage survey including a meta-analysis of NB and CAD genome-wide association data, prioritization of NB causal variants, and validation in an independent cohort of affected individuals-control subjects. The lead SNP, rs13337397 at the 16q23.1 locus, associated with both diseases in the meta-analysis and with NB in the validation study. All the SNPs in linkage disequilibrium with rs13337397 were annotated using the H3K27ac epigenetic marker of neural crest cells (NCC) and NB cell lines. Indeed, we identified the functional SNP rs13337017, mapping within an enhancer of NCCs and NB cell lines and showing long-range interactions with CFDP1 by Hi-C analysis. Luciferase assays indicated that the risk allele of rs13337017 increased CFDP1 expression in NB cell lines. Of note, CFDP1 high expression associated with unfavorable prognostic markers in an analysis including 498 NB transcriptomes. Moreover, depletion of CFDP1 markedly decreased viability and migration and increased apoptotic rates in NB cell lines. Finally, transcriptome and qPCR analyses revealed that the depletion of CFDP1 may affect noradrenergic neuron differentiation by downregulating master regulators of sympathetic noradrenergic identity, including PHOX2B, HAND2, and GATA3. Our data strongly suggest that CFDP1 acts as oncogene in NB. In addition, we provide evidence that genetic predisposition to NB can be mediated by the alteration of noradrenergic lineage-specific gene expression.

A pathogenic variant in the FLCN gene presenting with pure dementia: is autophagy at the intersection between neurodegeneration and cancer?

Introduction: Folliculin, encoded by FLCN gene, plays a role in the mTORC1 autophagy cascade and its alterations are responsible for the Birt-Hogg-Dubé (BHD) syndrome, characterized by follicle hamartomas, kidney tumors and pneumothorax. Patient and results: We report a 74-years-old woman diagnosed with dementia and carrying a FLCN alteration in absence of any sign of BHD. She also carried an alteration of MAT1A gene, which is also implicated in the regulation of mTORC1. Discussion: The MAT1A variant could have prevented the development of a FLCN-related oncological phenotype. Conversely, our patient presented with dementia that, to date, has yet to be documented in BHD. Folliculin belongs to the DENN family proteins, which includes C9orf72 whose alteration has been associated to neurodegeneration. The folliculin perturbation could affect the C9orf72 activity and our patient could represent the first human model of a relationship between FLCN and C9orf72 across the path of autophagy.

The melatonin receptor 1A (MTNR1A) gene is associated with recurrent and idiopathic calcium nephrolithiasis.

BACKGROUND: Experimental evidence indicate that melatonin regulates some renal tubular functions via specific melatonin receptors (MTNRs) located in the kidney of several avian and mammalian species, including humans. We hypothesized that single nucleotide polymorphisms (SNPs) in the melatonin receptor 1A gene (MTNR1A) might influence the risk of calcium nephrolithiasis. METHODS: We performed a systematic analysis of the MTNR1A gene in 246 recurrent calcium stone formers (136 men, 110 women; mean age 40.2 ± 12.0 years; body mass index 25.8 ± 4.5 kg/m2) and 269 healthy controls comparable for age and gender without a history of nephrolithiasis. RESULTS: Two SNPs in Intron 1 of MTNR1A were significantly associated with calcium nephrolithiasis: rs13140012 (P = 0.0004) and rs6553010 (P = 0.009). The haplotypes resulting from the two SNPs were also differently distributed between stone formers and controls, the haplotype A-T being more represented among stone formers (P = 0.00001) and the haplotype T-C being more common in healthy controls (P = 0.00001). Preliminary functional studies showed that the SNP rs13140012 could modify the binding sites for transcription factors. CONCLUSION: The results of this case-control study indicate a strong association between allelic variants of MTNR1A and recurrent calcium nephrolithiasis.

Common variants in the regulative regions of GRIA1 and GRIA3 receptor genes are associated with migraine susceptibility.

BACKGROUND: Glutamate is the principal excitatory neurotransmitter in the central nervous system which acts by the activation of either ionotropic (AMPA, NMDA and kainate receptors) or G-protein coupled metabotropic receptors. Glutamate is widely accepted to play a major role in the path physiology of migraine as implicated by data from animal and human studies. Genes involved in synthesis, metabolism and regulation of both glutamate and its receptors could be, therefore, considered as potential candidates for causing/predisposing to migraine when mutated. METHODS: The association of polymorphic variants of GRIA1-GRIA4 genes which encode for the four subunits (GluR1-GluR4) of the alpha-amino-3- hydroxy-5-methyl-4-isoxazole-propionic acid (AMPA) receptor for glutamate was tested in migraineurs with and without aura (MA and MO) and healthy controls. RESULTS: Two variants in the regulative regions of GRIA1 (rs2195450) and GRIA3 (rs3761555) genes resulted strongly associated with MA (P = 0.00002 and P = 0.0001, respectively), but not associated with MO, suggesting their role in cortical spreading depression. Whereas the rs548294 variant in GRIA1 gene showed association primarily with MO phenotype, supporting the hypothesis that MA and MO phenotypes could be genetically related. These variants modify binding sites for transcription factors altering the expression of GRIA1 and GRIA3 genes in different conditions. CONCLUSIONS: This study represents the first genetic evidence of a link between glutamate receptors and migraine.

Dualistic Role of BARD1 in Cancer.

BRCA1 Associated RING Domain 1 (BARD1) encodes a protein which interacts with the N-terminal region of BRCA1 in vivo and in vitro. The full length (FL) BARD1 mRNA includes 11 exons and encodes a protein comprising of six domains (N-terminal RING-finger domain, three Ankyrin repeats and two C-terminal BRCT domains) with different functions. Emerging data suggest that BARD1 can have both tumor-suppressor gene and oncogene functions in tumor initiation and progression. Indeed, whereas FL BARD1 protein acts as tumor-suppressor with and without BRCA1 interactions, aberrant splice variants of BARD1 have been detected in various cancers and have been shown to play an oncogenic role. Further evidence for a dualistic role came with the identification of BARD1 as a neuroblastoma predisposition gene in our genome wide association study which has demonstrated that single nucleotide polymorphisms in BARD1 can correlate with risk or can protect against cancer based on their association with the expression of FL and splice variants of BARD1. This review is an overview of how BARD1 functions in tumorigenesis with opposite effects in various types of cancer.

Exome and deep sequencing of clinically aggressive neuroblastoma reveal somatic mutations that affect key pathways involved in cancer progression.

The spectrum of somatic mutation of the most aggressive forms of neuroblastoma is not completely determined. We sought to identify potential cancer drivers in clinically aggressive neuroblastoma.Whole exome sequencing was conducted on 17 germline and tumor DNA samples from high-risk patients with adverse events within 36 months from diagnosis (HR-Event3) to identify somatic mutations and deep targeted sequencing of 134 genes selected from the initial screening in additional 48 germline and tumor pairs (62.5% HR-Event3 and high-risk patients), 17 HR-Event3 tumors and 17 human-derived neuroblastoma cell lines.We revealed 22 significantly mutated genes, many of which implicated in cancer progression. Fifteen genes (68.2%) were highly expressed in neuroblastoma supporting their involvement in the disease. CHD9, a cancer driver gene, was the most significantly altered (4.0% of cases) after ALK.Other genes (PTK2, NAV3, NAV1, FZD1 and ATRX), expressed in neuroblastoma and involved in cell invasion and migration were mutated at frequency ranged from 4% to 2%.Focal adhesion and regulation of actin cytoskeleton pathways, were frequently disrupted (14.1% of cases) thus suggesting potential novel therapeutic strategies to prevent disease progression.Notably BARD1, CHEK2 and AXIN2 were enriched in rare, potentially pathogenic, germline variants.In summary, whole exome and deep targeted sequencing identified novel cancer genes of clinically aggressive neuroblastoma. Our analyses show pathway-level implications of infrequently mutated genes in leading neuroblastoma progression.

Digenic mutational inheritance of the integrin alpha 7 and the myosin heavy chain 7B genes causes congenital myopathy with left ventricular non-compact cardiomyopathy.

BACKGROUND: We report an Italian family in which the proband showed a severe phenotype characterized by the association of congenital fiber type disproportion (CFTD) with a left ventricular non-compaction cardiomyopathy (LVNC). This study was focused on the identification of the responsible gene/s. METHODS AND RESULTS: Using the whole-exome sequencing approach, we identified the proband homozygous missense mutations in two genes, the myosin heavy chain 7B (MYH7B) and the integrin alpha 7 (ITGA7). Both genes are expressed in heart and muscle tissues, and both mutations were predicted to be deleterious and were not found in the healthy population.The R890C mutation in the MYH7B gene segregated with the LVNC phenotype in the examined family. It was also found in one unrelated patient affected by LVNC, confirming a causative role in cardiomyopathy.The E882K mutation in the ITGA7 gene, a key component of the basal lamina of muscle fibers, was found only in the proband, suggesting a role in CFTD. CONCLUSIONS: This study identifies two novel disease genes. Mutation in MYH7B causes a classical LVNC phenotype, whereas mutation in ITGA7 causes CFTD. Both phenotypes represent alterations of skeletal and cardiac muscle maturation and are usually not severe. The severe phenotype of the proband is most likely due to a synergic effect of these two mutations.This study provides new insights into the genetics underlying Mendelian traits and demonstrates a role for digenic inheritance in complex phenotypes.

Giant cell tumor occurring in familial Paget's disease of bone: report of clinical characteristics and linkage analysis of a large pedigree.

Neoplastic degeneration represents a rare but serious complication of Paget's disease of bone (PDB). Although osteosarcomas have been described in up to 1% of PDB cases, giant cell tumors are less frequent and mainly occur in patients with polyostotic disease. We recently characterized a large pedigree with 14 affected members of whom four developed giant cell tumors at pagetic sites. The high number of affected subjects across multiple generations allowed us to better characterize the clinical phenotype and look for possible susceptibility loci. Of interest, all the affected members had polyostotic PDB, but subjects developing giant cell tumors showed an increased disease severity with a reduced clinical response to bisphosphonate treatment and an increased prevalence of bone pain, deformities, and fractures. Together with an increased occurrence of common pagetic complications, affected patients of this pedigree also evidenced a fivefold higher prevalence of coronary artery disease with respect to either the unaffected family members or a comparative cohort of 150 unrelated PDB cases from the same geographical area. This association was further enhanced in the four cases with PDB and giant cell tumors, all of them developing coronary artery disease before 60 years of age. Despite the early onset and the severe phenotype, PDB patients from this pedigree were negative for the presence of SQSTM1 or TNFRSF11A mutations, previously associated with enhanced disease severity. Genome-wide linkage analysis identified six possible candidate regions on chromosomes 1, 5, 6, 8, 10, and 20. Because the chromosome 8 and 10 loci were next to the TNFRSF11B and OPTN genes, we extended the genetic screening to these two genes, but we failed to identify any causative mutation at both the genomic and transcription level, suggesting that a different genetic defect is associated with PDB and potentially giant cell tumor of bone in this pedigree.

Distinct disease phenotypes linked to different combinations of GAA mutations in a large late-onset GSDII sibship.

BACKGROUND: Glycogenosis type II (GSDII or Pompe disease) is an autosomal recessive disease, often characterized by a progressive accumulation of glycogen within lysosomes caused by a deficiency of α-1,4-glucosidase (GAA; acid maltase), a key enzyme of the glycogen degradation pathway. To date, more than 326 different mutations in the GAA gene have been identified in patients with GSDII but the course of the disease is difficult to be predicted on the basis of molecular genetic changes. Studies on large informative families are advisable to better define how genetics and non genetics factors like exercise and diet may influence the clinical phenotype. METHODS AND RESULTS: In this study, we report on clinical, instrumental, and pathological features as well as on molecular analysis of a family with 10 out of 13 siblings affected by late-onset Pompe disease. Three mutations segregated in the family, two of which are novel mutations. Siblings showing a more severe phenotype were compound heterozygous for c.118C > T [p.R40X] and c.2647-7G > A [p.N882fs] on GAA, whereas, two patients showing a mild phenotype were compound heterozygous c.2647-7G > A [p.N882fs] and c.2276G > C [p.G759A] mutations. Quantitative expression analysis showed, in the patients carrying p.R40X/ p.N882fs, a significant (p 0.01) correlation between the levels of expression of the mutated allele and the age at onset of the disease. CONCLUSIONS: As far as we know, this is the largest informative family with late-onset Pompe disease described in the literature showing a peculiar complex set of mutations of GAA gene that may partially elucidate the clinical heterogeneity of this family.

A nonsynonymous TNFRSF11A variation increases NFκB activity and the severity of Paget's disease.

Mutations in the SQSTM1 gene were identified as a common cause of Paget's disease of bone (PDB) but experimental evidence demonstrated that SQSTM1 mutation is not sufficient to induce PDB in vivo. Here, we identified two nonsynonymous single nucleotide polymorphisms (SNPs) (C421T, H141Y and T575C, V192A) in the TNFRSF11A gene, associated with PDB and with the severity of phenotype in a large population of 654 unrelated patients that were previously screened for SQSTM1 gene mutations. The largest effect was found for the T575C variant, yielding an odds ratio of 1.29 (p = 0.003), with the C allele as the risk allele. Moreover, an even more significant p-value (p = 0.0002) was observed in the subgroup of patients with SQSTM1 mutation, with an odds ratio of 1.71. Interestingly, patients with the C allele also showed an increased prevalence of polyostotic disease (68%, 53%, and 51% in patients with CC, CT, and TT genotypes, respectively; p = 0.01), as well as an increased number of affected skeletal sites (2.9, 2.5, and 2.0 in patients with CC, CT, and TT genotypes, respectively, p = 0.008). These differences increased when analyses were restricted to cases with SQSTM1 mutation. In human cell lines, cotrasfection with mutated SQSTM1 and TNFRSF11A(A192) produced a level of activation of NFκB signaling greater than cotrasfection with wild-type SQSTM1 and TNFRSF11A(V192), confirming genetics and clinical evidences. These results provide the first evidence that genetic variation within the OPG/RANK/RANKL system influences the severity of PBD in synergistic action with SQSTM1 gene mutations.

piR_015520 belongs to Piwi-associated RNAs regulates expression of the human melatonin receptor 1A gene.

Piwi-associated RNAs (piRNAs) are a distinct class of 24- to 30-nucleotide-long RNAs produced by a Dicer-independent mechanism, and are associated with Piwi-class Argonaute proteins. In contrast to the several hundred species of microRNAs (miRNAs) identified thus far, piRNAs consist of more than 30,000 different species in humans. Studies in flies, fish and mice implicate these piRNAs in regulating germ line development, the silencing of selfish DNA elements, and maintaining germ line DNA integrity. Most piRNAs map to unique sites in the human genome, including intergenic, intronic, and exonic sequences. However, the role of piRNAs in humans remains to be elucidated. Here, we uncover an unexpected function of the piRNA pathway in humans. We show for the first time, that the piRNA_015520, located in intron 1 of the human Melatonin receptor 1A (MTNR1A) gene, is expressed in adult human tissues (testes and brain) and in the human cell line HEK 293. Although the role of piR_015520 expression in brain tissue remains unknown, the testes-specific expression is consistent with previous findings in several species. Surprisingly, in contrast to the mechanism known for miRNA-mediated modulation of gene expression, piRNA_015520 negatively regulates MTNR1A gene expression by binding to its genomic region. This finding suggests that changes in individual piRNA levels could influence both autoregulatory gene expression and the expression of the gene in which the piRNA is located. These findings offer a new perspective for piRNAs functioning as gene regulators in humans.