GSVA: gene set variation analysis for microarray and RNA-seq data

Gene set enrichment (GSE) analysis is a popular framework for condensing information from gene expression profiles into a pathway or signature summary. The strengths of this approach over single gene analysis include noise and dimension reduction, as well as greater biological interpretability. As molecular profiling experiments move beyond simple case-control studies, robust and flexible GSE methodologies are needed that can model pathway activity within highly heterogeneous data sets. To address this challenge, we introduce Gene Set Variation Analysis (GSVA), a GSE method that estimates variation of pathway activity over a sample population in an unsupervised manner. We demonstrate the robustness of GSVA in a comparison with current state of the art sample-wise enrichment methods. Further, we provide examples of its utility in differential pathway activity and survival analysis. Lastly, we show how GSVA works analogously with data from both microarray and RNA-seq experiments. GSVA provides increased power to detect subtle pathway activity changes over a sample population in comparison to corresponding methods. While GSE methods are generally regarded as end points of a bioinformatic analysis, GSVA constitutes a starting point to build pathway-centric models of biology. Moreover, GSVA contributes to the current need of GSE methods for RNA-seq data. GSVA is an open source software package for R which forms part of the Bioconductor project and can be downloaded at http://www.bioconductor.org.

Normal variation in fronto-occipital circuitry and cerebellar structure with an autism-associated polymorphism of CNTNAP2.

Recent genetic studies have implicated a number of candidate genes in the pathogenesis of Autism Spectrum Disorder (ASD). Polymorphisms of CNTNAP2 (contactin-associated like protein-2), a member of the neurexin family, have already been implicated as a susceptibility gene for autism by at least 3 separate studies. We investigated variation in white and grey matter morphology using structural MRI and diffusion tensor imaging. We compared volumetric differences in white and grey matter and fractional anisotropy values in control subjects characterised by genotype at rs7794745, a single nucleotide polymorphism in CNTNAP2. Homozygotes for the risk allele showed significant reductions in grey and white matter volume and fractional anisotropy in several regions that have already been implicated in ASD, including the cerebellum, fusiform gyrus, occipital and frontal cortices. Male homozygotes for the risk alleles showed greater reductions in grey matter in the right frontal pole and in FA in the right rostral fronto-occipital fasciculus compared to their female counterparts who showed greater reductions in FA of the anterior thalamic radiation. Thus a risk allele for autism results in significant cerebral morphological variation, despite the absence of overt symptoms or behavioural abnormalities. The results are consistent with accumulating evidence of CNTNAP2's function in neuronal development. The finding suggests the possibility that the heterogeneous manifestations of ASD can be aetiologically characterised into distinct subtypes through genetic-morphological analysis.

Ip-10 Is Associated With Il28b Variation And Predicts The First Phase Decline Of Hcv Rna And Outcome Of Therapy In Chronic Hepatitis C

