Shotgun library construction for DNA sequencing

Mitochondrial genome analysis, phylogeny and divergence time evaluation of Strixaluco (Aves, Strigiformes, Strigidae)

Background Prior research has shown that the European peninsulas were the main sources of Strixaluco colonisation of Northern Europe during the late glacial period. However, the phylogenetic relationship and the divergence time between S.aluco from Leigong Mountain Nature Reserve, Guizhou Province, China and the Strigiformes from overseas remains unclear. The mitochondrial genome structure of birds is a covalent double-chain loop structure that is highly conserved and, thus, suitable for phylogenetic analysis. This study examined the phylogenetic relationship and divergence time of Strix using the whole mitochondrial genome of S.aluco. New information In this study, the complete mitochondrial genome of Strixaluco, with a total length of 18,632 bp, is reported for the first time. A total of 37 genes were found, including 22 tRNAs, two rRNAs, 13 protein-coding genes and two non-coding control regions. Certain species of Tytoninae were used as out-group and PhyloSuite software was applied to build the ML-tree and BI-tree of Strigiformes. Finally, the divergence time tree was constructed using BEAST 2.6.7 software and the age of Miosurniadiurna fossil-bearing sediments (6.0-9.5 Ma) was set as internal correction point. The common ancestor of Strix was confirmed to have diverged during the Pleistocene (2.58-0.01 Ma). The combined action of the dramatic uplift of the Qinling Mountains in the Middle Pleistocene and the climate oscillation of the Pleistocene caused Strix divergence between the northern and southern parts of mainland China. The isolation of glacial-interglacial rotation and glacier refuge was the main reason for the divergence of Strixuralensis and S.aluco from their common ancestor during this period. This study provides a reference for the evolutionary history of S.aluco.

Polydimethylsiloxane microstructure-induced acoustic streaming for enhanced ultrasonic DNA fragmentation on a microfluidic chip

Next-generation sequencing (NGS) is an essential technology for DNA identification in genomic research. DNA fragmentation is a critical step for NGS and doing this on-chip is of great interest for future integrated genomic solutions. Here we demonstrate fast acoustofluidic DNA fragmentation via ultrasound-actuated elastic polydimethylsiloxane (PDMS) microstructures that induce acoustic streaming and associated shear forces when placed in the field of an ultrasonic transducer. Indeed, acoustic streaming locally generates high tensile stresses that can mechanically stretch and break DNA molecule chains. The improvement in efficiency of the on-chip DNA fragmentation is due to the synergistic effect of these tensile stresses and ultrasonic cavitation phenomena. We tested these microstructure-induced effects in a DNA-containing microfluidic channel both experimentally and by simulation. The DNA fragmentation process was evaluated by measuring the change in the DNA fragment size over time. The chip works well with both long and short DNA chains; in particular, purified lambda (λ) DNA was cut from 48.5 kbp to 3 kbp in one minute with selected microstructures and further down to 300 bp within two and a half minutes. The fragment size of mouse genomic DNA was reduced from 1.4 kbp to 400 bp in one minute and then to 200 bp in two and a half minutes. The DNA fragmentation efficiency of the chip equipped with the PDMS microstructures was twice that of the chip without the microstructures. Exhaustive comparison shows that the on-chip fragmentation performance reaches the level of high-end professional standards. Recently, DNA fragmentation was shown to be enhanced using vibrating air microbubbles when the chip was placed in an acoustic field. We think the microbubble-free microstructure-based device we present is easier to operate and more reliable, as it avoids microbubble preparation and maintenance, while showing high DNA fragmentation performance.

The enhancement of DNA fragmentation in a bench top ultrasonic water bath with needle-induced air bubbles: Simulation and experimental investigation

Shearing DNA to a certain size is the first step in many medical and biological applications, especially in next-generation gene sequencing technology. In this article, we introduced a highly efficient ultrasonic DNA fragmentation method enhanced by needle-induced air bubbles, which is easy to operate with high throughput. The principle of the bubble-enhanced sonication system is introduced and verified by flow field and acoustic simulations and experiments. Lambda DNA long chains and mouse genomic DNA short chains are used in the experiments for testing the performance of the bubble-enhanced ultrasonic DNA fragmentation system. Air bubbles are an effective enhancement agent for ultrasonic DNA fragmentation; they can obviously improve the sound pressure level in the whole solution, thus, achieving better absorption of ultrasound energy. Growing bubbles also have a stretched function on DNA molecule chains and form a huge pressure gradient in the solution, which is beneficial to DNA fragmentation. Purified λDNA is cut from 48.5 to 2 kbp in 5 min and cut to 300 bp in 30 min. Mouse genomic DNA (≈1400 bp) decreases to 400 bp in 5 min and then reduces to 200 bp in 30 min. This bubble-enhanced ultrasonic method enables widespread access to genomic DNA fragmentation in a standard ultrasonic water bath for many virus sequencing demands even without good medical facilities.

