SEQUENCE AND ANALYSIS OF THE

Phymatotrichopsis omnivora, a soilborne fungus confined to the southwestern regions of the United States, causes root rot in more than 2,000 species of dicotyledonous plants. To investigate the gene repertoire in this plant pathogenic fungus, a cDNA library containing expressed genes from the three distinct morphological stages in its life cycle and on exposure to three different nutrient conditions and a whole genome shotgun library was sequenced using massive parallel pyrosequencing and analyzed.Unique expressed sequence tags (ESTs) obtained by sequencing the ends of the cDNA transcripts from each library were examined for homologs in GenBank, KEGG, KOG and COGEME using the Blast alignment program and categorized into groups based on biological function assignment. Normalization and comparison of the metabolic profile of the ESTs from each library revealed stage specific gene expression.Analysis of the draft genomic sequence indicated that the ~74 Mbp assembly was approximately twice the size estimated by electrophoretic gel karyotyping from P. omnivora protoplasts supporting the hypothesis that this fungus is an obligate heterokaryon with several heterokaryotic nuclei. Approximately 22,000 genes were predicted with 9000 having homologs in GenBank and a further 12,000 sharing domains with proteins in the Pfam database. Annotation of the P. omnivora predicted proteins based on their biological function and metabolic pathways by comparison against the KOG, KEGG and COGEME databases revealed a comparable number of proteins involved in metabolism and cellular processes as found in the well studied filamentous fungi N. crassa and M. grisea. Moreover,P. omnivora was found to encode slightly higher numbers of ABC-type transporters and calcium and heavy metal transporting P-type ATPases than N. crassa and M. grisea. These ATP-dependent proteins likely are involved in the survival of the pathogen in calcareous heavy metal containing soils and on exposure to plant toxins and fungicides. One such protein of particular interest is the homolog of Bcmfs1, a multidrug transporter involved in protection against natural toxins and fungicides. Since P. omnivora also is a filamentous fungus that lacks both functional conidia and an active sexual cycle, this is the first study to provide details of the life style and metabolic profile of fungi with a parasexual cycle

Cloning, characterization and expression analysis of porcine microRNAs

Background: MicroRNAs (miRNAs) are small ~22-nt regulatory RNAs that can silence target genes, by blocking their protein production or degrading the mRNAs. Pig is an important animal in the agriculture industry because of its utility in the meat production. Besides, pig has tremendous biomedical importance as a model organism because of its closer proximity to humans than the mouse model. Several hundreds of miRNAs have been identified from mammals, humans, mice and rats, but little is known about the miRNA component in the pig genome. Here, we adopted an experimental approach to identify conserved and unique miRNAs and characterize their expression patterns in diverse tissues of pig.Results: By sequencing a small RNA library generated using pooled RNA from the pig heart, liver and thymus; we identified a total of 120 conserved miRNA homologs in pig. Expression analysis of conserved miRNAs in 14 different tissue types revealed heart-specific expression of miR-499 and miR-208 and liver-specific expression of miR-122. Additionally, miR-1 and miR-133 in the heart, miR-181a and miR-142-3p in the thymus, miR-194 in the liver, and miR-143 in the stomach showed the highest levels of expression. miR-22, miR-26b, miR-29c and miR-30c showed ubiquitous expression in diverse tissues. The expression patterns of pig-specific miRNAs also varied among the tissues examined.Conclusion: Identification of 120 miRNAs and determination of the spatial expression patterns of a sub-set of these in the pig is a valuable resource for molecular biologists, breeders, and biomedical investigators interested in post-transcriptional gene regulation in pig and in related mammals, including humans.Peer reviewedBiochemistry and Molecular BiologyAnimal Scienc

Microcollinearity between autopolyploid sugarcane and diploid sorghum genomes

