Novel mutations in ACVR1 result in atypical features in two fibrodysplasia ossificans progressiva patients

Fibrodysplasia Ossificans Progressiva (FOP) is a rare, heritable condition typified by progression of extensive ossification within skeletal muscle, ligament and tendon together with defects in skeletal development. The condition is easily diagnosed by the presence of shortened great toes and there is severe advancement of disability with age. FOP has been shown to result from a point mutation (c.617G>A) in the ACVR1 gene in almost all patients reported. Very recently two other mutations have been described in three FOP patients. We present here evidence for two further unique mutations (c.605G>T and c.983G>A) in this gene in two FOP patients with some atypical digit abnormalities and other clinical features. The observation of disparate missense mutations mapped to the GS and kinase domains of the protein supports the disease model of mild kinase activation and provides a potential rationale for phenotypic variation. © 2009 Petrie et al.

Susceptibility to ankylosing spondylitis in twins: The role of genes, HLA, and the environment

Objective To determine the relative effects of genetic and environmental factors in susceptibility to ankylosing spondylitis (AS). Methods Twins with AS were identified from the Royal National Hospital for Rheumatic Diseases database. Clinical and radiographic examinations were performed to establish diagnoses, and disease severity was assessed using a combination of validated scoring systems. HLA typing for HLA-B27, HLA-B60, and HLA-DR1 was performed by polymerase chain reaction with sequence- specific primers, and zygosity was assessed using microsatellite markers. Genetic and environmental variance components were assessed with the program Mx, using data from this and previous studies of twins with AS. Results Six of 8 monozygotic (MZ) twin pairs were disease concordant, compared with 4 of 15 B27-positive dizygotic (DZ) twin pairs (27%) and 4 of 32 DZ twin pairs overall (12.5%). Nonsignificant increases in similarity with regard to age at disease onset and all of the disease severity scores assessed were noted in disease-concordant MZ twins compared with concordant DZ twins. HLA-B27 and B60 were associated with the disease in probands, and the rate of disease concordance was significantly increased among DZ twin pairs in which the co- twin was positive for both B27 and DR1. Additive genetic effects were estimated to contribute 97% of the population variance. Conclusion Susceptibility to AS is largely genetically determined, and the environmental trigger for the disease is probably ubiquitous. HLA-B27 accounts for a minority of the overall genetic susceptibility to AS.

Genetic studies in osteoporosis: The end of the beginning

Osteoporosis and disorders of bone fragility are highly heritable, but despite much effort the identities of few of the genes involved has been established. Recent developments in genetics such as genome-wide association studies are revolutionizing research in this field, and it is likely that further contributions will be made through application of next-generation sequencing technologies, analysis of copy number variation polymorphisms, and high-throughput mouse mutagenesis programs. This article outlines what we know about osteoporosis genetics to date and the probable future directions of research in this field.

Mapping genes for osteoporosis-Old dogs and new tricks

In stark contrast to its horticultural origins, modern genetics is an extremely technology-driven field. Almost all the major advances in the field over the past 20 years have followed technological developments that have permitted change in study designs. The development of PCR in the 1980s led to RFLP mapping of monogenic diseases. The development of fluorescent-tagged genotyping methods led to linkage mapping approaches for common diseases that dominated the 1990s. The development of microarray SNP genotyping has led to the genome-wide association study era of the new millennium. And now the development of next-generation sequencing technologies is about to open up a new era of gene-mapping, enabling many potential new study designs. This review aims to present the strengths and weaknesses of the current approaches, and present some new ideas about gene-mapping approaches that are likely to advance our knowledge of the genes involved in heritable bone traits such as bone mineral density (BMD) and fracture.

Progress in the genetics of ankylosing spondylitis

Ankylosing spondylitis (AS) is a common, highly heritable, inflammatory arthropathy. In addition to being strongly associated with HLA-B27, a further 13 genes have been robustly associated with the disease. These genes highlight the involvement of the IL-23 pathway in disease pathogenesis, and indicate overlaps between the pathogenesis of AS, and of inflammatory bowel disease. Genetic associations in B27-positive and -negative disease are similar, with the main exception of association with ERAP1, which is restricted in association to B27-positive cases. This restriction, and the known function of ERAP1 in peptide trimming prior to HLA Class I presentation, indicates that HLA-B27 is likely to operate in AS by a mechanism involving aberrant peptide handling. These advances point to several potential novel therapeutic approaches in AS.

PTHR1 polymorphisms influence BMD variation through effects on the growing skeleton

