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.

Language variation in source texts and their translations: the case of L3 in film translation

In addition to the two languages essentially involved in translation, that of the source text (L1) and that of the target text (L2), we propose a third language (L3) to refer to any other language(s) found in the text. L3 may appear in the source text (ST) or the target text (TT), actually appearing more frequently inSTs in our case studies. We present a range of combinations for the convergence and divergence of L1, L2 and L3, for the case of feature films and their translations using examples from dubbed and subtitled versions of films, but we are hopeful that our tentative conclusions may be relevant to other modalities of translation, audiovisual and otherwise. When L3 appears in an audiovisual ST,we find a variety of solutions whereby L3 is deleted from or adapted to the TT.In the latter case, L3 might be rendered in a number of ways, depending on factors such as the audience’s familiarity with L3, and the possibility that L3 inthe ST is an invented language.

Variation of grain nitrogen content in relation with grain yield in old and modern Spanish wheats grown under a wide range of agronomic conditions in a Mediterranean region

Wheat yield and grain nitrogen concentration (GNC; mg N/g grain) are frequently negatively correlated. In most growing conditions, this is mainly due to a feedback process between GNC and the number of grains/m2. In Mediterranean conditions, breeders may have produced cultivars with conservative grain set. The present study aimed at clarifying the main physiological determinants of grain nitrogen accumulation (GNA) in Mediterranean wheat and to analyse how breeding has affected them. Five field experiments were carried out in north-eastern Spain in the 2005/06 and 2006/ 07 growing seasons with three cultivars released at different times and an advanced line. Depending on the experiment, source-sink ratios during grain filling were altered by reducing grain number/m2 either through pre-anthesis shading (unshaded control or 0.75 shading only between jointing and anthesis) or by directly trimming the spikes after anthesis and before the onset of the effective grain filling period (un-trimmed control or spikes halved 7–10 days after anthesis). Grain nitrogen content (GN content ; mg N/grain) decreased with the year of release of the genotypes. As the number of grains/m2 was also increased by breeding there was a clear dilution effect on the amount of nitrogen allocated to each grain. However, the increase in GN content in old genotypes did not compensate for the loss in grain nitrogen yield (GNY) due to the lower number of grains/m2. GN content of all genotypes increased (increases ranged from 0.13 to 0.40 mg N/grain, depending on experiment and genotype) in response to the post-anthesis spike trimming or pre-anthesis shading. The degree of source-limitation for GNA increased with the year of release of the genotypes (and thus with increases in grain number/m2) from 0.22 (mean of the four manipulative experiments) in the oldest cultivar to 0.51 (mean of the four manipulative experiments) in the most modern line. It was found that final GN content depended strongly on the source-sink ratio established at anthesis between the number of grains set and the amount of nitrogen absorbed at this stage. Thus, Mediterranean wheat breeding that improved yield through increases in grain number/m2 reduced the GN content by diluting a rather limited source of nitrogen into more grains. This dilution effect produced by breeding was further confirmed by the reversal effect produced by grain number/m2 reductions due to either pre-anthesis shading or post-anthesis spike trimming.

Genetic variation for carbon isotope composition in Juglans regia L.: relationships with growth, phenology and climate of origin

Among the traits of breeding interest for the common walnut tree Juglans regia L., characteristics such as timing of budbreak and leaf fall, water-use efficiency and growth performance are regarded as being of utmost relevance in Mediterranean conditions. The authors evaluated intraspecific variation in $\delta$13C (carbon isotope composition, surrogate of intrinsic water-use efficiency, WUE$_{\rm i}$) for 22 J. regia families grown in a progeny test under supplementary irrigation, and investigated whether such variation correlated with climatic indicators of native habitats. The genetic relationships between $\delta$13C, growth and phenology were also assessed during two consecutive years. Overall, the most water-use-efficient families (i.e. with higher $\delta$13C), which originated mainly from drought-prone provenance regions which have a high vapour pressure deficit and low rainfall, exhibited less height growth and smaller DBH. Using a stepwise regression procedure, $\delta$13C was included as the main explanatory variable of genotypic variation in growth traits, together with growing season duration (for DBH in both years) and flushing (for height in 2007). It was concluded that WUE$_{\rm i}$ is largely unconnected to phenology effects in the explanation of growth performance for J. regia, therefore suggesting the opportunity of simultaneously selecting for low WUE$_{\rm i}$ and extended growing period to maximise productivity in non-water-limited environments.

Copy number variation in the porcine genome inferred from a 60 k SNP BeadChip

