High throughput single-cell detection of multiplex CRISPR-edited gene modifications.

CRISPR-Cas9 gene editing has transformed our ability to rapidly interrogate the functional impact of somatic mutations in human cancers. Droplet-based technology enables the analysis of Cas9-introduced gene edits in thousands of single cells. Using this technology, we analyze Ba/F3 cells engineered to express single or multiplexed loss-of-function mutations recurrent in chronic lymphocytic leukemia. Our approach reliably quantifies mutational co-occurrences, zygosity status, and the occurrence of Cas9 edits at single-cell resolution.

High-throughput single-cell DNA sequencing of acute myeloid leukemia tumors with droplet microfluidics.

To enable the characterization of genetic heterogeneity in tumor cell populations, we developed a novel microfluidic approach that barcodes amplified genomic DNA from thousands of individual cancer cells confined to droplets. The barcodes are then used to reassemble the genetic profiles of cells from next-generation sequencing data. By using this approach, we sequenced longitudinally collected acute myeloid leukemia (AML) tumor populations from two patients and genotyped up to 62 disease relevant loci across more than 16,000 individual cells. Targeted single-cell sequencing was able to sensitively identify cells harboring pathogenic mutations during complete remission and uncovered complex clonal evolution within AML tumors that was not observable with bulk sequencing. We anticipate that this approach will make feasible the routine analysis of AML heterogeneity, leading to improved stratification and therapy selection for the disease.

Synchrotron micro-scale measurement of metal distributions in Phragmites australis and Typha latifolia root tissue from an urban brownfield site.

Liberty State Park in New Jersey, USA, is a "brownfield" site containing various levels of contaminants. To investigate metal uptake and distributions in plants on the brownfield site, Phragmites australis and Typha latifolia were collected in Liberty State Park during the growing season (May-September) in 2011 at two sites with the high and low metal loads, respectively. The objective of this study was to understand the metal (Fe, Mn, Cu, Pb and Zn) concentration and spatial distributions in P. australis and T. latifolia root systems with micro-meter scale resolution using synchrotron X-ray microfluorescence (µXRF) and synchrotron X-ray computed microtomography (µCMT) techniques. The root structure measurement by synchrotron µCMT showed that high X-ray attenuation substance appeared in the epidermis. Synchrotron µXRF measurement showed that metal concentrations and distributions in the root cross-section between epidermis and vascular tissue were statistically different. Significant correlations were found between metals (Cu, Mn, Pb and Zn) and Fe in the epidermis, implying that metals were scavenged by Fe oxides. The results from this study suggest that the expression of metal transport and accumulation within the root systems may be element specific. The information derived from this study can improve our current knowledge of the wetland plant ecological function in brownfield remediation.

Synchrotron study of metal localization in Typha latifolia L. root sections.

Understanding mechanisms that control plant root metal assimilation in soil is critical to the sustainable management of metal-contaminated land. With the assistance of the synchrotron X-ray fluorescence technique, this study investigated possible mechanisms that control the localization of Fe, Cu, Mn, Pb and Zn in the root tissues of Typha latifolia L. collected from a contaminated wetland. Metal localizations especially in the case of Fe and Pb in the dermal tissue and the vascular bundles were different. Cluster analysis was performed to divide the dermal tissue into iron-plaque-enriched dermal tissue and regular dermal tissue based on the spatial distribution of Pb and Fe. Factor analysis showed that Cu and Zn were closely correlated to each other in the dermal tissues. The association of Cu, Zn and Mn with Fe was strong in both regular dermal tissue and iron-plaque-enriched dermal tissue, while significant (p < 0.05) correlation of Fe with Pb was only observed in tissues enriched with iron plaque. In the vascular bundles, Zn, Mn and Cu showed strong association, suggesting that the localization of these three elements was controlled by a similar mechanism. Iron plaque in the peripheral dermal tissues acted as a barrier for Pb and a buffer for Zn, Cu and Mn. The Casparian strip regulated the transportation of metals from dermal tissues to the vascular bundles. The results suggested that the mechanisms controlling metal localization in root tissues varied with both tissue types and metals.

