Post-genomic update on a classical candidate gene for coronary artery disease: ESR1.

BACKGROUND: After age, sex is the most important risk factor for coronary artery disease (CAD). The mechanism through which women are protected from CAD is still largely unknown, but the observed sex difference suggests the involvement of the reproductive steroid hormone signaling system. Genetic association studies of the gene-encoding Estrogen Receptor α (ESR1) have shown conflicting results, although only a limited range of variation in the gene has been investigated. METHODS AND RESULTS: We exploited information made available by advanced new methods and resources in complex disease genetics to revisit the question of ESR1's role in risk of CAD. We performed a meta-analysis of 14 genome-wide association studies (CARDIoGRAM discovery analysis, N=≈87,000) to search for population-wide and sex-specific associations between CAD risk and common genetic variants throughout the coding, noncoding, and flanking regions of ESR1. In addition to samples from the MIGen (N=≈6000), WTCCC (N=≈7400), and Framingham (N=≈3700) studies, we extended this search to a larger number of common and uncommon variants by imputation into a panel of haplotypes constructed using data from the 1000 Genomes Project. Despite the widespread expression of ERα in vascular tissues, we found no evidence for involvement of common or low-frequency genetic variation throughout the ESR1 gene in modifying risk of CAD, either in the general population or as a function of sex. CONCLUSIONS: We suggest that future research on the genetic basis of sex-related differences in CAD risk should initially prioritize other genes in the reproductive steroid hormone biosynthesis system.

Genetic lineages in the yellow fever mosquito Aedes (Stegomyia) aegypti (Diptera: Culicidae) from Peru

The yellow fever mosquito Aedes aegypti was introduced in Peru in 1852 and was considered to be eradicated in 1958. In 2001, Ae. aegypti had been recorded in 15 out of 24 Peruvian Departments. Peru has great ecological differences between the east and west sides of Andes. Because of this, we consider that Ae. aegypti populations of both east and west sides can have a genetically distinct population structure. In this study we examined genetic variability and genealogical relationships among three Ae. aegypti Peruvian populations: Lima, Piura (west Andes), and Iquitos (east Andes) using a fragment of the ND4 gene of the mitochondrial genome. Three haplotypes were detected among 55 samples. Lima and Iquitos showed the same haplotype (Haplotype I), whereas Piura has two haplotypes (Haplotype II and III). Haplotype II is four mutational steps apart from Haplotype I, while Haplotype III is 13 mutational steps apart from Haplotype I in the network. The analysis of molecular variation showed that mostly of the detected genetic variation occurs at interpopulational level. The significant value phist suggests that Piura population is structured in relation to Lima and Iquitos populations and the gene flow of the ND4 is restricted in Piura when compared to Lima and Iquitos. Genetic relationship between haplotype I and haplotype II suggests introduction of the same mtDNA lineage into those localities. However the existence of a genetically distant haplotype III also suggests introduction of at least two Ae. aegypti lineages in Peru.

Genetic diversity of Triatoma infestans (Hemiptera: Reduviidae) in Chuquisaca, Bolivia based on the mitochondrial cytochrome b gene

Partial cytochrome b DNA sequences for 62 Triatoma infestans were analyzed to determine the degree of genetic variation present in populations of this insect in the northwest region of Chuquisaca, Bolivia. A total of seven haplotypes were detected in the localities sampled. The phylogenetic relationship and population genetic structure of the haplotypes found in this region, indicate that there is greater variation in this relatively small region of Bolivia than what has been previously reported by studies using the same gene fragment, for more distant geographic areas of this country. In addition, a comparison of rural and peri-urban localities, indicate that there is no difference in the genetic variation of T. infestans between these two environments.

Population genetic structure of the dengue mosquito Aedes aegypti in Venezuela

The mosquito Aedes aegypti is the main vector of dengue in Venezuela. The genetic structure of this vector was investigated in 24 samples collected from eight geographic regions separated by up to 1160 km. We examined the distribution of a 359-basepair region of the NADH dehydrogenase subunit 4 mitochondrial gene among 1144 Ae. aegypti from eight collections. This gene was amplified by the polymerase chain reaction and tested for variation using single strand conformation polymorphism analysis. Seven haplotypes were detected throughout Venezuela and these were sorted into two clades. Significant differentiation was detected among collections and these were genetically isolated by distance.

