Menu
September 22, 2019  |  

Single-cell mRNA isoform diversity in the mouse brain.

Alternative mRNA isoform usage is an important source of protein diversity in mammalian cells. This phenomenon has been extensively studied in bulk tissues, however, it remains unclear how this diversity is reflected in single cells.Here we use long-read sequencing technology combined with unique molecular identifiers (UMIs) to reveal patterns of alternative full-length isoform expression in single cells from the mouse brain. We found a surprising amount of isoform diversity, even after applying a conservative definition of what constitutes an isoform. Genes tend to have one or a few isoforms highly expressed and a larger number of isoforms expressed at a low level. However, for many genes, nearly every sequenced mRNA molecule was unique, and many events affected coding regions suggesting previously unknown protein diversity in single cells. Exon junctions in coding regions were less prone to splicing errors than those in non-coding regions, indicating purifying selection on splice donor and acceptor efficiency.Our findings indicate that mRNA isoform diversity is an important source of biological variability also in single cells.


September 22, 2019  |  

Clonal distribution of BCR-ABL1 mutations and splice isoforms by single-molecule long-read RNA sequencing.

The evolution of mutations in the BCR-ABL1 fusion gene transcript renders CML patients resistant to tyrosine kinase inhibitor (TKI) based therapy. Thus screening for BCR-ABL1 mutations is recommended particularly in patients experiencing poor response to treatment. Herein we describe a novel approach for the detection and surveillance of BCR-ABL1 mutations in CML patients.To detect mutations in the BCR-ABL1 transcript we developed an assay based on the Pacific Biosciences (PacBio) sequencing technology, which allows for single-molecule long-read sequencing of BCR-ABL1 fusion transcript molecules. Samples from six patients with poor response to therapy were analyzed both at diagnosis and follow-up. cDNA was generated from total RNA and a 1,6 kb fragment encompassing the BCR-ABL1 transcript was amplified using long range PCR. To estimate the sensitivity of the assay, a serial dilution experiment was performed.Over 10,000 full-length BCR-ABL1 sequences were obtained for all samples studied. Through the serial dilution analysis, mutations in CML patient samples could be detected down to a level of at least 1%. Notably, the assay was determined to be sufficiently sensitive even in patients harboring a low abundance of BCR-ABL1 levels. The PacBio sequencing successfully identified all mutations seen by standard methods. Importantly, we identified several mutations that escaped detection by the clinical routine analysis. Resistance mutations were found in all but one of the patients. Due to the long reads afforded by PacBio sequencing, compound mutations present in the same molecule were readily distinguished from independent alterations arising in different molecules. Moreover, several transcript isoforms of the BCR-ABL1 transcript were identified in two of the CML patients. Finally, our assay allowed for a quick turn around time allowing samples to be reported upon within 2 days.In summary the PacBio sequencing assay can be applied to detect BCR-ABL1 resistance mutations in both diagnostic and follow-up CML patient samples using a simple protocol applicable to routine diagnosis. The method besides its sensitivity, gives a complete view of the clonal distribution of mutations, which is of importance when making therapy decisions.


September 22, 2019  |  

Great differences in performance and outcome of high-throughput sequencing data analysis platforms for fungal metabarcoding.

Along with recent developments in high-throughput sequencing (HTS) technologies and thus fast accumulation of HTS data, there has been a growing need and interest for developing tools for HTS data processing and communication. In particular, a number of bioinformatics tools have been designed for analysing metabarcoding data, each with specific features, assumptions and outputs. To evaluate the potential effect of the application of different bioinformatics workflow on the results, we compared the performance of different analysis platforms on two contrasting high-throughput sequencing data sets. Our analysis revealed that the computation time, quality of error filtering and hence output of specific bioinformatics process largely depends on the platform used. Our results show that none of the bioinformatics workflows appears to perfectly filter out the accumulated errors and generate Operational Taxonomic Units, although PipeCraft, LotuS and PIPITS perform better than QIIME2 and Galaxy for the tested fungal amplicon dataset. We conclude that the output of each platform requires manual validation of the OTUs by examining the taxonomy assignment values.


