Search

High-grade serous ovarian cancers (HGSOCs) with homologous recombination deficiency (HRD) are initially responsive to poly (ADP-ribose) polymerase inhibitors (PARPi), but resistance ultimately emerges. HGSOC with CCNE1 amplification (CCNE1 amp) are associated with resistance to PARPi and platinum treatments. High replication stress in HRD and CCNE1 amp HGSOC leads to increased reliance on checkpoint kinase 1 (CHK1), a key regulator of cell cycle progression and the replication stress response. Here, we investigated the anti-tumor activity of the potent, highly selective, orally bioavailable CHK1 inhibitor (CHK1i), SRA737, in both acquired PARPi-resistant BRCA1/2 mutant and CCNE1 amp HGSOC models. We demonstrated that SRA737 increased replication stress and induced subsequent cell death in vitro. SRA737 monotherapy in vivo prolonged survival in CCNE1 amp models, suggesting a potential biomarker for CHK1i therapy. Combination SRA737 and PARPi therapy increased tumor regression in both PARPi-resistant and CCNE1 amp patient-derived xenograft models, warranting further study in these HGSOC subgroups.

UNLABELLED: Ovarian carcinosarcoma (OCS) is an aggressive and rare tumor type with limited treatment options. OCS is hypothesized to develop via the combination theory, with a single progenitor resulting in carcinomatous and sarcomatous components, or alternatively via the conversion theory, with the sarcomatous component developing from the carcinomatous component through epithelial-to-mesenchymal transition (EMT). In this study, we analyzed DNA variants from isolated carcinoma and sarcoma components to show that OCS from 18 women is monoclonal. RNA sequencing indicated that the carcinoma components were more mesenchymal when compared with pure epithelial ovarian carcinomas, supporting the conversion theory and suggesting that EMT is important in the formation of these tumors. Preclinical OCS models were used to test the efficacy of microtubule-targeting drugs, including eribulin, which has previously been shown to reverse EMT characteristics in breast cancers and induce differentiation in sarcomas. Vinorelbine and eribulin more effectively inhibited OCS growth than standard-of-care platinum-based chemotherapy, and treatment with eribulin reduced mesenchymal characteristics and N-MYC expression in OCS patient-derived xenografts. Eribulin treatment resulted in an accumulation of intracellular cholesterol in OCS cells, which triggered a downregulation of the mevalonate pathway and prevented further cholesterol biosynthesis. Finally, eribulin increased expression of genes related to immune activation and increased the intratumoral accumulation of CD8+ T cells, supporting exploration of immunotherapy combinations in the clinic. Together, these data indicate that EMT plays a key role in OCS tumorigenesis and support the conversion theory for OCS histogenesis. Targeting EMT using eribulin could help improve OCS patient outcomes. SIGNIFICANCE: Genomic analyses and preclinical models of ovarian carcinosarcoma support the conversion theory for disease development and indica

Targeting cell division by chemotherapy is a highly effective strategy to treat a wide range of cancers. However, there are limitations of many standard-of-care chemotherapies: undesirable drug toxicity, side-effects, resistance and high cost. New small molecules which kill a wide range of cancer subtypes, with good therapeutic window in vivo, have the potential to complement the current arsenal of anti-cancer agents and deliver improved safety profiles for cancer patients. We describe results with a new anti-cancer small molecule, WEHI-7326, which causes cell cycle arrest in G2/M, cell death in vitro, and displays efficacious anti-tumor activity in vivo. WEHI-7326 induces cell death in a broad range of cancer cell lines, including taxane-resistant cells, and inhibits growth of human colon, brain, lung, prostate and breast tumors in mice xenografts. Importantly, the compound elicits tumor responses as a single agent in patient-derived xenografts of clinically aggressive, treatment-refractory neuroblastoma, breast, lung and ovarian cancer. In combination with standard-of-care, WEHI-7326 induces a remarkable complete response in a mouse model of high-risk neuroblastoma. WEHI-7326 is mechanistically distinct from known microtubule-targeting agents and blocks cells early in mitosis to inhibit cell division, ultimately leading to apoptotic cell death. The compound is simple to produce and possesses favorable pharmacokinetic and toxicity profiles in rodents. It represents a novel class of anti-cancer therapeutics with excellent potential for further development due to the ease of synthesis, simple formulation, moderate side effects and potent in vivo activity. WEHI-7326 has the potential to complement current frontline anti-cancer drugs and to overcome drug resistance in a wide range of cancers.

