Glycophorin A and glycophorin B are red blood cell surface proteins and are both receptors for the parasite Plasmodium falciparum, which is the principal cause of malaria in sub-Saharan Africa. DUP4 is a complex structural genomic variant that carries extra copies of a glycophorin A-glycophorin B fusion gene and has a dramatic effect on malaria risk by reducing the risk of severe malaria by up to 40%. Using fiber-FISH and Illumina sequencing, we validate the structural arrangement of the glycophorin locus in the DUP4 variant and reveal somatic variation in copy number of the glycophorin B-glycophorin A fusion gene. By developing a simple, specific, PCR-based assay for DUP4, we show that the DUP4 variant reaches a frequency of 13% in the population of a malaria-endemic village in south-eastern Tanzania. We genotype a substantial proportion of that village and demonstrate an association of DUP4 genotype with hemoglobin levels, a phenotype related to malaria, using a family-based association test. Taken together, we show that DUP4 is a complex structural variant that may be susceptible to somatic variation and show that DUP4 is associated with a malarial-related phenotype in a longitudinally followed population.

Background: Salivary amylase in humans is encoded by the copy variable gene AMY1 in the amylase gene cluster on chromosome 1. Although the role of salivary amylase is well established, the consequences of the copy number variation (CNV) at AMY1 on salivary amylase protein production are less well understood. The amylase gene cluster is highly structured with a fundamental difference between odd and even AMY1 copy number haplotypes. In this study, we aimed to explore, in samples from 119 unrelated individuals, not only the effects of AMY1 CNV on salivary amylase protein expression and amylase enzyme activity but also whether there is any evidence for underlying difference between the common haplotypes containing odd numbers of AMY1 and even copy number haplotypes.Results: AMY1 copy number was significantly correlated with the variation observed in salivary amylase production (11.7% of variance, P < 0.0005) and enzyme activity (13.6% of variance, P < 0.0005) but did not explain the majority of observed variation between individuals. AMY1-odd and AMY1-even haplotypes showed a different relationship between copy number and expression levels, but the difference was not statistically significant (P = 0.052).Conclusions: Production of salivary amylase is correlated with AMY1 CNV, but the majority of interindividual variation comes from other sources. Long-range haplotype structure may affect expression, but this was not significant in our data.

The human salivary amylase genes display extensive copy number variation (CNV), and recent work has implicated this variation in adaptation to starch-rich diets, and in association with body mass index. In this work, we use paralogue ratio tests, microsatellite analysis, read depth and fibre-FISH to demonstrate that human amylase CNV is not a smooth continuum, but is instead partitioned into distinct haplotype classes. There is a fundamental structural distinction between haplotypes containing odd or even numbers of AMY1 gene units, in turn coupled to CNV in pancreatic amylase genes AMY2A and AMY2B. Most haplotypes have one copy each of AMY2A and AMY2B and contain an odd number of copies of AMY1; consequently, most individuals have an even total number of AMY1. In contrast, haplotypes carrying an even number of AMY1 genes have rearrangements leading to CNVs of AMY2A/AMY2B. Read-depth and experimental data show that different populations harbour different proportions of these basic haplotype classes. In Europeans, the copy numbers of AMY1 and AMY2A are correlated, so that phenotypic associations caused by variation in pancreatic amylase copy number could be detected indirectly as weak association with AMY1 copy number. We show that the quantitative polymerase chain reaction (qPCR) assay previously applied to the high-throughput measurement of AMY1 copy number is less accurate than the measures we use and that qPCR data in other studies have been further compromised by systematic miscalibration. Our results uncover new patterns in human amylase variation and imply a potential role for AMY2 CNV in functional associations.

Tuberculosis is a major infectious disease and functional studies have provided evidence that both the chemokine MIP-1α and its receptor CCR5 play a role in susceptibility to TB. Thus by measuring copy number variation of CCL3L1, one of the genes that encode MIP-1α, and genotyping a functional promoter polymorphism -2459A > G in CCR5 (rs1799987) we investigate the influence of MIP-1α and CCR5, independently and combined, in susceptibility to clinically active TB in three populations, a Peruvian population (n = 1132), a !Xhosa population (n = 605) and a South African Coloured population (n = 221). The three populations include patients with clinically diagnosed pulmonary TB, as well as other, less prevalent forms of extrapulmonary TB.

