Presence Of Population Substructure Research Articles

Development of novel microsatellite markers for population differentiation and detection of natural selection in wild populations of butter catfish, Ompok bimaculatus (Bloch, 1794).

Butter catfish (Ompok bimaculatus) is a preferred species in South East Asia, with huge aquaculture potential. However, there is limited information about genetic stock composition due to insufficient markers. The goal of this study was to develop de novo microsatellite markers. For sequencing, genomic SMRT bell libraries (1.5 Kbp size) were prepared for O. bimaculatus. A total of 114 SSR containing sequences were used for primer designing. Polymorphic loci were validated by genotyping 83 individuals from four distant riverine populations, viz., Brahmaputra, Bichiya, Gomti and Kaveri. A total of 30 microsatellite loci were polymorphic, of which five were found to be associated with functional genes and eight (four positive and four negative) loci were found to be under selection pressure. A total of 115 alleles were detected in all loci and PIC ranged from 0.539 to 0.927 and pair-wise FST values from 0.1267 to 0.26002 (p < 0.001), with an overall FST value of 0.17047, indicating the presence of population sub-structure. Cross-species transferability of 29 loci (96.67%) was successful in congener species, Ompok pabda. The novel SSR markers developed in this study would facilitate stock characterization of natural populations, to be used in future selection breeding programs and planning conservation strategies in these species. Identified non-neutral markers will give insights into the effect of local adaptation on genetic differentiation in the natural population of this species.

Estimating SNP heritability in presence of population substructure in biobank-scale datasets.

Single nucleotide polymorphism heritability of a trait is measured as the proportion of total variance explained by the additive effects of genome-wide single nucleotide polymorphisms. Linear mixed models are routinely used to estimate single nucleotide polymorphism heritability for many complex traits, which requires estimation of a genetic relationship matrix among individuals. Heritability is usually estimated by the restricted maximum likelihood or method of moments approaches such as Haseman-Elston regression. The common practice of accounting for such population substructure is to adjust for the top few principal components of the genetic relationship matrix as covariates in the linear mixed model. This can get computationally very intensive on large biobank-scale datasets. Here, we propose a method of moments approach for estimating single nucleotide polymorphism heritability in presence of population substructure. Our proposed method is computationally scalable on biobank datasets and gives an asymptotically unbiased estimate of heritability in presence of discrete substructures. It introduces the adjustments for population stratification in a second-order estimating equation. It allows these substructures to vary in their single nucleotide polymorphism allele frequencies and in their trait distributions (means and variances) while the heritability is assumed to be the same across these substructures. Through extensive simulation studies and the application on 7 quantitative traits in the UK Biobank cohort, we demonstrate that our proposed method performs well in the presence of population substructure and much more computationally efficient than existing approaches.

On the Recombination Rate Estimation in the Presence of Population Substructure.

As recombination events are not uniformly distributed along the human genome, the estimation of fine-scale recombination maps, e.g. HapMap Project, has been one of the major research endeavors over the last couple of years. For simulation studies, these estimates provide realistic reference scenarios to design future study and to develop novel methodology. To achieve a feasible framework for the estimation of such recombination maps, existing methodology uses sample probabilities for a two-locus model with recombination, with recent advances allowing for computationally fast implementations. In this work, we extend the existing theoretical framework for the recombination rate estimation to the presence of population substructure. We show under which assumptions the existing methodology can still be applied. We illustrate our extension of the methodology by an extensive simulation study.

Family-based association analysis: a fast and efficient method of multivariate association analysis with multiple variants.

BackgroundMany disease phenotypes are outcomes of the complicated interplay between multiple genes, and multiple phenotypes are affected by a single or multiple genotypes. Therefore, joint analysis of multiple phenotypes and multiple markers has been considered as an efficient strategy for genome-wide association analysis, and in this work we propose an omnibus family-based association test for the joint analysis of multiple genotypes and multiple phenotypes.ResultsThe proposed test can be applied for both quantitative and dichotomous phenotypes, and it is robust under the presence of population substructure, as long as large-scale genomic data is available. Using simulated data, we showed that our method is statistically more efficient than the existing methods, and the practical relevance is illustrated by application of the approach to obesity-related phenotypes.ConclusionsThe proposed method may be more statistically efficient than the existing methods. The application was developed in C++ and is available at the following URL: http://healthstat.snu.ac.kr/software/mfqls/.

Impact of Population Stratification on Family-Based Association in an Admixed Population.

Population substructure is a well-known confounder in population-based case-control genetic studies, but its impact in family-based studies is unclear. We performed population substructure analysis using extended families of admixed population to evaluate power and Type I error in an association study framework. Our analysis shows that power was improved by 1.5% after principal components adjustment. Type I error was also reduced by 2.2% after adjusting for family substratification. The presence of population substructure was underscored by discriminant analysis, in which over 92% of individuals were correctly assigned to their actual family using only 100 principal components. This study demonstrates the importance of adjusting for population substructure in family-based studies of admixed populations.

Robust analysis with related samples under the presence of population substructure and its application to body mass index

