Quality control and integration of genotypes from two calling pipelines for whole genome sequence data in the Alzheimer's disease sequencing project

The Alzheimer's Disease Sequencing Project (ADSP) performed whole genome sequencing (WGS) of 584 subjects from 111 multiplex families at three sequencing centers. Genotype calling of single nucleotide variants (SNVs) and insertion-deletion variants (indels) was performed centrally using GATK-HaplotypeCaller and Atlas V2. The ADSP Quality Control (QC) Working Group applied QC protocols to project-level variant call format files (VCFs) from each pipeline, and developed and implemented a novel protocol, termed “consensus calling,” to combine genotype calls from both pipelines into a single high-quality set. QC was applied to autosomal bi-allelic SNVs and indels, and included pipeline-recommended QC filters, variant-level QC, and sample-level QC. Low-quality variants or genotypes were excluded, and sample outliers were noted. Quality was assessed by examining Mendelian inconsistencies (MIs) among 67 parent-offspring pairs, and MIs were used to establish additional genotype-specific filters for GATK calls. After QC, 578 subjects remained. Pipeline-specific QC excluded ~12.0% of GATK and 14.5% of Atlas SNVs. Between pipelines, ~91% of SNV genotypes across all QCed variants were concordant; 4.23% and 4.56% of genotypes were exclusive to Atlas or GATK, respectively; the remaining ~0.01% of discordant genotypes were excluded. For indels, variant-level QC excluded ~36.8% of GATK and 35.3% of Atlas indels. Between pipelines, ~55.6% of indel genotypes were concordant; while 10.3% and 28.3% were exclusive to Atlas or GATK, respectively; and ~0.29% of discordant genotypes were. The final WGS consensus dataset contains 27,896,774 SNVs and 3,133,926 indels and is publicly available.     

遗传距离的错误观点和病毒进化踪迹的探索

通过对病毒基因序列的遗传距离(K)的演算,发现公式中的p和q应理解为转换率和颠换率;计算所得的K值是一个百分率,更适合于显示生物进化的时间概念。因此建议将遗传距离改名为进化率。并从进化率计算中的规律性变化,探索了病毒进化的信息。     

Kinetics,Structure, and Mechanism of 8-Oxo-7,8-dihydro-2′-deoxyguanosine Bypass by Human DNA Polymerase η

DNA damage incurred by a multitude of endogenous and exogenous factors constitutes an inevitable challenge for the replication machinery. Cells rely on various mechanisms to either remove lesions or bypass them in a more or less error-prone fashion. The latter pathway involves the Y-family polymerases that catalyze trans-lesion synthesis across sites of damaged DNA. 7,8-Dihydro-8-oxo-2′-deoxyguanosine (8-oxoG) is a major lesion that is a consequence of oxidative stress and is associated with cancer, aging, hepatitis, and infertility. We have used steady-state and transient-state kinetics in conjunction with mass spectrometry to analyze in vitro bypass of 8-oxoG by human DNA polymerase η (hpol η). Unlike the high fidelity polymerases that show preferential insertion of A opposite 8-oxoG, hpol η is capable of bypassing 8-oxoG in a mostly error-free fashion, thus preventing GC→AT transversion mutations. Crystal structures of ternary hpol η-DNA complexes and incoming dCTP, dATP, or dGTP opposite 8-oxoG reveal that an arginine from the finger domain assumes a key role in avoiding formation of the nascent 8-oxoG:A pair. That hpol η discriminates against dATP exclusively at the insertion stage is confirmed by structures of ternary complexes that allow visualization of the extension step. These structures with G:dCTP following either 8-oxoG:C or 8-oxoG:A pairs exhibit virtually identical active site conformations. Our combined data provide a detailed understanding of hpol η bypass of the most common oxidative DNA lesion.     

Are transversion mutations better? A Mutagenesis Assistant Program analysis on P450 BM-3 heme domain

