Utilization of Probabilistic Models in Short Read Assembly from Second-Generation Sequencing

With the advent of cheaper and faster DNA sequencing technologies, assembly methods have greatly changed. Instead of outputting reads that are thousands of base pairs long, new sequencers parallelize the task by producing read lengths between 35 and 400 base pairs. Reconstructing an organism’s genome from these millions of reads is a computationally expensive task. Our algorithm solves this problem by organizing and indexing the reads using n-grams, which are short, fixed-length DNA sequences of length n. These n-grams are used to efficiently locate putative read joins, thereby eliminating the need to perform an exhaustive search over all possible read pairs. Our goal was develop a novel n-gram method for the assembly of genomes from next-generation sequencers. Specifically, a probabilistic, iterative approach was utilized to determine the most likely reads to join through development of a new metric that models the probability of any two arbitrary reads being joined together. Tests were run using simulated short read data based on randomly created genomes ranging in lengths from 10,000 to 100,000 nucleotides with 16 to 20x coverage. We were able to successfully re-assemble entire genomes up to 100,000 nucleotides in length.

SNP array profiling of childhood adrenocortical tumors reveals distinct pathways of tumorigenesis and highlights candidate driver genes

Childhood adrenocortical tumors (ACT) are rare malignancies, except in southern Brazil, where a higher incidence rate is associated to a high frequency of the founder R337H TP53 mutation. To date, copy number alterations in these tumors have only been analyzed by low-resolution comparative genomic hybridization.

Generation and analysis of a mouse intestinal metatranscriptome through Illumina based RNA-sequencing

With the advent of high through-put sequencing (HTS), the emerging science of metagenomics is transforming our understanding of the relationships of microbial communities with their environments. While metagenomics aims to catalogue the genes present in a sample through assessing which genes are actively expressed, metatranscriptomics can provide a mechanistic understanding of community inter-relationships. To achieve these goals, several challenges need to be addressed from sample preparation to sequence processing, statistical analysis and functional annotation. Here we use an inbred non-obese diabetic (NOD) mouse model in which germ-free animals were colonized with a defined mixture of eight commensal bacteria, to explore methods of RNA extraction and to develop a pipeline for the generation and analysis of metatranscriptomic data. Applying the Illumina HTS platform, we sequenced 12 NOD cecal samples prepared using multiple RNA-extraction protocols. The absence of a complete set of reference genomes necessitated a peptide-based search strategy. Up to 16% of sequence reads could be matched to a known bacterial gene. Phylogenetic analysis of the mapped ORFs revealed a distribution consistent with ribosomal RNA, the majority from Bacteroides or Clostridium species. To place these HTS data within a systems context, we mapped the relative abundance of corresponding Escherichia coli homologs onto metabolic and protein-protein interaction networks. These maps identified bacterial processes with components that were well-represented in the datasets. In summary this study highlights the potential of exploiting the economy of HTS platforms for metatranscriptomics.

Comparison of the gene expression profiles from normal and Fgfrl1 deficient mouse kidneys reveals downstream targets of Fgfrl1 signaling

Fgfrl1 (fibroblast growth factor receptor-like 1) is a transmembrane receptor that is essential for the development of the metanephric kidney. It is expressed in all nascent nephrogenic structures and in the ureteric bud. Fgfrl1 null mice fail to develop the metanephric kidneys. Mutant kidney rudiments show a dramatic reduction of ureteric branching and a lack of mesenchymal-to-epithelial transition. Here, we compared the expression profiles of wildtype and Fgfrl1 mutant kidneys to identify genes that act downstream of Fgfrl1 signaling during the early steps of nephron formation. We detected 56 differentially expressed transcripts with 2-fold or greater reduction, among them many genes involved in Fgf, Wnt, Bmp, Notch, and Six/Eya/Dach signaling. We validated the microarray data by qPCR and whole-mount in situ hybridization and showed the expression pattern of candidate genes in normal kidneys. Some of these genes might play an important role during early nephron formation. Our study should help to define the minimal set of genes that is required to form a functional nephron.

