New Features of CoDaPack. An Userfriendly Compositional Data Package

The statistical analysis of compositional data is commonly used in geological studies.As is well-known, compositions should be treated using logratios of parts, which aredifficult to use correctly in standard statistical packages. In this paper we describe thenew features of our freeware package, named CoDaPack, which implements most of thebasic statistical methods suitable for compositional data. An example using real data ispresented to illustrate the use of the package

Complete longitudinal analyses of the randomized, placebo-controlled, phase III trial of sunitinib in patients with gastrointestinal stromal tumor following imatinib failure.

PURPOSE: To analyze final long-term survival and clinical outcomes from the randomized phase III study of sunitinib in gastrointestinal stromal tumor patients after imatinib failure; to assess correlative angiogenesis biomarkers with patient outcomes. EXPERIMENTAL DESIGN: Blinded sunitinib or placebo was given daily on a 4-week-on/2-week-off treatment schedule. Placebo-assigned patients could cross over to sunitinib at disease progression/study unblinding. Overall survival (OS) was analyzed using conventional statistical methods and the rank-preserving structural failure time (RPSFT) method to explore cross-over impact. Circulating levels of angiogenesis biomarkers were analyzed. RESULTS: In total, 243 patients were randomized to receive sunitinib and 118 to placebo, 103 of whom crossed over to open-label sunitinib. Conventional statistical analysis showed that OS converged in the sunitinib and placebo arms (median 72.7 vs. 64.9 weeks; HR, 0.876; P = 0.306) as expected, given the cross-over design. RPSFT analysis estimated median OS for placebo of 39.0 weeks (HR, 0.505, 95% CI, 0.262-1.134; P = 0.306). No new safety concerns emerged with extended sunitinib treatment. No consistent associations were found between the pharmacodynamics of angiogenesis-related plasma proteins during sunitinib treatment and clinical outcome. CONCLUSIONS: The cross-over design provided evidence of sunitinib clinical benefit based on prolonged time to tumor progression during the double-blind phase of this trial. As expected, following cross-over, there was no statistical difference in OS. RPSFT analysis modeled the absence of cross-over, estimating a substantial sunitinib OS benefit relative to placebo. Long-term sunitinib treatment was tolerated without new adverse events.

CODAPACK3D. A new version of Compositional Data Package

The statistical analysis of compositional data should be treated using logratios of parts,which are difficult to use correctly in standard statistical packages. For this reason afreeware package, named CoDaPack was created. This software implements most of thebasic statistical methods suitable for compositional data.In this paper we describe the new version of the package that now is calledCoDaPack3D. It is developed in Visual Basic for applications (associated with Excel©),Visual Basic and Open GL, and it is oriented towards users with a minimum knowledgeof computers with the aim at being simple and easy to use.This new version includes new graphical output in 2D and 3D. These outputs could bezoomed and, in 3D, rotated. Also a customization menu is included and outputs couldbe saved in jpeg format. Also this new version includes an interactive help and alldialog windows have been improved in order to facilitate its use.To use CoDaPack one has to access Excel© and introduce the data in a standardspreadsheet. These should be organized as a matrix where Excel© rows correspond tothe observations and columns to the parts. The user executes macros that returnnumerical or graphical results. There are two kinds of numerical results: new variablesand descriptive statistics, and both appear on the same sheet. Graphical output appearsin independent windows. In the present version there are 8 menus, with a total of 38submenus which, after some dialogue, directly call the corresponding macro. Thedialogues ask the user to input variables and further parameters needed, as well aswhere to put these results. The web site http://ima.udg.es/CoDaPack contains thisfreeware package and only Microsoft Excel© under Microsoft Windows© is required torun the software.Kew words: Compositional data Analysis, Software

Probability distribution of the radionuclide half-lives

Power law distributions, a well-known model in the theory of real random variables, characterize a wide variety of natural and man made phenomena. The intensity of earthquakes, the word frequencies, the solar ares and the sizes of power outages are distributed according to a power law distribution. Recently, given the usage of power laws in the scientific community, several articles have been published criticizing the statistical methods used to estimate the power law behaviour and establishing new techniques to their estimation with proven reliability. The main object of the present study is to go in deep understanding of this kind of distribution and its analysis, and introduce the half-lives of the radioactive isotopes as a new candidate in the nature following a power law distribution, as well as a \canonical laboratory" to test statistical methods appropriate for long-tailed distributions.

