The microbiome in preadolescent acne: Assessment and prospective analysis of the influence of benzoyl peroxide.

BACKGROUND/OBJECTIVES: The pathogenesis of preadolescent acne has not been well studied, and it is uncertain if Cutibacterium acnes is a predominant organism in the microbiome in this age group. The aim of this study was to analyze the microbiome of preadolescent females and to assess whether benzoyl peroxide impacts the microbiome. METHODS: The study enrolled girls, aged 7-12 years, with evidence of at least six acne lesions who had not been previously treated. Participants' skin surface of forehead, cheeks, nose, chin, left retroauricular crease, and extruded contents of a comedonal lesion were sampled at baseline. Participants used benzoyl peroxide 4% wash for 6-8 weeks and returned for skin surface sampling and extraction collection. Microbiome analysis was performed using 16S ribosomal RNA gene amplicon sequencing on all swab and lesional extraction samples. RESULTS: Fifty-one participants were enrolled with a median IGA score of 2 (mild). Changes in microbiome diversity were associated with increasing age and number of acne lesions (P = 0.001). C. acnes had higher abundances on forehead and nose, as opposed to cheeks and chin (P = 0.009). Bacterial diversity (alpha diversity) of the skin microbiome was comparable between preadolescent at baseline and after treatment with benzoyl peroxide. CONCLUSION: This is the first large assessment characterizing female acne microbiome in early and late preadolescence. Results show that preadolescent acne can vary in its microbial profile, reflecting surrounding changes associated with the onset of puberty. Although benzoyl peroxide use was associated with decreased acne counts, its effect on microbial diversity was not demonstrated in our study.

In vitro, long-range sequence information for de novo genome assembly via transposase contiguity.

We describe a method that exploits contiguity preserving transposase sequencing (CPT-seq) to facilitate the scaffolding of de novo genome assemblies. CPT-seq is an entirely in vitro means of generating libraries comprised of 9216 indexed pools, each of which contains thousands of sparsely sequenced long fragments ranging from 5 kilobases to > 1 megabase. These pools are "subhaploid," in that the lengths of fragments contained in each pool sums to ∼5% to 10% of the full genome. The scaffolding approach described here, termed fragScaff, leverages coincidences between the content of different pools as a source of contiguity information. Specifically, CPT-seq data is mapped to a de novo genome assembly, followed by the identification of pairs of contigs or scaffolds whose ends disproportionately co-occur in the same indexed pools, consistent with true adjacency in the genome. Such candidate "joins" are used to construct a graph, which is then resolved by a minimum spanning tree. As a proof-of-concept, we apply CPT-seq and fragScaff to substantially boost the contiguity of de novo assemblies of the human, mouse, and fly genomes, increasing the scaffold N50 of de novo assemblies by eight- to 57-fold with high accuracy. We also demonstrate that fragScaff is complementary to Hi-C-based contact probability maps, providing midrange contiguity to support robust, accurate chromosome-scale de novo genome assemblies without the need for laborious in vivo cloning steps. Finally, we demonstrate CPT-seq as a means of anchoring unplaced novel human contigs to the reference genome as well as for detecting misassembled sequences.

Fitness landscape of antibiotic tolerance in Pseudomonas aeruginosa biofilms.

Bacteria in biofilms have higher antibiotic tolerance than their planktonic counterparts. A major outstanding question is the degree to which the biofilm-specific cellular state and its constituent genetic determinants contribute to this hyper-tolerant phenotype. Here, we used genome-wide functional profiling of a complex, heterogeneous mutant population of Pseudomonas aeruginosa MPAO1 in biofilm and planktonic growth conditions with and without tobramycin to systematically quantify the contribution of each locus to antibiotic tolerance under these two states. We identified large sets of mutations that contribute to antibiotic tolerance predominantly in the biofilm or planktonic setting only, offering global insights into the differences and similarities between biofilm and planktonic antibiotic tolerance. Our mixed population-based experimental design recapitulated the complexity of natural biofilms and, unlike previous studies, revealed clinically observed behaviors including the emergence of quorum sensing-deficient mutants. Our study revealed a substantial contribution of the cellular state to the antibiotic tolerance of biofilms, providing a rational foundation for the development of novel therapeutics against P. aeruginosa biofilm-associated infections.

