Association Study of Common Genetic Variants and HIV-1 Acquisition in 6,300 Infected Cases and 7,200 Controls

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⁻¹¹). 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.     

Passive sexual transmission of human immunodeficiency virus type 1 variants and adaptation in new hosts

Human immunodeficiency virus type 1 (HIV-1) genetic diversity is a major obstacle for the design of a successful vaccine. Certain viral polymorphisms encode human leukocyte antigen (HLA)-associated immune escape, potentially overcoming limited vaccine protection. Although transmission of immune escape variants has been reported, the overall extent to which this phenomenon occurs in populations and the degree to which it contributes to HIV-1 viral evolution are unknown. Selection on the HIV-1 env gene at transmission favors neutralization-sensitive variants, but it is not known to what degree selection acts on the internal HIV-1 proteins to restrict or enhance the transmission of immune escape variants. Studies have suggested that HLA class I may determine susceptibility to HIV-1 infection, but a definitive role for HLA at transmission remains unproven. Comparing populations of acute seroconverters and chronically infected patients, we found no evidence of selection acting to restrict transmission of HIV-1 variants. We found that statistical associations previously reported in chronic infection between viral polymorphisms and HLA class I alleles are not present in acute infection, suggesting that the majority of viral polymorphisms in these patients are the result of transmission rather than de novo adaptation. Using four episodes of HIV-1 transmission in which the donors and recipients were both sampled very close to the time of infection we found that, despite a transmission bottleneck, genetic variants of HIV-1 infection are transmitted in a frequency-dependent manner. As HIV-1 infections are seeded by unique donor-adapted viral variants, each episode is a highly individual antigenic challenge. Host-specific, idiosyncratic HIV-1 antigenic diversity will seriously tax the efficacy of immunization based on consensus sequences.     

A gene co-expression network in whole blood of schizophrenia patients is independent of antipsychotic-use and enriched for brain-expressed genes

Despite large-scale genome-wide association studies (GWAS), the underlying genes for schizophrenia are largely unknown. Additional approaches are therefore required to identify the genetic background of this disorder. Here we report findings from a large gene expression study in peripheral blood of schizophrenia patients and controls. We applied a systems biology approach to genome-wide expression data from whole blood of 92 medicated and 29 antipsychotic-free schizophrenia patients and 118 healthy controls. We show that gene expression profiling in whole blood can identify twelve large gene co-expression modules associated with schizophrenia. Several of these disease related modules are likely to reflect expression changes due to antipsychotic medication. However, two of the disease modules could be replicated in an independent second data set involving antipsychotic-free patients and controls. One of these robustly defined disease modules is significantly enriched with brain-expressed genes and with genetic variants that were implicated in a GWAS study, which could imply a causal role in schizophrenia etiology. The most highly connected intramodular hub gene in this module (ABCF1), is located in, and regulated by the major histocompatibility (MHC) complex, which is intriguing in light of the fact that common allelic variants from the MHC region have been implicated in schizophrenia. This suggests that the MHC increases schizophrenia susceptibility via altered gene expression of regulatory genes in this network.     

Recombination and mutation shape variations in the major histocompatibility complex

The major histocompatibility complex (MHC) is closely associated with numerous diseases, but its high degree of polymorphism complicates the discovery of disease-associated variants. In principle, recombination and de novo mutations are two critical factors responsible for MHC polymorphisms. However, direct evidence for this hypothesis is lacking. Here, we report the generation of fine-scale MHC recombination and de novo mutation maps of ～5 Mb by deep sequencing (> 100×) of the MHC genome for 17 MHC recombination and 30 non-recombination Han Chinese families (a total of 190 individuals). Recombination hotspots and Han-specific breakpoints are located in close proximity at haplotype block boundaries. The average MHC de novo mutation rate is higher than the genome-wide de novo mutation rate, particularly in MHC recombinant individuals. Notably, mutation and recombination generated polymorphisms are located within and outside linkage disequilibrium regions of the MHC, respectively, and evolution of the MHC locus was mainly controlled by positive selection. These findings provide insights on the evolutionary causes of the MHC diversity and may facilitate the identification of disease-associated genetic variants.     

Use of a combined ex vivo/in vivo population approach for screening of human genes involved in the human immunodeficiency virus type 1 life cycle for variants influencing disease progression

