FAVOR: functional annotation of variants online resource and annotator for variation across the human genome.

Large biobank-scale whole genome sequencing (WGS) studies are rapidly identifying a multitude of coding and non-coding variants. They provide an unprecedented resource for illuminating the genetic basis of human diseases. Variant functional annotations play a critical role in WGS analysis, result interpretation, and prioritization of disease- or trait-associated causal variants. Existing functional annotation databases have limited scope to perform online queries and functionally annotate the genotype data of large biobank-scale WGS studies. We develop the Functional Annotation of Variants Online Resources (FAVOR) to meet these pressing needs. FAVOR provides a comprehensive multi-faceted variant functional annotation online portal that summarizes and visualizes findings of all possible nine billion single nucleotide variants (SNVs) across the genome. It allows for rapid variant-, gene- and region-level queries of variant functional annotations. FAVOR integrates variant functional information from multiple sources to describe the functional characteristics of variants and facilitates prioritizing plausible causal variants influencing human phenotypes. Furthermore, we provide a scalable annotation tool, FAVORannotator, to functionally annotate large-scale WGS studies and efficiently store the genotype and their variant functional annotation data in a single file using the annotated Genomic Data Structure (aGDS) format, making downstream analysis more convenient. FAVOR and FAVORannotator are available at https://favor.genohub.org.

Pan-ancestry exome-wide association analyses of COVID-19 outcomes in 586,157 individuals.

Severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2) causes coronavirus disease 2019 (COVID-19), a respiratory illness that can result in hospitalization or death. We used exome sequence data to investigate associations between rare genetic variants and seven COVID-19 outcomes in 586,157 individuals, including 20,952 with COVID-19. After accounting for multiple testing, we did not identify any clear associations with rare variants either exome wide or when specifically focusing on (1) 13 interferon pathway genes in which rare deleterious variants have been reported in individuals with severe COVID-19, (2) 281 genes located in susceptibility loci identified by the COVID-19 Host Genetics Initiative, or (3) 32 additional genes of immunologic relevance and/or therapeutic potential. Our analyses indicate there are no significant associations with rare protein-coding variants with detectable effect sizes at our current sample sizes. Analyses will be updated as additional data become available, and results are publicly available through the Regeneron Genetics Center COVID-19 Results Browser.

Whole-genome sequencing association analysis of quantitative red blood cell phenotypes: The NHLBI TOPMed program.

Whole-genome sequencing (WGS), a powerful tool for detecting novel coding and non-coding disease-causing variants, has largely been applied to clinical diagnosis of inherited disorders. Here we leveraged WGS data in up to 62,653 ethnically diverse participants from the NHLBI Trans-Omics for Precision Medicine (TOPMed) program and assessed statistical association of variants with seven red blood cell (RBC) quantitative traits. We discovered 14 single variant-RBC trait associations at 12 genomic loci, which have not been reported previously. Several of the RBC trait-variant associations (RPN1, ELL2, MIDN, HBB, HBA1, PIEZO1, and G6PD) were replicated in independent GWAS datasets imputed to the TOPMed reference panel. Most of these discovered variants are rare/low frequency, and several are observed disproportionately among non-European Ancestry (African, Hispanic/Latino, or East Asian) populations. We identified a 3 bp indel p.Lys2169del (g.88717175_88717177TCT[4]) (common only in the Ashkenazi Jewish population) of PIEZO1, a gene responsible for the Mendelian red cell disorder hereditary xerocytosis (MIM: 194380), associated with higher mean corpuscular hemoglobin concentration (MCHC). In stepwise conditional analysis and in gene-based rare variant aggregated association analysis, we identified several of the variants in HBB, HBA1, TMPRSS6, and G6PD that represent the carrier state for known coding, promoter, or splice site loss-of-function variants that cause inherited RBC disorders. Finally, we applied base and nuclease editing to demonstrate that the sentinel variant rs112097551 (nearest gene RPN1) acts through a cis-regulatory element that exerts long-range control of the gene RUVBL1 which is essential for hematopoiesis. Together, these results demonstrate the utility of WGS in ethnically diverse population-based samples and gene editing for expanding knowledge of the genetic architecture of quantitative hematologic traits and suggest a continuum between complex trait and Mendelian red cell disorders.

