Transferring QTL technology to the pig breeding industry (PigQTech) - a demonstration project

The aim of the project has been to demonstrate how the farm animal breeding industry can utilise gene mapping technology to accelerate genetic improvement. Previous theoretical studies had suggested that the use of marker assisted selection could potentially increase the annual improvement for quantitative traits like backfat with about 10% and for more difficult traits such as meat quality and reproduction by as much as 40-60% compared with existing technology. The work has comprised two major tasks: 1. Commercially relevant populations have been screened for segregation at QTLs identified in experimental populations. The aim has been to establish optimal strategies for QTL detection in commercial pig populations and the extent to which QTLs explaining major phenotypic differences between divergent lines used in experimental studies also explain quantitative variation within commercial lines. The results are important for specifying future strategies for finding economically valuable QTLs. 2. Marker assisted backcrossing has been used to demonstrate how a QTL allele can be introgressed from one breed to another. The work has focused on the major fatness QTL on pig chromosome 4 previously identified in a wild pig/Large White intercross. The end result was not designed to be a commercially viable product in its own right, but the process has validated a number of points of major importance for the exploitation of QTLs in livestock.

Alternative ways to estimate change points in multinomial sequences. An application to an authorship attribution problem

The statistical analysis of literary style is the part of stylometry that compares measurable characteristicsin a text that are rarely controlled by the author, with those in other texts. When thegoal is to settle authorship questions, these characteristics should relate to the author’s style andnot to the genre, epoch or editor, and they should be such that their variation between authors islarger than the variation within comparable texts from the same author.For an overview of the literature on stylometry and some of the techniques involved, see for exampleMosteller and Wallace (1964, 82), Herdan (1964), Morton (1978), Holmes (1985), Oakes (1998) orLebart, Salem and Berry (1998).Tirant lo Blanc, a chivalry book, is the main work in catalan literature and it was hailed to be“the best book of its kind in the world” by Cervantes in Don Quixote. Considered by writterslike Vargas Llosa or Damaso Alonso to be the first modern novel in Europe, it has been translatedseveral times into Spanish, Italian and French, with modern English translations by Rosenthal(1996) and La Fontaine (1993). The main body of this book was written between 1460 and 1465,but it was not printed until 1490.There is an intense and long lasting debate around its authorship sprouting from its first edition,where its introduction states that the whole book is the work of Martorell (1413?-1468), while atthe end it is stated that the last one fourth of the book is by Galba (?-1490), after the death ofMartorell. Some of the authors that support the theory of single authorship are Riquer (1990),Chiner (1993) and Badia (1993), while some of those supporting the double authorship are Riquer(1947), Coromines (1956) and Ferrando (1995). For an overview of this debate, see Riquer (1990).Neither of the two candidate authors left any text comparable to the one under study, and thereforediscriminant analysis can not be used to help classify chapters by author. By using sample textsencompassing about ten percent of the book, and looking at word length and at the use of 44conjunctions, prepositions and articles, Ginebra and Cabos (1998) detect heterogeneities that mightindicate the existence of two authors. By analyzing the diversity of the vocabulary, Riba andGinebra (2000) estimates that stylistic boundary to be near chapter 383.Following the lead of the extensive literature, this paper looks into word length, the use of the mostfrequent words and into the use of vowels in each chapter of the book. Given that the featuresselected are categorical, that leads to three contingency tables of ordered rows and therefore tothree sequences of multinomial observations.Section 2 explores these sequences graphically, observing a clear shift in their distribution. Section 3describes the problem of the estimation of a suden change-point in those sequences, in the followingsections we propose various ways to estimate change-points in multinomial sequences; the methodin section 4 involves fitting models for polytomous data, the one in Section 5 fits gamma modelsonto the sequence of Chi-square distances between each row profiles and the average profile, theone in Section 6 fits models onto the sequence of values taken by the first component of thecorrespondence analysis as well as onto sequences of other summary measures like the averageword length. In Section 7 we fit models onto the marginal binomial sequences to identify thefeatures that distinguish the chapters before and after that boundary. Most methods rely heavilyon the use of generalized linear models

