Sunir Malla

Uncovering the genome-wide transcriptional responses of the filamentous fungus Aspergillus niger to lignocellulose using RNA sequencing.

A key challenge in the production of second generation biofuels is the conversion of lignocellulosic substrates into fermentable sugars. Enzymes, particularly those from fungi, are a central part of this process, and many have been isolated and characterised. However, relatively little is known of how fungi respond to lignocellulose and produce the enzymes necessary for dis-assembly of plant biomass. We studied the physiological response of the fungus Aspergillus niger when exposed to wheat straw as a model lignocellulosic substrate. Using RNA sequencing we showed that, 24 hours after exposure to straw, gene expression of known and presumptive plant cell wall-degrading enzymes represents a huge investment for the cells (about 20% of the total mRNA). Our results also uncovered new esterases and surface interacting proteins that might form part of the fungal arsenal of enzymes for the degradation of plant biomass. Using transcription factor deletion mutants (xlnR and creA) to study the response to both lignocellulosic substrates and low carbon source concentrations, we showed that a subset of genes coding for degradative enzymes is induced by starvation. Our data support a model whereby this subset of enzymes plays a scouting role under starvation conditions, testing for available complex polysaccharides and liberating inducing sugars, that triggers the subsequent induction of the majority of hydrolases. We also showed that antisense transcripts are abundant and that their expression can be regulated by growth conditions.

Nanopore sequencing and assembly of a human genome with ultra-long reads

We report the sequencing and assembly of a reference genome for the human GM12878 Utah/Ceph cell line using the MinION (Oxford Nanopore Technologies) nanopore sequencer. 91.2 Gb of sequence data, representing ∼30× theoretical coverage, were produced. Reference-based alignment enabled detection of large structural variants and epigenetic modifications. De novo assembly of nanopore reads alone yielded a contiguous assembly (NG50 ∼3 Mb). We developed a protocol to generate ultra-long reads (N50 > 100 kb, read lengths up to 882 kb). Incorporating an additional 5× coverage of these ultra-long reads more than doubled the assembly contiguity (NG50 ∼6.4 Mb). The final assembled genome was 2,867 million bases in size, covering 85.8% of the reference. Assembly accuracy, after incorporating complementary short-read sequencing data, exceeded 99.8%. Ultra-long reads enabled assembly and phasing of the 4-Mb major histocompatibility complex (MHC) locus in its entirety, measurement of telomere repeat length, and closure of gaps in the reference human genome assembly GRCh38.

Defining the molecular profile of planarian pluripotent stem cells using a combinatorial RNA-seq, RNA interference and irradiation approach

BackgroundPlanarian stem cells, or neoblasts, drive the almost unlimited regeneration capacities of freshwater planarians. Neoblasts are traditionally described by their morphological features and by the fact that they are the only proliferative cell type in asexual planarians. Therefore, they can be specifically eliminated by irradiation. Irradiation, however, is likely to induce transcriptome-wide changes in gene expression that are not associated with neoblast ablation. This has affected the accurate description of their specific transcriptomic profile.ResultsWe introduce the use of Smed-histone-2B RNA interference (RNAi) for genetic ablation of neoblast cells in Schmidtea mediterranea as an alternative to irradiation. We characterize the rapid, neoblast-specific phenotype induced by Smed-histone-2B RNAi, resulting in neoblast ablation. We compare and triangulate RNA-seq data after using both irradiation and Smed-histone-2B RNAi over a time course as means of neoblast ablation. Our analyses show that Smed-histone-2B RNAi eliminates neoblast gene expression with high specificity and discrimination from gene expression in other cellular compartments. We compile a high confidence list of genes downregulated by both irradiation and Smed-histone-2B RNAi and validate their expression in neoblast cells. Lastly, we analyze the overall expression profile of neoblast cells.ConclusionsOur list of neoblast genes parallels their morphological features and is highly enriched for nuclear components, chromatin remodeling factors, RNA splicing factors, RNA granule components and the machinery of cell division. Our data reveal that the regulation of planarian stem cells relies on posttranscriptional regulatory mechanisms and suggest that planarians are an ideal model for this understudied aspect of stem cell biology.