Acquired PARP inhibitor (PARPi) resistance in BRCA1- or BRCA2-mutant ovarian cancer often results from secondary mutations that restore expression of functional protein. RAD51C is a less commonly studied ovarian cancer susceptibility gene whose promoter is sometimes methylated, leading to homologous recombination (HR) deficiency and PARPi sensitivity. For this study, the PARPi-sensitive patient-derived ovarian cancer xenograft PH039, which lacks HR gene mutations but harbors RAD51C promoter methylation, was selected for PARPi resistance by cyclical niraparib treatment in vivo. PH039 acquired PARPi resistance by the third treatment cycle and grew through subsequent treatment with either niraparib or rucaparib. Transcriptional profiling throughout the course of resistance development showed widespread pathway level changes along with a marked increase in RAD51C mRNA, which reflected loss of RAD51C promoter methylation. Analysis of ovarian cancer samples from the ARIEL2 Part 1 clinical trial of rucaparib monotherapy likewise indicated an association between loss of RAD51C methylation prior to on-study biopsy and limited response. Interestingly, the PARPi resistant PH039 model remained platinum sensitive. Collectively, these results not only indicate that PARPi treatment pressure can reverse RAD51C methylation and restore RAD51C expression, but also provide a model for studying the clinical observation that PARPi and platinum sensitivity are sometimes dissociated.

OBJECTIVE: Mucinous ovarian carcinoma (MOC) is an uncommon ovarian cancer histotype that responds poorly to conventional chemotherapy regimens. Although long overall survival outcomes can occur with early detection and optimal surgical resection, recurrent and advanced disease are associated with extremely poor survival. There are no current guidelines specifically for the systemic management of recurrent MOC. We analyzed data from a large cohort of women with MOC to evaluate the potential for clinical utility from a range of systemic agents. METHODS: We analyzed gene copy number (n = 191) and DNA sequencing data (n = 184) from primary MOC to evaluate signatures of mismatch repair deficiency and homologous recombination deficiency, and other genetic events. Immunohistochemistry data were collated for ER, CK7, CK20, CDX2, HER2, PAX8 and p16 (n = 117-166). RESULTS: Molecular aberrations noted in MOC that suggest a match with current targeted therapies include amplification of ERBB2 (26.7%) and BRAF mutation (9%). Observed genetic events that suggest potential efficacy for agents currently in clinical trials include: KRAS/NRAS mutations (66%), TP53 missense mutation (49%), RNF43 mutation (11%), ARID1A mutation (10%), and PIK3CA/PTEN mutation (9%). Therapies exploiting homologous recombination deficiency (HRD) may not be effective in MOC, as only 1/191 had a high HRD score. Mismatch repair deficiency was similarly rare (1/184). CONCLUSIONS: Although genetically diverse, MOC has several potential therapeutic targets. Importantly, the lack of response to platinum-based therapy observed clinically corresponds to the lack of a genomic signature associated with HRD, and MOC are thus also unlikely to respond to PARP inhibition.

Epithelial ovarian cancer (EOC) is a complex disease comprising discrete histological and molecular subtypes, for which survival rates remain unacceptably low. Tailored approaches for this deadly heterogeneous disease are urgently needed. Efflux pumps belonging to the ATP-binding cassette (ABC) family of transporters are known for roles in both drug resistance and cancer biology and are also highly targetable. Here we have investigated the association of ABCC4/MRP4 expression to clinical outcome and its biological function in endometrioid and serous tumors, common histological subtypes of EOC. We found high expression of ABCC4/MRP4, previously shown to be directly regulated by c-Myc/N-Myc, was associated with poor prognosis in endometrioid EOC (P = .001) as well as in a subset of serous EOC with a "high-MYCN" profile (C5/proliferative; P = .019). Transient siRNA-mediated suppression of MRP4 in EOC cells led to reduced growth, migration and invasion, with the effects being most pronounced in endometrioid and C5-like serous cells compared to non-C5 serous EOC cells. Sustained knockdown of MRP4 also sensitized endometrioid cells to MRP4 substrate drugs. Furthermore, suppression of MRP4 decreased the growth of patient-derived EOC cells in vivo. Together, our findings provide the first evidence that MRP4 plays an important role in the biology of Myc-associated ovarian tumors and highlight this transporter as a potential therapeutic target for EOC.