A multilevel analysis on the determinants of regional health care expenditure. A note

We apply a multilevel hierarchical model to explore whether anaggregation fallacy exists in estimating the income elasticity of healthexpenditure by ignoring the regional composition of national healthexpenditure figures. We use data for 110 regions in eight OECD countriesin 1997: Australia, Canada, France, Germany, Italy, Spain, Sweden andUnited Kingdom. In doing this we have tried to identify two sources ofrandom variation: within countries and between-countries. Our resultsshow that: 1- Variability between countries amounts to (SD) 0.5433, andjust 13% of that can be attributed to income elasticity and the remaining87% to autonomous health expenditure; 2- Within countries, variabilityamounts to (SD) 1.0249; and 3- The intra-class correlation is 0.5300. Weconclude that we have to take into account the degree of fiscaldecentralisation within countries in estimating income elasticity ofhealth expenditure. Two reasons lie behind this: a) where there isdecentralisation to the regions, policies aimed at emulating diversitytend to increase national health care expenditure; and b) without fiscaldecentralisation, central monitoring of finance tends to reduce regionaldiversity and therefore decrease national health expenditure.

Oscillatory magnetic tweezers based on ferromagnetic beads and simple coaxial coils

We report the design and validation of simple magnetic tweezers for oscillating ferromagnetic beads in the piconewton and nanometer scales. The system is based on a single pair of coaxial coils operating in two sequential modes: permanent magnetization of the beads through a large and brief pulse of magnetic field and generation of magnetic gradients to produce uniaxial oscillatory forces. By using this two step method, the magnetic moment of the beads remains constant during measurements. Therefore, the applied force can be computed and varies linearly with the driving signal. No feedback control is required to produce well defined force oscillations over a wide bandwidth. The design of the coils was optimized to obtain high magnetic fields (280 mT) and gradients (2 T/m) with high homogeneity (5% variation) within the sample. The magnetic tweezers were implemented in an inverted optical microscope with a videomicroscopy-based multiparticle tracking system. The apparatus was validated with 4.5 ¿m magnetite beads obtaining forces up to ~2 pN and subnanometer resolution. The applicability of the device includes microrheology of biopolymer and cell cytoplasm, molecular mechanics, and mechanotransduction in living cells.

Dugesia sicula (Platyhelminthes, Tricladida): the colonizing success of an asexual Planarian

Background Dugesia sicula is the only species of its genus not presenting an endemic or restricted distribution within the Mediterranean area. It mostly comprises fissiparous populations (asexual reproduction by body division and regeneration), most likely sexually sterile, and characterized by an extremely low genetic diversity interpreted as the consequence of a recent anthropic expansion. However, its fissiparous reproduction can result in an apparent lack of diversity within the species, since genetic variation within individuals can be as large as between them because most individuals within a population are clones. We have estimated haplotype and nucleotide diversity of cytochrome oxidase I within and among individuals along the species distribution of a broad sample of D. sicula, including asexual and the two only sexual populations known today; and predicted its potential distribution based on climatic variables. Our aim was to determine the centre of colonisation origin, whether the populations are recent, and whether the species is expanding. Results The species presents 3 most frequent haplotypes, differing in a maximum of 11 base pairs. As expected from their fissiparous mode of reproduction, in half of all the analysed localities many individuals have multiple heteroplasmic haplotypes. The distribution of haplotypes is not geographically structured; however, the distribution of haplotypes and heteroplasmic populations shows higher diversity in the central Mediterranean region. The potential distribution predicted by climatic variables based modelling shows a preference for coastal areas and fits well with the observed data. Conclusions The distribution and frequency of the most frequent haplotypes and the presence of heteroplasmic individuals allow us to gain an understanding of the recent history of the species, together with previous knowledge on its phylogenetic relationships and age: The species most probably originated in Africa and dispersed through the central Mediterranean. After one or multiple populations became triploid and fissiparous, the species colonized the Mediterranean basin, likely both by its own means and helped by human activities. Its present distribution practically fulfils its potential distribution as modelled with climatic variables. Its prevalence in coastal regions with higher water temperatures predicts a likely future expansion to northern and more interior areas following the increase in temperatures due to climate change.