Background: Recent studies in pigs have detected copy number variants (CNVs) using the Comparative Genomic Hybridization technique in arrays designed to cover specific porcine chromosomes. The goal of this study was to identify CNV regions (CNVRs) in swine species based on whole genome SNP genotyping chips. Results: We used predictions from three different programs (cnvPartition, PennCNV and GADA) to analyze data from the Porcine SNP60 BeadChip. A total of 49 CNVRs were identified in 55 animals from an Iberian x Landrace cross (IBMAP) according to three criteria: detected in at least two animals, contained three or more consecutive SNPs and recalled by at least two programs. Mendelian inheritance of CNVRs was confirmed in animals belonging to several generations of the IBMAP cross. Subsequently, a segregation analysis of these CNVRs was performed in 372 additional animals from the IBMAP cross and its distribution was studied in 133 unrelated pig samples from different geographical origins. Five out of seven analyzed CNVRs were validated by real time quantitative PCR, some of which coincide with well known examples of CNVs conserved across mammalian species. Conclusions: Our results illustrate the usefulness of Porcine SNP60 BeadChip to detect CNVRs and show that structural variants can not be neglected when studying the genetic variability in this species.

Gene flow and gene flux shape evolutionary patterns of variation in Drosophila subobscura

Gene flow (defined as allele exchange between populations) and gene flux (defined as allele exchange during meiosis in heterokaryotypic females) are important factors decreasing genetic differentiation between populations and inversions. Many chromosomal inversions are under strong selection and their role in recombination reduction enhances the maintenance of their genetic distinctness. Here we analyze levels and patterns of nucleotide diversity, selection and demographic history, using 37 individuals of Drosophila subobscura from Mount Parnes (Greece) and Barcelona (Spain). Our sampling focused on two frequent O-chromosome arrangements that differ by two overlapping inversions (OST and O3+4), which are differentially adapted to the environment as observed by their opposing latitudinal clines in inversion frequencies. The six analyzed genes (Pif1A, Abi, Sqd, Yrt, Atpa and Fmr1) were selected for their location across the O-chromosome and their implication in thermal adaptation. Despite the extensive gene flux detected outside the inverted region, significant genetic differentiation between both arrangements was found inside it. However, high levels of gene flow were detected for all six genes when comparing the same arrangement among populations. These results suggest that the adaptive value of inversions is maintained, regardless of the lack of genetic differentiation within arrangements from different populations, and thus favors the Local Adaptation hypothesis over the Coadapted Genome hypothesis as the basis of the selection acting on inversions in these populations.

Nuclear DNA content variation in life history phases of the Bonnemaisoniaceae (Rhodophyta)

Nuclear DNA content in gametophytes and sporophytes or the prostrate phases of the following species of Bonnemaisoniaceae (Asparagopsis armata, Asparagopsis taxiformis, Bonnemaisonia asparagoides, Bonnemaisonia clavata and Bonnemaisonia hamifera) were estimated by image analysis and static microspectrophotometry using the DNA-localizing fluorochrome DAPI (4′, 6-diamidino-2-phenylindole, dilactate) and the chicken erythrocytes standard. These estimates expand on the Kew database of DNA nuclear content. DNA content values for 1C nuclei in the gametophytes (spermatia and vegetative cells) range from 0.5 pg to 0.8 pg, and for 2C nuclei in the sporophytes or the prostrate phases range from 1.15-1.7 pg. Although only the 2C and 4C values were observed in the sporophyte or the prostrate phase, in the vegetative cells of the gametophyte the values oscillated from 1C to 4C, showing the possible start of endopolyploidy. The results confirm the alternation of nuclear phases in these Bonnemaisoniaceae species, in those that have tetrasporogenesis, as well as those that have somatic meiosis. The availability of a consensus phylogenetic tree for Bonnemaisoniaceae has opened the way to determine evolutionary trends in DNA contents. Both the estimated genome sizes and the published chromosome numbers for Bonnemaisoniaceae suggest a narrow range of values consistent with the conservation of an ancestral genome.

Genetic Variation in the TP53 Pathway and Bladder Cancer Risk. A Comprehensive Analysis

