Pathogenic RAB34 variants impair primary cilium assembly and cause a novel oral-facial-digital syndrome.

Oral-facial-digital syndromes (OFDS) are a group of clinically and genetically heterogeneous disorders characterized by defects in the development of the face and oral cavity along with digit anomalies. Pathogenic variants in over 20 genes encoding ciliary proteins have been found to cause OFDS through deleterious structural or functional impacts on primary cilia. We identified by exome sequencing bi-allelic missense variants in a novel disease-causing ciliary gene RAB34 in four individuals from three unrelated families. Affected individuals presented a novel form of OFDS (OFDS-RAB34) accompanied by cardiac, cerebral, skeletal and anorectal defects. RAB34 encodes a member of the Rab GTPase superfamily and was recently identified as a key mediator of ciliary membrane formation. Unlike many genes required for cilium assembly, RAB34 acts selectively in cell types that use the intracellular ciliogenesis pathway, in which nascent cilia begin to form in the cytoplasm. We find that the protein products of these pathogenic variants, which are clustered near the RAB34 C-terminus, exhibit a strong loss of function. Although some variants retain the ability to be recruited to the mother centriole, cells expressing mutant RAB34 exhibit a significant defect in cilium assembly. While many Rab proteins have been previously linked to ciliogenesis, our studies establish RAB34 as the first small GTPase involved in OFDS and reveal the distinct clinical manifestations caused by impairment of intracellular ciliogenesis.

Evolutionary and functional analyses of LRP5 in archaic and extant modern humans.

BACKGROUND: The human lineage has undergone a postcranial skeleton gracilization (i.e. lower bone mass and strength relative to body size) compared to other primates and archaic populations such as the Neanderthals. This gracilization has been traditionally explained by differences in the mechanical load that our ancestors exercised. However, there is growing evidence that gracilization could also be genetically influenced. RESULTS: We have analyzed the LRP5 gene, which is known to be associated with high bone mineral density conditions, from an evolutionary and functional point of view. Taking advantage of the published genomes of archaic Homo populations, our results suggest that this gene has a complex evolutionary history both between archaic and living humans and within living human populations. In particular, we identified the presence of different selective pressures in archaics and extant modern humans, as well as evidence of positive selection in the African and South East Asian populations from the 1000 Genomes Project. Furthermore, we observed a very limited evidence of archaic introgression in this gene (only at three haplotypes of East Asian ancestry out of the 1000 Genomes), compatible with a general erasing of the fingerprint of archaic introgression due to functional differences in archaics compared to extant modern humans. In agreement with this hypothesis, we observed private mutations in the archaic genomes that we experimentally validated as putatively increasing bone mineral density. In particular, four of five archaic missense mutations affecting the first ß-propeller of LRP5 displayed enhanced Wnt pathway activation, of which two also displayed reduced negative regulation. CONCLUSIONS: In summary, these data suggest a genetic component contributing to the understanding of skeletal differences between extant modern humans and archaic Homo populations.

Rare predicted deleterious FEZF2 variants are associated with a neurodevelopmental phenotype.

FEZF2 encodes a transcription factor critical to neurodevelopment that regulates other neurodevelopment genes. Rare variants in FEZF2 have previously been suggested to play a role in autism, and cases of 3p14 microdeletions that include FEZF2 share a neurodevelopmental phenotype including mild dysmorphic features and intellectual disability. We identified seven heterozygous predicted deleterious variants in FEZF2 (three frameshifts, one recurrent missense in two independent cases, one nonsense, and one complete gene deletion) in unrelated individuals with neurodevelopmental disorders including developmental delay/intellectual disability, autism, and/or attention-deficit/hyperactivity. Variants were confirmed to be de novo in five of seven cases and paternally inherited from an affected father in one. Predicted deleterious variants in FEZF2 may affect the expression of genes that are involved in fate choice pathways in developing neurons, and thus contribute to the neurodevelopmental phenotype. Future studies are needed to clarify the mechanism by which FEZF2 leads to this neurodevelopmental disorder.

Genetics and Genomics of SOST: Functional Analysis of Variants and Genomic Regulation in Osteoblasts.

SOST encodes the sclerostin protein, which acts as a key extracellular inhibitor of the canonical Wnt pathway in bone, playing a crucial role in skeletal development and bone homeostasis. The objective of this work was to assess the functionality of two variants previously identified (the rare variant rs570754792 and the missense variant p.Val10Ile) and to investigate the physical interactors of the SOST proximal promoter region in bone cells. Through a promoter luciferase reporter assay we show that the minor allele of rs570754792, a variant located in the extended TATA box motif, displays a significant decrease in promoter activity. Likewise, through western blot studies of extracellular and intracellular sclerostin, we observe a reduced expression of the p.Val10Ile mutant protein. Finally, using a circular chromosome conformation capture assay (4C-seq) in 3 bone cell types (MSC, hFOB, Saos-2), we have detected physical interactions between the SOST proximal promoter and the ECR5 enhancer, several additional enhancers located between EVT4 and MEOX1 and a distant region containing exon 18 of DHX8. In conclusion, SOST presents functional regulatory and missense variants that affect its expression and displays physical contacts with far reaching genomic sequences, which may play a role in its regulation within bone cells.

