Search

Melanoma is a deadly form of skin cancer that accounts for a disproportionally large proportion of cancer-related deaths in younger people. Compared with most other skin cancers, a feature of melanoma is its high metastatic capacity, although the mechanisms that confer this are not well understood. The Hippo pathway is a key regulator of organ growth and cell fate that is deregulated in many cancers. To analyse the Hippo pathway in cutaneous melanoma, we generated a transcriptional signature of melanoma cells that overexpressed YAP, the key downstream Hippo pathway oncoprotein. YAP-mediated transcriptional activity varied in melanoma cell lines but did not cluster with known genetic drivers of melanomagenesis such as BRAF and NRAS mutations. Instead, it correlated strongly with published gene expression profiles linked to melanoma cell invasiveness and varied throughout the metastatic cascade in melanoma patient tumours. Consistent with this, YAP was both necessary and sufficient for melanoma cell invasion in vitro. In vivo, YAP promoted spontaneous melanoma metastasis, whilst the growth of YAP-expressing primary tumours was impeded. Finally, we identified the YAP target genes AXL, THBS1 and CYR61 as key mediators of YAP-induced melanoma cell invasion. These data suggest that YAP is a critical regulator of melanoma metastasis.

Background Prostate cancer is caused by genomic aberrations in normal epithelial cells, however clinical translation of findings from analyses of cancer cells alone has been very limited. A deeper understanding of the tumour microenvironment is needed to identify the key drivers of disease progression and reveal novel therapeutic opportunities. Results In this study, the experimental enrichment of selected cell-types, the development of a Bayesian inference model for continuous differential transcript abundance, and multiplex immunohistochemistry permitted us to define the transcriptional landscape of the prostate cancer microenvironment along the disease progression axis. An important role of monocytes and macrophages in prostate cancer progression and disease recurrence was uncovered, supported by both transcriptional landscape findings and by differential tissue composition analyses. These findings were corroborated and validated by spatial analyses at the single-cell level using multiplex immunohistochemistry. Conclusions This study advances our knowledge concerning the role of monocyte-derived recruitment in primary prostate cancer, and supports their key role in disease progression, patient survival and prostate microenvironment immune modulation.

Background The human microbiome plays an important role in cancer. Accumulating evidence indicates that commensal microbiome-derived DNA may be represented in minute quantities in the cell-free DNA of human blood and could possibly be harnessed as a new cancer biomarker. However, there has been limited use of rigorous experimental controls to account for contamination, which invariably affects low-biomass microbiome studies. Results We apply a combination of 16S-rRNA-gene sequencing and droplet digital PCR to determine if the specific detection of cell-free microbial DNA (cfmDNA) is possible in metastatic melanoma patients. Compared to matched stool and saliva samples, the absolute concentration of cfmDNA is low but significantly above the levels detected from negative controls. The microbial community of plasma is strongly influenced by laboratory and reagent contaminants introduced during the DNA extraction and sequencing processes. Through the application of an in silico decontamination strategy including the filtering of amplicon sequence variants (ASVs) with batch dependent abundances and those with a higher prevalence in negative controls, we identify known gut commensal bacteria, such as Faecalibacterium, Bacteroides and Ruminococcus, and also other uncharacterised ASVs. We analyse additional plasma samples, highlighting the potential of this framework to identify differences in cfmDNA between healthy and cancer patients. Conclusions Together, these observations indicate that plasma can harbour a low yet detectable level of cfmDNA. The results highlight the importance of accounting for contamination and provide an analytical decontamination framework to allow the accurate detection of cfmDNA for future biomarker studies in cancer and other diseases.

Accurate and sensitive somatic structural variant calling is essential to the comprehensive characterisation of genomic rearrangements in cancer. Current methods restrict their reporting of structural variants to breakpoint calls only, discarding critical information about genomic rearrangements involving low mappability sequences. Here we present GRIDSS2, the first structural variant caller to explicitly report single breakends - breakpoints in which only one side can be unambiguously determined. We show that treating single breakends as a fundamental genomic rearrangement signal on par with breakpoints enables GRIDSS2 to identify rearrangements in regions currently considered inaccessible to short read sequencing. GRIDSS2 outperforms recent state-of-the-art tools on cell line, validated patient samples, and large cohort data sets. Using a cohort of 3,782 deeply sequenced metastatic cancers, we show that GRIDSS2 achieves a false negative rate 3.1%, with 7.0% of calls being single breakend variants to low-mappability sequence. GRIDSS2 retains high sensitivity and precision even for small events. Our identification of a novel small (32-100bp) duplication overlaps the 50bp threshold typically used for SV calling. Demonstrating the power single breakend calling has in low mappability regions, we find that 47% of somatic centromeric copy number changes can be explained by breaks to non-centromeric sequence, with chromosome 1 exhibiting a unique centromeric rearrangement signature. Finally, we show that somatic structural variants are highly clustered, with GRIDSS2 able to phase 16% of somatic structural variants in our metastasis cohort from short read sequencing alone.

Malaria remains a major public health problem in many countries. Unlike influenza and HIV, where diversity in immunodominant surface antigens is understood geographically to inform disease surveillance, relatively little is known about the global population structure of PfEMP1, the major variant surface antigen of the malaria parasitePlasmodium falciparum. The complexity of thevarmultigene family that encodes PfEMP1 and that diversifies by recombination, has so far precluded its use in malaria surveillance. Recent studies have demonstrated that cost-effective deep sequencing of the region ofvargenes encoding the PfEMP1 DBLα domain and subsequent classification of within host sequences at 96% identity to define unique DBLα types, can reveal structure and strain dynamics within countries. However, to date there has not been a comprehensive comparison of these DBLα types between countries. By leveraging a bioinformatic approach (jumping hidden Markov model) designed specifically for the analysis of recombination withinvargenes and applying it to a dataset of DBLα types from 10 countries, we are able to describe population structure of DBLα types at the global scale. The sensitivity of the approach allows for the comparison of the global dataset to ape samples ofPlasmodium Laveraniaspecies. Our analyses show that the evolution of the parasite population emerging out of Africa underlies current patterns of DBLα type diversity. Most importantly, we can distinguish geographic population structure within Africa between Gabon and Ghana in West Africa and Uganda in East Africa. Our evolutionary findings have translational implications in the context of globalization. Firstly, DBLα type diversity can provide a simple diagnostic framework for geographic surveillance of the rapidly evolving transmission dynamics ofP.falciparum. It can also inform efforts to understand the presence or absence of global, regional and local population immunity to major surface antigen variants. Additionally, we identify a number of highly conserved DBLα types that are present globally that may be of biological significance and warrant further characterization.

