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/).

The gut microbial community of midas cichlid fish in repeatedly evolved limnetic-benthic species pairs

Gut bacterial communities are now known to influence a range of fitness related aspects of organisms. But how different the microbial community is in closely related species, and if these differences can be interpreted as adaptive is still unclear. In this study we compared microbial communities in two sets of closely related sympatric crater lake cichlid fish species pairs that show similar adaptations along the limnetic-benthic axis. The gut microbial community composition differs in the species pair inhabiting the older of two crater lakes. One major difference, relative to other fish, is that in these cichlids that live in hypersaline crater lakes, the microbial community is largely made up of Oceanospirillales (52.28%) which are halotolerant or halophilic bacteria. This analysis opens up further avenues to identify candidate symbiotic or co-evolved bacteria playing a role in adaptation to similar diets and life-styles or even have a role in speciation. Future functional and phylosymbiotic analyses might help to address these issues.

Bat species comparisons based on external morphology: A test of traditional versus geometric morphometric approaches

External morphology is commonly used to identify bats as well as to investigate flight and foraging behavior, typically relying on simple length and area measures or ratios. However, geometric morphometrics is increasingly used in the biological sciences to analyse variation in shape and discriminate among species and populations. Here we compare the ability of traditional versus geometric morphometric methods in discriminating between closely related bat species – in this case European horseshoe bats (Rhinolophidae, Chiroptera) – based on morphology of the wing, body and tail. In addition to comparing morphometric methods, we used geometric morphometrics to detect interspecies differences as shape changes. Geometric morphometrics yielded improved species discrimination relative to traditional methods. The predicted shape for the variation along the between group principal components revealed that the largest differences between species lay in the extent to which the wing reaches in the direction of the head. This strong trend in interspecific shape variation is associated with size, which we interpret as an evolutionary allometry pattern.

Genome-wide association study using extreme truncate selection identifies novel genes affecting bone mineral density and fracture risk

Osteoporotic fracture is a major cause of morbidity and mortality worldwide. Low bone mineral density (BMD) is a major predisposing factor to fracture and is known to be highly heritable. Site-, gender-, and age-specific genetic effects on BMD are thought to be significant, but have largely not been considered in the design of genome-wide association studies (GWAS) of BMD to date. We report here a GWAS using a novel study design focusing on women of a specific age (postmenopausal women, age 55-85 years), with either extreme high or low hip BMD (age- and gender-adjusted BMD z-scores of +1.5 to +4.0, n = 1055, or -4.0 to -1.5, n = 900), with replication in cohorts of women drawn from the general population (n = 20,898). The study replicates 21 of 26 known BMD-associated genes. Additionally, we report suggestive association of a further six new genetic associations in or around the genes CLCN7, GALNT3, IBSP, LTBP3, RSPO3, and SOX4, with replication in two independent datasets. A novel mouse model with a loss-of-function mutation in GALNT3 is also reported, which has high bone mass, supporting the involvement of this gene in BMD determination. In addition to identifying further genes associated with BMD, this study confirms the efficiency of extreme-truncate selection designs for quantitative trait association studies. © 2011 Duncan et al.

WNT16 influences bone mineral density, cortical bone thickness, bone strength, and osteoporotic fracture risk

We aimed to identify genetic variants associated with cortical bone thickness (CBT) and bone mineral density (BMD) by performing two separate genome-wide association study (GWAS) meta-analyses for CBT in 3 cohorts comprising 5,878 European subjects and for BMD in 5 cohorts comprising 5,672 individuals. We then assessed selected single-nucleotide polymorphisms (SNPs) for osteoporotic fracture in 2,023 cases and 3,740 controls. Association with CBT and forearm BMD was tested for ~2.5 million SNPs in each cohort separately, and results were meta-analyzed using fixed effect meta-analysis. We identified a missense SNP (Thr>Ile; rs2707466) located in the WNT16 gene (7q31), associated with CBT (effect size of -0.11 standard deviations [SD] per C allele, P = 6.2×10-9). This SNP, as well as another nonsynonymous SNP rs2908004 (Gly>Arg), also had genome-wide significant association with forearm BMD (-0.14 SD per C allele, P = 2.3×10-12, and -0.16 SD per G allele, P = 1.2×10-15, respectively). Four genome-wide significant SNPs arising from BMD meta-analysis were tested for association with forearm fracture. SNP rs7776725 in FAM3C, a gene adjacent to WNT16, was associated with a genome-wide significant increased risk of forearm fracture (OR = 1.33, P = 7.3×10-9), with genome-wide suggestive signals from the two missense variants in WNT16 (rs2908004: OR = 1.22, P = 4.9×10-6 and rs2707466: OR = 1.22, P = 7.2×10-6). We next generated a homozygous mouse with targeted disruption of Wnt16. Female Wnt16-/- mice had 27% (P<0.001) thinner cortical bones at the femur midshaft, and bone strength measures were reduced between 43%-61% (6.5×10-13<P<5.9×10-4) at both femur and tibia, compared with their wild-type littermates. Natural variation in humans and targeted disruption in mice demonstrate that WNT16 is an important determinant of CBT, BMD, bone strength, and risk of fracture. © 2012 Zheng et al.