Quantitative measurement of allele-specific protein expression in a diploid yeast hybrid by LC-MS.

Understanding the genetic basis of gene regulatory variation is a key goal of evolutionary and medical genetics. Regulatory variation can act in an allele-specific manner (cis-acting) or it can affect both alleles of a gene (trans-acting). Differential allele-specific expression (ASE), in which the expression of one allele differs from another in a diploid, implies the presence of cis-acting regulatory variation. While microarrays and high-throughput sequencing have enabled genome-wide measurements of transcriptional ASE, methods for measurement of protein ASE (pASE) have lagged far behind. We describe a flexible, accurate, and scalable strategy for measurement of pASE by liquid chromatography-coupled mass spectrometry (LC-MS). We apply this approach to a hybrid between the yeast species Saccharomyces cerevisiae and Saccharomyces bayanus. Our results provide the first analysis of the relative contribution of cis-acting and trans-acting regulatory differences to protein expression divergence between yeast species.

Newly identified genetic variations in common Escherichia coli MG1655 stock cultures.

We have recently identified seven mutations in commonly used stocks of the sequenced Escherichia coli strain MG1655 which do not appear in the reference sequence. The mutations are likely to cause loss of function of the glpR and crl genes, which may have serious implications for physiological experiments using the affected strains.

Genetic dissection of an exogenously induced biofilm in laboratory and clinical isolates of E. coli.

Microbial biofilms are a dominant feature of many human infections. However, developing effective strategies for controlling biofilms requires an understanding of the underlying biology well beyond what currently exists. Using a novel strategy, we have induced formation of a robust biofilm in Escherichia coli by utilizing an exogenous source of poly-N-acetylglucosamine (PNAG) polymer, a major virulence factor of many pathogens. Through microarray profiling of competitive selections, carried out in both transposon insertion and over-expression libraries, we have revealed the genetic basis of PNAG-based biofilm formation. Our observations reveal the dominance of electrostatic interactions between PNAG and surface structures such as lipopolysaccharides. We show that regulatory modulation of these surface structures has significant impact on biofilm formation behavior of the cell. Furthermore, the majority of clinical isolates which produced PNAG also showed the capacity to respond to the exogenously produced version of the polymer.

Regulatory and metabolic rewiring during laboratory evolution of ethanol tolerance in E. coli.

Understanding the genetic basis of adaptation is a central problem in biology. However, revealing the underlying molecular mechanisms has been challenging as changes in fitness may result from perturbations to many pathways, any of which may contribute relatively little. We have developed a combined experimental/computational framework to address this problem and used it to understand the genetic basis of ethanol tolerance in Escherichia coli. We used fitness profiling to measure the consequences of single-locus perturbations in the context of ethanol exposure. A module-level computational analysis was then used to reveal the organization of the contributing loci into cellular processes and regulatory pathways (e.g. osmoregulation and cell-wall biogenesis) whose modifications significantly affect ethanol tolerance. Strikingly, we discovered that a dominant component of adaptation involves metabolic rewiring that boosts intracellular ethanol degradation and assimilation. Through phenotypic and metabolomic analysis of laboratory-evolved ethanol-tolerant strains, we investigated naturally accessible pathways of ethanol tolerance. Remarkably, these laboratory-evolved strains, by and large, follow the same adaptive paths as inferred from our coarse-grained search of the fitness landscape.

Influence of genotype and nutrition on survival and metabolism of starving yeast.