Acquired resistance to PARP inhibitors (PARPi) is a major challenge for the clinical management of high grade serous ovarian cancer (HGSOC). Here, we demonstrate CX-5461, the first-in-class inhibitor of RNA polymerase I transcription of ribosomal RNA genes (rDNA), induces replication stress and activates the DNA damage response. CX-5461 co-operates with PARPi in exacerbating replication stress and enhances therapeutic efficacy against homologous recombination (HR) DNA repair-deficient HGSOC-patient-derived xenograft (PDX) in vivo. We demonstrate CX-5461 has a different sensitivity spectrum to PARPi involving MRE11-dependent degradation of replication forks. Importantly, CX-5461 exhibits in vivo single agent efficacy in a HGSOC-PDX with reduced sensitivity to PARPi by overcoming replication fork protection. Further, we identify CX-5461-sensitivity gene expression signatures in primary and relapsed HGSOC. We propose CX-5461 is a promising therapy in combination with PARPi in HR-deficient HGSOC and also as a single agent for the treatment of relapsed disease.

Gene editing using clustered regularly interspaced short palindromic repeats/Cas9 has great potential for improving the compatibility of porcine organs with human recipients. However, the risk of detrimental off-target mutations in gene-edited pigs remains largely undefined. We have previously generated GGTA1 knock-in pigs for xenotransplantation using FokI-dCas9, a variant of Cas9 that is reported to reduce the frequency of off-target mutagenesis. In this study, we used whole genome sequencing (WGS) and optimized bioinformatic analysis to assess the fidelity of FokI-dCas9 editing in the generation of these pigs. Genomic DNA was isolated from porcine cells before and after gene editing and sequenced by WGS. The genomic sequences were analyzed using GRIDSS variant-calling software to detect putative structural variations (SVs), which were validated by PCR of DNA from knock-in and wild-type pigs. Platypus variant-calling software was used to detect single-nucleotide variations (SNVs) and small insertions/deletions (indels). GRIDSS analysis confirmed the precise integration of one copy of the knock-in construct in the gene-edited cells. Three additional SVs were detected by GRIDSS: deletions in intergenic regions in chromosome 6 and the X chromosome and a duplication of part of the CALD1 gene on chromosome 18. These mutations were not associated with plausible off-target sites, and were not detected in a second line of knock-in pigs generated using the same pair of guide RNAs, suggesting that they were the result of background mutation rather than off-target activity. Platypus identified 1375 SNVs/indels after quality filtering, but none of these were located in proximity to potential off-target sites, indicating that they were probably also spontaneous mutations. This is the first WGS analysis of pigs generated from FokI-dCas9-edited cells. Our results demonstrate that FokI-dCas9 is capable of high-fidelity gene editing with negligible off-target or undesired on-targe

High-grade mucinous ovarian cancer (HGMOC) is often a misnomer as the majority of cases are metastatic disease with a gastro-intestinal origin. The standard platinum-based ovarian cancer (OC) chemotherapy regimens are often ineffective, and there are insufficient data to support the use of colorectal cancer (CRC) chemotherapy regimens due to the rarity of HGMOC. We described a cohort of four consecutive suspected HGMOC cases treated at the Royal Women’s Hospital, Melbourne in 2012. Two cases were treated as primary MOC, whereas the other two were considered to be metastatic CRC based on histopathological and clinical evidence. From the RNAseq analysis, we identified two cases of HGMOC whose gene expression profiles were consistent with mucinous epithelial OC, one case that was treated as metastatic CRC with gene expression profile correlated with CRC and one case with neuroendocrine (NET) gene expression features. Interestingly, glucagon was over-expressed in this tumor that was subsequently confirmed by immunohistochemistry. These findings suggest a rare glucagonoma-like NET appendiceal tumor that had metastasized to the surface of ovary and were unresponsive to CRC chemotherapy regimens. In summary, a carefully curated panel of expression markers and selected functional genomics could provide diagnosis and treatment guidance for patients with possible HGMOC.

