What Can Genome-Wide Association Studies Tell Us about the Genetics of Common Disease

The success of genome-wide association studies relies on much of the risk of common diseases being due to common genetic variants; but evidence for this is inconclusive. The results of published genome-wide association studies are examined to see what can be learnt about the distribution of disease-associated variants and how this might influence future study design. Although replicated disease-associated variants tend to be very common and frequency is inversely correlated with estimated effect size, our simulations suggest that such observations are the result of power. We find that for studies conducted to date, the frequency and effect size of significantly associated alleles are likely to be similar to those of the underlying disease alleles that they represent. Little of the genetic variation of disease has been explained so far, but current studies are only adequately powered to detect very common alleles unless they greatly increase disease risk. Thus, although the truth of the common disease / common variant hypothesis remains undecided, recent successes suggest that there are many more common genetic disease-associated variants, requiring larger studies to be identified.     

Annotating functional effects of non-coding variants in neuropsychiatric cell types by deep transfer learning

Genomewide association studies (GWAS) have identified a large number of loci associated with neuropsychiatric traits, however, understanding the molecular mechanisms underlying these loci remains difficult. To help prioritize causal variants and interpret their functions, computational methods have been developed to predict regulatory effects of non-coding variants. An emerging approach to variant annotation is deep learning models that predict regulatory functions from DNA sequences alone. While such models have been trained on large publicly available dataset such as ENCODE, neuropsychiatric trait-related cell types are under-represented in these datasets, thus there is an urgent need of better tools and resources to annotate variant functions in such cellular contexts. To fill this gap, we collected a large collection of neurodevelopment-related cell/tissue types, and trained deep Convolutional Neural Networks (ResNet) using such data. Furthermore, our model, called MetaChrom, borrows information from public epigenomic consortium to improve the accuracy via transfer learning. We show that MetaChrom is substantially better in predicting experimentally determined chromatin accessibility variants than popular variant annotation tools such as CADD and delta-SVM. By combining GWAS data with MetaChrom predictions, we prioritized 31 SNPs for Schizophrenia, suggesting potential risk genes and the biological contexts where they act. In summary, MetaChrom provides functional annotations of any DNA variants in the neuro-development context and the general method of MetaChrom can also be extended to other disease-related cell or tissue types.     

High-throughput SNP detection by using DNA pooling and denaturing high performance liquid chromatography (DHPLC)

One of the critical steps in the positional cloning of a complex disease gene involves association analysis between a phenotype and a set of densely spaced diallelic markers, typically single nucleotide repeats (SNPs), covering the region of interest. However, the effort and cost of detecting sufficient numbers of SNPs across relatively large physical distances represents a significant rate-limiting step. We have explored DNA pooling, in conjunction with denaturing high performance liquid chromatography (DHPLC), as a possible strategy for augmenting the efficiency, economy, and throughput of SNP detection. DHPLC is traditionally used to detect variants in polymerase chain reaction products containing both allelic forms of a polymorphism (e.g., heterozygotes or a 1:1 mix of both alleles) via heteroduplex separation and thereby requires separate analyses of multiple individual test samples. We have adapted this technology to identify variants in pooled DNA. To evaluate the utility and sensitivity of this approach, we constructed DNA pools comprised of 20 previously genotyped individuals with a frequency representation of 0%-50% for the variant allele. Mutation detection was performed by using temperature-modulated heteroduplex formation/DHPLC and dye-terminator sequencing. Using DHPLC, we could consistently detect SNPs at lower than 5% frequency, corresponding to the detection of one variant allele in a pool of 20 alleles. In contrast, fluorescent sequencing detected variants in the same pools only if the frequency of the less common allele was at least 10%. We conclude that DNA pooling of samples for DHPLC analysis is an effective way to increase throughput efficiency of SNP detection.     

Construction of a viable SV40 variant containing two functional origins of DNA replication.

Viable variants of simian virus 40 (SV40) have been constructed which contain two functional origins of DNA replication (Or). The variants were made by introducing, at 0.175 on the SV40 map, a segment of DNA containing the viral Or. Two types of experiments demonstrate that the second Or is functional. First, the distribution of radioactivity in pulse-labeled SV40 (I) DNA is dramatically altered in the variants when compared with the parental virus. Second, electron microscopic examination of viral replicative intermediates indicates that while there is one initiation site for DNA synthesis in the parental genome, there are two sites in the variant. It was possible to introduce a deletion which inactivated the original Or at 0.67 map units in this variant. The resulting mutant could be propagated, and its DNA replication originated at the site of the newly inserted Or.     

