A powerful approach to identify replicable variants in genome-wide association studies.

Replicability is the cornerstone of modern scientific research. Reliable identifications of genotype-phenotype associations that are significant in multiple genome-wide association studies (GWASs) provide stronger evidence for the findings. Current replicability analysis relies on the independence assumption among single-nucleotide polymorphisms (SNPs) and ignores the linkage disequilibrium (LD) structure. We show that such a strategy may produce either overly liberal or overly conservative results in practice. We develop an efficient method, ReAD, to detect replicable SNPs associated with the phenotype from two GWASs accounting for the LD structure. The local dependence structure of SNPs across two heterogeneous studies is captured by a four-state hidden Markov model (HMM) built on two sequences of p values. By incorporating information from adjacent locations via the HMM, our approach provides more accurate SNP significance rankings. ReAD is scalable, platform independent, and more powerful than existing replicability analysis methods with effective false discovery rate control. Through analysis of datasets from two asthma GWASs and two ulcerative colitis GWASs, we show that ReAD can identify replicable genetic loci that existing methods might otherwise miss.

mBAT-combo: A more powerful test to detect gene-trait associations from GWAS data.

Gene-based association tests aggregate multiple SNP-trait associations into sets defined by gene boundaries and are widely used in post-GWAS analysis. A common approach for gene-based tests is to combine SNPs associations by computing the sum of χ2 statistics. However, this strategy ignores the directions of SNP effects, which could result in a loss of power for SNPs with masking effects, e.g., when the product of two SNP effects and the linkage disequilibrium (LD) correlation is negative. Here, we introduce "mBAT-combo," a set-based test that is better powered than other methods to detect multi-SNP associations in the context of masking effects. We validate the method through simulations and applications to real data. We find that of 35 blood and urine biomarker traits in the UK Biobank, 34 traits show evidence for masking effects in a total of 4,273 gene-trait pairs, indicating that masking effects is common in complex traits. We further validate the improved power of our method in height, body mass index, and schizophrenia with different GWAS sample sizes and show that on average 95.7% of the genes detected only by mBAT-combo with smaller sample sizes can be identified by the single-SNP approach with a 1.7-fold increase in sample sizes. Eleven genes significant only in mBAT-combo for schizophrenia are confirmed by functionally informed fine-mapping or Mendelian randomization integrating gene expression data. The framework of mBAT-combo can be applied to any set of SNPs to refine trait-association signals hidden in genomic regions with complex LD structures.

Evidence of epistasis in regions of long-range linkage disequilibrium across five complex diseases in the UK Biobank and eMERGE datasets.

