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

Arabidopsis thaliana population analysis reveals high plasticity of the genomic region spanning MSH2, AT3G18530 and AT3G18535 genes and provides evidence for NAHR-driven recurrent CNV events occurring in this location

Background

Intraspecies copy number variations (CNVs), defined as unbalanced structural variations of specific genomic loci, ≥1 kb in size, are present in the genomes of animals and plants. A growing number of examples indicate that CNVs may have functional significance and contribute to phenotypic diversity. In the model plant Arabidopsis thaliana at least several hundred protein-coding genes might display CNV; however, locus-specific genotyping studies in this plant have not been conducted.

Results

We analyzed the natural CNVs in the region overlapping MSH2 gene that encodes the DNA mismatch repair protein, and AT3G18530 and AT3G18535 genes that encode poorly characterized proteins. By applying multiplex ligation-dependent probe amplification and droplet digital PCR we genotyped those genes in 189 A. thaliana accessions. We found that AT3G18530 and AT3G18535 were duplicated (2–14 times) in 20 and deleted in 101 accessions. MSH2 was duplicated in 12 accessions (up to 12-14 copies) but never deleted. In all but one case, the MSH2 duplications were associated with those of AT3G18530 and AT3G18535. Considering the structure of the CNVs, we distinguished 5 genotypes for this region, determined their frequency and geographical distribution. We defined the CNV breakpoints in 35 accessions with AT3G18530 and AT3G18535 deletions and tandem duplications and showed that they were reciprocal events, resulting from non-allelic homologous recombination between 99 %-identical sequences flanking these genes. The widespread geographical distribution of the deletions supported by the SNP and linkage disequilibrium analyses of the genomic sequence confirmed the recurrent nature of this CNV.

Conclusions

We characterized in detail for the first time the complex multiallelic CNV in Arabidopsis genome. The region encoding MSH2, AT3G18530 and AT3G18535 genes shows enormous variation of copy numbers among natural ecotypes, being a remarkable example of high Arabidopsis genome plasticity. We provided the molecular insight into the mechanism underlying the recurrent nature of AT3G18530-AT3G18535 duplications/deletions. We also performed the first direct comparison of the two leading experimental methods, suitable for assessing the DNA copy number status. Our comprehensive case study provides foundation information for further analyses of CNV evolution in Arabidopsis and other plants, and their possible use in plant breeding.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-016-3221-1) contains supplementary material, which is available to authorized users.

SRBreak: A Read-Depth and Split-Read Framework to Identify Breakpoints of Different Events Inside Simple Copy-Number Variable Regions

Copy-number variation (CNV) has been associated with increased risk of complex diseases. High-throughput sequencing (HTS) technologies facilitate the detection of copy-number variable regions (CNVRs) and their breakpoints. This helps in understanding genome structure as well as their evolution process. Various approaches have been proposed for detecting CNV breakpoints, but currently it is still challenging for tools based on a single analysis method to identify breakpoints of CNVs. It has been shown, however, that pipelines which integrate multiple approaches are able to report more reliable breakpoints. Here, based on HTS data, we have developed a pipeline to identify approximate breakpoints (±10 bp) relating to different ancestral events within a specific CNVR. The pipeline combines read-depth and split-read information to infer breakpoints, using information from multiple samples to allow an imputation approach to be taken. The main steps involve using a normal mixture model to cluster samples into different groups, followed by simple kernel-based approaches to maximize information obtained from read-depth and split-read approaches, after which common breakpoints of groups are inferred. The pipeline uses split-read information directly from CIGAR strings of BAM files, without using a re-alignment step. On simulated data sets, it was able to report breakpoints for very low-coverage samples including those for which only single-end reads were available. When applied to three loci from existing human resequencing data sets (NEGR1, LCE3, IRGM) the pipeline obtained good concordance with results from the 1000 Genomes Project (92, 100, and 82%, respectively). The package is available at https://github.com/hoangtn/SRBreak, and also as a docker-based application at https://registry.hub.docker.com/u/hoangtn/srbreak/.

A Genome-Wide Association Study to Identify Potential Germline Copy Number Variants for Sporadic Breast Cancer Susceptibility

Breast cancer (BC) predisposition in populations arises from both genetic and nongenetic risk factors. Structural variations such as copy number variations (CNVs) are heritable determinants for disease susceptibility. The primary objectives of this study are (1) to identify CNVs associated with sporadic BC using a genome-wide association study (GWAS) design; (2) to utilize 2 distinct CNV calling algorithms to identify concordant CNVs as a strategy to reduce false positive associations in the hypothesis-generating GWAS discovery phase, and (3) to identify potential candidate CNVs for follow-up replication studies. We used Affymetrix SNP Array 6.0 data profiled on Caucasian subjects (422 cases/348 controls) to call CNVs using algorithms implemented in Nexus Copy Number and Partek Genomics Suite software. Nexus algorithm identified CNVs associated with BC (731 autosomal CNVs with >5% frequency in the total sample and Q < 0.05). Thirteen CNVs were identified when Partek algorithm-called CNVs were overlapped with Nexus-identified CNVs; these CNVs showed concordances for frequency, effect size, and direction. Coding genes present within BC-associated CNVs were known to play a role in disease etiology and prognosis. Long noncoding RNAs identified within CNVs showed tissue-specific expression, indicating potential functional relevance of the findings. The identified candidate CNVs warrant independent replication.

Copy number variation of genes involved in the hepatitis C virus-human interactome

Copy number variation (CNV) is a newly discovered form of intra-species genetic polymorphism that is defined as deletions or duplications of genome segments ranging from 1 kbp to several Mbp. CNV accounts for the majority of the genetic variation observed in humans (CNV regions cover more than 10% of the human genome); therefore, it may significantly influence both the phenotype and susceptibility to various diseases. Unfortunately, the impact of CNV on a number of diseases, including hepatitis C virus (HCV) infection, remains largely unexplored. Here, we analyzed 421 human genes encoding proteins that have been shown to interact with HCV proteins or genomic RNA (proteins from the HCV-human interactome). We found that 19 of the 421 candidate genes are located in putative CNV regions. For all of these genes, copy numbers were determined for European, Asiatic and African populations using the multiplex ligation-dependent amplification (MLPA) method. As a result, we identified 4 genes, IGLL1, MLLT4, PDPK1, PPP1R13L, for which the CN-genotype ranged from 1 to 6. All of these genes are involved in host-virus interaction; thus, their polymorphism has a potential impact on the development of HCV infection and/or therapy outcome.

MGP-HMM: Detecting genome-wide CNVs using an HMM for modeling mate pair insertion sizes and read counts

MOTIVATION

Association of Copy Number Variation (CNV) with schizophrenia, autism, developmental disabilities and fatal diseases such as cancer is verified. Recent developments in Next Generation Sequencing (NGS) have facilitated the CNV studies. However, many of the current CNV detection tools are not capable of discriminating tandem duplication from non-tandem duplications.

RESULTS

In this study, we propose MGP-HMM as a tool which besides detecting genome-wide deletions discriminates tandem duplications from non-tandem duplications. MGP-HMM takes mate pair abnormalities into account and predicts the digitized number of tandem or non-tandem copies. Abnormalities in the mate pair directions and insertion sizes, after being mapped to the reference genome, are elucidated using a Hidden Markov Model (HMM). For this purpose, a Mixture Gaussian density with time-dependent parameters is applied for emitting mate pair insertion sizes from HMM states. Indeed, depending on observed abnormalities in mate pair insertion size or its orientation, each component in the mixture density will have different parameters. MGP-HMM also applies a Poisson distribution for modeling read depth data. This parametric modeling of the mate pair reads enables us to estimate the length of CNVs precisely, which is an advantage over methods which rely only on read depth approach for the CNV detection. Hidden state of the proposed HMM is the digitized copy number of a genomic segment and states correspond to the multipliers of the mixture Gaussian components. The accuracy of our model is validated on a set of next generation sequencing real and simulated data and is compared to other tools.

High-resolution arrays reveal burden of copy number variations on Parkinson disease genes associated with increased disease risk in random cohorts

BACKGROUND

Parkinson disease (PD) is a neurological disease responsible for a considerable rate of mortality and morbidity in the society. Since the symptoms of the disease appear much later than the actual onset of neuron degeneration, a majority of cases remain undiagnosed until the manifestation of the symptoms.