CoFeSe2 @DMSA@FA Nanocatalyst for Amplification of Oxidative Stress to Achieve Multimodal Tumor Therapy.

Nanomedicine has significantly advanced precise tumor therapy, providing essential technical blessing for active drug accumulation, targeted consignment, and mitigation of noxious side effects. To enhance anti-tumor efficacy, the integration of multiple therapeutic modalities has garnered significant attention. Here, we designed an innovative CoFeSe2 @DMSA@FA nanocatalyst with Se vacancies (abbreviated as CFSDF), which exhibits synergistic chemodynamic therapy (CDT) and photothermal therapy (PTT), leading to amplified tumor oxidative stress and enhanced photothermal effects. The multifunctional CFSDF nanocatalyst exhibits the remarkable ability to catalyze the Fenton reaction within the acidic tumor microenvironment, efficiently converting hydrogen peroxide (H2 O2 ) into highly harmful hydroxyl radicals (⋅OH). Moreover, the nanocatalyst effectively diminishes GSH levels and ameliorates intracellular oxidative stress. The incorporation of FA modification enables CFSDF to evade immune detection and selectively target tumor tissues. Numerous in vitro and in vivo investigations have consistently demonstrated that CFSDF optimizes its individual advantages and significantly enhances therapeutic efficiency through synergistic effects of multiple therapeutic modalities, offering a valuable and effective approach to cancer treatment.

Impact of Rare and Common Genetic Variants on Diabetes Diagnosis by Hemoglobin A1c in Multi-Ancestry Cohorts: The Trans-Omics for Precision Medicine Program.

Hemoglobin A1c (HbA1c) is widely used to diagnose diabetes and assess glycemic control in individuals with diabetes. However, nonglycemic determinants, including genetic variation, may influence how accurately HbA1c reflects underlying glycemia. Analyzing the NHLBI Trans-Omics for Precision Medicine (TOPMed) sequence data in 10,338 individuals from five studies and four ancestries (6,158 Europeans, 3,123 African-Americans, 650 Hispanics, and 407 East Asians), we confirmed five regions associated with HbA1c (GCK in Europeans and African-Americans, HK1 in Europeans and Hispanics, FN3K and/or FN3KRP in Europeans, and G6PD in African-Americans and Hispanics) and we identified an African-ancestry-specific low-frequency variant (rs1039215 in HBG2 and HBE1, minor allele frequency (MAF) = 0.03). The most associated G6PD variant (rs1050828-T, p.Val98Met, MAF = 12% in African-Americans, MAF = 2% in Hispanics) lowered HbA1c (-0.88% in hemizygous males, -0.34% in heterozygous females) and explained 23% of HbA1c variance in African-Americans and 4% in Hispanics. Additionally, we identified a rare distinct G6PD coding variant (rs76723693, p.Leu353Pro, MAF = 0.5%; -0.98% in hemizygous males, -0.46% in heterozygous females) and detected significant association with HbA1c when aggregating rare missense variants in G6PD. We observed similar magnitude and direction of effects for rs1039215 (HBG2) and rs76723693 (G6PD) in the two largest TOPMed African American cohorts, and we replicated the rs76723693 association in the UK Biobank African-ancestry participants. These variants in G6PD and HBG2 were monomorphic in the European and Asian samples. African or Hispanic ancestry individuals carrying G6PD variants may be underdiagnosed for diabetes when screened with HbA1c. Thus, assessment of these variants should be considered for incorporation into precision medicine approaches for diabetes diagnosis.

SAIGEgds-an efficient statistical tool for large-scale PheWAS with mixed models.

