Haptoglobin (HP) and Haptoglobin-related protein (HPR) copy number variation, natural selection, and trypanosomiasis

Long-read sequencing and structural variant characterization in 1,019 samples from the 1000 Genomes Project

Structural variants (SVs) contribute significantly to human genetic diversity and disease 1-4 . Previously, SVs have remained incompletely resolved by population genomics, with short-read sequencing facing limitations in capturing the whole spectrum of SVs at nucleotide resolution 5-7 . Here we leveraged nanopore sequencing 8 to construct an intermediate coverage resource of 1,019 long-read genomes sampled within 26 human populations from the 1000 Genomes Project. By integrating linear and graph-based approaches for SV analysis via pangenome graph-augmentation, we uncover 167,291 sequence-resolved SVs in these samples, considerably advancing SV characterization compared to population-wide short-read sequencing studies 3,4 . Our analysis details diverse SV classes-deletions, duplications, insertions, and inversions-at population-scale. LINE-1 and SVA retrotransposition activities frequently mediate transductions 9,10 of unique sequences, with both mobile element classes transducing sequences at either the 3'- or 5'-end, depending on the source element locus. Furthermore, analyses of SV breakpoint junctions suggest a continuum of homology-mediated rearrangement processes are integral to SV formation, and highlight evidence for SV recurrence involving repeat sequences. Our open-access dataset underscores the transformative impact of long-read sequencing in advancing the characterisation of polymorphic genomic architectures, and provides a resource for guiding variant prioritisation in future long-read sequencing-based disease studies.

Genomic structural variation: A complex but important driver of human evolution

Structural variants (SVs)-including duplications, deletions, and inversions of DNA-can have significant genomic and functional impacts but are technically difficult to identify and assay compared with single-nucleotide variants. With the aid of new genomic technologies, it has become clear that SVs account for significant differences across and within species. This phenomenon is particularly well-documented for humans and other primates due to the wealth of sequence data available. In great apes, SVs affect a larger number of nucleotides than single-nucleotide variants, with many identified SVs exhibiting population and species specificity. In this review, we highlight the importance of SVs in human evolution by (1) how they have shaped great ape genomes resulting in sensitized regions associated with traits and diseases, (2) their impact on gene functions and regulation, which subsequently has played a role in natural selection, and (3) the role of gene duplications in human brain evolution. We further discuss how to incorporate SVs in research, including the strengths and limitations of various genomic approaches. Finally, we propose future considerations in integrating existing data and biospecimens with the ever-expanding SV compendium propelled by biotechnology advancements.

The Legacy of Infectious Disease Exposure on the Genomic Diversity of Indigenous Southern Mexicans

To characterize host risk factors for infectious disease in Mesoamerican populations, we interrogated 857,481 SNPs assayed using the Affymetrix 6.0 genotyping array for signatures of natural selection in immune response genes. We applied three statistical tests to identify signatures of natural selection: locus-specific branch length (LSBL), the cross-population extended haplotype homozygosity (XP-EHH), and the integrated haplotype score (iHS). Each of the haplotype tests (XP-EHH and iHS) were paired with LSBL and significance was determined at the 1% level. For the paired analyses, we identified 95 statistically significant windows for XP-EHH/LSBL and 63 statistically significant windows for iHS/LSBL. Among our top immune response loci, we found evidence of recent directional selection associated with the major histocompatibility complex (MHC) and the peroxisome proliferator-activated receptor gamma (PPAR-γ) signaling pathway. These findings illustrate that Mesoamerican populations' immunity has been shaped by exposure to infectious disease. As targets of selection, these variants are likely to encode phenotypes that manifest themselves physiologically and therefore may contribute to population-level variation in immune response. Our results shed light on past selective events influencing the host response to modern diseases, both pathogenic infection as well as autoimmune disorders.

Genome structural variation in human evolution

Structural variation (SV) is a large difference (typically >100 bp) in the genomic structure of two genomes and includes both copy number variation and variation that does not change copy number of a genomic region, such as an inversion. Improved reference genomes, combined with widespread genome sequencing using short-read sequencing technology, and increasingly using long-read sequencing, have reignited interest in SV. Recent large-scale studies and functional focused analyses have highlighted the role of SV in human evolution. In this review, we highlight human-specific SVs involved in changes in the brain, population-specific SVs that affect response to the environment, including adaptation to diet and infectious diseases, and summarise the contribution of archaic hominin admixture to present-day human SV.