Tool evaluation for the detection of variably sized indels from next generation whole genome and targeted sequencing data

Insertions and deletions (indels) in human genomes are associated with a wide range of phenotypes, including various clinical disorders. High-throughput, next generation sequencing (NGS) technologies enable the detection of short genetic variants, such as single nucleotide variants (SNVs) and indels. However, the variant calling accuracy for indels remains considerably lower than for SNVs. Here we present a comparative study of the performance of variant calling tools for indel calling, evaluated with a wide repertoire of NGS datasets. While there is no single optimal tool to suit all circumstances, our results demonstrate that the choice of variant calling tool greatly impacts the precision and recall of indel calling. Furthermore, to reliably detect indels, it is essential to choose NGS technologies that offer a long read length and high coverage coupled with specific variant calling tools.     

Linkage analysis in the next-generation sequencing era

Linkage analysis was developed to detect excess co-segregation of the putative alleles underlying a phenotype with the alleles at a marker locus in family data. Many different variations of this analysis and corresponding study design have been developed to detect this co-segregation. Linkage studies have been shown to have high power to detect loci that have alleles (or variants) with a large effect size, i.e. alleles that make large contributions to the risk of a disease or to the variation of a quantitative trait. However, alleles with a large effect size tend to be rare in the population. In contrast, association studies are designed to have high power to detect common alleles which tend to have a small effect size for most diseases or traits. Although genome-wide association studies have been successful in detecting many new loci with common alleles of small effect for many complex traits, these common variants often do not explain a large proportion of disease risk or variation of the trait. In the past, linkage studies were successful in detecting regions of the genome that were likely to harbor rare variants with large effect for many simple Mendelian diseases and for many complex traits. However, identifying the actual sequence variant(s) responsible for these linkage signals was challenging because of difficulties in sequencing the large regions implicated by each linkage peak. Current 'next-generation' DNA sequencing techniques have made it economically feasible to sequence all exons or the whole genomes of a reasonably large number of individuals. Studies have shown that rare variants are quite common in the general population, and it is now possible to combine these new DNA sequencing methods with linkage studies to identify rare causal variants with a large effect size. A brief review of linkage methods is presented here with examples of their relevance and usefulness for the interpretation of whole-exome and whole-genome sequence data.     

DNA sequence variants in the G gamma-, A gamma-, delta- and beta-globin genes of man

DNA prepared from 60 unrelated individuals was cleaved with one of eight different restriction endonucleases and the resulting DNA fragments were separated by agarose gel electrophoresis. DNA fragments containing G gamma-, A gamma-, delta- or beta-globin genes were detected by Southern blot hybridization, using as probe either a 32P-labeled cloned DNA copy of rabbit beta-globin messenger RNA or labeled human beta- and G gamma- globin cDNA plasmids. Three types of variant restriction enzyme patterns of globin DNA fragments were detected in otherwise normal individuals. One variant pattern, found in only one person, was caused by an additional restriction endonuclease Pst I cleavage site in the center of the delta- globin gene intervening sequence; the subject was heterozygous for the presence of this cleavage site and was shown to have inherited it from her mother. Another variant pattern resulted from the appearance of an endonuclease Hind III cleavage site in the intervening sequence of the A gamma-globin gene; this variant is polymorphic, with a gene frequency for the presence of the intragenic Hind III site of 0.23. This Hind III cleavage site polymorphism is also found in the G gamma-globin gene intervening sequence and thus the polymorphism itself appears to be duplicated over the pair of gamma-globin loci. These variants can be used to derive an approximate estimate of the total number of different DNA sequence variants in man.     

Efficient and cost effective population resequencing by pooling and in-solution hybridization

