Submegabase copy number variations arise during cerebral cortical neurogenesis as revealed by single-cell whole-genome sequencing.

Somatic copy number variations (CNVs) exist in the brain, but their genesis, prevalence, forms, and biological impact remain unclear, even within experimentally tractable animal models. We combined a transposase-based amplification (TbA) methodology for single-cell whole-genome sequencing with a bioinformatic approach for filtering unreliable CNVs (FUnC), developed from machine learning trained on lymphocyte V(D)J recombination. TbA-FUnC offered superior genomic coverage and removed >90% of false-positive CNV calls, allowing extensive examination of submegabase CNVs from over 500 cells throughout the neurogenic period of cerebral cortical development in Mus musculus Thousands of previously undocumented CNVs were identified. Half were less than 1 Mb in size, with deletions 4× more common than amplification events, and were randomly distributed throughout the genome. However, CNV prevalence during embryonic cortical development was nonrandom, peaking at midneurogenesis with levels triple those found at younger ages before falling to intermediate quantities. These data identify pervasive small and large CNVs as early contributors to neural genomic mosaicism, producing genomically diverse cellular building blocks that form the highly organized, mature brain.

Characterizing transcriptional heterogeneity through pathway and gene set overdispersion analysis.

The transcriptional state of a cell reflects a variety of biological factors, from cell-type-specific features to transient processes such as the cell cycle, all of which may be of interest. However, identifying such aspects from noisy single-cell RNA-seq data remains challenging. We developed pathway and gene set overdispersion analysis (PAGODA) to resolve multiple, potentially overlapping aspects of transcriptional heterogeneity by testing gene sets for coordinated variability among measured cells.

Altered levels of mitochondrial DNA are associated with female age, aneuploidy, and provide an independent measure of embryonic implantation potential.

Mitochondria play a vital role in embryo development. They are the principal site of energy production and have various other critical cellular functions. Despite the importance of this organelle, little is known about the extent of variation in mitochondrial DNA (mtDNA) between individual human embryos prior to implantation. This study investigated the biological and clinical relevance of the quantity of mtDNA in 379 embryos. These were examined via a combination of microarray comparative genomic hybridisation (aCGH), quantitative PCR and next generation sequencing (NGS), providing information on chromosomal status, amount of mtDNA, and presence of mutations in the mitochondrial genome. The quantity of mtDNA was significantly higher in embryos from older women (P=0.003). Additionally, mtDNA levels were elevated in aneuploid embryos, independent of age (P=0.025). Assessment of clinical outcomes after transfer of euploid embryos to the uterus revealed that blastocysts that successfully implanted tended to contain lower mtDNA quantities than those failing to implant (P=0.007). Importantly, an mtDNA quantity threshold was established, above which implantation was never observed. Subsequently, the predictive value of this threshold was confirmed in an independent blinded prospective study, indicating that abnormal mtDNA levels are present in 30% of non-implanting euploid embryos, but are not seen in embryos forming a viable pregnancy. NGS did not reveal any increase in mutation in blastocysts with elevated mtDNA levels. The results of this study suggest that increased mtDNA may be related to elevated metabolism and are associated with reduced viability, a possibility consistent with the 'quiet embryo' hypothesis. Importantly, the findings suggest a potential role for mitochondria in female reproductive aging and the genesis of aneuploidy. Of clinical significance, we propose that mtDNA content represents a novel biomarker with potential value for in vitro fertilisation (IVF) treatment, revealing chromosomally normal blastocysts incapable of producing a viable pregnancy.

Whole-genome haplotyping by dilution, amplification, and sequencing.

