Molecularly defined and spatially resolved cell atlas of the whole mouse brain.

In mammalian brains, millions to billions of cells form complex interaction networks to enable a wide range of functions. The enormous diversity and intricate organization of cells have impeded our understanding of the molecular and cellular basis of brain function. Recent advances in spatially resolved single-cell transcriptomics have enabled systematic mapping of the spatial organization of molecularly defined cell types in complex tissues1-3, including several brain regions (for example, refs. 1-11). However, a comprehensive cell atlas of the whole brain is still missing. Here we imaged a panel of more than 1,100 genes in approximately 10 million cells across the entire adult mouse brains using multiplexed error-robust fluorescence in situ hybridization12 and performed spatially resolved, single-cell expression profiling at the whole-transcriptome scale by integrating multiplexed error-robust fluorescence in situ hybridization and single-cell RNA sequencing data. Using this approach, we generated a comprehensive cell atlas of more than 5,000 transcriptionally distinct cell clusters, belonging to more than 300 major cell types, in the whole mouse brain with high molecular and spatial resolution. Registration of this atlas to the mouse brain common coordinate framework allowed systematic quantifications of the cell-type composition and organization in individual brain regions. We further identified spatial modules characterized by distinct cell-type compositions and spatial gradients featuring gradual changes of cells. Finally, this high-resolution spatial map of cells, each with a transcriptome-wide expression profile, allowed us to infer cell-type-specific interactions between hundreds of cell-type pairs and predict molecular (ligand-receptor) basis and functional implications of these cell-cell interactions. These results provide rich insights into the molecular and cellular architecture of the brain and a foundation for functional investigations of neural circuits and their dysfunction in health and disease.

A prebiotic basis for ATP as the universal energy currency.

ATP is universally conserved as the principal energy currency in cells, driving metabolism through phosphorylation and condensation reactions. Such deep conservation suggests that ATP arose at an early stage of biochemical evolution. Yet purine synthesis requires 6 phosphorylation steps linked to ATP hydrolysis. This autocatalytic requirement for ATP to synthesize ATP implies the need for an earlier prebiotic ATP equivalent, which could drive protometabolism before purine synthesis. Why this early phosphorylating agent was replaced, and specifically with ATP rather than other nucleoside triphosphates, remains a mystery. Here, we show that the deep conservation of ATP might reflect its prebiotic chemistry in relation to another universally conserved intermediate, acetyl phosphate (AcP), which bridges between thioester and phosphate metabolism by linking acetyl CoA to the substrate-level phosphorylation of ADP. We confirm earlier results showing that AcP can phosphorylate ADP to ATP at nearly 20% yield in water in the presence of Fe3+ ions. We then show that Fe3+ and AcP are surprisingly favoured. A wide range of prebiotically relevant ions and minerals failed to catalyse ADP phosphorylation. From a panel of prebiotic phosphorylating agents, only AcP, and to a lesser extent carbamoyl phosphate, showed any significant phosphorylating potential. Critically, AcP did not phosphorylate any other nucleoside diphosphate. We use these data, reaction kinetics, and molecular dynamic simulations to infer a possible mechanism. Our findings might suggest that the reason ATP is universally conserved across life is that its formation is chemically favoured in aqueous solution under mild prebiotic conditions.

The Giardia ventrolateral flange is a lamellar membrane protrusion that supports attachment.

