Mutational and transcriptional landscape of pediatric B-cell precursor lymphoblastic lymphoma.

ABSTRACT: Pediatric B-cell precursor (BCP) lymphoblastic malignancies are neoplasms with manifestation either in the bone marrow or blood (BCP acute lymphoblastic leukemia [BCP-ALL]) or are less common in extramedullary tissue (BCP lymphoblastic lymphoma [BCP-LBL]). Although both presentations are similar in morphology and immunophenotype, molecular studies have been virtually restricted to BCP-ALL so far. The lack of molecular studies on BCP-LBL is due to its rarity and restriction on small, mostly formalin-fixed paraffin-embedded (FFPE) tissues. Here, to our knowledge, we present the first comprehensive mutational and transcriptional analysis of what we consider the largest BCP-LBL cohort described to date (n = 97). Whole-exome sequencing indicated a mutational spectrum of BCP-LBL, strikingly similar to that found in BCP-ALL. However, epigenetic modifiers were more frequently mutated in BCP-LBL, whereas BCP-ALL was more frequently affected by mutation in genes involved in B-cell development. Integrating copy number alterations, somatic mutations, and gene expression by RNA sequencing revealed that virtually all molecular subtypes originally defined in BCP-ALL are present in BCP-LBL, with only 7% of lymphomas that were not assigned to a subtype. Similar to BCP-ALL, the most frequent subtypes of BCP-LBL were high hyperdiploidy and ETV6::RUNX1. Tyrosine kinase/cytokine receptor rearrangements were detected in 7% of BCP-LBL. These results indicate that genetic subtypes can be identified in BCP-LBL using next-generation sequencing, even in FFPE tissue, and may be relevant to guide treatment.

Phenotypic and genomic dissection of colour pattern variation in a reef fish radiation.

Coral reefs rank among the most diverse species assemblages on Earth. A particularly striking aspect of coral reef communities is the variety of colour patterns displayed by reef fishes. Colour pattern is known to play a central role in the ecology and evolution of reef fishes through, for example, signalling or camouflage. Nevertheless, colour pattern is a complex trait in reef fishes-actually a collection of traits-that is difficult to analyse in a quantitative and standardized way. This is the challenge that we address in this study using the hamlets (Hypoplectrus spp., Serranidae) as a model system. Our approach involves a custom underwater camera system to take orientation- and size-standardized photographs in situ, colour correction, alignment of the fish images with a combination of landmarks and Bézier curves, and principal component analysis on the colour value of each pixel of each aligned fish. This approach identifies the major colour pattern elements that contribute to phenotypic variation in the group. Furthermore, we complement the image analysis with whole-genome sequencing to run a multivariate genome-wide association study for colour pattern variation. This second layer of analysis reveals sharp association peaks along the hamlet genome for each colour pattern element and allows to characterize the phenotypic effect of the single nucleotide polymorphisms that are most strongly associated with colour pattern variation at each association peak. Our results suggest that the diversity of colour patterns displayed by the hamlets is generated by a modular genomic and phenotypic architecture.

Two different epigenetic information channels in wild three-spined sticklebacks are involved in salinity adaptation.

Epigenetic inheritance has been proposed to contribute to adaptation and acclimation via two information channels: (i) inducible epigenetic marks that enable transgenerational plasticity and (ii) noninducible epigenetic marks resulting from random epimutations shaped by selection. We studied both postulated channels by sequencing methylomes and genomes of Baltic three-spined sticklebacks (Gasterosteus aculeatus) along a salinity cline. Wild populations differing in salinity tolerance revealed differential methylation (pop-DMS) at genes enriched for osmoregulatory processes. A two-generation experiment demonstrated that 62% of these pop-DMS were noninducible by salinity manipulation, suggesting that they are the result of either direct selection or associated genomic divergence at cis- or trans-regulatory sites. Two-thirds of the remaining inducible pop-DMS increased in similarity to patterns detected in wild populations from corresponding salinities. The level of similarity accentuated over consecutive generations, indicating a mechanism of transgenerational plasticity. While we can attribute natural DNA methylation patterns to the two information channels, their interplay with genomic variation in salinity adaptation is still unresolved.

