Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “sequence analysis”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,513 records · Page 84

Functional importance of cardiac enhancer-associated noncoding RNAs in heart development and disease

We report here that the key information processing units within gene regulatory networks are enhancers. Enhancer activity is associated with the production of tissue-specific noncoding RNAs, yet the existence of such transcripts during cardiac development has not been established. Using an integrated genomic approach, we demonstrate that fetal cardiac enhancers generate long noncoding RNAs (lncRNAs) during cardiac differentiation and morphogenesis. Enhancer expression correlates with the emergence of active enhancer chromatin states, the initiation of RNA polymerase II at enhancer loci and expression of target genes. Orthologous human sequences are also transcribed in fetal human hearts and cardiac progenitor cells. Through a systematic bioinformatic analysis, we identified and characterized, for the first time, a catalog of lncRNAs that are expressed during embryonic stem cell differentiation into cardiomyocytes and associated with active cardiac enhancer sequences. RNA-sequencing demonstrates that many of these transcripts are polyadenylated, multi-exonic long noncoding RNAs. Moreover, knockdown of two enhancer-associated lncRNAs resulted in the specific downregulation of their predicted target genes. Interestingly, the reactivation of the fetal gene program, a hallmark of the stress response in the adult heart, is accompanied by increased expression of fetal cardiac enhancer transcripts. Altogether, these findings demonstrate that the activity of cardiac enhancers and expression of their target genes are associated with the production of enhancer-derived lncRNAs.

59 BASIC BIOLOGICAL SCIENCES↗

Internal Fragments Generated by Electron Ionization Dissociation Enhance Protein Top-Down Mass Spectrometry

Top-down proteomics by mass spectrometry (MS) involves the mass measurement of an intact protein followed by subsequent activation of the protein to generate product ions. Electron-based fragmentation methods like electron capture dissociation and electron transfer dissociation are widely used for these types of analyses. Recently, electron ionization dissociation (EID), which utilizes higher energy electrons (>20 eV) has been suggested to be more efficient for top-down protein fragmentation compared to other electron-based dissociation methods. In this paper, we demonstrate that the use of EID enhances protein fragmentation and subsequent detection of protein fragments. Protein product ions can form by either single cleavage events, resulting in terminal fragments containing the C-terminus or N-terminus of the protein, or by multiple cleavage events to give rise to internal fragments that include neither the C-terminus nor the N-terminus of the protein. Conventionally, internal fragments have been disregarded, as reliable assignments of these fragments were limited. Here, we demonstrate that internal fragments generated by EID can account for ~20–40% of the mass spectral signals detected by top-down EID-MS experiments. By including internal fragments, the extent of the protein sequence that can be explained from a single tandem mass spectrum increases from ~50 to ~99% for 29 kDa carbonic anhydrase II and 8.6 kDa ubiquitin. When searching for internal fragments during data analysis, previously unassigned peaks can be readily and accurately assigned to confirm a given protein sequence and to enhance the utility of top-down protein sequencing experiments.

59 BASIC BIOLOGICAL SCIENCES↗

A sequence of transformations related to the formation of M{sub 3}X{sub 2}-type superstructures

A symmetry analysis of monoclinic, orthorhombic, and trigonal M{sub 3}X{sub 2}-type superstructures that can be formed in strongly stoichiometric MX{sub y} compounds with B1 structure is carried out. Channels of order-disorder transitions MX{sub y} → M{sub 3}X{sub 2} are determined. It is shown that, as temperature decreases, two physically admissible sequences of transformations associated with the formation of M{sub 3}X{sub 2} superstructures are possible in nonstoichiometric MX{sub y} compounds of group IV transition metals. By an example of vanadium carbide VC{sub y}, it is demonstrated that orthorhombic or monoclinic V{sub 3}C{sub 2} superstructures can be obtained with the formation of a nanostructure.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Xavier: High-Performance X-Drop Adaptive Banded Pairwise Alignment (Xavier) v1.0

