Joint study of two genome-wide association meta-analyses identified 20p12.1 and 20q13.33 for bone mineral density

The GWAS candidate far upstream element binding protein 3 (FUBP3) is required for normal skeletal growth, and adult bone mass and strength in mice

Bone mineral density (BMD) and height are highly heritable traits for which hundreds of genetic loci have been linked through genome wide association studies (GWAS). FUBP3 is a DNA and RNA binding protein best characterised as a transcriptional regulator of c-Myc, but little is known about its role in vivo. Single nucleotide polymorphisms in FUBP3 at the 9q34.11 locus have been associated with BMD, fracture and height in multiple GWAS, but FUBP3 has no previously established role in the skeleton. We analysed Fubp3-deficient mice to determine the consequence of FUBP3 deficiency in vivo. Mice lacking Fubp3 had reduced survival to adulthood and impaired skeletal growth. Bone mass was decreased, most strikingly in the vertebrae, with altered trabecular micro-architecture. Fubp3 deficient bones were also weak. These data provide the first functional demonstration that Fubp3 is required for normal skeletal growth and development and maintenance of adult bone structure and strength, indicating that FUBP3 contributes to the GWAS association of 9q34.11 with variation in height, BMD and fracture.

Regulation of WNT16 in bone may involve upstream enhancers within CPED1

WNT16 stands up as an essential gene for bone homeostasis. Here, we present new evidence of the functional role of a particular region within WNT16. Performing 4 C chromatin conformation analysis in three osteoblast-related cells (the human fetal osteoblast hFOB 1.19 cell line, Saos 2 osteosarcoma cell line and mesenchymal Stem Cells -MSC-), we identify physical interactions between the proximal part of WNT16 intron 2, shown here to be an active promoter in Saos 2 osteosarcoma cells, and several putative regulatory regions within CPED1. Analysis of previously published RNA-seq data from hFOB cells disclosed low expression of a region located downstream of this promoter. Our results suggest a novel regulatory mechanism of WNT16 in bone, mediated by physical interaction with various enhancer regions within CPED1.

The genetic architecture of age at menarche and its causal effects on other traits

Age at menarche (AAM) is a sign of puberty of females. It is a heritable trait associated with various adult diseases. However, the genetic mechanism that determines AAM and links it to disease risk is poorly understood. Aiming to uncover the genetic basis for AAM, we conducted a joint association study in up to 438,089 women from 3 genome-wide association studies of European and East Asian ancestries. A series of bioinformatical analyses and causal inference were then followed to explore in-depth annotations at the associated loci and infer the causal relationship between AAM and other complex traits/diseases. This largest meta-analysis identified a total of 21 novel AAM associated loci at the genome wide significance level (P < 5.0  × 10-8), 4 of which were European ancestry-specific loci. Functional annotations prioritized 33 candidate genes at newly identified loci. Significant genetic correlations were observed between AAM and 67 complex traits. Further causal inference demonstrated the effects of AAM on 13 traits, including forced vital capacity (FVC), high blood pressure, age at first live birth, etc, indicating that earlier AAM causes lower FVC, worse lung function, hypertension and earlier age at first (last) live birth. Enrichment analysis identified 5 enriched tissues, including the hypothalamus middle, hypothalamo hypophyseal system, neurosecretory systems, hypothalamus and retina. Our findings may provide useful insights that elucidate the mechanisms determining AAM and the genetic interplay between AAM and some traits of women.

The osteoporosis susceptibility SNP rs188303909 at 2q14.2 regulates EN1 expression by modulating DNA methylation and E2F6 binding

