Joint model with latent state for longitudinal and multistate data

In many chronic diseases, the patient's health status is followed up by quantitative markers. The evolution is often characterized by a 2-phase degradation process, that is, a normal phase followed by a pathological degradation phase preceding the disease diagnosis. We propose a joint multistate model with latent state for the joint modeling of repeated measures of a quantitative marker, time-to-illness and time-to-death. Using data from the PAQUID cohort on cognitive aging, we jointly studied cognitive decline, dementia risk, and death risk. We estimated the mean evolution of cognitive scores given age at dementia for subjects alive and demented, the mean evolution of cognitive scores for subjects alive and nondemented, in addition to age at acceleration of cognitive decline and duration of the pre-dementia phase.     

Association between IgM anti-herpes simplex virus and plasma amyloid-beta levels

Objective

Herpes simplex virus (HSV) reactivation has been identified as a possible risk factor for Alzheimer''s disease (AD) and plasma amyloid-beta (Aβ) levels might be considered as possible biomarkers of the risk of AD. The aim of our study was to investigate the association between anti-HSV antibodies and plasma Aβ levels.

Methods

The study sample consisted of 1222 subjects (73.9 y in mean) from the Three-City cohort. IgM and IgG anti-HSV antibodies were quantified using an ELISA kit, and plasma levels of Aβ_1–40 and Aβ_1–42 were measured using an xMAP-based assay technology. Cross-sectional analyses of the associations between anti-HSV antibodies and plasma Aβ levels were performed by multi-linear regression.

Results

After adjustment for study center, age, sex, education, and apolipoprotein E-e4 polymorphism, plasma Aβ_1–42 and Aβ_1–40 levels were specifically inversely associated with anti-HSV IgM levels (β = −20.7, P = 0.001 and β = −92.4, P = 0.007, respectively). In a sub-sample with information on CLU- and CR1-linked SNPs genotyping (n = 754), additional adjustment for CR1 or CLU markers did not modify these associations (adjustment for CR1 rs6656401, β = −25.6, P = 0.002 for Aβ_1–42 and β = −132.7, P = 0.002 for Aβ_1–40; adjustment for CLU rs2279590, β = −25.6, P = 0.002 for Aβ_1–42 and β = −134.8, P = 0.002 for Aβ_1–40). No association between the plasma Aβ_1–42-to-Aβ_1–40 ratio and anti-HSV IgM or IgG were evidenced.

Conclusion

High anti-HSV IgM levels, markers of HSV reactivation, are associated with lower plasma Aβ_1–40 and Aβ_1–42 levels, which suggest a possible involvement of the virus in the alterations of the APP processing and potentially in the pathogenesis of AD in human.     

Chytrid epidemics may increase genetic diversity of a diatom spring-bloom

Contrary to expectation, populations of clonal organisms are often genetically highly diverse. In phytoplankton, this diversity is maintained throughout periods of high population growth (that is, blooms), even though competitive exclusion among genotypes should hypothetically lead to the dominance of a few superior genotypes. Genotype-specific parasitism may be one mechanism that helps maintain such high-genotypic diversity of clonal organisms. Here, we present a comparison of population genetic similarity by estimating the beta-dispersion among genotypes of early and peak bloom populations of the diatom Asterionella formosa for three spring-blooms under high or low parasite pressure. The Asterionella population showed greater beta-dispersion at peak bloom than early bloom in the 2 years with high parasite pressure, whereas the within group dispersion did not change under low parasite pressure. Our findings support that high prevalence parasitism can promote genetic diversification of natural populations of clonal hosts.     

Expression and localization of trefoil factor family genes in rat submandibular glands

The trefoil factor (TFF) family, which comprises TFF1, TFF2 and TFF3, plays an essential role in epithelial regeneration within the gastrointestinal tract. All three TFFs are present in human saliva; TFF3 is the predominant trefoil peptide. Little is known about the expression and tissue distribution of TFFs in rats, which are commonly used as a model system for human studies. We investigated the localization of the TFF genes that encode secretory peptides in rat submandibular glands (SMG). All three TFFs were expressed in rat SMG, although their location varied. Substantial amounts of TFF1 were detected only in the cytoplasm of epithelial cells in the SMG granular convoluted tubules (GCT), while TFF2 and TFF3 were widely distributed in the cytoplasm of epithelial cells of intercalated ducts (ID), striated ducts (SD) and interlobular ducts (ILD). The three TFFs also were detected especially in the lumens of the SD and ILD. Semi-quantitative RT-PCR and in situ hybridization experiments confirmed TFF1, TFF2 and TFF3 mRNA expressions in the SMG. Greater expression of TFF peptides and mRNA was observed in male rats than in females. The broad expression of TFFs in rat SMG cells and lumens suggests that TFFs function in this organ by their secretion into the duct lumens. We also found differences in TFF expression profiles between rat and human SMG; therefore, caution should be exercised when using rats as a model for human TFF studies.     

tigaR: integrative significance analysis of temporal differential gene expression induced by genomic abnormalities

Background

To determine which changes in the host cell genome are crucial for cervical carcinogenesis, a longitudinal in vitro model system of HPV-transformed keratinocytes was profiled in a genome-wide manner. Four cell lines affected with either HPV16 or HPV18 were assayed at 8 sequential time points for gene expression (mRNA) and gene copy number (DNA) using high-resolution microarrays. Available methods for temporal differential expression analysis are not designed for integrative genomic studies.

Results