SNPs in IL4 and IFNG show no protective associations with human African trypanosomiasis in the Democratic Republic of the Congo: a case-control study

Background: Human African trypanosomiasis (HAT) is a protozoal disease transmitted by tsetse flies. Infection with trypanosomes can lead directly to active HAT or latent infection with no detectable parasites, which may progress to active HAT or to spontaneous self-cure. Genetic variation could explain these differences in the outcome of infection. To test this hypothesis, polymorphisms in 17 candidate genes were tested ( APOL1 [ G1 and G2], CFH, HLA-A, HPR, HP, IL1B, IL12B, IL12RB1, IL10, IL4R, MIF, TNFA , IL6, IL4, IL8, IFNG, and HLA-G). Methods: Samples were collected in Democratic Republic of the Congo. 233 samples were genotyped: 100 active HAT cases, 33 from subjects with latent infections and 100 negative controls. Commercial service providers genotyped polymorphisms at 96 single nucleotide polymorphisms (SNPs) on 17 genes. Data were analyzed using Plink V1.9 software and R. Loci, with suggestive associations (uncorrected p < 0.05) validated using an additional 594 individuals, including 164 cases and 430 controls. Results: After quality control, 87 SNPs remained in the analysis. Two SNPs in IL4 and two in IFNG were suggestively associated (uncorrected p<0.05) with a differential risk of developing a Trypanosoma brucei gambiense infection in the Congolese population. The IFNG minor allele (rs2430561, rs2069718) SNPs were protective in comparison between latent infections and controls. Carriers of the rs2243258_T and rs2243279_A alleles of IL4 and the rs2069728_T allele of IFNG had a reduced risk of developing illness or latent infection, respectively. None of these associations were significant after Bonferroni correction for multiple testing. A validation study using more samples was run to determine if the absence of significant association was due to lack of power. Conclusions: This study showed no evidence of an association of HAT with IL4 and IFNG SNPs or with APOL1 G1 and G2 alleles, which have been found to be protective in other studies.

Association of rs2000999 in the haptoglobin gene with total cholesterol, HDL-C, and LDL-C levels in Mexican type 2 diabetes patients

Recently, studies have shown significant association between the rs2000999 polymorphism in the haptoglobin-encoding gene (HP) and low-density lipoprotein cholesterol (LDL-C) and total cholesterol (TC) levels, which are important risk factors for cardiovascular diseases. However, the association of rs2000999 with serum lipids in Latin American diabetic populations is still uncharacterized. Here, we analyzed the association of rs2000999 with TC, high-density lipoprotein cholesterol (HDL-C), and LDL-C levels in 546 Mexican adults with type 2 diabetes (T2D) and in 654 controls without T2D. In this observational case-control study we included adults from 4 centers of the Mexican Social Security Institute in Mexico City recruited from 2012 to 2015. TC, HDL-C, LDL-C, triglycerides (TG), and glucose levels were measured by an enzymatic colorimetric method. The variant rs2000999 was genotyped using TaqMan real time polymerase chain reaction. The percentage of Native-American ancestry showed a negative association with the rs2000999 A allele. In contrast, the rs2000999 A allele had a strong positive association with European ancestry, and to a lesser extent, with African ancestry. Linear regression was used to estimate the association between the variant rs2000999 and lipid concentrations, using different genetic models. Under codominant and recessive models, rs2000999 was significantly associated with TC and LDL-C levels in the T2D group and in controls without T2D. In addition, the group with T2D showed a significant association between the variant and HDL-C levels. In summary, the rs2000999 A allele in Mexican population is positively associated with the percentage of European and negatively associated with Native American ancestry. Carriers of the A allele have increased levels of TC and LDL-C, independently of T2D diagnosis, and also increased concentrations of HDL-C in the T2D sample.

Association of APOL1 renal disease risk alleles with Trypanosoma brucei rhodesiense infection outcomes in the northern part of Malawi