EN1 encodes a homeodomain-containing transcription factor and is a determinant of bone density and fracture. Previous powerful genome-wide association studies (GWASs) have identified multiple single-nucleotide polymorphisms (SNPs) near EN1 at 2q14.2 locus for osteoporosis, but the causal SNPs and functional mechanisms underlying these associations are poorly understood. The target genes regulated by the transcription factor EN1 are also unclear. In this study, we identified rs188303909, a functional CpG-SNP, as a causal SNP for osteoporosis at 2q14.2 through the integration of functional and epigenomic analyses. Functional experiments demonstrated that unmethylated rs188303909 acted as a strong allele-specific distal enhancer to regulate EN1 expression by modifying the binding of transcription factor E2F6, but rs188303909 methylation attenuated the active effect of E2F6 on EN1 expression. Importantly, transcription factor EN1 could differentially bind osteoporosis GWAS lead SNPs rs4869739-T and rs4355801-G to upregulate CCDC170 and COLEC10 expression, thus promoting bone formation. Our study provided a mechanistic insight into expression regulation of the osteoporosis susceptibility gene EN1, which could be a potential therapeutic target for osteoporosis precision medicine. KEY MESSAGES: CpG-SNP rs188303909 is a causal SNP at the osteoporosis susceptibility locus 2q14.2. Rs188303909 distally regulates EN1 expression by modulating DNA methylation and E2F6 binding. EN1 upregulates CCDC170 and COLEC10 expression through osteoporosis GWAS lead SNPs rs4869739 and rs4355801.

A common variant SNP rs1937810 in the MPP7 gene contributes to the susceptibility of breast cancer in the Chinese Han population

BACKGROUND

Breast cancer (BC) is common cancer caused by environmental factors and genetic ones. Previous evidence has linked gene MAGUK P55 Scaffold Protein 7 (MPP7) to BC, despite that there has been no research evaluating the relationship between MPP7 genetic polymorphisms and BC susceptibility. We aimed to investigate the potential association of the MPP7 gene with the susceptibility to BC in Han Chinese individuals.

METHODS

In total, 1390 patients with BC and 2480 controls were enrolled. For genotyping, 20 tag SNPs were chosen. The serum levels of protein MPP7 were measured in all subjects using an enzyme-linked immunosorbent assay. Genetic association analysis was performed in both genotypic and allelic modes, and the relationship between BC patients' clinical features and genotypes of relevant SNPs was examined. The functional implications of significant markers were also evaluated.

RESULTS

After adjusting for Bonferroni correction, SNP rs1937810 was found to be significantly associated with the risk of BC (p = 1.19 × 10-4 ). The odds ratio of CC genotypes in BC patients was 49% higher than in controls (1.49 [1.23-1.81]). Serum MPP7 protein levels were significantly higher in BC patients than in controls (p < 0.001). The protein level of the CC genotype was the highest, and that of the CT and TT genotypes decreased in turn (both p < 0.001).

CONCLUSIONS

Our results linked SNP rs1937810 to the susceptibility of BC and the clinical features of BC patients. This SNP is also proved to be significantly related to the serum level of protein MPP7 in both BC patients and controls.

Methyltransferase Setdb1 Promotes Osteoblast Proliferation by Epigenetically Silencing Macrod2 with the Assistance of Atf7ip

Bone loss caused by mechanical unloading is a threat to prolonged space flight and human health. Epigenetic modifications play a crucial role in varied biological processes, but the mechanism of histone modification on unloading-induced bone loss has rarely been studied. Here, we discovered for the first time that the methyltransferase Setdb1 was downregulated under the mechanical unloading both in vitro and in vivo so as to attenuate osteoblast proliferation. Furthermore, we found these interesting processes depended on the repression of Macrod2 expression triggered by Setdb1 catalyzing the formation of H3K9me3 in the promoter region. Mechanically, we revealed that Macrod2 was upregulated under mechanical unloading and suppressed osteoblast proliferation through the GSK-3β/β-catenin signaling pathway. Moreover, Atf7ip cooperatively contributed to osteoblast proliferation by changing the localization of Setdb1 under mechanical loading. In summary, this research elucidated the role of the Atf7ip/Setdb1/Macrod2 axis in osteoblast proliferation under mechanical unloading for the first time, which can be a potential protective strategy against unloading-induced bone loss.

Identification of PDXDC1 as a novel pleiotropic susceptibility locus shared between lumbar spine bone mineral density and birth weight

