Crowdsourcing human mutations

The first Human Variome microattribution review shows that data citation and publication credit can work as incentives for systematic curation of gene variant and phenotype data. Analysis of the formal assertions in both databases and journal articles argues for better separation of data structures from narrative so that they can better support one another to communicate meaning.     

Estimating coverage and power for genetic association studies using near-complete variation data

Although studies suggest that SNPs derived from HapMap provide promising coverage and power for association studies, the lack of alternative variation datasets limits independent analysis. Using near-complete variation data for 76 genes resequenced in HapMap samples, we find that coverage of common variation by commercial genotyping arrays is substantially lower compared to the HapMap-based estimates. We quantify the power offered by these arrays for a range of disease models.     

Analysis of the human protein interactome and comparison with yeast, worm and fly interaction datasets

We present the first analysis of the human proteome with regard to interactions between proteins. We also compare the human interactome with the available interaction datasets from yeast (Saccharomyces cerevisiae), worm (Caenorhabditis elegans) and fly (Drosophila melanogaster). Of >70,000 binary interactions, only 42 were common to human, worm and fly, and only 16 were common to all four datasets. An additional 36 interactions were common to fly and worm but were not observed in humans, although a coimmunoprecipitation assay showed that 9 of the interactions do occur in humans. A re-examination of the connectivity of essential genes in yeast and humans indicated that the available data do not support the presumption that the number of interaction partners can accurately predict whether a gene is essential. Finally, we found that proteins encoded by genes mutated in inherited genetic disorders are likely to interact with proteins known to cause similar disorders, suggesting the existence of disease subnetworks. The human interaction map constructed from our analysis should facilitate an integrative systems biology approach to elucidating the cellular networks that contribute to health and disease states.     

Enzymatic production of RNAi libraries from cDNAs

RNA interference (RNAi) induced by small interfering (siRNA) or short hairpin RNA (shRNA) is an important research approach in mammalian genetics. Here we describe a technology called enzymatic production of RNAi library (EPRIL) by which cDNAs are converted by a sequence of enzymatic treatments into an RNAi library consisting of a vast array of different shRNA expression constructs. We applied EPRIL to a single cDNA source and prepared an RNAi library consisting of shRNA constructs with various RNAi efficiencies. High-throughput screening allowed us to rapidly identify the best shRNA constructs from the library. We also describe a new selection scheme using the thymidine kinase gene for obtaining efficient shRNA constructs. Furthermore, we show that EPRIL can be applied to constructing an RNAi library from a cDNA library, providing a basis for future whole-genome phenotypic screening of genes.     

Interaction between ERAP1 and HLA-B27 in ankylosing spondylitis implicates peptide handling in the mechanism for HLA-B27 in disease susceptibility

Ankylosing spondylitis is a common form of inflammatory arthritis predominantly affecting the spine and pelvis that occurs in approximately 5 out of 1,000 adults of European descent. Here we report the identification of three variants in the RUNX3, LTBR-TNFRSF1A and IL12B regions convincingly associated with ankylosing spondylitis (P < 5 × 10(-8) in the combined discovery and replication datasets) and a further four loci at PTGER4, TBKBP1, ANTXR2 and CARD9 that show strong association across all our datasets (P < 5 × 10(-6) overall, with support in each of the three datasets studied). We also show that polymorphisms of ERAP1, which encodes an endoplasmic reticulum aminopeptidase involved in peptide trimming before HLA class I presentation, only affect ankylosing spondylitis risk in HLA-B27-positive individuals. These findings provide strong evidence that HLA-B27 operates in ankylosing spondylitis through a mechanism involving aberrant processing of antigenic peptides.     

SUMO modification is required for in vivo Hox gene regulation by the Caenorhabditis elegans Polycomb group protein SOP-2

Post-translational modification of proteins by the ubiquitin-like molecule SUMO (sumoylation) regulates their subcellular localization and affects their functional properties in vitro, but the physiological function of sumoylation in multicellular organisms is largely unknown. Here, we show that the C. elegans Polycomb group (PcG) protein SOP-2 interacts with the SUMO-conjugating enzyme UBC-9 through its evolutionarily conserved SAM domain. Sumoylation of SOP-2 is required for its localization to nuclear bodies in vivo and for its physiological repression of Hox genes. Global disruption of sumoylation phenocopies a sop-2 mutation by causing ectopic Hox gene expression and homeotic transformations. Chimeric constructs in which the SOP-2 SAM domain is replaced with that derived from fruit fly or mammalian PcG proteins, but not those in which the SOP-2 SAM domain is replaced with the SAM domains of non-PcG proteins, confer appropriate in vivo nuclear localization and Hox gene repression. These observations indicate that sumoylation of PcG proteins, modulated by their evolutionarily conserved SAM domain, is essential to their physiological repression of Hox genes.     

A framework for variation discovery and genotyping using next-generation DNA sequencing data

