Reconstructing cancer genomes from paired-end sequencing data

<p>Abstract</p> <p>Background</p> <p>A cancer genome is derived from the germline genome through a series of somatic mutations. Somatic structural variants - including duplications, deletions, inversions, translocations, and other rearrangements - result in a cancer genome that is a scrambling of intervals, or "blocks" of the germline genome sequence. We present an efficient algorithm for reconstructing the block organization of a cancer genome from paired-end DNA sequencing data.</p> <p>Results</p> <p>By aligning paired reads from a cancer genome - and a matched germline genome, if available - to the human reference genome, we derive: (i) a partition of the reference genome into intervals; (ii) adjacencies between these intervals in the cancer genome; (iii) an estimated copy number for each interval. We formulate the Copy Number and Adjacency Genome Reconstruction Problem of determining the cancer genome as a sequence of the derived intervals that is consistent with the measured adjacencies and copy numbers. We design an efficient algorithm, called Paired-end Reconstruction of Genome Organization (PREGO), to solve this problem by reducing it to an optimization problem on an interval-adjacency graph constructed from the data. The solution to the optimization problem results in an Eulerian graph, containing an alternating Eulerian tour that corresponds to a cancer genome that is consistent with the sequencing data. We apply our algorithm to five ovarian cancer genomes that were sequenced as part of The Cancer Genome Atlas. We identify numerous rearrangements, or structural variants, in these genomes, analyze reciprocal vs. non-reciprocal rearrangements, and identify rearrangements consistent with known mechanisms of duplication such as tandem duplications and breakage/fusion/bridge (B/F/B) cycles.</p> <p>Conclusions</p> <p>We demonstrate that PREGO efficiently identifies complex and biologically relevant rearrangements in cancer genome sequencing data. An implementation of the PREGO algorithm is available at <url>http://compbio.cs.brown.edu/software/</url>.</p

Site Specific Modification of Adeno-Associated Virus Enables Both Fluorescent Imaging of Viral Particles and Characterization of the Capsid Interactome

Adeno-associated viruses (AAVs) are attractive gene therapy vectors due to their low toxicity, high stability, and rare integration into the host genome. Expressing ligands on the viral capsid can re-target AAVs to new cell types, but limited sites have been identified on the capsid that tolerate a peptide insertion. Here, we incorporated a site-specific tetracysteine sequence into the AAV serotype 9 (AAV9) capsid, to permit labelling of viral particles with either a fluorescent dye or biotin. We demonstrate that fluorescently labelled particles are detectable in vitro, and explore the utility of the method in vivo in mice with time-lapse imaging. We exploit the biotinylated viral particles to generate two distinct AAV interactomes, and identify several functional classes of proteins that are highly represented: actin/cytoskeletal protein binding, RNA binding, RNA splicing/processing, chromatin modifying, intracellular trafficking and RNA transport proteins. To examine the biological relevance of the capsid interactome, we modulated the expression of two proteins from the interactomes prior to AAV transduction. Blocking integrin αVβ6 receptor function reduced AAV9 transduction, while reducing histone deacetylase 4 (HDAC4) expression enhanced AAV transduction. Our method demonstrates a strategy for inserting motifs into the AAV capsid without compromising viral titer or infectivity

Evolution of a supergene that regulates a trans-species social polymorphism

Supergenes are clusters of linked genetic loci that jointly affect the expression of complex phenotypes, such as social organization. Little is known about the origin and evolution of these intriguing genomic elements. Here we analyse whole-genome sequences of males from native populations of six fire ant species and show that variation in social organization is under the control of a novel supergene haplotype (termed Sb), which evolved by sequential incorporation of three inversions spanning half of a 'social chromosome'. Two of the inversions interrupt protein-coding genes, resulting in the increased expression of one gene and modest truncation in the primary protein structure of another. All six socially polymorphic species studied harbour the same three inversions, with the single origin of the supergene in their common ancestor inferred by phylogenomic analyses to have occurred half a million years ago. The persistence of Sb along with the ancestral SB haplotype through multiple speciation events provides a striking example of a functionally important trans-species social polymorphism presumably maintained by balancing selection. We found that while recombination between the Sb and SB haplotypes is severely restricted in all species, a low level of gene flux between the haplotypes has occurred following the appearance of the inversions, potentially mitigating the evolutionary degeneration expected at genomic regions that cannot freely recombine. These results provide a detailed picture of the structural genomic innovations involved in the formation of a supergene controlling a complex social phenotype