Humans differ substantially with respect to susceptibility to human immunodeficiency virus type 1 (HIV-1). We evaluated variants of nine host genes participating in the viral life cycle for their role in modulating HIV-1 infection. Alleles were assessed ex vivo for their impact on viral replication in purified CD4 T cells from healthy blood donors (n = 128). Thereafter, candidate alleles were assessed in vivo in a cohort of HIV-1-infected individuals (n = 851) not receiving potent antiretroviral therapy. As a benchmark test, we tested 12 previously reported host genetic variants influencing HIV-1 infection as well as single nucleotide polymorphisms in the nine candidate genes. This led to the proposition of three alleles of PML, TSG101, and PPIA as potentially associated with differences in progression of HIV-1 disease. In a model considering the combined effects of new and previously reported gene variants, we estimated that their effect might be responsible for lengthening or shortening by up to 2.8 years the period from 500 CD4 T cells/mul to <200 CD4 T cells/mul.     

Genome-Wide Association Scan for Variants Associated with Early-Onset Prostate Cancer

Prostate cancer is the most common non-skin cancer and the second leading cause of cancer related mortality for men in the United States. There is strong empirical and epidemiological evidence supporting a stronger role of genetics in early-onset prostate cancer. We performed a genome-wide association scan for early-onset prostate cancer. Novel aspects of this study include the focus on early-onset disease (defined as men with prostate cancer diagnosed before age 56 years) and use of publically available control genotype data from previous genome-wide association studies. We found genome-wide significant (p<5×10⁻⁸) evidence for variants at 8q24 and 11p15 and strong supportive evidence for a number of previously reported loci. We found little evidence for individual or systematic inflated association findings resulting from using public controls, demonstrating the utility of using public control data in large-scale genetic association studies of common variants. Taken together, these results demonstrate the importance of established common genetic variants for early-onset prostate cancer and the power of including early-onset prostate cancer cases in genetic association studies.     

A comprehensive analysis of shared loci between systemic lupus erythematosus (SLE) and sixteen autoimmune diseases reveals limited genetic overlap

In spite of the well-known clustering of multiple autoimmune disorders in families, analyses of specific shared genes and polymorphisms between systemic lupus erythematosus (SLE) and other autoimmune diseases (ADs) have been limited. Therefore, we comprehensively tested autoimmune variants for association with SLE, aiming to identify pleiotropic genetic associations between these diseases. We compiled a list of 446 non–Major Histocompatibility Complex (MHC) variants identified in genome-wide association studies (GWAS) of populations of European ancestry across 17 ADs. We then tested these variants in our combined Caucasian SLE cohorts of 1,500 cases and 5,706 controls. We tested a subset of these polymorphisms in an independent Caucasian replication cohort of 2,085 SLE cases and 2,854 controls, allowing the computation of a meta-analysis between all cohorts. We have uncovered novel shared SLE loci that passed multiple comparisons adjustment, including the VTCN1 (rs12046117, P = 2.02×10⁻⁰⁶) region. We observed that the loci shared among the most ADs include IL23R, OLIG3/TNFAIP3, and IL2RA. Given the lack of a universal autoimmune risk locus outside of the MHC and variable specificities for different diseases, our data suggests partial pleiotropy among ADs. Hierarchical clustering of ADs suggested that the most genetically related ADs appear to be type 1 diabetes with rheumatoid arthritis and Crohn''s disease with ulcerative colitis. These findings support a relatively distinct genetic susceptibility for SLE. For many of the shared GWAS autoimmune loci, we found no evidence for association with SLE, including IL23R. Also, several established SLE loci are apparently not associated with other ADs, including the ITGAM-ITGAX and TNFSF4 regions. This study represents the most comprehensive evaluation of shared autoimmune loci to date, supports a relatively distinct non–MHC genetic susceptibility for SLE, provides further evidence for previously and newly identified shared genes in SLE, and highlights the value of studies of potentially pleiotropic genes in autoimmune diseases.     

The −1030/−862-linked TNF promoter single-nucleotide polymorphisms are associated with the inability to control HIV-1 viremia

Control of HIV-1 viremia and progression to AIDS has been associated with specific HLA genes. The tumor necrosis factor (TNF) and the non-classical major histocompatibility (MHC) class I chain-related A (MICA) genes are located in the genomic segment between the HLA class I and II genes and variants of both genes have been identified. We thus analyzed TNF promoter and MICA variants in a well-characterized group of HIV-1 infected individuals with different abilities to control HIV-1 viremia. In our cohort, the –1030/–862-linked TNF promoter single-nucleotide polymorphisms (SNPs), but not MICA variants, are significantly associated with lack of control of HIV-1 viremia (P=0.03). This association is independent of those HLA-B35 alleles associated with HIV-1 disease progression with which the –862 TNF SNP has previously been independently associated. Thus, non-randomly associated genes near the TNF locus are likely involved in control of HIV-1 viremia.     