Synchrotron micro-scale study of trace metal transport and distribution in Spartina alterniflora root system in Yangtze River intertidal zone.

This study is focused on micro-scale measurement of metal (Ca, Cl, Fe, K, Mn, Cu, Pb, and Zn) distributions in Spartina alterniflora root system. The root samples were collected in the Yangtze River intertidal zone in July 2013. Synchrotron X-ray fluorescence (XRF), computed microtomography (CMT), and X-ray absorption near-edge structure (XANES) techniques, which provide micro-meter scale analytical resolution, were applied to this study. Although it was found that the metals of interest were distributed in both epidermis and vascular tissue with the varying concentrations, the results showed that Fe plaque was mainly distributed in the root epidermis. Other metals (e.g., Cu, Mn, Pb, and Zn) were correlated with Fe in the epidermis possibly due to scavenge by Fe plaque. Relatively high metal concentrations were observed in the root hair tip. This micro-scale investigation provides insights of understanding the metal uptake and spatial distribution as well as the function of Fe plaque governing metal transport in the root system.

DNA sequencing versus standard prenatal aneuploidy screening.

BACKGROUND: In high-risk pregnant women, noninvasive prenatal testing with the use of massively parallel sequencing of maternal plasma cell-free DNA (cfDNA testing) accurately detects fetal autosomal aneuploidy. Its performance in low-risk women is unclear. METHODS: At 21 centers in the United States, we collected blood samples from women with singleton pregnancies who were undergoing standard aneuploidy screening (serum biochemical assays with or without nuchal translucency measurement). We performed massively parallel sequencing in a blinded fashion to determine the chromosome dosage for each sample. The primary end point was a comparison of the false positive rates of detection of fetal trisomies 21 and 18 with the use of standard screening and cfDNA testing. Birth outcomes or karyotypes were the reference standard. RESULTS: The primary series included 1914 women (mean age, 29.6 years) with an eligible sample, a singleton fetus without aneuploidy, results from cfDNA testing, and a risk classification based on standard screening. For trisomies 21 and 18, the false positive rates with cfDNA testing were significantly lower than those with standard screening (0.3% vs. 3.6% for trisomy 21, P<0.001; and 0.2% vs. 0.6% for trisomy 18, P=0.03). The use of cfDNA testing detected all cases of aneuploidy (5 for trisomy 21, 2 for trisomy 18, and 1 for trisomy 13; negative predictive value, 100% [95% confidence interval, 99.8 to 100]). The positive predictive values for cfDNA testing versus standard screening were 45.5% versus 4.2% for trisomy 21 and 40.0% versus 8.3% for trisomy 18. CONCLUSIONS: In a general obstetrical population, prenatal testing with the use of cfDNA had significantly lower false positive rates and higher positive predictive values for detection of trisomies 21 and 18 than standard screening. (Funded by Illumina; ClinicalTrials.gov number, NCT01663350.).

The curious case of miRNAs in circulation: potential diagnostic biomarkers?

The pervasive occurrence of cell-free miRNAs in circulation suggests that these species play an emerging role as regulatory molecules in the secretory environment. Are these molecules released fortuitously with no clear biological intent? Or do they constitute a regulatory architecture that has evolved to modulate gene expression using the highways and byways of the circulatory system? The study of circulating miRNAs continues to increase our understanding of the regulation of genomes. The diversity of acellular miRNAs from a functional perspective is discussed, and in particular we explore their utility in a clinical setting as blood-based biomarkers for diseases.

Genome-wide maps of circulating miRNA biomarkers for ulcerative colitis.