Our understanding of genomic heterogeneity in lung cancer is largely based on the analysis of early-stage surgical specimens. Here we used endoscopic sampling of paired primary and intrathoracic metastatic tumors from 11 lung cancer patients to map genomic heterogeneity inoperable lung cancer with deep whole-genome sequencing. Intra-patient heterogeneity in driver or targetable mutations was predominantly in the form of copy number gain. Private mutation signatures, including patterns consistent with defects in homologous recombination, were highly variable both within and between patients. Irrespective of histotype, we observed a smaller than expected number of private mutations, suggesting that ancestral clones accumulated large mutation burdens immediately prior to metastasis. Single-region whole-genome sequencing of from 20 patients showed that tumors in ever-smokers with the strongest tobacco signatures were associated with germline variants in genes implicated in the repair of cigarette-induced DNA damage. Our results suggest that lung cancer precursors in ever-smokers accumulate large numbers of mutations prior to the formation of frank malignancy followed by rapid metastatic spread. In advanced lung cancer, germline variants in DNA repair genes may interact with the airway environment to influence the pattern of founder mutations, whereas similar interactions with the tumor microenvironment may play a role in the acquisition of mutations following metastasis.

Tumors of distinct tissues of origin and genetic makeup display common hallmark cellular phenotypes, including sustained proliferation, suppression of cell death, and altered metabolism. These phenotypic commonalities have been proposed to stem from disruption of conserved regulatory mechanisms evolved during the transition to multicellularity to control fundamental cellular processes such as growth and replication. Dating the evolutionary emergence of human genes through phylostratigraphy uncovered close association between gene age and expression level in RNA sequencing data from The Cancer Genome Atlas for seven solid cancers. Genes conserved with unicellular organisms were strongly up-regulated, whereas genes of metazoan origin were primarily inactivated. These patterns were most consistent for processes known to be important in cancer, implicating both selection and active regulation during malignant transformation. The coordinated expression of strongly interacting multicellularity and unicellularity processes was lost in tumors. This separation of unicellular and multicellular functions appeared to be mediated by 12 highly connected genes, marking them as important general drivers of tumorigenesis. Our findings suggest common principles closely tied to the evolutionary history of genes underlie convergent changes at the cellular process level across a range of solid cancers. We propose altered activity of genes at the interfaces between multicellular and unicellular regions of human gene regulatory networks activate primitive transcriptional programs, driving common hallmark features of cancer. Manipulation of cross-talk between biological processes of different evolutionary origins may thus present powerful and broadly applicable treatment strategies for cancer.

The parasitic mite Sarcoptes scabiei is an economically highly significant parasite of the skin of humans and animals worldwide. In humans, this mite causes a neglected tropical disease (NTD), called scabies. This disease results in major morbidity, disability, stigma and poverty globally and is often associated with secondary bacterial infections. Currently, anti-scabies treatments are not sufficiently effective, resistance to them is emerging and no vaccine is available. Here, we report the first high-quality genome and transcriptomic data for S. scabiei. The genome is 56.6 Mb in size, has a a repeat content of 10.6% and codes for 9,174 proteins. We explored key molecules involved in development, reproduction, host-parasite interactions, immunity and disease. The enhanced ‘omic data sets for S. scabiei represent comprehensive and critical resources for genetic, functional genomic, metabolomic, phylogenetic, ecological and/or epidemiological investigations, and will underpin the design and development of new treatments, vaccines and/or diagnostic tests., Author summary Scabies is a highly significant parasitic disease caused by the mite S. scabiei. This NTD has a major adverse impact in disadvantaged communities around the world, particularly when associated with secondary bacterial infections and clinical complications. Here we report the first high-quality genome and transcriptomic data for S. scabiei and explore molecular aspects of S. scabiei/scabies. This genome (56.6 Mb, encoding ~ 9,200 proteins) provides a solid foundation for fundamental investigations of the molecular biology of the mite, host-parasite interactions and disease processes as well as for translational research to develop new treatments, vaccines and diagnostic tests.

Background The malaria parasite Plasmodium falciparum has an unusually euchromatic genome with poorly conserved positioning of nucleosomes in intergenic sequences and poorly understood mechanisms of gene regulation. Variant histones and histone modifications determine nucleosome stability and recruit trans factors, but their combinatorial contribution to gene regulation is unclear. Results Here, we show that the histone H3 acetylations H3K18ac and H3K27ac and the variant histone Pf H2A.Z are enriched together at regulatory sites upstream of genes. H3K18ac and H3K27ac together dynamically mark regulatory regions of genes expressed during the asexual life cycle. In contrast, H3K4me1 is depleted in intergenic sequence and dynamically depleted upstream of expressed genes. The temporal pattern of H3K27ac and H3K18ac enrichment indicates that they accumulate during S phase and mitosis and are retained at regulatory sequences until at least G1 phase and after cessation of expression of the cognate genes. We integrated our ChIPseq data with existing datasets to show that in schizont stages H3K18ac, H3K27ac and Pf H2A.Z colocalise with the transcription factor PfAP2-I and the bromodomain protein PfBDP1 and are enriched at stably positioned nucleosomes within regions of exposed DNA at active transcriptional start sites. Using transient transfections we showed that sequences enriched with colocalised H3K18ac, H3K27ac and Pf H2A.Z possess promoter activity in schizont stages, but no enhancer-like activity. Conclusions The dynamic H3 acetylations define P. falciparum regulatory sequences and contribute to gene activation. These findings expand the knowledge of the chromatin landscape that regulates gene expression in P. falciparum.