Trypanosoma brucei (T.b.) rhodesiense is the cause of the acute form of human African trypanosomiasis (HAT) in eastern and southern African countries. There is some evidence that there is diversity in the disease progression of T.b. rhodesiense in different countries. HAT in Malawi is associated with a chronic haemo-lymphatic stage infection compared to other countries, such as Uganda, where the disease is acute with more marked neurological impairment. This has raised the question of the role of host genetic factors in infection outcomes. A candidate gene association study was conducted in the northern region of Malawi. This was a case-control study involving 202 subjects, 70 cases and 132 controls. All individuals were from one area; born in the area and had been exposed to the risk of infection since birth. Ninety-six markers were genotyped from 17 genes: IL10, IL8, IL4, HLA-G, TNFA, IL6, IFNG, MIF, APOL, HLA-A, IL1B, IL4R, IL12B, IL12R, HP, HPR, and CFH. There was a strong significant association with APOL1 G2 allele (p = 0.0000105, OR = 0.14, CI95 = [0.05-0.41], BONF = 0.00068) indicating that carriers of the G2 allele were protected against T.b. rhodesiense HAT. SNP rs2069845 in IL6 had raw p < 0.05, but did not remain significant after Bonferroni correction. There were no associations found with the other 15 candidate genes. Our finding confirms results from other studies that the G2 variant of APOL1 is associated with protection against T.b. rhodesiense HAT.

Association between IL1 gene polymorphism and human African trypanosomiasis in populations of sleeping sickness foci of southern Cameroon

BACKGROUND

Human African Trypanosomiasis (HAT) is a neglected tropical disease caused by infections due to Trypanosoma brucei subspecies. In addition to the well-established environmental and behavioural risks of becoming infected, there is evidence for a genetic component to the response to trypanosome infection. We undertook a candidate gene case-control study to investigate genetic associations further.

METHODOLOGY

We genotyped one polymorphism in each of seven genes (IL1A, IL1RN, IL4RN, IL6, HP, HPR, and HLA-G) in 73 cases and 250 controls collected from 19 ethno-linguistic subgroups stratified into three major ethno-linguistic groups, 2 pooled ethno-linguistic groups and 11 ethno-linguistic subgroups from three Cameroonian HAT foci. The seven polymorphic loci tested consisted of three SNPs, three variable numbers of tandem repeat (VNTR) and one INDEL.

RESULTS

We found that the genotype (TT) and minor allele (T) of IL1A gene as well as the genotype 1A3A of IL1RN were associated with an increased risk of getting Trypanosoma brucei gambiense and develop HAT when all data were analysed together and also when stratified by the three major ethno-linguistic groups, 2 pooled ethno-linguistic subgroups and 11 ethno-linguistic subgroups.

CONCLUSION

This study revealed that one SNP rs1800794 of IL1A and one VNTR rs2234663 of IL1RN were associated with the increased risk to be infected by Trypanosoma brucei gambiense and develop sleeping sickness in southern Cameroon. The minor allele T and the genotype TT of SNP rs1800794 in IL1A as well as the genotype 1A3A of IL1RN rs2234663 VNTR seem to increase the risk of getting Trypanosoma brucei gambiense infections and develop sleeping sickness in southern Cameroon.

Association of the genetic variant rs2000999 with haptoglobin and diabetic macrovascular diseases in Chinese patients with type 2 diabetes

AIMS

The common copy number variant (CNV) in the haptoglobin (Hp) gene may influence the susceptibility to diabetic macrovascular diseases. We aimed to investigate the relationship of the genetic variant rs2000999, located in the haptoglobin-related protein (HPR) gene, with serum Hp levels and diabetic macrovascular diseases in Chinese type 2 diabetes patients.

METHODS

The Hp CNV and rs2000999 were genotyped in a group of 5457 Chinese patients with type 2 diabetes. Associations of rs2000999 with the common Hp CNV, susceptibility to diabetic macrovascular diseases and related metabolic traits were analysed. Furthermore, 886 patients were selected to detect serum Hp levels and to evaluate the correlation between rs2000999 and serum Hp levels.

RESULTS

The genetic variant rs2000999 was not associated with diabetic macrovascular diseases (P = 0.6109), while subjects carrying the A allele had higher levels of low-density lipoprotein cholesterol (P = 0.0578) and a smaller inter-adventitial diameter of the common carotid artery (P = 0.0266). Additionally, rs2000999 exhibited strong association with serum Hp levels (P = 2.03 × 10^-21).

CONCLUSIONS