Introduction: Germline variants in TP63 have been consistently associated with several tumors, including bladder cancer, indicating the importance of TP53 pathway in cancer genetic susceptibility. However, variants in other related genes, including TP53 rs1042522 (Arg72Pro), still present controversial results. We carried out an in depth assessment of associations between common germline variants in the TP53 pathway and bladder cancer risk. Material and Methods: We investigated 184 tagSNPs from 18 genes in 1,058 cases and 1,138 controls from the Spanish Bladder Cancer/EPICURO Study. Cases were newly-diagnosed bladder cancer patients during 1998–2001. Hospital controls were age-gender, and area matched to cases. SNPs were genotyped in blood DNA using Illumina Golden Gate and TaqMan assays. Cases were subphenotyped according to stage/grade and tumor p53 expression. We applied classical tests to assess individual SNP associations and the Least Absolute Shrinkage and Selection Operator (LASSO)-penalized logistic regression analysis to assess multiple SNPs simultaneously. Results: Based on classical analyses, SNPs in BAK1 (1), IGF1R (5), P53AIP1 (1), PMAIP1 (2), SERINPB5 (3), TP63 (3), and TP73 (1) showed significant associations at p-value#0.05. However, no evidence of association, either with overall risk or with specific disease subtypes, was observed after correction for multiple testing (p-value$0.8). LASSO selected the SNP rs6567355 in SERPINB5 with 83% of reproducibility. This SNP provided an OR = 1.21, 95%CI 1.05–1.38, p-value = 0.006, and a corrected p-value = 0.5 when controlling for over-estimation. Discussion: We found no strong evidence that common variants in the TP53 pathway are associated with bladder cancer susceptibility. Our study suggests that it is unlikely that TP53 Arg72Pro is implicated in the UCB in white Europeans. SERPINB5 and TP63 variation deserve further exploration in extended studies.

Genetic variation in the European hake, Merluccius merluccius. Description of protein loci and genetic divergence between Atlantic and Mediterranean populations = Anàlisi de la variabilitat genètica en el lluç europeu, Merluccius merluccius. Descripció dels loci electroforètics i divergència genètica entre les poblacions atlàntica i mediterrània

We examined the genetic population structure of the european hake (Merluccius merluccius) using electrophoretically detectable population markers in 35 protein loci. Samples were collected from 7 locations in the Atlantic Ocean and Mediterranean Sea. Six loci were polymorphic using the 0.05 criterion of polymorphism. Sample heterozigosities ranged from 0.052 to 0.072 and averaged 0.0625. In this study, significant allele frequency differences were detected between Atlantic and Mediterranean populations in three polymorphic loci: GAPDH-1*, GPI-2* and SOD-1*. Two major genetic groups were considered: a North-Atlantic stock and the Mediterranean stock. The Nei genetic distance, D, (based on 33 loci) between samples from these two groups ranged from 0.002 to 0.006. Genetic differenciation between these areas appears to reflect the barrier effect of Strait of Gibraltar. On average over loci, 96.92 % of the total gene diversity was contained within samples, 0.23 % expressed differences among locations within areas, and 2.64 % differences between regions. A review of morphological variation together with the genetic data presented here suggest that the populations of hake from these areas are subdivided into two different stocks: the North-Atlantic stock and the Mediterranean stock. The most conservative approach to the management of these stocks is to consider the Atlantic and Mediterranean stocks independently from oneanother

Geographic patterns of genetic variation in a broadly distributed marine vertebrate: new insights into loggerhead turtle stock structure from expanded mitochondrial DNA sequences

Previous genetic studies have demonstrated that natal homing shapes the stock structure of marine turtle nesting populations. However, widespread sharing of common haplotypes based on short segments of the mitochondrial control region often limits resolution of the demographic connectivity of populations. Recent studies employing longer control region sequences to resolve haplotype sharing have focused on regional assessments of genetic structure and phylogeography. Here we synthesize available control region sequences for loggerhead turtles from the Mediterranean Sea, Atlantic, and western Indian Ocean basins. These data represent six of the nine globally significant regional management units (RMUs) for the species and include novel sequence data from Brazil, Cape Verde, South Africa and Oman. Genetic tests of differentiation among 42 rookeries represented by short sequences (380 bp haplotypes from 3,486 samples) and 40 rookeries represented by long sequences (~800 bp haplotypes from 3,434 samples) supported the distinction of the six RMUs analyzed as well as recognition of at least 18 demographically independent management units (MUs) with respect to female natal homing. A total of 59 haplotypes were resolved. These haplotypes belonged to two highly divergent global lineages, with haplogroup I represented primarily by CC-A1, CC-A4, and CC-A11 variants and haplogroup II represented by CC-A2 and derived variants. Geographic distribution patterns of haplogroup II haplotypes and the nested position of CC-A11.6 from Oman among the Atlantic haplotypes invoke recent colonization of the Indian Ocean from the Atlantic for both global lineages. The haplotypes we confirmed for western Indian Ocean RMUs allow reinterpretation of previous mixed stock analysis and further suggest that contemporary migratory connectivity between the Indian and Atlantic Oceans occurs on a broader scale than previously hypothesized. This study represents a valuable model for conducting comprehensive international cooperative data management and research in marine ecology.

Sexual dimorphism of sonic apparatus and extreme intersexual variation of sounds in Ophidion rochei (Ophidiidae): first evidence of a tight relationship between morphology and sound characteristics in Ophidiidae