SUMMARY: Phenome-wide association studies (PheWASs) are known to be a powerful tool in discovery and replication of genetic association studies. To reduce the computational burden of PheWAS in the large cohorts, such as the UK Biobank, the SAIGE method has been proposed to control for case-control imbalance and sample relatedness in a tractable manner. However, SAIGE is still computationally intensive when deployed in analyzing the associations of thousands of ICD10-coded phenotypes with whole-genome imputed genotype data. Here, we present a new high-performance statistical R package (SAIGEgds) for large-scale PheWAS using generalized linear mixed models. The package implements the SAIGE method in optimized C++ codes, taking advantage of sparse genotype dosages and integrating the efficient genomic data structure file format. Benchmarks using the UK Biobank White British genotype data (N ≈ 430 K) with coronary heart disease and simulated cases show that the implementation in SAIGEgds is 5-6 times faster than the SAIGE R package. When used in conjunction with high-performance computing clusters, SAIGEgds provides an efficient analysis pipeline for biobank-scale PheWAS. AVAILABILITY AND IMPLEMENTATION: https://bioconductor.org/packages/SAIGEgds; vignettes included. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

A fully adjusted two-stage procedure for rank-normalization in genetic association studies.

When testing genotype-phenotype associations using linear regression, departure of the trait distribution from normality can impact both Type I error rate control and statistical power, with worse consequences for rarer variants. Because genotypes are expected to have small effects (if any) investigators now routinely use a two-stage method, in which they first regress the trait on covariates, obtain residuals, rank-normalize them, and then use the rank-normalized residuals in association analysis with the genotypes. Potential confounding signals are assumed to be removed at the first stage, so in practice, no further adjustment is done in the second stage. Here, we show that this widely used approach can lead to tests with undesirable statistical properties, due to both combination of a mis-specified mean-variance relationship and remaining covariate associations between the rank-normalized residuals and genotypes. We demonstrate these properties theoretically, and also in applications to genome-wide and whole-genome sequencing association studies. We further propose and evaluate an alternative fully adjusted two-stage approach that adjusts for covariates both when residuals are obtained and in the subsequent association test. This method can reduce excess Type I errors and improve statistical power.

A high-resolution HLA imputation system for the Taiwanese population: a study of the Taiwan Biobank.

An imputation algorithm for human leukocyte antigen (HLA) is helpful for exploring novel disease associations. However, population-specific HLA imputation references are essential for achieving high imputation accuracy. In this study, a subset of 1012 individuals from the Taiwan Biobank (TWB) who underwent both whole-genome SNP array and NGS-based HLA typing were used to establish Taiwanese HLA imputation references. The HIBAG package was used to generate the imputation references for eight HLA loci at a two- and three-field resolution. Internal validation was carried out to evaluate the call threshold and accuracy for each HLA gene. HLA class II genes found to be associated with rheumatoid arthritis (RA) were validated in this study by the imputed HLA alleles. Our Taiwanese population-specific references achieved average HLA imputation accuracies of 98.11% for two-field and 98.08% for three-field resolution. The frequency distribution of imputed HLA alleles among 23,972 TWB subjects were comparable with PCR-based HLA alleles in general Taiwanese reported in the allele frequency net database. We replicated four common HLA alleles (HLA-DRB1*03:01, DRB1*04:05, DQA1*03:03, and DQB1*04:01) significantly associated with RA. The population-specific references provide an informative tool to investigate the associations of HLA variants and human diseases in large-scale population-based studies.

Genome-wide minor histocompatibility matching as related to the risk of graft-versus-host disease.

