Statistical Methodological Issues in Handling of Fatty Acid Data: Percentage or Concentration, Imputation and Indices

Basic aspects in the handling of fatty acid-data have remained largely underexposed. Of these, we aimed to address three statistical methodological issues, by quantitatively exemplifying their imminent confounding impact on analytical outcomes: (1) presenting results as relative percentages or absolute concentrations, (2) handling of missing/non-detectable values, and (3) using structural indices for data-reduction. Therefore, we reanalyzed an example dataset containing erythrocyte fatty acid-concentrations of 137 recurrently depressed patients and 73 controls. First, correlations between data presented as percentages and concentrations varied for different fatty acids, depending on their correlation with the total fatty acid-concentration. Second, multiple imputation of non-detects resulted in differences in significance compared to zero-substitution or omission of non-detects. Third, patients’ chain length-, unsaturation-, and peroxidation-indices were significantly lower compared to controls, which corresponded with patterns interpreted from individual fatty acid tests. In conclusion, results from our example dataset show that statistical methodological choices can have a significant influence on outcomes of fatty acid analysis, which emphasizes the relevance of: (1) hypothesis-based fatty acid-presentation (percentages or concentrations), (2) multiple imputation, preventing bias introduced by non-detects; and (3) the possibility of using (structural) indices, to delineate fatty acid-patterns thereby preventing multiple testing

A Targeted Constitutive Mutation in the Apc Tumor Suppressor Gene Underlies Mammary But Not Intestinal Tumorigenesis

Germline mutations in the adenomatous polyposis coli (APC) gene are responsible for familial adenomatous polyposis (FAP), an autosomal dominant hereditary predisposition to the development of multiple colorectal adenomas and of a broad spectrum of extra-intestinal tumors. Moreover, somatic APC mutations play a rate-limiting and initiating role in the majority of sporadic colorectal cancers. Notwithstanding its multifunctional nature, the main tumor suppressing activity of the APC gene resides in its ability to regulate Wnt/β-catenin signaling. Notably, genotype–phenotype correlations have been established at the APC gene between the length and stability of the truncated proteins encoded by different mutant alleles, the corresponding levels of Wnt/β-catenin signaling activity they encode for, and the incidence and distribution of intestinal and extra-intestinal tumors. Here, we report a novel mouse model, Apc1572T, obtained by targeting a truncated mutation at codon 1572 in the endogenous Apc gene. This hypomorphic mutant allele results in intermediate levels of Wnt/β-catenin signaling activation when compared with other Apc mutations associated with multifocal intestinal tumors. Notwithstanding the constitutive nature of the mutation, Apc+/1572T mice have no predisposition to intestinal cancer but develop multifocal mammary adenocarcinomas and subsequent pulmonary metastases in both genders. The histology of the Apc1572T primary mammary tumours is highly heterogeneous with luminal, myoepithelial, and squamous lineages and is reminiscent of metaplastic carcinoma of the breast in humans. The striking phenotype of Apc+/1572T mice suggests that specific dosages of Wnt/β-catenin signaling activity differentially affect tissue homeostasis and initiate tumorigenesis in an organ-specific fashion

Genetic parameters for social effects on survival in cannibalistic layers: Combining survival analysis and a linear animal model

<p>Abstract</p> <p>Background</p> <p>Mortality due to cannibalism in laying hens is a difficult trait to improve genetically, because censoring is high (animals still alive at the end of the testing period) and it may depend on both the individual itself and the behaviour of its group members, so-called associative effects (social interactions). To analyse survival data, survival analysis can be used. However, it is not possible to include associative effects in the current software for survival analysis. A solution could be to combine survival analysis and a linear animal model including associative effects. This paper presents a two-step approach (2STEP), combining survival analysis and a linear animal model including associative effects (LAM).</p> <p>Methods</p> <p>Data of three purebred White Leghorn layer lines from Institut de Sélection Animale B.V., a Hendrix Genetics company, were used in this study. For the statistical analysis, survival data on 16,780 hens kept in four-bird cages with intact beaks were used. Genetic parameters for direct and associative effects on survival time were estimated using 2STEP. Cross validation was used to compare 2STEP with LAM. LAM was applied directly to estimate genetic parameters for social effects on observed survival days.</p> <p>Results</p> <p>Using 2STEP, total heritable variance, including both direct and associative genetic effects, expressed as the proportion of phenotypic variance, ranged from 32% to 64%. These results were substantially larger than when using LAM. However, cross validation showed that 2STEP gave approximately the same survival curves and rank correlations as LAM. Furthermore, cross validation showed that selection based on both direct and associative genetic effects, using either 2STEP or LAM, gave the best prediction of survival time.</p> <p>Conclusion</p> <p>It can be concluded that 2STEP can be used to estimate genetic parameters for direct and associative effects on survival time in laying hens. Using 2STEP increased the heritable variance in survival time. Cross validation showed that social genetic effects contribute to a large difference in survival days between two extreme groups. Genetic selection targeting both direct and associative effects is expected to reduce mortality due to cannibalism in laying hens.</p

