Yunchen Gong

Next-generation mapping of Arabidopsis genes

Next-generation genomic sequencing technologies have made it possible to directly map mutations responsible for phenotypes of interest via direct sequencing. However, most mapping strategies proposed to date require some prior genetic analysis, which can be very time-consuming even in genetically tractable organisms. Here we present a de novo method for rapidly and robustly mapping the physical location of EMS mutations by sequencing a small pooled F₂ population. This method, called Next Generation Mapping (NGM), uses a chastity statistic to quantify the relative contribution of the parental mutant and mapping lines to each SNP in the pooled F₂ population. It then uses this information to objectively localize the candidate mutation based on its exclusive segregation with the mutant parental line. A user-friendly, web-based tool for performing NGM analysis is available at http://bar.utoronto.ca/NGM. We used NGM to identify three genes involved in cell-wall biology in Arabidopsis thaliana, and, in a power analysis, demonstrate success in test mappings using as few as ten F₂ lines and a single channel of Illumina Genome Analyzer data. This strategy can easily be applied to other model organisms, and we expect that it will also have utility in crops and any other eukaryote with a completed genome sequence.

Lung microbiota across age and disease stage in cystic fibrosis

Understanding the significance of bacterial species that colonize and persist in cystic fibrosis (CF) airways requires a detailed examination of bacterial community structure across a broad range of age and disease stage. We used 16S ribosomal RNA sequencing to characterize the lung microbiota in 269 CF patients spanning a 60 year age range, including 76 pediatric samples from patients of age 4–17, and a broad cross-section of disease status to identify features of bacterial community structure and their relationship to disease stage and age. The CF lung microbiota shows significant inter-individual variability in community structure, composition and diversity. The core microbiota consists of five genera - Streptococcus, Prevotella, Rothia, Veillonella and Actinomyces. CF-associated pathogens such as Pseudomonas, Burkholderia, Stenotrophomonas and Achromobacter are less prevalent than core genera, but have a strong tendency to dominate the bacterial community when present. Community diversity and lung function are greatest in patients less than 10 years of age and lower in older age groups, plateauing at approximately age 25. Lower community diversity correlates with worse lung function in a multivariate regression model. Infection by Pseudomonas correlates with age-associated trends in community diversity and lung function.

Transcriptomic Responses of the Softwood-Degrading White-Rot Fungus Phanerochaete carnosa during Growth on Coniferous and Deciduous Wood

ABSTRACT To identify enzymes that could be developed to reduce the recalcitrance of softwood resources, the transcriptomes of the softwood-degrading white-rot fungus Phanerochaete carnosa were evaluated after growth on lodgepole pine, white spruce, balsam fir, and sugar maple and compared to the transcriptome of P. carnosa after growth on liquid nutrient medium. One hundred fifty-two million paired-end reads were obtained, and 63% of these reads were mapped to 10,257 gene models from P. carnosa. Five-hundred thirty-three of these genes had transcripts that were at least four times more abundant during growth on at least one wood medium than on nutrient medium. The 30 transcripts that were on average over 100 times more abundant during growth on wood than on nutrient medium included 6 manganese peroxidases, 5 cellulases, 2 hemicellulases, a lignin peroxidase, glyoxal oxidase, and a P450 monooxygenase. Notably, among the genes encoding putative cellulases, one encoding a glycosyl hydrolase family 61 protein had the highest relative transcript abundance during growth on wood. Overall, transcripts predicted to encode lignin-degrading activities were more abundant than those predicted to encode carbohydrate-active enzymes. Transcripts predicted to encode three MnPs represented the most highly abundant transcripts in wood-grown cultivations compared to nutrient medium cultivations. Gene set enrichment analyses did not distinguish transcriptomes resulting from softwood and hardwood cultivations, suggesting that similar sets of enzyme activities are elicited by P. carnosa grown on different wood substrates, albeit to different expression levels.

Comparative genomics of the white-rot fungi, Phanerochaete carnosa and P. chrysosporium, to elucidate the genetic basis of the distinct wood types they colonize

