Genome sequencing and analysis of the Tasmanian devil and its transmissible cancer.

The Tasmanian devil (Sarcophilus harrisii), the largest marsupial carnivore, is endangered due to a transmissible facial cancer spread by direct transfer of living cancer cells through biting. Here we describe the sequencing, assembly, and annotation of the Tasmanian devil genome and whole-genome sequences for two geographically distant subclones of the cancer. Genomic analysis suggests that the cancer first arose from a female Tasmanian devil and that the clone has subsequently genetically diverged during its spread across Tasmania. The devil cancer genome contains more than 17,000 somatic base substitution mutations and bears the imprint of a distinct mutational process. Genotyping of somatic mutations in 104 geographically and temporally distributed Tasmanian devil tumors reveals the pattern of evolution and spread of this parasitic clonal lineage, with evidence of a selective sweep in one geographical area and persistence of parallel lineages in other populations.

BAK α6 permits activation by BH3-only proteins and homooligomerization via the canonical hydrophobic groove.

BAK and BAX are the essential effectors of apoptosis because without them a cell is resistant to most apoptotic stimuli. BAK and BAX undergo conformation changes to homooligomerize then permeabilize the mitochondrial outer membrane during apoptosis. How BCL-2 homology 3 (BH3)-only proteins bind to activate BAK and BAX is unclear. We report that BH3-only proteins bind inactive full-length BAK at mitochondria and then dissociate following exposure of the BAK BH3 and BH4 domains before BAK homodimerization. Using a functional obstructive labeling approach, we show that activation of BAK involves important interactions of BH3-only proteins with both the canonical hydrophobic binding groove (α2-5) and α6 at the rear of BAK, with interaction at α6 promoting an open groove to receive a BH3-only protein. Once activated, how BAK homodimers multimerize to form the putative apoptotic pore is unknown. Obstructive labeling of BAK beyond the BH3 domain and hydrophobic groove did not inhibit multimerization and mitochondrial damage, indicating that critical protein-protein interfaces in BAK self-association are limited to the α2-5 homodimerization domain.

Genome analysis of the platypus reveals unique signatures of evolution.

We present a draft genome sequence of the platypus, Ornithorhynchus anatinus. This monotreme exhibits a fascinating combination of reptilian and mammalian characters. For example, platypuses have a coat of fur adapted to an aquatic lifestyle; platypus females lactate, yet lay eggs; and males are equipped with venom similar to that of reptiles. Analysis of the first monotreme genome aligned these features with genetic innovations. We find that reptile and platypus venom proteins have been co-opted independently from the same gene families; milk protein genes are conserved despite platypuses laying eggs; and immune gene family expansions are directly related to platypus biology. Expansions of protein, non-protein-coding RNA and microRNA families, as well as repeat elements, are identified. Sequencing of this genome now provides a valuable resource for deep mammalian comparative analyses, as well as for monotreme biology and conservation.

Reconstruction of the ancestral marsupial karyotype from comparative gene maps.

BACKGROUND: The increasing number of assembled mammalian genomes makes it possible to compare genome organisation across mammalian lineages and reconstruct chromosomes of the ancestral marsupial and therian (marsupial and eutherian) mammals. However, the reconstruction of ancestral genomes requires genome assemblies to be anchored to chromosomes. The recently sequenced tammar wallaby (Macropus eugenii) genome was assembled into over 300,000 contigs. We previously devised an efficient strategy for mapping large evolutionarily conserved blocks in non-model mammals, and applied this to determine the arrangement of conserved blocks on all wallaby chromosomes, thereby permitting comparative maps to be constructed and resolve the long debated issue between a 2n = 14 and 2n = 22 ancestral marsupial karyotype. RESULTS: We identified large blocks of genes conserved between human and opossum, and mapped genes corresponding to the ends of these blocks by fluorescence in situ hybridization (FISH). A total of 242 genes was assigned to wallaby chromosomes in the present study, bringing the total number of genes mapped to 554 and making it the most densely cytogenetically mapped marsupial genome. We used these gene assignments to construct comparative maps between wallaby and opossum, which uncovered many intrachromosomal rearrangements, particularly for genes found on wallaby chromosomes X and 3. Expanding comparisons to include chicken and human permitted the putative ancestral marsupial (2n = 14) and therian mammal (2n = 19) karyotypes to be reconstructed. CONCLUSIONS: Our physical mapping data for the tammar wallaby has uncovered the events shaping marsupial genomes and enabled us to predict the ancestral marsupial karyotype, supporting a 2n = 14 ancestor. Futhermore, our predicted therian ancestral karyotype has helped to understand the evolution of the ancestral eutherian genome.

