Prospects for estimating nucleotide divergence with RAPDs

The technique of random amplification of polymorphic DNA (RAPD), which is simply polymerase chain reaction (PCR) amplification of genomic DNA by a single short oligonucleotide primer, produces complex patterns of anonymous polymorphic DNA fragments. The information provided by these banding patterns has proved to be of great utility for mapping and for verification of identity of bacterial strains. Here we consider whether the degree of similarity of the banding patterns can be used to estimate nucleotide diversity and nucleotide divergence. With haploid data, fragments generated by RAPD-PCR can be treated in a fashion very similar to that for restriction-fragment data. Amplification of diploid samples, on the other hand, requires consideration of the fact that presence of a band is dominant to absence of the band. After describing a method for estimating nucleotide divergence on the basis of diploid samples, we summarize the restrictions and criteria that must be met when RAPD data are used for estimating population genetic parameters.      

Microsatellite variation in Drosophila melanogaster and Drosophila simulans: a reciprocal test of the ascertainment bias hypothesis

Interspecific comparisons of microsatellite loci have repeatedly shown that the loci are longer and more variable in the species from which they are derived (the focal species) than are homologous loci in other (nonfocal) species. There is debate as to whether this is due to directional evolution or to an ascertainment bias during the cloning and locus selection processes. This study tests these hypotheses by performing a reciprocal study. Eighteen perfect dinucleotide microsatellite loci identified from a Drosophila simulans library screen and 18 previously identified in an identical Drosophila melanogaster library screen were used to survey natural populations of each species. No difference between focal and nonfocal species was observed for mean PCR fragment length. However, heterozygosity and number of alleles were significantly higher in the focal species than in the nonfocal species. The most common allele in the Zimbabwe population of both species was sequenced for 31 of the 36 loci. The length of the longest stretch of perfect repeat units is, on average, longer in the focal species than in the non-focal species. There is a positive correlation between the length of the longest stretch of perfect repeats and heterozygosity. The difference in heterozygosity can thus be explained by a reduction in the length of the longest stretch of perfect repeats in the nonfocal species. Furthermore, flanking-sequence length difference was noted between the two species at 58% of the loci sequenced. These data do not support the predictions of the directional-evolution hypothesis; however, consistent with the ascertainment bias hypothesis, the lower variability in nonfocal species is an artifact of the microsatellite cloning and isolation process. Our results also suggest that the magnitude of ascertainment bias for repeat unit length is a function of the microsatellite size distribution in the genomes of different species.      

Evaluation of coagulation activation after Rhinovirus infection in patients with asthma and healthy control subjects: an observational study

Background

Asthma exacerbations are frequently triggered by rhinovirus infections. Both asthma and respiratory tract infection can activate haemostasis. Therefore we hypothesized that experimental rhinovirus-16 infection and asthmatic airway inflammation act in synergy on the haemostatic balance.

Methods

28 patients (14 patients with mild allergic asthma and 14 healthy non-allergic controls) were infected with low-dose rhinovirus type 16. Venous plasma and bronchoalveolar lavage fluid (BAL fluid) were obtained before and 6 days after infection to evaluate markers of coagulation activation, thrombin-antithrombin complexes, von Willebrand factor, plasmin-antiplasmin complexes, plasminogen activator inhibitor type-1, endogenous thrombin potential and tissue factor-exposing microparticles by fibrin generation test, in plasma and/or BAL fluid. Data were analysed by nonparametric tests (Wilcoxon, Mann Whitney and Spearman correlation).

Results

13 patients with mild asthma (6 females, 19-29 y) and 11 healthy controls (10 females, 19-31 y) had a documented Rhinovirus-16 infection. Rhinovirus-16 challenge resulted in a shortening of the fibrin generation test in BAL fluid of asthma patients (t = -1: 706 s vs. t = 6: 498 s; p = 0.02), but not of controls (t = -1: 693 s vs. t = 6: 636 s; p = 0.65). The fold change in tissue factor-exposing microparticles in BAL fluid inversely correlated with the fold changes in eosinophil cationic protein and myeloperoxidase in BAL fluid after virus infection (r = -0.517 and -0.528 resp., both p = 0.01).Rhinovirus-16 challenge led to increased plasminogen activator inhibitor type-1 levels in plasma in patients with asthma (26.0 ng/mL vs. 11.5 ng/mL in healthy controls, p = 0.04). Rhinovirus-16 load in BAL showed a linear correlation with the fold change in endogenous thrombin potential, plasmin-antiplasmin complexes and plasminogen activator inhibitor type-1.