Inflammatory Bowel Disease--comprised of Crohn's Disease and Ulcerative Colitis (UC)--is a complex, multi-factorial inflammatory disorder of the gastrointestinal tract. In this study we have explored the utility of naturally occurring circulating miRNAs as potential blood-based biomarkers for non-invasive prediction of UC incidences. Whole genome maps of circulating miRNAs in micro-vesicles, Peripheral Blood Mononuclear Cells and platelets have been constructed from a cohort of 20 UC patients and 20 normal individuals. Through Significance Analysis of Microarrays, a signature of 31 differentially expressed platelet-derived miRNAs has been identified and biomarker performance estimated through a non-probabilistic binary linear classification using Support Vector Machines. Through this approach, classifier measurements reveal a predictive score of 92.8% accuracy, 96.2% specificity and 89.5% sensitivity in distinguishing UC patients from normal individuals. Additionally, the platelet-derived biomarker signature can be validated at 88% accuracy through qPCR assays, and a majority of the miRNAs in this panel can be demonstrated to sub-stratify into 4 highly correlated intensity based clusters. Analysis of predicted targets of these biomarkers reveal an enrichment of pathways associated with cytoskeleton assembly, transport, membrane permeability and regulation of transcription factors engaged in a variety of regulatory cascades that are consistent with a cell-mediated immune response model of intestinal inflammation. Interestingly, comparison of the miRNA biomarker panel and genetic loci implicated in IBD through genome-wide association studies identifies a physical linkage between hsa-miR-941 and a UC susceptibility loci located on Chr 20. Taken together, analysis of these expression maps outlines a promising catalog of novel platelet-derived miRNA biomarkers of clinical utility and provides insight into the potential biological function of these candidates in disease pathogenesis.

Impact of cellular miRNAs on circulating miRNA biomarker signatures.

Effective diagnosis and surveillance of complex multi-factorial disorders such as cancer can be improved by screening of easily accessible biomarkers. Highly stable cell free Circulating Nucleic Acids (CNA) present as both RNA and DNA species have been discovered in the blood and plasma of humans. Correlations between tumor-associated genomic/epigenetic/transcriptional changes and alterations in CNA levels are strong predictors of the utility of this biomarker class as promising clinical indicators. Towards this goal microRNAs (miRNAs) representing a class of naturally occurring small non-coding RNAs of 19-25 nt in length have emerged as an important set of markers that can associate their specific expression profiles with cancer development. In this study we investigate some of the pre-analytic considerations for isolating plasma fractions for the study of miRNA biomarkers. We find that measurement of circulating miRNA levels are frequently confounded by varying levels of cellular miRNAs of different hematopoietic origins. In order to assess the relative proportions of this cell-derived class, we have fractionated whole blood into plasma and its ensuing sub-fractions. Cellular miRNA signatures in cohorts of normal individuals are catalogued and the abundance and gender specific expression of bona fide circulating markers explored after calibrating the signal for this interfering class. A map of differentially expressed profiles is presented and the intrinsic variability of circulating miRNA species investigated in subsets of healthy males and females.

Population-genetic nature of copy number variations in the human genome.

Copy number variations (CNVs) are universal genetic variations, and their association with disease has been increasingly recognized. We designed high-density microarrays for CNVs, and detected 3000-4000 CNVs (4-6% of the genomic sequence) per population that included CNVs previously missed because of smaller sizes and residing in segmental duplications. The patterns of CNVs across individuals were surprisingly simple at the kilo-base scale, suggesting the applicability of a simple genetic analysis for these genetic loci. We utilized the probabilistic theory to determine integer copy numbers of CNVs and employed a recently developed phasing tool to estimate the population frequencies of integer copy number alleles and CNV-SNP haplotypes. The results showed a tendency toward a lower frequency of CNV alleles and that most of our CNVs were explained only by zero-, one- and two-copy alleles. Using the estimated population frequencies, we found several CNV regions with exceptionally high population differentiation. Investigation of CNV-SNP linkage disequilibrium (LD) for 500-900 bi- and multi-allelic CNVs per population revealed that previous conflicting reports on bi-allelic LD were unexpectedly consistent and explained by an LD increase correlated with deletion-allele frequencies. Typically, the bi-allelic LD was lower than SNP-SNP LD, whereas the multi-allelic LD was somewhat stronger than the bi-allelic LD. After further investigation of tag SNPs for CNVs, we conclude that the customary tagging strategy for disease association studies can be applicable for common deletion CNVs, but direct interrogation is needed for other types of CNVs.

