Evaluation of next generation sequencing platforms for population targeted sequencing studies

Background

Next generation sequencing (NGS) platforms are currently being utilized for targeted sequencing of candidate genes or genomic intervals to perform sequence-based association studies. To evaluate these platforms for this application, we analyzed human sequence generated by the Roche 454, Illumina GA, and the ABI SOLiD technologies for the same 260 kb in four individuals.

Results

Local sequence characteristics contribute to systematic variability in sequence coverage (>100-fold difference in per-base coverage), resulting in patterns for each NGS technology that are highly correlated between samples. A comparison of the base calls to 88 kb of overlapping ABI 3730xL Sanger sequence generated for the same samples showed that the NGS platforms all have high sensitivity, identifying >95% of variant sites. At high coverage, depth base calling errors are systematic, resulting from local sequence contexts; as the coverage is lowered additional 'random sampling' errors in base calling occur.

Conclusions

Our study provides important insights into systematic biases and data variability that need to be considered when utilizing NGS platforms for population targeted sequencing studies.     

Acidophilic microbial communities catalyzing sludge bioleaching monitored by fluorescent in situ hybridization

Biological autotrophic sulfur oxidation processes have been proposed to remove heavy metals from wastewater treatment sludge by bioleaching. We made a characterization of the microbial population in batch and continuous sludge bioleaching reactors using fluorescent in situ hybridization of fluorescently-labeled oligonucleotidic probes targeting rRNA in a ‚top to bottom approach’. Batch incubations of sludge with 0.2% (w/v) elemental sulfur resulted in a pH value of 5. Alpha-Proteobacteria hybridizing with probe ALF1b were dominant in this incubation. Members of the Acidophilium-group (hybridizing with probe Acdp821) of Nitrospira/Leptospirillum phylum (Ntspa712 probe) and from the archaeal domain (ARCH915) were also detected. When sludge was incubated with 1% elemental sulfur in batch or continuous reactor experiments, final pH values were always below 2. Active microbial communities consisted almost exclusively of gamma-Proteobacteria (hybridizing with probe GAM42a). However, further hybridization experiments with probe Thio820 targeting Acidithiobacillus ferroxidans and Acidithiobacillus thioxidans gave negative results. A new probe, named THIO181, encompassing all known members of the genus was designed. Hybridization perfomed with THIO181 and GAM42a showed a perfect co-localization of the hybridization signals. Further hybridization experiments with probe THIO181 and THC642, specific for the species Acidithiobacillus caldus, confirmed that this bacteria was largely responsible for the sulfur oxidation reaction in our acidophilic sludge bioleaching reactors.     

Short-term temperature impact on soil heterotrophic respiration in limed agricultural soil samples

This study sought to investigate the hourly and daily timescale responses of soil CO₂ fluxes to temperature in a limed agricultural soil. Observations from different incubation experiments were compared with the results of a model combining biotic (heterotrophic respiration) and abiotic (carbonate weathering) components. Several samples were pre-incubated for 8–9 days at three temperatures (5, 15 and 25 °C) and then submitted to short-term temperature (STT) cycles (where the temperature was increased from 5 to 35 °C in 10 °C stages, with each stage being 3 h long). During the temperature cycles (hourly timescale), the soil CO₂ fluxes increased significantly with temperature under all pre-incubation temperature (PIT) treatments. A hysteresis effect and negative fluxes during cooling phases were also systematically observed. At a given hourly timescale temperature, there was a negative relationship of the CO₂ fluxes with the PIT. Using the combined model allowed the experimental results to be clearly described, including the negative fluxes and the hysteresis effect, showing the potentially large contribution of abiotic fluxes to total fluxes in limed soils, after STT changes. The fairly good agreement between the measured and simulated flux results also suggested that the biotic flux temperature sensitivity was probably unaffected by timescale (hourly or daily) or PIT. The negative relationship of the CO₂ fluxes with the PIT probably derived from very labile soil carbon depletion, as shown in the simulations. This was not, however, confirmed by soil carbon measurements, which leaves open the possibility of adaptation within the microbial community.     

Evaluation of prognostic factors following flow-cytometric DNA analysis after cytokeratin labelling: I. Breast cancer.