Unusual voltage sag event detection in power systems

Three multivariate statistical tools (principal component analysis, factor analysis, analysis discriminant) have been tested to characterize and model the sags registered in distribution substations. Those models use several features to represent the magnitude, duration and unbalanced grade of sags. They have been obtained from voltage and current waveforms. The techniques are tested and compared using 69 registers of sags. The advantages and drawbacks of each technique are listed

Clustering compositional data trajectories

Our essay aims at studying suitable statistical methods for the clustering ofcompositional data in situations where observations are constituted by trajectories ofcompositional data, that is, by sequences of composition measurements along a domain.Observed trajectories are known as “functional data” and several methods have beenproposed for their analysis.In particular, methods for clustering functional data, known as Functional ClusterAnalysis (FCA), have been applied by practitioners and scientists in many fields. To ourknowledge, FCA techniques have not been extended to cope with the problem ofclustering compositional data trajectories. In order to extend FCA techniques to theanalysis of compositional data, FCA clustering techniques have to be adapted by using asuitable compositional algebra.The present work centres on the following question: given a sample of compositionaldata trajectories, how can we formulate a segmentation procedure giving homogeneousclasses? To address this problem we follow the steps described below.First of all we adapt the well-known spline smoothing techniques in order to cope withthe smoothing of compositional data trajectories. In fact, an observed curve can bethought of as the sum of a smooth part plus some noise due to measurement errors.Spline smoothing techniques are used to isolate the smooth part of the trajectory:clustering algorithms are then applied to these smooth curves.The second step consists in building suitable metrics for measuring the dissimilaritybetween trajectories: we propose a metric that accounts for difference in both shape andlevel, and a metric accounting for differences in shape only.A simulation study is performed in order to evaluate the proposed methodologies, usingboth hierarchical and partitional clustering algorithm. The quality of the obtained resultsis assessed by means of several indices

Networks of PhD Students and Academic Performance: A Comparison across Countries

In this article we compare regression models obtained to predict PhD students’ academic performance in the universities of Girona (Spain) and Slovenia. Explanatory variables are characteristics of PhD student’s research group understood as an egocentered social network, background and attitudinal characteristics of the PhD students and some characteristics of the supervisors. Academic performance was measured by the weighted number of publications. Two web questionnaires were designed, one for PhD students and one for their supervisors and other research group members. Most of the variables were easily comparable across universities due to the careful translation procedure and pre-tests. When direct comparison was notpossible we created comparable indicators. We used a regression model in which the country was introduced as a dummy coded variable including all possible interaction effects. The optimal transformations of the main and interaction variables are discussed. Some differences between Slovenian and Girona universities emerge. Some variables like supervisor’s performance and motivation for autonomy prior to starting the PhD have the same positive effect on the PhD student’s performance in both countries. On the other hand, variables like too close supervision by the supervisor and having children have a negative influence in both countries. However, we find differences between countries when we observe the motivation for research prior to starting the PhD which increases performance in Slovenia but not in Girona. As regards network variables, frequency of supervisor advice increases performance in Slovenia and decreases it in Girona. The negative effect in Girona could be explained by the fact that additional contacts of the PhD student with his/her supervisor might indicate a higher workload in addition to or instead of a better advice about the dissertation. The number of external student’s advice relationships and social support mean contact intensity are not significant in Girona, but they have a negative effect in Slovenia. We might explain the negative effect of external advice relationships in Slovenia by saying that a lot of external advice may actually result from a lack of the more relevant internal advice

Genetics of cadmium hyperaccumulation variance in Thlaspi caerulescens natural populations from Switzerland