Hyaluronic acid enhances the mechanical properties of tissue-engineered cartilage constructs

There is a need for materials that are well suited for cartilage tissue engineering. Hydrogels have emerged as promising biomaterials for cartilage repair, since, like cartilage, they have high water content, and they allow cells to be encapsulated within the material in a genuinely three-dimensional microenvironment. In this study, we investigated the mechanical properties of tissue-engineered cartilage constructs using in vitro culture models incorporating human chondrocytes from osteoarthritis patients. We evaluated hydrogels formed from mixtures of photocrosslinkable gelatin-methacrylamide (Gel-MA) and varying concentrations (0–2%) of hyaluronic acid methacrylate (HA-MA). Initially, only small differences in the stiffness of each hydrogel existed. After 4 weeks of culture, and to a greater extent 8 weeks of culture, HA-MA had striking and concentration dependent impact on the changes in mechanical properties. For example, the initial compressive moduli of cell-laden constructs with 0 and 1% HA-MA were 29 and 41 kPa, respectively. After 8 weeks of culture, the moduli of these constructs had increased to 66 and 147 kPa respectively, representing a net improvement of 69 kPa for gels with 1% HA-MA. Similarly the equilibrium modulus, dynamic modulus, failure strength and failure strain were all improved in constructs containing HA-MA. Differences in mechanical properties did not correlate with glycosaminoglycan content, which did not vary greatly between groups, yet there were clear differences in aggrecan intensity and distribution as assessed using immunostaining. Based on the functional development with time in culture using human chondrocytes, mixtures of Gel-MA and HA-MA are promising candidates for cartilage tissue-engineering applications.

Spatiotemporal pattern of bacillary dysentery in China from 1990 to 2009: What is the driver behind?

BACKGROUND Little is known about the spatiotemporal pattern of bacillary dysentery (BD) in China. This study assessed the geographic distribution and seasonality of BD in China over the past two decades. METHODS Data on monthly BD cases in 31 provinces of China from January 1990 to December 2009 obtained from Chinese Center for Disease Control and Prevention, and data on demographic and geographic factors, as well as climatic factors, were compiled. The spatial distributions of BD in the four periods across different provinces were mapped, and heat maps were created to present the seasonality of BD by geography. A cosinor function combined with Poisson regression was used to quantify the seasonal parameters of BD, and a regression analysis was conducted to identify the potential drivers of morbidity and seasonality of BD. RESULTS Although most regions of China have experienced considerable declines in BD morbidity over the past two decades, Beijing and Ningxia still had high BD morbidity in 2009. BD morbidity decreased more slowly in North-west China than other regions. BD in China mainly peaked from July to September, with heterogeneity in peak time between regions. Relative humidity was associated with BD morbidity and peak time, and latitude was the major predictor of BD amplitude. CONCLUSIONS The transmission of BD was heterogeneous in China. Improved sanitation and hygiene in North-west China, and better access to clean water and food in the big floating population in some metropolises could be the focus of future preventive interventions against BD. BD control efforts should put more emphasis on those dry areas in summer.

Single-cell-precision microplasma-induced cancer cell apoptosis

The issue of single-cell control has recently attracted enormous interest. However, in spite of the presently achievable intracellular-level physiological probing through bio-photonics, nano-probe-based, and some other techniques, the issue of inducing selective, single-cell-precision apoptosis, without affecting neighbouring cells remains essentially open. Here we resolve this issue and report on the effective single-cell-precision cancer cell treatment using the reactive chemistry of the localized corona-type plasma discharge around a needle-like electrode with the spot size ∼1 µm. When the electrode is positioned with the micrometer precision against a selected cell, a focused and highly-localized micro-plasma discharge induces apoptosis in the selected individual HepG2 and HeLa cancer cells only, without affecting any surrounding cells, even in small cell clusters. This is confirmed by the real-time monitoring of the morphological and structural changes at the cellular and cell nucleus levels after the plasma exposure.

Combining transcriptome assemblies from multiple de novo assemblers in the allo-tetraploid plant Nicotiana benthamiana