Introduction: PARP inhibitors (PARPi) have revolutionized disease management of patients with homologous recombination (HR) DNA repair-deficient high-grade serous ovarian cancer (HGSOC). However, acquired resistance to PARPi is a major challenge in the clinic. The specific inhibitor of RNA polymerase I (Pol I) transcription of ribosomal RNA genes (rDNA) has demonstrated single-agent antitumor activity in p53 wild-type and p53-mutant hematologic malignancies (first-in-human trial, dose escalation study of CX-5461 at Peter MacCallum Cancer Centre) (Khot et al., Cancer Discov 2019). CX-5461 has also been reported to exhibit synthetic lethality with BRCA1/2 deficiency through stabilization of G-quadruplex DNA (GQ) structures. Here, we investigate the efficacy of CX-5461 in treating HGSOC. Experimental Design: The mechanisms by which CX-5461 induces DNA damage response (DDR) and displays synthetic lethality in HR-deficient HGSOC cells are explored. We present in vivo data of mice bearing two functionally and genomically profiled HGSOC-patient-derived xenograft (PDX)s treated with CX-5461 and olaparib, alone and in combination. We also investigate CX-5461-sensitivity gene expression signatures in primary and relapsed HGSOC. Results: Utilizing ovarian cancer cell lines, we demonstrate that sensitivity to CX-5461 is associated with “BRCA1 mutation” and “MYC targets” gene expression signatures. In addition, sensitivity to CX-5461 is associated with high basal rates of Pol I transcription. Importantly, we demonstrate a novel mechanism for CX-5461 synthetic lethal interaction with HR deficiency mediated through the induction of replication stress at rDNA repeats. Our data reveal CX-5461-mediated DDR in HR-deficient cells does not involve stabilization of GQ structures as previously proposed. On the contrary, we show definitively that CX-5461 inhibits Pol I recruitment leading to rDNA chromatin defects including stabilization of R-loops, single-stranded DNA, and replic

Adrenal chromaffin cells and sympathetic neurons synthesize and release catecholamines, and both cell types are derived from neural crest precursors. However, they have different developmental histories, with sympathetic neurons derived directly from neural crest precursors while adrenal chromaffin cells arise from neural crest-derived cells that express Schwann cell markers. We have sought to identify the genes, including imprinted genes, which regulate the development of the two cell types in mice. We developed a method of separating the two cell types as early as E12.5, using differences in expression of enhanced yellow fluorescent protein driven from the tyrosine hydroxylase gene, and then used RNA sequencing to confirm the characteristic molecular signatures of the two cell types. We identified genes differentially expressed by adrenal chromaffin cells and sympathetic neurons. Deletion of a gene highly expressed by adrenal chromaffin cells, NIK-related kinase, a gene on the X-chromosome, results in reduced expression of adrenaline-synthesizing enzyme, phenyl-N-methyl transferase, by adrenal chromaffin cells and changes in cell cycle dynamics. Finally, many imprinted genes are up-regulated in chromaffin cells and may play key roles in their development.

Mucinous ovarian carcinoma (MOC) is a unique subtype of ovarian cancer with an uncertain etiology, including whether it genuinely arises at the ovary or is metastatic disease from other organs. In addition, the molecular drivers of invasive progression, high-grade and metastatic disease are poorly defined. We perform genetic analysis of MOC across all histological grades, including benign and borderline mucinous ovarian tumors, and compare these to tumors from other potential extra-ovarian sites of origin. Here we show that MOC is distinct from tumors from other sites and supports a progressive model of evolution from borderline precursors to high-grade invasive MOC. Key drivers of progression identified are TP53 mutation and copy number aberrations, including a notable amplicon on 9p13. High copy number aberration burden is associated with worse prognosis in MOC. Our data conclusively demonstrate that MOC arise from benign and borderline precursors at the ovary and are not extra-ovarian metastases.

Accurately identifying patients with high-grade serous ovarian carcinoma (HGSOC) who respond to poly(ADP-ribose) polymerase inhibitor (PARPi) therapy is of great clinical importance. Here we show that quantitative BRCA1 methylation analysis provides new insight into PARPi response in preclinical models and ovarian cancer patients. The response of 12 HGSOC patient-derived xenografts (PDX) to the PARPi rucaparib was assessed, with variable dose-dependent responses observed in chemo-naive BRCA1/2-mutated PDX, and no responses in PDX lacking DNA repair pathway defects. Among BRCA1-methylated PDX, silencing of all BRCA1 copies predicts rucaparib response, whilst heterozygous methylation is associated with resistance. Analysis of 21 BRCA1-methylated platinum-sensitive recurrent HGSOC (ARIEL2 Part 1 trial) confirmed that homozygous or hemizygous BRCA1 methylation predicts rucaparib clinical response, and that methylation loss can occur after exposure to chemotherapy. Accordingly, quantitative BRCA1 methylation analysis in a pre-treatment biopsy could allow identification of patients most likely to benefit, and facilitate tailoring of PARPi therapy.

