Relatedness and HLA-DRB1 typing may discriminate the magnitude of the genetic susceptibility to tuberculosis using a household contact model

Background Tuberculosis clusters in families may be due to increased household exposure, shared genetic factors, or both. Household contact studies are useful to control exposure because socioeconomic and environmental conditions are similar to all subjects, allowing the evaluation of the contribution of relatedness to disease development. Methods In this study, the familial aggregation of tuberculosis using relatedness and a specific inherited marker (HLA-DRB1) was evaluated. Fifty families, which had at least two cases of tuberculosis diagnosed within the past 5 years, were selected from a cohort of tuberculosis carried out in Recife, Brazil. The first case diagnosed was considered to be a primary case. The secondary attack rate of tuberculosis in household contacts was estimated according to the degree of relatedness. The relative risk of having tuberculosis based on the degree of relatedness household and the population attributable fraction to relatedness were also estimated. HLA-DRB1 typing and attributable etiologic/preventive fractions were calculated among sick and healthy household contacts. Results Compared to unrelated contacts, the relative risk for tuberculosis adjusted for age was 1.38 (95% CI 0.86 to 2.21). Relatedness contributed 23% to the development of tuberculosis at the population levels. The HLA-DRB1*04 allele group (OR = 2.44; p =0.0324; etiologic fraction =0.15) was overrepresented and the DRB1*15 allele group (OR=0.48; p=0.0488; protective fraction=0.19) was underrepresented among household contacts exhibiting tuberculosis. The presence of DRB1 shared alleles between primary cases and their contacts was a risk factor for tuberculosis (p=0.0281). Conclusion This household contact model together with the utilisation of two genetic variables permitted the evaluation of genetic factors contributing towards tuberculosis development.

Chaotic model and memory in single calcium-activated potassium channel kinetics

Ion channels are pores formed by proteins and responsible for carrying ion fluxes through cellular membranes. The ion channels can assume conformational states thereby controlling ion flow. Physically, the conformational transitions from one state to another are associated with energy barriers between them and are dependent on stimulus, such as, electrical field, ligands, second messengers, etc. Several models have been proposed to describe the kinetics of ion channels. The classical Markovian model assumes that a future transition is independent of the time that the ion channel stayed in a previous state. Others models as the fractal and the chaotic assume that the rate of transitions between the states depend on the time that the ionic channel stayed in a previous state. For the calcium activated potassium channels of Leydig cells the R/S Hurst analysis has indicated that the channels are long-term correlated with a Hurst coefficient H around 0.7, showing a persistent memory in this kinetic. Here, we applied the R/S analysis to the opening and closing dwell time series obtained from simulated data from a chaotic model proposed by L. Liebovitch and T. Toth [J. Theor. Biol. 148, 243 (1991)] and we show that this chaotic model or any model that treats the set of channel openings and closings as independent events is inadequate to describe the long-term correlation (memory) already described for the experimental data. (C) 2008 American Institute of Physics.

Ultra-Deep Sequencing Reveals the microRNA Expression Pattern of the Human Stomach

Background: While microRNAs (miRNAs) play important roles in tissue differentiation and in maintaining basal physiology, little is known about the miRNA expression levels in stomach tissue. Alterations in the miRNA profile can lead to cell deregulation, which can induce neoplasia. Methodology/Principal Findings: A small RNA library of stomach tissue was sequenced using high-throughput SOLiD sequencing technology. We obtained 261,274 quality reads with perfect matches to the human miRnome, and 42% of known miRNAs were identified. Digital Gene Expression profiling (DGE) was performed based on read abundance and showed that fifteen miRNAs were highly expressed in gastric tissue. Subsequently, the expression of these miRNAs was validated in 10 healthy individuals by RT-PCR showed a significant correlation of 83.97% (P<0.05). Six miRNAs showed a low variable pattern of expression (miR-29b, miR-29c, miR-19b, miR-31, miR-148a, miR-451) and could be considered part of the expression pattern of the healthy gastric tissue. Conclusions/Significance: This study aimed to validate normal miRNA profiles of human gastric tissue to establish a reference profile for healthy individuals. Determining the regulatory processes acting in the stomach will be important in the fight against gastric cancer, which is the second-leading cause of cancer mortality worldwide.

Tumor slices as a model to evaluate doxorubicin in vitro treatment and expression of trios of genes PRSS11, MTSS1, CLPTM1 and PRSS11, MTSS1, SMYD2 in canine mammary gland cancer