The risk of acute graft-versus-host disease (GVHD) is higher after allogeneic hematopoietic cell transplantation (HCT) from unrelated donors as compared with related donors. This difference has been explained by increased recipient mismatching for major histocompatibility antigens or minor histocompatibility antigens. In the current study, we used genome-wide arrays to enumerate single nucleotide polymorphisms (SNPs) that produce graft-versus-host (GVH) amino acid coding differences between recipients and donors. We then tested the hypothesis that higher degrees of genome-wide recipient GVH mismatching correlate with higher risks of GVHD after allogeneic HCT. In HLA-genotypically matched sibling recipients, the average recipient mismatching of coding SNPs was 9.35%. Each 1% increase in genome-wide recipient mismatching was associated with an estimated 20% increase in the hazard of grades III-IV GVHD (hazard ratio [HR], 1.20; 95% confidence interval [CI], 1.05-1.37; P = .007) and an estimated 22% increase in the hazard of stage 2-4 acute gut GVHD (HR, 1.22; 95% CI, 1.02-1.45; P = .03). In HLA-A, B, C, DRB1, DQA1, DQB1, DPA1, DPB1-phenotypically matched unrelated recipients, the average recipient mismatching of coding SNPs was 17.3%. The estimated risks of GVHD-related outcomes in HLA-phenotypically matched unrelated recipients were low, relative to the large difference in genome-wide mismatching between the 2 groups. In contrast, the risks of GVHD-related outcomes were higher in HLA-DP GVH-mismatched unrelated recipients than in HLA-matched sibling recipients. Taken together, these results suggest that the increased GVHD risk after unrelated HCT is predominantly an effect of HLA-mismatching.

A functional FePt@MOFs (MIL-101(Fe)) nano-platform for high efficient colorimetric determination of H2O2.

The rapid and sensitive detection of H2O2 is of great importance for a series of industries. In this work, surface-modified FePt nanoparticles were loaded into metal-organic frameworks (MOFs, MIL-101(Fe)) in a facile in situ way and the resulting material (FePt@MOFs NCs) acted as a high efficient colorimetric detection agent. The as-prepared FePt@MOFs NCs possess a well-defined octahedron shape and high stability in water. Via the oxidation reaction of the chromogenic substrate 4,4'-bi-2,6-xylidine, the obtained FePt@MOFs NCs were proved to possess prominent peroxidase-mimic activity. Under the most favorable conditions, the linear range for H2O2 detection was 40 µM to 800 µM and the detection limit was 18.9 µM. Furthermore, FePt@MOFs NCs were successfully applied to the detection of H2O2 in real samples such as fruit juice. The developed colorimetric sensing platform is expected to have potential for the exact detection of H2O2 for the food, mining and bio-pharmaceutical industries, due to its high stability, sensitivity, rapidness and easy preparation.

Capture and separation of circulating tumor cells using functionalized magnetic nanocomposites with simultaneous in situ chemotherapy.

Circulating tumor cells (CTCs) are a type of rare cell that are firstly shed from solid tumors and then exist in the bloodstream. The effective capture and separation of CTCs has significant meaning in cancer diagnosis and prognosis. In this study, novel Fe3O4-FePt magnetic nanocomposites (Fe3O4-FePt MNCs) were constructed by integrating face centered cubic (fcc) FePt nanoparticles (NPs) onto the surface of the Fe3O4@SiO2 core. After further modification with NH2-PEG-COOH and the tumor-targeting molecule tLyP-1, the acquired Fe3O4-FePt MNCs possesses excellent biocompatibility and stability and could efficiently target and capture tLyP-1 receptor-positive CTCs. Based on the acidic microenvironment within cancer cells, the FePt layer could rapidly release active Fe2+ ions, which could catalyze H2O2 into reactive oxygen species (ROS) and further induce in situ apoptosis in cancer cells while having no distinct cytotoxicity to normal cells. Moreover, the Fe3O4@SiO2 core with its intrinsic magnetism has huge potential for the bioseparation of CTCs. The in vitro ROS fluorescence imaging experiments and cell capture and separation experiments indicated that the Fe3O4-FePt MNCs could specifically capture and separate cancer cells in the CTCs model and further induce in situ apoptosis. Therefore, the Fe3O4-FePt MNCs could serve as a promising multifunctional nanoseparator for efficiently capturing CTCs and simultaneously inducing in situ chemotherapy.

Characterizing the noncovalent binding behavior of tartrazine to lysozyme: A combined spectroscopic and computational analysis.

