Resolving the complexity of the human genome using single-molecule sequencing

The human genome is arguably the most complete mammalian reference assembly, yet more than 160 euchromatic gaps remain and aspects of its structural variation remain poorly understood ten years after its completion. To identify missing sequence and genetic variation, here we sequence and analyse a haploid human genome (CHM1) using single-molecule, real-time DNA sequencing. We close or extend 55% of the remaining interstitial gaps in the human GRCh37 reference genome - 78% of which carried long runs of degenerate short tandem repeats, often several kilobases in length, embedded within (G+C)-rich genomic regions. We resolve the complete sequence of 26,079 euchromatic structural variants at the base-pair level, including inversions, complex insertions and long tracts of tandem repeats. Most have not been previously reported, with the greatest increases in sensitivity occurring for events less than 5 kilobases in size. Compared to the human reference, we find a significant insertional bias (3:1) in regions corresponding to complex insertions and long short tandem repeats. Our results suggest a greater complexity of the human genome in the form of variation of longer and more complex repetitive DNA that can now be largely resolved with the application of this longer-read sequencing technology

Human Vav1 Expression in Hematopoietic and Cancer Cell Lines Is Regulated by c-Myb and by CpG Methylation

Vav1 is a signal transducer protein that functions as a guanine nucleotide exchange factor for the Rho/Rac GTPases in the hematopoietic system where it is exclusively expressed. Recently, Vav1 was shown to be involved in several human malignancies including neuroblastoma, lung cancer, and pancreatic ductal adenocarcinoma (PDA). Although some factors that affect vav1 expression are known, neither the physiological nor pathological regulation of vav1 expression is completely understood. We demonstrate herein that mutations in putative transcription factor binding sites at the vav1 promoter affect its transcription in cells of different histological origin. Among these sites is a consensus site for c-Myb, a hematopoietic-specific transcription factor that is also found in Vav1-expressing lung cancer cell lines. Depletion of c-Myb using siRNA led to a dramatic reduction in vav1 expression in these cells. Consistent with this, co-transfection of c-Myb activated transcription of a vav1 promoter-luciferase reporter gene construct in lung cancer cells devoid of Vav1 expression. Together, these results indicate that c-Myb is involved in vav1 expression in lung cancer cells. We also explored the methylation status of the vav1 promoter. Bisulfite sequencing revealed that the vav1 promoter was completely unmethylated in human lymphocytes, but methylated to various degrees in tissues that do not normally express vav1. The vav1 promoter does not contain CpG islands in proximity to the transcription start site; however, we demonstrated that methylation of a CpG dinucleotide at a consensus Sp1 binding site in the vav1 promoter interferes with protein binding in vitro. Our data identify two regulatory mechanisms for vav1 expression: binding of c-Myb and CpG methylation of 5′ regulatory sequences. Mutation of other putative transcription factor binding sites suggests that additional factors regulate vav1 expression as well

The Influence of cis-Regulatory Elements on DNA Methylation Fidelity

It is now established that, as compared to normal cells, the cancer cell genome has an overall inverse distribution of DNA methylation (“methylome”), i.e., predominant hypomethylation and localized hypermethylation, within “CpG islands” (CGIs). Moreover, although cancer cells have reduced methylation “fidelity” and genomic instability, accurate maintenance of aberrant methylomes that underlie malignant phenotypes remains necessary. However, the mechanism(s) of cancer methylome maintenance remains largely unknown. Here, we assessed CGI methylation patterns propagated over 1, 3, and 5 divisions of A2780 ovarian cancer cells, concurrent with exposure to the DNA cross-linking chemotherapeutic cisplatin, and observed cell generation-successive increases in total hyper- and hypo-methylated CGIs. Empirical Bayesian modeling revealed five distinct modes of methylation propagation: (1) heritable (i.e., unchanged) high- methylation (1186 probe loci in CGI microarray); (2) heritable (i.e., unchanged) low-methylation (286 loci); (3) stochastic hypermethylation (i.e., progressively increased, 243 loci); (4) stochastic hypomethylation (i.e., progressively decreased, 247 loci); and (5) considerable “random” methylation (582 loci). These results support a “stochastic model” of DNA methylation equilibrium deriving from the efficiency of two distinct processes, methylation maintenance and de novo methylation. A role for cis-regulatory elements in methylation fidelity was also demonstrated by highly significant (p<2.2×10−5) enrichment of transcription factor binding sites in CGI probe loci showing heritably high (118 elements) and low (47 elements) methylation, and also in loci demonstrating stochastic hyper-(30 elements) and hypo-(31 elements) methylation. Notably, loci having “random” methylation heritability displayed nearly no enrichment. These results demonstrate an influence of cis-regulatory elements on the nonrandom propagation of both strictly heritable and stochastically heritable CGIs

