Next generation sequencing of pooled samples reveals new SNRNP200 mutations associated with retinitis pigmentosa.

The gene SNRNP200 is composed of 45 exons and encodes a protein essential for pre-mRNA splicing, the 200 kDa helicase hBrr2. Two mutations in SNRNP200 have recently been associated with autosomal dominant retinitis pigmentosa (adRP), a retinal degenerative disease, in two families from China. In this work we analyzed the entire 35-Kb SNRNP200 genomic region in a cohort of 96 unrelated North American patients with adRP. To complete this large-scale sequencing project, we performed ultra high-throughput sequencing of pooled, untagged PCR products. We then validated the detected DNA changes by Sanger sequencing of individual samples from this cohort and from an additional one of 95 patients. One of the two previously known mutations (p.S1087L) was identified in 3 patients, while 4 new missense changes (p.R681C, p.R681H, p.V683L, p.Y689C) affecting highly conserved codons were identified in 6 unrelated individuals, indicating that the prevalence of SNRNP200-associated adRP is relatively high. We also took advantage of this research to evaluate the pool-and-sequence method, especially with respect to the generation of false positive and negative results. We conclude that, although this strategy can be adopted for rapid discovery of new disease-associated variants, it still requires extensive validation to be used in routine DNA screenings. © 2011 Wiley-Liss, Inc.

Mutation identification by direct comparison of whole-genome sequencing data from mutant and wild-type individuals using k-mers.

Genes underlying mutant phenotypes can be isolated by combining marker discovery, genetic mapping and resequencing, but a more straightforward strategy for mapping mutations would be the direct comparison of mutant and wild-type genomes. Applying such an approach, however, is hampered by the need for reference sequences and by mutational loads that confound the unambiguous identification of causal mutations. Here we introduce NIKS (needle in the k-stack), a reference-free algorithm based on comparing k-mers in whole-genome sequencing data for precise discovery of homozygous mutations. We applied NIKS to eight mutants induced in nonreference rice cultivars and to two mutants of the nonmodel species Arabis alpina. In both species, comparing pooled F2 individuals selected for mutant phenotypes revealed small sets of mutations including the causal changes. Moreover, comparing M3 seedlings of two allelic mutants unambiguously identified the causal gene. Thus, for any species amenable to mutagenesis, NIKS enables forward genetics without requiring segregating populations, genetic maps and reference sequences.

Use of constant denaturant capillary electrophoresis of pooled blood samples to identify single-nucleotide polymorphisms in the genes (Scnn1a and Scnn1b) encoding the alpha and beta subunits of the epithelial sodium channel.

BACKGROUND: The epithelial sodium channel (ENaC) is composed of three homologous subunits: alpha, beta, and gamma. Mutations in the Scnn1b and Scnn1g genes, which encode the beta and the gamma subunits of ENaC, cause a severe form of hypertension (Liddle syndrome). The contribution of genetic variants within the Scnn1a gene, which codes for the alpha subunit, has not been investigated. METHODS: We screened for mutations in the COOH termini of the alpha and beta subunits of ENaC. Blood from 184 individuals from 31 families participating in a study on the genetics of hypertension were analyzed. Exons 13 of Scnn1a and Scnn1b, which encode the second transmembrane segment and the COOH termini of alpha- and beta-ENaC, respectively, were amplified from pooled DNA samples of members of each family by PCR. Constant denaturant capillary electrophoresis (CDCE) was used to detect mutations in PCR products of the pooled DNA samples. RESULTS: The detection limit of CDCE for ENaC variants was 1%, indicating that all members of any family or up to 100 individuals can be analyzed in one CDCE run. CDCE profiles of the COOH terminus of alpha-ENaC in pooled family members showed that the 31 families belonged to four groups and identified families with genetic variants. Using this approach, we analyzed 31 rather than 184 samples. Individual CDCE analysis of members from families with different pooled CDCE profiles revealed five genotypes containing 1853G-->T and 1987A-->G polymorphisms. The presence of the mutations was confirmed by DNA sequencing. For the COOH terminus of beta-ENaC, only one family showed a different CDCE profile. Two members of this family (n = 5) were heterozygous at 1781C-->T (T594M). CONCLUSION: CDCE rapidly detects point mutations in these candidate disease genes.

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.

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.

Coffee and tea intake and risk of head and neck cancer: pooled analysis in the International Head and Neck Cancer Epidemiology Consortium.