Tartrazine is a stable water-soluble azo dye widely used as a food additive, which could pose potential threats to humans and the environment. In this paper, we evaluated the response mechanism between tartrazine and lysozyme under simulated conditions by means of biophysical methods, including multiple spectroscopic techniques, isothermal titration calorimetry (ITC), and molecular docking studies. From the multispectroscopic analysis, we found that tartrazine could effectively quench the intrinsic fluorescence of lysozyme to form a complex and lead to the conformational and microenvironmental changes of the enzyme. The ITC measurements suggested that the electrostatic forces played a major role in the binding of tartrazine to lysozyme with two binding sites. Finally, the molecular docking indicated that tartrazine had specific interactions with the residues of Trp108. The study provides an important insight within the binding mechanism of tartrazine to lysozyme in vitro.

A facile preparation of FePt-loaded few-layer MoS2 nanosheets nanocomposites (F-MoS2-FePt NCs) and their application for colorimetric detection of H2O2 in living cells.

BACKGROUND: Rapid and sensitive detection of H2O2 especially endogenous H2O2 is of great importance for series of industries including disease diagnosis and therapy. In this work, uniform FePt nanoparticles are successfully anchored onto Few-layer molybdenum disulfide nanosheets (F-MoS2 NSs). The powder X-ray diffraction, transmission electron microscopy, UV-Vis spectra and atomic force microscopy were employed to confirm the structure of the obtained nanocomposites (F-MoS2-FePt NCs). The prepared nanocomposites show efficient peroxidase-like catalytic activities verified by catalyzing the peroxidation substrate 4,4'-diamino-3,3',5,5'-tetramethylbiphenyl (TMB) with the existence of H2O2. RESULTS: The optimal conditions of the constructed colorimetric sensing platform is proved as 35 °C and pH 4.2. Under optimal catalytic conditions, the detection limit for H2O2 detection reaches 2.24 µM and the linear ranger is 8 µM to 300 µM. Furthermore, the proposed colorimetric sensing platform was successfully utilized to detect the intracellular H2O2 of cancer cells (MCF-7). CONCLUSIONS: These findings indicated that the F-MoS2-FePt-TMB-H2O2 system provides a potential sensing platform for hydrogen peroxide monitoring in living cells.

HLA-DPB1 and HLA class I confer risk of and protection from narcolepsy.

Type 1 narcolepsy, a disorder caused by a lack of hypocretin (orexin), is so strongly associated with human leukocyte antigen (HLA) class II HLA-DQA1(∗)01:02-DQB1(∗)06:02 (DQ0602) that very few non-DQ0602 cases have been reported. A known triggering factor for narcolepsy is pandemic 2009 influenza H1N1, suggesting autoimmunity triggered by upper-airway infections. Additional effects of other HLA-DQ alleles have been reported consistently across multiple ethnic groups. Using over 3,000 case and 10,000 control individuals of European and Chinese background, we examined the effects of other HLA loci. After careful matching of HLA-DR and HLA-DQ in case and control individuals, we found strong protective effects of HLA-DPA1(∗)01:03-DPB1(∗)04:02 (DP0402; odds ratio [OR] = 0.51 [0.38-0.67], p = 1.01 × 10(-6)) and HLA-DPA1(∗)01:03-DPB1(∗)04:01 (DP0401; OR = 0.61 [0.47-0.80], p = 2.07 × 10(-4)) and predisposing effects of HLA-DPB1(∗)05:01 in Asians (OR = 1.76 [1.34-2.31], p = 4.71 × 10(-05)). Similar effects were found by conditional analysis controlling for HLA-DR and HLA-DQ with DP0402 (OR = 0.45 [0.38-0.55] p = 8.99 × 10(-17)) and DP0501 (OR = 1.38 [1.18-1.61], p = 7.11 × 10(-5)). HLA-class-II-independent associations with HLA-A(∗)11:01 (OR = 1.32 [1.13-1.54], p = 4.92 × 10(-4)), HLA-B(∗)35:03 (OR = 1.96 [1.41-2.70], p = 5.14 × 10(-5)), and HLA-B(∗)51:01 (OR = 1.49 [1.25-1.78], p = 1.09 × 10(-5)) were also seen across ethnic groups in the HLA class I region. These effects might reflect modulation of autoimmunity or indirect effects of HLA class I and HLA-DP alleles on response to viral infections such as that of influenza.