Many Ophidiidae are active in dark environments and display complex sonic apparatus morphologies. However, sound recordings are scarce and little is known about acoustic communication in this family. This paper focuses on Ophidion rochei which is known to display an important sexual dimorphism in swimbladder and anterior skeleton. The aims of this study were to compare the sound producing morphology, and the resulting sounds in juveniles, females and males of O. rochei.Results: Males, females, and juveniles possessed different morphotypes. Females and juveniles contrasted with males because they possessed dramatic differences in morphology of their sonic muscles, swimbladder, supraoccipital crest, and first vertebrae and associated ribs. Further, they lacked the ‘rocker bone’ typically found in males. Sounds from each morphotype were highly divergent. Males generally produced non harmonic, multiple-pulsed sounds that lasted for several seconds (3.5 ± 1.3 s) with a pulse period of ca. 100 ms. Juvenile and female sounds were recorded for the first time in ophidiids. Female sounds were harmonic, had shorter pulse period (±3.7 ms), and never exceeded a few dozen milliseconds (18 ± 11 ms). Moreover, unlike male sounds, female sounds did not have alternating long and short pulse periods. Juvenile sounds were weaker but appear to be similar to female sounds.Conclusions: Although it is not possible to distinguish externally male from female in O. rochei, they show a sonic apparatus and sounds that are dramatically different. This difference is likely due to their nocturnal habits that may have favored the evolution of internal secondary sexual characters that help to distinguish males from females and that could facilitate mate choice by females. Moreover, the comparison of different morphotypes in this study shows that these morphological differences result from a peramorphosis that takes place during the development of the gonads

Improved understanding of drought controls on seasonal variation in Mediterranean forest canopy CO2 and water fluxes through combined in situ measurements and ecosystem modelling

Water stress is a defining characteristic of Mediterranean ecosystems, and is likely to become more severe in the coming decades. Simulation models are key tools for making predictions, but our current understanding of how soil moisture controls ecosystem functioning is not sufficient to adequately constrain parameterisations. Canopy-scale flux data from four forest ecosystems with Mediterranean-type climates were used in order to analyse the physiological controls on carbon and water flues through the year. Significant non-stomatal limitations on photosynthesis were detected, along with lesser changes in the conductance-assimilation relationship. New model parameterisations were derived and implemented in two contrasting modelling approaches. The effectiveness of two models, one a dynamic global vegetation model ('ORCHIDEE'), and the other a forest growth model particularly developed for Mediterranean simulations ('GOTILWA+'), was assessed and modelled canopy responses to seasonal changes in soil moisture were analysed in comparison with in situ flux measurements. In contrast to commonly held assumptions, we find that changing the ratio of conductance to assimilation under natural, seasonally-developing, soil moisture stress is not sufficient to reproduce forest canopy CO2 and water fluxes. However, accurate predictions of both CO2 and water fluxes under all soil moisture levels encountered in the field are obtained if photosynthetic capacity is assumed to vary with soil moisture. This new parameterisation has important consequences for simulated responses of carbon and water fluxes to seasonal soil moisture stress, and should greatly improve our ability to anticipate future impacts of climate changes on the functioning of ecosystems in Mediterranean-type climates.

Environmental factors, spatial variation, and specific requirements of Chironomidae in Mediterranean reference streams

Chironomidae spatial distribution was investigated at 63 near-pristine sites in 22 catchments of the Iberian Mediterranean coast. We used partial redundancy analysis to study Chironomidae community responses to a number of environmental factors acting at several spatial scales. The percentage of variation explained by local factors (23.3%) was higher than that explained by geographical (8.5%) or regional factors(8%). Catchment area, longitude, pH, % siliceous rocks in the catchment, and altitude were the best predictors of Chironomidae assemblages. We used a k-means cluster analysis to classified sites into 3 major groups based on Chironomidae assemblages. These groups were explained mainly by longitudinal zonation and geographical position, and were defined as 1) siliceous headwater streams, 2) mid-altitude streams with small catchment areas, and 3) medium-sized calcareous streams. Distinct species assemblages with associated indicator taxa were established for each stream category using IndVal analysis. Species responses to previously identified key environmental variables were determined, and optima and tolerances were established by weighted average regression. Distinct ecological requirements were observed among genera and among species of the same genus. Some genera were restricted to headwater systems (e.g., Diamesa), whereas others (e.g., Eukiefferiella) had wider ecological preferences but with distinct distributions among congenerics. In the present period of climate change, optima and tolerances of species might be a useful tool to predict responses of different species to changes in significant environmental variables, such as temperature and hydrology.