Functional Analyses of Four CYP1A1 Missense Mutations Present in Patients with Atypical Femoral Fractures.

Osteoporosis is the most common metabolic bone disorder and nitrogen-containing bisphosphonates (BP) are a first line treatment for it. Yet, atypical femoral fractures (AFF), a rare adverse effect, may appear after prolonged BP administration. Given the low incidence of AFF, an underlying genetic cause that increases the susceptibility to these fractures is suspected. Previous studies uncovered rare CYP1A1 mutations in osteoporosis patients who suffered AFF after long-term BP treatment. CYP1A1 is involved in drug metabolism and steroid catabolism, making it an interesting candidate. However, a functional validation for the AFF-associated CYP1A1 mutations was lacking. Here we tested the enzymatic activity of four such CYP1A1 variants, by transfecting them into Saos-2 cells. We also tested the effect of commonly used BPs on the enzymatic activity of the CYP1A1 forms. We demonstrated that the p.Arg98Trp and p.Arg136His CYP1A1 variants have a significant negative effect on enzymatic activity. Moreover, all the BP treatments decreased CYP1A1 activity, although no specific interaction with CYP1A1 variants was found. Our results provide functional support to the hypothesis that an additive effect between CYP1A1 heterozygous mutations p.Arg98Trp and p.Arg136His, other rare mutations and long-term BP exposure might generate susceptibility to AFF.

Chromatin assembly factor subunit CHAF1A as a monogenic cause for oculo-auriculo-vertebral spectrum.

Oculo-auriculo-vertebral spectrum (OAVS) is characterized by abnormal development of the 1st and 2nd branchial arches. Despite arguments against a monogenic condition, a few genes have been involved in a minority of cases. We now report heterozygous, presumably loss-of function variants in the CHAF1A gene in 8 individuals, including 3 members of the same family. Four cases fulfill stringent diagnostic criteria for OAVS, including asymmetric ear dysplasia, preauricular tags, mandibular asymmetry +/- vertebral malformations. Two patients also presented with kidney malformations. CHAF1A encodes a subunit of CAF-1 (chromatin assembly factor-1), a heterotrimeric protein complex responsible for the deposition of newly synthesized histones H3-H4 onto the newly synthetized DNA strand during replication. The identification of loss-of-unction variants in CHAF1A is consistent with the hypothesis of early developmental genes dysregulation driving OAVS and other associations recently lumped under the acronym Recurrent Constellations of Embryonic Malformations (RCEM).

The human ciliopathy protein RSG1 links the CPLANE complex to transition zone architecture.

Cilia are essential organelles and variants in genes governing ciliary function result in ciliopathic diseases. The Ciliogenesis and PLANar polarity Effectors (CPLANE) protein complex is essential for ciliogenesis in animals models but remains poorly defined. Notably, all but one subunit of the CPLANE complex have been implicated in human ciliopathy. Here, we identify three families in which variants in the remaining CPLANE subunit CPLANE2/RSG1 also cause ciliopathy. These patients display cleft palate, tongue lobulations and polydactyly, phenotypes characteristic of Oral-Facial-Digital Syndrome. We further show that these alleles disrupt two vital steps of ciliogenesis, basal body docking and recruitment of intraflagellar transport proteins. Moreover, APMS reveals that Rsg1 binds the CPLANE and also the transition zone protein Fam92 in a GTP-dependent manner. Finally, we show that CPLANE is generally required for normal transition zone architecture. Our work demonstrates that CPLANE2/RSG1 is a causative gene for human ciliopathy and also sheds new light on the mechanisms of ciliary transition zone assembly.

Wnt Pathway Extracellular Components and Their Essential Roles in Bone Homeostasis.

The Wnt pathway is involved in several processes essential for bone development and homeostasis. For proper functioning, the Wnt pathway is tightly regulated by numerous extracellular elements that act by both activating and inhibiting the pathway at different moments. This review aims to describe, summarize and update the findings regarding the extracellular modulators of the Wnt pathway, including co-receptors, ligands and inhibitors, in relation to bone homeostasis, with an emphasis on the animal models generated, the diseases associated with each gene and the bone processes in which each member is involved. The precise knowledge of all these elements will help us to identify possible targets that can be used as a therapeutic target for the treatment of bone diseases such as osteoporosis.

Clinical description and genetic analysis of a novel familial skeletal dysplasia characterized by high bone mass and lucent bone lesions.

