Estudios de asociación de genoma completo (GWAS) versus validación funcional: reto de la era post-GWAS / Genome-wide association studies (GWAS) versus functional validation: challenge of the post-GWAS era

En los últimos años se han dedicado muchos esfuerzos a la determinación de variantes y genes que pueden ser impor-tantes en la determinación de la densidad mineral ósea (DMO) y, a su vez, en diversas patologías óseas. Para conseguiresto, la aproximación que ha presentado mayores éxitos ha sido la de los estudios de asociación de genoma completo(GWAS). En particular, en la investigación sobre la biología ósea, se han publicado más de 50 grandes GWAS o metaa-nálisis de GWAS identificando más de 500 loci genéticos asociados con diferentes parámetros óseos como son la DMO,la resistencia ósea y el riesgo de fractura. Si bien el descubrimiento de las variantes asociadas es un aspecto esencial,es igualmente importante la validación funcional de dichas variantes para dilucidar su efecto y la relación causal quetienen con la enfermedad genética. Al tratarse de un aspecto mucho más lento y tedioso, se ha convertido en el nuevoreto de esta era post-GWAS. Entre los genes que ya se han abordado se incluyen varios de la vía de WNT y en especialel gen SOST, que juega un papel muy importante tanto en la determinación de la DMO poblacional como en enferme-dades monogénicas con elevada masa ósea y que ha dado lugar a un nuevo tratamiento contra la osteoporosis. En estarevisión recogemos los principales estudios GWAS con relación a fenotipos del hueso, así como algunos ejemplos devalidaciones funcionales para analizar las asociaciones encontradas en los mismos.(AU)

Case report of a child bearing a novel deleterious splicing variant in PIGT.

RATIONALE: Trio family-based whole exome sequencing (WES) is a powerful tool in the diagnosis of rare neurodevelopmental diseases, even in patients with the unclear diagnosis. There have been previous reports of variants in the phosphatidylinositol glycan anchor biosynthesis class T (PIGT) gene associated with multiple congenital anomalies, with a total of 14 affected individuals across 8 families. PATIENT CONCERNS: An 18-month-old boy of Greek ancestry presented with global developmental delay, generalized tonic-clonic seizures, hypotonia, renal cysts, esotropia, bilateral undescended testes, bilateral vesicoureteric reflux, marked cardiac dextroposition, bilateral talipes equinovarus, and dysmorphic features. DIAGNOSIS: WES revealed 2 compound heterozygous variants in the PIGT gene, c.[494-2A>G]; [547A>C]/p.[Asp122Glyfs*35]; [Thr183Pro]. The splicing mutation was demonstrated to lead to the skipping of exon 4. INTERVENTIONS: Seizures, infections, and other main symptoms were treated. OUTCOMES: The patient died at 2 years of age before the molecular diagnosis was achieved. Genetic counseling has been offered to the family. LESSONS: Most of the clinical features of the patient are in agreement with the previously described PIGT cases corroborating the usefulness of WES as a diagnostic tool.

Pro-osteoporotic miR-320a impairs osteoblast function and induces oxidative stress.

MicroRNAs (miRNAs) are important regulators of many cellular processes, including the differentiation and activity of osteoblasts, and therefore, of bone turnover. MiR-320a is overexpressed in osteoporotic bone tissue but its role in osteoblast function is unknown. In the present study, functional assays were performed with the aim to elucidate the mechanism of miR-320a action in osteoblastic cells. MiR-320a was either overexpressed or inhibited in human primary osteoblasts (hOB) and gene expression changes were evaluated through microarray analysis. In addition, the effect of miR-320a on cell proliferation, viability, and oxidative stress in hOB was evaluated. Finally, matrix mineralization and alkaline phosphatase activity were assessed in order to evaluate osteoblast functionality. Microarray results showed miR-320a regulation of a number of key osteoblast genes and of genes involved in oxidative stress. Regulation of osteoblast differentiation and ossification appeared as the best significant biological processes (PANTHER P value = 3.74E-05; and P value = 3.06E-04, respectively). The other enriched pathway was that of the cellular response to cadmium and zinc ions, mostly by the overexpression of metallothioneins. In hOBs, overexpression of miR-320a increased cell proliferation and oxidative stress levels whereas mineralization capacity was reduced. In conclusion, overexpression of miR-320a increased stress oxidation levels and was associated with reduced osteoblast differentiation and functionality, which could trigger an osteoporotic phenotype.

Functional Characterization of a GGPPS Variant Identified in Atypical Femoral Fracture Patients and Delineation of the Role of GGPPS in Bone-Relevant Cell Types.

