Modifying the m6A brain methylome by ALKBH5-mediated demethylation: a new contender for synaptic tagging.

Synaptic plasticity processes, which underlie learning and memory formation, require RNA to be translated local to synapses. The synaptic tagging hypothesis has previously been proposed to explain how mRNAs are available at specific activated synapses. However how RNA is regulated, and which transcripts are silenced or processed as part of the tagging process is still unknown. Modification of RNA by N6-methyladenosine (m6A/m) influences the cellular fate of mRNA. Here, by advanced microscopy, we showed that m6A demethylation by the eraser protein ALKBH5 occurs at active synaptic ribosomes and at synapses during short term plasticity. We demonstrated that at activated glutamatergic post-synaptic sites, both the YTHDF1 and YTHDF3 reader and the ALKBH5 eraser proteins increase in co-localisation to m6A-modified RNAs; but only the readers showed high co-localisation to modified RNAs during late-stage plasticity. The YTHDF1 and YTHFDF3 readers also exhibited differential roles during synaptic maturation suggesting that temporal and subcellular abundance may determine specific function. m6A-sequencing of human parahippocampus brain tissue revealed distinct white and grey matter m6A methylome profiles indicating that cellular context is a fundamental factor dictating regulated pathways. However, in both neuronal and glial cell-rich tissue, m6A effector proteins are themselves modified and m6A epitranscriptional and posttranslational modification processes coregulate protein cascades. We hypothesise that the availability m6A effector protein machinery in conjunction with RNA modification, may be important in the formation of condensed synaptic nanodomain assemblies through liquid-liquid phase separation. Our findings support that m6A demethylation by ALKBH5 is an intrinsic component of the synaptic tagging hypothesis and a molecular switch which leads to alterations in the RNA methylome, synaptic dysfunction and potentially reversible disease states.

TCRs with Distinct Specificity Profiles Use Different Binding Modes to Engage an Identical Peptide-HLA Complex.

The molecular rules driving TCR cross-reactivity are poorly understood and, consequently, it is unclear the extent to which TCRs targeting the same Ag recognize the same off-target peptides. We determined TCR-peptide-HLA crystal structures and, using a single-chain peptide-HLA phage library, we generated peptide specificity profiles for three newly identified human TCRs specific for the cancer testis Ag NY-ESO-1157-165-HLA-A2. Two TCRs engaged the same central peptide feature, although were more permissive at peripheral peptide positions and, accordingly, possessed partially overlapping peptide specificity profiles. The third TCR engaged a flipped peptide conformation, leading to the recognition of off-target peptides sharing little similarity with the cognate peptide. These data show that TCRs specific for a cognate peptide recognize discrete peptide repertoires and reconciles how an individual's limited TCR repertoire following negative selection in the thymus is able to recognize a vastly larger antigenic pool.

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.

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.

Transcriptomic responses of mixed cultures of ascomycete fungi to lignocellulose using dual RNA-seq reveal inter-species antagonism and limited beneficial effects on CAZyme expression.

Gaining new knowledge through fungal monoculture responses to lignocellulose is a widely used approach that can lead to better cocktails for lignocellulose saccharification (the enzymatic release of sugars which are subsequently used to make biofuels). However, responses in lignocellulose mixed cultures are rarely studied in the same detail even though in nature fungi often degrade lignocellulose as mixed communities. Using a dual RNA-seq approach, we describe the first study of the transcriptional responses of wild-type strains of Aspergillus niger, Trichoderma reesei and Penicillium chrysogenum in two and three mixed species shake-flask cultures with wheat straw. Based on quantification of species-specific rRNA, a set of conditions was identified where mixed cultures could be sampled so as to obtain sufficient RNA-seq reads for analysis from each species. The number of differentially-expressed genes varied from a couple of thousand to fewer than one hundred. The proportion of carbohydrate active enzyme (CAZy) encoding transcripts was lower in the majority of the mixed cultures compared to the respective straw monocultures. A small subset of P. chrysogenum CAZy genes showed five to ten-fold significantly increased transcript abundance in a two-species mixed culture with T. reesei. However, a substantial number of T. reesei CAZy transcripts showed reduced abundance in mixed cultures. The highly induced genes in mixed cultures indicated that fungal antagonism was a major part of the mixed cultures. In line with this, secondary metabolite producing gene clusters showed increased transcript abundance in mixed cultures and also mixed cultures with T. reesei led to a decrease in the mycelial biomass of A. niger. Significantly higher monomeric sugar release from straw was only measured using a minority of the mixed culture filtrates and there was no overall improvement. This study demonstrates fungal interaction with changes in transcripts, enzyme activities and biomass in the mixed cultures and whilst there were minor beneficial effects for CAZy transcripts and activities, the competitive interaction between T. reesei and the other fungi was the most prominent feature of this study.

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.

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.

