Metastatic site influences driver gene function in pancreatic cancer

Driver gene mutations can increase the metastatic potential of the primary tumor1-3, but their role in sustaining tumor growth at metastatic sites is poorly understood. A paradigm of such mutations is inactivation of SMAD4 - a transcriptional effector of TGFβ signaling - which is a hallmark of multiple gastrointestinal malignancies4,5. SMAD4 inactivation mediates TGFβ's remarkable anti- to pro-tumorigenic switch during cancer progression and can thus influence both tumor initiation and metastasis6-14. To determine whether metastatic tumors remain dependent on SMAD4 inactivation, we developed a mouse model of pancreatic ductal adenocarcinoma (PDAC) that enables Smad4 depletion in the pre-malignant pancreas and subsequent Smad4 reactivation in established metastases. As expected, Smad4 inactivation facilitated the formation of primary tumors that eventually colonized the liver and lungs. By contrast, Smad4 reactivation in metastatic disease had strikingly opposite effects depending on the tumor's organ of residence: suppression of liver metastases and promotion of lung metastases. Integrative multiomic analysis revealed organ-specific differences in the tumor cells' epigenomic state, whereby the liver and lungs harbored chromatin programs respectively dominated by the KLF and RUNX developmental transcription factors, with Klf4 depletion being sufficient to reverse Smad4's tumor-suppressive activity in liver metastases. Our results show how epigenetic states favored by the organ of residence can influence the function of driver genes in metastatic tumors. This organ-specific gene-chromatin interplay invites consideration of anatomical site in the interpretation of tumor genetics, with implications for the therapeutic targeting of metastatic disease.

MLL3 regulates the CDKN2A tumor suppressor locus in liver cancer

Mutations in genes encoding components of chromatin modifying and remodeling complexes are among the most frequently observed somatic events in human cancers. For example, missense and nonsense mutations targeting the mixed lineage leukemia family member 3 (MLL3, encoded by KMT2C) histone methyltransferase occur in a range of solid tumors, and heterozygous deletions encompassing KMT2C occur in a subset of aggressive leukemias. Although MLL3 loss can promote tumorigenesis in mice, the molecular targets and biological processes by which MLL3 suppresses tumorigenesis remain poorly characterized. Here we combined genetic, epigenomic, and animal modeling approaches to demonstrate that one of the mechanisms by which MLL3 links chromatin remodeling to tumor suppression is by co-activating the Cdkn2a tumor suppressor locus. Disruption of Kmt2c cooperates with Myc overexpression in the development of murine hepatocellular carcinoma (HCC), in which MLL3 binding to the Cdkn2a locus is blunted, resulting in reduced H3K4 methylation and low expression levels of the locus-encoded tumor suppressors p16/Ink4a and p19/Arf. Conversely, elevated KMT2C expression increases its binding to the CDKN2A locus and co-activates gene transcription. Endogenous Kmt2c restoration reverses these chromatin and transcriptional effects and triggers Ink4a/Arf-dependent apoptosis. Underscoring the human relevance of this epistasis, we found that genomic alterations in KMT2C and CDKN2A were associated with similar transcriptional profiles in human HCC samples. These results collectively point to a new mechanism for disrupting CDKN2A activity during cancer development and, in doing so, link MLL3 to an established tumor suppressor network.

Abstract 3512: Organ-specific effects of Smad4 in metastatic pancreatic cancer

Primary and metastatic tumors typically share the same driver gene mutations, but it is unclear if these mutations are functionally relevant across different anatomical sites. Among such mutations, inactivation of the tumor suppressor gene SMAD4 is a hallmark of pancreatic and other gastrointestinal cancers and has been associated with metastatic disease. Here, we develop a mouse model of pancreatic ductal adenocarcinoma that enables Smad4 depletion in the pre-malignant Kras-mutant pancreas and subsequent Smad4 reactivation in late-stage metastatic tumors. Whereas early Smad4 deficiency facilitated tumor formation, later Smad4 restoration had unexpected organ-specific outcomes: no effect on primary tumor growth, suppression of liver metastases, and promotion of lung metastases. Integrative multiomic analysis revealed a near-universal genomic deletion of the Cdkn2a/b locus and organ-specific changes in the tumor cells’ epigenomic state. In particular, the liver and lung differentially favored KLF vs. RUNX dominated chromatin programs, which were confirmed to underpin the divergent effects of Smad4 restoration in functional studies. Our results show how organ-dependent epigenomic changes lead to altered driver gene function in metastatic disease. This organ-specific gene-chromatin interplay may be a generalizable principle in cancer biology and invites a revised paradigm for precision oncology that considers anatomical site in the interpretation of tumor genetics.