Recent advances in sequencing technology make it possible to comprehensively catalog genetic variation in population samples, creating a foundation for understanding human disease, ancestry and evolution. The amounts of raw data produced are prodigious, and many computational steps are required to translate this output into high-quality variant calls. We present a unified analytic framework to discover and genotype variation among multiple samples simultaneously that achieves sensitive and specific results across five sequencing technologies and three distinct, canonical experimental designs. Our process includes (i) initial read mapping; (ii) local realignment around indels; (iii) base quality score recalibration; (iv) SNP discovery and genotyping to find all potential variants; and (v) machine learning to separate true segregating variation from machine artifacts common to next-generation sequencing technologies. We here discuss the application of these tools, instantiated in the Genome Analysis Toolkit, to deep whole-genome, whole-exome capture and multi-sample low-pass (～4×) 1000 Genomes Project datasets.     

Bayesian inference analyses of the polygenic architecture of rheumatoid arthritis

The genetic architectures of common, complex diseases are largely uncharacterized. We modeled the genetic architecture underlying genome-wide association study (GWAS) data for rheumatoid arthritis and developed a new method using polygenic risk-score analyses to infer the total liability-scale variance explained by associated GWAS SNPs. Using this method, we estimated that, together, thousands of SNPs from rheumatoid arthritis GWAS explain an additional 20% of disease risk (excluding known associated loci). We further tested this method on datasets for three additional diseases and obtained comparable estimates for celiac disease (43% excluding the major histocompatibility complex), myocardial infarction and coronary artery disease (48%) and type 2 diabetes (49%). Our results are consistent with simulated genetic models in which hundreds of associated loci harbor common causal variants and a smaller number of loci harbor multiple rare causal variants. These analyses suggest that GWAS will continue to be highly productive for the discovery of additional susceptibility loci for common diseases.     

Common variants at ABCA7, MS4A6A/MS4A4E, EPHA1, CD33 and CD2AP are associated with Alzheimer's disease

We sought to identify new susceptibility loci for Alzheimer's disease through a staged association study (GERAD+) and by testing suggestive loci reported by the Alzheimer's Disease Genetic Consortium (ADGC) in a companion paper. We undertook a combined analysis of four genome-wide association datasets (stage 1) and identified ten newly associated variants with P ≤ 1 × 10(-5). We tested these variants for association in an independent sample (stage 2). Three SNPs at two loci replicated and showed evidence for association in a further sample (stage 3). Meta-analyses of all data provided compelling evidence that ABCA7 (rs3764650, meta P = 4.5 × 10(-17); including ADGC data, meta P = 5.0 × 10(-21)) and the MS4A gene cluster (rs610932, meta P = 1.8 × 10(-14); including ADGC data, meta P = 1.2 × 10(-16)) are new Alzheimer's disease susceptibility loci. We also found independent evidence for association for three loci reported by the ADGC, which, when combined, showed genome-wide significance: CD2AP (GERAD+, P = 8.0 × 10(-4); including ADGC data, meta P = 8.6 × 10(-9)), CD33 (GERAD+, P = 2.2 × 10(-4); including ADGC data, meta P = 1.6 × 10(-9)) and EPHA1 (GERAD+, P = 3.4 × 10(-4); including ADGC data, meta P = 6.0 × 10(-10)).     

Meta-analysis identifies 29 additional ulcerative colitis risk loci, increasing the number of confirmed associations to 47

Genome-wide association studies and candidate gene studies in ulcerative colitis have identified 18 susceptibility loci. We conducted a meta-analysis of six ulcerative colitis genome-wide association study datasets, comprising 6,687 cases and 19,718 controls, and followed up the top association signals in 9,628 cases and 12,917 controls. We identified 29 additional risk loci (P < 5 × 10(-8)), increasing the number of ulcerative colitis-associated loci to 47. After annotating associated regions using GRAIL, expression quantitative trait loci data and correlations with non-synonymous SNPs, we identified many candidate genes that provide potentially important insights into disease pathogenesis, including IL1R2, IL8RA-IL8RB, IL7R, IL12B, DAP, PRDM1, JAK2, IRF5, GNA12 and LSP1. The total number of confirmed inflammatory bowel disease risk loci is now 99, including a minimum of 28 shared association signals between Crohn's disease and ulcerative colitis.     

Towards genetic genome projects: genomic library screening and gene-targeting vector construction in a single step.

We have developed technologies that simplify genomic library construction and screening, substantially reducing both the time and the cost associated with traditional library screening methods and facilitating the generation of gene-targeting constructs. By taking advantage of homologous recombination in Escherichia coli, we were able to use as little as 80 bp of total sequence homology to screen for a specific gene from a genomic library in plasmid or phage form. This method, called recombination cloning (REC), takes only a few days instead of the several weeks required for traditional plaque-lift methods. In addition, because every clone in the mouse genomic library we have constructed has a negative selection marker adjacent to the genomic insert, REC screening can generate gene-targeting vectors in one step, from library screening to finished construct. Conditional targeting constructs can be generated easily with minimal additional manipulation.     