September 22, 2019  |  

Soil microclimate changes affect soil fungal communities in a Mediterranean pine forest.

Soil microclimate is a potentially important regulator of the composition of plant-associated fungal communities in climates with significant drought periods. Here, we investigated the spatio-temporal dynamics of soil fungal communities in a Mediterranean Pinus pinaster forest in relation to soil moisture and temperature. Fungal communities in 336 soil samples collected monthly over 1 year from 28 long-term experimental plots were assessed by PacBio sequencing of ITS2 amplicons. Total fungal biomass was estimated by analysing ergosterol. Community changes were analysed in the context of functional traits. Soil fungal biomass was lowest during summer and late winter and highest during autumn, concurrent with a greater relative abundance of mycorrhizal species. Intra-annual spatio-temporal changes in community composition correlated significantly with soil moisture and temperature. Mycorrhizal fungi were less affected by summer drought than free-living fungi. In particular, mycorrhizal species of the short-distance exploration type increased in relative abundance under dry conditions, whereas species of the long-distance exploration type were more abundant under wetter conditions. Our observations demonstrate a potential for compositional and functional shifts in fungal communities in response to changing climatic conditions. Free-living fungi and mycorrhizal species with extensive mycelia may be negatively affected by increasing drought periods in Mediterranean forest ecosystems.© 2018 The Authors. New Phytologist © 2018 New Phytologist Trust.


September 22, 2019  |  

Human and rhesus macaque KIR haplotypes defined by their transcriptomes.

The killer-cell Ig-like receptors (KIRs) play a central role in the immune recognition in infection, pregnancy, and transplantation through their interactions with MHC class I molecules. KIR genes display abundant copy number variation as well as high levels of polymorphism. As a result, it is challenging to characterize this structurally dynamic region. KIR haplotypes have been analyzed in different species using conventional characterization methods, such as Sanger sequencing and Roche/454 pyrosequencing. However, these methods are time-consuming and often failed to define complete haplotypes, or do not reach allele-level resolution. In addition, most analyses were performed on genomic DNA, and thus were lacking substantial information about transcription and its corresponding modifications. In this paper, we present a single-molecule real-time sequencing approach, using Pacific Biosciences Sequel platform to characterize the KIR transcriptomes in human and rhesus macaque (Macaca mulatta) families. This high-resolution approach allowed the identification of novel Mamu-KIR alleles, the extension of reported allele sequences, and the determination of human and macaque KIR haplotypes. In addition, multiple recombinant KIR genes were discovered, all located on contracted haplotypes, which were likely the result of chromosomal rearrangements. The relatively high number of contracted haplotypes discovered might be indicative of selection on small KIR repertoires and/or novel fusion gene products. This next-generation method provides an improved high-resolution characterization of the KIR cluster in humans and macaques, which eventually may aid in a better understanding and interpretation of KIR allele-associated diseases, as well as the immune response in transplantation and reproduction. Copyright © 2018 by The American Association of Immunologists, Inc.


September 22, 2019  |  

Below-ground organic matter accumulation along a boreal forest fertility gradient relates to guild interaction within fungal communities.

Plant-soil interactions link ecosystem fertility and organic matter accumulation below ground. Soil microorganisms play a central role as mediators of these interactions, but mechanistic understanding is still largely lacking. Correlative data from a coniferous forest ecosystem support the hypothesis that interactions between fungal guilds play a central role in regulating organic matter accumulation in relation to fertility. With increasing ecosystem fertility, the proportion of saprotrophic basidiomycetes increased in deeper organic layers, at the expense of ectomycorrhizal fungal species. Saprotrophs correlated positively with the activity of oxidative enzymes, which in turn favoured organic matter turnover and nitrogen recycling to plants. Combined, our findings are consistent with a fungus-mediated feedback loop, which results in a negative correlation between ecosystem fertility and below-ground carbon storage. These findings call for a shift in focus from plant litter traits to fungal traits in explaining organic matter dynamics and ecosystem fertility in boreal forests.© 2017 John Wiley & Sons Ltd/CNRS.