Citation Format: Kaloyan M. Tsanov, Francisco M. Barriga, Yu-Jui Ho, Direna Alonso-Curbelo, Geulah Livshits, Richard Koche, Timour Baslan, Alexandra N. Wuest, Janelle Simon, Sha Tian, Wei Luan, Scott W. Lowe. Organ-specific effects of Smad4 in metastatic pancreatic cancer [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2023; Part 1 (Regular and Invited Abstracts); 2023 Apr 14-19; Orlando, FL. Philadelphia (PA): AACR; Cancer Res 2023;83(7_Suppl):Abstract nr 3512.

MACHETE identifies interferon-encompassing chromosome 9p21.3 deletions as mediators of immune evasion and metastasis

The most prominent homozygous deletions in cancer affect chromosome 9p21.3 and eliminate CDKN2A/B tumor suppressors, disabling a cell-intrinsic barrier to tumorigenesis. Half of 9p21.3 deletions, however, also encompass a type I interferon (IFN) gene cluster; the consequences of this co-deletion remain unexplored. To functionally dissect 9p21.3 and other large genomic deletions, we developed a flexible deletion engineering strategy, MACHETE (molecular alteration of chromosomes with engineered tandem elements). Applying MACHETE to a syngeneic mouse model of pancreatic cancer, we found that co-deletion of the IFN cluster promoted immune evasion, metastasis and immunotherapy resistance. Mechanistically, IFN co-deletion disrupted type I IFN signaling in the tumor microenvironment, leading to marked changes in infiltrating immune cells and escape from CD8+ T-cell surveillance, effects largely driven by the poorly understood interferon epsilon. These results reveal a chromosomal deletion that disables both cell-intrinsic and cell-extrinsic tumor suppression and provide a framework for interrogating large deletions in cancer and beyond.

A gene-environment-induced epigenetic program initiates tumorigenesis

Tissue damage increases the risk of cancer through poorly understood mechanisms1. In mouse models of pancreatic cancer, pancreatitis associated with tissue injury collaborates with activating mutations in the Kras oncogene to markedly accelerate the formation of early neoplastic lesions and, ultimately, adenocarcinoma2,3. Here, by integrating genomics, single-cell chromatin assays and spatiotemporally controlled functional perturbations in autochthonous mouse models, we show that the combination of Kras mutation and tissue damage promotes a unique chromatin state in the pancreatic epithelium that distinguishes neoplastic transformation from normal regeneration and is selected for throughout malignant evolution. This cancer-associated epigenetic state emerges within 48 hours of pancreatic injury, and involves an 'acinar-to-neoplasia' chromatin switch that contributes to the early dysregulation of genes that define human pancreatic cancer. Among the factors that are most rapidly activated after tissue damage in the pre-malignant pancreatic epithelium is the alarmin cytokine interleukin 33, which recapitulates the effects of injury in cooperating with mutant Kras to unleash the epigenetic remodelling program of early neoplasia and neoplastic transformation. Collectively, our study demonstrates how gene-environment interactions can rapidly produce gene-regulatory programs that dictate early neoplastic commitment, and provides a molecular framework for understanding the interplay between genetic and environmental cues in the initiation of cancer.

An in-depth study of the associations between osteoarthritis- and osteoporosis-related phenotypes at different skeletal locations

The relationship between OA and osteoporosis characteristics remains controversial. This study revealed that age-adjusted hand OA is associated with lower hand/arm BMD levels. Wrist fracture occurrence is associated with increased OA hand scores and low arm BMD. Conversely, age-adjusted knee and spine OA is associated with high spine, hip, and total BMDs.

INTRODUCTION

Osteoarthritis (OA) and osteoporosis are two common musculoskeletal diseases which contribute a high burden of disability, yet assessments of their relationship remains controversial. The aim of this study was to clarify the association between bone mineral densities (BMD) of the hand, arm, spine, hip, and total body, and OA of the hand and knee and lumbar disc degeneration in two different ethnic groups.

METHODS

Radiographic assessments of the hand, knee, and spine were collected and coded for joint space narrowing, osteophytes, and the Kellgren-Lawrence score from Chuvashian (n = 1504) and British (n = 2280) individuals. BMD measurements of standard skeletal sites were estimated by dual X-ray absorptiometry. Age- and familial-adjusted regression analyses were conducted to determine associations.

RESULTS

Knee OA affection was positively associated with elevated hip, spine, and total body BMD levels (p < 0.001). Additionally, disc degeneration phenotypes showed significant positive associations with the hip, spine, and total BMD (p < 0.001). However, increased hand OA scores was significantly negatively correlated with arm and hand BMD measurements in males and females in both samples (p < 0.001). Additionally, higher hand OA scores were significantly associated with wrist fracture.

CONCLUSIONS

We discovered a clear pattern of association between hand OA and low hand and arm BMD, with increased risk of wrist fracture, as well as reproducing previous associations between knee and spine OA and elevated spine, hip, and total body BMD. It appears that hand OA manifests differently in comparison to hip and knee OA.

A Gain-of-Function p53-Mutant Oncogene Promotes Cell Fate Plasticity and Myeloid Leukemia through the Pluripotency Factor FOXH1

Mutations in the TP53 tumor suppressor gene are common in many cancer types, including the acute myeloid leukemia (AML) subtype known as complex karyotype AML (CK-AML). Here, we identify a gain-of-function (GOF) Trp53 mutation that accelerates CK-AML initiation beyond p53 loss and, surprisingly, is required for disease maintenance. The Trp53^R172H mutation (TP53^R175H in humans) exhibits a neomorphic function by promoting aberrant self-renewal in leukemic cells, a phenotype that is present in hematopoietic stem and progenitor cells (HSPC) even prior to their transformation. We identify FOXH1 as a critical mediator of mutant p53 function that binds to and regulates stem cell-associated genes and transcriptional programs. Our results identify a context where mutant p53 acts as a bona fide oncogene that contributes to the pathogenesis of CK-AML and suggests a common biological theme for TP53 GOF in cancer. SIGNIFICANCE: Our study demonstrates how a GOF p53 mutant can hijack an embryonic transcription factor to promote aberrant self-renewal. In this context, mutant Trp53 functions as an oncogene to both initiate and sustain myeloid leukemia and suggests a potential convergent activity of mutant Trp53 across cancer types.This article is highlighted in the In This Issue feature, p. 813.

Traditional Approaches for Company Valuation Are Flawed for Valuing In Vivo Gene Therapy Companies

The era of gene therapy has begun. In recent years, potentially breakthrough datasets and rapidly expanding company pipelines have begun to overshadow the unfulfilled promise characteristic of the gene therapy sector in decades prior. One barometer for progress in the space can be seen in stock markets, where NASDAQ-listed in vivo gene therapy companies we follow have increased from 4 companies with $1.9 billion in market capitalization on January 31, 2014, to 24 companies with $30.5 billion in market capitalization on October 31, 2018. For many in the financial community, a tangible signal for the emergence of the broader gene therapy space is the recent notable mergers and acquisitions activity, a signal that previously heralded the arrival of blockbuster biotechnologies like monoclonal antibodies. Notably, Novartis' $8.7 billion acquisition of in vivo adeno-associated virus 9-based gene therapy player, AveXis, earlier this year has focused many on looking for new investment opportunities in the space, thereby increasing interest in the valuation of gene therapy companies. This perspective discusses the theoretical underpinnings of company valuation and explains why traditional approaches have limitations when valuing in vivo gene therapy companies, which produce single treatments that may achieve durable or curative benefits. We use the AveXis case study to illustrate certain points on the valuation of breakthrough innovation that we think have broader applicability throughout the in vivo gene therapy space. This publication is the first in a three-part series. Future discussions in this series on in vivo gene therapy companies will explore real-world approaches and considerations that have already proven successful in mitigating the limitations of traditional valuation approaches as well as those that may soon emerge.

Arid1a restrains Kras-dependent changes in acinar cell identity

Mutations in members of the SWI/SNF chromatin remodeling family are common events in cancer, but the mechanisms whereby disruption of SWI/SNF components alters tumorigenesis remain poorly understood. To model the effect of loss of function mutations in the SWI/SNF subunit Arid1a in pancreatic ductal adenocarcinoma (PDAC) initiation, we directed shRNA triggered, inducible and reversible suppression of Arid1a to the mouse pancreas in the setting of oncogenic Kras^G12D. Arid1a cooperates with Kras in the adult pancreas as postnatal silencing of Arid1a following sustained Kras^G12D expression induces rapid and irreversible reprogramming of acinar cells into mucinous PDAC precursor lesions. In contrast, Arid1a silencing during embryogenesis, concurrent with Kras^G12D activation, leads to retention of acinar cell fate. Together, our results demonstrate Arid1a as a critical modulator of Kras-dependent changes in acinar cell identity, and underscore an unanticipated influence of timing and genetic context on the effects of SWI/SNF complex alterations in epithelial tumorigenesis.

The pancreas produces many different hormones, as well as several substances important for digestion. To perform these roles, the pancreas contains different types of cells; for example, acinar cells make digestive enzymes that help to break down food. But, like other cells in the body, pancreatic cells can accumulate mutations in their DNA that cause them to divide, acquire an altered identity and form a cancerous tumor.

The DNA of cells is packed into a structure called chromatin. While the DNA sequence is essentially the same across all normal cells of a given individual, chromatin can be more or less compacted in the different cell types that comprise our body tissues. A collection of proteins called the SWI/SNF complex can reorganize the chromatin to change how tightly the DNA is packed. This determines which genes in the DNA are accessible and can be activated, and which ones cannot.

Around 25% of pancreatic cancers contain mutations in genes that produce proteins of the SWI/SNF complex. These mutations normally occur with an additional mutation that over-activates the gene that produces a potentially cancer-causing protein called Kras. Livshits et al. have now genetically engineered mice to investigate how one such SWI/SNF complex protein, called Arid1a, affects how pancreatic cancer develops using a genetic approach that made possible to temporarily halt the production of Arid1a in acinar cells by feeding these mice an antibiotic. The gene that produces Kras was also over-activated in the pancreases of the mice, making them more likely to develop cancer.

Within just two weeks of stopping the production of Arid1a, the acinar cells stopped producing digestive enzymes and started making other proteins that are typically found in cancerous cells, indicating that Arid1a is involved in maintaining the normal identity and activity of these cells. Restoring the ability of altered acinar cells to produce normal levels of Arid1a (by removing the mice from the antibiotic diet) did not reverse these changes. Biochemical experiments showed that acinar cells with reduced levels of Arid1a have altered chromatin. In particular, the genes that produce digestive enzymes, which are normally active in healthy pancreases, were less accessible in mice who had over-active Kras and reduced levels of Arid1a.

The results presented by Livshits et al. provide the first evidence of how alterations to Arid1a can lead to irreversible changes in the identity and activity of pancreatic acinar cells. These results will need to be carefully considered by researchers who are developing treatments for cancer patients with mutations in Arid1a and other SWI/SNF proteins. In particular, methods that attempt to restore the functions of absent SWI/SNF proteins to cancer cells are unlikely to treat the cancer successfully.

Hierarchical, imbalanced pro-inflammatory cytokine networks govern the pathogenesis of chronic arthropathies

BACKGROUND

Chronic inflammatory arthropathies, such as rheumatoid arthritis (RA), spondyloarthritis, including psoriatic arthritis (PsA), ankylosing spondyloarthritis (AS), osteoarthritis (OA), and intervertebral disc degenerative disease (DDD) constitute major public health problems that are anticipated to grow significantly as the human population ages. However, many aspects concerning the molecular mechanisms underlying their onset and progression remain unclear.

DESIGN

This narrative review critically analyzes the molecular mechanisms underlying the inflammation-associated pathogenesis of the aforementioned joint diseases. This includes, in particular, the major role played by several key soluble factors (such as cytokines and the associated signaling pathways, designated as "fragile nodes") produced by local cells and recruited to the joints' immune cells, whose elimination by specific drugs has dramatically improved the diseases' symptomatology and outcome in human clinical trials or in rodent arthritis models.

HYPOTHESIS AND THE AIM OF THIS REVIEW

We hypothesize that the pathogenesis of chronic inflammatory arthropathies is governed by hierarchical, imbalanced pro-inflammatory cytokine networks (HIPICNs) (comprising a combination of fragile nodes) that are created during the development of both autoimmune (RA, PsA, and AS) and non-autoimmune (OA and DDD) disorders. The main aim of this review is to provide evidence that despite substantial pathobiological differences between these arthropathies, the HIPICNs created are quite common, thus justifying the merging of these disorders mechanistically and suggesting that these common mechanisms exist in the onset and progression of different joint diseases.

Transplantation of engineered organoids enables rapid generation of metastatic mouse models of colorectal cancer

Colorectal cancer (CRC) is a leading cause of death in the developed world, yet facile preclinical models that mimic the natural stages of CRC progression are lacking. Through the orthotopic engraftment of colon organoids we describe a broadly usable immunocompetent CRC model that recapitulates the entire adenoma-adenocarcinoma-metastasis axis in vivo. The engraftment procedure takes less than 5 minutes, shows efficient tumor engraftment in two-thirds of mice, and can be achieved using organoids derived from genetically engineered mouse models (GEMMs), wild-type organoids engineered ex vivo, or from patient-derived human CRC organoids. In this model, we describe the genotype and time-dependent progression of CRCs from adenocarcinoma (6 weeks), to local disseminated disease (11-12 weeks), and spontaneous metastasis (>20 weeks). Further, we use the system to show that loss of dysregulated Wnt signaling is critical for the progression of disseminated CRCs. Thus, our approach provides a fast and flexible means to produce tailored CRC mouse models for genetic studies and pre-clinical investigation.

Using CRISPR/Cas to study gene function and model disease in vivo

The recent discovery of the CRISPR/Cas system and repurposing of this technology to edit a variety of different genomes have revolutionized an array of scientific fields, from genetics and translational research, to agriculture and bioproduction. In particular, the prospect of rapid and precise genome editing in laboratory animals by CRISPR/Cas has generated an immense interest in the scientific community. Here we review current in vivo applications of CRISPR/Cas and how this technology can improve our knowledge of gene function and our understanding of biological processes in animal models.

Quantitative genetics of circulating Dickkopf-related protein 1 (DKK1) in community-based sample of UK twins

Dickkopf-related protein 1 (DKK1) is a major inhibitor of Wnt signalling pathway but also plays an important role in bone formation. Its circulating levels appear to correlate significantly with plasma levels of inflammatory factors, fractalkine and IL-6. This study, using a large sample of UK twins, showed that the variation of each of these factors and correlation between them was explained by the genetic factors, and indicated possible association with DKK1 gene variants.

INTRODUCTION

DKK1 is involved in the development of several inflammatory conditions related to bone and joint degradation. Our objectives were to explore the genetic contribution (heritability) to circulating DKK1 variation and its correlation with other inflammatory cytokines, interleukin 6 (IL-6) and fractalkine, and to test whether the DKK1 heritability could be attributable to single nucleotide polymorphisms (SNPs) mapped to DKK1, IL-6 and FRCT genes.

METHODS

The study included a large community-based sample of 4939 women drawn from the general UK population. Plasma samples were analysed for circulating levels of DKK1, IL-6 and fractalkine (FRCT); 65 SNPs of DKK1, IL-6 and FRCT candidate genes, with MAF >0.1, were examined. We applied variance component analysis to evaluate contribution of putative genetic (including above SNPs) and environmental factors to variation of DKK1, and its correlation with IL-6 and FRCT.

RESULTS

Putative genetic factors explained 42.2 ± 2 % of the total variation of circulating DKK1 levels, and were also significant for fractalkine and IL-6 variations. Most importantly, we report significant phenotypic (0.208 ± 0.006-0.459 ± 0.007) and genetic (0.338 ± 0.069-0.617 ± 0.033) correlations between these molecules. We found evidence suggestive of association between the DKK1 and its structural genes variants.

CONCLUSIONS

Circulating DKK1 levels correlated significantly with levels of IL-6 and FRCT, known risk factors for several inflammatory processes suggesting a potential role of DKK1 in inflammation and tissue injury. Our results suggest the contribution of genetic factors in inter-individual variation of DKK1 levels in human population. However, further studies are required to determine genetic polymorphisms affecting DKK1 variation and its correlation with IL-6 and FRCT.

Late successful weight reduction and maintenance among overweight and obese adults--a two-year retrospective study

AIMS

Our purpose was to identify those factors associated with successful late weight reduction and maintenance among overweight and obese adults who failed to achieve initial weight reduction success.

METHODS

Medical computerized files of 5254 participants, who failed to achieve ≥ 5% weight reduction after an initial 6-month period, were retrospectively analyzed to identify predictors associated with late successful weight reduction and maintenance (≥ 5% during the first and second years, respectively). Over 40 independent variables were analyzed. The main outcome was the percentage of weight change.

RESULTS

Significant predictors of late success in weight reduction were as follows: more visits to a dietitian, higher baseline BMI, and any initial weight reduction (0-5%) (OR=3.69, compared with participants who initially gained weight). The use of insulin (OR=0.499) and the presence of hypertension (OR=0.75) were significantly correlated with failure to reduce weight. Predictors of late maintenance were as follows: more visits to a dietitian, higher baseline BMI, any initial weight reduction, a younger age, not being treated with insulin (OR=0.316), and more weighings (OR=1.68).

CONCLUSIONS

A substantial sub-group of obese and overweight patients was able to reduce their weight at a slower rate than the defined successful time of 6 months. Significant specific predictors were identified. Diabetic and hypertensive patients are at a significantly higher risk of failure to reduce and maintain weight. Using regression models, we calculated the probability of successful late weight reduction. This calculation could serve as a clinical tool for a professional team.

Osteocyte control of bone remodeling: is sclerostin a key molecular coordinator of the balanced bone resorption-formation cycles?

Osteocytes, entrapped within a newly mineralized bone matrix, possess a unique cellular identity due to a specialized morphology and a molecular signature. These features endow them to serve as a bone response mechanism for mechanical stress in their microenvironment. Sclerostin, a primarily osteocyte product, is widely considered as a mechanotranduction key molecule whose expression is suppressed by mechanical loading, or it is induced by unloading. This review presents a model suggesting that sclerostin is major mediator for integrating mechanical, local, and hormonal signals, sensed by the osteocytes, in controlling the remodeling apparatus. This central role is achieved through interplay between two opposing mechanisms: (1) unloading-induced high sclerostin levels, which antagonize Wnt-canonical-β-catenin signaling in osteocytes and osteoblasts, permitting simultaneously Wnt-noncanonical and/or other pathways in osteocytes and osteoclasts, directed at bone resorption; (2) mechanical loading results in low sclerostin levels, activation of Wnt-canonical signaling, and bone formation. Therefore, adaptive bone remodeling occurring at a distinct bone compartment is orchestrated by altered sclerostin levels, which regulate the expression of the other osteocyte-specific proteins, such as RANKL, OPG, and proteins encoded by "mineralization-related genes" (DMP1, PHEX, and probably FGF23). For example, under specific terms, sclerostin regulates differential RANKL and OPG production, and creates a dynamic RANKL/OPG ratio, leading either to bone formation or resorption. It also controls the expression of PHEX, DMP1, and most likely FGF23, leading to either bone matrix mineralization or its inhibition. Such opposing up- or down-regulation of remodeling phases allows osteocytes to function as an "external unit", ensuring transition from bone resorption to bone formation.Mini Abstract: The osteocyte network plays a central role in directing bone response either to mechanical loading, or to unloading, leading correspondingly to bone formation or resorption. This review shows a key role of the osteocyte-produced sclerostin as a major mediator of the molecular mechanisms involved in the process of adaptive bone remodeling.

Accelerating cancer modeling with RNAi and nongermline genetically engineered mouse models

For more than two decades, genetically engineered mouse models have been key to our mechanistic understanding of tumorigenesis and cancer progression. Recently, the massive quantity of data emerging from cancer genomics studies has demanded a corresponding increase in the efficiency and throughput of in vivo models for functional testing of putative cancer genes. Already a mainstay of cancer research, recent innovations in RNA interference (RNAi) technology have extended its utility for studying gene function and genetic interactions, enabling tissue-specific, inducible and reversible gene silencing in vivo. Concurrent advances in embryonic stem cell (ESC) culture and genome engineering have accelerated several steps of genetically engineered mouse model production and have facilitated the incorporation of RNAi technology into these models. Here, we review the current state of these technologies and examine how their integration has the potential to dramatically enhance the throughput and capabilities of animal models for cancer.

Association between cartilage and bone biomarkers and incidence of radiographic knee osteoarthritis (RKOA) in UK females: a prospective study

OBJECTIVE

There is a need to find biochemical markers that would identify people with increased risk of developing radiographic knee osteoarthritis (RKOA). The aim of this study was to evaluate the ability of cartilage and bone biomarkers (cartilage oligomeric matrix protein (COMP), aggrecan, cellular inhibitor of apoptosis protein (cIAP), N-telopeptide-to-helix (NTx)) to predict RKOA incidence in a 10-year follow-up of UK females from the Chingford community study.

METHOD

Joint space narrowing (JSN), osteophytes (OSP) and Kellgren-Lawrence (K/L) grades were scored from radiographs of both knees at study baseline and 10 years later in 1,003 women aged 45-64. Circulating levels of biomarkers and demographic variables were measured at baseline. Statistical association analysis was conducted between the potential predictor factors measured at baseline and documentation of RKOA at 10-year follow-up.

RESULTS

Age and body mass index (BMI), were significant predictors of incidence of RKOA as assessed by K/L and OSP. Considering biomarkers, independent significant association was found between COMP circulating levels and K/L scores (Odd Ratio (OR) = 2.87, 95% Confidence Interval (CI) = 1.19-6.89, P = 0.018). Significant negative association was detected between aggrecan plasma concentrations and JSN, with OR = 0.37 (95% CI 0.15-0.89), P = 0.026.

CONCLUSIONS

Aggrecan and COMP circulating levels contribute to identification of phenotype-specific RKOA incidence. These data suggest potentially protective role of aggrecan in cartilage loss, as measured by JSN. High COMP levels are risk factors for development of RKOA, as assessed by K/L scores.

Is interaction between age-dependent decline in mechanical stimulation and osteocyte-estrogen receptor levels the culprit for postmenopausal-impaired bone formation?

Declining estrogen levels during menopause are widely considered to be a major cause of age-dependent bone loss, which is primarily manifested by increased bone resorption by osteoclasts. We present accumulating evidence supporting another aspect of metabolic bone loss, suggesting that the combined interaction between age-dependent factors, namely, estrogen deficiency and reduced day-by-day activity/mechanical stimulation, directly leads to a reduction in anabolic processes. Such decreased bone formation results in diminished bone strength and failure to maintain the load-bearing competence of a healthy skeleton and to postmenopausal osteoporosis disorder. Estrogen receptors (ERs), as mediators of estrogenic actions, are essential components of bone osteocyte and osteoblast mechano-adaptive responses. ER expression appears to be upregulated by adequate circulating estrogen levels. ERα signaling pathways participate in the mechanotransduction response through obligatory "non-genomic" actions that occur independently of estrogen binding to ER and by a potentially "genomic", estrogen-dependent mode. The experimental data indicate that cross talk between the ERα-"non-genomic" and Wnt/β-catenin signaling pathways constitutes the major regulatory mechanism. This interaction uses mechanically and ER-induced prostaglandin E2 as a mediator for the downregulation of osteocyte production of sclerostin. Sclerostin suppression, in turn, is a central prerequisite for load-induced formation and mineralization of the bone matrix. It is therefore plausible that future strategies for preventing and treating postmenopausal osteoporosis may use estrogenic compounds (such as selective estrogen receptor modulators or phytoestrogens) with physical activity, to complement antiresorptive therapy, aimed at stopping further bone loss and possibly even reversing it by stimulation of bone gain.

A significant association exists between receptor tyrosine kinase-like orphan receptor 2 gene variants and the OPG/RANKL ratio in human plasma

SUMMARY

There is a paucity of studies investigating association between ROR2 gene variants and osteoporosis and osteoarthritis-related phenotypes. The published literature suggests that osteoprotegerin (OPG) and receptor activator of nuclear factor-kB ligand (RANKL) are essential for bone metabolism and correlate with osteoarthritis manifestation and progression. The present study provides evidence of the significant association between ROR2 variants and the OPG/RANKL ratio in human plasma. The present results also suggest significant association between ROR2 polymorphisms and severity of radiographic hand osteoarthritis.

INTRODUCTION

Despite the importance of the ROR-2 in skeletal physiology, there is a paucity of studies investigating the potential association of ROR2 gene variants with phenotypes relevant to osteoporosis and osteoarthritis. On the other hand, there is a considerable body of literature suggesting that OPG and RANKL and their ratio (OPG/RANKL) are essential for regulating bone resorption. This is also correlated with osteoarthritis manifestation and progression. The present study therefore examines whether ROR2 polymorphisms may be associated with the OPG/RANKL ratio and hand osteoarthritis (HOA).

METHODS

The study was conducted in a family-based sample of 1,515 Caucasian individuals, assessed for radiographic hand osteoarthritis, using the Kellgren/Lawrence score. Of these, 865 individuals were genotyped for 19 SNPs, relatively equally covering the ROR2 locus, and their plasma levels of OPG and RANKL were assayed. The association between the selected SNPs and OPG, along with the OPG/RANKL ratio and HOA, was explored using the pedigree disequilibrium test.

RESULTS

Of the total of 57 tests, 16 nominally significant results (p < 0.05) were obtained, which is considerably more than the three normally expected for type I error. The significant association signals for all three phenotypes were mapped to the intron 1 region. The most significant results were detected between OPG/RANKL and rs7048756 (p < 0.0005) and between adjacent rs4744107 and Kellgren/Lawrence score (p = 0.006).

CONCLUSIONS

The present study provides evidence of the significant association between ROR2 variants and the OPG/RANKL ratio in human plasma and also suggests ROR2 association with HOA.

Quantitative genetic study of the circulating osteopontin in community-selected families

The study assessed contribution of genetic factors to variability of osteopontin (OPN) levels. Evidence of association of OPN levels with polymorphisms in its structural gene and integrin-binding sialoprotein gene loci was obtained. The results motivate research of OPN-related proteins and genes with respect to biomineralization and other biological processes.

INTRODUCTION

OPN is a major phosphoprotein in bone, which plays key role in regulation of bone mineralization process. It is considered as a promising biomarker for osteoarthritis and osteoporosis, and various other pathological conditions. However, the contribution of genetics and other confounding factors to OPN circulating levels variation in general population has never been specifically determined. The main aims of the present study included (1) evaluation of the putative genetic and familial factors' effect on OPN variability and (2) testing the hypothesis that OPN plasma levels are associated with the genetic polymorphisms in its structural gene locus (SPP1) and in integrin-binding sialoprotein gene locus (IBSP).

METHODS

To address these questions, we used a family-based sample of 925 apparently healthy Caucasian individuals. Association of OPN levels with three SNPs in each of the two selected gene loci was explored using pedigree disequilibrium tests.

RESULTS

Some 58% and 13% of the OPN levels variability were attributable to genetic factors and common spouse environment, respectively. Three SNPs showed nominally significant association with OPN (p < 0.05). Of these, rs2616262 linked to IBSP promoter region remained significant after correction for multiple testing (p = 0.003). Significant association of this SNP and rs10516799 (distal segment of SPP1) with OPN was confirmed in several statistical tests. Using a special modification of variance component analysis, we examined gene-gene and gene-sex interaction effects, but found non-significant confirmation for these hypotheses.

CONCLUSIONS

Further studies are required to confirm the observed results and to explore the underlying molecular and physiological mechanisms.

Quantitative genetic study of amphiregulin and fractalkine circulating levels--potential markers of arthropathies

OBJECTIVE

Amphiregulin (AREG) and Fractalkine (FRACT), are involved in a variety of normal and pathological processes, and are both suggested to be relevant to joint degeneration. The aims of the present study included (1) testing association between circulating levels of these biomarkers and joint pathologies, (2) evaluation of the putative genetic and familial factors' effect on AREG and FRACT variability.

DESIGN

The study was conducted in the family-based sample of 923 Caucasian individuals. Variance component analysis was used to assess contribution of genetic and environmental factors to variability of AREG and FRACT concentration.

RESULTS

The mean levels of FRACT were significantly higher in the affected group with arthropathies (synovial joints osteoarthritis (OA) and disc degenerative disease, DDD) then in the control group (P<0.0004). Circulating AREG levels were higher in DDD (P=0.0272). Genetic factors constituted the main source of the interindividual differences of the AREG and FRACT levels in our sample, and explained 29.68% and 41.68% of the total variation, respectively. The phenotypic correlation between AREG and FRACT was substantial (r=0.55, P=0.0001) and was associated with both common genetic and environmental factors. Specifically, 30% of the phenotypic correlation between AREG and FRACT was due to common genetic effects.

CONCLUSIONS

Further studies are required to assess relevancy of FRACT to clinical diagnosis and prognosis of arthropathies, to investigate the mechanisms behind the observed phenotypic and genetic covariation among the studied biomarkers, and to explore specific genetic polymorphisms affecting AREG and FRACT variation.

Rapid functional dissection of genetic networks via tissue-specific transduction and RNAi in mouse embryos

Using ultrasound-guided in utero infections of fluorescently traceable lentiviruses carrying RNAi or Cre recombinase into mouse embryos, we have demonstrated noninvasive, highly efficient selective transduction of surface epithelium, in which progenitors stably incorporate and propagate the desired genetic alterations. We achieved epidermal-specific infection using small generic promoters of existing lentiviral short hairpin RNA libraries, thus enabling rapid assessment of gene function as well as complex genetic interactions in skin morphogenesis and disease in vivo. We adapted this technology to devise a new quantitative method for ascertaining whether a gene confers a growth advantage or disadvantage in skin tumorigenesis. Using alpha1-catenin as a model, we uncover new insights into its role as a widely expressed tumor suppressor and reveal physiological interactions between Ctnna1 and the Hras1-Mapk3 and Trp53 gene pathways in regulating skin cell proliferation and apoptosis. Our study illustrates the strategy and its broad applicability for investigations of tissue morphogenesis, lineage specification and cancers.

Morphological and biochemical features of obesity are associated with mineralization genes' polymorphisms

BACKGROUND

Ectonucleotide pyrophosphatase/phosphodiesterase 1 (ENPP1) was recently extensively studied as a candidate gene for obesity phenotypes. As the human homologue of the mouse progressive ankylosis (ANKH) and alkaline phosphatase (ALPL) are known functional partners of ENPP1 in bone mineralization, we hypothesized that these genes may also be jointly involved in determining obesity features.

AIM

To examine the effects of the three genes, possible gene-sex and gene-gene interactions on variability of four obesity phenotypes: the body mass index (BMI), the waist-hip ratio (WHR), the epidermal growth factor receptor (EGFR), and leptin.

SUBJECTS AND METHODS

In all, 962 healthy individuals from 230 families were genotyped for 45 single nucleotide polymorphisms (SNPs). The association analysis was performed using two family based association tests (family based association test and pedigree disequilibrium test). The combined P-values of the two tests were estimated by Monte-Carlo simulations. Relative magnitude of the genetic and familial effects, gene-sex and gene-gene interactions were assessed using variance component models.

RESULTS

Associations were observed between ENPP1 polymorphisms and BMI (P=0.0037) and leptin (P=0.0068). ALPL markers were associated with WHR (P=0.0026) and EGFR (P=0.0001). The ANKH gene was associated with all four studied obesity-related traits (P<0.0184), and its effects were modulated by sex. Gene-gene interactions were not detected.

CONCLUSION

The observed pattern of association signals indicates that ANKH may have a generalized effect on adipose tissue physiology, whereas ENPP1 and ALPL affect distinct obesity features. The joint analysis of related genes and integration of the results obtained by different methods used in this research should benefit other studies of similar design.

A genome-wide association study suggests that a locus within the ataxin 2 binding protein 1 gene is associated with hand osteoarthritis: the Treat-OA consortium

To identify the susceptibility gene in hand osteoarthritis (OA) the authors used a two-stage approach genome-wide association study using two discovery samples (the TwinsUK cohort and the Rotterdam discovery subset; a total of 1804 subjects) and four replication samples (the Chingford Study, the Chuvasha Skeletal Aging Study, the Rotterdam replication subset and the Genetics, Arthrosis, and Progression (GARP) Study; a total of 3266 people). Five single-nucleotide polymorphisms (SNPs) had a likelihood of association with hand OA in the discovery stage and one of them (rs716508), was successfully confirmed in the replication stage (meta-analysis p = 1.81×10⁻⁵). The C allele conferred a reduced risk of 33% to 41% using a case–control definition. The SNP is located in intron 1 of the A2BP1 gene. This study also found that the same allele of the SNP significantly reduced bone density at both the hip and spine (p<0.01), suggesting the potential mechanism of the gene in hand OA might be via effects on subchondral bone. The authors' findings provide a potential new insight into genetic mechanisms in the development of hand OA.

Epiphyseal expansion in hand bones: association with age, sex, and hand osteoarthritis

OBJECTIVES

Previous studies showed increased femoral, neck, and tibial plateau with age in individuals with and without osteoarthritis (OA) at adjacent joints. However, the question whether epiphyseal bone enlargement is a natural phenomenon of aging or associated with OA remains open. The aim of the present study was to evaluate age- and sex-related changes in the relative size of epiphyses of long hand bones and their association with radiographic OA.

DESIGN

The data were collected from a population-based European sample in 1994 (557 individuals) and in 2002 (513 individuals). The latter sample included 253 individuals who were previously investigated in 1994. The epiphyseal index (EI), reflecting the relative size of bone epiphyses and hand OA, was evaluated from hand radiographs. Statistical analyses included multiple regression analyses and a maximum likelihood-based model-fitting technique.

RESULTS

Hand bone epiphyses increased with age and with OA. In males, the EI gradually increased during their entire life span. In females, the EI remained almost unchanged up to the age of 40, after which, it increased more rapidly than in males. Individuals with OA had higher values of EI at any age. In both sexes, epiphyseal enlargement is a predisposing factor for hand OA progression in adjacent joints. This was clearly seen in males, where old individuals with high EI values had much higher OA scores in comparison with age-matched individuals.

CONCLUSIONS

Enlargement of long bone epiphyses with age appears to be a general tendency in the human skeleton. Our study shows that the enlargement of epiphyses may also be related to OA.

Family-based association study of polymorphisms in the RUNX2 locus with hand bone length and hand BMD

Osteoporosis is characterized by reduced bone strength. Bone size and bone mineral density (BMD) are major bone strength determinants. Identification of genes affecting the variability of these traits should improve prognosis and management of osteoporosis. This research was aimed to test the hypothesis of association of radiographic hand bone length (BL) and BMD with polymorphisms in the RUNX2 locus. Four SNPs linked to the two RUNX2 promoters were genotyped in 212 nuclear Caucasian families. These SNPs and four pairwise haplotypes were tested for association with eight BL and BMD traits, adjusted for covariates. We observed significant associations between polymorphisms linked to the RUNX2 P1 promoter and BL mean values for three studied bone groups: all 18 bones, proximal and medial bones (p = 0.0118, 0.0085, and 0.0056, respectively). Mean BMD values for all 18 bones, proximal and medial bones were associated with polymorphisms linked to the RUNX2 P2 promoter (p = 0.0032, 0.0077, 0.0007, respectively). Associations with BL and BMD mean values for medial and proximal bones remained significant even after correction for multiple testing. This study provides evidence of the association between polymorphisms linked to the two RUNX2 promoters and variability of hand BL and BMD. The results suggest independent roles for the two RUNX2 promoters in the determination of the traits studied.

Search for hand osteoarthritis susceptibility locus on chromosome 6p12.3-p12.1

The existence of osteoarthritis susceptibility loci on chromosome 6 for individuals suffering from hip and knee osteoarthritis has been suggested. We determined whether radiographic hand osteoarthritis in a demographically homogeneous population of European origin can be linked to loci on chromosome 6p12.3-p12.1. Nine single nucleotide polymorphisms (SNPs) were genotyped in 764 individuals (members of 189 nuclear and more complex two- or three-generation families). Radiographic hand osteoarthritis was characterized by two traits: (1) the total individual osteoarthritis score (PC1-OA) and (2) the osteophytes score (PC1-OS), obtained from the principal components analysis of sums of the Kellgren and Lawrence grade and of the osteophyte grades, respectively, for 14 joints on each hand. The contribution of genetic and environmental factors and of covariates such as age and body mass index to hand osteoarthritis was evaluated by variance components analysis. The association between the studied traits and selected DNA markers was evaluated by three types of transmission disequilibrium tests. The parent-offspring and sib-sib correlations were statistically significant for all studied traits. The additive genetic effects for PC1-OA and PC1-OS were estimated to be 43% and 37.9%, respectively. Transmission disequilibrium tests consistently revealed a statistically significant association (p values ranged from 0.017 to 0.030) between SNP rs1508632 and PC1-OS. In the tested cohort the putative genetic factors are influential enough to determine interindividual differences regarding the extent of hand osteoarthritis. SNP rs1508632 lies in immediate proximity to the TINAG gene, implicating it as a possible hand osteoarthritis susceptibility gene.

Family-based association study of ROR2 polymorphisms with an array of radiographic hand bone strength phenotypes

For the first time the study provides evidence of association of radiographic hand bone length (BL) and bone mineral density (BMD) with polymorphisms in ROR2 gene that plays important role in skeletal development. This contributes to better understanding of bone physiology and may have application in clinical practice.

INTRODUCTION AND HYPOTHESIS

Bone size and bone mineral density (BMD) are major determinants of bone strength. Identification of genes affecting these traits' variability is important for better understanding of normal and pathological bone physiology and identification of the individuals at risk for bone fracture. This study tested the hypothesis of association of radiographic hand bone length (BL) and BMD with polymorphisms in ROR2 gene that is important in skeletal development.

METHODS

Nineteen ROR2 SNPs were genotyped in 705 individuals, belonging to 212 nuclear families. The four tagging SNPs (tSNPs) and the pairwise haplotypes between adjacent tSNPs were tested for association with series of hand BL and BMD measurements, adjusted for covariates, using family-based association tests.

RESULTS

We observed significant associations with BL and BMD mean values for all 18 studied hand bones (p = 0.0080, 0.0030), mean BL and BMD for proximal phalanges (p = 0.0218, 0.0060) and metacarpal bones (p = 0.0014, 0.0004). In the latter, the association remained significant after correction for multiple testing.

CONCLUSIONS

The region of the first through the second ROR2 introns is most likely to contain the functional polymorphism/s responsible for the observed associations. Further studies are required to identify the ROR2 functional polymorphism/s affecting bone size and BMD variation.

Lumbar disc disease shows linkage to chromosome 19 overlapping with a QTL for hand OA

OBJECTIVE

Cervical and lumbar degenerative disc disease (CDD and LDD, respectively) form part of the spine osteoarthritis (OA) phenotype and are known to be influenced by genetic factors. A genome-wide linkage analysis was performed to identify new chromosomal regions of interest.

METHODS

Dizygotic healthy female twin volunteers (n = 348) from the TwinsUK register who had magnetic resonance imaging scans 10 years ago coded for degenerative disease, were identified. Multipoint genome-wide linkage analysis was conducted using 737 highly polymorphic markers of approximate spacing 10 cM.

RESULTS

The mean age of the twins was 52 years. Significant linkage peaks (log of the odds (LOD) >3) were identified for LDD at three chromosomal regions. These included chromosome 1 (position 285 cM), chromosome 5 (position 175 cM) and chromosome 19 (position 80 cM). The peak on chromosome 19 had LOD = 4.06, and the empirical p = 6.7x10(-4) confirmed reliability of the linkage signal. It lies close to a linkage peak previously obtained by our group for hand OA.

CONCLUSIONS

This genome-wide linkage study of CDD and LDD shows evidence of linkage for LDD on chromosome 19. The region of interest is likely to harbour genes that are common to LDD and hand OA.

In vivo imaging of germinal centres reveals a dynamic open structure

Germinal centres are specialized structures wherein B lymphocytes undergo clonal expansion, class switch recombination, antibody gene diversification and affinity maturation. Three to four antigen-specific B cells colonize a follicle to establish a germinal centre and become rapidly dividing germinal-centre centroblasts that give rise to dark zones. Centroblasts produce non-proliferating centrocytes that are thought to migrate to the light zone of the germinal centre, which is rich in antigen-trapping follicular dendritic cells and CD4+ T cells. It has been proposed that centrocytes are selected in the light zone on the basis of their ability to bind cognate antigen. However, there have been no studies of germinal-centre dynamics or the migratory behaviour of germinal-centre cells in vivo. Here we report the direct visualization of B cells in lymph node germinal centres by two-photon laser-scanning microscopy in mice. Nearly all antigen-specific B cells participating in a germinal-centre reaction were motile and physically restricted to the germinal centre but migrated bi-directionally between dark and light zones. Notably, follicular B cells were frequent visitors to the germinal-centre compartment, suggesting that all B cells scan antigen trapped in germinal centres. Consistent with this observation, we found that high-affinity antigen-specific B cells can be recruited to an ongoing germinal-centre reaction. We conclude that the open structure of germinal centres enhances competition and ensures that rare high-affinity B cells can participate in antibody responses.

Genetic and environmental determinants of circulating resistin level in a community-based sample

OBJECTIVE

Resistin is a hormone secreted by adipose tissue, monocytes, bone marrow, and other tissues. It was also proclaimed as an important link between obesity and diabetes. The main objective of this study was to elucidate the contribution of a number of endogenous factors, such as sex, age, obesity characteristics, and genetic effects to the production of resistin in apparently healthy individuals. We also tested the possible relationships between circulating levels of resistin and other adipokines (leptin, interleukin-6 (IL-6), and tumor necrosis factor-alpha (TNF-alpha)).

MEASUREMENTS

The plasma levels of studied adipokines were determined by enzyme-linked immunoassay in pedigree-based sample (n = 616), and subjected to model-based quantitative genetic analysis.

RESULTS

Resistin levels were significantly higher in women than in men (3.60 +/- 2.53 vs 3.15 +/- 2.48 ng/ml, P < 0.001), and varied independently of age in either sex. Statistical-genetic analysis revealed significant familial correlations (P < 0.01) for resistin. Adjusted for covariates, 66.38 +/- 10.28% of the resistin variation was attributable to putative genetic factors. A relatively small portion of the resistin variation (11.54 +/- 5.77%) was attributable to sharing a common household environment. The remaining variation, 22.12 +/- 17.69% was due to random environmental (i.e., unmeasured non-additive genetic) effects. The results of our analysis showed modest significant correlation of resistin with TNF-alpha and IL-6, and only in some groups; thus, while resistin was correlated with TNF-alpha in men, the correlation with IL-6 was significant only in the post-menopausal women group.

CONCLUSIONS

Our observations indicate that resistin is strongly influenced by genetic factors. The high heritability estimates for resistin concentrations clearly suggest the continuing need for further molecular genetic investigations.

Osteoprotegerin plasma levels are strongly associated with polymorphisms in human homologue of the mouse progressive ankylosis (ANKH) gene

Osteoprotegerin inhibits osteoclastogenesis and plays an important role in the control of bone resorption. However, the genetic mechanisms underlying regulation of OPG levels are currently not fully elucidated. The aim of the present study was to determine whether the ANKH gene, which plays a central role in bone mineralization, contributes to the genetic regulation of OPG levels. A family-based association study used a sample of 159 ethnically homogeneous nuclear families, comprising 556 apparently healthy individuals. Statistical analyses included family aggregation analysis of OPG variation and four types of transmission disequilibrium tests. Each individual was genotyped for 11 SNPs in the ANKH gene. Four TDTs consistently showed a highly significant association between OPG levels and the intronic SNP rs875525 located between exons 6 and 7. The combined p-value for four tests to reject the null hypothesis of no association was 0.0003. Furthermore, haplotypes generated between rs875525 and two additional neighbouring SNPs (rs2291943 and rs2288474) also revealed a significant association with OPG plasma levels (p < 10(-4)-10(-3)). ANKH genetic polymorphisms in the area between SNP rs2291943 and rs2288474 are strongly associated with OPG plasma levels. The molecular mechanism underlying this association is not obvious, and therefore these results should be regarded cautiously until they are confirmed in independent studies.

Applications of pooled DNA samples to the assessment of population affinities: short tandem repeats

Pooled DNA samples have been used in association studies of Mendelian disease genes. This method involves combining equal quantities of DNA from patients and control subjects into separate pools and comparing the pools for distributions of genetic markers. In this study identical quantities of DNA from 300 individuals representing 6 populations were pooled and amplified for 296 loci using the touchdown polymerase chain reaction (PCR) method. The purpose of this study is to test the efficacy of pooled DNA markers in the reconstruction of the genetic structure of human populations. The populations sampled included Chuvash, Buryats, Kizhi, Native Americans, South Africans, and New York City whites. To test the accuracy of the allele-frequency distributions, we genotyped the Buryats and New York samples individually for six microsatellite markers and compared their frequencies to the allele frequencies derived from the electropherogram peak heights for the pooled DNA, producing a correlation of 0.9811 with a variance of less than 0.04. Two-dimensional scaling of genetic distances among the six populations produced clusters that reflected known historical relationships. A distance matrix was created using all 296 loci, and matrices based on individual chromosomes were correlated against the total matrix. As expected, the largest chromosomes had the highest correlations with the total matrix, whereas one of the smallest chromosomes, chromosome 22, had the lowest correlation and differed most from the combined STR distance matrix.

The cyclin-dependent kinase (CDK) family member PNQALRE/CCRK supports cell proliferation but has no intrinsic CDK-activating kinase (CAK) activity

The cyclin-dependent kinases (CDKs) that drive the eukaryotic cell cycle must be phosphorylated within the activation segment (T-loop) by a CDK-activating kinase (CAK) to achieve full activity. Although a requirement for CDK-activating phosphorylation is conserved throughout eukaryotic evolution, CAK itself has diverged between metazoans and budding yeast, and fission yeast has two CAKs, raising the possibility that additional mammalian enzymes remain to be identified. We report here the characterization of PNQALRE (also known as CCRK or p42), a member of the mammalian CDK family most similar to the cell-cycle effectors Cdk1 and Cdk2 and to the CAK, Cdk7. Although PNQALRE/CCRK was recently proposed to activate Cdk2, we show that the monomeric protein has no intrinsic CAK activity. Depletion of PNQALRE by >80% due to RNA interference (RNAi) impairs cell proliferation, but fails to arrest the cell cycle at a discrete point. Instead, both the fraction of cells with a sub-G(1) DNA content and cleavage of poly(ADP-ribose) polymerase (PARP) increase. PNQALRE knockdown did not diminish Cdk2 T-loop phosphorylation in vivo or decrease CAK activity of a cell extract. In contrast, depletion of Cdk7 by RNAi causes a proportional decrease in the ability of an extract to activate recombinant Cdk2. Our data do not support the proposed function of PNQALRE/CCRK in activating CDKs, but instead reinforce the notion of Cdk7 as the major, and to date the only, CAK in mammalian cells.

Quantitative genetics of circulating molecules associated with bone metabolism: a review

This paper reviews recent advances in the studies of various biochemical factors (biomarkers) involved in bone metabolism and remodeling. The collected data in this area suggest the existence of complex and multilevel relationships between calciotropic hormones, various cytokines and growth factors. The paper summarizes the data on the magnitude of the familial and genetic effects on the interindividual variation in circulating levels of many of these biomarkers. The majority of the cited heritability estimates are well above 20%, reaching up to 80% for some cytokines (e.g., TNFalpha and VEGF). These estimates point to potential targets for the identification of novel quantitative trait loci involved in the control of the respective molecules variation. This information is of particular importance, because the available data on the association between specific genes/polymorphisms and the respective circulating molecules variation is still very limited. The paper also provides recent findings on the genetics of co-variation between the circulating levels of various biomarkers. It shows that only in a few instances, such as for example, between IGF-I and IGFBP-3, and IGFBP-1 and leptin, significant and substantial genetic (and environmental) correlations were found. It appears that despite the prominent strong genetic effects on variation of each of the numerous biomarkers, the pleiotropic effects are rather limited. We consider briefly some important new data obtained using the gene expression approach and microarray technique. The data, for instance, indicate that the genetic effects on bone metabolism appear to be an open system, which can be activated or modulated by external factors such as drugs, e.g., PTH. Extensive molecular genetic studies in this area are both timely and imperative to detect the specific genes affecting variation (and co-variation) of the circulating factors associated with bone metabolism.

Contribution of the familial and genetic factors on monocyte chemoattractant protein-1 variation in healthy human pedigrees

Monocyte chemoattractant protein-1 (MCP-1) is a chemokine whose circulating levels have been detected in the lesions of several diseases such as pulmonary fibrosis, rheumatoid arthritis and atherosclerosis. However, the factors involved in the regulation of its production remain largely unknown. The main aim of the present paper was to ascertain the contribution of the familial/genetic factors on the production of MCP-1 in apparently healthy individuals. We also tested the possible relationships between the plasma levels of MCP-1 and other cytokines involved in bone metabolism (receptor activator NF-kB ligand (RANKL), osteoprotegerin (OPG), interleukin-6, macrophage-colony stimulating factor, tumor necrosis factor-alpha). Using ELISA assays the cytokine levels were measured in 570 apparently healthy individuals belonging to ethnically homogeneous Caucasian families. We found that MCP-1 levels were significantly (P<0.01) correlated with RANKL (in both sexes) and with OPG only in women. The study showed that adjusted for potential covariates, 72% of the MCP-1 variance, was attributable to familial effects. About 49% was due to potential genetic factors and the rest was explained by common environmental sources shared by spouses within each family. In conclusion, our data provide reliable evidence for the substantial role of genetic factors in the determination of the phenotypic variability of MCP-1 plasma levels. The association between the osteoclastogenic cytokines and MCP-1 levels in healthy pedigrees is of special interest and might shed light on MCP-1 involvement in bone remodeling.

Heritability of circulating growth factors involved in the angiogenesis in healthy human population

The present study examined the extent of genetic and environmental influences on the populational variation of circulating growth factors (VEGF, EGF) involved in angiogenesis in healthy and ethnically homogeneous Caucasian families. The plasma levels of each of the studied biochemical indices were determined by enzyme-linked immunoassay in 478 healthy individuals aged 18-75 years. Quantitative genetic analysis showed that the VEGF and EGF variation was appreciably attributable to genetic effects, with heritability estimates of 79.9% and 48.4%, respectively. Yet, common environmental factors, shared by members of the same household, also played a significant role (P < 0.01) and explained between 20.1% and 32.6% of the variation. The present study additionally examined the covariations between these molecules and either transforming growth factor-beta 1 (TGF-beta 1) or tissue inhibitors of matrix metalloproteinases 1 (TIMP-1), likewise relevant for angiogenesis. Bivariate analysis revealed significant phenotypic correlations (P < 0.002) between all pairs of variables, thus indicating the possible existence of common genetic and environmental factors. The analysis suggested that the pleiotropic genetic effects were consistently the primary (or even the sole) source of correlation between all pairs of studied molecules. The results of our study affirm the existence of specific and common genetic pathways that commonly determine the greater part of the circulating variation of these molecules.

Association of leptin levels with obesity and blood pressure: possible common genetic variation

OBJECTIVE

To ascertain the extent to which relationships between obesity (OB) and blood pressure (BP) can be explained by an individual's leptin plasma levels.

DESIGN

Pedigree-based cross-sectional study in an apparently healthy population of European origin.

SUBJECTS

The study sample is comprised of 90 nuclear and more complex families totaling 210 male and 213 female subjects aged 18-75 y, randomly recruited in Bashkorstan Autonomic region, Russia.

MEASUREMENTS

Various fatness and fat distribution traits (including nine circumferences (CRCs), and eight skinfolds (CKFs) by anthropometry), blood pressure, and plasma leptin levels (by ELISA kits).

RESULTS

Adjustment for circulating leptin led to attenuation of the magnitude of correlations between OB and BP, regardless of trait pair and sex cohort. Some of these correlations became statistically nonsignificant. All familial effects were gone, and heritability estimates became virtually zero after adjustment of each of the OB traits and systolic blood pressure (SBP) in offspring for leptin values in parents.

CONCLUSION

BP and OB covariation is substantially mediated by circulating leptin levels. As a result, body fat has only a weak independent effect on BP variation after adjustment for leptin levels. Our findings also strongly suggest that genetic variation in body mass index, SKFs, and even body CRCs, as well as of SBP is due to genetic variation of leptin. Genetic variation of diastolic blood pressure in the present sample, however, shared very little with that of leptin.

The contribution of familial resemblance to variation in circulatory levels of tissue inhibitors of metalloproteinases and transforming growth factor-beta1

This study attempted to elucidate the genetic and environmental factors influencing interindividual variation of circulating TIMP-1, TIMP-2 and TGF-beta1 and to clarify the relationship between the latter biochemical indices and hand osteoarthritis in an ethnically homogeneous sample. Plasma levels of each of the above biochemical indices were measured in 401 healthy individuals (aged 18-75 years) belonging to 90 nuclear and more complex families. Variance component analysis showed that a major part of the interindividual differences in TGF-beta1, TIMP-1 and TIMP-2 levels was credibly attributable to genetic and familial factors. Adjusted for significant covariates, the putative genetic effects on the above three amounted to 0.40 +/- 0.10, 0.47 +/- 0.11 and 0.72 +/- 0.10, respectively. Common environmental factors, shared by members of the same household, also contributed significantly ( P < 0.01) to variation of each of the biochemical indices and explained between 27.6% (TIMP-2) and 38.7% (TGF-beta1) of their variation. A bivariate analysis revealed a strong and highly significant correlation between TIMP-1 and TGF-beta1 (r = 0.58, P < 0.001), which was due to common genetic and environmental sources (r(G) = 0.62 +/- 0.09, r(E) = 0.31 +/- 0.11, both P < 0.001). The analysis also detected modest but significant genetic correlation between TIMP-1 and TIMP-2 (r(G) = -0.307 +/- 0.108, P < 0.01). The present study evinces a strong genetic dependence for the plasma levels of both TIMPs and TGF-beta1 and provides a basis for the further analysis of genetic variation affecting and regulating the circulatory concentrations of TIMPs and TGF-beta1 in healthy humans.

Genetic effects of estrogen receptor alpha and collagen IA1 genes on the relationships of parathyroid hormone and 25 hydroxyvitamin D with bone mineral density in Caucasian women

There is a growing body of evidence that estrogen receptor alpha (ERalpha) and collagen IA1 (COLIA1) genes may affect bone mineral density (BMD) levels in postmenopausal women. In a recent study we found that the Px haplotype of the ERalpha gene (resulting from combined PvuII and XbaI restriction fragment-length polymorphisms [RFLPs] in intron 1) was associated with low radiographic phalangeal hand BMD in elderly women (62.7 +/- 6.5 years of age), of European origin. The combination of the Px haplotype and "s" allele of the COLIA1 gene (MscI RFLP in Sp1 locus) decreased BMD in these women. The major aim of the present study was to investigate whether the genetic effects of these genotypes on cancellous and cortical hand BMD, in the same elderly women (N = 122), are possibly mediated through circulating levels of parathyroid hormone (PTH) and/or 25 hydroxyvitamin D [25(OH)D], and may be related to biochemical markers of bone turnover (propeptide of type I procollagen [PICP] and osteocalcin). Multiple regression analyses of age-adjusted cancellous BMD revealed that ERalpha polymorphism and circulating levels of PTH were independent predictors of about 12.9% of its variation. Some 17.9% of cortical BMD variations were attributable to the combined effects of ERalpha polymorphism and plasma concentrations of 25(OH)D, estradiol, and PTH. The significant inverse association between PTH and BMD of both types was further confirmed by association analysis according to categorical subgroups of BMD values, as well by haplotype status. The mean difference in PTH concentrations between subjects carrying the Px haplotype (higher mean) and those lacking it (lower mean) reached 0.59 SD (P =.01). The difference in PTH levels further increased when explored in the 4 subgroups formed by combinations of polymorphic ERalpha and COLIA1 genotypes. Mean PTH of subjects carrying both the Px haplotype and "s" allele was higher by 1.52 SD (P =.001) than in subjects lacking both the Px haplotype and "s" allele. Those carrying both Px haplotype and "s" allele were also characterized by highest mean value of PICP and lowest means of 25(OH)D and BMD (both tissue types). We conclude that in the studied elderly women, the Px haplotype may be involved in causing the phenotypic expression of higher circulating levels of PTH and higher bone turnover, which, in turn, may lead to bone loss.

Search for linkage between hand osteoarthritis and 11q 12-13 chromosomal segment

OBJECTIVE

The aims of the present study were: (1) to evaluate the extent and mode of inheritance of hand osteoarthritis by using a large sample of ethnically homogeneous pedigrees of Caucasian origin; (2) to examine whether the synthetic measure of osteoarthritis according to Kellgren and Lawrence (K-L) and the more specific measure, namely, the extent of osteophytes development, have a similar putative genetic determination and pattern of biological inheritance and (3) to test the hypothesis that hand osteoarthritis dependent phenotypes are linked to the 11q 12-13 chromosomal region.

METHODS

The population of the present study comprised 1190 Chuvashians (Russian Federation) belonging to 295 nuclear families. Segregation analysis was carried out on a total sample. Sub-sample of 571 individuals was used to conduct Transmission/disequilibrium test (TDT) and model-based linkage analysis.

RESULTS

Adjusted for age, sex and other covariates, both OA phenotypes showed significant familial aggregation. The model fitting analysis strongly supported the hypothesis of a major gene effect on study traits. The inferred major gene explained about 52% of the osteophyte score (OPS) and 49% of the K-L score variation adjusted for confounding variables. The series of model-based linkage analyses and TDTs provided inconclusive evidence on possible linkage of both phenotypes to the 11q 12-13 chromosomal region.

CONCLUSIONS

We support the hypothesis of a major gene effect in heritability of hand osteoarthritis in both phenotypes. Despite the fact that some DNA markers showed statistically significant association to studied primary phenotypes, we find only weak evidence of linkage disequilibrium between hand osteoarthritis and the proximal part of the 11q 12-13 chromosomal segment (D11S1983 for K-L score and D11S1313 for OPS). The subject, however, a merit requires further investigation.

Association and linkage disequilibrium analyses suggest genetic effects of estrogen receptor alpha and collagen IA1 genes on bone mineral density in Caucasian women

Estrogen receptor alpha (ER alpha) and collagen IA1 (COLIA1) genes have been suggested as possibly implicated in reduced bone mineral density (BMD). The present study investigated the occurrence of association and linkage disequilibrium between radiographic hand BMD and polymorphic alleles of ER alpha and COLIA1 genes, in human pedigrees of a Chuvasha population in Russia. The study sample included 463 members of 113 pedigrees, mostly nuclear families. We performed association and transmission disequilibrium test (TDT) analyses of the combined PvuII and XbaI RFLPs alleles on the same chromosome (haplotype) of the ER alpha gene with BMD Z scores of cancellous or cortical bone in the hand phalanges. The association analyses were performed separately for both genders in the parental generation, i.e., 'fathers' (n = 114; average age 64.2 y) and 'mothers' (n = 122; average age 62.7 y). The Px haplotype was associated significantly with lower BMD Z scores in 'mothers' only. The difference between subjects who carried one or two copies of the Px haplotype and those lacking it was 0.68 Z scores, P = 0.003 and 0.51 Z scores, P = 0.025 for cancellous and cortical bone, respectively. Multiple linear regression model with age, height, weight, and Px haplotype status as predictors explained 26.7% and 28.3% of the total observed variance in BMD with Px haplotype as independent predictor explaining 5.9%; P = 0.002 and 3%; P = 0.028 (cancellous and cortical bone, respectively). Results of t-TDT for triads of two parents and just one of their female offspring (but not male offspring) suggested the existence of linkage disequilibrium between the two loci of Px haplotype and BMD trait (P = 0.047). No association was found between polymorphic alleles of COLIA1 gene and BMD, but 'mothers' with combined genotypes of Px haplotype of ER alpha gene and "s" allele of COLIA1 gene had the lowest mean Z scores (-0.944 and -0.788 for cancellous and cortical bone, respectively). We conclude that the Px haplotype of the ER alpha gene is associated with low BMD values in females, as the phenotype is gender dependent (the association was not observed in males), and the "s" allele of COLIA1 gene in combination with this haplotype contributes to reduced BMD.

Relationship between parameters of early growth in Israeli infants

A sample of 1931 Israeli infants was measured for body weight (WT), length (HT) and head circumference (HC) for approximately 2 years. The Count model with 3 parameters was chosen as the best fitting and most parsimonious function to approximate growth of all 3 studied traits. In the model parameter a relates to birth indices, b--to velocity of growth, and c--to rapid early childhood growth, or acceleration. Assuming a difference in growth patterns in the periods of different length, the whole sample was divided into 3 groups: 1) infants with last measurement around the age of 12 months; 2) infants with last measurement around the age of 18 months, and 3) infants with last measurement around the age of 24 months. The individuals measured up to 12 months were presented in all three groups. 27 curve fitting parameters, corresponding to 3 different follow-up intervals for WT, HT and HC were computed for each individual. A high correlation was detected between the a parameters regardless of time interval for 3 measured traits. A negative correlation was found between b and c parameters within the same time interval. A consistent positive correlation was indicated between a and b parameters, especially for body length and head circumference. A principal component analysis extracted five independent factors explaining 88.1% of the total variance. Three first factors retained parameters b and c, describing growth rate and pattern of each trait separately, namely, F1 was responsible for head circumference, F2 was a body length factor, F3 was a body weight factor. F4 extracted all birth indices, observed (HC0, HT0 and WT0) and expected (parameters a). The composition of principal factors allowed us to assume that there might be a strong involvement of a pleiotropic genetic source in determination of birth size traits and an independent genetic source controlling the pattern of growth for each trait separately.