Background: In women with breast cancer submitted to neoadjuvant chemotherapy based in doxorubicin, tumor expression of groups of three genes (PRSS11, MTSS1, CLPTM1 and PRSS11, MTSS1, SMYD2) have classified them as responsive or resistant. We have investigated whether expression of these trios of genes could predict mammary carcinoma response in dogs and whether tumor slices, which maintain epithelial-mesenchymal interactions, could be used to evaluate drug response in vitro. Methods: Tumors from 38 dogs were sliced and cultured with or without doxorubicin 1 mu M for 24 h. Tumor cells were counted by two observers to establish a percentage variation in cell number, between slices. Based on these results, a reduction in cell number between treated and control samples >= 21.7%, arbitrarily classified samples, as drug responsive. Tumor expression of PRSS11, MTSS1, CLPTM1 and SMYD2, was evaluated by real time PCR. Relative expression results were then transformed to their natural logarithm values, which were spatially disposed according to the expression of trios of genes, comprising PRSS11, MTSS1, CLPTM1 and PRSS11, MTSS1, SMYD2. Fisher linear discrimination test was used to generate a separation plane between responsive and non-responsive tumors. Results: Culture of tumor slices for 24 h was feasible. Nine samples were considered responsive and 29 non-responsive to doxorubicin, considering the pre-established cut-off value of cell number reduction = 21.7%, between doxorubicin treated and control samples. Relative gene expression was evaluated and tumor samples were then spatially distributed according to the expression of the trios of genes: PRSS11, MTSS1, CLPTM1 and PRSS11, MTSS1, SMYD2. A separation plane was generated. However, no clear separation between responsive and non-responsive samples could be observed. Conclusion: Three-dimensional distribution of samples according to the expression of the trios of genes PRSS11, MTSS1, CLPTM1 and PRSS11, MTSS1, SMYD2 could not predict doxorubicin in vitro responsiveness. Short term culture of mammary gland cancer slices may be an interesting model to evaluate chemotherapy activity.

A New Mouse Model for Marfan Syndrome Presents Phenotypic Variability Associated with the Genetic Background and Overall Levels of Fbn1 Expression

Marfan syndrome is an autosomal dominant disease of connective tissue caused by mutations in the fibrillin-1 encoding gene FBN1. Patients present cardiovascular, ocular and skeletal manifestations, and although being fully penetrant, MFS is characterized by a wide clinical variability both within and between families. Here we describe a new mouse model of MFS that recapitulates the clinical heterogeneity of the syndrome in humans. Heterozygotes for the mutant Fbn1 allele mg Delta(loxPneo), carrying the same internal deletion of exons 19-24 as the mg Delta mouse model, present defective microfibrillar deposition, emphysema, deterioration of aortic wall and kyphosis. However, the onset of a clinical phenotypes is earlier in the 129/Sv than in C57BL/6 background, indicating the existence of genetic modifiers of MFS between these two mouse strains. In addition, we characterized a wide clinical variability within the 129/Sv congenic heterozygotes, suggesting involvement of epigenetic factors in disease severity. Finally, we show a strong negative correlation between overall levels of Fbn1 expression and the severity of the phenotypes, corroborating the suggested protective role of normal fibrillin-1 in MFS pathogenesis, and supporting the development of therapies based on increasing Fbn1 expression.

Simplified Mathematical Model for an Anaerobic Sequencing Batch Biofilm Reactor Treating Lipid-Rich Wastewater Subject to Rising Organic Loading Rates