Recent investigations such as a more powerful quasi-likelihoods score test (MQLS) statistic have enabled the efficient association analysis with related samples. Although those approaches are robust against the mis-specified phenotypic distribution and covariance structure, it has been shown that MQLS statistic becomes violated under the presence of the population substructure if the level of population substructure depends on the genomic location. In this report, we propose a new statistical method which combines EIGENSTRAT approach and MQLS-statistic. The proposed method was evaluated with simulation data under various scenarios and we found that proposed method performs better than the traditional methods such as transmission disequilibrium test. The proposed method was applied to genetic association analysis for body mass index with Framingham heart study, and we found that rs1121980 and rs9940128 in the linkage block in FTO gene are associated with the body mass index.

FamLBL: detecting rare haplotype disease association based on common SNPs using case-parent triads.

In recent years, there has been an increasing interest in using common single-nucleotide polymorphisms (SNPs) amassed in genome-wide association studies to investigate rare haplotype effects on complex diseases. Evidence has suggested that rare haplotypes may tag rare causal single-nucleotide variants, making SNP-based rare haplotype analysis not only cost effective, but also more valuable for detecting causal variants. Although a number of methods for detecting rare haplotype association have been proposed in recent years, they are population based and thus susceptible to population stratification. We propose family-triad-based logistic Bayesian Lasso (famLBL) for estimating effects of haplotypes on complex diseases using SNP data. By choosing appropriate prior distribution, effect sizes of unassociated haplotypes can be shrunk toward zero, allowing for more precise estimation of associated haplotypes, especially those that are rare, thereby achieving greater detection power. We evaluate famLBL using simulation to gauge its type I error and power. Compared with its population counterpart, LBL, highlights famLBL's robustness property in the presence of population substructure. Further investigation by comparing famLBL with Family-Based Association Test (FBAT) reveals its advantage for detecting rare haplotype association. famLBL is implemented as an R-package available at http://www.stat.osu.edu/∼statgen/SOFTWARE/LBL/.

Impact of Population Substructure on Trend Tests for Genetic Case–Control Association Studies

Hidden population substructure in case-control data has the potential to distort the performance of Cochran-Armitage trend tests (CATTs) for genetic associations. Three possible scenarios that may arise are investigated here: (i) heterogeneity of genotype frequencies across unidentified subpopulations (PSI), (ii) heterogeneity of genotype frequencies and disease risk across unidentified subpopulations (PSII), and (iii) cryptic correlations within unidentified subpopulations. A unified approach is presented for deriving the bias and variance distortion under the three scenarios for any CATT in a general family. Using these analytical formulas, we evaluate the excess type I errors of the CATTs numerically in the presence of population substructure. Our results provide insight into the properties of some proposed corrections for bias and variance distortion and show why they may not fully correct for the effects of population substructure.

Gene-environment interaction tests for family studies with quantitative phenotypes: A review and extension to longitudinal measures

Longitudinal studies are an important tool for analysing traits that change over time, depending on individual characteristics and environmental exposures. Complex quantitative traits, such as lung function, may change over time and appear to depend on genetic and environmental factors, as well as on potential gene-environment interactions. There is a growing interest in modelling both marginal genetic effects and gene-environment interactions. In an admixed population, the use of traditional statistical models may fail to adjust for confounding by ethnicity, leading to bias in the genetic effect estimates. A variety of methods have been developed to account for the genetic substructure of human populations. Family-based designs provide an important resource for avoiding confounding due to admixture. To date, however, most genetic analyses have been applied to cross-sectional designs. In this paper, we propose a methodology which aims to improve the assessment of main genetic effect and gene-environment interaction effects by combining the advantages of both longitudinal studies for continuous phenotypes, and the family-based designs. This approach is based on an extension of ordinary linear mixed models for quantitative phenotypes, which incorporates information from a case-parent design. Our results indicate that use of this method allows both main genetic and gene-environment interaction effects to be estimated without bias, even in the presence of population substructure.

Correcting for cryptic relatedness by a regression-based genomic control method

BackgroundGenomic control (GC) method is a useful tool to correct for the cryptic relatedness in population-based association studies. It was originally proposed for correcting for the variance inflation of Cochran-Armitage's additive trend test by using information from unlinked null markers, and was later generalized to be applicable to other tests with the additional requirement that the null markers are matched with the candidate marker in allele frequencies. However, matching allele frequencies limits the number of available null markers and thus limits the applicability of the GC method. On the other hand, errors in genotype/allele frequencies may cause further bias and variance inflation and thereby aggravate the effect of GC correction.ResultsIn this paper, we propose a regression-based GC method using null markers that are not necessarily matched in allele frequencies with the candidate marker. Variation of allele frequencies of the null markers is adjusted by a regression method.ConclusionThe proposed method can be readily applied to the Cochran-Armitage's trend tests other than the additive trend test, the Pearson's chi-square test and other robust efficiency tests. Simulation results show that the proposed method is effective in controlling type I error in the presence of population substructure.

The expected performance of single nucleotide polymorphism loci in paternity testing

We discuss the utility of single nucleotide polymorphism loci for full trio and mother-unavailable paternity testing cases, in the presence of population substructure and relatedness of putative and actual fathers. We focus primarily on the expected number of loci required to gain specified probabilities of mismatches, and report the expected proportion of paternity indices greater than three threshold values for these loci.

