The impact of telomere length on prostate cancer aggressiveness, genomic instability and health disparities.

The telomere repetitive TTAGGG motif at the ends of chromosomes, serves to preserve genomic integrity and chromosomal stability. In turn, genomic instability is a hallmark of cancer-implicating telomere disturbance. Prostate cancer (PCa) shows significant ancestral disparities, with men of African ancestry at the greatest risk for aggressive disease and associated genomic instability. Yet, no study has explored the role of telomere length (TL) with respect to ancestrally driven PCa health disparities. Patient- and technically-matched tumour-blood whole genome sequencing data for 179 ancestrally defined treatment naïve PCa patients (117 African, 62 European), we assessed for TL (blood and tumour) associations. We found shortened tumour TL to be associated with aggressive PCa presentation and elevated genomic instabilities, including percentage of genome alteration and copy number gains, in men of African ancestry. For European patients, tumour TL showed significant associations with PCa driver genes PTEN, TP53, MSH2, SETBP1 and DDX11L1, while shorter blood TL (< 3200 base pairs) and tumour TL (< 2861 base pairs) were correlated with higher risk for biochemical recurrence. Concurring with previous studies linking TL to PCa diagnosis and/or prognosis, for the first time we correlated TL differences with patient ancestry with important implications for future treatments targeting telomere dysfunction.

ANO7 African-ancestral genomic diversity and advanced prostate cancer.

BACKGROUND: Prostate cancer (PCa) is a significant health burden for African men, with mortality rates more than double global averages. The prostate specific Anoctamin 7 (ANO7) gene linked with poor patient outcomes has recently been identified as the target for an African-specific protein-truncating PCa-risk allele. METHODS: Here we determined the role of ANO7 in a study of 889 men from southern Africa, leveraging exomic genotyping array PCa case-control data (n = 780, 17 ANO7 alleles) and deep sequenced whole genome data for germline and tumour ANO7 interrogation (n = 109), while providing clinicopathologically matched European-derived sequence data comparative analyses (n = 57). Associated predicted deleterious variants (PDVs) were further assessed for impact using computational protein structure analysis. RESULTS: Notably rare in European patients, we found the common African PDV p.Ile740Leu (rs74804606) to be associated with PCa risk in our case-control analysis (Wilcoxon rank-sum test, false discovery rate/FDR = 0.03), while sequencing revealed co-occurrence with the recently reported African-specific deleterious risk variant p.Ser914* (rs60985508). Additional findings included a novel protein-truncating African-specific frameshift variant p.Asp789Leu, African-relevant PDVs associated with altered protein structure at Ca2+ binding sites, early-onset PCa associated with PDVs and germline structural variants in Africans (Linear regression models, -6.42 years, 95% CI = -10.68 to -2.16, P-value = 0.003) and ANO7 as an inter-chromosomal PCa-related gene fusion partner in African derived tumours. CONCLUSIONS: Here we provide not only validation for ANO7 as an African-relevant protein-altering PCa-risk locus, but additional evidence for a role of inherited and acquired ANO7 variance in the observed phenotypic heterogeneity and African-ancestral health disparity.

Alterations in the Epigenetic Machinery Associated with Prostate Cancer Health Disparities.

Prostate cancer is driven by acquired genetic alterations, including those impacting the epigenetic machinery. With African ancestry as a significant risk factor for aggressive disease, we hypothesize that dysregulation among the roughly 656 epigenetic genes may contribute to prostate cancer health disparities. Investigating prostate tumor genomic data from 109 men of southern African and 56 men of European Australian ancestry, we found that African-derived tumors present with a longer tail of epigenetic driver gene candidates (72 versus 10). Biased towards African-specific drivers (63 versus 9 shared), many are novel to prostate cancer (18/63), including several putative therapeutic targets (CHD7, DPF3, POLR1B, SETD1B, UBTF, and VPS72). Through clustering of all variant types and copy number alterations, we describe two epigenetic PCa taxonomies capable of differentiating patients by ancestry and predicted clinical outcomes. We identified the top genes in African- and European-derived tumors representing a multifunctional "generic machinery", the alteration of which may be instrumental in epigenetic dysregulation and prostate tumorigenesis. In conclusion, numerous somatic alterations in the epigenetic machinery drive prostate carcinogenesis, but African-derived tumors appear to achieve this state with greater diversity among such alterations. The greater novelty observed in African-derived tumors illustrates the significant clinical benefit to be derived from a much needed African-tailored approach to prostate cancer healthcare aimed at reducing prostate cancer health disparities.