Directed evolution represents a versatile tool to tailor enzyme properties to needs in industrial applications and to understand structure-function relationships. Genetic diversity is commonly generated using error-prone PCR. Exploration of sequence space by random mutagenesis strongly favors transitions when enzyme-based mutagenesis methods are employed (Wong, T. S., Zhurina, D., Schwaneberg, U., Comb. Chem. High Throughput Screen. 2006, 9, 271-288). The genetic code has been organized in a manner that limits chemical diversity when a single transition mutation occurs in a codon (Wong, T. S., Roccatano, D., Schwaneberg, U., Biocatal. Biotransformation 2006, in press). Are transitions more beneficial than transversions for adapting biocatalysts to non-natural process conditions? In a statistical analysis performed with the Mutagenesis Assistant Program (MAP), we compared the consequences of transition and transversion bias on amino acid substitution patterns of the P450 BM-3 heme domain. For the analysis, we used a recently introduced benchmarking system consisting of a protein structure indicator, an amino acid diversity indicator with a codon diversity coefficient, and a chemical diversity indicator. A detailed analysis for the P450 BM-3 heme domain showed that an ideal transversion bias generates more diverse amino acid substitution patterns with a significantly different chemical composition than an ideal transition bias. Emphasis is given on the theoretical analysis with a brief discussion on potential implication of transition and transversion bias in directed evolution experiments.     

The study of neighboring nucleotide composition and transition/transversion bias

Base substitution is one of the raw fuels that produce genetic variation and drive evolution. Recent studies have shown that the genome components affect mutation patterns to some extent. In order to infer the correlation between the Transition/Transversion ratio (Ts/Tv) and the number of immediately adjacent A&T nucleotides, we investigated 3611007 Oryza sativa SNPs (including 45462 coding SNPs, and 242811 intronic SNPs) and 32019 Arabidopsis SNPs. The results show that Ts/Tv is negatively correlated with the number of immediately adjacent A&T in O. Sativa and Arabidopsis. We further calculated AT₂ (the number of SNPs whose immediately adjacent nucleotides are either A or T) and AT₀ (the number of SNPs whose immediately adjacent nucleotides are either C or G) for all 6 types of SNPs. C/G SNP of O. sativa and Arabidopsis has the highest AT₂/AT₀, which denotes C/G SNP may be influenced by the adjacent A&T nucleotides mostly. For SNPs in O. sativa, the neighboring effect of A&T nucleotides is limited to 2 nucleotides on both sides; for SNPs in Arabidopsis, the effect extends no more than 4 nucleotides on both sides.     

Plastid Sequence Evolution: A New Pattern of Nucleotide Substitutions in the Cucurbitaceae

Nucleotide substitutions (i.e., point mutations) are the primary driving force in generating DNA variation upon which selection can act. Substitutions called transitions, which entail exchanges between purines (A=adenine, G=guanine) or pyrimidines (C=cytosine, T=thymine), typically outnumber transversions (e.g., exchanges between a purine and a pyrimidine) in a DNA strand. With an increasing number of plant studies revealing a transversion rather than transition bias, we chose to perform a detailed substitution analysis for the plant family Cucurbitaceae using data from several short plastid DNA sequences. We generated a phylogenetic tree for 19 taxa of the tribe Benincaseae and related genera and then scored conservative substitution changes (e.g., those not exhibiting homoplasy or reversals) from the unambiguous branches of the tree. Neither the transition nor (A+T)/(G+C) biases found in previous studies were supported by our overall data. More importantly, we found a novel and symmetrical substitution bias in which Gs had been preferentially replaced by A, As by C, Cs by T, and Ts by G, resulting in the GACTG substitution series. Understanding this pattern will lead to new hypotheses concerning plastid evolution, which in turn will affect the choices of substitution models and other tree-building algorithms for phylogenetic analyses based on nucleotide data.     

Enhanced Synonymous Site Divergence in Positively Selected VertebrateAntimicrobial Peptide Genes