Background Nicotiana benthamiana is an allo-tetraploid plant, which can be challenging for de novo transcriptome assemblies due to homeologous and duplicated gene copies. Transcripts generated from such genes can be distinct yet highly similar in sequence, with markedly differing expression levels. This can lead to unassembled, partially assembled or mis-assembled contigs. Due to the different properties of de novo assemblers, no one assembler with any one given parameter space can re-assemble all possible transcripts from a transcriptome. Results In an effort to maximise the diversity and completeness of de novo assembled transcripts, we utilised four de novo transcriptome assemblers, TransAbyss, Trinity, SOAPdenovo-Trans, and Oases, using a range of k-mer sizes and different input RNA-seq read counts. We complemented the parameter space biologically by using RNA from 10 plant tissues. We then combined the output of all assemblies into a large super-set of sequences. Using a method from the EvidentialGene pipeline, the combined assembly was reduced from 9.9 million de novo assembled transcripts to about 235,000 of which about 50,000 were classified as primary. Metrics such as average bit-scores, feature response curves and the ability to distinguish paralogous or homeologous transcripts, indicated that the EvidentialGene processed assembly was of high quality. Of 35 RNA silencing gene transcripts, 34 were identified as assembled to full length, whereas in a previous assembly using only one assembler, 9 of these were partially assembled. Conclusions To achieve a high quality transcriptome, it is advantageous to implement and combine the output from as many different de novo assemblers as possible. We have in essence taking the ‘best’ output from each assembler while minimising sequence redundancy. We have also shown that simultaneous assessment of a variety of metrics, not just focused on contig length, is necessary to gauge the quality of assemblies.

Spatiotemporal transmission dynamics of hemorrhagic fever with renal syndrome in China, 2005-2012

Background Hemorrhagic fever with renal syndrome (HFRS) is a rodent-borne disease caused by many serotypes of hantaviruses. In China, HFRS has been recognized as a severe public health problem with 90% of the total reported cases in the world. This study describes the spatiotemporal dynamics of HFRS cases in China and identifies the regions, time, and populations at highest risk, which could help the planning and implementation of key preventative measures. Methods Data on all reported HFRS cases at the county level from January 2005 to December 2012 were collected from Chinese Center for Disease Control and Prevention. Geographic Information System-based spatiotemporal analyses including Local Indicators of Spatial Association and Kulldorff's space-time scan statistic were performed to detect local high-risk space-time clusters of HFRS in China. In addition, cases from high-risk and low-risk counties were compared to identify significant demographic differences. Results A total of 100,868 cases were reported during 2005–2012 in mainland China. There were significant variations in the spatiotemporal dynamics of HFRS. HFRS cases occurred most frequently in June, November, and December. There was a significant positive spatial autocorrelation of HFRS incidence during the study periods, with Moran's I values ranging from 0.46 to 0.56 (P<0.05). Several distinct HFRS cluster areas were identified, mainly concentrated in northeastern, central, and eastern of China. Compared with cases from low-risk areas, a higher proportion of cases were younger, non-farmer, and floating residents in high-risk counties. Conclusions This study identified significant space-time clusters of HFRS in China during 2005–2012 indicating that preventative strategies for HFRS should be particularly focused on the northeastern, central, and eastern of China to achieve the most cost-effective 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.

Reproducibility of fatmax and fat oxidation rates during exercise in recreationally trained males

Aerobic exercise training performed at the intensity eliciting maximal fat oxidation (Fatmax) has been shown to improve the metabolic profile of obese patients. However, limited information is available on the reproducibility of Fatmax and related physiological measures. The aim of this study was to assess the intra-individual variability of: a) Fatmax measurements determined using three different data analysis approaches and b) fat and carbohydrate oxidation rates at rest and at each stage of an individualized graded test. Fifteen healthy males [body mass index 23.1±0.6 kg/m2, maximal oxygen consumption () 52.0±2.0 ml/kg/min] completed a maximal test and two identical submaximal incremental tests on ergocycle (30-min rest followed by 5-min stages with increments of 7.5% of the maximal power output). Fat and carbohydrate oxidation rates were determined using indirect calorimetry. Fatmax was determined with three approaches: the sine model (SIN), measured values (MV) and 3rd polynomial curve (P3). Intra-individual coefficients of variation (CVs) and limits of agreement were calculated. CV for Fatmax determined with SIN was 16.4% and tended to be lower than with P3 and MV (18.6% and 20.8%, respectively). Limits of agreement for Fatmax were −2±27% of with SIN, −4±32 with P3 and −4±28 with MV. CVs of oxygen uptake, carbon dioxide production and respiratory exchange rate were <10% at rest and <5% during exercise. Conversely, CVs of fat oxidation rates (20% at rest and 24–49% during exercise) and carbohydrate oxidation rates (33.5% at rest, 8.5–12.9% during exercise) were higher. The intra-individual variability of Fatmax and fat oxidation rates was high (CV>15%), regardless of the data analysis approach employed. Further research on the determinants of the variability of Fatmax and fat oxidation rates is required.