Planarian MBD2/3 is required for adult stem cell pluripotency independently of DNA methylation.

Planarian adult stem cells (pASCs) or neoblasts represent an ideal system to study the evolution of stem cells and pluripotency as they underpin an unrivaled capacity for regeneration. We wish to understand the control of differentiation and pluripotency in pASCs and to understand how conserved, convergent or divergent these mechanisms are across the Bilateria. Here we show the planarian methyl-CpG Binding Domain 2/3 (mbd2/3) gene is required for pASC differentiation during regeneration and tissue homeostasis. The genome does not have detectable levels of 5-methylcytosine (5(m)C) and we find no role for a potential DNA methylase. We conclude that MBD proteins may have had an ancient role in broadly controlling animal stem cell pluripotency, but that DNA methylation is not involved in planarian stem cell differentiation.

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.

Mapping the functional impact of non-coding regulatory elements in primary T cells through single-cell CRISPR screens.

BACKGROUND: Drug targets with genetic evidence are expected to increase clinical success by at least twofold. Yet, translating disease-associated genetic variants into functional knowledge remains a fundamental challenge of drug discovery. A key issue is that the vast majority of complex disease associations cannot be cleanly mapped to a gene. Immune disease-associated variants are enriched within regulatory elements found in T-cell-specific open chromatin regions. RESULTS: To identify genes and molecular programs modulated by these regulatory elements, we develop a CRISPRi-based single-cell functional screening approach in primary human T cells. Our pipeline enables the interrogation of transcriptomic changes induced by the perturbation of regulatory elements at scale. We first optimize an efficient CRISPRi protocol in primary CD4+ T cells via CROPseq vectors. Subsequently, we perform a screen targeting 45 non-coding regulatory elements and 35 transcription start sites and profile approximately 250,000 T -cell single-cell transcriptomes. We develop a bespoke analytical pipeline for element-to-gene (E2G) mapping and demonstrate that our method can identify both previously annotated and novel E2G links. Lastly, we integrate genetic association data for immune-related traits and demonstrate how our platform can aid in the identification of effector genes for GWAS loci. CONCLUSIONS: We describe "primary T cell crisprQTL" - a scalable, single-cell functional genomics approach for mapping regulatory elements to genes in primary human T cells. We show how this framework can facilitate the interrogation of immune disease GWAS hits and propose that the combination of experimental and QTL-based techniques is likely to address the variant-to-function problem.

Genome-wide transcriptional response of Trichoderma reesei to lignocellulose using RNA sequencing and comparison with Aspergillus niger.

BACKGROUND: A major part of second generation biofuel production is the enzymatic saccharification of lignocellulosic biomass into fermentable sugars. Many fungi produce enzymes that can saccarify lignocellulose and cocktails from several fungi, including well-studied species such as Trichoderma reesei and Aspergillus niger, are available commercially for this process. Such commercially-available enzyme cocktails are not necessarily representative of the array of enzymes used by the fungi themselves when faced with a complex lignocellulosic material. The global induction of genes in response to exposure of T. reesei to wheat straw was explored using RNA-seq and compared to published RNA-seq data and model of how A. niger senses and responds to wheat straw. RESULTS: In T. reesei, levels of transcript that encode known and predicted cell-wall degrading enzymes were very high after 24h exposure to straw (approximately 13% of the total mRNA) but were less than recorded in A. niger (approximately 19% of the total mRNA). Closer analysis revealed that enzymes from the same glycoside hydrolase families but different carbohydrate esterase and polysaccharide lyase families were up-regulated in both organisms. Accessory proteins which have been hypothesised to possibly have a role in enhancing carbohydrate deconstruction in A. niger were also uncovered in T. reesei and categories of enzymes induced were in general similar to those in A. niger. Similarly to A. niger, antisense transcripts are present in T. reesei and their expression is regulated by the growth condition. CONCLUSIONS: T. reesei uses a similar array of enzymes, for the deconstruction of a solid lignocellulosic substrate, to A. niger. This suggests a conserved strategy towards lignocellulose degradation in both saprobic fungi. This study provides a basis for further analysis and characterisation of genes shown to be highly induced in the presence of a lignocellulosic substrate. The data will help to elucidate the mechanism of solid substrate recognition and subsequent degradation by T. reesei and provide information which could prove useful for efficient production of second generation biofuels.

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

BACKGROUND: Comparative 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. RESULTS: Here, 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. CONCLUSIONS: We suggest that the genome of S. arboricolus will be useful in future comparative genomics analysis of the Saccharomyces sensu stricto yeasts.