Integrated detection and population-genetic analysis of SNPs and copy number variation.

Dissecting the genetic basis of disease risk requires measuring all forms of genetic variation, including SNPs and copy number variants (CNVs), and is enabled by accurate maps of their locations, frequencies and population-genetic properties. We designed a hybrid genotyping array (Affymetrix SNP 6.0) to simultaneously measure 906,600 SNPs and copy number at 1.8 million genomic locations. By characterizing 270 HapMap samples, we developed a map of human CNV (at 2-kb breakpoint resolution) informed by integer genotypes for 1,320 copy number polymorphisms (CNPs) that segregate at an allele frequency >1%. More than 80% of the sequence in previously reported CNV regions fell outside our estimated CNV boundaries, indicating that large (>100 kb) CNVs affect much less of the genome than initially reported. Approximately 80% of observed copy number differences between pairs of individuals were due to common CNPs with an allele frequency >5%, and more than 99% derived from inheritance rather than new mutation. Most common, diallelic CNPs were in strong linkage disequilibrium with SNPs, and most low-frequency CNVs segregated on specific SNP haplotypes.

Improved detection of global copy number variation using high density, non-polymorphic oligonucleotide probes.

BACKGROUND: DNA sequence diversity within the human genome may be more greatly affected by copy number variations (CNVs) than single nucleotide polymorphisms (SNPs). Although the importance of CNVs in genome wide association studies (GWAS) is becoming widely accepted, the optimal methods for identifying these variants are still under evaluation. We have previously reported a comprehensive view of CNVs in the HapMap DNA collection using high density 500 K EA (Early Access) SNP genotyping arrays which revealed greater than 1,000 CNVs ranging in size from 1 kb to over 3 Mb. Although the arrays used most commonly for GWAS predominantly interrogate SNPs, CNV identification and detection does not necessarily require the use of DNA probes centered on polymorphic nucleotides and may even be hindered by the dependence on a successful SNP genotyping assay. RESULTS: In this study, we have designed and evaluated a high density array predicated on the use of non-polymorphic oligonucleotide probes for CNV detection. This approach effectively uncouples copy number detection from SNP genotyping and thus has the potential to significantly improve probe coverage for genome-wide CNV identification. This array, in conjunction with PCR-based, complexity-reduced DNA target, queries over 1.3 M independent NspI restriction enzyme fragments in the 200 bp to 1100 bp size range, which is a several fold increase in marker density as compared to the 500 K EA array. In addition, a novel algorithm was developed and validated to extract CNV regions and boundaries. CONCLUSION: Using a well-characterized pair of DNA samples, close to 200 CNVs were identified, of which nearly 50% appear novel yet were independently validated using quantitative PCR. The results indicate that non-polymorphic probes provide a robust approach for CNV identification, and the increasing precision of CNV boundary delineation should allow a more complete analysis of their genomic organization.

High-density oligonucleotide array with sub-kilobase resolution reveals breakpoint information of submicroscopic deletions in nevoid basal cell carcinoma syndrome.

Small submicroscopic genomic deletions and duplications constitute up to 15% of all mutations underlying human monogenic diseases. In this study, we used newly designed high-resolution oligonucleotide microarrays with a median distance between the probes of 776 bp (average probe interval 2,271 bp) to detect gene deletions in nevoid basal cell carcinoma syndrome (NBCCS) patients. NBCCS, also called Gorlin syndrome, is characterized by developmental defects and tumorigenesis such as medulloblastomas and basal cell carcinomas, caused by mutations of the human patched-1 (PTCH1) gene. Two out of three deletions could not be detected by a conventional chromosomal analysis. A submicroscopic deletion as small as 165 kb was detected affecting only PTCH1, whereas the other two deletions were much larger (5 and 11 Mb). We demonstrated not only the exact number of genes involved in the deletion but also rapidly determined the junction sequences after pinpointing the breakpoint regions in all individuals analyzed. This report of an array-based determination of junction sequences of long deletions circumvented a labor-intensive analysis such as Southern blotting or FISH. Alu-mediated recombination in one case and non-homologous end joining in the other two were probably implicated in the generation of deletions. This method will contribute to the understanding of molecular pathogenesis of gene deletions as well as rapid genetic testing.

