Clonality of circulating tumor cells in breast cancer brain metastasis patients.

BACKGROUND: The incidence of brain metastases in breast cancer (BCBM) patients is increasing. These patients have a very poor prognosis, and therefore, identification of blood-based biomarkers, such as circulating tumor cells (CTCs), and understanding the genomic heterogeneity could help to personalize treatment options. METHODS: Both EpCAM-dependent (CellSearch® System) and EpCAM-independent Ficoll-based density centrifugation methods were used to detect CTCs from 57 BCBM patients. DNA from individual CTCs and corresponding primary tumors and brain metastases were analyzed by next-generation sequencing (NGS) in order to evaluate copy number aberrations and single nucleotide variations (SNVs). RESULTS: CTCs were detected after EpCAM-dependent enrichment in 47.7% of the patients (≥ 5 CTCs/7.5 ml blood in 20.5%). The CTC count was associated with ERBB2 status (p = 0.029) of the primary tumor as well as with the prevalence of bone metastases (p = 0.021). EpCAM-independent enrichment revealed CTCs in 32.6% of the patients, especially among triple-negative breast cancer (TNBC) patients (70.0%). A positive CTC status after enrichment of either method was significantly associated with decreased overall survival time (p < 0.05). Combining the results of both enrichment methods, 63.6% of the patients were classified as CTC positive. In three patients, the matched tumor tissue and single CTCs were analyzed by NGS showing chromosomal aberrations with a high genomic clonality and mutations in pathways potentially important in brain metastasis formation. CONCLUSION: The detection of CTCs, regardless of the enrichment method, is of prognostic relevance in BCBM patients and in combination with molecular analysis of CTCs can help defining patients with higher risk of early relapse and suitability for targeted treatment.

Corrigendum: Common genetic variation drives molecular heterogeneity in human iPSCs.

This corrects the article DOI: 10.1038/nature22403.

Common genetic variation drives molecular heterogeneity in human iPSCs.

Technology utilizing human induced pluripotent stem cells (iPS cells) has enormous potential to provide improved cellular models of human disease. However, variable genetic and phenotypic characterization of many existing iPS cell lines limits their potential use for research and therapy. Here we describe the systematic generation, genotyping and phenotyping of 711 iPS cell lines derived from 301 healthy individuals by the Human Induced Pluripotent Stem Cells Initiative. Our study outlines the major sources of genetic and phenotypic variation in iPS cells and establishes their suitability as models of complex human traits and cancer. Through genome-wide profiling we find that 5-46% of the variation in different iPS cell phenotypes, including differentiation capacity and cellular morphology, arises from differences between individuals. Additionally, we assess the phenotypic consequences of genomic copy-number alterations that are repeatedly observed in iPS cells. In addition, we present a comprehensive map of common regulatory variants affecting the transcriptome of human pluripotent cells.

Tyrosine kinase 2 is not limiting human antiviral type III interferon responses.

Tyrosine kinase 2 (TYK2) associates with interferon (IFN) alpha receptor, IL-10 receptor (IL-10R) beta and other cytokine receptor subunits for signal transduction, in response to various cytokines, including type-I and type-III IFNs, IL-6, IL-10, IL-12 and IL-23. Data on TYK2 dependence on cytokine responses and in vivo consequences of TYK2 deficiency are inconsistent. We investigated a TYK2 deficient patient, presenting with eczema, skin abscesses, respiratory infections and IgE levels >1000 U/mL, without viral or mycobacterial infections and a corresponding cellular model to analyze the role of TYK2 in type-III IFN mediated responses and NK-cell function. We established a novel simple diagnostic monocyte assay to show that the mutation completely abolishes the IFN-α mediated antiviral response. It also partly reduces IL-10 but not IL-6 mediated signaling associated with reduced IL-10Rß expression. However, we found almost normal type-III IFN signaling associated with minimal impairment of virus control in a TYK2 deficient human cell line. Contrary to observations in TYK2 deficient mice, NK-cell phenotype and function, including IL-12/IL-18 mediated responses, were normal in the patient. Thus, preserved type-III IFN responses and normal NK-cell function may contribute to antiviral protection in TYK2 deficiency leading to a surprisingly mild human phenotype.

TTC25 Deficiency Results in Defects of the Outer Dynein Arm Docking Machinery and Primary Ciliary Dyskinesia with Left-Right Body Asymmetry Randomization.