September 22, 2019  |  

Fungal community shifts underpin declining mycelial production and turnover across a Pinus sylvestris chronosequence

Fungi play critical roles in ecosystem processes such as decomposition and nutrient cycling, but have also been highlighted as significant contributors to organic matter build-up in boreal forest soils. Ectomycorrhizal (ECM) mycelial biomass and necromass dynamics have recently been highlighted as essential for regulating build-up of soil organic matter. Understanding the extent to which shifts in mycelial community or growth trait composition cause changes in mycelial production and turnover over ecological gradients would aid a mechanistic understanding of these important processes at an ecosystem scale. Here, we test the hypotheses that shifting species and mycelial trait (exploration type) composition within the mycelial community underpin changes in biomass turnover with increasing forest age. We quantified mycelial turnover and assessed fungal community composition in a chronosequence of eight, 12- to 158-year-old, managed Pinus sylvestris forests. Turnover was estimated by determining mycelial biomass (ergosterol) in a sequence of ingrowth mesh bags and applying mathematical models. Fungal communities in the bags were identified using Pacific Biosciences sequencing of fungal ITS2 amplicons. To evaluate the accuracy of this method to represent all ECM fungi, community composition in bags was followed over time and compared with communities in soil. Mycelial communities changed with stand age, but we found no evidence that there were concurrent shifts in mycelial exploration types. Forest age and turnover were significantly correlated with ECM mycelial community composition and collectively explained 39.4% of total variation. The similarity between fungal communities in mesh bags and in soil was strongly forest age dependent, with communities in mesh bags diverging from soil communities in stands older than 60 years. However, in all stands, when bag incubation time exceeded 75 days, communities became more similar to soil communities. Synthesis. Our results support the idea that shifts in fungal community composition underpin the forest age-related decrease in mycelial turnover; however, since ingrowth mesh bags exclude some mycorrhizal species in older forests, it remains a possibility that turnover estimates were not reflecting the entire community. While we found no evidence that mycelial exploration types of fungi changed systematically with forest age, we suggest that other traits that relate to biomass turnover and necromass degradation require further study, as they may explain the extent to which ectomycorrhizal fungi regulate and contribute to soil organic matter accumulation.


September 22, 2019  |  

Retention of seed trees fails to lifeboat ectomycorrhizal fungal diversity in harvested Scots pine forests.

Fennoscandian forestry has in the past decades changed from natural regeneration of forests towards replantation of clear-cuts, which negatively impacts ectomycorrhizal fungal (EMF) diversity. Retention of trees during harvesting enables EMF survival, and we therefore expected EMF communities to be more similar to those in old natural stands after forest regeneration using seed trees compared to full clear-cutting and replanting. We sequenced fungal internal transcribed spacer 2 (ITS2) amplicons to assess EMF communities in 10- to 60-year-old Scots pine stands regenerated either using seed trees or through replanting of clear-cuts with old natural stands as reference. We also investigated local EMF communities around retained old trees. We found that retention of seed trees failed to mitigate the impact of harvesting on EMF community composition and diversity. With increasing stand age, EMF communities became increasingly similar to those in old natural stands and permanently retained trees maintained EMF locally. From our observations, we conclude that EMF communities, at least common species, post-harvest are more influenced by environmental filtering, resulting from environmental changes induced by harvest, than by the continuity of trees. These results suggest that retention of intact forest patches is a more efficient way to conserve EMF diversity than retaining dispersed single trees.© FEMS 2017. All rights reserved. For permissions, please e-mail: journals.permissions@oup.com.


September 22, 2019  |  

Improved metagenome assemblies and taxonomic binning using long-read circular consensus sequence data.