Genome-wide detection of human copy number variations using high-density DNA oligonucleotide arrays.

Recent reports indicate that copy number variations (CNVs) within the human genome contribute to nucleotide diversity to a larger extent than single nucleotide polymorphisms (SNPs). In addition, the contribution of CNVs to human disease susceptibility may be greater than previously expected, although a complete understanding of the phenotypic consequences of CNVs is incomplete. We have recently reported a comprehensive view of CNVs among 270 HapMap samples using high-density SNP genotyping arrays and BAC array CGH. In this report, we describe a novel algorithm using Affymetrix GeneChip Human Mapping 500K Early Access (500K EA) arrays that identified 1203 CNVs ranging in size from 960 bp to 3.4 Mb. The algorithm consists of three steps: (1) Intensity pre-processing to improve the resolution between pairwise comparisons by directly estimating the allele-specific affinity as well as to reduce signal noise by incorporating probe and target sequence characteristics via an improved version of the Genomic Imbalance Map (GIM) algorithm; (2) CNV extraction using an adapted SW-ARRAY procedure to automatically and robustly detect candidate CNV regions; and (3) copy number inference in which all pairwise comparisons are summarized to more precisely define CNV boundaries and accurately estimate CNV copy number. Independent testing of a subset of CNVs by quantitative PCR and mass spectrometry demonstrated a >90% verification rate. The use of high-resolution oligonucleotide arrays relative to other methods may allow more precise boundary information to be extracted, thereby enabling a more accurate analysis of the relationship between CNVs and other genomic features.

Global variation in copy number in the human genome.

Copy number variation (CNV) of DNA sequences is functionally significant but has yet to be fully ascertained. We have constructed a first-generation CNV map of the human genome through the study of 270 individuals from four populations with ancestry in Europe, Africa or Asia (the HapMap collection). DNA from these individuals was screened for CNV using two complementary technologies: single-nucleotide polymorphism (SNP) genotyping arrays, and clone-based comparative genomic hybridization. A total of 1,447 copy number variable regions (CNVRs), which can encompass overlapping or adjacent gains or losses, covering 360 megabases (12% of the genome) were identified in these populations. These CNVRs contained hundreds of genes, disease loci, functional elements and segmental duplications. Notably, the CNVRs encompassed more nucleotide content per genome than SNPs, underscoring the importance of CNV in genetic diversity and evolution. The data obtained delineate linkage disequilibrium patterns for many CNVs, and reveal marked variation in copy number among populations. We also demonstrate the utility of this resource for genetic disease studies.

Copy number variation: new insights in genome diversity.

DNA copy number variation has long been associated with specific chromosomal rearrangements and genomic disorders, but its ubiquity in mammalian genomes was not fully realized until recently. Although our understanding of the extent of this variation is still developing, it seems likely that, at least in humans, copy number variants (CNVs) account for a substantial amount of genetic variation. Since many CNVs include genes that result in differential levels of gene expression, CNVs may account for a significant proportion of normal phenotypic variation. Current efforts are directed toward a more comprehensive cataloging and characterization of CNVs that will provide the basis for determining how genomic diversity impacts biological function, evolution, and common human diseases.

Noise reduction from genotyping microarrays using probe level information.