Cancer is an evolutionary and ecological process in which complex interactions between tumour cells and their environment share many similarities with organismal evolution. Tumour cells with highest adaptive potential have a selective advantage over less fit cells. Naturally occurring transmissible cancers provide an ideal model system for investigating the evolutionary arms race between cancer cells and their surrounding micro-environment and macro-environment. However, the evolutionary landscapes in which contagious cancers reside have not been subjected to comprehensive investigation. Here, we provide a multifocal analysis of transmissible tumour progression and discuss the selection forces that shape it. We demonstrate that transmissible cancers adapt to both their micro-environment and macro-environment, and evolutionary theories applied to organisms are also relevant to these unique diseases. The three naturally occurring transmissible cancers, canine transmissible venereal tumour (CTVT) and Tasmanian devil facial tumour disease (DFTD) and the recently discovered clam leukaemia, exhibit different evolutionary phases: (i) CTVT, the oldest naturally occurring cell line is remarkably stable; (ii) DFTD exhibits the signs of stepwise cancer evolution; and (iii) clam leukaemia shows genetic instability. While all three contagious cancers carry the signature of ongoing and fairly recent adaptations to selective forces, CTVT appears to have reached an evolutionary stalemate with its host, while DFTD and the clam leukaemia appear to be still at a more dynamic phase of their evolution. Parallel investigation of contagious cancer genomes and transcriptomes and of their micro-environment and macro-environment could shed light on the selective forces shaping tumour development at different time points: during the progressive phase and at the endpoint. A greater understanding of transmissible cancers from an evolutionary ecology perspective will provide novel avenues for the prevention and treatment of both contagious and non-communicable cancers.

Background The suboptimal sensitivity and specificity of available diagnostic methods for scabies hampers clinical management, trials of new therapies and epidemiologic studies. Additionally, parasitologic diagnosis by microscopic examination of skin scrapings requires sample collection with a sharp scalpel blade, causing discomfort to patients and difficulty in children. Polymerase chain reaction (PCR)-based diagnostic assays, combined with non-invasive sampling methods, represent an attractive approach. In this study, we aimed to develop a real-time probe-based PCR test for scabies, test a non-invasive sampling method and evaluate its diagnostic performance in two clinical settings. Methodology/Principal findings High copy-number repetitive DNA elements were identified in draft Sarcoptes scabiei genome sequences and used as assay targets for diagnostic PCR. Two suitable repetitive DNA sequences, a 375 base pair microsatellite (SSR5) and a 606 base pair long tandem repeat (SSR6), were identified. Diagnostic sensitivity and specificity were tested using relevant positive and negative control materials and compared to a published assay targeting the mitochondrial cox1 gene. Both assays were positive at a 1:100 dilution of DNA from a single mite; no amplification was observed in DNA from samples from 19 patients with other skin conditions nor from house dust, sheep or dog mites, head and body lice or from six common skin bacterial and fungal species. Moderate sensitivity of the assays was achieved in a pilot study, detecting 5/7 (71.4% [95% CI: 29.0% - 96.3%]) of clinically diagnosed untreated scabies patients). Greater sensitivity was observed in samples collected by FLOQ swabs compared to skin scrapings. Conclusions/Significance This newly developed qPCR assay, combined with the use of an alternative non-invasive swab sampling technique offers the possibility of enhanced diagnosis of scabies. Further studies will be required to better define the diagnostic performance of these tests., Author summary As scabies control efforts continue to grow, scarcity of diagnostic options hinders success of elimination efforts in endemic areas. Efficiency in large-scale monitoring is further obstructed by invasive sample collection techniques, which are often uncomfortable for patients, and lack sensitivity. We have developed two PCR-based diagnostic assays targeting repetitive DNA elements. These were identified using new data on the S. scabiei genome. Targeting these elements by PCR improved the detection of scabies DNA. Enhanced sensitivity was demonstrated when tested against routine microscopy and a published PCR-based diagnostic assay. When combined with a non-invasive, effective FLOQ swab sampling method, the developed qPCR-based assays may provide a useful complementary tool for diagnosis of scabies, and its application will likely improve scabies control in target populations.

In recent years, many software packages for identifying structural variants (SVs) using whole-genome sequencing data have been released. When published, a new method is commonly compared with those already available, but this tends to be selective and incomplete. The lack of comprehensive benchmarking of methods presents challenges for users in selecting methods and for developers in understanding algorithm behaviours and limitations. Here we report the comprehensive evaluation of 10 SV callers, selected following a rigorous process and spanning the breadth of detection approaches, using high-quality reference cell lines, as well as simulations. Due to the nature of available truth sets, our focus is on general-purpose rather than somatic callers. We characterise the impact on performance of event size and type, sequencing characteristics, and genomic context, and analyse the efficacy of ensemble calling and calibration of variant quality scores. Finally, we provide recommendations for both users and methods developers., A number of computational methods have been developed for calling structural variants (SVs) using short read sequencing data. Here, the authors perform a comprehensive benchmarking analysis comparing 10 general-purpose callers and provide recommendations for both users and methods developers.

Background Scabies is worldwide one of the most common, yet neglected, parasitic skin infections, affecting a wide range of mammals including humans. Limited treatment options and evidence of emerging mite resistance against the currently used drugs drive our research to explore new therapeutic candidates. Previously, we discovered a multicopy family of genes encoding cysteine proteases with their catalytic sites inactivated by mutation (SMIPP-Cs). This protein family is unique in parasitic scabies mites and is absent in related non-burrowing mites. We postulated that the SMIPP-Cs have evolved as an adaptation to the parasitic lifestyle of the scabies mite. To formulate testable hypotheses for their functions and to propose possible strategies for translational research we investigated whether the SMIPP-Cs are common to all scabies mite varieties and where within the mite body as well as when throughout the parasitic life-cycle they are expressed. Results SMIPP-C sequences from human, pig and dog mites were analysed bioinformatically and the phylogenetic relationships between the SMIPP-C multi-copy gene families of human, pig and dog mites were established. Results suggest that amplification of the SMIPP-C genes occurred in a common ancestor and individual genes evolved independently in the different mite varieties. Recombinant human mite SMIPP-C proteins were produced and used for murine polyclonal antibody production. Immunohistology on skin sections from human patients localised the SMIPP-Cs in the mite gut and in mite faeces within in the epidermal skin burrows. SMIPP-C transcription into mRNA in different life stages was assessed in human and pig mites by reverse transcription followed by droplet digital PCR (ddPCR). High transcription levels of SMIPP-C genes were detected in the adult female life stage in comparison to all other life stages. Conclusions The fact that the SMIPP-Cs are unique to three Sarcoptes varieties, present in all burrowing life stages and highly expressed in the digestive system of the infective adult female life stage may highlight an essential role in parasitism. As they are excreted from the gut in scybala they presumably are able to interact or interfere with host proteins present in the epidermis. Electronic supplementary material The online version of this article (10.1186/s13071-018-2862-0) contains supplementary material, which is available to authorized users.