Here, we present a method that allows for the identification of differential gene expression associated with DNA copy number changes over time. The temporal variation in gene expression is described by a generalized linear mixed model employing low-rank thin-plate splines. Model parameters are estimated with an empirical Bayes procedure, which exploits integrated nested Laplace approximation for fast computation. Iteratively, posteriors of hyperparameters and model parameters are estimated. The empirical Bayes procedure shrinks multiple dispersion-related parameters. Shrinkage leads to more stable estimates of the model parameters, better control of false positives and improvement of reproducibility. In addition, to make estimates of the DNA copy number more stable, model parameters are also estimated in a multivariate way using triplets of features, imposing a spatial prior for the copy number effect.

Conclusion

With the proposed method for analysis of time-course multilevel molecular data, more profound insight may be gained through the identification of temporal differential expression induced by DNA copy number abnormalities. In particular, in the analysis of an integrative oncogenomics study with a time-course set-up our method finds genes previously reported to be involved in cervical carcinogenesis. Furthermore, the proposed method yields improvements in sensitivity, specificity and reproducibility compared to existing methods. Finally, the proposed method is able to handle count (RNAseq) data from time course experiments as is shown on a real data set.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2105-15-327) contains supplementary material, which is available to authorized users.     

Do leucocytes reflect condition in nestling burrowing parrots Cyanoliseus patagonus in the wild?

The different leucocyte types are an important part of the immune system. Thus, they have been used in ecological studies to assess immune function and physiological stress in wild birds. It is generally assumed that increased stress and decreased condition are associated with an increase in the ratio of heterophils to lymphocytes, the H/L ratio. We studied leucocyte profiles in relation to body condition in nestling Burrowing Parrots (Cyanoliseus patagonus) in North-eastern Patagonia, Argentina. As in other wild parrots, heterophils were the most numerous leucocyte type, suggesting strong investment into innate immunity. Leucocyte profiles did not change with the age, while nestlings in better body condition increased the number of heterophils. Because the number of lymphocytes was independent of body condition, as a result we observed a positive correlation between body condition and the H/L ratio. The total number of leucocytes relative to erythrocytes increased in nestlings in better body condition, indicating a larger overall investment into immune function in well-nourished nestlings. The observed heterophilic profiles of nestling Burrowing Parrots together with the positive relationship between H/L ratio and body condition may indicate a favoured investment in a robust innate immunity that reduces the risk of infection taking hold in these long-lived birds.     

Is Intracranial Atherosclerosis an Independent Risk Factor for Cerebral Atrophy? A Retrospective Evaluation

Background

Our purpose was to study the association between the intracranial atherosclerosis as measured by cavernous carotid artery calcification (ICAC) observed on head CT and atrophic changes of supra-tentorial brain demonstrated by MRI.

Methods

Institutional review board approval was obtained for this retrospective study incorporating 65 consecutive patients presenting acutely who had both head CT and MRI. Arterial calcifications of the intracranial cavernous carotids (ICAC) were assigned a number (1 to 4) in the bone window images from CT scans. These 4 groups were then combined into high (grades 3 and 4) and low calcium (grades 1 and 2) subgroups. Brain MRI was independently evaluated to identify cortical and central atrophy. Demographics and cardiovascular risk factors were evaluated in subjects with high and low ICAC. Relationship between CT demonstrated ICAC and brain atrophy patterns were evaluated both without and with adjustment for cerebral ischemic scores and cardiovascular risk factors.

Results

Forty-six of the 65 (71%) patients had high ICAC on head CT. Subjects with high ICAC were older, and had higher prevalence of hypertension, diabetes, coronary artery disease (CAD), atrial fibrillation and history of previous stroke (CVA) compared to those with low ICAC. Age demonstrated strong correlation with both supratentorial atrophy patterns. There was no correlation between ICAC and cortical atrophy. There was correlation however between central atrophy and ICAC. This persisted even after adjustment for age.

Conclusion

Age is the most important determinant of atrophic cerebral changes. However, high ICAC demonstrated age independent association with central atrophy.

     

Specific localization of scallop gill epithelial calmodulin in cilia

Calmodulin has been isolated and characterized from the gill of the bay scallop aequipecten irradians. Quantitative electrophoretic analysis of epithelial cell fractions show most of the calmodulin to be localized in the cilia, specifically in the detergent- solubilized membrane-matrix fraction. Calmodulin represents 2.2 +/- 0.3 percent of the membrane-matrix protein or 0.41 +/- 0.5 percent of the total ciliary protein. Its concentration is at least 10(-4) M if distributed uniformly within the matrix. Extraction in the presence of calcium suggests that the calmodulin is not bound to the axoneme proper. The ciliary protein is identified as a calmodulin on the basis of its calcium- dependent binding to a fluphenazine-sepharose affinity column and its comigration with bovine brain calmodulin on alkaline-urea and SDS polyacrylamide gels in both the presence and absence of calcium. Scallop ciliary calmodulin activates bovine brain phosphodiesterase to the same extent as bovine brain and chicken gizzard calmodulins. Containing trimethyllysine and lacking cysteine and tryptophan, the amino acid composition of gill calmodulin is typical of known calmodulins, except that it is relatively high in serine and low in methionine. Its composition is less acidic than other calmodulins, in agreement with an observed isoelectric point approximately 0.2 units higher than that of bovine brain. Comparative tryptic peptide mapping of scallop gill ciliary and bovine brain calmodulins indicates coincidence of over 75 percent of the major peptides, but at least two major peptides in each show no near-equivalency. Preliminary results using ATP-reactivated gill cell models show no effect of calcium at micromolar levels on ciliary beat or directionality of the lateral cilia, the cilia which constitute the vast majority of those isolated. However, ciliary arrest will occur at calcium levels more than 150 muM. Because calmodulin usually functions in the micromolar range, its role in this system is unclear. Scallop gill ciliary calmodulin may be involved in the direct regulation of dyneintubule sliding, or it may serve some coupled calcium transport function. At the concentration in which it is found, it must also at least act as a calcium buffer.