Cyclin E facilitates dysplastic hepatocytes to bypass G1/S checkpoint in hepatocarcinogenesis.

BACKGROUND AND AIM: By array-comparative genomic hybridization, we demonstrated cyclin E as one of seven genes associated with hepatocellular carcinoma (HCC) development in Ku70 DNA repair-deficient mice. We therefore explored the hypothesis that during hepatocarcinogenesis, cyclin E kinase can overcome the inhibitory effects of p53 and establish whether abnormal miRNA(mi-R)-34, a co-regulator of cyclin E and p53, can account for their interactions as "drivers" of HCC. METHODS: Dysplastic hepatocytes (DNs) and HCCs were generated from diethylnitrosamine (DEN)-injected C57BL/6 male mice at 3-12 months. RESULTS: Cyclin E/cdk2 was barely expressed in normal liver, but was readily detected in dysplastic hepatocytes, localizing to glutathione-S transferase pi-form positive cells dissected by laser-dissection. Cyclin E kinase activity preceded cyclin D1, proliferating cell nuclear antigen expression in DNs and HCCs despite maximal p53 and p21 expression. We confirmed that cyclin E, rather than cyclin D1, is the proliferative driver in hepatocarcinogenesis by immunoprecipitation experiments demonstrating preferential binding of p21 to cyclin D1, allowing cyclin E-mediated "escape" from G1/S checkpoint. We then showed cyclin E was responsible for regulating wild-type p53 by knockdown experiments in primary HCC cells; cyclin E-knockdown increased p53 and p21, diminished anti-apoptotic Bcl-XL and reduced cell viability. Conversely, blocking p53 augmented cyclin E, Bcl-XL expression and increased proliferation. Physiological interactions between cyclin E/p53/p21 were confirmed in primary hepatocytes. miR-34a,c were upregulated in dysplastic murine, human liver and HCCs compared with normal liver, and appeared to be linked to cyclin E/p53. CONCLUSION: Upregulation of functionally active cyclin E via miR34 with loss of p53 function is associated with cell-cycle checkpoint failure increasing proliferative drive that favors hepatocarcinogenesis.

Defining the target specificity of ABT-737 and synergistic antitumor activities in combination with histone deacetylase inhibitors.

The apoptotic and therapeutic activities of the histone deacetylase inhibitor (HDACi) vorinostat are blocked by overexpression of Bcl-2 or Bcl-X(L). Herein, we used the small molecule inhibitor ABT-737 to restore sensitivity of Emu-myc lymphomas overexpressing Bcl-2 or Bcl-X(L) to vorinostat and valproic acid (VPA). Combining low-dose ABT-737 with vorinostat or VPA resulted in synergistic apoptosis of these cells. ABT-737 was ineffective against Emu-myc/Mcl-1 and Emu-myc/A1 cells either as a single agent or in combination with HDACi. However, in contrast to the reported binding specificity data, Emu-myc/Bcl-w lymphomas were insensitive to ABT-737 used alone or in combination with HDACi, indicating that the regulatory activity of ABT-737 is restricted to Bcl-2 and Bcl-X(L). Emu-myc lymphomas that expressed Bcl-2 throughout the tumorigenesis process were especially sensitive to ABT-737, while those forced to overexpress Mcl-1 were not. This supports the notion that tumor cells "addicted" to ABT-737 target proteins (ie, Bcl-2 or Bcl-X(L)) are likely to be the most sensitive target cell population. Our studies provide important preclinical data on the binding specificity of ABT-737 and its usefulness against primary hematologic malignancies when used as a single agent and in combination with HDACi.

A second-generation anchored genetic linkage map of the tammar wallaby (Macropus eugenii).