Comparative phylogeography in the Atlantic forest and Brazilian savannas: pleistocene fluctuations and dispersal shape spatial patterns in two bumblebees

Background: Bombus morio and B. pauloensis are sympatric widespread bumblebee species that occupy two major Brazilian biomes, the Atlantic forest and the savannas of the Cerrado. Differences in dispersion capacity, which is greater in B. morio, likely influence their phylogeographic patterns. This study asks which processes best explain the patterns of genetic variation observed in B. morio and B. pauloensis, shedding light on the phenomena that shaped the range of local populations and the spatial distribution of intra-specific lineages. Results: Results suggest that Pleistocene climatic oscillations directly influenced the population structure of both species. Correlative species distribution models predict that the warmer conditions of the Last Interglacial contributed to population contraction, while demographic expansion happened during the Last Glacial Maximum. These results are consistent with physiological data suggesting that bumblebees are well adapted to colder conditions. Intra-specific mitochondrial genealogies are not congruent between the two species, which may be explained by their documented differences in dispersal ability. Conclusions: While populations of the high-dispersal B. morio are morphologically and genetically homogeneous across the species range, B. pauloensis encompasses multiple (three) mitochondrial lineages, and show clear genetic, geographic, and morphological differences. Because the lineages of B. pauloensis are currently exposed to distinct climatic conditions (and elevations), parapatric diversification may occur within this taxon. The eastern portion of the state of São Paulo, the most urbanized area in Brazil, represents the center of genetic diversity for B. pauloensis

Keratan sulphate in the tumour environment

Keratan sulphate (KS) is a bioactive glycosaminoglycan (GAG) of some complexity composed of the repeat disaccharide D-galactose β1→4 glycosidically linked to N-acetyl glucosamine. During the biosynthesis of KS, a family of glycosyltransferase and sulphotransferase enzymes act sequentially and in a coordinated fashion to add D-galactose (D-Gal) then N-acetyl glucosamine (GlcNAc) to a GlcNAc acceptor residue at the reducing terminus of a nascent KS chain to effect chain elongation. D-Gal and GlcNAc can both undergo sulphation at C6 but this occurs more frequently on GlcNAc than D-Gal. Sulphation along the developing KS chain is not uniform and contains regions of variable length where no sulphation occurs, regions which are monosulphated mainly on GlcNAc and further regions of high sulphation where both of the repeat disaccharides are sulphated. Each of these respective regions in the KS chain can be of variable length leading to KS complexity in terms of chain length and charge localization along the KS chain. Like other GAGs, it is these variably sulphated regions in KS which define its interactive properties with ligands such as growth factors, morphogens and cytokines and which determine the functional properties of tissues containing KS. Further adding to KS complexity is the identification of three different linkage structures in KS to asparagine (N-linked) or to threonine or serine residues (O-linked) in proteoglycan core proteins which has allowed the categorization of KS into three types, namely KS-I (corneal KS, N-linked), KS-II (skeletal KS, O-linked) or KS-III (brain KS, O-linked). KS-I to -III are also subject to variable addition of L-fucose and sialic acid groups. Furthermore, the GlcNAc residues of some members of the mucin-like glycoprotein family can also act as acceptor molecules for the addition of D-Gal and GlcNAc residues which can also be sulphated leading to small low sulphation glycoforms of KS. These differ from the more heavily sulphated KS chains found on proteoglycans. Like other GAGs, KS has evolved molecular recognition and information transfer properties over hundreds of millions of years of vertebrate and invertebrate evolution which equips them with cell mediatory properties in normal cellular processes and in aberrant pathological situations such as in tumourogenesis. Two KS-proteoglycans in particular, podocalyxin and lumican, are cell membrane, intracellular or stromal tissue–associated components with roles in the promotion or regulation of tumour development, mucin-like KS glycoproteins may also contribute to tumourogenesis. A greater understanding of the biology of KS may allow better methodology to be developed to more effectively combat tumourogenic processes