In gynecologic oncology valid prognostic factors are necessary to estimate the course of disease and to define biologically similar subgroups for analysis of therapeutic efficacy. The presented study is a prospective study concerning prognostic significance of DNA ploidy and S-phase fraction in breast cancer following enrichment of tumor cells by cytokeratin labelling. Epithelial cells were labeled by FITC-conjugated cytokeratin antibody (CK 5, 6, 8, and CK 17) prior to flow cytometric cell cycle analysis in 327 fresh specimens of primary breast cancer. Univariate analysis in breast cancer detected the prognostic significance of DNA-ploidy, S-phase fraction and CV (coefficient of variation) of G(0)G(1)-peak of tumor cells for clinical outcome, especially for nodal-negative patients. Multivariate analysis could not confirm prognostic evidence of DNA-ploidy and S-phase fraction.In conclusion, in breast cancer no clinical significance for determination of DNA-parameters was found.     

Cellular lithium and transepithelial transport across toad urinary bladder

Summary Toad urinary bladders were exposed on either their mucosal or serosal surfaces, or on both surfaces, to medium in which sodium was replaced completely by lithium. With mucosal lithium Ringer's, serosal sodium Ringer's, short-circuit current (SCC) declined by about 50 percent over the first 60 min and was then maintained over a further 180 min. Cellular lithium content was comparable to the sodium transport pool. With lithium Ringer's serosa, SCC was abolished over 60 to 120 min whether the mucosal cation was sodium or lithium. Measurements of cellular ionic composition revealed that the epithelial cells gained lithium from both the mucosal and serosal media. With lithium Ringer's mucosa and serosa, cells lost potassium and gained lithium and a little chloride and water, but these changes in cellular ions could not account for the current flow across the tissue under these conditions, which must, therefore, have been carried by a transepithelial movement of lithium itself. The inhibition by serosal lithium of SCC was overcome by exposure of the mucosal surface of the bladders to amphotericin B. Thus it reflected, predominantly, an inhibition of lithium entry to the cells across the apical membrane. It is suggested that this inhibition is a consequence of cellular lithium accumulation.     

The N-glycosylation of classical swine fever virus E2 glycoprotein extracellular domain expressed in the milk of goat

Classical swine fever virus (CSFV) outer surface E2 glycoprotein represents an important target to induce protective immunization during infection but the influence of N-glycosylation pattern in antigenicity is yet unclear. In the present work, the N-glycosylation of the E2-CSFV extracellular domain expressed in goat milk was determined. Enzymatic N-glycans releasing, 2-aminobenzamide (2AB) labeling, weak anion-exchange and normal-phase HPLC combined with exoglycosidase digestions and mass spectrometry of 2AB-labeled and unlabeled N-glycans showed a heterogenic population of oligomannoside, hybrid and complex-type structures. The detection of two Man₈GlcNAc₂ isomers indicates an alternative active pathway in addition to the classical endoplasmic reticulum processing. N-acetyl or N-glycolyl monosialylated species predominate over neutral complex-type N-glycans. Asn207 site-specific micro-heterogeneity of the E2 most relevant antigenic and virulence site was determined by HPLC-mass spectrometry of glycopeptides. The differences in N-glycosylation with respect to the native E2 may not disturb the main antigenic domains when expressed in goat milk.     

The death domain kinase RIP1 links the immunoregulatory CD40 receptor to apoptotic signaling in carcinomas

CD40, a tumor necrosis factor (TNF) receptor family member, is widely recognized for its prominent role in the antitumor immune response. The immunostimulatory effects of CD40 ligation on malignant cells can be switched to apoptosis upon disruption of survival signals transduced by the binding of the adaptor protein TRAF6 to CD40. Apoptosis induction requires a TRAF2-interacting CD40 motif but is initiated within a cytosolic death-inducing signaling complex after mobilization of receptor-bound TRAF2 to the cytoplasm. We demonstrate that receptor-interacting protein 1 (RIP1) is an integral component of this complex and is required for CD40 ligand-induced caspase-8 activation and tumor cell killing. Degradation of the RIP1 K63 ubiquitin ligases cIAP1/2 amplifies the CD40-mediated cytotoxic effect, whereas inhibition of CYLD, a RIP1 K63 deubiquitinating enzyme, reduces it. This two-step mechanism of apoptosis induction expands our appreciation of commonalities in apoptosis regulatory pathways across the TNF receptor superfamily and provides a telling example of how TNF family receptors usurp alternative programs to fulfill distinct cellular functions.     