Evaluating Germline Testing Panels in Southern African Males With Advanced Prostate Cancer.

BACKGROUND: Germline testing for prostate cancer is on the increase, with clinical implications for risk assessment, treatment, and management. Regardless of family history, NCCN recommends germline testing for patients with metastatic, regional, very-high-risk localized, and high-risk localized prostate cancer. Although African ancestry is a significant risk factor for aggressive prostate cancer, due to a lack of available data no testing criteria have been established for ethnic minorities. PATIENTS AND METHODS: Through deep sequencing, we interrogated the 20 most common germline testing panel genes in 113 Black South African males presenting with largely advanced prostate cancer. Bioinformatic tools were then used to identify the pathogenicity of the variants. RESULTS: After we identified 39 predicted deleterious variants (16 genes), further computational annotation classified 17 variants as potentially oncogenic (12 genes; 17.7% of patients). Rare pathogenic variants included CHEK2 Arg95Ter, BRCA2 Trp31Arg, ATM Arg3047Ter (2 patients), and TP53 Arg282Trp. Notable oncogenic variants of unknown pathogenicity included novel BRCA2 Leu3038Ile in a patient with early-onset disease, whereas patients with FANCA Arg504Cys and RAD51C Arg260Gln reported a family history of prostate cancer. Overall, rare pathogenic and early-onset or familial-associated oncogenic variants were identified in 6.9% (5/72) and 9.2% (8/87) of patients presenting with a Gleason score ≥8 or ≥4 + 3 prostate cancer, respectively. CONCLUSIONS: In this first-of-its-kind study of southern African males, we provide support of African inclusion for advanced, early-onset, and familial prostate cancer genetic testing, indicating clinical value for 30% of current gene panels. Recognizing current panel limitations highlights an urgent need to establish testing guidelines for men of African ancestry. We provide a rationale for considering lowering the pathologic diagnostic inclusion criteria and call for further genome-wide interrogation to ensure the best possible African-relevant prostate cancer gene panel.

African-specific molecular taxonomy of prostate cancer.

Prostate cancer is characterized by considerable geo-ethnic disparity. African ancestry is a significant risk factor, with mortality rates across sub-Saharan Africa of 2.7-fold higher than global averages1. The contributing genetic and non-genetic factors, and associated mutational processes, are unknown2,3. Here, through whole-genome sequencing of treatment-naive prostate cancer samples from 183 ancestrally (African versus European) and globally distinct patients, we generate a large cancer genomics resource for sub-Saharan Africa, identifying around 2 million somatic variants. Significant African-ancestry-specific findings include an elevated tumour mutational burden, increased percentage of genome alteration, a greater number of predicted damaging mutations and a higher total of mutational signatures, and the driver genes NCOA2, STK19, DDX11L1, PCAT1 and SETBP1. Examining all somatic mutational types, we describe a molecular taxonomy for prostate cancer differentiated by ancestry and defined as global mutational subtypes (GMS). By further including Chinese Asian data, we confirm that GMS-B (copy-number gain) and GMS-D (mutationally noisy) are specific to African populations, GMS-A (mutationally quiet) is universal (all ethnicities) and the African-European-restricted subtype GMS-C (copy-number losses) predicts poor clinical outcomes. In addition to the clinical benefit of including individuals of African ancestry, our GMS subtypes reveal different evolutionary trajectories and mutational processes suggesting that both common genetic and environmental factors contribute to the disparity between ethnicities. Analogous to gene-environment interaction-defined here as a different effect of an environmental surrounding in people with different ancestries or vice versa-we anticipate that GMS subtypes act as a proxy for intrinsic and extrinsic mutational processes in cancers, promoting global inclusion in landmark studies.

Genome-wide interrogation of structural variation reveals novel African-specific prostate cancer oncogenic drivers.

