Cluster-based network model for time-course gene expression data

We propose a model-based approach to unify clustering and network modeling using time-course gene expression data. Specifically, our approach uses a mixture model to cluster genes. Genes within the same cluster share a similar expression profile. The network is built over cluster-specific expression profiles using state-space models. We discuss the application of our model to simulated data as well as to time-course gene expression data arising from animal models on prostate cancer progression. The latter application shows that with a combined statistical/bioinformatics analyses, we are able to extract gene-to-gene relationships supported by the literature as well as new plausible relationships.

Time course of foreign gene integration and expression in transgenic fish embryos

Using a nuclear transplantation approach, the integration and expression of the green fluorescent protein (GFP) gene in the embryogenesis of transgenic leach (Misgurnus anguillicaudatus Cantor) have been studied. The GFP gene expression is first observed at the gastrula stage, which is consistent with the initiation of cell differentiation of fish embryos. The time course of the foreign gene expression is correlated with the regulatory sequences. The expression efficiency also depends on the gene configuration: the expression of pre-integrating circular plasmid at early embryos is higher than that of the linear plasmid. The integration of the GFP gene is first detected at the blastula stage and lasts for quite a long period. When two types of different plasmids are co-injected into fertilized eggs, the behavior of their integration and expression is not identical.

Time course of training-induced microcirculatory changes and of vegf expression in skeletal muscles of spontaneously hypertensive female rats

Exercise-induced vessel changes modulate arterial pressure (AP) in male spontaneously hypertensive rats (SHR). Vascular endothelial growth factor (VEGF) is important for angiogenesis of skeletal muscle. The present study evaluated the time course of VEGF and angiogenesis after short- and long-term exercise training of female SHR and Wistar Kyoto (WKY) rats, 8-9 weeks (200-250 g). Rats were allocated to daily training or remained sedentary for 3 days (N = 23) or 13 weeks (N = 23). After training, the carotid artery was catheterized for AP measurements. Locomotor (tibialis anterior and gracilis) and non-locomotor skeletal muscles (temporalis) were harvested and prepared for histologic and protein expression analyses. Training increased treadmill performance by all groups (SHR = 28%, WKY = 64%, 3 days) and (SHR = 141%, WKY = 122%, 13 weeks). SHR had higher values of AP than WKY (174 ± 4 vs 111 ± 2 mmHg) that were not altered by training. Three days of running increased VEGF expression (SHR = 28%, WKY = 36%) simultaneously with an increase in capillary-to-fiber ratio in gracilis muscle (SHR = 19%, WKY = 15%). In contrast, 13 weeks of training increased gracilis capillary-to-fiber ratio (SHR = 18%, WKY = 19%), without simultaneous changes in VEGF expression. Training did not change VEGF expression and capillarity of temporalis muscle. We conclude that training stimulates time- and tissue-dependent VEGF protein expression, independent of pressure levels. VEGF triggers angiogenesis in locomotor skeletal muscle shortly after the exercise starts, but is not involved in the maintenance of capillarity after long-term exercise in female rats.

Time course-dependent changes in the transcriptome of human skeletal muscle during recovery from endurance exercise: From inflammation to adaptive remodeling

Re-programming of gene expression is fundamental for skeletal muscle adaptations in response to endurance exercise. This study investigated the time-course dependent changes in the muscular transcriptome following an endurance exercise trial consisting of 1 h of intense cycling immediately followed by 1 h of intense running. Skeletal muscle samples were taken at baseline, 3 h, 48 h, and 96 h post-exercise from eight healthy, endurance-trained, male individuals. RNA was extracted from muscle. Differential gene expression was evaluated using Illumina microarrays and validated with qPCR. Gene set enrichment analysis identified enriched molecular signatures chosen from the Molecular Signatures Database. Three h post-exercise, 102 gene sets were up-regulated [family wise error rate (FWER), P < 0.05]; including groups of genes related with leukocyte migration, immune and chaperone activation, and cyclic AMP responsive element binding protein (CREB) 1-signaling. Forty-eight h post-exercise, among 19 enriched gene sets (FWER, P < 0.05), two gene sets related to actin cytoskeleton remodeling were up-regulated. Ninety-six h post-exercise, 83 gene sets were enriched (FWER, P < 0.05), 80 of which were up-regulated; including gene groups related to chemokine signaling, cell stress management, and extracellular matrix remodeling. These data provide comprehensive insights into the molecular pathways involved in acute stress, recovery, and adaptive muscular responses to endurance exercise. The novel 96 h post-exercise transcriptome indicates substantial transcriptional activity, potentially associated with the prolonged presence of leukocytes in the muscles. This suggests that muscular recovery, from a transcriptional perspective, is incomplete 96 h after endurance exercise involving muscle damage.