Attachment to the intestinal epithelium is critical to the lifestyle of the ubiquitous parasite Giardia lamblia. The ventrolateral flange is a sheet-like membrane protrusion at the interface between parasites and attached surfaces. This structure has been implicated in attachment, but its role has been poorly defined. Here, we identified a novel actin associated protein with putative WH2-like actin binding domains we named Flangin. Flangin complexes with Giardia actin (GlActin) and is enriched in the ventrolateral flange making it a valuable marker for studying the flanges' role in Giardia biology. Live imaging revealed that the flange grows to around 1 µm in width after cytokinesis, then remains uniform in size during interphase, grows in mitosis, and is resorbed during cytokinesis. A flangin truncation mutant stabilizes the flange and blocks cytokinesis, indicating that flange disassembly is necessary for rapid myosin-independent cytokinesis in Giardia. Rho family GTPases are important regulators of membrane protrusions and GlRac, the sole Rho family GTPase in Giardia, was localized to the flange. Knockdown of Flangin, GlActin, and GlRac result in flange formation defects. This indicates a conserved role for GlRac and GlActin in forming membrane protrusions, despite the absence of canonical actin binding proteins that link Rho GTPase signaling to lamellipodia formation. Flangin-depleted parasites had reduced surface contact and when challenged with fluid shear force in flow chambers they had a reduced ability to remain attached, confirming a role for the flange in attachment. This secondary attachment mechanism complements the microtubule based adhesive ventral disc, a feature that may be particularly important during mitosis when the parental ventral disc disassembles in preparation for cytokinesis. This work supports the emerging view that Giardia's unconventional actin cytoskeleton has an important role in supporting parasite attachment.

Multiplexed single-cell profiling of chromatin states at genomic loci by expansion microscopy.

Proper regulation of genome architecture and activity is essential for the development and function of multicellular organisms. Histone modifications, acting in combination, specify these activity states at individual genomic loci. However, the methods used to study these modifications often require either a large number of cells or are limited to targeting one histone mark at a time. Here, we developed a new method called Single Cell Evaluation of Post-TRanslational Epigenetic Encoding (SCEPTRE) that uses Expansion Microscopy (ExM) to visualize and quantify multiple histone modifications at non-repetitive genomic regions in single cells at a spatial resolution of ∼75 nm. Using SCEPTRE, we distinguished multiple histone modifications at a single housekeeping gene, quantified histone modification levels at multiple developmentally-regulated genes in individual cells, and evaluated the relationship between histone modifications and RNA polymerase II loading at individual loci. We find extensive variability in epigenetic states between individual gene loci hidden from current population-averaged measurements. These findings establish SCEPTRE as a new technique for multiplexed detection of combinatorial chromatin states at single genomic loci in single cells.

Strelka2: fast and accurate calling of germline and somatic variants.

