Detecting fitness epistasis in recently admixed populations with genome-wide data.

BACKGROUND: Fitness epistasis, the interaction effect of genes at different loci on fitness, makes an important contribution to adaptive evolution. Although fitness interaction evidence has been observed in model organisms, it is more difficult to detect and remains poorly understood in human populations as a result of limited statistical power and experimental constraints. Fitness epistasis is inferred from non-independence between unlinked loci. We previously observed ancestral block correlation between chromosomes 4 and 6 in African Americans. The same approach fails when examining ancestral blocks on the same chromosome due to the strong confounding effect observed in a recently admixed population. RESULTS: We developed a novel approach to eliminate the bias caused by admixture linkage disequilibrium when searching for fitness epistasis on the same chromosome. We applied this approach in 16,252 unrelated African Americans and identified significant ancestral correlations in two pairs of genomic regions (P-value< 8.11 × 10- 7) on chromosomes 1 and 10. The ancestral correlations were not explained by population admixture. Historical African-European crossover events are reduced between pairs of epistatic regions. We observed multiple pairs of co-expressed genes shared by the two regions on each chromosome, including ADAR being co-expressed with IFI44 in almost all tissues and DARC being co-expressed with VCAM1, S1PR1 and ELTD1 in multiple tissues in the Genotype-Tissue Expression (GTEx) data. Moreover, the co-expressed gene pairs are associated with the same diseases/traits in the GWAS Catalog, such as white blood cell count, blood pressure, lung function, inflammatory bowel disease and educational attainment. CONCLUSIONS: Our analyses revealed two instances of fitness epistasis on chromosomes 1 and 10, and the findings suggest a potential approach to improving our understanding of adaptive evolution.

Genome-wide survey in African Americans demonstrates potential epistasis of fitness in the human genome.

The role played by epistasis between alleles at unlinked loci in shaping population fitness has been debated for many years and the existing evidence has been mainly accumulated from model organisms. In model organisms, fitness epistasis can be systematically inferred by detecting nonindependence of genotypic values between loci in a population and confirmed through examining the number of offspring produced in two-locus genotype groups. No systematic study has been conducted to detect epistasis of fitness in humans owing to experimental constraints. In this study, we developed a novel method to detect fitness epistasis by testing the correlation between local ancestries on different chromosomes in an admixed population. We inferred local ancestry across the genome in 16,252 unrelated African Americans and systematically examined the pairwise correlations between the genomic regions on different chromosomes. Our analysis revealed a pair of genomic regions on chromosomes 4 and 6 that show significant local ancestry correlation (P-value = 4.01 × 10-8 ) that can be potentially attributed to fitness epistasis. However, we also observed substantial local ancestry correlation that cannot be explained by systemic ancestry inference bias. To our knowledge, this study is the first to systematically examine evidence of fitness epistasis across the human genome.

Whole exome analyses to examine the impact of rare variants on left ventricular traits in African American participants from the HyperGEN and GENOA studies.

Left ventricular (LV) hypertrophy, highest in prevalence among African Americans, is an established risk factor heart failure. Several genome wide association studies have identified common variants associated with LV-related quantitative-traits in African Americans. To date, however, the effect of rare variants on these traits has not been extensively studied, especially in minority groups. We therefore investigated the association between rare variants and LV traits among 1,934 African Americans using exome chip data from the Hypertension Genetic Epidemiology Network (HyperGEN) study, with replication in 1,090 African American from the Genetic Epidemiology Network of Arteriopathy (GENOA) study. We used single-variant analyses and gene-based tests to investigate the association between 86,927 variants and six structural and functional LV traits including LV mass, LV internal dimension-diastole, relative wall thickness, left atrial dimension (LAD), fractional shortening (FS), and the ratio of LV early-to-late transmitral velocity (E/A ratio). Only rare variants (MAF <1% and <5%) were considered in gene-based analyses. In gene-based analyses, we found a statistically significant association between potassium voltage-gated channel subfamily H member 4 (KCNH4) and E/A ratio (P=8.7*10-8 using a burden test). Endonuclease G (ENDOG) was associated with LAD using the Madsen Browning weighted burden (MB) test (P=1.4*10-7). Neither gene result was replicated in GENOA, but the direction of effect of single variants in common was comparable. G protein-coupled receptor 55 (GPR55) was marginally associated with LAD in HyperGEN (P=3.2*10-5 using the MB test) and E/A ratio in GENOA, but with opposing directions of association for variants in common (P=0.03 for the MB test). No single variant was statistically significantly associated with any trait after correcting for multiple testing. The findings in this study highlight the potential cumulative contributions of rare variants to LV traits which, if validated, could improve our understanding of heart failure in African Americans.

