Semi-automated assembly of high-quality diploid human reference genomes.

The current human reference genome, GRCh38, represents over 20 years of effort to generate a high-quality assembly, which has benefitted society1,2. However, it still has many gaps and errors, and does not represent a biological genome as it is a blend of multiple individuals3,4. Recently, a high-quality telomere-to-telomere reference, CHM13, was generated with the latest long-read technologies, but it was derived from a hydatidiform mole cell line with a nearly homozygous genome5. To address these limitations, the Human Pangenome Reference Consortium formed with the goal of creating high-quality, cost-effective, diploid genome assemblies for a pangenome reference that represents human genetic diversity6. Here, in our first scientific report, we determined which combination of current genome sequencing and assembly approaches yield the most complete and accurate diploid genome assembly with minimal manual curation. Approaches that used highly accurate long reads and parent-child data with graph-based haplotype phasing during assembly outperformed those that did not. Developing a combination of the top-performing methods, we generated our first high-quality diploid reference assembly, containing only approximately four gaps per chromosome on average, with most chromosomes within ±1% of the length of CHM13. Nearly 48% of protein-coding genes have non-synonymous amino acid changes between haplotypes, and centromeric regions showed the highest diversity. Our findings serve as a foundation for assembling near-complete diploid human genomes at scale for a pangenome reference to capture global genetic variation from single nucleotides to structural rearrangements.

Long-read genome sequencing identifies causal structural variation in a Mendelian disease.

PurposeCurrent clinical genomics assays primarily utilize short-read sequencing (SRS), but SRS has limited ability to evaluate repetitive regions and structural variants. Long-read sequencing (LRS) has complementary strengths, and we aimed to determine whether LRS could offer a means to identify overlooked genetic variation in patients undiagnosed by SRS.MethodsWe performed low-coverage genome LRS to identify structural variants in a patient who presented with multiple neoplasia and cardiac myxomata, in whom the results of targeted clinical testing and genome SRS were negative.ResultsThis LRS approach yielded 6,971 deletions and 6,821 insertions > 50 bp. Filtering for variants that are absent in an unrelated control and overlap a disease gene coding exon identified three deletions and three insertions. One of these, a heterozygous 2,184 bp deletion, overlaps the first coding exon of PRKAR1A, which is implicated in autosomal dominant Carney complex. RNA sequencing demonstrated decreased PRKAR1A expression. The deletion was classified as pathogenic based on guidelines for interpretation of sequence variants.ConclusionThis first successful application of genome LRS to identify a pathogenic variant in a patient suggests that LRS has significant potential for the identification of disease-causing structural variation. Larger studies will ultimately be required to evaluate the potential clinical utility of LRS.

Genome and Methylome Variation in Helicobacter pylori With a cag Pathogenicity Island During Early Stages of Human Infection.

BACKGROUND & AIMS: Helicobacter pylori is remarkable for its genetic variation; yet, little is known about its genetic changes during early stages of human infection, as the bacteria adapt to their new environment. We analyzed genome and methylome variations in a fully virulent strain of H pylori during experimental infection. METHODS: We performed a randomized Phase I/II, observer-blind, placebo-controlled study of 12 healthy, H pylori-negative adults in Germany from October 2008 through March 2010. The volunteers were given a prophylactic vaccine candidate (n = 7) or placebo (n = 5) and then challenged with H pylori strain BCM-300. Biopsy samples were collected and H pylori were isolated. Genomes of the challenge strain and 12 reisolates, obtained 12 weeks after (or in 1 case, 62 weeks after) infection were sequenced by single-molecule, real-time technology, which, in parallel, permitted determination of genome-wide methylation patterns for all strains. Functional effects of genetic changes observed in H pylori strains during human infection were assessed by measuring release of interleukin 8 from AGS cells (to detect cag pathogenicity island function), neutral red uptake (to detect vacuolating cytotoxin activity), and adhesion assays. RESULTS: The observed mutation rate was in agreement with rates previously determined from patients with chronic H pylori infections, without evidence of a mutation burst. A loss of cag pathogenicity island function was observed in 3 reisolates. In addition, 3 reisolates from the vaccine group acquired mutations in the vacuolating cytotoxin gene vacA, resulting in loss of vacuolization activity. We observed interstrain variation in methylomes due to phase variation in genes encoding methyltransferases. CONCLUSIONS: We analyzed adaptation of a fully virulent strain of H pylori to 12 different volunteers to obtain a robust estimate of the frequency of genetic and epigenetic changes in the absence of interstrain recombination. Our findings indicate that the large amount of genetic variation in H pylori poses a challenge to vaccine development. ClinicalTrials.gov no: NCT00736476.

