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.     

Hip joint replacement surgery for idiopathic osteoarthritis aggregates in families

In order to determine whether there is a genetic component to hip or knee joint failure due to idiopathic osteoarthritis (OA), we invited patients (probands) undergoing hip or knee arthroplasty for management of idiopathic OA to provide detailed family histories regarding the prevalence of idiopathic OA requiring joint replacement in their siblings. We also invited their spouses to provide detailed family histories about their siblings to serve as a control group. In the probands, we confirmed the diagnosis of idiopathic OA using American College of Rheumatology criteria. The cohorts included the siblings of 635 probands undergoing total hip replacement, the siblings of 486 probands undergoing total knee replacement, and the siblings of 787 spouses. We compared the prevalence of arthroplasty for idiopathic OA among the siblings of the probands with that among the siblings of the spouses, and we used logistic regression to identify independent risk factors for hip and knee arthroplasty in the siblings. Familial aggregation for hip arthroplasty, but not for knee arthroplasty, was observed after controlling for age and sex, suggesting a genetic contribution to end-stage hip OA but not to end-stage knee OA. We conclude that attempts to identify genes that predispose to idiopathic OA resulting in joint failure are more likely to be successful in patients with hip OA than in those with knee OA.     

Biosynthesis of high density lipoprotein by chicken liver: nature of nascent intracellular high density lipoprotein

Young chickens were administered L-[(3)H]leucine and after 10 or 30 min the livers were removed and fractioned into rough (RER) and smooth (SER) endoplasmic reticulum fractions and into light, intermediate, and heavy golgo cell fractions. The labeled high density lipoprotein (HDL), contained within these intracellular organelles was isolated either by immunoprecipitation using rabbit antiserum to rooster HDL, or by ultracentrifugal glotation between densities 1.063 and 1.21 g/ml. The radioactive apoproteins of nascent HDL were analyzed by SDS PAGE and detected by fluorography. Analyses of radioactive apoproteins obtained by immunoprecipitation from the contents of the RER, the SER, and the three golgi complex fractions revealed only one apoprotein, A1. The C peptide present in serum HDL was not detected intracellularly. The radioactive apoprotein A1 which is present within the cisternae of the RER and the SER fractions failed to float, whereas apoprotein A1, present within the golgi apparatus, readily floated between densities 1.063 and 1.21 g/ml. The HDL particles, isolated by flotation from the golgi apparatus content, were further characterized by lipid and protein analyses and by electron microscopy. Golgi HDL particles have the same density as serum HDL. On a percentage basis, golgi HDL contains less protein and more phospholipids than does serum HDL. Morphologically, golgi HDL is different in appearance from serum HDL. It is more heterogeneous in size, with most of the particles ranging 8.3-25 nm in diameter. The spherical particles contain small membrane tails. Occasionally, a few disk-shaped bilayer structures are also found within the golgi apparatus. These studies show that the newly synthesized apoprotein A1, present within the RER and the SER cell fractions, is not fully complexed with lipid and that apoprotein A1 does not acquire sufficient lipid to float at the proper HDL density until it enters the golgi apparatus. The difference in chemical composition and the heterogeneous size of golgi HDL may be attributed to the different stages of HDL maturation.     

Energy homeostasis and gastrointestinal endocrine differentiation do not require the anorectic hormone peptide YY

The gastrointestinal hormone peptide YY is a potent inhibitor of food intake and is expressed early during differentiation of intestinal and pancreatic endocrine cells. In order to better understand the role of peptide YY in energy homeostasis and development, we created mice with a targeted deletion of the peptide YY gene. All intestinal and pancreatic endocrine cells developed normally in the absence of peptide YY with the exception of pancreatic polypeptide (PP) cells, indicating that peptide YY expression was not required for terminal differentiation. We used recombination-based cell lineage trace to determine if peptide YY cells were progenitors for gastrointestinal endocrine cells. Peptide YY(+) cells gave rise to all L-type enteroendocrine cells and to islet partial differential and PP cells. In the pancreas, approximately 40% of pancreatic alpha and rare beta cells arose from peptide YY(+) cells, suggesting that most beta cells and surprisingly the majority of alpha cells are not descendants of peptide YY(+)/glucagon-positive/insulin-positive cells that appear during early pancreagenesis. Despite the anorectic effects of exogenous peptide YY(3-36) following intraperitoneal administration, mice lacking peptide YY showed normal growth, food intake, energy expenditure, and responsiveness to peptide YY(3-36). These observations suggest that targeted disruption of the peptide YY gene does not perturb terminal endocrine cell differentiation or the control of food intake and energy homeostasis.     

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.     