The transcriptome of embryonic and adult rat lateral ventricular choroid plexus, using a combination of RNA-Sequencing and microarray data, was analyzed by functional groups of influx transporters, particularly solute carrier (SLC) transporters. RNA-Seq was performed at embryonic day (E) 15 and adult with additional data obtained at intermediate ages from microarray analysis. The largest represented functional group in the embryo was amino acid transporters (twelve) with expression levels 2-98 times greater than in the adult. In contrast, in the adult only six amino acid transporters were up-regulated compared to the embryo and at more modest enrichment levels (<5-fold enrichment above E15). In E15 plexus five glucose transporters, in particular Glut-1, and only one monocarboxylate transporter were enriched compared to the adult, whereas only two glucose transporters but six monocarboxylate transporters in the adult plexus were expressed at higher levels than in embryos. These results are compared with earlier published physiological studies of amino acid and monocarboxylate transport in developing rodents. This comparison shows correlation of high expression of some transporters in the developing brain with higher amino acid transport activity reported previously. Data for divalent metal transporters are also considered. Immunohistochemistry of several transporters (e.g., Slc16a10, a thyroid hormone transporter) gene products was carried out to confirm translational activity and to define cellular distribution of the proteins. Overall the results show that there is substantial expression of numerous influx transporters in the embryonic choroid plexus, many at higher levels than in the adult. This, together with immunohistochemical evidence and data from published physiological transport studies suggests that the choroid plexus in embryonic brain plays a major role in supplying the developing brain with essential nutrients.

We provide comprehensive identification of embryonic (E15) and adult rat lateral ventricular choroid plexus transcriptome, with focus on junction-associated proteins, ionic influx transporters and channels. Additionally, these data are related to new structural and previously published permeability studies. Results reveal that most genes associated with intercellular junctions are expressed at similar levels at both ages. In total, 32 molecules known to be associated with brain barrier interfaces were identified. Nine claudins showed unaltered expression, while two claudins (6 and 8) were expressed at higher levels in the embryo. Expression levels for most cytoplasmic/regulatory adaptors (10 of 12) were similar at the two ages. A few junctional genes displayed lower expression in embryos, including 5 claudins, occludin and one junctional adhesion molecule. Three gap junction genes were enriched in the embryo. The functional effectiveness of these junctions was assessed using blood-delivered water-soluble tracers at both the light and electron microscopic level: embryo and adult junctions halted movement of both 286Da and 3kDa molecules into the cerebrospinal fluid (CSF). The molecular identities of many ion channel and transporter genes previously reported as important for CSF formation and secretion in the adult were demonstrated in the embryonic choroid plexus (and validated with immunohistochemistry of protein products), but with some major age-related differences in expression. In addition, a large number of previously unidentified ion channel and transporter genes were identified for the first time in plexus epithelium. These results, in addition to data obtained from electron microscopical and physiological permeability experiments in immature brains, indicate that exchange between blood and CSF is mainly transcellular, as well-formed tight junctions restrict movement of small water-soluble molecules from early in development. These data strongly indicate the

The phosphate-specific transport operon, pstSCAB-phoU, of Gram-negative bacteria is an essential part of the Pho regulon. Its key roles are to encode a high-affinity inorganic phosphate transport system and to prevent activation of PhoB in phosphate-rich environments. In general, mutations in pstSCAB-phoU lead to the constitutive expression of the Pho regulon. Previously, we constructed a pstCA deletion mutant of Citrobacter rodentium and found it to be attenuated for virulence in mice, its natural host. This attenuation was dependent on PhoB or PhoB-regulated gene(s) because a phoB mutation restored virulence for mice to the pstCA mutant. To investigate how downstream genes may contribute to the virulence of C. rodentium, we used microarray analysis to investigate global gene expression of C. rodentium strain ICC169 and its isogenic pstCA mutant when grown in phosphate-rich medium. Overall 323 genes of the pstCA mutant were differentially expressed by at least 1.5-fold compared to the wild-type C. rodentium. Of these 145 were up-regulated and 178 were down-regulated. Differentially expressed genes included some involved in phosphate homoeostasis, cellular metabolism and protein metabolism. A large number of genes involved in stress responses and of unknown function were also differentially expressed, as were some virulence-associated genes. Up-regulated virulence-associated genes in the pstCA mutant included that for DegP, a serine protease, which appeared to be directly regulated by PhoB. Down-regulated genes included those for the production of the urease, flagella, NleG8 (a type III-secreted protein) and the tad focus (which encodes type IVb pili in Yersinia enterocolitica). Infection studies using C57/BL6 mice showed that DegP and NleG8 play a role in bacterial virulence. Overall, our study provides evidence that Pho is a global regulator of gene expression in C. rodentium and indicates the presence of at least two previously unrecognized virulence determinants

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. © 2011 Renfree et al.; licensee BioMed Central Ltd.

