Plant species diversity in remnant grasslands on arable soils in the cropping Pampa

In order to provide a baseline for management decisions, we evaluated vegetation heterogeneity among 38 remnant grasslands located on high agronomic potential soils along a five degree-east–west transect in the cropping Pampa. This extensive region has been cultivated for well over a century, small-scale traditional forms of agriculture were progressively replaced with large-scale industrial forms and its cropped area has more than doubled in the last 25 years, threatening the maintenance of biodiversity associated to productive soils. We characterized the diversity of this grassland patches, identified its main environmental drivers and searched for factors that favor their invasion by exotic plant species. Vegetation samples were analysed through Canonical Correspondence Analysis, proportional distribution of functional groups was compared and regression models were estimated to identify variables associated with exotic species importance. Water availability, determined by climate and soil characteristics, was the main driving factor of floristic heterogeneity in these grasslands. Exotic richness, mostly of Eurasian species, increases toward the east, related to proximity to the entrance port. The functional group of C₄ native grasses seems to be responsible for biotic resistance to invasion. More than 80% of the species recorded during the first vegetation registers in 1930 are still present, while the number of exotics was seven-fold increased. Despite covering small areas as isolated stands, these are hot spots of native biodiversity and must be the object of increasing attention for conservation and restoration in order to protect the biodiversity of an extensive agricultural landscape.     

A GWAS assessment of the contribution of genomic imprinting to the variation of body mass index in mice

Background

Genomic imprinting is an epigenetic mechanism that can lead to differential gene expression depending on the parent-of-origin of a received allele. While most studies on imprinting address its underlying molecular mechanisms or attempt at discovering genomic regions that might be subject to imprinting, few have focused on the amount of phenotypic variation contributed by such epigenetic process. In this report, we give a brief review of a one-locus imprinting model in a quantitative genetics framework, and provide a decomposition of the genetic variance according to this model. Analytical deductions from the proposed imprinting model indicated a non-negligible contribution of imprinting to genetic variation of complex traits. Also, we performed a whole-genome scan analysis on mouse body mass index (BMI) aiming at revealing potential consequences when existing imprinting effects are ignored in genetic analysis.

Results

10,021 SNP markers were used to perform a whole-genome single marker regression on mouse BMI using an additive and an imprinting model. Markers significant for imprinting indicated that BMI is subject to imprinting. Marked variance changed from 1.218 ×10⁻⁴ to 1.842 ×10⁻⁴ when imprinting was considered in the analysis, implying that one third of marked variance would be lost if existing imprinting effects were not accounted for. When both marker and pedigree information were used, estimated heritability increased from 0.176 to 0.195 when imprinting was considered.

Conclusions

When a complex trait is subject to imprinting, using an additive model that ignores this phenomenon may result in an underestimate of additive variability, potentially leading to wrong inferences about the underlying genetic architecture of that trait. This could be a possible factor explaining part of the missing heritability commonly observed in genome-wide association studies (GWAS).     

Defining window-boundaries for genomic analyses using smoothing spline techniques

Background

High-density genomic data is often analyzed by combining information over windows of adjacent markers. Interpretation of data grouped in windows versus at individual locations may increase statistical power, simplify computation, reduce sampling noise, and reduce the total number of tests performed. However, use of adjacent marker information can result in over- or under-smoothing, undesirable window boundary specifications, or highly correlated test statistics. We introduce a method for defining windows based on statistically guided breakpoints in the data, as a foundation for the analysis of multiple adjacent data points. This method involves first fitting a cubic smoothing spline to the data and then identifying the inflection points of the fitted spline, which serve as the boundaries of adjacent windows. This technique does not require prior knowledge of linkage disequilibrium, and therefore can be applied to data collected from individual or pooled sequencing experiments. Moreover, in contrast to existing methods, an arbitrary choice of window size is not necessary, since these are determined empirically and allowed to vary along the genome.

Results