Genomic rearrangements are common in cancer, with demonstrated links to disease progression and treatment response. These rearrangements can be complex, resulting in fusions of multiple chromosomal fragments and generation of derivative chromosomes. Although methods exist for detecting individual fusions, they are generally unable to reconstruct complex chained events. To overcome these limitations, we adopted a new optical mapping approach, allowing megabase-length genome maps to be reconstructed and rearranged genomes to be visualized without loss of integrity. Whole-genome mapping (Bionano Genomics) of a well-studied highly rearranged liposarcoma cell line resulted in 3338 assembled consensus genome maps, including 72 fusion maps. These fusion maps represent 112.3 Mb of highly rearranged genomic regions, illuminating the complex architecture of chained fusions, including content, order, orientation, and size. Spanning the junction of 147 chromosomal translocations, we found a total of 28 Mb of interspersed sequences that could not be aligned to the reference genome. Traversing these interspersed sequences using short-read sequencing breakpoint calls, we were able to identify and place 399 sequencing fragments within the optical mapping gaps, thus illustrating the complementary nature of optical mapping and short-read sequencing. We demonstrate that optical mapping provides a powerful new approach for capturing a higher level of complex genomic architecture, creating a scaffold for renewed interpretation of sequencing data of particular relevance to human cancer.

Within the human host, the malaria parasite Plasmodium falciparum is exposed to multiple selection pressures. The host environment changes dramatically in severe malaria, but the extent to which the parasite responds to—or is selected by—this environment remains unclear. From previous studies, the parasites that cause severe malaria appear to increase expression of a restricted but poorly defined subset of the PfEMP1 variant, surface antigens. PfEMP1s are major targets of protective immunity. Here, we used RNA sequencing (RNAseq) to analyse gene expression in 44 parasite isolates that caused severe and uncomplicated malaria in Papuan patients. The transcriptomes of 19 parasite isolates associated with severe malaria indicated that these parasites had decreased glycolysis without activation of compensatory pathways; altered chromatin structure and probably transcriptional regulation through decreased histone methylation; reduced surface expression of PfEMP1; and down-regulated expression of multiple chaperone proteins. Our RNAseq also identified novel associations between disease severity and PfEMP1 transcripts, domains, and smaller sequence segments and also confirmed all previously reported associations between expressed PfEMP1 sequences and severe disease. These findings will inform efforts to identify vaccine targets for severe malaria and also indicate how parasites adapt to—or are selected by—the host environment in severe malaria., Author summary Infection by Plasmodium falciparum—the parasite responsible for malaria in humans—can result in a severe disease that can be fatal or in an uncomplicated disease that can be resolved by the host immune system. However, whether the parasites causing severe disease differ from those causing uncomplicated disease is unknown. Several strands of evidence have suggested that parasites causing severe disease may express a restricted set of the Plasmodium falciparum Erythrocyte Membrane Protein 1 (PfEMP1) proteins. PfEMP1 proteins are expressed on the surface of the infected red blood cells and elicit protective immunity. We compared the transcriptomes of parasites causing severe and uncomplicated malaria to determine whether these parasites differed in the genes they expressed. We found that the parasites causing severe malaria had altered expression of genes involved in basic metabolism, nuclear processes, and surface expression of PfEMP1. The parasites causing severe malaria had up-regulated expression of a set of PfEMP1 proteins. Some of these PfEMP1s had been previously implicated in severe malaria, lending support to our data. Multiple associations identified between severe malaria and expressed PfEMP1 sequences were novel. These novel, severe disease–associated PfEMP1 sequences could be useful for informing design of vaccines targeting severe malaria disease.

Cellular dedifferentiation is the regression of a cell from a specialized state to a more multipotent state and is implicated in cancer. However, the transcriptional network that prevents differentiated cells from reacquiring stem cell fate is so far unclear. Neuroblasts (NBs), the Drosophila neural stem cells, are a model for the regulation of stem cell self-renewal and differentiation. Here we show that the Drosophila zinc finger transcription factor Nervous fingers 1 (Nerfin-1) locks neurons into differentiation, preventing their reversion into NBs. Following Prospero-dependent neuronal specification in the ganglion mother cell (GMC), a Nerfin-1-specific transcriptional program maintains differentiation in the post-mitotic neurons. The loss of Nerfin-1 causes reversion to multipotency and results in tumors in several neural lineages. Both the onset and rate of neuronal dedifferentiation in nerfin-1 mutant lineages are dependent on Myc- and target of rapamycin (Tor)-mediated cellular growth. In addition, Nerfin-1 is required for NB differentiation at the end of neurogenesis. RNA sequencing (RNA-seq) and chromatin immunoprecipitation (ChIP) analysis show that Nerfin-1 administers its function by repression of self-renewing-specific and activation of differentiation-specific genes. Our findings support the model of bidirectional interconvertibility between neural stem cells and their post-mitotic progeny and highlight the importance of the Nerfin-1-regulated transcriptional program in neuronal maintenance.

The identification of genomic rearrangements, particularly in cancers, with high sensitivity and specificity using massively parallel sequencing remains a major challenge. Here, we describe the Genome Rearrangement IDentification Software Suite (GRIDSS), a high-speed structural variant (SV) caller that performs efficient genome-wide break-end assembly prior to variant calling using a novel positional de Bruijn graph assembler. By combining assembly, split read and read pair evidence using a probabilistic scoring, GRIDSS achieves high sensitivity and specificity on simulated, cell line and patient tumour data, recently winning SV sub-challenge #5 of the ICGC-TCGA DREAM Somatic Mutation Calling Challenge. On human cell line data, GRIDSS halves the false discovery rate compared to other recent methods. GRIDSS identifies non-template sequence insertions, micro-homologies and large imperfect homologies, and supports multi-sample analysis. GRIDSS is freely available athttps://github.com/PapenfussLab/gridss.