The genetic variant rs2000999 was not associated with diabetic macrovascular diseases but showed an association with metabolic traits and serum Hp levels in Chinese patients with type 2 diabetes.

The Diversity of Mammalian Hemoproteins and Microbial Heme Scavengers Is Shaped by an Arms Race for Iron Piracy

Iron is an essential micronutrient for most living species. In mammals, hemoglobin (Hb) stores more than two thirds of the body's iron content. In the bloodstream, haptoglobin (Hp) and hemopexin (Hpx) sequester free Hb or heme. Pathogenic microorganisms usually acquire iron from their hosts and have evolved complex systems of iron piracy to circumvent nutritional immunity. Herein, we performed an evolutionary analysis of genes coding for mammalian heme-binding proteins and heme-scavengers in pathogen species. The underlying hypothesis is that these molecules are engaged in a molecular arms race. We show that positive selection drove the evolution of mammalian Hb and Hpx. Positively selected sites in Hb are located at the interaction surface with Neisseria meningitidis heme scavenger HpuA and with Staphylococcus aureus iron-regulated surface determinant B (IsdB). In turn, positively selected sites in HpuA and IsdB are located in the flexible protein regions that contact Hb. A residue in Hb (S45H) was also selected on the Caprinae branch. This site stabilizes the interaction with Trypanosoma brucei hemoglobin-haptoglobin (HbHp) receptor (TbHpHbR), a molecule that also mediates trypanosome lytic factor (TLF) entry. In TbHpHbR, positive selection drove the evolution of a variant (L210S) which allows evasion from TLF but reduces affinity for HbHp. Finally, selected sites in Hpx are located at the interaction surface with the Haemophilus influenzae hemophore HxuA, which in turn displays fast evolving sites at the Hpx-binding interface. These results shed light into host-pathogens conflicts and establish the importance of nutritional immunity as an evolutionary force.

Pan-Genomic Study of Mycobacterium tuberculosis Reflecting the Primary/Secondary Genes, Generality/Individuality, and the Interconversion Through Copy Number Variations

Tuberculosis (TB) has surpassed HIV as the leading infectious disease killer worldwide since 2014. The main pathogen, Mycobacterium tuberculosis (Mtb), contains ~4,000 genes that account for ~90% of the genome. However, it is still unclear which of these genes are primary/secondary, which are responsible for generality/individuality, and which interconvert during evolution. Here we utilized a pan-genomic analysis of 36 Mtb genomes to address these questions. We identified 3,679 Mtb core (i.e., primary) genes, determining their phenotypic generality (e.g., virulence, slow growth, dormancy). We also observed 1,122 dispensable and 964 strain-specific secondary genes, reflecting partially shared and lineage-/strain-specific individualities. Among which, five L2 lineage-specific genes might be related to the increased virulence of the L2 lineage. Notably, we discovered 28 Mtb “Super Core Genes” (SCGs: more than a copy in at least 90% strains), which might be of increased importance, and reflected the “super phenotype generality.” Most SCGs encode PE/PPE, virulence factors, antigens, and transposases, and have been verified as playing crucial roles in Mtb pathogenicity. Further investigation of the 28 SCGs demonstrated the interconversion among SCGs, single-copy core, dispensable, and strain-specific genes through copy number variations (CNVs) during evolution; different mutations on different copies highlight the delicate adaptive-evolution regulation amongst Mtb lineages. This reflects that the importance of genes varied through CNVs, which might be driven by selective pressure from environment/host-adaptation. In addition, compared with Mycobacterium bovis (Mbo), Mtb possesses 48 specific single core genes that partially reflect the differences between Mtb and Mbo individuality.

Do Cryptic Reservoirs Threaten Gambiense-Sleeping Sickness Elimination?

Trypanosoma brucei gambiense causes human African trypanosomiasis (HAT). Between 1990 and 2015, almost 440000 cases were reported. Large-scale screening of populations at risk, drug donations, and efforts by national and international stakeholders have brought the epidemic under control with <2200 cases in 2016. The World Health Organization (WHO) has set the goals of gambiense-HAT elimination as a public health problem for 2020, and of interruption of transmission to humans for 2030. Latent human infections and possible animal reservoirs may challenge these goals. It remains largely unknown whether, and to what extend, they have an impact on gambiense-HAT transmission. We argue that a better understanding of the contribution of human and putative animal reservoirs to gambiense-HAT epidemiology is mandatory to inform elimination strategies.