Atypical femoral fractures (AFFs) are a rare but potentially devastating event, often but not always linked to bisphosphonate (BP) therapy. The pathogenic mechanisms underlying AFFs remain obscure, and there are no tests available that might assist in identifying those at high risk of AFF. We previously used exome sequencing to explore the genetic background of three sisters with AFFs and three additional unrelated AFF cases, all previously treated with BPs. We detected 37 rare mutations (in 34 genes) shared by the three sisters. Notably, we found a p.Asp188Tyr mutation in the enzyme geranylgeranyl pyrophosphate synthase, a component of the mevalonate pathway, which is critical to osteoclast function and is inhibited by N-BPs. In addition, the CYP1A1 gene, responsible for the hydroxylation of 17ß-estradiol, estrone, and vitamin D, was also mutated in all three sisters and one unrelated patient. Here we present a detailed list of the variants found and report functional analyses of the GGPS1 p.Asp188Tyr mutation, which showed a severe reduction in enzyme activity together with oligomerization defects. Unlike BP treatment, this genetic mutation will affect all cells in the carriers. RNAi knockdown of GGPS1 in osteoblasts produced a strong mineralization reduction and a reduced expression of osteocalcin, osterix, and RANKL, whereas in osteoclasts, it led to a lower resorption activity. Taken together, the impact of the mutated GGPPS and the relevance of the downstream effects in bone cells make it a strong candidate for AFF susceptibility. We speculate that other genes such as CYP1A1 might be involved in AFF pathogenesis, which remains to be functionally proved. The identification of the genetic background for AFFs provides new insights for future development of novel risk assessment tools. © 2018 American Society for Bone and Mineral Research.

The ASXL1 mutation p.Gly646Trpfs*12 found in a Turkish boy with Bohring-Opitz Syndrome.

In line with a recent study showing that ASXL1 mutations found in the common population cannot be ruled out as pathogenic, we have identified the ASXL1 p.Gly646Trpfs*12 mutation-present in 132 individuals in ExAC-as a very probable cause of the disease in a Bohring-Opitz syndrome patient.

Common and rare variants of WNT16, DKK1 and SOST and their relationship with bone mineral density.