OBJECTIVES

In order to investigate the existence of such susceptibility in the population, we analyzed Copy Number Variation (CNV) influences on PD genes in 1715 individuals from 12 different populations.

RESULTS

Overall, 16 CNV-PD genes, 3 known to be causal and 13 associated, were found to be significantly enriched. PARK2, was under heavy burden with ~1% of the population containing CNV in the exonic region. The impact of these genes on the genome and disease pathway was analyzed using several genome analysis tools. Protein interaction network of CNV-PD genes revealed a complex interaction of molecules forming a major hub by the α-Synuclein, whose direct interactors, LRRK2, PARK2 and ATP13A2 are under CNV influence.

CONCLUSIONS

We hypothesize that CNVs may not be the initiating event in the pathogenesis of PD and remain latent until additional secondary hits are acquired and also propose novel genes that may fall under the PD pathway which contribute in pathogenesis.

New quality measure for SNP array based CNV detection

MOTIVATION

Only a few large systematic studies have evaluated the impact of copy number variants (CNVs) on common diseases. Several million individuals have been genotyped on single nucleotide variation arrays, which could be used for genome-wide CNVs association studies. However, CNV calls remain prone to false positives and only empirical filtering strategies exist in the literature. To overcome this issue, we defined a new quality score (QS) estimating the probability of a CNV called by PennCNV to be confirmed by other software.

RESULTS

Out-of-sample comparison showed that the correlation between the consensus CNV status and the QS is twice as high as it is for any previously proposed CNV filters. ROC curves displayed an AUC higher than 0.8 and simulations showed an increase up to 20% in statistical power when using QS in comparison to other filtering strategies. Superior performance was confirmed also for alternative consensus CNV definition and through improving known CNV-trait associations.

AVAILABILITY AND IMPLEMENTATION

http://goo.gl/T6yuFM CONTACT: zoltan.kutalik@unil.ch or aurelien@mace@unil.chSupplementary information: Supplementary data are available at Bioinformatics online.

Impact of genetic variations in C-C chemokine receptors and ligands on infectious diseases

Chemokine receptors and ligands are crucial for extensive immune response against infectious diseases such as malaria, leishmaniasis, HIV and tuberculosis and a wide variety of other diseases. Role of chemokines are evidenced in the activation and regulation of immune cell migration which is important for immune response against diseases. Outcome of disease is determined by complex interaction among pathogen, host genetic variability and surrounding milieu. Variation in expression or function of chemokines caused by genetic polymorphisms could be associated with attenuated immune responses. Exploration of chemokine genetic polymorphisms in therapeutic response, gene regulation and disease outcome is important. Infectious agents in human host alter the expression of chemokines via epigenetic alterations and thus contribute to disease pathogenesis. Although some fragmentary data are available on chemokine genetic variations and their contribution in diseases, no unequivocal conclusion has been arrived as yet. We therefore, aim to investigate the association of CCR5-CCL5 and CCR2-CCL2 genetic polymorphisms with different infectious diseases, transcriptional regulation of gene, disease severity and response to therapy. Furthermore, the role of epigenetics in genes related to chemokines and infectious disease are also discussed.

Copy number variants in the sheep genome detected using multiple approaches

Background

Copy number variants (CNVs) are a type of polymorphism found to underlie phenotypic variation, both in humans and livestock. Most surveys of CNV in livestock have been conducted in the cattle genome, and often utilise only a single approach for the detection of copy number differences. Here we performed a study of CNV in sheep, using multiple methods to identify and characterise copy number changes. Comprehensive information from small pedigrees (trios) was collected using multiple platforms (array CGH, SNP chip and whole genome sequence data), with these data then analysed via multiple approaches to identify and verify CNVs.

Results

In total, 3,488 autosomal CNV regions (CNVRs) were identified in this study, which substantially builds on an initial survey of the sheep genome that identified 135 CNVRs. The average length of the identified CNVRs was 19 kb (range of 1 kb to 3.6 Mb), with shorter CNVRs being more frequent than longer CNVRs. The total length of all CNVRs was 67.6Mbps, which equates to 2.7 % of the sheep autosomes. For individuals this value ranged from 0.24 to 0.55 %, and the majority of CNVRs were identified in single animals. Rather than being uniformly distributed throughout the genome, CNVRs tended to be clustered. Application of three independent approaches for CNVR detection facilitated a comparison of validation rates. CNVs identified on the Roche-NimbleGen 2.1M CGH array generally had low validation rates with lower density arrays, while whole genome sequence data had the highest validation rate (>60 %).

Conclusions

This study represents the first comprehensive survey of the distribution, prevalence and characteristics of CNVR in sheep. Multiple approaches were used to detect CNV regions and it appears that the best method for verifying CNVR on a large scale involves using a combination of detection methodologies. The characteristics of the 3,488 autosomal CNV regions identified in this study are comparable to other CNV regions reported in the literature and provide a valuable and sizeable addition to the small subset of published sheep CNVs.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-016-2754-7) contains supplementary material, which is available to authorized users.

Copy Number Variation of UGT 2B Genes in Indian Families Using Whole Genome Scans

Background and Objectives. Uridine diphospho-glucuronosyltransferase 2B (UGT2B) is a family of genes involved in metabolizing steroid hormones and several other xenobiotics. These UGT2B genes are highly polymorphic in nature and have distinct polymorphisms associated with specific regions around the globe. Copy number variations (CNVs) status of UGT2B17 in Indian population is not known and their disease associations have been inconclusive. It was therefore of interest to investigate the CNV profile of UGT2B genes. Methods. We investigated the presence of CNVs in UGT2B genes in 31 members from eight Indian families using Affymetrix Genome-Wide Human SNP Array 6.0 chip. Results. Our data revealed >50% of the study members carried CNVs in UGT2B genes, of which 76% showed deletion polymorphism. CNVs were observed more in UGT2B17 (76.4%) than in UGT2B15 (17.6%). Molecular network and pathway analysis found enrichment related to steroid metabolic process, carboxylesterase activity, and sequence specific DNA binding. Interpretation and Conclusion. We report the presence of UGT2B gene deletion and duplication polymorphisms in Indian families. Network analysis indicates the substitutive role of other possible genes in the UGT activity. The CNVs of UGT2B genes are very common in individuals indicating that the effect is neutral in causing any suspected diseases.

Structure, function, and genetics of lipoprotein (a)

Lipoprotein (a) [Lp(a)] has attracted the interest of researchers and physicians due to its intriguing properties, including an intragenic multiallelic copy number variation in the LPA gene and the strong association with coronary heart disease (CHD). This review summarizes present knowledge of the structure, function, and genetics of Lp(a) with emphasis on the molecular and population genetics of the Lp(a)/LPA trait, as well as aspects of genetic epidemiology. It highlights the role of genetics in establishing Lp(a) as a risk factor for CHD, but also discusses uncertainties, controversies, and lack of knowledge on several aspects of the genetic Lp(a) trait, not least its function.

A Novel Role for the Receptor of the Complement Cleavage Fragment C5a, C5aR1, in CCR5-Mediated Entry of HIV into Macrophages