An increasing number of epidemiological studies have suggested that birth weight (BW) may be a determinant of bone health later in life, although the underlying genetic mechanism remains unclear. Here, we applied a pleiotropic conditional false discovery rate (cFDR) approach to the genome-wide association study (GWAS) summary statistics for lumbar spine bone mineral density (LS BMD) and BW, aiming to identify novel susceptibility variants shared between these two traits. We detected 5 novel potential pleiotropic loci which are located at or near 7 different genes (NTAN1, PDXDC1, CACNA1G, JAG1, FAT1P1, CCDC170, ESR1), among which PDXDC1 and FAT1P1 have not previously been linked to these phenotypes. To partially validate the findings, we demonstrated that the expression of PDXDC1 was dramatically reduced in ovariectomized (OVX) mice in comparison with sham-operated (SHAM) mice in both the growth plate and trabecula bone. Furthermore, immunohistochemistry assay with serial sections showed that both osteoclasts and osteoblasts express PDXDC1, supporting its potential role in bone metabolism. In conclusion, our study provides insights into some shared genetic mechanisms for BMD and BW as well as a novel potential therapeutic target for the prevention of OP in the early stages of the disease development. KEY MESSAGES : We investigated pleiotropy-informed enrichment between LS BMD and BW. We identified genetic variants related to both LS BMD and BW by utilizing a cFDR approach. PDXDC1 is a novel pleiotropic gene which may be related to both LS BMD and BW. Elevated expression of PDXDC1 is related to higher BMD and lower ratio n-6/n-3 PUFA indicating a bone protective effect of PDXDC1.

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.

Biobank data for studying the genetic architecture of osteoporosis and developing genetic risk scores

Osteoporosis is a chronic systemic disease of the skeleton, characterized by a decrease in bone mass and an impairment of bone microarchitecture, which can lead to a decrease in bone strength and an increase in the risk of minor trauma fractures. Osteoporosis is diagnosed on the basis of bone mineral density (BMD). BMD is characterized by high heritability that ranges according to various sources from 50 to 85%. As in the case of other complex traits, the most common approach to searching for genetic variants that affect BMD is a genome-wide association study. The lower effect size or frequency of a variant is, the larger the sample size is required to achieve statistically significant data on associations. Therefore, the studies involving hundreds of thousands of participants based on biobank data can identify the largest number of variants associated with BMD. In addition, biobank data are used in the development of genetic risk scores for osteoporosis that can be used both in combination with existing prognosis algorithms and independently of them. The aim of this review was to present the most significant studies of osteoporosis genetics, including those based on biobank data and genome-wide association studies, as well as studies on the genetic risk scores and the contribution of rare variants.

The osteoporosis risk variant rs9820407 at 3p22.1 acts as an allele-specific enhancer to regulate CTNNB1 expression by long-range chromatin loop formation

Previous powerful genome-wide association studies (GWASs) and whole-genome sequencing have identified multiple single-nucleotide polymorphisms (SNPs) located over 69 kb upstream of CTNNB1 at 3p22.1 locus associated with osteoporosis. The CTNNB1 gene encodes β-catenin that is an integral part of adherens junctions and the primary mediator of the canonical Wnt signaling pathway. The causal variants and underlying molecular mechanisms of the osteoporosis susceptibility locus 3p22.1 remains unknown. Through comprehensive computational analyses, including expression quantitative trait locus (eQTL), high-throughput chromatin interaction (Hi-C), epigenomic and functional annotation, four enhancer SNPs (rs9820407, rs9878224, rs454690 and rs9832204) were prioritized as potential causal SNPs at 3p22.1 for osteoporosis. Rs9820407 displayed the strongest enhancer activity in dual-luciferase assays. Specifically, the minor rs9820407-A can preferentially bind transcription factor FOXC1, elevate the enhancer activity and increase CTNNB1 expression. The architectural protein CTCF was presumably involved in long-range chromatin interaction between rs9820407 and CTNNB1. Our study provided a mechanistic insight into how noncoding enhancer SNP rs9820407 distally regulates CTNNB1 expression and modulates osteoporosis risk.

Gut microbiota is associated with bone mineral density : an observational and genome-wide environmental interaction analysis in the UK Biobank cohort

AIMS

Despite the interest in the association of gut microbiota with bone health, limited population-based studies of gut microbiota and bone mineral density (BMD) have been made. Our aim is to explore the possible association between gut microbiota and BMD.