The evolution of microendemism in a reef fish (Hypoplectrus maya).

Marine species tend to have extensive distributions, which are commonly attributed to the dispersal potential provided by planktonic larvae and the rarity of absolute barriers to dispersal in the ocean. Under this paradigm, the occurrence of marine microendemism without geographic isolation in species with planktonic larvae poses a dilemma. The recently described Maya hamlet (Hypoplectrus maya, Serranidae) is exactly such a case, being endemic to a 50-km segment of the Mesoamerican Barrier Reef System (MBRS). We use whole-genome analysis to infer the demographic history of the Maya hamlet and contrast it with the sympatric and pan-Caribbean black (H. nigricans), barred (H. puella) and butter (H. unicolor) hamlets, as well as the allopatric but phenotypically similar blue hamlet (H. gemma). We show that H. maya is indeed a distinct evolutionary lineage, with genomic signatures of inbreeding and a unique demographic history of continuous decrease in effective population size since it diverged from congeners just ~3,000 generations ago. We suggest that this case of microendemism may be driven by the combination of a narrow ecological niche and restrictive oceanographic conditions in the southern MBRS, which is consistent with the occurrence of an unusually high number of marine microendemics in this region. The restricted distribution of the Maya hamlet, its decline in both census and effective population sizes, and the degradation of its habitat place it at risk of extinction. We conclude that the evolution of marine microendemism can be a fast and dynamic process, with extinction possibly occurring before speciation is complete.

Inter-chromosomal coupling between vision and pigmentation genes during genomic divergence.

Recombination between loci underlying mate choice and ecological traits is a major evolutionary force acting against speciation with gene flow. The evolution of linkage disequilibrium between such loci is therefore a fundamental step in the origin of species. Here, we show that this process can take place in the absence of physical linkage in hamlets-a group of closely related reef fishes from the wider Caribbean that differ essentially in colour pattern and are reproductively isolated through strong visually-based assortative mating. Using full-genome analysis, we identify four narrow genomic intervals that are consistently differentiated among sympatric species in a backdrop of extremely low genomic divergence. These four intervals include genes involved in pigmentation (sox10), axial patterning (hoxc13a), photoreceptor development (casz1) and visual sensitivity (SWS and LWS opsins) that develop islands of long-distance and inter-chromosomal linkage disequilibrium as species diverge. The relatively simple genomic architecture of species differences facilitates the evolution of linkage disequilibrium in the presence of gene flow.

Serendipitous Meta-Transcriptomics: The Fungal Community of Norway Spruce (Picea abies).

After performing de novo transcript assembly of >1 billion RNA-Sequencing reads obtained from 22 samples of different Norway spruce (Picea abies) tissues that were not surface sterilized, we found that assembled sequences captured a mix of plant, lichen, and fungal transcripts. The latter were likely expressed by endophytic and epiphytic symbionts, indicating that these organisms were present, alive, and metabolically active. Here, we show that these serendipitously sequenced transcripts need not be considered merely as contamination, as is common, but that they provide insight into the plant's phyllosphere. Notably, we could classify these transcripts as originating predominantly from Dothideomycetes and Leotiomycetes species, with functional annotation of gene families indicating active growth and metabolism, with particular regards to glucose intake and processing, as well as gene regulation.

A universal genomic coordinate translator for comparative genomics.