BackgroundSoftwood is the predominant form of land plant biomass in the Northern hemisphere, and is among the most recalcitrant biomass resources to bioprocess technologies. The white rot fungus, Phanerochaete carnosa, has been isolated almost exclusively from softwoods, while most other known white-rot species, including Phanerochaete chrysosporium, were mainly isolated from hardwoods. Accordingly, it is anticipated that P. carnosa encodes a distinct set of enzymes and proteins that promote softwood decomposition. To elucidate the genetic basis of softwood bioconversion by a white-rot fungus, the present study reports the P. carnosa genome sequence and its comparative analysis with the previously reported P. chrysosporium genome.ResultsP. carnosa encodes a complete set of lignocellulose-active enzymes. Comparative genomic analysis revealed that P. carnosa is enriched with genes encoding manganese peroxidase, and that the most divergent glycoside hydrolase families were predicted to encode hemicellulases and glycoprotein degrading enzymes. Most remarkably, P. carnosa possesses one of the largest P450 contingents (266 P450s) among the sequenced and annotated wood-rotting basidiomycetes, nearly double that of P. chrysosporium. Along with metabolic pathway modeling, comparative growth studies on model compounds and chemical analyses of decomposed wood components showed greater tolerance of P. carnosa to various substrates including coniferous heartwood.ConclusionsThe P. carnosa genome is enriched with genes that encode P450 monooxygenases that can participate in extractives degradation, and manganese peroxidases involved in lignin degradation. The significant expansion of P450s in P. carnosa, along with differences in carbohydrate- and lignin-degrading enzymes, could be correlated to the utilization of heartwood and sapwood preparations from both coniferous and hardwood species.

Tri6 Is a Global Transcription Regulator in the Phytopathogen Fusarium graminearum

In F. graminearum, the transcriptional regulator Tri6 is encoded within the trichothecene gene cluster and regulates genes involved in the biosynthesis of the secondary metabolite deoxynivalenol (DON). The Tri6 protein with its Cys2His2 zinc-finger may also conform to the class of global transcription regulators. This class of global transcriptional regulators mediate various environmental cues and generally responds to the demands of cellular metabolism. To address this issue directly, we sought to find gene targets of Tri6 in F. graminearum grown in optimal nutrient conditions. Chromatin immunoprecipitation followed by Illumina sequencing (ChIP-Seq) revealed that in addition to identifying six genes within the trichothecene gene cluster, Tri1, Tri3, Tri6, Tri7, Tri12 and Tri14, the ChIP-Seq also identified 192 additional targets potentially regulated by Tri6. Functional classification revealed that, among the annotated genes, ∼40% are associated with cellular metabolism and transport and the rest of the target genes fall into the category of signal transduction and gene expression regulation. ChIP-Seq data also revealed Tri6 has the highest affinity toward its own promoter, suggesting that this gene could be subject to self-regulation. Electro mobility shift assays (EMSA) performed on the promoter of Tri6 with purified Tri6 protein identified a minimum binding motif of GTGA repeats as a consensus sequence. Finally, expression profiling of F. graminearum grown under nitrogen-limiting conditions revealed that 49 out of 198 target genes are differentially regulated by Tri6. The identification of potential new targets together with deciphering novel binding sites for Tri6, casts new light into the role of this transcriptional regulator in the overall growth and development of F. graminearum.

Quantitative Interactor Screening with next- generation Sequencing (QIS-Seq) identifies Arabidopsis thaliana MLO2 as a target of the Pseudomonas syringae type III effector HopZ2

BackgroundIdentification of protein-protein interactions is a fundamental aspect of understanding protein function. A commonly used method for identifying protein interactions is the yeast two-hybrid system.ResultsHere we describe the application of next-generation sequencing to yeast two-hybrid interaction screens and develop Quantitative Interactor Screen Sequencing (QIS-Seq). QIS-Seq provides a quantitative measurement of enrichment for each interactor relative to its frequency in the library as well as its general stickiness (non-specific binding). The QIS-Seq approach is scalable and can be used with any yeast two-hybrid screen and with any next-generation sequencing platform. The quantitative nature of QIS-Seq data make it amenable to statistical evaluation, and importantly, facilitates the standardization of experimental design, data collection, and data analysis. We applied QIS-Seq to identify the Arabidopsis thaliana MLO2 protein as a target of the Pseudomonas syringae type III secreted effector protein HopZ2. We validate the interaction between HopZ2 and MLO2 in planta and show that the interaction is required for HopZ2-associated virulence.ConclusionsWe demonstrate that QIS-Seq is a high-throughput quantitative interactor screen and validate MLO2 as an interactor and novel virulence target of the P. syringae type III secreted effector HopZ2.

Metatranscriptome of an Anaerobic Benzene-Degrading, Nitrate-Reducing Enrichment Culture Reveals Involvement of Carboxylation in Benzene Ring Activation

