The Growing Importance of CNVs: New Insights for Detection and Clinical Interpretation.

Differences between genomes can be due to single nucleotide variants, translocations, inversions, and copy number variants (CNVs, gain or loss of DNA). The latter can range from sub-microscopic events to complete chromosomal aneuploidies. Small CNVs are often benign but those larger than 500 kb are strongly associated with morbid consequences such as developmental disorders and cancer. Detecting CNVs within and between populations is essential to better understand the plasticity of our genome and to elucidate its possible contribution to disease. Hence there is a need for better-tailored and more robust tools for the detection and genome-wide analyses of CNVs. While a link between a given CNV and a disease may have often been established, the relative CNV contribution to disease progression and impact on drug response is not necessarily understood. In this review we discuss the progress, challenges, and limitations that occur at different stages of CNV analysis from the detection (using DNA microarrays and next-generation sequencing) and identification of recurrent CNVs to the association with phenotypes. We emphasize the importance of germline CNVs and propose strategies to aid clinicians to better interpret structural variations and assess their clinical implications.

Decreased PM10 exposure attenuates age-related lung function decline: genetic variants in p53, p21, and CCND1 modify this effect.

BACKGROUND: Decreasing exposure to airborne particulates was previously associated with reduced age-related decline in lung function. However, whether the benefit from improved air quality depends on genetic background is not known. Recent evidence points to the involvement of the genes p53 and p21 and of the cell cycle control gene cyclin D1 (CCND1) in the response of bronchial cells to air pollution. OBJECTIVE: We determined in 4,326 participants of the Swiss Cohort Study on Air Pollution and Lung and Heart Diseases in Adults (SAPALDIA) whether four single-nucleotide polymorphisms in three genes [CCND1 (rs9344 [P242P], rs667515), p53 (rs1042522 [R72P]), and p21 (rs1801270 [S31R])] modified the previously observed attenuation of the decline in the forced expiratory flow between 25% and 75% of the forced vital capacity (FEF(25-75)) associated with improved air quality. METHODS: Subjects of the prospective population-based SAPALDIA cohort were assessed in 1991 and 2002 by spirometry, questionnaires, and biological sample collection for genotyping. We assigned spatially resolved concentrations of particulate matter with aerodynamic diameter < or = 10 microm (PM(10)) to each participant's residential history 12 months before the baseline and follow-up assessments. RESULTS: The effect of diminishing PM(10) exposure on FEF(25-75) decline appeared to be modified by p53 R72P, CCND1 P242P, and CCND1 rs667515. For example, a 10-microg/m(3) decline in average PM(10) exposure over an 11-year period attenuated the average annual decline in FEF(25-75) by 21.33 mL/year (95% confidence interval, 10.57-32.08) among participants homozygous for the CCND1 (P242P) GG genotype, by 13.72 mL/year (5.38-22.06) among GA genotypes, and by 6.00 mL/year (-4.54 to 16.54) among AA genotypes. CONCLUSIONS: Our results suggest that cell cycle control genes may modify the degree to which improved air quality may benefit respiratory function in adults.

A large study of androgen receptor germline variants and their relation to sex hormone levels and prostate cancer risk. Results from the National Cancer Institute Breast and Prostate Cancer Cohort Consortium

Background: Androgens are key regulators of prostate gland maintenance and prostate cancer growth, and androgen deprivation therapy has been the mainstay of treatment for advanced prostate cancer for many years. A long-standing hypothesis has been that inherited variation in the androgen receptor (AR) gene plays a role in prostate cancer initiation. However, studies to date have been inconclusive and often suffered from small sample sizes. Objective and Methods: We investigated the association of AR sequence variants with circulating sex hormone levels and prostate cancer risk in 6058 prostate cancer cases and 6725 controls of Caucasian origin within the Breast and Prostate Cancer Cohort Consortium. We genotyped a highly polymorphic CAG microsatellite in exon 1 and six haplotype tagging single nucleotide polymorphisms and tested each genetic variant for association with prostate cancer risk and with sex steroid levels. Results: We observed no association between AR genetic variants and prostate cancer risk. However, there was a strong association between longer CAG repeats and higher levels of testosterone (P = 4.73 × 10−5) and estradiol (P = 0.0002), although the amount of variance explained was small (0.4 and 0.7%, respectively). Conclusions: This study is the largest to date investigating AR sequence variants, sex steroid levels, and prostate cancer risk. Although we observed no association between AR sequence variants and prostate cancer risk, our results support earlier findings of a relation between the number of CAG repeats and circulating levels of testosterone and estradiol.

