Predicting survival from microarray data--a comparative study

MOTIVATION: Survival prediction from gene expression data and other high-dimensional genomic data has been subject to much research during the last years. These kinds of data are associated with the methodological problem of having many more gene expression values than individuals. In addition, the responses are censored survival times. Most of the proposed methods handle this by using Cox's proportional hazards model and obtain parameter estimates by some dimension reduction or parameter shrinkage estimation technique. Using three well-known microarray gene expression data sets, we compare the prediction performance of seven such methods: univariate selection, forward stepwise selection, principal components regression (PCR), supervised principal components regression, partial least squares regression (PLS), ridge regression and the lasso. RESULTS: Statistical learning from subsets should be repeated several times in order to get a fair comparison between methods. Methods using coefficient shrinkage or linear combinations of the gene expression values have much better performance than the simple variable selection methods. For our data sets, ridge regression has the overall best performance. AVAILABILITY: Matlab and R code for the prediction methods are available at http://www.med.uio.no/imb/stat/bmms/software/microsurv/.     

Variations of the candidate SEZ6L2 gene on Chromosome 16p11.2 in patients with autism spectrum disorders and in human populations

Background

Autism spectrum disorders (ASD) are a group of severe childhood neurodevelopmental disorders with still unknown etiology. One of the most frequently reported associations is the presence of recurrent de novo or inherited microdeletions and microduplications on chromosome 16p11.2. The analysis of rare variations of 8 candidate genes among the 27 genes located in this region suggested SEZ6L2 as a compelling candidate.

Methodology/Principal Findings

We further explored the role of SEZ6L2 variations by screening its coding part in a group of 452 individuals, including 170 patients with ASD and 282 individuals from different ethnic backgrounds of the Human Genome Diversity Panel (HGDP), complementing the previously reported screening. We detected 7 previously unidentified non-synonymous variations of SEZ6L2 in ASD patients. We also identified 6 non-synonymous variations present only in HGDP. When we merged our results with the previously published, no enrichment of non-synonymous variation in SEZ6L2 was observed in the ASD group compared with controls.

Conclusions/Significance

Our results provide an extensive ascertainment of the genetic variability of SEZ6L2 in human populations and do not support a major role for SEZ6L2 sequence variations in the susceptibility to ASD.     

N-CAM exhibits a regulatory function in pathological angiogenesis in oxygen induced retinopathy

Background

Diabetic retinopathy and retinopathy of prematurity are diseases caused by pathological angiogenesis in the retina as a consequence of local hypoxia. The underlying mechanism for epiretinal neovascularization (tuft formation), which contributes to blindness, has yet to be identified. Neural cell adhesion molecule (N-CAM) is expressed by Müller cells and astrocytes, which are in close contact with the retinal vasculature, during normal developmental angiogenesis.

Methodology/Principal Findings

Notably, during oxygen induced retinopathy (OIR) N-CAM accumulated on astrocytes surrounding the epiretinal tufts. Here, we show that N-CAM ablation results in reduced vascular tuft formation due to reduced endothelial cell proliferation despite an elevation in VEGFA mRNA expression, whereas retinal developmental angiogenesis was unaffected.

Conclusion/Significance

We conclude that N-CAM exhibits a regulatory function in pathological angiogenesis in OIR. This is a novel finding that can be of clinical relevance in diseases associated with proliferative vasculopathy.     

Greenhouse gas emissions from a constructed wetland in southern Sweden

This paper investigates the greenhouse gas emissions from a Swedish wetland, constructed to decrease nutrient content in sewage treatment water. To evaluate the effect of the construction in terms of greenhouse gas emissions we carried out ecosystem-atmosphere flux measurements of CO₂, CH₄ and N₂O using a closed chamber technique. To evaluate the importance of vascular plant species composition to gas emissions we distributed the measurement plots over the three dominating plant species at the field site, i.e., Typha latifolia, Phragmites australis and Juncus effusus. The fluxes of CO₂ (total respiration), CH₄ and N₂O from vegetated plots ranged from 1.39 to 77.5 (g m⁻² day⁻¹), −377 to 1387 and −13.9 to 31.5 (mg m⁻² day⁻¹) for CO₂, CH₄ and N₂O, respectively. Presence of vascular plants lead as expected to significantly higher total respiration rates compared with un-vegetated control plots. Furthermore, we found that the emission rates of N₂O and CH₄ was affected by presence of vascular plants and tended to be species-specific. We assessed the integrated greenhouse warming effect of the emissions using a Global Warming Potential over a 100-year horizon (GWP₁₀₀) and it corresponded to 431 kg CO₂ equivalents m⁻² day⁻¹. Assuming a 7-month season with conditions similar to the study period this is equal to 90 tonnes of CO₂ equivalents annually. N₂O emissions were responsible for one third of the estimated total greenhouse forcing. Furthermore, we estimated that the emission from the forested bog that was the precursor land to Magle constructed wetland amounted to 18.6 tonnes of CO₂ equivalents annually. Hence, the constructed wetland has increased annual greenhouse gas emissions by 71.4 tonnes of CO₂ equivalents for the whole area. Our findings indicate that management processes in relation to wetland construction projects must consider the primary function of the wetland in decreasing eutrophication, in relation to other positive aspects on for instance plant and animal life and recreation as well as possible negative climatic aspects of increased emissions of CH₄ and N₂O.     

