Resampling-based approaches to study variation in morphological modularity

Modularity has been suggested to be connected to evolvability because a higher degree of independence among parts allows them to evolve as separate units. Recently, the Escoufier RV coefficient has been proposed as a measure of the degree of integration between modules in multivariate morphometric datasets. However, it has been shown, using randomly simulated datasets, that the value of the RV coefficient depends on sample size. Also, so far there is no statistical test for the difference in the RV coefficient between a priori defined groups of observations. Here, we (1), using a rarefaction analysis, show that the value of the RV coefficient depends on sample size also in real geometric morphometric datasets; (2) propose a permutation procedure to test for the difference in the RV coefficient between a priori defined groups of observations; (3) show, through simulations, that such a permutation procedure has an appropriate Type I error; (4) suggest that a rarefaction procedure could be used to obtain sample-size-corrected values of the RV coefficient; and (5) propose a nearest-neighbor procedure that could be used when studying the variation of modularity in geographic space. The approaches outlined here, readily extendable to non-morphometric datasets, allow study of the variation in the degree of integration between a priori defined modules. A Java application – that will allow performance of the proposed test using a software with graphical user interface – has also been developed and is available at the Morphometrics at Stony Brook Web page (http://life.bio.sunysb.edu/morph/).

Geographical and morphological variation within and between colour phases in Coris julis (L. 1758), a protogynous marine fish

The possible differences between sexes in patterns of morphological variation in geographical space have been explored only in gonochorist freshwater species. We explored patterns of body shape variation in geographical space in a marine sequential hermaphrodite species, Coris julis (L. 1758), analyzing variation both within and between colour phases, through the use of geometric morphometrics and spatially-explicit statistical analyses. We also tested for the association of body shape with two environmental variables: temperature and chlorophyll a concentration, as obtained from time-series of satellite-derived data. Both colour phases showed a significant morphological variation in geographical space and patterns of variation divergent between phases. Although the morphological variation was qualitatively similar, individuals in the initial colour phase showed a more marked variation than individuals in the terminal phase. Body shape showed a weak but significant correlation with environmental variables, which was more pronounced in primary specimens.

Atlantic-Mediterranean and within-Mediterranean molecular variation in Coris julis (L. 1758) (Teleostei, Labridae)

Sequence variation in the mitochondrial control region was studied in the Mediterranean rainbow wrasse (Coris julis), a species with pronounced pelagic larval phase inhabiting the Mediterranean Sea and the adjacent coastal eastern Atlantic Ocean. A total of 309 specimens from 19 sampling sites were analysed with the aim of elucidating patterns of molecular variation between the Atlantic and the Mediterranean as well as within the Mediterranean Sea. Phylogeographic analyses revealed a pronounced structuring into a Mediterranean and an Atlantic group. Samples from a site at the Moroccan Mediterranean coast in the Alboran Sea showed intermediate frequencies of “Mediterranean” and “Atlantic” haplotypes. We recognised a departure from molecular neutrality and a star-like genealogy for samples from the Mediterranean Sea, which we propose to have happened due to a recent demographic expansion. The results are discussed in the light of previous studies on molecular variation in fish species between the Atlantic and the Mediterranean and within the Mediterranean.

Body shape variation and colour change during growth in a protogynous fish

Protogynous sequential hermaphroditism is very common in marine fish. Despite a large number of studies on various aspects of sequential hermaphroditism in fish, the relationship between body shape and colour during growth in dichromatic species has not been assessed. Using geometric morphometrics, the present study explores the relationship between growth, body shape and colouration in Coris julis (L. 1758), a small protogynous labrid species with distinct colour phases. Results show that body shape change during growth is independent of change in colour phase, a result which can be explained by the biology of the species and by the social control of sex change. Also, during growth the body grows deeper and the head has a steeper profile. It is hypothesized that a deeper body and a steeper profile might have a function in agonistic interactions between terminal phase males and that the marked chromatic difference between colour phases allows the lack of strict interdependence of body shape and colour during growth.

Traditional and geometric morphometrics detect morphological variation of lower pharyngeal jaw inCoris julis(Teleostei, Labridae)

In the present study, variation in the morphology of the lower pharyngeal element between two Sicilian populations of the rainbow wrasse Coris julis has been explored by the means of traditional morphometrics for size and geometric morphometrics for shape. Despite close geographical distance and probable high genetic flow between the populations, statistically significant differences have been found both for size and shape. In fact, one population shows a larger lower pharyngeal element that has a larger central tooth. Compared to the other population, this population also has medially enlarged lower pharyngeal jaws with a more pronounced convexity of the medial-posterior margin. The results are discussed in the light of a possible more pronounced durophagy of this population.

Common polygenic variation contributes to risk of migraine in the Norfolk Island population

Migraine has been defined as a common disabling primary headache disorder. Epidemiology studies have provided with the undeniable evidence of genetic components as active players in the development of the disease under a polygenic model in which multiple risk alleles exert modest individual effects. Our objective was to test the contribution of a polygenic effect to migraine risk in the Norfolk Island population using a panel of SNPs reported to be disease associated in published migraine GWAS. We also investigated whether individual SNPs were associated with gene expression levels measured in whole-blood. Polygenic scores were calculated in a total of 285 related individuals (74 cases, 211 controls) from the Norfolk Island using 51 SNPs previously reported to be associated with migraine in published GWAS. The association between polygenic score and migraine case-control status was tested using logistic regression. Results indicate that a migraine polygenic risk score was associated with migraine case-control status in this population (P=0.016). This supports the hypothesis that multiple SNPs with weak effects collectively contribute to migraine risk in this population. Amongst the SNPs included in the polygenic model, 4 were associated with the expression of the USMG5 gene, including rs171251 (P = 0.012). Results from this study provide evidence for a polygenic contribution to migraine risk in an isolated population and highlight specific SNPs that regulate the expression of USMG5, a gene critical for mitochondrial function.

Matters Arising: Zhu et al, "A novel gene variation of TNFα associated with ankylosing spondylitis: A reconfirmed study"

In a recently published study1 involving 79 ankylosing spondylitis (AS) patients and 132 unrelated healthy blood donors, Zhu et al report association of a tumour necrosis factor α (TNFα) single nucleotide polymorphism (SNP), rs1799724, with AS.1 I have many concerns with this paper and its conclusions...

Expectation and variation with a virtual die

By the time students reach the middle years they have experienced many chance activities based on dice. Common among these are rolling one die to explore the relationship of frequency and theoretical probability, and rolling two dice and summing the outcomes to consider their probabilities. Although dice may be considered overused by some, the advantage they offer is a familiar context within which to explore much more complex concepts. If the basic chance mechanism of the device is understood, it is possible to enter quickly into an arena of more complex concepts. This is what happened with a two hour activity engaged in by four classes of Grade 6 students in the same school. The activity targeted the concepts of variation and expectation. The teachers held extended discussions with their classes on variation and expectation at the beginning of the activity, with students contributing examples of the two concepts from their own experience. These notions are quite sophisticated for Grade 6, but the underlying concepts describe phenomena that students encounter every day. For example, time varies continuously; sporting results vary from game to game; the maximum temperature varies from day to day. However, there is an expectation about tomorrow’s maximum temperature based on the expert advice from the weather bureau. There may also be an expectation about a sporting result based on the participants’ previous results. It is this juxtaposition that makes life interesting. Variation hence describes the differences we see in phenomena around us. In a scenario displaying variation, expectation describes the effort to characterise or summarise the variation and perhaps make a prediction about the message arising from the scenario. The explicit purpose of the activity described here was to use the familiar scenario of rolling a die to expose these two concepts. Because the students had previously experienced rolling physical dice they knew instinctively about the variation that occurs across many rolls and about the theoretical expectation that each side should “come up” one-sixth of the time. They had observed the instances of the concepts in action, but had not consolidated the underlying terminology to describe it. As the two concepts are so fundamental to understanding statistics, we felt it would be useful to begin building in the familiar environment of rolling a die. Because hand-held dice limit the explorations students can undertake, the classes used the soft-ware TinkerPlots (Konold & Miller, 2011) to simulate rolling a die multiple times.

Regional and seasonal variation in airborne grass pollen levels between cities of Australia and New Zealand

Although grass pollen is widely regarded as the major outdoor aeroallergen source in Australia and New Zealand (NZ), no assemblage of airborne pollen data for the region has been previously compiled. Grass pollen count data collected at 14 urban sites in Australia and NZ over periods ranging from 1 to 17 years were acquired, assembled and compared, revealing considerable spatiotemporal variability. Although direct comparison between these data is problematic due to methodological differences between monitoring sites, the following patterns are apparent. Grass pollen seasons tended to have more than one peak from tropics to latitudes of 37°S and single peaks at sites south of this latitude. A longer grass pollen season was therefore found at sites below 37°S, driven by later seasonal end dates for grass growth and flowering. Daily pollen counts increased with latitude; subtropical regions had seasons of both high intensity and long duration. At higher latitude sites, the single springtime grass pollen peak is potentially due to a cooler growing season and a predominance of pollen from C3 grasses. The multiple peaks at lower latitude sites may be due to a warmer season and the predominance of pollen from C4 grasses. Prevalence and duration of seasonal allergies may reflect the differing pollen seasons across Australia and NZ. It must be emphasized that these findings are tentative due to limitations in the available data, reinforcing the need to implement standardized pollen-monitoring methods across Australasia. Furthermore, spatiotemporal differences in grass pollen counts indicate that local, current, standardized pollen monitoring would assist with the management of pollen allergen exposure for patients at risk of allergic rhinitis and asthma. © 2015 Springer Science+Business Media Dordrecht

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.

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.

Use of observed within-person variation of cardiac troponin in emergency department patients for determination of biological variation and percentage and absolute reference change values

BACKGROUND Many patients presenting to the emergency department (ED) for assessment of possible acute coronary syndrome (ACS) have low cardiac troponin concentrations that change very little on repeat blood draw. It is unclear if a lack of change in cardiac troponin concentration can be used to identify acutely presenting patients at low risk of ACS. METHODS We used the hs-cTnI assay from Abbott Diagnostics, which can detect cTnI in the blood of nearly all people. We identified a population of ED patients being assessed for ACS with repeat cTnI measurement who ultimately were proven to have no acute cardiac disease at the time of presentation. We used data from the repeat sampling to calculate total within-person CV (CV(T)) and, knowing the assay analytical CV (CV(A)), we could calculate within-person biological variation (CV(i)), reference change values (RCVs), and absolute RCV delta cTnI concentrations. RESULTS We had data sets on 283 patients. Men and women had similar CV(i) values of approximately 14%, which was similar at all concentrations <40 ng/L. The biological variation was not dependent on the time interval between sample collections (t = 1.5-17 h). The absolute delta critical reference change value was similar no matter what the initial cTnI concentration was. More than 90% of subjects had a critical reference change value <5 ng/L, and 97% had values of <10 ng/L. CONCLUSIONS With this hs-cTnI assay, delta cTnI seems to be a useful tool for rapidly identifying ED patients at low risk for possible ACS.

Global variation in the effects of ambient temperature on mortality: A systematic evaluation

Background: Studies have examined the effects of temperature on mortality in a single city, country, or region. However, less evidence is available on the variation in the associations between temperature and mortality in multiple countries, analyzed simultaneously. Methods: We obtained daily data on temperature and mortality in 306 communities from 12 countries/regions (Australia, Brazil, Thailand, China, Taiwan, Korea, Japan, Italy, Spain, United Kingdom, United States, and Canada). Two-stage analyses were used to assess the nonlinear and delayed relation between temperature and mortality. In the first stage, a Poisson regression allowing overdispersion with distributed lag nonlinear model was used to estimate the community-specific temperature-mortality relation. In the second stage, a multivariate meta-analysis was used to pool the nonlinear and delayed effects of ambient temperature at the national level, in each country. Results: The temperatures associated with the lowest mortality were around the 75th percentile of temperature in all the countries/regions, ranging from 66th (Taiwan) to 80th (UK) percentiles. The estimated effects of cold and hot temperatures on mortality varied by community and country. Meta-analysis results show that both cold and hot temperatures increased the risk of mortality in all the countries/regions. Cold effects were delayed and lasted for many days, whereas heat effects appeared quickly and did not last long. Conclusions: People have some ability to adapt to their local climate type, but both cold and hot temperatures are still associated with increased risk of mortality. Public health strategies to alleviate the impact of ambient temperatures are important, in particular in the context of climate change.

Temperature variation between neighboring days and mortality: a distributed lag non-linear analysis

Objectives To investigate whether a sudden temperature change between neighboring days has significant impact on mortality. Methods A Poisson generalized linear regression model combined with a distributed lag non-linear models was used to estimate the association of temperature change between neighboring days with mortality in a subtropical Chinese city during 2008–2012. Temperature change was calculated as the current day’s temperature minus the previous day’s temperature. Results A significant effect of temperature change between neighboring days on mortality was observed. Temperature increase was significantly associated with elevated mortality from non-accidental and cardiovascular diseases, while temperature decrease had a protective effect on non-accidental mortality and cardiovascular mortality. Males and people aged 65 years or older appeared to be more vulnerable to the impact of temperature change. Conclusions Temperature increase between neighboring days has a significant adverse impact on mortality. Further health mitigation strategies as a response to climate change should take into account temperature variation between neighboring days.

Polymorphisms in a desaturase 2 ortholog associate with cuticular hydrocarbon and male mating success variation in a natural population of Drosophila serrata

Elucidating the nature of genetic variation underlying both sexually selected traits and the fitness components of sexual selection is essential to understanding the broader consequences of sexual selection as an evolutionary process. To date, there have been relatively few attempts to connect the genetic variance in sexually selected traits with segregating DNA sequence polymorphisms. We set out to address this in a well-characterized sexual selection system - the cuticular hydrocarbons (CHCs) of Drosophila serrata - using an indirect association study design that allowed simultaneous estimation of the genetic variance in CHCs, sexual fitness and single nucleotide polymorphism (SNP) effects in an outbred population. We cloned and sequenced an ortholog of the D. melanogaster desaturase 2 gene, previously shown to affect CHC biosynthesis in D. melanogaster, and associated 36 SNPs with minor allele frequencies > 0.02 with variance in CHCs and sexual fitness. Three SNPs had significant multivariate associations with CHC phenotype (q-value < 0.05). At these loci, minor alleles had multivariate effects on CHCs that were weakly associated with the multivariate direction of sexual selection operating on these traits. Two of these SNPs had pleiotropic associations with male mating success, suggesting these variants may underlie responses to sexual selection due to this locus. There were 15 significant male mating success associations (q-value < 0.1), and interestingly, we detected a nonrandom pattern in the relationship between allele frequency and direction of effect on male mating success. The minor-frequency allele usually reduced male mating success, suggesting a positive association between male mating success and total fitness at this locus.