BACKGROUND: Only a few studies have explored the relation between coffee and tea intake and head and neck cancers, with inconsistent results. METHODS: We pooled individual-level data from nine case-control studies of head and neck cancers, including 5,139 cases and 9,028 controls. Logistic regression was used to estimate odds ratios (OR) and 95% confidence intervals (95% CI), adjusting for potential confounders. RESULTS: Caffeinated coffee intake was inversely related with the risk of cancer of the oral cavity and pharynx: the ORs were 0.96 (95% CI, 0.94-0.98) for an increment of 1 cup per day and 0.61 (95% CI, 0.47-0.80) in drinkers of >4 cups per day versus nondrinkers. This latter estimate was consistent for different anatomic sites (OR, 0.46; 95% CI, 0.30-0.71 for oral cavity; OR, 0.58; 95% CI, 0.41-0.82 for oropharynx/hypopharynx; and OR, 0.61; 95% CI, 0.37-1.01 for oral cavity/pharynx not otherwise specified) and across strata of selected covariates. No association of caffeinated coffee drinking was found with laryngeal cancer (OR, 0.96; 95% CI, 0.64-1.45 in drinkers of >4 cups per day versus nondrinkers). Data on decaffeinated coffee were too sparse for detailed analysis, but indicated no increased risk. Tea intake was not associated with head and neck cancer risk (OR, 0.99; 95% CI, 0.89-1.11 for drinkers versus nondrinkers). CONCLUSIONS: This pooled analysis of case-control studies supports the hypothesis of an inverse association between caffeinated coffee drinking and risk of cancer of the oral cavity and pharynx. IMPACT: Given widespread use of coffee and the relatively high incidence and low survival of head and neck cancers, the observed inverse association may have appreciable public health relevance.

An atlas of human gene expression from massively parallel signature sequencing (MPSS).

We have used massively parallel signature sequencing (MPSS) to sample the transcriptomes of 32 normal human tissues to an unprecedented depth, thus documenting the patterns of expression of almost 20,000 genes with high sensitivity and specificity. The data confirm the widely held belief that differences in gene expression between cell and tissue types are largely determined by transcripts derived from a limited number of tissue-specific genes, rather than by combinations of more promiscuously expressed genes. Expression of a little more than half of all known human genes seems to account for both the common requirements and the specific functions of the tissues sampled. A classification of tissues based on patterns of gene expression largely reproduces classifications based on anatomical and biochemical properties. The unbiased sampling of the human transcriptome achieved by MPSS supports the idea that most human genes have been mapped, if not functionally characterized. This data set should prove useful for the identification of tissue-specific genes, for the study of global changes induced by pathological conditions, and for the definition of a minimal set of genes necessary for basic cell maintenance. The data are available on the Web at http://mpss.licr.org and http://sgb.lynxgen.com.

Superior outcome of women with stage I/II cutaneous melanoma: pooled analysis of four European Organisation for Research and Treatment of Cancer phase III trials.

PURPOSE: Several studies observed a female advantage in the prognosis of cutaneous melanoma, for which behavioral factors or an underlying biologic mechanism might be responsible. Using complete and reliable follow-up data from four phase III trials of the European Organisation for Research and Treatment of Cancer (EORTC) Melanoma Group, we explored the female advantage across multiple end points and in relation to other important prognostic indicators. PATIENTS AND METHODS: Patients diagnosed with localized melanoma were included in EORTC adjuvant treatment trials 18832, 18871, 18952, and 18961 and randomly assigned during the period of 1984 to 2005. Cox proportional hazard models were used to calculate hazard ratios (HRs) and 95% CIs for women compared with men, adjusted for age, Breslow thickness, body site, ulceration, performed lymph node dissection, and treatment. RESULTS: A total of 2,672 patients with stage I/II melanoma were included. Women had a highly consistent and independent advantage in overall survival (adjusted HR, 0.70; 95% CI, 0.59 to 0.83), disease-specific survival (adjusted HR, 0.74; 95% CI, 0.62 to 0.88), time to lymph node metastasis (adjusted HR, 0.70; 95% CI, 0.51 to 0.96), and time to distant metastasis (adjusted HR, 0.69; 95% CI, 0.59 to 0.81). Subgroup analysis showed that the female advantage was consistent across all prognostic subgroups (with the possible exception of head and neck melanomas) and in pre- and postmenopausal age groups. CONCLUSION: Women have a consistent and independent relative advantage in all aspects of the progression of localized melanoma of approximately 30%, most likely caused by an underlying biologic sex difference.

