A Pseudolikelihood Approach for Simultaneous Analysis of Array Comparative Genomic Hybridizations (aCGH)

DNA sequence copy number has been shown to be associated with cancer development and progression. Array-based Comparative Genomic Hybridization (aCGH) is a recent development that seeks to identify the copy number ratio at large numbers of markers across the genome. Due to experimental and biological variations across chromosomes and across hybridizations, current methods are limited to analyses of single chromosomes. We propose a more powerful approach that borrows strength across chromosomes and across hybridizations. We assume a Gaussian mixture model, with a hidden Markov dependence structure, and with random effects to allow for intertumoral variation, as well as intratumoral clonal variation. For ease of computation, we base estimation on a pseudolikelihood function. The method produces quantitative assessments of the likelihood of genetic alterations at each clone, along with a graphical display for simple visual interpretation. We assess the characteristics of the method through simulation studies and through analysis of a brain tumor aCGH data set. We show that the pseudolikelihood approach is superior to existing methods both in detecting small regions of copy number alteration and in accurately classifying regions of change when intratumoral clonal variation is present.

Large genomic fibrillin-1 (FBN1) gene deletions provide evidence for true haploinsufficiency in Marfan syndrome

Mutations in the FBN1 gene are the major cause of Marfan syndrome (MFS), an autosomal dominant connective tissue disorder, which displays variable manifestations in the cardiovascular, ocular, and skeletal systems. Current molecular genetic testing of FBN1 may miss mutations in the promoter region or in other noncoding sequences as well as partial or complete gene deletions and duplications. In this study, we tested for copy number variations by successively applying multiplex ligation-dependent probe amplification (MLPA) and the Affymetrix Human Mapping 500 K Array Set, which contains probes for approximately 500,000 single-nucleotide polymorphisms (SNPs) across the genome. By analyzing genomic DNA of 101 unrelated individuals with MFS or related phenotypes in whom standard genetic testing detected no mutation, we identified FBN1 deletions in two patients with MFS. Our high-resolution approach narrowed down the deletion breakpoints. Subsequent sequencing of the junctional fragments revealed the deletion sizes of 26,887 and 302,580 bp, respectively. Surprisingly, both deletions affect the putative regulatory and promoter region of the FBN1 gene, strongly indicating that they abolish transcription of the deleted allele. This expectation of complete loss of function of one allele, i.e. true haploinsufficiency, was confirmed by transcript analyses. Our findings not only emphasize the importance of screening for large genomic rearrangements in comprehensive genetic testing of FBN1 but, importantly, also extend the molecular etiology of MFS by providing hitherto unreported evidence that true haploinsufficiency is sufficient to cause MFS.

PHYLOGENETIC ANALYSES, GENE DOSAGE EFFECT AND SELECTION IN MEMBERS OF THE LATERAL BOUNDARY DOMAIN (LBD) GENES TRANSCRIPTION FACTOR FAMILY

Wood formation is an economically and environmentally important process and has played a significant role in the evolution of terrestrial plants. Despite its significance, the molecular underpinnings of the process are still poorly understood. We have previously shown that four Lateral Boundary Domain (LBD) transcription factors have important roles in the regulation of wood formation with two (LBD1 and LBD4) involved in secondary phloem and ray cell development and two (LBD15 and LBD18) in secondary xylem formation. Here, we used comparative phylogenetic analyses to test potential roles of the four LBD genes in the evolution of woodiness. We studied the copy number and variation in DNA and amino acid sequences of the four LBDs in a wide range of woody and herbaceous plant taxa with fully sequenced and annotated genomes. LBD1 showed the highest gene copy number across the studied species, and LBD1 gene copy number was strongly and significantly correlated with the level of ray seriation. The lianas, cucumber and grape, with multiseriate ray cells showed the highest gene copy number (12 and 11, respectively). Because lianas’ growth habit requires significant twisting and bending, the less lignified ray parenchyma cells likely facilitate stem flexibility and maintenance of xylem conductivity. We further demonstrate conservation of amino acids in the LBD18 protein sequences that are specific to woody taxa. Neutrality tests showed evidence for strong purifying selection on these gene regions across various orders, indicating adaptive convergent evolution of LBD18. Structural modeling demonstrates that the conserved amino acids have a significant impact on the tertiary protein structure and thus are likely of significant functional importance.

