The influence of CCL3L1 gene-containing segmental duplications on HIV-1/AIDS susceptibility

Copy number variations of TBK1 in Australian patients with primary open-angle glaucoma

PURPOSE

To investigate the presence of TBK1 copy number variations in a large, well-characterized Australian cohort of patients with glaucoma comprising both normal-tension glaucoma and high-tension glaucoma cases.

DESIGN

A retrospective cohort study.

METHODS

DNA samples from patients with normal-tension glaucoma and high-tension glaucoma and unaffected controls were screened for TBK1 copy number variations using real-time quantitative polymerase chain reaction. Samples with additional copies of the TBK1 gene were further tested using custom comparative genomic hybridization arrays.

RESULTS

Four out of 334 normal-tension glaucoma cases (1.2%) were found to carry TBK1 copy number variations using quantitative polymerase chain reaction. One extra dose of the TBK1 gene (duplication) was detected in 3 normal-tension glaucoma patients, while 2 extra doses of the gene (triplication) were detected in a fourth normal-tension glaucoma patient. The results were further confirmed by custom comparative genomic hybridization arrays. Further, the TBK1 copy number variation segregated with normal-tension glaucoma in the family members of the probands, showing an autosomal dominant pattern of inheritance. No TBK1 copy number variations were detected in 1045 Australian patients with high-tension glaucoma or in 254 unaffected controls.

CONCLUSION

We report the presence of TBK1 copy number variations in our Australian normal-tension glaucoma cohort, including the first example of more than 1 extra copy of this gene in glaucoma patients (gene triplication). These results confirm TBK1 to be an important cause of normal-tension glaucoma, but do not suggest common involvement in high-tension glaucoma.

Impact of copy number variations burden on coding genome in humans using integrated high resolution arrays

Copy number variations (CNVs) alter the transcriptional and translational levels of genes by disrupting the coding structure and this burden of CNVs seems to be a significant contributor to phenotypic variations. Therefore it was necessary to assess the complexities of CNV burden on the coding genome. A total of 1715 individuals from 12 populations were used for CNV analysis in the present investigation. Analysis was performed using Affymetrix Genome-Wide Human SNP Array 6·0 chip and CytoScan High-Density arrays. CNVs were more frequently observed in the coding region than in the non-coding region. CNVs were observed vastly more frequently in the coding region than the non-coding region. CNVs were found to be enriched in the regions containing functional genes (83-96%) compared with the regions containing pseudogenes (4-17%). CNVs across the genome of an individual showed multiple hits across many genes, whose proteins interact physically and function under the same pathway. We identified varying numbers of proteins and degrees of interactions within protein complexes of single individual genomes. This study represents the first draft of a population-specific CNV genes map as well as a cross-populational map. The complex relationship of CNVs on genes and their physically interacting partners unravels many complexities involved in phenotype expression. This study identifies four mechanisms contributing to the complexities caused by the presence of multiple CNVs across many genes in the coding part of the genome.

Biological and Structural Characterization of Rotamers of C-C Chemokine Receptor Type 5 (CCR5) Inhibitor GSK214096

We recently reported the discovery of preclinical CCR5 inhibitor GSK214096, 1 (J. Med. Chem. 2011, 54, 756). Detailed characterization of 1 revealed that it exists as a mixture of four separable atropisomers A-D. The two slow-interconverting pairs of rotamers A + B and C + D were separated and further characterized. HIV and CCR5-mediated chemotaxis data strongly suggest that the antiviral potency of 1 is due to rotamers A + B and not C + D. Furthermore, integrated UV, vibrational circular dichroism VCD and computational approach allowed to determine the M chirality in C + D (and P chirality in A + B). These findings imply additional avenues to be pursued toward new CCR5 antagonists.

Copy number variations and human genetic disease

PURPOSE OF REVIEW

Recent studies clearly demonstrate that copy number variations (CNVs) are widespread in our genome and play an important role in human genetic variation, accounting for both human population diversity and human genetic disease. This review will discuss the most current knowledge regarding our understanding of the biology of CNVs in relation to human genetic disease.

RECENT FINDINGS

CNVs associated with human genetic disease can be either recurrent, with a common size and breakpoint clustering, or nonrecurrent, with different sizes and variable breakpoints. Two types of recurrent CNVs have been distinguished, including the syndromic forms in which the phenotypic features are relatively consistent, and those in which the same recurrent CNV can be associated with a diverse set of diagnoses. Recently, the 'Two-hit model' was used to explain the phenotypic variability associated with the latter group of recurrent CNVs. Nonrecurrent CNVs, on the contrary, occur at a relatively lower frequency at the individual locus level but collectively they are as common as recurrent CNVs. Finally, the study of CNV burden in different diseases demonstrated a clear trend of an increasing CNV burden in diseases with more severe phenotypes.

SUMMARY

In spite of the advances in the study of the CNV landscape associated with human genetic disease, there still remain many unexplored questions especially regarding the role of CNVs in the pathogenesis of complex human genetic diseases.

New genes contribute to genetic and phenotypic novelties in human evolution

New genes in human genomes have been found relevant in evolution and biology of humans. It was conservatively estimated that the human genome encodes more than 300 human-specific genes and 1000 primate-specific genes. These new arrivals appear to be implicated in brain function and male reproduction. Surprisingly, increasing evidence indicates that they may also bring negative pleiotropic effects, while assuming various possible biological functions as sources of phenotypic novelties, suggesting a non-progressive route for functional evolution. Similar to these fixed new genes, polymorphic new genes were found to contribute to functional evolution within species, for example, with respect to digestion or disease resistance, revealing that new genes can acquire new or diverged functions in its initial stage as prototypic genes. These progresses have provided new opportunities to explore the genetic basis of human biology and human evolutionary history in a new dimension.

Frequency distribution of autoimmunity associated FCGR3B gene copy number in Indian population

Amongst several human genome variations, copy number variations (CNVs) are considered as an important source of variability contributing to susceptibility to wide range of diseases. Although CNV is scattered for genes throughout the human genome, several of autoimmunity related genes have CN variation and therefore play an important role in susceptibility to autoimmune diseases. The association of the Fc gamma receptor 3B (FCGR3B) gene copy number in autoimmunity is well characterized in various populations studied. The Fc gamma receptor is a low affinity, glycosylphosphatidylinositol-linked receptor for IgG molecule predominantly expressed on human neutrophils. The variable gene copy number of FCGR3B is found to be involved in the impaired clearance of immune complexes, which significantly contribute to the pathogenesis of several autoimmune diseases such as systemic lupus erythematosus (SLE), rheumatoid arthritis (RA), type-1 diabetes and others. The FCGR3B copy number ranged from 0 to ≥ 2 copies per diploid genome in other populations, but yet not explored in Indian population. Hence, this study aims to evaluate the variation in the frequency distribution of FCGR3B CNV in Indian population. FCGR3B gene copy number varied significantly when compared to other population of the world. This observation will help us in exploring the potential role of CNV in FCGR3B gene and its association to autoimmune disorders in Indian population.

Cellular differentiation regulator BLIMP1 induces Epstein-Barr virus lytic reactivation in epithelial and B cells by activating transcription from both the R and Z promoters

Epstein-Barr virus (EBV) maintains a lifelong latent infection within a subset of its host's memory B cells, while lytic EBV replication takes place in plasma cells and differentiated epithelial cells. Therefore, cellular transcription factors, such as BLIMP1, that are key mediators of differentiation likely contribute to the EBV latent-to-lytic switch. Previous reports showed that ectopic BLIMP1 expression induces reactivation in some EBV-positive (EBV(+)) B-cell lines and transcription from Zp, with all Z(+) cells in oral hairy leukoplakia being BLIMP1(+). Here, we examined BLIMP1's role in inducing EBV lytic gene expression in numerous EBV(+) epithelial and B-cell lines and activating transcription from Rp. BLIMP1 addition was sufficient to induce reactivation in latently infected epithelial cells derived from gastric cancers, nasopharyngeal carcinomas, and normal oral keratinocytes (NOK) as well as some, but not all B-cell lines. BLIMP1 strongly induced transcription from Rp as well as Zp, with there being three or more synergistically acting BLIMP1-responsive elements (BRE) within Rp. BLIMP1's DNA-binding domain was required for reactivation, but BLIMP1 did not directly bind the nucleotide (nt) -660 Rp BRE. siRNA knockdown of BLIMP1 inhibited 12-O-tetradecanoyl-phorbol-13-acetate (TPA)-induced lytic reactivation in NOK-Akata cells, cells that can be reactivated by R, but not Z. Thus, we conclude that BLIMP1 expression is both necessary and sufficient to induce EBV lytic replication in many (possibly all) EBV(+) epithelial-cell types, but in only a subset of EBV(+) B-cell types; it does so, at least in part, by strongly activating expression of both EBV immediately early genes, BZLF1 and BRLF1.