DNA assembly is a core methodological step in metagenomic pipelines used to study the structure and function within microbial communities. Here we investigate the utility of Pacific Biosciences long and high accuracy circular consensus sequencing (CCS) reads for metagenomic projects. We compared the application and performance of both PacBio CCS and Illumina HiSeq data with assembly and taxonomic binning algorithms using metagenomic samples representing a complex microbial community. Eight SMRT cells produced approximately 94 Mb of CCS reads from a biogas reactor microbiome sample that averaged 1319 nt in length and 99.7% accuracy. CCS data assembly generated a comparative number of large contigs greater than 1?kb, to those assembled from a ~190x larger HiSeq dataset (~18 Gb) produced from the same sample (i.e approximately 62% of total contigs). Hybrid assemblies using PacBio CCS and HiSeq contigs produced improvements in assembly statistics, including an increase in the average contig length and number of large contigs. The incorporation of CCS data produced significant enhancements in taxonomic binning and genome reconstruction of two dominant phylotypes, which assembled and binned poorly using HiSeq data alone. Collectively these results illustrate the value of PacBio CCS reads in certain metagenomics applications.


September 22, 2019  |  

Genome characterization of oleaginous Aspergillus oryzae BCC7051: A potential fungal-based platform for lipid production.

The selected robust fungus, Aspergillus oryzae strain BCC7051 is of interest for biotechnological production of lipid-derived products due to its capability to accumulate high amount of intracellular lipids using various sugars and agro-industrial substrates. Here, we report the genome sequence of the oleaginous A. oryzae BCC7051. The obtained reads were de novo assembled into 25 scaffolds spanning of 38,550,958 bps with predicted 11,456 protein-coding genes. By synteny mapping, a large rearrangement was found in two scaffolds of A. oryzae BCC7051 as compared to the reference RIB40 strain. The genetic relationship between BCC7051 and other strains of A. oryzae in terms of aflatoxin production was investigated, indicating that the A. oryzae BCC7051 was categorized into group 2 nonaflatoxin-producing strain. Moreover, a comparative analysis of the structural genes focusing on the involvement in lipid metabolism among oleaginous yeast and fungi revealed the presence of multiple isoforms of metabolic enzymes responsible for fatty acid synthesis in BCC7051. The alternative routes of acetyl-CoA generation as oleaginous features and malate/citrate/pyruvate shuttle were also identified in this A. oryzae strain. The genome sequence generated in this work is a dedicated resource for expanding genome-wide study of microbial lipids at systems level, and developing the fungal-based platform for production of diversified lipids with commercial relevance.


September 22, 2019  |  

Packaging of Dinoroseobacter shibae DNA into gene transfer agent particles is not random.

Gene transfer agents (GTAs) are phage-like particles which contain a fragment of genomic DNA of the bacterial or archaeal producer and deliver this to a recipient cell. GTA gene clusters are present in the genomes of almost all marine Rhodobacteraceae (Roseobacters) and might be important contributors to horizontal gene transfer in the world’s oceans. For all organisms studied so far, no obvious evidence of sequence specificity or other nonrandom process responsible for packaging genomic DNA into GTAs has been found. Here, we show that knock-out of an autoinducer synthase gene of Dinoroseobacter shibae resulted in overproduction and release of functional GTA particles (DsGTA). Next-generation sequencing of the 4.2-kb DNA fragments isolated from DsGTAs revealed that packaging was not random. DNA from low-GC conjugative plasmids but not from high-GC chromids was excluded from packaging. Seven chromosomal regions were strongly overrepresented in DNA isolated from DsGTA. These packaging peaks lacked identifiable conserved sequence motifs that might represent recognition sites for the GTA terminase complex. Low-GC regions of the chromosome, including the origin and terminus of replication, were underrepresented in DNA isolated from DsGTAs. DNA methylation reduced packaging frequency while the level of gene expression had no influence. Chromosomal regions found to be over- and underrepresented in DsGTA-DNA were regularly spaced. We propose that a “headful” type of packaging is initiated at the sites of coverage peaks and, after linearization of the chromosomal DNA, proceeds in both directions from the initiation site. GC-content, DNA-modifications, and chromatin structure might influence at which sides GTA packaging can be initiated.© The Author(s) 2018. Published by Oxford University Press on behalf of the Society for Molecular Biology and Evolution.