Psoriasis patients are enriched for genetic variants that protect against HIV-1 disease

An important paradigm in evolutionary genetics is that of a delicate balance between genetic variants that favorably boost host control of infection but which may unfavorably increase susceptibility to autoimmune disease. Here, we investigated whether patients with psoriasis, a common immune-mediated disease of the skin, are enriched for genetic variants that limit the ability of HIV-1 virus to replicate after infection. We analyzed the HLA class I and class II alleles of 1,727 Caucasian psoriasis cases and 3,581 controls and found that psoriasis patients are significantly more likely than controls to have gene variants that are protective against HIV-1 disease. This includes several HLA class I alleles associated with HIV-1 control; amino acid residues at HLA-B positions 67, 70, and 97 that mediate HIV-1 peptide binding; and the deletion polymorphism rs67384697 associated with high surface expression of HLA-C. We also found that the compound genotype KIR3DS1 plus HLA-B Bw4-80I, which respectively encode a natural killer cell activating receptor and its putative ligand, significantly increased psoriasis susceptibility. This compound genotype has also been associated with delay of progression to AIDS. Together, our results suggest that genetic variants that contribute to anti-viral immunity may predispose to the development of psoriasis.     

HLA Diversity in the 1000 Genomes Dataset

The 1000 Genomes Project aims to provide a deep characterization of human genome sequence variation by sequencing at a level that should allow the genome-wide detection of most variants with frequencies as low as 1%. However, in the major histocompatibility complex (MHC), only the top 10 most frequent haplotypes are in the 1% frequency range whereas thousands of haplotypes are present at lower frequencies. Given the limitation of both the coverage and the read length of the sequences generated by the 1000 Genomes Project, the highly variable positions that define HLA alleles may be difficult to identify. We used classical Sanger sequencing techniques to type the HLA-A, HLA-B, HLA-C, HLA-DRB1 and HLA-DQB1 genes in the available 1000 Genomes samples and combined the results with the 103,310 variants in the MHC region genotyped by the 1000 Genomes Project. Using pairwise identity-by-descent distances between individuals and principal component analysis, we established the relationship between ancestry and genetic diversity in the MHC region. As expected, both the MHC variants and the HLA phenotype can identify the major ancestry lineage, informed mainly by the most frequent HLA haplotypes. To some extent, regions of the genome with similar genetic or similar recombination rate have similar properties. An MHC-centric analysis underlines departures between the ancestral background of the MHC and the genome-wide picture. Our analysis of linkage disequilibrium (LD) decay in these samples suggests that overestimation of pairwise LD occurs due to a limited sampling of the MHC diversity. This collection of HLA-specific MHC variants, available on the dbMHC portal, is a valuable resource for future analyses of the role of MHC in population and disease studies.     

Genome-wide association study of relative telomere length

Telomere function is essential to maintaining the physical integrity of linear chromosomes and healthy human aging. The probability of forming proper telomere structures depends on the length of the telomeric DNA tract. We attempted to identify common genetic variants associated with log relative telomere length using genome-wide genotyping data on 3,554 individuals from the Nurses'' Health Study and the Prostate, Lung, Colorectal, and Ovarian Cancer Screening Trial that took part in the National Cancer Institute Cancer Genetic Markers of Susceptibility initiative for breast and prostate cancer. After genotyping 64 independent SNPs selected for replication in additional Nurses'' Health Study and Women''s Genome Health Study participants, we did not identify genome-wide significant loci; however, we replicated the inverse association of log relative telomere length with the minor allele variant [C] of rs16847897 at the TERC locus (per allele β = −0.03, P = 0.003) identified by a previous genome-wide association study. We did not find evidence for an association with variants at the OBFC1 locus or other loci reported to be associated with telomere length. With this sample size we had >80% power to detect β estimates as small as ±0.10 for SNPs with minor allele frequencies of ≥0.15 at genome-wide significance. However, power is greatly reduced for β estimates smaller than ±0.10, such as those for variants at the TERC locus. In general, common genetic variants associated with telomere length homeostasis have been difficult to detect. Potential biological and technical issues are discussed.     

Association of variations in TPH1 and HTR2B with gestational weight gain and measures of obesity

