Modeling associations between genetic markers using Bayesian networks

Motivation: Understanding the patterns of association between polymorphisms at different loci in a population ( linkage disequilibrium, LD) is of fundamental importance in various genetic studies. Many coefficients were proposed for measuring the degree of LD, but they provide only a static view of the current LD structure. Generative models (GMs) were proposed to go beyond these measures, giving not only a description of the actual LD structure but also a tool to help understanding the process that generated such structure. GMs based in coalescent theory have been the most appealing because they link LD to evolutionary factors. Nevertheless, the inference and parameter estimation of such models is still computationally challenging. Results: We present a more practical method to build GM that describe LD. The method is based on learning weighted Bayesian network structures from haplotype data, extracting equivalence structure classes and using them to model LD. The results obtained in public data from the HapMap database showed that the method is a promising tool for modeling LD. The associations represented by the learned models are correlated with the traditional measure of LD D`. The method was able to represent LD blocks found by standard tools. The granularity of the association blocks and the readability of the models can be controlled in the method. The results suggest that the causality information gained by our method can be useful to tell about the conservability of the genetic markers and to guide the selection of subset of representative markers.

Human leukocyte antigen (HLA) and single nucleotide polymorphisms (SNPs) tumor necrosis factor (TNF)-alpha-238 and-308 as genetic markers of susceptibility to psoriasis and severity of the disease in a long-term follow-up Brazilian study

Background The strongest genetic marker for psoriasis is Cw*06. Polymorphisms in the tumor necrosis factor (TNF)-alpha promoter region, especially replacement of guanine with adenine in positions -238 and -308 are related to higher TNF-alpha production and higher risk for psoriasis in Caucasoid populations, not found in Asians. We performed a case-control study of 69 patients with psoriasis type I and 70 controls, characterized clinical progression along 10-years of follow-up in mild or severe disease and determined HLA class I, II, and TNF single nucleotide polymorphisms (SNPs) -238 and -308 polymorphisms to demonstrate whether these polymorphisms may be genetic risk for susceptibility to psoriasis or severity of the disease in Brazilians. Methods Polymorphisms were identified using PCR/SSP. Alleles, genotypes, and haplotypes frequencies were compared using Fisher`s test. Results More severe disease was found in male patients. It may be suggested that alleles B*37, Cw*06, Cw*12, and DRB1*07 were associated with severe disease course, while B*57 with mild disease. No statistical difference was found between the patients and controls regarding polymorphisms frequencies in TNF SNPs. This study pointed to a higher TNF-238 G/G genotype frequency (OR: 3.21; CI: 1.06-9.71; P = 0.04) in the group with severe disease. Conclusions Polymorphisms in the TNF-alpha SNPs do not seem to be a more important genetic risk factor for psoriasis than the already known Cw*06 in Brazilian patients, but these markers may be related to clinical manifestations.

Systemic lupus erythematosus: Association with KIR and SLC11A1 polymorphisms, ethnic predisposition and influence in clinical manifestations at onset revealed by ancestry genetic markers in an urban Brazilian population