Transcriptional control is dependent on a vast network of epigenetic modifications. One epigenetic mark of particular interest is tri-methylation of lysine 27 on histone H3 (H3K27me3), which is catalysed and maintained by Polycomb Repressive Complex 2 (PRC2). Although this histone mark is studied widely, the precise relationship between its local pattern of enrichment and regulation of gene expression is currently unclear. We have used ChIP-seq to generate genome-wide maps of H3K27me3 enrichment, and have identified three enrichment profiles with distinct regulatory consequences. First, a broad domain of H3K27me3 enrichment across the body of genes corresponds to the canonical view of H3K27me3 as inhibitory to transcription. Second, a peak of enrichment around the transcription start site (TSS) is commonly associated with 'bivalent' genes, where H3K4me3 also marks the TSS. Finally and most surprisingly, we identified an enrichment profile with a peak in the promoter of genes that is associated with active transcription. Genes with each of these three profiles were found in different proportions in each of the cell types studied. The data analysis techniques developed here will be useful for the identification of common enrichment profiles for other histone modifications that have important consequences for transcriptional regulation.

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 e

We present GOseq, an application for performing Gene Ontology (GO) analysis on RNA-seq data. GO analysis is widely used to reduce complexity and highlight biological processes in genome-wide expression studies, but standard methods give biased results on RNA-seq data due to over-detection of differential expression for long and highly expressed transcripts. Application of GOseq to a prostate cancer data set shows that GOseq dramatically changes the results, highlighting categories more consistent with the known biology.

BACKGROUND: Several recent studies have demonstrated the effectiveness of deep sequencing for transcriptome analysis (RNA-seq) in mammals. As RNA-seq becomes more affordable, whole genome transcriptional profiling is likely to become the platform of choice for species with good genomic sequences. As yet, a rigorous analysis methodology has not been developed and we are still in the stages of exploring the features of the data. RESULTS: We investigated the effect of transcript length bias in RNA-seq data using three different published data sets. For standard analyses using aggregated tag counts for each gene, the ability to call differentially expressed genes between samples is strongly associated with the length of the transcript. CONCLUSION: Transcript length bias for calling differentially expressed genes is a general feature of current protocols for RNA-seq technology. This has implications for the ranking of differentially expressed genes, and in particular may introduce bias in gene set testing for pathway analysis and other multi-gene systems biology analyses. REVIEWERS: This article was reviewed by Rohan Williams (nominated by Gavin Huttley), Nicole Cloonan (nominated by Mark Ragan) and James Bullard (nominated by Sandrine Dudoit).

We have implemented in Python the COmparative GENomic Toolkit, a fully integrated and thoroughly tested framework for novel probabilistic analyses of biological sequences, devising workflows, and generating publication quality graphics. PyCogent includes connectors to remote databases, built-in generalized probabilistic techniques for working with biological sequences, and controllers for third-party applications. The toolkit takes advantage of parallel architectures and runs on a range of hardware and operating systems, and is available under the general public license from http://sourceforge.net/projects/pycogent.

A report on 'The Biology of Genomes' meeting, Cold Spring Harbor, USA, 10-14 May 2006.

The first sequenced marsupial genome promises to reveal unparalleled insights into mammalian evolution. We have used the Monodelphis domestica (gray short-tailed opossum) sequence to construct the first map of a marsupial major histocompatibility complex (MHC). The MHC is the most gene-dense region of the mammalian genome and is critical to immunity and reproductive success. The marsupial MHC bridges the phylogenetic gap between the complex MHC of eutherian mammals and the minimal essential MHC of birds. Here we show that the opossum MHC is gene dense and complex, as in humans, but shares more organizational features with non-mammals. The Class I genes have amplified within the Class II region, resulting in a unique Class I/II region. We present a model of the organization of the MHC in ancestral mammals and its elaboration during mammalian evolution. The opossum genome, together with other extant genomes, reveals the existence of an ancestral "immune supercomplex" that contained genes of both types of natural killer receptors together with antigen processing genes and MHC genes.