We present SVclone, a computational method for inferring the cancer cell fraction of structural variant breakpoints from whole-genome sequencing data. We validate our approach using simulated and real tumour samples, and demonstrate its utility on 2,778 whole-genome sequenced tumours. We find a subset of liver, breast and ovarian cancer cases with decreased overall survival that have subclonally enriched copy-number neutral rearrangements, an observation that could not be discovered with currently available methods.

Genomic rearrangements are common in cancer, with demonstrated links to disease progression and treatment response. These rearrangements can be complex, resulting in fusions of multiple chromosomal fragments and generation of derivative chromosomes. While methods exist for detecting individual fusions, they are generally unable to reconstruct complex chained events. To overcome these limitations, we adopted a new optical mapping approach, allowing for megabase length DNA to be captured, and in turn rearranged genomes to be visualized without loss of integrity. Whole genome mapping (Bionano Genomics) of a well-studied highly rearranged liposarcoma cell line, resulted in 3,338 assembled haploid genome maps, including 101 fusion maps. These fusion maps represent 175 Mb of highly rearranged genomic regions, illuminating the complex architecture of chained fusions, including content, order, orientation, and size. Spanning the junction of 151 chromosomal translocations, we found a total of 32 Mb of novel interspersed sequences that were not detected from short-read sequencing. We demonstrate that optical mapping provides a powerful new approach for capturing a higher level of complex genomic architecture, creating a scaffold for renewed interpretation of sequencing data of particular relevance to human cancer.

SummaryWe isolated and analyzed, at single-nucleotide resolution, cancer-associated neochromosomes from well- and/or dedifferentiated liposarcomas. Neochromosomes, which can exceed 600 Mb in size, initially arise as circular structures following chromothripsis involving chromosome 12. The core of the neochromosome is amplified, rearranged, and corroded through hundreds of breakage-fusion-bridge cycles. Under selective pressure, amplified oncogenes are overexpressed, while coamplified passenger genes may be silenced epigenetically. New material may be captured during punctuated chromothriptic events. Centromeric corrosion leads to crisis, which is resolved through neocentromere formation or native centromere capture. Finally, amplification terminates, and the neochromosome core is stabilized in linear form by telomere capture. This study investigates the dynamic mutational processes underlying the life history of a special form of cancer mutation.

Background The Major Histocompatibility Complex (MHC) class I family of genes encode for molecules that have well-conserved structures, but have evolved to perform diverse functions. The availability of the gray, short-tailed opossum, Monodelphis domestica whole genome sequence has allowed for analysis of MHC class I gene content in this marsupial. Utilization of a novel method to search for MHC related domain structures revealed a previously unknown family of MHC class I-related genes. These genes, named UT1-17, are clustered on chromosome 1 in the opossum, unlinked to the MHC region. UT genes are only found in marsupial and monotreme genomes, consistent with being ancient in mammals yet lost in eutherian mammals. This study investigates the expression and polymorphism of the UT loci in the opossum to gain insight into their possible function. Results Of the 17 opossum UT genes, most have restricted tissue transcription patterns, with the thymus and skin being the most common sites. Full-length structure of 11 UT transcripts revealed genes varying between five and eight exons, typical for class I family members. There were only two alternative splice variants found. The UT genes also have limited polymorphism and little evidence of positive selection. One locus, UT8, was chosen for further analysis due to its conservation amongst marsupials and generic characteristics. UT8 transcription is limited to developing αβ thymocytes, and is absent from mature αβ T cells in peripheral lymphoid tissues. Conclusion The overall characteristics and features of UT genes including low polymorphism and restricted tissue expression make it likely that the molecules encoded by UT genes perform roles other than antigenic peptide presentation. Electronic supplementary material The online version of this article (doi:10.1186/s12865-016-0181-9) contains supplementary material, which is available to authorized users.

Sherlock Holmes was renowned for his ability to discern subtle clues imprinted on people’s history to solve a case. A study recently published in Nature identified breast cancer culprits, scrutinising the genome of a large number of patients.1 The authors found specific genomic ‘prints’ or ‘signatures’ left by mutational processes in the DNA, providing valuable insights into the enigmatic genesis of breast cancer. The mutational landscape of 560 breast cancer patients was determined using whole-genome sequencing, RNAseq, miRNAseq and DNA methylation arrays.1 Malignant tissues, mostly obtained from drug-naive patients, were compared with non-malignant tissues (lymphocytes, adjacent normal tissue or skin). The power of this study is that the resulting catalogue of somatic mutations was used not only to reveal driver mutations that confer clonal advantage, but also to empirically identify the mutational processes that generate these somatic mutations.

Summary 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. PaperClip, Graphical Abstract Highlights ► Whole-genome sequences of the Tasmanian devil and two distant cancer subclones ► The Tasmanian devil cancer lineage originated recently in a female devil ► The devil cancer genome is relatively stable despite ongoing evolution ► Clonal divergence and geographic spread elucidated through patterns of mutation, Whole-genome sequences of the Tasmanian devil and two devil cancer subclones suggest that the cancer first arose from a female devil and that the clone has subsequently genetically diverged during its spread across Tasmania.

Gene duplication followed by adaptive selection is believed to be the primary driver of venom evolution. However, to date, no studies have evaluated the importance of gene duplications for venom evolution using a genomic approach. The availability of a sequenced genome and a venom gland transcriptome for the enigmatic platypus provides a unique opportunity to explore the role that gene duplication plays in venom evolution. Here, we identify gene duplication events and correlate them with expressed transcripts in an in-season venom gland. Gene duplicates (1,508) were identified. These duplicated pairs (421), including genes that have undergone multiple rounds of gene duplications, were expressed in the venom gland. The majority of these genes are involved in metabolism and protein synthesis not toxin functions. Twelve secretory genes including serine proteases, metalloproteinases, and protease inhibitors likely to produce symptoms of envenomation such as vasodilation and pain were detected. Only 16 of 107 platypus genes with high similarity to known toxins evolved through gene duplication. Platypus venom C-type natriuretic peptides and nerve growth factor do not possess lineage-specific gene duplicates. Extensive duplications, believed to increase the potency of toxic content and promote toxin diversification, were not found. This is the first study to take a genome-wide approach in order to examine the impact of gene duplication on venom evolution. Our findings support the idea that adaptive selection acts on gene duplicates to drive the independent evolution and functional diversification of similar venom genes in venomous species. However, gene duplications alone do not explain the “venome” of the platypus. Other mechanisms, such as alternative splicing and mutation, may be important in venom innovation.

