Differential expression analysis for sequence count data via mixtures of negative binomials

The recent advent of Next-generation sequencing technologies has revolutionized the way of analyzing the genome. This innovation allows to get deeper information at a lower cost and in less time, and provides data that are discrete measurements. One of the most important applications with these data is the differential analysis, that is investigating if one gene exhibit a different expression level in correspondence of two (or more) biological conditions (such as disease states, treatments received and so on). As for the statistical analysis, the final aim will be statistical testing and for modeling these data the Negative Binomial distribution is considered the most adequate one especially because it allows for "over dispersion". However, the estimation of the dispersion parameter is a very delicate issue because few information are usually available for estimating it. Many strategies have been proposed, but they often result in procedures based on plug-in estimates, and in this thesis we show that this discrepancy between the estimation and the testing framework can lead to uncontrolled first-type errors. We propose a mixture model that allows each gene to share information with other genes that exhibit similar variability. Afterwards, three consistent statistical tests are developed for differential expression analysis. We show that the proposed method improves the sensitivity of detecting differentially expressed genes with respect to the common procedures, since it is the best one in reaching the nominal value for the first-type error, while keeping elevate power. The method is finally illustrated on prostate cancer RNA-seq data.

High variability and non-neutral evolution of the mammalian avpr1a gene

BACKGROUND: The arginine-vasopressin 1a receptor has been identified as a key determinant for social behaviour in Microtus voles, humans and other mammals. Nevertheless, the genetic bases of complex phenotypic traits like differences in social and mating behaviour among species and individuals remain largely unknown. Contrary to previous studies focusing on differences in the promotor region of the gene, we investigate here the level of functional variation in the coding region (exon 1) of this locus. RESULTS: We detected high sequence diversity between higher mammalian taxa as well as between species of the genus Microtus. This includes length variation and radical amino acid changes, as well as the presence of distinct protein variants within individuals. Additionally, negative selection prevails on most parts of the first exon of the arginine-vasopressin receptor 1a (avpr1a) gene but it contains regions with higher rates of change that harbour positively selected sites. Synonymous and non-synonymous substitution rates in the avpr1a gene are not exceptional compared to other genes, but they exceed those found in related hormone receptors with similar functions. DISCUSSION: These results stress the importance of considering variation in the coding sequence of avpr1a in regards to associations with life history traits (e.g. social behaviour, mating system, habitat requirements) of voles, other mammals and humans in particular.

1200 years of decadal-scale variability of Mediterranean vegetation and climate at Pantelleria Island, Italy

A new sedimentary sequence from Lago di Venere on Pantelleria Island, located in the Strait of Sicily between Tunisia and Sicily was recovered. The lake is located in the coastal infra-Mediterranean vegetation belt at 2 m a.s.l. Pollen, charcoal and sedimentological analyses are used to explore linkages among vegetation, fire and climate at a decadal scale over the past 1200 years. A dry period from ad 800 to 1000 that corresponds to the ‘Medieval Warm Period’ (WMP) is inferred from sedimentological analysis. The high content of carbonate recorded in this period suggests a dry phase, when the ratio of evaporation/precipitation was high. During this period the island was dominated by thermophilous and drought-tolerant taxa, such as Quercus ilex, Olea, Pistacia and Juniperus. A marked shift in the sediment properties is recorded at ad 1000, when carbonate content became very low suggesting wetter conditions until ad 1850–1900. Broadly, this period coincides with the ‘Little Ice Age’ (LIA), which was characterized by wetter and colder conditions in Europe. During this time rather mesic conifers (i.e. Pinus pinaster), shrubs and herbs (e.g. Erica arborea and Selaginella denticulata) expanded, whereas more drought-adapted species (e.g. Q. ilex) declined. Charcoal data suggest enhanced fire activity during the LIA probably as a consequence of anthropogenic burning and/or more flammable fuel (e.g. resinous Pinus biomass). The last century was characterized by a shift to high carbonate content, indicating a change towards drier conditions, and re-expansion of Q. ilex and Olea. The post-LIA warming is in agreement with historical documents and meteorological time series. Vegetation dynamics were co-determined by agricultural activities on the island. Anthropogenic indicators (e.g. Cerealia-type, Sporormiella) reveal the importance of crops and grazing on the island. Our pollen data suggest that extensive logging caused the local extinction of deciduous Quercus pubescens around ad1750.

Sedimentological processes and environmental variability at Lake Ohrid (Macedonia, Albania) between 637 ka and the present

