Twenty five years after the discovery of HIV as the cause of AIDS, there is still no effective vaccine and no cure for this disease. HIV susceptibility shows a substantial degree of individual heterogeneity, much of which can be conferred by host genetic variation. In an effort to discover host factors required for HIV replication, identify critical pathogenic pathways, and reveal the full armament of host defenses, there has been a shift from candidate gene studies to unbiased genome wide genetic and functional studies. However, the number of securely identified host factors involved in HIV disease remains small, explaining only ~15–20% of the observed heterogeneity – most of which is attributable to HLA. Multidisciplinary approaches integrating genetic epidemiology to systems biology will be required to fully understand viral-host interactions to effectively combat HIV/AIDS.
Nucleoside analogue reverse transcriptase inhibitors are an integral component of combination antiretroviral treatment regimens. However, their ability to inhibit polymerase-γ has been associated with several mitochondrial toxicities, including potentially life-threatening lactic acidosis. A total of 650 antiretroviral-naive adults (69% female) initiated combination antiretroviral therapy (cART) and were intensively screened for toxicities including lactic acidosis as part of a 3-year clinical trial in Botswana. Patients were categorized as no lactic acidosis symptoms, minor symptoms but lactate <4.4 mmol/liter, and symptoms with lactate ≥4.4 mmol/liter [moderate to severe symptomatic hyperlactatemia (SH) or lactic acidosis (LA)]. Of 650 participants 111 (17.1%) developed symptoms and/or laboratory results suggestive of lactic acidosis and had a serum lactate drawn; 97 (87.4%) of these were female. There were 20 events, 13 having SH and 7 with LA; all 20 (100%) were female (p<0.001). Cox proportional hazard analysis limited to the 451 females revealed that having a higher baseline BMI was predictive for the development of SH/LA [aHR=1.17 per one-unit increase (1.08–1.25), p<0.0001]. Ordered logistic regression performed among all 650 patients revealed that having a lower baseline hemoglobin [aOR=1.28 per one-unit decrease (1.1–1.49), p=0.002] and being randomized to d4T/3TC-based cART [aOR=1.76 relative to ZDV/3TC (1.03–3.01), p=0.04] were predictive of the symptoms and/or the development of SH/LA. cART-treated women in sub-Saharan Africa, especially those having higher body mass indices, should receive additional monitoring for SH/LA. Women presently receiving d4T/3TC-based cART in such settings also warrant more intensive monitoring.
Multiple genome-wide association studies (GWAS) have been performed in HIV-1 infected individuals, identifying common genetic influences on viral control and disease course. Similarly, common genetic correlates of acquisition of HIV-1 after exposure have been interrogated using GWAS, although in generally small samples. Under the auspices of the International Collaboration for the Genomics of HIV, we have combined the genome-wide single nucleotide polymorphism (SNP) data collected by 25 cohorts, studies, or institutions on HIV-1 infected individuals and compared them to carefully matched population-level data sets (a list of all collaborators appears in Note S1 in Text S1). After imputation using the 1,000 Genomes Project reference panel, we tested approximately 8 million common DNA variants (SNPs and indels) for association with HIV-1 acquisition in 6,334 infected patients and 7,247 population samples of European ancestry. Initial association testing identified the SNP rs4418214, the C allele of which is known to tag the HLA-B*57:01 and B*27:05 alleles, as genome-wide significant (p = 3.6×10−11). However, restricting analysis to individuals with a known date of seroconversion suggested that this association was due to the frailty bias in studies of lethal diseases. Further analyses including testing recessive genetic models, testing for bulk effects of non-genome-wide significant variants, stratifying by sexual or parenteral transmission risk and testing previously reported associations showed no evidence for genetic influence on HIV-1 acquisition (with the exception of CCR5Δ32 homozygosity). Thus, these data suggest that genetic influences on HIV acquisition are either rare or have smaller effects than can be detected by this sample size.
Comparing the frequency differences between common DNA variants in disease-affected cases and in unaffected controls has been successful in uncovering the genetic component of multiple diseases. This approach is most effective when large samples of cases and controls are available. Here we combine information from multiple studies of HIV infected patients, including more than 6,300 HIV+ individuals, with data from 7,200 general population samples of European ancestry to test nearly 8 million common DNA variants for an impact on HIV acquisition. With this large sample we did not observe any single common genetic variant that significantly associated with HIV acquisition. We further tested 22 variants previously identified by smaller studies as influencing HIV acquisition. With the exception of a deletion polymorphism in the CCR5 gene (CCR5Δ32) we found no convincing evidence to support these previous associations. Taken together these data suggest that genetic influences on HIV acquisition are either rare or have smaller effects than can be detected by this sample size.
Despite intensive anti-hypertensive therapy there was a high incidence of renal end-points in participants of the African American Study of Kidney Disease and Hypertension (AASK) cohort. To better understand this, coding variants in the apolipoprotein L1 (APOL1) and the non-muscle myosin heavy chain 9 (MYH9) genes were evaluated for an association with hypertension-attributed nephropathy and clinical outcomes in a case-control study. Clinical data and DNA were available for 675 AASK participant cases and 618 African American non-nephropathy control individuals. APOL1 G1 and G2, and MYH9 E1 variants along with 44 ancestry informative markers were genotyped with allele frequency differences between cases and controls analyzed by logistic regression multivariable models adjusting for ancestry, age, and gender. In recessive models, APOL1 risk variants were significantly associated with kidney disease in all cases compared to controls with an odds ratio of 2.57. In AASK cases with more advanced disease, such as a baseline urine protein to creatinine ratio over 0.6 g/g or a serum creatinine over 3 mg/dL during follow-up, the association was strengthened with odds ratios of 6.29 and 4.61, respectively. APOL1 risk variants were consistently associated with renal disease progression across medication classes and blood pressure targets. Thus, kidney disease in AASK participants was strongly associated with APOL1 renal risk variants.
The human APOBEC3 family of cytidine deaminases provides intrinsic immunity to retroviral infection. A naturally occurring 29.5-kb deletion removes the entire APOBEC3B gene. We examined the impact of the APOBEC3B gene deletion in >4000 individuals from five HIV-1 natural history cohorts. The hemizygous genotype had no effect on either infection or progression. However, the homozygous deletion was significantly associated with unfavorable outcomes for HIV-1 acquisition (OR=7.37, P=0.024), progression to AIDS (RH = 4.01, P=0.03), and viral set-point (P=0.04). These findings suggest that the loss of APOBEC3B may increase host susceptibility to HIV-1/AIDS and warrant further study.
Background. To date, only mutations in CCR5 have been shown to confer resistance to human immunodeficiency virus type 1 (HIV-1) infection, and these explain only a small fraction of the observed variability in HIV susceptibility.
Methods. We performed a meta-analysis between 2 independent European genomewide association studies, each comparing HIV-1 seropositive cases with normal population controls known to be HIV uninfected, to identify single-nucleotide polymorphisms (SNPs) associated with the HIV-1 acquisition phenotype. SNPs exhibiting P < 10−5 in this first stage underwent second-stage analysis in 2 independent US cohorts of European descent.
Results. After the first stage, a single highly significant association was revealed for the chromosome 8 rs6996198 with HIV-1 acquisition and was replicated in both second-stage cohorts. Across the 4 groups, the rs6996198-T allele was consistently associated with a significant reduced risk of HIV-1 infection, and the global meta-analysis reached genomewide significance: Pcombined = 7.76 × 10−8.
Conclusions. We provide strong evidence of association for a common variant with HIV-1 acquisition in populations of European ancestry. This protective signal against HIV-1 infection is the first identified outside the CCR5 nexus. First clues point to a potential functional role for a nearby candidate gene, CYP7B1, but this locus warrants further investigation.
Recently, an association was found between non-diabetic kidney disease in African Americans and two independent sequence variants in the APOL1 gene, encoding apolipoprotein L1. In this study we determined the frequency of APOL1 risk variants in patients with biopsy-proven HIV-associated nephropathy (HIVAN) and distinctive pathological characteristics potentially driven by those risk variants. Among 76 patients with HIVAN, 60 were successfully genotyped for APOL1 G1 and G2 polymorphisms. In this cohort, 37 had two risk alleles, 18 were heterozygous and 5 had neither risk variant. There were no differences in the pathological findings of HIVAN and the number of APOL1 risk alleles. Further, the progression to end stage kidney disease or death did not differ by the number of risk alleles. Median renal survival was 9.3 months in patients with none or one risk allele compared to 11.7 months in patients with two APOL1 risk alleles. Thus, our study suggests that although the majority of African American patients with HIVAN have two APOL1 risk alleles, other as yet unknown factors in the host including genetic risk variants and environmental or viral factors may influence the development of this disorder in those with none or one APOL1 risk allele.
HIV-1 derives from multiple independent transfers of simian immunodeficiency virus (SIV) strains from chimpanzees to human populations. We hypothesized that human populations in west central Africa may have been exposed to SIV prior to the pandemic, and that previous outbreaks may have selected for genetic resistance to immunodeficiency viruses. To test this hypothesis, we examined the genomes of Biaka Western Pygmies, who historically resided in communities within the geographic range of the central African chimpanzee subspecies (Pan troglodytes troglodytes) that carries strains of SIV ancestral to HIV-1.
SNP genotypes of the Biaka were compared to those of African human populations who historically resided outside the range of P. t. troglodytes, including the Mbuti Eastern Pygmies. Genomic regions showing signatures of selection were compared to the genomic locations of genes reported to be associated with HIV infection or pathogenesis. In the Biaka, a strong signal of selection was detected at CUL5, which codes for a component of the vif-mediated APOBEC3 degradation pathway. A CUL5 allele protective against AIDS progression was fixed in the Biaka. A signal of selection was detected at TRIM5, which codes for an HIV post-entry restriction factor. A protective mis-sense mutation in TRIM5 had the highest frequency in Biaka compared to other African populations, as did a protective allele for APOBEC3G, which codes for an anti-HIV-1 restriction factor. Alleles protective against HIV-1 for APOBEC3H, CXCR6 and HLA-C were at higher frequencies in the Biaka than in the Mbuti. Biaka genomes showed a strong signal of selection at TSG101, an inhibitor of HIV-1 viral budding.
We found protective alleles or evidence for selection in the Biaka at a number of genes associated with HIV-1 infection or progression. Pygmies have also been reported to carry genotypes protective against HIV-1 for the genes CCR5 and CCL3L1. Our hypothesis that HIV-1 may have shaped the genomes of some human populations in West Central Africa appears to merit further investigation.
HIV dependency factors; Single nucleotide polymorphisms; Biaka pygmies; Mbuti pygmies
Nasopharyngeal carcinoma (NPC) is an epithelial malignancy facilitated by Epstein-Barr Virus infection. Here we resolve the major genetic influences for NPC incidence using a genome-wide association study (GWAS), independent cohort replication, and high-resolution molecular HLA class I gene typing including 4,055 study participants from the Guangxi Zhuang Autonomous Region and Guangdong province of southern China. We detect and replicate strong association signals involving SNPs, HLA alleles, and amino acid (aa) variants across the major histocompatibility complex-HLA-A, HLA –B, and HLA -C class I genes (PHLA-A-aa-site-62 = 7.4×10−29; P HLA-B-aa-site-116 = 6.5×10−19; P HLA-C-aa-site-156 = 6.8×10−8 respectively). Over 250 NPC-HLA associated variants within HLA were analyzed in concert to resolve separate and largely independent HLA-A, -B, and -C gene influences. Multivariate logistical regression analysis collapsed significant associations in adjacent genes spanning 500 kb (OR2H1, GABBR1, HLA-F, and HCG9) as proxies for peptide binding motifs carried by HLA- A*11:01. A similar analysis resolved an independent association signal driven by HLA-B*13:01, B*38:02, and B*55:02 alleles together. NPC resistance alleles carrying the strongly associated amino acid variants implicate specific class I peptide recognition motifs in HLA-A and -B peptide binding groove as conferring strong genetic influence on the development of NPC in China.
NPC is a deadly throat cancer in China that is dependent on EBV infection. Here, we performed a 1 M SNP genome-wide association study using a large cohort of Chinese study participants at risk for NPC. Although several putative gene regions show significant associations, the strongest statistical signals involved scores of variants within the HLA region on chromosome 6. HLA poses a formidable association-genetics challenge because of extensive linkage disequilibrium, rather low allele frequencies, and multiple physically close interacting genes of diverse function. We examined over 250 NPC-HLA associated variants detected with sequence-based nucleotide alleles and amino acid variants. The multiple associations were collapsed to implicate causal signals by multivariate logistical regression to resolve allele association interaction. One operative variant was identified as the HLA-A*11:01 allele motif, specifically in the peptide binding groove, which recognizes invading antigens; a second involved two aa sites with HLA-B tracking B*13:01 and B*55:02 alleles. We synthesize these new and previous discoveries to help resolve the important gene influences on this disease.
In the absence of HIV-1 Vif, cellular cytosine deaminases such as APOBEC3G, inhibit the virus by inducing hypermutations on viral DNA, among other mechanisms of action. We investigated the association of APOBEC3G mRNA levels and APOBEC3G genetic variants on HIV-1 susceptibility, and early disease pathogenesis using viral load and CD4+ T cell counts as outcomes.
Study subjects were 250 South African females at high risk for HIV-1C infection. We used quantitative real-time PCR to measure the expression of APOBEC3G in HIV−ve and HIV+ve primary infection samples. APOBEC3G variants were identified by DNA re-sequencing and TaqMan genotyping.
We found no correlation between APOBEC3G expression levels and plasma viral loads (r=0.053, p=0.596) or CD4+ T cell counts (r=0.030, p=0.762) in 32 seroconverters. However, APOBEC3G expression levels were significantly higher in HIV−ve individuals compared to HIV+ve individuals (p<0.0001), including matched pre- and post infection samples from the same individuals (n=13, p<0.0001). 25 single nucleotide polymorphisms (SNPs), nine of which were novel, were identified within APOBEC3G by re-sequencing followed by genotyping of 168 individuals. The H186R mutation, a codon changing variant in exon 4, was associated with high viral loads (p=0.0097) and decreased CD4+ T cell levels (p=0.0081).
These data suggest that APOBEC3G transcription is rapidly downregulated upon HIV-1 infection. During primary infection, APOBEC3G expression levels in PBMCs do not correlate with viral loads or CD4+ T cell counts. However, structural variation of APOBEC3G may significantly affect early HIV-1 pathogenesis, although the mechanism remains unclear and warrants further investigation.
APOBEC3G; HIV-1 C; Primary Infection; mRNA Expression; Polymorphisms; Host proteins
LEDGF/p75, encoded by the PSIP1 gene, interacts with HIV-1 integrase and targets HIV-1 integration into active genes. We investigated the influence of polymorphisms in PSIP1 on HIV-1 acquisition and disease progression in black South Africans.
Integrase binding domain (IBD) of LEDGF/p75 was sequenced in 126 participants. Four haplotype tagging SNPs, SNP1-SNP4 and one exonic SNP, SNP5 were genotyped in 195 HIV-1 seronegative, 52 primary and 403 chronically infected individuals using TaqMan assays. LEDGF/p75 expression was quantified by real-time RT-PCR. The impact of Q472L mutation on the interaction with HIV-1 IN was measured by AlphaScreen.
rs2277191A was more frequent among seropositives (p=0.06, Fisher's exact test), and among individuals followed longitudinally trended towards association with higher likelihood of HIV-1 acquisition (RH=2.21, p=0.08; Cox model) and it was also associated with rapid disease progression (RH=5.98, p=0.04; Cox model). rs12339417C was associated with slower decline of CD4+ T cell (p=0.02) and lower levels of LEDGF/p75 (p<0.01). Seroconverters had higher preinfection levels of LEDGF/p75 (p<0.01) but levels decreased after HIV infection (p=0.02).
Genetic variants of PSIP1 may affect HIV-1 outcomes. Further studies are needed to confirm the effect of genetic variation of PSIP1 on HIV-1 pathogenesis in different cohorts.
Seven genome-wide association studies (GWAS) have been published in AIDS and only associations in the HLA region on chromosome 6 and CXCR6 have passed genome-wide significance.
We reanalyzed the data from three previously published GWAS, targeting specifically low frequency SNPs (minor allele frequency (MAF)<5%). Two groups composed of 365 slow progressors (SP) and 147 rapid progressors (RP) from Europe and the US were compared with a control group of 1394 seronegative individuals using Eigenstrat corrections.
Of the 8584 SNPs with MAF<5% in cases and controls (Bonferroni threshold=5.8×10−6), four SNPs showed statistical evidence of association with the SP phenotype. The best result was for HCP5 rs2395029 (p=8.54×10−15, OR=3.41) in the HLA locus, in partial linkage disequilibrium with two additional chromosome 6 associations in C6orf48 (p=3.03×10−10, OR=2.9) and NOTCH4 (9.08×10−07, OR=2.32). The fourth association corresponded to rs2072255 located in RICH2 (p=3.30×10−06, OR=0.43) in chromosome 17. Using HCP5 rs2395029 as a covariate, the C6orf48 and NOTCH4 signals disappeared, but the RICH2 signal still remained significant.
Besides the already known chromosome 6 associations, the analysis of low frequency SNPs brought up a new association in the RICH2 gene. Interestingly, RICH2 interacts with BST-2 known to be a major restriction factor for HIV-1 infection. Our study has thus identified a new candidate gene for AIDS molecular etiology and confirms the interest of singling out low frequency SNPs in order to exploit GWAS data.
AIDS; HIV-1; genome-wide association study; SNP; disease progression; RICH2
Background. Host genetic variation influences human immunodeficiency virus (HIV) infection and progression to AIDS. Here we used clinically well-characterized subjects from 5 pretreatment HIV/AIDS cohorts for a genome-wide association study to identify gene associations with rate of AIDS progression.
Methods. European American HIV seroconverters (n = 755) were interrogated for single-nucleotide polymorphisms (SNPs) (n = 700,022) associated with progression to AIDS 1987 (Cox proportional hazards regression analysis, co-dominant model).
Results. Association with slower progression was observed for SNPs in the gene PARD3B. One of these, rs11884476, reached genome-wide significance (relative hazard = 0.3; P =3. 370 × 10−9) after statistical correction for 700,022 SNPs and contributes 4.52% of the overall variance in AIDS progression in this study. Nine of the top-ranked SNPs define a PARD3B haplotype that also displays significant association with progression to AIDS (hazard ratio, 0.3; P = 3.220 × 10−8). One of these SNPs, rs10185378, is a predicted exonic splicing enhancer; significant alteration in the expression profile of PARD3B splicing transcripts was observed in B cell lines with alternate rs10185378 genotypes. This SNP was typed in European cohorts of rapid progressors and was found to be protective for AIDS 1993 definition (odds ratio, 0.43, P = .025).
Conclusions. These observations suggest a potential unsuspected pathway of host genetic influence on the dynamics of AIDS progression.
Diabetic nephropathy (DN) is a leading cause of mortality and morbidity in patients with type 1 and type 2 diabetes. The multicenter FIND consortium aims to identify genes for DN and its associated quantitative traits, e.g. the urine albumin:creatinine ratio (ACR). Herein, the results of whole-genome linkage analysis and a sparse association scan for ACR and a dichotomous DN phenotype are reported in diabetic individuals.
A genomewide scan comprising more than 5,500 autosomal single nucleotide polymorphism markers (average spacing of 0.6 cM) was performed on 1,235 nuclear and extended pedigrees (3,972 diabetic participants) ascertained for DN from African-American (AA), American-Indian (AI), European-American (EA) and Mexican-American (MA) populations.
Strong evidence for linkage to DN was detected on chromosome 6p (p = 8.0 × 10−5, LOD = 3.09) in EA families as well as suggestive evidence for linkage to chromosome 7p in AI families. Regions on chromosomes 3p in AA, 7q in EA, 16q in AA and 22q in MA displayed suggestive evidence of linkage for urine ACR. The linkage peak on chromosome 22q overlaps the MYH9/APOL1 gene region, previously implicated in AA diabetic and nondiabetic nephropathies.
These results strengthen the evidence for previously identified genomic regions and implicate several novel loci potentially involved in the pathogenesis of DN.
Albuminuria; Diabetes mellitus; Renal failure; End-stage renal disease; Linkage; Allelic association
A recent genome-wide study showed that the single nucleotide polymorphisms (SNPs) in the HLA-DP region were associated with chronic hepatitis B virus (HBV) infection in Japanese and Thai persons. We tested the effects of HLA-DP SNPs for all major HBV outcomes in Han Chinese (n = 1742): HBV resistance, clearance, chronic infection, cirrhosis, and hepatocellular carcinoma. HLA - DPA1 rs3077 T was strongly associated with decreased risk of chronic HBV infection (odds ratio, .62; P = .001), consistent with the previous report. We showed for the first time to our knowledge that it is a predictor for HBV clearance (odds ratio, 2.41; P < .001). However, rs3077 was not associated with the development of cirrhosis or hepatocellular carcinoma.
Chromosome 3p21–22 harbors two clusters of chemokine receptor genes, several of which serve as major or minor coreceptors of HIV-1. Although the genetic association of CCR5 and CCR2 variants with HIV-1 pathogenesis is well known, the role of variation in other nearby chemokine receptor genes remain unresolved. We genotyped exonic single nucleotide polymorphisms (SNPs) in chemokine receptor genes: CCR3, CCRL2, and CXCR6 (at 3p21) and CCR8 and CX3CR1 (at 3p22), the majority of which were non-synonymous. The individual SNPs were tested for their effects on disease progression and outcomes in five treatment-naïve HIV-1/AIDS natural history cohorts. In addition to the known CCR5 and CCR2 associations, significant associations were identified for CCR3, CCR8, and CCRL2 on progression to AIDS. A multivariate survival analysis pointed to a previously undetected association of a non-conservative amino acid change F167Y in CCRL2 with AIDS progression: 167F is associated with accelerated progression to AIDS (RH = 1.90, P = 0.002, corrected). Further analysis indicated that CCRL2-167F was specifically associated with more rapid development of pneumocystis pneumonia (PCP) (RH = 2.84, 95% CI 1.28–6.31) among four major AIDS–defining conditions. Considering the newly defined role of CCRL2 in lung dendritic cell trafficking, this atypical chemokine receptor may affect PCP through immune regulation and inducing inflammation.
Human chemokine receptors are cell surface proteins that may be utilized by HIV-1 for entry into host cells. DNA variation in the HIV-1 major coreceptor CCR5 affects HIV-1 infection and progression. This study comprehensively assesses the role of genetic variation of multiple chemokine receptor genes clustered in the chromosome 3p21 and 3p22 on HIV-1 disease outcomes in HIV-1 natural history cohorts. The multivariate survival analyses identified functional variants that altered disease progression rate in CCRL2, CCR3, and CCR8. CCRL2-F167Y affects the rate to AIDS development through a specific protection against pneumocystis pneumonia (PCP), a common AIDS–defining condition. Our study identified this atypical chemokine receptor CCRL2 as a key factor involved in PCP, possibly through inducing inflammation in the lung.
Background. The compilation of previous genomewide association studies of AIDS shows a major polymorphism in the HCP5 gene associated with both control of the viral load and long-term nonprogression (LTNP) to AIDS.
Methods. To look for genetic variants that affect LTNP without necessary control of the viral load, we reanalyzed the genomewide data of the unique LTNP Genomics of Resistance to Immunodeficiency Virus (GRIV) cohort by excluding “elite controller” patients, who were controlling the viral load at very low levels (<100 copies/mL).
Results. The rs2234358 polymorphism in the CXCR6 gene was the strongest signal (P = 2.5 × 10−7; odds ratio, 1.85) obtained for the genomewide association study comparing the 186 GRIV LTNPs who were not elite controllers with 697 uninfected control subjects. This association was replicated in 3 additional independent European studies, reaching genomewide significance of Pcombined = 9.7 × 10−10. This association with LTNP is independent of the combined CCR2-CCR5 locus and the HCP5 polymorphisms.
Conclusion. The statistical significance, the replication, and the magnitude of the association demonstrate that CXCR6 is likely involved in the molecular etiology of AIDS and, in particular, in LTNP, emphasizing the power of extreme-phenotype cohorts. CXCR6 is a chemokine receptor that is known as a minor coreceptor in human immunodeficiency virus type 1 infection but could participate in disease progression through its role as a mediator of inflammation.
African Americans (AA) disproportionately develop lupus nephritis (LN) relative to European Americans and familial clustering supports causative genes. Since MYH9 underlies approximately 40% of end-stage renal disease (ESRD) in AA, we tested for genetic association with LN.
Seven MYH9 single nucleotide polymorphisms (SNPs) and the E1 risk haplotype were tested for association with LN in three cohorts of AA.
A preliminary analysis revealed that the MYH9 E1 risk haplotype was associated with ESRD in 25 cases with presumed systemic lupus erythematosus (SLE)-associated ESRD, compared to 735 non-SLE controls (odds ratio 3.1; p = 0.010 recessive). Replication analyses were performed in 583 AA with SLE in the PROFILE cohort (318 with LN; 265 with SLE but without nephropathy) and 60 AA from the NIH (39 with LN; 21 with SLE but without nephropathy). Analysis of the NIH and larger PROFILE cohorts, as well as a combined analysis, did not support this association.
These results suggest that AA with ESRD and coincident SLE who were recruited from dialysis clinics more likely have kidney diseases in the MYH9-associated spectrum of focal segmental glomerulosclerosis. PROFILE and NIH participants, recruited from rheumatology practices, demonstrate that MYH9 does not contribute substantially to the development of LN in AA.
African Americans; Genetics; Lupus nephritis; Kidney; MYH9; Systemic lupus erythematosus
Single nucleotide polymorphisms (SNPs) in MYH9 and APOL1 on chromosome 22 (c22) are powerfully associated with non-diabetic end-stage renal disease (ESRD) in African Americans (AAs). Many AAs diagnosed with type 2 diabetic nephropathy (T2DN) have non-diabetic kidney disease, potentially masking detection of DN genes. Therefore, genome-wide association analyses were performed using the Affymetrix SNP Array 6.0 in 966 AA with T2DN and 1,032 non-diabetic, non-nephropathy (NDNN) controls, with and without adjustment for c22 nephropathy risk variants. No associations were seen between FRMD3 SNPs and T2DN before adjusting for c22 variants. However, logistic regression analysis revealed seven FRMD3 SNPs significantly interacting with MYH9—a finding replicated in 640 additional AA T2DN cases and 683 NDNN controls. Contrasting all 1,592 T2DN cases with all 1,671 NDNN controls, FRMD3 SNPs appeared to interact with the MYH9 E1 haplotype (e.g., rs942280 interaction p-value = 9.3E−7 additive; odds ratio [OR] 0.67). FRMD3 alleles were associated with increased risk of T2DN only in subjects lacking two MYH9 E1 risk haplotypes (rs942280 OR = 1.28), not in MYH9 E1 risk allele homozygotes (rs942280 OR = 0.80; homogeneity p-value = 4.3E−4). Effects were weaker stratifying on APOL1. FRMD3 SNPS were associated with T2DN, not type 2 diabetes per se, comparing AAs with T2DN to those with diabetes lacking nephropathy. T2DN-associated FRMD3 SNPs were detectable in AAs only after accounting for MYH9, with differential effects for APOL1. These analyses reveal a role for FRMD3 in AA T2DN susceptibility and accounting for c22 nephropathy risk variants can assist in detecting DN susceptibility genes.
African Americans have high rates of kidney disease attributed to type 2 diabetes mellitus. However, approximately 25% of patients are misclassified and have non-diabetic kidney disease on renal biopsy. The APOL1-MYH9 gene region on chromosome 22 is powerfully associated with non-diabetic kidney diseases in African Americans. Therefore, we tested for interactions between single nucleotide polymorphisms across the genome with APOL1 and MYH9 non-diabetic nephropathy risk variants in African Americans with presumed diabetic nephropathy. Markers in FRMD3, a gene associated with type 1 diabetic nephropathy in Caucasians, appeared to interact with MYH9; however, increased nephropathy risk was seen in diabetic cases lacking two MYH9 risk haplotypes, and protective effects were seen in those with two MYH9 risk haplotypes. Stratified analyses based on the chromosome 22 nephropathy risk haplotypes demonstrated that FRMD3 variants were associated with diabetic nephropathy risk in cases without two MYH9 (or APOL1) risk haplotypes. It appears that African Americans with diabetes and kidney disease who are not chromosome 22 nephropathy risk variant homozygotes are enriched for the presence of diabetic nephropathy and FRMD3 risk alleles. This genetic dissection ultimately allowed for detection of the FRMD3 diabetic nephropathy gene association in a subset of cases enriched for this disorder.
Genetic variation in MYH9, encoding nonmuscle myosin IIA heavy chain, has been associated recently with increased risk for kidney disease. Previously, MYH9 missense mutations have been shown to cause the autosomal-dominant MYH9 (ADM9) spectrum, characterized by large platelets, leukocyte Döhle bodies, and, variably, sensorineural deafness, cataracts, and glomerulopathy. Genetic testing is indicated for familial and sporadic cases that fit this spectrum. By contrast, the MYH9 kidney risk variant is characterized by multiple intronic single nucleotide polymorphisms, but the causative variant has not been identified. Disease associations include human immunodeficiency virus-associated collapsing glomerulopathy, focal segmental glomerulosclerosis, hypertension-attributed end-stage kidney disease, and diabetes-attributed end-stage kidney disease. One plausible hypothesis is that the MYH9 kidney risk variant confers a fragile podocyte phenotype. In the case of hypertension-attributed kidney disease, it remains unclear if the hypertension is a contributing cause or a consequence of glomerular injury. The MYH9 kidney risk variant is strikingly more common among individuals of African descent, but only some will develop clinical kidney disease in their lifetime. Thus, it is likely that additional genes and/or environmental factors interact with the MYH9 kidney risk variant to trigger glomerular injury. A preliminary genetic risk stratification scheme, using 10 single nucleotide polymorphisms, may estimate lifetime risk for kidney disease. Nevertheless, at present, no role has been established for genetic testing as part of personalized medicine, but testing should be considered in clinical studies of glomerular diseases among populations of African descent. Such studies will address critical questions pertaining to MYH9-associated kidney disease, including mechanism, course, and response to therapy.
Focal segmental glomerulosclerosis; HIV-associated nephropathy; hypertensive nephrosclerosis; chronic kidney disease; end-stage kidney disease; African American
Until recently knowledge of genetic causes of glomerular disease was limited to certain rare or uncommon inherited diseases, and to a genes, either rare or with small effect, identified in candidate gene studies. These genetic factors accounted for only a very small fraction of kidney disease. However, the striking differences in frequency of many forms of kidney disease between African Americans and European Americans, which could not be completely explained by cultural or economic factors, pointed to a large unidentified genetic influence. Since FSGS and HIV-associated collapsing glomerulopathy (HVAN) have striking racial disparities, we performed an admixture mapping study to identify contributing genetic factors. Admixture mapping identified genetic variants in the non-muscle myosin gene MYH9 as having an extreme influence on both FSGS and HIVAN, with odds ratios from 4 to 8 and attributable fractions of 70–100%. Previously identified, rare inherited MYH9 disorders point to a mechanism by which MYH9 variation disrupts the actin-myosin filaments responsible for maintaining the structure of podocytes, the cells that provide one of three filtration barriers in the glomeruli. MYH9 variation has a smaller but still highly significant effect on non-diabetic kidney disease, and a weaker but significant effect on diabetic kidney disease; it is unclear whether underlying cryptic FSGS is responsible for the MYH9 association with these diseases. The strong predicted power of MYH9 variation for disease indicates a clear role for genetic testing for these variants in personalized medicine, for assessment of genetic risk, and potentially for diagnosis.
As we enter an era when testing millions of SNPs in a single gene association study will become the standard, consideration of multiple comparisons is an essential part of determining statistical significance. Bonferroni adjustments can be made but are conservative due to the preponderance of linkage disequilibrium (LD) between genetic markers, and permutation testing is not always a viable option. Three major classes of corrections have been proposed to correct the dependent nature of genetic data in Bonferroni adjustments: permutation testing and related alternatives, principal components analysis (PCA), and analysis of blocks of LD across the genome. We consider seven implementations of these commonly used methods using data from 1514 European American participants genotyped for 700,078 SNPs in a GWAS for AIDS.
A Bonferroni correction using the number of LD blocks found by the three algorithms implemented by Haploview resulted in an insufficiently conservative threshold, corresponding to a genome-wide significance level of α = 0.15 - 0.20. We observed a moderate increase in power when using PRESTO, SLIDE, and simpleℳ when compared with traditional Bonferroni methods for population data genotyped on the Affymetrix 6.0 platform in European Americans (α = 0.05 thresholds between 1 × 10-7 and 7 × 10-8).
Correcting for the number of LD blocks resulted in an anti-conservative Bonferroni adjustment. SLIDE and simpleℳ are particularly useful when using a statistical test not handled in optimized permutation testing packages, and genome-wide corrected p-values using SLIDE, are much easier to interpret for consumers of GWAS studies.
African-Americans have higher rates of kidney disease than European-Americans. Here we show that in African-Americans, focal segmental glomerulosclerosis (FSGS) and hypertension-attributed end-stage kidney disease (H-ESKD) are associated with two independent sequence variants in the APOL1 gene on chromosome 22 [FSGS odds ratio = 10.5 (95% CI 6.0–18.4); H-ESKD odds ratio = 7.3 (95% CI 5.6–9.5)]. The two APOL1 variants are common in African chromosomes but absent from European chromosomes and both reside within haplotypes that harbor signatures of positive selection. ApoL1 is a serum factor that lyses trypanosomes. In vitro assays revealed that only the kidney disease-associated ApoL1 variants lysed Trypanosoma brucei rhodesiense. We speculate that evolution of a critical survival factor in Africa may have contributed to the high rates of renal disease in African-Americans.
Background. Although MYH9 is strongly associated with biopsy-proven idiopathic and HIV-associated focal segmental glomerulosclerosis (FSGS) and clinically diagnosed ‘hypertension-associated’ end-stage renal disease (ESRD) in African Americans, its role in type 2 diabetes mellitus (T2DM)-associated ESRD is unclear.
Methods. To assess whether MYH9 was associated with T2DM-ESRD, 751 African Americans with T2DM-ESRD, 227 with T2DM lacking nephropathy and 925 non-diabetic non-nephropathy controls were genotyped for 14 MYH9 SNPs. Association analyses used SNPGWA and Dandelion.
Results. Comparing T2DM-ESRD cases with non-diabetic controls, single SNP associations were detected with 8 of 14 SNPs, gender- and admixture-adjusted P-values 0.047–0.005 [recessive model, odds ratio (OR) range 1.30–1.55]. The previously associated MYH9 E1 and L1 haplotypes were associated with T2DM-ESRD (E1: OR 1.27, 95% CI 1.04–1.56, P = 0.021 recessive and L1: OR 1.43, 95% CI 1.09–1.87, P = 0.009 dominant). Contrasting the 751 T2DM-ESRD cases with 227 T2DM non-nephropathy controls revealed that E1 haplotype SNPs rs4821480, rs2032487 and rs4821481 were associated with kidney failure (OR 1.38–1.40 recessive, all P < 0.048). Among E1 and L1 risk homozygotes, respectively, mean (SD) diabetes duration prior to renal replacement therapy was 16.6 (9.7) and 16.4 (10.0) years, and 65% had diabetic retinopathy.
Conclusions. Genetic dissection of T2DM-associated ESRD reveals that MYH9 underlies a portion of this clinically diagnosed disorder in African Americans. It is likely that a subset of African Americans with T2DM and coincident nephropathy have primary MYH9-related kidney disease (e.g. FSGS or global glomerulosclerosis), although renal biopsy studies need to be performed.
African American; diabetic nephropathy; kidney; MYH9; type 2 diabetes mellitus