Restriction site-associated DNA sequencing, genotyping error estimation and de novo assembly optimization for population genetic inference.

Restriction site-associated DNA sequencing (RADseq) provides researchers with the ability to record genetic polymorphism across thousands of loci for nonmodel organisms, potentially revolutionizing the field of molecular ecology. However, as with other genotyping methods, RADseq is prone to a number of sources of error that may have consequential effects for population genetic inferences, and these have received only limited attention in terms of the estimation and reporting of genotyping error rates. Here we use individual sample replicates, under the expectation of identical genotypes, to quantify genotyping error in the absence of a reference genome. We then use sample replicates to (i) optimize de novo assembly parameters within the program Stacks, by minimizing error and maximizing the retrieval of informative loci; and (ii) quantify error rates for loci, alleles and single-nucleotide polymorphisms. As an empirical example, we use a double-digest RAD data set of a nonmodel plant species, Berberis alpina, collected from high-altitude mountains in Mexico.

Clinical and genetic investigation of a large Tunisian family with complete achromatopsia: identification of a new nonsense mutation in GNAT2 gene.

Complete achromatopsia is a rare autosomal recessive disease associated with CNGA3, CNGB3, GNAT2 and PDE6C mutations. This retinal disorder is characterized by complete loss of color discrimination due to the absence or alteration of the cones function. The purpose of the present study was the clinical and the genetic characterization of achromatopsia in a large consanguineous Tunisian family. Ophthalmic evaluation included a full clinical examination, color vision testing and electroretinography. Linkage analysis using microsatellite markers flanking CNGA3, CNGB3, GNAT2 and PDE6C genes was performed. Mutations were screened by direct sequencing. A total of 12 individuals were diagnosed with congenital complete achromatopsia. They are members of six nuclear consanguineous families belonging to the same large consanguineous family. Linkage analysis revealed linkage to GNAT2. Mutational screening of GNAT2 revealed three intronic variations c.119-69G>C, c.161+66A>T and c.875-31G>C that co-segregated with a novel mutation p.R313X. An identical GNAT2 haplotype segregating with this mutation was identified, indicating a founder mutation. All patients were homozygous for the p.R313X mutation. This is the first report of the clinical and genetic investigation of complete achromatopsia in North Africa and the largest family with recessive achromatopsia involving GNAT2; thus, providing a unique opportunity for genotype-phenotype correlation for this extremely rare condition.

Mutation screening of the glutamate cysteine ligase modifier (GCLM) gene in patients with schizophrenia.

BACKGROUND: Experimental evidences show that glutathione and its rate-limiting synthesizing enzyme, the glutamate-cysteine ligase (GCL), are involved in the pathogenesis of schizophrenia. Furthermore, genetic association has been previously reported between two single nucleotide polymorphisms lying in noncoding regions of glutamate cysteine ligase modifier (GCLM) gene, which specifies for the modifier subunit of GCL and schizophrenia. OBJECTIVE: We wanted to investigate the presence of GCLM true functional mutations, likely in linkage disequilibrium with the previously identified single nucleotide polymorphism alleles, in the same set of cases that allowed the detection of the original association signal. METHODS: We screened all the coding regions of GCLM and their intronic flanking vicinities in 353 patients with schizophrenia by direct DNA sequencing. RESULTS: Ten sequence variations were identified, five of which were not previously described. None of these DNA changes was within the GCLM coding sequence and in-silico analysis failed to indicate functional impairment induced by these variations. Furthermore, screening of normal controls and downstream statistical analyses revealed no significant relationship of any of these DNA variants with schizophrenia. CONCLUSION: It is unlikely that functional mutations in the GCLM gene could play a major role in genetic predisposition to schizophrenia and further studies will be required to assess its etiological function in the disease.

Behind the Scenes of an Ant Genome Project

Dramatic improvements in DNA sequencing technologies have led to amore than 1,000-fold reduction in sequencing costs over the past five years.Genome-wide research approaches can thus now be applied beyond medicallyrelevant questions to examine the molecular-genetic basis of behavior,development and unique life histories in almost any organism. A first step foran emerging model organism is usually establishing a reference genomesequence. I offer insight gained from the fire ant genome project. First, I detailhow the project came to be and how sequencing, assembly and annotationstrategies were chosen. Subsequently, I describe some of the issues linked toworking with data from recently sequenced genomes. Finally, I discuss anapproach undertaken in a follow-up project based on the fire ant genomesequence.