Genomic copy number change is one of the important phenomenon observed in cancer and other genetic disorders. Recently oligonucleotide microarrays have been used to analyze changes in the copy number. Although high density microarrays provide genome wide useful data on copy number, they are often associated with substantial amount of experimental noise that could affect the performance of the analyses. We used the high density oligonucleotide genotyping microarrays in our experiments that uses redundant probe tiling approach for individual SNPs. We found that the noise in the genotyping microarray data is associated with several experimental steps during target preparation and devised an algorithm that takes into account those experimental parameters. Additionally, defective probes that do not hybridize well to the target and therefore could not be modified inherently were detected and omitted automatically by using the algorithm. When we applied the algorithm to actual datasets, we could reduce the noise substantially without compressing the dynamic range. Additionally, combinatorial use of our noise reduction algorithm and conventional breakpoint detection algorithm successfully detected a microamplification of c-myc which was overlooked in the raw data. The algorithm described here is freely available with the software upon request to all non-profit researchers.

CARAT: a novel method for allelic detection of DNA copy number changes using high density oligonucleotide arrays.

BACKGROUND: DNA copy number alterations are one of the main characteristics of the cancer cell karyotype and can contribute to the complex phenotype of these cells. These alterations can lead to gains in cellular oncogenes as well as losses in tumor suppressor genes and can span small intervals as well as involve entire chromosomes. The ability to accurately detect these changes is central to understanding how they impact the biology of the cell. RESULTS: We describe a novel algorithm called CARAT (Copy Number Analysis with Regression And Tree) that uses probe intensity information to infer copy number in an allele-specific manner from high density DNA oligonuceotide arrays designed to genotype over 100,000 SNPs. Total and allele-specific copy number estimations using CARAT are independently evaluated for a subset of SNPs using quantitative PCR and allelic TaqMan reactions with several human breast cancer cell lines. The sensitivity and specificity of the algorithm are characterized using DNA samples containing differing numbers of X chromosomes as well as a test set of normal individuals. Results from the algorithm show a high degree of agreement with results from independent verification methods. CONCLUSION: Overall, CARAT automatically detects regions with copy number variations and assigns a significance score to each alteration as well as generating allele-specific output. When coupled with SNP genotype calls from the same array, CARAT provides additional detail into the structure of genome wide alterations that can contribute to allelic imbalance.

Evaluating potential for whole-genome studies in Kosrae, an isolated population in Micronesia.

Whole-genome association studies are predicted to be especially powerful in isolated populations owing to increased linkage disequilibrium (LD) and decreased allelic diversity, but this possibility has not been empirically tested. We compared genome-wide data on 113,240 SNPs typed on 30 trios from the Pacific island of Kosrae to the same markers typed in the 270 samples from the International HapMap Project. The extent of LD is longer and haplotype diversity is lower in Kosrae than in the HapMap populations. More than 98% of Kosraen haplotypes are present in HapMap populations, indicating that HapMap will be useful for genetic studies on Kosrae. The long-range LD around common alleles and limited diversity result in improved efficiency in genetic studies in this population and augments the power to detect association of 'hidden SNPs'.

Large-scale SNP analysis reveals clustered and continuous patterns of human genetic variation.

Understanding the distribution of human genetic variation is an important foundation for research into the genetics of common diseases. Some of the alleles that modify common disease risk are themselves likely to be common and, thus, amenable to identification using gene-association methods. A problem with this approach is that the large sample sizes required for sufficient statistical power to detect alleles with moderate effect make gene-association studies susceptible to false-positive findings as the result of population stratification. Such type I errors can be eliminated by using either family-based association tests or methods that sufficiently adjust for population stratification. These methods require the availability of genetic markers that can detect and, thus, control for sources of genetic stratification among populations. In an effort to investigate population stratification and identify appropriate marker panels, we have analysed 11,555 single nucleotide polymorphisms in 203 individuals from 12 diverse human populations. Individuals in each population cluster to the exclusion of individuals from other populations using two clustering methods. Higher-order branching and clustering of the populations are consistent with the geographic origins of populations and with previously published genetic analyses. These data provide a valuable resource for the definition of marker panels to detect and control for population stratification in population-based gene identification studies. Using three US resident populations (European-American, African-American and Puerto Rican), we demonstrate how such studies can proceed, quantifying proportional ancestry levels and detecting significant admixture structure in each of these populations.