More than 25 loci have been linked to type 1 diabetes (T1D) in the nonobese diabetic (NOD) mouse, but identification of the underlying genes remains challenging. We describe here the positional cloning of a T1D susceptibility locus, Idd11, located on mouse chromosome 4. Sequence analysis of a series of congenic NOD mouse strains over a critical 6.9-kb interval in these mice and in 25 inbred strains identified several haplotypes, including a unique NOD haplotype, associated with varying levels of T1D susceptibility. Haplotype diversity within this interval between congenic NOD mouse strains was due to a recombination hotspot that generated four crossover breakpoints, including one with a complex conversion tract. The Idd11 haplotype and recombination hotspot are located within a predicted gene of unknown function, which exhibits decreased expression in relevant tissues of NOD mice. Notably, it was the recombination hotspot that aided our mapping of Idd11 and confirms that recombination hotspots can create genetic variation affecting a common polygenic disease. This finding has implications for human genetic association studies, which may be affected by the approximately 33,000 estimated hotspots in the genome.

Techniques for targeted genetic disruption in Plasmodium, the causative agent of malaria, are currently intractable for those genes that are essential for blood stage development. The ability to use RNA interference (RNAi) to silence gene expression would provide a powerful means to gain valuable insight into the pathogenic blood stages but its functionality in Plasmodium remains controversial. Here we have used various RNA-based gene silencing approaches to test the utility of RNAi in malaria parasites and have undertaken an extensive comparative genomics search using profile hidden Markov models to clarify whether RNAi machinery exists in malaria. These investigative approaches revealed that Plasmodium lacks the enzymology required for RNAi-based ablation of gene expression and indeed no experimental evidence for RNAi was observed. In its absence, the most likely explanations for previously reported RNAi-mediated knockdown are either the general toxicity of introduced RNA (with global down-regulation of gene expression) or a specific antisense effect mechanistically distinct from RNAi, which will need systematic analysis if it is to be of use as a molecular genetic tool for malaria parasites.

The platypus, as an egg-laying mammal, displays an unusual mixture of reptilian and mammalian characteristics. It is also venomous, and further investigations into its little-studied venom may lead to the development of novel pharmaceuticals and drug targets and provide insights into the origins of mammalian venom. Here we investigate the expression patterns of antimicrobial genes called defensins, and also the venom peptides called defensin-like peptides (OvDLPs). We show, in the first expression study on any platypus venom gene, that the OvDLPs are expressed in a greater range of tissues than would be expected for genes with specific venom function, and thus that they may have a wider role than previously suspected.

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.

The scabies mite, Sarcoptes scabiei, is an obligate parasite of the skin that infects humans and other animal species, causing scabies, a contagious disease characterized by extreme itching. Scabies infections are a major health problem, particularly in remote Indigenous communities in Australia, where co-infection of epidermal scabies lesions by Group A Streptococci or Staphylococcus aureus is thought to be responsible for the high rate of rheumatic heart disease and chronic kidney disease. We collected and separately sequenced mite DNA from several pools of thousands of whole mites from a porcine model of scabies (S. scabiei var. suis) and two human patients (S. scabiei var. hominis) living in different regions of northern Australia. Our sequencing samples the mite and its metagenome, including the mite gut flora and the wound micro-environment. Here, we describe the mitochondrial genome of the scabies mite. We developed a new de novo assembly pipeline based on a bait-and-reassemble strategy, which produced a 14 kilobase mitochondrial genome sequence assembly. We also annotated 35 genes and have compared these to other Acari mites. We identified single nucleotide polymorphisms (SNPs) and used these to infer the presence of six haplogroups in our samples, Remarkably, these fall into two closely-related clades with one clade including both human and pig varieties. This supports earlier findings that only limited genetic differences may separate some human and animal varieties, and raises the possibility of cross-host infections. Finally, we used these mitochondrial haplotypes to show that the genetic diversity of individual infections is typically small with 1–3 distinct haplotypes per infestation., Author Summary The scabies mite is a skin parasite that infects humans and other animal species, causing scabies, a contagious disease characterized by extreme itching. Scabies infections are a major health problem in developing countries and in indigenous Australian populations, where scabies is associated with pyoderma (skin sores) and linked to rheumatic fever and rheumatic heart disease. Little is known about the genetics of the scabies mite. We have assembled the mitochondrial genome of scabies mites obtained from human patients in Australia and from a pig model. While investigating the genetic diversity of each infestation, we found that mitochrondial genomes clustered into two broad clades and showed limited genetic diversity within each infestation. Remarkably, one closely related clade included both human and pig mites, suggesting that mite transmission from pig to human may be possible. This could have major implications in the management of porcine mange and human scabies.

Tasmanian devil (Sarcophilus harrisii) pouch young, like other marsupials, are born underdeveloped and immunologically naïve, and are unable to mount an adaptive immune response. The mother’s milk provides nutrients for growth and development as well as providing passive immunity. To better understand immune response in this endangered species, we set out to characterise the genes involved in passive immunity by sequencing and annotating the transcriptome of a devil milk sample collected during mid-lactation. At mid-lactation we expect the young to have heightened immune responses, as they have emerged from the pouch, encountering new pathogens. A total of 233,660 transcripts were identified, including approximately 17,827 unique protein-coding genes and 846 immune genes. The most highly expressed transcripts were dominated by milk protein genes such as those encoding early lactation protein, late lactation proteins,α-lactalbumin,α-casein andβ-casein. There were numerous highly expressed immune genes including lysozyme, whey acidic protein, ferritin and major histocompatibility complex I and II. Genes encoding immunoglobulins, antimicrobial peptides, chemokines and immune cell receptors were also identified. The array of immune genes identified in this study reflects the importance of the milk in providing immune protection to Tasmanian devil young and provides the first insight into Tasmanian devil milk.