Quantitative genetic study of radiographic hand bone size and geometry

Despite the obvious epidemiological significance of bone size (BS) and geometry (BG) traits as risk factors for osteoporotic fracture, very little is still known concerning the extent of their genetic determination. In the present paper we report the results of quantitative genetic analysis of a number of BG and BS indices, as well as of BMD measurements, obtained on a large pedigree-based sample (296 nuclear families, 1208 individuals) of plain hand radiographs. The families studied were all ethnically Caucasians (Chuvasha) living in small villages along the Volga River (Russia). The sample consisted of 636 men and 572 women, aged 18-91 years. To assess hand bone size we used the outcome of principal component analysis conducted on 48 measurements of metacarpal bones and proximal phalanges (PC-BS). Two BG indices, average metacarpal cortical index and breaking bending resistance index (BBRI), also measured on metacarpal and proximal phalanges were used. Again the outcome of the principal component PC-BBRI was examined in the genetic analysis. PC-BS measurements strongly correlated with body length (r = 0.75, P < 0.001) and weight (r = 0.39, P < 0.001), suggesting that they indeed reflected hand skeleton size. Familial correlations for all studied traits, adjusted for covariates (sex, age, etc.), were all highly significant statistically. For example, parent/offspring correlations ranged between 0.248 (P < 0.001) for phalangeal BMD and 0.385 (P < 0.001) for PC-BBRI. Maximum likelihood estimates of the variance component analysis confirmed these results, indicating that approximately 58 to 66% of the residual variance of the studied traits was attributable to genetic effects. Bivariate analysis clearly revealed that while genetic variation of the phalangeal BMD was independent of the genetic effects influencing hand BS and BG, the latter two were strongly interrelated. A substantial proportion of PC-BS and PC-BBRI variation was due to shared genetic (r(G) = 0.468 +/- 0.063) and environmental (r(E) = 0.704 +/- 0.052) factors.