The complement system is an ancient pattern recognition system that becomes activated during all stages of HIV infection. Previous studies have shown that C5a can enhance the infection of monocyte-derived macrophages and T cells indirectly through the production of interleukin (IL)-6 and tumor necrosis factor (TNF)-α and the attraction of dendritic cells. C5a exerts its multiple biologic functions mainly through activation of C5a receptor 1 (C5aR1). Here, we assessed the role of C5aR1 as an enhancer of CCR5-mediated HIV infection. We determined CCR5 and C5aR1 heterodimer formation in myeloid cells and the impact of C5aR1 blockade on HIV entry and genomic integration. C5aR1/CCR5 heterodimer formation was identified by immunoprecipitation and western blotting. THP-1 cells and monocyte-derived macrophages (MDM) were infected by R5 laboratory strains or HIV pseudotyped for the vesicular stomatitis virus (VSV) envelope. Levels of integrated HIV were measured by quantitative PCR after targeting of C5aR1 by a C5aR antagonist, neutralizing C5aR1 monoclonal antibody (mAb) or hC5a. C5aR1 was also silenced by specific siRNA prior to viral entry. We found that C5aR1 forms heterodimers with the HIV coreceptor CCR5 in myeloid cells. Targeting C5aR1 significantly decreased integration by R5 viruses but not by VSV-pseudotyped viruses, suggesting that C5aR1 is critical for viral entry. The level of inhibition achieved with C5aR1-blocking reagents was comparable to that of CCR5 antagonists. Mechanistically, C5aR1 targeting decreased CCR5 expression. MDM from CCR5Δ32 homozygous subjects expressed levels of C5aR1 similar to CCR5 WT individuals, suggesting that mere C5aR1 expression is not sufficient for HIV infection. HIV appeared to preferentially enter THP-1 cells expressing high levels of both C5aR1 and CCR5. Targeted reduction of C5aR1 expression in such cells reduced HIV infection by ~50%. Our data thus suggest that C5aR1 acts as an enhancer of CCR5-mediated HIV entry into macrophages, the targeting of which may prove useful to reduce HIV infection by R5 strains.

HIV-associated neurocognitive disorder--pathogenesis and prospects for treatment

In the past two decades, several advancements have improved the care of HIV-infected individuals. Most importantly, the development and deployment of combination antiretroviral therapy (CART) has resulted in a dramatic decline in the rate of deaths from AIDS, so that people living with HIV today have nearly normal life expectancies if treated with CART. The term HIV-associated neurocognitive disorder (HAND) has been used to describe the spectrum of neurocognitive dysfunction associated with HIV infection. HIV can enter the CNS during early stages of infection, and persistent CNS HIV infection and inflammation probably contribute to the development of HAND. The brain can subsequently serve as a sanctuary for ongoing HIV replication, even when systemic viral suppression has been achieved. HAND can remain in patients treated with CART, and its effects on survival, quality of life and everyday functioning make it an important unresolved issue. In this Review, we describe the epidemiology of HAND, the evolving concepts of its neuropathogenesis, novel insights from animal models, and new approaches to treatment. We also discuss how inflammation is sustained in chronic HIV infection. Moreover, we suggest that adjunctive therapies--treatments targeting CNS inflammation and other metabolic processes, including glutamate homeostasis, lipid and energy metabolism--are needed to reverse or improve HAND-related neurological dysfunction.

Defensin gene variation and HIV/AIDS: a comprehensive perspective needed

Both α- and β-defensins have anti-human immunodeficiency virus activity. These defensins achieve human immunodeficiency virus inhibition through a variety of mechanisms, including direct binding with virions, binding to and modulation of host cell-surface receptors with disruption of intracellular signaling, and functioning as chemokines or cytokines to augment and alter adaptive immune responses. Polymorphisms in the defensin genes have been associated with susceptibility to human immunodeficiency virus infection and disease progression. However, the roles that these defensins and their genetic polymorphisms have in influencing human immunodeficiency virus/acquired immunodeficiency syndrome outcomes are not straightforward and, at times, appear contradictory. Differences in populations, study designs, and techniques for genotyping defensin gene polymorphisms may have contributed to this lack of clarity. In addition, a comprehensive approach, where both subfamilies of defensins and their all-inclusive genetic polymorphism profiles are analyzed, is lacking. Such an approach may reveal whether the human immunodeficiency virus inhibitory activities of α- and β-defensins are based on parallel or divergent mechanisms and may provide further insights into how the genetic predisposition for susceptibility or resistance to human immunodeficiency virus/acquired immunodeficiency syndrome is orchestrated between these molecules.

The expression of the chemokine receptor CCR5 in tick-borne encephalitis

Background

Chemokine receptor 5 (CCR5) is hypothesized to drive the lymphocyte migration to central nervous system in flavivirus encephalitis, and the non-functional CCR5Δ32 genetic variant was identified as a risk factor of a West Nile virus infection and of tick-borne encephalitis (TBE). We have attempted to investigate how CCR5 expression corresponds to the clinical course and severity of TBE.

Methods

We have repeatedly studied CCR5 expression in 76 patients during encephalitic and convalescent TBE phase, analyzing its association with clinical features, cerebrospinal fluid (csf) pleocytosis, and concentrations of CCR5 ligands (chemokines CCL3, CCL4, and CCL5) and CCR5 genotype. Fifteen patients with neuroborreliosis, 7 with aseptic meningitis, 17 in whom meningitis/encephalitis had been excluded, and 18 healthy blood donors were studied as controls. Expression of CCR5 was measured cytometrically in blood and csf-activated Th lymphocytes (CD3+CD4+CD45RO+). Concentrations of chemokines in serum and csf were measured immunoenzymatically, and CCR5Δ32 was detected with sequence-specific primers. Data were analyzed with non-parametric tests, and p < 0.05 was considered significant.

Results

The blood expression of CCR5 did neither differ between the groups nor change in the course of TBE. The CCR5 expression in the inflammatory csf was several-fold increased in comparison with blood but lower in TBE than in neuroborreliosis. The csf concentration of CCL5 was increased in TBE, the highest in the most severe presentation (meningoencephalomyelitis) and correlated with pleocytosis. The CCR5Δ32/wt genotype present in 7 TBE patients was associated with a decreased CCR5 expression, but enrichment of csf Th population in CCR5-positive cells and the intrathecal inflammatory response were preserved, without a compensatory increase of CCL5 expression.

Conclusions

We infer CCR5 and CCL5 participate in the response to TBE virus, as well as to other neurotropic pathogens. The intrathecal response to TBE is not hampered in the bearers of a single copy of CCR5Δ32 allele, suggesting that the association of CCR5Δ32 with TBE may be mediated in the periphery at the earlier stage of the infection. Otherwise, a variability of the CCR5 expression in the peripheral blood lymphocytes seems not to be associated with a variable susceptibility to TBE.

Theoretical Aspects of GPCR-Ligand Complex Pharmacology

Over the past 50 years in pharmacology, an understanding of seven transmembrane (7TMR) function has been gained from the comparison of experimental data to receptor models. These models have been constructed from building blocks composed of systems consisting of series and parallel mass action binding reactions. Basic functions such as the the isomerization of receptors upon ligand binding, the sequential binding of receptors to membrane coupling proteins, and the selection of multiple receptor conformations have been combined in various ways to build receptor systems such as the ternary complex, extended ternary complex, and cubic ternary complex models for 7TMR function. Separately, the Black/Leff operational model has furnished an extremely valuable method of quantifying drug agonism. In the past few years, incorporation of the basic allosteric nature of 7TMRs has led to additional useful models of functional receptor allosteric mechanisms; these models yield valuable methods for quantifying allosteric effects. Finally, molecular dynamics has provided yet another new set of models describing the probability of formation of multiple receptor states; these radically new models are extremely useful in the prediction of functionally selective drug effects.

Different Facets of Copy Number Changes: Permanent, Transient, and Adaptive

Chromosomal copy number changes are frequently associated with harmful consequences and are thought of as an underlying mechanism for the development of diseases. However, changes in copy number are observed during development and occur during normal biological processes. In this review, we highlight the causes and consequences of copy number changes in normal physiologic processes as well as cover their associations with cancer and acquired drug resistance. We discuss the permanent and transient nature of copy number gains and relate these observations to a new mechanism driving transient site-specific copy gains (TSSGs). Finally, we discuss implications of TSSGs in generating intratumoral heterogeneity and tumor evolution and how TSSGs can influence the therapeutic response in cancer.

CD26/dipeptidylpeptidase IV-chemokine interactions: double-edged regulation of inflammation and tumor biology

Review of how chemokine processing by CD26/DPP IV regulates leukocyte trafficking.