BACKGROUND: African ancestry is a significant risk factor for advanced prostate cancer (PCa). Mortality rates in sub-Saharan Africa are 2.5-fold greater than global averages. However, the region has largely been excluded from the benefits of whole genome interrogation studies. Additionally, while structural variation (SV) is highly prevalent, PCa genomic studies are still biased towards small variant interrogation. METHODS: Using whole genome sequencing and best practice workflows, we performed a comprehensive analysis of SVs for 180 (predominantly Gleason score ≥ 8) prostate tumours derived from 115 African, 61 European and four ancestrally admixed patients. We investigated the landscape and relationship of somatic SVs in driving ethnic disparity (African versus European), with a focus on African men from southern Africa. RESULTS: Duplication events showed the greatest ethnic disparity, with a 1.6- (relative frequency) to 2.5-fold (count) increase in African-derived tumours. Furthermore, we found duplication events to be associated with CDK12 inactivation and MYC copy number gain, and deletion events associated with SPOP mutation. Overall, African-derived tumours were 2-fold more likely to present with a hyper-SV subtype. In addition to hyper-duplication and deletion subtypes, we describe a new hyper-translocation subtype. While we confirm a lower TMPRSS2-ERG fusion-positive rate in tumours from African cases (10% versus 33%), novel African-specific PCa ETS family member and TMPRSS2 fusion partners were identified, including LINC01525, FBXO7, GTF3C2, NTNG1 and YPEL5. Notably, we found 74 somatic SV hotspots impacting 18 new candidate driver genes, with CADM2, LSAMP, PTPRD, PDE4D and PACRG having therapeutic implications for African patients. CONCLUSIONS: In this first African-inclusive SV study for high-risk PCa, we demonstrate the power of SV interrogation for the identification of novel subtypes, oncogenic drivers and therapeutic targets. Identifying a novel spectrum of SVs in tumours derived from African patients provides a mechanism that may contribute, at least in part, to the observed ethnic disparity in advanced PCa presentation in men of African ancestry.

DNA methylation profiles unique to Kalahari KhoeSan individuals.

Genomes of KhoeSan individuals of the Kalahari Desert provide the greatest understanding of single nucleotide diversity in the human genome. Compared with individuals in industrialized environments, the KhoeSan have a unique foraging and hunting lifestyle. Given these dramatic environmental differences, and the responsiveness of the methylome to environmental exposures of many types, we hypothesized that DNA methylation patterns would differ between KhoeSan and neighbouring agropastoral and/or industrial Bantu. We analysed Illumina HumanMethylation 450 k array data generated from blood samples from 38 KhoeSan and 42 Bantu, and 6 Europeans. After removing CpG positions associated with annotated and novel polymorphisms and controlling for white blood cell composition, sex, age and technical variation we identified 816 differentially methylated CpG loci, out of which 133 had an absolute beta-value difference of at least 0.05. Notably SLC39A4/ZIP4, which plays a role in zinc transport, was one of the most differentially methylated loci. Although the chronological ages of the KhoeSan are not formally recorded, we compared historically estimated ages to methylation-based calculations. This study demonstrates that the epigenetic profile of KhoeSan individuals reveals differences from other populations, and along with extensive genetic diversity, this community brings increased accessibility and understanding to the diversity of the human genome.

The Impact of Whole Genome Data on Therapeutic Decision-Making in Metastatic Prostate Cancer: A Retrospective Analysis.

BACKGROUND: While critical insights have been gained from evaluating the genomic landscape of metastatic prostate cancer, utilizing this information to inform personalized treatment is in its infancy. We performed a retrospective pilot study to assess the current impact of precision medicine for locally advanced and metastatic prostate adenocarcinoma and evaluate how genomic data could be harnessed to individualize treatment. METHODS: Deep whole genome-sequencing was performed on 16 tumour-blood pairs from 13 prostate cancer patients; whole genome optical mapping was performed in a subset of 9 patients to further identify large structural variants. Tumour samples were derived from prostate, lymph nodes, bone and brain. RESULTS: Most samples had acquired genomic alterations in multiple therapeutically relevant pathways, including DNA damage response (11/13 cases), PI3K (7/13), MAPK (10/13) and Wnt (9/13). Five patients had somatic copy number losses in genes that may indicate sensitivity to immunotherapy (LRP1B, CDK12, MLH1) and one patient had germline and somatic BRCA2 alterations. CONCLUSIONS: Most cases, whether primary or metastatic, harboured therapeutically relevant alterations, including those associated with PARP inhibitor sensitivity, immunotherapy sensitivity and resistance to androgen pathway targeting agents. The observed intra-patient heterogeneity and presence of genomic alterations in multiple growth pathways in individual cases suggests that a precision medicine model in prostate cancer needs to simultaneously incorporate multiple pathway-targeting agents. Our whole genome approach allowed for structural variant assessment in addition to the ability to rapidly reassess an individual's molecular landscape as knowledge of relevant biomarkers evolve. This retrospective oncological assessment highlights the genomic complexity of prostate cancer and the potential impact of assessing genomic data for an individual at any stage of the disease.