A genome-wide association study reveals variants in ARL15 that influence adiponectin levels.

The adipocyte-derived protein adiponectin is highly heritable and inversely associated with risk of type 2 diabetes mellitus (T2D) and coronary heart disease (CHD). We meta-analyzed 3 genome-wide association studies for circulating adiponectin levels (n = 8,531) and sought validation of the lead single nucleotide polymorphisms (SNPs) in 5 additional cohorts (n = 6,202). Five SNPs were genome-wide significant in their relationship with adiponectin (P< or =5x10(-8)). We then tested whether these 5 SNPs were associated with risk of T2D and CHD using a Bonferroni-corrected threshold of P< or =0.011 to declare statistical significance for these disease associations. SNPs at the adiponectin-encoding ADIPOQ locus demonstrated the strongest associations with adiponectin levels (P-combined = 9.2x10(-19) for lead SNP, rs266717, n = 14,733). A novel variant in the ARL15 (ADP-ribosylation factor-like 15) gene was associated with lower circulating levels of adiponectin (rs4311394-G, P-combined = 2.9x10(-8), n = 14,733). This same risk allele at ARL15 was also associated with a higher risk of CHD (odds ratio [OR] = 1.12, P = 8.5x10(-6), n = 22,421) more nominally, an increased risk of T2D (OR = 1.11, P = 3.2x10(-3), n = 10,128), and several metabolic traits. Expression studies in humans indicated that ARL15 is well-expressed in skeletal muscle. These findings identify a novel protein, ARL15, which influences circulating adiponectin levels and may impact upon CHD risk.

Candidate gene analysis of ocular toxoplasmosis in Brazil: evidence for a role for toll-like receptor 9 (TLR9)

Toxoplasma gondii infection is an important mediator of ocular disease in Brazil more frequently than reported from elsewhere. Infection and pathology are characterized by a strong proinflammatory response which in mice is triggered by interaction of the parasite with the toll-like receptor (TLR)/MyD88 pathway. A powerful way to identify the role of TLRs in humans is to determine whether polymorphisms at these loci influence susceptibility to T. gondii-mediated pathologies. Here we report on a small family-based study (60 families; 68 affected offspring) undertaken in Brazil which was powered for large effect sizes using single nucleotide polymorphisms with minor alleles frequencies > 0.3. Of markers in TLR2, TLR5 and TLR9 that met these criteria, we found an association Family Based Association Tests [(FBAT) Z score = 4.232; p = 1.5 x 10-5; p corrected = 1.2 x 10-4] between the C allele (frequency = 0.424; odds ratio = 7; 95% confidence interval 1.6-30.8) of rs352140 at TLR9 and toxoplasmic retinochoroiditis in Brazil. This supports the hypothesis that direct interaction between T. gondii and TLR9 may trigger proinflammatory responses that lead to severe pathologies such as the ocular disease that is associated with this infection in Brazil.

High prevalence and low E6 genetic variability of human papillomavirus 58 in women with cervical cancer and precursor lesions in Southeast Mexico