Familial history, age and smoking are important risk factors for disc degeneration disease in Arabic pedigrees

The present study used computed tomography imaging to evaluate the extent and pattern of the intergenerational transmission of spinal disc degeneration disease (DDD) in complex pedigrees. Contribution of a number of the potential covariates was also studied using univariate and multivariate logistic regression analysis, as well as two types of complex segregation analysis models. Among 161 individuals studied, DDD was diagnosed in 60 individuals. The number of protruded discs varied from 1 to 4, mostly in lumbar or lumbosacral regions. The average age at onset of the disease was similar for both women (36.0 years) and men (34.8 years). The proportion of the individuals affected by the DDD status of their parents ranged from 10% in families of two healthy parents to 55.5% of two affected parents (p < 0.01). The results of the logistic regression analyses and complex segregation analysis were qualitatively the same: DDD status of parents, age and smoking were the main risk factors for disc herniation in the Arabic families we examined. All analyses showed a predominating role of the family history as a risk factor for DDD in offsprings. It showed, for example, four times higher risk at age 50 for individuals with two affected parents vs. those who have two non-affected parents. However, the results of models-fitting genetic analysis, did not confirm a monogenic Mendelian pattern of inheritance.

Modelling of age-related bone loss using cross-sectional data

BACKGROUND