SeqArray-a storage-efficient high-performance data format for WGS variant calls.

MOTIVATION: Whole-genome sequencing (WGS) data are being generated at an unprecedented rate. Analysis of WGS data requires a flexible data format to store the different types of DNA variation. Variant call format (VCF) is a general text-based format developed to store variant genotypes and their annotations. However, VCF files are large and data retrieval is relatively slow. Here we introduce a new WGS variant data format implemented in the R/Bioconductor package 'SeqArray' for storing variant calls in an array-oriented manner which provides the same capabilities as VCF, but with multiple high compression options and data access using high-performance parallel computing. RESULTS: Benchmarks using 1000 Genomes Phase 3 data show file sizes are 14.0 Gb (VCF), 12.3 Gb (BCF, binary VCF), 3.5 Gb (BGT) and 2.6 Gb (SeqArray) respectively. Reading genotypes in the SeqArray package are two to three times faster compared with the htslib C library using BCF files. For the allele frequency calculation, the implementation in the SeqArray package is over 5 times faster than PLINK v1.9 with VCF and BCF files, and over 16 times faster than vcftools. When used in conjunction with R/Bioconductor packages, the SeqArray package provides users a flexible, feature-rich, high-performance programming environment for analysis of WGS variant data. AVAILABILITY AND IMPLEMENTATION: http://www.bioconductor.org/packages/SeqArray. CONTACT: zhengx@u.washington.edu. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

pH-Responsive, Self-Sacrificial Nanotheranostic Agent for Potential In Vivo and In Vitro Dual Modal MRI/CT Imaging, Real-Time, and In Situ Monitoring of Cancer Therapy.

Multifunctional nanotheranostic agents have been highly commended due to the application to image-guided cancer therapy. Herein, based on the chemically disordered face centered cubic (fcc) FePt nanoparticles (NPs) and graphene oxide (GO), we develop a pH-responsive FePt-based multifunctional theranostic agent for potential in vivo and in vitro dual modal MRI/CT imaging and in situ cancer inhibition. The fcc-FePt will release highly active Fe ions due to the low pH in tumor cells, which would catalyze H2O2 decomposition into reactive oxygen species (ROS) within the cells and further induce cancer cell apoptosis. Conjugated with folic acid (FA), the iron platinum-dimercaptosuccinnic acid/PEGylated graphene oxide-folic acid (FePt-DMSA/GO-PEG-FA) composite nanoassemblies (FePt/GO CNs) could effectively target and show significant toxicity to FA receptor-positive tumor cells, but no obvious toxicity to FA receptor-negative normal cells, which was evaluated by WST-1 assay. The FePt-based multifunctional nanoparticles allow real-time monitoring of Fe release by T2-weighted MRI, and the selective contrast enhancement in CT could be estimated in vivo after injection. The results showed that FePt-based NPs displayed excellent biocompatibility and favorable MRI/CT imaging ability in vivo and in vitro. Meanwhile, the decomposition of FePt will dramatically decrease the T2-weighted MRI signal and increase the ROS signal, which enables real-time and in situ visualized monitoring of Fe release in tumor cells. In addition, the self-sacrificial decomposition of fcc-FePt will be propitious to the self-clearance of the as-prepared FePt-based nanocomposite in vivo. Therefore, the FePt/GO CNs could serve as a potential multifunctional theranostic nanoplatform of MRI/CT imaging guided cancer diagnosis and therapy in the clinic.

Phenome-wide association study maps new diseases to the human major histocompatibility complex region.