We describe Strelka2 ( https://github.com/Illumina/strelka ), an open-source small-variant-calling method for research and clinical germline and somatic sequencing applications. Strelka2 introduces a novel mixture-model-based estimation of insertion/deletion error parameters from each sample, an efficient tiered haplotype-modeling strategy, and a normal sample contamination model to improve liquid tumor analysis. For both germline and somatic calling, Strelka2 substantially outperformed the current leading tools in terms of both variant-calling accuracy and computing cost.

A reference data set of 5.4 million phased human variants validated by genetic inheritance from sequencing a three-generation 17-member pedigree.

Improvement of variant calling in next-generation sequence data requires a comprehensive, genome-wide catalog of high-confidence variants called in a set of genomes for use as a benchmark. We generated deep, whole-genome sequence data of 17 individuals in a three-generation pedigree and called variants in each genome using a range of currently available algorithms. We used haplotype transmission information to create a phased "Platinum" variant catalog of 4.7 million single-nucleotide variants (SNVs) plus 0.7 million small (1-50 bp) insertions and deletions (indels) that are consistent with the pattern of inheritance in the parents and 11 children of this pedigree. Platinum genotypes are highly concordant with the current catalog of the National Institute of Standards and Technology for both SNVs (>99.99%) and indels (99.92%) and add a validated truth catalog that has 26% more SNVs and 45% more indels. Analysis of 334,652 SNVs that were consistent between informatics pipelines yet inconsistent with haplotype transmission ("nonplatinum") revealed that the majority of these variants are de novo and cell-line mutations or reside within previously unidentified duplications and deletions. The reference materials from this study are a resource for objective assessment of the accuracy of variant calls throughout genomes.

Spinal muscular atrophy diagnosis and carrier screening from genome sequencing data.

PURPOSE: Spinal muscular atrophy (SMA), caused by loss of the SMN1 gene, is a leading cause of early childhood death. Due to the near identical sequences of SMN1 and SMN2, analysis of this region is challenging. Population-wide SMA screening to quantify the SMN1 copy number (CN) is recommended by the American College of Medical Genetics and Genomics. METHODS: We developed a method that accurately identifies the CN of SMN1 and SMN2 using genome sequencing (GS) data by analyzing read depth and eight informative reference genome differences between SMN1/2. RESULTS: We characterized SMN1/2 in 12,747 genomes, identified 1568 samples with SMN1 gains or losses and 6615 samples with SMN2 gains or losses, and calculated a pan-ethnic carrier frequency of 2%, consistent with previous studies. Additionally, 99.8% of our SMN1 and 99.7% of SMN2 CN calls agreed with orthogonal methods, with a recall of 100% for SMA and 97.8% for carriers, and a precision of 100% for both SMA and carriers. CONCLUSION: This SMN copy-number caller can be used to identify both carrier and affected status of SMA, enabling SMA testing to be offered as a comprehensive test in neonatal care and an accurate carrier screening tool in GS sequencing projects.

Expansion microscopy with conventional antibodies and fluorescent proteins.

Expansion microscopy is a technique in which fluorophores on fixed specimens are linked to a swellable polymer that is physically expanded to enable super-resolution microscopy with ordinary microscopes. We have developed and characterized new methods for linking fluorophores to the polymer that now enable expansion microscopy with conventional fluorescently labeled antibodies and fluorescent proteins. Our methods simplify the procedure and expand the palette of compatible labels, allowing rapid dissemination of the technique.

Extended-Depth 3D Super-Resolution Imaging Using Probe-Refresh STORM.

Single-molecule localization microscopy methods for super-resolution fluorescence microscopy such as STORM (stochastic optical reconstruction microscopy) are generally limited to thin three-dimensional (3D) sections (≤600 nm) because of photobleaching of molecules outside the focal plane. Although multiple focal planes may be imaged before photobleaching by focusing progressively deeper within the sample, image quality is compromised in this approach because the total number of measurable localizations is divided between detection planes. Here, we solve this problem on fixed samples by developing an imaging method that we call probe-refresh STORM (prSTORM), which allows bleached fluorophores to be straightforwardly replaced with nonbleached fluorophores. We accomplish this by immunostaining the sample with DNA-conjugated antibodies and then reading out their distribution using fluorescently-labeled DNA-reporter oligonucleotides that can be fully replaced in successive rounds of imaging. We demonstrate that prSTORM can acquire 3D images over extended depths without sacrificing the density of localizations at any given plane. We also show that prSTORM can be adapted to obtain high-quality, 3D multichannel images with extended depth that would be challenging or impossible to achieve using established probe methods.

Accurate whole-genome sequencing and haplotyping from 10 to 20 human cells.

Recent advances in whole-genome sequencing have brought the vision of personal genomics and genomic medicine closer to reality. However, current methods lack clinical accuracy and the ability to describe the context (haplotypes) in which genome variants co-occur in a cost-effective manner. Here we describe a low-cost DNA sequencing and haplotyping process, long fragment read (LFR) technology, which is similar to sequencing long single DNA molecules without cloning or separation of metaphase chromosomes. In this study, ten LFR libraries were made using only ∼100 picograms of human DNA per sample. Up to 97% of the heterozygous single nucleotide variants were assembled into long haplotype contigs. Removal of false positive single nucleotide variants not phased by multiple LFR haplotypes resulted in a final genome error rate of 1 in 10 megabases. Cost-effective and accurate genome sequencing and haplotyping from 10-20 human cells, as demonstrated here, will enable comprehensive genetic studies and diverse clinical applications.

Evaluation of somatic copy number variation detection by NGS technologies and bioinformatics tools on a hyper-diploid cancer genome.

BACKGROUND: Copy number variation (CNV) is a key genetic characteristic for cancer diagnostics and can be used as a biomarker for the selection of therapeutic treatments. Using data sets established in our previous study, we benchmark the performance of cancer CNV calling by six most recent and commonly used software tools on their detection accuracy, sensitivity, and reproducibility. In comparison to other orthogonal methods, such as microarray and Bionano, we also explore the consistency of CNV calling across different technologies on a challenging genome. RESULTS: While consistent results are observed for copy gain, loss, and loss of heterozygosity (LOH) calls across sequencing centers, CNV callers, and different technologies, variation of CNV calls are mostly affected by the determination of genome ploidy. Using consensus results from six CNV callers and confirmation from three orthogonal methods, we establish a high confident CNV call set for the reference cancer cell line (HCC1395). CONCLUSIONS: NGS technologies and current bioinformatics tools can offer reliable results for detection of copy gain, loss, and LOH. However, when working with a hyper-diploid genome, some software tools can call excessive copy gain or loss due to inaccurate assessment of genome ploidy. With performance matrices on various experimental conditions, this study raises awareness within the cancer research community for the selection of sequencing platforms, sample preparation, sequencing coverage, and the choice of CNV detection tools.

Electrodeposition of polydopamine thin films for DNA patterning and microarrays.

The controlled electrodeposition of functional polydopamine (PDA) thin films from aqueous dopamine solutions is demonstrated with a combination of electrochemistry, atomic force microscopy (AFM), and surface plasmon resonance (SPR) measurements. PDA micropatterns are then fabricated by electrodeposition on micrometer length scale gold electrodes and used for attaching amino-modified single-stranded DNA (ssDNA). After hybridization with fluorescently labeled ssDNA, the fluorescence microscopy characterization reveals that: (i) PDA can be toposelectively deposited at the microscale and (ii) electrochemically deposited PDA can be functionalized with amino-terminated ssDNA using the same chemistry as that for spontaneously deposited PDA. Finally, the application of electrodeposited PDA thin films to fabricate ssDNA microarrays is reported using SPR imaging (SPRI) measurements for the detection of DNA and DNA-modified gold nanoparticles.

Fabrication of DNA microarrays on polydopamine-modified gold thin films for SPR imaging measurements.

Polydopamine (PDA) films were fabricated on thin film gold substrates in a single-step polymerization-deposition process from dopamine solutions and then employed in the construction of robust DNA microarrays for the ultrasensitive detection of biomolecules with nanoparticle-enhanced surface plasmon resonance (SPR) imaging. PDA multilayers with thicknesses varying from 1 to 5 nm were characterized with a combination of scanning angle SPR and AFM experiments, and 1.3 ± 0.2 nm PDA multilayers were chosen as an optimal thickness for the SPR imaging measurements. DNA microarrays were then fabricated by the reaction of amine-functionalized single-stranded DNA (ssDNA) oligonucleotides with PDA-modified gold thin film microarray elements, and were subsequently employed in SPR imaging measurements of DNA hybridization adsorption and protein-DNA binding. Concurrent control experiments with non-complementary ssDNA sequences demonstrated that the adhesive PDA multilayer was also able to provide good resistance to the nonspecific binding of biomolecules. Finally, a series of SPR imaging measurements of the hybridization adsorption of DNA-modified gold nanoparticles onto mixed sequence DNA microarrays were used to confirm that the use of PDA multilayer films is a simple, rapid, and versatile method for fabricating DNA microarrays for ultrasensitive nanoparticle-enhanced SPR imaging biosensing.

Characterization of Prochlorococcus clades from iron-depleted oceanic regions.

Prochlorococcus describes a diverse and abundant genus of marine photosynthetic microbes. It is primarily found in oligotrophic waters across the globe and plays a crucial role in energy and nutrient cycling in the ocean ecosystem. The abundance, global distribution, and availability of isolates make Prochlorococcus a model system for understanding marine microbial diversity and biogeochemical cycling. Analysis of 73 metagenomic samples from the Global Ocean Sampling expedition acquired in the Atlantic, Pacific, and Indian Oceans revealed the presence of two uncharacterized Prochlorococcus clades. A phylogenetic analysis using six different genetic markers places the clades close to known lineages adapted to high-light environments. The two uncharacterized clades consistently cooccur and dominate the surface waters of high-temperature, macronutrient-replete, and low-iron regions of the Eastern Equatorial Pacific upwelling and the tropical Indian Ocean. They are genetically distinct from each other and other high-light Prochlorococcus isolates and likely define a previously unrecognized ecotype. Our detailed genomic analysis indicates that these clades comprise organisms that are adapted to iron-depleted environments by reducing their iron quota through the loss of several iron-containing proteins that likely function as electron sinks in the photosynthetic pathway in other Prochlorococcus clades from high-light environments. The presence and inferred physiology of these clades may explain why Prochlorococcus populations from iron-depleted regions do not respond to iron fertilization experiments and further expand our understanding of how phytoplankton adapt to variations in nutrient availability in the ocean.

Massive parallelism, randomness and genomic advances.

In reviewing the past decade, it is clear that genomics was, and still is, driven by innovative technologies, perhaps more so than any other scientific area in recent memory. From the outset, computing, mathematics and new automated laboratory techniques have been key components in allowing the field to move forward rapidly. We highlight some key innovations that have come together to nurture the explosive growth that makes a new era of genomics a reality. We also document how these new approaches have fueled further innovations and discoveries.

Biophysical Interactions Underpin the Emergence of Information in the Genetic Code.

The genetic code conceals a 'code within the codons', which hints at biophysical interactions between amino acids and their cognate nucleotides. Yet, research over decades has failed to corroborate systematic biophysical interactions across the code. Using molecular dynamics simulations and NMR, we have analysed interactions between the 20 standard proteinogenic amino acids and 4 RNA mononucleotides in 3 charge states. Our simulations show that 50% of amino acids bind best with their anticodonic middle base in the -1 charge state common to the backbone of RNA, while 95% of amino acids interact most strongly with at least 1 of their codonic or anticodonic bases. Preference for the cognate anticodonic middle base was greater than 99% of randomised assignments. We verify a selection of our results using NMR, and highlight challenges with both techniques for interrogating large numbers of weak interactions. Finally, we extend our simulations to a range of amino acids and dinucleotides, and corroborate similar preferences for cognate nucleotides. Despite some discrepancies between the predicted patterns and those observed in biology, the existence of weak stereochemical interactions means that random RNA sequences could template non-random peptides. This offers a compelling explanation for the emergence of genetic information in biology.

A molecularly defined and spatially resolved cell atlas of the whole mouse brain.

In mammalian brains, tens of millions to billions of cells form complex interaction networks to enable a wide range of functions. The enormous diversity and intricate organization of cells in the brain have so far hindered our understanding of the molecular and cellular basis of its functions. Recent advances in spatially resolved single-cell transcriptomics have allowed systematic mapping of the spatial organization of molecularly defined cell types in complex tissues1-3. However, these approaches have only been applied to a few brain regions1-11 and a comprehensive cell atlas of the whole brain is still missing. Here, we imaged a panel of >1,100 genes in ~8 million cells across the entire adult mouse brain using multiplexed error-robust fluorescence in situ hybridization (MERFISH)12 and performed spatially resolved, single-cell expression profiling at the whole-transcriptome scale by integrating MERFISH and single-cell RNA-sequencing (scRNA-seq) data. Using this approach, we generated a comprehensive cell atlas of >5,000 transcriptionally distinct cell clusters, belonging to ~300 major cell types, in the whole mouse brain with high molecular and spatial resolution. Registration of the MERFISH images to the common coordinate framework (CCF) of the mouse brain further allowed systematic quantifications of the cell composition and organization in individual brain regions defined in the CCF. We further identified spatial modules characterized by distinct cell-type compositions and spatial gradients featuring gradual changes in the gene-expression profiles of cells. Finally, this high-resolution spatial map of cells, with a transcriptome-wide expression profile associated with each cell, allowed us to infer cell-type-specific interactions between several hundred pairs of molecularly defined cell types and predict potential molecular (ligand-receptor) basis and functional implications of these cell-cell interactions. These results provide rich insights into the molecular and cellular architecture of the brain and a valuable resource for future functional investigations of neural circuits and their dysfunction in diseases.

On-chip synthesis of protein microarrays from DNA microarrays via coupled in vitro transcription and translation for surface plasmon resonance imaging biosensor applications.

Protein microarrays are fabricated from double-stranded DNA (dsDNA) microarrays by a one-step, multiplexed enzymatic synthesis in an on-chip microfluidic format and then employed for antibody biosensing measurements with surface plasmon resonance imaging (SPRI). A microarray of dsDNA elements (denoted as generator elements) that encode either a His-tagged green fluorescent protein (GFP) or a His-tagged luciferase protein is utilized to create multiple copies of mRNA (mRNA) in a surface RNA polymerase reaction; the mRNA transcripts are then translated into proteins by cell-free protein synthesis in a microfluidic format. The His-tagged proteins diffuse to adjacent Cu(II)-NTA microarray elements (denoted as detector elements) and are specifically adsorbed. The net result is the on-chip, cell-free synthesis of a protein microarray that can be used immediately for SPRI protein biosensing. The dual element format greatly reduces any interference from the nonspecific adsorption of enzyme or proteins. SPRI measurements for the detection of the antibodies anti-GFP and antiluciferase were used to verify the formation of the protein microarray. This convenient on-chip protein microarray fabrication method can be implemented for multiplexed SPRI biosensing measurements in both clinical and research applications.

Wafer-scale fabrication of nanofluidic arrays and networks using nanoimprint lithography and lithographically patterned nanowire electrodeposition gold nanowire masters.

Wafer scale (cm(2)) arrays and networks of nanochannels were created in polydimethylsiloxane (PDMS) from a surface pattern of electrodeposited gold nanowires in a master-replica process and characterized with scanning electron microscopy (SEM), atomic force microscopy (AFM), and fluorescence imaging measurements. Patterns of gold nanowires with cross-sectional dimensions as small as 50 nm in height and 100 nm in width were prepared on silica substrates using the process of lithographically patterned nanowire electrodeposition (LPNE). These nanowire patterns were then employed as masters for the fabrication of inverse replica nanochannels in a special formulation of PDMS. SEM and AFM measurements verified a linear correlation between the widths and heights of the nanowires and nanochannels over a range of 50 to 500 nm. The PDMS replica was then oxygen plasma-bonded to a glass substrate in order to create a linear array of nanofluidic channels (up to 1 mm in length) filled with solutions of either fluorescent dye or 20 nm diameter fluorescent polymer nanoparticles. Nanochannel continuity and a 99% fill success rate was determined from the fluorescence imaging measurements, and the electrophoretic injection of both dye and nanoparticles in the nanochannel arrays was also demonstrated. Employing a double LPNE fabrication method, this master-replica process was also used to create a large two-dimensional network of crossed nanofluidic channels.

Near infrared surface plasmon resonance phase imaging and nanoparticle-enhanced surface plasmon resonance phase imaging for ultrasensitive protein and DNA biosensing with oligonucleotide and aptamer microarrays.

The techniques of surface plasmon resonance-phase imaging (SPR-PI) and nanoparticle-enhanced SPR-PI have been implemented for the multiplexed bioaffinity detection of proteins and nucleic acids. The SPR-PI experiments utilized a near-infrared 860 nm light emitting diode (LED) light source and a wedge depolarizer to create a phase grating on a four-element single-stranded DNA (ssDNA) microarray; bioaffinity adsorption onto the various microarray elements was detected via multiplexed real time phase shift measurements. In a first set of demonstration experiments, an ssDNA aptamer microarray was used to directly detect thrombin at concentrations down to 100 pM with SPR-PI. Two different ssDNA aptamers were used in these experiments with two different Langmuir adsorption coefficients, K(A1) = 4.4 × 10(8) M(-1) and K(A2) = 1.2 × 10(8) M(-1). At concentrations below 1 nM, the equilibrium phase shifts observed upon thrombin adsorption vary linearly with concentration with a slope that is proportional to the appropriate Langmuir adsorption coefficient. The observed detection limit of 100 pM is approximately 20 times more sensitive than that observed previously with SPRI. In a second set of experiments, two short ssDNA oligonucleotides (38mers) were simultaneously detected at concentrations down to 25 fM using a three-sequence hybridization format that employed 120 nm DNA-modified silica nanoparticles to enhance the SPR-PI signal. In this first demonstration of nanoparticle-enhanced SPR-PI, the adsorbed silica nanoparticles provided a greatly enhanced phase shift upon bioaffinity adsorption due to a large increase in the real component of the interfacial refractive index from the adsorbed nanoparticle. As in the case of SPR-PI, the detection limit of 25 fM for nanoparticle-enhanced SPR-PI is approximately 20 times more sensitive than that observed previously with nanoparticle-enhanced SPRI.