High-throughput sequencing of targeted genomic loci in large populations is an effective approach for evaluating the contribution of rare variants to disease risk. We evaluated the feasibility of using in-solution hybridization-based target capture on pooled DNA samples to enable cost-efficient population sequencing studies. For this, we performed pooled sequencing of 100 HapMap samples across ～ 600 kb of DNA sequence using the Illumina GAIIx. Using our accurate variant calling method for pooled sequence data, we were able to not only identify single nucleotide variants with a low false discovery rate (<1%) but also accurately detect short insertion/deletion variants. In addition, with sufficient coverage per individual in each pool (30-fold) we detected 97.2% of the total variants and 93.6% of variants below 5% in frequency. Finally, allele frequencies for single nucleotide variants (SNVs) estimated from the pooled data and the HapMap genotype data were tightly correlated (correlation coefficient > = 0.995).     

Development of molecular approaches to estimating germinal mutation rates. I. Detection of insertion/deletion/rearrangement variants in the human genome

DNA from 130 individuals was studied with up to 18 (primarily cDNA) probes for the frequency of variants in this initial experiment to determine the feasibility of this approach to screening for germinal gene mutations. This approach, a modification of the usual restriction enzyme mapping strategy, focuses on the detection of insertion/deletion/rearrangement (I/D/R) variants, because the DNA is digested with only two restriction enzymes before transfer to membranes and hybridization with an extensive series of unrelated probes. Some 4000 noncontiguous, independent DNA fragments ("loci"), functional loci, pseudogenes or anonymous fragments, (a total of approximately 77,400 kb) were screened. 19 different classes and 31 copies of presumably I/D/R variants were detected while 4 different classes and 24 individuals exhibiting base substitution variants were observed. 18 of the 19 I/D/R classes were rare variants, that is, each were observed at a frequency, within this population, of less than 0.01; 3 of the base substitution classes existed at polymorphic frequencies and only 1 was a rare variant. 10 of the I/D/R classes, occurring in a total of 18 individuals, were detected with probes which are not known to be associated with repetitive elements. This is a variant frequency for I/D/R variants without known repetitive elements of 0.15 classes and 0.23 copies for each 1000 kb screened; this would extrapolate to 1600 such variant sites in the genome of each individual. Within the context of a mutation screening program, the rare variants, either with or without repetitive elements, would have a higher probability of being de novo mutations than would polymorphic variants; this former group would be the focus of family studies to test for the heritability of the allele (fragment pattern). Sufficient DNA probes are available to screen a significant portion of the human genome for genetic variation and de novo mutations of this type.     

Physicochemical properties of temperate phage DNA and its possible effect on the variability of Mycobacterium lacticolum

The temperate phage 104 S was isolated from the S variant of Mycobacterium lacticolum, strain 104, and some of its characteristics were studied. The content of GC pairs in the phage DNA was 77 mole% as was calculated from the melting profile or 65 mole% as was calculated from the value of buoyant density in CsCl. The DNA was shown to be composed of 18,000 nucleotide pairs. DNA restriction fragments of M. lacticolum R, S and M variants were subjected for the first time to molecular hybridization with [32P]DNA of the temperate phage. The genome of the three M. lacticolum variants and the genome of a non-dissociating S variant clone were shown to contain sequences homologous to the DNA sequence of phage 104 S. Differences are found among the variants in the hybridizing DNA fragments. These data indicate that the phage DNA may actively be involved in the variability of the culture. Its participation can be realized by the different mode of prophage incorporation into the genome of the variants.     

DNA Sequence Variants in the Five Prime Untranslated Region of the Cyclooxygenase-2 Gene Are Commonly Found in Healthy Dogs and Gray Wolves

The aim of this study was to investigate the frequency of regional DNA variants upstream to the translation initiation site of the canine Cyclooxygenase-2 (Cox-2) gene in healthy dogs. Cox-2 plays a role in various disease conditions such as acute and chronic inflammation, osteoarthritis and malignancy. A role for Cox-2 DNA variants in genetic predisposition to canine renal dysplasia has been proposed and dog breeders have been encouraged to select against these DNA variants. We sequenced 272–422 bases in 152 dogs unaffected by renal dysplasia and found 19 different haplotypes including 11 genetic variants which had not been described previously. We genotyped 7 gray wolves to ascertain the wildtype variant and found that the wolves we analyzed had predominantly the second most common DNA variant found in dogs. Our results demonstrate an elevated level of regional polymorphism that appears to be a feature of healthy domesticated dogs.     

Combinatorial redesign of the DNA binding specificity of a prokaryotic helix-turn-helix repressor