We report a high-quality draft of the genome sequence of the grey, short-tailed opossum (Monodelphis domestica). As the first metatherian ('marsupial') species to be sequenced, the opossum provides a unique perspective on the organization and evolution of mammalian genomes. Distinctive features of the opossum chromosomes provide support for recent theories about genome evolution and function, including a strong influence of biased gene conversion on nucleotide sequence composition, and a relationship between chromosomal characteristics and X chromosome inactivation. Comparison of opossum and eutherian genomes also reveals a sharp difference in evolutionary innovation between protein-coding and non-coding functional elements. True innovation in protein-coding genes seems to be relatively rare, with lineage-specific differences being largely due to diversification and rapid turnover in gene families involved in environmental interactions. In contrast, about 20% of eutherian conserved non-coding elements (CNEs) are recent inventions that postdate the divergence of Eutheria and Metatheria. A substantial proportion of these eutherian-specific CNEs arose from sequence inserted by transposable elements, pointing to transposons as a major creative force in the evolution of mammalian gene regulation.

Neochromosomes are a little-studied class of chromosome-scale mutations that drive some cancers. By sequencing isolated neochromosomes from liposarcomas, we recently defined their structure at single-nucleotide resolution and proposed a model for their life history. Here, we summarize that work, highlighting significant aspects and providing historical context and insight into the discovery process.

To follow the fate of CD8+ T cells responsive to Plasmodium berghei ANKA (PbA) infection, we generated an MHC I-restricted TCR transgenic mouse line against this pathogen. T cells from this line, termed PbT-I T cells, were able to respond to blood-stage infection by PbA and two other rodent malaria species, P. yoelii XNL and P. chabaudi AS. These PbT-I T cells were also able to respond to sporozoites and to protect mice from liver-stage infection. Examination of the requirements for priming after intravenous administration of irradiated sporozoites, an effective vaccination approach, showed that the spleen rather than the liver was the main site of priming and that responses depended on CD8α+ dendritic cells. Importantly, sequential exposure to irradiated sporozoites followed two days later by blood-stage infection led to augmented PbT-I T cell expansion. These findings indicate that PbT-I T cells are a highly versatile tool for studying multiple stages and species of rodent malaria and suggest that cross-stage reactive CD8+ T cells may be utilized in liver-stage vaccine design to enable boosting by blood-stage infections., Author Summary Malaria is a disease caused by Plasmodium species, which have a highly complex life cycle involving both liver and blood stages of mammalian infection. To prevent disease, one strategy has been to induce CD8+ T cells against liver-stage parasites, usually by immunization with stage-specific antigens. Here we describe a T cell receptor specificity that recognizes an antigen expressed in both the liver and blood stages of several rodent Plasmodium species. We generated a T cell receptor transgenic mouse with this specificity and showed that T cells from this line could protect against liver-stage infection. We used this novel tool to identify the site and cell-type involved in priming to a recently developed intravenous attenuated sporozoite vaccine shown to have efficacy in humans. We showed that CD8+ T cells with this specificity could protect against liver-stage infection while causing pathology to the blood stage. Finally, we provided evidence that T cells with cross-stage specificity can be primed and boosted on alternative stages, raising the possibility that antigens expressed in multiple stages might be ideal vaccine candidates for generating strong immunity to liver-stage parasites.

T cell repertoire diversity and clonotype follow-up in vaccination, cancer, infectious and immune diseases represent a major challenge owing to the enormous complexity of the data generated. Here we describe a next generation methodology, which combines 5′RACE PCR, 454 sequencing and, for analysis, IMGT, the international ImMunoGeneTics information system (IMGT), IMGT/HighV-QUEST web portal and IMGT-ONTOLOGY concepts. The approach is validated in a human case study of T cell receptor beta (TRB) repertoire, by chronologically tracking the effects of influenza vaccination on conventional and regulatory T cell subpopulations. The IMGT/HighV-QUEST paradigm defines standards for genotype/haplotype analysis and characterization of IMGT clonotypes for clonal diversity and expression and achieves a degree of resolution for next generation sequencing verifiable by the user at the sequence level, while providing a normalized reference immunoprofile for human TRB., Dynamic changes in T cell repertoire underlie immune responses during infection, allergy, autoimmunity and cancer. Here, Li et al. present a workflow for high throughput sequencing and analysis of T cell receptor sequences, and use it to monitor the T cell response to influenza vaccination in a human patient.

We sought to explore the etiology, frequency and clinicopathologic associations of clustered noncoding hotspot mutations in cutaneous malignancies (melanoma, basal cell carcinoma (BCC) and Merkel cell carcinoma (MCC)). We did so by applying a sliding window approach which identified clusters in melanoma whole exome and genome sequencing data. These regions were then evaluated for recurrence in other skin malignancies by using whole exome data from BCC and MCC samples; motifs were extracted using a maximum expectation algorithm. We then sequenced 170 annotated melanoma clinical samples using a targeted amplicon panel which permitted analysis of clinicopathological correlations. Numerous recurrent clustered hotspot mutations were identified in melanoma whole exome and genome data sets in proximal promoters, with some gene-associated clusters present in up to 25% of cases. Mutations were clustered around canonical ETS and SP1 transcription factor binding sites and had a clear UV signature. These promoter mutations correlated with overall mutation load and NF1 mutation status, and showed significant enrichment in specific melanoma subtypes (desmoplastic and lentigo maligna melanoma), sun exposed anatomical locations, male gender, ulceration and increasing age. Identical clusters were identified in BCC and MCC samples, including several previously annotated in melanoma. In conclusion, we identify a novel class of mutation in cutaneous malignancies which rather than being defined by a particular genomic position are related to a position in a transcription factor binding motif. These mutations are also independent of cell of origin. Clinicopathologic correlation and mutation analysis support an etiologic role for chronic UV irradiation. Citation Format: Andrew J. Colebatch, Leon di Stefano, Stephen Q. Wong, Anthony T. Papenfuss, Grant A. McArthur. Clustered somatic mutations are frequent in transcription factor binding motifs in melanoma and other cutaneous malignancies. [abstract]. In: Proceedings of the 107th Annual Meeting of the American Association for Cancer Research; 2016 Apr 16-20; New Orleans, LA. Philadelphia (PA): AACR; Cancer Res 2016;76(14 Suppl):Abstract nr 87.