A New Large Deletion in the DFNB1 Locus Causes Nonsyndromic Hearing Loss.

Mutations in the GJB2 gene encoding the gap junction protein connexin 26 are responsible for up to 30% of all cases of autosomal recessive nonsyndromic hearing impairment (HI) with prelingual onset in most populations. The corresponding locus DFNB1, located on chromosome 13q11-q12, is also affected by three distinct deletions. These deletions extended distally to GJB2, which remains intact. We report a novel large deletion in DFNB1 observed in a patient presenting profound prelingual HI. This deletion was observed in trans to a GJB2 mutated allele carrying the p.Val84Met (V84M) mutation and was shown to be associated with hearing loss. The deletion caused a false homozygosity of V84M in the proband. Quantification of alleles by quantitative fluorescent multiplex PCR (QFM-PCR) enabled us to study the breakpoints of the deletion. The deleted segment extended through at least 920kb and removed the three connexin genes GJA3, GJB2 and GJB6. The distal breakpoint inside intron 2 of CRYL1 gene differed from the breakpoints of the known DFNB1 deletions. This case highlights the importance of screening for large deletions in molecular studies of GJB2.

Multiple refugia and barriers explain the phylogeography of the Valais shrew, Sorex antinorii (Mammalia: Soricomorpha), Sorex antinorii (Mammalia: Soricomorpha)

The aim of the present study was to investigate the genetic structure of the Valais shrew (Sorex antinorii) by a combined phylogeographical and landscape genetic approach, and thereby to infer the locations of glacial refugia and establish the influence of geographical barriers. We sequenced part of the mitochondrial cytochrome b (cyt b) gene of 179 individuals of S. antinorii sampled across the entire species' range. Six specimens attributed to S. arunchi were included in the analysis. The phylogeographical pattern was assessed by Bayesian molecular phylogenetic reconstruction, population genetic analyses, and a species distribution modelling (SDM)-based hindcasting approach. We also used landscape genetics (including isolation-by-resistance) to infer the determinants of current intra-specific genetic structure. The phylogeographical analysis revealed shallow divergence among haplotypes and no clear substructure within S. antinorii. The starlike structure of the median-joining network is consistent with population expansion from a single refugium, probably located in the Apennines. Long branches observed on the same network also suggest that another refugium may have existed in the north-eastern part of Italy. This result is consistent with SDM, which also suggests several habitable areas for S. antinorii in the Italian peninsula during the LGM. Therefore S. antinorii appears to have occupied disconnected glacial refugia in the Italian peninsula, supporting previous data for other species showing multiple refugia within southern refugial areas. By coupling genetic analyses and SDM, we were able to infer how past climatic suitability contributed to genetic divergence of populations. The genetic differentiation shown in the present study does not support the specific status of S. arunchi.

Genetic variability of Aedes aegypti in the Americas using a mitochondrial gene: evidence of multiple introductions

To analyze the genetic relatedness and phylogeographic structure of Aedes aegypti, we collected samples from 36 localities throughout the Americas (Brazil, Peru, Venezuela, Guatemala, US), three from Africa (Guinea, Senegal, Uganda), and three from Asia (Singapore, Cambodia, Tahiti). Amplification and sequencing of a fragment of the mitochondrial NADH dehydrogenase subunit 4 gene identified 20 distinct haplotypes, of which 14 are exclusive to the Americas, four to African/Asian countries, one is common to the Americas and Africa, and one to the Americas and Asia. Nested clade analysis (NCA), pairwise distribution, statistical parsimony, and maximum parsimony analyses were used to infer evolutionary and historic processes, and to estimate phylogenetic relationships among haplotypes. Two clusters were found in all the analyses. Haplotypes clustered in the two clades were separated by eight mutational steps. Phylogeographic structure detected by the NCA was consistent with distant colonization within one clade and fragmentation followed by range expansion via long distance dispersal in the other. Three percent of nucleotide divergence between these two clades is suggestive of a gene pool division that may support the hypothesis of occurrence of two subspecies of Ae. aegypti in the Americas.

Population structure of the malaria vector Anopheles darlingi in Rondônia, Brazilian Amazon, based on mitochondrial DNA