BACKGROUND: The tammar wallaby, Macropus eugenii, a small kangaroo used for decades for studies of reproduction and metabolism, is the model Australian marsupial for genome sequencing and genetic investigations. The production of a more comprehensive cytogenetically-anchored genetic linkage map will significantly contribute to the deciphering of the tammar wallaby genome. It has great value as a resource to identify novel genes and for comparative studies, and is vital for the ongoing genome sequence assembly and gene ordering in this species. RESULTS: A second-generation anchored tammar wallaby genetic linkage map has been constructed based on a total of 148 loci. The linkage map contains the original 64 loci included in the first-generation map, plus an additional 84 microsatellite loci that were chosen specifically to increase coverage and assist with the anchoring and orientation of linkage groups to chromosomes. These additional loci were derived from (a) sequenced BAC clones that had been previously mapped to tammar wallaby chromosomes by fluorescence in situ hybridization (FISH), (b) End sequence from BACs subsequently FISH-mapped to tammar wallaby chromosomes, and (c) tammar wallaby genes orthologous to opossum genes predicted to fill gaps in the tammar wallaby linkage map as well as three X-linked markers from a published study. Based on these 148 loci, eight linkage groups were formed. These linkage groups were assigned (via FISH-mapped markers) to all seven autosomes and the X chromosome. The sex-pooled map size is 1402.4 cM, which is estimated to provide 82.6% total coverage of the genome, with an average interval distance of 10.9 cM between adjacent markers. The overall ratio of female/male map length is 0.84, which is comparable to the ratio of 0.78 obtained for the first-generation map. CONCLUSIONS: Construction of this second-generation genetic linkage map is a significant step towards complete coverage of the tammar wallaby genome and considerably extends that of the first-generation map. It will be a valuable resource for ongoing tammar wallaby genetic research and assembling the genome sequence. The sex-pooled map is available online at http://compldb.angis.org.au/.

Retrotransposon silencing by DNA methylation can drive mammalian genomic imprinting.

Among mammals, only eutherians and marsupials are viviparous and have genomic imprinting that leads to parent-of-origin-specific differential gene expression. We used comparative analysis to investigate the origin of genomic imprinting in mammals. PEG10 (paternally expressed 10) is a retrotransposon-derived imprinted gene that has an essential role for the formation of the placenta of the mouse. Here, we show that an orthologue of PEG10 exists in another therian mammal, the marsupial tammar wallaby (Macropus eugenii), but not in a prototherian mammal, the egg-laying platypus (Ornithorhynchus anatinus), suggesting its close relationship to the origin of placentation in therian mammals. We have discovered a hitherto missing link of the imprinting mechanism between eutherians and marsupials because tammar PEG10 is the first example of a differentially methylated region (DMR) associated with genomic imprinting in marsupials. Surprisingly, the marsupial DMR was strictly limited to the 5' region of PEG10, unlike the eutherian DMR, which covers the promoter regions of both PEG10 and the adjacent imprinted gene SGCE. These results not only demonstrate a common origin of the DMR-associated imprinting mechanism in therian mammals but provide the first demonstration that DMR-associated genomic imprinting in eutherians can originate from the repression of exogenous DNA sequences and/or retrotransposons by DNA methylation.

Homozygous deletions may be markers of nearby heterozygous mutations: The complex deletion at FRA16D in the HCT116 colon cancer cell line removes exons of WWOX.

Homozygous deletions in cancer cells have been thought to harbor tumor suppressor genes. We show that the 25 and 50 kb homozygous deletions in WWOX in the colon cancer cell line HCT116 result from a complex set of heterozygous deletions, some of which overlap to give homozygous loss. One of the heterozygous deletions has removed exons 6-8 of one allele of WWOX, and there is also a third copy of the distal region of WWOX in an unbalanced translocation. The exon 6-8 deletion results in allele-specific expression of a deleted transcript, which seems likely to be the main biological consequence of the deletions, since similar transcripts are found in other tumors. We show that such a complex set of deletions could form in a single exchange event between two homologous chromosomes, so that the selective advantage of such rearrangements need not be within the homozygous deletion. We conclude that homozygous deletions can be markers of complex rearrangements that have targets outside the homozygous deletion itself and that the target of deletions in the FRA16D region is indeed WWOX, the common outcome being the removal of particular WWOX exons. This article contains supplementary material available at http://www.interscience.wiley.com/jpages/1045-2257/suppmat.