Serotonin is involved in appetite regulation and energy homeostasis. Recently, it has been reported that 5-hydroxytryptamine receptor 2B (Htr2b) and tryptophan hydroxylase 1 (Tph1) play major role in β-cell proliferation in mouse during pregnancy. We investigated the genetic association of HTR2B and TPH1 with risk of gestational diabetes mellitus (GDM) and measures of obesity, in 869 Korean GDM women and carefully selected 632 nondiabetic control subjects. Six single-nucleotide polymorphisms (SNPs) in HTR2B and ten SNPs in TPH1 were selected for genotyping according to their tagging status. Genetic variants in HTR2B and TPH1 were not associated with the risk of GDM. In GDM women, SNPs of TPH1 were significantly associated with weight gain during pregnancy. In nondiabetic controls, SNPs of TPH1 were associated with waist circumference and BMI. We also found that a variant of TPH1 (rs623580) was associated with BMI in a genome-wide association study comprised of 8,842 subjects. Although genetic variants in HTR2B and TPH1 were not associated with risk of GDM, we found significant association of these variants with measures of obesity. However, further replication studies in a different population are required to confirm our findings.     

Defining the spectrum of alleles that contribute to blood lipid concentrations in humans

PURPOSE OF REVIEW: Recently, genome-wide genetic screening of common DNA sequence variants has proven a successful approach to identify novel genetic contributors to complex traits. This review summarizes recent genome-wide association studies for lipid phenotypes, and evaluates the next steps needed to obtain a full picture of genotype-phenotype correlation and apply these findings to inform clinical practice. RECENT FINDINGS: So far, genome-wide association studies have defined at least 19 genomic regions that contain common DNA single nucleotide polymorphisms associated with LDL cholesterol, HDL cholesterol and/or triglycerides. Of these, eight represent novel loci in humans, whereas 11 genes have been previously implicated in lipoprotein metabolism. Many of the same loci with common variants have already been shown to lead to monogenic lipid disorders in humans and/or mice, suggesting that a spectrum of common and rare alleles at each validated locus contributes to blood lipid concentrations. SUMMARY: At least 19 loci harbor common variations that contribute to blood lipid concentrations in humans. Larger scale genome-wide association studies should identify additional loci, and sequencing of these loci should pinpoint all relevant alleles. With a full catalog of DNA polymorphisms in hand, a panel of lipid-related variants can be studied to provide clinical risk stratification and targeting of therapeutic interventions.     

A Genome-Wide Investigation of SNPs and CNVs in Schizophrenia

We report a genome-wide assessment of single nucleotide polymorphisms (SNPs) and copy number variants (CNVs) in schizophrenia. We investigated SNPs using 871 patients and 863 controls, following up the top hits in four independent cohorts comprising 1,460 patients and 12,995 controls, all of European origin. We found no genome-wide significant associations, nor could we provide support for any previously reported candidate gene or genome-wide associations. We went on to examine CNVs using a subset of 1,013 cases and 1,084 controls of European ancestry, and a further set of 60 cases and 64 controls of African ancestry. We found that eight cases and zero controls carried deletions greater than 2 Mb, of which two, at 8p22 and 16p13.11-p12.4, are newly reported here. A further evaluation of 1,378 controls identified no deletions greater than 2 Mb, suggesting a high prior probability of disease involvement when such deletions are observed in cases. We also provide further evidence for some smaller, previously reported, schizophrenia-associated CNVs, such as those in NRXN1 and APBA2. We could not provide strong support for the hypothesis that schizophrenia patients have a significantly greater “load” of large (>100 kb), rare CNVs, nor could we find common CNVs that associate with schizophrenia. Finally, we did not provide support for the suggestion that schizophrenia-associated CNVs may preferentially disrupt genes in neurodevelopmental pathways. Collectively, these analyses provide the first integrated study of SNPs and CNVs in schizophrenia and support the emerging view that rare deleterious variants may be more important in schizophrenia predisposition than common polymorphisms. While our analyses do not suggest that implicated CNVs impinge on particular key pathways, we do support the contribution of specific genomic regions in schizophrenia, presumably due to recurrent mutation. On balance, these data suggest that very few schizophrenia patients share identical genomic causation, potentially complicating efforts to personalize treatment regimens.     

Distinct Genetic Loci Control Plasma HIV-RNA and Cellular HIV-DNA Levels in HIV-1 Infection: The ANRS Genome Wide Association 01 Study