Species differential regulation of COX2 can be described by an NFκB-dependent logic AND gate

Cyclooxygenase 2 (COX2), a key regulatory enzyme of the prostaglandin/eicosanoid pathway, is an important target for anti-inflammatory therapy. It is highly induced by pro-inflammatory cytokines in a Nuclear factor kappa B (NFκB)-dependent manner. However, the mechanisms determining the amplitude and dynamics of this important pro-inflammatory event are poorly understood. Furthermore, there is significant difference between human and mouse COX2 expression in response to the inflammatory stimulus tumor necrosis factor alpha (TNFα). Here, we report the presence of a molecular logic AND gate composed of two NFκB response elements (NREs) which controls the expression of human COX2 in a switch-like manner. Combining quantitative kinetic modeling and thermostatistical analysis followed by experimental validation in iterative cycles, we show that the human COX2 expression machinery regulated by NFκB displays features of a logic AND gate. We propose that this provides a digital, noise-filtering mechanism for a tighter control of expression in response to TNFα, such that a threshold level of NFκB activation is required before the promoter becomes active and initiates transcription. This NFκB-regulated AND gate is absent in the mouse COX2 promoter, most likely contributing to its differential graded response in promoter activity and protein expression to TNFα. Our data suggest that the NFκB-regulated AND gate acts as a novel mechanism for controlling the expression of human COX2 to TNFα, and its absence in the mouse COX2 provides the foundation for further studies on understanding species-specific differential gene regulation

Antidepressant augmentation using the N-methyl-D-aspartate antagonist memantine: a randomized, double-blind, placebo-controlled trial

OBJECTIVE: Intravenous N-methyl-d-aspartate (NMDA) antagonists have shown promising results in rapidly ameliorating depression symptoms, but placebo-controlled trials of oral NMDA antagonists as monotherapy have not observed efficacy. We conducted a randomized, double-blind, placebo-controlled trial of the NMDA antagonist memantine as an augmentation treatment for patients with DSM-IV major depressive disorder. METHOD: Adult outpatients with major depressive disorder and partial response or nonresponse to their current antidepressant (as indicated by a 17-item Hamilton Depression Rating Scale score of \u3e/= 16 at baseline) were randomized (from July 2006-December 2011) to add memantine (flexible dose 5-20 mg/d, with all memantine group participants reaching the dose of 20 mg/d) (n = 15) or placebo (n = 16) to their existing treatment for 8 weeks. The primary outcome, change in Montgomery-Asberg Depression Rating Score (MADRS), was evaluated with repeated-measures mixed effects models using last-observation-carried-forward methods. Secondary outcomes included other depression and anxiety rating scales, suicidal and delusional ideation, and other adverse effects. RESULTS: 84% of participants completed the trial, including 93% of participants receiving memantine. Participants receiving memantine did not show a statistically or clinically significant change in MADRS scores compared to placebo, either over the entire study (beta = 0.133, favoring placebo, P = .74) or at study completion (week 8 mean [SD] MADRS score change = -7.13 [6.61] [memantine]; -7.25 [11.14] [placebo]; P = .97). A minimal to small effect size (comparing change to baseline variability) favoring placebo was observed (Cohen d = 0.19). Similarly, no substantial effect sizes favoring memantine nor statistically significant between-group differences were observed on secondary efficacy outcomes. CONCLUSIONS: This trial did not detect significant statistical or effect size differences between memantine and placebo augmentation among nonresponders or poor responders to conventional antidepressants. While the small number of participants is a limitation, this study suggests memantine lacks substantial efficacy as an augmentation treatment for major depressive disorder. TRIAL REGISTRATION: ClinicalTrials.gov identifier: NCT00344682

Highly accurate long-read HiFi sequencing data for five complex genomes