Infection with some genotypes of human papillomavirus (HPV) is the most important risk factor associated with cervical cancer (CC). Throughout the world, HPV type 58 prevalence varies from one region to another; it is higher in women from certain countries in Asia and Latin America, such as China and Mexico. Although intratypic variants have been reported on a few occasions, our knowledge about HPV 58 genetic variation remains limited. Therefore, this work aims to (i) determine the prevalence of HPV type 58 amongst Mexican women with invasive CC or precursor lesions and (ii) identify HPV 58 sequence variants. One hundred and forty five colposcopy clinic patients were studied. Genotyping of HPV 16, 18 and 58 was determined by specific nested PCR and HPV 58 variants were detected by direct sequencing. The general prevalence of HPV was 51.7% (75/145). HPV 16 was found in 30.6% (23/75) and HPV 58 in 24% (18/75) of the patients. HPV 18 was not identified in patients with cervical intraepithelial neoplasia (CIN) grade I; it was only found in those with CIN II, with a prevalence of 6.8% (3/44). In patients with CC, the prevalence of HPV 16 and 58 was 78.9%. Regarding HPV 58 variants, 94.4% of the HPV 58 sequences were identical to the prototype strain, whereas one sample showed changes at a single nucleotide. This study demonstrates a high prevalence of HPV 58 and a low genetic variability of E6 sequences amongst Mexican colposcopy patients.

Additive effects of LPL, APOA5 and APOE variant combinations on triglyceride levels and hypertriglyceridemia: results of the ICARIA genetic sub-study.

BACKGROUND Hypertriglyceridemia (HTG) is a well-established independent risk factor for cardiovascular disease and the influence of several genetic variants in genes related with triglyceride (TG) metabolism has been described, including LPL, APOA5 and APOE. The combined analysis of these polymorphisms could produce clinically meaningful complementary information. METHODS A subgroup of the ICARIA study comprising 1825 Spanish subjects (80% men, mean age 36 years) was genotyped for the LPL-HindIII (rs320), S447X (rs328), D9N (rs1801177) and N291S (rs268) polymorphisms, the APOA5-S19W (rs3135506) and -1131T/C (rs662799) variants, and the APOE polymorphism (rs429358; rs7412) using PCR and restriction analysis and TaqMan assays. We used regression analyses to examine their combined effects on TG levels (with the log-transformed variable) and the association of variant combinations with TG levels and hypertriglyceridemia (TG > or = 1.69 mmol/L), including the covariates: gender, age, waist circumference, blood glucose, blood pressure, smoking and alcohol consumption. RESULTS We found a significant lowering effect of the LPL-HindIII and S447X polymorphisms (p < 0.0001). In addition, the D9N, N291S, S19W and -1131T/C variants and the APOE-epsilon4 allele were significantly associated with an independent additive TG-raising effect (p < 0.05, p < 0.01, p < 0.001, p < 0.0001 and p < 0.001, respectively). Grouping individuals according to the presence of TG-lowering or TG-raising polymorphisms showed significant differences in TG levels (p < 0.0001), with the lowest levels exhibited by carriers of two lowering variants (10.2% reduction in TG geometric mean with respect to individuals who were homozygous for the frequent alleles of all the variants), and the highest levels in carriers of raising combinations (25.1% mean TG increase). Thus, carrying two lowering variants was protective against HTG (OR = 0.62; 95% CI, 0.39-0.98; p = 0.042) and having one single raising polymorphism (OR = 1.20; 95% CI, 1.39-2.87; p < 0.001) or more (2 or 3 raising variants; OR = 2.90; 95% CI, 1.56-5.41; p < 0.001) were associated with HTG. CONCLUSION Our results showed a significant independent additive effect on TG levels of the LPL polymorphisms HindIII, S447X, D9N and N291S; the S19W and -1131T/C variants of APOA5, and the epsilon4 allele of APOE in our study population. Moreover, some of the variant combinations studied were significantly associated with the absence or the presence of hypertriglyceridemia.

Sequence analysis of coding DNA fragments of pfcrt and pfmdr-1 genes in Plasmodium falciparum isolates from Odisha, India