Post‐translational modification of chemokines is an essential regulatory mechanism to enhance or dampen the inflammatory response. CD26/dipeptidylpeptidase IV, ubiquitously expressed in tissues and blood, removes NH₂‐terminal dipeptides from proteins with a penultimate Pro or Ala. A large number of human chemokines, including CXCL2, CXCL6, CXCL9, CXCL10, CXCL11, CXCL12, CCL3L1, CCL4, CCL5, CCL11, CCL14, and CCL22, are cleaved by CD26; however, the efficiency is clearly influenced by the amino acids surrounding the cleavage site and although not yet proven, potentially affected by the chemokine concentration and interactions with third molecules. NH₂‐terminal cleavage of chemokines by CD26 has prominent effects on their receptor binding, signaling, and hence, in vitro and in vivo biologic activities. However, rather than having a similar result, the outcome of NH₂‐terminal truncation is highly diverse. Either no difference in activity or drastic alterations in receptor recognition/specificity and hence, chemotactic activity are observed. Analogously, chemokine‐dependent inhibition of HIV infection is enhanced (for CCL3L1 and CCL5) or decreased (for CXCL12) by CD26 cleavage. The occurrence of CD26‐processed chemokine isoforms in plasma underscores the importance of the in vitro‐observed CD26 cleavages. Through modulation of chemokine activity, CD26 regulates leukocyte/tumor cell migration and progenitor cell release from the bone marrow, as shown by use of mice treated with CD26 inhibitors or CD26 knockout mice. As chemokine processing by CD26 has a significant impact on physiologic and pathologic processes, application of CD26 inhibitors to affect chemokine function is currently explored, e.g., as add‐on therapy in viral infection and cancer.

Characterization of copy number variants for CCL3L1 gene in rheumatoid arthritis for French trio families and Tunisian cases and controls

Analyses of copy number variants (CNVs) for candidate genes in complex diseases are currently a promising research field. CNVs of C-C chemokine ligand 3-like 1 (CCL3L1) gene are candidate genomic factors in rheumatoid arthritis (RA). We investigated CCL3L1 CNVs association with a case-control study in Tunisians and a transmission analysis in French trio families. Relative copy number (rCN) of CCL3L1 gene was quantified by droplet digital PCR (ddPCR) in 100 French trio families (RA patients and their two parents) and in 166 RA cases and 102 healthy controls from Tunisia. We calculated odds ratio (OR) to investigate association risk for CCL3L1 CNVs in RA. rCN identified varied from 0 to 4 in the French population and from 0 to 7 in the Tunisian population. A significant difference was observed in the distribution of these rCNs between the two populations (p = 2.34 × 10(-10)), as when rCN from French and Tunisian RA patients were compared (p = 2.83 × 10(-5)). CNVs transmission in French RA trios allowed the characterization of genotypes with the presence of tandem duplication and triplication on the same chromosome. RA association tests highlighted a protective effect of rCN = 5 for CCL3L1 gene in the Tunisian population (OR = 0.056; CI 95 % [0.01-0.46]). Characterization of CCL3L1 CNVs with ddPCR methodology highlighted specific CN genotypes in a French family sample. A copy number polymorphism of a RA candidate gene was quantified, and its significant association with RA was revealed in a Tunisian sample.

Rapid Increase in frequency of gene copy-number variants during experimental evolution in Caenorhabditis elegans

BACKGROUND

Gene copy-number variation (CNVs), which provides the raw material for the evolution of novel genes, is widespread in natural populations. We investigated whether CNVs constitute a common mechanism of genetic change during adaptation in experimental Caenorhabditis elegans populations. Outcrossing C. elegans populations with low fitness were evolved for >200 generations. The frequencies of CNVs in these populations were analyzed by oligonucleotide array comparative genome hybridization, quantitative PCR, PCR, DNA sequencing across breakpoints, and single-worm PCR.

RESULTS

Multiple duplications and deletions rose to intermediate or high frequencies in independent populations. Several lines of evidence suggest that these changes were adaptive: (i) copy-number changes reached high frequency or were fixed in a short time, (ii) many independent populations harbored CNVs spanning the same genes, and (iii) larger average size of CNVs in adapting populations relative to spontaneous CNVs. The latter is expected if larger CNVs are more likely to encompass genes under selection for a change in gene dosage. Several convergent CNVs originated in populations descended from different low fitness ancestors as well as high fitness controls.

CONCLUSIONS

We show that gene copy-number changes are a common class of adaptive genetic change. Due to the high rates of origin of spontaneous duplications and deletions, copy-number changes containing the same genes arose readily in independent populations. Duplications that reached high frequencies in these adapting populations were significantly larger in span. Many convergent CNVs may be general adaptations to laboratory conditions. These results demonstrate the great potential borne by CNVs for evolutionary adaptation.

Multilevel analysis of neuropathogenesis of neurocognitive impairment in HIV

The neuropathogenesis of HIV-associated neurocognitive disorders (HAND) remains puzzling. We interrogated several levels of data (host genetic, histopathology, brain viral load, and neurocognitive) to identify histopathological changes most relevant to HAND. The design of the study is a clinicopathological study employing genetic association analyses. Data and brain tissue from 80 HIV-infected adults were used. Markers in monocyte chemoattractant protein-1 (MCP-1), interleukin 1-alpha (IL1-α), macrophage inflammatory protein 1-alpha (MIP1-α), DRD3, DRD2, and apolipoprotein E (ApoE) were genotyped. Microtubule associated protein 2 (MAP2), synaptophysin (SYP), human leukocyte antigen-DR (HLA-DR), glial fibrillary acidic protein (GFAP), amyloid beta (A-Beta), and ionized calcium-binding adaptor molecule-1 (Iba-1) immunoreactivity were quantified in the frontal cortex, putamen, and hippocampus. A composite score for each marker (mean of the three brain regions) was used. Neurocognitive functioning and other clinical variables were determined within 1 year of death. Brain HIV RNA viral load was available for a subset of cases. MAP2 and SYP proved most relevant to neurocognitive functioning. Immunoreactivity of these markers, as well as A-Beta and Iba-1, was correlated with brain HIV RNA viral load. Several genetic markers in combination with other factors predicted histopathology: HIV blood viral load, MIP1-α genotype, and DRD3 genotype predicted Iba-1 immunoreactivity; the duration of infection and IL1-α genotype predicted GFAP immunoreactivity; ApoE genotype and age at death predicted A-Beta immunoreactivity. These data indicate that HIV replication in the brain is the primary driving force leading to neuroinflammation and dysfunctional protein clearance, as reflected by A-Beta and Iba-1. Downstream to these changes are synaptodendritic degeneration, which is the immediate histopathological substrate of the neurocognitive impairment characteristic of HAND. These intermediate histopathological phenotypes are influenced by host genetic polymorphisms in genes encoding cytokines/chemokines, neuronal protein clearance pathways, and dopaminergic factors.

Identifying rhesus macaque gene orthologs using heterospecific human CNV probes

We used the Affymetrix(®) Genome-Wide Human SNP Array 6.0 to identify heterospecific markers and compare copy number and structural genomic variation between humans and rhesus macaques. Over 200,000 human copy number variation (CNV) probes were mapped to a Chinese and an Indian rhesus macaque sample. Observed genomic rearrangements and synteny were in agreement with the results of a previously published genomic comparison between humans and rhesus macaques. Comparisons between each of the two rhesus macaques and humans yielded 206 regions with copy numbers that differed by at least two fold in the Indian rhesus macaque and human, 32 in the Chinese rhesus macaque and human, and 147 in both rhesus macaques. The detailed genomic map and preliminary CNV data are useful for better understanding genetic variation in rhesus macaques, identifying derived changes in human CNVs that may have evolved by selection, and determining the suitability of rhesus macaques as human models for particular biomedical studies.

Population Structure Shapes Copy Number Variation in Malaria Parasites

If copy number variants (CNVs) are predominantly deleterious, we would expect them to be more efficiently purged from populations with a large effective population size (Ne) than from populations with a small Ne. Malaria parasites (Plasmodium falciparum) provide an excellent organism to examine this prediction, because this protozoan shows a broad spectrum of population structures within a single species, with large, stable, outbred populations in Africa, small unstable inbred populations in South America and with intermediate population characteristics in South East Asia. We characterized 122 single-clone parasites, without prior laboratory culture, from malaria-infected patients in seven countries in Africa, South East Asia and South America using a high-density single-nucleotide polymorphism/CNV microarray. We scored 134 high-confidence CNVs across the parasite exome, including 33 deletions and 102 amplifications, which ranged in size from <500 bp to 59 kb, as well as 10,107 flanking, biallelic single-nucleotide polymorphisms. Overall, CNVs were rare, small, and skewed toward low frequency variants, consistent with the deleterious model. Relative to African and South East Asian populations, CNVs were significantly more common in South America, showed significantly less skew in allele frequencies, and were significantly larger. On this background of low frequency CNV, we also identified several high-frequency CNVs under putative positive selection using an FST outlier analysis. These included known adaptive CNVs containing rh2b and pfmdr1, and several other CNVs (e.g., DNA helicase and three conserved proteins) that require further investigation. Our data are consistent with a significant impact of genetic structure on CNV burden in an important human pathogen.