IMPORTANCE

This study is the first one to show that the cellular transcription factor BLIMP1, a key player in both epithelial and B-cell differentiation, induces reactivation of the oncogenic herpesvirus Epstein-Barr virus (EBV) out of latency into lytic replication in a variety of cancerous epithelial cell types as well as in some, but not all, B-cell types that contain this virus in a dormant state. The mechanism by which BLIMP1 does so involves strongly turning on expression of both of the immediate early genes of the virus, probably by directly acting upon the promoters as part of protein complexes or indirectly by altering the expression or activities of some cellular transcription factors and signaling pathways. The fact that EBV(+) cancers usually contain mostly undifferentiated cells may be due in part to these cells dying from lytic EBV infection when they differentiate and express wild-type BLIMP1.

Host genetic factors associated with symptomatic primary HIV infection and disease progression among Argentinean seroconverters

BACKGROUND

Variants in HIV-coreceptor C-C chemokine receptor type 5 (CCR5) and Human leukocyte antigen (HLA) genes are the most important host genetic factors associated with HIV infection and disease progression. Our aim was to analyze the association of these genetic factors in the presence of clinical symptoms during Primary HIV Infection (PHI) and disease progression within the first year.

METHODS

Seventy subjects diagnosed during PHI were studied (55 symptomatic and 15 asymptomatic). Viral load (VL) and CD4 T-cell count were evaluated. HIV progression was defined by presence of B or C events and/or CD4 T-cell counts <350 cell/mm3. CCR5 haplotypes were characterized by polymerase chain reaction and SDM-PCR-RFLP. HLA-I characterization was performed by Sequencing.

RESULTS