September 22, 2019  |  

Analysis of the hybrid genomes of two field isolates of the soil-borne fungal species Verticillium longisporum.

Brassica plant species are attacked by a number of pathogens; among them, the ones with a soil-borne lifestyle have become increasingly important. Verticillium stem stripe caused by Verticillium longisporum is one example. This fungal species is thought to be of a hybrid origin, having a genome composed of combinations of lineages denominated A and D. In this study we report the draft genomes of 2 V. longisporum field isolates sequenced using the Illumina technology. Genomic characterization and lineage composition, followed by selected gene analysis to facilitate the comprehension of its genomic features and potential effector categories were performed.The draft genomes of 2 Verticillium longisporum single spore isolates (VL1 and VL2) have an estimated ungapped size of about 70 Mb. The total number of protein encoding genes identified in VL1 was 20,793, whereas 21,072 gene models were predicted in VL2. The predicted genome size, gene contents, including the gene families coding for carbohydrate active enzymes were almost double the numbers found in V. dahliae and V. albo-atrum. Single nucleotide polymorphisms (SNPs) were frequently distributed in the two genomes but the distribution of heterozygosity and depth was not independent. Further analysis of potential parental lineages suggests that the V. longisporum genome is composed of two parts, A1 and D1, where A1 is more ancient than the parental lineage genome D1, the latter being more closer related to V. dahliae. Presence of the mating-type genes MAT1-1-1 and MAT1-2-1 in the V. longisporum genomes were confirmed. However, the MAT genes in V. dahliae, V. albo-atrum and V. longisporum have experienced extensive nucleotide changes at least partly explaining the present asexual nature of these fungal species.The established draft genome of V. longisporum is comparatively large compared to other studied ascomycete fungi. Consequently, high numbers of genes were predicted in the two V. longisporum genomes, among them many secreted proteins and carbohydrate active enzyme (CAZy) encoding genes. The genome is composed of two parts, where one lineage is more ancient than the part being more closely related to V. dahliae. Dissimilar mating-type sequences were identified indicating possible ancient hybridization events.


September 22, 2019  |  

Functional metagenomics reveals a novel carbapenem-hydrolyzing mobile beta-lactamase from Indian river sediments contaminated with antibiotic production waste.

Evolution has provided environmental bacteria with a plethora of genes that give resistance to antibiotic compounds. Under anthropogenic selection pressures, some of these genes are believed to be recruited over time into pathogens by horizontal gene transfer. River sediment polluted with fluoroquinolones and other drugs discharged from bulk drug production in India constitute an environment with unprecedented, long-term antibiotic selection pressures. It is therefore plausible that previously unknown resistance genes have evolved and/or are promoted here. In order to search for novel resistance genes, we therefore analyzed such river sediments by a functional metagenomics approach. DNA fragments providing resistance to different antibiotics in E. coli were sequenced using Sanger and PacBio RSII platforms. We recaptured the majority of known antibiotic resistance genes previously identified by open shot-gun metagenomics sequencing of the same samples. In addition, seven novel resistance gene candidates (six beta-lactamases and one amikacin resistance gene) were identified. Two class A beta-lactamases, blaRSA1 and blaRSA2, were phylogenetically close to clinically important ESBLs like blaGES, blaBEL and blaL2, and were further characterized for their substrate spectra. The blaRSA1 protein, encoded as an integron gene cassette, efficiently hydrolysed penicillins, first generation cephalosporins and cefotaxime, while blaRSA2 was an inducible class A beta-lactamase, capable of hydrolyzing carbapenems albeit with limited efficiency, similar to the L2 beta-lactamase from Stenotrophomonas maltophilia. All detected novel genes were associated with plasmid mobilization proteins, integrons, and/or other resistance genes, suggesting a potential for mobility. This study provides insight into a resistome shaped by an exceptionally strong and long-term antibiotic selection pressure. An improved knowledge of mobilized resistance factors in the external environment may make us better prepared for the resistance challenges that we may face in clinics in the future. Copyright © 2017 Elsevier Ltd. All rights reserved.