Time course analysis of hypoxia, granulation tissue and blood vessel growth, and remodeling in healing rat cutaneous incisional primary intention wounds

Hypoxia and the development and remodeling of blood vessels and connective tissue in granulation tissue that forms in a wound gap following full-thickness skin incision in the rat were examined as a function of time. A 1.5 cm-long incisional wound was created in rat groin skin and the opposed edges sutured together. Wounds were harvested between 3 days and 16 weeks and hypoxia, percent vascular volume, cell proliferation and apoptosis, α-smooth muscle actin, vascular endothelial growth factor-A, vascular endothelial growth factor receptor-2, and transforming growth factor-β 1 expression in granulation tissue were then assessed. Hypoxia was evident between 3 and 7 days while maximal cell proliferation at 3 days (123.6 ± 22.2 cells/mm 2, p < 0.001 when compared with normal skin) preceded the peak percent vascular volume that occurred at 7 days (15.83 ± 1.10%, p < 0.001 when compared with normal skin). The peak in cell apoptosis occurred at 3 weeks (12.1 ± 1.3 cells/mm 2, p < 0.001 when compared with normal skin). Intense α-smooth muscle actin labeling in myofibroblasts was evident at 7 and 10 days. Vascular endothelial growth factor receptor-2 and vascular endothelial growth factor-A were detectable until 2 and 3 weeks, respectively, while transforming growth factor-β 1 protein was detectable in endothelial cells and myofibroblasts until 3-4 weeks and in the extracellular matrix for 16 weeks. Incisional wound granulation tissue largely developed within 3-7 days in the presence of hypoxia. Remodeling, marked by a decline in the percent vascular volume and increased cellular apoptosis, occurred largely in the absence of detectable hypoxia. The expression of vascular endothelial growth factor-A, vascular endothelial growth factor receptor-2, and transforming growth factor-β 1 is evident prior, during, and after the peak of vascular volume reflecting multiple roles for these factors during wound healing.

Robust test method for time-course microarray experiments.

BACKGROUND: In a time-course microarray experiment, the expression level for each gene is observed across a number of time-points in order to characterize the temporal trajectories of the gene-expression profiles. For many of these experiments, the scientific aim is the identification of genes for which the trajectories depend on an experimental or phenotypic factor. There is an extensive recent body of literature on statistical methodology for addressing this analytical problem. Most of the existing methods are based on estimating the time-course trajectories using parametric or non-parametric mean regression methods. The sensitivity of these regression methods to outliers, an issue that is well documented in the statistical literature, should be of concern when analyzing microarray data. RESULTS: In this paper, we propose a robust testing method for identifying genes whose expression time profiles depend on a factor. Furthermore, we propose a multiple testing procedure to adjust for multiplicity. CONCLUSIONS: Through an extensive simulation study, we will illustrate the performance of our method. Finally, we will report the results from applying our method to a case study and discussing potential extensions.

Time-dependent expression of annexin 1 in a model of chronic granulomatous inflammation

Objective and design: To determine the expression pattern and distribution of the glucocorticoid-inducible protein annexin 1 (ANXA1) in a murine model of chronic granulomatous inflammation.Materials or subjects: TO Mouse.Treatment: Chronic granulomatous inflammation was induced by injecting into dorsal sub-cutaneous air-pouches in mice, a mixture of croton oil and Freund's complete adjuvant (CO/FCA).Methods: Western and northern analysis, corticosterone assay, and immunohistochemistry. Statistical analysis was performed using ANOVA followed by Tukey's pair-wise comparisons or Dunnett's multiple comparisons.Results: ANXA1 protein levels changed significantly throughout the 4-week time course, with an initial peak at day 7 and a later elevation at 28 days. ANXA1 mRNA levels peaked at days 1 and 3, with a significant decline at day 7 followed by an upward trend to day 28. Plasma corticosterone measurements taken throughout the time course revealed an increase from 14 days onward, suggesting that corticosterone does not influence ANXA1 expression during the initial stages of the model. Immunogold staining revealed that ANXA1 expression in the inflamed tissue was mainly in extravasated neutrophils, with intact protein (37 kDa) being predominantly observed on the cell membrane.Conclusions: the pattern of ANXA1 expression indicates that infiltrated neutrophils are responsible for the majority of ANXA1 present both at early and later stages of this model of granulomatous inflammation.

Time-course morphological and functional disorders of the kidney induced by long-term high-fat diet intake in female rats

