The validity and heritability of self-report osteoarthritis in an Australian older twin sample

In order to investigate the genetic and environmental antecedents of osteoarthritis (CA), self-report measures of joint pain, stiffness and swelling were obtained from a population-based sample of 1242 twin pairs over 50 years of age. In order to provide validation for these self-report measures, a subsample of 118 twin pairs were examined according to the American College of Rheumatology clinical and radiographic criteria for the classification of osteoarthritis. A variety of statistical methods were employed to identify the model derived from self-report variables which would provide optimal prediction of these standardised assessments, and structural equation modelling was used to determine the relative influences of genetic and environmental influences on the development of osteoarthritis. Significant genetic effects were found to contribute to osteoarthritis of the hands, hips and knees in women, with heritability estimates ranging from 30-46% depending on the site. In addition, the additive genetic effects contributing to osteoarthritis in various parts of the body were confirmed to be the same. Statistically significant familial aggregation of osteoarthritis in men was also observed, but it was not possible to determine whether this was due to genetic or shared environmental effects.

Genetics of Serum Dehydroepiandrosterone Sulfate and Its Relationship to Insulin in a Population-Based Cohort of Twin Subjects

Previous studies have shown a significant effect of insulin administration on serum dehydroepiandrosterone sulfate (DHEA-S) concentration and its metabolic rate, with evidence for the effect in men, but not in women. This could lead to differences in the sources of variation in serum DHEA-S between men and women and in its covariation with insulin concentration. This study aimed to test whether these hypotheses were supported in a sample of healthy adult twins. Serum DHEA-S (n=2287) and plasma insulin (n=2436) were measured in samples from adult male and female twins recruited through the Australian Twin Registry. Models of genetic and environmental sources of variation and covariation were tested against the data. DHEA-S showed substantial genetic effects in both men and women after adjustment for covariates, including sex, age, body mass index, and time since the last meal. There was no significant phenotypic or genetic correlation between DHEA-S and insulin in either men or women. Despite the experimental evidence for insulin infusion producing a reduction in serum DHEA-S and some effect of meals on the observed DHEA-S concentration, there were no associations between insulin and DHEA-S at the population level. Variations in DHEA-S are due to age, sex, obesity, and substantial polygenic genetic influences.

Geographical variation in the penetrance of CDKN2A mutations for melanoma

Background: Germline mutations in the CDKN2A gene, which encodes two proteins (p16INK4A and p14ARF), are the most common cause of inherited susceptibility to melanoma. We examined the penetrance of such mutations using data from eight groups from Europe, Australia and the United States that are part of The Melanoma Genetics Consortium Methods: We analyzed 80 families with documented CDKN2A mutations and multiple cases of cutaneous melanoma. We modeled penetrance for melanoma using a logistic regression model incorporating survival analysis. Hypothesis testing was based on likelihood ratio tests. Covariates included gender, alterations in p14APF protein, and population melanoma incidence rates. All statistical tests were two-sided. Results: The 80 analyzed families contained 402 melanoma patients, 320 of whom were tested for mutations and 291 were mutation carriers. We also tested 713 unaffected family members for mutations and 194 were carriers. Overall, CDKN2A mutation penetrance was estimated to be 0.30 (95% confidence interval (CI) = 0.12 to 0.62) by age 50 years and 0.67 (95% CI = 0.31 to 0.96) by age 80 years. Penetrance was not statistically significantly modified by gender or by whether the CDKN2A mutation altered p14ARF protein. However, there was a statistically significant effect of residing in a location with a high population incidence rate of melanoma (P = .003). By age 50 years CDKN2A mutation penetrance reached 0.13 in Europe, 0.50 in the United States, and 0.32 in Australia; by age 80 years it was 0.58 in Europe, 0.76 in the United States, and 0.91 in Australia. Conclusions: This study, which gives the most informed estimates of CDKN2A mutation penetrance available, indicates that the penetrance varies with melanoma population incidence rates. Thus, the same factors that affect population incidence of melanoma may also mediate CDKN2A penetrance.

