Abstract:The dimension of the population genetics data produced by next-generation sequencing platforms is extremely high. However, the "intrinsic dimensionality" of sequence data, which determines the structure of populations, is much lower. This motivates us to use locally linear embedding (LLE) which projects high dimensional genomic data into low dimensional, neighborhood preserving embedding, as a general framework for population structure and historical inference. To facilitate application of the LLE to population genetic analysis, we systematically investigate several important properties of the LLE and reveal the connection between the LLE and principal component analysis (PCA). Identifying a set of markers and genomic regions which could be used for population structure analysis will provide invaluable information for population genetics and association studies. In addition to identifying the LLE-correlated or PCA-correlated structure informative marker, we have developed a new statistic that integrates genomic information content in a genomic region for collectively studying its association with the population structure and LASSO algorithm to search such regions across the genomes. We applied the developed methodologies to a low coverage pilot dataset in the 1000 Genomes Project and a PHASE III Mexico dataset of the HapMap. We observed that 25.1%, 44.9% and 21.4% of the common variants and 89.2%, 92.4% and 75.1% of the rare variants were the LLE-correlated markers in CEU, YRI and ASI, respectively. This showed that rare variants, which are often private to specific populations, have much higher power to identify population substructure than common variants. The preliminary results demonstrated that next generation sequencing offers a rich resources and LLE provide a powerful tool for population structure analysis.

LEA: An R package for landscape and ecological association studies

Regulatory landscape enrichment analysis (RLEA): a computational toolkit for non-coding variant enrichment and cell type prioritization

Evaluation of redundancy analysis to identify signatures of local adaptation

Testing for Associations between Loci and Environmental Gradients Using Latent Factor Mixed Models

RDAforest: identifying environmental drivers of polygenic adaptation.

LFMM 2: Fast and Accurate Inference of Gene-Environment Associations in Genome-Wide Studies

Netreg: Network-Regularized Linear Models for Biological Association Studies

PEAS V1.0: a Package for Elementary Analysis of SNP Data

Gjls2: an R Package for Generalized Joint Location and Scale Analysis in X-inclusive Genome-Wide Association Studies

Using Environmental Correlations to Identify Loci Underlying Local Adaptation

BPEC: An R Package for Bayesian Phylogeographic and Ecological Clustering

LAVAA: a lightweight association viewer across ailments

Combined Linkage Disequilibrium and Linkage Analysis (cLDLA): implementation of a powerful approach to identify the genetic basis of complex traits in a bioinformatics workflow.

A Unifying Model for the Analysis of Phenotypic, Genetic and Geographic Data

Leveraging massively parallel reporter assays for evolutionary questions

Landscape Genetics: From Classic Molecular Markers to Genomics

DISENTANGLING THE EFFECTS OF GEOGRAPHIC AND ECOLOGICAL ISOLATION ON GENETIC DIFFERENTIATION

An integrated approach to reduce the impact of minor allele frequency and linkage disequilibrium on variable importance measures for genome-wide data

Manifold Learning for Human Population Structure Studies

slimr: An R package for tailor-made integrations of data in population genomic simulations over space and time

BoundaryStats: An R package to calculate boundary overlap statistics