Abstract\ud \ud \ud \ud Background\ud \ud Sugarcane (Saccharum spp.) has become an increasingly important crop for its leading role in biofuel production. The high sugar content species S. officinarum is an octoploid without known diploid or tetraploid progenitors. Commercial sugarcane cultivars are hybrids between S. officinarum and wild species S. spontaneum with ploidy at ~12×. The complex autopolyploid sugarcane genome has not been characterized at the DNA sequence level.\ud \ud \ud \ud Results\ud \ud The microsynteny between sugarcane and sorghum was assessed by comparing 454 pyrosequences of 20 sugarcane bacterial artificial chromosomes (BACs) with sorghum sequences. These 20 BACs were selected by hybridization of 1961 single copy sorghum overgo probes to the sugarcane BAC library with one sugarcane BAC corresponding to each of the 20 sorghum chromosome arms. The genic regions of the sugarcane BACs shared an average of 95.2% sequence identity with sorghum, and the sorghum genome was used as a template to order sequence contigs covering 78.2% of the 20 BAC sequences. About 53.1% of the sugarcane BAC sequences are aligned with sorghum sequence. The unaligned regions contain non-coding and repetitive sequences. Within the aligned sequences, 209 genes were annotated in sugarcane and 202 in sorghum. Seventeen genes appeared to be sugarcane-specific and all validated by sugarcane ESTs, while 12 appeared sorghum-specific but only one validated by sorghum ESTs. Twelve of the 17 sugarcane-specific genes have no match in the non-redundant protein database in GenBank, perhaps encoding proteins for sugarcane-specific processes. The sorghum orthologous regions appeared to have expanded relative to sugarcane, mostly by the increase of retrotransposons.\ud \ud \ud \ud Conclusions\ud \ud The sugarcane and sorghum genomes are mostly collinear in the genic regions, and the sorghum genome can be used as a template for assembling much of the genic DNA of the autopolyploid sugarcane genome. The comparable gene density between sugarcane BACs and corresponding sorghum sequences defied the notion that polyploidy species might have faster pace of gene loss due to the redundancy of multiple alleles at each locus.We acknowledge our colleagues at the University of Oklahomas Advanced Center for Genome Technology, Chunmei Qu and Ping Wang for their assistance with 454 GSFLX sequencing sample preparation and Steve Kenton for his help with deconvoluting the pooled BACs and their subsequent assembly. We also thank Eric Tang for assistance on sequencing two BACs using Sanger sequencers. This project is supported by startup funds from the University of Illinois to RM and a grant from the Energy Bioscience Institute (EBI) to SPM, MEH, RM, and DSR.We acknowledge our colleagues at the University of Oklahoma's Advanced Center for Genome Technology, Chunmei Qu and Ping Wang for their assistance with 454 GS-FLX sequencing sample preparation and Steve Kenton for his help with deconvoluting the pooled BACs and their subsequent assembly. We also thank Eric Tang for assistance on sequencing two BACs using Sanger sequencers. This project is supported by start-up funds from the University of Illinois to RM and a grant from the Energy Bioscience Institute (EBI) to SPM, MEH, RM, and DSR

Improved genome assembly and evidence-based global gene model set for the chordate Ciona intestinalis: new insight into intron and operon populations

An improved assembly of the Ciona intestinalis genome reveals that it contains non-canonical introns and that about 20% of Ciona genes reside in operons

Genome-Wide DNA Methylation Maps in Follicular Lymphoma Cells Determined by Methylation-Enriched Bisulfite Sequencing

BACKGROUND: Follicular lymphoma (FL) is a form of non-Hodgkin's lymphoma (NHL) that arises from germinal center (GC) B-cells. Despite the significant advances in immunotherapy, FL is still not curable. Beyond transcriptional profiling and genomics datasets, there currently is no epigenome-scale dataset or integrative biology approach that can adequately model this disease and therefore identify novel mechanisms and targets for successful prevention and treatment of FL. METHODOLOGY/PRINCIPAL FINDINGS: We performed methylation-enriched genome-wide bisulfite sequencing of FL cells and normal CD19(+) B-cells using 454 sequencing technology. The methylated DNA fragments were enriched with methyl-binding proteins, treated with bisulfite, and sequenced using the Roche-454 GS FLX sequencer. The total number of bases covered in the human genome was 18.2 and 49.3 million including 726,003 and 1.3 million CpGs in FL and CD19(+) B-cells, respectively. 11,971 and 7,882 methylated regions of interest (MRIs) were identified respectively. The genome-wide distribution of these MRIs displayed significant differences between FL and normal B-cells. A reverse trend in the distribution of MRIs between the promoter and the gene body was observed in FL and CD19(+) B-cells. The MRIs identified in FL cells also correlated well with transcriptomic data and ChIP-on-Chip analyses of genome-wide histone modifications such as tri-methyl-H3K27, and tri-methyl-H3K4, indicating a concerted epigenetic alteration in FL cells. CONCLUSIONS/SIGNIFICANCE: This study is the first to provide a large scale and comprehensive analysis of the DNA methylation sequence composition and distribution in the FL epigenome. These integrated approaches have led to the discovery of novel and frequent targets of aberrant epigenetic alterations. The genome-wide bisulfite sequencing approach developed here can be a useful tool for profiling DNA methylation in clinical samples

