Identifying Signatures of Selection in Genetic Time Series

Both genetic drift and natural selection cause the frequencies of alleles in a population to vary over time. Discriminating between these two evolutionary forces, based on a time series of samples from a population, remains an outstanding problem with increasing relevance to modern data sets. Even in the idealized situation when the sampled locus is independent of all other loci, this problem is difficult to solve, especially when the size of the population from which the samples are drawn is unknown. A standard χ²-based likelihood-ratio test was previously proposed to address this problem. Here we show that the χ²-test of selection substantially underestimates the probability of type I error, leading to more false positives than indicated by its P-value, especially at stringent P-values. We introduce two methods to correct this bias. The empirical likelihood-ratio test (ELRT) rejects neutrality when the likelihood-ratio statistic falls in the tail of the empirical distribution obtained under the most likely neutral population size. The frequency increment test (FIT) rejects neutrality if the distribution of normalized allele-frequency increments exhibits a mean that deviates significantly from zero. We characterize the statistical power of these two tests for selection, and we apply them to three experimental data sets. We demonstrate that both ELRT and FIT have power to detect selection in practical parameter regimes, such as those encountered in microbial evolution experiments. Our analysis applies to a single diallelic locus, assumed independent of all other loci, which is most relevant to full-genome selection scans in sexual organisms, and also to evolution experiments in asexual organisms as long as clonal interference is weak. Different techniques will be required to detect selection in time series of cosegregating linked loci.     

Automated Classification and Cluster Visualization of Genotypes Derived from High Resolution Melt Curves

Introduction

High Resolution Melting (HRM) following PCR has been used to identify DNA genotypes. Fluorescent dyes bounded to double strand DNA lose their fluorescence with increasing temperature, yielding different signatures for different genotypes. Recent software tools have been made available to aid in the distinction of different genotypes, but they are not fully automated, used only for research purposes, or require some level of interaction or confirmation from an analyst.

Materials and Methods

We describe a fully automated machine learning software algorithm that classifies unknown genotypes. Dynamic melt curves are transformed to multidimensional clusters of points whereby a training set is used to establish the distribution of genotype clusters. Subsequently, probabilistic and statistical methods were used to classify the genotypes of unknown DNA samples on 4 different assays (40 VKORC1, CYP2C9*2, CYP2C9*3 samples in triplicate, and 49 MTHFR c.665C>T samples in triplicate) run on the Roche LC480. Melt curves of each of the triplicates were genotyped separately.

Results

Automated genotyping called 100% of VKORC1, CYP2C9*3 and MTHFR c.665C>T samples correctly. 97.5% of CYP2C9*2 melt curves were genotyped correctly with the remaining 2.5% given a no call due to the inability to decipher 3 melt curves in close proximity as either homozygous mutant or wild-type with greater than 99.5% posterior probability.

Conclusions

We demonstrate the ability to fully automate DNA genotyping from HRM curves systematically and accurately without requiring any user interpretation or interaction with the data. Visualization of genotype clusters and quantification of the expected misclassification rate is also available to provide feedback to assay scientists and engineers as changes are made to the assay or instrument.     

Applicability of an Arrhenius Model for the Combined Effect of Temperature and CO2 Packaging on the Spoilage Microflora of Fish

The temperature behavior of the natural microflora on the Mediterranean fish red mullet (Mullus barbatus) was examined as a case study. The growth of the spoilage bacteria Pseudomonas spp., Shewanella putrefaciens, Brochothrix thermosphacta, and lactic acid bacteria was modeled as a function of temperature and the concentration of carbon dioxide in modified atmosphere packaging. Combined models were developed and comparatively assessed based on polynomial, Belehradek, and Arrhenius equations. The activation energy parameter of the Arrhenius model, E_A, was independent of the packaging atmosphere and ranged from 75 to 85 kJ/mol for the different bacteria, whereas the preexponential constant decreased exponentially with the packaging CO₂ concentration. We evaluated the applicability of the models developed by using experimental bacterial growth rates obtained from 42 independent experiments performed with three Mediterranean fish species and growth rates predicted from the models under the same temperature and packaging conditions. The accuracy factor and bias factor were used as statistical tools for evaluation, and the developed Arrhenius model and the Belehradek model were judged satisfactory overall.     