Body mass index, cigarette smoking, and alcohol consumption and cancers of the oral cavity, pharynx, and larynx: modeling odds ratios in pooled case-control data.

Odds ratios for head and neck cancer increase with greater cigarette and alcohol use and lower body mass index (BMI; weight (kg)/height(2) (m(2))). Using data from the International Head and Neck Cancer Epidemiology Consortium, the authors conducted a formal analysis of BMI as a modifier of smoking- and alcohol-related effects. Analysis of never and current smokers included 6,333 cases, while analysis of never drinkers and consumers of < or =10 drinks/day included 8,452 cases. There were 8,000 or more controls, depending on the analysis. Odds ratios for all sites increased with lower BMI, greater smoking, and greater drinking. In polytomous regression, odds ratios for BMI (P = 0.65), smoking (P = 0.52), and drinking (P = 0.73) were homogeneous for oral cavity and pharyngeal cancers. Odds ratios for BMI and drinking were greater for oral cavity/pharyngeal cancer (P < 0.01), while smoking odds ratios were greater for laryngeal cancer (P < 0.01). Lower BMI enhanced smoking- and drinking-related odds ratios for oral cavity/pharyngeal cancer (P < 0.01), while BMI did not modify smoking and drinking odds ratios for laryngeal cancer. The increased odds ratios for all sites with low BMI may suggest related carcinogenic mechanisms; however, BMI modification of smoking and drinking odds ratios for cancer of the oral cavity/pharynx but not larynx cancer suggests additional factors specific to oral cavity/pharynx cancer.

Exome sequencing extends the phenotypic spectrum for ABHD12 mutations: from syndromic to nonsyndromic retinal degeneration.

OBJECTIVE: To identify the genetic causes underlying autosomal recessive retinitis pigmentosa (arRP) and to describe the associated phenotype. DESIGN: Case series. PARTICIPANTS: Three hundred forty-seven unrelated families affected by arRP and 33 unrelated families affected by retinitis pigmentosa (RP) plus noncongenital and progressive hearing loss, ataxia, or both, respectively. METHODS: A whole exome sequencing (WES) analysis was performed in 2 families segregating arRP. A mutational screening was performed in 378 additional unrelated families for the exon-intron boundaries of the ABHD12 gene. To establish a genotype-phenotype correlation, individuals who were homozygous or compound heterozygotes of mutations in ABHD12 underwent exhaustive clinical examinations by ophthalmologists, neurologists, and otologists. MAIN OUTCOME MEASURES: DNA sequence variants, best-corrected visual acuity, visual field assessments, electroretinogram responses, magnetic resonance imaging, and audiography. RESULTS: After a WES analysis, we identified 4 new mutations (p.Arg107Glufs*8, p.Trp159*, p.Arg186Pro, and p.Thr202Ile) in ABHD12 in 2 families (RP-1292 and W08-1833) previously diagnosed with nonsyndromic arRP, which cosegregated with the disease among the family members. Another homozygous mutation (p.His372Gln) was detected in 1 affected individual (RP-1487) from a cohort of 378 unrelated arRP and syndromic RP patients. After exhaustive clinical examinations by neurologists and otologists, the 4 affected members of the RP-1292 had no polyneuropathy or ataxia, and the sensorineural hearing loss and cataract were attributed to age or the normal course of the RP, whereas the affected members of the families W08-1833 and RP-1487 showed clearly symptoms associated with polyneuropathy, hearing loss, cerebellar ataxia, RP, and early-onset cataract (PHARC) syndrome. CONCLUSIONS: Null mutations in the ABHD12 gene lead to PHARC syndrome, a neurodegenerative disease including polyneuropathy, hearing loss, cerebellar ataxia, RP, and early-onset cataract. Our study allowed us to report 5 new mutations in ABHD12. This is the first time missense mutations have been described for this gene. Furthermore, these findings are expanding the spectrum of phenotypes associated with ABHD12 mutations ranging from PHARC syndrome to a nonsyndromic form of retinal degeneration.

Risk factors for head and neck cancer in young adults : a pooled analysis in the INHANCE consortium

