Identification of a common risk haplotype for canine idiopathic epilepsy in the ADAM23 gene

Background

Idiopathic epilepsy is a common neurological disease in human and domestic dogs but relatively few risk genes have been identified to date. The seizure characteristics, including focal and generalised seizures, are similar between the two species, with gene discovery facilitated by the reduced genetic heterogeneity of purebred dogs. We have recently identified a risk locus for idiopathic epilepsy in the Belgian Shepherd breed on a 4.4 megabase region on CFA37.

Results

We have expanded a previous study replicating the association with a combined analysis of 157 cases and 179 controls in three additional breeds: Schipperke, Finnish Spitz and Beagle (p_c = 2.9e–07, p_GWAS = 1.74E-02). A targeted resequencing of the 4.4 megabase region in twelve Belgian Shepherd cases and twelve controls with opposite haplotypes identified 37 case-specific variants within the ADAM23 gene. Twenty-seven variants were validated in 285 cases and 355 controls from four breeds, resulting in a strong replication of the ADAM23 locus (p_raw = 2.76e–15) and the identification of a common 28 kb-risk haplotype in all four breeds. Risk haplotype was present in frequencies of 0.49–0.7 in the breeds, suggesting that ADAM23 is a low penetrance risk gene for canine epilepsy.

Conclusions

These results implicate ADAM23 in common canine idiopathic epilepsy, although the causative variant remains yet to be identified. ADAM23 plays a role in synaptic transmission and interacts with known epilepsy genes, LGI1 and LGI2, and should be considered as a candidate gene for human epilepsies.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1651-9) contains supplementary material, which is available to authorized users.     

A bi-filtering method for processing single nucleotide polymorphism array data improves the quality of genetic map and accuracy of quantitative trait locus mapping in doubled haploid populations of polyploid Brassica napus

Background

Single nucleotide polymorphism (SNP) markers have a wide range of applications in crop genetics and genomics. Due to their polyploidy nature, many important crops, such as wheat, cotton and rapeseed contain a large amount of repeat and homoeologous sequences in their genomes, which imposes a huge challenge in high-throughput genotyping with sequencing and/or array technologies. Allotetraploid Brassica napus (AACC, 2n = 4x = 38) comprises of two highly homoeologous sub-genomes derived from its progenitor species B. rapa (AA, 2n = 2x = 20) and B. oleracea (CC, 2n = 2x = 18), and is an ideal species to exploit methods for reducing the interference of extensive inter-homoeologue polymorphisms (mHemi-SNPs and Pseudo-simple SNPs) between closely related sub-genomes.

Results

Based on a recent B. napus 6K SNP array, we developed a bi-filtering procedure to identify unauthentic lines in a DH population, and mHemi-SNPs and Pseudo-simple SNPs in an array data matrix. The procedure utilized both monomorphic and polymorphic SNPs in the DH population and could effectively distinguish the mHemi-SNPs and Pseudo-simple SNPs that resulted from superposition of the signals from multiple SNPs. Compared with conventional procedure for array data processing, the bi-filtering method could minimize the pseudo linkage relationship caused by the mHemi-SNPs and Pseudo-simple SNPs, thus improving the quality of SNP genetic map. Furthermore, the improved genetic map could increase the accuracies of mapping of QTLs as demonstrated by the ability to eliminate non-real QTLs in the mapping population.

Conclusions

The bi-filtering analysis of the SNP array data represents a novel approach to effectively assigning the multi-loci SNP genotypes in polyploid B. napus and may find wide applications to SNP analyses in polyploid crops.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1559-4) contains supplementary material, which is available to authorized users.     

Impaired macrophage phagocytosis of bacteria in severe asthma

Background

Bacteria are frequently cultured from sputum samples of severe asthma patients suggesting a defect in bacterial clearance from the airway. We measured the capacity of macrophages from patients with asthma to phagocytose bacteria.

Methods

Phagocytosis of fluorescently-labelled polystyrene beads, Haemophilus influenzae or Staphylococcus aureus by broncholaveolar lavage alveolar macrophages (AM) and by monocyte-derived macrophages (MDM) from non-asthmatics, mild-moderate and severe asthmatic patients was assessed using fluorimetry.

Results