The first sequenced marsupial genome promises to reveal unparalleled insights into mammalian evolution. We have used the Monodelphis domestica (gray short-tailed opossum) sequence to construct the first map of a marsupial major histocompatibility complex (MHC). The MHC is the most gene-dense region of the mammalian genome and is critical to immunity and reproductive success. The marsupial MHC bridges the phylogenetic gap between the complex MHC of eutherian mammals and the minimal essential MHC of birds. Here we show that the opossum MHC is gene dense and complex, as in humans, but shares more organizational features with non-mammals. The Class I genes have amplified within the Class II region, resulting in a unique Class I/II region. We present a model of the organization of the MHC in ancestral mammals and its elaboration during mammalian evolution. The opossum genome, together with other extant genomes, reveals the existence of an ancestral "immune supercomplex" that contained genes of both types of natural killer receptors together with antigen processing genes and MHC genes. Copyright: © 2006 Belov et al.

A recent landmark paper demonstrates the unique contribution of marsupials and monotremes to comparative genome analysis, filling an evolutionary gap between the eutherian mammals (including humans) and more distant vertebrate species.

BACKGROUND: Phylogenetic footprinting is the identification of functional regions of DNA by their evolutionary conservation. This is achieved by comparing orthologous regions from multiple species and identifying the DNA regions that have diverged less than neutral DNA. Vestige is a phylogenetic footprinting package built on the PyEvolve toolkit that uses probabilistic molecular evolutionary modelling to represent aspects of sequence evolution, including the conventional divergence measure employed by other footprinting approaches. In addition to measuring the divergence, Vestige allows the expansion of the definition of a phylogenetic footprint to include variation in the distribution of any molecular evolutionary processes. This is achieved by displaying the distribution of model parameters that represent partitions of molecular evolutionary substitutions. Examination of the spatial incidence of these effects across regions of the genome can identify DNA segments that differ in the nature of the evolutionary process. RESULTS: Vestige was applied to a reference dataset of the SCL locus from four species and provided clear identification of the known conserved regions in this dataset. To demonstrate the flexibility to use diverse models of molecular evolution and dissect the nature of the evolutionary process Vestige was used to footprint the Ka/Ks ratio in primate BRCA1 with a codon model of evolution. Two regions of putative adaptive evolution were identified illustrating the ability of Vestige to represent the spatial distribution of distinct molecular evolutionary processes. CONCLUSION: Vestige provides a flexible, open platform for phylogenetic footprinting. Underpinned by the PyEvolve toolkit, Vestige provides a framework for visualising the signatures of evolutionary processes across the genome of numerous organisms simultaneously. By exploiting the maximum-likelihood statistical framework, the complex interplay between mutational processes, DNA repair and se

BACKGROUND: Examining the distribution of variation has proven an extremely profitable technique in the effort to identify sequences of biological significance. Most approaches in the field, however, evaluate only the conserved portions of sequences - ignoring the biological significance of sequence differences. A suite of sophisticated likelihood based statistical models from the field of molecular evolution provides the basis for extracting the information from the full distribution of sequence variation. The number of different problems to which phylogeny-based maximum likelihood calculations can be applied is extensive. Available software packages that can perform likelihood calculations suffer from a lack of flexibility and scalability, or employ error-prone approaches to model parameterisation. RESULTS: Here we describe the implementation of PyEvolve, a toolkit for the application of existing, and development of new, statistical methods for molecular evolution. We present the object architecture and design schema of PyEvolve, which includes an adaptable multi-level parallelisation schema. The approach for defining new methods is illustrated by implementing a novel dinucleotide model of substitution that includes a parameter for mutation of methylated CpG's, which required 8 lines of standard Python code to define. Benchmarking was performed using either a dinucleotide or codon substitution model applied to an alignment of BRCA1 sequences from 20 mammals, or a 10 species subset. Up to five-fold parallel performance gains over serial were recorded. Compared to leading alternative software, PyEvolve exhibited significantly better real world performance for parameter rich models with a large data set, reducing the time required for optimisation from approximately 10 days to approximately 6 hours. CONCLUSION: PyEvolve provides flexible functionality that can be used either for statistical modelling of molecular evolution, or the development of new methods in the fi