The global emergence and spread of malaria parasites resistant to antimalarial drugs is the major problem in malaria control. The genetic basis of the parasite's resistance to the antimalarial drug chloroquine (CQ) is well-documented, allowing for the analysis of field isolates of malaria parasites to address evolutionary questions concerning the origin and spread of CQ-resistance. Here, we present DNA sequence analyses of both the second exon of the Plasmodium falciparum CQ-resistance transporter (pfcrt) gene and the 5' end of the P. falciparum multidrug-resistance 1 (pfmdr-1) gene in 40 P. falciparum field isolates collected from eight different localities of Odisha, India. First, we genotyped the samples for the pfcrt K76T and pfmdr-1 N86Y mutations in these two genes, which are the mutations primarily implicated in CQ-resistance. We further analyzed amino acid changes in codons 72-76 of the pfcrt haplotypes. Interestingly, both the K76T and N86Y mutations were found to co-exist in 32 out of the total 40 isolates, which were of either the CVIET or SVMNT haplotype, while the remaining eight isolates were of the CVMNK haplotype. In total, eight nonsynonymous single nucleotide polymorphisms (SNPs) were observed, six in the pfcrt gene and two in the pfmdr-1 gene. One poorly studied SNP in the pfcrt gene (A97T) was found at a high frequency in many P. falciparum samples. Using population genetics to analyze these two gene fragments, we revealed comparatively higher nucleotide diversity in the pfcrt gene than in the pfmdr-1 gene. Furthermore, linkage disequilibrium was found to be tight between closely spaced SNPs of the pfcrt gene. Finally, both the pfcrt and the pfmdr-1 genes were found to evolve under the standard neutral model of molecular evolution.

Change in mutation patterns of Plasmodium vivax dihydrofolate reductase (Pvdhfr) and dihydropteroate synthase (Pvdhps) in P. vivax isolates from malaria endemic areas of Thailand

Malaria is the most important public health problem in several countries. In Thailand, co-infections of Plasmodium vivax and Plasmodium falciparum are common. We examined the prevalence and patterns of mutations in P. vivax dihydrofolate reductase (Pvdhfr) and P. vivax dihydropteroate synthase (Pvdhps) in 103 blood samples collected from patients with P. vivax infection who had attended the malaria clinic in Mae Sot, Tak Province during 2009 and 2010. Using nested polymerase chain reaction-restriction fragment length polymorfism, we examined single nucleotide polymorphisms-haplotypes at amino acid positions 13, 33, 57, 58, 61, 117 and 173 of Pvdhfr and 383 and 553 of Pvdhps. All parasite isolates carried mutant Pvdhfr alleles, of which the most common alleles were triple mutants (99%). Eight different types of Pvdhfr and combination alleles were found, as follows: 57I/58R/117T, 57I/58R/117T, 57I/58R/117T/N, 57L/58R/117T, 57L/58R/117T, 58R/61M/117N, 58R/61M/117N and 13L/57L/58R/117T. The most common Pvdhfr alleles were 57I/58R/117T (77.7%), 57I/58R/117T/N (1%), 57L/58R/117T (5.8%) and 58R/61M/117N (14.5%). The most common Pvdhfr alleles were 57I/58R/117T (77.7%), 57I/58R/117T/N (1%), 57L/58R/117T (5.8%) and 58R/61M/117N (14.5%). Additionally, we recovered one isolate of a carrying a quadruple mutant allele, 13L/57L/58R/117T. The most prevalent Pvdhps allele was a single mutation in amino acid 383 (82.5%), followed by the wild-type A383/A553 (17.5%) allele. Results suggest that all P. vivax isolates in Thailand carry some combination of mutations in Pvdhfr and Pvdhps. Our findings demonstrate that development of new antifolate drugs effective against sulfadoxine-pyrimethamine-resistant P. vivax is required.

Rapid tests for the detection of the Mycobacterium abscessus subsp. bolletii strain responsible for an epidemic of surgical-site infections in Brazil