BackgroundEvidence is emerging that highlights the far-reaching consequences of a high-fat diet (HFD) on kidney morphology and function disorders.MethodsThe present study was performed on 3-, 5-, 7- and 9-week-old HFD female rats compared with the appropriate gender and age-matched animals. We evaluated the kidney expression of angiotensin type II receptor and fibrotic and epithelial-to-mesenchymal transition (EMT) markers, by immunoblotting and immunohistochemical and histological techniques, in parallel with kidney function.ResultsIn the current study, the time-course HFD-treated group showed, by immunoblotting and immunohistochemical analysis, an early time-course increase in the expression of transforming growth factor β-1 (TGFβ-1) in the entire kidney of HFD-treated rats, compared with that observed in the control group. Simultaneously, the study shows a transient increase in the expression of ZEB2 in the HFD whole kidney accompanied by a fall in the E-cadherin expression and increased collagen and fibronectin deposition. A pronounced decrease in fractional urinary sodium excretion was also demonstrated in the long-term HFD-treated rats. The decreased FENa + was accompanied by a fall in FEPNa + and FEPP Na +, which occurred in association with significantly decreased CCr and, certainly on the sodium-filtered load. The reduction in the glomerular filtration rate (GFR) occurred in parallel to proteinuria and glomerular desmin overexpression.ConclusionsThe results of the current study suggest that podocyte injury in parallel with observed proteinuria and evidence of EMT transformation are associated with long-term loss of kidney function and renal sodium and water retention. © 2013 The Author. All rights reserved.

Time course involvement of matrix metalloproteinases in the vascular alterations of renovascular hypertension

Increased vascular matrix metalloproteinases (MMPs) levels play a role in late phases of hypertensive vascular remodeling. However, no previous study has examined the time course of MMPs in the various phases of two-kidney, one-clip hypertension (2K1C). We examined structural vascular changes, collagen and elastin content, vascular oxidative stress, and MMPs levels/activities during the development of 2K1C hypertension. Plasma angiotensin converting enzyme (ACE) activity was measured to assess renin-angiotensin system activation. Sham or 2K1C hypertensive rats were studied after 2, 4, 6, and 10 weeks of hypertension. Systolic blood pressure (SBP) was monitored weekly. Morphometry of structural changes in the aortic wall was studied in hematoxylin/eosin, orcein and picrosirius red sections. Aortic NADPH activity and superoxide production was evaluated. Aortic gelatinolytic activity was determined by in situ zymography, and MMP-2, MMP-14, and tissue inhibitor of MMPs (TIMP)-2 levels were determined by gelatin zymography, immunofluorescence and immunohistochemistry. 2K1C hypertension was associated with increased ACE activity, which decreased to normal after 10 weeks. We found increased aortic collagen and elastin content in the early phase of hypertension, which were associated with vascular hypertrophy, increased vascular MMP-2 and MMP-14 (but not TIMP-2) levels, and increased gelatinolytic activity, possibly as a result of increased vascular NADPH oxidase activity and oxidative stress. These results indicate that vascular remodeling of renovascular hypertension is an early process associated with early increases in MMPs activities, enhanced matrix deposition and oxidative stress. Using antioxidants or MMPs inhibitors in the early phase of hypertension may prevent the vascular alterations of hypertension. (C) 2012 Elsevier B.V. All rights reserved.

Time course of hydrogen peroxide-thioredoxin balance and its influence on the intracellular signalling in myocardial infarction

We investigated the myocardial thioredoxin-1 and hydrogen peroxide concentrations and their association with some prosurvival and pro-apoptotic proteins, during the transition from myocardial infarction (MI) to heart failure in rats. Male Wistar rats were divided into the following six groups: three sham-operated groups and three MI groups, each at at 2, 7 and 28 days postsurgery. Cardiac function was analysed by echocardiography; the concentration of H2O2 and the ratio of reduced to oxidized glutathione were measured spectrophotometrically, while the myocardial immunocontent of thioredoxin-1, angiotensin II, angiotensin II type 1 and type 2 receptors, p-JNK/JNK, p-ERK/ERK, p-Akt/Akt, p-mTOR/mTOR and p-GSK3 beta/GSK3 beta was evaluated by Western blot. Our results show that thioredoxin-1 appears to make an important contribution to the reduced H2O2 concentration. It was associated with lower JNK expression in the early period post-MI (2 days). However, thioredoxin-1 decreased, while reninangiotensin system markers and levels of H2O2 increased, over 28 days post-MI, in parallel with some signalling proteins involved in maladaptative cardiac remodelling and ventricular dysfunction. These findings provide insight into the time course profile of endogenous antioxidant adaptation to ischaemic injury, which may be useful for the design of therapeutical strategies targeting oxidative stress post-MI.

Time course of temperature effects on cardiovascular mortality in Brisbane, Australia