Redesign of the bacteriophage 434 Cro repressor was accomplished by using an in vivo genetic screening system to identify new variants that specifically bound previously unrecognized DNA sequences. Site-directed, combinatorial mutagenesis of the 434 Cro helix-turn-helix (HTH) motif generated libraries of new variants which were screened for binding to new target sequences. Multiple mutations of 434 Cro that functionally converted wild-type (wt) 434 Cro DNA binding-sequence specificity to that of a lambda bacteriophage-specific repressor were identified. The libraries contained variations within the HTH sequence at only three positions. In vivo and in vitro analysis of several of the identified 434 Cro variants showed that the relatively few changes in the recognition helix of the HTH motif of 434 Cro resulted in specific and tight binding of the target DNA sequences. For the best 434 Cro variant identified, an apparent K(d) for lambda O(R)3 of 1 nM was observed. In competition experiments, this Cro variant was observed to be highly selective. We conclude that functional 434 Cro repressor variants with new DNA binding specificities can be generated from wt 434 Cro by mutating just the recognition helix. Important characteristics of the screening system responsible for the successful identifications are discussed. Application of the techniques presented here may allow the identification of DNA binding protein variants that functionally affect DNA regulatory sequences important in disease and industrial and biotechnological processes.     

Analysis of archived residual newborn screening blood spots after whole genome amplification

Background

Deidentified newborn screening bloodspot samples (NBS) represent a valuable potential resource for genomic research if impediments to whole exome sequencing of NBS deoxyribonucleic acid (DNA), including the small amount of genomic DNA in NBS material, can be overcome. For instance, genomic analysis of NBS could be used to define allele frequencies of disease-associated variants in local populations, or to conduct prospective or retrospective studies relating genomic variation to disease emergence in pediatric populations over time. In this study, we compared the recovery of variant calls from exome sequences of amplified NBS genomic DNA to variant calls from exome sequencing of non-amplified NBS DNA from the same individuals.

Results

Using a standard alignment-based Genome Analysis Toolkit (GATK), we find 62,000–76,000 additional variants in amplified samples. After application of a unique kmer enumeration and variant detection method (RUFUS), only 38,000–47,000 additional variants are observed in amplified gDNA. This result suggests that roughly half of the amplification-introduced variants identified using GATK may be the result of mapping errors and read misalignment.

Conclusions

Our results show that it is possible to obtain informative, high-quality data from exome analysis of whole genome amplified NBS with the important caveat that different data generation and analysis methods can affect variant detection accuracy, and the concordance of variant calls in whole-genome amplified and non-amplified exomes.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1747-2) contains supplementary material, which is available to authorized users.     

Recombination in simian virus 40-infected cells. Structure of naturally arising variants ev-2114, ev-2102, and ev-1110

The complete nucleotide sequence has been determined for three newly cloned evolutionary variants from two different independently generated evolutionary series (1100 and 2100 series) of simian virus 40 (SV40). These naturally arising variants, designated ev-1110, ev-2102, and ev-2114, were isolated after five high multiplicity serial passages. The structure of the variants consists of a monomeric unit tandemly repeated four times (ev-2102 and ev-2114) or six times (ev-1110) in the variant genome; the variants have four or six copies, respectively, of the viral origin signal for DNA replication. The DNA content in the three variants is vastly different in that the genome of variant ev-2114 contains only rearranged viral sequences, while variant ev-2102 contains a substitution with monkey DNA sequences consisting of a nearly complete dimeric unit of Alu family sequences as well as less repetitive sequences and variant ev-1110 contains monkey DNA sequences derived solely from repetitive alpha-component DNA. Recombination events, cellular sequences, and structural features of these and other naturally arising SV40 variants are compared.     

Germline EPHB2 receptor variants in familial colorectal cancer