Standard whole-genome genotyping technologies are unable to determine haplotypes. Here we describe a method for rapid and cost-effective long-range haplotyping. Genomic DNA is diluted and distributed into multiple aliquots such that each aliquot receives a fraction of a haploid copy. The DNA template in each aliquot is amplified by multiple displacement amplification, converted into barcoded sequencing libraries using Nextera technology, and sequenced in multiplexed pools. To assess the performance of our method, we combined two male genomic DNA samples at equal ratios, resulting in a sample with diploid X chromosomes with known haplotypes. Pools of the multiplexed sequencing libraries were subjected to targeted pull-down of a 1-Mb contiguous region of the X-chromosome Duchenne muscular dystrophy gene. We were able to phase the Duchenne muscular dystrophy region into two contiguous haplotype blocks with a mean length of 494 kb. The haplotypes showed 99% agreement with the consensus base calls made by sequencing the individual DNAs. We subsequently used the strategy to haplotype two human genomes. Standard genomic sequencing to identify all heterozygous SNPs in the sample was combined with dilution-amplification-based sequencing data to resolve the phase of identified heterozygous SNPs. Using this procedure, we were able to phase >95% of the heterozygous SNPs from the diploid sequence data. The N50 for a Yoruba male DNA was 702 kb whereas the N50 for a European female DNA was 358 kb. Therefore, the strategy described here is suitable for haplotyping of a set of targeted regions as well as of the entire genome.

Quantitative phenotyping via deep barcode sequencing.

Next-generation DNA sequencing technologies have revolutionized diverse genomics applications, including de novo genome sequencing, SNP detection, chromatin immunoprecipitation, and transcriptome analysis. Here we apply deep sequencing to genome-scale fitness profiling to evaluate yeast strain collections in parallel. This method, Barcode analysis by Sequencing, or "Bar-seq," outperforms the current benchmark barcode microarray assay in terms of both dynamic range and throughput. When applied to a complex chemogenomic assay, Bar-seq quantitatively identifies drug targets, with performance superior to the benchmark microarray assay. We also show that Bar-seq is well-suited for a multiplex format. We completely re-sequenced and re-annotated the yeast deletion collection using deep sequencing, found that approximately 20% of the barcodes and common priming sequences varied from expectation, and used this revised list of barcode sequences to improve data quality. Together, this new assay and analysis routine provide a deep-sequencing-based toolkit for identifying gene-environment interactions on a genome-wide scale.

Yeast Barcoders: a chemogenomic application of a universal donor-strain collection carrying bar-code identifiers.

The ability to perform complex bioassays in parallel enables experiments that are otherwise impossible because of throughput and cost constraints. For example, highly parallel chemical-genetic screens using pooled collections of thousands of defined Saccharomyces cerevisiae gene deletion strains are feasible because each strain is bar-coded with unique DNA sequences. It is, however, time-consuming and expensive to individually bar-code individual strains. To provide a simple and general method of barcoding yeast collections, we built a set of donor strains, called Barcoders, with unique bar codes that can be systematically transferred to any S. cerevisiae collection. We applied this technology by generating a collection of bar-coded 'decreased abundance by mRNA perturbation' (DAmP) loss-of-function strains comprising 87.1% of all essential yeast genes. These experiments validate both the Barcoders and the DAmP strain collection as useful tools for genome-wide chemical-genetic assays.

Circulating transcripts in maternal blood reflect a molecular signature of early-onset preeclampsia.

Circulating RNA (C-RNA) is continually released into the bloodstream from tissues throughout the body, offering an opportunity to noninvasively monitor all aspects of pregnancy health from conception to birth. We asked whether C-RNA analysis could robustly detect aberrations in patients diagnosed with preeclampsia (PE), a prevalent and potentially fatal pregnancy complication. As an initial examination, we sequenced the circulating transcriptome from 40 pregnancies at the time of severe, early-onset PE diagnosis and 73 gestational age-matched controls. Differential expression analysis identified 30 transcripts with gene ontology annotations and tissue expression patterns consistent with the placental dysfunction, impaired fetal development, and maternal immune and cardiovascular system dysregulation characteristic of PE. Furthermore, machine learning identified combinations of 49 C-RNA transcripts that classified an independent cohort of patients (early-onset PE, n = 12; control, n = 12) with 85 to 89% accuracy. C-RNA may thus hold promise for improving the diagnosis and identification of at-risk pregnancies.