Genome-wide detection of allele specific copy number variation associated with insulin resistance in African Americans from the HyperGEN study.

African Americans have been understudied in genome wide association studies of diabetes and related traits. In the current study, we examined the joint association of single nucleotide polymorphisms (SNPs) and copy number variants (CNVs) with fasting insulin and an index of insulin resistance (HOMA-IR) in the HyperGEN study, a family based study with proband ascertainment for hypertension. This analysis is restricted to 1,040 African Americans without diabetes. We generated allele specific CNV genotypes at 872,243 autosomal loci using Birdsuite, a freely available multi-stage program. Joint tests of association for SNPs and CNVs were performed using linear mixed models adjusting for covariates and familial relationships. Our results highlight SNPs associated with fasting insulin and HOMA-IR (rs6576507 and rs8026527, 3.7*10(-7)≤P≤1.1*10(-5)) near ATPase, class V, type 10A (ATP10A), and the L Type voltage dependent calcium channel (CACNA1D, rs1401492, P≤5.2*10(-6)). ATP10A belongs to a family of aminophospholipid-transporting ATPases and has been associated with type 2 diabetes in mice. CACNA1D has been linked to pancreatic beta cell generation in mice. The two most significant copy variable markers (rs10277702 and rs361367; P<2.0*10(-4)) were in the beta variable region of the T-cell receptor gene (TCRVB). Human and mouse TCR has been shown to mimic insulin and its receptor and could contribute to insulin resistance. Our findings differ from genome wide association studies of fasting insulin and other diabetes related traits in European populations, highlighting the continued need to investigate unique genetic influences for understudied populations such as African Americans.

Bivariate genetic association of KIAA1797 with heart rate in American Indians: the Strong Heart Family Study.

Heart rate (HR) has been identified as a risk factor for cardiovascular disease (CVD), yet little is known regarding genetic factors influencing this phenotype. Previous research in American Indians (AIs) from the Strong Heart Family Study (SHFS) identified a significant quantitative trait locus (QTL) for HR on chromosome 9p21. Genetic association on HR was conducted in the SHFS. HR was measured from electrocardiogram (ECG) and echocardiograph (Echo) Doppler recordings. We examined 2248 single-nucleotide polymorphisms (SNPs) on chromosome 9p21 for association using a gene-centric statistical test. We replicated the aforementioned QTL [logarithm of odds (LOD) = 4.83; genome-wide P= 0.0003] on chromosome 9p21 in one SHFS population using joint linkage of ECG and Echo HR. After correcting for effective number of SNPs using a gene-centric test, six SNPs (rs7875153, rs7848524, rs4446809, rs10964759, rs1125488 and rs7853123) remained significant. We applied a novel bivariate association method, which was a joint test of association of a single locus to two traits using a standard additive genetic model. The SNP, rs7875153, provided the strongest evidence for association (P = 7.14 x 10(-6)). This SNP (rs7875153) is rare (minor allele frequency = 0.02) in AIs and is located within intron 9 of the gene KIAA1797. To support this association, we applied lymphocyte RNA expression data from the San Antonio Family Heart Study, a longitudinal study of CVD in Mexican Americans. Expression levels of KIAA1797 were significantly associated (P = 0.012) with HR. These findings in independent populations support that KIAA1797 genetic variation may be associated with HR but elucidation of a functional relationship requires additional study.