Multiprotein complexes referred to as outer dynein arms (ODAs) develop the main mechanical force to generate the ciliary and flagellar beat. ODA defects are the most common cause of primary ciliary dyskinesia (PCD), a congenital disorder of ciliary beating, characterized by recurrent infections of the upper and lower airways, as well as by progressive lung failure and randomization of left-right body asymmetry. Using a whole-exome sequencing approach, we identified recessive loss-of-function mutations within TTC25 in three individuals from two unrelated families affected by PCD. Mice generated by CRISPR/Cas9 technology and carrying a deletion of exons 2 and 3 in Ttc25 presented with laterality defects. Consistently, we observed immotile nodal cilia and missing leftward flow via particle image velocimetry. Furthermore, transmission electron microscopy (TEM) analysis in TTC25-deficient mice revealed an absence of ODAs. Consistent with our findings in mice, we were able to show loss of the ciliary ODAs in humans via TEM and immunofluorescence (IF) analyses. Additionally, IF analyses revealed an absence of the ODA docking complex (ODA-DC), along with its known components CCDC114, CCDC151, and ARMC4. Co-immunoprecipitation revealed interaction between the ODA-DC component CCDC114 and TTC25. Thus, here we report TTC25 as a new member of the ODA-DC machinery in humans and mice.

Deficient methylation and formylation of mt-tRNA(Met) wobble cytosine in a patient carrying mutations in NSUN3.

Epitranscriptome modifications are required for structure and function of RNA and defects in these pathways have been associated with human disease. Here we identify the RNA target for the previously uncharacterized 5-methylcytosine (m(5)C) methyltransferase NSun3 and link m(5)C RNA modifications with energy metabolism. Using whole-exome sequencing, we identified loss-of-function mutations in NSUN3 in a patient presenting with combined mitochondrial respiratory chain complex deficiency. Patient-derived fibroblasts exhibit severe defects in mitochondrial translation that can be rescued by exogenous expression of NSun3. We show that NSun3 is required for deposition of m(5)C at the anticodon loop in the mitochondrially encoded transfer RNA methionine (mt-tRNA(Met)). Further, we demonstrate that m(5)C deficiency in mt-tRNA(Met) results in the lack of 5-formylcytosine (f(5)C) at the same tRNA position. Our findings demonstrate that NSUN3 is necessary for efficient mitochondrial translation and reveal that f(5)C in human mitochondrial RNA is generated by oxidative processing of m(5)C.

Whole-exome sequencing in an isolated population from the Dalmatian island of Vis.

We have whole-exome sequenced 176 individuals from the isolated population of the island of Vis in Croatia in order to describe exonic variation architecture. We found 290 577 single nucleotide variants (SNVs), 65% of which are singletons, low frequency or rare variants. A total of 25 430 (9%) SNVs are novel, previously not catalogued in NHLBI GO Exome Sequencing Project, UK10K-Generation Scotland, 1000Genomes Project, ExAC or NCBI Reference Assembly dbSNP. The majority of these variants (76%) are singletons. Comparable to data obtained from UK10K-Generation Scotland that were sequenced and analysed using the same protocols, we detected an enrichment of potentially damaging variants (non-synonymous and loss-of-function) in the low frequency and common variant categories. On average 115 (range 93-140) genotypes with loss-of-function variants, 23 (15-34) of which were homozygous, were identified per person. The landscape of loss-of-function variants across an exome revealed that variants mainly accumulated in genes on the xenobiotic-related pathways, of which majority coded for enzymes. The frequency of loss-of-function variants was additionally increased in Vis runs of homozygosity regions where variants mainly affected signalling pathways. This work confirms the isolate status of Vis population by means of whole-exome sequence and reveals the pattern of loss-of-function mutations, which resembles the trails of adaptive evolution that were found in other species. By cataloguing the exomic variants and describing the allelic structure of the Vis population, this study will serve as a valuable resource for future genetic studies of human diseases, population genetics and evolution in this population.

DNAH11 Localization in the Proximal Region of Respiratory Cilia Defines Distinct Outer Dynein Arm Complexes.