This study proposes a simplified mathematical model to describe the processes occurring in an anaerobic sequencing batch biofilm reactor (ASBBR) treating lipid-rich wastewater. The reactor, subjected to rising organic loading rates, contained biomass immobilized cubic polyurethane foam matrices, and was operated at 32 degrees C +/- 2 degrees C, using 24-h batch cycles. In the adaptation period, the reactor was fed with synthetic substrate for 46 days and was operated without agitation. Whereas agitation was raised to 500 rpm, the organic loading rate (OLR) rose from 0.3 g chemical oxygen demand (COD) . L(-1) . day(-1) to 1.2 g COD . L(-1) . day(-1). The ASBBR was fed fat-rich wastewater (dairy wastewater), in an operation period lasting for 116 days, during which four operational conditions (OCs) were tested: 1.1 +/- 0.2 g COD . L(-1) . day(-1) (OC1), 4.5 +/- 0.4 g COD . L(-1) . day(-1) (OC2), 8.0 +/- 0.8 g COD . L(-1) . day(-1) (OC3), and 12.1 +/- 2.4 g COD . L(-1) . day(-1) (OC4). The bicarbonate alkalinity (BA)/COD supplementation ratio was 1:1 at OC1, 1:2 at OC2, and 1:3 at OC3 and OC4. Total COD removal efficiencies were higher than 90%, with a constant production of bicarbonate alkalinity, in all OCs tested. After the process reached stability, temporal profiles of substrate consumption were obtained. Based on these experimental data a simplified first-order model was fit, making possible the inference of kinetic parameters. A simplified mathematical model correlating soluble COD with volatile fatty acids (VFA) was also proposed, and through it the consumption rates of intermediate products as propionic and acetic acid were inferred. Results showed that the microbial consortium worked properly and high efficiencies were obtained, even with high initial substrate concentrations, which led to the accumulation of intermediate metabolites and caused low specific consumption rates.

The Coupled Aerosol and Tracer Transport model to the Brazilian developments on the Regional Atmospheric Modeling System (CATT-BRAMS) - Part 1: Model description and evaluation

We introduce the Coupled Aerosol and Tracer Transport model to the Brazilian developments on the Regional Atmospheric Modeling System (CATT-BRAMS). CATT-BRAMS is an on-line transport model fully consistent with the simulated atmospheric dynamics. Emission sources from biomass burning and urban-industrial-vehicular activities for trace gases and from biomass burning aerosol particles are obtained from several published datasets and remote sensing information. The tracer and aerosol mass concentration prognostics include the effects of sub-grid scale turbulence in the planetary boundary layer, convective transport by shallow and deep moist convection, wet and dry deposition, and plume rise associated with vegetation fires in addition to the grid scale transport. The radiation parameterization takes into account the interaction between the simulated biomass burning aerosol particles and short and long wave radiation. The atmospheric model BRAMS is based on the Regional Atmospheric Modeling System (RAMS), with several improvements associated with cumulus convection representation, soil moisture initialization and surface scheme tuned for the tropics, among others. In this paper the CATT-BRAMS model is used to simulate carbon monoxide and particulate material (PM(2.5)) surface fluxes and atmospheric transport during the 2002 LBA field campaigns, conducted during the transition from the dry to wet season in the southwest Amazon Basin. Model evaluation is addressed with comparisons between model results and near surface, radiosondes and airborne measurements performed during the field campaign, as well as remote sensing derived products. We show the matching of emissions strengths to observed carbon monoxide in the LBA campaign. A relatively good comparison to the MOPITT data, in spite of the fact that MOPITT a priori assumptions imply several difficulties, is also obtained.

Turbulent magnetic pumping in a Babcock-Leighton solar dynamo model

Context. The turbulent pumping effect corresponds to the transport of magnetic flux due to the presence of density and turbulence gradients in convectively unstable layers. In the induction equation it appears as an advective term and for this reason it is expected to be important in the solar and stellar dynamo processes. Aims. We explore the effects of turbulent pumping in a flux-dominated Babcock-Leighton solar dynamo model with a solar-like rotation law. Methods. As a first step, only vertical pumping has been considered through the inclusion of a radial diamagnetic term in the induction equation. In the second step, a latitudinal pumping term was included and then, a near-surface shear was included. Results. The results reveal the importance of the pumping mechanism in solving current limitations in mean field dynamo modeling, such as the storage of the magnetic flux and the latitudinal distribution of the sunspots. If a meridional flow is assumed to be present only in the upper part of the convective zone, it is the full turbulent pumping that regulates both the period of the solar cycle and the latitudinal distribution of the sunspot activity. In models that consider shear near the surface, a second shell of toroidal field is generated above r = 0.95 R(circle dot) at all latitudes. If the full pumping is also included, the polar toroidal fields are efficiently advected inwards, and the toroidal magnetic activity survives only at the observed latitudes near the equator. With regard to the parity of the magnetic field, only models that combine turbulent pumping with near-surface shear always converge to the dipolar parity. Conclusions. This result suggests that, under the Babcock-Leighton approach, the equartorward motion of the observed magnetic activity is governed by the latitudinal pumping of the toroidal magnetic field rather than by a large scale coherent meridional flow. Our results support the idea that the parity problem is related to the quadrupolar imprint of the meridional flow on the poloidal component of the magnetic field and the turbulent pumping positively contributes to wash out this imprint.