Systemic lupus erythematosus (SLE) is an autoimmune disorder of the connective tissue with a wide and heterogeneous spectrum of manifestations, with renal and neurological involvement usually related to worse prognosis. SLE more frequently affects females of reproductive age, and a high prevalence and renal manifestation seem to be associated with non-European ethnicity. The present study aims to investigate candidate loci to SLE predisposition and evaluate the influence of ethnic ancestry in the disease risk and clinical phenotypic heterogeneity of lupus at onset. Samples represented by 111 patients and 345 controls, originated from the city of Belem, located in the Northern Region of Brazil, were investigated for polymorphisms in HLA-G, HLA-C, SLC11A1, MTHFR, CASP8 and 15 KIR genes, in addition to 89 Amerindian samples genotyped for SLC11A1. We also investigated 48 insertion/deletion ancestry markers to characterize individual African, European and Amerindian ancestry proportions in the samples. Predisposition to SLE was associated with GTGT deletion at the SLC11A1 3`UTR, presence of KIR2DS2 +/KIR2DS5 +/KIR3DS1 + profile, increased number of stimulatory KIR genes, and European and Amerindian ancestries. The ancestry analysis ruled out ethnic differences between controls and patients as the source of the observed associations. Moreover, the African ancestry was associated with renal manifestations. Lupus (2011) 20, 265-273.

Testing the applicability of molecular genetic markers to population analyses of scleratinian corals

The abundance of coral reefs worldwide is in decline, and despite the ecological importance of reefs, only a limited number of DNA markers have been identified for scleractinian coral genetic studies. This paper addresses the search for new coral molecular markers and investigates the applicability of the cytochrome c oxidase subunit I (COI), the internal transcribed spacer region 1 (ITS1), and the pocilloporin gene to the question of intraspecific variation in the scleractinian coral Pocillopora verrucosa along the southeast African coastline. The COI fragment was 710 bp long and was identical for P. verrucosa (n = 10) and P. damicornis (n = 3). Only two different ITS1 sequences were found (differing by 13 bp insertion), but more importantly, 24% of the sequences were heterogenous indicating that different multiple copies of the sequence exist. Pocilloporin is an intronless gene that was absolutely conserved throughout all P. verrucosa populations (n = 50). Thus, the three DNA regions studied appear unsuitable for the population genetic analyses of P. verrucosa.

Genetic markers from Biomphalaria tenagophila (Gastropoda: Pulmonata: Planorbidae) obtained by the double stringency polymerase chain reaction technique

Biomphalaria tenagophila, one of the intermediate hosts of the trematoda Schistosoma mansoni, is a simultaneous hermafrodite snail species. In order to analyse the genetic structure of these populations, we performed a double-stringency PCR technique to obtain genetic markers with microsatellites and arbitrary primers in a single reaction.

Genetic markers between Biomphalaria glabrata snails susceptible and resistant to Schistosoma mansoni infection

The analysis of the genetic variability related to susceptibility to Schistosoma mansoni infection in the vector of the genus Biomphalaria is important in terms of a better understanding of the epidemiology of schistosomiasis itself, the possible pathological implications of this interaction in vertebrate hosts, and the formulation of new strategies and approaches for disease control. In the present study, the genetic variability of B. glabrata strains found to be resistant or susceptible to S. mansoni infection was investigated using DNA amplification by random amplified polymorphic DNA-polymerase chain reaction (RAPD-PCR). The amplification products were analyzed on 8% polyacrylamide gel and stained with silver. We selected 10 primers, since they have previously been useful to detect polymorphism among B. glabrata and/or B. tenagophila. The results showed polymorphisms with 5 primers. Polymorphic bands observed only in the susceptible strain. The RAPD-PCR methodology represents an adequate approach for the analysis of genetic polymorphisms. The understanding of the genetic polymorphisms associated to resistance may contribute to the future identification of genomic sequences related to the resistance/susceptibility of Biomphalaria to the larval forms of S. mansoni and to the development of new strategies for the control of schistosomiasis.

The association of genetic markers and malaria infection in the Brazilian Western Amazonian region

Almost all individuals (182) belonging to an Amazonian riverine population (Portuchuelo, RO, Brazil) were investigated for ascertaining data on epidemiological aspects of malaria. Thirteen genetic blood polymorphisms were investigated (ABO, MNSs, Rh, Kell, and Duffy systems, haptoglobins, hemoglobins, and the enzymes glucose-6-phosphate dehydrogenase, glyoxalase, phosphoglucomutase, carbonic anhydrase, red cell acid phosphatase, and esterase D). The results indicated that the Duffy system is associated with susceptibility to malaria, as observed in other endemic areas. Moreover, suggestions also arose indicating that the EsD and Rh loci may be significantly associated with resistance to malaria. If statistical type II errors and sample stratification could be ruled out, hypotheses on the existence of a causal mechanism or an unknown closely linked locus involved in susceptibility to malaria infection may explain the present findings.

Tests for sex-biased dispersal using bi-parentally inherited genetic markers.

Understanding why dispersal is sex-biased in many taxa is still a major concern in evolutionary ecology. Dispersal tends to be male-biased in mammals and female-biased in birds, but counter-examples exist and little is known about sex bias in other taxa. Obtaining accurate measures of dispersal in the field remains a problem. Here we describe and compare several methods for detecting sex-biased dispersal using bi-parentally inherited, codominant genetic markers. If gene flow is restricted among populations, then the genotype of an individual tells something about its origin. Provided that dispersal occurs at the juvenile stage and that sampling is carried out on adults, genotypes sampled from the dispersing sex should on average be less likely (compared to genotypes from the philopatric sex) in the population in which they were sampled. The dispersing sex should be less genetically structured and should present a larger heterozygote deficit. In this study we use computer simulations and a permutation test on four statistics to investigate the conditions under which sex-biased dispersal can be detected. Two tests emerge as fairly powerful. We present results concerning the optimal sampling strategy (varying number of samples, individuals, loci per individual and level of polymorphism) under different amounts of dispersal for each sex. These tests for biases in dispersal are also appropriate for any attribute (e.g. size, colour, status) suspected to influence the probability of dispersal. A windows program carrying out these tests can be freely downloaded from http://www.unil.ch/izea/softwares/fstat.html

Distinct genetic structure in populations of Chrysoperla externa (Hagen) (Neuroptera, Chrysopidae) shown by genetic markers ISSR and COI gene

Distinct genetic structure in populations of Chrysoperla externa (Hagen) (Neuroptera, Chrysopidae) shown by genetic markers ISSR and COI gene. Green lacewings are generalist predators, and the species Chrysoperla externa presents a great potential for use in biological control of agricultural pests due to its high predation and reproduction capacities, as well as its easy mass rearing in the laboratory. The adaptive success of a species is related to genetic variability, so that population genetic studies are extremely important in order to maximize success of the biological control. Thus, the present study used nuclear (Inter Simple Sequence Repeat - ISSR) and mitochondrial (Cytochrome Oxidase I - COI) molecular markers to estimate the genetic variability of 12 populations in the São Paulo State, Brazil, as well as the genetic relationships between populations. High levels of genetic diversity were observed for both markers, and the highest values of genetic diversity appear associated with municipalities that have the greatest areas of native vegetation. There was high haplotype sharing, and there was no correlation between the markers and the geographic distribution of the populations. The AMOVA indicated absence of genetic structure for the COI gene, suggesting that the sampled areas formed a single population unit. However, the great genetic differentiation among populations showed by ISSR demonstrates that these have been under differentiation after their expansion or may also reflect distinct dispersal behavior between males and females.

Sex-linked genetic markers in the shrews of the Sorex araneus group : from free gene flow to speciation

SUMMARY: The shrews of the Sorex araneus group are morphologically .very similar, but have undergone a spectacular chromosomal evolution. Altogether, the shrews of this group present a complete array of every possible level of chromosomal and genetic differentiation. In South-Western Europe, four species are recognised: S. antiriorii, S. araneus, S. coronatus and S. granarius, which differ essentially by the amount and the composition of Robertsonian metacentric chromosomés. Additionally, several chromosome races of S. araneus are also present in the same region (i.e. Bretolet, Carlit, Cordon, Jura and Vaud). The objective of this thesis was to examine the genetic relationships between populations, races and /or species of the Sorex araneus group with a special emphasis onsex-specific markers (mtDNA and Y chromosome). We first investigate the evolutionary history of the shrews of the Sorex araneus group distributed in the South-Western Europe. The results of. these analyses confirmed the difficulty to draw a single dichotomic tree within this group. Incongruent mtDNA and Y chromosome phylogenies suggest further that genetic and chromosomal evolution are in this group partially independent processes and that the evolutionary history of the south-western European populations of the S. araneus group can only be understood if we consider secondary contacts between taxa, after their divergence (with genetic exchanges by means of hybridization and / or introgression). Using one male-inherited, one female inherited and eight biparentally inherited markers, we investigate the population genetic structure of the Valais shrew (Sorex antinorii). Overall there results suggest that two already well-differentiated genetic lineages colonized the Swiss Alps after the last glacial period and came into contact in the Rhône Valley. After the Valais shrew (Sorex antinorii) reached the Swiss Alps, it came into contact with the common shrew (Sorex araneus). When two species come into contact and hybridize, endogenous counter-selection of hybrids is usually first expressed as a reduced fertility or viability in hybrids of the heterogametic sex, a mechanism know as Haldane's rule (Haldane 1922). We first evaluated the extent of introgression for Y chromosome, mtDNA and autosomal markers in a hybrid zone between S. antinoriii and S. araneus. The overall level of genetic and karyotypic differentiation between the two species must be strong .enough to allow the detection asymmetric introgression. Secondly, we compared the levels of gene flow between chromosome common to both species and chromosome differently rearranged in each of them. We detected a significantly stronger genetic structure in rearranged chromosomes. Over a 10-year period, we even observed a decrease of genetic structure for common chromosomes. These results strongly support the role of chromosomal rearrangements in the reproductive barrier between S. araneus and S. anfinorii. Overall, this thesis underlines the need to use different inherited (paternally, maternally and / or biparentally) and chromosomally located (on common vs. on rearranged chromosomes) markers to obtain more accurate pictures of genetic relationships between populations or species. RÉSUMÉ: Les musaraignes du groupe Sorex araneus sont morphologiquement très proches, mais ont connu une spectaculaire évolution chromosomique. Prises dans leur ensemble, les musaraignes de ce groupe présentent tous les nivaux possibles de différenciation génétique et chromosomique. Dans le sud-ouest de l'Europe, quatre espèces appartenant à ce groupe sont présentes : S. antinorii, S. araneus, S. coronatus et S. granarius. Celles-ci diffèrent essentiellement par leur caryotype dont la variabilité est principalement due à des fusions Robertsoniennes. De plus, plusieurs races chromosomiques appartenant à S. araneus sont aussi présentes dans la même région (i.e. les races Bretolet, Carlit, Cordon, Jura et Vaud). L'objectif de cette thèse était d'examiner les relations génétiques entre populations, races et/ou espèces du groupe S. araneus, en utilisant particulièrement des marqueurs liés aux sexes (ADN mitochondrial et Chromosome Y). Nous avons dans un premier temps retracé l'histoire évolutive des musaraignes de ce groupe dans le sud-ouest de l'Europe. Les résultats dé ces analyses confirment qu'il est difficile de tracer un simple arbre dichotomique au sein de ce groupe. Les arbres phylogénétiques obtenus sur l'ADN mitochondrial et le chromosome Y sont incongruents et suggèrent de plus que l'évolution génétique et chromosomique sont des processus indépendants. L'histoire évolutive -des populations de ce groupe ne peut. être comprise qu'en considérant des contacts secondaires entre taxa postérieure à leur divergence et induisant des échanges génétiques par hybridation et/ou introgression. Par la suite, nous avons examiné la structure génétique des populations de la musaraigne du Valais, S. antinorii, en utilisant un marqueur transmis par les mâles, un marqueur transmis par les femelles et huit marqueurs transmis par les 2 sexes. Nos résultats suggèrent que deux lignées génétiquement bien différenciées aient colonisé les Alpes Suisses, après les dernières glaciations et entrent en contact dans là Vallée du Rhône. Après avoir franchi les Alpes Suisses, la musaraigne du Valais est entrée en contact avec là musaraigne commune (S. araneus). Lorsque deux espèces entrent en contact et s'hybrident, la sélection contre les hybrides implique habituellement une baisse de fertilité ou de viabilité des hybrides du sexe hétérogamétique (i.e. les mâles XY chez les mammifères). Ce mécanisme est connu sous le nom de règle de Haldane (Haldane 1922) et implique une plus forte structuration génétique de marqueurs males - spécifiques que des marqueurs femelles spécifiques. Nous avons donc évalué le degré d'introgression des marqueurs situés sur le chromosome Y, sur l'ADN mitochondrial et sur des autosomes dans une zone hybride entre S. araneus et S. antinorii. Le niveau de différenciation chromosomique et génétique entre les 2 espèces doit être suffisamment fort pour ne pas permettre la détection d'une introgression asymétrique entre les sexes. Dans un second temps, nous avons comparé les niveaux de flux de gênes mesurés à l'échelle du chromosome, pour des chromosomes communs aux deux espèces et pour des chromosomes différemment arrangées dans chacune des deux espèces. Nous avons détecté une structure génétique significativement plus forte sur les chromosomes réarrangés et comme la zone hybride a été étudiée à dix années d'intervalle, nous observons même une diminution de la structure génétique pour les chromosomes communs au cours du temps.. Ces résultats soutiennent fortement l'hypothèse d'un rôle des réarrangements chromosomiques dans l'établissement d'une barrière reproductive entre S. araneus et S. antinorii. Ainsi cette thèse souligne l'utilité d'utiliser des marqueurs génétiques avec différents modes de transmission. (par les mâles, par les femelles et/ou par les 2 sexes) ou localisés au niveau du chromosome (chromosomes communs vs chromosomes réarrangés) afin d'obtenir une image plus juste ou du moins plus complète des relations génétiques entre populations ou espèces.

Association study of lipoprotein(a) genetic markers, traditional risk factors, and coronary heart disease in HIV-1-infected patients.

Objectives: General population studies have shown associations between copy number variation (CNV) of the LPA gene Kringle-IV type-2 (KIV-2) coding region, single-nucleotide polymorphism (SNP) rs6415084 in LPA and coronary heart disease (CHD). Because risk factors for HIV-infected patients may differ from the general population, we aimed to assess whether these potential associations also occur in HIV-infected patients. Methods: A unicenter, retrospective, case-control (1:3) study. Eighteen HIV-patients with confirmed diagnosis of acute myocardial infarction (AMI) were adjusted for age, gender, and time since HIV diagnosis to 54 HIV-patients without CHD. After gDNA extraction from frozen blood, both CNV and SNP genotyping were performed using real-time quantitative PCR. All genetic and non-genetic variables for AMI were assessed in a logistic regression analysis. Results: Our results did not confirm any association in terms of lipoprotein(a) LPA structural genetic variants when comparing KIV-2 CNV (p = 0.67) and SNP genotypes (p = 0.44) between AMI cases and controls. However, traditional risk factors such as diabetes mellitus, hypertension, and CD4(+) T cell count showed association (p < 0.05) with CHD. Conclusion: Although significant associations of AMI with diabetes, hypertension and CD4(+) T cell count in HIV-patients were found, this study could not confirm the feasibility neither of KIV-2 CNV nor rs6415084 in LPA as genetic markers of CHD in HIV-infected patients.Highlights:● Individuals with HIV infection are at higher risk of coronary heart disease (CHD) than the non-infected population.● Our results showed no evidence of LPA structural genetic variants associated with CHD in HIV-1-infected patients.● Associations were found between diabetes mellitus, arterial hypertension, CD4(+) T cell count, and CHD.● The clinical usefulness of these biomarkers to predict CHD in HIV-1-infected population remains unproven.● Further studies are needed to assess the contribution of common genetic variations to CHD in HIV-infected individuals.

Genetic markers of the Va/Ba balanced lethal strain of Drosophila subobscura

The Va/Ba strain, constructed by Sperlich et al. (1977), is the only balanced lethal strain in D. subobscura. It allows the production of homozygous O chromosomes and has been a useful tool not only to analyse chromosomal viabilities but also to obtain homokaryotypic lines (Mestres and Serra, 2008). Besides the morphological dominant mutations Va (Varicose) and Ba (Bare), other genetic markers have been characterized in this strain, some of them by our group and not described previously. Here we present a list of these markers.

Clinical and genetic markers in disease progression – a study among subjects with obstructive lung diseases

Asthma, COPD, and asthma and COPD overlap syndrome (ACOS) are chronic pulmonary diseases with an obstructive component. In COPD, the obstruction is irreversible and the disease is progressive. The aim of the study was to define and analyze factors that affected disease progression and patients’ well-being, prognosis and mortality in Chronic Airway Disease (CAD) cohort. The main focus was on COPD and ACOS patients. Retrospective data from medical records was combined with genetic and prospective follow-up data. Smoking is the biggest risk factor for COPD and even after the diagnosis of the disease, smoking plays an important role in disease development and patient’s prognosis. Sixty percent of the COPD patients had succeeded in smoking cessation. Patients who had managed to quit smoking had lower mortality rates and less psychiatric diseases and alcohol abuse although they were older and had more cardiovascular diseases than patients who continued smoking. Genetic polymorphism rs1051730 in the nicotinic acethylcholine receptor gene (CHRNA3/5) associated with heavy smoking, cancer prevalence and mortality in two Finnish independent cohorts consisting of COPD patients and male smokers. Challenges in smoking cessation and higher mortality rates may be partly due to individual patient’s genetic composition. Approximately 50% of COPD patients are physically inactive and the proportion was higher among current smokers. Physically active and inactive patients didn’t differ from each other in regard to age, gender or comorbidities. Bronchial obstruction explained inactivity only in severe disease. Subjective sensation of dyspnea, however, had very strong association to inactivity and was also associated to low health related quality of life (HRQoL). ACOS patients had a significantly lower HRQoL than either the patients with asthma or with COPD even though they were younger than COPD patients, had better lung functions and smaller tobacco exposure.