Response to Bakker et al.

Robinson and colleagues respond to the points raised about their paper by Bakker et al.

Genome-wide diversity in the California condor tracks its prehistoric abundance and decline.

Due to their small population sizes, threatened and endangered species frequently suffer from a lack of genetic diversity, potentially leading to inbreeding depression and reduced adaptability.1 During the latter half of the twentieth century, North America's largest soaring bird,2 the California condor (Gymnogyps californianus; Critically Endangered3), briefly went extinct in the wild. Though condors once ranged throughout North America, by 1982 only 22 individuals remained. Following decades of captive breeding and release efforts, there are now >300 free-flying wild condors and ∼200 in captivity. The condor's recent near-extinction from lead poisoning, poaching, and loss of habitat is well documented,4 but much about its history remains obscure. To fill this gap and aid future management of the species, we produced a high-quality chromosome-length genome assembly for the California condor and analyzed its genome-wide diversity. For comparison, we also examined the genomes of two close relatives: the Andean condor (Vultur gryphus; Vulnerable3) and the turkey vulture (Cathartes aura; Least Concern3). The genomes of all three species show evidence of historic population declines. Interestingly, the California condor genome retains a high degree of variation, which our analyses reveal is a legacy of its historically high abundance. Correlations between genome-wide diversity and recombination rate further suggest a history of purifying selection against linked deleterious alleles, boding well for future restoration. We show how both long-term evolutionary forces and recent inbreeding have shaped the genome of the California condor, and provide crucial genomic resources to enable future research and conservation.

A genome wide study of genetic adaptation to high altitude in feral Andean Horses of the páramo.

BACKGROUND: Life at high altitude results in physiological and metabolic challenges that put strong evolutionary pressure on performance due to oxidative stress, UV radiation and other factors dependent on the natural history of the species. To look for genes involved in altitude adaptation in a large herbivore, this study explored genome differentiation between a feral population of Andean horses introduced by the Spanish in the 1500s to the high Andes and their Iberian breed relatives. RESULTS: Using allelic genetic models and Fst analyses of ~50 K single nucleotide polymorphisms (SNPs) across the horse genome, 131 candidate genes for altitude adaptation were revealed (Bonferoni of p ≤ 2 × 10(-7)). Significant signals included the EPAS1 in the hypoxia-induction-pathway (HIF) that was previously discovered in human studies (p = 9.27 × 10(-8)); validating the approach and emphasizing the importance of this gene to hypoxia adaptation. Strong signals in the cytochrome P450 3A gene family (p = 1.5 ×10(-8)) indicate that other factors, such as highly endemic vegetation in altitude environments are also important in adaptation. Signals in tenuerin 2 (TENM2, p = 7.9 × 10(-14)) along with several other genes in the nervous system (gene categories representation p = 5.1 × 10(-5)) indicate the nervous system is important in altitude adaptation. CONCLUSIONS: In this study of a large introduced herbivore, it becomes apparent that some gene pathways, such as the HIF pathway are universally important for high altitude adaptation in mammals, but several others may be selected upon based on the natural history of a species and the unique ecology of the altitude environment.

Host genomic influences on HIV/AIDS.

The AIDS era has seen multiple advances in the power of genetics research; scores of host genetic protective factors have been nominated and several have translated to the bedside. We discuss how genomics may inform HIV/AIDS prevention, treatment and eradication.