First stars XIV. Sulfur abundances in extremely metal-poor stars

Context. Precise S abundances are important in the study of the early chemical evolution of the Galaxy. In particular the site of the formation remains uncertain because, at low metallicity, the trend of this alpha-element versus [Fe/H] remains unclear. Moreover, although sulfur is not bound significantly in dust grains in the ISM, it seems to behave differently in DLAs and old metal-poor stars. Aims. We attempt a precise measurement of the S abundance in a sample of extremely metal-poor stars observed with the ESO VLT equipped with UVES, taking into account NLTE and 3D effects. Methods. The NLTE profiles of the lines of multiplet 1 of S I were computed with a version of the program MULTI, including opacity sources from ATLAS9 and based on a new model atom for S. These profiles were fitted to the observed spectra. Results. We find that sulfur in EMP stars behaves like the other alpha-elements, with [S/Fe] remaining approximately constant below [Fe/H] = -3. However, [S/Mg] seems to decrease slightly with increasing [Mg/H]. The overall abundance patterns of O, Na, Mg, Al, S, and K are most closely matched by the SN model yields by Heger & Woosley. The [S/Zn] ratio in EMP stars is solar, as also found in DLAs. We derive an upper limit to the sulfur abundance [S/Fe] < +0.5 for the ultra metal-poor star CS 22949-037. This, along with a previously reported measurement of zinc, argues against the conjecture that the light-element abundance pattern of this star (and by analogy, the hyper iron-poor stars HE 0107-5240 and HE 1327-2326) would be due to dust depletion.

First stars XIII. Two extremely metal-poor RR Lyrae stars

Context. The chemical composition of extremely metal-poor stars (EMP stars; [Fe/H] < similar to -3) is a unique tracer of early nucleosynthesis in the Galaxy. As such stars are rare, we wish to find classes of luminous stars which can be studied at high spectral resolution. Aims. We aim to determine the detailed chemical composition of the two EMP stars CS 30317-056 and CS 22881-039, originally thought to be red horizontal-branch (RHB) stars, and compare it to earlier results for EMP stars as well as to nucleosynthesis yields from various supernova (SN) models. In the analysis, we discovered that our targets are in fact the two most metal-poor RR Lyrae stars known. Methods. Our detailed abundance analysis, taking into account the variability of the stars, is based on VLT/UVES spectra (R similar or equal to 43 000) and 1D LTE OSMARCS model atmospheres and synthetic spectra. For comparison with SN models we also estimate NLTE corrections for a number of elements. Results. We derive LTE abundances for the 16 elements O, Na, Mg, Al, Si, S, Ca, Sc, Ti, Cr, Mn, Fe, Co, Ni, Sr and Ba, in good agreement with earlier values for EMP dwarf, giant and RHB stars. Li and C are not detected in either star. NLTE abundance corrections are newly calculated for O and Mg and taken from the literature for other elements. The resulting abundance pattern is best matched by model yields for supernova explosions with high energy and/or significant asphericity effects. Conclusions. Our results indicate that, except for Li and C, the surface composition of EMP RR Lyr stars is not significantly affected by mass loss, mixing or diffusion processes; hence, EMP RR Lyr stars should also be useful tracers of the chemical evolution of the early Galactic halo. The observed abundance ratios indicate that these stars were born from an ISM polluted by energetic, massive (25-40 M(circle dot)) and/or aspherical supernovae, but the NLTE corrections for Sc and certain other elements do play a role in the choice of model.

uvby-beta photometry of solar twins The solar colors, model atmospheres, and the T(eff) and metallicity scales