Symptoms during PHI were significantly associated with lower frequency of CCR5-CF1 (1.8% vs. 26.7%, p = 0.006). Rapid progression was significantly associated with higher frequency of CCR5-CF2 (16.7% vs. 0%, p = 0.024) and HLA-A*11 (16.7% vs. 1.2%, p = 0.003) and lower frequency of HLA-C*3 (2.8% vs. 17.5%, p = 0.035). Higher baseline VL was significantly associated with presence of HLA-A*11, HLA-A*24, and absence of HLA-A*31 and HLA-B*57. Higher 6-month VL was significantly associated with presence of CCR5-HHE, HLA-A*24, HLA-B*53, and absence of HLA-A*31 and CCR5-CF1. Lower baseline CD4 T-cell count was significantly associated with presence of HLA-A*24/*33, HLA-B*53, CCR5-CF2 and absence of HLA-A*01/*23 and CCR5-HHA. Lower 6-month CD4 T-cell count was associated with presence of HLA-A*24 and HLA-B*53, and absence of HLA-A*01 and HLA-B*07/*39. Moreover, lower 12-month CD4 T-cell count was significantly associated with presence of HLA-A*33, HLA-B*14, HLA-C*08, CCR5-CF2, and absence of HLA-B*07 and HLA-C*07.

CONCLUSION

Several host factors were significantly associated with disease progression in PHI subjects. Most results agree with previous studies performed in other groups. However, some genetic factor associations are being described for the first time, highlighting the importance of genetic studies at a local level.

A genome-wide detection of copy number variation using SNP genotyping arrays in Beijing-You chickens

Copy number variation (CNV) has been recently examined in many species and is recognized as being a source of genetic variability, especially for disease-related phenotypes. In this study, the PennCNV software, a genome-wide CNV detection system based on the 60 K SNP BeadChip was used on a total sample size of 1,310 Beijing-You chickens (a Chinese local breed). After quality control, 137 high confidence CNVRs covering 27.31 Mb of the chicken genome and corresponding to 2.61 % of the whole chicken genome. Within these regions, 131 known genes or coding sequences were involved. Q-PCR was applied to verify some of the genes related to disease development. Results showed that copy number of genes such as, phosphatidylinositol-5-phosphate 4-kinase II alpha, PHD finger protein 14, RHACD8 (a CD8α- like messenger RNA), MHC B-G, zinc finger protein, sarcosine dehydrogenase and ficolin 2 varied between individual chickens, which also supports the reliability of chip-detection of the CNVs. As one source of genomic variation, CNVs may provide new insight into the relationship between the genome and phenotypic characteristics.

Copy number deletion has little impact on gene expression levels in racehorses

Copy number variations (CNVs), important genetic factors for study of human diseases, may have as large of an effect on phenotype as do single nucleotide polymorphisms. Indeed, it is widely accepted that CNVs are associated with differential disease susceptibility. However, the relationships between CNVs and gene expression have not been characterized in the horse. In this study, we investigated the effects of copy number deletion in the blood and muscle transcriptomes of Thoroughbred racing horses. We identified a total of 1,246 CNVs of deletion polymorphisms using DNA re-sequencing data from 18 Thoroughbred racing horses. To discover the tendencies between CNV status and gene expression levels, we extracted CNVs of four Thoroughbred racing horses of which RNA sequencing was available. We found that 252 pairs of CNVs and genes were associated in the four horse samples. We did not observe a clear and consistent relationship between the deletion status of CNVs and gene expression levels before and after exercise in blood and muscle. However, we found some pairs of CNVs and associated genes that indicated relationships with gene expression levels: a positive relationship with genes responsible for membrane structure or cytoskeleton and a negative relationship with genes involved in disease. This study will lead to conceptual advances in understanding the relationship between CNVs and global gene expression in the horse.

Genetics of Sub-Saharan African Human Population: Implications for HIV/AIDS, Tuberculosis, and Malaria

Sub-Saharan Africa has continued leading in prevalence and incidence of major infectious disease killers such as HIV/AIDS, tuberculosis, and malaria. Epidemiological triad of infectious diseases includes susceptible host, pathogen, and environment. It is imperative that all aspects of vertices of the infectious disease triad are analysed to better understand why this is so. Studies done to address this intriguing reality though have mainly addressed pathogen and environmental components of the triad. Africa is the most genetically diverse region of the world as well as being the origin of modern humans. Malaria is relatively an ancient infection in this region as compared to TB and HIV/AIDS; from the evolutionary perspective, we would draw lessons that this ancestrally unique population now under three important infectious diseases both ancient and exotic will be skewed into increased genetic diversity; moreover, other evolutionary forces are also still at play. Host genetic diversity resulting from many years of malaria infection has been well documented in this population; we are yet to account for genetic diversity from the trio of these infections. Effect of host genetics on treatment outcome has been documented. Host genetics of sub-Saharan African population and its implication to infectious diseases are an important aspect that this review seeks to address.

Whole-genome sequencing of individuals from a founder population identifies candidate genes for asthma

Asthma is a complex genetic disease caused by a combination of genetic and environmental risk factors. We sought to test classes of genetic variants largely missed by genome-wide association studies (GWAS), including copy number variants (CNVs) and low-frequency variants, by performing whole-genome sequencing (WGS) on 16 individuals from asthma-enriched and asthma-depleted families. The samples were obtained from an extended 13-generation Hutterite pedigree with reduced genetic heterogeneity due to a small founding gene pool and reduced environmental heterogeneity as a result of a communal lifestyle. We sequenced each individual to an average depth of 13-fold, generated a comprehensive catalog of genetic variants, and tested the most severe mutations for association with asthma. We identified and validated 1960 CNVs, 19 nonsense or splice-site single nucleotide variants (SNVs), and 18 insertions or deletions that were out of frame. As follow-up, we performed targeted sequencing of 16 genes in 837 cases and 540 controls of Puerto Rican ancestry and found that controls carry a significantly higher burden of mutations in IL27RA (2.0% of controls; 0.23% of cases; nominal p = 0.004; Bonferroni p = 0.21). We also genotyped 593 CNVs in 1199 Hutterite individuals. We identified a nominally significant association (p = 0.03; Odds ratio (OR) = 3.13) between a 6 kbp deletion in an intron of NEDD4L and increased risk of asthma. We genotyped this deletion in an additional 4787 non-Hutterite individuals (nominal p = 0.056; OR = 1.69). NEDD4L is expressed in bronchial epithelial cells, and conditional knockout of this gene in the lung in mice leads to severe inflammation and mucus accumulation. Our study represents one of the early instances of applying WGS to complex disease with a large environmental component and demonstrates how WGS can identify risk variants, including CNVs and low-frequency variants, largely untested in GWAS.

Combined analysis with copy number variation identifies risk loci in lung cancer

Background. Lung cancer is the most important cause of cancer mortality worldwide, but the underlying mechanisms of this disease are not fully understood. Copy number variations (CNVs) are promising genetic variations to study because of their potential effects on cancer. Methodology/Principal Findings. Here we conducted a pilot study in which we systematically analyzed the association of CNVs in two lung cancer datasets: the Environment And Genetics in Lung cancer Etiology (EAGLE) and the Prostate, Lung, Colorectal and Ovarian (PLCO) Cancer Screening Trial datasets. We used a preestablished association method to test the datasets separately and conducted a combined analysis to test the association accordance between the two datasets. Finally, we identified 167 risk SNP loci and 22 CNVs associated with lung cancer and linked them with recombination hotspots. Functional annotation and biological relevance analyses implied that some of our predicted risk loci were supported by other studies and might be potential candidate loci for lung cancer studies. Conclusions/Significance. Our results further emphasized the importance of copy number variations in cancer and might be a valuable complement to current genome-wide association studies on cancer.

The CNVrd2 package: measurement of copy number at complex loci using high-throughput sequencing data

Recent advances in high-throughout sequencing technologies have made it possible to accurately assign copy number (CN) at CN variable loci. However, current analytic methods often perform poorly in regions in which complex CN variation is observed. Here we report the development of a read depth-based approach, CNVrd2, for investigation of CN variation using high-throughput sequencing data. This methodology was developed using data from the 1000 Genomes Project from the CCL3L1 locus, and tested using data from the DEFB103A locus. In both cases, samples were selected for which paralog ratio test data were also available for comparison. The CNVrd2 method first uses observed read-count ratios to refine segmentation results in one population. Then a linear regression model is applied to adjust the results across multiple populations, in combination with a Bayesian normal mixture model to cluster segmentation scores into groups for individual CN counts. The performance of CNVrd2 was compared to that of two other read depth-based methods (CNVnator, cn.mops) at the CCL3L1 and DEFB103A loci. The highest concordance with the paralog ratio test method was observed for CNVrd2 (77.8/90.4% for CNVrd2, 36.7/4.8% for cn.mops and 7.2/1% for CNVnator at CCL3L1 and DEF103A). CNVrd2 is available as an R package as part of the Bioconductor project: http://www.bioconductor.org/packages/release/bioc/html/CNVrd2.html.

Inferring mechanisms of copy number change from haplotype structures at the human DEFA1A3 locus

BACKGROUND

The determination of structural haplotypes at copy number variable regions can indicate the mechanisms responsible for changes in copy number, as well as explain the relationship between gene copy number and expression. However, obtaining spatial information at regions displaying extensive copy number variation, such as the DEFA1A3 locus, is complex, because of the difficulty in the phasing and assembly of these regions. The DEFA1A3 locus is intriguing in that it falls within a region of high linkage disequilibrium, despite its high variability in copy number (n = 3-16); hence, the mechanisms responsible for changes in copy number at this locus are unclear.

RESULTS

In this study, a region flanking the DEFA1A3 locus was sequenced across 120 independent haplotypes with European ancestry, identifying five common classes of DEFA1A3 haplotype. Assigning DEFA1A3 class to haplotypes within the 1000 Genomes project highlights a significant difference in DEFA1A3 class frequencies between populations with different ancestry. The features of each DEFA1A3 class, for example, the associated DEFA1A3 copy numbers, were initially assessed in a European cohort (n = 599) and replicated in the 1000 Genomes samples, showing within-class similarity, but between-class and between-population differences in the features of the DEFA1A3 locus. Emulsion haplotype fusion-PCR was used to generate 61 structural haplotypes at the DEFA1A3 locus, showing a high within-class similarity in structure.

CONCLUSIONS

Structural haplotypes across the DEFA1A3 locus indicate that intra-allelic rearrangement is the predominant mechanism responsible for changes in DEFA1A3 copy number, explaining the conservation of linkage disequilibrium across the locus. The identification of common structural haplotypes at the DEFA1A3 locus could aid studies into how DEFA1A3 copy number influences expression, which is currently unclear.

Novel population specific autosomal copy number variation and its functional analysis amongst Negritos from Peninsular Malaysia

Copy number variation (CNV) has been recognized as a major contributor to human genome diversity. It plays an important role in determining phenotypes and has been associated with a number of common and complex diseases. However CNV data from diverse populations is still limited. Here we report the first investigation of CNV in the indigenous populations from Peninsular Malaysia. We genotyped 34 Negrito genomes from Peninsular Malaysia using the Affymetrix SNP 6.0 microarray and identified 48 putative novel CNVs, consisting of 24 gains and 24 losses, of which 5 were identified in at least 2 unrelated samples. These CNVs appear unique to the Negrito population and were absent in the DGV, HapMap3 and Singapore Genome Variation Project (SGVP) datasets. Analysis of gene ontology revealed that genes within these CNVs were enriched in the immune system (GO:0002376), response to stimulus mechanisms (GO:0050896), the metabolic pathways (GO:0001852), as well as regulation of transcription (GO:0006355). Copy number gains in CNV regions (CNVRs) enriched with genes were significantly higher than the losses (P value <0.001). In view of the small population size, relative isolation and semi-nomadic lifestyles of this community, we speculate that these CNVs may be attributed to recent local adaptation of Negritos from Peninsular Malaysia.

Digital genotyping of macrosatellites and multicopy genes reveals novel biological functions associated with copy number variation of large tandem repeats

Tandem repeats are common in eukaryotic genomes, but due to difficulties in assaying them remain poorly studied. Here, we demonstrate the utility of Nanostring technology as a targeted approach to perform accurate measurement of tandem repeats even at extremely high copy number, and apply this technology to genotype 165 HapMap samples from three different populations and five species of non-human primates. We observed extreme variability in copy number of tandemly repeated genes, with many loci showing 5–10 fold variation in copy number among humans. Many of these loci show hallmarks of genome assembly errors, and the true copy number of many large tandem repeats is significantly under-represented even in the high quality ‘finished’ human reference assembly. Importantly, we demonstrate that most large tandem repeat variations are not tagged by nearby SNPs, and are therefore essentially invisible to SNP-based GWAS approaches. Using association analysis we identify many cis correlations of large tandem repeat variants with nearby gene expression and DNA methylation levels, indicating that variations of tandem repeat length are associated with functional effects on the local genomic environment. This includes an example where expansion of a macrosatellite repeat is associated with increased DNA methylation and suppression of nearby gene expression, suggesting a mechanism termed “repeat induced gene silencing”, which has previously been observed only in transgenic organisms. We also observed multiple signatures consistent with altered selective pressures at tandemly repeated loci, suggesting important biological functions. Our studies show that tandemly repeated loci represent a highly variable fraction of the genome that have been systematically ignored by most previous studies, copy number variation of which can exert functionally significant effects. We suggest that future studies of tandem repeat loci will lead to many novel insights into their role in modulating both genomic and phenotypic diversity.

Here we utilize Nanostring digital assays and show their utility for estimating copy number of 186 multicopy genes and tandem repeats. By analyzing patterns of single nucleotide variation around these variants, we show that copy number variation at the vast majority of tandem repeat variations is not effectively tagged by nearby SNPs, and thus standard genome-wide association studies that focus on SNPs provide little or no information about such variants. By comparing patterns of tandem repeat copy number with variation in local gene expression and DNA methylation, we also identify extensive functional effects on local genome function. This includes an example of a non-coding macrosatellite repeat, expansion of which exerts a repressive effect on a nearby gene accompanied by accumulations of local DNA methylation. Finally, comparison of diverse human populations with a number of primate genomes shows that many of these sequences have undergone extreme changes in copy number during recent human and primate evolution, and show signatures that suggest possible selective effects. Overall, we conclude that multicopy genes and macrosatellites represent a highly variable fraction of the genome with important functional effects that has been systematically ignored by previous studies.

Exploring the function of genetic variants in the non-coding genomic regions: approaches for identifying human regulatory variants affecting gene expression

Understanding the genetic basis of human traits/diseases and the underlying mechanisms of how these traits/diseases are affected by genetic variations is critical for public health. Current genome-wide functional genomics data uncovered a large number of functional elements in the noncoding regions of human genome, providing new opportunities to study regulatory variants (RVs). RVs play important roles in transcription factor bindings, chromatin states and epigenetic modifications. Here, we systematically review an array of methods currently used to map RVs as well as the computational approaches in annotating and interpreting their regulatory effects, with emphasis on regulatory single-nucleotide polymorphism. We also briefly introduce experimental methods to validate these functional RVs.

Human gene copy number variation and infectious disease

Variability in the susceptibility to infectious disease and its clinical manifestation can be determined by variation in the environment and by genetic variation in the pathogen and the host. Despite several successes based on candidate gene studies, defining the host variation affecting infectious disease has not been as successful as for other multifactorial diseases. Both single nucleotide variation and copy number variation (CNV) of the host contribute to the host's susceptibility to infectious disease. In this review we focus on CNV, particularly on complex multiallelic CNV that is often not well characterised either directly by hybridisation methods or indirectly by analysis of genotypes and flanking single nucleotide variants. We summarise the well-known examples, such as α-globin deletion and susceptibility to severe malaria, as well as more recent controversies, such as the extensive CNV of the chemokine gene CCL3L1 and HIV infection. We discuss the potential biological mechanisms that could underly any genetic association and reflect on the extensive complexity and functional variation generated by a combination of CNV and sequence variation, as illustrated by the Fc gamma receptor genes FCGR3A, FCGR3B and FCGR2C. We also highlight some understudied areas that might prove fruitful areas for further research.

Transcriptional and posttranscriptional regulation of cytokine gene expression in HIV-1 antigen-specific CD8+ T cells that mediate virus inhibition

The ability of CD8+ T cells to effectively limit HIV-1 replication and block HIV-1 acquisition is determined by the capacity to rapidly respond to HIV-1 antigens. Understanding both the functional properties and regulation of an effective CD8+ response would enable better evaluation of T cell-directed vaccine strategies and may inform the design of new therapies. We assessed the antigen specificity, cytokine signature, and mechanisms that regulate antiviral gene expression in CD8+ T cells from a cohort of HIV-1-infected virus controllers (VCs) (<5,000 HIV-1 RNA copies/ml and CD4+ lymphocyte counts of >400 cells/μl) capable of soluble inhibition of HIV-1. Gag p24 and Nef CD8+ T cell-specific soluble virus inhibition was common among the VCs and correlated with substantial increases in the abundance of mRNAs encoding the antiviral cytokines macrophage inflammatory proteins MIP-1α, MIP-1αP (CCL3L1), and MIP-1β; granulocyte-macrophage colony-stimulating factor (GM-CSF); lymphotactin (XCL1); tumor necrosis factor receptor superfamily member 9 (TNFRSF9); and gamma interferon (IFN-γ). The induction of several of these mRNAs was driven through a coordinated response of both increased transcription and stabilization of mRNA, which together accounted for the observed increase in mRNA abundance. This coordinated response allows rapid and robust induction of mRNA messages that can enhance the CD8+ T cells' ability to inhibit virus upon antigen encounter.

IMPORTANCE

We show that mRNA stability, in addition to transcription, is key in regulating the direct anti-HIV-1 function of antigen-specific memory CD8+ T cells. Regulation at the level of RNA helps enable rapid recall of memory CD8+ T cell effector functions for HIV-1 inhibition. By uncovering and understanding the mechanisms employed by CD8+ T cell subsets with antigen-specific anti-HIV-1 activity, we can identify new strategies for comprehensive identification of other important antiviral genes. This will, in turn, enhance our ability to inhibit virus replication by informing both cure strategies and HIV-1 vaccine designs that aim to reduce transmission and can aid in blocking HIV-1 acquisition.

Copy number variants and selective sweeps in natural populations of the house mouse (Mus musculus domesticus)

Copy-number variants (CNVs) may play an important role in early adaptations, potentially facilitating rapid divergence of populations. We describe an approach to study this question by investigating CNVs present in natural populations of mice in the early stages of divergence and their involvement in selective sweeps. We have analyzed individuals from two recently diverged natural populations of the house mouse (Mus musculus domesticus) from Germany and France using custom, high-density, comparative genome hybridization arrays (CGH) that covered almost 164 Mb and 2444 genes. One thousand eight hundred and sixty one of those genes we previously identified as differentially expressed between these populations, while the expression of the remaining genes was invariant. In total, we identified 1868 CNVs across all 10 samples, 200 bp to 600 kb in size and affecting 424 genic regions. Roughly two thirds of all CNVs found were deletions. We found no enrichment of CNVs among the differentially expressed genes between the populations compared to the invariant ones, nor any meaningful correlation between CNVs and gene expression changes. Among the CNV genes, we found cellular component gene ontology categories of the synapse overrepresented among all the 2444 genes tested. To investigate potential adaptive significance of the CNV regions, we selected six that showed large differences in frequency of CNVs between the two populations and analyzed variation in at least two microsatellites surrounding the loci in a sample of 46 unrelated animals from the same populations collected in field trappings. We identified two loci with large differences in microsatellite heterozygosity (Sfi1 and Glo1/Dnahc8 regions) and one locus with low variation across the populations (Cmah), thus suggesting that these genomic regions might have recently undergone selective sweeps. Interestingly, the Glo1 CNV has previously been implicated in anxiety-like behavior in mice, suggesting a differential evolution of a behavioral trait.

Detection of copy number variation by SNP-allelotyping

Charcot-Marie-Tooth disease type 1A (CMT1A) is caused by an abnormal copy number variation (CNV) with a trisomy of chromosome 17p12. The increase of the DNA-segment copy number is expected to alter the allele frequency of single nucleotide polymorphism (SNP) within the duplicated region. We tested whether SNP allele frequency determined by a Sequenom MassArray can be used to detect the CMT1A mutation. Our results revealed distinct patterns of SNP allele frequency distribution, which reliably differentiated CMT1A patients from controls. This finding suggests that this technique may serve as an alternative approach to identifying CNV in certain diseases, including CMT1A.

An MLPA-based approach for high-resolution genotyping of disease-related multi-allelic CNVs

Copy number variation has recently been recognized as an important type of genetic variation that modifies human phenotypes. Copy number variants (CNVs) are being increasingly associated with various human phenotypes and diseases. However, the lack of an appropriate method that allows fast, inexpensive and, most importantly, accurate CNVs genotyping significantly hampers CNV analysis. This limitation especially affects the analysis of multi-allelic CNVs that frequently modify various phenotypes. Recently, we developed a multiplex ligation-dependent probe amplification (MLPA)-based strategy for multiplex copy number genotyping and the validation of candidate CNV-miRNAs. Here we present the adaptation and optimization of this recently developed method for high-resolution genotyping of individual disease-related multi-allelic CNVs. We developed appropriate assays for three well-known and extensively studied CNVs: CNV-CCL3L1, CNV-DEFB, and CNV-UGT2B17, which have been associated with various human phenotypes including inflammation-related and infectious diseases. With the use of these assays we identified several general factors that allow to increase the resolution of the copy number genotyping. Performed experiments confirmed the high reproducibility and accuracy of the obtained genotyping results. The reliability of the results and relatively low per-genotype cost makes this strategy an attractive method for large-scale experiments such as genotype-phenotype association studies.

African ancestry influences CCR5 -2459G>A genotype-associated virologic success of highly active antiretroviral therapy

INTRODUCTION

In a North American, HIV-positive, highly active antiretroviral therapy (HAART)-treated, adherent cohort of self-identified white and black patients, we previously observed that chemokine (C-C motif) receptor 5 (CCR5) -2459G>A genotype had a strong association with time to achieve virologic success (TVLS) in black but not in white patients.

METHODS

Using 128 genome-wide ancestry informative markers, we performed a quantitative assessment of ancestry in these patients (n = 310) to determine (1) whether CCR5 -2459G>A genotype is still associated with TVLS of HAART when ancestry, not self-identified race, is considered and (2) whether this association is influenced by varying African ancestry.

RESULTS

We found that the interaction between CCR5 -2459G>A genotype and African ancestry (≤ 0.125 vs. ≥ 0.425 and <0.71 vs. ≥ 0.71) was significantly associated with TVLS (GG compared with AA, P = 0.044 and 0.018, respectively). Furthermore, the association between CCR5 -2459G>A genotype and TVLS was stronger in patients with African ancestry ≥ 0.71 than in patients with African ancestry ≥ 0.452, in both Kaplan-Meier (log-rank P = 0.039 and 0.057, respectively, for AA, GA, and GG) and Cox proportional hazards regression (relative hazard for GG compared with AA 2.59 [95% confidence interval: 1.27 to 5.22; P = 0.01] and 2.26 [95% confidence interval: 1.18 to 4.32; P = 0.01], respectively) analyses.

CONCLUSIONS

We observed that the association between CCR5 -2459G>A genotype and TVLS of HAART increased with stronger African ancestry. Understanding the genomic mechanisms by which African ancestry influences this association is critical and requires further studies.

Rare autosomal copy number variations in early-onset familial Alzheimer's disease

Over 200 rare and fully penetrant pathogenic mutations in amyloid precursor protein (APP), presenilin 1 and 2 (PSEN1 and PSEN2) cause a subset of early-onset familial Alzheimer's disease (EO-FAD). Of these, 21 cases of EO-FAD families carrying unique APP locus duplications remain the only pathogenic copy number variations (CNVs) identified to date in Alzheimer's disease (AD). Using high-density DNA microarrays, we performed a comprehensive genome-wide analysis for the presence of rare CNVs in 261 EO-FAD and early/mixed-onset pedigrees. Our analysis revealed 10 novel private CNVs in 10 EO-FAD families overlapping a set of genes that includes: A2BP1, ABAT, CDH2, CRMP1, DMRT1, EPHA5, EPHA6, ERMP1, EVC, EVC2, FLJ35024 and VLDLR. In addition, CNVs encompassing two known frontotemporal dementia genes, CHMP2B and MAPT were found. To our knowledge, this is the first study reporting rare gene-rich CNVs in EO-FAD and early/mixed-onset AD that are likely to underlie pathogenicity in familial AD and perhaps related dementias.

A genome-wide assessment of the role of untagged copy number variants in type 1 diabetes

Genome-wide association studies (GWAS) for type 1 diabetes (T1D) have successfully identified more than 40 independent T1D associated tagging single nucleotide polymorphisms (SNPs). However, owing to technical limitations of copy number variants (CNVs) genotyping assays, the assessment of the role of CNVs has been limited to the subset of these in high linkage disequilibrium with tag SNPs. The contribution of untagged CNVs, often multi-allelic and difficult to genotype using existing assays, to the heritability of T1D remains an open question. To investigate this issue, we designed a custom comparative genetic hybridization array (aCGH) specifically designed to assay untagged CNV loci identified from a variety of sources. To overcome the technical limitations of the case control design for this class of CNVs, we genotyped the Type 1 Diabetes Genetics Consortium (T1DGC) family resource (representing 3,903 transmissions from parents to affected offspring) and used an association testing strategy that does not necessitate obtaining discrete genotypes. Our design targeted 4,309 CNVs, of which 3,410 passed stringent quality control filters. As a positive control, the scan confirmed the known T1D association at the INS locus by direct typing of the 5′ variable number of tandem repeat (VNTR) locus. Our results clarify the fact that the disease association is indistinguishable from the two main polymorphic allele classes of the INS VNTR, class I-and class III. We also identified novel technical artifacts resulting into spurious associations at the somatically rearranging loci, T cell receptor, TCRA/TCRD and TCRB, and Immunoglobulin heavy chain, IGH, loci on chromosomes 14q11.2, 7q34 and 14q32.33, respectively. However, our data did not identify novel T1D loci. Our results do not support a major role of untagged CNVs in T1D heritability.

For many complex traits, and in particular type 1 diabetes (T1D), the genome-wide association study (GWAS) design has been successful at detecting a large number of loci that contribute disease risk. However, in the case of T1D as well as almost all other traits, the sum of these loci does not fully explain the heritability estimated from familial studies. This observation raises the possibility that additional variants exist but have not yet been found because they have not effectively been targeted by the GWAS design. Here, we focus on a specific class of large deletions/duplications called copy number variants (CNVs), and more precisely to the subset of these loci that mutate rapidly, which are highly polymorphic. A consequence of this high level of polymorphism is that these variants have typically not been captured by previous GWAS studies. We use a family based design that is optimized to capture these previously untested variants. We then perform a genome-wide scan to assess their contribution to T1D. Our scan was technically successful but did not identify novel associations. This suggests that little was missed by the GWAS strategy, and that the remaining heritability of T1D is most likely driven by a large number of variants, either rare of common, but with a small individual contribution to disease risk.

Viral and cellular factors underlying neuropathogenesis in HIV associated neurocognitive disorders (HAND)

As the HIV-1 epidemic enters its fourth decade, HIV-1 associated neurological disorders (HAND) continue to be a major concern in the infected population, despite the widespread use of anti-retroviral therapy. Advancing age and increased life expectancy of the HIV-1 infected population have been shown to increase the risk of cognitive dysfunction. Over the past 10 years, there has been a significant progress in our understanding of the mechanisms and the risk factors involved in the development of HAND. Key events that lead up to neuronal damage in HIV-1 infected individuals can be categorized based on the interaction of HIV-1 with the various cell types, including but not limited to macrophages, brain endothelial cells, microglia, astrocytes and the neurons. This review attempts to decipher these interactions, beginning with HIV-1 infection of macrophages and ultimately resulting in the release of neurotoxic viral and host products. These include: interaction with endothelial cells, resulting in the impairment of the blood brain barrier; interaction with the astrocytes, leading to metabolic and neurotransmitter imbalance; interactions with resident immune cells in the brain, leading to release of toxic cytokines and chemokines. We also review the mechanisms underlying neuronal damage caused by the factors mentioned above. We have attempted to bring together recent findings in these areas to help appreciate the viral and host factors that bring about neurological dysfunction. In addition, we review host factors and viral genotypic differences that affect phenotypic pathological outcomes, as well as recent advances in treatment options to specifically address the neurotoxic mechanisms in play.

Association between copy number variation losses and alcohol dependence across African American and European American ethnic groups

BACKGROUND

Copy number variations (CNVs) are structural genetic mutations consisting of segmental gains or losses in DNA sequence. Although CNVs contribute substantially to genomic variation, few genetic and imaging studies report association of CNVs with alcohol dependence (AD). Our purpose is to find evidence of this association across ethnic populations and genders. This work is the first AD-CNV study across ethnic groups and the first to include the African American (AA) population.

METHODS

This study considers 2 CNV data sets, one for discovery (2,345 samples) and the other for validation (239 samples), both including subjects with AD and healthy controls of European and African ancestry. Our analysis assesses the association between AD and CNV losses across ethnic groups and gender by examining the effect of overall losses across the whole genome, collective losses within individual cytogenetic bands, and specific losses in CNV regions.

RESULTS

Results from the discovery data set showed an association between CNV losses within 16q12.2 and AD diagnosis (p = 4.53 × 10(-3) ). An overlapping CNV region from the validation data set exhibited the same direction of effect with respect to AD (p = 0.051). This CNV region affects the genes CES1p1 and CES1, which are members of the carboxylesterase (CES) family. The enzyme encoded by CES1 is a major liver enzyme that typically catalyzes the decomposition of ester into alcohol and carboxylic acid and is involved in drug or xenobiotics, fatty acid, and cholesterol metabolisms. In addition, the most significantly associated CNV region was located at 9p21.2 (p = 1.9 × 10(-3) ) in our discovery data set. Although not observed in the validation data set, probably due to small sample size, this result might hold potential connection to AD given its connection with neuronal death. In contrast, we did not find any association between AD and the overall total losses or the collective losses within individual cytogenetic bands.

CONCLUSIONS

Overall, our study provides evidence that the specific CNVs at 16q12.2 contribute to the development of alcoholism in AA and European American populations.

Technical considerations for genotyping multi-allelic copy number variation (CNV), in regions of segmental duplication

BACKGROUND

Intrachromosomal segmental duplications provide the substrate for non-allelic homologous recombination, facilitating extensive copy number variation in the human genome. Many multi-copy gene families are embedded within genomic regions with high levels of sequence identity (>95%) and therefore pose considerable analytical challenges. In some cases, the complexity involved in analyzing such regions is largely underestimated. Rapid, cost effective analysis of multi-copy gene regions have typically implemented quantitative approaches, however quantitative data are not an absolute means of certainty. Therefore any technique prone to degrees of measurement error can produce ambiguous results that may lead to spurious associations with complex disease.

RESULTS

In this study we have focused on testing the accuracy and reproducibility of quantitative analysis techniques. With reference to the C-C Chemokine Ligand-3-like-1 (CCL3L1) gene, we performed analysis using real-time Quantitative PCR (QPCR), Multiplex Ligation-dependent Probe Amplification (MLPA) and Paralogue Ratio Test (PRT). After controlling for potential outside variables on assay performance, including DNA concentration, quality, preparation and storage conditions, we find that real-time QPCR produces data that does not cluster tightly around copy number integer values, with variation substantially greater than that of the MLPA or PRT systems. We find that the method of rounding real-time QPCR measurements can potentially lead to mis-scoring of copy number genotypes and suggest caution should be exercised in interpreting QPCR data.

CONCLUSIONS

We conclude that real-time QPCR is inherently prone to measurement error, even under conditions that would seem favorable for association studies. Our results indicate that potential variability in the physicochemical properties of the DNA samples cannot solely explain the poor performance exhibited by the real-time QPCR systems. We recommend that more robust approaches such as PRT or MLPA should be used to genotype multi-allelic copy number variation in disease association studies and suggest several approaches which can be implemented to ensure the quality of the copy number typing using quantitative methods.

Gene-rich large deletions are overrepresented in POAG patients of Indian and Caucasian origins

PURPOSE

Large copy number variations (CNV) can contribute to increased burden for neurodegenerative diseases. In this study, we analyzed the genome-wide burden of large CNVs > 100 kb in primary open angle glaucoma (POAG), a neurodegenerative disease of the eye that is the largest cause of irreversible blindness.

METHODS

Genome-wide analysis of CNVs > 100 kb were analyzed in a total of 1720 individuals, including an Indian cohort (347 POAG cases and 345 controls) and a Caucasian cohort (624 cases and 404 controls). All the CNV data were obtained from experiments performed on Illumina 660W-Quad (infinium) arrays.

RESULTS

We observed that for both the populations CNVs > 1 Mb was significantly enriched for gene-rich regions unique to the POAG cases (P < 10(-11)). In the Indian cohort CNVs > 1 Mb (39 calls) in patients influenced 125 genes while in controls 31 such CNVs influenced only 5 genes with no overlap. In both cohorts we observed 1.9-fold gene enrichment in patients for deletions compared to duplications, while such a bias was not observed in controls (0.3-fold). Overall duplications > 1 Mb were more than deletions (Del/Dup = 0.82) confirming that the enrichment of gene-rich deletions in patients was associated with the disease. Of the 39 CNVs > 1 Mb from Indian patients, 28 (72%) also were implicated in other neurodegenerative disorders, like autism, schizophrenia, sensorineural hearing loss, and so forth. We found one large duplication encompassing CNTN4 gene in Indian and Caucasian POAG patients that was absent in the controls.

CONCLUSIONS

To our knowledge, our study is the first report on large CNV bias for gene-rich regions in glaucomatous neurodegeneration, implicating its impact across populations of contrasting ethnicities. We identified CNTN4 as a novel candidate gene for POAG.

Genetic, transcriptomic, and epigenetic studies of HIV-associated neurocognitive disorder

The Human Genome Project, coupled with rapidly evolving high-throughput technologies, has opened the possibility of identifying heretofore unknown biological processes underlying human disease. Because of the opaque nature of HIV-associated neurocognitive disorder (HAND) neuropathogenesis, the utility of such methods has gained notice among NeuroAIDS researchers. Furthermore, the merging of genetics with other research areas has also allowed for application of relatively nascent fields, such as neuroimaging genomics, and pharmacogenetics, to the context of HAND. In this review, we detail the development of genetic, transcriptomic, and epigenetic studies of HAND, beginning with early candidate gene association studies and culminating in current "omics" approaches that incorporate methods from systems biology to interpret data from multiple levels of biological functioning. Challenges with this line of investigation are discussed, including the difficulty of defining a valid phenotype for HAND. We propose that leveraging known associations between biology and pathology across multiple levels will lead to a more reliable and valid phenotype. We also discuss the difficulties of interpreting the massive and multitiered mountains of data produced by current high-throughput omics assays and explore the utility of systems biology approaches in this regard.

Comparison of a quantitative Real-Time PCR assay and droplet digital PCR for copy number analysis of the CCL4L genes

The controversy surrounding the findings that copy number variation, of the CCL3 encoding genes, influences HIV-1 infection and disease progression has been in part attributed to the variable results obtained from methods used for copy number evaluation. Like CCL3, the genes encoding the CC chemokine CCL4, also a natural ligand of the CCR5 receptor, are found to occur in population-specific multiple copy number and have been shown to play a protective role against HIV-1. This study evaluated the standard method of quantitative Real-Time PCR (qPCR) and droplet digital PCR (ddPCR) for CCL4L gene copy number determination. The CCL4 encoding genes are CCL4, occurring in two copies per diploid genome (pdg), and the non-allelic CCL4L genes, comprised of CCL4L1 and CCL4L2, which are both found in multiple copies pdg. Copy number of CCL4L, CCL4L1 and CCL4L2 was determined in a cohort of HIV-1-uninfected individuals from the South African Black (n=23) and Caucasian (n=32) population groups using qPCR and ddPCR. A stronger correlation between the number of CCL4L copies and the sum of CCL4L1 and CCL4L2 copies generated by ddPCR (r=0.99, p<0.0001) compared to qPCR (r=0.87, p<0.0001) was observed. Real-Time qPCR exhibited greater inaccuracy at higher copy numbers which is particularly relevant to our cohort of Black individuals who have a higher range of CCL4L copies (3-6) compared to Caucasians (0-4) and a higher population median (4 and 2, respectively). Medians and ranges of CCL4L1 (Black: 2, 0-4, Caucasian: 0, 0-2) and CCL4L2 (Black: 2, 1-5, Caucasian: 2, 0-3) were also higher in the Black population. Droplet digital PCR was shown to be a far superior method to qPCR for assessment of CCL4 gene copy number variation, the accuracy of which is essential for studies of the contribution of variable gene copy number to phenotypic outcomes of host infection and disease course.

Simultaneous copy number losses within multiple subtelomeric regions in early-onset type 2 diabetes mellitus

Genetic factors play very important roles in the onset and progression of type 2 diabetes mellitus (T2DM). However, the genetic factors correlating with T2DM onset have not as yet been fully clarified. We previously found that copy number losses in the subtelomeric region on chromosome 4p16.3 were detected in early-onset Japanese T2DM patients (onset age <35 years) at a high frequency. Herein, we additionally found two novel copy number losses within the subtelomeric regions on chromosomes 16q24.2-3 and 22q13.31-33, which have significant associations with early-onset Japanese T2DM. The associations were statistically significant by Fisher's exact tests with P values of 5.19×10⁻³ and 1.81×10⁻³ and odds ratios of 5.7 and 4.4 for 16q24.2-3 and 22q13.31-33, respectively. Furthermore, copy number variation (CNV) analysis of the whole genome using the CNV BeadChip system verified simultaneous copy number losses in all three subtelomeric regions in 11 of our 100 T2DM subjects, while none of 100 non-diabetic controls showed the copy number losses in all three regions. Our results suggest that the mechanism underlying induction of CNVs is involved in the pathogenesis of early-onset T2DM. Thus, copy number losses within multiple subtelomeric regions are strongly associated with early-onset T2DM and examination of simultaneous CNVs in these three regions may lead to the development of an accurate and selective procedure for detecting genetic susceptibility to T2DM.

Analytical pharmacology and allosterism: the importance of quantifying drug parameters in drug discovery

Allosteric ligands bind to receptors at sites separate from those binding endogenous ligands; this allows for a wide range of effects from antagonism to potentiation to direct agonism. This paper discusses techniques to quantify complex allosteric behaviors that yield parameters to characterize direct effects (tB, the efficacy of allosteric agonists), a and b (the effect of the allosteric ligand on endogenous ligand affinity and efficacy respectively). These parameters are independent of the system used to determine them and thus can be used to predict allosteric effects in all systems.

Low copy number of the salivary amylase gene predisposes to obesity

Common multi-allelic copy number variants (CNVs) appear enriched for phenotypic associations compared to their biallelic counterparts. Here we investigated the influence of gene dosage effects on adiposity through a CNV association study of gene expression levels in adipose tissue. We identified significant association of a multi-allelic CNV encompassing the salivary amylase gene (AMY1) with body mass index (BMI) and obesity, and we replicated this finding in 6,200 subjects. Increased AMY1 copy number was positively associated with both amylase gene expression (P = 2.31 × 10(-14)) and serum enzyme levels (P < 2.20 × 10(-16)), whereas reduced AMY1 copy number was associated with increased BMI (change in BMI per estimated copy = -0.15 (0.02) kg/m(2); P = 6.93 × 10(-10)) and obesity risk (odds ratio (OR) per estimated copy = 1.19, 95% confidence interval (CI) = 1.13-1.26; P = 1.46 × 10(-10)). The OR value of 1.19 per copy of AMY1 translates into about an eightfold difference in risk of obesity between subjects in the top (copy number > 9) and bottom (copy number < 4) 10% of the copy number distribution. Our study provides a first genetic link between carbohydrate metabolism and BMI and demonstrates the power of integrated genomic approaches beyond genome-wide association studies.

Piecewise-constant and low-rank approximation for identification of recurrent copy number variations

MOTIVATION

The post-genome era sees urgent need for more novel approaches to extracting useful information from the huge amount of genetic data. The identification of recurrent copy number variations (CNVs) from array-based comparative genomic hybridization (aCGH) data can help understand complex diseases, such as cancer. Most of the previous computational methods focused on single-sample analysis or statistical testing based on the results of single-sample analysis. Finding recurrent CNVs from multi-sample data remains a challenging topic worth further study.

RESULTS

We present a general and robust method to identify recurrent CNVs from multi-sample aCGH profiles. We express the raw dataset as a matrix and demonstrate that recurrent CNVs will form a low-rank matrix. Hence, we formulate the problem as a matrix recovering problem, where we aim to find a piecewise-constant and low-rank approximation (PLA) to the input matrix. We propose a convex formulation for matrix recovery and an efficient algorithm to globally solve the problem. We demonstrate the advantages of PLA compared with alternative methods using synthesized datasets and two breast cancer datasets. The experimental results show that PLA can successfully reconstruct the recurrent CNV patterns from raw data and achieve better performance compared with alternative methods under a wide range of scenarios.

AVAILABILITY AND IMPLEMENTATION

The MATLAB code is available at http://bioinformatics.ust.hk/pla.zip.

Adaptive gene amplification as an intermediate step in the expansion of virus host range

The majority of recently emerging infectious diseases in humans is due to cross-species pathogen transmissions from animals. To establish a productive infection in new host species, viruses must overcome barriers to replication mediated by diverse and rapidly evolving host restriction factors such as protein kinase R (PKR). Many viral antagonists of these restriction factors are species specific. For example, the rhesus cytomegalovirus PKR antagonist, RhTRS1, inhibits PKR in some African green monkey (AGM) cells, but does not inhibit human or rhesus macaque PKR. To model the evolutionary changes necessary for cross-species transmission, we generated a recombinant vaccinia virus that expresses RhTRS1 in a strain that lacks PKR inhibitors E3L and K3L (VVΔEΔK+RhTRS1). Serially passaging VVΔEΔK+RhTRS1 in minimally-permissive AGM cells increased viral replication 10- to 100-fold. Notably, adaptation in these AGM cells also improved virus replication 1000- to 10,000-fold in human and rhesus cells. Genetic analyses including deep sequencing revealed amplification of the rhtrs1 locus in the adapted viruses. Supplying additional rhtrs1 in trans confirmed that amplification alone was sufficient to improve VVΔEΔK+RhTRS1 replication. Viruses with amplified rhtrs1 completely blocked AGM PKR, but only partially blocked human PKR, consistent with the replication properties of these viruses in AGM and human cells. Finally, in contrast to AGM-adapted viruses, which could be serially propagated in human cells, VVΔEΔK+RhTRS1 yielded no progeny virus after only three passages in human cells. Thus, rhtrs1 amplification in a minimally permissive intermediate host was a necessary step, enabling expansion of the virus range to previously nonpermissive hosts. These data support the hypothesis that amplification of a weak viral antagonist may be a general evolutionary mechanism to permit replication in otherwise resistant host species, providing a molecular foothold that could enable further adaptations necessary for efficient replication in the new host.

The spread of microbes from animals to humans has been responsible for most recently emerging human infectious diseases, including AIDS, bird flu, and SARS. Therefore, understanding the evolutionary and molecular mechanisms underlying cross-species transmission is of critical importance for public health. After entering a new host cell, the success of a virus depends on its ability to overcome antiviral factors in the cell, such as protein kinase R (PKR). To investigate the process of virus transmission between species, we employed a recombinant vaccinia virus (VVΔEΔK+RhTRS1) expressing the rhesus cytomegalovirus PKR antagonist RhTRS1. This protein inhibits some African green monkey (AGM) PKRs; however, it does not inhibit human or rhesus variants of PKR. Serial passaging VVΔEΔK+RhTRS1 in RhTRS1-resistant AGM cells resulted in rhtrs1 duplication in the viral genome, which improved VVΔEΔK+RhTRS1 replication in AGM cells. Remarkably, rhtrs1 duplication also enhanced virus replication in human and rhesus cells. In contrast, passage of VVΔEΔK+RhTRS1 in human cells, without prior adaptation in AGM cells, did not improve VVΔEΔK+RhTRS1 replication. These results support the hypothesis that amplification of a weak viral antagonist of a host defense protein in one species may enable cross-species transmission into new hosts that are nonpermissive to the initial virus.

T cell interleukin-15 surface expression in chimpanzees infected with human immunodeficiency virus

Interleukin-15 (IL-15) contributes to natural killer cell development and immune regulation. However, IL-15 and interferon-gamma (IFN-γ) production are significantly reduced during progression to AIDS. We have previously reported that HIV infected chimpanzees (Pan troglodytes) express CD3-CD8+ IFN-γ+ natural killer (NK) cells with an inverse correlation to plasma HIV viral load. To expand on our initial study, we examined a larger population of HIV infected chimpanzees (n=10). Whole blood flow cytometry analyses showed that recombinant gp120 (rgp120) or recombinant IL-15 induces specific CD3-CD8+ IFN-γ+ NK cells at higher levels than CD3+CD8+ IFN-γ+ T cells in HIV infected specimens. Interestingly, peripheral blood T cells exhibited 0.5-3% IL-15 surface Tcell/NKT cell phenotypes, and rIL-15 stimulation significantly (P<0.007) up-regulated CD4+CD25+ T cell expression. Importantly, these data demonstrate novel T cell interleukin-15 expression and indicate a plausible regulatory mechanism for this cell-type during viral infection.

Copy number variations burden on miRNA genes reveals layers of complexities involved in the regulation of pathways and phenotypic expression

MicroRNAs are involved in post-transcriptional down-regulation of gene expression. Variations in miRNA genes can severely affect downstream-regulated genes and their pathways. However, population-specific burden of CNVs on miRNA genes and the complexities created towards the phenotype is not known. From a total of 44109 CNVs investigated from 1715 individuals across 12 populations using high-throughput arrays, 4007 miRNA-CNVs (∼ 9%) consisting 6542 (∼ 5%) miRNA genes with a total of 333 (∼ 5%) singleton miRNA genes were identified. We found miRNA-CNVs across the genomes of individuals showing multiple hits in many targets, co-regulated under the same pathway. This study proposes four mechanisms unraveling the many complexities in miRNA genes, targets and co-regulated miRNA genes towards establishment of phenotypic diversity.

The application of genome-wide SNP genotyping methods in studies on livestock genomes

Animal genomics is currently undergoing dynamic development, which is driven by the flourishing of high-throughput genome analysis methods. Recently, a large number of animals has been genotyped with the use of whole-genome genotyping assays in the course of genomic selection programmes. The results of such genotyping can also be used for studies on different aspects of livestock genome functioning and diversity. In this article, we review the recent literature concentrating on various aspects of animal genomics, including studies on linkage disequilibrium, runs of homozygosity, selection signatures, copy number variation and genetic differentiation of animal populations. Our work is aimed at providing insight into certain achievements of animal genomics and to arouse interest in basic research on the complexity and structure of the genomes of livestock.

Evaluation of copy number variation detection for a SNP array platform

Background

Copy Number Variations (CNVs) are usually inferred from Single Nucleotide Polymorphism (SNP) arrays by use of some software packages based on given algorithms. However, there is no clear understanding of the performance of these software packages; it is therefore difficult to select one or several software packages for CNV detection based on the SNP array platform.

We selected four publicly available software packages designed for CNV calling from an Affymetrix SNP array, including Birdsuite, dChip, Genotyping Console (GTC) and PennCNV. The publicly available dataset generated by Array-based Comparative Genomic Hybridization (CGH), with a resolution of 24 million probes per sample, was considered to be the “gold standard”. Compared with the CGH-based dataset, the success rate, average stability rate, sensitivity, consistence and reproducibility of these four software packages were assessed compared with the “gold standard”. Specially, we also compared the efficiency of detecting CNVs simultaneously by two, three and all of the software packages with that by a single software package.

Results

Simply from the quantity of the detected CNVs, Birdsuite detected the most while GTC detected the least. We found that Birdsuite and dChip had obvious detecting bias. And GTC seemed to be inferior because of the least amount of CNVs it detected. Thereafter we investigated the detection consistency produced by one certain software package and the rest three software suits. We found that the consistency of dChip was the lowest while GTC was the highest. Compared with the CNVs detecting result of CGH, in the matching group, GTC called the most matching CNVs, PennCNV-Affy ranked second. In the non-overlapping group, GTC called the least CNVs. With regards to the reproducibility of CNV calling, larger CNVs were usually replicated better. PennCNV-Affy shows the best consistency while Birdsuite shows the poorest.

Conclusion

We found that PennCNV outperformed the other three packages in the sensitivity and specificity of CNV calling. Obviously, each calling method had its own limitations and advantages for different data analysis. Therefore, the optimized calling methods might be identified using multiple algorithms to evaluate the concordance and discordance of SNP array-based CNV calling.

Copy number variations of HLA-DRB5 is associated with systemic lupus erythematosus risk in Chinese Han population

Systemic lupus erythematosus (SLE) is a polygenic, systemic, autoimmune disease. Copy number variants (CNVS) have been discovered to be associated with a number of complex disorders. We undertook the current study to explore the potential associations between genomic CNVS and SLE in Chinese Han population. In the discovery stage, seven SLE patients were examined with the high-density comparative genomic hybridization microarrays in the screening test for SLE associated CNVS. Then, in the validation stage, 135 SLE patients and 219 matched healthy subjects were investigated for the CNVS of gene HLA-DRB5 by AccuCopy™ technology. Quantitative polymerase chain reaction was carried out to determine the copy number (CN) and mRNA level of HLA-DRB5 in SLE patients. Although the mRNA level of HLA-DRB5 between the CN deletion group and the CN normal group in SLE patients was not statistically positive (P = 0.46), our results still showed more CN of HLA-DRB5 in SLE patients than in healthy controls (P = 3.98 × 10(-6)). Odds ratio for CN deletion was 0.38 (95% confidence interval (CI), 0.23-0.61, P = 7.79 × 10(-5)) and for CN duplication was 1.89 (95% CI, 0.56-7.66, P = 0.37), respectively. These findings indicated that CNVS of HLA-DRB5 was associated with the risk of SLE, and CN deletion appeared to be protective for SLE.

Gene duplication and phenotypic changes in the evolution of mammalian metabolic networks

Metabolic networks attempt to describe the complete suite of biochemical reactions available to an organism. One notable feature of these networks in mammals is the large number of distinct proteins that catalyze the same reaction. While the existence of these isoenzymes has long been known, their evolutionary significance is still unclear. Using a phylogenetically-aware comparative genomics approach, we infer enzyme orthology networks for sixteen mammals as well as for their common ancestors. We find that the pattern of isoenzymes copy-number alterations (CNAs) in these networks is suggestive of natural selection acting on the retention of certain gene duplications. When further analyzing these data with a machine-learning approach, we found that that the pattern of CNAs is also predictive of several important phenotypic traits, including milk composition and geographic range. Integrating tools from network analyses, phylogenetics and comparative genomics both allows the prediction of phenotypes from genetic data and represents a means of unifying distinct biological disciplines.

Identification of medium-sized copy number alterations in whole-genome sequencing

The genome-wide discoveries such as detection of copy number alterations (CNA) from high-throughput whole-genome sequencing data enabled new developments in personalized medicine. The CNAs have been reported to be associated with various diseases and cancers including acute myeloid leukemia. However, there are multiple challenges to the use of current CNA detection tools that lead to high false-positive rates and thus impede widespread use of such tools in cancer research. In this paper, we discuss these issues and propose possible solutions. First, since the entire genome cannot be mapped due to some regions lacking sequence uniqueness, current methods cannot be appropriately adjusted to handle these regions in the analyses. Thus, detection of medium-sized CNAs is also being directly affected by these mappability problems. The requirement for matching control samples is also an important limitation because acquiring matching controls might not be possible or might not be cost efficient. Here we present an approach that addresses these issues and detects medium-sized CNAs in cancer genomes by (1) masking unmappable regions during the initial CNA detection phase, (2) using pool of a few normal samples as control, and (3) employing median filtering to adjust CNA ratios to its surrounding coverage and eliminate false positives.

Copy number polymorphism in plant genomes

Copy number variants (CNVs) are genomic rearrangements resulting from gains or losses of DNA segments. Typically, the term refers to rearrangements of sequences larger than 1 kb. This type of polymorphism has recently been shown to be a key contributor to intra-species genetic variation, along with single-nucleotide polymorphisms and short insertion-deletion polymorphisms. Over the last decade, a growing number of studies have highlighted the importance of copy number variation (CNV) as a factor affecting human phenotype and individual CNVs have been linked to risks for severe diseases. In plants, the exploration of the extent and role of CNV is still just beginning. Initial genomic analyses indicate that CNVs are prevalent in plants and have greatly affected plant genome evolution. Many CNV events have been observed in outcrossing and autogamous species. CNVs are usually found on all chromosomes, with CNV hotspots interspersed with regions of very low genetic variation. Although CNV is mainly associated with intergenic regions, many CNVs encompass protein-coding genes. The collected data suggest that CNV mainly affects the members of large families of functionally redundant genes. Thus, the effects of individual CNV events on phenotype are usually modest. Nevertheless, there are many cases in which CNVs for specific genes have been linked to important traits such as flowering time, plant height and resistance to biotic and abiotic stress. Recent reports suggest that CNVs may form rapidly in response to stress.

Emerging paradigms in GPCR allostery: implications for drug discovery

Allosteric ligands bind to G protein-coupled receptors (GPCRs; also known as seven-transmembrane receptors) at sites that are distinct from the sites to which endogenous ligands bind. The existence of allosteric ligands has enriched the ways in which the functions of GPCRs can be manipulated for potential therapeutic benefit, yet the complexity of their actions provides both challenges and opportunities for drug screening and development. Converging avenues of research in areas such as biased signalling by allosteric ligands and the mechanisms by which allosteric ligands modulate the effects of diverse endogenous ligands have provided new insights into how interactions between allosteric ligands and GPCRs could be exploited for drug discovery. These new findings have the potential to alter how screening for allosteric drugs is performed and may increase the chances of success in the development of allosteric modulators as clinical lead compounds.