Selection of appropriate reference genes for gene expression studies by quantitative real-time polymerase chain reaction in cucumber

Quantitative real-time polymerase chain reaction (QRT-PCR) has become one of the most widely used methods for gene expression analysis. However, the expression profile of a target gene may be misinterpreted due to unstable expression of the reference genes under different experimental conditions. Thus, a systematic evaluation of these reference genes is necessary before experiments are performed. In this study, 10 putative reference genes were chosen for identifying expression stability using geNorm, NormFinder, and BestKeeper statistical algorithms in 12 different cucumber sample pools, including those from different plant tissues and from plants treated with hormones and abiotic stresses. EF1α and UBI-ep exhibited the most stable expression across all of the tested cucumber samples. In different tissues, in addition to expression of EF1α and UBI-ep, the expression of TUA was also stable and was considered as an appropriate reference gene. Evaluation of samples treated with different hormones revealed that TUA and UBI-ep were the most stably expressed genes. However, for abiotic stress treatments, only EF1α showed a relatively stable expression level. In conclusion, TUA, UBI-ep, and EF1α will be particularly helpful for reliable QRT-PCR data normalization in these types of samples. This study also provides guidelines for selecting different reference genes under different conditions.     

Evaluation of appropriate reference genes for gene expression studies in pepper by quantitative real-time PCR

Quantitative real-time polymerase chain reaction (qRT-PCR) has been extensively used in several plant species as an accurate technique for gene expression analysis. However, the expression level of a target gene may be misconstrued due to unstable expression of the reference genes under different experimental conditions. Therefore, it is necessary to systematically evaluate these reference genes before experiments are conducted. Recently, more and more studies have focused on gene expression in pepper (Capsicum annuum L.). In this study, ten putative reference genes were chosen to identify expression stability by using geNorm and NormFinder statistical algorithms in ten different pepper sample pools, including those from different plant tissues (root, stem, leaf and flower) and from plants treated with hormones (salicylic acid and gibberellic acid) and abiotic stresses (cold, heat, salt and drought). EF1?? and UEP exhibited the most stable expression across all of the tested pepper samples. For abiotic stress or different hormone treatment, the ranking of candidate reference genes was not completely consistent, except for EF1?? which showed a relatively stable expression level. For different tissues, the expression of Actin1 was stable and it was considered an appropriate reference gene. It is concluded that EF1??, UEP and Actin1 are suitable reference genes for reliable qRT-PCR data normalization for the tissues and experimental conditions used in this experiment.     

Kinetic analysis of P-700 photoconversion: Effect of secondary electron donation and plastocyanin inhibition

The kinetics of P-700 photoconversion under weak continuous actinic illumination were quantitatively analyzed to provide information on the relative absorption cross-section σ_PSI of the light-harvesting pigments associated with photosystem I and on the number of electrons stored between the two photosystems in dark-adapted chloroplasts. The theory of chemical kinetics for a system of monomolecular consecutive first-order reactions is reviewed briefly to provide support for the experimental approach taken. A complete inhibition of plastocyanin by cyanide eliminated all secondary electron donation to P-700⁺ and allowed the registration of the exponential (monomolecular) P-700 photoconversion at room temperature. The rate constant K_p-700 of the exponential kinetics was independent of the ionic (± Mg²⁺) and osmotic (± sucrose) strength of the chloroplast suspension medium, and of the oxidation-reduction state of photosystem II. The extent of plastocyanin inhibition in partially inhibited samples was greater under low ionic and low osmotic conditions. In dark-adapted chloroplast samples that were not cyanide treated, the number of electrons stored between the two photosystems was 3.9 ± 0.2 and independent of divalent cations. It is concluded that plastocyanin inhibition by cyanide is favored under low ionic and low osmotic conditions. The Mg²⁺ ion and redox state of photosystem II-independent photoconversion of P-700 does not support significant changes in the spillover of excitation from photosystem II to photosystem I in isolated chloroplasts.     