Stable S/MAR-based episomal vectors are regulated at the chromatin level

Episomal vectors assembled from defined genetic components are a promising alternative to traditional gene therapy vectors that integrate in the host genome and may cause insertional mutations. The vector pEPI-eGFP is stably retained in the episomal state in cultured mammalian cells at low copy number for many generations without integration into the host genome. Although pEPI-eGFP is a fully engineered vector, little is known about how it interacts with the host genome and about the molecular mechanisms that are responsible for its transcriptional activity. We have analyzed the expression of the episomal reporter gene eGFP under conditions that affect the chromatin state of the genome. We have also constructed pEPI derivatives carrying a tandem array of lac operator sequences, which allows in vivo visualization and manipulation of the chromatin state of the episome. We show that changes in chromatin state of both the host and pEPI-eGFP induces changes in episomal gene activity and influences the episome’s nuclear distributions. We conclude that episomal genes are subject to control systems of the host, similarly to their counterparts in the host genome

Tree height integrated into pantropical forest biomass estimates

Copyright © 2012 European Geosciences Union. This is the published version available at http://www.biogeosciences.net/9/3381/2012/bg-9-3381-2012.htmlAboveground tropical tree biomass and carbon storage estimates commonly ignore tree height (H). We estimate the effect of incorporating H on tropics-wide forest biomass estimates in 327 plots across four continents using 42 656 H and diameter measurements and harvested trees from 20 sites to answer the following questions: 1. What is the best H-model form and geographic unit to include in biomass models to minimise site-level uncertainty in estimates of destructive biomass? 2. To what extent does including H estimates derived in (1) reduce uncertainty in biomass estimates across all 327 plots? 3. What effect does accounting for H have on plot- and continental-scale forest biomass estimates? The mean relative error in biomass estimates of destructively harvested trees when including H (mean 0.06), was half that when excluding H (mean 0.13). Power- and Weibull-H models provided the greatest reduction in uncertainty, with regional Weibull-H models preferred because they reduce uncertainty in smaller-diameter classes (≤40 cm D) that store about one-third of biomass per hectare in most forests. Propagating the relationships from destructively harvested tree biomass to each of the 327 plots from across the tropics shows that including H reduces errors from 41.8 Mg ha−1 (range 6.6 to 112.4) to 8.0 Mg ha−1 (−2.5 to 23.0). For all plots, aboveground live biomass was −52.2 Mg ha−1 (−82.0 to −20.3 bootstrapped 95% CI), or 13%, lower when including H estimates, with the greatest relative reductions in estimated biomass in forests of the Brazilian Shield, east Africa, and Australia, and relatively little change in the Guiana Shield, central Africa and southeast Asia. Appreciably different stand structure was observed among regions across the tropical continents, with some storing significantly more biomass in small diameter stems, which affects selection of the best height models to reduce uncertainty and biomass reductions due to H. After accounting for variation in H, total biomass per hectare is greatest in Australia, the Guiana Shield, Asia, central and east Africa, and lowest in east-central Amazonia, W. Africa, W. Amazonia, and the Brazilian Shield (descending order). Thus, if tropical forests span 1668 million km2 and store 285 Pg C (estimate including H), then applying our regional relationships implies that carbon storage is overestimated by 35 Pg C (31–39 bootstrapped 95% CI) if H is ignored, assuming that the sampled plots are an unbiased statistical representation of all tropical forest in terms of biomass and height factors. Our results show that tree H is an important allometric factor that needs to be included in future forest biomass estimates to reduce error in estimates of tropical carbon stocks and emissions due to deforestation

Meta-analysis of genome-wide association studies for cattle stature identifies common genes that regulate body size in mammals