Pairwise sequence alignment is one of the most compute-intense challenges for the analysis of genomics data and it can account for more than 90% of the runtime depending on the application. Despite the effort to reduce the computational cost introducing heuristics such as banded and x-drop alignment, the performance of state-of-the-art software is often unsatisfactory for long sequences. Recent advances in single-molecule sequencing technologies keep increasing the length of the sequences to align at the cost of higher error rates b ranging from 5% to 15% b, emphasizing the need for fast and accurate software. Our long-read alignment software, named Xavier, computes high-performance x-drop adaptive banded pairwise alignment exploiting instruction-level parallelism via Single-Instruction-Multiple-Data (SIMD) operations. Xavier uses a narrow bandwidth that appreciably improves performance reducing the search space for the optimal alignment. In our experiment, the adaptive mechanism enables Xavier to find the optimal alignment between related sequences even when the error rate is as high as 50%. Xavier achieves speed-up up to 5X over competitor software when enabling the x-drop termination and up to 415X when disabling it, with negligible accuracy loss.

Guidi, Giulia↗

Persistent minimal sequences of SARS-CoV-2

Abstract Motivation Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has caused more than 14 million cases and more than half million deaths. Given the absence of implemented therapies, new analysis, diagnosis and therapeutics are of great importance. Results Analysis of SARS-CoV-2 genomes from the current outbreak reveals the presence of short persistent DNA/RNA sequences that are absent from the human genome and transcriptome (PmRAWs). For the PmRAWs with length 12, only four exist at the same location in all SARS-CoV-2. At the gene level, we found one PmRAW of size 13 at the Spike glycoprotein coding sequence. This protein is fundamental for binding in human ACE2 and further use as an entry receptor to invade target cells. Applying protein structural prediction, we localized this PmRAW at the surface of the Spike protein, providing a potential targeted vector for diagnostics and therapeutics. In addition, we show a new pattern of relative absent words (RAWs), characterized by the progressive increase of GC content (Guanine and Cytosine) according to the decrease of RAWs length, contrarily to the virus and host genome distributions. New analysis shows the same property during the Ebola virus outbreak. At a computational level, we improved the alignment-free method to identify pathogen-specific signatures in balance with GC measures and removed previous size limitations. Availability and implementation https://github.com/cobilab/eagle. Supplementary information Supplementary data are available at Bioinformatics online.

Pratas, Diogo↗

Chicago Area Waterway System Microbiome Research

This report summarizes the research conducted by Argonne National Laboratory (Argonne) for the Metropolitan Water Reclamation District of Greater Chicago (MWRD) for the first five years (2013-2017) of a seven-year study with emphasis on the impact of disinfection (2016-2017) on the microbial community of the Chicago Area Waterway System (CAWS). The study examines microbial communities following disinfection of secondary treated effluents at the T.J. O’Brien (UV) and Calumet (chlorine and dechlorination) Water Reclamation Plants (WRPs) and the phased implementation of the Tunnel and Reservoir Plan (TARP). The Thornton Composite Reservoir (TCR) was completed in 2015. It provides 7.9 billion gallons of storage and since its completion has captured more than 11.0 billion gallons of combined stormwater and sewage from Calumet WRP that would otherwise overflow into CAWS in rainy weather. In addition to the effects of TARP and disinfection, we present an analysis of potential sources of the microbes found in the CAWS and the results of 16S rRNA amplicon gene sequencing which was used to characterize microbial community variability as a function of location, season, and environmental conditions. A total of 2,077 samples were collected from 2013 to 2017 from WRPs and CAWS and analyzed using high-throughput sequencing. The whole microbial community composition using all the genetic material in the environmental sample was characterized using 16S rRNA gene-based analysis, while functional subsystem attributes of these microbial communities were characterized using shotgun metagenomics. We used a Bayesian statistical tool, SourceTracker 2.0 to determine the likely sources of microbes found in CAWS by sampling location and sampling period.

54 ENVIRONMENTAL SCIENCES↗

Complementary Metagenomic Approaches Improve Reconstruction of Microbial Diversity in a Forest Soil