Anopheles darlingi is the most important Brazilian malaria vector, with a widespread distribution in the Amazon forest. Effective strategies for vector control could be better developed through knowledge of its genetic structure and gene flow among populations, to assess the vector diversity and competence in transmitting Plasmodium. The aim of this study was to assess the genetic diversity of An. darlingi collected at four locations in Porto Velho, by sequencing a fragment of the ND4 mitochondrial gene. From 218 individual mosquitoes, we obtained 20 different haplotypes with a diversity index of 0.756, equivalent to that found in other neotropical anophelines. The analysis did not demonstrate significant population structure. However, haplotype diversity within some populations seems to be over-represented, suggesting the presence of sub-populations, but the presence of highly represented haplotypes complicates this analysis. There was no clear correlation among genetic and geographical distance and there were differences in relation to seasonality, which is important for malarial epidemiology.

Polymorphism in the promoter region of the mannose-binding lectin gene among human T-cell lymphotropic virus infected subjects

The present study investigated the frequency of the mutations at positions -550 and -221 of the mannose-binding lectin (MBL) gene in a sample of 75 human T-cell lymphotropic virus (HTLV) infected patients and 96 HTLV seronegative controls, in order to evaluate the occurrence of a possible association between the polymorphism and HTLV infection. A sequence specific primer-polymerase chain reaction was used for discrimination of the polymorphism. The analysis of allele frequencies at position -550 did not show any significant differences between HTLV infected group and controls, but there was a significant difference at position -221. The comparative analysis of haplotypes frequencies were not significant, but the genotype frequencies between the two groups, revealed a higher prevalence of genotype LYLX (25.3%), associated with medium and low MBL serum levels among HTLV infected subjects. The odds ratio estimation demonstrated that the presence of genotype LYLX was associated with an increased risk of HTLV infection (p = 0.0096; 1.38 < IC95% < 7.7605). There was no association between proviral load and the promoter polymorphism, but when promoter and exon 1 mutations were matched, it was possible to identify a significant higher proviral load among HTLV infected individuals carrying haplotypes correlated to low serum levels of MBL. The present study shows that the polymorphism in the promoter region of the MBL gene may be a genetic marker associated with HTLV infection, and emphasizes the need for further studies to determinate if the present polymorphism have any impact on diseases linked to HTLV infection.

Characterization of polymorphisms in the mannose-binding lectin gene promoter among human immunodeficiency virus 1 infected subjects

The present study investigated the prevalence of mutations in the -550 (H/L) and -221 (X/Y) mannose-binding lectin (MBL) gene promoter regions and their impact on infection by human immunodeficiency virus 1 (HIV-1) in a population of 128 HIV-1 seropositive and 97 seronegative patients. The allele identification was performed through the sequence-specific primer polymerase chain reaction method, using primer sequences specific to each polymorphism. The evolution of the infection was evaluated through CD4+ T-lymphocyte counts and plasma viral load. The allele and haplotype frequencies among HIV-1-infected patients and seronegative healthy control patients did not show significant differences. CD4+ T-lymphocyte counts showed lower levels among seropositive patients carrying haplotypes LY, LX and HX, as compared to those carrying the HY haplotype. Mean plasma viral load was higher among seropositive patients with haplotypes LY, LX and HX than among those carrying the HY haplotype. When promoter and exon 1 mutations were matched, it was possible to identify a significantly higher viral load among HIV-1 infected individuals carrying haplotypes correlated to low serum levels of MBL. The current study shows that haplotypes related to medium and low MBL serum levels might directly influence the evolution of viral progression in patients. Therefore, it is suggested that the identification of haplotypes within the promoter region of the MBL gene among HIV-1 infected persons should be further evaluated as a prognostic tool for AIDS progression.

Population structure and diversity of the pathogenic fungus Aspergillus fumigatus isolated from different sources and geographic origins

Fifty-five clinical and environmental Aspergillus fumigatus isolates from Mexico, Argentina, France and Peru were analyzed to determine their genetic variability, reproductive system and level of differentiation using amplified fragment length polymorphism markers. The level of genetic variability was assessed by measuring the percentage of polymorphic loci, number of effective alleles, expected heterozygocity and by performing an association index test (I A). The degree of genetic differentiation and variation was determined using analysis of molecular variance at three levels. Using the paired genetic distances, a dendrogram was built to detect the genetic relationship among alleles. Finally, a network of haplotypes was constructed to determine the geographic relationship among them. The results indicate that the clinical isolates have greater genetic variability than the environmental isolates. The I A of the clinical and environmental isolates suggests a recombining population structure. The genetic differentiation among isolates and the dendrogram suggest that the groups of isolates are different. The network of haplotypes demonstrates that the majority of the isolates are grouped according to geographic origin.