BACKGROUND:Multi-allelic copy number variants include examples of extensive variation between individuals in the copy number of important genes, most notably genes involved in immune function. The definition of this variation, and analysis of its impact on function, has been hampered by the technical difficulty of large-scale but accurate typing of genomic copy number. The copy-variable alpha-defensin locus DEFA1A3 on human chromosome 8 commonly varies between 4 and 10 copies per diploid genome, and presents considerable challenges for accurate high-throughput typing.RESULTS:In this study, we developed two paralogue ratio tests and three allelic ratio measurements that, in combination, provide an accurate and scalable method for measurement of DEFA1A3 gene number. We combined information from different measurements in a maximum-likelihood framework which suggests that most samples can be assigned to an integer copy number with high confidence, and applied it to typing 589 unrelated European DNA samples. Typing the members of three-generation pedigrees provided further reassurance that correct integer copy numbers had been assigned. Our results have allowed us to discover that the SNP rs4300027 is strongly associated with DEFA1A3 gene copy number in European samples.CONCLUSIONS:We have developed an accurate and robust method for measurement of DEFA1A3 copy number. Interrogation of rs4300027 and associated SNPs in Genome-Wide Association Study SNP data provides no evidence that alpha-defensin copy number is a strong risk factor for phenotypes such as Crohn's disease, type I diabetes, HIV progression and multiple sclerosis.

Copy number variation (CNV) is becoming increasingly important as a feature of human variation in disease susceptibility studies. However, the consequences of CNV are not so well understood. Here, we present data exploring the functional consequences of CNV of CCL3L1 in 55 independent UK samples with no known clinical phenotypes. The copy number of CCL3L1 was determined by the paralogue ratio test, and expression levels of macrophage inflammatory protein-1α (MIP-1α) and mRNA from stimulated monocytes were measured and analysed. The data show no statistically significant association of MIP-1α protein levels with copy number. However, there was a significant correlation between copy number and CCL3L1:CCL3 mRNA ratio. The data also provide evidence that expression of CCL3 predominates in both protein and mRNA, and therefore the observed variation of CCL3 is potentially more important biologically than that of CNV of CCL3L1.

Genotyping of multiallelic copy number variants (CNVs) is technically difficult and can lead to inaccurate conclusions. This is reflected by inconsistent results published for the CNV C-C chemokine ligand 3-like 1 (CCL3L1) and its contribution to rheumatoid arthritis (RA) susceptibility. In order to draw robust conclusions about CCL3L1 involvement in RA, we have performed association analysis (CNVtools) using genotyping by the paralogue ratio test of a Norwegian RA case–control material (N=1877). We also analyzed the associations after stratification for anti-citrullinated peptide antibody (ACPA) status. Clear clusters representing specific copy number classes were evident, but significant differential bias was observed resulting in a systematic trend toward slightly higher apparent copy number for cases relative to controls. Controlling for bias revealed no significant differences in copy number distribution either between all patients and controls, or after ACPA stratification. Our results do not support involvement of the CCL3L1 CNV in RA susceptibility.

Copy number variation can contribute to the variation observed in susceptibility to complex diseases. Here we present the first study to investigate copy number variation of the chemokine gene CCL3L1 with susceptibility to malaria. We present a family-based genetic analysis of a Tanzanian population (n = 922), using parasite load, mean number of clinical infections of malaria and haemoglobin levels as phenotypes. Copy number of CCL3L1 was measured using the paralogue ratio test (PRT) and the dataset exhibited copy numbers ranging between 1 and 10 copies per diploid genome (pdg). Association between copy number and phenotypes was assessed. Furthermore, we were able to identify copy number haplotypes in some families, using microsatellites within the copy variable region, for transmission disequilibrium testing. We identified a high level of copy number haplotype diversity and find some evidence for an association of low CCL3L1 copy number with protection from anaemia.

Background Copy number variation (CNV) contributes to the variation observed between individuals and can influence human disease progression, but the accurate measurement of individual copy numbers is technically challenging. In the work presented here we describe a modification to a previously described paralogue ratio test (PRT) method for genotyping the CCL3L1/CCL4L1 copy variable region, which we use to ascertain CCL3L1/CCL4L1 copy number in 1581 European samples. As the products of CCL3L1 and CCL4L1 potentially play a role in autoimmunity we performed case control association studies with Crohn's disease, rheumatoid arthritis and psoriasis clinical cohorts. Results We evaluate the PRT methodology used, paying particular attention to accuracy and precision, and highlight the problems of differential bias in copy number measurements. Our PRT methods for measuring copy number were of sufficient precision to detect very slight but systematic differential bias between results from case and control DNA samples in one study. We find no evidence for an association between CCL3L1 copy number and Crohn's disease, rheumatoid arthritis or psoriasis. Conclusions Differential bias of this small magnitude, but applied systematically across large numbers of samples, would create a serious risk of false positive associations in copy number, if measured using methods of lower precision, or methods relying on single uncorroborated measurements. In this study the small differential bias detected by PRT in one sample set was resolved by a simple pre-treatment by restriction enzyme digestion.