SynteBase/SynteView: a tool to visualize gene order conservation in prokaryotic genomes

Background

Most microarray studies are made using labelling with one or two dyes which allows the hybridization of one or two samples on the same slide. In such experiments, the most frequently used dyes areCy3 andCy5. Recent improvements in the technology (dye-labelling, scanner and, image analysis) allow hybridization up to four samples simultaneously. The two additional dyes areAlexa488 andAlexa494. The triple-target or four-target technology is very promising, since it allows more flexibility in the design of experiments, an increase in the statistical power when comparing gene expressions induced by different conditions and a scaled down number of slides. However, there have been few methods proposed for statistical analysis of such data. Moreover the lowess correction of the global dye effect is available for only two-color experiments, and even if its application can be derived, it does not allow simultaneous correction of the raw data.

Results

We propose a two-step normalization procedure for triple-target experiments. First the dye bleeding is evaluated and corrected if necessary. Then the signal in each channel is normalized using a generalized lowess procedure to correct a global dye bias. The normalization procedure is validated using triple-self experiments and by comparing the results of triple-target and two-color experiments. Although the focus is on triple-target microarrays, the proposed method can be used to normalizepdifferently labelled targets co-hybridized on a same array, for any value ofpgreater than 2.

Conclusion

The proposed normalization procedure is effective: the technical biases are reduced, the number of false positives is under control in the analysis of differentially expressed genes, and the triple-target experiments are more powerful than the corresponding two-color experiments. There is room for improving the microarray experiments by simultaneously hybridizing more than two samples.     

Metal ions-nucleobases interactions: preparation and crystal structures of trichloroadeninium zinc(ii)(form ii) and a similar zinc complex of arprinocid [6-amino-9-(2-chloro-6-fluorobenzyl)purine]

Two zinc complexes—trichloroadeninium zinc(II)(Form 11), C₅H₆N₅Cl₃Zn [structure(I)] and a similar complex of Arprinocid, (6-amino-9-(2-chloro-6-fluorobenzyl)purine], C₁₂H₁₀N₅FCl₄Zn [structure(II)]—have been prepared Structure(I) crystallizes in the space group P2₁/c with a = 8.223(1)Å, b = 6.755(1) Å, c = 18.698(3) Å, β = 96.10(2)°,and Z = 4. Structure(II) crystallizes in the space group P2₁/c with a = 8.209(2) Å, b = 6.421(8) Å, c = 31.794(8) Å, β = 90.76(2)°, and Z = 4. Both of these structures were solved by the heavy atom method using diffractometric data and refined to R = 0.028 [structure(I)] and 0.038 [structure(II)]. Zinc with a distorted tetrahedral coordination having three chlorines and N(7) as ligators, protonation of the adenine moiety at N(1), dissymmetry of exocyclic angles at N(7), and an interligand hydrogen bond (“indirect chelation”) involving one of the three chlorines, coordinated to zinc and a proton of the exocylic amino group are the striking features common to both structures. Similar types of indirect chelation as observed in the different complexes of purines have been discussed. The zinc ion deviates from the imidazole plane by 0.412 Å in structure(I) and 0.524 Å in Structure(II). The imidazol and pyrimidine planes fold about the C(4)-C(5) bond by 2.4° in strctur(I) and 3.8° in structure(II). In structure(I), inversion related molecules are paired through N(9)-H…N(3) hydrogen bonds. N-H…Cl hydrogen bonds and C(8)-H…Cl interactions have been observed in both structures.     

Co-segregation of trichorhinophalangeal syndrome with a t(8;13)(q23.3;q21.31) familial translocation that appears to increase TRPS1 gene expression