Quantitative sequence analysis of FBN1 premature termination codons provides evidence for incomplete NMD in leukocytes

We improved, evaluated, and used Sanger sequencing for quantification of single nucleotide polymorphism (SNP) variants in transcripts and gDNA samples. This improved assay resulted in highly reproducible relative allele frequencies (e.g., for a heterozygous gDNA 50.0+/-1.4%, and for a missense mutation-bearing transcript 46.9+/-3.7%) with a lower detection limit of 3-9%. It provided excellent accuracy and linear correlation between expected and observed relative allele frequencies. This sequencing assay, which can also be used for the quantification of copy number variations (CNVs), methylations, mosaicisms, and DNA pools, enabled us to analyze transcripts of the FBN1 gene in fibroblasts and blood samples of patients with suspected Marfan syndrome not only qualitatively but also quantitatively. We report a total of 18 novel and 19 known FBN1 sequence variants leading to a premature termination codon (PTC), 26 of which we analyzed by quantitative sequencing both at gDNA and cDNA levels. The relative amounts of PTC-containing FBN1 transcripts in fresh and PAXgene-stabilized blood samples were significantly higher (33.0+/-3.9% to 80.0+/-7.2%) than those detected in affected fibroblasts with inhibition of nonsense-mediated mRNA decay (NMD) (11.0+/-2.1% to 25.0+/-1.8%), whereas in fibroblasts without NMD inhibition no mutant alleles could be detected. These results provide evidence for incomplete NMD in leukocytes and have particular importance for RNA-based analyses not only in FBN1 but also in other genes.

RNA-dependent genome processing during nuclear differentiation: the model systems of stichotrichous ciliates

We introduce ciliated protozoa, and more specifically the stichotrichous ciliates Oxytricha and Stylonychia, as biological model systems for the analysis of programmed DNA-reorganization processes during nuclear differentiation. These include DNA excision, DNA elimination, reordering of gene segments and specific gene amplification. We show that small nuclear RNAs specify DNA sequences to be excised or retained, but also discuss the need for a RNA template molecule derived from the parental nucleus for these processes. This RNA template guides reordering of gene segments to become functional genes and determines gene copy number in the differentiated nucleus. Since the template is derived from the parental macronucleus, gene reordering and DNA amplification are inherited in a non-Mendelian epigenetic manner.

Trypanosomal TAC40 constitutes a novel subclass of mitochondrial -barrel proteins specialized in mitochondrial genome inheritance

Mitochondria cannot form de novo but require mechanisms allowing their inheritance to daughter cells. In contrast to most other eukaryotes Trypanosoma brucei has a single mitochondrion whose single-unit genome is physically connected to the flagellum. Here we identify a β-barrel mitochondrial outer membrane protein, termed tripartite attachment complex 40 (TAC40), that localizes to this connection. TAC40 is essential for mitochondrial DNA inheritance and belongs to the mitochondrial porin protein family. However, it is not specifically related to any of the three subclasses of mitochondrial porins represented by the metabolite transporter voltage-dependent anion channel (VDAC), the protein translocator of the outer membrane 40 (TOM40), or the fungi-specific MDM10, a component of the endoplasmic reticulum–mitochondria encounter structure (ERMES). MDM10 and TAC40 mediate cellular architecture and participate in transmembrane complexes that are essential for mitochondrial DNA inheritance. In yeast MDM10, in the context of the ERMES, is postulated to connect the mitochondrial genomes to actin filaments, whereas in trypanosomes TAC40 mediates the linkage of the mitochondrial DNA to the basal body of the flagellum. However, TAC40 does not colocalize with trypanosomal orthologs of ERMES components and, unlike MDM10, it regulates neither mitochondrial morphology nor the assembly of the protein translocase. TAC40 therefore defines a novel subclass of mitochondrial porins that is distinct from VDAC, TOM40, and MDM10. However, whereas the architecture of the TAC40-containing complex in trypanosomes and the MDM10-containing ERMES in yeast is very different, both are organized around a β-barrel protein of the mitochondrial porin family that mediates a DNA–cytoskeleton linkage that is essential for mitochondrial DNA inheritance.