Nonrandom patterns associated with adaptively evolving genes can shed light on how selection and mutation produce rapid changes in sequences. I examine such patterns in two independent families of antimicrobial peptide genes: those in frogs, which are known to have evolved under positive selection, and those in flatfishes, which I show have also evolved under positive selection. I address two recently proposed hypotheses about the molecular evolution of antimicrobial peptide genes. The first is that the mature peptide region is replicated by an error-prone polymerase that increases the mutation rate and the transversion/transition ratio compared to the signal sequence of the same genes. The second is that mature peptides evolve in a coordinated fashion with their propieces, such that a change in net charge in one molecular region prompts an opposite change in charge in the other region. I test these hypotheses using alternative methods that minimize alignment errors, correct for phylogenetic nonindependence, reduce sequence saturation, and account for differing selection pressures on different regions of the gene. In both gene families I show that divergence at both synonymous and nonsynonymous sites within the mature peptide region is enhanced. However, in neither gene family is there evidence of an increased mutational transversion/transition ratio or coordinated evolution. My observations are consistent with either an elevated mutation rate in an adaptively evolving gene region or widespread selection on “silent” sites. These hypotheses challenge the assumption that mutations are random and can be measured by the synonymous substitution rate. [Reviewing Editor: Dr. Willie J. Swanson]     

Synonymous substitution rates in Drosophila: Mitochondrial versus nuclear genes

Synonymous substitution rates in mitochondrial and nuclear genes of Drosophila were compared. To make accurate comparisons, we considered the following: (1) relative synonymous rates, which do not require divergence time estimates, should be used; (2) methods estimating divergence should take into account base composition; (3) only very closely related species should be used to avoid effects of saturation; (4) the heterogeneity of rates should be examined. We modified the methods estimating synonymous substitution numbers to account for base composition bias. By using these methods, we found that mitochondrial genes have 1.7–3.4 times higher synonymous substitution rates than the fastest nuclear genes or 4.5–9.0 times higher rates than the average nuclear genes. The average rate of synonymous transversions was 2.7 (estimated from the melanogaster species subgroup) or 2.9 (estimated from the obscura group) times higher in mitochondrial genes than in nuclear genes. Synonymous transversions in mitochondrial genes occurred at an approximately equivalent rate to those in the fastest nuclear genes. This last result is not consistent with the hypothesis that the difference in turnover rates between mitochondrial and nuclear genomes is the major factor determining higher synonymous substitution rates in mtDNA. We conclude that the difference in synonymous substitution rates is due to a combination of two factors: a higher transitional mutation rate in mtDNA and constraints on nuclear genes due to selection for codon usage. Received: 27 November 1996 / Accepted: 8 May 1997     

Accelerated Evolution of Cytochrome b in Simian Primates: Adaptive Evolution in Concert with Other Mitochondrial Proteins?

We have sequenced the cytochrome b gene of Horsfield's tarsier, Tarsius bancanus, to complete a data set of sequences for this gene from representatives of each primate infraorder. These primate cytochrome b sequences were combined with those from representatives of three other mammalian orders (cat, whale, and rat) in an analysis of relative evolutionary rates. The nonsynonymous nucleotide substitution rate of the cytochrome b gene has increased approximately twofold along lineages leading to simian primates compared to that of the tarsier and other primate and nonprimate mammalian species. However, the rate of transversional substitutions at fourfold degenerate sites has remained uniform among all lineages. This increase in the evolutionary rate of cytochrome b is similar in character and magnitude to that described previously for the cytochrome c oxidase subunit II gene. We propose that the evolutionary rate increase observed for cytochrome b and cytochrome c oxidase subunit II may underlie an episode of coadaptive evolution of these two proteins in the mitochondria of simian primates. Received: 15 December 1997 / Accepted: 24 February 1998     

Unbiased estimation of the rates of synonymous and nonsynonymous substitution

Summary The current convention in estimating the number of substitutions per synonymous site (K _S ) and per nonsynonymous site (K _A ) between two protein-coding genes is to count each twofold degenerate site as one-third synonymous and two-thirds nonsynonymous because one of the three possible changes at such a site is synonymous and the other two are nonsynonymous. This counting rule can considerably overestimate theK _S value because transitional mutations tend to occur more often than transversional mutations and because most transitional mutations at twofold degenerate sites are synonymous. A new method that gives unbiased estimates is proposed. An application of the new and the old method to 14 pairs of mouse and rat genes shows that the new method gives aK _S value very close to the number of substitutions per fourfold degenerate site whereas the old method gives a value 30% higher. Both methods give aK _A value close to the number of substitutions per nondegenerate site.