Numerous GWAS and candidate gene studies have highlighted the role of the Wnt pathway in bone biology. Our objective has been to study in detail the allelic architecture of three Wnt pathway genes: WNT16, DKK1 and SOST, in the context of osteoporosis. We have resequenced the coding and some regulatory regions of these three genes in two groups with extreme bone mineral density (BMD) (n = ∼50, each) from the BARCOS cohort. No interesting novel variants were identified. Thirteen predicted functional variants have been genotyped in the full cohort (n = 1490), and for ten of them (with MAF > 0.01), the association with BMD has been studied. We have found six variants nominally associated with BMD, of which 2 WNT16 variants predicted to be eQTLs for FAM3C (rs55710688, in the Kozak sequence and rs142005327, within a putative enhancer) withstood multiple-testing correction. In addition, two rare variants in functional regions (rs190011371 in WNT16b 3'UTR and rs570754792 in the SOST TATA box) were found only present in three women each, all with BMD below the mean of the cohort. Our results reinforce the higher importance of regulatory versus coding variants in these Wnt pathway genes and open new ways for functional studies of the relevant variants.

A De Novo FOXP1 Truncating Mutation in a Patient Originally Diagnosed as C Syndrome.

De novo FOXP1 mutations have been associated with intellectual disability (ID), motor delay, autistic features and a wide spectrum of speech difficulties. C syndrome (Opitz C trigonocephaly syndrome) is a rare and genetically heterogeneous condition, characterized by trigonocephaly, craniofacial anomalies and ID. Several different chromosome deletions and and point mutations in distinct genes have been associated with the disease in patients originally diagnosed as Opitz C. By whole exome sequencing we identified a de novo splicing mutation in FOXP1 in a patient, initially diagnosed as C syndrome, who suffers from syndromic intellectual disability with trigonocephaly. The mutation (c.1428 + 1 G > A) promotes the skipping of exon 16, a frameshift and a premature STOP codon (p.Ala450GLyfs*13), as assessed by a minigene strategy. The patient reported here shares speech difficulties, intellectual disability and autistic features with other FOXP1 syndrome patients, and thus the diagnosis for this patient should be changed. Finally, since trigonocephaly has not been previously reported in FOXP1 syndrome, it remains to be proved whether it may be associated with the FOXP1 mutation.

Expression profiling of microRNAs in human bone tissue from postmenopausal women.

Bone tissue is composed of several cell types, which express their own microRNAs (miRNAs) that will play a role in cell function. The set of total miRNAs expressed in all cell types configures the specific signature of the bone tissue in one physiological condition. The aim of this study was to explore the miRNA expression profile of bone tissue from postmenopausal women. Tissue was obtained from trabecular bone and was analyzed in fresh conditions (n = 6). Primary osteoblasts were also obtained from trabecular bone (n = 4) and human osteoclasts were obtained from monocyte precursors after in vitro differentiation (n = 5). MicroRNA expression profiling was obtained for each sample by microarray and a global miRNA analysis was performed combining the data acquired in all the microarray experiments. From the 641 miRNAs detected in bone tissue samples, 346 (54%) were present in osteoblasts and/or osteoclasts. The other 46% were not identified in any of the bone cells analyzed. Intersection of osteoblast and osteoclast arrays identified 101 miRNAs shared by both cell types, which accounts for 30-40% of miRNAs detected in these cells. In osteoblasts, 266 miRNAs were detected, of which 243 (91%) were also present in the total bone array, representing 38% of all bone miRNAs. In osteoclasts, 340 miRNAs were detected, of which 196 (58%) were also present in the bone tissue array, representing 31% of all miRNAs detected in total bone. These analyses provide an overview of miRNAs expressed in bone tissue, broadening our knowledge in the microRNA field.

A De Novo Nonsense Mutation in MAGEL2 in a Patient Initially Diagnosed as Opitz-C: Similarities Between Schaaf-Yang and Opitz-C Syndromes.

Opitz trigonocephaly C syndrome (OTCS) is a rare genetic disorder characterized by craniofacial anomalies, variable intellectual and psychomotor disability, and variable cardiac defects with a high mortality rate. Different patterns of inheritance and genetic heterogeneity are known in this syndrome. Whole exome and genome sequencing of a 19-year-old girl (P7), initially diagnosed with OTCS, revealed a de novo nonsense mutation, p.Q638*, in the MAGEL2 gene. MAGEL2 is an imprinted, maternally silenced, gene located at 15q11-13, within the Prader-Willi region. Patient P7 carried the mutation in the paternal chromosome. Recently, mutations in MAGEL2 have been described in Schaaf-Yang syndrome (SHFYNG) and in severe arthrogryposis. Patient P7 bears resemblances with SHFYNG cases but has other findings not described in this syndrome and common in OTCS. We sequenced MAGEL2 in nine additional OTCS patients and no mutations were found. This study provides the first clear molecular genetic basis for an OTCS case, indicates that there is overlap between OTCS and SHFYNG syndromes, and confirms that OTCS is genetically heterogeneous. Genes encoding MAGEL2 partners, either in the retrograde transport or in the ubiquitination-deubiquitination complexes, are promising candidates as OTCS disease-causing genes.

Screening of CD96 and ASXL1 in 11 patients with Opitz C or Bohring-Opitz syndromes.

Opitz C trigonocephaly (or Opitz C syndrome, OTCS) and Bohring-Opitz syndrome (BOS or C-like syndrome) are two rare genetic disorders with phenotypic overlap. The genetic causes of these diseases are not understood. However, two genes have been associated with OTCS or BOS with dominantly inherited de novo mutations. Whereas CD96 has been related to OTCS (one case) and to BOS (one case), ASXL1 has been related to BOS only (several cases). In this study we analyze CD96 and ASXL1 in a group of 11 affected individuals, including 2 sibs, 10 of them were diagnosed with OTCS, and one had a BOS phenotype. Exome sequences were available on six patients with OTCS and three parent pairs. Thus, we could analyze the CD96 and ASXL1 sequences in these patients bioinformatically. Sanger sequencing of all exons of CD96 and ASXL1 was carried out in the remaining patients. Detailed scrutiny of the sequences and assessment of variants allowed us to exclude putative pathogenic and private mutations in all but one of the patients. In this patient (with BOS) we identified a de novo mutation in ASXL1 (c.2100dupT). By nature and location within the gene, this mutation resembles those previously described in other BOS patients and we conclude that it may be responsible for the condition. Our results indicate that in 10 of 11, the disease (OTCS or BOS) cannot be explained by small changes in CD96 or ASXL1. However, the cohort is too small to make generalizations about the genetic etiology of these diseases.

MiRNA profiling of whole trabecular bone: identification of osteoporosis-related changes in MiRNAs in human hip bones.

BACKGROUND: MicroRNAs (miRNAs) are important regulators of gene expression, with documented roles in bone metabolism and osteoporosis, suggesting potential therapeutic targets. Our aim was to identify miRNAs differentially expressed in fractured vs nonfractured bones. Additionally, we performed a miRNA profiling of primary osteoblasts to assess the origin of these differentially expressed miRNAs. METHODS: Total RNA was extracted from (a) fresh femoral neck trabecular bone from women undergoing hip replacement due to either osteoporotic fracture (OP group, n = 6) or osteoarthritis in the absence of osteoporosis (Control group, n = 6), matching the two groups by age and body mass index, and (b) primary osteoblasts obtained from knee replacement due to osteoarthritis (n = 4). Samples were hybridized to a microRNA array containing more than 1900 miRNAs. Principal component analysis (PCA) plots and heat map hierarchical clustering were performed. For comparison of expression levels, the threshold was set at log fold change > 1.5 and a p-value < 0.05 (corrected for multiple testing). RESULTS: Both PCA and heat map analyses showed that the samples clustered according to the presence or absence of fracture. Overall, 790 and 315 different miRNAs were detected in fresh bone samples and in primary osteoblasts, respectively, 293 of which were common to both groups. A subset of 82 miRNAs was differentially expressed (p < 0.05) between osteoporotic and control osteoarthritic samples. The eight miRNAs with the lowest p-values (and for which a validated miRNA qPCR assay was available) were assayed, and two were confirmed: miR-320a and miR-483-5p. Both were over-expressed in the osteoporotic samples and expressed in primary osteoblasts. miR-320a is known to target CTNNB1 and predicted to regulate RUNX2 and LEPR, while miR-483-5p down-regulates IGF2. We observed a reduction trend for this target gene in the osteoporotic bone. CONCLUSIONS: We identified two osteoblast miRNAs over-expressed in osteoporotic fractures, which opens novel prospects for research and therapy.

Alelo doble mutante en el gen EXT1 no informado previamente en una adolescente con exostosis múltiple hereditaria / Double mutant alleles in the EXT1 gene not previously reported in a teenager with hereditary multiple exostoses

Las formas hereditarias de exostosis múltiple, actualmente denominada EXT1 / EXT2-CDG dentro de los desórdenes congénitos de la glicosilación, son los tumores óseos benignos más comunes y se caracterizan por la formación de lesiones óseas cubiertas de cartílago, localizadas en yuxtaposición a epífisis de huesos largos, aunque, en los casos graves, pueden presentar una amplia distribución. El inicio es variable desde los 2-3 años hasta los 13-15 y presenta una incidencia estimada que va de 1/18 000 a 1/50 000 casos en los países europeos. Se presenta el caso de un doble alelo mutante en el gen EXT1 no informado previamente en una adolescente y su familia con exostosis múltiple hereditaria.

Hereditary forms of multiple exostoses, now called EXT1/EXT2-CDG within Congenital Disorders of Glycosylation, are the most common benign bone tumors in humans and clinical description consists of the formation of several cartilage-capped bone tumors, usually benign and localized in the juxta-epiphyseal region of long bones, although wide body dissemination in severe cases is not uncommon. Onset of the Alelo doble mutante en el gen EXT1 no informado previamente en una adolescente con exostosis múltiple hereditariaDouble mutant alleles in the EXT1 gene not previously reported in a teenager with hereditary multiple exostosesdisease is variable ranging from 2-3 years up to 13-15 years with an estimated incidence ranging from 1/18 000 to 1/50 000 cases in European countries. We present a double mutant alleles in the EXT1 gene not previously reported in a teenager and her family with hereditary multiple exostoses

Alelo doble mutante en el gen EXT1 no informado previamente en una adolescente con exostosis múltiple hereditaria / Double mutant alleles in the EXT1 gene not previously reported in a teenager with hereditary multiple exostoses

Las formas hereditarias de exostosis múltiple, actualmente denominada EXT1 / EXT2-CDG dentro de los desórdenes congénitos de la glicosilación, son los tumores óseos benignos más comunes y se caracterizan por la formación de lesiones óseas cubiertas de cartílago, localizadas en yuxtaposición a epífisis de huesos largos, aunque, en los casos graves, pueden presentar una amplia distribución. El inicio es variable desde los 2-3 años hasta los 13-15 y presenta una incidencia estimada que va de 1/18 000 a 1/50 000 casos en los países europeos. Se presenta el caso de un doble alelo mutante en el gen EXT1 no informado previamente en una adolescente y su familia con exostosis múltiple hereditaria.(AU)

[Double mutant alleles in the EXT1 gene not previously reported in a teenager with hereditary multiple exostoses]. / Alelo doble mutante en el gen EXT1 no informado previamente en una adolescente con exostosis múltiple hereditaria.

Hereditary forms of multiple exostoses, now called EXT1/EXT2-CDG within Congenital Disorders of Glycosylation, are the most common benign bone tumors in humans and clinical description consists of the formation of several cartilage-capped bone tumors, usually benign and localized in the juxta-epiphyseal region of long bones, although wide body dissemination in severe cases is not uncommon. Onset of the disease is variable ranging from 2-3 years up to 13-15 years with an estimated incidence ranging from 1/18,000 to 1/50,000 cases in European countries. We present a double mutant alleles in the EXT1 gene not previously reported in a teenager and her family with hereditary multiple exostoses.

Genetic analysis of high bone mass cases from the BARCOS cohort of Spanish postmenopausal women.

The aims of the study were to establish the prevalence of high bone mass (HBM) in a cohort of Spanish postmenopausal women (BARCOS) and to assess the contribution of LRP5 and DKK1 mutations and of common bone mineral density (BMD) variants to a HBM phenotype. Furthermore, we describe the expression of several osteoblast-specific and Wnt-pathway genes in primary osteoblasts from two HBM cases. A 0.6% of individuals (10/1600) displayed Z-scores in the HBM range (sum Z-score >4). While no mutation in the relevant exons of LRP5 was detected, a rare missense change in DKK1 was found (p.Y74F), which cosegregated with the phenotype in a small pedigree. Fifty-five BMD SNPs from Estrada et al. [NatGenet 44:491-501,2012] were genotyped in the HBM cases to obtain risk scores for each individual. In this small group of samples, Z-scores were found inversely related to risk scores, suggestive of a polygenic etiology. There was a single exception, which may be explained by a rare penetrant genetic variant, counterbalancing the additive effect of the risk alleles. The expression analysis in primary osteoblasts from two HBM cases and five controls suggested that IL6R, DLX3, TWIST1 and PPARG are negatively related to Z-score. One HBM case presented with high levels of RUNX2, while the other displayed very low SOX6. In conclusion, we provide evidence of lack of LRP5 mutations and of a putative HBM-causing mutation in DKK1. Additionally, we present SNP genotyping and expression results that suggest additive effects of several genes for HBM.

Description of extreme longevity in the Balearic Islands: Exploring a potential Blue Zone in Menorca, Spain.

AIM: We aimed to determine whether there was a Blue Zone, an area characterized by extreme longevity, in Menorca, Spain. METHODS: We explored official statistics of the Balearic Islands, Spain, and calculated life expectancy from 1991 to 2009, by sex and island, among other demographic estimators. RESULTS: The life expectancy at birth in Menorca reached a peak in 2007 with 82.3 years, but since 2001 in females and 2007 in males, it plateaued and then descended. The percentage of centenarians in Es Migjorn Gran was 0.17% in 2001, only because of the 0.33% male contribution, and it returned to within the average Balearic rates. There is no record of any other surviving centenarian after 2006 in Es Migjorn Gran. Results obtained with other demographic indices are confirmatory. CONCLUSIONS: It appears there was no extreme longevity or Blue Zone in Es Migjorn Gran and the surrounding villages in Menorca.

Genetic determinants of aromatase inhibitor-related arthralgia: the B-ABLE cohort study.

A major side effect of aromatase inhibitor (AI) therapy is AI-related arthralgia (AIA), which often leads to therapy discontinuation. We aimed to identify genetic variants associated with AIA and therapy discontinuation in the first year of AI treatment. Our prospective cohort study included 343 postmenopausal women with early breast cancer starting AI therapy. Single nucleotide polymorphisms (SNPs) in candidate genes involved in estrogen and vitamin D signaling were selected. Univariate and multivariate linear/logistic regressions were fitted in order to asses the association between studied SNPs and AIA intensity (visual analogic scale score) at 3 and 12 months of follow-up, worsening pain, and therapy discontinuation. We also tested for a priori-defined interactions by introducing multiplicative terms in the regression equations. SNPs in CYP17A1 and VDR genes appeared significantly associated with AIA (P = 0.003, P = 0.012, respectively). One SNP in CYP27B1 gene was related to therapy discontinuation [P = 0.02; OR 0.29 (0.09-0.99)]. We revealed interactions between CYP27B1 and both CYP17A1 (P = 0.01) and VDR SNPs (P = 0.06). Furthermore, an additive effect on pain intensity was shown for unfavorable alleles, with two points higher mean absolute pain increase and up to 5.3-fold higher risk of worsening pain compared to favorable genotypes. SNPs in CYP17A1, VDR, and CYP27B1 genes predict the risk of AIA. Their determination would be useful to trigger the monitoring strategies in women at risk of therapy discontinuation.