The Draft Assembly of the Radically Organized Stylonychia lemnae Macronuclear Genome.

Stylonychia lemnae is a classical model single-celled eukaryote, and a quintessential ciliate typified by dimorphic nuclei: A small, germline micronucleus and a massive, vegetative macronucleus. The genome within Stylonychia's macronucleus has a very unusual architecture, comprised variably and highly amplified "nanochromosomes," each usually encoding a single gene with a minimal amount of surrounding noncoding DNA. As only a tiny fraction of the Stylonychia genes has been sequenced, and to promote research using this organism, we sequenced its macronuclear genome. We report the analysis of the 50.2-Mb draft S. lemnae macronuclear genome assembly, containing in excess of 16,000 complete nanochromosomes, assembled as less than 20,000 contigs. We found considerable conservation of fundamental genomic properties between S. lemnae and its close relative, Oxytricha trifallax, including nanochromosomal gene synteny, alternative fragmentation, and copy number. Protein domain searches in Stylonychia revealed two new telomere-binding protein homologs and the presence of linker histones. Among the diverse histone variants of S. lemnae and O. trifallax, we found divergent, coexpressed variants corresponding to four of the five core nucleosomal proteins (H1.2, H2A.6, H2B.4, and H3.7) suggesting that these ciliates may possess specialized nucleosomes involved in genome processing during nuclear differentiation. The assembly of the S. lemnae macronuclear genome demonstrates that largely complete, well-assembled highly fragmented genomes of similar size and complexity may be produced from one library and lane of Illumina HiSeq 2000 shotgun sequencing. The provision of the S. lemnae macronuclear genome sets the stage for future detailed experimental studies of chromatin-mediated, RNA-guided developmental genome rearrangements.

A conserved mitochondrial outer membrane protein mediates kDNA maintenace in Trypanosoma brucei

Kinetoplastids are defined by the unique organization of their mitochondrial DNA (kDNA). It forms a highly concatenated DNA network that is linked to the basal body of the flagellum by the tripartite attachment complex (TAC). The TAC encompasses intra and extramitochondrial filaments and a highly differentiated region of the two mitochondrial membranes. Here we identify and characterize a mitochondrial outer membrane protein of Trypanosoma brucei that is predominantly localized in the TAC. The protein is essential for growth in both life cycle stages. Immunofluorescence shows that ablation of the protein does not affect kDNA replication but abolishes the segregation of the replicated kDNA network causing rapid loss of kDNA. Besides its role in kDNA maintenance in vivo and in vitro experiments show that the protein is involved in mitochondrial protein import and that it interacts with a recently discovered protein import factor. RNAi experiments in a T. brucei cell line in which the kDNA is dispensable suggest that the essential function is linked to kDNA maintenance. Bioinformatic analysis shows that the studied outer membrane protein has beta-barrel topology and that it belongs to the mitochondrial porin family comprising VDAC, Tom40 and Mdm10. Interestingly, Mdm10 has sofar only been found in yeast. Ist function in protein import and mitochondrial DNA maintenance suggests that the protein in our study is the functional homologue of Mdm10. Thus, the TAC – a defining structure of Kinetoplastids – contains a conserved protein which in yeast and trypanosomes performs the same function. Our study therefore provides an example that trypanosomal biology, rather than being unique, often simply represents a more extreme manifestation of a conserved biological concept.

A conserved mitochondrial outer membrane protein mediates kDNA maintenace in Trypanosoma brucei