High systemic levels of IP-10 at onset of combination therapy for chronic hepatitis C mirror intrahepatic mRNA levels and predict a slower first phase decline in HCV RNA as well as poor outcome. Recently several genome wide association studies have revealed that single nucleotide polymorphisms (SNPs) on chromosome19 within proximity of IL28B predict spontaneous clearance of HCV infection and as therapeutic outcome among patients infected with HCV genotype 1, with three such SNPs being highly predictive: rs12979860, rs12980275, and rs8099917. In the present study, we correlated genetic variations in these SNPs from 253 Caucasian patients with pretreatment plasma levels of IP-10 and HCV RNA throughout therapy within a phase III treatment trial (HCV-DITTO). The favorable genetic variations in all three SNPs (CC, AA, and TT respectively) was significantly associated with lower baseline IP-10 (CC vs. CT/TT at rs12979860: median 189 vs. 258 pg/mL, P=0.02, AA vs. AG/GG at rs12980275: median 189 vs. 258 pg/mL, P=0.01, TT vs. TG/GG at rs8099917: median 224 vs. 288 pg/mL, P=0.04), were significantly less common among HCV genotype 1 infected patients than genotype 2/3 (P<0.0001, P<0.0001, and P=0.01 respectively) and had significantly higher baseline viral load than carriers of the SNP genotypes (6.3 vs. 5.9 log 10 IU/mL, P=0.0012, 6.3 vs. 6.0 log 10 IU/mL, P=0.026, and 6.3 vs. 5.8 log 10 IU/mL, P=0.0003 respectively). Among HCV genotype 1 infected homozygous or heterogeneous carriers of the favorable C, A, and T genotypes, lower baseline IP-10 was significantly associated with greater decline in HCV-RNA day 0-4, which translated into increased rates of achieving SVR among homozygous patients with baseline IP-10 below 150 pg/mL (85%, 75%, and 75% respectively). In a multivariate analysis among genotype 1 infected patients, both baseline IP-10 and the SNPs were significant independent predictors of SVR. Conclusion: Baseline plasma IP-10 is significantly associated with IL28B variations, and augments the predictiveness of the first phase decline in HCV RNA and final treatment outcome.

Structural variation and its effect on expression.

Structural variation, whether it is caused by copy number variants or present in a balanced form, such as reciprocal translocations and inversions, can have a profound and dramatic effect on the expression of genes mapping within and close to the rearrangement, as well as affecting others genome wide. These effects can be caused by altering the copy number of one or more genes or regulatory elements (dosage effect) or from physical disruption of links between regulatory elements and their associated gene or genes, resulting in perturbation of expression. Similarly, large-scale structural variants can result in genome-wide expression changes by altering the positions that chromosomes occupy within the nucleus, potentially disrupting not only local cis interactions, but also trans interactions that occur throughout the genome. Structural variation is, therefore, a significant factor in the study of gene expression and is discussed here in more detail.

A guide to understanding the variation in premiums in rural health insurance marketplaces (RUPRI Center for Rural Health Policy Analysis rural policy brief)

This brief discusses several important factors that should be considered when comparing health insurance plans in the health insurance marketplaces across geographic areas.

A hyperactive quantitative trait locus allele of Arabidopsis BRX contributes to natural variation in root growth vigor.

Quantitative trait loci analysis of natural Arabidopsis thaliana accessions is increasingly exploited for gene isolation. However, to date this has mostly revealed deleterious mutations. Among them, a loss-of-function allele identified the root growth regulator BREVIS RADIX (BRX). Here we present evidence that BRX and the paralogous BRX-LIKE (BRXL) genes are under selective constraint in monocotyledons as well as dicotyledons. Unexpectedly, however, whereas none of the Arabidopsis orthologs except AtBRXL1 could complement brx null mutants when expressed constitutively, nearly all monocotyledon BRXLs tested could. Thus, BRXL proteins seem to be more diversified in dicotyledons than in monocotyledons. This functional diversification was correlated with accelerated rates of sequence divergence in the N-terminal regions. Population genetic analyses of 30 haplotypes are suggestive of an adaptive role of AtBRX and AtBRXL1. In two accessions, Lc-0 and Lov-5, seven amino acids are deleted in the variable region between the highly conserved C-terminal, so-called BRX domains. Genotyping of 42 additional accessions also found this deletion in Kz-1, Pu2-7, and Ws-0. In segregating recombinant inbred lines, the Lc-0 allele (AtBRX(Lc-0)) conferred significantly enhanced root growth. Moreover, when constitutively expressed in the same regulatory context, AtBRX(Lc-0) complemented brx mutants more efficiently than an allele without deletion. The same was observed for AtBRXL1, which compared with AtBRX carries a 13 amino acid deletion that encompasses the deletion found in AtBRX(Lc-0). Thus, the AtBRX(Lc-0) allele seems to contribute to natural variation in root growth vigor and provides a rare example of an experimentally confirmed, hyperactive allelic variant.