Interspecific variation in plumage colour among birds: species recognition or light environment?

The traditional explanation for interspecific plumage colour variation in birds is that colour differences between species are adaptations to minimize the risk of hybridization. Under this explanation, colour differences between closely related species of birds represent reproductive character displacement. An alternative explanation is that interspecific variation in plumage colour is an adaptive response to variation in light environments across habitats. Under this explanation, differences in colour between closely related species are a product of selection on signal efficiency. We use a comparative approach to examine these two hypotheses, testing the effects of sympatry and habitat use, respectively, on divergence in male plumage colour. Contrary to the prediction of the Species Isolation Hypothesis, we find no evidence that sympatric pairs of species are consistently more divergent in coloration than are allopatric pairs of species. However, in agreement with the Light Environment Hypothesis, we find significant associations between plumage coloration and habitat use. All of these results remain qualitatively unchanged irrespective of the statistical methodology used to compare reflectance spectra, the body regions used in the analyses, or the exclusion of areas of plumage not used in sexual displays. Our results suggest that, in general, interspecific variation in plumage colour among birds is more strongly influenced by the signalling environment than by the risk of hybridization.

Categorization and characterization of transcript-confirmed constitutively and alternatively spliced introns and exons from human

By spliced alignment of human DNA and transcript sequence data we constructed a data set of transcript-confirmed exons and introns from 2793 genes, 796 of which (28%) were seen to have multiple isoforms. We find that over one-third of human exons can translate in more than one frame, and that this is highly correlated with G+C content. Introns containing adenosine at donor site position +3 (A3), rather than guanosine (G3), are more common in low G+C regions, while the converse is true in high G+C regions. These two classes of introns are shown to have distinct lengths, consensus sequences and correlations among splice signals, leading to the hypothesis that A3 donor sites are associated with exon definition, and G3 donor sites with intron definition. Minor classes of introns, including GC-AG, U12-type GT-AG, weak, and putative AG-dependant introns are identified and characterized. Cassette exons are more prevalent in low G+C regions, while exon isoforms are more prevalent in high G+C regions. Cassette exon events outnumber other alternative events, while exon isoform events involve truncation twice as often as extension, and occur at acceptor sites twice as often as at donor sites. Alternative splicing is usually associated with weak splice signals, and in a majority of cases, preserves the coding frame. The reported characteristics of constitutive and alternative splice signals, and the hypotheses offered regarding alternative splicing and genome organization, have important implications for experimental research into RNA processing. The 'AltExtron' data sets are available at http://www.bit.uq.edu.au/altExtron/ and http://www.ebi.ac.uk/similar tothanaraj/altExtron/.

Bridging the gap between different statistical approaches: An integrated framework for modelling

This paper proposes a template for modelling complex datasets that integrates traditional statistical modelling approaches with more recent advances in statistics and modelling through an exploratory framework. Our approach builds on the well-known and long standing traditional idea of 'good practice in statistics' by establishing a comprehensive framework for modelling that focuses on exploration, prediction, interpretation and reliability assessment, a relatively new idea that allows individual assessment of predictions. The integrated framework we present comprises two stages. The first involves the use of exploratory methods to help visually understand the data and identify a parsimonious set of explanatory variables. The second encompasses a two step modelling process, where the use of non-parametric methods such as decision trees and generalized additive models are promoted to identify important variables and their modelling relationship with the response before a final predictive model is considered. We focus on fitting the predictive model using parametric, non-parametric and Bayesian approaches. This paper is motivated by a medical problem where interest focuses on developing a risk stratification system for morbidity of 1,710 cardiac patients given a suite of demographic, clinical and preoperative variables. Although the methods we use are applied specifically to this case study, these methods can be applied across any field, irrespective of the type of response.

Oceanic interchange and nonequilibrium population structure in the estuarine dependent Indo-Pacific tasselfish, Polynemus sheridani