Copy number variations in the genome of the Qatari population

Background

The populations of the Arabian Peninsula remain the least represented in public genetic databases, both in terms of single nucleotide variants and of larger genomic mutations. We present the first high-resolution copy number variation (CNV) map for a Gulf Arab population, using a hybrid approach that integrates array genotyping intensity data and next-generation sequencing reads to call CNVs in the Qatari population.

Methods

CNVs were detected in 97 unrelated Qatari individuals by running two calling algorithms on each of two primary datasets: high-resolution genotyping (Illumina Omni 2.5M) and high depth whole-genome sequencing (Illumina PE 100bp). The four call-sets were integrated to identify high confidence CNV regions, which were subsequently annotated for putative functional effect and compared to public databases of CNVs in other populations. The availability of genome sequence was leveraged to identify tagging SNPs in high LD with common deletions in this population, enabling their imputation from genotyping experiments in the future.

Results

Genotyping intensities and genome sequencing data from 97 Qataris were analyzed with four different algorithms and integrated to discover 16,660 high confidence CNV regions (CNVRs) in the total population, affecting ~28 Mb in the median Qatari genome. Up to 40 % of all CNVs affected genes, including novel CNVs affecting Mendelian disease genes, segregating at different frequencies in the 3 major Qatari subpopulations, including those with Bedouin, Persian/South Asian, and African ancestry. Consistent with high consanguinity levels in the Bedouin subpopulation, we found an increased burden for homozygous deletions in this group. In comparison to known CNVs in the comprehensive Database of Genomic Variants, we found that 5 % of all CNVRs in Qataris were completely novel, with an enrichment of CNVs affecting several known chromosomal disorder loci and genes known to regulate sugar metabolism and type 2 diabetes in the Qatari cohort. Finally, we leveraged the availability of genome sequence to find suitable tagging SNPs for common deletions in this population.

Conclusion

We combine four independently generated datasets from 97 individuals to study CNVs for the first time at high-resolution in a Gulf Arab population.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1991-5) contains supplementary material, which is available to authorized users.

CCR2, CCR5, and CXCL12 variation and HIV/AIDS in Papua New Guinea

Polymorphisms in chemokine receptors, serving as HIV co-receptors, and their ligands are among the well-known host genetic factors associated with susceptibility to HIV infection and/or disease progression. Papua New Guinea (PNG) has one of the highest adult HIV prevalences in the Asia-Pacific region. However, information regarding the distribution of polymorphisms in chemokine receptor (CCR5, CCR2) and chemokine (CXCL12) genes in PNG is very limited. In this study, we genotyped a total of nine CCR2-CCR5 polymorphisms, including CCR2 190G >A, CCR5 -2459G >A and Δ32, and CXCL12 801G >A in PNG (n=258), North America (n=184), and five countries in West Africa (n=178). Using this data, we determined previously characterized CCR5 haplotypes. In addition, based on the previously reported associations of CCR2 190, CCR5 -2459, CCR5 open reading frame, and CXCL12 801 genotypes with HIV acquisition and/or disease progression, we calculated composite full risk scores, considering both protective as well as susceptibility effects of the CXCL12 801 AA genotype. We observed a very high frequency of the CCR5 -2459A allele (0.98) in the PNG population, which together with the absence of Δ32 resulted in a very high frequency of the HHE haplotype (0.92). These frequencies were significantly higher than in any other population (all P-values<0.001). Regardless of whether we considered the CXCL12 801 AA genotype protective or susceptible, the risk scores were significantly higher in the PNG population compared with any other population (all P-values<0.001). The results of this study provide new insights regarding CCR5 variation in the PNG population, and suggest that the collective variation in CCR2, CCR5, and CXCL12 may increase the risk of HIV/AIDS in a large majority of Papua New Guineans.

Epigenetic mechanisms, T-cell activation, and CCR5 genetics interact to regulate T-cell expression of CCR5, the major HIV-1 coreceptor

T-cell expression levels of CC chemokine receptor 5 (CCR5) are a critical determinant of HIV/AIDS susceptibility, and manifest wide variations (i) between T-cell subsets and among individuals and (ii) in T-cell activation-induced increases in expression levels. We demonstrate that a unifying mechanism for this variation is differences in constitutive and T-cell activation-induced DNA methylation status of CCR5 cis-regulatory regions (cis-regions). Commencing at an evolutionarily conserved CpG (CpG -41), CCR5 cis-regions manifest lower vs. higher methylation in T cells with higher vs. lower CCR5 levels (memory vs. naïve T cells) and in memory T cells with higher vs. lower CCR5 levels. HIV-related and in vitro induced T-cell activation is associated with demethylation of these cis-regions. CCR5 haplotypes associated with increased vs. decreased gene/surface expression levels and HIV/AIDS susceptibility magnify vs. dampen T-cell activation-associated demethylation. Methylation status of CCR5 intron 2 explains a larger proportion of the variation in CCR5 levels than genotype or T-cell activation. The ancestral, protective CCR5-HHA haplotype bears a polymorphism at CpG -41 that is (i) specific to southern Africa, (ii) abrogates binding of the transcription factor CREB1 to this cis-region, and (iii) exhibits a trend for overrepresentation in persons with reduced susceptibility to HIV and disease progression. Genotypes lacking the CCR5-Δ32 mutation but with hypermethylated cis-regions have CCR5 levels similar to genotypes heterozygous for CCR5-Δ32. In HIV-infected individuals, CCR5 cis-regions remain demethylated, despite restoration of CD4+ counts (≥800 cells per mm(3)) with antiretroviral therapy. Thus, methylation content of CCR5 cis-regions is a central epigenetic determinant of T-cell CCR5 levels, and possibly HIV-related outcomes.

Multiply to conquer: Copy number variations at Ppd-B1 and Vrn-A1 facilitate global adaptation in wheat

Background

Copy number variation was found to be a frequent type of DNA polymorphism in the human genome often associated with diseases but its importance in crops and the effects on agronomic traits are still largely unknown.

Results

Here, we employed a large worldwide panel of 1110 winter wheat varieties to assess the frequency and the geographic distribution of copy number variants at the Photoperiod-B1 (Ppd-B1) and the Vernalization-A1 (Vrn-A1) loci as well as their effects on flowering time under field conditions. We identified a novel four copy variant of Vrn-A1 and based on the phylogenetic relationships among the lines show that the higher copy variants at both loci are likely to have arisen independently multiple times. In addition, we found that the frequency of the different copy number variants at both loci reflects the environmental conditions in the varieties’ region of origin and based on multi-location field trials show that Ppd-B1 copy number has a substantial effect on the fine-tuning of flowering time.

Conclusions

In conclusion, our results show the importance of copy number variation at Ppd-B1 and Vrn-A1 for the global adaptation of wheat making it a key factor for wheat success in a broad range of environments and in a wider context substantiate the significant role of copy number variation in crops.

Electronic supplementary material

The online version of this article (doi:10.1186/s12863-015-0258-0) contains supplementary material, which is available to authorized users.

Experimental Evolution Identifies Vaccinia Virus Mutations in A24R and A35R That Antagonize the Protein Kinase R Pathway and Accompany Collapse of an Extragenic Gene Amplification