Neoadjuvant chemoradiotherapy with or without panitumumab in patients with wild-type KRAS, locally advanced rectal cancer (LARC): a randomized, multicenter, phase II trial SAKK 41/07.

BACKGROUND: We conducted a randomized, phase II, multicenter study to evaluate the anti-epidermal growth factor receptor (EGFR) mAb panitumumab (P) in combination with chemoradiotherapy (CRT) with standard-dose capecitabine as neoadjuvant treatment for wild-type KRAS locally advanced rectal cancer (LARC). PATIENTS AND METHODS: Patients with wild-type KRAS, T3-4 and/or N+ LARC were randomly assigned to receive CRT with or without P (6 mg/kg). The primary end-point was pathological near-complete or complete tumor response (pNC/CR), defined as grade 3 (pNCR) or 4 (pCR) histological regression by Dworak classification (DC). RESULTS: Forty of 68 patients were randomly assigned to P + CRT and 28 to CRT. pNC/CR was achieved in 21 patients (53%) treated with P + CRT [95% confidence interval (CI) 36%-69%] versus 9 patients (32%) treated with CRT alone (95% CI: 16%-52%). pCR was achieved in 4 (10%) and 5 (18%) patients, and pNCR in 17 (43%) and 4 (14%) patients. In immunohistochemical analysis, most DC 3 cells were not apoptotic. The most common grade ≥3 toxic effects in the P + CRT/CRT arm were diarrhea (10%/6%) and anastomotic leakage (15%/4%). CONCLUSIONS: The addition of panitumumab to neoadjuvant CRT in patients with KRAS wild-type LARC resulted in a high pNC/CR rate, mostly grade 3 DC. The results of both treatment arms exceeded prespecified thresholds. The addition of panitumumab increased toxicity.

Methylation profile of TP53 regulatory pathway and mtDNA alterations in breast cancer patients lacking TP53 mutations.

The present study investigated promoter hypermethylation of TP53 regulatory pathways providing a potential link between epigenetic changes and mitochondrial DNA (mtDNA) alterations in breast cancer patients lacking a TP53 mutation. The possibility of using the cancer-specific alterations in serum samples as a blood-based test was also explored. Triple-matched samples (cancerous tissues, matched adjacent normal tissues and serum samples) from breast cancer patients were screened for TP53 mutations, and the promoter methylation profile of P14(ARF), MDM2, TP53 and PTEN genes was analyzed as well as mtDNA alterations, including D-loop mutations and mtDNA content. In the studied cohort, no mutation was found in TP53 (DNA-binding domain). Comparison of P14(ARF) and PTEN methylation patterns showed significant hypermethylation levels in tumor tissues (P < 0.05 and <0.01, respectively) whereas the TP53 tumor suppressor gene was not hypermethylated (P < 0.511). The proportion of PTEN methylation was significantly higher in serum than in the normal tissues and it has a significant correlation to tumor tissues (P < 0.05). mtDNA analysis revealed 36.36% somatic and 90.91% germline mutations in the D-loop region and also significant mtDNA depletion in tumor tissues (P < 0.01). In addition, the mtDNA content in matched serum was significantly lower than in the normal tissues (P < 0.05). These data can provide an insight into the management of a therapeutic approach based on the reversal of epigenetic silencing of the crucial genes involved in regulatory pathways of the tumor suppressor TP53. Additionally, release of significant aberrant methylated PTEN in matched serum samples might represent a promising biomarker for breast cancer.

Variation in novel exons (RACEfrags) of the MECP2 gene in Rett syndrome patients and controls.

The study of transcription using genomic tiling arrays has lead to the identification of numerous additional exons. One example is the MECP2 gene on the X chromosome; using 5'RACE and RT-PCR in human tissues and cell lines, we have found more than 70 novel exons (RACEfrags) connecting to at least one annotated exon.. We sequenced all MECP2-connected exons and flanking sequences in 3 groups: 46 patients with the Rett syndrome and without mutations in the currently annotated exons of the MECP2 and CDKL5 genes; 32 patients with the Rett syndrome and identified mutations in the MECP2 gene; 100 control individuals from the same geoethnic group. Approximately 13 kb were sequenced per sample, (2.4 Mb of DNA resequencing). A total of 75 individuals had novel rare variants (mostly private variants) but no statistically significant difference was found among the 3 groups. These results suggest that variants in the newly discovered exons may not contribute to Rett syndrome. Interestingly however, there are about twice more variants in the novel exons than in the flanking sequences (44 vs. 21 for approximately 1.3 Mb sequenced for each class of sequences, p=0.0025). Thus the evolutionary forces that shape these novel exons may be different than those of neighboring sequences.