MethylToSNP: identifying SNPs in Illumina DNA methylation array data.

BACKGROUND: Current array-based methods for the measurement of DNA methylation rely on the process of sodium bisulfite conversion to differentiate between methylated and unmethylated cytosine bases in DNA. In the absence of genotype data this process can lead to ambiguity in data interpretation when a sample has polymorphisms at a methylation probe site. A common way to minimize this problem is to exclude such potentially problematic sites, with some methods removing as much as 60% of array probes from consideration before data analysis. RESULTS: Here, we present an algorithm implemented in an R Bioconductor package, MethylToSNP, which detects a characteristic data pattern to infer sites likely to be confounded by polymorphisms. Additionally, the tool provides a stringent reliability score to allow thresholding on SNP predictions. We calibrated parameters and thresholds used by the algorithm on simulated and real methylation data sets. We illustrate findings using methylation data from YRI (Yoruba in Ibadan, Nigeria), CEPH (European descent) and KhoeSan (southern African) populations. Our polymorphism predictions made using MethylToSNP have been validated through SNP databases and bisulfite and genomic sequencing. CONCLUSIONS: The benefits of this method are threefold. First, it prevents extensive data loss by considering only SNPs specific to the individuals in the study. Second, it offers the possibility to identify new polymorphisms in samples for which there is little known about the genetic landscape. Third, it identifies variants as they exist in functional regions of a genome, such as in CTCF (transcriptional repressor) sites and enhancers, that may be common alleles or personal mutations with potential to deleteriously affect genomic regulatory activities. We demonstrate that MethylToSNP is applicable to the Illumina 450K and Illumina 850K EPIC array data and is also backwards compatible to the 27K methylation arrays. Going forward, this kind of nuanced approach can increase the amount of information derived from precious data sets by considering samples of the project individually to enable more informed decisions about data cleaning.

Metagenomic analysis reveals a rich bacterial content in high-risk prostate tumors from African men.

BACKGROUND: Inflammation is a hallmark of prostate cancer (PCa), yet no pathogenic agent has been identified. Men from Africa are at increased risk for both aggressive prostate disease and infection. We hypothesize that pathogenic microbes may be contributing, at least in part, to high-risk PCa presentation within Africa and in turn the observed ethnic disparity. METHODS: Here we reveal through metagenomic analysis of host-derived whole-genome sequencing data, the microbial content within prostate tumor tissue from 22 men. What is unique about this study is that patients were separated by ethnicity, African vs European, and environments, Africa vs Australia. RESULTS: We identified 23 common bacterial genera between the African, Australian, and Chinese prostate tumor samples, while nonbacterial microbes were notably absent. While the most abundant genera across all samples included: Escherichia, Propionibacterium, and Pseudomonas, the core prostate tumor microbiota was enriched for Proteobacteria. We observed a significant increase in the richness of the bacterial communities within the African vs Australian samples (t = 4.6-5.5; P = .0004-.001), largely driven by eight predominant genera. Considering core human gut microbiota, African prostate tissue samples appear enriched for Escherichia and Acidovorax, with an abundance of Eubacterium associated with host tumor hypermutation. CONCLUSIONS: Our study provides suggestive evidence for the presence of a core, bacteria-rich, prostate microbiome. While unable to exclude for fecal contamination, the observed increased bacterial content and richness within the African vs non-African samples, together with elevated tumor mutational burden, suggests the possibility that bacterially-driven oncogenic transformation within the prostate microenvironment may be contributing to aggressive disease presentation in Africa.

African KhoeSan ancestry linked to high-risk prostate cancer.