BACKGROUND: Over 160 disease phenotypes have been mapped to the major histocompatibility complex (MHC) region on chromosome 6 by genome-wide association study (GWAS), suggesting that the MHC region as a whole may be involved in the aetiology of many phenotypes, including unstudied diseases. The phenome-wide association study (PheWAS), a powerful and complementary approach to GWAS, has demonstrated its ability to discover and rediscover genetic associations. The objective of this study is to comprehensively investigate the MHC region by PheWAS to identify new phenotypes mapped to this genetically important region. METHODS: In the current study, we systematically explored the MHC region using PheWAS to associate 2692 MHC-linked variants (minor allele frequency ≥0.01) with 6221 phenotypes in a cohort of 7481 subjects from the Marshfield Clinic Personalized Medicine Research Project. RESULTS: Findings showed that expected associations previously identified by GWAS could be identified by PheWAS (eg, psoriasis, ankylosing spondylitis, type I diabetes and coeliac disease) with some having strong cross-phenotype associations potentially driven by pleiotropic effects. Importantly, novel associations with eight diseases not previously assessed by GWAS (eg, lichen planus) were also identified and replicated in an independent population. Many of these associated diseases appear to be immune-related disorders. Further assessment of these diseases in 16 484 Marshfield Clinic twins suggests that some of these diseases, including lichen planus, may have genetic aetiologies. CONCLUSIONS: These results demonstrate that the PheWAS approach is a powerful and novel method to discover SNP-disease associations, and is ideal when characterising cross-phenotype associations, and further emphasise the importance of the MHC region in human health and disease.

Eigenanalysis of SNP data with an identity by descent interpretation.

Principal component analysis (PCA) is widely used in genome-wide association studies (GWAS), and the principal component axes often represent perpendicular gradients in geographic space. The explanation of PCA results is of major interest for geneticists to understand fundamental demographic parameters. Here, we provide an interpretation of PCA based on relatedness measures, which are described by the probability that sets of genes are identical-by-descent (IBD). An approximately linear transformation between ancestral proportions (AP) of individuals with multiple ancestries and their projections onto the principal components is found. In addition, a new method of eigenanalysis "EIGMIX" is proposed to estimate individual ancestries. EIGMIX is a method of moments with computational efficiency suitable for millions of SNP data, and it is not subject to the assumption of linkage equilibrium. With the assumptions of multiple ancestries and their surrogate ancestral samples, EIGMIX is able to infer ancestral proportions (APs) of individuals. The methods were applied to the SNP data from the HapMap Phase 3 project and the Human Genome Diversity Panel. The APs of individuals inferred by EIGMIX are consistent with the findings of the program ADMIXTURE. In conclusion, EIGMIX can be used to detect population structure and estimate genome-wide ancestral proportions with a relatively high accuracy.

Genetic variation associated with circulating monocyte count in the eMERGE Network.

With white blood cell count emerging as an important risk factor for chronic inflammatory diseases, genetic associations of differential leukocyte types, specifically monocyte count, are providing novel candidate genes and pathways to further investigate. Circulating monocytes play a critical role in vascular diseases such as in the formation of atherosclerotic plaque. We performed a joint and ancestry-stratified genome-wide association analyses to identify variants specifically associated with monocyte count in 11 014 subjects in the electronic Medical Records and Genomics Network. In the joint and European ancestry samples, we identified novel associations in the chromosome 16 interferon regulatory factor 8 (IRF8) gene (P-value = 2.78×10(-16), ß = -0.22). Other monocyte associations include novel missense variants in the chemokine-binding protein 2 (CCBP2) gene (P-value = 1.88×10(-7), ß = 0.30) and a region of replication found in ribophorin I (RPN1) (P-value = 2.63×10(-16), ß = -0.23) on chromosome 3. The CCBP2 and RPN1 region is located near GATA binding protein2 gene that has been previously shown to be associated with coronary heart disease. On chromosome 9, we found a novel association in the prostaglandin reductase 1 gene (P-value = 2.29×10(-7), ß = 0.16), which is downstream from lysophosphatidic acid receptor 1. This region has previously been shown to be associated with monocyte count. We also replicated monocyte associations of genome-wide significance (P-value = 5.68×10(-17), ß = -0.23) at the integrin, alpha 4 gene on chromosome 2. The novel IRF8 results and further replications provide supporting evidence of genetic regions associated with monocyte count.