Accelerated growth in the absence of DNA replication origins

DNA replication initiates at defined sites called origins, which serve as binding sites for initiator proteins that recruit the replicative machinery. Origins differ in number and structure across the three domains of life and their properties determine the dynamics of chromosome replication. Bacteria and some archaea replicate from single origins, whereas most archaea and all eukaryotes replicate using multiple origins. Initiation mechanisms that rely on homologous recombination operate in some viruses. Here we show that such mechanisms also operate in archaea. We use deep sequencing to study replication in Haloferax volcanii and identify four chromosomal origins of differing activity. Deletion of individual origins results in perturbed replication dynamics and reduced growth. However, a strain lacking all origins has no apparent defects and grows significantly faster than wild type. Origin-less cells initiate replication at dispersed sites rather than at discrete origins and have an absolute requirement for the recombinase RadA, unlike strains lacking individual origins. Our results demonstrate that homologous recombination alone can efficiently initiate the replication of an entire cellular genome. This raises the question of what purpose replication origins serve and why they have evolved.

Transient Exposure to Low Levels of Insecticide Affects Metabolic Networks of Honeybee Larvae

The survival of a species depends on its capacity to adjust to changing environmental conditions, and new stressors. Such new, anthropogenic stressors include the neonicotinoid class of crop-protecting agents, which have been implicated in the population declines of pollinating insects, including honeybees (Apis mellifera). The low-dose effects of these compounds on larval development and physiological responses have remained largely unknown. Over a period of 15 days, we provided syrup tainted with low levels (2 µg/L−1) of the neonicotinoid insecticide imidacloprid to beehives located in the field. We measured transcript levels by RNA sequencing and established lipid profiles using liquid chromatography coupled with mass spectrometry from worker-bee larvae of imidacloprid-exposed (IE) and unexposed, control (C) hives. Within a catalogue of 300 differentially expressed transcripts in larvae from IE hives, we detect significant enrichment of genes functioning in lipid-carbohydrate-mitochondrial metabolic networks. Myc-involved transcriptional response to exposure of this neonicotinoid is indicated by overrepresentation of E-box elements in the promoter regions of genes with altered expression. RNA levels for a cluster of genes encoding detoxifying P450 enzymes are elevated, with coordinated downregulation of genes in glycolytic and sugar-metabolising pathways. Expression of the environmentally responsive Hsp90 gene is also reduced, suggesting diminished buffering and stability of the developmental program. The multifaceted, physiological response described here may be of importance to our general understanding of pollinator health. Muscles, for instance, work at high glycolytic rates and flight performance could be impacted should low levels of this evolutionarily novel stressor likewise induce downregulation of energy metabolising genes in adult pollinators.

High quality de novo sequencing and assembly of the Saccharomyces arboricolus genome

BackgroundComparative genomics is a formidable tool to identify functional elements throughout a genome. In the past ten years, studies in the budding yeast Saccharomyces cerevisiae and a set of closely related species have been instrumental in showing the benefit of analyzing patterns of sequence conservation. Increasing the number of closely related genome sequences makes the comparative genomics approach more powerful and accurate.ResultsHere, we report the genome sequence and analysis of Saccharomyces arboricolus, a yeast species recently isolated in China, that is closely related to S. cerevisiae. We obtained high quality de novo sequence and assemblies using a combination of next generation sequencing technologies, established the phylogenetic position of this species and considered its phenotypic profile under multiple environmental conditions in the light of its gene content and phylogeny.ConclusionsWe suggest that the genome of S. arboricolus will be useful in future comparative genomics analysis of the Saccharomyces sensu stricto yeasts.

The role of carbon starvation in the induction of enzymes that degrade plant-derived carbohydrates in Aspergillus niger