BACKGROUNDS: Genetic diversity is greatest within Africa, in particular the KhoeSan click-speaking peoples of southern Africa. South African populations represent admixture fractions including differing degrees of African, African-KhoeSan and non-African genetic ancestries. Within the United States, African ancestry has been linked to prostate cancer presentation and mortality. Together with environmental contributions, genetics is a significant risk factor for high-risk prostate cancer, defined by a pathological Gleason score ≥ 8. METHODS: Using genotype array data merged with ancestry informative reference data, we investigate the contribution of African ancestral fractions to high-risk prostate cancer. Our study includes 152 South African men of African (Black) or African-admixed (Coloured) ancestries, in which 40% showed high-risk prostate cancer. RESULTS: Genetic fractions were determined for averaging an equal African to non-African genetic ancestral contribution in the Coloured; we found African ancestry to be linked to high-risk prostate cancer (P-value = 0.0477). Adjusting for age, the associated African ancestral fraction was driven by a significant KhoeSan over Bantu contribution, defined by Gleason score ≥ 8 (P-value = 0.02329) or prostate specific antigen levels ≥20 ng/ml (P-value = 0.03713). Additionally, we observed the mean overall KhoeSan contribution to be increased in Black patients with high-risk (11.8%) over low-risk (10.9%) disease. Linking for the first time KhoeSan ancestry to a common modern disease, namely high-risk prostate cancer, we tested in this small study the validity of using KhoeSan ancestry as a surrogate for identifying potential high-risk prostate cancer risk loci. As such, we identified four loci within chromosomal regions 2p11.2, 3p14, 8q23 and 22q13.2 (P-value = all age-adjusted < 0.01), two of which have previously been associated with high-risk prostate cancer. CONCLUSIONS: Our study suggests that ancient KhoeSan ancestry may be linked to common modern diseases, specifically those of late onset and therefore unlikely to have undergone exclusive selective pressure. As such we show within a uniquely admixed South African population a link between KhoeSan ancestry and high-risk prostate cancer, which may explain the 2-fold increase in presentation in Black South Africans compared with African Americans.

TMPRSS2-ERG fusions linked to prostate cancer racial health disparities: A focus on Africa.

BACKGROUND: The androgen-regulated gene TMPRSS2 to the ETS transcription factor gene ERG fusion is the most common genomic alteration acquired during prostate tumorigenesis and biased toward men of European ancestry. In contrast, African American men present with more advanced disease, yet their tumors are less likely to acquire TMPRSS2-ERG. Data for Africa is scarce. METHODS: RNA was made available for genomic analyses from 181 prostate tissue biopsy cores from Black South African men, 94 with and 87 without pathological evidence for prostate cancer. Reverse transcription polymerase chain reaction was used to screen for the TMPRSS2-ERG fusion, while transcript junction coordinates and isoform frequencies, including novel gene fusions, were determined using targeted RNA sequencing. RESULTS: Here we report a frequency of 13% for TMPRSS2-ERG in tumors from Black South Africans. Present in 12/94 positive versus 1/87 cancer negative prostate tissue cores, this suggests a 92.62% predictivity for a positive cancer diagnosis (P = 0.0031). At a frequency of almost half that reported for African Americans and roughly a quarter of that reported for men of European ancestry, acquisition of TMPRSS2-ERG appears to be inversely associated with aggressive prostate cancer. Further support was provided by linking the presence of TMPRSS2-ERG to low-grade disease in younger patients (P = 0.0466), with higher expressing distal ERG fusion junction coordinates. CONCLUSIONS: Only the second study of its kind for the African continent, we support a link between TMPRSS2-ERG status and prostate cancer racial health disparity beyond the borders of the United States. We call for urgent evaluation of androgen deprivation therapy within Africa.

Whole-Genome Sequencing Reveals Elevated Tumor Mutational Burden and Initiating Driver Mutations in African Men with Treatment-Naïve, High-Risk Prostate Cancer.