We investigated whether polymorphisms in PTHR1 are associated with bone mineral density (BMD), to determine whether the association of this gene with BMD was due to effects on attainment of peak bone mass or effects on subsequent bone loss. The PTHR1 gene, including its 14 exons, their exon-intron boundaries, and 1,500 bp of its promoter region, was screened for polymorphisms by denaturing high-performance liquid chromatography (dHPLC) and sequencing in 36 osteoporotic cases. Eleven single-nucleotide polymorphisms (SNPs), one tetranucleotide repeat, and one tetranucleotide deletion were identified. A cohort of 634 families, including 1,236 men (39%) and 1,926 women (61%) ascertained with probands with low BMD (Z< -2.0) and the Children in Focus subset of the Avon Longitudinal Study of Parents and Children (ALSPAC) cohort (785 unrelated individuals, mean age 118 months), were genotyped for the five most informative SNPs (minor allele frequency >5%) and the tetranucleotide repeat. In our osteoporosis families, association was noted between lumbar spine BMD and alleles of a known functional tetranucleotide repeat (U4) in the PTHR1 promoter region (P = 0.042) and between two and three marker haplotypes of PTHR1 polymorphisms with lumbar spine, femoral neck, and total hip BMD (P = 0.021-0.047). This association was restricted to the youngest tertile of the population (age 16-39 years, P = 0.013-0.048). A similar association was found for the ALSPAC cohort: two marker haplotypes of SNPs A48609T and C52813T were associated with height (P = 0.006) and total body less head BMD (P = 0.02), corrected for age and gender, confirming the family findings. These findings suggest a role for PTHR1 variation in determining peak BMD.

Evidence of genetic association between TNFRSF1A encoding the p55 tumour necrosis factor receptor, and ankylosing spondylitis in UK Caucasians

Objectives: To replicate the possible genetic association between ankylosing spondylitis (AS) and TNFRSF1A. Methods: TNFRSF1A was re-sequenced in 48 individuals with AS to identify novel polymorphisms. Nine single nucleotide polymorphisms (SNPs) in TNFRSF1A and 5 SNPs in the neighbouring gene SCNN1A were genotyped in 1604 UK Caucasian individuals with AS and 1019 matched controls. An extended study was implemented using additional genotype data on 8 of these SNPs from 1400 historical controls from the 1958 British Birth Cohort. A meta-analysis of previously published results was also undertaken. Results: One novel variant in intron 6 was identified but no new coding variants. No definite associations were seen in the initial study but in the extended study there were weak associations with rs4149576 (p=0.04) and rs4149577 (p=0.007). In the metaanalysis consistent, somewhat stronger associations were seen with rs4149577 (p=0.002) and rs4149578 (p=0.006). Conclusions: These studies confirm the weak genetic associations between AS and TNFRSF1A. In view of the previously reported associations of TNFRSF1A with AS, in Caucasians and Chinese, and the biological plausibility of this candidate gene, replication of this finding in well powered studies is clearly indicated.

Cuff tear arthropathy: Evidence of functional variation in pyrophosphate metabolism genes

We investigated the role of two genes, ANKH and TNAP, in patients with cuff tear arthropathy. These genes encode proteins which regulate the extracellular concentration of inorganic pyrophosphate, fluctuations of which can lead to calcium crystal formation. Variants were detected by direct sequencing of DNA and their frequencies compared with healthy controls. The effect of variants on protein function was further studied by in vitro approaches. Variant genotypes were observed more frequently in the cases when compared with controls in ANKH (45% and 20%) and TNAP (32% and 9%). Variants in ANKH altered inorganic pyrophosphate (PPi) concentrations in transfected human chondrocytes. There was a higher mean serum concentration of TNAP detected in female patients compared with normal ranges. Cuff tear arthropathy is associated with variants in ANKH and TNAP that alter extracellular inorganic pyrophosphate concentrations causing calcium crystal deposition. This supports a theory that genetic variants predispose patients to primary crystal deposition which when combined with a massive rotator cuff tear leads to the development of arthritis.

Quantification of stability in an agility drill using linear and nonlinear measures of variability

This study implemented linear and nonlinear methods of measuring variability to determine differences in stability of two groups of skilled (n = 10) and unskilled (n = 10) participants performing 3m forward/backward shuttle agility drill. We also determined whether stability measures differed between the forward and backward segments of the drill. Finally, we sought to investigate whether local dynamic stability, measured using largest finite-time Lyapunov exponents, changed from distal to proximal lower extremity segments. Three-dimensional coordinates of five lower extremity markers data were recorded. Results revealed that the Lyapunov exponents were lower (P < 0.05) for skilled participants at all joint markers indicative of higher levels of local dynamic stability. Additionally, stability of motion did not differ between forward and backward segments of the drill (P > 0.05), signifying that almost the same control strategy was used in forward and backward directions by all participants, regardless of skill level. Furthermore, local dynamic stability increased from distal to proximal joints (P < 0.05) indicating that stability of proximal segments are prioritized by the neuromuscular control system. Finally, skilled participants displayed greater foot placement standard deviation values (P < 0.05), indicative of adaptation to task constraints. The results of this study provide new methods for sport scientists, coaches to characterize stability in agility drill performance.

