Genomic, transcriptomic and epigenomic sequencing data of the B-cell leukemia cell line REH.

OBJECTIVES: The aim of this data paper is to describe a collection of 33 genomic, transcriptomic and epigenomic sequencing datasets of the B-cell acute lymphoblastic leukemia (ALL) cell line REH. REH is one of the most frequently used cell lines for functional studies of pediatric ALL, and these data provide a multi-faceted characterization of its molecular features. The datasets described herein, generated with short- and long-read sequencing technologies, can both provide insights into the complex aberrant karyotype of REH, and be used as reference datasets for sequencing data quality assessment or for methods development. DATA DESCRIPTION: This paper describes 33 datasets corresponding to 867 gigabases of raw sequencing data generated from the REH cell line. These datasets include five different approaches for whole genome sequencing (WGS) on four sequencing platforms, two RNA sequencing (RNA-seq) techniques on two different sequencing platforms, DNA methylation sequencing, and single-cell ATAC-sequencing.

Inversions maintain differences between migratory phenotypes of a songbird.

Structural rearrangements have been shown to be important in local adaptation and speciation, but have been difficult to reliably identify and characterize in non-model species. Here we combine long reads, linked reads and optical mapping to characterize three divergent chromosome regions in the willow warbler Phylloscopus trochilus, of which two are associated with differences in migration and one with an environmental gradient. We show that there are inversions (0.4-13 Mb) in each of the regions and that the divergence times between inverted and non-inverted haplotypes are similar across the regions (~1.2 Myrs), which is compatible with a scenario where inversions arose in either of two allopatric populations that subsequently hybridized. The improved genomes allow us to detect additional functional differences in the divergent regions, providing candidate genes for migration and adaptations to environmental gradients.

Hybrid sequencing resolves two germline ultra-complex chromosomal rearrangements consisting of 137 breakpoint junctions in a single carrier.

Chromoanagenesis is a genomic event responsible for the formation of complex structural chromosomal rearrangements (CCRs). Germline chromoanagenesis is rare and the majority of reported cases are associated with an affected phenotype. Here, we report a healthy female carrying two de novo CCRs involving chromosomes 4, 19, 21 and X and chromosomes 7 and 11, respectively, with a total of 137 breakpoint junctions (BPJs). We characterized the CCRs using a hybrid-sequencing approach, combining short-read sequencing, nanopore sequencing, and optical mapping. The results were validated using multiple cytogenetic methods, including fluorescence in situ hybridization, spectral karyotyping, and Sanger sequencing. We identified 137 BPJs, which to our knowledge is the highest number of reported breakpoint junctions in germline chromoanagenesis. We also performed a statistical assessment of the positioning of the breakpoints, revealing a significant enrichment of BPJ-affecting genes (96 intragenic BPJs, 26 genes, p < 0.0001), indicating that the CCRs formed during active transcription of these genes. In addition, we find that the DNA fragments are unevenly and non-randomly distributed across the derivative chromosomes indicating a multistep process of scattering and re-joining of DNA fragments. In summary, we report a new maximum number of BPJs (137) in germline chromoanagenesis. We also show that a hybrid sequencing approach is necessary for the correct characterization of complex CCRs. Through in-depth statistical assessment, it was found that the CCRs most likely was formed through an event resembling chromoplexy-a catastrophic event caused by erroneous transcription factor binding.

Evaluation of Single-Molecule Sequencing Technologies for Structural Variant Detection in Two Swedish Human Genomes.

Long-read single molecule sequencing is increasingly used in human genomics research, as it allows to accurately detect large-scale DNA rearrangements such as structural variations (SVs) at high resolution. However, few studies have evaluated the performance of different single molecule sequencing platforms for SV detection in human samples. Here we performed Oxford Nanopore Technologies (ONT) whole-genome sequencing of two Swedish human samples (average 32× coverage) and compared the results to previously generated Pacific Biosciences (PacBio) data for the same individuals (average 66× coverage). Our analysis inferred an average of 17k and 23k SVs from the ONT and PacBio data, respectively, with a majority of them overlapping with an available multi-platform SV dataset. When comparing the SV calls in the two Swedish individuals, we find a higher concordance between ONT and PacBio SVs detected in the same individual as compared to SVs detected by the same technology in different individuals. Downsampling of PacBio reads, performed to obtain similar coverage levels for all datasets, resulted in 17k SVs per individual and improved overlap with the ONT SVs. Our results suggest that ONT and PacBio have a similar performance for SV detection in human whole genome sequencing data, and that both technologies are feasible for population-scale studies.

