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.

     

Sex differences in provisioning rules and honest signalling of need in Manx shearwaters, Puffinus puffinus

Sex differences in provisioning rules have been found in a number of bird species, but the reasons remain unclear. Studies of begging in species with single-chick broods exclude the influence of nestling competition and may provide especially useful models for the study of signalling during parent-offspring conflict. We tested whether sex differences in provisioning rules occur in a species with an obligate brood size of one, the Manx shearwater. We found that chicks conveyed information about their body condition through begging, but male and female parents responded differently to that information. Females varied meal sizes according to the begging intensity of the chick and adjusted subsequent trip duration according to the chick's body condition after feeding, but males did not alter meal size or adjust trip duration. We discuss these findings in the context of recorded differences in the contributions to food provisioning by male and female parents, and we discuss why females may respond more sensitively than males to changes in chicks' nutritional requirements.     

Genetic monogamy in burrowing parrots Cyanoliseus patagonus?

We report on a first DNA fingerprinting study of paternity in a Psittaciform bird, the burrowing parrot Cyanoliseus patagonus . In two consecutive breeding seasons, a total of 49 families was sampled, of which 11 breeding pairs were investigated each of two years. Extra-pair paternity was not encountered suggesting that burrowing parrots are socially as well as genetically monogamous. Strict genetic monogamy is comparatively rare in birds and occurs predominantly in some groups of non-passeriformes all of which exhibit long reproductive lifespans and essential paternal care. Psittaciformes fit this pattern. We conclude that paternal care plays a crucial role in the evolution and maintenance of genetic monogamy in the study species. Cases of intraspecific brood parasitism are rarely observed.     

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.     

Early events in the cellular formation of proparathyroid hormone

Early events in the cellular synthesis and subsequent transfer into membrane-limited compartments of pre-proparathyroid hormone (pre-proPTH) and proparathyroid hormone (proPTH) were investigated by electrophoretic analyses of newly synthesized proteins in subcellular fractions of parthyroid gland slices pulse-labeled for 0.5-5 min with [(35)S] methionine. During these short times of incubation, both pre-proPTH and proPTH were confined to the microsomal fraction. Labeled pre-proPTH and proPTH were detected in a 30-s interval between 0.5 and 1.0 min of incubation. The radioactivity in proPTH became relatively constant between 3 and 5 min, whereas the radioactivity in ProPTH increased markedly over this period. When corrected for the known content of methionine in the prohormone and the prohormone, we found four times as much radiolabeled prohormone as prehormone between 0.5 and 1.0 min of synthesis. Sequestration of labeled prohomrone into endoplasmic reticulum compartments was shown by treatment of the microsomal fraction with chymotrypsin and trypsin, which resulted in the degradation of the prehormone but not of the prohormones. Approximately 50 percent of pre-prohormone and 25 percent of prohormone were released from the microsomes by their extraction with 1.0 M KCl, whereas 80-90 percent of both was released by treatment with Triton X-100. These results in intact cells support the signal hypothesis proposed by Blobel and his co-workers in studies utilizing cell-free systems, inasmuch as the results indicate transfer of prohormone into the cisternal space of the rough endoplasmic reticulum concomitant with the growth of the nascent polypeptide chain. Appearance of membrane-sequestered proPTH takes place without entry of pre-proPTH into the cisternal space, suggesting that proteolytic removal of the leader peptide occurs during transfer of the polypeptide through the lipid bilayer. Further evidence in support of this process is that pre-proPTH is only partly extracted from the microsomes by treatment with 1.0 M KCl, suggesting that a substantial fraction of the nascent pre-proPTH is integrally inserted into the membranes before it is cleaved to form proPTH.     

Phylogenetic relationships within the Alcidae (Charadriiformes: Aves) inferred from total molecular evidence

The Alcidae is a unique assemblage of Northern Hemisphere seabirds that forage by "flying" underwater. Despite obvious affinities among the species, their evolutionary relationships are unclear. We analyzed nucleotide sequences of 1,045 base pairs of the mitochondrial cytochrome b gene and allelic profiles for 37 allozyme loci in all 22 extant species. Trees were constructed on independent and combined data sets using maximum parsimony and distance methods that correct for superimposed changes. Alternative methods of analysis produced only minor differences in relationships that were supported strongly by bootstrapping or standard error tests. Combining sequence and allozyme data into a single analysis provided the greatest number of relationships receiving strong support. Addition of published morphological and ecological data did not improve support for any additional relationship. All analyses grouped species into six distinct lineages: (1) the dovekie (Alle alle) and auks, (2) guillemots, (3) brachyramphine murrelets, (4) synthliboramphine murrelets, (5) true auklets, and (6) the rhinoceros auklet (Cerorhinca monocerata) and puffins. The two murres (genus Uria) were sister taxa, and the black guillemot (Cepphus grylle) was basal to the other guillemots. The Asian subspecies of the marbled murrelet (Brachyramphus marmoratus perdix) was the most divergent brachyramphine murrelet, and two distinct lineages occurred within the synthliboramphine murrelets. Cassin's auklet (Ptychoramphus aleuticus) and the rhinoceros auklet were basal to the other auklets and puffins, respectively, and the Atlantic (Fratercula arctica) and horned (Fratercula corniculata) puffins were sister taxa. Several relationships among tribes, among the dovekie and auks, and among the auklets could not be resolved but resembled "star" phylogenies indicative of adaptive radiations at different depths within the trees.      

Evidence for para Dechlorination of Polychlorobiphenyls by Methanogenic Bacteria

When microorganisms eluted from upper Hudson River sediment were cultured without any substrate except polychlorobiphenyl (PCB)-free Hudson River sediment, methane formation was the terminal step of the anaerobic food chain. In sediments containing Aroclor 1242, addition of eubacterium-inhibiting antibiotics, which should have directly inhibited fermentative bacteria and thereby should have indirectly inhibited methanogens, resulted in no dechlorination activity or methane production. However, when substrates for methanogenic bacteria were provided along with the antibiotics (to free the methanogens from dependence on eubacteria), concomitant methane production and dechlorination of PCBs were observed. The dechlorination of Aroclor 1242 was from the para positions, a pattern distinctly different from, and more limited than, the pattern observed with untreated or pasteurized inocula. Both methane production and dechlorination in cultures amended with antibiotics plus methanogenic substrates were inhibited by 2-bromoethanesulfonic acid. These results suggest that the methanogenic bacteria are among the physiological groups capable of anaerobic dechlorination of PCBs, but that the dechlorination observed with methanogenic bacteria is less extensive than the dechlorination observed with more complex anaerobic consortia.