Mcl-1 and Bcl-xL sequestration of Bak confers differential resistance to BH3-only proteins.

The prosurvival Bcl-2 family proteins Mcl-1 and Bcl-xL inhibit apoptosis by sequestering BH3-only proteins such as Bid and Bim (MODE 1) or the effector proteins Bak and Bax (MODE 2). To better understand the contributions of MODE 1 and MODE 2 in blocking cell death, and thus how to bypass resistance to cell death, we examined prescribed mixtures of Bcl-2 family proteins. In a Bim and Bak mixture, Bcl-xL and Mcl-1 each sequestered not only Bim but also Bak as it became activated by Bim. In contrast, in a Bid and Bak mixture, Bcl-xL preferentially sequestered Bid while Mcl-1 preferentially sequestered Bak. Notably, Bcl-xL could sequester Bak in response to the BH3 mimetic ABT-737, despite this molecule targeting Bcl-xL. These findings highlight the importance of Bak sequestration in resistance to anti-cancer treatments, including BH3 mimetics.

Disordered clusters of Bak dimers rupture mitochondria during apoptosis.

During apoptosis, Bak and Bax undergo major conformational change and form symmetric dimers that coalesce to perforate the mitochondrial outer membrane via an unknown mechanism. We have employed cysteine labelling and linkage analysis to the full length of Bak in mitochondria. This comprehensive survey showed that in each Bak dimer the N-termini are fully solvent-exposed and mobile, the core is highly structured, and the C-termini are flexible but restrained by their contact with the membrane. Dimer-dimer interactions were more labile than the BH3:groove interaction within dimers, suggesting there is no extensive protein interface between dimers. In addition, linkage in the mobile Bak N-terminus (V61C) specifically quantified association between dimers, allowing mathematical simulations of dimer arrangement. Together, our data show that Bak dimers form disordered clusters to generate lipidic pores. These findings provide a molecular explanation for the observed structural heterogeneity of the apoptotic pore.

Epigenetic control of mitochondrial cell death through PACS1-mediated regulation of BAX/BAK oligomerization.

PCAF and ADA3 associate within the same macromolecular complexes to control the transcription of many genes, including some that regulate apoptosis. Here we show that PCAF and ADA3 regulate the expression of PACS1, whose protein product is a key component of the machinery that sorts proteins among the trans-Golgi network and the endosomal compartment. We describe a novel role for PACS1 as a regulator of the intrinsic pathway of apoptosis and mitochondrial outer membrane permeabilization. Cells with decreased PACS1 expression were refractory to cell death mediated by a variety of stimuli that operate through the mitochondrial pathway, including human granzyme B, staurosporine, ultraviolet radiation and etoposide, but remained sensitive to TRAIL receptor ligation. The mitochondria of protected cells failed to release cytochrome c as a result of perturbed oligomerization of BAX and BAK. We conclude that PCAF and ADA3 transcriptionally regulate PACS1 and that PACS1 is a key regulator of BAX/BAK oligomerization and the intrinsic (mitochondrial) pathway to apoptosis.

Distribution of breakpoints on chromosome 18 in breast, colorectal, and pancreatic carcinoma cell lines.

Chromosome 18 is frequently rearranged in carcinomas. We explored the distribution of breakpoints affecting chromosome 18 by mapping 56 breakpoints in 26 carcinoma cell lines by fluorescence in situ hybridization (FISH) using bacterial artificial chromosomes (BACs) and band paints. The distribution of breaks among 18 intervals of chromosome 18 was significantly nonrandom. The interval spanning the centromere contained the greatest number of breaks and had the highest average copy number of any interval. There was a high density of breaks close to the centromere as well as actually within the centromere. A cluster of breaks encompassing SMAD4 was associated with the minimum average copy number, consistent with SMAD4 being a tumor suppressor gene. There may be another cluster of breaks around 18q12. We offer two interpretations of the concentration of breaks near the centromere. It may reflect selection for an oncogene near the centromere, or there may be an underlying bias of breakage toward the centromere. We show that the latter is predicted by a simple model that invokes random breakage following anchorage of some random point on the chromosome, or selection of breaks proximal to one of several tumor suppressor genes.