Disruption of chromosome 11 in canine fibrosarcomas highlights an unusual variability of CDKN2B in dogs

Background  

In dogs in the western world neoplasia constitutes the most frequently diagnosed cause of death. Although there appear to be similarities between canine and human cancers, rather little is known about the cytogenetic and molecular alterations in canine tumours. Different dog breeds are susceptible to different types of cancer, but the genetic basis of the great majority of these predispositions has yet to be discovered. In some retriever breeds there is a high incidence of soft tissue sarcomas and we have previously reported alterations of chromosomes 11 and 30 in two poorly differentiated fibrosarcomas. Here we extend our observations and present a case report on detail rearrangements on chromosome 11 as well as genetic variations in a tumour suppressor gene in normal dogs.     

Low sclerostin levels: a predictive marker of persistent inflammation in ankylosing spondylitis during anti-tumor necrosis factor therapy?

Introduction

Sclerostin levels have been reported to be low in ankylosing spondylitis (AS), but there is no data regarding the possible role of this Wnt inhibitor during anti-tumor necrosis factor (TNF) therapy. The present study longitudinally evaluated sclerostin levels, inflammatory markers and bone mineral density (BMD) in AS patients under anti-TNF therapy.

Methods

Thirty active AS patients were assessed at baseline, 6 and 12 months after anti-TNF therapy regarding clinical parameters, inflammatory markers, BMD and baseline radiographic damage (mSASSS). Thirty age- and sex-matched healthy individuals comprised the control group. Patients'' sclerostin levels, sclerostin binding low-density lipoprotein receptor-related protein 6 (LRP6) and BMD were evaluated at the same time points and compared to controls.

Results

At baseline, AS patients had lower sclerostin levels (60.5 ± 32.7 vs. 96.7 ± 52.9 pmol/L, P = 0.002) and comparable sclerostin binding to LRP6 (P = 0.387) than controls. Improvement of Bath Ankylosing Spondylitis Disease Activity Index (BASDAI), Bath Ankylosing Spondylitis Functional Index (BASFI), Bath Ankylosing Spondylitis Metrology Index (BASMI), Ankylosing Spondylitis quality of life (ASQoL) was observed at baseline vs. 6 vs. 12 months (P < 0.01). Concomitantly, a gradual increase in spine BMD (P < 0.001) and a positive correlation between baseline mSASSS and spine BMD was found (r = 0.468, P < 0.01). Inflammatory parameters reduction was observed comparing baseline vs. 6 vs. 12 months (P <0.01). Sclerostin levels progressively increased [baseline (60.5 ± 32.7) vs. 6 months (67.1 ± 31.9) vs. 12 months (72.7 ± 32.3) pmol/L, P <0.001]. At 12 months, the sclerostin levels remained significantly lower in patients compared to controls (72.7 ± 32.3 vs. 96.70 ± 52.85 pmol/L, P = 0.038). Moreover, sclerostin serum levels at 12 months were lower in the 10 patients with high C reactive protein (CRP) (≥ 5 mg/l) compared to the other 20 patients with normal CRP (P = 0.004). Of note, these 10 patients with persistent inflammation also had lower sclerostin serum levels at baseline compared to the other patients (P = 0.023). Univariate logistic regression analysis demonstrated that AS patients with lower sclerostin serum levels had an increased risk to have high CRP at 12 months (odds ratio = 7.43, 95% CI 1.23 to 45.01, P = 0.020) than those with higher sclerostin values.