peer-reviewedH.D.D., A.J.C., P.J.B. and B.J.H. would like to acknowledge the Dairy Futures Cooperative Research Centre for funding. H.P. and R.F. acknowledge funding from the German Federal Ministry of Education and Research (BMBF) within the AgroClustEr ‘Synbreed—Synergistic Plant and Animal Breeding’ (grant 0315527B). H.P., R.F., R.E. and K.-U.G. acknowledge the Arbeitsgemeinschaft Süddeutscher Rinderzüchter, the Arbeitsgemeinschaft Österreichischer Fleckviehzüchter and ZuchtData EDV Dienstleistungen for providing genotype data. A. Bagnato acknowledges the European Union (EU) Collaborative Project LowInputBreeds (grant agreement 222623) for providing Brown Swiss genotypes. Braunvieh Schweiz is acknowledged for providing Brown Swiss phenotypes. H.P. and R.F. acknowledge the German Holstein Association (DHV) and the Confederación de Asociaciones de Frisona Española (CONCAFE) for sharing genotype data. H.P. was financially supported by a postdoctoral fellowship from the Deutsche Forschungsgemeinschaft (DFG) (grant PA 2789/1-1). D.B. and D.C.P. acknowledge funding from the Research Stimulus Fund (11/S/112) and Science Foundation Ireland (14/IA/2576). M.S. and F.S.S. acknowledge the Canadian Dairy Network (CDN) for providing the Holstein genotypes. P.S. acknowledges funding from the Genome Canada project entitled ‘Whole Genome Selection through Genome Wide Imputation in Beef Cattle’ and acknowledges WestGrid and Compute/Calcul Canada for providing computing resources. J.F.T. was supported by the National Institute of Food and Agriculture, US Department of Agriculture, under awards 2013-68004-20364 and 2015-67015-23183. A. Bagnato, F.P., M.D. and J.W. acknowledge EU Collaborative Project Quantomics (grant 516 agreement 222664) for providing Brown Swiss and Finnish Ayrshire sequences and genotypes. A.C.B. and R.F.V. acknowledge funding from the public–private partnership ‘Breed4Food’ (code BO-22.04-011- 001-ASG-LR) and EU FP7 IRSES SEQSEL (grant 317697). A.C.B. and R.F.V. acknowledge CRV (Arnhem, the Netherlands) for providing data on Dutch and New Zealand Holstein and Jersey bulls.Stature is affected by many polymorphisms of small effect in humans1. In contrast, variation in dogs, even within breeds, has been suggested to be largely due to variants in a small number of genes2,3. Here we use data from cattle to compare the genetic architecture of stature to those in humans and dogs. We conducted a meta-analysis for stature using 58,265 cattle from 17 populations with 25.4 million imputed whole-genome sequence variants. Results showed that the genetic architecture of stature in cattle is similar to that in humans, as the lead variants in 163 significantly associated genomic regions (P < 5 × 10−8) explained at most 13.8% of the phenotypic variance. Most of these variants were noncoding, including variants that were also expression quantitative trait loci (eQTLs) and in ChIP–seq peaks. There was significant overlap in loci for stature with humans and dogs, suggesting that a set of common genes regulates body size in mammals

Genome-wide association analysis implicates dysregulation of immunity genes in chronic lymphocytic leukaemia

Several chronic lymphocytic leukaemia (CLL) susceptibility loci have been reported; however, much of the heritable risk remains unidentified. Here we perform a meta-analysis of six genome-wide association studies, imputed using a merged reference panel of 1,000 Genomes and UK10K data, totalling 6,200 cases and 17,598 controls after replication. We identify nine risk loci at 1p36.11 (rs34676223, P=5.04 × 10−13), 1q42.13 (rs41271473, P=1.06 × 10−10), 4q24 (rs71597109, P=1.37 × 10−10), 4q35.1 (rs57214277, P=3.69 × 10−8), 6p21.31 (rs3800461, P=1.97 × 10−8), 11q23.2 (rs61904987, P=2.64 × 10−11), 18q21.1 (rs1036935, P=3.27 × 10−8), 19p13.3 (rs7254272, P=4.67 × 10−8) and 22q13.33 (rs140522, P=2.70 × 10−9). These new and established risk loci map to areas of active chromatin and show an over-representation of transcription factor binding for the key determinants of B-cell development and immune response