The yeast Wickerhamomyces anomalus AS1 secretes a multifunctional exo-ß-1,3-glucanase with implications for winemaking.

A multifunctional exo-ß-1,3-glucanase (WaExg2) was purified from the culture supernatant of the yeast Wickerhamomyces anomalus AS1. The enzyme was identified by mass spectroscopic analysis of tryptic peptide fragments and the encoding gene WaEXG2 was sequenced. The latter codes for a protein of 427 amino acids, beginning with a probable signal peptide (17 aa) for secretion. The mature protein has a molecular mass of 47 456 Da with a calculated pI of 4.84. The somewhat higher mass of the protein in SDS-PAGE might be due to bound carbohydrates. Presumptive disulphide bridges confer a high compactness to the molecule. This explains the apparent smaller molecular mass (35 kDa) of the native enzyme determined by electrophoresis, whereas the unfolded form is consistent with the theoretical mass. Enzymatic hydrolysis of selected glycosides and glycans by WaExg2 was proved by TLC analysis of cleavage products. Glucose was detected as the sole hydrolysis product from laminarin, underlining that the enzyme acts as an exoglucanase. In addition, the enzyme efficiently hydrolysed small ß-linked glycosides (arbutin, esculin, polydatin, salicin) and disaccharides (cellobiose, gentiobiose). WaExg2 was active under typical wine-related conditions, such as low pH (3.5-4.0), high sugar concentrations (up to 20% w/v), high ethanol concentrations (10-15% v/v), presence of sulphites (up to 2 mm) and various cations. Therefore, the characterized enzyme might have multiple uses in winemaking, to increase concentrations of sensory and bioactive compounds by splitting glycosylated precursors or to reduce viscosity by hydrolysis of glycan slimes.

Introgressive hybridization between anciently diverged lineages of Silene (Caryophyllaceae).

Hybridization has played a major role during the evolution of angiosperms, mediating both gene flow between already distinct species and the formation of new species. Newly formed hybrids between distantly related taxa are often sterile. For this reason, interspecific crosses resulting in fertile hybrids have rarely been described to take place after more than a few million years after divergence. We describe here the traces of a reproductively successful hybrid between two ancestral species of Silene, diverged for about six million years prior to hybridization. No extant hybrids between the two parental lineages are currently known, but introgression of the RNA polymerase gene NRPA2 provides clear evidence of a temporary and fertile hybrid. Parsimony reconciliation between gene trees and the species tree, as well as consideration of clade ages, help exclude gene paralogy and lineage sorting as alternative hypotheses. This may represent one of the most extreme cases of divergence between species prior to introgressive hybridization discovered yet, notably at a homoploid level. Although species boundaries are generally believed to be stable after millions of years of divergence, we believe that this finding may indicate that gene flow between distantly related species is merely largely undetected at present.

Use of allele-specific sequencing primers is an efficient alternative to PCR subcloning of low-copy nuclear genes.

Direct Sanger sequencing of polymerase chain reaction (PCR)-amplified nuclear genes leads to polymorphic sequences when allelic variation is present. To overcome this problem, most researchers subclone the PCR products to separate alleles. An alternative is to directly sequence the separate alleles using allele-specific primers. We tested two methods to enhance the specificity of allele-specific primers for use in direct sequencing: using short primers and amplification refractory mutation system (ARMS) technique. By shortening the allele-specific primer to 15-13 nucleotides, the single mismatch in the ultimate base of the primer is enough to hinder the amplification of the nontarget allele in direct sequencing and recover only the targeted allele at high accuracy. The deliberate addition of a second mismatch, as implemented in the ARMS technique, was less successful and seems better suited for allele-specific amplification in regular PCR rather than in direct sequencing.

Metaxa: a software tool for automated detection and discrimination among ribosomal small subunit (12S/16S/18S) sequences of archaea, bacteria, eukaryotes, mitochondria, and chloroplasts in metagenomes and environmental sequencing datasets.