Conclusions

Persistent low sclerostin levels may underlie continuous inflammation in AS patients under anti-TNF therapy.     

Gene expression profiling of early intervertebral disc degeneration reveals a down-regulation of canonical Wnt signaling and caveolin-1 expression: implications for development of regenerative strategies

Introduction

Early degeneration of the intervertebral disc (IVD) involves a change in cellular differentiation from notochordal cells (NCs) in the nucleus pulposus (NP) to chondrocyte-like cells (CLCs). The purpose of this study was to investigate the gene expression profiles involved in this process using NP tissue from non-chondrodystrophic and chondrodystrophic dogs, a species with naturally occurring IVD degeneration.

Methods

Dual channel DNA microarrays were used to compare 1) healthy NP tissue containing only NCs (NC-rich), 2) NP tissue with a mixed population of NCs and CLCs (Mixed), and 3) NP tissue containing solely CLCs (CLC-rich) in both non-chondrodystrophic and chondrodystrophic dogs. Based on previous reports and the findings of the microarray analyses, canonical Wnt signaling was further evaluated using qPCR of relevant Wnt target genes. We hypothesized that caveolin-1, a regulator of Wnt signaling that showed significant changes in gene expression in the microarray analyses, played a significant role in early IVD degeneration. Caveolin-1 expression was investigated in IVD tissue sections and in cultured NCs. To investigate the significance of Caveolin-1 in IVD health and degeneration, the NP of 3-month-old Caveolin-1 knock-out mice was histopathologically evaluated and compared with the NP of wild-type mice of the same age.

Results

Early IVD degeneration involved significant changes in numerous pathways, including Wnt/β-catenin signaling. With regard to Wnt/β-catenin signaling, axin2 gene expression was significantly higher in chondrodystrophic dogs compared with non-chondrodystrophic dogs. IVD degeneration involved significant down-regulation of axin2 gene expression. IVD degeneration involved significant down-regulation in Caveolin-1 gene and protein expression. NCs showed abundant caveolin-1 expression in vivo and in vitro, whereas CLCs did not. The NP of wild-type mice was rich in viable NCs, whereas the NP of Caveolin-1 knock-out mice contained chondroid-like matrix with mainly apoptotic, small, rounded cells.

Conclusions

Early IVD degeneration involves down-regulation of canonical Wnt signaling and Caveolin-1 expression, which appears to be essential to the physiology and preservation of NCs. Therefore, Caveolin-1 may be regarded an exciting target for developing strategies for IVD regeneration.     

Inactivation of the sodium channel. I. Sodium current experiments

Inactivation of sodium conductance has been studied in squid axons with voltage clamp techniques and with the enzyme pronase which selectively destroys inactivation. Comparison of the sodium current before and after pronase treatment shows a lag of several hundred microseconds in the onset of inactivation after depolarization. This lag can of several hundred microseconds in the onset of inactivation after polarization. This lag can also be demonstrated with double-pulse experiments. When the membrane potential is hyperpolarized to -140 mV before depolarization, both activation and inactivation are delayed. These findings suggest that inactivation occurs only after activation are delayed. These findings suggest that inactivation occurs only after activation; i.e. that the channels must open before they can inactivate. The time constant of inactivation measured with two pulses (τ(c)) is the same as the one measured from the decay of the sodium current during a single pulse (τ(h)). For large depolarizations, steady-state inactivation becomes more incomplete as voltage increases; but it is relatively complete and appears independent of voltage when determined with a two- pulse method. This result confirms the existence of a second open state for Na channels, as proposed by Chandler and Meves (1970. J. Physiol. [Lond.]. 211:653-678). The time constant of recovery from inactivation is voltage dependent and decreases as the membrane potential is made more negative. A model for Na channels is presented which has voltage-dependent transitions between the closed and open states, and a voltage-independent transition between the open and the inactivated state. In this model the voltage dependence of inactivation is a consequence of coupling to the activation process.