METHODS

A total of 3,321 independent loci of gut microbiota were used to calculate the individual polygenic risk score (PRS) for 114 gut microbiota-related traits. The individual genotype data were obtained from UK Biobank cohort. Linear regressions were then conducted to evaluate the possible association of gut microbiota with L1-L4 BMD (n = 4,070), total BMD (n = 4,056), and femur total BMD (n = 4,054), respectively. PLINK 2.0 was used to detect the single-nucleotide polymorphism (SNP) × gut microbiota interaction effect on the risks of L1-L4 BMD, total BMD, and femur total BMD, respectively.

RESULTS

We detected five, three, and seven candidate gut microbiota-related traits for L1-L4 BMD, total BMD, and femur BMD, respectively, such as genus Dialister (p = 0.004) for L1-L4 BMD, and genus Eisenbergiella (p = 0.046) for total BMD. We also detected two common gut microbiota-related traits shared by L1-L4 BMD, total BMD, and femur total BMD, including genus Escherichia Shigella and genus Lactococcus. Interaction analysis of BMD detected several genes that interacted with gut microbiota, such as phospholipase D1 (PLD1) and endomucin (EMCN) interacting with genus Dialister in total BMD, and COL12A1 and Discs Large MAGUK Scaffold Protein 2 (DLG2) interacting with genus Lactococcus in femur BMD.

CONCLUSION

Our results suggest associations between gut microbiota and BMD, which will be helpful to further explore the regulation mechanism and intervention gut microbiota of BMD. Cite this article: Bone Joint Res 2021;10(11):734-741.

Twelve years of GWAS discoveries for osteoporosis and related traits: advances, challenges and applications

Osteoporosis is a common skeletal disease, affecting ~200 million people around the world. As a complex disease, osteoporosis is influenced by many factors, including diet (e.g. calcium and protein intake), physical activity, endocrine status, coexisting diseases and genetic factors. In this review, we first summarize the discovery from genome-wide association studies (GWASs) in the bone field in the last 12 years. To date, GWASs and meta-analyses have discovered hundreds of loci that are associated with bone mineral density (BMD), osteoporosis, and osteoporotic fractures. However, the GWAS approach has sometimes been criticized because of the small effect size of the discovered variants and the mystery of missing heritability, these two questions could be partially explained by the newly raised conceptual models, such as omnigenic model and natural selection. Finally, we introduce the clinical use of GWAS findings in the bone field, such as the identification of causal clinical risk factors, the development of drug targets and disease prediction. Despite the fruitful GWAS discoveries in the bone field, most of these GWAS participants were of European descent, and more genetic studies should be carried out in other ethnic populations to benefit disease prediction in the corresponding population.

Assessing causal relationship from gut microbiota to heel bone mineral density

Recent studies have demonstrated the important role played by gut microbiota in regulating bone development, but the evidence of such causal relationship is still sparse in human population. The aim of this study is to assess the causal relationship from gut microbiota to bone development and to identify specific causal bacteria taxa via a Mendelian randomization (MR) approach. A genome-wide association study (GWAS) summary statistic based two-sample MR analysis was performed. Summary statistics of microbiome GWAS (MGWAS) in 1126 twin pairs of the TwinsUK study was used as discovery sample, and the MGWAS in 984 Dutch participants from the LifeLines-DEEP cohort was used as replication sample. Estimated heel bone mineral density (eBMD) GWAS in 426,824 participants from the UK biobank (UKB) cohort was used as outcome. Bacteria were grouped into taxa features at both order and family levels. In the discovery sample, a total of 25 bacteria features including 9 orders and 16 families were analyzed. Fourteen features (5 orders + 9 families) were nominally significant, including 5 orders (Bacteroidales, Clostridiales, Lactobacillales, Pasteurellales and Verrucomicrobiales) and 9 families (Bacteroidaceae, Clostridiaceae, Lachnospiraceae, Mogibacteriaceae, Pasteurellaceae, Porphyromonadaceae, Streptococcaceae, Verrucomicrobiaceae and Veillonellaceae). One order Clostridiales and its child taxon, family Lachnospiraceae, were successfully replicated in the replication sample (Clostridiales P_discovery = 3.32 × 10^-3P_replication = 7.29 × 10^-3; Lachnospiraceae P_discovery = 0.03 P_replication = 7.29 × 10^-3). Our findings provided evidence of causal relationship from microbiota to bone development, as well as identified specific bacteria taxa that regulated bone mass variation, thus providing new insights into the microbiota mediated bone development mechanism.