: African-American men are more likely than any other racial group to die from prostate cancer. The contribution of acquired genomic variation to this racial disparity is largely unknown, as genomic from Africa is lacking. Here, we performed the first tumor-normal paired deep whole-genome sequencing for Africa. A direct study-matched comparison between African- and European-derived, treatment-naïve, high-risk prostate tumors for 15 cases allowed for further comparative analyses of existing data. Excluding a single hypermutated tumor with 55 mutations per megabase, we observed a 1.8-fold increase in small somatic variants in African- versus European-derived tumors (P = 1.02e-04), rising to 4-fold when compared with published tumor-matched data. Furthermore, we observed an increase in oncogenic driver mutations in African tumors (P = 2.92e-03); roughly 30% of impacted genes were novel to prostate cancer, and 79% of recurrent driver mutations appeared early in tumorigenesis. Although complex genomic rearrangements were less frequent in African tumors, we describe a uniquely hyperduplicated tumor affecting 149 transposable elements. Comparable with African Americans, ERG fusions and PIK3CA mutations were absent and PTEN loss less frequent. CCND1 and MYC were frequently gained, with somatic copy-number changes more likely to occur late in tumorigenesis. In addition to traditional prostate cancer gene pathways, genes regulating calcium ion-ATPase signal transduction were disrupted in African tumors. Although preliminary, our results suggest that further validation and investigation into the potential implications for elevated tumor mutational burden and tumor-initiating mutations in clinically unfavorable prostate cancer can improve patient outcomes in Africa. SIGNIFICANCE: The first whole-genome sequencing study for high-risk prostate cancer in African men allows a simultaneous comparison of ethnic differences relative to European populations and of the influences of the environment relative to African-American men. GRAPHICAL ABSTRACT: http://cancerres.aacrjournals.org/content/canres/78/24/6736/F1.large.jpg.See related commentary by Huang, p. 6726.

Altered mitochondrial genome content signals worse pathology and prognosis in prostate cancer.

BACKGROUND: Mitochondrial genome (mtDNA) content is depleted in many cancers. In prostate cancer, there is intra-glandular as well as inter-patient mtDNA copy number variation. In this study, we determine if mtDNA content can be used as a predictor for prostate cancer staging and outcomes. METHODS: Fresh prostate cancer biopsies from 115 patients were obtained at time of surgery. All cores underwent pathological review, followed by isolation of cancer and normal tissue. DNA was extracted and qPCR performed to quantify the total amount of mtDNA as a ratio to genomic DNA. Differences in mtDNA content were compared for prostate cancer pathology features and disease outcomes. RESULTS: We showed a significantly reduced mtDNA content in prostate cancer compared with normal adjacent prostate tissue (mean difference 1.73-fold, P-value <0.001). Prostate cancer with increased mtDNA content showed unfavorable pathologic characteristics including, higher disease stage (PT2 vs PT3 P-value = 0.018), extracapsular extension (P-value = 0.02) and a trend toward an increased Gleason score (P-value = 0.064). No significant association was observed between changes in mtDNA content and biochemical recurrence (median follow up of 107 months). CONCLUSIONS: Contrary to other cancer types, prostate cancer tissue shows no universally depleted mtDNA content. Rather, the change in mtDNA content is highly variable, mirroring known prostate cancer genome heterogeneity. Patients with high mtDNA content have an unfavorable pathology, while a high mtDNA content in normal adjacent prostate tissue is associated with worse prognosis.

Next generation mapping reveals novel large genomic rearrangements in prostate cancer.

Complex genomic rearrangements are common molecular events driving prostate carcinogenesis. Clinical significance, however, has yet to be fully elucidated. Detecting the full range and subtypes of large structural variants (SVs), greater than one kilobase in length, is challenging using clinically feasible next generation sequencing (NGS) technologies. Next generation mapping (NGM) is a new technology that allows for the interrogation of megabase length DNA molecules outside the detection range of single-base resolution NGS. In this study, we sought to determine the feasibility of using the Irys (Bionano Genomics Inc.) nanochannel NGM technology to generate whole genome maps of a primary prostate tumor and matched blood from a Gleason score 7 (4 + 3), ETS-fusion negative prostate cancer patient. With an effective mapped coverage of 35X and sequence coverage of 60X, and an estimated 43% tumor purity, we identified 85 large somatic structural rearrangements and 6,172 smaller somatic variants, respectively. The vast majority of the large SVs (89%), of which 73% are insertions, were not detectable ab initio using high-coverage short-read NGS. However, guided manual inspection of single NGS reads and de novo assembled scaffolds of NGM-derived candidate regions allowed for confirmation of 94% of these large SVs, with over a third impacting genes with oncogenic potential. From this single-patient study, the first cancer study to integrate NGS and NGM data, we hypothesise that there exists a novel spectrum of large genomic rearrangements in prostate cancer, that these large genomic rearrangements are likely early events in tumorigenesis, and they have potential to enhance taxonomy.