Lake Ohrid (Macedonia, Albania) is thought to be more than 1.2 million years old and host more than 300 endemic species. As a target of the International Continental scientific Drilling Program (ICDP), a successful deep drilling campaign was carried out within the scope of the Scientific Collaboration on Past Speciation Conditions in Lake Ohrid (SCOPSCO) project in 2013. Here, we present lithological, sedimentological, and (bio-)geochemical data from the upper 247.8 m composite depth of the overall 569 m long DEEP site sediment succession from the central part of the lake. According to an age model, which is based on 11 tephra layers (first-order tie points) and on tuning of bio-geochemical proxy data to orbital parameters (second-order tie points), the analyzed sediment sequence covers the last 637 kyr. The DEEP site sediment succession consists of hemipelagic sediments, which are interspersed by several tephra layers and infrequent, thin (< 5 cm) mass wasting deposits. The hemipelagic sediments can be classified into three different lithotypes. Lithotype 1 and 2 deposits comprise calcareous and slightly calcareous silty clay and are predominantly attributed to interglacial periods with high primary productivity in the lake during summer and reduced mixing during winter. The data suggest that high ion and nutrient concentrations in the lake water promoted calcite precipitation and diatom growth in the epilimnion during MIS15, 13, and 5. Following a strong primary productivity, highest interglacial temperatures can be reported for marine isotope stages (MIS) 11 and 5, whereas MIS15, 13, 9, and 7 were comparably cooler. Lithotype 3 deposits consist of clastic, silty clayey material and predominantly represent glacial periods with low primary productivity during summer and longer and intensified mixing during winter. The data imply that the most severe glacial conditions at Lake Ohrid persisted during MIS16, 12, 10, and 6, whereas somewhat warmer temperatures can be inferred for MIS14, 8, 4, and 2. Interglacial-like conditions occurred during parts of MIS14 and 8.

Distribution and variability of deformed wing virus of honeybees (Apis mellifera) in the Middle East and North Africa.

Three hundred eleven honeybee samples from twelve countries in the Middle East and North Africa (MENA) (Jordan, Lebanon, Syria, Iraq, Egypt, Libya, Tunisia, Algeria, Morocco, Yemen, Palestine and Sudan) were analyzed for the presence of deformed wing virus (DWV). The prevalence of DWV throughout the MENA region was pervasive, but variable. The highest prevalence was found in Lebanon and Syria, with prevalence dropping in Palestine, Jordan and Egypt before increasing slightly moving westwards to Algeria and Morocco Phylogenetic analysis of a 194 nucleotide section of the DWV Lp gene did not identify any significant phylogenetic resolution among the samples, although the sequences did show consistent regional clustering, including an interesting geographic gradient from Morocco through North Africa to Jordan and Syria. The sequences revealed several clear variability hotspots in the deduced amino acid sequence, that furthermore showed some patterns of regional identity. Furthermore, the sequence variants from the Middle East and North Africa appear more numerous and diverse than those from Europe. This article is protected by copyright. All rights reserved.

Middle to late Miocene oxygen isotope stratigraphy of ODP site 1085 (SE Atlantic): new constrains on Miocene climate variability and sea-level fluctuations

The middle Miocene delta18O increase represents a fundamental change in earth's climate system due to a major expansion and permanent establishment of the East Antarctic Ice Sheet accompanied by some effect of deepwater cooling. The long-term cooling trend in the middle to late Miocene was superimposed by several punctuated periods of glaciations (Mi-Events) characterized by oxygen isotopic shifts that have been related to the waxing and waning of the Antarctic ice-sheet and bottom water cooling. Here, we present a high-resolution benthic stable oxygen isotope record from ODP Site 1085 located at the southwestern African continental margin that provides a detailed chronology for the middle to late Miocene (13.9-7.3 Ma) climate transition in the eastern South Atlantic. A composite Fe intensity record obtained by XRF core scanning ODP Sites 1085 and 1087 was used to construct an astronomically calibrated chronology based on orbital tuning. The oxygen isotope data exhibit four distinct delta18O excursions, which have astronomical ages of 13.8, 13.2, 11.7, and 10.4 Ma and correspond to the Mi3, Mi4, Mi5, and Mi6 events. A global climate record was extracted from the oxygen isotopic composition. Both long- and short-term variabilities in the climate record are discussed in terms of sea-level and deep-water temperature changes. The oxygen isotope data support a causal link between sequence boundaries traced from the shelf and glacioeustatic changes due to ice-sheet growth. Spectral analysis of the benthic delta18O record shows strong power in the 400-kyr and 100-kyr bands documenting a paleoceanographic response to eccentricity-modulated variations in precession. A spectral peak around 180-kyr might be related to the asymmetry of the obliquity cycle indicating that the response of the dominantly unipolar Antarctic ice-sheet to obliquityinduced variations probably controlled the middle to late Miocene climate system. Maxima in the delta18O record, interpreted as glacial periods, correspond to minima in 100-kyr eccentricity cycle and minima in the 174-kyr obliquity modulation. Strong middle to late Miocene glacial events are associated with 400-kyr eccentricity minima and obliquity modulation minima. Thus, fluctuations in the amplitude of obliquity and eccentricity seem to be the driving force for the middle to late Miocene climate variability.