Methylomic and phenotypic analysis of the ModH5 phasevarion of Helicobacter pylori.

The Helicobacter pylori phase variable gene modH, typified by gene HP1522 in strain 26695, encodes a N6-adenosine type III DNA methyltransferase. Our previous studies identified multiple strain-specific modH variants (modH1 - modH19) and showed that phase variation of modH5 in H. pylori P12 influenced expression of motility-associated genes and outer membrane protein gene hopG. However, the ModH5 DNA recognition motif and the mechanism by which ModH5 controls gene expression were unknown. Here, using comparative single molecule real-time sequencing, we identify the DNA site methylated by ModH5 as 5'-Gm6ACC-3'. This motif is vastly underrepresented in H. pylori genomes, but overrepresented in a number of virulence genes, including motility-associated genes, and outer membrane protein genes. Motility and the number of flagella of H. pylori P12 wild-type were significantly higher than that of isogenic modH5 OFF or ΔmodH5 mutants, indicating that phase variable switching of modH5 expression plays a role in regulating H. pylori motility phenotypes. Using the flagellin A (flaA) gene as a model, we show that ModH5 modulates flaA promoter activity in a GACC methylation-dependent manner. These findings provide novel insights into the role of ModH5 in gene regulation and how it mediates epigenetic regulation of H. pylori motility.

Identification of a Pseudomonas aeruginosa PAO1 DNA Methyltransferase, Its Targets, and Physiological Roles.

DNA methylation is widespread among prokaryotes, and most DNA methylation reactions are catalyzed by adenine DNA methyltransferases, which are part of restriction-modification (R-M) systems. R-M systems are known for their role in the defense against foreign DNA; however, DNA methyltransferases also play functional roles in gene regulation. In this study, we used single-molecule real-time (SMRT) sequencing to uncover the genome-wide DNA methylation pattern in the opportunistic pathogen Pseudomonas aeruginosa PAO1. We identified a conserved sequence motif targeted by an adenine methyltransferase of a type I R-M system and quantified the presence of N6-methyladenine using liquid chromatography-tandem mass spectrometry (LC-MS/MS). Changes in the PAO1 methylation status were dependent on growth conditions and affected P. aeruginosa pathogenicity in a Galleria mellonella infection model. Furthermore, we found that methylated motifs in promoter regions led to shifts in sense and antisense gene expression, emphasizing the role of enzymatic DNA methylation as an epigenetic control of phenotypic traits in P. aeruginosa Since the DNA methylation enzymes are not encoded in the core genome, our findings illustrate how the acquisition of accessory genes can shape the global P. aeruginosa transcriptome and thus may facilitate adaptation to new and challenging habitats.IMPORTANCE With the introduction of advanced technologies, epigenetic regulation by DNA methyltransferases in bacteria has become a subject of intense studies. Here we identified an adenosine DNA methyltransferase in the opportunistic pathogen Pseudomonas aeruginosa PAO1, which is responsible for DNA methylation of a conserved sequence motif. The methylation level of all target sequences throughout the PAO1 genome was approximated to be in the range of 65 to 85% and was dependent on growth conditions. Inactivation of the methyltransferase revealed an attenuated-virulence phenotype in the Galleria mellonella infection model. Furthermore, differential expression of more than 90 genes was detected, including the small regulatory RNA prrF1, which contributes to a global iron-sparing response via the repression of a set of gene targets. Our finding of a methylation-dependent repression of the antisense transcript of the prrF1 small regulatory RNA significantly expands our understanding of the regulatory mechanisms underlying active DNA methylation in bacteria.

Enzymatic Degradation of Phenazines Can Generate Energy and Protect Sensitive Organisms from Toxicity.

