Data from the 1000 Genomes project is quite often used as a reference for human genomic analysis. However, its accuracy needs to be assessed to understand the quality of predictions made using this ...reference. We present here an assessment of the genotyping, phasing, and imputation accuracy data in the 1000 Genomes project. We compare the phased haplotype calls from the 1000 Genomes project to experimentally phased haplotypes for 28 of the same individuals sequenced using the 10X Genomics platform.
We observe that phasing and imputation for rare variants are unreliable, which likely reflects the limited sample size of the 1000 Genomes project data. Further, it appears that using a population specific reference panel does not improve the accuracy of imputation over using the entire 1000 Genomes data set as a reference panel. We also note that the error rates and trends depend on the choice of definition of error, and hence any error reporting needs to take these definitions into account.
The quality of the 1000 Genomes data needs to be considered while using this database for further studies. This work presents an analysis that can be used for these assessments.
Full text
Available for:
DOBA, IZUM, KILJ, NUK, PILJ, PNG, SAZU, SIK, UILJ, UKNU, UL, UM, UPUK
Personalized medical care focuses on prediction of disease risk and response to medications. To build the risk models, access to both large-scale genomic resources and human genetic studies is ...required. The Taiwan Biobank (TWB) has generated high-coverage, whole-genome sequencing data from 1492 individuals and genome-wide SNP data from 103,106 individuals of Han Chinese ancestry using custom SNP arrays. Principal components analysis of the genotyping data showed that the full range of Han Chinese genetic variation was found in the cohort. The arrays also include thousands of known functional variants, allowing for simultaneous ascertainment of Mendelian disease-causing mutations and variants that affect drug metabolism. We found that 21.2% of the population are mutation carriers of autosomal recessive diseases, 3.1% have mutations in cancer-predisposing genes, and 87.3% carry variants that affect drug response. We highlight how TWB data provide insight into both population history and disease burden, while showing how widespread genetic testing can be used to improve clinical care.
Simulation is a key tool in population genetics for both methods development and empirical research, but producing simulations that recapitulate the main features of genomic datasets remains a major ...obstacle. Today, more realistic simulations are possible thanks to large increases in the quantity and quality of available genetic data, and the sophistication of inference and simulation software. However, implementing these simulations still requires substantial time and specialized knowledge. These challenges are especially pronounced for simulating genomes for species that are not well-studied, since it is not always clear what information is required to produce simulations with a level of realism sufficient to confidently answer a given question. The community-developed framework stdpopsim seeks to lower this barrier by facilitating the simulation of complex population genetic models using up-to-date information. The initial version of stdpopsim focused on establishing this framework using six well-characterized model species (Adrion et al., 2020). Here, we report on major improvements made in the new release of stdpopsim (version 0.2), which includes a significant expansion of the species catalog and substantial additions to simulation capabilities. Features added to improve the realism of the simulated genomes include non-crossover recombination and provision of species-specific genomic annotations. Through community-driven efforts, we expanded the number of species in the catalog more than threefold and broadened coverage across the tree of life. During the process of expanding the catalog, we have identified common sticking points and developed the best practices for setting up genome-scale simulations. We describe the input data required for generating a realistic simulation, suggest good practices for obtaining the relevant information from the literature, and discuss common pitfalls and major considerations. These improvements to stdpopsim aim to further promote the use of realistic whole-genome population genetic simulations, especially in non-model organisms, making them available, transparent, and accessible to everyone.
We analyze the role of solvation for enzymatic catalysis in two distinct, artificially designed Kemp Eliminases, KE07 and KE70, and mutated variants that were optimized by laboratory directed ...evolution. Using a spatially resolved analysis of hydration patterns, intermolecular vibrations, and local solvent entropies, we identify distinct classes of hydration water and follow their changes upon substrate binding and transition state formation for the designed KE07 and KE70 enzymes and their evolved variants. We observe that differences in hydration of the enzymatic systems are concentrated in the active site and undergo significant changes during substrate recruitment. For KE07, directed evolution reduces variations in the hydration of the polar catalytic center upon substrate binding, preserving strong protein–water interactions, while the evolved enzyme variant of KE70 features a more hydrophobic reaction center for which the expulsion of low-entropy water molecules upon substrate binding is substantially enhanced. While our analysis indicates a system-dependent role of solvation for the substrate binding process, we identify more subtle changes in solvation for the transition state formation, which are less affected by directed evolution.
Full text
Available for:
IJS, KILJ, NUK, PNG, UL, UM
We have used the AMOEBA model to simulate the THz spectra of two zwitterionic amino acids in aqueous solution, which is compared to the results on these same systems using ab initio molecular ...dynamics (AIMD) simulations. Overall we find that the polarizable force field shows promising agreement with AIMD data for both glycine and valine in water. This includes the THz spectral assignments and the mode-specific spectral decomposition into intramolecular solute motions as well as distinct solute-water cross-correlation modes some of which cannot be captured by non-polarizable force fields that rely on fixed partial charges. This bodes well for future studies for simulating and decomposing the THz spectra for larger solutes such as proteins or polymers for which AIMD studies are presently intractable. Furthermore, we believe that the current study on rather simple aqueous solutions offers a way to systematically investigate the importance of charge transfer, nuclear quantum effects, and the validity of computationally practical density functionals, all of which are needed to fully quantitatively capture complex dynamical motions in the condensed phase.
This work reports on the role of solvation for enzymatic catalysis in two distinct, artificially designed Kemp Eliminases, KE07 and KE70, and mutated variants that were optimized by laboratory ...directed evolution. Using a spatially resolved analysis of hydration patterns, intermolecular vibrations, and local solvent entropies, we identify distinct classes of hydration water and follow their changes upon substrate binding and transition state formation for the designed KE07 and KE70 enzymes and their evolved variants. We observe that differences in hydration of the enzymatic systems are concentrated in the active site and undergo significant changes during substrate recruitment. For KE07, directed evolution reduces variations in the hydration of the polar catalytic center upon substrate binding, preserving strong protein-water interactions, while the evolved enzyme variant of KE70 features a more hydrophobic reaction center for which the expulsion of low-entropy water molecules upon substrate binding is substantially enhanced. While our analysis suggests a system-dependent role of solvation for the substrate binding process, we identify more subtle changes in solvation for the transition state formation, which are less affected by directed evolution.
Full text
Available for:
IJS, KILJ, NUK, PNG, UL, UM
We have used the AMOEBA model to simulate the THz spectra of two zwitterionic amino acids in aqueous solution, which is compared to the results on these same systems using
ab initio
molecular ...dynamics (AIMD) simulations. Overall we find that the polarizable force field shows promising agreement with AIMD data for both glycine and valine in water. This includes the THz spectral assignments and the mode-specific spectral decomposition into intramolecular solute motions as well as distinct solute-water cross-correlation modes some of which cannot be captured by non-polarizable force fields that rely on fixed partial charges. This bodes well for future studies for simulating and decomposing the THz spectra for larger solutes such as proteins or polymers for which AIMD studies are presently intractable. Furthermore, we believe that the current study on rather simple aqueous solutions offers a way to systematically investigate the importance of charge transfer, nuclear quantum effects, and the validity of computationally practical density functionals, all of which are needed to fully quantitatively capture complex dynamical motions in the condensed phase.
We have used the AMOEBA model to simulate the THz spectra of two zwitterionic amino acids in aqueous solution, which is compared to the results on these same systems using
ab initio
molecular dynamics (AIMD) simulations.
Baboons (genus
) are broadly studied in the wild and in captivity. They are widely used as a nonhuman primate model for biomedical studies, and the Southwest National Primate Research Center (SNPRC) ...at Texas Biomedical Research Institute has maintained a large captive baboon colony for more than 50 yr. Unlike other model organisms, however, the genomic resources for baboons are severely lacking. This has hindered the progress of studies using baboons as a model for basic biology or human disease. Here, we describe a data set of 100 high-coverage whole-genome sequences obtained from the mixed colony of olive (
) and yellow (
) baboons housed at the SNPRC. These data provide a comprehensive catalog of common genetic variation in baboons, as well as a fine-scale genetic map. We show how the data can be used to learn about ancestry and admixture and to correct errors in the colony records. Finally, we investigated the consequences of inbreeding within the SNPRC colony and found clear evidence for increased rates of infant mortality and increased homozygosity of putatively deleterious alleles in inbred individuals.
PDF