Bubble-enhanced ultrasonic microfluidic chip for rapid DNA fragmentation

DNA fragmentation is an essential process in developing genetic sequencing strategies, genetic research, as well as for the diagnosis of diseases with a genetic signature like cancer. Efficient on-chip DNA fragmentation protocols would be beneficial to process integration and open new opportunities for microfluidics in genetic applications. Here we present an acoustic microfluidic chip comprising an array of ultrasound-actuated microbubbles located at dedicated positions adjacent to a channel containing the DNA sample solution. The efficiency of the on-chip DNA fragmentation process arises mainly from tensile forces generated by acoustic streaming near the oscillating bubble interfaces, as well as a synergistic effect of streaming stress and ultrasonic cavitation. Acoustic microstreaming and the pressure distribution in the DNA channel were assessed by finite element simulation. We characterized the bubble-enhanced effect by measuring gene fragment size distributions with respect to different ultrasound parameters. For optimized on-chip conditions, purified lambda (λ) DNA (48.5 kbp) could be disrupted to fragments with an average size of 2 kbp after 30 s and down to 300 bp after 90 s. Mouse genomic DNA (1.4 kbp) fragmentation size decreased to 500 bp in 30 s and reduced further to 250 bp in 90 s. Bubble-induced fragmentation was more than 3 times faster than without bubbles. On-chip performance and process yield were found to be comparable to a sophisticated high-end commercial system. In this view, our new bubble-enhanced microfluidic approach is a promising tool for current and next generation sequencing platforms with high efficiency and good capacity. Moreover, the availability of an efficient on-chip DNA fragmentation process opens perspectives for implementing full molecular protocols on a single microfluidic platform.

Cavitation Enhancing Nanodroplets Mediate Efficient DNA Fragmentation in a Bench Top Ultrasonic Water Bath

A perfluorocarbon nanodroplet formulation is shown to be an effective cavitation enhancement agent, enabling rapid and consistent fragmentation of genomic DNA in a standard ultrasonic water bath. This nanodroplet-enhanced method produces genomic DNA libraries and next-generation sequencing results indistinguishable from DNA samples fragmented in dedicated commercial acoustic sonication equipment, and with higher throughput. This technique thus enables widespread access to fast bench-top genomic DNA fragmentation.

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 ~12×. 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.

A device for automated hydrodynamic shearing of genomic DNA

We describe a device for automated fragmentation of genomic DNA by hydrodynamic shearing using a filter screen with uniform pores. Human genomic DNA can be fragmented reproducibly to a range of 2000-12,000 bp by using various fluid flow rates and screens with 0.5-10-microm pores. The utilization of disposable screens eliminates sample cross-contamination and the tendency of device clogging, commonly encountered in single-orifice shearing devices.

A database of expressed genes from Cochliomyia hominivorax (Diptera: Calliphoridae)

We used an expressed sequence tag and 454 pyrosequencing approach to initiate a study of the genome of the screwworm, Cochliomyia hominivorax (Coquerel) (Diptera: Calliphoridae). Two normalized cDNA libraries were constructed from RNA isolated from embryos and second instar larvae from the Panama 95 strain. Approximately 5,400 clones from each library were sequenced from both the 5' and 3' directions using the Sanger method. In addition, double-stranded cDNA was prepared from random-primed polyA RNA purified from embryos, second-instar larvae, adult males, and adult females. These four cDNA samples were used for 454 pyrosequencing that produced approximately 300,000 independent sequences. Sequences were assembled into a database of assembled contigs and singletons and used to search public protein databases and annotate the sequences. The full database consists of 6,076 contigs and 58,221 singletons assembled from both the traditional expressed sequence tag (EST) and 454 sequences. Annotation of the data led to the identification of several gene coding regions with possible roles in sex determination in the screwworm. This database will facilitate the design of microarray and other experiments to study screwworm gene expression on a larger scale than previously possible.

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.

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.

Effects of Muclin (Dmbt1) deficiency on the gastrointestinal system