September 22, 2019  |  

Screening and genomic characterization of filamentous hemagglutinin-deficient Bordetella pertussis.

Despite high vaccine coverage, pertussis cases in the United States have increased over the last decade. Growing evidence suggests that disease resurgence results, in part, from genetic divergence of circulating strain populations away from vaccine references. The United States employs acellular vaccines exclusively, and current Bordetella pertussis isolates are predominantly deficient in at least one immunogen, pertactin (Prn). First detected in the United States retrospectively in a 1994 isolate, the rapid spread of Prn deficiency is likely vaccine driven, raising concerns about whether other acellular vaccine immunogens experience similar pressures, as further antigenic changes could potentially threaten vaccine efficacy. We developed an electrochemiluminescent antibody capture assay to monitor the production of the acellular vaccine immunogen filamentous hemagglutinin (Fha). Screening 722 U.S. surveillance isolates collected from 2010 to 2016 identified two that were both Prn and Fha deficient. Three additional Fha-deficient laboratory strains were also identified from a historic collection of 65 isolates dating back to 1935. Whole-genome sequencing of deficient isolates revealed putative, underlying genetic changes. Only four isolates harbored mutations to known genes involved in Fha production, highlighting the complexity of its regulation. The chromosomes of two Fha-deficient isolates included unexpected structural variation that did not appear to influence Fha production. Furthermore, insertion sequence disruption of fhaB was also detected in a previously identified pertussis toxin-deficient isolate that still produced normal levels of Fha. These results demonstrate the genetic potential for additional vaccine immunogen deficiency and underscore the importance of continued surveillance of circulating B. pertussis evolution in response to vaccine pressure. Copyright © 2018 American Society for Microbiology.


September 22, 2019  |  

Characterization of plasmids harboring blaCTX-M and blaCMY genes in E. coli from French broilers.

Resistance to extended-spectrum cephalosporins (ESC) is a global health issue. The aim of this study was to analyze and compare plasmids coding for resistance to ESC isolated from 16 avian commensal and 17 avian pathogenic Escherichia coli (APEC) strains obtained respectively at slaughterhouse or from diseased broilers in 2010-2012. Plasmid DNA was used to transform E. coli DH5alpha, and the resistances of the transformants were determined. The sequences of the ESC-resistance plasmids prepared from transformants were obtained by Illumina (33 plasmids) or PacBio (1 plasmid). Results showed that 29 of these plasmids contained the blaCTX-M-1 gene and belonged to the IncI1/ST3 type, with 27 and 20 of them carrying the sul2 or tet(A) genes respectively. Despite their diverse origins, several plasmids showed very high percentages of identity. None of the blaCTX-M-1-containing plasmid contained APEC virulence genes, although some of them were detected in the parental strains. Three plasmids had the blaCMY-2 gene, but no other resistance gene. They belonged to IncB/O/K/Z-like or IncFIA/FIB replicon types. The blaCMY-2 IncFIA/FIB plasmid was obtained from a strain isolated from a diseased broiler and also containing a blaCTX-M-1 IncI1/ST3 plasmid. Importantly APEC virulence genes (sitA-D, iucA-D, iutA, hlyF, ompT, etsA-C, iss, iroB-E, iroN, cvaA-C and cvi) were detected on the blaCMY-2 plasmid. In conclusion, our results show the dominance and high similarity of blaCTX-M-1 IncI1/ST3 plasmids, and the worrying presence of APEC virulence genes on a blaCMY-2 plasmid.


Talk with an expert

If you have a question, need to check the status of an order, or are interested in purchasing an instrument, we're here to help.