There were no differences in phagocytosis of polystyrene beads by AMs or MDMs from any of the subject groups. There was reduced phagocytosis of Haemophilus influenzae and Staphylococcus aureus in MDMs from patients with severe asthma compared to non-severe asthma (p < 0.05 and p < 0.01, respectively) and healthy subjects (p < 0.01and p < 0.001, respectively). Phagocytosis of Haemophilus influenzae and Staphylococcus aureus by AM was also reduced in severe asthma compared to normal subjects (p < 0.05). Dexamethasone and formoterol did not suppress phagocytosis of bacteria by MDMs from any of the groups.

Conclusions

Persistence of bacteria in the lower airways may result partly from a reduced phagocytic capacity of macrophages for bacteria. This may contribute to increased exacerbations, airway colonization and persistence of inflammation.     

The prevalence and identity of Chlamydia-specific IgE in children with asthma and other chronic respiratory symptoms

Background

Recent studies have confirmed the presence of viable Chlamydia in the bronchoalveolar lavage (BAL) fluid of pediatric patients with airway hyperresponsiveness. While specific IgG and IgM responses to C. pneumoniae are well described, the response and potential contribution of Ag-specific IgE are not known. The current study sought to determine if infection with Chlamydia triggers the production of pathogen-specific IgE in children with chronic respiratory diseases which might contribute to inflammation and pathology.

Methods

We obtained BAL fluid and serum from pediatric respiratory disease patients who were generally unresponsive to corticosteroid treatment as well as sera from age-matched control patients who saw their doctor for wellness checkups. Chlamydia-specific IgE was isolated from BAL and serum samples and their specificity determined by Western blot techniques. The presence of Chlamydia was confirmed by species-specific PCR and BAL culture assays.

Results

Chlamydial DNA was detected in the BAL fluid of 134/197 (68%) patients. Total IgE increased with age until 15 years old and then decreased. Chlamydia-specific IgE was detected in the serum and/or BAL of 107/197 (54%) patients suffering from chronic respiratory disease, but in none of the 35 healthy control sera (p < 0.0001). Of the 74 BAL culture-positive patients, 68 (91.9%, p = 0.0001) tested positive for Chlamydia-specific IgE. Asthmatic patients had significantly higher IgE levels compared to non-asthmatics (p = 0.0001). Patients who were positive for Chlamydia DNA or culture had significantly higher levels of serum IgE compared to negative patients (p = 0.0071 and p = 0.0001 respectively). Only 6 chlamydial antigens induced Chlamydia-specific IgE and patients with C. pneumoniae-specific IgE had significantly greater levels of total IgE compared to C. pneumoniae-specific IgE negative ones (p = 0.0001).

Conclusions

IgE antibodies play a central role in allergic inflammation; therefore production of Chlamydia-specific IgE may prove significant in the exacerbation of chronic, allergic airway diseases, thus highlighting a direct role for Chlamydia in asthma pathogenesis.     

De novo mutations in ARID1B associated with both syndromic and non-syndromic short stature

Background

Human height is a complex trait with a strong genetic basis. Recently, a significant association between rare copy number variations (CNVs) and short stature has been identified, and candidate genes in these rare CNVs are being explored. This study aims to evaluate the association between mutations in ARID1B gene and short stature, both the syndromic and non-syndromic form.

Results

Based on a case-control study of whole genome chromosome microarray analysis (CMA), three overlapping CNVs were identified in patients with developmental disorders who exhibited short stature. ARID1B, a causal gene for Coffin Siris syndrome, is the only gene encompassed by all three CNVs. A following retrospective genotype-phenotype analysis based on a literature review confirmed that short stature is a frequent feature in those Coffin-Siris syndrome patients with ARID1B mutations. Mutation screening of ARID1B coding regions was further conducted in a cohort of 48 non-syndromic short stature patients,andfour novel missense variants including two de novo mutations were found.

Conclusion

These results suggest that haploinsufficient mutations of ARID1B are associated with syndromic short stature including Coffin-Siris syndrome and intellectual disability, while rare missense variants in ARID1B are associated with non-syndromic short stature. This study supports the notion that mutations in genes related to syndromic short stature may exert milder effect and contribute to short stature in the general population.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1898-1) contains supplementary material, which is available to authorized users.     