No evidence for association between APOL1 kidney disease risk alleles and Human African Trypanosomiasis in two Ugandan populations

BACKGROUND

Human African trypanosomiasis (HAT) manifests as an acute form caused by Trypanosoma brucei rhodesiense (Tbr) and a chronic form caused by Trypanosoma brucei gambiense (Tbg). Previous studies have suggested a host genetic role in infection outcomes, particularly for APOL1. We have undertaken candidate gene association studies (CGAS) in a Ugandan Tbr and a Tbg HAT endemic area, to determine whether polymorphisms in IL10, IL8, IL4, HLAG, TNFA, TNX4LB, IL6, IFNG, MIF, APOL1, HLAA, IL1B, IL4R, IL12B, IL12R, HP, HPR, and CFH have a role in HAT.

METHODOLOGY AND RESULTS

We included 238 and 202 participants from the Busoga Tbr and Northwest Uganda Tbg endemic areas respectively. Single Nucleotide Polymorphism (SNP) genotype data were analysed in the CGAS. The study was powered to find odds ratios > 2 but association testing of the SNPs with HAT yielded no positive associations i.e. none significant after correction for multiple testing. However there was strong evidence for no association with Tbr HAT and APOL1 G2 of the size previously reported in the Kabermaido district of Uganda.

CONCLUSIONS/SIGNIFICANCE

A recent study in the Soroti and Kaberamaido focus in Central Uganda found that the APOL1 G2 allele was strongly associated with protection against Tbr HAT (odds ratio = 0.2, 95% CI: 0.07 to 0.48, p = 0.0001). However, in our study no effect of G2 on Tbr HAT was found, despite being well powered to find a similar sized effect (OR = 0.9281, 95% CI: 0.482 to 1.788, p = 0.8035). It is possible that the G2 allele is protective from Tbr in the Soroti/Kabermaido focus but not in the Iganga district of Busoga, which differ in ethnicity and infection history. Mechanisms underlying HAT infection outcome and virulence are complex and might differ between populations, and likely involve several host, parasite or even environmental factors.

Potential role of ferric hemoglobin in MS pathogenesis: Effects of oxidative stress and extracellular methemoglobin or its degradation products on myelin components

There is a well-documented relationship between cerebral vasculature and multiple sclerosis (MS) lesions: abnormal accumulations of iron have been found in the walls of the dilated veins in cerebral MS plaques. The source of this iron is unknown, but could be related to the recognized phenomenon of capillary and venous hemorrhages leading to blood extravasation. In turn, hemorrhaging leading to hemolysis results in extracellular release of hemoglobin, a reactive molecule that could induce local oxidative stress, inflammation, and tissue damage. Our previous studies with a reduced form of hemoglobin (oxyHb) have demonstrated its ability to cause extensive lipid and protein oxidation in vitro, which would result in membrane destabilization. Here, we investigated in further detail the mechanism by which the more abundant oxidized form of extracellular hemoglobin (metHb), and dissociated hemin, cause direct oxidative damage to myelin components, specifically membrane-mimetic lipid vesicles and myelin basic protein (MBP), a highly-abundant protein in the CNS. Oxidation of lipids was assessed by the formation of conjugated diene/triene and malondialdehyde, and oxidation of MBP was demonstrated by the bityrosine formation and by the change in protein mass. Our results show that metHb causes oxidative damage to MBP and myelin lipids, partly by transferring its hemin moiety to protein and lipid, but mostly as an intact protein possibly via formation of a ferryl radical. These results elucidating the mechanism of extracellular hemoglobin-induced oxidative damage to myelin components support the need for further research into vascular pathology in MS pathogenesis, to gain insight into the role of iron deposits and/or in stimulation of different comorbidities associated with the disease.

A polymorphism in the haptoglobin, haptoglobin related protein locus is associated with risk of human sleeping sickness within Cameroonian populations

Background

Human African Trypanosomiasis (HAT) is a neglected disease targeted for elimination as a public health problem by 2020. Elimination requires a better understanding of the epidemiology and clinical evolution of HAT. In addition to the classical clinical evolution of HAT, asymptomatic carriers and spontaneous cure have been reported in West Africa. A genetic component to human susceptibility to HAT has been suggested to explain these newly observed responses to infection. In order to test for genetic associations with infection response, genetic polymorphism in 17 genes were tested (APOL1, IL1B, IL4, IL4R, IL6, IL8, IL12B, IL12RB1, IL10, TNFA, INFG, MIF, HLA-G, HLA-A, HP, HPR and CFH).

