Identification of candidate gene regions in the rat by co-localization of QTLs for bone density, size, structure and strength

The genetics of bone mass and susceptibility to bone diseases

Osteoporosis is characterized by low bone mass and an increased risk of fracture. Genetic factors, environmental factors and gene-environment interactions all contribute to a person's lifetime risk of developing an osteoporotic fracture. This Review summarizes key advances in understanding of the genetics of bone traits and their role in osteoporosis. Candidate-gene approaches dominated this field 20 years ago, but clinical and preclinical genetic studies published in the past 5 years generally utilize more-sophisticated and better-powered genome-wide association studies (GWAS). High-throughput DNA sequencing, large genomic databases and improved methods of data analysis have greatly accelerated the gene-discovery process. Linkage analyses of single-gene traits that segregate in families with extreme phenotypes have led to the elucidation of critical pathways controlling bone mass. For example, components of the Wnt-β-catenin signalling pathway have been validated (in both GWAS and functional studies) as contributing to various bone phenotypes. These notable advances in gene discovery suggest that the next decade will witness cataloguing of the hundreds of genes that influence bone mass and osteoporosis, which in turn will provide a roadmap for the development of new drugs that target diseases of low bone mass, including osteoporosis.

Genetic regulation of bone strength: a review of animal model studies

Population- and family-based studies have established that fragility fracture risk is heritable; yet, the genome-wide association studies published to date have only accounted for a small fraction of the known variation for fracture risk of either the femur or the lumbar spine. Much work has been carried out using animal models toward finding genetic loci that are associated with bone strength. Studies using animal models overcome some of the issues associated with using patient data, but caution is needed when interpreting the results. In this review, we examine the types of tests that have been used for forward genetics mapping in animal models to identify loci and/or genes that regulate bone strength and discuss the limitations of these test methods. In addition, we present a summary of the quantitative trait loci that have been mapped for bone strength in mice, rats and chickens. The majority of these loci co-map with loci for bone size and/or geometry and thus likely dictate strength via modulating bone size. Differences in bone matrix composition have been demonstrated when comparing inbred strains of mice, and these matrix differences may be associated with differences in bone strength. However, additional work is needed to identify loci that act on bone strength at the materials level.

Widespread differential maternal and paternal genome effects on fetal bone phenotype at mid-gestation

Parent-of-origin-dependent (epi)genetic factors are important determinants of prenatal development that program adult phenotype. However, data on magnitude and specificity of maternal and paternal genome effects on fetal bone are lacking. We used an outbred bovine model to dissect and quantify effects of parental genomes, fetal sex, and nongenetic maternal effects on the fetal skeleton and analyzed phenotypic and molecular relationships between fetal muscle and bone. Analysis of 51 bone morphometric and weight parameters from 72 fetuses recovered at day 153 gestation (54% term) identified six principal components (PC1-6) that explained 80% of the variation in skeletal parameters. Parental genomes accounted for most of the variation in bone wet weight (PC1, 72.1%), limb ossification (PC2, 99.8%), flat bone size (PC4, 99.7%), and axial skeletal growth (PC5, 96.9%). Limb length showed lesser effects of parental genomes (PC3, 40.8%) and a significant nongenetic maternal effect (gestational weight gain, 29%). Fetal sex affected bone wet weight (PC1, p < 0.0001) and limb length (PC3, p < 0.05). Partitioning of variation explained by parental genomes revealed strong maternal genome effects on bone wet weight (74.1%, p < 0.0001) and axial skeletal growth (93.5%, p < 0.001), whereas paternal genome controlled limb ossification (95.1%, p < 0.0001). Histomorphometric data revealed strong maternal genome effects on growth plate height (98.6%, p < 0.0001) and trabecular thickness (85.5%, p < 0.0001) in distal femur. Parental genome effects on fetal bone were mirrored by maternal genome effects on fetal serum 25-hydroxyvitamin D (96.9%, p < 0.001) and paternal genome effects on alkaline phosphatase (90.0%, p < 0.001) and their correlations with maternally controlled bone wet weight and paternally controlled limb ossification, respectively. Bone wet weight and flat bone size correlated positively with muscle weight (r = 0.84 and 0.77, p < 0.0001) and negatively with muscle H19 expression (r = -0.34 and -0.31, p < 0.01). Because imprinted maternally expressed H19 regulates growth factors by miRNA interference, this suggests muscle-bone interaction via epigenetic factors.

A genome-wide association study demonstrates significant genetic variation for fracture risk in Thoroughbred racehorses

Background

Thoroughbred racehorses are subject to non-traumatic distal limb bone fractures that occur during racing and exercise. Susceptibility to fracture may be due to underlying disturbances in bone metabolism which have a genetic cause. Fracture risk has been shown to be heritable in several species but this study is the first genetic analysis of fracture risk in the horse.

Results

Fracture cases (n = 269) were horses that sustained catastrophic distal limb fractures while racing on UK racecourses, necessitating euthanasia. Control horses (n = 253) were over 4 years of age, were racing during the same time period as the cases, and had no history of fracture at the time the study was carried out. The horses sampled were bred for both flat and National Hunt (NH) jump racing. 43,417 SNPs were employed to perform a genome-wide association analysis and to estimate the proportion of genetic variance attributable to the SNPs on each chromosome using restricted maximum likelihood (REML). Significant genetic variation associated with fracture risk was found on chromosomes 9, 18, 22 and 31. Three SNPs on chromosome 18 (62.05 Mb – 62.15 Mb) and one SNP on chromosome 1 (14.17 Mb) reached genome-wide significance (p < 0.05) in a genome-wide association study (GWAS). Two of the SNPs on ECA 18 were located in a haplotype block containing the gene zinc finger protein 804A (ZNF804A). One haplotype within this block has a protective effect (controls at 1.95 times less risk of fracture than cases, p = 1 × 10^-4), while a second haplotype increases fracture risk (cases at 3.39 times higher risk of fracture than controls, p = 0.042).

Conclusions

Fracture risk in the Thoroughbred horse is a complex condition with an underlying genetic basis. Multiple genomic regions contribute to susceptibility to fracture risk. This suggests there is the potential to develop SNP-based estimators for genetic risk of fracture in the Thoroughbred racehorse, using methods pioneered in livestock genetics such as genomic selection. This information would be useful to racehorse breeders and owners, enabling them to reduce the risk of injury in their horses.