Pfam is a widely used database of protein families and domains. This article describes a set of major updates that we have implemented in the latest release (version 24.0). The most important change is that we now use HMMER3, the latest version of the popular profile hidden Markov model package. This software is approximately 100 times faster than HMMER2 and is more sensitive due to the routine use of the forward algorithm. The move to HMMER3 has necessitated numerous changes to Pfam that are described in detail. Pfam release 24.0 contains 11,912 families, of which a large number have been significantly updated during the past two years. Pfam is available via servers in the UK (http://pfam.sanger.ac.uk/), the USA (http://pfam.janelia.org/) and Sweden (http://pfam.sbc.su.se/).

/pdf/the-pfam-protein-families-database-3t3lguk949.pdf

The Pfam protein families database

The 1000 Genomes Project set out to provide a comprehensive description of common human genetic variation by applying whole-genome sequencing to a diverse set of individuals from multiple populations. Here we report completion of the project, having reconstructed the genomes of 2,504 individuals from 26 populations using a combination of low-coverage whole-genome sequencing, deep exome sequencing, and dense microarray genotyping. We characterized a broad spectrum of genetic variation, in total over 88 million variants (84.7 million single nucleotide polymorphisms (SNPs), 3.6 million short insertions/deletions (indels), and 60,000 structural variants), all phased onto high-quality haplotypes. This resource includes >99% of SNP variants with a frequency of >1% for a variety of ancestries. We describe the distribution of genetic variation across the global sample, and discuss the implications for common disease studies.

https://cdr.lib.unc.edu/downloads/v692t8920

A global reference for human genetic variation.

A 2.91-billion base pair (bp) consensus sequence of the euchromatic portion of the human genome was generated by the whole-genome shotgun sequencing method. The 14.8-billion bp DNA sequence was generated over 9 months from 27,271,853 high-quality sequence reads (5.11-fold coverage of the genome) from both ends of plasmid clones made from the DNA of five individuals. Two assembly strategies-a whole-genome assembly and a regional chromosome assembly-were used, each combining sequence data from Celera and the publicly funded genome effort. The public data were shredded into 550-bp segments to create a 2.9-fold coverage of those genome regions that had been sequenced, without including biases inherent in the cloning and assembly procedure used by the publicly funded group. This brought the effective coverage in the assemblies to eightfold, reducing the number and size of gaps in the final assembly over what would be obtained with 5.11-fold coverage. The two assembly strategies yielded very similar results that largely agree with independent mapping data. The assemblies effectively cover the euchromatic regions of the human chromosomes. More than 90% of the genome is in scaffold assemblies of 100,000 bp or more, and 25% of the genome is in scaffolds of 10 million bp or larger. Analysis of the genome sequence revealed 26,588 protein-encoding transcripts for which there was strong corroborating evidence and an additional approximately 12,000 computationally derived genes with mouse matches or other weak supporting evidence. Although gene-dense clusters are obvious, almost half the genes are dispersed in low G+C sequence separated by large tracts of apparently noncoding sequence. Only 1.1% of the genome is spanned by exons, whereas 24% is in introns, with 75% of the genome being intergenic DNA. Duplications of segmental blocks, ranging in size up to chromosomal lengths, are abundant throughout the genome and reveal a complex evolutionary history. Comparative genomic analysis indicates vertebrate expansions of genes associated with neuronal function, with tissue-specific developmental regulation, and with the hemostasis and immune systems. DNA sequence comparisons between the consensus sequence and publicly funded genome data provided locations of 2.1 million single-nucleotide polymorphisms (SNPs). A random pair of human haploid genomes differed at a rate of 1 bp per 1250 on average, but there was marked heterogeneity in the level of polymorphism across the genome. Less than 1% of all SNPs resulted in variation in proteins, but the task of determining which SNPs have functional consequences remains an open challenge.

http://science.sciencemag.org/content/sci/291/5507/1304.full.pdf

The sequence of the human genome.

Conserved seed pairing, often flanked by adenosines, indicates that thousands of human genes are microRNA targets

As vertebrate genome sequences near completion and research refocuses to their analysis, the issue of effective genome annotation display becomes critical. A mature web tool for rapid and reliable display of any requested portion of the genome at any scale, together with several dozen aligned annotation tracks, is provided at http://genome.ucsc.edu. This browser displays assembly contigs and gaps, mRNA and expressed sequence tag alignments, multiple gene predictions, cross-species homologies, single nucleotide polymorphisms, sequence-tagged sites, radiation hybrid data, transposon repeats, and more as a stack of coregistered tracks. Text and sequence-based searches provide quick and precise access to any region of specific interest. Secondary links from individual features lead to sequence details and supplementary off-site databases. One-half of the annotation tracks are computed at the University of California, Santa Cruz from publicly available sequence data; collaborators worldwide provide the rest. Users can stably add their own custom tracks to the browser for educational or research purposes. The conceptual and technical framework of the browser, its underlying MYSQL database, and overall use are described. The web site currently serves over 50,000 pages per day to over 3000 different users.

/pdf/the-human-genome-browser-at-ucsc-4zclgthyly.pdf

The Human Genome Browser at UCSC

The Center for Disease Control estimates 128,000 people in the U.S. are hospitalized annually due to food borne illnesses. This has created a demand for food safety testing targeting the detection of pathogenic mold and bacteria on agricultural products. This risk extends to medical Cannabis and is of particular concern with inhaled, vaporized and even concentrated Cannabis products. As a result, third party microbial testing has become a regulatory requirement in the medical and recreational Cannabis markets, yet knowledge of the Cannabis microbiome is limited. Here we describe the first next generation sequencing survey of the microbial communities found in dispensary based Cannabis flowers and demonstrate the limitations in the culture based regulations that are being superimposed from the food industry.

/pdf/cannabis-microbiome-sequencing-reveals-several-mycotoxic-50slat85jz.pdf

Cannabis microbiome sequencing reveals several mycotoxic fungi native to dispensary grade cannabis flowers

Rare genetic variants in the core endocannabinoid system genes CNR1, CNR2, DAGLA, MGLL and FAAH were identified in molecular testing data from up to 6.032 patients with a broad spectrum of neurological disorders. The variants were evaluated for association with phenotypes similar to those observed in the orthologous gene knockouts in mice. Heterozygous rare coding variants in CNR1, which encodes the type 1 cannabinoid receptor (CB1), were found to be significantly associated with pain sensitivity (especially migraine), sleep and memory disorders - alone or in combination with anxiety - compared to a set of controls without such CNR1 variants. Similarly, heterozygous rare variants in DAGLA, which encodes diacylglycerol lipase alpha, were found to be significantly associated with seizures and developmental disorders, including abnormalities of brain morphology, compared to controls. Rare variants in MGLL, FAAH and CNR2 were not associated with any neurological phenotypes in the patients tested. Diacylglycerol lipase alpha synthesizes the endocannabinoid 2-AG in the brain, which interacts with CB1 receptors. The phenotypes associated with rare CNR1 variants are reminiscent of those implicated in the theory of clinical endocannabinoid deficiency syndrome. The severe phenotypes associated with rare DAGLA variants underscore the critical role of rapid 2-AG synthesis and the endocannabinoid system in regulating neurological function and development. Mapping of the variants to the 3D structure of the type 1 cannabinoid receptor, or primary structure of diacylglycerol lipase alpha, reveals clustering of variants in certain structural regions and is consistent with impacts to function.

/pdf/rare-genetic-variants-in-the-endocannabinoid-system-genes-1b2f36zhwz.pdf

Rare genetic variants in the endocannabinoid system genes CNR1 and DAGLA are associated with neurological phenotypes in humans

Publisher Summary This chapter describes specific methods for generating yeast artificial chromosome (YAC) clones and for amplifying the artificial chromosomes they contain. Yeast artificial chromosome (YAC) cloning is a useful method for propagating very large DNA segments, from 100 kilobases (kb) to over 1 megabase (Mb) in size. However, artificial chromosomes are normally present in yeast at only one copy per haploid genome, hence, a number of routine applications, (such as colony-screening and preparative-scale DNA isolation) are either inconvenient to perform or are plagued with problems of low sensitivity. This chapter discusses the system that permits copy number amplification of YACs by introducing certain cis-acting elements into the cloning vector, pCGS966.

Copy number amplification of yeast artificial chromosomes.

The present teachings relate to improved methods, kits, and compositions for making nucleic acid libraries and sequencing nucleic acids. In some embodiments, directionally defined concatamers are generated, facilitating sequencing efforts.

Directed assembly of amplicons to enhance read pairing signature with massively parallel short read sequencers

Purpose We aimed to gain insight into frequencies of genetic variants in genes implicated in neurodevelopmental disorder with epilepsy (NDD+E) by investigating large cohorts of patients in a diagnostic setting. Methods We analyzed variants in NDD+E using epilepsy gene panel sequencing performed between 2013 and 2017 by two large diagnostic companies. We compared variant frequencies in 6,994 panels to other 8,588 recently published panels as well as exome-wide de novo variants in 1,942 individuals with NDD+E and 10,937 controls. Results Genes with highest frequencies of ultra-rare variants in NDD+E comprised SCN1A, KCNQ2, SCN2A, CDKL5, SCN8A and STXBP1, concordant with the two other epilepsy cohorts we investigated. Only 46% of the analysed 262 dominant and X-linked panel genes contained ultra-rare variants in patients. Among genes with contradictory evidence of association with epilepsy CACNB4, CLCN2, EFHC1, GABRD, MAGI2 and SRPX2 showed equal frequencies in cases and controls. Conclusion We show that improvement of panel design increased diagnostic yield over time, but panels still display genes with low or no diagnostic yield. With our data, we hope to improve current diagnostic NDD+E panel design and provide a resource of ultra-rare variants in individuals with NDD+E to the community.

/pdf/targeted-gene-sequencing-in-6994-individuals-with-4yxrzs6dc8.pdf

Douglas R. Smith

Papers

Cannabis microbiome sequencing reveals several mycotoxic fungi native to dispensary grade cannabis flowers

Rare genetic variants in the endocannabinoid system genes CNR1 and DAGLA are associated with neurological phenotypes in humans

Copy number amplification of yeast artificial chromosomes.

Directed assembly of amplicons to enhance read pairing signature with massively parallel short read sequencers

Targeted gene sequencing in 6994 individuals with neurodevelopmental disorder with epilepsy