Association of 25-Hydroxyvitamin D status and genetic variation in the vitamin D metabolic pathway with FEV1 in the Framingham Heart Study

Background

Vitamin D is associated with lung function in cross-sectional studies, and vitamin D inadequacy is hypothesized to play a role in the pathogenesis of chronic obstructive pulmonary disease. Further data are needed to clarify the relation between vitamin D status, genetic variation in vitamin D metabolic genes, and cross-sectional and longitudinal changes in lung function in healthy adults.

Methods

We estimated the association between serum 25-hydroxyvitamin D [25(OH)D] and cross-sectional forced expiratory volume in the first second (FEV₁) in Framingham Heart Study (FHS) Offspring and Third Generation participants and the association between serum 25(OH)D and longitudinal change in FEV₁ in Third Generation participants using linear mixed-effects models. Using a gene-based approach, we investigated the association between 241 SNPs in 6 select vitamin D metabolic genes in relation to longitudinal change in FEV₁ in Offspring participants and pursued replication of these findings in a meta-analyzed set of 4 independent cohorts.

Results

We found a positive cross-sectional association between 25(OH)D and FEV₁ in FHS Offspring and Third Generation participants (P = 0.004). There was little or no association between 25(OH)D and longitudinal change in FEV₁ in Third Generation participants (P = 0.97). In Offspring participants, the CYP2R1 gene, hypothesized to influence usual serum 25(OH)D status, was associated with longitudinal change in FEV₁ (gene-based P < 0.05). The most significantly associated SNP from CYP2R1 had a consistent direction of association with FEV₁ in the meta-analyzed set of replication cohorts, but the association did not reach statistical significance thresholds (P = 0.09).

Conclusions

Serum 25(OH)D status was associated with cross-sectional FEV₁, but not longitudinal change in FEV₁. The inconsistent associations may be driven by differences in the groups studied. CYP2R1 demonstrated a gene-based association with longitudinal change in FEV₁ and is a promising candidate gene for further studies.

Electronic supplementary material

The online version of this article (doi:10.1186/s12931-015-0238-y) contains supplementary material, which is available to authorized users.     

Mapping in the era of sequencing: high density genotyping and its application for mapping TYLCV resistance in Solanum pimpinellifolium

Background

A RIL population between Solanum lycopersicum cv. Moneymaker and S. pimpinellifolium G1.1554 was genotyped with a custom made SNP array. Additionally, a subset of the lines was genotyped by sequencing (GBS).

Results

A total of 1974 polymorphic SNPs were selected to develop a linkage map of 715 unique genetic loci. We generated plots for visualizing the recombination patterns of the population relating physical and genetic positions along the genome.This linkage map was used to identify two QTLs for TYLCV resistance which contained favourable alleles derived from S. pimpinellifolium. Further GBS was used to saturate regions of interest, and the mapping resolution of the two QTLs was improved. The analysis showed highest significance on Chromosome 11 close to the region of 51.3 Mb (qTy-p11) and another on Chromosome 3 near 46.5 Mb (qTy-p3). Furthermore, we explored the population using untargeted metabolic profiling, and the most significant differences between susceptible and resistant plants were mainly associated with sucrose and flavonoid glycosides.

Conclusions

The SNP information obtained from an array allowed a first QTL screening of our RIL population. With additional SNP data of a RILs subset, obtained through GBS, we were able to perform an in silico mapping improvement to further confirm regions associated with our trait of interest. With the combination of different ~ omics platforms we provide valuable insight into the genetics of S. pimpinellifolium-derived TYLCV resistance.

Electronic supplementary material

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

DNAH5 is associated with total lung capacity in chronic obstructive pulmonary disease

Background

Chronic obstructive pulmonary disease (COPD) is characterized by expiratory flow limitation, causing air trapping and lung hyperinflation. Hyperinflation leads to reduced exercise tolerance and poor quality of life in COPD patients. Total lung capacity (TLC) is an indicator of hyperinflation particularly in subjects with moderate-to-severe airflow obstruction. The aim of our study was to identify genetic variants associated with TLC in COPD.

Methods