ABSTRACT Soil ecosystems harbor diverse microorganisms and yet remain only partially characterized as neither single-cell sequencing nor whole-community sequencing offers a complete picture of these complex communities. Thus, the genetic and metabolic potential of this “uncultivated majority” remains underexplored. To address these challenges, we applied a pooled-cell-sorting-based mini-metagenomics approach and compared the results to bulk metagenomics. Informatic binning of these data produced 200 mini-metagenome assembled genomes (sorted-MAGs) and 29 bulk metagenome assembled genomes (MAGs). The sorted and bulk MAGs increased the known phylogenetic diversity of soil taxa by 7.2% with respect to the Joint Genome Institute IMG/M database and showed clade-specific sequence recruitment patterns across diverse terrestrial soil metagenomes. Additionally, sorted-MAGs expanded the rare biosphere not captured through MAGs from bulk sequences, exemplified through phylogenetic and functional analyses of members of the phylum Bacteroidetes . Analysis of 67 Bacteroidetes sorted-MAGs showed conserved patterns of carbon metabolism across four clades. These results indicate that mini-metagenomics enables genome-resolved investigation of predicted metabolism and demonstrates the utility of combining metagenomics methods to tap into the diversity of heterogeneous microbial assemblages. IMPORTANCE Microbial ecologists have historically used cultivation-based approaches as well as amplicon sequencing and shotgun metagenomics to characterize microbial diversity in soil. However, challenges persist in the study of microbial diversity, including the recalcitrance of the majority of microorganisms to laboratory cultivation and limited sequence assembly from highly complex samples. The uncultivated majority thus remains a reservoir of untapped genetic diversity. To address some of the challenges associated with bulk metagenomics as well as low throughput of single-cell genomics, we applied flow cytometry-enabled mini-metagenomics to capture expanded microbial diversity from forest soil and compare it to soil bulk metagenomics. Our resulting data from this pooled-cell sorting approach combined with bulk metagenomics revealed increased phylogenetic diversity through novel soil taxa and rare biosphere members. In-depth analysis of genomes within the highly represented Bacteroidetes phylum provided insights into conserved and clade-specific patterns of carbon metabolism.

59 BASIC BIOLOGICAL SCIENCES↗

Differential DNA methylation profile of key genes in malignant prostate epithelial cells transformed by inorganic arsenic or cadmium

Previous work shows altered methylation patterns in inorganic arsenic (iAs)- or cadmium (Cd)-transformed epithelial cells. Here, the methylation status near the transcriptional start site was assessed in the normal human prostate epithelial cell line (RWPE-1) that was malignantly transformed by 10 μM Cd for 11 weeks (CTPE) or 5 μM iAs for 29 weeks (CAsE-PE), at which time cells showed multiple markers of acquired cancer phenotype. Next generation sequencing of the transcriptome of CAsE-PE cells identified multiple dysregulated genes. Of the most highly dysregulated genes, five genes that can be relevant to the carcinogenic process (S100P, HYAL1, NTM, NES, ALDH1A1) were chosen for an in-depth analysis of the DNA methylation profile. DNA was isolated, bisulfite converted, and combined bisulfite restriction analysis was used to identify differentially methylated CpG sites, which was confirmed with bisulfite sequencing. Four of the five genes showed differential methylation in transformants relative to control cells that was inversely related to altered gene expression. Increased expression of HYAL1 (> 25-fold) and S100P (> 40-fold) in transformants was correlated with hypomethylation near the transcriptional start site. Decreased expression of NES (> 15-fold) and NTM (> 1000-fold) in transformants was correlated with hypermethylation near the transcriptional start site. ALDH1A1 expression was differentially expressed in transformed cells but was not differentially methylated relative to control. In conclusion, altered gene expression observed in Cd and iAs transformed cells may result from altered DNA methylation status. - Highlights: • Cd and iAs are known human carcinogens, yet neither appears directly mutagenic. • Prior data suggest epigenetic modification plays a role in Cd or iAs induced cancer. • Altered methylation of four misregulated genes was found in Cd or iAs transformants. • The resulting altered gene expression may be relevant to cellular transformation.

60 APPLIED LIFE SCIENCES↗

MATEY: multiscale adaptive transformer models for spatiotemporal physical systems

Accurate representation of the multiscale features in spatiotemporal physical systems using vision transformer architectures requires extremely long, computationally prohibitive token sequences. To address this issue, we propose two novel adaptive tokenization schemes that dynamically adjust patch sizes based on local features: one ensures convergent behavior to uniform patch refinement, while the other offers better computational efficiency. Moreover, we present a set of spatiotemporal attention schemes, where the temporal or axial spatial dimensions are decoupled, to evaluate their baseline computational and data efficiencies and to determine whether adaptive tokenization can improve this performance. We assess the performance of the proposed multiscale adaptive model, MATEY, in a sequence of experiments. Compared to a full spatiotemporal attention scheme or a scheme that decouples only the temporal dimension, we find that fully decoupled axial attention is less efficient and expressive, requiring more training time and model parameters to achieve the same accuracy. The experiments on the adaptive tokenization schemes show that, compared to a uniformly refined model, the proposed schemes achieve comparable or improved accuracy at a much lower cost in the tested two-dimensional settings. While the asymptotic analysis suggests the potential for favorable scaling, empirical validation at substantially longer sequence lengths remains to be performed in future work. Finally, we demonstrate in two fine-tuning tasks featuring different physics that models pretrained on PDEBench data outperform the ones trained from scratch, especially in the low data regime with frozen attention.