The ribosomal small subunit (SSU) rRNA gene has emerged as an important genetic marker for taxonomic identification in environmental sequencing datasets. In addition to being present in the nucleus of eukaryotes and the core genome of prokaryotes, the gene is also found in the mitochondria of eukaryotes and in the chloroplasts of photosynthetic eukaryotes. These three sets of genes are conceptually paralogous and should in most situations not be aligned and analyzed jointly. To identify the origin of SSU sequences in complex sequence datasets has hitherto been a time-consuming and largely manual undertaking. However, the present study introduces Metaxa ( http://microbiology.se/software/metaxa/ ), an automated software tool to extract full-length and partial SSU sequences from larger sequence datasets and assign them to an archaeal, bacterial, nuclear eukaryote, mitochondrial, or chloroplast origin. Using data from reference databases and from full-length organelle and organism genomes, we show that Metaxa detects and scores SSU sequences for origin with very low proportions of false positives and negatives. We believe that this tool will be useful in microbial and evolutionary ecology as well as in metagenomics.

Molecular and biochemical properties of the S-layer protein from the wine bacterium Lactobacillus hilgardii B706.

Different strains of the genus Lactobacillus can be regularly isolated from must and wine samples. By various physiological activities, they can improve or reduce the wine quality. Lactobacillus hilgardii that is known to survive under harsh wine conditions is classified as a spoilage bacterium, e.g. due to the production of histamine. Many lactobacilli form an S-layer as the outermost cell wall component which has been found to facilitate the colonization of special ecological niches. A detailed understanding of the properties related to their S-layer proteins is necessary to improve the knowledge of the interactions between different bacterial cells and with the surrounding environments. The S-layer protein from the wine-related L. hilgardii strain B706 has been isolated and its gene sequence determined. The deduced amino acid sequence corresponds to a 41 kDa protein with an isoelectric point of 9.6 without additional posttranslational modifications after splitting off the leader peptide. The complete protein is organized in a 32 amino acids signal sequence for membrane translocation, a positively charged N-terminal domain that binds to the cell wall and a negatively charged C-terminal domain. When the S-layer was removed, the corresponding L. hilgardii B706 cells became more sensitive to bacteriolytic enzymes and some wine-related stress conditions. From a practical point of view, the S-layer may be considered as a target for the inhibition of food-spoiling lactobacilli.

Study of the distribution of autotrophic CO2 fixation cycles in Crenarchaeota.

Two new autotrophic carbon fixation cycles have been recently described in Crenarchaeota. The 3-hydroxypropionate/4-hydroxybutyrate cycle using acetyl-coenzyme A (CoA)/propionyl-CoA carboxylase as the carboxylating enzyme has been identified for (micro)aerobic members of the Sulfolobales. The dicarboxylate/4-hydroxybutyrate cycle using oxygen-sensitive pyruvate synthase and phosphoenolpyruvate carboxylase as carboxylating enzymes has been found in members of the anaerobic Desulfurococcales and Thermoproteales. However, Sulfolobales include anaerobic and Desulfurococcales aerobic autotrophic representatives, raising the question of which of the two cycles they use. We studied the mechanisms of autotrophic CO(2) fixation in the strictly anaerobic Stygiolobus azoricus (Sulfolobales) and in the facultatively aerobic Pyrolobus fumarii (Desulfurococcales). The activities of all enzymes of the 3-hydroxypropionate/4-hydroxybutyrate cycle were found in the anaerobic S. azoricus. In contrast, the aerobic or denitrifying P. fumarii possesses all enzyme activities of the dicarboxylate/4-hydroxybutyrate cycle. We conclude that autotrophic Crenarchaeota use one of the two cycles, and that their distribution correlates with the 16S rRNA-based phylogeny of this group, rather than with the aerobic or anaerobic lifestyle.

Inferring polyploid phylogenies from multiply-labeled gene trees.

BACKGROUND: Gene trees that arise in the context of reconstructing the evolutionary history of polyploid species are often multiply-labeled, that is, the same leaf label can occur several times in a single tree. This property considerably complicates the task of forming a consensus of a collection of such trees compared to usual phylogenetic trees. RESULTS: We present a method for computing a consensus tree of multiply-labeled trees. As with the well-known greedy consensus tree approach for phylogenetic trees, our method first breaks the given collection of gene trees into a set of clusters. It then aims to insert these clusters one at a time into a tree, starting with the clusters that are supported by most of the gene trees. As the problem to decide whether a cluster can be inserted into a multiply-labeled tree is computationally hard, we have developed a heuristic method for solving this problem. CONCLUSION: We illustrate the applicability of our method using two collections of trees for plants of the genus Silene, that involve several allopolyploids at different levels.