Objective To quantify the lagged effects of mean temperature on deaths from cardiovascular diseases in Brisbane, Australia. Design Polynomial distributed lag models were used to assess the percentage increase in mortality up to 30 days associated with an increase (or decrease) of 1°C above (or below) the threshold temperature. Setting Brisbane, Australia. Patients 22 805 cardiovascular deaths registered between 1996 and 2004. Main outcome measures Deaths from cardiovascular diseases. Results The results show a longer lagged effect in cold days and a shorter lagged effect in hot days. For the hot effect, a statistically significant association was observed only for lag 0–1 days. The percentage increase in mortality was found to be 3.7% (95% CI 0.4% to 7.1%) for people aged ≥65 years and 3.5% (95% CI 0.4% to 6.7%) for all ages associated with an increase of 1°C above the threshold temperature of 24°C. For the cold effect, a significant effect of temperature was found for 10–15 lag days. The percentage estimates for older people and all ages were 3.1% (95% CI 0.7% to 5.7%) and 2.8% (95% CI 0.5% to 5.1%), respectively, with a decrease of 1°C below the threshold temperature of 24°C. Conclusions The lagged effects lasted longer for cold temperatures but were apparently shorter for hot temperatures. There was no substantial difference in the lag effect of temperature on mortality between all ages and those aged ≥65 years in Brisbane, Australia.

Event-related brain potentials dissociate the developmental time-course of automatic numerical magnitude analysis and cognitive control functions during the first three years of primary school

In this study we set out to dissociate the developmental time course of automatic symbolic number processing and cognitive control functions in grade 1-3 British primary school children. Event-related potential (ERP) and behavioral data were collected in a physical size discrimination numerical Stroop task. Task-irrelevant numerical information was processed automatically already in grade 1. Weakening interference and strengthening facilitation indicated the parallel development of general cognitive control and automatic number processing. Relationships among ERP and behavioral effects suggest that control functions play a larger role in younger children and that automaticity of number processing increases from grade 1 to 3.

The time course of in vivo recovery of transverse strain in high-stress tendons following exercise

Objective To evaluate the time course of the recovery of transverse strain in the Achilles and patellar tendon following a bout of resistance exercise. Methods Seventeen healthy adults underwent sonographic examination of the right patellar (n=9) and Achilles (n=8) tendons immediately prior to and following 90 repetitions of weight-bearing quadriceps and gastrocnemius-resistance exercise performed against an effective resistance of 175% and 250% body weight, respectively. Sagittal tendon thickness was determined 20 mm from the enthesis and transverse strain, as defined by the stretch ratio, was repeatedly monitored over a 24 h recovery period. Results Resistance exercise resulted in an immediate decrease in Achilles (t7=10.6, p<0.01) and patellar (t8=8.9, p<0.01) tendon thickness, resulting in an average transverse stretch ratio of 0.86±0.04 and 0.82±0.05, which was not significantly different between tendons. The magnitude of the immediate transverse strain response, however, was reduced with advancing age (r=0.63, p<0.01). Recovery in transverse strain was prolonged compared with the duration of loading and exponential in nature. The average primary recovery time was not significantly different between the Achilles (6.5±3.2 h) and patellar (7.1±3.2 h) tendons. Body weight accounted for 62% and 64% of the variation in recovery time, respectively. Conclusions Despite structural and biochemical differences between the Achilles and patellar tendon, the mechanisms underlying transverse creep recovery in vivo appear similar and are highly time dependent. These novel findings have important implications concerning the time required for the mechanical recovery of high-stress tendons following an acute bout of exercise.

Early time course of Akt phosphorylation after endurance and resistance exercise

Purpose The aim of this study was to determine the early time course of exercise-induced signaling after divergent contractile activity associated with resistance and endurance exercise. Methods Sixteen male subjects were randomly assigned to either a cycling (CYC; n = 8, 60 min, 70% V?O2peak) or resistance (REX; n = 8, 8×5 leg extension, 80% one-repetition maximum, 3-min recovery) exercise group. Serial muscle biopsies were obtained from vastus lateralis at rest before, immediately after, and after 15, 30, and 60 min of passive recovery to determine early signaling responses after exercise. Results There were comparable increases from rest in AktThr308/Ser473 and mTORSer2448 phosphorylation during the postexercise time course that peaked 30-60 min after both CYC and REX (P<0.05). There were also similar patterns in p70S6K Thr389 and 4E-BP1Thr37/46 phosphorylation, but a greater magnitude of effect was observed for REX and CYC, respectively (P<0.05). However, AMPKThr172 phosphorylation was only significantly elevated after CYC (P<0.05), and we observed divergent responses for glycogen synthaseSer641 and AS160 phosphorylation that were enhanced after CYC but not REX (P<0.05). Conclusions We show a similar time course for Akt-mTOR-S6K phosphorylation during the initial 60-min recovery period after divergent contractile stimuli. Conversely, enhanced phosphorylation status of proteins that promote glucose transport and glycogen synthesis only occurred after endurance exercise. Our results indicate that endurance and resistance exercise initiate translational signaling, but high-load, low-repetition contractile activity failed to promote phosphorylation of pathways regulating glucose metabolism.