adaptive tokenization↗

Spectral energy distribution analysis of class I and class II FU Orionis stars

FU Orionis stars (FUors) are eruptive pre-main sequence objects thought to represent quasi-periodic or recurring stages of enhanced accretion during the low-mass star-forming process. We characterize the sample of known and candidate FUors in a homogeneous and consistent way, deriving stellar and circumstellar parameters for each object. We emphasize the analysis in those parameters that are supposed to vary during the FUor stage. We modeled the spectral energy distributions of 24 of the 26 currently known FUors, using the radiative transfer code of Whitney et al. We compare our models with those obtained by Robitaille et al. for Taurus class II and I sources in quiescence periods by calculating the cumulative distribution of the different parameters. FUors have more massive disks: we find that ∼80% of the disks in FUors are more massive than any Taurus class II and I sources in the sample. Median values for the disk mass accretion rates are ∼10{sup –7} M {sub ☉} yr{sup –1} versus ∼10{sup –5} M {sub ☉} yr{sup –1} for standard young stellar objects (YSOs) and FUors, respectively. While the distributions of envelope mass accretion rates for class I FUors and standard class I objects are similar, FUors, on average, have higher envelope mass accretion rates than standard class II and class I sources. Most FUors (∼70%) have envelope mass accretion rates above 10{sup –7} M {sub ☉} yr{sup –1}. In contrast, 60% of the classical YSO sample has an accretion rate below this value. Our results support the current scenario in which changes experimented by the circumstellar disk explain the observed properties of these stars. However, the increase in the disk mass accretion rate is smaller than theoretically predicted, although in good agreement with previous determinations.

79 ASTRONOMY AND ASTROPHYSICS↗

Milky Way Tomography IV: Dissecting Dust

We use SDSS photometry of 73 million stars to simultaneously obtain best-fit main-sequence stellar energy distribution (SED) and amount of dust extinction along the line of sight towards each star. Using a subsample of 23 million stars with 2MASS photometry, whose addition enables more robust results, we show that SDSS photometry alone is sufficient to break degeneracies between intrinsic stellar color and dust amount when the shape of extinction curve is fixed. When using both SDSS and 2MASS photometry, the ratio of the total to selective absorption, R{sub V}, can be determined with an uncertainty of about 0.1 for most stars in high-extinction regions. These fits enable detailed studies of the dust properties and its spatial distribution, and of the stellar spatial distribution at low Galactic latitudes (|b| < 30{sup o}). Our results are in good agreement with the extinction normalization given by the Schlegel et al. (1998, SFD) dust maps at high northern Galactic latitudes, but indicate that the SFD extinction map appears to be consistently overestimated by about 20% in the southern sky, in agreement with recent study by Schlafly et al. (2010). The constraints on the shape of the dust extinction curve across the SDSS and 2MASS bandpasses disfavor the reddening law of O'Donnell (1994), but support the models by Fitzpatrick (1999) and Cardelli et al. (1989). For the latter, we find a ratio of the total to selective absorption to be R{sub V} = 3.0 {+-} 0.1(random) {+-} 0.1 (systematic) over most of the high-latitude sky. At low Galactic latitudes (|b| < 5{sup o}), we demonstrate that the SFD map cannot be reliably used to correct for extinction because most stars are embedded in dust, rather than behind it, as is the case at high Galactic latitudes. We analyze three-dimensional maps of the best-fit R{sub V} and find that R{sub V} = 3.1 cannot be ruled out in any of the ten SEGUE stripes at a precision level of {approx} 0.1 - 0.2. Our best estimate for the intrinsic scatter of R{sub V} in the regions probed by SEGUE stripes is {approx} 0.2. We introduce a method for efficient selection of candidate red giant stars in the disk, dubbed 'dusty parallax relation', which utilizes a correlation between distance and the extinction along the line of sight. We make these best-fit parameters, as well as all the input SDSS and 2MASS data, publicly available in a user-friendly format. These data can be used for studies of stellar number density distribution, the distribution of dust properties, for selecting sources whose SED differs from SEDs for high-latitude main sequence stars, and for estimating distances to dust clouds and, in turn, to molecular gas clouds.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Fosmid library end sequencing reveals a rarely known genome structure of marine shrimp Penaeus monodon

