GeneEvolve: a fast and memory efficient forward-time simulator of realistic whole-genome sequence and SNP data.

MOTIVATION: Computer simulations are excellent tools for understanding the evolutionary and genetic consequences of complex processes that cannot be analytically predicted and for creating realistic genetic data. There are many software packages that simulate genetic data, but they are typically not fast or memory efficient enough to simulate realistic, individual-level genome-wide SNP/sequence data. RESULTS: GeneEvolve is a user-friendly and efficient population genetics simulator that handles complex evolutionary and life history scenarios and generates individual-level phenotypes and realistic whole-genome sequence or SNP data. GeneEvolve runs forward-in-time, which allows it to provide a wide range of scenarios for mating systems, selection, population size and structure, migration, recombination and environmental effects. The software is designed to use as input data from real or previously simulated phased haplotypes, allowing it to mimic very closely the properties of real genomic data. AVAILABILITY AND IMPLEMENTATION: GeneEvolve is freely available at https://github.com/rtahmasbi/GeneEvolve CONTACT: Rasool.Tahmasbi@Colorado.eduSupplementary information: Supplementary data are available at Bioinformatics online.

Narrow-sense heritability estimation of complex traits using identity-by-descent information.

Heritability is a fundamental parameter in genetics. Traditional estimates based on family or twin studies can be biased due to shared environmental or non-additive genetic variance. Alternatively, those based on genotyped or imputed variants typically underestimate narrow-sense heritability contributed by rare or otherwise poorly tagged causal variants. Identical-by-descent (IBD) segments of the genome share all variants between pairs of chromosomes except new mutations that have arisen since the last common ancestor. Therefore, relating phenotypic similarity to degree of IBD sharing among classically unrelated individuals is an appealing approach to estimating the near full additive genetic variance while possibly avoiding biases that can occur when modeling close relatives. We applied an IBD-based approach (GREML-IBD) to estimate heritability in unrelated individuals using phenotypic simulation with thousands of whole-genome sequences across a range of stratification, polygenicity levels, and the minor allele frequencies of causal variants (CVs). In simulations, the IBD-based approach produced unbiased heritability estimates, even when CVs were extremely rare, although precision was low. However, population stratification and non-genetic familial environmental effects shared across generations led to strong biases in IBD-based heritability. We used data on two traits in ~120,000 people from the UK Biobank to demonstrate that, depending on the trait and possible confounding environmental effects, GREML-IBD can be applied to very large genetic datasets to infer the contribution of very rare variants lost using other methods. However, we observed apparent biases in these real data, suggesting that more work may be required to understand and mitigate factors that influence IBD-based heritability estimates.

Multiomics Screening Identifies Molecular Biomarkers Causally Associated With the Risk of Coronary Artery Disease.

BACKGROUND: In this study, we aimed to investigate functional mechanisms underlying coronary artery disease (CAD) loci and find molecular biomarkers for CAD. METHODS: We devised a multiomics data analysis approach based on Mendelian randomization and utilized it to search for molecular biomarkers causally associated with the risk of CAD within genomic regions known to be associated with CAD. RESULTS: Through our CAD-centered multiomics data analysis approach, we identified 33 molecular biomarkers (probes) that were causally associated with the risk of CAD. The majority of these (N=19) were methylation probes; moreover, methylation was often behind the causal effect of expression/protein probes. We identified a number of novel loci that have a causal impact on CAD including C5orf38, SF3A3, DHX36, and MRPL33. Furthermore, by integrating the risk factors of CAD in our analysis, we were able to investigate the clinical pathways whereby several of our probes exert their effect. We found that the SELE protein level in the blood is under the trans-regulatory impact of methylation sites within the ABO gene and that SELE exerts its effect on CAD through immune, glycemic, and lipid metabolism, making it a candidate of interest for therapeutic interventions. We found the methylation site, cg05126514 within the BSN gene exert its effect on CAD through central nervous system-lifestyle risk factors. Finally, genes with a transcriptional regulatory role (SF3A3, ILF3, and N4BP2L2) exert their effect on CAD through height. CONCLUSIONS: We demonstrate that multiomics data analysis is a powerful approach to unravel the functional mechanisms underlying CAD loci and to identify novel molecular biomarkers. Our results indicate epigenetic modifications are important in the pathogenesis of CAD and identifying and targeting these sites is of potential therapeutic interest to address the detrimental effects of both environmental and genetic factors.

Comparison of methods that use whole genome data to estimate the heritability and genetic architecture of complex traits.

Multiple methods have been developed to estimate narrow-sense heritability, h2, using single nucleotide polymorphisms (SNPs) in unrelated individuals. However, a comprehensive evaluation of these methods has not yet been performed, leading to confusion and discrepancy in the literature. We present the most thorough and realistic comparison of these methods to date. We used thousands of real whole-genome sequences to simulate phenotypes under varying genetic architectures and confounding variables, and we used array, imputed, or whole genome sequence SNPs to obtain 'SNP-heritability' estimates. We show that SNP-heritability can be highly sensitive to assumptions about the frequencies, effect sizes, and levels of linkage disequilibrium of underlying causal variants, but that methods that bin SNPs according to minor allele frequency and linkage disequilibrium are less sensitive to these assumptions across a wide range of genetic architectures and possible confounding factors. These findings provide guidance for best practices and proper interpretation of published estimates.