CCL3LI and CCL4LI are chemokine genes, located on chromosome 17q12. They are copy number variable genes which share 95% sequence identity with their non-copy number variable paralogues CCL3 and CCL4. The copy number of these genes varies between populations and has been reported to be associated with phenotypes such as susceptibility to HIV infection, hepatitis C virus infection, Kawasaki disease and SLE. The aim of this study is to understand the evolutionary history of variation at the CCL3L1/CCL4LI cluster. To accomplish this goal, several approaches including typing microsatellites, single nucleotide polymorphisms (SNPs) and CCL3L 1/CCL4L 1 sequence haplotypes were used to investigate the association with CCL3L 1 and CCL4L 1 copy number. However, the results showed that there is no strong association between a single-copy marker and CCL3L 1 and CCL4LI copy number, but there is evidence of recombination. Therefore, this may suggest that CCL3L 1/CCL4L 1 is a complex region and one plausible hypothesis is that there is a high rate of recombination in this region. This study of the evolution of CCL3L 1/CCL4L 1 haplotypes showed that a major one-copy CCL3L 1/CCL4L I haplotype (about 70% haplotype frequency) identified in humans, represents the ancestral state, as inferred from comparison with chimpanzee.

Previous studies have established a genetic component for susceptibility to malaria. Here we use a pedigree based approach, and transmission disequilibrium testing (TDT), to identify immune response genes that influence susceptibility to Plasmodium falciparum malarial phenotypes (parasite density and frequency of clinical episodes) in a Tanzanian population. Evidence for association was observed between markers in the TNF gene cluster and both the malarial phenotypes. There was weaker evidence for associations between HLA-DRB1*04, HLA-DRB1*10, and loci in the TCRBV region with parasite density. There was no evidence for association with polymorphisms in the IL10 promoter, IL1 gene cluster, or from the IL4/IL13 region.

Malignant hyperthermia (MH) is an inherited pharmacogenetic disorder of skeletal muscle, characterised by an elevated calcium release from the skeletal muscle sarcoplasmic reticulum. The dihydropyridine receptor (DHPR) plays an essential role in excitation-contraction coupling and calcium homeostasis in skeletal muscle. This study focuses on the gene CACNA1S which encodes the α1 subunit of the DHPR, in order to establish whether CACNA1S plays a major role in MH susceptibility in the UK.

In this study we present 3 families with malignant hyperthermia (MH), all of Indian subcontinent descent. One individual from each of these families was fully sequenced for RYR1 and presented with the non‐synonymous change c.11315G>A/p.R3772Q. When present in the homozygous state c.11315*A is associated with myopathic symptoms. Muscle Nerve, 2009

This study builds upon the established genetic control of antimalarial immune responses and prior association studies by using a family-based approach, transmission disequilibrium testing, to identify immune response genes that influence antibody responses to Plasmodium falciparum infection in an endemic Tanzanian population. Candidate polymorphisms are within the interleukin-1 (IL-1) gene cluster, the IL-10 promoter, Major histocompatibility complex class II and III, the 5q31-q33 region, and the T-Cell Receptor beta variable region. There was a significant association between the IL1RN alleles and total IgE. Weak evidence for association was present between polymorphisms in the IL10 promoter region and both anti-P falciparum IgE and IgG4 antibodies.

The RYR1 gene encodes the skeletal muscle isoform ryanodine receptor and is fundamental to the process of excitation-contraction coupling and skeletal muscle calcium homeostasis. Mapping to chromosome 19q13.2, the gene comprises 106 exons and encodes a protein of 5,038 amino acids. Mutations in the gene have been found in association with several diseases: the pharmacogenetic disorder, malignant hyperthermia (MH); and three congenital myopathies, including central core disease (CCD), multiminicore disease (MmD), and in an isolated case of a congenital myopathy characterized on histology by cores and rods. The majority of gene mutations reported are missense changes identified in cases of MH and CCD. In vitro analysis has confirmed that alteration of normal calcium homeostasis is a functional consequence of some of these changes. Genotype-phenotype correlation studies performed using data from MH and CCD patients have also suggested that mutations may be associated with a range of disease severity phenotypes. This review aims to summarize the current understanding of RYR1 mutations reported in association with MH and CCD and the present viewpoint on the use of mutation data to aid clinical diagnosis of these conditions.

This study investigated polymorphisms of genes in two regions of the T-cell antigen receptor β-subunit (TCRB) locus, including BV9S2P, and BV6S7 in a 5′ linkage group, and BV8S3, BV24S1, BV25S1, BV18S1, BV2S1, BV15S1 and BV3S1 in a 3′ linkage group. These loci have been genotyped in individuals from five regions in Africa, including The Gambia, Nigeria, Cameroon, Tanzania, and Zambia, and in individuals from northern Britain, northern India, and Papua New Guinea (PNG). In the 3′ linkage group, 11 unique haplotypes were identified in the combined African populations; two equally frequent haplotypes represent the majority of African chromosomes. One haplotype was found in all four regions studied. This is the most frequent haplotype in the northern British, northern Indian and PNG populations. Although present, it is infrequent in the African populations. A North-South gradient in the frequency of a common African haplotype was observed. The distribution did not represent that of a known disease. Evidence suggests that malaria is not responsible for selection of these haplotypes. Overall, this study highlights large differences in the genetic constitution of the TCRB locus between Africans and other populations.