Primary ciliary dyskinesia (PCD) is a recessively inherited disease that leads to chronic respiratory disorders owing to impaired mucociliary clearance. Conventional transmission electron microscopy (TEM) is a diagnostic standard to identify ultrastructural defects in respiratory cilia but is not useful in approximately 30% of PCD cases, which have normal ciliary ultrastructure. DNAH11 mutations are a common cause of PCD with normal ciliary ultrastructure and hyperkinetic ciliary beating, but its pathophysiology remains poorly understood. We therefore characterized DNAH11 in human respiratory cilia by immunofluorescence microscopy (IFM) in the context of PCD. We used whole-exome and targeted next-generation sequence analysis as well as Sanger sequencing to identify and confirm eight novel loss-of-function DNAH11 mutations. We designed and validated a monoclonal antibody specific to DNAH11 and performed high-resolution IFM of both control and PCD-affected human respiratory cells, as well as samples from green fluorescent protein (GFP)-left-right dynein mice, to determine the ciliary localization of DNAH11. IFM analysis demonstrated native DNAH11 localization in only the proximal region of wild-type human respiratory cilia and loss of DNAH11 in individuals with PCD with certain loss-of-function DNAH11 mutations. GFP-left-right dynein mice confirmed proximal DNAH11 localization in tracheal cilia. DNAH11 retained proximal localization in respiratory cilia of individuals with PCD with distinct ultrastructural defects, such as the absence of outer dynein arms (ODAs). TEM tomography detected a partial reduction of ODAs in DNAH11-deficient cilia. DNAH11 mutations result in a subtle ODA defect in only the proximal region of respiratory cilia, which is detectable by IFM and TEM tomography.

Homozygous loss-of-function variants in European cosmopolitan and isolate populations.

Homozygous loss of function (HLOF) variants provide a valuable window on gene function in humans, as well as an inventory of the human genes that are not essential for survival and reproduction. All humans carry at least a few HLOF variants, but the exact number of inactivated genes that can be tolerated is currently unknown­as are the phenotypic effects of losing function for most human genes. Here, we make use of 1432 whole exome sequences from five European populations to expand the catalogue of known human HLOF mutations; after stringent filtering of variants in our dataset, we identify a total of 173 HLOF mutations, 76 (44%) of which have not been observed previously. We find that population isolates are particularly well suited to surveys of novel HLOF genes because individuals in such populations carry extensive runs of homozygosity, which we show are enriched for novel, rare HLOF variants. Further, we make use of extensive phenotypic data to show that most HLOFs, ascertained in population-based samples, appear to have little detectable effect on the phenotype. On the contrary, we document several genes directly implicated in disease that seem to tolerate HLOF variants. Overall HLOF genes are enriched for olfactory receptor function and are expressed in testes more often than expected, consistent with reduced purifying selection and incipient pseudogenisation.

Deficiency of ECHS1 causes mitochondrial encephalopathy with cardiac involvement.

OBJECTIVE: Short-chain enoyl-CoA hydratase (ECHS1) is a multifunctional mitochondrial matrix enzyme that is involved in the oxidation of fatty acids and essential amino acids such as valine. Here, we describe the broad phenotypic spectrum and pathobiochemistry of individuals with autosomal-recessive ECHS1 deficiency. METHODS: Using exome sequencing, we identified ten unrelated individuals carrying compound heterozygous or homozygous mutations in ECHS1. Functional investigations in patient-derived fibroblast cell lines included immunoblotting, enzyme activity measurement, and a palmitate loading assay. RESULTS: Patients showed a heterogeneous phenotype with disease onset in the first year of life and course ranging from neonatal death to survival into adulthood. The most prominent clinical features were encephalopathy (10/10), deafness (9/9), epilepsy (6/9), optic atrophy (6/10), and cardiomyopathy (4/10). Serum lactate was elevated and brain magnetic resonance imaging showed white matter changes or a Leigh-like pattern resembling disorders of mitochondrial energy metabolism. Analysis of patients' fibroblast cell lines (6/10) provided further evidence for the pathogenicity of the respective mutations by showing reduced ECHS1 protein levels and reduced 2-enoyl-CoA hydratase activity. While serum acylcarnitine profiles were largely normal, in vitro palmitate loading of patient fibroblasts revealed increased butyrylcarnitine, unmasking the functional defect in mitochondrial ß-oxidation of short-chain fatty acids. Urinary excretion of 2-methyl-2,3-dihydroxybutyrate - a potential derivative of acryloyl-CoA in the valine catabolic pathway - was significantly increased, indicating impaired valine oxidation. INTERPRETATION: In conclusion, we define the phenotypic spectrum of a new syndrome caused by ECHS1 deficiency. We speculate that both the ß-oxidation defect and the block in l-valine metabolism, with accumulation of toxic methacrylyl-CoA and acryloyl-CoA, contribute to the disorder that may be amenable to metabolic treatment approaches.

Immunofluorescence Analysis and Diagnosis of Primary Ciliary Dyskinesia with Radial Spoke Defects.