Increasing functional variability in the preparatory phase of the takeoff improves elite springboard diving performance

PURPOSE: Previous research demonstrating that specific performance outcome goals can be achieved in different ways is functionally significant for springboard divers whose performance environment can vary extensively. This body of work raises questions about the traditional approach of balking (terminating the takeoff) by elite divers aiming to perform only identical, invariant movement patterns during practice. METHOD: A 12-week training program (2 times per day; 6.5 hr per day) was implemented with 4 elite female springboard divers to encourage them to adapt movement patterns under variable takeoff conditions and complete intended dives, rather than balk. RESULTS: Intraindividual analyses revealed small increases in variability in the board-work component of each diver's pretraining and posttraining program reverse-dive takeoffs. No topological differences were observed between movement patterns of dives completed pretraining and posttraining. Differences were noted in the amount of movement variability under different training conditions (evidenced by higher normalized root mean square error indexes posttraining). An increase in the number of completed dives (from 78.91%-86.84% to 95.59%-99.29%) and a decrease in the frequency of balked takeoffs (from 13.16%-19.41% to 0.63%-4.41%) showed that the elite athletes were able to adapt their behaviors during the training program. These findings coincided with greater consistency in the divers' performance during practice as scored by qualified judges. CONCLUSION: Results suggested that on completion of training, athletes were capable of successfully adapting their movement patterns under more varied takeoff conditions to achieve greater consistency and stability of performance outcomes.

Population of computational rabbit-specific ventricular action potential models for investigating sources of variability in cellular repolarisation

Variability is observed at all levels of cardiac electrophysiology. Yet, the underlying causes and importance of this variability are generally unknown, and difficult to investigate with current experimental techniques. The aim of the present study was to generate populations of computational ventricular action potential models that reproduce experimentally observed intercellular variability of repolarisation (represented by action potential duration) and to identify its potential causes. A systematic exploration of the effects of simultaneously varying the magnitude of six transmembrane current conductances (transient outward, rapid and slow delayed rectifier K(+), inward rectifying K(+), L-type Ca(2+), and Na(+)/K(+) pump currents) in two rabbit-specific ventricular action potential models (Shannon et al. and Mahajan et al.) at multiple cycle lengths (400, 600, 1,000 ms) was performed. This was accomplished with distributed computing software specialised for multi-dimensional parameter sweeps and grid execution. An initial population of 15,625 parameter sets was generated for both models at each cycle length. Action potential durations of these populations were compared to experimentally derived ranges for rabbit ventricular myocytes. 1,352 parameter sets for the Shannon model and 779 parameter sets for the Mahajan model yielded action potential duration within the experimental range, demonstrating that a wide array of ionic conductance values can be used to simulate a physiological rabbit ventricular action potential. Furthermore, by using clutter-based dimension reordering, a technique that allows visualisation of multi-dimensional spaces in two dimensions, the interaction of current conductances and their relative importance to the ventricular action potential at different cycle lengths were revealed. Overall, this work represents an important step towards a better understanding of the role that variability in current conductances may play in experimentally observed intercellular variability of rabbit ventricular action potential repolarisation.

Mixture models of nucleotide sequence evolution that account for heterogeneity in the substitution process across sites and across lineages