Fungi are an important source of enzymes for saccharification of plant polysaccharides and production of biofuels. Understanding of the regulation and induction of expression of genes encoding these enzymes is still incomplete. To explore the induction mechanism, we analysed the response of the industrially important fungus Aspergillus niger to wheat straw, with a focus on events occurring shortly after exposure to the substrate. RNA sequencing showed that the transcriptional response after 6h of exposure to wheat straw was very different from the response at 24h of exposure to the same substrate. For example, less than half of the genes encoding carbohydrate active enzymes that were induced after 24h of exposure to wheat straw, were also induced after 6h exposure. Importantly, over a third of the genes induced after 6h of exposure to wheat straw were also induced during 6h of carbon starvation, indicating that carbon starvation is probably an important factor in the early response to wheat straw. The up-regulation of the expression of a high number of genes encoding CAZymes that are active on plant-derived carbohydrates during early carbon starvation suggests that these enzymes could be involved in a scouting role during starvation, releasing inducing sugars from complex plant polysaccharides. We show, using proteomics, that carbon-starved cultures indeed release CAZymes with predicted activity on plant polysaccharides. Analysis of the enzymatic activity and the reaction products, indicates that these proteins are enzymes that can degrade various plant polysaccharides to generate both known, as well as potentially new, inducers of CAZymes.

Real-time selective sequencing using nanopore technology

The Oxford Nanopore Technologies MinION sequencer enables the selection of specific DNA molecules for sequencing by reversing the driving voltage across individual nanopores. To directly select molecules for sequencing, we used dynamic time warping to match reads to reference sequences. We demonstrate our open-source Read Until software in real-time selective sequencing of regions within small genomes, individual amplicon enrichment and normalization of an amplicon set.

The dynamics of genome replication using deep sequencing

Eukaryotic genomes are replicated from multiple DNA replication origins. We present complementary deep sequencing approaches to measure origin location and activity in Saccharomyces cerevisiae. Measuring the increase in DNA copy number during a synchronous S-phase allowed the precise determination of genome replication. To map origin locations, replication forks were stalled close to their initiation sites; therefore, copy number enrichment was limited to origins. Replication timing profiles were generated from asynchronous cultures using fluorescence-activated cell sorting. Applying this technique we show that the replication profiles of haploid and diploid cells are indistinguishable, indicating that both cell types use the same cohort of origins with the same activities. Finally, increasing sequencing depth allowed the direct measure of replication dynamics from an exponentially growing culture. This is the first time this approach, called marker frequency analysis, has been successfully applied to a eukaryote. These data provide a high-resolution resource and methodological framework for studying genome biology.

High through-put sequencing of the Parhyale hawaiensis mRNAs and microRNAs to aid comparative developmental studies.

Understanding the genetic and evolutionary basis of animal morphological diversity will require comparative developmental studies that use new model organisms. This necessitates development of tools for the study of genetics and also the generation of sequence information of the organism to be studied. The development of next generation sequencing technology has enabled quick and cost effective generation of sequence information. Parhyale hawaiensis has emerged as a model organism of choice due to the development of advanced molecular tools, thus P. hawaiensis genetic information will help drive functional studies in this organism. Here we present a transcriptome and miRNA collection generated using next generation sequencing platforms. We generated approximately 1.7 million reads from a P. hawaiensis cDNA library constructed from embryos up to the germ band stage. These reads were assembled into a dataset comprising 163,501 transcripts. Using the combined annotation of Annot8r and pfam2go, Gene Ontology classifications was assigned to 20,597 transcripts. Annot8r was used to provide KEGG orthology to our transcript dataset. A total of 25,292 KEGG pathway assignments were defined and further confirmed with reciprocal blast against the NCBI nr protein database. This has identified many P. hawaiensis gene orthologs of key conserved signalling pathways involved in development. We also generated small RNA sequences from P. hawaiensis, identifying 55 conserved miRNAs. Sequenced small RNAs that were not annotated by stringent comparison to mirBase were used to search the Daphnia pulex for possible novel miRNAs. Using a conservative approach, we have identified 51 possible miRNA candidates conserved in the Daphnia pulex genome, which could be potential crustacean/arthropod specific miRNAs. Our study presents gene and miRNA discovery in a new model organism that does not have a sequenced genome. The data provided by our work will be valuable for the P. hawaiensis community as well as the wider evolutionary developmental biology community.