Primary ciliary dyskinesia (PCD) is a genetically heterogeneous recessive disorder caused by several distinct defects in genes responsible for ciliary beating, leading to defective mucociliary clearance often associated with randomization of left/right body asymmetry. Individuals with PCD caused by defective radial spoke (RS) heads are difficult to diagnose owing to lack of gross ultrastructural defects and absence of situs inversus. Thus far, most mutations identified in human radial spoke genes (RSPH) are loss-of-function mutations, and missense variants have been rarely described. We studied the consequences of different RSPH9, RSPH4A, and RSPH1 mutations on the assembly of the RS complex to improve diagnostics in PCD. We report 21 individuals with PCD (16 families) with biallelic mutations in RSPH9, RSPH4A, and RSPH1, including seven novel mutations comprising missense variants, and performed high-resolution immunofluorescence analysis of human respiratory cilia. Missense variants are frequent genetic defects in PCD with RS defects. Absence of RSPH4A due to mutations in RSPH4A results in deficient axonemal assembly of the RS head components RSPH1 and RSPH9. RSPH1 mutant cilia, lacking RSPH1, fail to assemble RSPH9, whereas RSPH9 mutations result in axonemal absence of RSPH9, but do not affect the assembly of the other head proteins, RSPH1 and RSPH4A. Interestingly, our results were identical in individuals carrying loss-of-function mutations, missense variants, or one amino acid deletion. Immunofluorescence analysis can improve diagnosis of PCD in patients with loss-of-function mutations as well as missense variants. RSPH4A is the core protein of the RS head.

Insights into hominid evolution from the gorilla genome sequence.

Gorillas are humans' closest living relatives after chimpanzees, and are of comparable importance for the study of human origins and evolution. Here we present the assembly and analysis of a genome sequence for the western lowland gorilla, and compare the whole genomes of all extant great ape genera. We propose a synthesis of genetic and fossil evidence consistent with placing the human-chimpanzee and human-chimpanzee-gorilla speciation events at approximately 6 and 10 million years ago. In 30% of the genome, gorilla is closer to human or chimpanzee than the latter are to each other; this is rarer around coding genes, indicating pervasive selection throughout great ape evolution, and has functional consequences in gene expression. A comparison of protein coding genes reveals approximately 500 genes showing accelerated evolution on each of the gorilla, human and chimpanzee lineages, and evidence for parallel acceleration, particularly of genes involved in hearing. We also compare the western and eastern gorilla species, estimating an average sequence divergence time 1.75 million years ago, but with evidence for more recent genetic exchange and a population bottleneck in the eastern species. The use of the genome sequence in these and future analyses will promote a deeper understanding of great ape biology and evolution.

A Runx1-Smad6 rheostat controls Runx1 activity during embryonic hematopoiesis.

The oncogenic transcription factor Runx1 is required for the specification of definitive hematopoietic stem cells (HSC) in the developing embryo. The activity of this master regulator is tightly controlled during development. The transcription factors that upregulate the expression of Runx1 also upregulate the expression of Smad6, the inhibitory Smad, which controls Runx1 activity by targeting it to the proteasome. Here we show that Runx1, in conjunction with Fli1, Gata2, and Scl, directly regulates the expression of Smad6 in the aorta-gonad-mesonephros (AGM) region in the developing embryo, where HSCs originate. Runx1 regulates Smad6 activity via a novel upstream enhancer, and Runx1 null embryos show reduced Smad6 transcripts in the yolk-sac and c-Kit-positive fetal liver cells. By directly regulating the expression of Smad6, Runx1 sets up a functional rheostat to control its own activity. The perturbation of this rheostat, using a proteasomal inhibitor, results in an increase in Runx1 and Smad6 levels that can be directly attributed to increased Runx1 binding to tissue-specific regulatory elements of these genes. Taken together, we describe a scenario in which a key hematopoietic transcription factor controls its own expression levels by transcriptionally controlling its controller.

Endoglin expression in blood and endothelium is differentially regulated by modular assembly of the Ets/Gata hemangioblast code.

Endoglin is an accessory receptor for TGF-beta signaling and is required for normal hemangioblast, early hematopoietic, and vascular development. We have previously shown that an upstream enhancer, Eng -8, together with the promoter region, mediates robust endothelial expression yet is inactive in blood. To identify hematopoietic regulatory elements, we used array-based methods to determine chromatin accessibility across the entire locus. Subsequent transgenic analysis of candidate elements showed that an endothelial enhancer at Eng +9 when combined with an element at Eng +7 functions as a strong hemato-endothelial enhancer. Chromatin immunoprecipitation (ChIP)-chip analysis demonstrated specific binding of Ets factors to the promoter as well as to the -8, +7+9 enhancers in both blood and endothelial cells. By contrast Pu.1, an Ets factor specific to the blood lineage, and Gata2 binding was only detected in blood. Gata2 was bound only at +7 and GATA motifs were required for hematopoietic activity. This modular assembly of regulators gives blood and endothelial cells the regulatory freedom to independently fine-tune gene expression and emphasizes the role of regulatory divergence in driving functional divergence.