Modeling sequencing errors by combining Hidden Markov models.

Among the largest resources for biological sequence data is the large amount of expressed sequence tags (ESTs) available in public and proprietary databases. ESTs provide information on transcripts but for technical reasons they often contain sequencing errors. Therefore, when analyzing EST sequences computationally, such errors must be taken into account. Earlier attempts to model error prone coding regions have shown good performance in detecting and predicting these while correcting sequencing errors using codon usage frequencies. In the research presented here, we improve the detection of translation start and stop sites by integrating a more complex mRNA model with codon usage bias based error correction into one hidden Markov model (HMM), thus generalizing this error correction approach to more complex HMMs. We show that our method maintains the performance in detecting coding sequences.

Results of a collaborative study of the EDNAP group regarding mitochondrial DNA heteroplasmy and segregation in hair shafts.

A collaborative exercise was carried out by the European DNA Profiling Group (EDNAP) in order to evaluate the distribution of mitochondrial DNA (mtDNA) heteroplasmy amongst the hairs of an individual who displays point heteroplasmy in blood and buccal cells. A second aim of the exercise was to study reproducibility of mtDNA sequencing of hairs between laboratories using differing chemistries, further to the first mtDNA reproducibility study carried out by the EDNAP group. Laboratories were asked to type 2 sections from each of 10 hairs, such that each hair was typed by at least two laboratories. Ten laboratories participated in the study, and a total of 55 hairs were typed. The results showed that the C/T point heteroplasmy observed in blood and buccal cells at position 16234 segregated differentially between hairs, such that some hairs showed only C, others only T and the remainder, C/T heteroplasmy at varying ratios. Additionally, differential segregation of heteroplasmic variants was confirmed in independent extracts at positions 16093 and the poly(C) tract at 302-309, whilst a complete A-G transition was confirmed at position 16129 in one hair. Heteroplasmy was observed at position 16195 on both strands of a single extract from one hair segment, but was not observed in the extracts from any other segment of the same hair. Similarly, heteroplasmy at position 16304 was observed on both strands of a single extract from one hair. Additional variants at positions 73, 249 and the HVII poly(C) region were reported by one laboratory; as these were not confirmed in independent extracts, the possibility of contamination cannot be excluded. Additionally, the electrophoresis and detection equipment used by this laboratory was different to those of the other laboratories, and the discrepancies at position 249 and the HVII poly(C) region appear to be due to reading errors that may be associated with this technology. The results, and their implications for forensic mtDNA typing, are discussed in the light of the biology of hair formation.

High throughput sequencing and proteomics to identify immunogenic proteins of a new pathogen: the dirty genome approach.

With the availability of new generation sequencing technologies, bacterial genome projects have undergone a major boost. Still, chromosome completion needs a costly and time-consuming gap closure, especially when containing highly repetitive elements. However, incomplete genome data may be sufficiently informative to derive the pursued information. For emerging pathogens, i.e. newly identified pathogens, lack of release of genome data during gap closure stage is clearly medically counterproductive. We thus investigated the feasibility of a dirty genome approach, i.e. the release of unfinished genome sequences to develop serological diagnostic tools. We showed that almost the whole genome sequence of the emerging pathogen Parachlamydia acanthamoebae was retrieved even with relatively short reads from Genome Sequencer 20 and Solexa. The bacterial proteome was analyzed to select immunogenic proteins, which were then expressed and used to elaborate the first steps of an ELISA. This work constitutes the proof of principle for a dirty genome approach, i.e. the use of unfinished genome sequences of pathogenic bacteria, coupled with proteomics to rapidly identify new immunogenic proteins useful to develop in the future specific diagnostic tests such as ELISA, immunohistochemistry and direct antigen detection. Although applied here to an emerging pathogen, this combined dirty genome sequencing/proteomic approach may be used for any pathogen for which better diagnostics are needed. These genome sequences may also be very useful to develop DNA based diagnostic tests. All these diagnostic tools will allow further evaluations of the pathogenic potential of this obligate intracellular bacterium.