Mutations in the PI3K/PTEN/TSC2 pathway contribute to mammalian target of rapamycin activity and increased translation under hypoxic conditions.

Decreased oxygen causes a rapid inhibition of mRNA translation. An important regulatory mechanism of translational repression under hypoxic conditions involves inhibition of the mammalian target of rapamycin (mTOR). mTOR is a target of the phosphatase and tensin homologue detected on chromosome 10 (PTEN)/phosphatidylinositol 3-kinase/AKT/TSC2 pathway, a pathway that is frequently mutated in human cancers. Although hypoxia has been shown to inhibit mTOR activity, we show here that the hypoxia-induced inhibition of mTOR activity is attenuated in cells lacking TSC2 or PTEN, resulting in a higher translation rate even under hypoxic conditions. Comparison of mTOR inhibition by hypoxia alone or in combination with rapamycin showed that prolonged exposure to hypoxia was required to fully inhibit mTOR activity even in wild-type cells. Increased mTOR activity and protein synthesis did not translate into enhanced cell proliferation rates. However, lack of TSC2 resulted in a survival advantage when cells were exposed to hypoxia. Protection against hypoxia-induced cell death due to TSC2 deficiency is rapamycin-resistant, suggesting that TSC2 affects an apoptotic pathway. Tumors derived from TSC2 wild-type cells exhibited a growth delay compared with TSC2-deficient tumors, indicating that enhanced mTOR activity is advantageous in the initial phase of tumor growth. Therefore, failure to inhibit mTOR under oxygen-limiting conditions can be affected by upstream activating mutations and increases the survival and growth of hypoxic tumor cells.

Connective tissue growth factor-specific monoclonal antibody therapy inhibits pancreatic tumor growth and metastasis.

Pancreatic cancer is highly aggressive and refractory to most existing therapies. Past studies have shown that connective tissue growth factor (CTGF) expression is elevated in human pancreatic adenocarcinomas and some pancreatic cancer cell lines. To address whether and how CTGF influences tumor growth, we generated pancreatic tumor cell lines that overexpress different levels of human CTGF. The effect of CTGF overexpression on cell proliferation was measured in vitro in monolayer culture, suspension culture, or soft agar, and in vivo in tumor xenografts. Although there was no effect of CTGF expression on proliferation in two-dimensional cultures, anchorage-independent growth (AIG) was enhanced. The capacity of CTGF to enhance AIG in vitro was linked to enhanced pancreatic tumor growth in vivo when these cells were implanted s.c. in nude mice. Administration of a neutralizing CTGF-specific monoclonal antibody, FG-3019, had no effect on monolayer cell proliferation, but blocked AIG in soft agar. Consistent with this observation, anti-CTGF treatment of mice bearing established CTGF-expressing tumors abrogated CTGF-dependent tumor growth and inhibited lymph node metastases without any toxicity observed in normal tissue. Together, these studies implicate CTGF as a new target in pancreatic cancer and suggest that inhibition of CTGF with a human monoclonal antibody may control primary and metastatic tumor growth.

Regulation of hypoxia-inducible factor 1alpha expression and function by the mammalian target of rapamycin.