The black tiger shrimp (Penaeus monodon) is one of the most important aquaculture species in the world, representing the crustacean lineage which possesses the greatest species diversity among marine invertebrates. Yet, we barely know anything about their genomic structure. To understand the organization and evolution of the P. monodon genome, a fosmid library consisting of 288,000 colonies and was constructed, equivalent to 5.3-fold coverage of the 2.17 Gb genome. Approximately 11.1 Mb of fosmid end sequences (FESs) from 20,926 non-redundant reads representing 0.45% of the P. monodon genome were obtained for repetitive and protein-coding sequence analyses. We found that microsatellite sequences were highly abundant in the P. monodon genome, comprising 8.3% of the total length. The density and the average length of microsatellites were evidently higher in comparison to those of other taxa. AT-rich microsatellite motifs, especially poly (AT) and poly (AAT), were the most abundant. High abundance of microsatellite sequences were also found in the transcribed regions. Furthermore, via self-BlastN analysis we identified 103 novel repetitive element families which were categorized into four groups, i.e., 33 WSSV-like repeats, 14 retrotransposons, 5 gene-like repeats, and 51 unannotated repeats. Overall, various types of repeats comprise 51.18% of the P. monodon genome in length. Approximately 7.4% of the FESs contained protein-coding sequences, and the Inhibitor of Apoptosis Protein (IAP) gene and the Innexin 3 gene homologues appear to be present in high abundance in the P. monodon genome. The redundancy of various repeat types in the P. monodon genome illustrates its highly repetitive nature. In particular, long and dense microsatellite sequences as well as abundant WSSV-like sequences highlight the uniqueness of genome organization of penaeid shrimp from those of other taxa. These results provide substantial improvement to our current knowledge not only for shrimp but also for marine crustaceans of large genome size.

59 BASIC BIOLOGICAL SCIENCES↗

Multidimensional metrics for estimating phage abundance, distribution, gene density, and sequence coverage in metagenomes

Phages are the most abundant biological entities on Earth and play major ecological roles, yet the current sequenced phage genomes do not adequately represent their diversity, and little is known about the abundance and distribution of these sequenced genomes in nature. Although the study of phage ecology has benefited tremendously from the emergence of metagenomic sequencing, a systematic survey of phage genes and genomes in various ecosystems is still lacking, and fundamental questions about phage biology, lifestyle, and ecology remain unanswered. To address these questions and improve comparative analysis of phages in different metagenomes, we screened a core set of publicly available metagenomic samples for sequences related to completely sequenced phages using the web tool, Phage Eco-Locator. We then adopted and deployed an array of mathematical and statistical metrics for a multidimensional estimation of the abundance and distribution of phage genes and genomes in various ecosystems. Experiments using those metrics individually showed their usefulness in emphasizing the pervasive, yet uneven, distribution of known phage sequences in environmental metagenomes. Using these metrics in combination allowed us to resolve phage genomes into clusters that correlated with their genotypes and taxonomic classes as well as their ecological properties. By adding this set of metrics to current metaviromic analysis pipelines, where they can provide insight regarding phage mosaicism, habitat specificity, and evolution.

59 BASIC BIOLOGICAL SCIENCES↗

Structural insights into the thermal decomposition sequence of barium tetrahydrogenorthotellurate(VI), Ba[H{sub 4}TeO{sub 6}]