Development of ten new EST-derived microsatellites in Atlantic cod (Gadus morhua L.)

Ten polymorphic microsatellite markers were developed from approximately 1,300 expressed sequence tags (ESTs) of Atlantic cod (Gadus morhua L.). Thirty two primer pairs were designed for EST sequences containing perfect di- tri- tetra- and pentanucleotide motifs and characterised in 96 unrelated fish. Ten markers were successfully amplified with number of alleles from 2 to 13 per locus and observed and expected heterozygosity ranging from 0.03 to 0.69 and 0.03 to 0.74, respectively. Loci Gmo-C131, C132 and C136 deviated from Hardy-Weinberg equilibrium. Genetic linkage disequilibrium analysis between all pairs of the loci showed significant departure from the null hypothesis between loci Gmo-C131 and Gmo-C132 and C128 and Gmo-C133. The gene identity was determined at five of the loci, confirming the associated microsatellites as Type I markers. The new microsatellites reported in this work can be used for conservation and enhancement of wild stocks for commercial harvesting. Jon-Ivar Westgaard and Tekle Tafese have contributed equally to the work.     

The major outer membrane protein of Fusobacterium nucleatum (FomA) folds and inserts into lipid bilayers via parallel folding pathways

Membrane protein insertion and folding was studied for the major outer membrane protein of Fusobacterium nucleatum (FomA), which is a voltage-dependent general diffusion porin. The transmembrane domain of FomA forms a beta-barrel that is predicted to consist of 14 beta-strands. Here, unfolded FomA is shown to insert and fold spontaneously and quantitatively into phospholipid bilayers upon dilution of the denaturant urea, which was shown previously only for outer membrane protein A (OmpA) of Escherichia coli. Folding of FomA is demonstrated by circular dichroism and fluorescence spectroscopy, by SDS-polyacrylamide gel electrophoresis, and by single-channel recordings. Refolded FomA had a single-channel conductance of 1.1 nS at 1 M KCl, in agreement with the conductance of FomA isolated from membranes in native form. In contrast to OmpA, which forms a smaller eight-stranded beta-barrel domain, folding kinetics of the larger FomA were slower and provided evidence for parallel folding pathways of FomA into lipid bilayers. Two pathways were observed independent of membrane thickness with two different lipid bilayers, which were either composed of dicapryl phosphatidylcholine or dioleoyl phosphatidylcholine. This is the first observation of parallel membrane insertion and folding pathways of a beta-barrel membrane protein from an unfolded state in urea into lipid bilayers. The kinetics of both folding pathways depended on the chain length of the lipid and on temperature with estimated activation energies of 19 kJ/mol (dicapryl phosphatidylcholine) and 70 kJ/mol (dioleoyl phosphatidylcholine) for the faster pathways.     

Genetics of self/nonself recognition in Serpula lacrymans

This study provides an analysis of the vegetative incompatibility system in Serpula lacrymans (Basidiomycota), a genetic system used to recognize nonself in fungi. Seventy-five worldwide isolates could be grouped into eight vegetative compatibility (VC) types, some of them distributed on different continents. Mating studies combined with vegetative incompatibility analyses revealed that the vegetative incompatibility response between isolates mainly could be explained by two biallelic vegetative incompatibility (vic) loci. The frequency distributions of the interpreted vic alleles do not seem to support the idea of frequency-dependent or balancing selection acting on the vic loci. We find little genetic variation at the vic loci and in one of the loci there was a significant heterozyote deficiency among strains in the overall material. The results may be explained by a recent worldwide dispersal of a few S. lacrymans isolates and, correspondingly, only a few vic alleles are being maintained in these populations.     

Phosphorus and carbon availability regulate structural composition and complexity of AM fungal mycelium

The regulation of the structural composition and complexity of the mycelium of arbuscular mycorrhizal (AM) fungi is not well understood due to their obligate biotrophic nature. The aim of this study was to investigate the structure of extraradical mycelium at high and low availability of carbon (C) to the roots and phosphorus (P) to the fungus. We used monoxenic cultures of the AM fungus Rhizophagus irregularis (formerly Glomus intraradices) with transformed carrot roots as the host in a cultivation system including a root-free compartment into which the extraradical mycelium could grow. We found that high C availability increased hyphal length and spore production and anastomosis formation within individual mycelia. High P availability increased the formation of branched absorbing structures and reduced spore production and the overall length of runner hyphae. The complexity of the mycelium, as indicated by its fractal dimensions, increased with both high C and P availability. The results indicate that low P availability induces a growth pattern that reflects foraging for both P and C. Low C availability to AM roots could still support the explorative development of the mycelium when P availability was low. These findings help us to better understand the development of AM fungi in ecosystems with high P input and/or when plants are subjected to shading, grazing or any management practice that reduces the photosynthetic ability of the plant.