A functional variant in the Stearoyl-CoA desaturase gene promoter enhances fatty acid desaturation in pork

There is growing public concern about reducing saturated fat intake. Stearoyl-CoA desaturase (SCD) is the lipogenic enzyme responsible for the biosynthesis of oleic acid (18:1) by desaturating stearic acid (18:0). Here we describe a total of 18 mutations in the promoter and 3′ non-coding region of the pig SCD gene and provide evidence that allele T at AY487830:g.2228T>C in the promoter region enhances fat desaturation (the ratio 18:1/18:0 in muscle increases from 3.78 to 4.43 in opposite homozygotes) without affecting fat content (18:0+18:1, intramuscular fat content, and backfat thickness). No mutations that could affect the functionality of the protein were found in the coding region. First, we proved in a purebred Duroc line that the C-T-A haplotype of the 3 single nucleotide polymorphisms (SNPs) (g.2108C>T; g.2228T>C; g.2281A>G) of the promoter region was additively associated to enhanced 18:1/18:0 both in muscle and subcutaneous fat, but not in liver. We show that this association was consistent over a 10-year period of overlapping generations and, in line with these results, that the C-T-A haplotype displayed greater SCD mRNA expression in muscle. The effect of this haplotype was validated both internally, by comparing opposite homozygote siblings, and externally, by using experimental Duroc-based crossbreds. Second, the g.2281A>G and the g.2108C>T SNPs were excluded as causative mutations using new and previously published data, restricting the causality to g.2228T>C SNP, the last source of genetic variation within the haplotype. This mutation is positioned in the core sequence of several putative transcription factor binding sites, so that there are several plausible mechanisms by which allele T enhances 18:1/18:0 and, consequently, the proportion of monounsaturated to saturated fat.

Multi-scale effects of nestling diet on breeding performance in a terrestrial top predator inferred from stable isotope analysis

Inter-individual diet variation within populations is likely to have important ecological and evolutionary implications. The diet-fitness relationships at the individual level and the emerging population processes are, however, poorly understood for most avian predators inhabiting complex terrestrial ecosystems. In this study, we use an isotopic approach to assess the trophic ecology of nestlings in a long-lived raptor, the Bonelli"s eagle Aquila fasciata, and investigate whether nestling dietary breath and main prey consumption can affect the species" reproductive performance at two spatial scales: territories within populations and populations over a large geographic area. At the territory level, those breeding pairs whose nestlings consumed similar diets to the overall population (i.e. moderate consumption of preferred prey, but complemented by alternative prey categories) or those disproportionally consuming preferred prey were more likely to fledge two chicks. An increase in the diet diversity, however, related negatively with productivity. The age and replacements of breeding pair members had also an influence on productivity, with more fledglings associated to adult pairs with few replacements, as expected in long-lived species. At the population level, mean productivity was higher in those population-years with lower dietary breadth and higher diet similarity among territories, which was related to an overall higher consumption of preferred prey. Thus, we revealed a correspondence in diet-fitness relationships at two spatial scales: territories and populations. We suggest that stable isotope analyses may be a powerful tool to monitor the diet of terrestrial avian predators on large spatio-temporal scales, which could serve to detect potential changes in the availability of those prey on which predators depend for breeding. We encourage ecologists and evolutionary and conservation biologists concerned with the multi-scale fitness consequences of inter-individual variation in resource use to employ similar stable isotope-based approaches, which can be successfully applied to complex ecosystems such as the Mediterranean.

DnaSP Version 5. DNA Sequence Polymorphism

DnaSP, DNA Sequence Polymorphism, is a software package for the analysis of nucleotide polymorphism from aligned DNA sequence data. DnaSP can estimate several measures of DNA sequence variation within and between populations (in noncoding, synonymous or nonsynonymous sites, or in various sorts of codon positions), as well as linkage disequilibrium, recombination, gene flow and gene conversion parameters. DnaSP can also carry out several tests of neutrality: Hudson, Kreitman and Aguadé (1987), Tajima (1989), McDonald and Kreitman (1991), Fu and Li (1993), and Fu (1997) tests. Additionally, DnaSP can estimate the confidence intervals of some test-statistics by the coalescent. The results of the analyses are displayed on tabular and graphic form.