ABSTRACT The enzymes involved in the initial steps of anaerobic benzene catabolism are not known. To try to elucidate this critical step, a metatranscriptomic analysis was conducted to compare the genes transcribed during the metabolism of benzene and benzoate by an anaerobic benzene-degrading, nitrate-reducing enrichment culture. RNA was extracted from the mixed culture and sequenced without prior mRNA enrichment, allowing simultaneous examination of the active community composition and the differential gene expression between the two treatments. Ribosomal and mRNA sequences attributed to a member of the family Peptococcaceae from the order Clostridiales were essentially only detected in the benzene-amended culture samples, implicating this group in the initial catabolism of benzene. Genes similar to each of two subunits of a proposed benzene-carboxylating enzyme were transcribed when the culture was amended with benzene. Anaerobic benzoate degradation genes from strict anaerobes were transcribed only when the culture was amended with benzene. Genes for other benzoate catabolic enzymes and for nitrate respiration were transcribed in both samples, with those attributed to an Azoarcus species being most abundant. These findings indicate that the mineralization of benzene starts with its activation by a strict anaerobe belonging to the Peptococcaceae, involving a carboxylation step to form benzoate. These data confirm the previously hypothesized syntrophic association between a benzene-degrading Peptococcaceae strain and a benzoate-degrading denitrifying Azoarcus strain for the complete catabolism of benzene with nitrate as the terminal electron acceptor.

Analysis of the Cystic Fibrosis Lung Microbiota via Serial Illumina Sequencing of Bacterial 16S rRNA Hypervariable Regions

The characterization of bacterial communities using DNA sequencing has revolutionized our ability to study microbes in nature and discover the ways in which microbial communities affect ecosystem functioning and human health. Here we describe Serial Illumina Sequencing (SI-Seq): a method for deep sequencing of the bacterial 16S rRNA gene using next-generation sequencing technology. SI-Seq serially sequences portions of the V5, V6 and V7 hypervariable regions from barcoded 16S rRNA amplicons using an Illumina short-read genome analyzer. SI-Seq obtains taxonomic resolution similar to 454 pyrosequencing for a fraction of the cost, and can produce hundreds of thousands of reads per sample even with very high multiplexing. We validated SI-Seq using single species and mock community controls, and via a comparison to cystic fibrosis lung microbiota sequenced using 454 FLX Titanium. Our control runs show that SI-Seq has a dynamic range of at least five orders of magnitude, can classify >96% of sequences to the genus level, and performs just as well as 454 and paired-end Illumina methods in estimation of standard microbial ecology diversity measurements. We illustrate the utility of SI-Seq in a pilot sample of central airway secretion samples from cystic fibrosis patients.

Extensive remodeling of the Pseudomonas syringae pv. avellanae type III secretome associated with two independent host shifts onto hazelnut

BackgroundHazelnut (Corylus avellana) decline disease in Greece and Italy is caused by the convergent evolution of two distantly related lineages of Pseudomonas syringae pv. avellanae (Pav). We sequenced the genomes of three Pav isolates to determine if their convergent virulence phenotype had a common genetic basis due to either genetic exchange between lineages or parallel evolution.ResultsWe found little evidence for horizontal transfer (recombination) of genes between Pav lineages, but two large genomic islands (GIs) have been recently acquired by one of the lineages. Evolutionary analyses of the genes encoding type III secreted effectors (T3SEs) that are translocated into host cells and are important for both suppressing and eliciting defense responses show that the two Pav lineages have dramatically different T3SE profiles, with only two shared putatively functional T3SEs. One Pav lineage has undergone unprecedented secretome remodeling, including the acquisition of eleven new T3SEs and the loss or pseudogenization of 15, including five of the six core T3SE families that are present in the other Pav lineage. Molecular dating indicates that divergence within both of the Pav lineages predates their observation in the field. This suggest that both Pav lineages have been cryptically infecting hazelnut trees or wild relatives for many years, and that the emergence of hazelnut decline in the 1970s may have been due to changes in agricultural practice.ConclusionsThese data show that divergent lineages of P. syringae can converge on identical disease etiology on the same host plant using different virulence mechanisms and that dramatic shifts in the arsenal of T3SEs can accompany disease emergence.

Semi-Automatic In Silico Gap Closure Enabled De Novo Assembly of Two Dehalobacter Genomes from Metagenomic Data

Typically, the assembly and closure of a complete bacterial genome requires substantial additional effort spent in a wet lab for gap resolution and genome polishing. Assembly is further confounded by subspecies polymorphism when starting from metagenome sequence data. In this paper, we describe an in silico gap-resolution strategy that can substantially improve assembly. This strategy resolves assembly gaps in scaffolds using pre-assembled contigs, followed by verification with read mapping. It is capable of resolving assembly gaps caused by repetitive elements and subspecies polymorphisms. Using this strategy, we realized the de novo assembly of the first two Dehalobacter genomes from the metagenomes of two anaerobic mixed microbial cultures capable of reductive dechlorination of chlorinated ethanes and chloroform. Only four additional PCR reactions were required even though the initial assembly with Newbler v. 2.5 produced 101 contigs within 9 scaffolds belonging to two Dehalobacter strains. By applying this strategy to the re-assembly of a recently published genome of Bacteroides, we demonstrate its potential utility for other sequencing projects, both metagenomic and genomic.