A single strain of Mycobacterium abscessus subsp. bolletii, characterised by a particular rpoB sequevar and two highly related pulsed field gel electrophoresis patterns has been responsible for a nationwide outbreak of surgical infections in Brazil since 2004. In this study, we developed molecular tests based on polymerase chain reaction restriction-enzyme analysis (PRA) and sequencing for the rapid identification of this strain. Sequences of 15 DNA regions conserved in mycobacteria were retrieved from GenBank or sequenced and analysed in silico. Single nucleotide polymorphisms specific to the epidemic strain and located in enzyme recognition sites were detected in rpoB, the 3' region of the 16S rDNA and gyrB. The three tests that were developed, i.e., PRA-rpoB, PRA-16S and gyrB sequence analysis, showed 100%, 100% and 92.31% sensitivity and 93.06%, 90.28% and 100% specificity, respectively, for the discrimination of the surgical strain from other M. abscessus subsp. bolletii isolates, including 116 isolates from 95 patients, one environmental isolate and two type strains. The results of the three tests were stable, as shown by results obtained for different isolates from the same patient. In conclusion, due to the clinical and epidemiological importance of this strain, these tests could be implemented in reference laboratories for the rapid preliminary diagnosis and epidemiological surveillance of this epidemic strain.

Genotyping of Mycobacterium leprae present on Ziehl-Neelsen-stained microscopic slides and in skin biopsy samples from leprosy patients in different geographic regions of Brazil

We analysed 16 variable number tandem repeats (VNTR) and three single-nucleotide polymorphisms (SNP) in Mycobacterium leprae present on 115 Ziehl-Neelsen (Z-N)-stained slides and in 51 skin biopsy samples derived from leprosy patients from Ceará (n = 23), Pernambuco (n = 41), Rio de Janeiro (n = 22) and Rondônia (RO) (n = 78). All skin biopsies yielded SNP-based genotypes, while 48 of the samples (94.1%) yielded complete VNTR genotypes. We evaluated two procedures for extracting M. leprae DNA from Z-N-stained slides: the first including Chelex and the other combining proteinase and sodium dodecyl sulfate. Of the 76 samples processed using the first procedure, 30.2% were positive for 16 or 15 VNTRs, whereas of the 39 samples processed using the second procedure, 28.2% yielded genotypes defined by at least 10 VNTRs. Combined VNTR and SNP analysis revealed large variability in genotypes, but a high prevalence of SNP genotype 4 in the Northeast Region of Brazil. Our observation of two samples from RO with an identical genotype and seven groups with similar genotypes, including four derived from residents of the same state or region, suggest a tendency to form groups according to the origin of the isolates. This study demonstrates the existence of geographically related M. leprae genotypes and that Z-N-stained slides are an alternative source for M. leprae genotyping.

Examining ERBB2 as a candidate gene for susceptibility to leprosy (Hansen’s disease) in Brazil

Leprosy remains prevalent in Brazil. ErbB2 is a receptor for leprosy bacilli entering Schwann cells, which mediates Mycobacterium leprae-induced demyelination and the ERBB2 gene lies within a leprosy susceptibility locus on chromosome 17q11-q21. To determine whether polymorphisms at the ERBB2 locus contribute to this linkage peak, three haplotype tagging single nucleotide polymorphisms (tag-SNPs) (rs2517956, rs2952156, rs1058808) were genotyped in 72 families (208 cases; 372 individuals) from the state of Pará (PA). All three tag-SNPs were associated with leprosy per se [best SNP rs2517959 odds ratio (OR) = 2.22; 95% confidence interval (CI) 1.37-3.59; p = 0.001]. Lepromatous (LL) (OR = 3.25; 95% CI 1.37-7.70; p = 0.007) and tuberculoid (TT) (OR = 1.79; 95% CI 1.04-3.05; p = 0.034) leprosy both contributed to the association, which is consistent with the previous linkage to chromosome 17q11-q21 in the population from PA and supports the functional role of ErbB2 in disease pathogenesis. To attempt to replicate these findings, six SNPs (rs2517955, rs2517956, rs1810132, rs2952156, rs1801200, rs1058808) were genotyped in a population-based sample of 570 leprosy cases and 370 controls from the state of Rio Grande do Norte (RN) and the results were analysed using logistic regression analysis. However, none of the associations were replicated in the RN sample, whether analysed for leprosy per se, LL leprosy, TT leprosy, erythema nodosum leprosum or reversal reaction conditions. The role of polymorphisms at ERBB2 in controlling susceptibility to leprosy in Brazil therefore remains unclear.

Tumour necrosis factor -308 and -238 promoter polymorphisms are predictors of a null virological response in the treatment of Brazilian hepatitis C patients