Current applications of bone mineral density (BMD) data in age studies are not free of certain drawbacks. Since it is well established that age-related patterns of BMD changes involve three distinct periods (bone acquisition in youth, stabilization at maturity, and decrease with ageing), adjusting for age via an inappropriate mathematical function may lead to inconsistencies and wrong conclusions.

HYPOTHESIS

The piecewise model, which encompasses the above three periods, will accurately describe the BMD dependence on age.

OBJECTIVE

To examine age-related patterns of BMD changes using a number of possible mathematical functions and to find among them the best-fitting function. Next, to test whether the chosen function is universally applicable or if there are diverse population-specific functions.

MATERIAL AND METHODS

Thirteen ethnic samples from various regions of Europe and Asia, assigned into five ethnic-geographic groups, were examined. The total sample included 2430 males and 2515 females. Compact BMD of hand phalanges was measured by photodensitometry from plain radiographs of each individual studied. Statistical software was developed for the purposes of the present study; this software gave a maximum likelihood of the parameter estimates for various statistical models (functions).

RESULTS

In all samples of sufficient size and representative age range, a two-interval function was found as the best fitting and most parsimonious model to describe the BMD age-related changes. This two-interval function was characterized by age-related bone mass increase, positive slope beta(1s) in young age or a plateau (beta(1s) = 0, i.e. no age-related changes) until a sex-specific age threshold, T(0), after which annual bone loss ensued with a slope coefficient beta(2s). Threshold of BMD loss in women of different ethnic groups ranged between 37.85 and 47.77 years, and roughly coincided with perimenopausal age. In males, the age T(0) varied between 27.85 and 49.07 years. The ensuing cortical bone loss appeared to be linear in both sexes, averaging between 0.51% and 1.15% in men and between 0.74% and 1.77% per year of young age BMD value in females.

CONCLUSIONS

The change of phalangeal BMD with age may be best described by a two-interval function, regardless of sex and ethnic background. However, specific parameter estimates depend both on gender and ethnic affiliation. This study has yielded a well-fitted model of BMD dependence on age suitable for further use in population studies.