We performed genome-wide association studies (GWASs) in white subjects from three cohorts: the COPDGene Study; the Evaluation of COPD Longitudinally to Identify Predictive Surrogate Endpoints (ECLIPSE); and GenKOLS (Bergen, Norway). All subjects were current or ex-smokers with at least moderate airflow obstruction, defined by a ratio of forced expiratory volume in 1 second to forced vital capacity (FEV₁/FVC) <0.7 and FEV₁ < 80% predicted on post-bronchodilator spirometry. TLC was calculated by using volumetric computed tomography scans at full inspiration (TLC_CT). Genotyping in each cohort was completed, with statistical imputation of additional markers. To find genetic variants associated with TLC_CT, linear regression models were used, with adjustment for age, sex, pack-years of smoking, height, and principal components for genetic ancestry. Results were summarized using fixed-effect meta-analysis.

Results

Analysis of a total of 4,543 COPD subjects identified one genome-wide significant locus on chromosome 5p15.2 (rs114929486, β = 0.42L, P = 4.66 × 10⁻⁸).

Conclusions

In COPD, TLC_CT was associated with a SNP in dynein, axonemal, heavy chain 5 (DNAH5), a gene in which genetic variants can cause primary ciliary dyskinesia. DNAH5 could have an effect on hyperinflation in COPD.

Electronic supplementary material

The online version of this article (doi:10.1186/s12931-014-0097-y) contains supplementary material, which is available to authorized users.     

The influences of PRG-1 on the expression of small RNAs and mRNAs

Background

In metazoans, Piwi-related Argonaute proteins play important roles in maintaining germline integrity and fertility and have been linked to a class of germline-enriched small RNAs termed piRNAs. Caenorhabditis elegans encodes two Piwi family proteins called PRG-1 and PRG-2, and PRG-1 interacts with the C. elegans piRNAs (21U-RNAs). Previous studies found that mutation of prg-1 causes a marked reduction in the expression of 21U-RNAs, temperature-sensitive defects in fertility and other phenotypic defects.

Results

In this study, we wanted to systematically demonstrate the function of PRG-1 in the regulation of small RNAs and their targets. By analyzing small RNAs and mRNAs with and without a mutation in prg-1 during C. elegans development, we demonstrated that (1) mutation of prg-1 leads to a decrease in the expression of 21U-RNAs, and causes 35 ~ 40% of miRNAs to be down-regulated; (2) in C. elegans, approximately 3% (6% in L4) of protein-coding genes are differentially expressed after mutating prg-1, and 60 ~ 70% of these substantially altered protein-coding genes are up-regulated; (3) the target genes of the down-regulated miRNAs and the candidate target genes of the down-regulated 21U-RNAs are enriched in the up-regulated protein-coding genes; and (4) PRG-1 regulates protein-coding genes by down-regulating small RNAs (miRNAs and 21U-RNAs) that target genes that participate in the development of C. elegans.

Conclusions

In prg-1-mutated C. elegans, the expression of miRNAs and 21U-RNAs was reduced, and the protein-coding targets, which were associated with the development of C. elegans, were up-regulated. This may be the mechanism underlying PRG-1 function.

Electronic supplementary material

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

Assessment of pharmacological activities of two medicinal plant of Bangladesh: Launaea sarmentosa and Aegialitis rotundifolia roxb in the management of pain,pyrexia and inflammation

Background

The current study aims at evaluating the analgesic, anti-pyretic and anti-inflammatory properties of methanolic extract of the stem, bark and leaves of Launaea sarmentosa and Aegialitis rotundifolia roxb.

Results

The AELS and AEAR extract presented a significant (***p < 0.001) dose dependent increase in reaction time in writhing method and showed inhibition of 63.1% and 57.1% respectively at the doses of 400 mg/kg body weight while standard drug showed (P < 0.001) inhibition of 69.23%. In tail immersion method, AELS and AEAR showed maximum time of tail retention at 30 min in hot water i.e. 6.93 sec and 6.54 sec respectively at highest doses of 400 mg/kg body weight than lower dose while standard pentazocine showed reaction time of 7.62 sec. The AELS and AEAR extract also exhibited promising anti-inflammatory effect as demonstrated by statistically significant inhibition of paw volume by 32.48% and 26.75% respectively at the dose of 400 mg/kg body weight while the value at the dose of 200 mg/kg body weight were linear to higher dose at the 3^rd hour of study. On the other hand, Standard indomethacin inhibited 40.13% of inflammation (***P < 0.001). In Cotton-pellet granuloma method, AELS and AEAR extract at the dose of 400 mg/kg body weight exhibited inhibition of inflammation of 34.7% and 29.1% respectively while standard drug showed (P < 0.001) inhibition of 63.22%. Intraperitoneal administration of AELS and AEAR showed dose dependent decrease in body temperature in brewer’s yeast induced hyperthermia in rats at both doses. However, AELS significantly decreased body temperature (***p < 0.001) at 400 mg/kg compared to control.