Previous studies of the HIV-1 disease have shown that HLA and Chemokine receptor genetic variants influence disease progression and early viral load. We performed a Genome Wide Association study in a cohort of 605 HIV-1-infected seroconverters for detection of novel genetic factors that influence plasma HIV-RNA and cellular HIV-DNA levels. Most of the SNPs strongly associated with HIV-RNA levels were localised in the 6p21 major histocompatibility complex (MHC) region and were in the vicinity of class I and III genes. Moreover, protective alleles for four disease-associated SNPs in the MHC locus (rs2395029, rs13199524, rs12198173 and rs3093662) were strikingly over-represented among forty-five Long Term HIV controllers. Furthermore, we show that the HIV-DNA levels (reflecting the HIV reservoir) are associated with the same four SNPs, but also with two additional SNPs on chromosome 17 (rs6503919; intergenic region flanked by the DDX40 and YPEL2 genes) and chromosome 8 (rs2575735; within the Syndecan 2 gene). Our data provide evidence that the MHC controls both HIV replication and HIV reservoir. They also indicate that two additional genomic loci may influence the HIV reservoir.     

Identification of a Potential Regulatory Variant for Colorectal Cancer Risk Mapping to Chromosome 5q31.1: A Post-GWAS Study

Large-scale genome-wide association studies (GWAS) have established chromosome 5q31.1 as a susceptibility locus for colorectal cancer (CRC), which was still lack of causal genetic variants. We searched potentially regulatory single nucleotide polymorphisms (SNPs) in the overlap region between linkage disequilibrium (LD) block of 5q31.1 and regulatory elements predicted by histone modifications, then tested their association with CRC via a case-control study. Among three candidate common variants, we found rs17716310 conferred significantly (heterozygous model: OR = 1.273, 95% confidence interval (95%CI) = 1.016–1.595, P = 0.036) and marginally (dominant model: OR = 1.238, 95%CI = 1.000–1.532, P = 0.050) increase risk for CRC in a Chinese population including 695 cases and 709 controls. This variation was suggested to be regulatory altering the activity of enhancer that control PITX1 expression. Using epigenetic information such as chromatin immunoprecipitation-sequencing (ChIP-seq) data might help researchers to interpret the results of GWAS and locate causal variants for diseases in post-GWAS era.     

Genome-wide analysis of copy number variation in type 1 diabetes

Type 1 diabetes (T1D) tends to cluster in families, suggesting there may be a genetic component predisposing to disease. However, a recent large-scale genome-wide association study concluded that identified genetic factors, single nucleotide polymorphisms, do not account for overall familiality. Another class of genetic variation is the amplification or deletion of >1 kilobase segments of the genome, also termed copy number variations (CNVs). We performed genome-wide CNV analysis on a cohort of 20 unrelated adults with T1D and a control (Ctrl) cohort of 20 subjects using the Affymetrix SNP Array 6.0 in combination with the Birdsuite copy number calling software. We identified 39 CNVs as enriched or depleted in T1D versus Ctrl. Additionally, we performed CNV analysis in a group of 10 monozygotic twin pairs discordant for T1D. Eleven of these 39 CNVs were also respectively enriched or depleted in the Twin cohort, suggesting that these variants may be involved in the development of islet autoimmunity, as the presently unaffected twin is at high risk for developing islet autoimmunity and T1D in his or her lifetime. These CNVs include a deletion on chromosome 6p21, near an HLA-DQ allele. CNVs were found that were both enriched or depleted in patients with or at high risk for developing T1D. These regions may represent genetic variants contributing to development of islet autoimmunity in T1D.     

Genetic Association Studies: An Information Content Perspective

The availability of high-density single nucleotide polymorphisms (SNPs) data has made the human genetic association studies possible to identify common and rare variants underlying complex diseases in a genome-wide scale. A handful of novel genetic variants have been identified, which gives much hope and prospects for the future of genetic association studies. In this process, statistical and computational methods play key roles, among which information-based association tests have gained large popularity. This paper is intended to give a comprehensive review of the current literature in genetic association analysis casted in the framework of information theory. We focus our review on the following topics: (1) information theoretic approaches in genetic linkage and association studies; (2) entropy-based strategies for optimal SNP subset selection; and (3) the usage of theoretic information criteria in gene clustering and gene regulatory network construction.     

Prediction of individual genetic risk of complex disease

Most common diseases are caused by multiple genetic and environmental factors. In the last 2 years, genome-wide association studies (GWAS) have identified polymorphisms that are associated with risk to common disease, but the effect of any one risk allele is typically small. By combining information from many risk variants, will it be possible to predict accurately each individual person's genetic risk for a disease? In this review we consider the lessons from GWAS and the implications for genetic risk prediction to common disease. We conclude that with larger GWAS sample sizes or by combining studies, accurate prediction of genetic risk will be possible, even if the causal mutations or the mechanisms by which they affect susceptibility are unknown.