Trichorhinophalangeal syndrome type I (TRPS I) is a rare autosomal dominant syndrome caused by haploinsufficiency of TRPS1 due to point mutations or deletions. Here, we report the first familial TRPS I due to a t(8;13)(q23.3;q21.31) translocation breakpoint <100 kb from the 5′ end of TRPS1. Based on the additional abnormalities observed exclusively in the index patient that are mainly compatible with clinical features of TRPS, her phenotype was defined as expanded TRPS I including brain malformations and intellectual disability. Initial analyses did not reveal any genetic defect affecting TRPS1 or any genomic alteration within the breakpoint regions or elsewhere in the genome. The pathogenic chromosome 8q23.3 breakpoint is at position g.116,768,309_116,768,310 within a transposon type I element, 87 kb from the TRPS1 5′ end. The 13q21.31 breakpoint is within a tandem repeat region at position g.65,101,509_65,101,510 (genome assembly GRCh37/hg19). This breakpoint is flanked by protocadherin 9 (PCDH9) and protocadherin 20 (PCDH20). As an outcome of the translocation, an evolutionarily conserved non-coding VISTA enhancer element from 13q21.31 is placed within the TRPS1 5′ region, 1,294 bp from the breakpoint. The increased expression of TRPS1 found by three independent methods is most probably translocation allele derived and driven by the translocated enhancer element. The index patient’s expanded phenotype presumably involves the epithelial-to-mesenchymal transition pathway that may be due to TRPS1 overexpression. Together, these findings support that the reported translocation-associated phenotypes are “cis-ruption” and TRPS1 overexpression related, the latter most probably caused by the novel enhancer element in the TRPS1 5′ region.     

Accuracy of Urine Circulating Cathodic Antigen (CCA) Test for Schistosoma mansoni Diagnosis in Different Settings of C?te d'Ivoire

Background

Promising results have been reported for a urine circulating cathodic antigen (CCA) test for the diagnosis of Schistosoma mansoni. We assessed the accuracy of a commercially available CCA cassette test (designated CCA-A) and an experimental formulation (CCA-B) for S. mansoni diagnosis.

Methodology

We conducted a cross-sectional survey in three settings of Côte d''Ivoire: settings A and B are endemic for S. mansoni, whereas S. haematobium co-exists in setting C. Overall, 446 children, aged 8–12 years, submitted multiple stool and urine samples. For S. mansoni diagnosis, stool samples were examined with triplicate Kato-Katz, whereas urine samples were tested with CCA-A. The first stool and urine samples were additionally subjected to an ether-concentration technique and CCA-B, respectively. Urine samples were examined for S. haematobium using a filtration method, and for microhematuria using Hemastix dipsticks.

Principal Findings

Considering nine Kato-Katz as diagnostic ‘gold’ standard, the prevalence of S. mansoni in setting A, B and C was 32.9%, 53.1% and 91.8%, respectively. The sensitivity of triplicate Kato-Katz from the first stool and a single CCA-A test was 47.9% and 56.3% (setting A), 73.9% and 69.6% (setting B), and 94.2% and 89.6% (setting C). The respective sensitivity of a single CCA-B was 10.4%, 29.9% and 75.0%. The ether-concentration technique showed a low sensitivity for S. mansoni diagnosis (8.3–41.0%). The specificity of CCA-A was moderate (76.9–84.2%); CCA-B was high (96.7–100%). The likelihood of a CCA-A color reaction increased with higher S. mansoni fecal egg counts (odds ratio: 1.07, p<0.001). A concurrent S. haematobium infection or the presence of microhematuria did not influence the CCA-A test results for S. mansoni diagnosis.

Conclusion/Significance

CCA-A showed similar sensitivity than triplicate Kato-Katz for S. mansoni diagnosis with no cross-reactivity to S. haematobium and microhematuria. The low sensitivity of CCA-B in our study area precludes its use for S. mansoni diagnosis.     

Use of a Single, Triplicate Arbitrarily Primed-PCR Procedure for Molecular Fingerprinting of Lactic Acid Bacteria