Aims. Solar colors have been determined on the uvby-beta photometric system to test absolute solar fluxes, to examine colors predicted by model atmospheres as a function of stellar parameters (T(eff), log g, [Fe/H]), and to probe zero-points of T(eff) and metallicity scales. Methods. New uvby-beta photometry is presented for 73 solar-twin candidates. Most stars of our sample have also been observed spectroscopically to obtain accurate stellar parameters. Using the stars that most closely resemble the Sun, and complementing our data with photometry available in the literature, the solar colors on the uvby-beta system have been inferred. Our solar colors are compared with synthetic solar colors computed from absolute solar spectra and from the latest Kurucz (ATLAS9) and MARCS model atmospheres. The zero-points of different T(eff) and metallicity scales are verified and corrections are proposed. Results. Our solar colors are (b - y)(circle dot) = 0.4105 +/- 0.0015, m(1,circle dot) = 0.2122 +/- 0.0018, c(1,circle dot) = 0.3319 +/- 0.0054, and beta(circle dot) = 2.5915 +/- 0.0024. The (b - y)(circle dot) and m(1,circle dot) colors obtained from absolute spectrophotometry of the Sun agree within 3-sigma with the solar colors derived here when the photometric zero-points are determined from either the STIS HST observations of Vega or an ATLAS9 Vega model, but the c(1,circle dot) and beta(circle dot) synthetic colors inferred from absolute solar spectra agree with our solar colors only when the zero-points based on the ATLAS9 model are adopted. The Kurucz solar model provides a better fit to our observations than the MARCS model. For photometric values computed from the Kurucz models, (b - y)(circle dot) and m(1,circle dot) are in excellent agreement with our solar colors independently of the adopted zero-points, but for c(1,circle dot) and beta circle dot agreement is found only when adopting the ATLAS9 zero-points. The c(1,circle dot) color computed from both the Kurucz and MARCS models is the most discrepant, probably revealing problems either with the models or observations in the u band. The T(eff) calibration of Alonso and collaborators has the poorest performance (similar to 140 K off), while the relation of Casagrande and collaborators is the most accurate (within 10 K). We confirm that the Ramirez & Melendez uvby metallicity calibration, recommended by Arnadottir and collaborators to obtain [Fe/H] in F, G, and K dwarfs, needs a small (similar to 10%) zero-point correction to place the stars and the Sun on the same metallicity scale. Finally, we confirm that the c(1) index in solar analogs has a strong metallicity sensitivity.

Phylogeography of the common vampire bat (Desmodus rotundus): Marked population structure, Neotropical Pleistocene vicariance and incongruence between nuclear and mtDNA markers

Background: The common vampire bat Desmodus rotundus is an excellent model organism for studying ecological vicariance in the Neotropics due to its broad geographic range and its preference for forested areas as roosting sites. With the objective of testing for Pleistocene ecological vicariance, we sequenced a mitocondrial DNA (mtDNA) marker and two nuclear markers (RAG2 and DRB) to try to understand how Pleistocene glaciations affected the distribution of intraspecific lineages in this bat. Results: Five reciprocally monophyletic clades were evident in the mitochondrial gene tree, and in most cases with high bootstrap support: Central America (CA), Amazon and Cerrado (AMC), Pantanal (PAN), Northern Atlantic Forest (NAF) and Southern Atlantic Forest (SAF). The Atlantic forest clades formed a monophyletic clade with high bootstrap support, creating an east/west division for this species in South America. On the one hand, all coalescent and non-coalescent estimates point to a Pleistocene time of divergence between the clades. On the other hand, the nuclear markers showed extensive sharing of haplotypes between distant localities, a result compatible with male-biased gene flow. In order to test if the disparity between the mitochondrial and nuclear markers was due to the difference in mutation rate and effective size, we performed a coalescent simulation to examine the feasibility that, given the time of separation between the observed lineages, even with a gene flow rate close to zero, there would not be reciprocal monophyly for a neutral nuclear marker. We used the observed values of theta and an estimated mutation rate for the nuclear marker gene to perform 1000 iterations of the simulation. The results of this simulation were inconclusive: the number of iterations with and without reciprocal monophyly of one or more clades are similar. Conclusions: We therefore conclude that the pattern exhibited by the common vampire bat, with marked geographical structure for a mitochondrial marker and no phylogeographic structure for nuclear markers is compatible with a historical scenario of complete isolation of refuge-like populations during the Pleistocene. The results on demographic history on this species is compatible with the Carnaval-Moritz model of Pleistocene vicariance, with demographic expansions in the southern Atlantic forest.

Human Dental Pulp Cells: A New Source of Cell Therapy in a Mouse Model of Compressive Spinal Cord Injury