The gap between the concept and definitions in the Evolutionarily Significant Unit: the need to integrate neutral genetic variation and adaptive variation

The Evolutionarily Significant Unit (ESU) was conceptualized in 1986 as a conservation unit below the species level, theoretically applicable to a wide range of taxa. The concept has gained support, and various definitions or criteria, some of which are inconsistent with each other, have since been proposed. Recent critiques of the ESU have pointed out the dominance of definitions biased to the identification of long-term isolation or neutral genetic variation, which has largely ignored the adaptive components. We present here the validity of such claims and show how the ESU definitions have actually been applied in research. We surveyed scientific journals for original papers supporting ESU designations and determined who among the proponents of ESU definitions have gained wider support. Our results indicate that indeed there are inconsistencies with the original concept and with the existing definitions. Although the original concept recommended both ecological and genetic data as the basis for identification of ESUs, which reflect true evolutionary variation, recent definitions have become biased to either neutral genetic variation or adaptive variation. The definition which uses genetic data to assess neutral genetic variation (long-term isolation) has gained major support, and therefore validates the earlier claims. To bridge the gap between the original concept and the practical application, we propose the use of partial ESU and full ESU designations. The application of full ESU should be limited solely to when both information about neutral genetic variation and adaptive variation are available. In other cases, in which only a part of the variation is examined, we should use the term partial ESU (e.g., molecular-based ESU) and continue to investigate focal populations from other aspects of variations to designate full ESU.     

Functional Promiscuity of Gene Regulation by Serpentine Receptors in Dictyostelium discoideum