The Dmbt1 gene encodes alternatively spliced glycoproteins that are either membrane-associated or secreted epithelial products. Functions proposed for Dmbt1 include it being a tumor suppressor, having roles in innate immune defense and inflammation, and being a Golgi-sorting receptor in the exocrine pancreas. The heavily sulfated membrane glycoprotein mucin-like glycoprotein (Muclin) is a Dmbt1 product that is strongly expressed in organs of the gastrointestinal (GI) system. To explore Muclin's functions in the GI system, the Dmbt1 gene was targeted to produce Muclin-deficient mice. Muclin-deficient mice have normal body weight gain and are fertile. The Muclin-deficient mice did not develop GI tumors, even when crossed with mice lacking the known tumor suppressor p53. When colitis was induced by dextran sulfate sodium, there was no significant difference in disease severity in Muclin-deficient mice. Also, when acute pancreatitis was induced with supraphysiological caerulein, there was no difference in disease severity in the Muclin-deficient mice. Exocrine pancreatic function was impaired, as measured by attenuated neurohormonal-stimulated amylase release from Muclin-deficient acinar cells. Also, by [(35)S]Met/Cys pulse-chase analysis, traffic of newly synthesized protein to the stimulus-releasable pool was significantly retarded in Muclin-deficient cells compared with wild type. Thus Muclin deficiency impairs trafficking of regulated proteins to a stimulus-releasable pool in the exocrine pancreas.

Metagenomic analysis of the microbial community at Zodletone Spring (Oklahoma): insights into the genome of a member of the novel candidate division OD1

A metagenomic library was constructed from the anaerobic sediments of a mesophilic sulfur spring. Thirty-five bacterial 16S rRNA gene-containing clones were identified in this library. Analysis of a genomic fragment belonging to candidate division OD1 provided useful insights into the physiology and biochemistry of this novel, yet-uncultured candidate division.

Genomic comparison of plant pathogenic and nonpathogenic Serratia marcescens strains by suppressive subtractive hybridization

Cucurbit yellow vine disease (CYVD) is caused by disease-associated Serratia marcescens strains that have phenotypes significantly different from those of nonphytopathogenic strains. To identify the genetic differences responsible for pathogenicity-related phenotypes, we used a suppressive subtractive hybridization (SSH) strategy. S. marcescens strain Z01-A, isolated from CYVD-affected zucchini, was used as the tester, whereas rice endophytic S. marcescens strain R02-A (IRBG 502) was used as the driver. SSH revealed 48 sequences, ranging from 200 to 700 bp, that were present in Z01-A but absent in R02-A. Sequence analysis showed that a large proportion of these sequences resembled genes involved in synthesis of surface structures. By construction of a fosmid library, followed by colony hybridization, selection, and DNA sequencing, a phage gene cluster and a genome island containing a fimbrial-gene cluster were identified. Arrayed dot hybridization showed that the conservation of subtracted sequences among CYVD pathogenic and nonpathogenic S. marcescens strains varied. Thirty-four sequences were present only in pathogenic strains. Primers were designed based on one Z01-A-specific sequence, A79, and used in a multiplex PCR to discriminate between S. marcescens strains causing CYVD and those from other ecological niches.

Comparative analysis of the Borrelia garinii genome

Three members of the genus Borrelia (B.burgdorferi, B.garinii, B.afzelii) cause tick-borne borreliosis. Depending on the Borrelia species involved, the borreliosis differs in its clinical symptoms. Comparative genomics opens up a way to elucidate the underlying differences in Borrelia species. We analysed a low redundancy whole-genome shotgun (WGS) assembly of a B.garinii strain isolated from a patient with neuroborreliosis in comparison to the B.burgdorferi genome. This analysis reveals that most of the chromosome is conserved (92.7% identity on DNA as well as on amino acid level) in the two species, and no chromosomal rearrangement or larger insertions/deletions could be observed. Furthermore, two collinear plasmids (lp54 and cp26) seem to belong to the basic genome inventory of Borrelia species. These three collinear parts of the Borrelia genome encode 861 genes, which are orthologous in the two species examined. The majority of the genetic information of the other plasmids of B.burgdorferii is also present in B.garinii although orthology is not easy to define due to a high redundancy of the plasmid fraction. Yet, we did not find counterparts of the B.burgdorferi plasmids lp36 and lp38 or their respective gene repertoire in the B.garinii genome. Thus, phenotypic differences between the two species could be attributable to the presence or absence of these two plasmids as well as to the potentially positively selected genes.