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

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.

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.

Microcollinearity between autopolyploid sugarcane and diploid sorghum genomes.

BACKGROUND: 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 approximately 12x. The complex autopolyploid sugarcane genome has not been characterized at the DNA sequence level. RESULTS: 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. CONCLUSIONS: 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.

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 approximately 250,000 high-quality reads corresponding to 8790 genes, approximately 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 approximately 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 approximately 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, Ca(2+) homeostasis, and actin cytoskeleton.

Cloning and characterization of small RNAs from Medicago truncatula reveals four novel legume-specific microRNA families.

MicroRNAs (miRNAs) and small-interfering RNAs (siRNAs) have emerged as important regulators of gene expression in higher eukaryotes. Recent studies indicate that genomes in higher plants encode lineage-specific and species-specific miRNAs in addition to the well-conserved miRNAs. Leguminous plants are grown throughout the world for food and forage production. To date the lack of genomic sequence data has prevented systematic examination of small RNAs in leguminous plants. Medicago truncatula, a diploid plant with a near-completely sequenced genome has recently emerged as an important model legume. We sequenced a small RNA library generated from M. truncatula to identify not only conserved miRNAs but also novel small RNAs, if any. Eight novel small RNAs were identified, of which four (miR1507, miR2118, miR2119 and miR2199) are annotated as legume-specific miRNAs because these are conserved in related legumes. Three novel transcripts encoding TIR-NBS-LRR proteins are validated as targets for one of the novel miRNA, miR2118. Small RNA sequence analysis coupled with the small RNA blot analysis, confirmed the expression of around 20 conserved miRNA families in M. truncatula. Fifteen transcripts have been validated as targets for conserved miRNAs. We also characterized Tas3-siRNA biogenesis in M. truncatula and validated three auxin response factor (ARF) transcripts that are targeted by tasiRNAs. These findings indicate that M. truncatula and possibly other related legumes have complex mechanisms of gene regulation involving specific and common small RNAs operating post-transcriptionally.

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

The horn fly, 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 stage- and 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.

Methods for generating shotgun and mixed shotgun/paired-end libraries for the 454 DNA sequencer.

With the introduction of massively parallel, microminiature-based instrumentation for DNA sequencing, robust, reproducible, optimized methods are needed to prepare the target DNA for analysis using these high-throughput approaches because the cost per instrument run is orders of magnitude more than for typical Sanger dideoxynucleotide sequencing on fluorescence-based capillary systems. The methods provided by the manufacturer for genome sequencing using the 454/Roche GS-20 and GS-FLX instruments are robust. However, in an effort to streamline them for automation, we have incorporated several novel changes and deleted several extraneous steps. As a result of modifying these sample preparation protocols, the number of manual manipulations has also been minimized, and the overall yields have been improved for both shotgun and mixed shotgun/paired-end libraries.

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.

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.

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

BACKGROUND: The draft genome sequence of the ascidian Ciona intestinalis, along with associated gene models, has been a valuable research resource. However, recently accumulated expressed sequence tag (EST)/cDNA data have revealed numerous inconsistencies with the gene models due in part to intrinsic limitations in gene prediction programs and in part to the fragmented nature of the assembly. RESULTS: We have prepared a less-fragmented assembly on the basis of scaffold-joining guided by paired-end EST and bacterial artificial chromosome (BAC) sequences, and BAC chromosomal in situ hybridization data. The new assembly (115.2 Mb) is similar in length to the initial assembly (116.7 Mb) but contains 1,272 (approximately 50%) fewer scaffolds. The largest scaffold in the new assembly incorporates 95 initial-assembly scaffolds. In conjunction with the new assembly, we have prepared a greatly improved global gene model set strictly correlated with the extensive currently available EST data. The total gene number (15,254) is similar to that of the initial set (15,582), but the new set includes 3,330 models at genomic sites where none were present in the initial set, and 1,779 models that represent fusions of multiple previously incomplete models. In approximately half, 5'-ends were precisely mapped using 5'-full-length ESTs, an important refinement even in otherwise unchanged models. CONCLUSION: Using these new resources, we identify a population of non-canonical (non-GT-AG) introns and also find that approximately 20% of Ciona genes reside in operons and that operons contain a high proportion of single-exon genes. Thus, the present dataset provides an opportunity to analyze the Ciona genome much more precisely than ever.

Purification and characterization of chitinase from Alcaligenes xylosoxydans.

Extracellular chitinase from Alcaligenes xylosoxydans was purified to electrophoretic homogeneity using affinity and gel filtration chromatography. The molecular mass of chitinase was estimated to be 45 kDa and 44 kDa by SDS-PAGE and gel-filtration, respectively. The enzyme was optimally active at 50 degrees C (over 30 min) and pH 5. Activity staining after PAGE showed a single band. The Km for chitin was 3 g l-1. Cu2+ and Na+ at 5 mM inhibited chitinase activity to 25% while Ca2+, Mg2+ and Ba2+ had no effect at the same concentration. The purified enzyme degraded mycelia of Aspergillus niger.