Use of Bulk Segregant Analysis for Determining the Genetic Basis of Azole Resistance in the Opportunistic Pathogen Aspergillus fumigatus.

A sexual cycle was described in 2009 for the opportunistic fungal pathogen Aspergillus fumigatus, opening up for the first time the possibility of using techniques reliant on sexual crossing for genetic analysis. The present study was undertaken to evaluate whether the technique 'bulk segregant analysis' (BSA), which involves detection of differences between pools of progeny varying in a particular trait, could be applied in conjunction with next-generation sequencing to investigate the underlying basis of monogenic traits in A. fumigatus. Resistance to the azole antifungal itraconazole was chosen as a model, with a dedicated bioinformatic pipeline developed to allow identification of SNPs that differed between the resistant progeny pool and resistant parent compared to the sensitive progeny pool and parent. A clinical isolate exhibiting monogenic resistance to itraconazole of unknown basis was crossed to a sensitive parent and F1 progeny used in BSA. In addition, the use of backcrossing and increasing the number in progeny pools was evaluated as ways to enhance the efficiency of BSA. Use of F1 pools of 40 progeny led to the identification of 123 candidate genes with SNPs distributed over several contigs when aligned to an A1163 reference genome. Successive rounds of backcrossing enhanced the ability to identify specific genes and a genomic region, with BSA of progeny (using 40 per pool) from a third backcross identifying 46 genes with SNPs, and BSA of progeny from a sixth backcross identifying 20 genes with SNPs in a single 292 kb region of the genome. The use of an increased number of 80 progeny per pool also increased the resolution of BSA, with 29 genes demonstrating SNPs between the different sensitive and resistant groupings detected using progeny from just the second backcross with the majority of variants located on the same 292 kb region. Further bioinformatic analysis of the 292 kb region identified the presence of a cyp51A gene variant resulting in a methionine to lysine (M220K) change in the CYP51A protein, which was concluded to be the causal basis of the observed resistance to itraconazole. The future use of BSA in genetic analysis of A. fumigatus is discussed.

Site-specific recombination in Schizosaccharomyces pombe and systematic assembly of a 400kb transgene array in mammalian cells using the integrase of Streptomyces phage phiBT1.

We have established the integrase of the Streptomyces phage phiBT1 as a tool for eukaryotic genome manipulation. We show that the phiBT1 integrase promotes efficient reciprocal and conservative site-specific recombination in vertebrate cells and in Schizosaccharomyces pombe, thus establishing the utility of this protein for genome manipulation in a wide range of eukaryotes. We show that the phiBT1 integrase can be used in conjunction with Cre recombinase to promote the iterative integration of transgenic DNA. We describe five cycles of iterative integration of a candidate mouse centromeric sequence 80 kb in length into a human mini-chromosome within a human-Chinese hamster hybrid cell line. These results establish the generality of the iterative site-specific integration technique.

Gfi1aa and Gfi1b set the pace for primitive erythroblast differentiation from hemangioblasts in the zebrafish embryo.

The transcriptional repressors Gfi1(a) and Gfi1b are epigenetic regulators with unique and overlapping roles in hematopoiesis. In different contexts, Gfi1 and Gfi1b restrict or promote cell proliferation, prevent apoptosis, influence cell fate decisions, and are essential for terminal differentiation. Here, we show in primitive red blood cells (prRBCs) that they can also set the pace for cellular differentiation. In zebrafish, prRBCs express 2 of 3 zebrafish Gfi1/1b paralogs, Gfi1aa and Gfi1b. The recently identified zebrafish gfi1aa gene trap allele qmc551 drives erythroid green fluorescent protein (GFP) instead of Gfi1aa expression, yet homozygous carriers have normal prRBCs. prRBCs display a maturation defect only after splice morpholino-mediated knockdown of Gfi1b in gfi1aa qmc551 homozygous embryos. To study the transcriptome of the Gfi1aa/1b double-depleted cells, we performed an RNA-Seq experiment on GFP-positive prRBCs sorted from 20-hour-old embryos that were heterozygous or homozygous for gfi1aa qmc551 , as well as wt or morphant for gfi1b We subsequently confirmed and extended these data in whole-mount in situ hybridization experiments on newly generated single- and double-mutant embryos. Combined, the data showed that in the absence of Gfi1aa, the synchronously developing prRBCs were delayed in activating late erythroid differentiation, as they struggled to suppress early erythroid and endothelial transcription programs. The latter highlighted the bipotent nature of the progenitors from which prRBCs arise. In the absence of Gfi1aa, Gfi1b promoted erythroid differentiation as stepwise loss of wt gfi1b copies progressively delayed Gfi1aa-depleted prRBCs even further, showing that Gfi1aa and Gfi1b together set the pace for prRBC differentiation from hemangioblasts.

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.