Familial clustering of colorectal cancer occurs in 15-20% of cases, however recognized cancer syndromes explain only a small fraction of this disease. Thus, the genetic basis for the majority of hereditary colorectal cancer remains unknown. EPHB2 has recently been implicated as a candidate tumor suppressor gene in colorectal cancer. The aim of this study was to evaluate the contribution of EPHB2 to hereditary colorectal cancer. We screened for germline EPHB2 sequence variants in 116 population-based familial colorectal cancer cases by DNA sequencing. We then estimated the population frequencies and characterized the biological activities of the EPHB2 variants identified. Three novel nonsynonymous missense alterations were detected. Two of these variants (A438T and G787R) result in significant residue changes, while the third leads to a conservative substitution in the carboxy-terminal SAM domain (V945I). The former two variants were found once in the 116 cases, while the V945I variant was present in 2 cases. Genotyping of additional patients with colorectal cancer and control subjects revealed that A438T and G787R represent rare EPHB2 alleles. In vitro functional studies show that the G787R substitution, located in the kinase domain, causes impaired receptor kinase activity and is therefore pathogenic, whereas the A438T variant retains its receptor function and likely represents a neutral polymorphism. Tumor tissue from the G787R variant case manifested loss of heterozygosity, with loss of the wild-type allele, supporting a tumor suppressor role for EPHB2 in rare colorectal cancer cases. Rare germline EPHB2 variants may contribute to a small fraction of hereditary colorectal cancer.     

Generation and chracterization of variants of mouse hepatoma cells with defects in hepato-specific gene expression. I. Albumin synthesis variants

Clonal variants of mouse hepatoma cells that either fail to produce albumin (variant 19/2) or show significantly reduced levels (100-fold less) of albumin production (variant 1/c/1) were isolated from the parental line. Hepa la, after a single exposure to N-methyl-N'-nitrosoguanidine (MNNG). Intracellular levels of albumin in both variants were below detection by our assay. Analyses by cDNA-RNA reassociation kinetics indicate that there are approximately 3900 molecules of cytoplasmic albumin mRNA per cell in the parent and less than 10 molecules per cell in both variants. Southern blotting of the Eco RI restriction fragments of cellular DNA from the parent and variants did not indicate any major deletions in the albumin gene DNA sequences. We conclude that in the two variants studied, processes that regulate albumin production via alterations in the level of cytoplasmic albumin mRNA have been affected. Our analyses have also shown that alpha-fetoprotein (AFP) production is lacking in one variant (19/2) and is slightly reduced in the other (1/c/1). Transferrin secretion is lower than the parental line in both variants. Thus multiple nonlethal defects in hepatic gene expression can be obtained in Hepa la cells in culture that will be useful in determining the number and kinds of genes that control the expression of liver-specific loci.     

Fluorescent-labeled single-strand ATP aptamer DNA: chemo- and enantio-selectivity in sensing adenosine

One of the intriguing applications of aptamers is sensing molecules. In principle, an aptamer can specifically recognize and bind to a unique ligand, leading to a structural change of an aptamer. By acquiring information for the structural change, the detection of the ligand can be achieved. To design and explore an aptamer molecule to detect adenosine, we have synthesized some ATP aptamer variants labeled with donor and acceptor fluorophores. Although the fluorescent response of the aptamer variants was highly dependent on experimental temperature, we have found one of the variants showing suitable fluorescent response by titration with adenosine. The aptamer variant showed remarkable selectivity for adenosine over the other ribonucleosides. On the other hand, the enantio-specificity of the aptamer variant in the ligand recognition was not enough to selectively detect d-adenosine over l-adenosine.     

Restriction endonuclease DNA analysis of antigenic variants of leptospires selected by monoclonal antibodies

The genome of antigenic variant CV (CT3)-1 derived from Leptospira interrogans serovar canicola was compared by cleavage with restriction endonucleases with the parent and serovar bafani, to which the variant was serologically most closely related. No differences were observed between the parent and variant in DNA restriction endonuclease patterns using eight restriction endonucleases. Serovar bafani was different in the patterns from the parent and antigenic variant CV (CT3)-1. The two antigenic variants derived from serovar hebdomadis, HV (H16)-1 and HV (H19)-1 which belonged serologically to serovars jules and hebdomadis, respectively, were compared by restriction endonuclease DNA analysis with the parent and serovar jules. No differences were observed between the parent and variants in DNA restriction endonuclease patterns using the same enzymes. But some differences were observed in DNA restriction endonuclease patterns between HV (H16)-1 and serovar jules. Thus, the antigenic variant selected from the parent by the anti-parent monoclonal antibody and serologically different from the parent, being identified either as a new serovar or as a known one, was found to be similar to the parent by the restriction endonuclease DNA analysis.