Identification of an activation site in Bak and mitochondrial Bax triggered by antibodies.

During apoptosis, Bak and Bax are activated by BH3-only proteins binding to the α2-α5 hydrophobic groove; Bax is also activated via a rear pocket. Here we report that antibodies can directly activate Bak and mitochondrial Bax by binding to the α1-α2 loop. A monoclonal antibody (clone 7D10) binds close to α1 in non-activated Bak to induce conformational change, oligomerization, and cytochrome c release. Anti-FLAG antibodies also activate Bak containing a FLAG epitope close to α1. An antibody (clone 3C10) to the Bax α1-α2 loop activates mitochondrial Bax, but blocks translocation of cytosolic Bax. Tethers within Bak show that 7D10 binding directly extricates α1; a structural model of the 7D10 Fab bound to Bak reveals the formation of a cavity under α1. Our identification of the α1-α2 loop as an activation site in Bak paves the way to develop intrabodies or small molecules that directly and selectively regulate these proteins.

Dissociation of Bak α1 helix from the core and latch domains is required for apoptosis.

During apoptosis, Bak permeabilizes mitochondria after undergoing major conformational changes, including poorly defined N-terminal changes. Here, we characterize those changes using 11 antibodies that were epitope mapped using peptide arrays and mutagenesis. After Bak activation by Bid, epitopes throughout the α1 helix are exposed indicating complete dissociation of α1 from α2 in the core and from α6-α8 in the latch. Moreover, disulfide tethering of α1 to α2 or α6 blocks cytochrome c release, suggesting that α1 dissociation is required for further conformational changes during apoptosis. Assaying epitope exposure when α1 is tethered shows that Bid triggers α2 movement, followed by α1 dissociation. However, α2 reaches its final position only after α1 dissociates from the latch. Thus, α1 dissociation is a key step in unfolding Bak into three major components, the N terminus, the core (α2-α5) and the latch (α6-α8).

Fas ligand-mediated immune surveillance by T cells is essential for the control of spontaneous B cell lymphomas.

Loss of function of the tumor suppressor gene PRDM1 (also known as BLIMP1) or deregulated expression of the oncogene BCL6 occurs in a large proportion of diffuse large B cell lymphoma (DLBCL) cases. However, targeted mutation of either gene in mice leads to only slow and infrequent development of malignant lymphoma, and despite frequent mutation of BCL6 in activated B cells of healthy individuals, lymphoma development is rare. Here we show that T cells prevent the development of overt lymphoma in mice caused by Blimp1 deficiency or overexpression of Bcl6 in the B cell lineage. Impairment of T cell control results in rapid development of DLBCL-like disease, which can be eradicated by polyclonal CD8(+) T cells in a T cell receptor-, CD28- and Fas ligand-dependent manner. Thus, malignant transformation of mature B cells requires mutations that impair intrinsic differentiation processes and permit escape from T cell-mediated tumor surveillance.

Deciphering the molecular events necessary for synergistic tumor cell apoptosis mediated by the histone deacetylase inhibitor vorinostat and the BH3 mimetic ABT-737.

The concept of personalized anticancer therapy is based on the use of targeted therapeutics through in-depth knowledge of the molecular mechanisms of action of these agents when used alone and in combination. We have identified the apoptotic proteins and pathways necessary for synergistic tumor cell apoptosis and in vivo antitumor responses seen when the HDAC inhibitor vorinostat is combined with the BH3-mimetic ABT-737 in lymphomas overexpressing Bcl-2. Vorinostat "primes" tumors overexpressing Bcl-2 for rapid ABT-737-mediated apoptosis by inducing expression of the BH3-only gene bmf. Moreover, these synergistic effects of vorinostat/ABT-737 were blunted in cells with an inactive p53 pathway or in cells lacking expression of the p53 target gene, noxa. These studies show the important and complex functional interaction between specific proapoptotic BH3-only proteins and the BH3-mimetic compound ABT-737 and provide the most comprehensive functional link between tumor genotype and the apoptotic and therapeutic effects of HDACi combined with ABT-737.

Bird-like sex chromosomes of platypus imply recent origin of mammal sex chromosomes.