Serpentine receptors such as smoothened and frizzled play important roles in cell fate determination during animal development. In Dictyostelium discoideum, four serpentine cyclic AMP (cAMP) receptors (cARs) regulate expression of multiple classes of developmental genes. To understand their function, it is essential to know whether each cAR is coupled to a specific gene regulatory pathway or whether specificity results from the different developmental regulation of individual cARs. To distinguish between these possibilities, we measured gene induction in car1 car3 double mutant cell lines that express equal levels of either cAR1, cAR2, or cAR3 under a constitutive promoter. We found that all cARs efficiently mediate both aggregative gene induction by cAMP pulses and induction of postaggregative and prespore genes by persistent cAMP stimulation. Two exceptions to this functional promiscuity were observed. (i) Only cAR1 can mediate adenosine inhibition of cAMP-induced prespore gene expression, a phenomenon that was found earlier in wild-type cells. cAR1’s mediation of adenosine inhibition suggests that cAR1 normally mediates prespore gene induction. (ii) Only cAR2 allows entry into the prestalk pathway. Prestalk gene expression is induced by differentiation-inducing factor (DIF) but only after cells have been prestimulated with cAMP. We found that DIF-induced prestalk gene expression is 10 times higher in constitutive cAR2 expressors than in constitutive cAR1 or cAR3 expressors (which still have endogenous cAR2), suggesting that cAR2 mediates induction of DIF competence. Since in wild-type slugs cAR2 is expressed only in anterior cells, this could explain the so far puzzling observations that prestalk cells differentiate at the anterior region but that DIF levels are actually higher at the posterior region. After the initial induction of DIF competence, cAMP becomes a repressor of prestalk gene expression. This function can again be mediated by cAR1, cAR2, and cAR3.Recent years have seen the discovery of critical roles in animal development for serpentine receptors, which are usually coupled to heterotrimeric G proteins. The insect sigaling peptides hedgehog and wingless and their mammalian counterparts sonic hedgehog, desert hedgehog, and indian hedgehog and the wnt factors control a multitude of inductive events during all stages of embryogenesis. The hedgehog signal is detected by two different serpentine receptors, smoothened (1, 40) and patched (21, 38), whereas the wingless or wnt signal is detected by the serpentine receptor D-frizzled-2 (3). In the social amoeba Dictyostelium discoideum, serpentine cyclic AMP (cAMP) receptors (cARs) control induction of cell differentiation during the entire course of development. Starving cells secrete cAMP pulses that induce chemotaxis and expression of genes required for the aggregation process. Cells aggregate to form mounds, which ultimately transform into fruiting structures that consist of a globular spore mass supported by a column of stalk cells. cAMP induces entry into the spore differentiation pathway as well as synthesis of a lipophilic factor, differentiation-inducing factor (DIF), which induces entry into the stalk differentiation pathway (see reference 5). At an early stage of development cAMP synergizes with DIF to induce prestalk genes, but later it becomes an inhibitor of stalk gene expression (2). cARs were shown previously to mediate induction of aggregative genes by cAMP pulses (20) as well as cAMP induction of prespore genes and repression of prestalk genes (31, 37). Remarkably, the target for the latter critical step in cell fate determination is glycogen synthase kinase 3 (GSK-3), a zeste white-3 homolog, which is the target for the effects of wingless and wnt in insects and vertebrates, respectively (7, 34).Four cARs, showing 54 to 69% amino acid identity, are expressed in a stage- and cell-type-specific manner. cAR1 is predominantly expressed before and during aggregation (18). cAR3 is expressed at late aggregation, and expression is later restricted to the prespore cell population (13, 44). cAR2 and cAR4 are both expressed exclusively in the prestalk cell population after aggregation (19, 30). cAR knockout cell lines were generated to examine the role of the individual cARs in Dictyostelium development. car1 null cells neither aggregate nor express developmental genes but can be triggered to express aggregative and postaggregative genes by stimulation with cAMP (37, 39). car3 null cells aggregate and develop normally (13). car1 car3 double gene disruptants do not aggregate, and developmental gene expression cannot be restored with cAMP, indicating that cAR1 or cAR3 shows functional redundancy and that either one or the other has to be present for gene induction to occur (10, 36). car2 null cells are blocked in the mound stage, while car4 null cells show abnormal slug morphogenesis and culmination. Both lines show reduced expression of prestalk genes and enhanced expression of prespore genes (19, 29).To understand the function of the four cARs, it is essential to know whether each receptor is coupled to a specific signal transduction pathway that controls a specific cell differentiation event or whether each receptor can activate multiple cell differentiation pathways. In the latter case, it is not the presence of a specific receptor that determines whether a response occurs but the availability of the downstream signaling pathway. To determine whether individual receptors have unique functions in developmental gene expression, we examined gene regulation in cell lines that display about equal levels of cAR1, cAR2, and cAR3 in a car1 car3 mutant background. Our results show that with two exceptions, all three receptors can transduce both the excitation and adaptation components of the different cAMP-regulated gene induction events with almost equal levels of efficiency.     

Dry mass changes in germinating spores of Diplodia maydis

Summary The dry mass of two-celled Diplodia maydis spores was measured both before and after germination by quantitative interference microscopy. The dry mass of spores declined approximately 50% during germination. However, the dry mass of germinating spores plus the dry mass of their germ tubes was greater than the dry mass of spores before germination. We conclude that the germinating spores absorbed nutrients released from non-germinating spores.The dry mass of fungal spores can be estimated by weighing large numbers of spores and determining the mean from sample spore counts. Mumford and Pappelis(4) determined the total dry mass of individual spores of Fusarium roseum and the contained lipid bodies before and after spores germinated using quantitative interference microscopy. The mean spore dry mass before germination was 57 pg. Lipid bodies accounted for about 61% of that mass and decreased as spores germinated. The total dry mass of the spore and germ tube 24 hr later greatly exceeded that of the spore before germination. Quantitative interference microscopy has been used to measure the dry mass of various types of cells. Kulfinski and Pappelis (3) recently reviewed how this technique has been applied to plant cells. Technical aspects of interference microscopy have been described by Ross (6).The purpose of this study was to examine the dry mass changes in Diplodia maydis (Berk.) Sacc. with and without germ tubes through the use of interference microscopy.