The principal genetic determinants for nasopharyngeal carcinoma in China involve the HLA class I antigen recognition groove.

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 (P(HLA-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.

Risk factors for symptomatic hyperlactatemia and lactic acidosis among combination antiretroviral therapy-treated adults in Botswana: results from a clinical trial.

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.

Genome-wide association study implicates PARD3B-based AIDS restriction.

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.

Genetic variants in nuclear-encoded mitochondrial genes influence AIDS progression.

BACKGROUND: The human mitochondrial genome includes only 13 coding genes while nuclear-encoded genes account for 99% of proteins responsible for mitochondrial morphology, redox regulation, and energetics. Mitochondrial pathogenesis occurs in HIV patients and genetically, mitochondrial DNA haplogroups with presumed functional differences have been associated with differential AIDS progression. METHODOLOGY/PRINCIPAL FINDINGS: Here we explore whether single nucleotide polymorphisms (SNPs) within 904 of the estimated 1,500 genes that specify nuclear-encoded mitochondrial proteins (NEMPs) influence AIDS progression among HIV-1 infected patients. We examined NEMPs for association with the rate of AIDS progression using genotypes generated by an Affymetrix 6.0 genotyping array of 1,455 European American patients from five US AIDS cohorts. Successfully genotyped SNPs gave 50% or better haplotype coverage for 679 of known NEMP genes. With a Bonferroni adjustment for the number of genes and tests examined, multiple SNPs within two NEMP genes showed significant association with AIDS progression: acyl-CoA synthetase medium-chain family member 4 (ACSM4) on chromosome 12 and peroxisomal D3,D2-enoyl-CoA isomerase (PECI) on chromosome 6. CONCLUSIONS: Our previous studies on mitochondrial DNA showed that European haplogroups with presumed functional differences were associated with AIDS progression and HAART mediated adverse events. The modest influences of nuclear-encoded mitochondrial genes found in the current study add support to the idea that mitochondrial function plays a role in AIDS pathogenesis.

Effect of host genetics on the development of cytomegalovirus retinitis in patients with AIDS.

BACKGROUND: Cytomegalovirus (CMV) retinitis is a common opportunistic infection among patients with AIDS and still causes visual morbidity despite the wide spread usage of highly active antiretroviral therapy (HAART). The ubiquitous CMV pathogen contains a human interleukin-10 (IL-10) homolog in its genome and utilizes it to evade host immune reactions through an IL-10 receptor mediated immune-suppression pathway. METHODS: Effects of IL-10R1, IL-10 and previously described AIDS restriction gene variants are investigated on the development of CMV retinitis in the Longitudinal Study of the Ocular Complications of AIDS (LSOCA) cohort (N = 1284). RESULTS: In European Americans (n = 750), a haplotype carrying an amino acid changing variation in the cytoplasmic domain (S420L) of IL-10R1 can be protective (OR, 0.14; 95% CI, 0.02-0.94; P = .04) against, whereas another haplotype carrying an amino acid changing variation in the extracellular domain (I224V) of IL-10R1 can be more susceptible (OR, 6.21; 95% CI, 1.22- 31.54; P = .03) to CMV retinitis. In African Americans (n = 534), potential effects of IL-10 variants are observed. CONCLUSION: Host genetics may have a role in the occurrence of CMV retinitis in patients infected with HIV.

Effect of host genetics on incidence of HIV neuroretinal disorder in patients with AIDS.

Approximately 10%-15% of patients with AIDS but without ocular opportunistic infections will have a presumed neuroretinal disorder (HIV-NRD), manifested by reduced contrast sensitivity and abnormal visual fields. The loss of contrast sensitivity often is sufficient to impair reading speed. To evaluate the effect of host genetics on HIV-NRD, we explored validated AIDS restriction gene variants CCR5Delta32, CCR2-64I, CCR5 P1, SDF-3'A, IL-10-5'A, RANTES -403A, RANTES -28G, RANTES-In1.1C, CX3CR1-249I, CX3CR1-280M, IFNG-179T, MDR1-3435T, and MCP-1364G, each of which has been implicated previously to influence HIV-1 infection, AIDS progression, therapy response, and antiviral drug metabolism, and an IL-10 receptor gene, IL-10R1, in the Longitudinal Study of the Ocular Complications of AIDS cohort. In European Americans (cases = 55, controls = 290), IL-10-5'A variant and its promoter haplotype (hazard ratio = 2.09, confidence interval. 1.19 to 3.67, P = 0.01), in African Americans (cases = 54, controls = 180), RANTES-In1.1C and the associated haplotype (hazard ratio = 2.72, confidence interval.: 1.48 to 5.00, P = 0.001), showed increased HIV-NRD susceptibility. Although sample sizes are small and P values do not pass a strict Bonferroni correction, our results suggest that, in European Americans, an IL-10-related pathway, and, in African Americans, chemokine receptor ligand polymorphisms in RANTES are risk factors for HIV-NRD development. Clearly, further studies are warrented.

Mitochondrial haplogroups are associated with risk of neuroretinal disorder in HIV-positive patients.

Although highly active antiretroviral therapy has improved survivorship dramatically and decreased the incidence of cytomegalovirus retinitis among patients with AIDS, other ophthalmic complications continue to occur. One complication observed in approximately 12% of HIV-infected patients is a presumed neuroretinal disorder (NRD), manifested as decreased contrast sensitivity and associated with vague subjective complaints of hazy vision. Pathologically, patients with AIDS even without ocular opportunistic infections have loss of optic nerve axons, suggestive of mitochondrial dysfunction. We explored whether variation in mitochondrial DNA was associated with time to NRD in HIV-infected patients in the Longitudinal Study of Ocular Complications of AIDS cohort. Within the Western European, or "N", mitochondrial DNA macrohaplogroup, haplogroup J, was associated with 80% decrease in the risk of progression to NRD during the study (hazard ratio = 0.20, P = 0.039) and suggested an independent association with protection against NRD in a cross-section of all patients taken at enrollment (1.5% vs. 8.9% in patients with vs. without haplogroup J, respectively, P = 0.05). These data suggest that mitochondrial genotype may influence propensity to develop HIV-associated NRD in patients with AIDS.

Mitochondrial DNA haplogroups influence lipoatrophy after highly active antiretroviral therapy.

Although highly active antiretroviral therapy (HAART) has been extremely effective in lowering AIDS incidence among patients infected with HIV, certain drugs included in HAART can cause serious mitochondrial toxicities. One of the most frequent adverse events is lipoatrophy, which is the loss of subcutaneous fat in the face, arms, buttocks, and/or legs as an adverse reaction to nucleoside reverse transcriptase inhibitors. The clinical symptoms of lipoatrophy resemble those of inherited mitochondrial diseases, which suggest that host mitochondrial genotype may play a role in susceptibility. We analyzed the association between mitochondrial haplogroup and severity of lipoatrophy in HIV-infected European American patients on HAART in the Multicenter AIDS cohort Study and found that mitochondrial haplogroup H was strongly associated with increased atrophy [arms: P = 0.007, odds ratio (OR) = 1.77, 95% confidence interval (CI) = 1.17 to 2.69; legs: P = 0.037, OR = 1.54, 95% CI = 1.03 to 2.31; and buttocks: P = 0.10, OR = 1.41 95% CI = 0.94 to 2.12]. We also saw borderline significance for haplogroup T as protective against lipoatrophy (P = 0.05, OR = 0.52, 95% CI = 0.20 to 1.00). These data suggest that mitochondrial DNA haplogroup may influence the propensity for lipoatrophy in patients receiving nucleoside reverse transcriptase inhibitors.

Mitochondrial DNA haplogroups influence AIDS progression.

OBJECTIVE: Mitochondrial function plays a role in both AIDS progression and HAART toxicity; therefore, we sought to determine whether mitochondrial DNA variation revealed novel AIDS restriction genes, particularly as mitochondrial DNA single-nucleotide polymorphisms are known to influence regulation of oxidative phosphorylation, reactive oxygen species production, and apoptosis. DESIGN: This is a retrospective cohort study. METHODS: We performed an association study of mitochondrial DNA haplogroups among 1833 European American HIV-1 patients from five US cohorts: the Multicenter AIDS Cohort Study, the San Francisco City Clinic Study, Hemophilia Growth and Development Study, the Multicenter Hemophilia Cohort Study, and the AIDS Linked to Intravenous Experiences cohort to determine whether the mitochondrial DNA haplogroup correlated with AIDS progression rate. RESULTS: Mitochondrial DNA haplogroups J and U5a were elevated among HIV-1 infected people who display accelerated progression to AIDS and death. Haplogroups Uk, H3, and IWX appeared to be highly protective against AIDS progression. CONCLUSION: The associations found in our study appear to support a functional explanation by which mitochondrial DNA variation among haplogroups, influencing ATP production, reactive oxygen species generation, and apoptosis, is correlated to AIDS disease progression; however, repeating these results in cohorts with different ethnic backgrounds would be informative. These data suggest that mitochondrial genes are important indicators of AIDS disease progression in HIV-1 infected persons.

Host genetic influences on highly active antiretroviral therapy efficacy and AIDS-free survival.

We studied the influence of AIDS restriction genes (ARGs) CCR5-Delta32, CCR2-64I, SDF1-3'A, IL10-5'A, CX3CR1-V249I, CX3CR1-T280M, and MDR1-C3435T and haplotypes of the CCR5 P1 promoter and RANTES variants -403A, In1.1C, 3'222C, and -28G among HIV-1 infected patients on highly active antiretroviral therapy (HAART) in the Multicenter AIDS Cohort Study (MACS) and the Multicenter Hemophilia Cohort Study (MHCS). Our results indicate that several ARGs also influence therapy efficacy (ie, the success in viral suppression) and subsequent progression to AIDS while on HAART. CCR5-Delta32 decreased time to viral suppression (<200 HIV RNA copies/mL, relative hazard [RH]=1.40; P=0.008) and was protective against AIDS (RH=0.11; P=or<0.0001), whereas the CCR5 P1 haplotype was associated with delayed viral suppression (RNA<50 copies/mL, odds ratio [OR]=0.65; P=0.03) and accelerated time to AIDS (RH=2.68; P=0.02). SDF1-3'A reduced viral suppression (OR=0.61; P=0.02) and accelerated AIDS (RH=3.18; P=0.009). Accelerated AIDS progression was also observed with the RANTES haplotype carrying RANTES-IN1.1C and RANTES-3'222C (P=0.005 to 0.007). In contrast, the RANTES haplotype H1, which lacks suspected deleterious single-nucleotide polymorphisms, was protective against AIDS. CX3CR1-V249I seemed to accelerate viral suppression (RNA<50 copies/mL, OR=1.27; P=0.01). ARG influence after HAART suggests residual HIV-1 replication, and spread continues even in patients successfully suppressing detectable viral RNA.

NPHS2 variation in sporadic focal segmental glomerulosclerosis.

Mutations in NPHS2, the gene that encodes podocin, are well-established causes of both familial and sporadic steroid-resistant focal segmental glomerulosclerosis (FSGS) in the pediatric population, but have not been well-characterized in late-onset disease. To investigate the role of NPHS2 polymorphisms in sporadic cases of late-onset FSGS, we studied 377 biopsy-confirmed FSGS cases and 919 controls. We identified 18 single nucleotide polymorphisms (SNPs) by resequencing a subgroup of cases and controls, and subsequently genotyped African-American and European-American cases and controls for five missense SNPs, three SNPs within introns, and four SNPs in the 3' untranslated region. No homozygotes or compound heterozygotes were observed for any missense mutation. R138Q carriers were more frequent among FSGS cases relative to controls (OR = 4.9, P = 0.06), but heterozygosity for the other four missense mutations was equally distributed among FSGS cases and controls. Finally, a common haplotype of noncoding SNPs carried by 20% of African-Americans, but not observed in European-Americans, was strongly associated with a 50% reduction in risk for sporadic FSGS (OR = 0.5, P = 0.001). These results indicate that genetic variation or mutation of NPHS2 may play a role in late-onset sporadic FSGS.