Apes have culture but may not know that they do

The research leading to these results has received funding from the People Programme (Marie Curie Actions) and from the European Research Council under the European Union’s Seventh Framework Programme for research, technological development and demonstration under REA grant agreement n° 329197 awarded to Thibaud Gruber, ERC grant agreement n° 283871 awarded to Klaus Zuberbühler. Date of Acceptance: 16/01/2015There is good evidence that some ape behaviors can be transmitted socially and that this can lead to group-specific traditions. However, many consider animal traditions, including those in great apes, to be fundamentally different from human cultures, largely because of lack of evidence for cumulative processes and normative conformity, but perhaps also because current research on ape culture is usually restricted to behavioral comparisons. Here, we propose to analyze ape culture not only at the surface behavioral level but also at the underlying cognitive level. To this end, we integrate empirical findings in apes with theoretical frameworks developed in developmental psychology regarding the representation of tools and the development of metarepresentational abilities, to characterize the differences between ape and human cultures at the cognitive level. Current data are consistent with the notion of apes possessing mental representations of tools that can be accessed through re-representations: apes may reorganize their knowledge of tools in the form of categories or functional schemes. However, we find no evidence for metarepresentations of cultural knowledge: apes may not understand that they or others hold beliefs about their cultures. The resulting Jourdain Hypothesis, based on Moliere's character, argues that apes express their cultures without knowing that they are cultural beings because of cognitive limitations in their ability to represent knowledge, a determining feature of modern human cultures, allowing representing and modifying the current norms of the group. Differences in metarepresentational processes may thus explain fundamental differences between human and other animals' cultures, notably limitations in cumulative behavior and normative conformity. Future empirical work should focus on how animals mentally represent their cultural knowledge to conclusively determine the ways by which humans are unique in their cultural behavior.Publisher PDFPeer reviewe

Different alleles of the response regulator gene bldM arrest Streptomyces coelicolor development at distinct stages

International audiencewhiK was one of five new whi loci identified in a recent screen of NTG-induced whi mutants and was defined by three mutants, R273, R318 and R655. R273 and R318 produce long, tightly coiled aerial hyphae with frequent septation. In contrast, R655 shows a more severe phenotype; it produces straight, undifferentiated aerial hyphae with very rare short chains of spores. Subcloning and sequencing showed that whiK encodes a member of the FixJ subfamily of response regulators, with a C-terminal helix-turn-helix DNA-binding domain and an apparently typical N-terminal phosphorylation pocket. Unexpectedly, a constructed whiK null mutant failed to form aerial mycelium, showing that different alleles of this locus can arrest Streptomyces coelicolor development at very distinct stages. As a consequence of the null mutant phenotype, whiK was renamed bldM. The bldM null mutant fits into the extracellular signalling cascade proposed for S. coelicolor and is a member of the bldD extracellular complementation group. The three original NTG-induced mutations that defined the whiK/bldM locus each affected the putative phosphorylation pocket. The mutations in R273 and in R318 were the same, replacing a highly conserved glycine (G-62) with aspartate. The more severe mutant, R655, carried a C-7Y substitution adjacent to the highly conserved DD motif at positions 8-9. However, although bldM has all the highly conserved residues associated with the phosphorylation pocket of conventional response regulators, aspartate-54, the putative site of phosphorylation, is not required for bldM function. Constructed mutant alleles carrying either D-54N or D-54A substitutions complemented the bldM null mutant in single copy in trans, and strains carrying the D-54N or the D-54A substitution at the native chromosomal bldM locus sporulated normally. bldM was not phosphorylated in vitro with either of the small-molecule phosphodonors acetyl phosphate or carbamoyl phosphate under conditions in which a control response regulator protein, NtrC, was labelled efficiently