Simulations applying this method were performed to identify selection signatures from pooled sequencing F_ST data, for which allele frequencies were estimated from a pool of individuals. The relative ratio of true to false positives was twice that generated by existing techniques. A comparison of the approach to a previous study that involved pooled sequencing F_ST data from maize suggested that outlying windows were more clearly separated from their neighbors than when using a standard sliding window approach.

Conclusions

We have developed a novel technique to identify window boundaries for subsequent analysis protocols. When applied to selection studies based on F_ST data, this method provides a high discovery rate and minimizes false positives. The method is implemented in the R package GenWin, which is publicly available from CRAN.     

Deficiency in the mitochondrial apoptotic pathway reveals the toxic potential of autophagy under ER stress conditions

Endoplasmic reticulum (ER) stress-induced cell death is normally associated with activation of the mitochondrial apoptotic pathway, which is characterized by CYCS (cytochrome c, somatic) release, apoptosome formation, and caspase activation, resulting in cell death. In this study, we demonstrate that under conditions of ER stress cells devoid of CASP9/caspase-9 or BAX and BAK1, and therefore defective in the mitochondrial apoptotic pathway, still undergo a delayed form of cell death associated with the activation of caspases, therefore revealing the existence of an alternative stress-induced caspase activation pathway. We identified CASP8/caspase-8 as the apical protease in this caspase cascade, and found that knockdown of either of the key autophagic genes, ATG5 or ATG7, impacted on CASP8 activation and cell death induction, highlighting the crucial role of autophagy in the activation of this novel ER stress-induced death pathway. In line with this, we identified a protein complex composed of ATG5, FADD, and pro-CASP8 whose assembly coincides with caspase activation and cell death induction. Together, our results reveal the toxic potential of autophagy in cells undergoing ER stress that are defective in the mitochondrial apoptotic pathway, and suggest a model in which the autophagosome functions as a platform facilitating pro-CASP8 activation. Chemoresistance, a common problem in the treatment of cancer, is frequently caused by the downregulation of key mitochondrial death effector proteins. Alternate stress-induced apoptotic pathways, such as the one described here, may become of particular relevance for tackling the problem of chemoresistance in cancer cells.     

An improved formulation of the zirconyl hematoxylin stain for acidic mucins

Zirconyl hematoxylin stains acidic mucins darkly and specifically using a solution of 100 mg hematoxylin, 5 ml ethanol, 5 ml 0.5% sodium iodate, 400 mg zirconyl chloride octahydrate, and 30 ml 25% aqueous glycerol. The stain is especially advantageous for studying goblet cells and Paget cells.     

Rate of progression of CT-quantified emphysema in male current and ex-smokers: a follow-up study

Background

Little is known about the factors associated with CT-quantified emphysema progression in heavy smokers. The objective of this study was to investigate the effect of length of smoking cessation and clinical / demographical factors on the rate of emphysema progression and FEV₁-decline in male heavy smokers.

Methods

3,670 male smokers with mean (SD) 40.8 (17.9) packyears underwent chest CT scans and pulmonary function tests at baseline and after 1 and 3 years follow-up. Smoking status (quitted ≥5, ≥1-<5, <1 years or current smoker) was noted. Rate of progression of emphysema and FEV₁-decline after follow-up were assessed by analysis of variance adjusting for age, height, baseline pulmonary function and emphysema severity, packyears, years in study and respiratory symptoms. The quitted ≥5 group was used as reference.

Results

Median (Q1-Q3) emphysema severity,<-950 HU, was 8.8 (5.1 – 14.1) and mean (SD) FEV₁ was 3.4 (0.73) L or 98.5 (18.5) % of predicted. The group quitted ‘>5 years’ showed significantly lower rates of progression of emphysema compared to current smokers, 1.07% and 1.12% per year, respectively (p<0.001). Current smokers had a yearly FEV₁-decline of 69 ml, while subjects quit smoking >5 years had a yearly decline of 57.5 ml (p<0.001).

Conclusion

Quit smoking >5 years significantly slows the rate of emphysema progression and lung function decline.

Trial registration

Registered at http://www.trialregister.nl with trial number ISRCTN63545820.