UNLABELLED: Diverse bacteria, including several Pseudomonas species, produce a class of redox-active metabolites called phenazines that impact different cell types in nature and disease. Phenazines can affect microbial communities in both positive and negative ways, where their presence is correlated with decreased species richness and diversity. However, little is known about how the concentration of phenazines is modulated in situ and what this may mean for the fitness of members of the community. Through culturing of phenazine-degrading mycobacteria, genome sequencing, comparative genomics, and molecular analysis, we identified several conserved genes that are important for the degradation of three Pseudomonas-derived phenazines: phenazine-1-carboxylic acid (PCA), phenazine-1-carboxamide (PCN), and pyocyanin (PYO). PCA can be used as the sole carbon source for growth by these organisms. Deletion of several genes in Mycobacterium fortuitum abolishes the degradation phenotype, and expression of two genes in a heterologous host confers the ability to degrade PCN and PYO. In cocultures with phenazine producers, phenazine degraders alter the abundance of different phenazine types. Not only does degradation support mycobacterial catabolism, but also it provides protection to bacteria that would otherwise be inhibited by the toxicity of PYO. Collectively, these results serve as a reminder that microbial metabolites can be actively modified and degraded and that these turnover processes must be considered when the fate and impact of such compounds in any environment are being assessed. IMPORTANCE: Phenazine production by Pseudomonas spp. can shape microbial communities in a variety of environments ranging from the cystic fibrosis lung to the rhizosphere of dryland crops. For example, in the rhizosphere, phenazines can protect plants from infection by pathogenic fungi. The redox activity of phenazines underpins their antibiotic activity, as well as providing pseudomonads with important physiological benefits. Our discovery that soil mycobacteria can catabolize phenazines and thereby protect other organisms against phenazine toxicity suggests that phenazine degradation may influence turnover in situ. The identification of genes involved in the degradation of phenazines opens the door to monitoring turnover in diverse environments, an essential process to consider when one is attempting to understand or control communities influenced by phenazines.

Characterization of fusion genes and the significantly expressed fusion isoforms in breast cancer by hybrid sequencing.

We developed an innovative hybrid sequencing approach, IDP-fusion, to detect fusion genes, determine fusion sites and identify and quantify fusion isoforms. IDP-fusion is the first method to study gene fusion events by integrating Third Generation Sequencing long reads and Second Generation Sequencing short reads. We applied IDP-fusion to PacBio data and Illumina data from the MCF-7 breast cancer cells. Compared with the existing tools, IDP-fusion detects fusion genes at higher precision and a very low false positive rate. The results show that IDP-fusion will be useful for unraveling the complexity of multiple fusion splices and fusion isoforms within tumorigenesis-relevant fusion genes.

Single-molecule fluorescence imaging of processive myosin with enhanced background suppression using linear zero-mode waveguides (ZMWs) and convex lens induced confinement (CLIC).

Resolving single fluorescent molecules in the presence of high fluorophore concentrations remains a challenge in single-molecule biophysics that limits our understanding of weak molecular interactions. Total internal reflection fluorescence (TIRF) imaging, the workhorse of single-molecule fluorescence microscopy, enables experiments at concentrations up to about 100 nM, but many biological interactions have considerably weaker affinities, and thus require at least one species to be at micromolar or higher concentration. Current alternatives to TIRF often require three-dimensional confinement, and thus can be problematic for extended substrates, such as cytoskeletal filaments. To address this challenge, we have demonstrated and applied two new single-molecule fluorescence microscopy techniques, linear zero-mode waveguides (ZMWs) and convex lens induced confinement (CLIC), for imaging the processive motion of molecular motors myosin V and VI along actin filaments. Both technologies will allow imaging in the presence of higher fluorophore concentrations than TIRF microscopy. They will enable new biophysical measurements of a wide range of processive molecular motors that move along filamentous tracks, such as other myosins, dynein, and kinesin. A particularly salient application of these technologies will be to examine chemomechanical coupling by directly imaging fluorescent nucleotide molecules interacting with processive motors as they traverse their actin or microtubule tracks.

Enhanced 5-methylcytosine detection in single-molecule, real-time sequencing via Tet1 oxidation.