Phylotyping and Functional Analysis of Two Ancient Human Microbiomes

Background: The Human Microbiome Project (HMP) is one of the U.S. National Institutes of Health Roadmap for Medical Research. Primary interests of the HMP include the distinctiveness of different gut microbiomes, the factors influencing microbiome diversity, and the functional redundancies of the members of human microbiotas. In this present work, we contribute to these interests by characterizing two extinct human microbiotas. Methodology/Principal Findings: We examine two paleofecal samples originating from cave deposits in Durango Mexico and dating to approximately 1300 years ago. Contamination control is a serious issue in ancient DNA research; we use a novel approach to control contamination. After we determined that each sample originated from a different human, we generated 45 thousand shotgun DNA sequencing reads. The phylotyping and functional analysis of these reads reveals a signature consistent with the modern gut ecology. Interestingly, inter-individual variability for phenotypes but not functional pathways was observed. The two ancient samples have more similar functional profiles to each other than to a recently published profile for modern humans. This similarity could not be explained by a chance sampling of the databases. Conclusions/Significance: We conduct a phylotyping and functional analysis of ancient human microbiomes, while providing novel methods to control for DNA contamination and novel hypotheses about past microbiome biogeography. We postulate that natural selection has more of an influence on microbiome functional profiles than it does on the species represented in the microbial ecology. We propose that human microbiomes were more geographically structured during pre-Columbian times than today

BAC-pool sequencing and analysis of large segments of A12 and D12 homoeologous chromosomes in upland cotton.

Acknowledgments “Dedicated to Dr. Ramesh Kantety, a mentor, colleague and friend”. We would like to acknowledge the support offered by Padmini Sripathi during data analysis and submissions. Author Contributions Conceived and designed the experiments: RVK JZY. Performed the experiments: RB ZX SM GBW. Analyzed the data: RB. Contributed reagents/materials/analysis tools: RVK RB JZY RJK BAR. Wrote the manuscript: RB. Revised the manuscript: RB RVK JZY RGP BAR GCS. Advised the research: RVK JZY RGP BAR GCS.Author Contributions Conceived and designed the experiments: RVK JZY. Performed the experiments: RB ZX SM GBW. Analyzed the data: RB. Contributed reagents/materials/analysis tools: RVK RB JZY RJK BAR. Wrote the manuscript: RB. Revised the manuscript: RB RVK JZY RGP BAR GCS. Advised the research: RVK JZY RGP BAR GCS.Although new and emerging next-generation sequencing (NGS) technologies have reduced sequencing costs significantly, much work remains to implement them for de novo sequencing of complex and highly repetitive genomes such as the tetraploid genome of Upland cotton (Gossypium hirsutum L.). Herein we report the results from implementing a novel, hybrid Sanger/454-based BAC-pool sequencing strategy using minimum tiling path (MTP) BACs from Ctg-3301 and Ctg-465, two large genomic segments in A12 and D12 homoeologous chromosomes (Ctg). To enable generation of longer contig sequences in assembly, we implemented a hybrid assembly method to process ~35x data from 454 technology and 2.8-3x data from Sanger method. Hybrid assemblies offered higher sequence coverage and better sequence assemblies. Homology studies revealed the presence of retrotransposon regions like Copia and Gypsy elements in these contigs and also helped in identifying new genomic SSRs. Unigenes were anchored to the sequences in Ctg-3301 and Ctg-465 to support the physical map. Gene density, gene structure and protein sequence information derived from protein prediction programs were used to obtain the functional annotation of these genes. Comparative analysis of both contigs with Arabidopsis genome exhibited synteny and microcollinearity with a conserved gene order in both genomes. This study provides insight about use of MTP-based BAC-pool sequencing approach for sequencing complex polyploid genomes with limited constraints in generating better sequence assemblies to build reference scaffold sequences. Combining the utilities of MTP-based BAC-pool sequencing with current longer and short read NGS technologies in multiplexed format would provide a new direction to cost-effectively and precisely sequence complex plant genomes.Yeshttp://www.plosone.org/static/editorial#pee