Genetic variability and limited clonality of Mycoplasma hyorhinis in pig herds

Mycoplasma hyorhinis is a common inhabitant of the upper respiratory tract and tonsils of pigs. Its role as a possible pathogen remains controversial. In order to gain more insight into the epidemiology and population structure of M. hyorhinis we genetically characterized 60 isolates by multi locus sequence typing (MLST). The M. hyorhinis strains originated from Swiss and German pig herds with knowledge on the clinical background. The MLST scheme of Tocqueville et al. (J. Clin. Microbiol. 2014) was optimized, primers for the six MLST gene fragments were newly designed to allow amplification and sequencing with a single protocol. A total of 27 ST were observed with the 60 strains, 26 of those were previously unknown types. Generally identical genotypes were observed within a farm but they differed between farms. The identical genotype was also observed in three different Swiss farms. On the other Hand different genotypes within a farm were found with three German farms. The Swiss isolates formed a distinct cluster but otherwise there was no geographical nor a clinical association with specific Clusters observed. Data shows a high variability of M. hyorhinis comparable to what is observed for Mycoplasma hyopneumoniae. Similar to this pathogen the population structure of M. hyorhinis also shows some limited clonality with predominant genotypes within an animal and a single farm but different ones between farms. The comparable population structure of M. hyopneumoniae and M. hyorhinis could indicate a similar evolution of the two species in the common pig host.

Genetic variability of Tomato spotted wilt virus in Australia and validation of real time RT-PCR for its detection in single and bulked leaf samples

The potential for large-scale use of a sensitive real time reverse transcription polymerase chain reaction (RT-PCR) assay was evaluated for the detection of Tomato spotted wilt virus (TSWV) in single and bulked leaf samples by comparing its sensitivity with that of DAS-ELISA. Using total RNA extracted with RNeasy (R) or leaf soak methods, real time RT-PCR detected TSWV in all infected samples collected from 16 horticultural crop species (including flowers, herbs and vegetables), two arable crop species, and four weed species by both assays. In samples in which DAS-ELISA had previously detected TSWV, real time RT-PCR was effective at detecting it in leaf tissues of all 22 plant species tested at a wide range of concentrations. Bulk samples required more robust and extensive extraction methods with real time RT-PCR, but it generally detected one infected sample in 1000 uninfected ones. By contrast, ELISA was less sensitive when used to test bulked samples, once detecting up to I infected in 800 samples with pepper but never detecting more than I infected in 200 samples in tomato and lettuce. It was also less reliable than real time RT-PCR when used to test samples from parts of the leaf where the virus concentration was low. The genetic variability among Australian isolates of TSWV was small. Direct sequencing of a 587 bp region of the nucleoprotein gene (S RNA) of 29 isolates from diverse crops and geographical locations yielded a maximum of only 4.3% nucleotide sequence difference. Phylogenetic analysis revealed no obvious groupings of isolates according to geographic origin or host species. TSWV isolates, that break TSWV resistance genes in tomato or pepper did not differ significantly in the N gene region studied, indicating that a different region of the virus genome is responsible for this trait.

Nucleocapsid gene variability reveals two subgroups of Lettuce necrotic yellows virus

The complete nucleocapsid (N) genes of eight Australian isolates of Lettuce necrotic yellows virus (LNYV) were amplified by reverse transcription PCR, cloned and sequenced. Phylogenetic analyses of these sequences revealed two distinct subgroups of LNYV isolates. Nucleotide sequences within each subgroup were more than 96% identical but heterogeneity between groups was about 20% at the nucleotide sequence level. However, less than 4% heterogeneity was noted at the amino acid level, indicating mostly third nucleotide position changes and a strong conservation for N protein function. There was no obvious geographical or temporal separation of the subgroups in Australia.

Completion of the genome sequence of Lettuce necrotic yellows virus, type species of the genus Cytorhabdovirus