Most new human infectious diseases emerge from cross-species pathogen transmissions; however, it is not clear how viruses adapt to productively infect new hosts. Host restriction factors represent one species-specific barrier that viruses may initially have little ability to inhibit in new hosts. For example, viral antagonists of protein kinase R (PKR) vary in their ability to block PKR-mediated inhibition of viral replication, in part due to PKR allelic variation between species. We previously reported that amplification of a weak PKR antagonist encoded by rhesus cytomegalovirus, rhtrs1, improved replication of a recombinant poxvirus (VVΔEΔK+RhTRS1) in several resistant primate cells. To test whether amplification increases the opportunity for mutations that improve virus replication with only a single copy of rhtrs1 to evolve, we passaged rhtrs1-amplified viruses in semipermissive primate cells. After passage, we isolated two viruses that contained only a single copy of rhtrs1 yet replicated as well as the amplified virus. Surprisingly, rhtrs1 was not mutated in these viruses; instead, we identified mutations in two vaccinia virus (VACV) genes, A24R and A35R, either of which was sufficient to improve VVΔEΔK+RhTRS1 replication. Neither of these genes has previously been implicated in PKR antagonism. Furthermore, the mutation in A24R, but not A35R, increased resistance to the antipoxviral drug isatin-β-thiosemicarbazone, suggesting that these mutations employ different mechanisms to evade PKR. This study supports our hypothesis that gene amplification may provide a "molecular foothold," broadly improving replication to facilitate rapid adaptation, while subsequent mutations maintain this efficient replication in the new host without requiring gene amplification.

IMPORTANCE

Understanding how viruses adapt to a new host may help identify viruses poised to cross species barriers before an outbreak occurs. Amplification of rhtrs1, a weak viral antagonist of the host antiviral protein PKR, enabled a recombinant vaccinia virus to replicate in resistant cells from humans and other primates. After serial passage of rhtrs1-amplified viruses, there arose in two vaccinia virus genes mutations that improved viral replication without requiring rhtrs1 amplification. Neither of these genes has previously been associated with inhibition of the PKR pathway. These data suggest that gene amplification can improve viral replication in a resistant host species and facilitate the emergence of novel adaptations that maintain the foothold needed for continued replication and spread in the new host.

Detection of HIV-1-specific T-cell immune responses in highly HIV-exposed uninfected individuals by in-vitro dendritic cell co-culture

BACKGROUND

Although virus-specific responses are rarely detected by conventional approaches, we report here the detection of T-cell responses in HIV-exposed seronegative (HESN) patients by two distinct assays.

METHODS

HIV-specific T-cell responses were analyzed by enzyme-linked immunospot in peripheral blood mononuclear cells from HESN patients after a 48-h co-culture with boosted dendritic cells. Additionally, a boosted flow cytometry approach was used to capture antiviral T-cell responses. Host genetic factors and T-cell activation were also analyzed to assess their implication on HIV exposure.

RESULTS

Of the 45 HESN individuals tested, up to 11 (24.4%) showed at least one response to peptide pools covering HIV Gag and Nef. A positive correlation was observed between the intensity (P = 0.0022) and magnitude (P = 0.0174) of the response detected in the HESN, and the viral load of the HIV-positive partner. Moreover, the result from the boosted flow and cytomix analyses showed a dominant Th1-like response pattern against HIV antigens, especially in CD8 T-cell populations.

CONCLUSIONS

The combined use of our boosted dendritic cell technique with a boosted flow cytometric approach allows us both to detect specific HIV-positive responses in a higher percentage of HESN patients and to define specific effector function profiles. This study contributes to a better understanding of resistance to HIV infection.

Complex and multi-allelic copy number variation in human disease

Hundreds of copy number variants are complex and multi-allelic, in that they have many structural alleles and have rearranged multiple times in the ancestors who contributed chromosomes to current humans. Not only are the relationships of these multi-allelic CNVs (mCNVs) to phenotypes generally unknown, but many mCNVs have not yet been described at the basic levels—alleles, allele frequencies, structural features—that support genetic investigation. To date, most reported disease associations to these variants have been ascertained through candidate gene studies. However, only a few associations have reached the level of acceptance defined by durable replications in many cohorts. This likely stems from longstanding challenges in making precise molecular measurements of the alleles individuals have at these loci. However, approaches for mCNV analysis are improving quickly, and some of the unique characteristics of mCNVs may assist future association studies. Their various structural alleles are likely to have different magnitudes of effect, creating a natural allelic series of growing phenotypic impact and giving investigators a set of natural predictions and testable hypotheses about the extent to which each allele of an mCNV predisposes to a phenotype. Also, mCNVs’ low-to-modest correlation to individual single-nucleotide polymorphisms (SNPs) may make it easier to distinguish between mCNVs and nearby SNPs as the drivers of an association signal, and perhaps, make it possible to preliminarily screen candidate loci, or the entire genome, for the many mCNV–disease relationships that remain to be discovered.

Divergence patterns of genic copy number variation in natural populations of the house mouse (Mus musculus domesticus) reveal three conserved genes with major population-specific expansions

Copy number variation represents a major source of genetic divergence, yet the evolutionary dynamics of genic copy number variation in natural populations during differentiation and adaptation remain unclear. We applied a read depth approach to genome resequencing data to detect copy number variants (CNVs) ≥1 kb in wild-caught mice belonging to four populations of Mus musculus domesticus. We complemented the bioinformatics analyses with experimental validation using droplet digital PCR. The specific focus of our analysis is CNVs that include complete genes, as these CNVs could be expected to contribute most directly to evolutionary divergence. In total, 1863 transcription units appear to be completely encompassed within CNVs in at least one individual when compared to the reference assembly. Further, 179 of these CNVs show population-specific copy number differences, and 325 are subject to complete deletion in multiple individuals. Among the most copy-number variable genes are three highly conserved genes that encode the splicing factor CWC22, the spindle protein SFI1, and the Holliday junction recognition protein HJURP. These genes exhibit population-specific expansion patterns that suggest involvement in local adaptations. We found that genes that overlap with large segmental duplications are generally more copy-number variable. These genes encode proteins that are relevant for environmental and behavioral interactions, such as vomeronasal and olfactory receptors, as well as major urinary proteins and several proteins of unknown function. The overall analysis shows that genic CNVs contribute more to population differentiation in mice than in humans and may promote and speed up population divergence.

Adaptive potential of genomic structural variation in human and mammalian evolution

Because phenotypic innovations must be genetically heritable for biological evolution to proceed, it is natural to consider new mutation events as well as standing genetic variation as sources for their birth. Previous research has identified a number of single-nucleotide polymorphisms that underlie a subset of adaptive traits in organisms. However, another well-known class of variation, genomic structural variation, could have even greater potential to produce adaptive phenotypes, due to the variety of possible types of alterations (deletions, insertions, duplications, among others) at different genomic positions and with variable lengths. It is from these dramatic genomic alterations, and selection on their phenotypic consequences, that adaptations leading to biological diversification could be derived. In this review, using studies in humans and other mammals, we highlight examples of how phenotypic variation from structural variants might become adaptive in populations and potentially enable biological diversification. Phenotypic change arising from structural variants will be described according to their immediate effect on organismal metabolic processes, immunological response and physical features. Study of population dynamics of segregating structural variation can therefore provide a window into understanding current and historical biological diversification.

Function of cancer associated genes revealed by modern univariate and multivariate association tests

Copy number variation (CNV) plays a role in pathogenesis of many human diseases, especially cancer. Several whole genome CNV association studies have been performed for the purpose of identifying cancer associated CNVs. Here we undertook a novel approach to whole genome CNV analysis, with the goal being identification of associations between CNV of different genes (CNV-CNV) across 60 human cancer cell lines. We hypothesize that these associations point to the roles of the associated genes in cancer, and can be indicators of their position in gene networks of cancer-driving processes. Recent studies show that gene associations are often non-linear and non-monotone. In order to obtain a more complete picture of all CNV associations, we performed omnibus univariate analysis by utilizing dCov, MIC, and HHG association tests, which are capable of detecting any type of association, including non-monotone relationships. For comparison we used Spearman and Pearson association tests, which detect only linear or monotone relationships. Application of dCov, MIC and HHG tests resulted in identification of twice as many associations compared to those found by Spearman and Pearson alone. Interestingly, most of the new associations were detected by the HHG test. Next, we utilized dCov's and HHG's ability to perform multivariate analysis. We tested for association between genes of unknown function and known cancer-related pathways. Our results indicate that multivariate analysis is much more effective than univariate analysis for the purpose of ascribing biological roles to genes of unknown function. We conclude that a combination of multivariate and univariate omnibus association tests can reveal significant information about gene networks of disease-driving processes. These methods can be applied to any large gene or pathway dataset, allowing more comprehensive analysis of biological processes.