Molecular phylogenetic studies of homologous sequences of nucleotides often assume that the underlying evolutionary process was globally stationary, reversible, and homogeneous (SRH), and that a model of evolution with one or more site-specific and time-reversible rate matrices (e.g., the GTR rate matrix) is enough to accurately model the evolution of data over the whole tree. However, an increasing body of data suggests that evolution under these conditions is an exception, rather than the norm. To address this issue, several non-SRH models of molecular evolution have been proposed, but they either ignore heterogeneity in the substitution process across sites (HAS) or assume it can be modeled accurately using the distribution. As an alternative to these models of evolution, we introduce a family of mixture models that approximate HAS without the assumption of an underlying predefined statistical distribution. This family of mixture models is combined with non-SRH models of evolution that account for heterogeneity in the substitution process across lineages (HAL). We also present two algorithms for searching model space and identifying an optimal model of evolution that is less likely to over- or underparameterize the data. The performance of the two new algorithms was evaluated using alignments of nucleotides with 10 000 sites simulated under complex non-SRH conditions on a 25-tipped tree. The algorithms were found to be very successful, identifying the correct HAL model with a 75% success rate (the average success rate for assigning rate matrices to the tree's 48 edges was 99.25%) and, for the correct HAL model, identifying the correct HAS model with a 98% success rate. Finally, parameter estimates obtained under the correct HAL-HAS model were found to be accurate and precise. The merits of our new algorithms were illustrated with an analysis of 42 337 second codon sites extracted from a concatenation of 106 alignments of orthologous genes encoded by the nuclear genomes of Saccharomyces cerevisiae, S. paradoxus, S. mikatae, S. kudriavzevii, S. castellii, S. kluyveri, S. bayanus, and Candida albicans. Our results show that second codon sites in the ancestral genome of these species contained 49.1% invariable sites, 39.6% variable sites belonging to one rate category (V1), and 11.3% variable sites belonging to a second rate category (V2). The ancestral nucleotide content was found to differ markedly across these three sets of sites, and the evolutionary processes operating at the variable sites were found to be non-SRH and best modeled by a combination of eight edge-specific rate matrices (four for V1 and four for V2). The number of substitutions per site at the variable sites also differed markedly, with sites belonging to V1 evolving slower than those belonging to V2 along the lineages separating the seven species of Saccharomyces. Finally, sites belonging to V1 appeared to have ceased evolving along the lineages separating S. cerevisiae, S. paradoxus, S. mikatae, S. kudriavzevii, and S. bayanus, implying that they might have become so selectively constrained that they could be considered invariable sites in these species.

The koala immunological toolkit: Sequence identification and comparison of key markers of the koala (Phascolarctos cinereus) immune response

The koala (Phascolarctos cinereus) is an Australian marsupial that continues to experience significant population declines. Infectious diseases caused by pathogens such as Chlamydia are proposed to have a major role. Very few species-specific immunological reagents are available, severely hindering our ability to respond to the threat of infectious diseases in the koala. In this study, we utilise data from the sequencing of the koala transcriptome to identify key immunological markers of the koala adaptive immune response and cytokines known to be important in the host response to chlamydial infection in other species. This report describes the identification and preliminary sequence analysis of (1) T lymphocyte glycoprotein markers (CD4, CD8); (2) IL-4, a marker for the Th2 response; (3) cytokines such as IL-6, IL-12 and IL-1β, that have been shown to have a role in chlamydial clearance and pathology in other hosts; and (4) the sequences for the koala immunoglobulins, IgA, IgG, IgE and IgM. These sequences will enable the development of a range of immunological reagents for understanding the koala’s innate and adaptive immune responses, while also providing a resource that will enable continued investigations into the origin and evolution of the marsupial immune system.

Evaluation of the relationship between Chlamydia pecorum sequence types and disease using a species-specific multi-locus sequence typing scheme (MLST)

Chlamydia pecorum is globally associated with several ovine diseases including keratoconjunctivitis and polyarthritis. The exact relationship between the variety of C. pecorum strains reported and the diseases described in sheep remains unclear, challenging efforts to accurately diagnose and manage infected flocks. In the present study, we applied C. pecorum multi-locus sequence typing (MLST) to C. pecorum positive samples collected from sympatric flocks of Australian sheep presenting with conjunctivitis, conjunctivitis with polyarthritis, or polyarthritis only and with no clinical disease (NCD) in order to elucidate the exact relationships between the infecting strains and the range of diseases. Using Bayesian phylogenetic and cluster analyses on 62 C. pecorum positive ocular, vaginal and rectal swab samples from sheep presenting with a range of diseases and in a comparison to C. pecorum sequence types (STs) from other hosts, one ST (ST 23) was recognised as a globally distributed strain associated with ovine and bovine diseases such as polyarthritis and encephalomyelitis. A second ST (ST 69) presently only described in Australian animals, was detected in association with ovine as well as koala chlamydial infections. The majority of vaginal and rectal C. pecorum STs from animals with NCD and/or anatomical sites with no clinical signs of disease in diseased animals, clustered together in a separate group, by both analyses. Furthermore, 8/13 detected STs were novel. This study provides a platform for strain selection for further research into the pathogenic potential of C. pecorum in animals and highlights targets for potential strain-specific diagnostic test development.

The revolution in human monogenic disease mapping

The successful completion of the Human Genome Project (HGP) was an unprecedented scientific advance that has become an invaluable resource in the search for genes that cause monogenic and common (polygenic) diseases. Prior to the HGP, linkage analysis had successfully mapped many disease genes for monogenic disorders; however, the limitations of this approach were particularly evident for identifying causative genes in rare genetic disorders affecting lifespan and/or reproductive fitness, such as skeletal dysplasias. In this review, we illustrate the challenges of mapping disease genes in such conditions through the ultra-rare disorder fibrodysplasia ossificans progressiva (FOP) and we discuss the advances that are being made through current massively parallel (“next generation”) sequencing (MPS) technologies.