Estimating Genetic Variability in Non-Model Taxa: A General Procedure for Discriminating Sequence Errors from Actual Variation

Genetic variation is the driving force of evolution and as such is of central interest for biologists. However, inadequate discrimination of errors from true genetic variation could lead to incorrect estimates of gene copy number, population genetic parameters, phylogenetic relationships and the deposition of gene and protein sequences in databases that are not actually present in any organism. Misincorporation errors in multi-template PCR cloning methods, still commonly used for obtaining novel gene sequences in non-model species, are difficult to detect, as no previous information may be available about the number of expected copies of genes belonging to multi-gene families. However, studies employing these techniques rarely describe in any great detail how errors arising in the amplification process were detected and accounted for. Here, we estimated the rate of base misincorporation of a widely-used PCR-cloning method, using a single copy mitochondrial gene from a single individual to minimise variation in the template DNA, as 1.62×10−3 errors per site, or 9.26×10−5 per site per duplication. The distribution of errors among sequences closely matched that predicted by a binomial distribution function. The empirically estimated error rate was applied to data, obtained using the same methods, from the Phospholipase A2 toxin family from the pitviper Ovophis monticola. The distribution of differences detected closely matched the expected distribution of errors and we conclude that, when undertaking gene discovery or assessment of genetic diversity using this error-prone method, it will be informative to empirically determine the rate of base misincorporation

Investigation into the annotation of protocol sequencing steps in the sequence read archive

BACKGROUND: The workflow for the production of high-throughput sequencing data from nucleic acid samples is complex. There are a series of protocol steps to be followed in the preparation of samples for next-generation sequencing. The quantification of bias in a number of protocol steps, namely DNA fractionation, blunting, phosphorylation, adapter ligation and library enrichment, remains to be determined. RESULTS: We examined the experimental metadata of the public repository Sequence Read Archive (SRA) in order to ascertain the level of annotation of important sequencing steps in submissions to the database. Using SQL relational database queries (using the SRAdb SQLite database generated by the Bioconductor consortium) to search for keywords commonly occurring in key preparatory protocol steps partitioned over studies, we found that 7.10%, 5.84% and 7.57% of all records (fragmentation, ligation and enrichment, respectively), had at least one keyword corresponding to one of the three protocol steps. Only 4.06% of all records, partitioned over studies, had keywords for all three steps in the protocol (5.58% of all SRA records). CONCLUSIONS: The current level of annotation in the SRA inhibits systematic studies of bias due to these protocol steps. Downstream from this, meta-analyses and comparative studies based on these data will have a source of bias that cannot be quantified at present

The role of RAS oncogene in survival of patients with lung cancer: a systematic review of the literature with meta-analysis

The proto-oncogene RAS, coding for a 21 kDa protein (p21), is mutated in 20% of lung cancer. However, the literature remains controversial on its prognostic significance for survival in lung cancer. We performed a systematic review of the literature with meta-analysis to assess its possible prognostic value on survival. Published studies on lung cancer assessing prognostic value of RAS mutation or p21 overexpression on survival were identified by an electronic search. After a methodological assessment, we estimated individual hazard ratios (HR) estimating RAS protein alteration or RAS mutation effect on survival and combined them using meta-analytic methods. In total, 53 studies were found eligible, with 10 concerning the same cohorts of patients. Among the 43 remaining studies, the revelation method was immunohistochemistry (IHC) in nine and polymerase chain reaction (PCR) in 34. Results in terms of survival were significantly pejorative, significantly favourable, not significant and not conclusive in 9, 1, 31, 2, respectively. In total, 29 studies were evaluable for meta-analysis but we aggregated only the 28 dealing with non-small-cell lung cancer (NSCLC) and not the only one dealing with small-cell-lung cancer (SCLC). The quality scores were not statistically significantly different between studies with or without significant results in terms of survival, allowing us to perform a quantitative aggregation. The combined HR was 1.35 (95% CI: 1.16–1.56), showing a worse survival for NSCLC with KRAS2 mutations or p21 overexpression and, particularly, in adenocarcinomas (ADC) (HR 1.59; 95% CI 1.26–2.02) and in studies using PCR (HR 1.40; 95% CI 1.18–1.65) but not in studies using IHC (HR 1.08; 95% CI 0.86–1.34). RAS appears to be a pejorative prognostic factor in terms of survival in NSCLC globally, in ADC and when it is studied by PCR

Promoter methylation of RASSF1A and DAPK and mutations of K-ras, p53, and EGFR in lung tumors from smokers and never-smokers