Application of phage display to high throughput antibody generation and characterization.

We have created a high quality phage display library containing over 1010 human antibodies and describe its use in the generation of antibodies on an unprecedented scale. We have selected, screened and sequenced over 38,000 recombinant antibodies to 292 antigens, yielding over 7,200 unique clones. 4,400 antibodies were characterized by specificity testing and detailed sequence analysis and the data/clones are available online. Sensitive detection was demonstrated in a bead based flow cytometry assay. Furthermore, positive staining by immunohistochemistry on tissue microarrays was found for 37% (143/381) of antibodies. Thus, we have demonstrated the potential of and illuminated the issues associated with genome-wide monoclonal antibody generation.

Gata2, Fli1, and Scl form a recursively wired gene-regulatory circuit during early hematopoietic development.

Conservation of the vertebrate body plan has been attributed to the evolutionary stability of gene-regulatory networks (GRNs). We describe a regulatory circuit made up of Gata2, Fli1, and Scl/Tal1 and their enhancers, Gata2-3, Fli1+12, and Scl+19, that operates during specification of hematopoiesis in the mouse embryo. We show that the Fli1+12 enhancer, like the Gata2-3 and Scl+19 enhancers, targets hematopoietic stem cells (HSCs) and relies on a combination of Ets, Gata, and E-Box motifs. We show that the Gata2-3 enhancer also uses a similar cluster of motifs and that Gata2, Fli1, and Scl are expressed in embryonic day-11.5 dorsal aorta where HSCs originate and in fetal liver where they multiply. The three HSC enhancers in these tissues and in ES cell-derived hemangioblast equivalents are bound by each of these transcription factors (TFs) and form a fully connected triad that constitutes a previously undescribed example of both this network motif in mammalian development and a GRN kernel operating during the specification of a mammalian stem cell.

The systematic functional characterisation of Xq28 genes prioritises candidate disease genes.

BACKGROUND: Well known for its gene density and the large number of mapped diseases, the human sub-chromosomal region Xq28 has long been a focus of genome research. Over 40 of approximately 300 X-linked diseases map to this region, and systematic mapping, transcript identification, and mutation analysis has led to the identification of causative genes for 26 of these diseases, leaving another 17 diseases mapped to Xq28, where the causative gene is still unknown. To expedite disease gene identification, we have initiated the functional characterisation of all known Xq28 genes. RESULTS: By using a systematic approach, we describe the Xq28 genes by RNA in situ hybridisation and Northern blotting of the mouse orthologs, as well as subcellular localisation and data mining of the human genes. We have developed a relational web-accessible database with comprehensive query options integrating all experimental data. Using this database, we matched gene expression patterns with affected tissues for 16 of the 17 remaining Xq28 linked diseases, where the causative gene is unknown. CONCLUSION: By using this systematic approach, we have prioritised genes in linkage regions of Xq28-mapped diseases to an amenable number for mutational screens. Our database can be queried by any researcher performing highly specified searches including diseases not listed in OMIM or diseases that might be linked to Xq28 in the future.

Alternative pre-mRNA processing regulates cell-type specific expression of the IL4l1 and NUP62 genes.

BACKGROUND: Given the complexity of higher organisms, the number of genes encoded by their genomes is surprisingly small. Tissue specific regulation of expression and splicing are major factors enhancing the number of the encoded products. Commonly these mechanisms are intragenic and affect only one gene. RESULTS: Here we provide evidence that the IL4I1 gene is specifically transcribed from the apparent promoter of the upstream NUP62 gene, and that the first two exons of NUP62 are also contained in the novel IL4I1_2 variant. While expression of IL4I1 driven from its previously described promoter is found mostly in B cells, the expression driven by the NUP62 promoter is restricted to cells in testis (Sertoli cells) and in the brain (e.g., Purkinje cells). Since NUP62 is itself ubiquitously expressed, the IL4I1_2 variant likely derives from cell type specific alternative pre-mRNA processing. CONCLUSION: Comparative genomics suggest that the promoter upstream of the NUP62 gene originally belonged to the IL4I1 gene and was later acquired by NUP62 via insertion of a retroposon. Since both genes are apparently essential, the promoter had to serve two genes afterwards. Expression of the IL4I1 gene from the "NUP62" promoter and the tissue specific involvement of the pre-mRNA processing machinery to regulate expression of two unrelated proteins indicate a novel mechanism of gene regulation.