High bone mass (HBM) disorders are a clinically and genetically heterogeneous subgroup of rare skeletal dysplasias. Here we present a case of a previously unreported familial skeletal dysplasia characterized by HBM and lucent bone lesions that we aimed to clinically characterize and genetically investigate. For phenotyping, we reviewed past clinical records and imaging tests, and performed physical examination (PE), bone densitometry, and mineral panels in affected individuals, including a male proband, his son and daughter, in addition to unaffected controls, including the proband's wife and brother. Affected individuals also underwent impact microindentation (IMI). In an effort to elucidate the disorder's molecular etiology, whole exome sequencing (WES) was performed in all individuals to filter for rare variants present only in affected ones. The cases displayed a unique skeletal phenotype with a mix of sclerotic features and lucent bone lesions, and high IMI values. Bone mineral density was very elevated in the proband and his daughter. The proband's daughter also exhibited idiopathic scoliosis (IS), in addition to mild thrombocytopenia and mild structural thyroid abnormalities, which were the only extra-skeletal abnormalities identified. WES analysis yielded 5 rare putative pathogenic variants in affected members in genes that are associated with bone metabolism including: SEM4AD, TBX18, PTCH1, PTK7, and ADGRE5. The PTK7 variant appeared as possibly implicated in the development of IS while the TBX18 and SEMA4D variants stood out as the strongest candidates for the lucent bone lesions and HBM, respectively, given their high predicted pathogenicity and putative role in bone biology. Variant functionality should be addressed in the future to assess their implication in skeletal metabolism as it is the first time that mutations in TBX18 and SEMA4D have been associated to bone developmental lesions and mineral metabolism in a clinical setting.

Genetic Analysis in a Familial Case With High Bone Mineral Density Suggests Additive Effects at Two Loci.

Osteoporosis is the most common bone disease, characterized by a low bone mineral density (BMD) and increased risk of fracture. At the other end of the BMD spectrum, some individuals present strong, fracture-resistant, bones. Both osteoporosis and high BMD are heritable and their genetic architecture encompasses polygenic inheritance of common variants and some cases of monogenic highly penetrant variants in causal genes. We have investigated the genetics of high BMD in a family segregating this trait in an apparently Mendelian dominant pattern. We searched for rare causal variants by whole-exome sequencing in three affected and three nonaffected family members. Using this approach, we have identified 38 rare coding variants present in the proband and absent in the three individuals with normal BMD. Although we have found four variants shared by the three affected members of the family, we have not been able to relate any of these to the high-BMD phenotype. In contrast, we have identified missense variants in two genes, VAV3 and ADGRE5, each shared by two of out of three affected members, whose loss of function fits with the phenotype of the family. In particular, the proband, a woman displaying the highest BMD (sum Z-score = 7), carries both variants, whereas the other two affected members carry one each. VAV3 encodes a guanine-nucleotide-exchange factor with an important role in osteoclast activation and function. Although no previous cases of VAV3 mutations have been reported in humans, Vav3 knockout (KO) mice display dense bones, similarly to the high-BMD phenotype present in our family. The ADGRE5 gene encodes an adhesion G protein-coupled receptor expressed in osteoclasts whose KO mouse displays increased trabecular bone volume. Combined, these mouse and human data highlight VAV3 and ADGRE5 as novel putative high-BMD genes with additive effects, and potential therapeutic targets for osteoporosis. © 2022 The Authors. JBMR Plus published by Wiley Periodicals LLC on behalf of American Society for Bone and Mineral Research.

On the association between Chiari malformation type 1, bone mineral density and bone related genes.

Background: Chiari malformation type 1 (C1M) is a neurological disease characterized by herniation of the cerebellar tonsils below the foramen magnum. Cranial bone constriction is suspected to be its main cause. To date, genes related to bone development (e.g. DKK1 or COL1A2) have been associated with C1M, while some bone diseases (e.g. Paget) have been found to cosegregate with C1M. Nevertheless, the association between bone mineral density (BMD) and C1M has not been investigated, yet. Here, we systematically investigate the association between C1M and BMD, and between bone related genes and C1M. Methods: We have recruited a small cohort of C1M patients (12 unrelated patients) in whom we have performed targeted sequencing of an in-house bone-related gene panel and BMD determination through non-invasive DXA. Results: In the search for association between the bone related genes and C1M we have found variants in more than one C1M patient in WNT16, CRTAP, MYO7A and NOTCH2. These genes have been either associated with craniofacial development in different ways, or previously associated with C1M (MYO7A). Regarding the potential link between BMD and C1M, we have found three osteoporotic patients and one patient who had high BMD, very close to the HBM phenotype values, although most patients had normal BMD. Conclusions: Variants in bone related genes have been repeatedly found in some C1M cases. The relationship of bone genes with C1M deserves further study, to get a clearer estimate of their contribution to its etiology. No direct correlation between BMD and C1M was observed.