Hypoxia-inducible factor 1 (HIF-1) is a heterodimeric transcription factor containing an inducibly expressed HIF-1alpha subunit and a constititutively expressed HIF-1beta subunit. Under hypoxic conditions, the HIF-1alpha subunit accumulates due to a decrease in the rate of proteolytic degradation, and the resulting HIF-1alpha-HIF-1beta heterodimers undergo posttranslational modifications that promote transactivation. Recent studies suggest that amplified signaling through phosphoinositide 3-kinase, and its downstream target, mTOR, enhances HIF-1-dependent gene expression in certain cell types. In the present study, we have explored further the linkage between mTOR and HIF-1 in PC-3 prostate cancer cells treated with hypoxia or the hypoxia mimetic agent, CoCl(2). Pretreatment of PC-3 cells with the mTOR inhibitor, rapamycin, inhibited both the accumulation of HIF-1alpha and HIF-1-dependent transcription induced by hypoxia or CoCl(2). Transfection of these cells with wild-type mTOR enhanced HIF-1 activation by hypoxia or CoCl(2), while expression of a rapamycin-resistant mTOR mutant rendered both HIF-1alpha stabilization and HIF-1 transactivating function refractory to inhibition by rapamycin. Studies with GAL4-HIF-1alpha fusion proteins pinpointed the oxygen-dependent degradation domain as a critical target for the rapamycin-sensitive, mTOR-dependent signaling pathway leading to HIF-1alpha stabilization by CoCl(2). These studies position mTOR as an upstream activator of HIF-1 function in cancer cells and suggest that the antitumor activity of rapamycin is mediated, in part, through the inhibition of cellular responses to hypoxic stress.

Anoxia is necessary for tumor cell toxicity caused by a low-oxygen environment.

Cells exposed to oxygen deprivation in vitro have been shown to reduce proliferation and/or engage in programmed cell death. There is considerable controversy in the literature as to the role of hypoxia-inducible factor-1 (HIF-1) and HIF-1 target genes in initiating these responses. We therefore examined the oxygen dependence and the role of the hypoxia-responsive transcription factor HIF-1 in making the cellular death decision. Oxygen concentrations as low as 0.5% did not alter the growth of HIF-1-proficient or HIF-1-deficient murine fibroblasts, or human tumor cells, despite the appropriate induction of HIF-1 target genes. Severe hypoxia (<0.01% oxygen) did induced apoptosis, resulting in decreased colony formation, chromatin condensation, DNA fragmentation, and caspase activation but also independent of HIF1alpha status. Transcriptional induction of HIF-1-dependent genes putatively involved in cell death like BNip3 and BNip3L was therefore disassociated from hypoxia-dependent toxicity. Likewise, forced overexpression of a nondegradable form of HIF-1alpha in several human tumor cell lines was not sufficient to induce apoptosis under normoxic conditions. Taken together, these findings indicate that additional molecular events are triggered by anoxia in a HIF-1-independent manner, and these changes are necessary for cell death observed in low-oxygen environments.

Neuronal subtypes and diversity revealed by single-nucleus RNA sequencing of the human brain.

The human brain has enormously complex cellular diversity and connectivities fundamental to our neural functions, yet difficulties in interrogating individual neurons has impeded understanding of the underlying transcriptional landscape. We developed a scalable approach to sequence and quantify RNA molecules in isolated neuronal nuclei from a postmortem brain, generating 3227 sets of single-neuron data from six distinct regions of the cerebral cortex. Using an iterative clustering and classification approach, we identified 16 neuronal subtypes that were further annotated on the basis of known markers and cortical cytoarchitecture. These data demonstrate a robust and scalable method for identifying and categorizing single nuclear transcriptomes, revealing shared genes sufficient to distinguish previously unknown and orthologous neuronal subtypes as well as regional identity and transcriptomic heterogeneity within the human brain.

Effects of perillyl alcohol on melanoma in the TPras mouse model.

This study evaluates the chemopreventive effects of topically applied perillyl alcohol on the development of melanoma in TPras transgenic mice. Our strategy was to target critical pathways in the development of melanoma, in particular, the ras pathway. Ras has been shown in our experimental mouse model, as well as others, to be important in the development and maintenance of melanomas. Perillyl alcohol (POH), a naturally occurring monoterpene, inhibits the isoprenylation of small G protein, including Ras. POH (10 mM) was applied to the shaved dorsal skin of TPras mice starting 1 week before five treatments of dimethylbenz[a]anthracene (50 microg) and was continued for 38 weeks. We observed a delay in the appearance of tumors and a 25-35% reduction in melanoma incidence. POH treatment of melanoma cells in vitro reduced the levels of detectable Ras protein and inhibits the activation of downstream targets, mitogen-activated protein kinases and Akt. POH only minimally induced apoptosis in this system. Pretreatment but not post-treatment of the melanoma cells with POH, however, markedly reduced levels of UV-induced reactive oxygen species. These studies suggest that POH inhibition of the Ras signaling pathway may be an effective target for chemoprevention of melanoma.