<p>Abstract</p> <p>Background</p> <p>Epidemiological studies indicate that some characteristics of lung cancer among never-smokers significantly differ from those of smokers. Aberrant promoter methylation and mutations in some oncogenes and tumor suppressor genes are frequent in lung tumors from smokers but rare in those from never-smokers. In this study, we analyzed promoter methylation in the <it>ras-association domain isoform A (RASSF1A) </it>and the <it>death-associated protein kinase (DAPK) </it>genes in lung tumors from patients with primarily non-small cell lung cancer (NSCLC) from the Western Pennsylvania region. We compare the results with the smoking status of the patients and the mutation status of the K-<it>ras</it>, <it>p53</it>, and <it>EGFR </it>genes determined previously on these same lung tumors.</p> <p>Methods</p> <p>Promoter methylation of the <it>RASSF1A </it>and <it>DAPK </it>genes was analyzed by using a modified two-stage methylation-specific PCR. Data on mutations of K-<it>ras</it>, <it>p53</it>, and <it>EGFR </it>were obtained from our previous studies.</p> <p>Results</p> <p>The <it>RASSF1A </it>gene promoter methylation was found in tumors from 46.7% (57/122) of the patients and was not significantly different between smokers and never-smokers, but was associated significantly in multiple variable analysis with tumor histology (p = 0.031) and marginally with tumor stage (p = 0.063). The <it>DAPK </it>gene promoter methylation frequency in these tumors was 32.8% (40/122) and did not differ according to the patients' smoking status, tumor histology, or tumor stage. Multivariate analysis adjusted for age, gender, smoking status, tumor histology and stage showed that the frequency of promoter methylation of the <it>RASSF1A </it>or <it>DAPK </it>genes did not correlate with the frequency of mutations of the K<it>-ras, p53</it>, and <it>EGFR </it>gene.</p> <p>Conclusion</p> <p>Our results showed that <it>RASSF1A </it>and <it>DAPK </it>genes' promoter methylation occurred frequently in lung tumors, although the prevalence of this alteration in these genes was not associated with the smoking status of the patients or the occurrence of mutations in the K-<it>ras</it>, <it>p53 </it>and <it>EGFR </it>genes, suggesting each of these events may represent independent event in non-small lung tumorigenesis.</p

Qualitatively and quantitatively similar effects of active and passive maternal tobacco smoke exposure on in utero mutagenesis at the HPRT locus

BACKGROUND: Induced mutagenesis in utero is likely to have life-long repercussions for the exposed fetus, affecting survival, birth weight and susceptibility to both childhood and adult-onset diseases, such as cancer. In the general population, such exposures are likely to be a consequence of the lifestyle choices of the parents, with exposure to tobacco smoke one of the most pervasive and easily documented. Previous studies attempting to establish a direct link between active smoking and levels of somatic mutation have largely discounted the effects of passive or secondary exposure, and have produced contradictory results. METHODS: Data from three studies of possible smoking effects on in utero mutagenesis at the HPRT locus were compiled and reanalyzed, alone and in combination. Where possible, passive exposure to environmental tobacco smoke was considered as a separate category of exposure, rather than being included in the non-smoking controls. Molecular spectra from these studies were reanalyzed after adjustment for reported mutation frequencies from the individual studies and the entire data set. RESULTS: A series of related studies on mutation at the X-linked HPRT locus in human newborn cord blood samples has led to the novel conclusion that only passive maternal exposure to tobacco mutagens has a significant effect on the developing baby. We performed a pooled analysis of the complete data from these studies, at the levels of both induced mutation frequency and the resulting mutational spectrum. CONCLUSION: Our analysis reveals a more commonsensical, yet no less cautionary result: both active maternal smoking and secondary maternal exposure produce quantitatively and qualitatively indistinguishable increases in fetal HPRT mutation. Further, it appears that this effect is not perceptibly ameliorated if the mother adjusts her behavior (i.e. stops smoking) when pregnancy is confirmed, although this conclusion may also be affected by continued passive exposure

An integrative approach to discovering cryptic species within the Bemisia tabaci whitefly species complex

Bemisia tabaci is a cryptic whitefly-species complex that includes some of the most damaging pests and plant-virus vectors of a diverse range of food and fibre crops worldwide. We combine experimental evidence of: (i) differences in reproductive compatibility, (ii) hybrid verification using a specific nuclear DNA marker and hybrid fertility confirmation and (iii) high-throughput sequencing-derived mitogenomes, to show that the “Mediterranean” (MED) B. tabaci comprises at least two distinct biological species; the globally invasive MED from the Mediterranean Basin and the “African silver-leafing” (ASL) from sub-Saharan Africa, which has no associated invasion records. We demonstrate that, contrary to its common name, the “ASL” does not induce squash silver-leafing symptoms and show that species delimitation based on the widely applied 3.5% partial mtCOI gene sequence divergence threshold produces discordant results, depending on the mtCOI region selected. Of the 292 published mtCOI sequences from MED/ASL groups, 158 (54%) are low quality and/or potential pseudogenes. We demonstrate fundamental deficiencies in delimiting cryptic B. tabaci species, based solely on partial sequences of a mitochondrial barcoding gene. We advocate an integrative approach to reveal the true species richness within cryptic species complexes, which is integral to the deployment of effective pest and disease management strategies