Kinetoplastids are defined by the unique organization of their mitochondrial DNA (kDNA). It forms a highly concatenated DNA network that is linked to the basal body of the flagellum by the tripartite attachment complex (TAC). The TAC encompasses intra and extramitochondrial filaments and a highly differentiated region of the two mitochondrial membranes. Here we identify and characterize a mitochondrial outer membrane protein of Trypanosoma brucei that is predominantly localized in the TAC. The protein is essential for growth in both life cycle stages. Immunofluorescence shows that ablation of the protein does not affect kDNA replication but abolishes the segregation of the replicated kDNA network causing rapid loss of kDNA. Besides its role in kDNA maintenance in vivo and in vitro experiments show that the protein is involved in mitochondrial protein import and that it interacts with a recently discovered protein import factor. RNAi experiments in a T. brucei cell line in which the kDNA is dispensable suggest that the essential function is linked to kDNA maintenance. Bioinformatic analysis shows that the studied outer membrane protein has beta-barrel topology and that it belongs to the mitochondrial porin family comprising VDAC, Tom40 and Mdm10. Interestingly, Mdm10 has so far only been found in yeast. Its function in protein import and mitochondrial DNA maintenance suggests that the protein in our study is the functional homologue of Mdm10. Thus, the TAC – a defining structure of Kinetoplastids – contains a conserved protein which in yeast and trypanosomes performs the same function. Our study therefore provides an example that trypanosomal biology, rather than being unique, often simply represents a more extreme manifestation of a conserved biological concept.

Tumour genomic and microenvironmental heterogeneity for integrated prediction of 5-year biochemical recurrence of prostate cancer: a retrospective cohort study.

BACKGROUND Clinical prognostic groupings for localised prostate cancers are imprecise, with 30-50% of patients recurring after image-guided radiotherapy or radical prostatectomy. We aimed to test combined genomic and microenvironmental indices in prostate cancer to improve risk stratification and complement clinical prognostic factors. METHODS We used DNA-based indices alone or in combination with intra-prostatic hypoxia measurements to develop four prognostic indices in 126 low-risk to intermediate-risk patients (Toronto cohort) who will receive image-guided radiotherapy. We validated these indices in two independent cohorts of 154 (Memorial Sloan Kettering Cancer Center cohort [MSKCC] cohort) and 117 (Cambridge cohort) radical prostatectomy specimens from low-risk to high-risk patients. We applied unsupervised and supervised machine learning techniques to the copy-number profiles of 126 pre-image-guided radiotherapy diagnostic biopsies to develop prognostic signatures. Our primary endpoint was the development of a set of prognostic measures capable of stratifying patients for risk of biochemical relapse 5 years after primary treatment. FINDINGS Biochemical relapse was associated with indices of tumour hypoxia, genomic instability, and genomic subtypes based on multivariate analyses. We identified four genomic subtypes for prostate cancer, which had different 5-year biochemical relapse-free survival. Genomic instability is prognostic for relapse in both image-guided radiotherapy (multivariate analysis hazard ratio [HR] 4·5 [95% CI 2·1-9·8]; p=0·00013; area under the receiver operator curve [AUC] 0·70 [95% CI 0·65-0·76]) and radical prostatectomy (4·0 [1·6-9·7]; p=0·0024; AUC 0·57 [0·52-0·61]) patients with prostate cancer, and its effect is magnified by intratumoral hypoxia (3·8 [1·2-12]; p=0·019; AUC 0·67 [0·61-0·73]). A novel 100-loci DNA signature accurately classified treatment outcome in the MSKCC low-risk to intermediate-risk cohort (multivariate analysis HR 6·1 [95% CI 2·0-19]; p=0·0015; AUC 0·74 [95% CI 0·65-0·83]). In the independent MSKCC and Cambridge cohorts, this signature identified low-risk to high-risk patients who were most likely to fail treatment within 18 months (combined cohorts multivariate analysis HR 2·9 [95% CI 1·4-6·0]; p=0·0039; AUC 0·68 [95% CI 0·63-0·73]), and was better at predicting biochemical relapse than 23 previously published RNA signatures. INTERPRETATION This is the first study of cancer outcome to integrate DNA-based and microenvironment-based failure indices to predict patient outcome. Patients exhibiting these aggressive features after biopsy should be entered into treatment intensification trials. FUNDING Movember Foundation, Prostate Cancer Canada, Ontario Institute for Cancer Research, Canadian Institute for Health Research, NIHR Cambridge Biomedical Research Centre, The University of Cambridge, Cancer Research UK, Cambridge Cancer Charity, Prostate Cancer UK, Hutchison Whampoa Limited, Terry Fox Research Institute, Princess Margaret Cancer Centre Foundation, PMH-Radiation Medicine Program Academic Enrichment Fund, Motorcycle Ride for Dad (Durham), Canadian Cancer Society.