Parameters for reliable results in genetic association studies in common disease.

It is increasingly apparent that the identification of true genetic associations in common multifactorial disease will require studies comprising thousands rather than the hundreds of individuals employed to date. Using 2,873 families, we were unable to confirm a recently published association of the interleukin 12B gene in 422 type I diabetic families. These results emphasize the need for large datasets, small P values and independent replication if results are to be reliable.     

MAGIC, an in vivo genetic method for the rapid construction of recombinant DNA molecules

We describe a highly engineered in vivo cloning method, mating-assisted genetically integrated cloning (MAGIC), that facilitates the rapid construction of recombinant DNA molecules. MAGIC uses bacterial mating, in vivo site-specific endonuclease cleavage and homologous recombination to catalyze the transfer of a DNA fragment between a donor vector in one bacterial strain and a recipient plasmid in a separate bacterial strain. Recombination events are genetically selected and result in placement of the gene of interest under the control of new regulatory elements with high efficiency. MAGIC eliminates the need for restriction enzymes, DNA ligases, preparation of DNA and all in vitro manipulations required for subcloning and allows the rapid construction of multiple constructs with minimal effort. We show that MAGIC can generate constructs for expression in multiple organisms. As this new method requires only the simple mixing of bacterial strains, it represents a substantial advance in high-throughput recombinant DNA production that will save time, effort and expense in functional genomics studies.     

Rescue of a deficiency in ATP synthesis by transfer of MTATP6, a mitochondrial DNA-encoded gene, to the nucleus

A T-->G transversion at nt 8993 in mitochondrial DNA of MTATP6 (encoding ATPase 6 of complex V of the respiratory chain) causes impaired mitochondrial ATP synthesis in two related mitochondrial disorders: neuropathy, ataxia and retinitis pigmentosa and maternally inherited Leigh syndrome. To overcome the biochemical defect, we expressed wildtype ATPase 6 protein allotopically from nucleus-transfected constructs encoding an amino-terminal mitochondrial targeting signal appended to a recoded ATPase 6 gene (made compatible with the universal genetic code) that also contained a carboxy-terminal FLAG epitope tag. After transfection of human cells, the precursor polypeptide was expressed, imported into and processed within mitochondria, and incorporated into complex V. Allotopic expression of stably transfected constructs in cytoplasmic hybrids (cybrids) homoplasmic with respect to the 8993T-->G mutation showed a significantly improved recovery after growth in selective medium as well as a significant increase in ATP synthesis. This is the first successful demonstration of allotopic expression of an mtDNA-encoded polypeptide in mammalian cells and could form the basis of a genetic approach to treat a number of human mitochondrial disorders.     

A copy number variation morbidity map of developmental delay

To understand the genetic heterogeneity underlying developmental delay, we compared copy number variants (CNVs) in 15,767 children with intellectual disability and various congenital defects (cases) to CNVs in 8,329 unaffected adult controls. We estimate that ～14.2% of disease in these children is caused by CNVs >400 kb. We observed a greater enrichment of CNVs in individuals with craniofacial anomalies and cardiovascular defects compared to those with epilepsy or autism. We identified 59 pathogenic CNVs, including 14 new or previously weakly supported candidates, refined the critical interval for several genomic disorders, such as the 17q21.31 microdeletion syndrome, and identified 940 candidate dosage-sensitive genes. We also developed methods to opportunistically discover small, disruptive CNVs within the large and growing diagnostic array datasets. This evolving CNV morbidity map, combined with exome and genome sequencing, will be critical for deciphering the genetic basis of developmental delay, intellectual disability and autism spectrum disorders.     

Restriction enzyme-generated siRNA (REGS) vectors and libraries

Small interfering RNA (siRNA) technology facilitates the study of loss of gene function in mammalian cells and animal models, but generating multiple siRNA vectors using oligonucleotides is slow, inefficient and costly. Here we describe a new, enzyme-mediated method for generating numerous functional siRNA constructs from any gene of interest or pool of genes. To test our restriction enzyme-generated siRNA (REGS) system, we silenced a transgene and two endogenous genes and obtained the predicted phenotypes. REGS generated on average 34 unique siRNAs per kilobase of sequence. REGS enabled us to create enzymatically a complex siRNA library (>4 x 10(5) clones) from double-stranded cDNA encompassing known and unknown genes with 96% of the clones containing inserts of the appropriate size.     

Learning to think continentally

The US Department of Health and Human Services (DHHS) is proposing to enhance federal regulation intended to protect human research subjects, in particular to increase measures aimed at security of personal data. Since the ethical review process is partially based on respect for people and their autonomy, harmonization of these rules will be a process of convincing individuals and their states to accept uniform standards that give enough privacy but do not lock away personal data from either research participants or researchers.