Background The HOX gene clusters are thought to be highly conserved amongst mammals and other vertebrates, but the long non-coding RNAs have only been studied in detail in human and mouse. The sequencing of the kangaroo genome provides an opportunity to use comparative analyses to compare the HOX clusters of a mammal with a distinct body plan to those of other mammals. Results Here we report a comparative analysis of HOX gene clusters between an Australian marsupial of the kangaroo family and the eutherians. There was a strikingly high level of conservation of HOX gene sequence and structure and non-protein coding genes including the microRNAs miR-196a, miR-196b, miR-10a and miR-10b and the long non-coding RNAs HOTAIR, HOTAIRM1 and HOX A11AS that play critical roles in regulating gene expression and controlling development. By microRNA deep sequencing and comparative genomic analyses, two conserved microRNAs (miR-10a and miR-10b) were identified and one new candidate microRNA with typical hairpin precursor structure that is expressed in both fibroblasts and testes was found. The prediction of microRNA target analysis showed that several known microRNA targets, such as miR-10, miR-414 and miR-464, were found in the tammar HOX clusters. In addition, several novel and putative miRNAs were identified that originated from elsewhere in the tammar genome and that target the tammar HOXB and HOXD clusters. Conclusions This study confirms that the emergence of known long non-coding RNAs in the HOX clusters clearly predate the marsupial-eutherian divergence 160 Ma ago. It also identified a new potentially functional microRNA as well as conserved miRNAs. These non-coding RNAs may participate in the regulation of HOX genes to influence the body plan of this marsupial.

Background Bats are the natural reservoir host for a range of emerging and re-emerging viruses, including SARS-like coronaviruses, Ebola viruses, henipaviruses and Rabies viruses. However, the mechanisms responsible for the control of viral replication in bats are not understood and there is little information available on any aspect of antiviral immunity in bats. Massively parallel sequencing of the bat transcriptome provides the opportunity for rapid gene discovery. Although the genomes of one megabat and one microbat have now been sequenced to low coverage, no transcriptomic datasets have been reported from any bat species. In this study, we describe the immune transcriptome of the Australian flying fox, Pteropus alecto, providing an important resource for identification of genes involved in a range of activities including antiviral immunity. Results Towards understanding the adaptations that have allowed bats to coexist with viruses, we have de novo assembled transcriptome sequence from immune tissues and stimulated cells from P. alecto. We identified about 18,600 genes involved in a broad range of activities with the most highly expressed genes involved in cell growth and maintenance, enzyme activity, cellular components and metabolism and energy pathways. 3.5% of the bat transcribed genes corresponded to immune genes and a total of about 500 immune genes were identified, providing an overview of both innate and adaptive immunity. A small proportion of transcripts found no match with annotated sequences in any of the public databases and may represent bat-specific transcripts. Conclusions This study represents the first reported bat transcriptome dataset and provides a survey of expressed bat genes that complement existing bat genomic data. In addition, these data provide insight into genes relevant to the antiviral responses of bats, and form a basis for examining the roles of these molecules in immune response to viral infection.

Background Devil Facial Tumour Disease (DFTD) is a unique clonal cancer that threatens the world's largest carnivorous marsupial, the Tasmanian devil (Sarcophilus harrisii) with extinction. This transmissible cancer is passed between individual devils by cell implantation during social interactions. The tumour arose in a Schwann cell of a single devil over 15 years ago and since then has expanded clonally, without showing signs of replicative senescence; in stark contrast to a somatic cell that displays a finite capacity for replication, known as the “Hayflick limit”. Methodology/Principal Findings In the present study we investigate the role of telomere length, measured as Telomere Copy Number (TCN), and telomerase and shelterin gene expression, as well as telomerase activity in maintaining hyperproliferation of Devil Facial Tumour (DFT) cells. Our results show that DFT cells have short telomeres. DFTD TCN does not differ between geographic regions or between strains. However, TCN has increased over time. Unlimited cell proliferation is likely to have been achieved through the observed up-regulation of the catalytic subunit of telomerase (TERT) and concomitant activation of telomerase. Up-regulation of the central component of shelterin, the TRF1-intercating nuclear factor 2 (TINF2) provides DFT a mechanism for telomere length homeostasis. The higher expression of both TERT and TINF2 may also protect DFT cells from genomic instability and enhance tumour proliferation. Conclusions/Significance DFT cells appear to monitor and regulate the length of individual telomeres: i.e. shorter telomeres are elongated by up-regulation of telomerase-related genes; longer telomeres are protected from further elongation by members of the shelterin complex, which may explain the lack of spatial and strain variation in DFT telomere copy number. The observed longitudinal increase in gene expression in DFT tissue samples and telomerase activity in DFT cell lines might indicate a selection for more stable tumours with higher proliferative potential.

BACKGROUND: We present the genome sequence of the tammar wallaby, Macropus eugenii, which is a member of the kangaroo family and the first representative of the iconic hopping mammals that symbolize Australia to be sequenced. The tammar has many unusual biological characteristics, including the longest period of embryonic diapause of any mammal, extremely synchronized seasonal breeding and prolonged and sophisticated lactation within a well-defined pouch. Like other marsupials, it gives birth to highly altricial young, and has a small number of very large chromosomes, making it a valuable model for genomics, reproduction and development. RESULTS: The genome has been sequenced to 2 × coverage using Sanger sequencing, enhanced with additional next generation sequencing and the integration of extensive physical and linkage maps to build the genome assembly. We also sequenced the tammar transcriptome across many tissues and developmental time points. Our analyses of these data shed light on mammalian reproduction, development and genome evolution: there is innovation in reproductive and lactational genes, rapid evolution of germ cell genes, and incomplete, locus-specific X inactivation. We also observe novel retrotransposons and a highly rearranged major histocompatibility complex, with many class I genes located outside the complex. Novel microRNAs in the tammar HOX clusters uncover new potential mammalian HOX regulatory elements. CONCLUSIONS: Analyses of these resources enhance our understanding of marsupial gene evolution, identify marsupial-specific conserved non-coding elements and critical genes across a range of biological systems, including reproduction, development and immunity, and provide new insight into marsupial and mammalian biology and genome evolution.