Novel mitochondrial tRNA(Ile) m.4282A>G gene mutation leads to chronic progressive external ophthalmoplegia plus phenotype

BACKGROUND/AIM To investigate the underlying pathomechanism in a 33-year-old female Caucasian patient presenting with chronic progressive external ophthalmoplegia (CPEO) plus symptoms. METHODS Histochemical analysis of skeletal muscle and biochemical measurements of individual oxidative phosphorylation (OXPHOS) complexes. Genetic analysis of mitochondrial DNA in various tissues with subsequent investigation of single muscle fibres for correlation of mutational load. RESULTS The patient's skeletal muscle showed 20% of cytochrome c oxidase-negative fibres and 8% ragged-red fibres. Genetic analysis of the mitochondrial DNA revealed a novel point mutation in the mitochondrial tRNA(Ile) (MTTI) gene at position m.4282G>A. The heteroplasmy was determined in blood, buccal cells and muscle by restriction fragment length polymorphism (RFLP) combined with a last fluorescent cycle. The total mutational load was 38% in skeletal muscle, but was not detectable in blood or buccal cells of the patient. The phenotype segregated with the mutational load as determined by analysis of single cytochrome c oxidase-negative/positive fibres by laser capture microdissection and subsequent LFC-RFLP. CONCLUSIONS We describe a novel MTTI transition mutation at nucleotide position m.4282G>A associated with a CPEO plus phenotype. The novel variant at position m.4282G>A disrupts the middle bond of the D-stem of the tRNA(Ile) and is highly conserved. The conservation and phenotype-genotype segregation strongly suggest pathogenicity and is in good agreement with the MTTI gene being frequently associated with CPEO. This novel variant broadens the spectrum of MTTI mutations causing CPEO.

The TIM translocase in T. brucei contains Rhomboid-like proteins that are essential for mitochondrial protein import

The parasitic protozoon Trypanosoma brucei is often considered as one of the earliest branching eukaryotes that have mitochondria capable of oxidative phosphorylation. Its protein import systems are therefore of great interest. Recently, it was shown that the outer mitochondrial membrane protein translocase is of similar complexity yet different composition than in other eukaryotes (1). In the inner membrane however, only a single orthologue of the pore forming Tim17/22/23 protein family was identified and termed TbTim17. Based on this finding it has been suggested that, instead of separate TIM22 and TIM23 complexes as in other eukaryotes, trypanosomes may have a single multifunctional translocase of the inner mitochondrial membrane (TIM) of reduced complexity. To elucidate the composition of the trypanosomal TIM complex we performed co-immunoprecipitations (CoIP) of epitope-tagged TbTim17 in combination with SILAC-based quantitative mass spectrometry. This led to the identification of 22 highly enriched TbTim17-interacting proteins. We tagged two of the top-scoring proteins for reciprocal CoIP analyses and recovered a set of ten proteins that are highly enriched in all three CoIPs. These proteins are excellent candidates for core subunits of the trypanosomal TIM complex. Eight of them were present in the previously determined inner membrane proteome and four show homology to small Tim chaperones. Three candidates, a novel trypanosomatid-specific 42 kDa protein, termed Tim42, and two putative orthologues of probably inactive rhomboid proteases were chosen for further analysis. All three proteins are essential in both life cycle stages and in a cell line that can grow in the absence of mitochondrial DNA. Additionally, their ablation by RNAi results in a strong protein import defect both in vivo and in vitro. Blue native PAGE reveals that Tim42, like TbTim17 is present in a high molecular weight complex. Moreover, ablation of either Tim42 or TbTim17 leads to a destabilization of the complex containing the other protein, suggesting a tight interaction of the two proteins. In summary our study shows that unlike anticipated trypanosomes have a highly complex TIM translocase that has extensively been redesigned. We have characterized three novel TIM subunits that have never been associated with mitochondrial protein import before. Two of them belong to the rhomboid protease family, a member of which recently has been implicated in the ERAD translocation system. Our study provides insight into mitochondrial evolution over large phylogenetic distances and suggests an exciting analogy between protein translocation systems of mitochondria and the ER.

What can a unicellular parasite tell us about mitochondrial biogenesis?