Gene Network of Susceptibility to Atypical Femoral Fractures Related to Bisphosphonate Treatment.

Atypical femoral fractures (AFF) are rare fragility fractures in the subtrocantheric or diaphysis femoral region associated with long-term bisphosphonate (BP) treatment. The etiology of AFF is still unclear even though a genetic basis is suggested. We performed whole exome sequencing (WES) analysis of 12 patients receiving BPs for at least 5 years who sustained AFFs and 4 controls, also long-term treated with BPs but without any fracture. After filtration and prioritization of rare variants predicted to be damaging and present in genes shared among at least two patients, a total of 272 variants in 132 genes were identified. Twelve of these genes were known to be involved in bone metabolism and/or AFF, highlighting DAAM2 and LRP5, both involved in the Wnt pathway, as the most representative. Afterwards, we intersected all mutated genes with a list of 34 genes obtained from a previous study of three sisters with BP-related AFF, identifying nine genes. One of these (MEX3D) harbored damaging variants in two AFF patients from the present study and one shared among the three sisters. Gene interaction analysis using the AFFNET web suggested a complex network among bone-related genes as well as with other mutated genes. BinGO biological function analysis highlighted cytoskeleton and cilium organization. In conclusion, several genes and their interactions could provide genetic susceptibility to AFF, that along with BPs treatment and in some cases with glucocorticoids may trigger this so feared complication.

Functional Assessment of Coding and Regulatory Variants From the DKK1 Locus.

The DKK1 gene encodes an extracellular inhibitor of the Wnt pathway with an important role in bone tissue development, bone homeostasis, and different critical aspects of bone biology. Several BMD genome-wide association studies (GWASs) have consistently found association with SNPs in the DKK1 genomic region. For these reasons, it is important to assess the functionality of coding and regulatory variants in the gene. Here, we have studied the functionality of putative regulatory variants, previously found associated with BMD in different studies by others and ourselves, and also six missense variants present in the general population. Using a Wnt-pathway-specific luciferase reporter assay, we have determined that the variants p.Ala41Thr, p.Tyr74Phe, p.Arg120Leu, and p.Ser157Ile display a reduced DKK1 inhibitory capacity as compared with WT. This result agrees with the high-bone-mass (HBM) phenotype of two women from our cohort who carried mutations p.Tyr74Phe or p.Arg120Leu. On the other hand, by means of a circularized chromosome conformation capture- (4C-) sequencing experiment, we have detected that the region containing 24 BMD-GWA variants, located 350-kb downstream of DKK1, interacts both with DKK1 and the LNCAROD (LncRNA-activating regulator of DKK1, AKA LINC0148) in osteoblastic cells. In conclusion, we have shown that some rare coding variants are partial loss-of-function mutations that may lead to a HBM phenotype, whereas the common SNPs associated with BMD in GWASs belong to a putative long-range regulatory region, through a yet unknown mechanism involving LNCAROD. © 2020 The Authors. JBMR Plus published by Wiley Periodicals LLC on behalf of American Society for Bone and Mineral Research.

Functional characterization of the C7ORF76 genomic region, a prominent GWAS signal for osteoporosis in 7q21.3.

Genome-wide association studies (GWAS) have repeatedly identified genetic variants associated with bone mineral density (BMD) and osteoporotic fracture in non-coding regions of C7ORF76, a poorly studied gene of unknown function. The aim of the present study was to elucidate the causality and molecular mechanisms underlying the association. We re-sequenced the genomic region in two extreme BMD groups from the BARCOS cohort of postmenopausal women to search for functionally relevant variants. Eight selected variants were tested for association in the complete cohort and 2 of them (rs4342521 and rs10085588) were found significantly associated with lumbar spine BMD and nominally associated with osteoporotic fracture. cis-eQTL analyses of these 2 SNPs, together with SNP rs4727338 (GWAS lead SNP in Estrada et al., Nat Genet. 44:491-501, 2012), performed in human primary osteoblasts, disclosed a statistically significant influence on the expression of the proximal neighbouring gene SLC25A13 and a tendency on the distal SHFM1. We then studied the functionality of a putative upstream regulatory element (UPE), containing rs10085588. Luciferase reporter assays showed transactivation capability with a strong allele-dependent effect. Finally, 4C-seq experiments in osteoblastic cell lines showed that the UPE interacted with different tissue-specific enhancers and a lncRNA (LOC100506136) in the region. In summary, this study is the first one to analyse in depth the functionality of C7ORF76 genomic region. We provide functional regulatory evidence for the rs10085588, which may be a causal SNP within the 7q21.3 GWAS signal for osteoporosis.

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)

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.