Mutational load of the mitochondrial genome predicts pathological features and biochemical recurrence in prostate cancer.

Prostate cancer management is complicated by extreme disease heterogeneity, which is further limited by availability of prognostic biomarkers. Recognition of prostate cancer as a genetic disease has prompted a focus on the nuclear genome for biomarker discovery, with little attention given to the mitochondrial genome. While it is evident that mitochondrial DNA (mtDNA) mutations are acquired during prostate tumorigenesis, no study has evaluated the prognostic value of mtDNA variation. Here we used next-generation sequencing to interrogate the mitochondrial genomes from prostate tissue biopsies and matched blood of 115 men having undergone a radical prostatectomy for which there was a mean of 107 months clinical follow-up. We identified 74 unique prostate cancer specific somatic mtDNA variants in 50 patients, providing significant expansion to the growing catalog of prostate cancer mtDNA mutations. While no single variant or variant cluster showed recurrence across multiple patients, we observe a significant positive correlation between the total burden of acquired mtDNA variation and elevated Gleason Score at diagnosis and biochemical relapse. We add to accumulating evidence that total acquired genomic burden, rather than specific mtDNA mutations, has diagnostic value. This is the first study to demonstrate the prognostic potential of mtDNA mutational burden in prostate cancer.

Comparative genome analyses reveal distinct structure in the saltwater crocodile MHC.

The major histocompatibility complex (MHC) is a dynamic genome region with an essential role in the adaptive immunity of vertebrates, especially antigen presentation. The MHC is generally divided into subregions (classes I, II and III) containing genes of similar function across species, but with different gene number and organisation. Crocodylia (crocodilians) are widely distributed and represent an evolutionary distinct group among higher vertebrates, but the genomic organisation of MHC within this lineage has been largely unexplored. Here, we studied the MHC region of the saltwater crocodile (Crocodylus porosus) and compared it with that of other taxa. We characterised genomic clusters encompassing MHC class I and class II genes in the saltwater crocodile based on sequencing of bacterial artificial chromosomes. Six gene clusters spanning ∼452 kb were identified to contain nine MHC class I genes, six MHC class II genes, three TAP genes, and a TRIM gene. These MHC class I and class II genes were in separate scaffold regions and were greater in length (2-6 times longer) than their counterparts in well-studied fowl B loci, suggesting that the compaction of avian MHC occurred after the crocodilian-avian split. Comparative analyses between the saltwater crocodile MHC and that from the alligator and gharial showed large syntenic areas (>80% identity) with similar gene order. Comparisons with other vertebrates showed that the saltwater crocodile had MHC class I genes located along with TAP, consistent with birds studied. Linkage between MHC class I and TRIM39 observed in the saltwater crocodile resembled MHC in eutherians compared, but absent in avian MHC, suggesting that the saltwater crocodile MHC appears to have gene organisation intermediate between these two lineages. These observations suggest that the structure of the saltwater crocodile MHC, and other crocodilians, can help determine the MHC that was present in the ancestors of archosaurs.

Three crocodilian genomes reveal ancestral patterns of evolution among archosaurs.

To provide context for the diversification of archosaurs--the group that includes crocodilians, dinosaurs, and birds--we generated draft genomes of three crocodilians: Alligator mississippiensis (the American alligator), Crocodylus porosus (the saltwater crocodile), and Gavialis gangeticus (the Indian gharial). We observed an exceptionally slow rate of genome evolution within crocodilians at all levels, including nucleotide substitutions, indels, transposable element content and movement, gene family evolution, and chromosomal synteny. When placed within the context of related taxa including birds and turtles, this suggests that the common ancestor of all of these taxa also exhibited slow genome evolution and that the comparatively rapid evolution is derived in birds. The data also provided the opportunity to analyze heterozygosity in crocodilians, which indicates a likely reduction in population size for all three taxa through the Pleistocene. Finally, these data combined with newly published bird genomes allowed us to reconstruct the partial genome of the common ancestor of archosaurs, thereby providing a tool to investigate the genetic starting material of crocodilians, birds, and dinosaurs.