Leveraging linkage disequilibrium (LD) patterns as representative of population substructure enables the discovery of additive association signals in genome-wide association studies (GWASs). Standard GWASs are well-powered to interrogate additive models; however, new approaches are required for invesigating other modes of inheritance such as dominance and epistasis. Epistasis, or non-additive interaction between genes, exists across the genome but often goes undetected because of a lack of statistical power. Furthermore, the adoption of LD pruning as customary in standard GWASs excludes detection of sites that are in LD but might underlie the genetic architecture of complex traits. We hypothesize that uncovering long-range interactions between loci with strong LD due to epistatic selection can elucidate genetic mechanisms underlying common diseases. To investigate this hypothesis, we tested for associations between 23 common diseases and 5,625,845 epistatic SNP-SNP pairs (determined by Ohta's D statistics) in long-range LD (>0.25 cM). Across five disease phenotypes, we identified one significant and four near-significant associations that replicated in two large genotype-phenotype datasets (UK Biobank and eMERGE). The genes that were most likely involved in the replicated associations were (1) members of highly conserved gene families with complex roles in multiple pathways, (2) essential genes, and/or (3) genes that were associated in the literature with complex traits that display variable expressivity. These results support the highly pleiotropic and conserved nature of variants in long-range LD under epistatic selection. Our work supports the hypothesis that epistatic interactions regulate diverse clinical mechanisms and might especially be driving factors in conditions with a wide range of phenotypic outcomes.

LDER-GE estimates phenotypic variance component of gene-environment interactions in human complex traits accurately with GE interaction summary statistics and full LD information.

Gene-environment (GE) interactions are essential in understanding human complex traits. Identifying these interactions is necessary for deciphering the biological basis of such traits. In this study, we review state-of-art methods for estimating the proportion of phenotypic variance explained by genome-wide GE interactions and introduce a novel statistical method Linkage-Disequilibrium Eigenvalue Regression for Gene-Environment interactions (LDER-GE). LDER-GE improves the accuracy of estimating the phenotypic variance component explained by genome-wide GE interactions using large-scale biobank association summary statistics. LDER-GE leverages the complete Linkage Disequilibrium (LD) matrix, as opposed to only the diagonal squared LD matrix utilized by LDSC (Linkage Disequilibrium Score)-based methods. Our extensive simulation studies demonstrate that LDER-GE performs better than LDSC-based approaches by enhancing statistical efficiency by ~23%. This improvement is equivalent to a sample size increase of around 51%. Additionally, LDER-GE effectively controls type-I error rate and produces unbiased results. We conducted an analysis using UK Biobank data, comprising 307 259 unrelated European-Ancestry subjects and 966 766 variants, across 217 environmental covariate-phenotype (E-Y) pairs. LDER-GE identified 34 significant E-Y pairs while LDSC-based method only identified 23 significant E-Y pairs with 22 overlapped with LDER-GE. Furthermore, we employed LDER-GE to estimate the aggregated variance component attributed to multiple GE interactions, leading to an increase in the explained phenotypic variance with GE interactions compared to considering main genetic effects only. Our results suggest the importance of impacts of GE interactions on human complex traits.

Association Analysis Provides Insights into Plant Mitonuclear Interactions.

Cytonuclear interaction refers to the complex and ongoing process of coevolution between nuclear and organelle genomes, which are responsible for cellular respiration, photosynthesis, lipid metabolism, etc. and play a significant role in adaptation and speciation. There have been a large number of studies to detect signatures of cytonuclear interactions. However, identification of the specific nuclear and organelle genetic polymorphisms that are involved in these interactions within a species remains relatively rare. The recent surge in whole genome sequencing has provided us an opportunity to explore cytonuclear interaction from a population perspective. In this study, we analyzed a total of 3,439 genomes from 7 species to identify signals of cytonuclear interactions by association (linkage disequilibrium) analysis of variants in both the mitochondrial and nuclear genomes across flowering plants. We also investigated examples of nuclear loci identified based on these association signals using subcellular localization assays, gene editing, and transcriptome sequencing. Our study provides a novel perspective on the investigation of cytonuclear coevolution, thereby enriching our understanding of plant fitness and offspring sterility.

Correlated Allele Frequency Changes Reveal Clonal Structure and Selection in Temporal Genetic Data.

In evolving populations where the rate of beneficial mutations is large, subpopulations of individuals with competing beneficial mutations can be maintained over long times. Evolution with this kind of clonal structure is commonly observed in a wide range of microbial and viral populations. However, it can be difficult to completely resolve clonal dynamics in data. This is due to limited read lengths in high-throughput sequencing methods, which are often insufficient to directly measure linkage disequilibrium or determine clonal structure. Here, we develop a method to infer clonal structure using correlated allele frequency changes in time-series sequence data. Simulations show that our method recovers true, underlying clonal structures when they are known and accurately estimate linkage disequilibrium. This information can then be combined with other inference methods to improve estimates of the fitness effects of individual mutations. Applications to data suggest novel clonal structures in an E. coli long-term evolution experiment, and yield improved predictions of the effects of mutations on bacterial fitness and antibiotic resistance. Moreover, our method is computationally efficient, requiring orders of magnitude less run time for large data sets than existing methods. Overall, our method provides a powerful tool to infer clonal structures from data sets where only allele frequencies are available, which can also improve downstream analyses.

Inference of Host-Pathogen Interaction Matrices from Genome-Wide Polymorphism Data.

Host-pathogen coevolution is defined as the reciprocal evolutionary changes in both species due to genotype × genotype (G×G) interactions at the genetic level determining the outcome and severity of infection. While co-analyses of hosts and pathogen genomes (co-genome-wide association studies) allow us to pinpoint the interacting genes, these do not reveal which host genotype(s) is/are resistant to which pathogen genotype(s). The knowledge of this so-called infection matrix is important for agriculture and medicine. Building on established theories of host-pathogen interactions, we here derive four novel indices capturing the characteristics of the infection matrix. These indices can be computed from full genome polymorphism data of randomly sampled uninfected hosts, as well as infected hosts and their pathogen strains. We use these indices in an approximate Bayesian computation method to pinpoint loci with relevant G×G interactions and to infer their underlying interaction matrix. In a combined single nucleotide polymorphism dataset of 451 European humans and their infecting hepatitis C virus (HCV) strains and 503 uninfected individuals, we reveal a new human candidate gene for resistance to HCV and new virus mutations matching human genes. For two groups of significant human-HCV (G×G) associations, we infer a gene-for-gene infection matrix, which is commonly assumed to be typical of plant-pathogen interactions. Our model-based inference framework bridges theoretical models of G×G interactions with host and pathogen genomic data. It, therefore, paves the way for understanding the evolution of key G×G interactions underpinning HCV adaptation to the European human population after a recent expansion.

Partitioning gene-level contributions to complex-trait heritability by allele frequency identifies disease-relevant genes.

Recent works have shown that SNP heritability-which is dominated by low-effect common variants-may not be the most relevant quantity for localizing high-effect/critical disease genes. Here, we introduce methods to estimate the proportion of phenotypic variance explained by a given assignment of SNPs to a single gene ("gene-level heritability"). We partition gene-level heritability by minor allele frequency (MAF) to find genes whose gene-level heritability is explained exclusively by "low-frequency/rare" variants (0.5% ≤ MAF < 1%). Applying our method to ∼16K protein-coding genes and 25 quantitative traits in the UK Biobank (N = 290K "White British"), we find that, on average across traits, ∼2.5% of nonzero-heritability genes have a rare-variant component and only ∼0.8% (327 gene-trait pairs) have heritability exclusively from rare variants. Of these 327 gene-trait pairs, 114 (35%) were not detected by existing gene-level association testing methods. The additional genes we identify are significantly enriched for known disease genes, and we find several examples of genes that have been previously implicated in phenotypically related Mendelian disorders. Notably, the rare-variant component of gene-level heritability exhibits trends different from those of common-variant gene-level heritability. For example, while total gene-level heritability increases with gene length, the rare-variant component is significantly larger among shorter genes; the cumulative distributions of gene-level heritability also vary across traits and reveal differences in the relative contributions of rare/common variants to overall gene-level polygenicity. While nonzero gene-level heritability does not imply causality, if interpreted in the correct context, gene-level heritability can reveal useful insights into complex-trait genetic architecture.

Cost-effective genomic prediction of critical economic traits in sturgeons through low-coverage sequencing.

Low-coverage whole-genome sequencing (LCS) offers a cost-effective alternative for sturgeon breeding, especially given the lack of SNP chips and the high costs associated with whole-genome sequencing. In this study, the efficiency of LCS for genotype imputation and genomic prediction was assessed in 643 sequenced Russian sturgeons (∼13.68×). The results showed that using BaseVar+STITCH at a sequencing depth of 2× with a sample size larger than 300 resulted in the highest genotyping accuracy. In addition, when the sequencing depth reached 0.5× and SNP density was reduced to 50 K through linkage disequilibrium pruning, the prediction accuracy was comparable to that of whole sequencing depth. Furthermore, an incremental feature selection method has the potential to improve prediction accuracy. This study suggests that the combination of LCS and imputation can be a cost-effective strategy, contributing to the genetic improvement of economic traits and promoting genetic gains in aquaculture species.

Viral genomic variation and the severity of genital HSV-2 infection as quantified by shedding rate: a viral genome-wide association study.

BACKGROUND: The clinical severity of genital HSV-2 infection varies widely among infected persons with some experiencing frequent genital lesions while others are asymptomatic. The viral genital shedding rate is closely associated with and has been established as a surrogate marker of clinical severity. METHODS: To assess the relationship between viral genetics and shedding, we assembled a set of 145 persons who had the severity of their genital herpes quantified through determination of their HSV genital shedding rate. An HSV-2 sample from each person was sequenced and biallelic variants among these genomes were identified. RESULTS: We found no association between metrics of genome-wide variation in HSV-2 and shedding rate. A viral genome-wide association study (vGWAS) identified the minor alleles of three individual unlinked variants as significantly associated with higher shedding rate (p<8.4x10-5): C44973T (A512T), a non-synonymous variant in UL22 (glycoprotein H); A74534G, a synonymous variant in UL36 (large tegument protein); and T119283C, an intergenic variant. We also found an association between the total number of minor alleles for the significant variants and shedding rate (p=6.6x10-7). CONCLUSIONS: These results add to a growing body of literature for HSV suggesting a connection between viral genetic variation and clinically important phenotypes of infection.

Improved two-step testing of genome-wide gene-environment interactions.

Two-step tests for gene-environment ( G × E $G\times E$ ) interactions exploit marginal single-nucleotide polymorphism (SNP) effects to improve the power of a genome-wide interaction scan. They combine a screening step based on marginal effects used to "bin" SNPs for weighted hypothesis testing in the second step to deliver greater power over single-step tests while preserving the genome-wide Type I error. However, the presence of many SNPs with detectable marginal effects on the trait of interest can reduce power by "displacing" true interactions with weaker marginal effects and by adding to the number of tests that need to be corrected for multiple testing. We introduce a new significance-based allocation into bins for Step-2 G × E $G\times E$ testing that overcomes the displacement issue and propose a computationally efficient approach to account for multiple testing within bins. Simulation results demonstrate that these simple improvements can provide substantially greater power than current methods under several scenarios. An application to a multistudy collaboration for understanding colorectal cancer reveals a G × Sex interaction located near the SMAD7 gene.

Effective population size in field pea.

BACKGROUND: Effective population size (Ne) is a pivotal parameter in population genetics as it can provide information on the rate of inbreeding and the contemporary status of genetic diversity in breeding populations. The population with smaller Ne can lead to faster inbreeding, with little potential for genetic gain making selections ineffective. The importance of Ne has become increasingly recognized in plant breeding, which can help breeders monitor and enhance the genetic variability or redesign their selection protocols. Here, we present the first Ne estimates based on linkage disequilibrium (LD) in the pea genome. RESULTS: We calculated and compared Ne using SNP markers from North Dakota State University (NDSU) modern breeding lines and United States Department of Agriculture (USDA) diversity panel. The extent of LD was highly variable not only between populations but also among different regions and chromosomes of the genome. Overall, NDSU had a higher and longer-range LD than the USDA that could extend up to 500 Kb, with a genome-wide average r2 of 0.57 (vs 0.34), likely due to its lower recombination rates and the selection background. The estimated Ne for the USDA was nearly three-fold higher (Ne = 174) than NDSU (Ne = 64), which can be confounded by a high degree of population structure due to the selfing nature of pea. CONCLUSIONS: Our results provided insights into the genetic diversity of the germplasm studied, which can guide plant breeders to actively monitor Ne in successive cycles of breeding to sustain viability of the breeding efforts in the long term.

Genome-wide association scan reveals the reinforcing effect of nano-potassium in improving the yield and quality of salt-stressed barley via enhancing the antioxidant defense system.

Salinity is one of the major environmental factor that can greatly impact the growth, development, and productivity of barley. Our study aims to detect the natural phenotypic variation of morphological and physiological traits under both salinity and potassium nanoparticles (n-K) treatment. In addition to understanding the genetic basis of salt tolerance in barley is a critical aspect of plant breeding for stress resilience. Therefore, a foliar application of n-K was applied at the vegetative stage for 138 barley accessions to enhance salt stress resilience. Interestingly, barley accessions showed high significant increment under n-K treatment compared to saline soil. Based on genome-wide association studies (GWAS) analysis, causative alleles /reliable genomic regions were discovered underlying improved salt resilience through the application of potassium nanoparticles. On chromosome 2H, a highly significant QTN marker (A:C) was located at position 36,665,559 bp which is associated with APX, AsA, GSH, GS, WGS, and TKW under n-K treatment. Inside this region, our candidate gene is HORVU.MOREX.r3.2HG0111480 that annotated as NAC domain protein. Allelic variation detected that the accessions carrying C allele showed higher antioxidants (APX, AsA, and GSH) and barley yield traits (GS, WGS, and TKW) than the accessions carrying A allele, suggesting a positive selection of the accessions carrying C allele that could be used to develop barley varieties with improved salt stress resilience.

Rapid Evolution of the Fine-scale Recombination Landscape in Wild House Mouse (Mus musculus) Populations.

Meiotic recombination is an important evolutionary force and an essential meiotic process. In many species, recombination events concentrate into hotspots defined by the site-specific binding of PRMD9. Rapid evolution of Prdm9's zinc finger DNA-binding array leads to remarkably abrupt shifts in the genomic distribution of hotspots between species, but the question of how Prdm9 allelic variation shapes the landscape of recombination between populations remains less well understood. Wild house mice (Mus musculus) harbor exceptional Prdm9 diversity, with >150 alleles identified to date, and pose a particularly powerful system for addressing this open question. We employed a coalescent-based approach to construct broad- and fine-scale sex-averaged recombination maps from contemporary patterns of linkage disequilibrium in nine geographically isolated wild house mouse populations, including multiple populations from each of three subspecies. Comparing maps between wild mouse populations and subspecies reveals several themes. First, we report weak fine- and broad-scale recombination map conservation across subspecies and populations, with genetic divergence offering no clear prediction for recombination map divergence. Second, most hotspots are unique to one population, an outcome consistent with minimal sharing of Prdm9 alleles between surveyed populations. Finally, by contrasting aggregate hotspot activity on the X versus autosomes, we uncover evidence for population-specific differences in the degree and direction of sex dimorphism for recombination. Overall, our findings illuminate the variability of both the broad- and fine-scale recombination landscape in M. musculus and underscore the functional impact of Prdm9 allelic variation in wild mouse populations.

Evolutionary Insights from a Large-Scale Survey of Population-Genomic Variation.

The field of genomics has ushered in new methods for studying molecular-genetic variation in natural populations. However, most population-genomic studies still rely on small sample sizes (typically, <100 individuals) from single time points, leaving considerable uncertainties with respect to the behavior of relatively young (and rare) alleles and, owing to the large sampling variance of measures of variation, to the specific gene targets of unusually strong selection. Genomic sequences of ∼1,700 haplotypes distributed over a 10-year period from a natural population of the microcrustacean Daphnia pulex reveal evolutionary-genomic features at a refined scale, including previously hidden information on the behavior of rare alleles predicted by recent theory. Background selection, resulting from the recurrent introduction of deleterious alleles, appears to strongly influence the dynamics of neutral alleles, inducing indirect negative selection on rare variants and positive selection on common variants. Temporally fluctuating selection increases the persistence of nonsynonymous alleles with intermediate frequencies, while reducing standing levels of variation at linked silent sites. Combined with the results from an equally large metapopulation survey of the study species, classes of genes that are under strong positive selection can now be confidently identified in this key model organism. Most notable among rapidly evolving Daphnia genes are those associated with ribosomes, mitochondrial functions, sensory systems, and lifespan determination.

Haplotype analysis and linkage disequilibrium of ApoB gene polymorphisms and its relationship with hyperlipidemia in patients with acne vulgaris.

BACKGROUND: Acne vulgaris (AV) is a chronic, multifactorial inflammatory disease of the pilosebaceous unit brought on by hormonal imbalance, excessive sebum production, follicular hyperkeratinization, inflammation and Cutibacterium acne. Acne patients are characterized by alteration of the lipid profile. Apolipoprotein B gene (ApoB) plays an essential role in lipoprotein biosynthesis and multiple single-nucleotide polymorphisms (SNPs) in ApoB are associated with dyslipidemia. AIM: The aim of this study was to estimate the alteration of lipid profiles in AV, determine the genetic association with lipid profile alteration by studying the ApoB gene polymorphisms, and to identify the exact haplotypes associated with acne and lipid profile alteration. SUBJECTS AND METHODS: In a case-control study consisting of 63 non-obese acne patients and 43 healthy controls, all participants underwent biochemical, anthropological assessments, and genetic analysis for ApoB polymorphisms. RESULT: Our results indicate that serum ApoB and the lipid profile were higher in acne patients compared with healthy subject. The most common haplotypes in acne patients were rs562338 A/rs17240441 I/c.12669 A/rs1042034 G, whereas the most common haplotypes in healthy subjects were rs562338 G/rs17240441 D/c.12669 A/rs1042034 G. Patients with mild acne had higher serum ApoB levels p = 0.005. Also, the low-density lipoprotein cholesterol (LDL-C) level was higher in mild acne compared with other acne groups, with a highly significant variation of p ≤ 0.001. CONCLUSION: We found a significant variation between the acne group and healthy controls in serum ApoB, triglycerides, total cholesterol and LDL-C. The most common haplotypes in acne patients are rs562338 A/, rs17240441 I/, c.12669 A/ and rs1042034 G, and there is a linkage disequilibrium between the four selected SNPs.

Genome-wide association study identifies novel loci and candidate genes for rust resistance in wheat (Triticum aestivum L.).

BACKGROUND: Wheat rusts are important biotic stresses, development of rust resistant cultivars through molecular approaches is both economical and sustainable. Extensive phenotyping of large mapping populations under diverse production conditions and high-density genotyping would be the ideal strategy to identify major genomic regions for rust resistance in wheat. The genome-wide association study (GWAS) population of 280 genotypes was genotyped using a 35 K Axiom single nucleotide polymorphism (SNP) array and phenotyped at eight, 10, and, 10 environments, respectively for stem/black rust (SR), stripe/yellow rust (YR), and leaf/brown rust (LR). RESULTS: Forty-one Bonferroni corrected marker-trait associations (MTAs) were identified, including 17 for SR and 24 for YR. Ten stable MTAs and their best combinations were also identified. For YR, AX-94990952 on 1A + AX-95203560 on 4A + AX-94723806 on 3D + AX-95172478 on 1A showed the best combination with an average co-efficient of infection (ACI) score of 1.36. Similarly, for SR, AX-94883961 on 7B + AX-94843704 on 1B and AX-94883961 on 7B + AX-94580041 on 3D + AX-94843704 on 1B showed the best combination with an ACI score of around 9.0. The genotype PBW827 have the best MTA combinations for both YR and SR resistance. In silico study identifies key prospective candidate genes that are located within MTA regions. Further, the expression analysis revealed that 18 transcripts were upregulated to the tune of more than 1.5 folds including 19.36 folds (TraesCS3D02G519600) and 7.23 folds (TraesCS2D02G038900) under stress conditions compared to the control conditions. Furthermore, highly expressed genes in silico under stress conditions were analyzed to find out the potential links to the rust phenotype, and all four genes were found to be associated with the rust phenotype. CONCLUSION: The identified novel MTAs, particularly stable and highly expressed MTAs are valuable for further validation and subsequent application in wheat rust resistance breeding. The genotypes with favorable MTA combinations can be used as prospective donors to develop elite cultivars with YR and SR resistance.

Genome-wide association study of drought tolerance in wheat (Triticum aestivum L.) identifies SNP markers and candidate genes.

Drought stress poses a severe threat to global wheat production, necessitating an in-depth exploration of the genetic basis for drought tolerance associated traits. This study employed a 90 K SNP array to conduct a genome-wide association analysis, unravelling genetic determinants of key traits related to drought tolerance in wheat, namely plant height, root length, and root and shoot dry weight. Using the mixed linear model (MLM) method on 125 wheat accessions subjected to both well-watered and drought stress treatments, we identified 53 SNPs significantly associated with stress susceptibility (SSI) and tolerance indices (STI) for the targeted traits. Notably, chromosomes 2A and 3B stood out with ten and nine associated markers, respectively. Across 17 chromosomes, 44 unique candidate genes were pinpointed, predominantly located on the distal ends of 1A, 1B, 1D, 2A, 3A, 3B, 4A, 6A, 6B, 7A, 7B, and 7D chromosomes. These genes, implicated in diverse functions related to plant growth, development, and stress responses, offer a rich resource for future investigation. A clustering pattern emerged, notably with seven genes associated with SSI for plant height and four genes linked to both STI of plant height and shoot dry weight, converging on specific regions of chromosome arms of 2AS and 3BL. Additionally, shared genes encoding polygalacturonase, auxilin-related protein 1, peptide deformylase, and receptor-like kinase underscored the interconnectedness between plant height and shoot dry weight. In conclusion, our findings provide insights into the molecular mechanisms governing wheat drought tolerance, identifying promising genomic loci for further exploration and crop improvement strategies.

LmTag: functional-enrichment and imputation-aware tag SNP selection for population-specific genotyping arrays.

Despite the rapid development of sequencing technology, single-nucleotide polymorphism (SNP) arrays are still the most cost-effective genotyping solutions for large-scale genomic research and applications. Recent years have witnessed the rapid development of numerous genotyping platforms of different sizes and designs, but population-specific platforms are still lacking, especially for those in developing countries. SNP arrays designed for these countries should be cost-effective (small size), yet incorporate key information needed to associate genotypes with traits. A key design principle for most current platforms is to improve genome-wide imputation so that more SNPs not included in the array (imputed SNPs) can be predicted. However, current tag SNP selection methods mostly focus on imputation accuracy and coverage, but not the functional content of the array. It is those functional SNPs that are most likely associated with traits. Here, we propose LmTag, a novel method for tag SNP selection that not only improves imputation performance but also prioritizes highly functional SNP markers. We apply LmTag on a wide range of populations using both public and in-house whole-genome sequencing databases. Our results show that LmTag improved both functional marker prioritization and genome-wide imputation accuracy compared to existing methods. This novel approach could contribute to the next generation genotyping arrays that provide excellent imputation capability as well as facilitate array-based functional genetic studies. Such arrays are particularly suitable for under-represented populations in developing countries or non-model species, where little genomics data are available while investment in genome sequencing or high-density SNP arrays is limited. $\textrm{LmTag}$ is available at: https://github.com/datngu/LmTag.

Genetic insight into a polygenic trait using a novel genome-wide association approach in a wild amphibian population.

Body size variation is central in the evolution of life-history traits in amphibians, but the underlying genetic architecture of this complex trait is still largely unknown. Herein, we studied the genetic basis of body size and fecundity of the alternative morphotypes in a wild population of the Greek smooth newt (Lissotriton graecus). By combining a genome-wide association approach with linkage disequilibrium network analysis, we were able to identify clusters of highly correlated loci thus maximizing sequence data for downstream analysis. The putatively associated variants explained 12.8% to 44.5% of the total phenotypic variation in body size and were mapped to genes with functional roles in the regulation of gene expression and cell cycle processes. Our study is the first to provide insights into the genetic basis of complex traits in newts and provides a useful tool to identify loci potentially involved in fitness-related traits in small data sets from natural populations in non-model species.