Rearranging the centromere of the human Y chromosome with phiC31 integrase.

We have investigated the ability of the integrase from the Streptomyces phiC31 'phage to either delete or invert 1 Mb of DNA around the centromere of the human Y chromosome in chicken DT40 hybrid somatic cells. Reciprocal and conservative site-specific recombination was observed in 54% of cells expressing the integrase. The sites failed to recombine in the remaining cells because the sites had been damaged. The sequences of the damaged sites indicated that the damage arose as a result of repair of recombination intermediates by host cell pathways. The liability of recombination intermediates to damage is consistent with what is known about the mechanism of serine recombinase reactions. The structures of the products of the chromosome rearrangements were consistent with the published sequence of the Y chromosome indicating that the assembly of the highly repeated region between the sites is accurate to a resolution of about 50 kb. Mini-chromosomes lacking a centromere were not recovered which also suggested that neo-centromere formation occurs infrequently in vertebrate somatic cells. No ectopic recombination was observed between a phiC31 integrase attB site and the chicken genome.

Iterative in vivo assembly of large and complex transgenes by combining the activities of phiC31 integrase and Cre recombinase.

We have used the phiC31 integrase to introduce large DNA sequences into a vertebrate genome and measure the efficiency of integration of intact DNA as a function of insert size. Inserts of 110 kb and 140 kb in length may be integrated with about 25% and 10% efficiency respectively. In order to overcome the problems of constructing transgenes longer than approximately 150 kb we have established a method that we call; 'Iterative Site Specific Integration' (ISSI). ISSI combines the activities of phiC31 integrase and Cre recombinase to enable the iterative and serial integration of transgenic DNA sequences. In principle the procedure may be repeated an arbitrary number of times and thereby allow the integration of tracts of DNA many hundreds of kilobase pairs long. In practice it may be limited by the time needed to check the accuracy of integration at each step of the procedure. We describe two ISSI experiments, in one of which we have constructed a complex array of vertebrate centromeric sequences of 150 kb in size. The principle that underlies ISSI is applicable to transgenesis in all organisms. ISSI may thus facilitate the reconstitution of biosynthetic pathways encoded by many different genes in transgenic plants, the assembly of large vertebrate loci as transgenes and the synthesis of complete genomes in bacteria.

Chromosome engineering in DT40 cells and mammalian centromere function.

Chromosome engineering is the term given to procedures which modify the long range structure of a chromosome by homologous and site specific recombination or by telomere directed chromosome breakage. DT40 cells are uniquely powerful for chromosome engineering because mammalian chromosomes may be moved into them, efficiently modified and then moved back into a mammalian cell lines (Dieken et al., 1996). The high rate of sequence targeting seen in DT40 cells carrying human chromosomes is necessary but not sufficient for chromosome engineering. The ability to either delete or introduce long tracts of DNA subsequent to a sequence targeting reaction depends upon the use of site specific recombinases. We have made important progress in the development of this technology in the past few years and much of this review will be used to describe this work.

Development of Gene-Based SSR Markers in Winged Bean (Psophocarpus tetragonolobus (L.) DC.) for Diversity Assessment.

Winged bean (Psophocarpus tetragonolobus) is an herbaceous multipurpose legume grown in hot and humid countries as a pulse, vegetable (leaves and pods), or root tuber crop depending on local consumption preferences. In addition to its different nutrient-rich edible parts which could contribute to food and nutritional security, it is an efficient nitrogen fixer as a component of sustainable agricultural systems. Generating genetic resources and improved lines would help to accelerate the breeding improvement of this crop, as the lack of improved cultivars adapted to specific environments has been one of the limitations preventing wider use. A transcriptomic de novo assembly was constructed from four tissues: leaf, root, pod, and reproductive tissues from Malaysian accessions, comprising of 198,554 contigs with a N50 of 1462 bp. Of these, 138,958 (70.0%) could be annotated. Among 9682 genic simple sequence repeat (SSR) motifs identified (excluding monomer repeats), trinucleotide-repeats were the most abundant (4855), followed by di-nucleotide (4500) repeats. A total of 18 SSR markers targeting di- and tri-nucleotide repeats have been validated as polymorphic markers based on an initial assessment of nine genotypes originated from five countries. A cluster analysis revealed provisional clusters among this limited, yet diverse selection of germplasm. The developed assembly and validated genic SSRs in this study provide a foundation for a better understanding of the plant breeding system for the genetic improvement of winged bean.