In the unicellular parasite Trypanosoma brucei, as in other eukaryotes, more than 95% of all mitochondrial proteins are imported from the cytosol. The recently characterized multisubunit ATOM complex, the functional analogue of the TOM complex of yeast, mediates import of essentially all proteins across the outer mitochondrial membrane in T. brucei. Moreover, an additional protein termed pATOM36, which is loosely associated with the ATOM complex, has been implicated in the import of only a subset of mitochondrial proteins. Here we have investigated more precisely which role pATOM36 plays in mitochondrial protein import. RNAi mediated ablation of pATOM36 specifically depletes a subset of outer mitochondrial membrane proteins including ATOM complex subunits and as a consequence results in the collapse of the ATOM complex as shown by Blue native PAGE. In addition, a SILAC-based global proteomic analysis of uninduced and induced pATOM36 RNAi cells together with in vitro import experiments suggest that pATOM36 might be a novel protein import factor acting on a subset of alpha-helically anchored mitochondrial outer membrane proteins. Identification of pATOM36 interaction partners by co-immunoprecipitation together with immunofluorescence analysis shows that unexpectedly a fraction of the protein is associated with the tripartite attachment complex (TAC). This complex is essential for proper inheritance of the mitochondrial DNA in T. brucei. It forms a physical connection between the single unit mitochondrial DNA and the basal body of the flagellum that is stable throughout the cell cycle. Thus, pATOM36 simultaneously mediates ATOM assembly, and thus protein import, as well as mitochondrial DNA inheritance since it is an essential component of the TAC.

Spatial genomic heterogeneity within localized, multifocal prostate cancer

Herein we provide a detailed molecular analysis of the spatial heterogeneity of clinically localized, multifocal prostate cancer to delineate new oncogenes or tumor suppressors. We initially determined the copy number aberration (CNA) profiles of 74 patients with index tumors of Gleason score 7. Of these, 5 patients were subjected to whole-genome sequencing using DNA quantities achievable in diagnostic biopsies, with detailed spatial sampling of 23 distinct tumor regions to assess intraprostatic heterogeneity in focal genomics. Multifocal tumors are highly heterogeneous for single-nucleotide variants (SNVs), CNAs and genomic rearrangements. We identified and validated a new recurrent amplification of MYCL, which is associated with TP53 deletion and unique profiles of DNA damage and transcriptional dysregulation. Moreover, we demonstrate divergent tumor evolution in multifocal cancer and, in some cases, tumors of independent clonal origin. These data represent the first systematic relation of intraprostatic genomic heterogeneity to predicted clinical outcome and inform the development of novel biomarkers that reflect individual prognosis.

A certified plasmid reference material for the standardisation of BCR-ABL1 mRNA quantification by real-time quantitative PCR.

Serial quantification of BCR-ABL1 mRNA is an important therapeutic indicator in chronic myeloid leukaemia, but there is a substantial variation in results reported by different laboratories. To improve comparability, an internationally accepted plasmid certified reference material (CRM) was developed according to ISO Guide 34:2009. Fragments of BCR-ABL1 (e14a2 mRNA fusion), BCR and GUSB transcripts were amplified and cloned into pUC18 to yield plasmid pIRMM0099. Six different linearised plasmid solutions were produced with the following copy number concentrations, assigned by digital PCR, and expanded uncertainties: 1.08±0.13 × 10(6), 1.08±0.11 × 10(5), 1.03±0.10 × 10(4), 1.02±0.09 × 10(3), 1.04±0.10 × 10(2) and 10.0±1.5 copies/μl. The certification of the material for the number of specific DNA fragments per plasmid, copy number concentration of the plasmid solutions and the assessment of inter-unit heterogeneity and stability were performed according to ISO Guide 35:2006. Two suitability studies performed by 63 BCR-ABL1 testing laboratories demonstrated that this set of 6 plasmid CRMs can help to standardise a number of measured transcripts of e14a2 BCR-ABL1 and three control genes (ABL1, BCR and GUSB). The set of six plasmid CRMs is distributed worldwide by the Institute for Reference Materials and Measurements (Belgium) and its authorised distributors (https://ec.europa.eu/jrc/en/reference-materials/catalogue/; CRM code ERM-AD623a-f).