Calculating the exclusion probability and paternity index for X-chromosomal loci in the presence of substructure

There has been recent interest in the use of X-chromosomal loci for forensic and relatedness testing casework, with many authors developing new X-linked short tandem repeat (STR) loci suitable for forensic use. Here we present formulae for two key quantities in paternity testing, the average probability of exclusion and the paternity index, which are suitable for X-chromosomal loci in the presence of population substructure.

Family-based association studies.

Over the past decade, attention has turned from positional cloning of Mendelian disease genes to the dissection of complex diseases. Both theoretical and empirical studies have shown that traditional linkage studies may be inferior in power compared to studies that directly utilize allele status. Case-control association studies, as an alternative, are subject to bias due to population stratification. As a compromise between linkage studies and case-control studies, family-based association designs have received great attention recently due to their potentially higher power to identify complex disease genes and their robustness in the presence of population substructure. In this review, we first describe the basic family-based association design involving one affected offspring with its two parents, all genotyped for a biallelic genetic marker. Extensions of the original transmission disequilibrium tests to multiallelic markers, families with multiple siblings, families with incomplete parental genotypes, and general pedigree structures are discussed. Further developments of statistical methods to study quantitative traits, to analyse genes on the X chromosome, to incorporate multiple tightly linked markers, to identify imprinting genes, and to detect gene-environment interactions are also reviewed. Finally, we discuss the implications of the completion of the Human Genome Project and the identification of hundreds of thousands of genetic polymorphisms on employing family-based association designs to search for complex disease genes.

TDT statistics for mapping quantitative trait loci.

The original transmission disequilibrium test (TDT), was introduced to test for linkage between a marker and a disease-susceptibility locus (Spielman et al. 1993). Allison (1997) extended the TDT procedure to quantitative traits. Allison's test, however, is restrictive in that it requires family trios consisting of one heterozygous parent, one homozygous parent and one child, and considers only the situation where the marker locus is analogous to the quantitative trait locus itself. In this paper, we propose, investigate and apply a general TDT for quantitative traits that permits more than one child per family, does not require only one parent to be heterozygous, and allows for the fact that the various alleles at the marker and trait loci may be at varying degree of linkage disequilibrium. We also show that this TDT for quantitative traits is still a valid test of linkage in the presence of population substructure. To provide guidelines for study design, we develop analytic formulae for calculation of the power of the TDT for mapping quantitative trait loci and investigate the impact of various factors on the power. Power calculations show that the proposed TDT for quantitative traits is more powerful than Allison's basic test statistic and the extreme discordant sib pair linkage method. The proposed TDT statistic for quantitative traits is applied to systolic blood pressure variation in the Rochester Family Heart Study using an extremely discordant sibling pair design.

A test of allelic independence based on distributions of allele size differences at microsatellite loci.

Thousands of polymorphic microsatellite loci are now mapped in the human genome, most of which exhibit a large number of segregating alleles. For forensic, evolutionary and gene mapping applications, it is important to establish intralocus allelic independence at these loci. We develop a test for intralocus allelic independence based on repeat size differences between alleles within genotypes of individuals, and show its relationship with the intraclass correlation of allele sizes within loci. Applications of this test to 13 short tandem repeat loci in 4 subpopulations indicate that the distribution of size differences of alleles has adequate power to detect intralocus dependence of alleles caused by population substructure. In addition, size differences between randomly chosen alleles provide information regarding heterozygosity and contiguity of allele sizes at the locus, which by themselves do not necessarily indicate the presence of population substructure.

Potassium Perchlorate in Thyrotoxicosis

Abstract SNP heritability of a trait is measured by the proportion of total variance explained by the additive effects of genome-wide single nucleotide polymorphisms (SNPs). Linear mixed models are routinely used to estimate SNP heritability for many complex traits. The basic concept behind this approach is to model genetic contribution as a random effect, where the variance of this genetic contribution attributes to the heritability of the trait. This linear mixed model approach requires estimation of ‘relatedness’ among individuals in the sample, which is usually captured by estimating a genetic relationship matrix (GRM). Heritability is estimated by the restricted maximum likelihood (REML) or method of moments (MOM) approaches, and this estimation relies heavily on the GRM computed from the genetic data on individuals. Presence of population substructure in the data could significantly impact the GRM estimation and may introduce bias in heritability estimation. The common practice of accounting for such population substructure is to adjust for the top few principal components of the GRM as covariates in the linear mixed model. Here we propose an alternative way of estimating heritability in multi-ethnic studies. Our proposed approach is a MOM estimator derived from the Haseman-Elston regression and gives an asymptotically unbiased estimate of heritability in presence of population stratification. It introduces adjustments for the population stratification in a second-order estimating equation and allows for the total phenotypic variance vary by ethnicity. We study the performance of different MOM and REML approaches in presence of population stratification through extensive simulation studies. We estimate the heritability of height, weight and other anthropometric traits in the UK Biobank cohort to investigate the impact of subtle population substructure on SNP heritability estimation.