Starvation of yeast cultures limited by auxotrophic requirements results in glucose wasting and failure to achieve prompt cell-cycle arrest when compared with starvation for basic natural nutrients like phosphate or sulfate. We measured the survival of yeast auxotrophs upon starvation for different nutrients and found substantial differences: When deprived of leucine or uracil, viability declined exponentially with a half-life of <2 days, whereas when the same strains were deprived of phosphate or sulfate, the half-life was approximately 10 days. The survival rates of nongrowing auxotrophs deprived of uracil or leucine depended on the carbon source available during starvation, but were independent of the carbon source during prior growth. We performed an enrichment procedure for mutants that suppress lethality during auxotrophy starvation. We repeatedly found loss-of-function mutations in a number of functionally related genes. Mutations in PPM1, which methylates protein phosphatase 2A, and target of rapamycin (TOR1) were characterized further. Deletion of PPM1 almost completely suppressed the rapid lethality and substantially suppressed glucose wasting during starvation for leucine or uracil. Suppression by a deletion of TOR1 was less complete. We suggest that, similar to the Warburg effect observed in tumor cells, starving yeast auxotrophs wastes glucose as a consequence of the failure of conserved growth control pathways. Furthermore, we suggest that our results on condition-dependent chronological lifespan have important implications for the interpretation and design of studies on chronological aging.

Contiguity-Preserving Transposition Sequencing (CPT-Seq) for Genome-Wide Haplotyping, Assembly, and Single-Cell ATAC-Seq.

Most genomes to date have been sequenced without taking into account the diploid nature of the genome. However, the distribution of variants on each individual chromosome can (1) significantly impact gene regulation and protein function, (2) have important implications for analyses of population history and medical genetics, and (3) be of great value for accurate interpretation of medically relevant genetic variation. Here, we describe a comprehensive and detailed protocol for an ultra fast (<3 h library preparation), cost-effective, and scalable haplotyping method, named Contiguity Preserving Transposition sequencing or CPT-seq (Amini et al., Nat Genet 46(12):1343-1349, 2014). CPT-seq accurately phases >95 % of the whole human genome in Mb-scale phasing blocks. Additionally, the same workflow can be used to aid de novo assembly (Adey et al., Genome Res 24(12):2041-2049, 2014), detect structural variants, and perform single cell ATAC-seq analysis (Cusanovich et al., Science 348(6237):910-914, 2015).

A multi-enzyme model for Pyrosequencing.

Pyrosequencing is a DNA sequencing technique based on sequencing-by-synthesis enabling rapid real-time sequence determination. This technique employs four enzymatic reactions in a single tube to monitor DNA synthesis. Nucleotides are added iteratively to the reaction and in case of incorporation, pyrophosphate (PPi) is released. PPi triggers a series of reactions resulting in production of light, which is proportional to the amount of DNA and number of incorporated nucleotides. Generated light is detected and recorded by a detector system in the form of a peak signal, which reflects the activity of all four enzymes in the reaction. We have developed simulations to model the kinetics of the enzymes. These simulations provide a full model for the Pyrosequencing four-enzyme system, based on which the peak height and shape can be predicted depending on the concentrations of enzymes and substrates. Simulation results are shown to be compatible with experimental data. Based on these simulations, the rate-limiting steps in the chain can be determined, and K(M) and kcat of all four enzymes in Pyrosequencing can be calculated.

Haplotype-resolved whole-genome sequencing by contiguity-preserving transposition and combinatorial indexing.

Haplotype-resolved genome sequencing enables the accurate interpretation of medically relevant genetic variation, deep inferences regarding population history and non-invasive prediction of fetal genomes. We describe an approach for genome-wide haplotyping based on contiguity-preserving transposition (CPT-seq) and combinatorial indexing. Tn5 transposition is used to modify DNA with adaptor and index sequences while preserving contiguity. After DNA dilution and compartmentalization, the transposase is removed, resolving the DNA into individually indexed libraries. The libraries in each compartment, enriched for neighboring genomic elements, are further indexed via PCR. Combinatorial 96-plex indexing at both the transposition and PCR stage enables the construction of phased synthetic reads from each of the nearly 10,000 'virtual compartments'. We demonstrate the feasibility of this method by assembling >95% of the heterozygous variants in a human genome into long, accurate haplotype blocks (N50 = 1.4-2.3 Mb). The rapid, scalable and cost-effective workflow could enable haplotype resolution to become routine in human genome sequencing.