BACKGROUND: Genomic duplications constitute major events in the evolution of species, allowing paralogous copies of genes to take on fine-tuned biological roles. Unambiguously identifying the orthology relationship between copies across multiple genomes can be resolved by synteny, i.e. the conserved order of genomic sequences. However, a comprehensive analysis of duplication events and their contributions to evolution would require all-to-all genome alignments, which increases at N2 with the number of available genomes, N. RESULTS: Here, we introduce Kraken, software that omits the all-to-all requirement by recursively traversing a graph of pairwise alignments and dynamically re-computing orthology. Kraken scales linearly with the number of targeted genomes, N, which allows for including large numbers of genomes in analyses. We first evaluated the method on the set of 12 Drosophila genomes, finding that orthologous correspondence computed indirectly through a graph of multiple synteny maps comes at minimal cost in terms of sensitivity, but reduces overall computational runtime by an order of magnitude. We then used the method on three well-annotated mammalian genomes, human, mouse, and rat, and show that up to 93% of protein coding transcripts have unambiguous pairwise orthologous relationships across the genomes. On a nucleotide level, 70 to 83% of exons match exactly at both splice junctions, and up to 97% on at least one junction. We last applied Kraken to an RNA-sequencing dataset from multiple vertebrates and diverse tissues, where we confirmed that brain-specific gene family members, i.e. one-to-many or many-to-many homologs, are more highly correlated across species than single-copy (i.e. one-to-one homologous) genes. Not limited to protein coding genes, Kraken also identifies thousands of newly identified transcribed loci, likely non-coding RNAs that are consistently transcribed in human, chimpanzee and gorilla, and maintain significant correlation of expression levels across species. CONCLUSIONS: Kraken is a computational genome coordinate translator that facilitates cross-species comparisons, distinguishes orthologs from paralogs, and does not require costly all-to-all whole genome mappings. Kraken is freely available under LPGL from http://github.com/nedaz/kraken.

Modular and configurable optimal sequence alignment software: Cola.

BACKGROUND: The fundamental challenge in optimally aligning homologous sequences is to define a scoring scheme that best reflects the underlying biological processes. Maximising the overall number of matches in the alignment does not always reflect the patterns by which nucleotides mutate. Efficiently implemented algorithms that can be parameterised to accommodate more complex non-linear scoring schemes are thus desirable. RESULTS: We present Cola, alignment software that implements different optimal alignment algorithms, also allowing for scoring contiguous matches of nucleotides in a nonlinear manner. The latter places more emphasis on short, highly conserved motifs, and less on the surrounding nucleotides, which can be more diverged. To illustrate the differences, we report results from aligning 14,100 sequences from 3' untranslated regions of human genes to 25 of their mammalian counterparts, where we found that a nonlinear scoring scheme is more consistent than a linear scheme in detecting short, conserved motifs. CONCLUSIONS: Cola is freely available under LPGL from https://github.com/nedaz/cola.

Algal genomes reveal evolutionary mosaicism and the fate of nucleomorphs.

Cryptophyte and chlorarachniophyte algae are transitional forms in the widespread secondary endosymbiotic acquisition of photosynthesis by engulfment of eukaryotic algae. Unlike most secondary plastid-bearing algae, miniaturized versions of the endosymbiont nuclei (nucleomorphs) persist in cryptophytes and chlorarachniophytes. To determine why, and to address other fundamental questions about eukaryote-eukaryote endosymbiosis, we sequenced the nuclear genomes of the cryptophyte Guillardia theta and the chlorarachniophyte Bigelowiella natans. Both genomes have >21,000 protein genes and are intron rich, and B. natans exhibits unprecedented alternative splicing for a single-celled organism. Phylogenomic analyses and subcellular targeting predictions reveal extensive genetic and biochemical mosaicism, with both host- and endosymbiont-derived genes servicing the mitochondrion, the host cell cytosol, the plastid and the remnant endosymbiont cytosol of both algae. Mitochondrion-to-nucleus gene transfer still occurs in both organisms but plastid-to-nucleus and nucleomorph-to-nucleus transfers do not, which explains why a small residue of essential genes remains locked in each nucleomorph.

Artificially designed promoters: understanding the role of spatial features and canonical binding sites in transcription.

The promoter is a key element in gene transcription and regulation. We previously reported that artificial sequences rich in the dinucleotide CpG are sufficient to drive expression in vitro in mammalian cell lines, without requiring canonical binding sites for transcription factor proteins. Here, we report that introducing a promoter organization that alternates in CpGs and regions rich in A and T further increases expression strength, as well as how insertion of specific binding sites makes such sequences respond to induced levels of the transcription factor NFκB. Our findings further contribute to the mechanistic understanding of promoters, as well as how these sequences might be shaped by evolutionary pressure in living organisms.