Potential drought stress in a Swiss mountain catchment-Ensemble forecasting of high mountain soil moisture reveals a drastic decrease, despite major uncertainties

Climate change is expected to profoundly influence the hydrosphere of mountain ecosystems. The focus of current process-based research is centered on the reaction of glaciers and runoff to climate change; spatially explicit impacts on soil moisture remain widely neglected. We spatio-temporally analyzed the impact of the climate on soil moisture in a mesoscale high mountain catchment to facilitate the development of mitigation and adaptation strategies at the level of vegetation patterns. Two regional climate models were downscaled using three different approaches (statistical downscaling, delta change, and direct use) to drive a hydrological model (WaSiM-ETH) for reference and scenario period (1960–1990 and 2070–2100), resulting in an ensemble forecast of six members. For all ensembles members we found large changes in temperature, resulting in decreasing snow and ice storage and earlier runoff, but only small changes in evapotranspiration. The occurrence of downscaled dry spells was found to fluctuate greatly, causing soil moisture depletion and drought stress potential to show high variability in both space and time. In general, the choice of the downscaling approach had a stronger influence on the results than the applied regional climate model. All of the results indicate that summer soil moisture decreases, which leads to more frequent declines below a critical soil moisture level and an advanced evapotranspiration deficit. Forests up to an elevation of 1800 m a.s.l. are likely to be threatened the most, while alpine areas and most pastures remain nearly unaffected. Nevertheless, the ensemble variability was found to be extremely high and should be interpreted as a bandwidth of possible future drought stress situations.

Urinary metabolite profiling reveals CYP1A2-mediated metabolism of NSC686288 (aminoflavone)

NSC686288 [aminoflavone (AF)], a candidate chemotherapeutic agent, possesses a unique antiproliferative profile against tumor cells. Metabolic bioactivation of AF by drug-metabolizing enzymes, especially CYP1A monooxygenases, has been implicated as an underlying mechanism for its selective cytotoxicity in several cell culture-based studies. However, in vivo metabolism of AF has not been investigated in detail. In this study, the structural identities of 13 AF metabolites (12 of which are novel) in mouse urine or from microsomal incubations, including three monohydroxy-AFs, two dihydroxy-AFs and their sulfate and glucuronide conjugates, as well as one N-glucuronide, were determined by accurate mass measurements and liquid chromatography-tandem mass spectrometry fragmentation patterns, and a comprehensive map of the AF metabolic pathways was constructed. Significant differences between wild-type and Cyp1a2-null mice, within the relative composition of urinary metabolites of AF, demonstrated that CYP1A2-mediated regioselective oxidation was a major contributor to the metabolism of AF. Comparisons between wild-type and CYP1A2-humanized mice further revealed interspecies differences in CYP1A2-mediated catalytic activity. Incubation of AF with liver microsomes from all three mouse lines and with pooled human liver microsomes confirmed the observations from urinary metabolite profiling. Results from enzyme kinetic analysis further indicated that in addition to CYP1A P450s, CYP2C P450s may also play some role in the metabolism of AF.

Gene expression profiling reveals consistent differences between clinical samples of human leukaemias and their model cell lines

Microarray gene expression profiles of fresh clinical samples of chronic myeloid leukaemia in chronic phase, acute promyelocytic leukaemia and acute monocytic leukaemia were compared with profiles from cell lines representing the corresponding types of leukaemia (K562, NB4, HL60). In a hierarchical clustering analysis, all clinical samples clustered separately from the cell lines, regardless of leukaemic subtype. Gene ontology analysis showed that cell lines chiefly overexpressed genes related to macromolecular metabolism, whereas in clinical samples genes related to the immune response were abundantly expressed. These findings must be taken into consideration when conclusions from cell line-based studies are extrapolated to patients.

Real-time quantitative broad-range PCR assay for detection of the 16S rRNA gene followed by sequencing for species identification