BACKGROUND: Increasing incidence of head and neck cancer (HNC) in young adults has been reported. We aimed to compare the role of major risk factors and family history of cancer in HNC in young adults and older patients. METHODS: We pooled data from 25 case-control studies and conducted separate analyses for adults ≤45 years old ('young adults', 2010 cases and 4042 controls) and >45 years old ('older adults', 17 700 cases and 22 704 controls). Using logistic regression with studies treated as random effects, we estimated adjusted odds ratios (ORs) and 95% confidence intervals (CIs). RESULTS: The young group of cases had a higher proportion of oral tongue cancer (16.0% in women; 11.0% in men) and unspecified oral cavity / oropharynx cancer (16.2%; 11.1%) and a lower proportion of larynx cancer (12.1%; 16.6%) than older adult cases. The proportions of never smokers or never drinkers among female cases were higher than among male cases in both age groups. Positive associations with HNC and duration or pack-years of smoking and drinking were similar across age groups. However, the attributable fractions (AFs) for smoking and drinking were lower in young when compared with older adults (AFs for smoking in young women, older women, young men and older men, respectively, = 19.9% (95% CI = 9.8%, 27.9%), 48.9% (46.6%, 50.8%), 46.2% (38.5%, 52.5%), 64.3% (62.2%, 66.4%); AFs for drinking = 5.3% (-11.2%, 18.0%), 20.0% (14.5%, 25.0%), 21.5% (5.0%, 34.9%) and 50.4% (46.1%, 54.3%). A family history of early-onset cancer was associated with HNC risk in the young [OR = 2.27 (95% CI = 1.26, 4.10)], but not in the older adults [OR = 1.10 (0.91, 1.31)]. The attributable fraction for family history of early-onset cancer was 23.2% (8.60% to 31.4%) in young compared with 2.20% (-2.41%, 5.80%) in older adults. CONCLUSIONS: Differences in HNC aetiology according to age group may exist. The lower AF of cigarette smoking and alcohol drinking in young adults may be due to the reduced length of exposure due to the lower age. Other characteristics, such as those that are inherited, may play a more important role in HNC in young adults compared with older adults.

Global transcriptome sequencing identifies chlamydospore specific markers in Candida albicans and Candida dubliniensis.

Candida albicans and Candida dubliniensis are pathogenic fungi that are highly related but differ in virulence and in some phenotypic traits. During in vitro growth on certain nutrient-poor media, C. albicans and C. dubliniensis are the only yeast species which are able to produce chlamydospores, large thick-walled cells of unknown function. Interestingly, only C. dubliniensis forms pseudohyphae with abundant chlamydospores when grown on Staib medium, while C. albicans grows exclusively as a budding yeast. In order to further our understanding of chlamydospore development and assembly, we compared the global transcriptional profile of both species during growth in liquid Staib medium by RNA sequencing. We also included a C. albicans mutant in our study which lacks the morphogenetic transcriptional repressor Nrg1. This strain, which is characterized by its constitutive pseudohyphal growth, specifically produces masses of chlamydospores in Staib medium, similar to C. dubliniensis. This comparative approach identified a set of putatively chlamydospore-related genes. Two of the homologous C. albicans and C. dubliniensis genes (CSP1 and CSP2) which were most strongly upregulated during chlamydospore development were analysed in more detail. By use of the green fluorescent protein as a reporter, the encoded putative cell wall related proteins were found to exclusively localize to C. albicans and C. dubliniensis chlamydospores. Our findings uncover the first chlamydospore specific markers in Candida species and provide novel insights in the complex morphogenetic development of these important fungal pathogens.

Rapid bacterial genome sequencing: methods and applications in clinical microbiology.

The recent advances in sequencing technologies have given all microbiology laboratories access to whole genome sequencing. Providing that tools for the automated analysis of sequence data and databases for associated meta-data are developed, whole genome sequencing will become a routine tool for large clinical microbiology laboratories. Indeed, the continuing reduction in sequencing costs and the shortening of the 'time to result' makes it an attractive strategy in both research and diagnostics. Here, we review how high-throughput sequencing is revolutionizing clinical microbiology and the promise that it still holds. We discuss major applications, which include: (i) identification of target DNA sequences and antigens to rapidly develop diagnostic tools; (ii) precise strain identification for epidemiological typing and pathogen monitoring during outbreaks; and (iii) investigation of strain properties, such as the presence of antibiotic resistance or virulence factors. In addition, recent developments in comparative metagenomics and single-cell sequencing offer the prospect of a better understanding of complex microbial communities at the global and individual levels, providing a new perspective for understanding host-pathogen interactions. Being a high-resolution tool, high-throughput sequencing will increasingly influence diagnostics, epidemiology, risk management, and patient care.