We assayed mtDNA haplotype [300 base pairs (bp) control region] geography and genealogy in the Indo-Pacific tasselfish, Polynemus sheridani from its contiguous estuarine distribution across northern Australia (n = 169). Eight estuaries were sampled from three oceanographic regions (Timor Sea, Gulf of Carpentaria and the Coral Sea) to assess the impact of Pleistocene sea level changes on the historical connectivity among P. sheridani populations. Specifically, we investigated the genetic consequences of disruption to Indian-Pacific Ocean connectivity brought about by the closure of the Torres Strait. Overall there was significant population subdivision among estuaries (F-ST = 0.161, (Phi(ST) = 0.187). Despite a linear distribution, P. sheridani did not show isolation by distance over the entire sampled range because of genetic similarity of estuaries greater than 3000 km apart. However, significant isolation by distance was detected between estuaries separated by less than 3000 km of coastline. Unlike many genetic studies of Indo-Pacific marine species, there was no evidence for an historical division between eastern and western populations. Instead, phylogeographical patterns were dominated by a starlike intraspecific phylogeny coupled with evidence for population expansion in both the Gulf of Carpentaria and the Coral Sea but not the Timor Sea. This was interpreted as evidence for recent west to east recolonization across of northern Australia following the last postglacial marine advance. We argue that although sufficient time has elapsed postcolonization for populations to approach gene flow-drift equilibrium over smaller spatial scales (< 3000 km), the signal of historical colonization persists to obscure the expected equilibrium pattern of isolation by distance over large spatial scales (> 3000 km).

Effect of changing data collection parameters on statistical motor unit number estimates

The effect of number of samples and selection of data for analysis on the calculation of surface motor unit potential (SMUP) size in the statistical method of motor unit number estimates (MUNE) was determined in 10 normal subjects and 10 with amyotrophic lateral sclerosis (ALS). We recorded 500 sequential compound muscle action potentials (CMAPs) at three different stable stimulus intensities (10–50% of maximal CMAP). Estimated mean SMUP sizes were calculated using Poisson statistical assumptions from the variance of 500 sequential CMAP obtained at each stimulus intensity. The results with the 500 data points were compared with smaller subsets from the same data set. The results using a range of 50–80% of the 500 data points were compared with the full 500. The effect of restricting analysis to data between 5–20% of the CMAP and to standard deviation limits was also assessed. No differences in mean SMUP size were found with stimulus intensity or use of different ranges of data. Consistency was improved with a greater sample number. Data within 5% of CMAP size gave both increased consistency and reduced mean SMUP size in many subjects, but excluded valid responses present at that stimulus intensity. These changes were more prominent in ALS patients in whom the presence of isolated SMUP responses was a striking difference from normal subjects. Noise, spurious data, and large SMUP limited the Poisson assumptions. When these factors are considered, consistent statistical MUNE can be calculated from a continuous sequence of data points. A 2 to 2.5 SD or 10% window are reasonable methods of limiting data for analysis. Muscle Nerve 27: 320–331, 2003

The role of melanocortin-1 receptor polymorphism in skin cancer risk phenotypes

We have examined melanocortin-1 receptor (MC1R) variant allele frequencies in the general population and in a collection of adolescent dizygotic and monozygotic twins to determine statistical associations of pigmentation phenotypes with increased skin cancer risk. This included hair and skin color, freckling, mole count and sun exposed skin reflectance. Nine variants were studied and designated as either strong R (OR = 63; 95% CI 32-140) or weak r (OR = 5; 95% CI 3-11) red hair alleles. Penetrance of each MC1R variant allele was consistent with an allelic model where effects were multiplicative for red hair but additive for skin reflectance. To assess the interaction of the brown eye color gene BEY2/OCA2 on the phenotypic effects of variant MC1R alleles we imputed OCA2 genotype in the twin collection. A modifying effect of OCA2 on MC1R variant alleles was seen on constitutive skin color, freckling and mole count. In order to study the individual effects of these variants on pigmentation phenotype we have established a series of human primary melanocyte strains genotyped for the MC1R receptor. These include strains which are MC1R wild-type consensus, variant heterozygotes, and homozygotes for strong R alleles Arg151Cys and Arg160Trp. Ultrastructural analysis demonstrated that only consensus strains contained stage III and IV melanosomes in their terminal dendrites whereas Arg151Cys and Arg160Trp homozygous strains contained only immature stage I and II melanosomes. Such genetic association studies combined with the functional analysis of MC1R variant alleles in melanocytic cells should provide a link in understanding the association between pigmentary phototypes and skin cancer risk.

Genetics of melanoma predisposition

Predisposition to melanoma is genetically heterogeneous. Two high penetrance susceptibility genes, CDKN2A and CDK4, have so far been identified and mapping is ongoing to localize and identify others. With the advent of a catalogue of millions of potential DNA polymorphisms, attention is now also being focused on identification of genes that confer a more modest contribution to melanoma risk, such as those encoding proteins involved in pigmentation, DNA repair, cell growth and differentiation or detoxification of metabolites. One such pigmentation gene, MC1R, has not only been found to be a low penetrance melanoma gene but has also been shown to act as a genetic modifier of melanoma risk in individuals carrying CDKN2A mutations. Most recently, an environmental agent, ultraviolet radiation, has also been established as a modifier of melanoma risk in CDKN2A mutation carriers. Hence, melanoma is turning out to be an excellent paradigm for studying gene-gene and gene-environment interactions.

The genetics of coronary heart disease: the contribution of twin studies

Despite the decline in coronary heart disease in many European countries, the disease remains an enormous public health problem. Although we know a great deal about environmental risk factors for coronary heart disease, a heritable component was recognized a long time ago. The earliest and best known examples of how our genetic constitution may determine cardiovascular risk relate to lipoprotein(a), familial hypercholesterolaemia and apolipoprotein E. In the past 20 years a fair number of polymorphisms assessed singly have shown strong associations with the disease but most are subject to poor repeatability. Twins constitute a compelling natural experiment to establish the genetic contribution to coronary heart disease and its risk factors. GenomEUtwin, a recently funded Framework 5 Programme of the European Community, affords the opportunity of comparing the heritability of risk factors in different European Twin Registries. As an illustration we present the heritabilities of systolic and diastolic blood pressure, based on data from over 4000 twin pairs from six different European countries and Australia. Heritabilities for systolic blood pressure are between 52 and 66% and for diastolic blood pressure between 44 and 66%. There is no evidence of sex differences in heritability estimates and very little to no evidence for a significant contribution of shared family environment. A non-twin based prospective case/cohort study of coronary heart disease and stroke (MORGAM) will allow hypotheses relating to cardiovascular disease, generated in the twin cohorts, to be tested prospectively in adult populations. Twin studies have also contributed to our understanding of the life course hypothesis, and GenomEUtwin has the potential to add to this.

Heritability of adult body height: A comparative study of twin cohorts in eight countries

A major component of variation in body height is due to genetic differences, but environmental factors have a substantial contributory effect. In this study we aimed to analyse whether the genetic architecture of body height varies between affluent western societies. We analysed twin data from eight countries comprising 30,111 complete twin pairs by using the univariate genetic model of the Mx statistical package. Body height and zygosity were self-reported in seven populations and measured directly in one population. We found that there was substantial variation in mean body height between countries; body height was least in Italy (177 cm in men and 163 cm in women) and greatest in the Netherlands (184 cm and 171 cm, respectively). In men there was no corresponding variation in heritability of body height, heritability estimates ranging from 0.87 to 0.93 in populations under an additive genes/unique environment (AE) model. Among women the heritability estimates were generally lower than among men with greater variation between countries, ranging from 0.68 to 0.84 when an additive genes/shared environment/unique environment (ACE) model was used. In four populations where an AE model fit equally well or better, heritability ranged from 0.89 to 0.93. This difference between the sexes was mainly due to the effect of the shared environmental component of variance, which appears to be more important among women than among men in our study populations. Our results indicate that, in general, there are only minor differences in the genetic architecture of height between affluent Caucasian populations, especially among men.