Strategies aimed at improving spinal cord regeneration after trauma are still challenging neurologists and neuroscientists throughout the world. Many cell-based therapies have been tested, with limited success in terms of functional outcome. In this study, we investigated the effects of human dental pulp cells (HDPCs) in a mouse model of compressive spinal cord injury (SCI). These cells present some advantages, such as the ease of the extraction process, and expression of trophic factors and embryonic markers from both ecto-mesenchymal and mesenchymal components. Young adult female C57/BL6 mice were subjected to laminectomy at T9 and compression of the spinal cord with a vascular clip for 1 min. The cells were transplanted 7 days or 28 days after the lesion, in order to compare the recovery when treatment is applied in a subacute or chronic phase. We performed quantitative analyses of white-matter preservation, trophic-factor expression and quantification, and ultrastructural and functional analysis. Our results for the HDPC-transplanted animals showed better white-matter preservation than the DMEM groups, higher levels of trophic-factor expression in the tissue, better tissue organization, and the presence of many axons being myelinated by either Schwann cells or oligodendrocytes, in addition to the presence of some healthy-appearing intact neurons with synapse contacts on their cell bodies. We also demonstrated that HDPCs were able to express some glial markers such as GFAP and S-100. The functional analysis also showed locomotor improvement in these animals. Based on these findings, we propose that HDPCs may be feasible candidates for therapeutic intervention after SCI and central nervous system disorders in humans.

Renal denervation in an animal model of diabetes and hypertension: Impact on the autonomic nervous system and nephropathy

Background: The effects of renal denervation on cardiovascular reflexes and markers of nephropathy in diabetic-hypertensive rats have not yet been explored. Methods: Aim: To evaluate the effects of renal denervation on nephropathy development mechanisms (blood pressure, cardiovascular autonomic changes, renal GLUT2) in diabetic-hypertensive rats. Forty-one male spontaneously hypertensive rats (SHR) similar to 250 g were injected with STZ or not; 30 days later, surgical renal denervation (RD) or sham procedure was performed; 15 days later, glycemia and albuminuria (ELISA) were evaluated. Catheters were implanted into the femoral artery to evaluate arterial pressure (AP) and heart rate variability (spectral analysis) one day later in conscious animals. Animals were killed, kidneys removed, and cortical renal GLUT2 quantified (Western blotting). Results: Higher glycemia (p < 0.05) and lower mean AP were observed in diabetics vs. nondiabetics (p < 0.05). Heart rate was higher in renal-denervated hypertensive and lower in diabetic-hypertensive rats (384.8 +/- 37, 431.3 +/- 36, 316.2 +/- 5, 363.8 +/- 12 bpm in SHR, RD-SHR, STZ-SHR and RD-STZ-SHR, respectively). Heart rate variability was higher in renal-denervated diabetic-hypertensive rats (55.75 +/- 25.21, 73.40 +/- 53.30, 148.4 +/- 93 in RD-SHR, STZ-SHR-and RD-STZ-SHR, respectively, p < 0.05), as well as the LF component of AP variability (1.62 +/- 0.9, 2.12 +/- 0.9, 7.38 +/- 6.5 in RD-SHR, STZ-SHR and RD-STZ-SHR, respectively, p < 0.05). GLUT2 renal content was higher in all groups vs. SHR. Conclusions: Renal denervation in diabetic-hypertensive rats improved previously reduced heart rate variability. The GLUT2 equally overexpressed by diabetes and renal denervation may represent a maximal derangement effect of each condition.

Melatonin treatment decreases c-fos expression in a headache model induced by capsaicin

The aim of the present work was to analyze c-fos response within the trigeminal nucleus caudalis (TNC) of pinealectomized rats and animals that received intraperitoneal melatonin, after intracisternal infusion of capsaicin, used to induce intracranial trigeminovascular stimulation. Experimental groups consisted of animals that received vehicle solution (saline-ethanol-Tween 80, 8:1:1, diluted 1:50) only (VEI, n = 5); animals that received capsaicin solution (200 nM) only (CAP, n = 6); animals submitted to pinealectomy (PX, n = 5); sham-operated animals (SH, n = 5); animals submitted to pinealectomy followed by capsaicin stimulation (200 nM) after 15 days (PX + CAP, n = 7); and animals that received capsaicin solution (200 nM) and intraperitoneal melatonin (10 mg/kg) (CAP + MEL, n = 5). Control rats, receiving vehicle in the cisterna magna, showed a small number of c-fos-positive cells in the TNC (layer I/II) as well as the sham-operated and pinealectomized rats, when compared to animals stimulated by capsaicin. On the other hand, pinealectomized rats, which received capsaicin, presented the highest number of c-fos-positive cells. Animals receiving capsaicin and melatonin treatment had similar expression of the vehicle group. Our data provide experimental evidence to support the role of melatonin and pineal gland in the pathophysiology of neurovascular headaches.