Clinical implications of copy number variations in autoimmune disorders

Human genetic variation is represented by the genetic differences both within and among populations, and most genetic variants do not cause overt diseases but contribute to disease susceptibility and influence drug response. During the last century, various genetic variants, such as copy number variations (CNVs), have been associated with diverse human disorders. Here, we review studies on the associations between CNVs and autoimmune diseases to gain some insight. First, some CNV loci are commonly implicated in various autoimmune diseases, such as Fcγ receptors in patients with systemic lupus erythemoatosus or idiopathic thrombocytopenic purpura and β-defensin genes in patients with psoriasis or Crohn's disease. This means that when a CNV locus is associated with a particular autoimmune disease, we should examine its potential associations with other diseases. Second, interpopulation or interethnic differences in the effects of CNVs on phenotypes exist, including disease susceptibility, and evidence suggests that CNVs are important to understand susceptibility to and pathogenesis of autoimmune diseases. However, many findings need to be replicated in independent populations and different ethnic groups. The validity and reliability of detecting CNVs will improve quickly as genotyping technology advances, which will support the required replication.

The impact of human copy number variation on gene expression

Recent years have witnessed a flurry of important technological and methodological developments in the discovery and analysis of copy number variations (CNVs), which are increasingly enabling the systematic evaluation of their impact on a broad range of phenotypes from molecular-level (intermediate) traits to higher-order clinical phenotypes. Like single nucleotide variants in the human genome, CNVs have been linked to complex traits in humans, including disease and drug response. These recent developments underscore the importance of incorporating complex forms of genetic variation into disease mapping studies and promise to transform our understanding of genome function and the genetic basis of disease. Here we review some of the findings that have emerged from transcriptome studies of CNVs facilitated by the rapid advances in -omics technologies and corresponding methodologies.

A decade of structural variants: description, history and methods to detect structural variation

In the past decade, the view on genomic structural variation (SV) has been changed completely. SVs, previously considered rare events, are now recognized as the largest source of interindividual genetic variation affecting more bases than single nucleotide polymorphisms, variable number of tandem repeats and other small genetic variants. They have also been shown to play a role in phenotypic variation and in disease. In this review, the authors will provide an introduction to SV; a short historical perspective on the research of this source of genomic variation; a description of the types of structural variants, and on how they may have arisen; and an overview on methods of detecting structural variants, focusing on the analysis of high-throughput sequencing data.

A method for generating new datasets based on copy number for cancer analysis

New data sources for the analysis of cancer data are rapidly supplementing the large number of gene-expression markers used for current methods of analysis. Significant among these new sources are copy number variation (CNV) datasets, which typically enumerate several hundred thousand CNVs distributed throughout the genome. Several useful algorithms allow systems-level analyses of such datasets. However, these rich data sources have not yet been analyzed as deeply as gene-expression data. To address this issue, the extensive toolsets used for analyzing expression data in cancerous and noncancerous tissue (e.g., gene set enrichment analysis and phenotype prediction) could be redirected to extract a great deal of predictive information from CNV data, in particular those derived from cancers. Here we present a software package capable of preprocessing standard Agilent copy number datasets into a form to which essentially all expression analysis tools can be applied. We illustrate the use of this toolset in predicting the survival time of patients with ovarian cancer or glioblastoma multiforme and also provide an analysis of gene- and pathway-level deletions in these two types of cancer.

Haplotype phasing and inheritance of copy number variants in nuclear families

DNA copy number variants (CNVs) that alter the copy number of a particular DNA segment in the genome play an important role in human phenotypic variability and disease susceptibility. A number of CNVs overlapping with genes have been shown to confer risk to a variety of human diseases thus highlighting the relevance of addressing the variability of CNVs at a higher resolution. So far, it has not been possible to deterministically infer the allelic composition of different haplotypes present within the CNV regions. We have developed a novel computational method, called PiCNV, which enables to resolve the haplotype sequence composition within CNV regions in nuclear families based on SNP genotyping microarray data. The algorithm allows to i) phase normal and CNV-carrying haplotypes in the copy number variable regions, ii) resolve the allelic copies of rearranged DNA sequence within the haplotypes and iii) infer the heritability of identified haplotypes in trios or larger nuclear families. To our knowledge this is the first program available that can deterministically phase null, mono-, di-, tri- and tetraploid genotypes in CNV loci. We applied our method to study the composition and inheritance of haplotypes in CNV regions of 30 HapMap Yoruban trios and 34 Estonian families. For 93.6% of the CNV loci, PiCNV enabled to unambiguously phase normal and CNV-carrying haplotypes and follow their transmission in the corresponding families. Furthermore, allelic composition analysis identified the co-occurrence of alternative allelic copies within 66.7% of haplotypes carrying copy number gains. We also observed less frequent transmission of CNV-carrying haplotypes from parents to children compared to normal haplotypes and identified an emergence of several de novo deletions and duplications in the offspring.

Deletion of REXO1L1 locus in a patient with malabsorption syndrome, growth retardation, and dysmorphic features: a novel recognizable microdeletion syndrome?

Background

Copy number variations (CNVs) can contribute to genetic variation among individuals and/or have a significant influence in causing diseases. Many studies consider new CNVs’ effects on protein family evolution giving rise to gene duplicates or losses. “Unsuccessful” duplicates that remain in the genome as pseudogenes often exhibit functional roles. So, changes in gene and pseudogene number may contribute to development or act as susceptibility alleles of diseases.

Case presentation

We report a de novo heterozygous 271 Kb microdeletion at 8q21.2 region which includes the family of REXO1L genes and pseudogenes in a young man affected by global development delay, progeroid signs, and gastrointestinal anomalies. Molecular and cellular analysis showed that the REXO1L1 gene hemizygosity in a patient’s fibroblasts induces genetic instability and increased apoptosis after treatment with different DNA damage-induced agents.

Conclusions

The present results support the hypothesis that low copy gene number within REXO1L1 cluster could play a significant role in this complex clinical and cellular phenotype.

Effect of MPG gene rs2858056 polymorphism, copy number variation, and level of serum MPG protein on the risk for rheumatoid arthritis

OBJECTIVE

This study examined the role of SNP rs2858056 of the MPG gene on the incidence and severity of rheumatoid arthritis (RA).

METHODS

This cohort study enrolled 365 RA patients and 375 age- and gender-matched healthy controls, all of whom had Han Chinese ethnicity and were from Taiwan. Gene polymorphism of the SNP rs2858056 of MPG was determined from genomic DNA. Allelic frequencies and genotypes were compared among cases and controls. Quantitation of rs2858056 copy number variation (CNV) was determined. Serum samples from RA patients and controls were analyzed to determine serum levels of MPG. The relationship between rs2858056 polymorphism and clinical manifestations of RA was evaluated.

RESULTS

Our results indicated a statistically significant difference in genotype frequency distributions at rs2858056 for RA patients and controls (p = 0.05) and a significant difference in allelic frequency in patients and controls (p = 0.04). Furthermore, there was a significantly greater level of serum MPG protein in patients than controls (p < 0.001). However, the cases and controls had no significant differences in MPG CNV (p = 0.12). We also did not detect any association of the MPG rs2858056 with rheumatoid factor (RF), extraarticular involvement, or bone erosion in the RA patients.

CONCLUSION

Our study suggests that RA is associated with a polymorphism in the MPG gene (rs2858056) and increased serum level of the MPG protein.

p21(WAF1/CIP1) RNA expression in highly HIV-1 exposed, uninfected individuals