Here we determined the analytical sensitivities of broad-range real-time PCR-based assays employing one of three different genomic DNA extraction protocols in combination with one of three different primer pairs targeting the 16S rRNA gene to detect a panel of 22 bacterial species. DNA extraction protocol III, using lysozyme, lysostaphin, and proteinase K, followed by PCR with the primer pair Bak11W/Bak2, giving amplicons of 796 bp in length, showed the best overall sensitivity, detecting DNA of 82% of the strains investigated at concentrations of < or =10(2) CFU in water per reaction. DNA extraction protocols I and II, using less enzyme treatment, combined with other primer pairs giving shorter amplicons of 466 bp and 342 or 346 bp, respectively, were slightly more sensitive for the detection of gram-negative but less sensitive for the detection of gram-positive bacteria. The obstacle of detecting background DNA in blood samples spiked with bacteria was circumvented by introducing a broad-range hybridization probe, and this preserved the minimal detection limits observed in samples devoid of blood. Finally, sequencing of the amplicons generated using the primer pair Bak11W/Bak2 allowed species identification of the detected bacterial DNA. Thus, broad-spectrum PCR targeting the 16S rRNA gene in the quantitative real-time format can achieve an analytical sensitivity of 1 to 10 CFU per reaction in water, avoid detection of background DNA with the introduction of a broad-range probe, and generate amplicons that allow species identification of the detected bacterial DNA by sequencing. These prerequisites are important for its application to blood-containing patient samples.

Crystal structure of the platelet glycoprotein Ib(alpha) N-terminal domain reveals an unmasking mechanism for receptor activation

Glycoprotein Ib (GPIb) is a platelet receptor with a critical role in mediating the arrest of platelets at sites of vascular damage. GPIb binds to the A1 domain of von Willebrand factor (vWF-A1) at high blood shear, initiating platelet adhesion and contributing to the formation of a thrombus. To investigate the molecular basis of GPIb regulation and ligand binding, we have determined the structure of the N-terminal domain of the GPIb(alpha) chain (residues 1-279). This structure is the first determined from the cell adhesion/signaling class of leucine-rich repeat (LRR) proteins and reveals the topology of the characteristic disulfide-bonded flanking regions. The fold consists of an N-terminal beta-hairpin, eight leucine-rich repeats, a disulfide-bonded loop, and a C-terminal anionic region. The structure also demonstrates a novel LRR motif in the form of an M-shaped arrangement of three tandem beta-turns. Negatively charged binding surfaces on the LRR concave face and anionic region indicate two-step binding kinetics to vWF-A1, which can be regulated by an unmasking mechanism involving conformational change of a key loop. Using molecular docking of the GPIb and vWF-A1 crystal structures, we were also able to model the GPIb.vWF-A1 complex.

Combined voltage and calcium epifluorescence imaging in vitro and in vivo reveals subthreshold and suprathreshold dynamics of mouse barrel cortex

Cortical dynamics can be imaged at high spatiotemporal resolution with voltage-sensitive dyes (VSDs) and calcium-sensitive dyes (CaSDs). We combined these two imaging techniques using epifluorescence optics together with whole cell recordings to measure the spatiotemporal dynamics of activity in the mouse somatosensory barrel cortex in vitro and in the supragranular layers in vivo. The two optical signals reported distinct aspects of cortical function. VSD fluorescence varied linearly with membrane potential and was dominated by subthreshold postsynaptic potentials, whereas the CaSD signal predominantly reflected local action potential firing. Combining VSDs and CaSDs allowed us to monitor the synaptic drive and the spiking activity of a given area at the same time in the same preparation. The spatial extent of the two dye signals was different, with VSD signals spreading further than CaSD signals, reflecting broad subthreshold and narrow suprathreshold receptive fields. Importantly, the signals from the dyes were differentially affected by pharmacological manipulations, stimulation strength, and depth of isoflurane anesthesia. Combined VSD and CaSD measurements can therefore be used to specify the temporal and spatial relationships between subthreshold and suprathreshold activity of the neocortex.