The PacBio(®) HiFi sequencing method yields highly accurate long-read sequencing datasets with read lengths averaging 10-25 kb and accuracies greater than 99.5%. These accurate long reads can be used to improve results for complex applications such as single nucleotide and structural variant detection, genome assembly, assembly of difficult polyploid or highly repetitive genomes, and assembly of metagenomes. Currently, there is a need for sample data sets to both evaluate the benefits of these long accurate reads as well as for development of bioinformatic tools including genome assemblers, variant callers, and haplotyping algorithms. We present deep coverage HiFi datasets for five complex samples including the two inbred model genomes Mus musculus and Zea mays, as well as two complex genomes, octoploid Fragaria × ananassa and the diploid anuran Rana muscosa. Additionally, we release sequence data from a mock metagenome community. The datasets reported here can be used without restriction to develop new algorithms and explore complex genome structure and evolution. Data were generated on the PacBio Sequel II System

Highly accurate long-read HiFi sequencing data for five complex genomes.

The PacBio® HiFi sequencing method yields highly accurate long-read sequencing datasets with read lengths averaging 10-25 kb and accuracies greater than 99.5%. These accurate long reads can be used to improve results for complex applications such as single nucleotide and structural variant detection, genome assembly, assembly of difficult polyploid or highly repetitive genomes, and assembly of metagenomes. Currently, there is a need for sample data sets to both evaluate the benefits of these long accurate reads as well as for development of bioinformatic tools including genome assemblers, variant callers, and haplotyping algorithms. We present deep coverage HiFi datasets for five complex samples including the two inbred model genomes Mus musculus and Zea mays, as well as two complex genomes, octoploid Fragaria × ananassa and the diploid anuran Rana muscosa. Additionally, we release sequence data from a mock metagenome community. The datasets reported here can be used without restriction to develop new algorithms and explore complex genome structure and evolution. Data were generated on the PacBio Sequel II System

The transcriptional diversity of 25 Drosophila cell lines

Drosophila melanogaster cell lines are important resources for cell biologists. Here, we catalog the expression of exons, genes, and unannotated transcriptional signals for 25 lines. Unannotated transcription is substantial (typically 19% of euchromatic signal). Conservatively, we identify 1405 novel transcribed regions; 684 of these appear to be new exons of neighboring, often distant, genes. Sixty-four percent of genes are expressed detectably in at least one line, but only 21% are detected in all lines. Each cell line expresses, on average, 5885 genes, including a common set of 3109. Expression levels vary over several orders of magnitude. Major signaling pathways are well represented: most differentiation pathways are “off” and survival/growth pathways “on.” Roughly 50% of the genes expressed by each line are not part of the common set, and these show considerable individuality. Thirty-one percent are expressed at a higher level in at least one cell line than in any single developmental stage, suggesting that each line is enriched for genes characteristic of small sets of cells. Most remarkable is that imaginal disc-derived lines can generally be assigned, on the basis of expression, to small territories within developing discs. These mappings reveal unexpected stability of even fine-grained spatial determination. No two cell lines show identical transcription factor expression. We conclude that each line has retained features of an individual founder cell superimposed on a common “cell line“ gene expression pattern.Organismic and Evolutionary Biolog

Genome-wide analysis of promoter architecture in Drosophila melanogaster

Core promoters are critical regions for gene regulation in higher eukaryotes. However, the boundaries of promoter regions, the relative rates of initiation at the transcription start sites (TSSs) distributed within them, and the functional significance of promoter architecture remain poorly understood. We produced a high-resolution map of promoters active in the Drosophila melanogaster embryo by integrating data from three independent and complementary methods: 21 million cap analysis of gene expression (CAGE) tags, 1.2 million RNA ligase mediated rapid amplification of cDNA ends (RLM-RACE) reads, and 50,000 cap-trapped expressed sequence tags (ESTs). We defined 12,454 promoters of 8037 genes. Our analysis indicates that, due to non-promoter-associated RNA background signal, previous studies have likely overestimated the number of promoter-associated CAGE clusters by fivefold. We show that TSS distributions form a complex continuum of shapes, and that promoters active in the embryo and adult have highly similar shapes in 95% of cases. This suggests that these distributions are generally determined by static elements such as local DNA sequence and are not modulated by dynamic signals such as histone modifications. Transcription factor binding motifs are differentially enriched as a function of promoter shape, and peaked promoter shape is correlated with both temporal and spatial regulation of gene expression. Our results contribute to the emerging view that core promoters are functionally diverse and control patterning of gene expression in Drosophila and mammals