Case Study of Seasonal Variation in the Subgrade and Subbase Layers of Highway US 20, TR-516, 2008

Seasonal variations in ground temperature and moisture content influence the load carrying capacity of pavement subgrade layers. To improve pavement performance, pavement design guidelines require knowledge of environmental factors and subgrade stiffness relationships. As part of this study, in-ground instrumentation was installed in the pavement foundation layers of a newly constructed section along US Highway 20 near Fort Dodge, Iowa, to monitor the seasonal variations in temperature, frost depth, groundwater levels, and moisture regime. Dynamic cone penetrometer (DCP), nuclear gauge, and Clegg hammer tests were performed at 64 test points in a 6-ft x 6-ft grid pattern to characterize the subgrade stiffness properties (i.e., resilient modulus) prior to paving. The purpose of this paper is to present the field instrumentation results and the observed changes in soil properties due to seasonal environmental effects.

Variation in biomass related variables of reed canary grass

Selostus: Ruokohelven biomassan tuotantoon vaikuttavien ominaisuuksien vaihtelu

Natural genetic variation in Arabidopsis: tools, traits and prospects for evolutionary ecology.

BACKGROUND: The model plant Arabidopsis thaliana (Arabidopsis) shows a wide range of genetic and trait variation among wild accessions. Because of its unparalleled biological and genomic resources, the potential of Arabidopsis for molecular genetic analysis of this natural variation has increased dramatically in recent years. SCOPE: Advanced genomics has accelerated molecular phylogenetic analysis and gene identification by quantitative trait loci (QTL) mapping and/or association mapping in Arabidopsis. In particular, QTL mapping utilizing natural accessions is now becoming a major strategy of gene isolation, offering an alternative to artificial mutant lines. Furthermore, the genomic information is used by researchers to uncover the signature of natural selection acting on the genes that contribute to phenotypic variation. The evolutionary significance of such genes has been evaluated in traits such as disease resistance and flowering time. However, although molecular hallmarks of selection have been found for the genes in question, a corresponding ecological scenario of adaptive evolution has been difficult to prove. Ecological strategies, including reciprocal transplant experiments and competition experiments, and utilizing near-isogenic lines of alleles of interest will be a powerful tool to measure the relative fitness of phenotypic and/or allelic variants. CONCLUSIONS: As the plant model organism, Arabidopsis provides a wealth of molecular background information for evolutionary genetics. Because genetic diversity between and within Arabidopsis populations is much higher than anticipated, combining this background information with ecological approaches might well establish Arabidopsis as a model organism for plant evolutionary ecology.

Temporal variation in glucocorticoid levels during the resting phase is associated in opposite way with maternal and paternal melanic coloration.

Sex-dependent selection can help maintain sexual dimorphism. When the magnitude of selection exerted on a heritable sex trait differs between the sexes, it may prevent each sex to reach its phenotypic optimum. As a consequence, the benefit of expressing a sex trait to a given value may differ between males and females favouring sex-specific adaptations associated with different values of a sex trait. The level of metabolites regulated by genes that are under sex-dependent selection may therefore covary with the degree of ornamentation differently in the two sexes. We investigated this prediction in the barn owl, a species in which females display on average larger black spots on the plumage than males, a heritable ornament. This melanin-based colour trait is strongly selected in females and weakly counter-selected in males indicating sex-dependent selection. In nestling barn owls, we found that daily variation in baseline corticosterone levels, a key hormone that mediates life history trade-offs, covaries with spot diameter displayed by their biological parents. When their mother displayed larger spots, nestlings had lower corticosterone levels in the morning and higher levels in the evening, whereas the opposite pattern was found with the size of paternal spots. Our study suggests a link between daily regulation of glucocorticoids and sex-dependent selection exerted on sexually dimorphic melanin-based ornaments.

Genetic variation and environmental effects on beta-conglycinin and glycinin content in Brazilian soybean cultivars