We completed the genome sequence of Lettuce necrotic yellows virus (LNYV) by determining the nucleotide sequences of the 4a (putative phosphoprotein), 4b, M (matrix protein), G (glycoprotein) and L (polymerase) genes. The genome consists of 12,807 nucleotides and encodes six genes in the order 3' leader-N-4a(P)-4b-M-G-L-5' trailer. Sequences were derived from clones of a cDNA library from LNYV genomic RNA and from fragments amplified using reverse transcription-polymerase chain reaction. The 4a protein has a low isoelectric point characteristic for rhabdovirus phosphoproteins. The 4b protein has significant sequence similarities with the movement proteins of capillo- and trichoviruses and may be involved in cell-to-cell movement. The putative G protein sequence contains a predicted 25 amino acids signal peptide and endopeptidase cleavage site, three predicted glycosylation sites and a putative transmembrane domain. The deduced L protein sequence shows similarities with the L proteins of other plant rhabdoviruses and contains polymerase module motifs characteristic for RNA-dependent RNA polymerases of negative-strand RNA viruses. Phylogenetic analysis of this motif among rhabdoviruses placed LNYV in a group with other sequenced cytorhabdoviruses, most closely related to Strawberry crinkle virus. (c) 2005 Elsevier B.V. All rights reserved.

Effect of sequence variation in Plasmodium falciparum Histidine-Rich protein 2 on binding of specific monoclonal antibodies: Implications for rapid diagnostic tests for malaria

This Article Right arrow Full Text Right arrow Full Text (PDF) Right arrow Supplemental material Right arrow Alert me when this article is cited Right arrow Alert me if a correction is posted Services Right arrow Similar articles in this journal Right arrow Similar articles in PubMed Right arrow Alert me to new issues of the journal Right arrow Download to citation manager Right arrow Reprints and Permissions Right arrow Copyright Information Right arrow Books from ASM Press Right arrow MicrobeWorld Citing Articles Right arrow Citing Articles via HighWire Right arrow Citing Articles via Google Scholar Google Scholar Right arrow Articles by Lee, N. Right arrow Articles by McCarthy, J. Right arrow Search for Related Content PubMed Right arrow PubMed Citation Right arrow Articles by Lee, N. Right arrow Articles by McCarthy, J. Right arrow Pubmed/NCBI databases * Substance via MeSH Previous Article | Next Article Journal of Clinical Microbiology, August 2006, p. 2773-2778, Vol. 44, No. 8 0095-1137/06/$08.00+0 doi:10.1128/JCM.02557-05 Copyright © 2006, American Society for Microbiology. All Rights Reserved. Effect of Sequence Variation in Plasmodium falciparum Histidine- Rich Protein 2 on Binding of Specific Monoclonal Antibodies: Implications for Rapid Diagnostic Tests for Malaria{dagger} Nelson Lee,1,2 Joanne Baker,2 Kathy T. Andrews,1 Michelle L. Gatton,1,3 David Bell,4 Qin Cheng,2,3 and James McCarthy1* Australian Centre for International and Tropical Health and Nutrition, Queensland Institute of Medical Research and School of Population Health, University of Queensland, Queensland, Australia,1 Department of Drug Resistance and Diagnostics, Australian Army Malaria Institute, Brisbane, Australia,2 Malaria Drug Resistance and Chemotherapy, Queensland Institute of Medical Research, Queensland, Australia,3 World Health Organization, Regional Office for the Western Pacific, Manila, Philippines4 Received 8 December 2005/ Returned for modification 23 February 2006/ Accepted 26 May 2006 The ability to accurately diagnose malaria infections, particularly in settings where laboratory facilities are not well developed, is of key importance in the control of this disease. Rapid diagnostic tests (RDTs) offer great potential to address this need. Reports of significant variation in the field performance of RDTs based on the detection of Plasmodium falciparum histidine-rich protein 2 (HRP2) (PfHRP2) and of significant sequence polymorphism in PfHRP2 led us to evaluate the binding of four HRP2-specific monoclonal antibodies (MABs) to parasite proteins from geographically distinct P. falciparum isolates, define the epitopes recognized by these MABs, and relate the copy number of the epitopes to MAB reactivity. We observed a significant difference in the reactivity of the same MAB to different isolates and between different MABs tested with single isolates. When the target epitopes of three of the MABs were determined and mapped onto the peptide sequences of the field isolates, significant variability in the frequency of these epitopes was observed. These findings support the role of sequence variation as an explanation for variations in the performance of HRP2-based RDTs and point toward possible approaches to improve their diagnostic sensitivities