Multiple sclerosis risk variant HLA-DRB1*1501 associates with high expression of DRB1 gene in different human populations.

The human leukocyte antigen (HLA) DRB1*1501 has been consistently associated with multiple sclerosis (MS) in nearly all populations tested. This points to a specific antigen presentation as the pathogenic mechanism though this does not fully explain the disease association. The identification of expression quantitative trait loci (eQTL) for genes in the HLA locus poses the question of the role of gene expression in MS susceptibility. We analyzed the eQTLs in the HLA region with respect to MS-associated HLA-variants obtained from genome-wide association studies (GWAS). We found that the Tag of DRB1*1501, rs3135388 A allele, correlated with high expression of DRB1, DRB5 and DQB1 genes in a Caucasian population. In quantitative terms, the MS-risk AA genotype carriers of rs3135388 were associated with 15.7-, 5.2- and 8.3-fold higher expression of DQB1, DRB5 and DRB1, respectively, than the non-risk GG carriers. The haplotype analysis of expression-associated variants in a Spanish MS cohort revealed that high expression of DRB1 and DQB1 alone did not contribute to the disease. However, in Caucasian, Asian and African American populations, the DRB1*1501 allele was always highly expressed. In other immune related diseases such as type 1 diabetes, inflammatory bowel disease, ulcerative colitis, asthma and IgA deficiency, the best GWAS-associated HLA SNPs were also eQTLs for different HLA Class II genes. Our data suggest that the DR/DQ expression levels, together with specific structural properties of alleles, seem to be the causal effect in MS and in other immunopathologies rather than specific antigen presentation alone.

Characterisation of pvmdr1 and pvdhfr genes associated with chemoresistance in Brazilian Plasmodium vivax isolates

Plasmodium vivax control is now being hampered by drug resistance. Orthologous Plasmodium falciparum genes linked to chloroquine or sulfadoxine-pyrimethamine chemoresistance have been identified in P. vivax parasites, but few studies have been performed. The goal of the present work is to characterise pvmdr1 and pvdhfr genes in parasite isolates from a Brazilian endemic area where no molecular investigation had been previously conducted. The pvmdr1 analysis revealed the existence of single (85.7%) and double (14.3%) mutant haplotypes, while the pvdhfr examination showed the presence of double (57.2%) and triple (42.8%) mutant haplotypes. The implications of these findings are discussed.

Trypanosoma cruzi: ancestral genomes and population structure

Although the genome of Trypanosoma cruzi has been completely sequenced, little is known about its population structure and evolution. Since 1999, two major evolutionary lineages presenting distinct epidemiological characteristics have been recognised: T. cruzi I and T. cruzi II. We describe new and important aspects of the population structure of the parasite, and unequivocally characterise a third ancestral lineage that we propose to name T. cruzi III. Through a careful analysis of haplotypes (blocks of genes that are stably transmitted from generation to generation of the parasite), we inferred at least two hybridisation events between the parental lineages T. cruzi II and T. cruzi III. The strain CL Brener, whose genome was sequenced, is one such hybrid. Based on these results, we propose a simple evolutionary model based on three ancestral genomes, T. cruzi I, T. cruzi II and T. cruzi III. At least two hybridisation events produced evolutionarily viable progeny, and T. cruzi III was the cytoplasmic donor for the resulting offspring (as identified by the mitochondrial clade of the hybrid strains) in both events. This model should be useful to inform evolutionary and pathogenetic hypotheses regarding T. cruzi.

On the use of principal components in contemporany population genetics: a case study

In human Population Genetics, routine applications of principal component techniques are oftenrequired. Population biologists make widespread use of certain discrete classifications of humansamples into haplotypes, the monophyletic units of phylogenetic trees constructed from severalsingle nucleotide bimorphisms hierarchically ordered. Compositional frequencies of the haplotypesare recorded within the different samples. Principal component techniques are then required as adimension-reducing strategy to bring the dimension of the problem to a manageable level, say two,to allow for graphical analysis.Population biologists at large are not aware of the special features of compositional data and normally make use of the crude covariance of compositional relative frequencies to construct principalcomponents. In this short note we present our experience with using traditional linear principalcomponents or compositional principal components based on logratios, with reference to a specificdataset