Certain host single nucleotide polymorphisms (SNPs) affect the likelihood of a sustained virological response (SVR) to treatment in subjects infected with hepatitis C virus (HCV). SNPs in the promoters of interleukin (IL)-10 (-1082 A/G, rs1800896), myxovirus resistance protein 1 (-123 C/A, rs17000900 and -88 G/T, rs2071430) and tumour necrosis factor (TNF) (-308 G/A, rs1800629 and -238 G/A, rs361525) genes and the outcome of PEGylated α-interferon plus ribavirin therapy were investigated. This analysis was performed in 114 Brazilian, HCV genotype 1-infected patients who had a SVR and in 85 non-responders and 64 relapsers. A significantly increased risk of having a null virological response was observed in patients carrying at least one A allele at positions -308 [odds ratios (OR) = 2.58, 95% confidence intervals (CI) = 1.44-4.63, p = 0.001] or -238 (OR = 7.33, 95% CI = 3.59-14.93, p < 0.001) in the TNF promoter. The risk of relapsing was also elevated (-308: OR = 2.87, 95% CI = 1.51-5.44, p = 0.001; -238: OR = 4.20, 95% CI = 1.93-9.10, p < 0.001). Multiple logistic regression of TNF diplotypes showed that patients with at least two copies of the A allele had an even higher risk of having a null virological response (OR = 16.43, 95% CI = 5.70-47.34, p < 0.001) or relapsing (OR = 6.71, 95% CI = 2.18-20.66, p = 0.001). No statistically significant association was found between the other SNPs under study and anti-HCV therapy response.

LINGO1 and LINGO2 variants are associated with essential tremor and Parkinson disease.

Genetic variation in the leucine-rich repeat and Ig domain containing 1 gene (LINGO1) was recently associated with an increased risk of developing essential tremor (ET) and Parkinson disease (PD). Herein, we performed a comprehensive study of LINGO1 and its paralog LINGO2 in ET and PD by sequencing both genes in patients (ET, n=95; PD, n=96) and by examining haplotype-tagging single-nucleotide polymorphisms (tSNPs) in a multicenter North American series of patients (ET, n=1,247; PD, n= 633) and controls (n=642). The sequencing study identified six novel coding variants in LINGO1 (p.S4C, p.V107M, p.A277T, p.R423R, p.G537A, p.D610D) and three in LINGO2 (p.D135D, p.P217P, p.V565V), however segregation analysis did not support pathogenicity. The association study employed 16 tSNPs at the LINGO1 locus and 21 at the LINGO2 locus. One variant in LINGO1 (rs9652490) displayed evidence of an association with ET (odds ratio (OR) =0.63; P=0.026) and PD (OR=0.54; P=0.016). Additionally, four other tSNPs in LINGO1 and one in LINGO2 were associated with ET and one tSNP in LINGO2 associated with PD (P<0.05). Further analysis identified one tSNP in LINGO1 and two in LINGO2 which influenced age at onset of ET and two tSNPs in LINGO1 which altered age at onset of PD (P<0.05). Our results support a role for LINGO1 and LINGO2 in determining risk for and perhaps age at onset of ET and PD. Further studies are warranted to confirm these findings and to determine the pathogenic mechanisms involved.

Mitochondrial genome sequence diversity of Indian Plasmodium falciparum isolates

We have analysed the whole mitochondrial (mt) genome sequences (each ~6 kilo nucleotide base pairs in length) of four field isolates of the malaria parasite Plasmodium falciparum collected from different locations in India. Comparative genomic analyses of mt genome sequences revealed three novel India-specific single nucleotide polymorphisms. In general, high mt genome diversity was found in Indian P. falciparum, at a level comparable to African isolates. A population phylogenetic tree placed the presently sequenced Indian P. falciparum with the global isolates, while a previously sequenced Indian isolate was an outlier. Although this preliminary study is limited to a few numbers of isolates, the data have provided fundamental evidence of the mt genome diversity and evolutionary relationships of Indian P. falciparum with that of global isolates.