Mutation in ST6GALNAC5 identified in family with coronary artery disease.

We aimed to identify the genetic cause of coronary artery disease (CAD) in an Iranian pedigree. Genetic linkage analysis identified three loci with an LOD score of 2.2. Twelve sequence variations identified by exome sequencing were tested for segregation with disease. A p.Val99Met causing mutation in ST6GALNAC5 was considered the likely cause of CAD. ST6GALNAC5 encodes sialyltransferase 7e. The variation affects a highly conserved amino acid, was absent in 800 controls, and was predicted to damage protein function. ST6GALNAC5 is positioned within loci previously linked to CAD-associated parameters. While hypercholesterolemia was a prominent feature in the family, clinical and genetic data suggest that this condition is not caused by the mutation in ST6GALNAC5. Sequencing of ST6GALNAC5 in 160 Iranian patients revealed a candidate causative stop-loss mutation in two other patients. The p.Val99Met and stop-loss mutations both caused increased sialyltransferase activity. Sequence data from combined Iranian and US controls and CAD affected individuals provided evidence consistent with potential role of ST6GALNAC5 in CAD. We conclude that ST6GALNAC5 mutations can cause CAD. There is substantial literature suggesting a relation between sialyltransferase and sialic acid levels and coronary disease. Our findings provide strong evidence for the existence of this relation.

Mutation in CYP27A1 identified in family with coronary artery disease.

Coronary artery disease (CAD) is a leading cause of death worldwide. Myocardial infarction is the most severe outcome of CAD. Despite extensive efforts, the genetics of CAD is poorly understood. We aimed to identify the genetic cause of CAD in a pedigree with several affected individuals. Exome sequencing led to identification of a mutation in CYP27A1 that causes p.Arg225His in the encoded protein sterol 27-hydroxylase as the likely cause of CAD in the pedigree. The enzyme is multifunctional, and several of its functions including its functions in vitamin D metabolism and reverse cholesterol transport (RCT) are relevant to the CAD phenotype. Measurements of vitamin D levels suggested that the mutation does not affect CAD by affecting this parameter. We suggest that the mutation may cause CAD by affecting RCT. Screening of all coding regions of the CYP27A1 in 100 additional patients led to finding four variations (p.Arg14Gly, p.Arg26Lys, p.Ala27Arg, and p.Val86Met) in seven patients that may contribute to their CAD status. CYP27A1 is the known causative gene of cerebrotendinous xanthomatosis, a disorder which is sometimes accompanied by early onset atherosclerosis. This and the observation of potentially harmful variations in unrelated CAD patients provide additional evidence for the suggested causative role of the p.Arg225His mutation in CAD.

Antibiotics and the post-genome revolution.

The emergence of pathogenic bacteria resistant to multiple antimicrobial agents is turning into a major crisis in human and veterinary medicine. This necessitates a serious re-evaluation of our approaches toward antibacterial drug discovery and use. Concurrent advances in genomics including whole-genome sequencing, genotyping, and gene expression profiling have the potential to transform our basic understanding of antimicrobial pathways and lead to the discovery of novel targets and therapeutics.

Accurate proteome-wide protein quantification from high-resolution 15N mass spectra.

In quantitative mass spectrometry-based proteomics, the metabolic incorporation of a single source of 15N-labeled nitrogen has many advantages over using stable isotope-labeled amino acids. However, the lack of a robust computational framework for analyzing the resulting spectra has impeded wide use of this approach. We have addressed this challenge by introducing a new computational methodology for analyzing 15N spectra in which quantification is integrated with identification. Application of this method to an Escherichia coli growth transition reveals significant improvement in quantification accuracy over previous methods.