In therian mammals (placentals and marsupials), sex is determined by an XX female: XY male system, in which a gene (SRY) on the Y affects male determination. There is no equivalent in other amniotes, although some taxa (notably birds and snakes) have differentiated sex chromosomes. Birds have a ZW female: ZZ male system with no homology with mammal sex chromosomes, in which dosage of a Z-borne gene (possibly DMRT1) affects male determination. As the most basal mammal group, the egg-laying monotremes are ideal for determining how the therian XY system evolved. The platypus has an extraordinary sex chromosome complex, in which five X and five Y chromosomes pair in a translocation chain of alternating X and Y chromosomes. We used physical mapping to identify genes on the pairing regions between adjacent X and Y chromosomes. Most significantly, comparative mapping shows that, contrary to earlier reports, there is no homology between the platypus and therian X chromosomes. Orthologs of genes in the conserved region of the human X (including SOX3, the gene from which SRY evolved) all map to platypus chromosome 6, which therefore represents the ancestral autosome from which the therian X and Y pair derived. Rather, the platypus X chromosomes have substantial homology with the bird Z chromosome (including DMRT1) and to segments syntenic with this region in the human genome. Thus, platypus sex chromosomes have strong homology with bird, but not to therian sex chromosomes, implying that the therian X and Y chromosomes (and the SRY gene) evolved from an autosomal pair after the divergence of monotremes only 166 million years ago. Therefore, the therian X and Y are more than 145 million years younger than previously thought.

Defensins and the convergent evolution of platypus and reptile venom genes.

When the platypus (Ornithorhynchus anatinus) was first discovered, it was thought to be a taxidermist's hoax, as it has a blend of mammalian and reptilian features. It is a most remarkable mammal, not only because it lays eggs but also because it is venomous. Rather than delivering venom through a bite, as do snakes and shrews, male platypuses have venomous spurs on each hind leg. The platypus genome sequence provides a unique opportunity to unravel the evolutionary history of many of these interesting features. While searching the platypus genome for the sequences of antimicrobial defensin genes, we identified three Ornithorhynchus venom defensin-like peptide (OvDLP) genes, which produce the major components of platypus venom. We show that gene duplication and subsequent functional diversification of beta-defensins gave rise to these platypus OvDLPs. The OvDLP genes are located adjacent to the beta-defensins and share similar gene organization and peptide structures. Intriguingly, some species of snakes and lizards also produce venoms containing similar molecules called crotamines and crotamine-like peptides. This led us to trace the evolutionary origins of other components of platypus and reptile venom. Here we show that several venom components have evolved separately in the platypus and reptiles. Convergent evolution has repeatedly selected genes coding for proteins containing specific structural motifs as templates for venom molecules.

DMRT gene cluster analysis in the platypus: new insights into genomic organization and regulatory regions.

We isolated and characterized a cluster of platypus DMRT genes and compared their arrangement, location, and sequence across vertebrates. The DMRT gene cluster on human 9p24.3 harbors, in order, DMRT1, DMRT3, and DMRT2, which share a DM domain. DMRT1 is highly conserved and involved in sexual development in vertebrates, and deletions in this region cause sex reversal in humans. Sequence comparisons of DMRT genes between species have been valuable in identifying exons, control regions, and conserved nongenic regions (CNGs). The addition of platypus sequences is expected to be particularly valuable, since monotremes fill a gap in the vertebrate genome coverage. We therefore isolated and fully sequenced platypus BAC clones containing DMRT3 and DMRT2 as well as DMRT1 and then generated multispecies alignments and ran prediction programs followed by experimental verification to annotate this gene cluster. We found that the three genes have 58-66% identity to their human orthologues, lie in the same order as in other vertebrates, and colocate on 1 of the 10 platypus sex chromosomes, X5. We also predict that optimal annotation of the newly sequenced platypus genome will be challenging. The analysis of platypus sequence revealed differences in structure and sequence of the DMRT gene cluster. Multispecies comparison was particularly effective for detecting CNGs, revealing several novel potential regulatory regions within DMRT3 and DMRT2 as well as DMRT1. RT-PCR indicated that platypus DMRT1 and DMRT3 are expressed specifically in the adult testis (and not ovary), but DMRT2 has a wider expression profile, as it does for other mammals. The platypus DMRT1 expression pattern, and its location on an X chromosome, suggests an involvement in monotreme sexual development.