Methodology

A case-control study was performed on 180 blood samples collected from 56 cases and 124 controls from Cameroon. DNA was extracted from blood samples. After quality control, 25 samples (24 controls and 1 case) were eliminated. The genotyping undertaken on 155 individuals including 55 cases and 100 controls were investigated at 96 loci (88 SNPs and 8 indels) located on 17 genes. Associations between these loci and HAT were estimated via a case-control association test.

Results

Analyses of 64 SNPs and 4 indels out of 96 identified in the selected genes reveal that the minor allele (T) of rs8062041 in haptoglobin (HP) appeared to be protective against HAT (p = 0.0002395, OR 0.359 (CI₉₅ [0.204–0.6319])); indicating higher frequency in cases compared to controls. This minor allele with adjusted p value of 0.0163 is associated with a lower risk (protective effect) of developing sleeping sickness.

Conclusion

The haptoglobin related protein HPR and HP are tightly linked and both are duplicated in some people and may lead to higher activity. This increased production could be responsible of the protection associated with rs8062041 even though this SNP is within HP.

Human African trypanosomiasis (HAT) or sleeping sickness is a neglected tropical disease targeted for elimination by 2020. This elimination requires a better understanding of the epidemiology and clinical evolution of this disease. Beside the classical clinical evolution, asymptomatic carriers, seropositive and spontaneous cure of infected persons have been reported in West Africa. Arguments in favor of human genetic susceptibility to HAT have been raised to explain this variability in clinical presentation. This study investigated the genetic polymorphism of 17 genes between controls and sleeping sickness patients in Southern Cameroon in order to improve our knowledge of human susceptibility to trypanosome infections. We identified single nucleotide polymorphisms and indels in 17 selected genes involved in immune responses and carried out a case-control candidate gene association study and demonstrated differences between variants associated with the disease. From these genes, only haptoglobin (HP) at the SNP rs8062041 was found to have polymorphisms which were strongly associated with trypanosomiasis. The minor allele (T) at this SNP position appeared to be protective against HAT (p = 0.0002395, OR 0.359 (CI₉₅ [0.204–0.6319])) reducing the risk of developing disease approximately threefold. The haptoglobin related protein (HPR) is adjacent to HP and is a component of the Trypanolytic factor that kills trypanosomes. The HP and HPR locus is duplicated in some people. The rs8062041 variant may be associated with this duplication and it is possible that increased production of HPR is the cause of the protection associated with rs8062041. The results reported here will contribute to the knowledge of the role of human genetics in disease progression, and thus lead to the identification of novel biomarkers which could involve development of new diagnostics, treatments and intervention strategies.

Transcription and translation of APOL1 variants

It is highly important to document the molecular alterations existing in normal cells prior to the onset of any disease. Knowledge of genetic mutations and associated molecular mechanisms will be helpful for better diagnosis and management of disease. The major focus of this commentary on providing understanding about the apolipoprotein 1 (APOL1) gene, the protein encoded by this gene (apoL1) and the mechanistic details regarding the role of apoL1 in the lysis of Trypanosoma brucei. Information about APOL1 genetic variants, APOL1G1 and APOL1G2, is provided along with the association of these variants with hypertension-attributed end-stage renal disease (ESRD) and focal segmental glomerulosclerosis (FSGS). Moreover, this commentary presents a brief overview of how the authors of a recent Bioscience Reports article [Haque et al (2017) 37, BSR20160531, doi: 10.1042/BSR20160531] have evaluated the functional impact of G1 and G2 alleles on the transcription and translation of APOL1 mRNA.

Dosage-sensitive genes in evolution and disease

For a subset of genes in our genome a change in gene dosage, by duplication or deletion, causes a phenotypic effect. These dosage-sensitive genes may confer an advantage upon copy number change, but more typically they are associated with disease, including heart disease, cancers and neuropsychiatric disorders. This gene copy number sensitivity creates characteristic evolutionary constraints that can serve as a diagnostic to identify dosage-sensitive genes. Though the link between copy number change and disease is well-established, the mechanism of pathogenicity is usually opaque. We propose that gene expression level may provide a common basis for the pathogenic effects of many copy number variants.