SUMMARY Heavy metal presence in the environment is a serious concern since some of them can be toxic to plants, animals and humans once accumulated along the food chain. Cadmium (Cd) is one of the most toxic heavy metal. It is naturally present in soils at various levels and its concentration can be increased by human activities. Several plants however have naturally developed strategies allowing them to grow on heavy metal enriched soils. One of them consists in the accumulation and sequestration of heavy metals in the above-ground biomass. Some plants present in addition an extreme strategy by which they accumulate a limited number of heavy metals in their shoots in amounts 100 times superior to those expected for a non-accumulating plant in the same conditions. Understanding the genetic basis of the hyperaccumulation trait - particularly for Cd - remains an important challenge which may lead to biotechnological applications in the soil phytoremediation. In this thesis, Thlaspi caerulescens J. & C. Presl (Brassicaceae) was used as a model plant to study the Cd hyperaccumulation trait, owing to its physiological and genetic characteristics. Twenty-four wild populations were sampled in different regions of Switzerland. They were characterized for environmental and soil parameters as well as intrinsic characteristics of plants (i.e. metal concentrations in shoots). They were as well genetically characterized by AFLPs, plastid DNA polymorphism and genes markers (CAPS and microsatellites) mainly developed in this thesis. Some of the investigated genes were putatively linked to the Cd hyperaccumulation trait. Since the study of the Cd hyperaccumulation in the field is important as it allows the identification of patterns of selection, the present work offered a methodology to define the Cd hyperaccumulation capacity of populations from different habitats permitting thus their comparison in the field. We showed that Cd, Zn, Fe and Cu accumulations were linked and that populations with higher Cd hyperaccumulation capacity had higher shoot and reproductive fitness. Using our genetic data, statistical methods (Beaumont & Nichols's procedure, partial Mantel tests) were applied to identify genomic signatures of natural selection related to the Cd hyperaccumulation capacity. A significant genetic difference between populations related to their Cd hyperaccumulation capacity was revealed based on somè specific markers (AFLP and candidate genes). Polymorphism at the gene encoding IRTl (Iron-transporter also participating to the transport of Zn) was suggested as explaining part of the variation in Cd hyperaccumulation capacity of populations supporting previous physiological investigations. RÉSUMÉ La présence de métaux lourds dans l'environnement est un phénomène préoccupant. En effet, certains métaux lourds - comme le cadmium (Cd) -sont toxiques pour les plantes, les animaux et enfin, accumulés le long de la chaîne alimentaire, pour les hommes. Le Cd est naturellement présent dans le sol et sa concentration peut être accrue par différentes activités humaines. Certaines plantes ont cependant développé des stratégies leur permettant de pousser sur des sols contaminés en métaux lourds. Parmi elles, certaines accumulent et séquestrent les métaux lourds dans leurs parties aériennes. D`autres présentent une stratégie encore plus extrême. Elles accumulent un nombre limité de métaux lourds en quantités 100 fois supérieures à celles attendues pour des espèces non-accumulatrices sous de mêmes conditions. La compréhension des bases génétiques de l'hyperaccumulation -particulièrement celle du Cd - représente un défi important avec des applications concrètes en biotechnologies, tout particulièrement dans le but appliqué de la phytoremediation des sols contaminés. Dans cette thèse, Thlaspi caerulescens J. & C. Presl (Brassicaceae) a été utilisé comme modèle pour l'étude de l'hyperaccumulation du Cd de par ses caractéristiques physiologiques et génétiques. Vingt-quatre populations naturelles ont été échantillonnées en Suisse et pour chacune d'elles les paramètres environnementaux, pédologique et les caractéristiques intrinsèques aux plantes (concentrations en métaux lourds) ont été déterminés. Les populations ont été caractérisées génétiquement par des AFLP, des marqueurs chloroplastiques et des marqueurs de gènes spécifiques, particulièrement ceux potentiellement liés à l'hyperaccumulation du Cd (CAPS et microsatellites). La plupart ont été développés au cours de cette thèse. L'étude de l'hyperaccumulation du Cd en conditions naturelles est importante car elle permet d'identifier la marque, éventuelle de sélection naturelle. Ce travail offre ainsi une méthodologie pour définir et comparer la capacité des populations à hyperaccumuler le Cd dans différents habitats. Nous avons montré que les accumulations du Cd, Zn, Fe et Cu sont liées et que les populations ayant une grande capacité d'hyperaccumuler le Cd ont également une meilleure fitness végétative et reproductive. Des méthodes statistiques (l'approche de Beaumont & Nichols, tests de Martel partiels) ont été utilisées sur les données génétiques pour identifier la signature génomique de la sélection naturelle liée à la capacité d'hyperaccumuler le Cd. Une différenciation génétique des populations liée à leur capacité d'hyperaccumuler le Cd a été mise en évidence sur certains marqueurs spécifiques. En accord avec les études physiologiques connues, le polymorphisme au gène codant IRT1 (un transporteur de Fe impliqué dans le transport du Zn) pourrait expliquer une partie de la variance de la capacité des populations à hyperaccumuler le Cd.

Interregional compensatory mechanisms of motor functioning in progressing preclinical neurodegeneration.

Understanding brain reserve in preclinical stages of neurodegenerative disorders allows determination of which brain regions contribute to normal functioning despite accelerated neuronal loss. Besides the recruitment of additional regions, a reorganisation and shift of relevance between normally engaged regions are a suggested key mechanism. Thus, network analysis methods seem critical for investigation of changes in directed causal interactions between such candidate brain regions. To identify core compensatory regions, fifteen preclinical patients carrying the genetic mutation leading to Huntington's disease and twelve controls underwent fMRI scanning. They accomplished an auditory paced finger sequence tapping task, which challenged cognitive as well as executive aspects of motor functioning by varying speed and complexity of movements. To investigate causal interactions among brain regions a single Dynamic Causal Model (DCM) was constructed and fitted to the data from each subject. The DCM parameters were analysed using statistical methods to assess group differences in connectivity, and the relationship between connectivity patterns and predicted years to clinical onset was assessed in gene carriers. In preclinical patients, we found indications for neural reserve mechanisms predominantly driven by bilateral dorsal premotor cortex, which increasingly activated superior parietal cortices the closer individuals were to estimated clinical onset. This compensatory mechanism was restricted to complex movements characterised by high cognitive demand. Additionally, we identified task-induced connectivity changes in both groups of subjects towards pre- and caudal supplementary motor areas, which were linked to either faster or more complex task conditions. Interestingly, coupling of dorsal premotor cortex and supplementary motor area was more negative in controls compared to gene mutation carriers. Furthermore, changes in the connectivity pattern of gene carriers allowed prediction of the years to estimated disease onset in individuals. Our study characterises the connectivity pattern of core cortical regions maintaining motor function in relation to varying task demand. We identified connections of bilateral dorsal premotor cortex as critical for compensation as well as task-dependent recruitment of pre- and caudal supplementary motor area. The latter finding nicely mirrors a previously published general linear model-based analysis of the same data. Such knowledge about disease specific inter-regional effective connectivity may help identify foci for interventions based on transcranial magnetic stimulation designed to stimulate functioning and also to predict their impact on other regions in motor-associated networks.

The COSMIN checklist for evaluating the methodological quality of studies on measurement properties: a clarification of its content

Background: The COSMIN checklist (COnsensus-based Standards for the selection of health status Measurement INstruments) was developed in an international Delphi study to evaluate the methodological quality of studies on measurement properties of health-related patient reported outcomes (HR-PROs). In this paper, we explain our choices for the design requirements and preferred statistical methods for which no evidence is available in the literature or on which the Delphi panel members had substantial discussion. Methods: The issues described in this paper are a reflection of the Delphi process in which 43 panel members participated. Results: The topics discussed are internal consistency (relevance for reflective and formative models, and distinction with unidimensionality), content validity (judging relevance and comprehensiveness), hypotheses testing as an aspect of construct validity (specificity of hypotheses), criterion validity (relevance for PROs), and responsiveness (concept and relation to validity, and (in) appropriate measures).Conclusions: We expect that this paper will contribute to a better understanding of the rationale behind the items, thereby enhancing the acceptance and use of the COSMIN checklist.

Bloodstream infections related to totally implantable venous access port: What is the situation in our hospital?

Port-a-Cath© (PAC) are totally implantable devices that offer an easy and long term access to venous circulation. They have been extensively used for intravenous therapy administration and are particularly well suited for chemotherapy in oncologic patients. Previous comparative studies have shown that these devices have the lowest catheter-related bloodstream infection rates among all intravascular access systems. However, bloodstream infection (BSI) still remains a major issue of port use and epidemiology data for PAC-associated BSI (PABSI) rates differ strongly depending on studies. Also, current literature about PABSI risk factors is scarce and sometimes controversial. Such heterogeneity may depend on type of studied population and local factors. Therefore, the aim of this study was to describe local epidemiology and risk factors for PABSI in adult patients in our tertiary- care university hospital. We conducted a retrospective cohort study in order to describe local epidemiology. We also performed a nested case-control study to identify local risk factors of PABSI. We analyzed medical files of adult patients who had a PAC implanted between January 1st, 2008 and December 31st, 2009 and looked for PABSI occurrence before May 1st, 2011 to define cases. Thirty nine PABSI occurred in this population with an attack rate of 5.8%. We estimated an incidence rate of 0.08/1000 PAC-days using the case-control study. PABSI causative agents were mainly Gram positive cocci (62%). We identified three predictive factors of PABSI by multivariate statistical analysis: neutropenia on outcome date (Odds Ratio [OR]: 4.05; 95% confidence interval [CI]:1.05- 15.66; p=0.042), diabetes (OR: 11.53; 95% CI: 1.07-124.70; p=0.044) and having another infection than PABSI on outcome date (OR: 6.35; 95% CI: 1.50-26.86; p=0.012). Patients suffering from acute or renal failure (OR: 4.26; 95% CI: 0.94-19.21; p=0.059) or wearing another invasive device (OR: 2.99; 95%CI:0.96-9.31; p=0.059) did not have a statistically increased risk for developing a PABSI according to classical threshold (p<0.05) but nevertheless remained close to significance. Our study demonstrated that local epidemiology and microbiology of PABSI in our institution was similar to previous reports. A larger prospective study is required to confirm our results or to test preventive measures.

Analysis of susceptibility to human immunodeficiency virus 1 (HIV-1) by comparative genetics

Summary [résumé français voir ci-dessous] From the beginning of the 20th century the world population has been confronted with the human immune deficiency virus 1 (HIV-1). This virus has the particularity to mutate fast, and could thus evade and adapt to the human host. Our closest evolutionary related organisms, the non-human primates, are less susceptible to HIV-1. In a broader sense, primates are differentially susceptible to various retrovirus. Species specificity may be due to genetic differences among primates. In the present study we applied evolutionary and comparative genetic techniques to characterize the evolutionary pattern of host cellular determinants of HIV-1 pathogenesis. The study of the evolution of genes coding for proteins participating to the restriction or pathogenesis of HIV-1 may help understanding the genetic basis of modern human susceptibility to infection. To perform comparative genetics analysis, we constituted a collection of primate DNA and RNA to allow generation of de novo sequence of gene orthologs. More recently, release to the public domain of two new primate complete genomes (bornean orang-utan and common marmoset) in addition of the three previously available genomes (human, chimpanzee and Rhesus monkey) help scaling up the evolutionary and comparative genome analysis. Sequence analysis used phylogenetic and statistical methods for detecting molecular adaptation. We identified different selective pressures acting on host proteins involved in HIV-1 pathogenesis. Proteins with HIV-1 restriction properties in non-human primates were under strong positive selection, in particular in regions of interaction with viral proteins. These regions carried key residues for the antiviral activity. Proteins of the innate immunity presented an evolutionary pattern of conservation (purifying selection) but with signals of relaxed constrain if we compared them to the average profile of purifying selection of the primate genomes. Large scale analysis resulted in patterns of evolutionary pressures according to molecular function, biological process and cellular distribution. The data generated by various analyses served to guide the ancestral reconstruction of TRIM5a a potent antiviral host factor. The resurrected TRIM5a from the common ancestor of Old world monkeys was effective against HIV-1 and the recent resurrected hominoid variants were more effective against other retrovirus. Thus, as the result of trade-offs in the ability to restrict different retrovirus, human might have been exposed to HIV-1 at a time when TRIM5a lacked the appropriate specific restriction activity. The application of evolutionary and comparative genetic tools should be considered for the systematical assessment of host proteins relevant in viral pathogenesis, and to guide biological and functional studies. Résumé La population mondiale est confrontée depuis le début du vingtième siècle au virus de l'immunodéficience humaine 1 (VIH-1). Ce virus a un taux de mutation particulièrement élevé, il peut donc s'évader et s'adapter très efficacement à son hôte. Les organismes évolutivement le plus proches de l'homme les primates nonhumains sont moins susceptibles au VIH-1. De façon générale, les primates répondent différemment aux rétrovirus. Cette spécificité entre espèces doit résider dans les différences génétiques entre primates. Dans cette étude nous avons appliqué des techniques d'évolution et de génétique comparative pour caractériser le modèle évolutif des déterminants cellulaires impliqués dans la pathogenèse du VIH- 1. L'étude de l'évolution des gènes, codant pour des protéines impliquées dans la restriction ou la pathogenèse du VIH-1, aidera à la compréhension des bases génétiques ayant récemment rendu l'homme susceptible. Pour les analyses de génétique comparative, nous avons constitué une collection d'ADN et d'ARN de primates dans le but d'obtenir des nouvelles séquences de gènes orthologues. Récemment deux nouveaux génomes complets ont été publiés (l'orang-outan du Bornéo et Marmoset commun) en plus des trois génomes déjà disponibles (humain, chimpanzé, macaque rhésus). Ceci a permis d'améliorer considérablement l'étendue de l'analyse. Pour détecter l'adaptation moléculaire nous avons analysé les séquences à l'aide de méthodes phylogénétiques et statistiques. Nous avons identifié différentes pressions de sélection agissant sur les protéines impliquées dans la pathogenèse du VIH-1. Des protéines avec des propriétés de restriction du VIH-1 dans les primates non-humains présentent un taux particulièrement haut de remplacement d'acides aminés (sélection positive). En particulier dans les régions d'interaction avec les protéines virales. Ces régions incluent des acides aminés clé pour l'activité de restriction. Les protéines appartenant à l'immunité inné présentent un modèle d'évolution de conservation (sélection purifiante) mais avec des traces de "relaxation" comparé au profil général de sélection purifiante du génome des primates. Une analyse à grande échelle a permis de classifier les modèles de pression évolutive selon leur fonction moléculaire, processus biologique et distribution cellulaire. Les données générées par les différentes analyses ont permis la reconstruction ancestrale de TRIM5a, un puissant facteur antiretroviral. Le TRIM5a ressuscité, correspondant à l'ancêtre commun entre les grands singes et les groupe des catarrhiniens, est efficace contre le VIH-1 moderne. Les TRIM5a ressuscités plus récents, correspondant aux ancêtres des grands singes, sont plus efficaces contre d'autres rétrovirus. Ainsi, trouver un compromis dans la capacité de restreindre différents rétrovirus, l'homme aurait été exposé au VIH-1 à une période où TRIM5a manquait d'activité de restriction spécifique contre celui-ci. L'application de techniques d'évolution et de génétique comparative devraient être considérées pour l'évaluation systématique de protéines impliquées dans la pathogenèse virale, ainsi que pour guider des études biologiques et fonctionnelles

Morphological studies on castes of Protopolybia chartergoides (Hymenoptera, Vespidae, Epiponini) observed in colonies during male production stage

Social wasps are important for the study of social behavior evolution because their colonies present different degrees of caste differentiation, from morphologically similar to highly conspicuous, associated with the occasional presence of intermediate females, which bears developed ovaries but no insemination. In the Polistinae, depending on the taxon, such differentiation can be discreet or conspicuous. This work intended to study morphological and physiological differences between castes in Protopolybia chartergoides by using morphometric analyses associated with multivariate statistical analyses and physiological evaluations from females' ovarian development. Results evidence low morphological and physiological differences among the castes in P. chartergoides, indicating three groups of females: queens, workers and intermediates. In this way, it was possible to suggest that Protopolybia chartergoides presents post-imaginal caste differentiation (or a very subtle form of pre-imaginal determination).

Dinámica espacio temporal en la ecología trófica de la anchoveta (Engraulis ringens Jenyns, 1842; pisces: engraulide) en el mar peruano

El objetivo de este trabajo es estudiar las variaciones ontogénicas, espaciales y nictemerales en la dieta de la anchoveta peruana (Engraulis ringens). El trabajo se basa en el análisis de contenido estomacal de 13 430 anchovetas (peces de 4 a 19,5 cm) muestreados durante 22 cruceros de investigación realizados a bordo de las embarcaciones científicas del Instituto del Mar del Perú durante el período 1996-2003. Previo a los análisis estadísticos la contribución de cada grupo de presa fue convertido a carbonó. Los datos fueron procesados por medio de estadística uni y multivariada en función de la talla de la anchoveta, la latitud, la distancia de la costa, la distancia del borde de la plataforma continental y el período nictemeral. Los resultados enseñan que la anchoveta se alimenta mayormente de zooplancton cualquier sea su tamaño o su distribución espaciotemporal. La anchoveta cambio gradualmente su dieta a lo largo de su ciclo de vida. Así, los individuos mas pequeños han ingerido mayormente copépodos calanoideos (68.0%) y eufausidos (23.3%), mientras que los individuos mas grandes ingirieron mayormente eufausidos (84.9%) y copépodos calanoideos (11.5%). Latitudinalmente, la mayor contribución de los copépodos calanoideos se observó en las zonas de afloramiento de Chicama y San Juan, con el 52.3% y 47.5%, respectivamente. Hasta los 160 km de distancia de la costa, no se ha observado un patrón regular en la contribución de las principales presas, mientras que por fuera, se insinúa un patrón caracterizado por el predominio de los eufausidos (>84%) y la menor contribución de copépodos calanoideos (<14,1%). En la plataforma continental los eufausidos contribuyeron con el 55.6% y los copépodos calanoideos con el 36.9%; mientras que por fuera de esa zona, los eufausidos se incrementaron al 75.3% y los copépodos calanoideos disminuyeron al 20.1%. Al considerar el período nictemeral, el ítem presa predominante durante el día fueron los eufausidos (52.7%) seguido de los copépodos calanoideos (40.9%); mientras que durante la noche, los eufausidos se incrementaron al 80.2% y los copépodos calanoideos disminuyeron al 15.7%. Estos resultados contradicen lo que se conocía acerca de la variabilidad de la dieta de la anchoveta, desestimándose: (i) que en el norte es fitoplanctívora y en el sur es zooplanctívora; (ii) que la fracción zooplanctónica es mayor en la zona oceánica. Asimismo, se ha demostrado que: (i) la anchoveta es principalmente zooplanctófaga durante su ciclo de vida, con modificaciones que corresponden a la mayor ingesta de presas pequeñas cuando son juveniles y de presas grandes cuando son adultas, (ii) en la zona de influencia de la plataforma continental consume aproximadamente el doble de copépodos calanoideos que en la zona oceánica, y aproximadamente una tercera parte menos de eufausidos en la zona interna que en la zona oceánica, y (iii) durante la noche y el día ingiere principalmente eufausidos pero la contribución de los copépodos calanoideos sube al 40%.

Trophic ecology of jumbo squid and predatory fishes in the Northern Humboldt Current System

This work provides a contribution to a better understanding of the trophic ecology of important predators in the Northern Humboldt Current System, the jack mackerel (Trachurus murphyi), the chub mackerel (Scomber japonicus) and the jumbo squid (Dosidicus gigas) by the characterization of the highly variable feeding patterns of these species at different spatiotemporal scales. We provided new knowledge on the comparative trophic behaviour of these species, defined as opportunistic in previous investigations. For that purpose we applied a variety of statistical methods to an extensive dataset of 27,188 non-empty stomachs. We defined the spatial organization of the forage fauna of these predators and documented changes in prey composition according to predators’ size and spatiotemporal features of environment. Our results highligh the key role played by the dissolved oxygen. We also deciphered an important paradox on the jumbo squid diet: why do they hardly forage on the huge anchovy (Engraulis ringens) biomass distributed of coastal Peru? We showed that the shallow oxygen minimum zone present off coastal Peru could hamper the co-occurrence of jumbo squids and anchovies. In addition, we proposed a conceptual model on jumbo squid trophic ecology including the ontogenetic cycle, oxygen and prey availability. Moreover we showed that the trophic behaviour of jack mackerel and chub mackerel is adapted to forage on more accessible species such as for example the squat lobster Pleurocondes monodon and Zoea larvae. Besides, both predators present a trophic overlap. But jack mackerel was not as oracious as chub mackerel, contradictorily to what was observed by others authors. Fish diet presented a high spatiotemporal variability, and the shelf break appeared as a strong biogeographical frontier. Diet composition of our fish predators was not necessarily a consistent indicator of changes in prey biomass. El Niño events had a weak effect on the stomach fullness and diet composition of chub mackerel and jack mackerel. Moreover, decadal changes in diet diversity challenged the classic paradigm of positive correlation between species richness and temperature. Finally, the global patterns that we described in this work, illustrated the opportunistic foraging behaviour, life strategies and the high degree of plasticity of these species. Such behaviour allows adaptation to changes in the environment.