Arbitrarily primed (AP)-PCR can be used to generate characteristic DNA fingerprint patterns. However, small changes in reaction conditions can cause band irreproducibility. In this study, a single methodology encompassing triplicate reactions, which were intentionally exposed to three different annealing temperatures, enabled bands that were reproducibly generated to be recognized. A single triplicate AP-PCR (TAP-PCR) procedure, using an 18-mer primer, was developed and used to fingerprint representative isolates from the major genera of lactic acid bacteria and Bifidobacterium to the strain level.     

Statistical tests for differential expression in cDNA microarray experiments

Extracting biological information from microarray data requires appropriate statistical methods. The simplest statistical method for detecting differential expression is the t test, which can be used to compare two conditions when there is replication of samples. With more than two conditions, analysis of variance (ANOVA) can be used, and the mixed ANOVA model is a general and powerful approach for microarray experiments with multiple factors and/or several sources of variation.     

Equalizer reduces SNP bias in Affymetrix microarrays

Background

Gene expression microarrays measure the levels of messenger ribonucleic acid (mRNA) in a sample using probe sequences that hybridize with transcribed regions. These probe sequences are designed using a reference genome for the relevant species. However, most model organisms and all humans have genomes that deviate from their reference. These variations, which include single nucleotide polymorphisms, insertions of additional nucleotides, and nucleotide deletions, can affect the microarray’s performance. Genetic experiments comparing individuals bearing different population-associated single nucleotide polymorphisms that intersect microarray probes are therefore subject to systemic bias, as the reduction in binding efficiency due to a technical artifact is confounded with genetic differences between parental strains. This problem has been recognized for some time, and earlier methods of compensation have attempted to identify probes affected by genome variants using statistical models. These methods may require replicate microarray measurement of gene expression in the relevant tissue in inbred parental samples, which are not always available in model organisms and are never available in humans.

Results

By using sequence information for the genomes of organisms under investigation, potentially problematic probes can now be identified a priori. However, there is no published software tool that makes it easy to eliminate these probes from an annotation. I present equalizer, a software package that uses genome variant data to modify annotation files for the commonly used Affymetrix IVT and Gene/Exon platforms. These files can be used by any microarray normalization method for subsequent analysis. I demonstrate how use of equalizer on experiments mapping germline influence on gene expression in a genetic cross between two divergent mouse species and in human samples significantly reduces probe hybridization-induced bias, reducing false positive and false negative findings.

Conclusions

The equalizer package reduces probe hybridization bias from experiments performed on the Affymetrix microarray platform, allowing accurate assessment of germline influence on gene expression.     

Experiments and ANFIS modelling for the biodegradation of penicillin-G wastewater using anaerobic hybrid reactor

The performance of an anaerobic hybrid reactor (AHR) for treating penicillin-G wastewater was investigated at the ambient temperatures of 30-35 °C for 245 days in three phases. The experimental data were analysed by adopting an adaptive network-based fuzzy inference system (ANFIS) model, which combines the merits of both fuzzy systems and neural network technology. The statistical quality of the ANFIS model was significant due to its high correlation coefficient R² between experimental and simulated COD values. The R² was found to be 0.9718, 0.9268 and 0.9796 for the I, II and III phases, respectively. Furthermore, one to one correlation among the simulated and observed values was also observed. The results showed the proposed ANFIS model was well performed in predicting the performance of AHR.     

A Lab Assembled Microcontroller-Based Sensor Module for Continuous Oxygen Measurement in Portable Hypoxia Chambers

Background

Hypoxia-based cell culture experiments are routine and essential components of in vitro cancer research. Most laboratories use low-cost portable modular chambers to achieve hypoxic conditions for cell cultures, where the sealed chambers are purged with a gas mixture of preset O₂ concentration. Studies are conducted under the assumption that hypoxia remains unaltered throughout the 48 to 72 hour duration of such experiments. Since these chambers lack any sensor or detection system to monitor gas-phase O₂, the cell-based data tend to be non-uniform due to the ad hoc nature of the experimental setup.