Laser capture microdissection and multiplex-tandem PCR analysis of proximal tubular epithelial cell signaling in human kidney disease

Interstitial fibrosis, a histological process common to many kidney diseases, is the precursor state to end stage kidney disease, a devastating and costly outcome for the patient and the health system. Fibrosis is historically associated with chronic kidney disease (CKD) but emerging evidence is now linking many forms of acute kidney disease (AKD) with the development of CKD. Indeed, we and others have observed at least some degree of fibrosis in up to 50% of clinically defined cases of AKD. Epithelial cells of the proximal tubule (PTEC) are central in the development of kidney interstitial fibrosis. We combine the novel techniques of laser capture microdissection and multiplex-tandem PCR to identify and quantitate “real time��� gene transcription profiles of purified PTEC isolated from human kidney biopsies that describe signaling pathways associated with this pathological fibrotic process. Our results: (i) confirm previous in-vitro and animal model studies; kidney injury molecule-1 is up-regulated in patients with acute tubular injury, inflammation, neutrophil infiltration and a range of chronic disease diagnoses, (ii) provide data to inform treatment; complement component 3 expression correlates with inflammation and acute tubular injury, (iii) identify potential new biomarkers; proline 4-hydroxylase transcription is down-regulated and vimentin is up-regulated across kidney diseases, (iv) describe previously unrecognized feedback mechanisms within PTEC; Smad-3 is down-regulated in many kidney diseases suggesting a possible negative feedback loop for TGF-β in the disease state, whilst tight junction protein-1 is up-regulated in many kidney diseases, suggesting feedback interactions with vimentin expression. These data demonstrate that the combined techniques of laser capture microdissection and multiplex-tandem PCR have the power to study molecular signaling within single cell populations derived from clinically sourced tissue.

Stochastic dynamics of interacting haematopoietic stem cell niche lineages

Since we still know very little about stem cells in their natural environment, it is useful to explore their dynamics through modelling and simulation, as well as experimentally. Most models of stem cell systems are based on deterministic differential equations that ignore the natural heterogeneity of stem cell populations. This is not appropriate at the level of individual cells and niches, when randomness is more likely to affect dynamics. In this paper, we introduce a fast stochastic method for simulating a metapopulation of stem cell niche lineages, that is, many sub-populations that together form a heterogeneous metapopulation, over time. By selecting the common limiting timestep, our method ensures that the entire metapopulation is simulated synchronously. This is important, as it allows us to introduce interactions between separate niche lineages, which would otherwise be impossible. We expand our method to enable the coupling of many lineages into niche groups, where differentiated cells are pooled within each niche group. Using this method, we explore the dynamics of the haematopoietic system from a demand control system perspective. We find that coupling together niche lineages allows the organism to regulate blood cell numbers as closely as possible to the homeostatic optimum. Furthermore, coupled lineages respond better than uncoupled ones to random perturbations, here the loss of some myeloid cells. This could imply that it is advantageous for an organism to connect together its niche lineages into groups. Our results suggest that a potential fruitful empirical direction will be to understand how stem cell descendants communicate with the niche and how cancer may arise as a result of a failure of such communication.

Altered formalin-induced pain and Fos induction in the periaqueductal grey of preadolescent rats following neonatal LPS exposure

Animal and human studies have demonstrated that early pain experiences can produce alterations in the nociceptive systems later in life including increased sensitivity to mechanical, thermal, and chemical stimuli. However, less is known about the impact of neonatal immune challenge on future responses to noxious stimuli and the reactivity of neural substrates involved in analgesia. Here we demonstrate that rats exposed to Lipopolysaccharide (LPS; 0.05 mg/kg IP, Salmonella enteritidis) during postnatal day (PND) 3 and 5 displayed enhanced formalin-induced flinching but not licking following formalin injection at PND 22. This LPS-induced hyperalgesia was accompanied by distinct recruitment of supra-spinal regions involved in analgesia as indicated by significantly attenuated Fos-protein induction in the rostral dorsal periaqueductal grey (DPAG) as well as rostral and caudal axes of the ventrolateral PAG (VLPAG). Formalin injections were associated with increased Fos-protein labelling in lateral habenula (LHb) as compared to medial habenula (MHb), however the intensity of this labelling did not differ as a result of neonatal immune challenge. These data highlight the importance of neonatal immune priming in programming inflammatory pain sensitivity later in development and highlight the PAG as a possible mediator of this process