Some individuals remain HIV-1 antibody and PCR negative after repeated exposures to the virus, and are referred to as HIV-exposed seronegatives (HESN). However, the causes of resistance to HIV-1 infection in cases other than those with a homozygous CCR5Δ32 deletion are unclear. We hypothesized that human p21^WAF1/CIP1 (a cyclin-dependent kinase inhibitor) could play a role in resistance to HIV-1 infection in HESN, as p21 expression has been associated with suppression of HIV-1 in elite controllers and reported to block HIV-1 integration in cell culture. We measured p21 RNA expression in PBMC from 40 HESN and 40 low exposure HIV-1 seroconverters (LESC) prior to their infection using a real-time PCR assay. Comparing the 20 HESN with the highest exposure risk (median = 111 partners/2.5 years prior to the 20 LESC with the lowest exposure risk (median = 1 partner/2.5 years prior), p21 expression trended higher in HESN in only one of two experiments (P = 0.11 vs. P = 0.80). Additionally, comparison of p21 expression in the top 40 HESN (median = 73 partners/year) and lowest 40 LESC (median = 2 partners/year) showed no difference between the groups (P = 0.84). There was a weak linear trend between risk of infection after exposure and increasing p21 gene expression (R² = 0.02, P = 0.12), but again only in one experiment. Hence, if p21 expression contributes to the resistance to viral infection in HESN, it likely plays a minor role evident only in those with extremely high levels of exposure to HIV-1.

Increased Levels of Macrophage Inflammatory Proteins Result in Resistance to R5-Tropic HIV-1 in a Subset of Elite Controllers

Elite controllers (ECs) are a rare group of HIV seropositive individuals who are able to control viral replication without antiretroviral therapy. The mechanisms responsible for this phenotype, however, have not been fully elucidated. In this study, we examined CD4(+) T cell resistance to HIV in a cohort of elite controllers and explored transcriptional signatures associated with cellular resistance. We demonstrate that a subgroup of elite controllers possess CD4(+) T cells that are specifically resistant to R5-tropic HIV while remaining fully susceptible to X4-tropic and vesicular stomatitis virus G (VSV-G)-pseudotyped viruses. Transcriptome analysis revealed 17 genes that were differentially regulated in resistant elite controllers relative to healthy controls. Notably, the genes encoding macrophage inflammatory protein 1α (MIP-1α), CCL3 and CCL3L1, were found to be upregulated. The MIP-1α, MIP-1β, and RANTES chemokines are natural ligands of CCR5 and are known to interfere with HIV replication. For three elite controllers, we observed increased production of MIP-1α and/or MIP-1β at the protein level. The supernatant from resistant EC cells contained MIP-1α and MIP-1β and was sufficient to confer R5-tropic resistance to susceptible CD4(+) T cells. Additionally, this effect was reversed by using inhibitory anti-MIP antibodies. These results suggest that the T cells of these particular elite controllers may be naturally resistant to HIV infection by blocking R5-tropic viral entry.

IMPORTANCE

HIV is a pandemic health problem, and the majority of seropositive individuals will eventually progress to AIDS unless antiretroviral therapy (ART) is administered. However, rare patients, termed elite controllers, have a natural ability to control HIV infection in the absence of ART, but the mechanisms by which they achieve this phenotype have not been fully explored. This paper identifies one mechanism that may contribute to this natural resistance: some elite controllers have CD4(+) T cells that produce high levels of MIP chemokines, which block R5-tropic HIV entry. This mechanism could potentially be exploited to achieve a therapeutic effect in other HIV-seropositive individuals.

Multiple Break-Points Detection in Array CGH Data via the Cross-Entropy Method

Array comparative genome hybridization (aCGH) is a widely used methodology to detect copy number variations of a genome in high resolution. Knowing the number of break-points and their corresponding locations in genomic sequences serves different biological needs. Primarily, it helps to identify disease-causing genes that have functional importance in characterizing genome wide diseases. For human autosomes the normal copy number is two, whereas at the sites of oncogenes it increases (gain of DNA) and at the tumour suppressor genes it decreases (loss of DNA). The majority of the current detection methods are deterministic in their set-up and use dynamic programming or different smoothing techniques to obtain the estimates of copy number variations. These approaches limit the search space of the problem due to different assumptions considered in the methods and do not represent the true nature of the uncertainty associated with the unknown break-points in genomic sequences. We propose the Cross-Entropy method, which is a model-based stochastic optimization technique as an exact search method, to estimate both the number and locations of the break-points in aCGH data. We model the continuous scale log-ratio data obtained by the aCGH technique as a multiple break-point problem. The proposed methodology is compared with well established publicly available methods using both artificially generated data and real data. Results show that the proposed procedure is an effective way of estimating number and especially the locations of break-points with high level of precision. Availability: The methods described in this article are implemented in the new R package breakpoint and it is available from the Comprehensive R Archive Network at http://CRAN.R-project.org/package=breakpoint.

Genome-wide association study of copy number variations (CNVs) with opioid dependence

Single-nucleotide polymorphisms that have been associated with opioid dependence (OD) altogether account for only a small proportion of the known heritability. Most of the genetic risk factors are unknown. Some of the 'missing heritability' might be explained by copy number variations (CNVs) in the human genome. We used Illumina HumanOmni1 arrays to genotype 5152 African-American and European-American OD cases and screened controls and implemented combined CNV calling methods. After quality control measures were applied, a genome-wide association study (GWAS) of CNVs with OD was performed. For common CNVs, two deletions and one duplication were significantly associated with OD genome-wide (eg, P=2 × 10(-8) and OR (95% CI)=0.64 (0.54-0.74) for a chromosome 18q12.3 deletion). Several rare or unique CNVs showed suggestive or marginal significance with large effect sizes. This study is the first GWAS of OD using CNVs. Some identified CNVs harbor genes newly identified here to be of biological importance in addiction, whereas others affect genes previously known to contribute to substance dependence risk. Our findings augment our specific knowledge of the importance of genomic variation in addictive disorders, and provide an addiction CNV pool for further research. These findings require replication.

The use of a two-tiered testing strategy for the simultaneous detection of small EGFR mutations and EGFR amplification in lung cancer

Lung cancer is the leading cause of cancer-related death worldwide. Recent progress in lung cancer diagnosis and treatment has been achieved due to a better understanding the molecular mechanisms of the disease and the identification of biomarkers that allow more specific cancer treatments. One of the best known examples of personalized therapy is the use of tyrosine kinase inhibitors, such as gefitinib and erlotinib, for the successful treatment of non-small-cell lung cancer patients selected based on the specific EGFR mutations. Therefore, the reliable detection of mutations is critical for the application of appropriate therapy. In this study, we tested a two-tiered mutation detection strategy using real-time PCR assays as a well-validated high-sensitivity method and multiplex ligation-dependent probe amplification (MLPA)-based EGFRmut+ assay as a second-tier standard-sensitivity method. One additional advantage of the applied MLPA method is that it allows the simultaneous detection of EGFR mutations and copy-number alterations (i.e., amplifications) in EGFR, MET and ERBB2. Our analysis showed high concordance between these two methods. With the use of this two-tier strategy, we reliably determined the frequency of EGFR mutations and EGFR, MET and ERBB2 amplifications in over 200 lung cancer samples. Additionally, taking advantage of simultaneous copy number and small mutation analyses, we showed a very strong correlation between EGFR mutations and EGFR amplifications and a mutual exclusiveness of EGFR mutations/amplifications with MET and ERBB2 amplifications. Our results proved the reliability and usefulness of the two-tiered EGFR testing strategy.

A comparison of assays for accurate copy number measurement of the low-affinity Fc gamma receptor genes FCGR3A and FCGR3B

The FCGR3 locus encoding the low affinity activating receptor FcγRIII, plays a vital role in immunity triggered by cellular effector and regulatory functions. Copy number of the genes FCGR3A and FCGR3B has previously been reported to affect susceptibility to several autoimmune diseases and chronic inflammatory conditions. However, such genetic association studies often yield inconsistent results; hence require assays that are robust with low error rate. We investigated the accuracy and efficiency in estimating FCGR3 CNV by comparing Sequenom MassARRAY and paralogue ratio test-restriction enzyme digest variant ratio (PRT-REDVR). In addition, since many genetic association studies of FCGR3B CNV were carried out using real-time quantitative PCR, we have also included the evaluation of that method’s performance in estimating the multi-allelic CNV of FCGR3B. The qPCR assay exhibited a considerably broader distribution of signal intensity, potentially introducing error in estimation of copy number and higher false positive rates. Both Sequenom and PRT-REDVR showed lesser systematic bias, but Sequenom skewed towards copy number normal (CN = 2). The discrepancy between Sequenom and PRT-REDVR might be attributed either to batch effects noise in individual measurements. Our study suggests that PRT-REDVR is more robust and accurate in genotyping the CNV of FCGR3, but highlights the needs of multiple independent assays for extensive validation when performing a genetic association study with multi-allelic CNVs.