BACKGROUND: DNA methylation serves as an important epigenetic mark in both eukaryotic and prokaryotic organisms. In eukaryotes, the most common epigenetic mark is 5-methylcytosine, whereas prokaryotes can have 6-methyladenine, 4-methylcytosine, or 5-methylcytosine. Single-molecule, real-time sequencing is capable of directly detecting all three types of modified bases. However, the kinetic signature of 5-methylcytosine is subtle, which presents a challenge for detection. We investigated whether conversion of 5-methylcytosine to 5-carboxylcytosine using the enzyme Tet1 would enhance the kinetic signature, thereby improving detection. RESULTS: We characterized the kinetic signatures of various cytosine modifications, demonstrating that 5-carboxylcytosine has a larger impact on the local polymerase rate than 5-methylcytosine. Using Tet1-mediated conversion, we show improved detection of 5-methylcytosine using in vitro methylated templates and apply the method to the characterization of 5-methylcytosine sites in the genomes of Escherichia coli MG1655 and Bacillus halodurans C-125. CONCLUSIONS: We have developed a method for the enhancement of directly detecting 5-methylcytosine during single-molecule, real-time sequencing. Using Tet1 to convert 5-methylcytosine to 5-carboxylcytosine improves the detection rate of this important epigenetic marker, thereby complementing the set of readily detectable microbial base modifications, and enhancing the ability to interrogate eukaryotic epigenetic markers.

The birth of the Epitranscriptome: deciphering the function of RNA modifications.

Recent studies have found methyl-6-adenosine in thousands of mammalian genes, and this modification is most pronounced near the beginning of the 3' UTR. We present a perspective on current work and new single-molecule sequencing methods for detecting RNA base modifications.

Live-cell imaging of single receptor composition using zero-mode waveguide nanostructures.

We exploit the optical and spatial features of subwavelength nanostructures to examine individual receptors on the plasma membrane of living cells. Receptors were sequestered in portions of the membrane projected into zero-mode waveguides. Using single-step photobleaching of green fluorescent protein incorporated into individual subunits, the resulting spatial isolation was used to measure subunit stoichiometry in α4ß4 and α4ß2 nicotinic acetylcholine and P2X2 ATP receptors. We also show that nicotine and cytisine have differential effects on α4ß2 stoichiometry.

Direct detection and sequencing of damaged DNA bases.

Products of various forms of DNA damage have been implicated in a variety of important biological processes, such as aging, neurodegenerative diseases, and cancer. Therefore, there exists great interest to develop methods for interrogating damaged DNA in the context of sequencing. Here, we demonstrate that single-molecule, real-time (SMRT®) DNA sequencing can directly detect damaged DNA bases in the DNA template - as a by-product of the sequencing method - through an analysis of the DNA polymerase kinetics that are altered by the presence of a modified base. We demonstrate the sequencing of several DNA templates containing products of DNA damage, including 8-oxoguanine, 8-oxoadenine, O6-methylguanine, 1-methyladenine, O4-methylthymine, 5-hydroxycytosine, 5-hydroxyuracil, 5-hydroxymethyluracil, or thymine dimers, and show that these base modifications can be readily detected with single-modification resolution and DNA strand specificity. We characterize the distinct kinetic signatures generated by these DNA base modifications.

Selective aluminum passivation for targeted immobilization of single DNA polymerase molecules in zero-mode waveguide nanostructures.

Optical nanostructures have enabled the creation of subdiffraction detection volumes for single-molecule fluorescence microscopy. Their applicability is extended by the ability to place molecules in the confined observation volume without interfering with their biological function. Here, we demonstrate that processive DNA synthesis thousands of bases in length was carried out by individual DNA polymerase molecules immobilized in the observation volumes of zero-mode waveguides (ZMWs) in high-density arrays. Selective immobilization of polymerase to the fused silica floor of the ZMW was achieved by passivation of the metal cladding surface using polyphosphonate chemistry, producing enzyme density contrasts of glass over aluminum in excess of 400:1. Yields of single-molecule occupancies of approximately 30% were obtained for a range of ZMW diameters (70-100 nm). Results presented here support the application of immobilized single DNA polymerases in ZMW arrays for long-read-length DNA sequencing.

Conformational changes of calmodulin upon Ca2+ binding studied with a microfluidic mixer.

A microfluidic mixer is applied to study the kinetics of calmodulin conformational changes upon Ca2+ binding. The device facilitates rapid, uniform mixing by decoupling hydrodynamic focusing from diffusive mixing and accesses time scales of tens of microseconds. The mixer is used in conjunction with multiphoton microscopy to examine the fast Ca2+-induced transitions of acrylodan-labeled calmodulin. We find that the kinetic rates of the conformational changes in two homologous globular domains differ by more than an order of magnitude. The characteristic time constants are approximately 490 micros for the transitions in the C-terminal domain and approximately 20 ms for those in the N-terminal domain of the protein. We discuss possible mechanisms for the two distinct events and the biological role of the stable intermediate, half-saturated calmodulin.