The compounds Ba[H{sub 4}TeO{sub 6}] (I), Ba[H{sub 2}TeO{sub 5}] (II), Ba[Te{sub 2}O{sub 6}(OH){sub 2}] (III) and Ba[TeO{sub 4}] (IV) were prepared by application of a diffusion method (I), under hydrothermal conditions (II and III) and from solid state reactions (IV), respectively. Structure analysis on the basis of single crystal X-ray diffraction data revealed novel structure types for (I), (II) and (III) and isotypism of (IV) with PrSbO{sub 4} and LaSbO{sub 4}. Common feature of the four oxotellurate(VI) structures are [TeO{sub 6}] octahedra. Whereas in the crystal structure of (I) the octahedral units are isolated, they are condensed into chains via corner-sharing in (II) and via edge-sharing in (III) and (IV). The coordination numbers of the barium cations in the four structures range from seven to ten. Although hydrogen atom positions could not be located for the structures of (I) and (II), short interpolyhedral O···O contacts are evident for strong hydrogen bonding. The temperature behaviour of (I), (II) and (IV) was monitored by simultaneous thermal analysis (STA) measurements and in situ powder X-ray diffraction, revealing the decomposition sequence Ba[H{sub 4}TeO{sub 6}] → Ba[H{sub 2}TeO{sub 5}] → Ba[TeO{sub 4}]→ Ba[TeO{sub 3}] upon heating to temperatures up to 900 °C. - Graphical abstract: The crystal structures of the four oxotellurates(VI) were determined from single crystal data. The thermal decomposition of Ba[H{sub 4}TeO{sub 6}], monitored by temperature-dependent X-ray powder diffraction and simultaneous thermal analysis measurements, involves two condensation reactions according to Ba[H{sub 4}TeO{sub 6}]→Ba[H{sub 2}TeO{sub 5}]+H{sub 2}O(↑)→Ba[TeO{sub 4}]+ H{sub 2}O(↑). Display Omitted.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Systematic Analysis of Vulnerability to Intrusion

SAVI models physical protection systems for fixed sites using Adversary Sequence Diagrams (ASDs) and then uses the ASD to find most-vulnerable adversary paths through the ASD. The identified paths have the lowest Probability of Interruption among all the paths through the ASD.

Bouchard, Ann MArie↗

A Retrospective Analysis Reveals That the 2021 Outbreaks of African Swine Fever Virus in Ghana Were Caused by Two Distinct Genotypes

African swine fever virus (ASFV) is the causative agent of African swine fever (ASF), a highly infectious and lethal disease of domesticated swine. Outbreaks of ASF have been mostly restricted to the continent of Africa. The outbreaks that have occurred outside of Africa were controlled by extensive depopulation of the domesticated pig population. However, in 2007, an outbreak occurred in the country of Georgia, where ASFV infected wild pigs and quickly spread across eastern Europe. Since the reintroduction of ASF into Europe, variants of the current pandemic strain, ASFV Georgia 2007/01 (ASFV-G), which is classified as Genotype 2 based on p72 sequencing, have been reported in countries within western Europe, Asia, and the island of Hispaniola. Additionally, isolates collected in 2020 confirmed the presence of variants of ASFV-G in Nigeria. Recently, we reported similar variants of ASFV-G collected from domestic pigs suspected of dying of ASF in Ghana in 2022. Here, we retroactively report, based on full-length sequencing, that similar variants were present in Ghana in 2021. The SNP analysis revealed derivatives of ASFV with distinct genetic markers. Furthermore, we identified three full-length ASFV genomes as Genotype 1, indicating that there were two genotypes circulating in proximity during the 2021 ASF outbreaks in Ghana.

Virology↗

High-throughput single-cell transcriptomics of bacteria using combinatorial barcoding

Microbial split-pool ligation transcriptomics (microSPLiT) is a high-throughput single-cell RNA sequencing method for bacteria. With four combinatorial barcoding rounds, microSPLiT can profile transcriptional states in hundreds of thousands of Gram-negative and Gram-positive bacteria in a single experiment without specialized equipment. As bacterial samples are fixed and permeabilized before barcoding, they can be collected and stored ahead of time. During the first barcoding round, the fixed and permeabilized bacteria are distributed into a 96-well plate, where their transcripts are reverse transcribed into cDNA and labeled with the first well-specific barcode inside the cells. The cells are mixed and redistributed two more times into new 96-well plates, where the second and third barcodes are appended to the cDNA via in-cell ligation reactions. Finally, the cells are mixed and divided into aliquot sub-libraries, which can be stored until future use or prepared for sequencing with the addition of a fourth barcode. It takes 4 days to generate sequencing-ready libraries, including 1 day for collection and overnight fixation of samples. Here, the standard plate setup enables single-cell transcriptional profiling of up to 1 million bacterial cells and up to 96 samples in a single barcoding experiment, with the possibility of expansion by adding barcoding rounds. The protocol requires experience in basic molecular biology techniques, handling of bacterial samples and preparation of DNA libraries for next-generation sequencing. It can be performed by experienced undergraduate or graduate students. Data analysis requires access to computing resources, familiarity with Unix command line and basic experience with Python or R.

59 BASIC BIOLOGICAL SCIENCES↗