Conclusions

Experimental rhinovirus infection induces procoagulant changes in the airways of patients with asthma through increased activity of tissue factor-exposing microparticles. These microparticle-associated procoagulant changes are associated with both neutrophilic and eosinophilic inflammation. Systemic activation of haemostasis increases with Rhinoviral load.

Trial registration

This trial was registered at the Dutch trial registry (http://www.trialregister.nl): NTR1677.     

Dermatological conditions during TNF-α-blocking therapy in patients with rheumatoid arthritis: a prospective study

Various dermatological conditions have been reported during tumor necrosis factor (TNF)-α-blocking therapy, but until now no prospective studies have been focused on this aspect. The present study was set up to investigate the number and nature of clinically important dermatological conditions during TNF-α-blocking therapy in patients with rheumatoid arthritis (RA). RA patients starting on TNF-α-blocking therapy were prospectively followed up. The numbers and natures of dermatological events giving rise to a dermatological consultation were recorded. The patients with a dermatological event were compared with a group of prospectively followed up RA control patients, naive to TNF-α-blocking therapy and matched for follow-up period. 289 RA patients started TNF-α-blocking therapy. 128 dermatological events were recorded in 72 patients (25%) during 911 patient-years of follow-up. TNF-α-blocking therapy was stopped in 19 (26%) of these 72 patients because of the dermatological event. More of the RA patients given TNF-α-blocking therapy (25%) than of the anti-TNF-α-naive patients (13%) visited a dermatologist during follow-up (P < 0.0005). Events were recorded more often during active treatment (0.16 events per patient-year) than during the period of withdrawal of TNF-α-blocking therapy (0.09 events per patient-year, P < 0.0005). The events recorded most frequently were skin infections (n = 33), eczema (n = 20), and drug-related eruptions (n = 15). Other events with a possible relation to TNF-α-blocking therapy included vasculitis, psoriasis, drug-induced systemic lupus erythematosus, dermatomyositis, and a lymphomatoid-papulosis-like eruption. This study is the first large prospective study focusing on dermatological conditions during TNF-α-blocking therapy. It shows that dermatological conditions are a significant and clinically important problem in RA patients receiving TNF-α-blocking therapy.     

Initial steps towards a production platform for DNA sequence analysis on the grid

Background  

Bioinformatics is confronted with a new data explosion due to the availability of high throughput DNA sequencers. Data storage and analysis becomes a problem on local servers, and therefore it is needed to switch to other IT infrastructures. Grid and workflow technology can help to handle the data more efficiently, as well as facilitate collaborations. However, interfaces to grids are often unfriendly to novice users.     

Near-terminal creep damage does not substantially influence fatigue life under physiological loading

Cortical bone specimens were damaged using repeated blocks of tensile creep loading until a near-terminal amount of creep damage was generated (corresponding to a reduction in elastic modulus of 15%). One group of cortical bone specimens was submitted to the near-terminal damage protocol and subsequently underwent fatigue loading in tension with a maximum strain of 2000 με (Damage Fatigue, n=5). A second group was submitted to cyclic fatigue loading but was not pre-damaged (Control Fatigue, n=5). All but one specimen (a damaged specimen) reached run-out (10 million cycles, 7.7 days). No significant differences in microscopic cracks or other tissue damage were observed between the two groups or between either group and additional, completely unloaded specimens. Our results suggest that damage in cortical bone allograft that is not obvious or associated with a stress riser may not substantially affect its fatigue life under physiologic loading.     

Microdamage Caused by Fatigue Loading in Human Cancellous Bone: Relationship to Reductions in Bone Biomechanical Performance

Vertebral fractures associated with osteoporosis are often the result of tissue damage accumulated over time. Microscopic tissue damage (microdamage) generated in vivo is believed to be a mechanically relevant aspect of bone quality that may contribute to fracture risk. Although the presence of microdamage in bone tissue has been documented, the relationship between loading, microdamage accumulation and mechanical failure is not well understood. The aim of the current study was to determine how microdamage accumulates in human vertebral cancellous bone subjected to cyclic fatigue loading. Cancellous bone cores (n = 32) from the third lumbar vertebra of 16 donors (10 male, 6 female, age 76±8.8, mean ± SD) were subjected to compressive cyclic loading at σ/E₀ = 0.0035 (where σ is stress and E₀ is the initial Young’s modulus). Cyclic loading was suspended before failure at one of seven different amounts of loading and specimens were stained for microdamage using lead uranyl acetate. Damage volume fraction (DV/BV) varied from 0.8±0.5% (no loading) to 3.4±2.1% (fatigue-loaded to complete failure) and was linearly related to the reductions in Young’s modulus caused by fatigue loading (r² = 0.60, p<0.01). The relationship between reductions in Young’s modulus and proportion of fatigue life was nonlinear and suggests that most microdamage generation occurs late in fatigue loading, during the tertiary phase. Our results indicate that human vertebral cancellous bone tissue with a DV/BV of 1.5% is expected to have, on average, a Young’s modulus 31% lower than the same tissue without microdamage and is able to withstand 92% fewer cycles before failure than the same tissue without microdamage. Hence, even small amounts of microscopic tissue damage in human vertebral cancellous bone may have large effects on subsequent biomechanical performance.     

Number and mode of inheritance of QTL influencing backfat thickness on SSC2p in Sino-European pig pedigrees

Background

In the pig, multiple QTL associated with growth and fatness traits have been mapped to chromosome 2 (SSC2) and among these, at least one shows paternal expression due to the IGF2-intron3-G3072A substitution. Previously published results on the position and imprinting status of this QTL disagree between analyses from French and Dutch F2 crossbred pig populations obtained with the same breeds (Meishan crossed with Large White or Landrace).

Methods

To study the role of paternal and maternal alleles at the IGF2 locus and to test the hypothesis of a second QTL affecting backfat thickness on the short arm of SSC2 (SSC2p), a QTL mapping analysis was carried out on a combined pedigree including both the French and Dutch F2 populations, on the progeny of F1 males that were heterozygous (A/G) and homozygous (G/G) at the IGF2 locus. Simulations were performed to clarify the relations between the two QTL and to understand to what extent they can explain the discrepancies previously reported.

Results

The QTL analyses showed the segregation of at least two QTL on chromosome 2 in both pedigrees, i.e. the IGF2 locus and a second QTL segregating at least in the G/G F1 males and located between positions 30 and 51 cM. Statistical analyses highlighted that the maternally inherited allele at the IGF2 locus had a significant effect but simulation studies showed that this is probably a spurious effect due to the segregation of the second QTL.

Conclusions

Our results show that two QTL on SSC2p affect backfat thickness. Differences in the pedigree structures and in the number of heterozygous females at the IGF2 locus result in different imprinting statuses in the two pedigrees studied. The spurious effect observed when a maternally allele is present at the IGF2 locus, is in fact due to the presence of a second closely located QTL. This work confirms that pig chromosome 2 is a major region associated with fattening traits.     

A gene frequency model for QTL mapping using Bayesian inference

Background

Information for mapping of quantitative trait loci (QTL) comes from two sources: linkage disequilibrium (non-random association of allele states) and cosegregation (non-random association of allele origin). Information from LD can be captured by modeling conditional means and variances at the QTL given marker information. Similarly, information from cosegregation can be captured by modeling conditional covariances. Here, we consider a Bayesian model based on gene frequency (BGF) where both conditional means and variances are modeled as a function of the conditional gene frequencies at the QTL. The parameters in this model include these gene frequencies, additive effect of the QTL, its location, and the residual variance. Bayesian methodology was used to estimate these parameters. The priors used were: logit-normal for gene frequencies, normal for the additive effect, uniform for location, and inverse chi-square for the residual variance. Computer simulation was used to compare the power to detect and accuracy to map QTL by this method with those from least squares analysis using a regression model (LSR).

Results

To simplify the analysis, data from unrelated individuals in a purebred population were simulated, where only LD information contributes to map the QTL. LD was simulated in a chromosomal segment of 1 cM with one QTL by random mating in a population of size 500 for 1000 generations and in a population of size 100 for 50 generations. The comparison was studied under a range of conditions, which included SNP density of 0.1, 0.05 or 0.02 cM, sample size of 500 or 1000, and phenotypic variance explained by QTL of 2 or 5%. Both 1 and 2-SNP models were considered. Power to detect the QTL for the BGF, ranged from 0.4 to 0.99, and close or equal to the power of the regression using least squares (LSR). Precision to map QTL position of BGF, quantified by the mean absolute error, ranged from 0.11 to 0.21 cM for BGF, and was better than the precision of LSR, which ranged from 0.12 to 0.25 cM.

Conclusions

In conclusion given a high SNP density, the gene frequency model can be used to map QTL with considerable accuracy even within a 1 cM region.