The objective of this work was to determine genetic and environmental effects on beta-conglycinin and glycinin content in Brazilian soybean cultivars. The concentrations of these protein fractions were analyzed by scanning densitometry after electrophoresis, in 90 Brazilian soybean cultivars sown in Ponta Grossa, PR, in 2001. The effects of the sowing location were determined in the cultivar MG/BR 46 (Conquista), sown in 16 locations of Goiás and Minas Gerais states (Central Brazil), and in the cultivar IAS 5, sown in 12 locations of Paraná and São Paulo states (Southern Brazil), in 2002 soybean season. A significant variability for beta-conglycinin (7S) and glycinin (11S) protein fractions ratio was observed among the 90 Brazilian soybean cultivars. 'MS/BRS 169' (Bacuri) and 'BR-8' (Pelotas) presented the highest and the lowest 11S/7S ratios (2.76 and 1.17, respectively). Beta-conglycinin protein fractions presented more variability than glycinin protein fractions. Grouping test classified 7S proteins in seven groups, 11S proteins in four groups, and protein fraction ratios (11S/7S) in nine groups. Significant effect of sowing locations was also observed on protein fractions contents. There is a good possibility of breeding for individual protein fractions, and their subunits, without affecting protein content.

Genetic variation in GIPR influences the glucose and insulin responses to an oral glucose challenge.

Glucose levels 2 h after an oral glucose challenge are a clinical measure of glucose tolerance used in the diagnosis of type 2 diabetes. We report a meta-analysis of nine genome-wide association studies (n = 15,234 nondiabetic individuals) and a follow-up of 29 independent loci (n = 6,958-30,620). We identify variants at the GIPR locus associated with 2-h glucose level (rs10423928, beta (s.e.m.) = 0.09 (0.01) mmol/l per A allele, P = 2.0 x 10(-15)). The GIPR A-allele carriers also showed decreased insulin secretion (n = 22,492; insulinogenic index, P = 1.0 x 10(-17); ratio of insulin to glucose area under the curve, P = 1.3 x 10(-16)) and diminished incretin effect (n = 804; P = 4.3 x 10(-4)). We also identified variants at ADCY5 (rs2877716, P = 4.2 x 10(-16)), VPS13C (rs17271305, P = 4.1 x 10(-8)), GCKR (rs1260326, P = 7.1 x 10(-11)) and TCF7L2 (rs7903146, P = 4.2 x 10(-10)) associated with 2-h glucose. Of the three newly implicated loci (GIPR, ADCY5 and VPS13C), only ADCY5 was found to be associated with type 2 diabetes in collaborating studies (n = 35,869 cases, 89,798 controls, OR = 1.12, 95% CI 1.09-1.15, P = 4.8 x 10(-18)).

The branch-site test of positive selection is surprisingly robust but lacks power under synonymous substitution saturation and variation in GC

Positive selection is widely estimated from protein coding sequence alignments by the nonsynonymous-to-synonymous ratio omega. Increasingly elaborate codon models are used in a likelihood framework for this estimation. Although there is widespread concern about the robustness of the estimation of the omega ratio, more efforts are needed to estimate this robustness, especially in the context of complex models. Here, we focused on the branch-site codon model. We investigated its robustness on a large set of simulated data. First, we investigated the impact of sequence divergence. We found evidence of underestimation of the synonymous substitution rate for values as small as 0.5, with a slight increase in false positives for the branch-site test. When dS increases further, underestimation of dS is worse, but false positives decrease. Interestingly, the detection of true positives follows a similar distribution, with a maximum for intermediary values of dS. Thus, high dS is more of a concern for a loss of power (false negatives) than for false positives of the test. Second, we investigated the impact of GC content. We showed that there is no significant difference of false positives between high GC (up to similar to 80%) and low GC (similar to 30%) genes. Moreover, neither shifts of GC content on a specific branch nor major shifts in GC along the gene sequence generate many false positives. Our results confirm that the branch-site is a very conservative test.