A co-evolutionary arms race: trypanosomes shaping the human genome, humans shaping the trypanosome genome

Trypanosoma brucei is the causative agent of African sleeping sickness in humans and one of several pathogens that cause the related veterinary disease Nagana. A complex co-evolution has occurred between these parasites and primates that led to the emergence of trypanosome-specific defences and counter-measures. The first line of defence in humans and several other catarrhine primates is the trypanolytic protein apolipoprotein-L1 (APOL1) found within two serum protein complexes, trypanosome lytic factor 1 and 2 (TLF-1 and TLF-2). Two sub-species of T. brucei have evolved specific mechanisms to overcome this innate resistance, Trypanosoma brucei gambiense and Trypanosoma brucei rhodesiense. In T. b. rhodesiense, the presence of the serum resistance associated (SRA) gene, a truncated variable surface glycoprotein (VSG), is sufficient to confer resistance to lysis. The resistance mechanism of T. b. gambiense is more complex, involving multiple components: reduction in binding affinity of a receptor for TLF, increased cysteine protease activity and the presence of the truncated VSG, T. b. gambiense-specific glycoprotein (TgsGP). In a striking example of co-evolution, evidence is emerging that primates are responding to challenge by T. b. gambiense and T. b. rhodesiense, with several populations of humans and primates displaying resistance to infection by these two sub-species.

Correlation of geographic distributions of haptoglobin alleles with prevalence of multiple sclerosis (MS) - a narrative literature review

We have proposed that the myelin damage observed in multiple sclerosis (MS) may be partly mediated through the long-term release and degradation of extracellular hemoglobin (Hb) and the products of its oxidative degradation [Cellular and Molecular Life Sciences, 71, 1789-1798, 2014]. The protein haptoglobin (Hpt) binds extracellular Hb as a first line of defense, and can serve as a vascular antioxidant. Humans have two different Hpt alleles: Hpt1 and Hpt2, giving either homozygous Hpt1-1 or Hpt2-2 phenotypes, or a heterozygous Hpt1-2 phenotype. We questioned whether those geographic regions with higher frequency of the Hpt2 allele (conversely, lower frequency of Hpt1 allele) would correlate with an increased incidence of MS, because different Hpt phenotypes will have variable anti-oxidative potentials in protecting myelin from damage inflicted by extracellular Hb and its degradation products. To test this hypothesis, we undertook a systematic analysis of the literature on reported geographic distributions of Hpt alleles to compare them with data reported in the World Health Organization Atlas of worldwide MS prevalence. We found the frequency of the Hpt1 allele to be low in European and North American countries with a high prevalence of MS, consistent with our hypothesis. However, this correlation was not observed in China and India, countries with the lowest Hpt1 frequencies, yet low reported prevalence of MS. Nevertheless, this work shows the need for continued refinement of geographic patterns of MS prevalence, including data on ethnic or racial origin, and for new clinical studies to probe the observed correlation and evaluate Hpt phenotype as a predictor of disease variability and progression, severity, and/or comorbidity with cardiovascular disorders.

Analysis of Copy Number Variation Using the Paralogue Ratio Test (PRT)

Copy number variation (CNV), where a segment of DNA differs in copy number between different individuals, is an extensive and often underappreciated source of genetic variation within species. However, reliably determining copy number of a particular DNA sequence for a large number of samples can be challenging. Here, I describe and review the paralogue ratio test (PRT) in detail. PRT was developed to robustly type the CNV of the beta-defensin locus using small amounts of genomic DNA in a high-throughput manner, and has been applied successfully at many other loci. I discuss the strategies for designing successful PRT assays using both manual and bioinformatics methods, how to optimize experimental conditions, and approaches for analyzing the data. I discuss strengths and weaknesses of the approach, and how to troubleshoot results, as well as the range of problems to which PRT can be a potential solution.

Copy number variations in Saudi family with intellectual disability and epilepsy

Background

Epilepsy is genetically complex but common brain disorder of the world affecting millions of people with almost of all age groups. Novel Copy number variations (CNVs) are considered as important reason for the numerous neurodevelopmental disorders along with intellectual disability and epilepsy. DNA array based studies contribute to explain a more severe clinical presentation of the disease but interoperation of many detected CNVs are still challenging.