Conclusions

Present work propose that the methanolic extract of Launaea sarmentosa and Aegialitis rotundifolia roxb possesses dose dependent pharmacological action which supports its therapeutic use in folk medicine possibly mediated through the inhibition or blocking of release of prostaglandin and/or actions of vasoactive substances such as histamine, serotonin and kinins.     

Identification of genome-wide single nucleotide polymorphisms in allopolyploid crop Brassica napus

Background

Single nucleotide polymorphisms (SNPs) are the most common type of genetic variation. Identification of large numbers of SNPs is helpful for genetic diversity analysis, map-based cloning, genome-wide association analyses and marker-assisted breeding. Recently, identifying genome-wide SNPs in allopolyploid Brassica napus (rapeseed, canola) by resequencing many accessions has become feasible, due to the availability of reference genomes of Brassica rapa (2n = AA) and Brassica oleracea (2n = CC), which are the progenitor species of B. napus (2n = AACC). Although many SNPs in B. napus have been released, the objective in the present study was to produce a larger, more informative set of SNPs for large-scale and efficient genotypic screening. Hence, short-read genome sequencing was conducted on ten elite B. napus accessions for SNP discovery. A subset of these SNPs was randomly selected for sequence validation and for genotyping efficiency testing using the Illumina GoldenGate assay.

Results

A total of 892,536 bi-allelic SNPs were discovered throughout the B. napus genome. A total of 36,458 putative amino acid variants were located in 13,552 protein-coding genes, which were predicted to have enriched binding and catalytic activity as a result. Using the GoldenGate genotyping platform, 94 of 96 SNPs sampled could effectively distinguish genotypes of 130 lines from two mapping populations, with an average call rate of 92%.

Conclusions

Despite the polyploid nature of B. napus, nearly 900,000 simple SNPs were identified by whole genome resequencing. These SNPs were predicted to be effective in high-throughput genotyping assays (51% polymorphic SNPs, 92% average call rate using the GoldenGate assay, leading to an estimated >450 000 useful SNPs). Hence, the development of a much larger genotyping array of informative SNPs is feasible. SNPs identified in this study to cause non-synonymous amino acid substitutions can also be utilized to directly identify causal genes in association studies.     

Effects of crowding and sex on fecal cortisol levels of captive forest musk deer

Background

Restricted space and close contact with conspecifics in captivity may be stressful for musk deer, as they are highly territorial and solitary in the wild. So we tested the effects of crowding on stress of forest musk deer (Moschus berezovskii) in heterosexual groups, using fecal cortisol analysis as a non-invasive method. 32 healthy adults during non-breeding seasons were chose as our experimental objects. Group 1 was defined as higher crowding condition, with 10-15 m²/deer (6 enclosures, 10♀ and 6♂); group 2 was defined as lower crowding condition, with 23-33 m²/deer (6 enclosures, 10♀ and 6♂). Every enclosure contained 1 male and 3 female. These patterns had been existed for years.

Results

The results showed that females in lower crowding condition (217.1 ± 9.5 ug/g) had significantly higher fecal cortisol levels than those in higher crowding condition (177.2 ± 12.1 ug/g). Interestingly, crowding seemed have no effect on male fecal cortisol levels (148.1 ± 9.1 ug/g and 140.5 ± 13.3 ug/g, respectively). At both groups, cortisol was significantly lower in males than in females.

Conclusions

These results showed that chronic crowding may affect stress status of captive forest musk deer. The captive environment should consider the space need for musk deer.     

Decrease in quality of life predicts mortality in adult patients with pulmonary arterial hypertension due to congenital heart disease

Background

Decrease in quality of life (QoL) in left-sided heart failure precedes poor survival, which can be reversed with exercise training. We investigated whether QoL is associated with mortality in pulmonary arterial hypertension due to congenital heart disease (PAH-CHD) patients.