Noninvasive identification and monitoring of cancer mutations by targeted deep sequencing of plasma DNA.

Plasma of cancer patients contains cell-free tumor DNA that carries information on tumor mutations and tumor burden. Individual mutations have been probed using allele-specific assays, but sequencing of entire genes to detect cancer mutations in circulating DNA has not been demonstrated. We developed a method for tagged-amplicon deep sequencing (TAm-Seq) and screened 5995 genomic bases for low-frequency mutations. Using this method, we identified cancer mutations present in circulating DNA at allele frequencies as low as 2%, with sensitivity and specificity of >97%. We identified mutations throughout the tumor suppressor gene TP53 in circulating DNA from 46 plasma samples of advanced ovarian cancer patients. We demonstrated use of TAm-Seq to noninvasively identify the origin of metastatic relapse in a patient with multiple primary tumors. In another case, we identified in plasma an EGFR mutation not found in an initial ovarian biopsy. We further used TAm-Seq to monitor tumor dynamics, and tracked 10 concomitant mutations in plasma of a metastatic breast cancer patient over 16 months. This low-cost, high-throughput method could facilitate analysis of circulating DNA as a noninvasive "liquid biopsy" for personalized cancer genomics.

Highly parallel genome-wide expression analysis of single mammalian cells.

BACKGROUND: We have developed a high-throughput amplification method for generating robust gene expression profiles using single cell or low RNA inputs. METHODOLOGY/PRINCIPAL FINDINGS: The method uses tagged priming and template-switching, resulting in the incorporation of universal PCR priming sites at both ends of the synthesized cDNA for global PCR amplification. Coupled with a whole-genome gene expression microarray platform, we routinely obtain expression correlation values of R(2)~0.76-0.80 between individual cells and R(2)~0.69 between 50 pg total RNA replicates. Expression profiles generated from single cells or 50 pg total RNA correlate well with that generated with higher input (1 ng total RNA) (R(2)~0.80). Also, the assay is sufficiently sensitive to detect, in a single cell, approximately 63% of the number of genes detected with 1 ng input, with approximately 97% of the genes detected in the single-cell input also detected in the higher input. CONCLUSIONS/SIGNIFICANCE: In summary, our method facilitates whole-genome gene expression profiling in contexts where starting material is extremely limiting, particularly in areas such as the study of progenitor cells in early development and tumor stem cell biology.

Hop, an active Mutator-like element in the genome of the fungus Fusarium oxysporum.

A new type of active DNA transposon has been identified in the genome of Fusarium oxysporum by its transposition into the niaD target gene. Two insertions within the final exon, in opposite orientations at the same nucleotide site, have been characterized. These elements, called Hop, are 3,299 bp long, with perfect terminal inverted repeats (TIRs) of 99 bp. The sequencing of genomic copies reveals a 9-bp target site duplication and no apparent sequence specificity at the insertion sites. The sequencing of a cDNA indicates that Hop does not contain an intron and encodes a putative transposase of 836 amino acids. The structural features (length, TIRs size, and 9-bp duplication), together with the presence of conserved domains in the transposase, strongly suggest that Hop is a Mutator-like element (MULE). Hop is thus the first active member of this family found beyond plants. The high rate of excision observed indicates that Hop is very active and thus represents a promising efficient tagging system for the isolation of fungal genes. The distribution of Hop elements within the Fusarium genus revealed that they are present in different species, suggesting that related elements could be present in other fungal genomes. In fact, Hop-related sequences have been identified in the survey of the entire genome sequence of three other ascomycetes, Magnaporthe grisea, Neurospora crassa, and Aspergillus fumigatus.