Results

In order to study novel CNVs with epilepsy related genes in Saudi family with six affected and two normal individuals with several forms of epileptic seizures, intellectual disability (ID), and minor dysmorphism, we performed the high density whole genome Agilent sure print G3 Hmn CGH 2x 400 K array-CGH chips analysis. Our results showed de novo deletions, duplications and deletion plus duplication on differential chromosomal regions in the affected individuals that were not shown in the normal fathe and normal kids by using Agilent CytoGenomics 3.0.6.6 softwear. Copy number gain were observed in the chromosome 1, 16 and 22 with LCE3C, HPR, GSTT2, GSTTP2, DDT and DDTL genes respectively whereas the deletions observed in the chromosomal regions 8p23-p21 (4303127–4337759) and the potential gene in this region is CSMD1 (OMIM: 612279). Moreover, the array CGH results deletions and duplication were also validated by using primer design of deleted regions utilizing the flanked SNPs using simple PCR and also by using quantitative real time PCR.

Conclusions

We found some of the de novo deletions and duplication in our study in Saudi family with intellectual disability and epilepsy. Our results suggest that array-CGH should be used as a first line of genetic test for epilepsy except there is a strong indication for a monogenic syndrome. The advanced high through put array-CGH technique used in this study aim to collect the data base and to identify new mechanisms describing epileptic disorder, may help to improve the clinical management of individual cases in decreasing the burden of epilepsy in Saudi Arabia.

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/.

Recurring exon deletions in the HP (haptoglobin) gene contribute to lower blood cholesterol levels

Two exons of the human haptoglobin (HP) gene exhibit copy number variation that affects HP multimerization and underlies one of the first protein polymorphisms identified in humans. The evolutionary origins and medical significance of this polymorphism have been uncertain. Here we show that this variation has likely arisen from the recurring reversion of an ancient hominin-specific duplication of these exons. Though this polymorphism has been largely invisible to genome-wide genetic studies to date, we describe a way to analyze it by imputation from SNP haplotypes and find among 22,288 individuals that these HP exonic deletions associate with reduced LDL and total cholesterol levels. We show that these deletions, and a SNP that affects HP expression, are the likely drivers of the strong but complex association of cholesterol levels to SNPs near HP. Recurring exonic deletions in the haptoglobin gene likely enhance human health by lowering cholesterol levels in the blood.

Global diversity, population stratification, and selection of human copy-number variation

In order to explore the diversity and selective signatures of duplication and deletion human copy-number variants (CNVs), we sequenced 236 individuals from 125 distinct human populations. We observed that duplications exhibit fundamentally different population genetic and selective signatures than deletions and are more likely to be stratified between human populations. Through reconstruction of the ancestral human genome, we identify megabases of DNA lost in different human lineages and pinpoint large duplications that introgressed from the extinct Denisova lineage now found at high frequency exclusively in Oceanic populations. We find that the proportion of CNV base pairs to single-nucleotide-variant base pairs is greater among non-Africans than it is among African populations, but we conclude that this difference is likely due to unique aspects of non-African population history as opposed to differences in CNV load.

Host genetics and parasitic infections

Parasites still impose a high death and disability burden on human populations, and are therefore likely to act as selective factors for genetic adaptations. Genetic epidemiological investigation of parasitic diseases is aimed at disentangling the mechanisms underlying immunity and pathogenesis by looking for associations or linkages between loci and susceptibility phenotypes. Until recently, most studies used a candidate gene approach and were relatively underpowered, with few attempts at replicating findings in different populations. However, in the last 5 years, genome-wide and/or multicentre studies have been conducted for severe malaria, visceral leishmaniasis, and cardiac Chagas disease, providing some novel important insights. Furthermore, studies of helminth infections have repeatedly shown the involvement of common loci in regulating susceptibility to distinct diseases such as schistosomiasis, ascariasis, trichuriasis, and onchocherciasis. As more studies are conducted, evidence is increasing that at least some of the identified susceptibility loci are shared not only among parasitic diseases but also with immunological disorders such as allergy or autoimmune disease, suggesting that parasites may have played a role in driving the evolution of the immune system.

Human gene copy number variation and infectious disease

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