Economic crisis and recovery: changes in second birth rates within occupational classes and educational groups

This study assesses the decline in second birth rates for men and women across different skill levels in transitional Russia. Changes within educational groups and occupational classes are observed over three distinct time periods: the Soviet era, economic crisis, and economic recovery. The most remarkable finding is the similarity in the extent second birth rates declined within educational groups and occupational classes during the economic crisis. Although further decline occurred in the recovery period, more variation emerged across groups.

Seasonal patterns op biomass variation of Ruppia cirrhosa (Petagna) Grande and Potamageton pectinatus L. in a coastal logoon

Coastal lagoons where salinity varies within a wide range during the year are colonized by euryhaline macrophytes which can develop extensive beds. Seasonal changes in biomass of Ruppia cirrhosa and Potamogeton pectinatus were studied in Tancada Lagoon (Ebro Delta, NE Spain) in order to reveal the environmental factors controlling their population development. Ruppia cirrhosa occupy a larger area of the lagoon than Potarnogeton pectinatus. Their maximum above ground biomasses are also different (495 g m-2 and 351 g m-2 ash free dry weight, respectively). Below ground biomass of Ruppia cirrhosa is between 9 and 53 % of the above ground biomass, while it is 3-40 % for Potamogeton pectinatus. Chlorophyll a contents show fluctuations similar to biomass. Low salinity and high turbidity caused by freshwater inflows favour Potamogeton expansion, while Ruppia development is favoured by high salinity and transparent water.

Color differences among feral pigeons (Columba livia) are not attributable to sequence variation in the coding region of the melanocortin-1 receptor gene (MC1R)

Genetic variation at the melanocortin-1 receptor (MC1R) gene is correlated with melanin color variation in many birds. Feral pigeons (Columba livia) show two major melanin-based colorations: a red coloration due to pheomelanic pigment and a black coloration due to eumelanic pigment. Furthermore, within each color type, feral pigeons display continuous variation in the amount of melanin pigment present in the feathers, with individuals varying from pure white to a full dark melanic color. Coloration is highly heritable and it has been suggested that it is under natural or sexual selection, or both. Our objective was to investigate whether MC1R allelic variants are associated with plumage color in feral pigeons.We sequenced 888 bp of the coding sequence of MC1R among pigeons varying both in the type, eumelanin or pheomelanin, and the amount of melanin in their feathers. We detected 10 non-synonymous substitutions and 2 synonymous substitution but none of them were associated with a plumage type. It remains possible that non-synonymous substitutions that influence coloration are present in the short MC1R fragment that we did not sequence but this seems unlikely because we analyzed the entire functionally important region of the gene.Our results show that color differences among feral pigeons are probably not attributable to amino acid variation at the MC1R locus. Therefore, variation in regulatory regions of MC1R or variation in other genes may be responsible for the color polymorphism of feral pigeons.

Gorilla genome structural variation reveals evolutionary parallelisms with chimpanzee

Structural variation has played an important role in the evolutionary restructuring of human and great ape genomes. Recent analyses have suggested that the genomes of chimpanzee and human have been particularly enriched for this form of genetic variation. Here, we set out to assess the extent of structural variation in the gorilla lineage by generating 10-fold genomic sequence coverage from a western lowland gorilla and integrating these data into a physical and cytogenetic framework of structural variation. We discovered and validated over 7665 structural changes within the gorilla lineage, including sequence resolution of inversions, deletions, duplications, and mobile element insertions. A comparison with human and other ape genomes shows that the gorilla genome has been subjected to the highest rate of segmental duplication. We show that both the gorilla and chimpanzee genomes have experienced independent yet convergent patterns of structural mutation that have not occurred in humans, including the formation of subtelomeric heterochromatic caps, the hyperexpansion of segmental duplications, and bursts of retroviral integrations. Our analysis suggests that the chimpanzee and gorilla genomes are structurally more derived than either orangutan or human genomes.

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.

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.

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