Methods

In this observational study, PAH-CHD adults referred for PAH-specific therapy were included. QoL surveys (SF36) were recorded during 2 years of therapy. Based on shift in SF36 scores during this period, patients had either decreased or non-decreased QoL. Subsequently, the patients were followed for mortality.

Results

Thirty-nine PAH-CHD patients (mean age 42, 44 % male, 49 % Down’s syndrome) were analysed. Following PAH-specific therapy, SF36 physical component summary (PCS) decreased in 13 (35–31 points, p = 0.001) and showed no decrease in 26 patients (34–43 points, mean values, p < 0.001). Post-initiation phase, median follow-up was 4.5 years, during which 12 deaths occurred (31 %), 10 (56 %) in the decreased and 2 (10 %) in the non-decreased group (p = 0.002). Cox regression showed a decrease in SF36 PCS predicted mortality (HR 3.4, 95 % CI 1.03–11, p = 0.045).

Conclusions

In PAH-CHD patients, decrease in SF36 PCS following initiation of PAH-specific therapy is a determinant of mortality.

Electronic supplementary material

The online version of this article (doi:10.1007/s12471-015-0666-9) contains supplementary material, which is available to authorized users.     

Inferring mechanisms of copy number change from haplotype structures at the human DEFA1A3 locus

Background

The determination of structural haplotypes at copy number variable regions can indicate the mechanisms responsible for changes in copy number, as well as explain the relationship between gene copy number and expression. However, obtaining spatial information at regions displaying extensive copy number variation, such as the DEFA1A3 locus, is complex, because of the difficulty in the phasing and assembly of these regions. The DEFA1A3 locus is intriguing in that it falls within a region of high linkage disequilibrium, despite its high variability in copy number (n = 3–16); hence, the mechanisms responsible for changes in copy number at this locus are unclear.

Results

In this study, a region flanking the DEFA1A3 locus was sequenced across 120 independent haplotypes with European ancestry, identifying five common classes of DEFA1A3 haplotype. Assigning DEFA1A3 class to haplotypes within the 1000 Genomes project highlights a significant difference in DEFA1A3 class frequencies between populations with different ancestry. The features of each DEFA1A3 class, for example, the associated DEFA1A3 copy numbers, were initially assessed in a European cohort (n = 599) and replicated in the 1000 Genomes samples, showing within-class similarity, but between-class and between-population differences in the features of the DEFA1A3 locus. Emulsion haplotype fusion-PCR was used to generate 61 structural haplotypes at the DEFA1A3 locus, showing a high within-class similarity in structure.

Conclusions

Structural haplotypes across the DEFA1A3 locus indicate that intra-allelic rearrangement is the predominant mechanism responsible for changes in DEFA1A3 copy number, explaining the conservation of linkage disequilibrium across the locus. The identification of common structural haplotypes at the DEFA1A3 locus could aid studies into how DEFA1A3 copy number influences expression, which is currently unclear.

Electronic supplementary material

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

FEV1 inversely correlates with metalloproteinases 1, 7, 9 and CRP in COPD by biomass smoke exposure

Background

Matrix metalloproteinases (MMPs) and C-reactive protein (CRP) are involved in chronic obstructive pulmonary disease (COPD) pathogenesis. The aim of the present work was to determine plasma concentrations of MMPs and CRP in COPD associated to biomass combustion exposure (BE) and tobacco smoking (TS).

Methods

Pulmonary function tests, plasma levels of MMP-1, MMP-7, MMP-9, MMP-9/TIMP-1 and CRP were measured in COPD associated to BE (n = 40) and TS (n =40) patients, and healthy non-smoking (NS) healthy women (controls, n = 40).

Results

Plasma levels of MMP-1, MMP-7, MMP-9, and MMP-9/TIMP-1 and CRP were higher in BE and TS than in the NS healthy women (p <0.01). An inverse correlation between MMP-1, MMP-7, MMP-9, MMP-9/TIMP-1 and CRP plasma concentrations and FEV₁ was observed.

Conclusions

Increase of MMPs and CRP plasma concentrations in BE suggests a systemic inflammatory phenomenon similar to that observed in COPD associated to tobacco smoking, which may also play a role in COPD pathogenesis.