Three pleiotropic loci associated with bone mineral density and lean body mass

Both bone mineral density (BMD) and lean body mass (LBM) are important physiological measures with strong genetic determination. Besides, BMD and LBM might have common genetic factors. Aiming to identify pleiotropic genomic loci underlying BMD and LBM, we performed bivariate genome-wide association study meta-analyses of femoral neck bone mineral density and LBM at arms and legs, and replicated in the large-scale UK Biobank cohort sample. Combining the results from discovery meta-analysis and replication sample, we identified three genomic loci at the genome-wide significance level (p < 5.0 × 10-8): 2p23.2 (lead SNP rs4477866, discovery p = 3.47 × 10-8, replication p = 1.03 × 10-4), 16q12.2 (rs1421085, discovery p = 2.04 × 10-9, replication p = 6.47 × 10-14) and 18q21.32 (rs11152213, discovery p = 3.47 × 10-8, replication p = 6.69 × 10-6). Our findings not only provide useful insights into lean mass and bone mass development, but also enhance our understanding of the potential genetic correlation between BMD and LBM.

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.

Gene Expression and RNA Splicing Imputation Identifies Novel Candidate Genes Associated with Osteoporosis

CONTEXT

Though genome-wide association studies (GWASs) have identified hundreds of genetic variants associated with osteoporosis related traits, such as bone mineral density (BMD) and fracture, it remains a challenge to interpret their biological functions and underlying biological mechanisms.

OBJECTIVE

Integrate diverse expression quantitative trait loci and splicing quantitative trait loci data with several powerful GWAS datasets to identify novel candidate genes associated with osteoporosis.

DESIGN, SETTING, AND PARTICIPANTS

Here, we conducted a transcriptome-wide association study (TWAS) for total body BMD (TB-BMD) (n = 66 628 for discovery and 7697 for validation) and fracture (53 184 fracture cases and 373 611 controls for discovery and 37 857 cases and 227 116 controls for validation), respectively. We also conducted multi-SNP-based summarized mendelian randomization analysis to further validate our findings.

RESULTS

In total, we detected 88 genes significantly associated with TB-BMD or fracture through expression or ribonucleic acid splicing. Summarized mendelian randomization analysis revealed that 78 of the significant genes may have potential causal effects on TB-BMD or fracture in at least 1 specific tissue. Among them, 64 genes have been reported in previous GWASs or TWASs for osteoporosis, such as ING3, CPED1, and WNT16, as well as 14 novel genes, such as DBF4B, GRN, TMUB2, and UNC93B1.

CONCLUSIONS

Overall, our findings provide novel insights into the pathogenesis mechanisms of osteoporosis and highlight the power of a TWAS to identify and prioritize potential causal genes.

Pleiotropic genomic variants at 17q21.31 associated with bone mineral density and body fat mass: a bivariate genome-wide association analysis

Osteoporosis and obesity are two severe complex diseases threatening public health worldwide. Both diseases are under strong genetic determinants as well as genetically correlated. Aiming to identify pleiotropic genes underlying obesity and osteoporosis, we performed a bivariate genome-wide association (GWA) meta-analysis of hip bone mineral density (BMD) and total body fat mass (TBFM) in 12,981 participants from seven samples, and followed by in silico replication in the UK biobank (UKB) cohort sample (N = 217,822). Combining the results from discovery meta-analysis and replication sample, we identified one novel locus, 17q21.31 (lead SNP rs12150327, NC_000017.11:g.44956910G > A, discovery bivariate P = 4.83 × 10^-9, replication P = 5.75 × 10^-5) at the genome-wide significance level (ɑ = 5.0 × 10^-8), which may have pleiotropic effects to both hip BMD and TBFM. Functional annotations highlighted several candidate genes, including KIF18B, C1QL1, and PRPF19 that may exert pleiotropic effects to the development of both body mass and bone mass. Our findings can improve our understanding of the etiology of osteoporosis and obesity, as well as shed light on potential new therapies.