Methodology

With the availability of low-cost open-source microcontroller-based electronic project kits, it is now possible for researchers to program these with easy-to-use software, link them to sensors, and place them in basic scientific apparatus to monitor and record experimental parameters. We report here the design and construction of a small-footprint kit for continuous measurement and recording of O₂ concentration in modular hypoxia chambers. The low-cost assembly (US$135) consists of an Arduino-based microcontroller, data-logging freeware, and a factory pre-calibrated miniature O₂ sensor. A small, intuitive software program was written by the authors to control the data input and output. The basic nature of the kit will enable any student in biology with minimal experience in hobby-electronics to assemble the system and edit the program parameters to suit individual experimental conditions.

Results/Conclusions

We show the kit’s utility and stability of data output via a series of hypoxia experiments. The studies also demonstrated the critical need to monitor and adjust gas-phase O₂ concentration during hypoxia-based experiments to prevent experimental errors or failure due to partial loss of hypoxia. Thus, incorporating the sensor-microcontroller module to a portable hypoxia chamber provides a researcher a capability that was previously available only to labs with access to sophisticated (and expensive) cell culture incubators.     

Isolation of bacteriophage MX4, a generalized transducing phage for Myxococcus xanthus.

The structure of α-chitin has been determined by X-ray diffraction, based on the intensity data from deproteinized lobster tendon. Least-squares refinement shows that adjacent chains have alternating sense (i.e. are antiparallel). In addition, there is a statistical distribution of side-chain orientations, such that all the hydroxyl groups form hydrogen bonds. The unit cell is orthorhombic with dimensions a = 0.474 ± 0.001 nm, b = 1.886 ± 0.002 nm and c = 1.032 ± 0.002 nm (fiber axis); the space group is P2₁2₁2₁ and the cell contains disaccharide sections of the two chains passing through the center and corner of the ab projection. The chains form hydrogen-bonded sheets linked by CO…HN bonds approximately parallel to the a axis, and each chain has an O-3′H…O.5 intramolecular hydrogen bond, similar to that in cellulose. Adjacent chains along the ab diagonal have different conformations for the CH₂OH groups: on one chain these groups form O.6H…O.6′ intermolecular hydrogen bonds to the CH₂OH group on the adjacent chain along the ab diagonal. The latter group is oriented to form an intramolecular O.6′H…O.7 bond to the carboxyl oxygen on the next residue. The results indicate that a statistical mixture of CH₂OH orientations is present, equivalent to half oxygens on each residue, each forming inter- and intramolecular hydrogen bonds. As a result the structure contains two types of amide groups, which differ in their hydrogen bonding, and account for the splitting of the amide I band in the infrared spectrum. The Inability of this chitin polymorph to swell on soaking in water is explained by the extensive intermolecular hydrogen bonding.     

The DNA sequence on bacteriophage G4 recognized by the Escherichia coli B restriction enzyme.

Bacteriophage G4 possesses a single EcoB site located in the overlap between restriction fragments HinfI-12 and HaeIII-6. The sequence 5′-T-G-A … 8N … T-G-C-T occurs once in this segment and nowhere else in the DNA sequence of G4. Four independent G4 mutants that were not restricted by Escherichia coli B possessed the sequence 5′-T-G-A … 8N … T-G-C-C. The common sequence shared by the previously mapped EcoB sites on φXsB1, simian virus 40, f1, and fd DNAs is 5′-T-G-A … 8N … T-G-C-T … 9N … T. However, the sequence in the region of the G4 EcoB site contains an A instead of the final T conserved in these other examples. When the G4 EcoB site is aligned with the other EcoB sites, there are no conserved residues within 50 bases of the common sequence, 5′-T-G-A … 8N … T-G-C-T, except for those seven residues. The analysis of the EcoB site on G4 provides further evidence that only those seven bases are recognized by the E. coli B restriction enzyme.