Microarray analysis of femaleand larval-specific gene expression in the horn fly (Diptera: Muscidae)

The horn My, Haematobia irritans L., is an obligate blood-feeding parasite of cattle, and control of this pest is a continuing problem because the fly is becoming resistant to pesticides. Dominant conditional lethal gene systems are being studied as population control technologies against agricultural pests. One of the components of these systems is a female-specific gene promoter that drives expression of a lethality-inducing gene. To identify candidate genes to supply this promoter, microarrays were designed from a horn fly expressed sequence tag (EST) database and probed to identify female-specific and larval-specific gene expression. Analysis of dye swap experiments found 432 and 417 transcripts whose expression levels were higher or lower in adult female flies, respectively, compared with adult male flies. Additionally, 419 and 871 transcripts were identified whose expression levels were higher or lower in first-instar larvae compared with adult flies, respectively. Three transcripts were expressed more highly in adult females flies compared with adult males and also higher in the first-instar larval lifestage compared with adult flies. One of these transcripts, a putative nanos ortholog, has a high female-to-male expression ratio, a moderate expression level in first-instar larvae, and has been well characterized in Drosophila. melanogaster (Meigen). In conclusion, we used microarray technology, verified by reverse transcriptase-polymerase chain reaction and massively parallel pyrosequencing, to study life stageand sex-specific gene expression in the horn fly and identified three gene candidates for detailed evaluation as a gene promoter source for the development of a female-specific conditional lethality system

High-throughput sequence analysis of Ciona intestinalis SL trans-spliced mRNAs: Alternative expression modes and gene function correlates

Pre-mRNA 5′ spliced-leader (SL) trans-splicing occurs in some metazoan groups but not in others. Genome-wide characterization of the trans-spliced mRNA subpopulation has not yet been reported for any metazoan. We carried out a high-throughput analysis of the SL trans-spliced mRNA population of the ascidian tunicate Ciona intestinalis by 454 Life Sciences (Roche) pyrosequencing of SL-PCR-amplified random-primed reverse transcripts of tailbud embryo RNA. We obtained ∼250,000 high-quality reads corresponding to 8790 genes, ∼58% of the Ciona total gene number. The great depth of this data revealed new aspects of trans-splicing, including the existence of a significant class of “infrequently trans-spliced” genes, accounting for ∼28% of represented genes, that generate largely non-trans-spliced mRNAs, but also produce trans-spliced mRNAs, in part through alternative promoter use. Thus, the conventional qualitative dichotomy of trans-spliced versus non-trans-spliced genes should be supplanted by a more accurate quantitative view recognizing frequently and infrequently trans-spliced gene categories. Our data include reads representing ∼80% of Ciona frequently trans-spliced genes. Our analysis also revealed significant use of closely spaced alternative trans-splice acceptor sites which further underscores the mechanistic similarity of cis- and trans-splicing and indicates that the prevalence of ±3-nt alternative splicing events at tandem acceptor sites, NAGNAG, is driven by spliceosomal mechanisms, and not nonsense-mediated decay, or selection at the protein level. The breadth of gene representation data enabled us to find new correlations between trans-splicing status and gene function, namely the overrepresentation in the frequently trans-spliced gene class of genes associated with plasma/endomembrane system, Ca2+ homeostasis, and actin cytoskeleton