Twelve New Genomic Loci Associated With Bone Mineral Density

Aiming to identify more genomic loci associated with bone mineral density (BMD), we conducted a joint association analysis of 2 genome-wide association study (GWAS) by the integrative association method multi-trait analysis of GWAS (MTAG). The first one is the single GWAS of estimated heel BMD (eBMD) in the UK biobank (UKB) cohort (N = 426,824), and the second one is the GWAS meta-analysis of total body BMD (TB-BMD) in 66,628 participants from 30 studies. Approximate conditional association analysis was performed in the identified novel loci to identify secondary association signal. Statistical fine-mapping was conducted to prioritize plausible credible risk variants (CRVs). Candidate genes were prioritized based on the analyses of cis- expression quantitative trait locus (cis-eQTL) and cis-protein QTL (cis-pQTL) information as well as the functional category of the SNP. By integrating the information carried in over 490,000 participants, this largest joint analysis of BMD GWAS identified 12 novel genomic loci at the genome-wide significance level (GWS, p = 5.0 × 10^-8), nine of which were for eBMD and four were for TB-BMD, explaining an additional 0.11 and 0.23% heritability for the two traits, respectively. These loci include 1p33 (lead SNP rs10493130, p_eBMD = 3.19 × 10^-8), 5q13.2 (rs4703589, p_eBMD = 4.78 × 10^-8), 5q31.3 (rs9324887, p_TB-BMD = 1.36 × 10^-9), 6p21.32 (rs6905837, p_eBMD = 3.32 × 10^-8), 6q14.1 (rs10806234, p_eBMD = 2.63 × 10^-8), 7q21.11 (rs10806234, p_TB-BMD = 3.37 × 10^-8), 8q24.12 (rs11995866, p_eBMD = 6.72 × 10^-9), 12p13.31 (rs1639122, p_eBMD = 4.43 × 10^-8), 12p12.1 (rs58489179, p_eBMD = 4.74 × 10^-8), 12q24.23 (rs75499226, p_eBMD = 1.44 × 10^-8), 19q13.31 (rs7255083, p_TB-BMD = 2.18 × 10^-8) and 22q11.23 (rs13056137, p_TB-BMD = 2.54 × 10^-8). All lead SNPs in these 12 loci are nominally significant in both original studies as well as consistent in effect direction between them, providing solid evidence of replication. Approximate conditional analysis identified one secondary signal in 5q13.2 (rs11738874, p_conditional = 5.06 × 10^-9). Statistical fine-mapping analysis prioritized 269 CRVs. A total of 65 candidate genes were prioritized, including those with known biological function to bone development (such as FGF1, COL11A2 and DEPTOR). Our findings provide novel insights into a better understanding of the genetic mechanism underlying bone development as well as candidate genes for future functional investigation.

Split hand/foot malformation associated with 20p12.1 deletion: A case report

Split hand/foot malformation (SHFM) or ectrodactyly is a rare congenital disorder affecting limb development characterized by clinical and genetic heterogeneity. SHFM is usually inherited as an autosomal dominant trait with incomplete penetrance. Isolated and syndromic forms are described. The extent of associated malformations is highly variable and multiple syndromes with clinical and genetic overlap have been described. We report here a 28 year-old man presenting with SHFM, sparse hair and widespread freckles. Array-CGH identified a 450 kb de novo 20p12.1 microdeletion encompassing three exons (exon 6 to 8) of MACROD2. Although MACROD2 mutations have not been associated with limb malformation until now, it is located next to KIF16B, which is involved in fibroblast growth factor receptor (FGFR) signaling. Additionally, the deletion encompassed a histone modification H3K27ac mark, known as a provider of quantitative readout of promoter and enhancer activity during human limb development. Altogether, these findings suggest that the 20p12.1 CNV is causative of SHFM in the present case through disturbance of regulatory elements functioning.