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.

Molecular markers and genetic diversity of Plasmodium vivax

Enhanced understanding of the transmission dynamics and population genetics for Plasmodium vivax is crucial in predicting the emergence and spread of novel parasite phenotypes with major public health implications, such as new relapsing patterns, drug resistance and increased virulence. Suitable molecular markers are required for these population genetic studies. Here, we focus on two groups of molecular markers that are commonly used to analyse natural populations of P. vivax. We use markers under selective pressure, for instance, antigen-coding polymorphic genes, and markers that are not under strong natural selection, such as most minisatellite and microsatellite loci. First, we review data obtained using genes encoding for P. vivax antigens: circumsporozoite protein, merozoite surface proteins 1 and 3α, apical membrane antigen 1 and Duffy binding antigen. We next address neutral or nearly neutral molecular markers, especially microsatellite loci, providing a complete list of markers that have already been used in P. vivax populations studies. We also analyse the microsatellite loci identified in the P. vivax genome project. Finally, we discuss some practical uses for P. vivax genotyping, for example, detecting multiple-clone infections and tracking the geographic origin of isolates.

Nuclear and mitochondrial DNA markers in traceability of retail beef samples

Several characteristics are important in a traceability system of animal products, such as age at slaughter, breed composition, besides information of the productive chain. In general, the certification agent records information about the animals and the system which it came from, although cannot guarantee that the slaughtering, meat processing and distribution are error proof. Besides, there is a differential price, at least at the international market, based on sex and breed composition of the animals. Genetic markers allow identification of characteristics controlled in the beef cattle traceability program, as sex and breed composition, in order to correctly identify and appraise the final product for the consumer. The hypothesis of this study was that the majority beef samples retailed in the local market originate from female with a great participation of zebu breeds. Therefore, the objective of this work was to characterize retail beef samples with DNA markers that identify cattle sex and breed composition. Within 10 beef shops localized in Pirassununga, SP, Brazil, 61 samples were collected, all were genotyped as harboring Bos taurus mitochondrial DNA and 18 were positive for the Y chromosome amplification (male). For the marker sat1711b-Msp I the frequency of the allele A was 0.278 and for the marker Lhr-Hha I the frequency of the allele T was 0.417. The results of sat1711b-Msp I and Lhr-Hha I allelic frequencies are suggestive that the proportion of indicus genome compared with the taurine genome in the market meat is smaller than the observed in the Nellore breed. The procedure described in this study identified sex and subspecies characteristics of beef meat samples, with potential application in meat products certification in special as an auxiliary tool in beef cattle traceability programs.

PRACTICAL USE OF MOLECULAR MARKERS IN BEEF AND DAIRY ZEBU CATTLE IN BRAZIL

Selection of zebu (Bos indicus) beef and dairy cattle in Brazil and the validation process of genetic markers for growth, carcass and meat quality traits and also for milk production, fat and protein milk content are discussed as concerned to the concepts and details of their use as auxiliary tools in selection processes. It is highlighted, also, the importance of right selection of ova donor cows for production of embryos to be transferred.

Migration in Afro-Brazilian Rural Communities: Crossing Demographic and Genetic Data

Many studies have used genetic markers to understand global migration patterns of our species. However, there are only few studies of human migration on a local scale. We, therefore, researched migration dynamics in three Afro-Brazilian rural communities, using demographic data and ten Ancestry Informative Markers. In addition to the description of migration and marriage structures, we carried out genetic comparisons between the three populations, as well as between locals and migrants from each community. Genetic admixture analyses were conducted according to the gene-identity method, with Sub-Saharan Africans, Amerindians, and Europeans as parental populations. The three analyzed Afro-Brazilian rural communities consisted of 16% to 30% of migrants, most of them women. The age pyramid revealed a gap in the segment of men aged between 20 to 30 yrs. While endogamous marriages predominated, exogamous marriages were mainly patrilocal. Migration dynamics are apparently associated with matrimonial customs and other social practices of such communities. The impact of migration upon the populations` genetic composition was low but showed an increase in European alleles with a concomitant decrease in the Amerindian contribution. Admixture analysis evidenced a higher African contribution to the gene pool of the studied populations, followed by the contribution of Europeans and Amerindians, respectively.

Genetic diversity among Toxoplasma gondii isolates of small ruminants from Brazil: Novel genotypes revealed

Recent studies indicated that Toxoplasma gondii isolates of many domestic animal hosts from Brazil are genetically and biologically different from those in USA and Europe. Despite of high pathogenicity of this parasite to small ruminants, the epidemiology and genetic diversity of T. gondii in these animals are not well understood in Brazil. In this study, a total of 28 T. gondii samples (16 isolates from sheep in Sao Paulo state, and 12 isolates from goats in the states of Sao Paulo and Rio Grande do Norte) were genotyped using genetic markers SAG1, SAG2, SAG3, BTUB, CRAG, c22-8, c29-2, L358, PK1, Apico and CS3. Eleven genotypes were identified from these T. gondii isolates. Eight isolates (4 from sheep and 4 from goats) were grouped into the common clonal type Brl lineage. One sheep isolate was grouped to the type BrIII lineage. Five isolates grouped to three previously identified genotypes in Brazil, and 13 isolates grouped to six novel genotypes. Mixed genotype was found in one isolate from goat in Sao Paulo. No classical clonal Type I. II or III isolates were found, confirming previous reports that these clonal lineages are rare in Brazil. The allele types at the CS3 locus are strongly linked to mouse virulence of the parasite. The results of this study indicate that even though a large number of T. gondii genotypes have been identified from a variety of animal hosts in Brazil, high percentage of new genotypes are continuously identified from different animal species, suggesting extremely high diversity of T. gondii in the population. (C) 2010 Elsevier B.V. All rights reserved.

Craniometric Similarities Within and Between Human Populations in Comparison with Neutral Genetic Data

The statement that pairs of individuals from different populations are often more genetically similar than pairs from the same population is a widespread idea inside and outside the scientific community. Witherspoon et al. [""Genetic similarities within and between human populations,"" Genetics 176:351-359 (2007)] proposed an index called the dissimilarity fraction (omega) to access in a quantitative way the validity of this statement for genetic systems. Witherspoon demonstrated that, as the number of loci increases, omega decreases to a point where, when enough sampling is available, the statement is false. In this study, we applied the dissimilarity fraction to Howells`s craniometric database to establish whether or not similar results are obtained for cranial morphological traits. Although in genetic studies thousands of loci are available, Howells`s database provides no more than 55 metric traits, making the contribution of each variable important. To cope with this limitation, we developed a routine that takes this effect into consideration when calculating. omega Contrary to what was observed for the genetic data, our results show that cranial morphology asymptotically approaches a mean omega of 0.3 and therefore supports the initial statement-that is, that individuals from the same geographic region do not form clear and discrete clusters-further questioning the idea of the existence of discrete biological clusters in the human species. Finally, by assuming that cranial morphology is under an additive polygenetic model, we can say that the population history signal of human craniometric traits presents the same resolution as a neutral genetic system dependent on no more than 20 loci.

Endothelial Nitric Oxide Genotypes and Haplotypes Are Not Associated with End-Stage Renal Disease

The identification of genetic markers associated with chronic kidney disease (CKD) may help to predict its development. Because reduced nitric oxide (NO) bioavailability and endothelial dysfunction are involved in CKD, genetic polymorphisms in the gene encoding the enzyme involved in NO synthesis (endothelial NO synthase [eNos]) may affect the susceptibility to CKD and the development of end-stage renal disease (ESRD). We compared genotype and haplotype distributions of three relevant eNOS polymorphisms (T(-786) C in the promoter region, Glu298Asp in exon 7, and 4b/4a in intron 4) in 110 healthy control subjects and 127 ESRD patients. Genotypes for the T(-786) C and Glu298Asp polymorphisms were determined by TaqMan (R) Allele Discrimination assay and real-time polymerase chain reaction. Genotypes for the intron 4 polymorphism were determined by polymerase chain reaction and fragment separation by electrophoresis. The software program PHASE 2.1 was used to estimate the haplotypes frequencies. We considered significant a probability value of p < 0.05/number of haplotypes (p < 0.05/8 = 0.0063). We found no significant differences between groups with respect to age, ethnicity, and gender. CKD patients had higher blood pressure, total cholesterol, and creatinine levels than healthy control subjects (all p < 0.05). Genotype and allele distributions for the three eNOS polymorphisms were similar in both groups (p > 0.05). We found no significant differences in haplotype distribution between groups (p > 0.05). The lack of significant associations between eNOS polymorphisms and ESRD suggests that eNOS polymorphisms may not be relevant to the genetic component of CKD that leads to ESRD.

The role of major histocompatibility complex Alleles in the susceptibility of Brazilian patients to develop the myogenic type of Graves' orbitopathy

Objective: To assess the frequency of the genetic markers HLA-DRB1 and DQB1 in patients with Graves' orbitopathy ( GO) with and without extraocular muscle involvement. Design: The frequencies of class II HLA-DRB1 and DQB1 allele groups were determined for 81 Brazilian patients with GO and 161 normal subjects. The patients were divided into myogenic and nonmyogenic groups based on the clinical characteristics of the orbitopathy and quantitative computed tomography analysis of the extraocular muscle ( EOM) dimensions. Main outcome: Compared to the frequency obtained for samples of normal subjects of the Brazilian population, HLA-DRB1*16 (p(c)= 0.008) was overrepresented in myogenic and HLA-DRB1*03 (p(c)= 0.02) in nonmyogenic patients. Conclusions: The association between the HLA-DRB1* 16 and the myogenic subtype of GO suggests that EOM involvement in GO may be genetically predisposed.

CONTRASTING PHYLOGEOGRAPHIC PATTERNS IN MITOCHONDRIAL DNA AND MICROSATELLITES: EVIDENCE OF FEMALE PHILOPATRY AND MALE-BIASED GENE FLOW AMONG REGIONAL POPULATIONS OF THE BLUE-AND-YELLOW MACAW (PSITTACIFORMES: ARA ARARAUNA) IN BRAZIL

Comparing the patterns of population differentiation among genetic markers with different modes of inheritance call provide insights into patterns of sex-biased dispersal and gene flow. The blue-and-yellow Macaw (Ara ararauna) is a Neotropical parrot with a broad geographic distribution ill South America. However, little is known about the natural history and current status Of remaining wild populations, including levels of genetic variability. The progressive decline and possible fragmentation of populations may endanger this species in the near future. We analyzed mitochondrial DNA (mtDNA) control-region sequences and six microsatellite 106 Of Blue-and-yellow Macaws sampled throughout their geographic range ill Brazil to describe population genetic Structure, to make inferences about historical demography and dispersal behavior, and to provide insight for conservation efforts. Analyses of population genetic structure based on mtDNA showed evidence of two major populations ill western and eastern Brazil that share a few low-frequency haplotypes. This phylogeographic pattern seems to have originated by the historical isolation of Blue-and-yellow Macaw populations similar to 374,000 years ago and has been maintained by restricted gene flow and female philopatry. By contrast, variation ill biparentally inherited microsatellites was not structured geographically, Male-biased dispersal and female philopatry best explain the different patterns observed in these two markers. Because females disperse less than males, the two regional populations with well-differentiated mtDNA haplogroups should be considered two different management units for conservation purposes. Received 4 November 2007 accepted 10 December 2008.

Induction of mitotic crossing-over in diploid strains of Aspergillus nidulans using low-dose X-rays

As a contribution towards detecting the genetic effects of low doses of genotoxic physical agents, this paper deals with the consequences of low-dose X-rays in the Aspergillus nidulans genome. The irradiation doses studied were those commonly used in dental clinics (1-5 cGy). Even very low doses promoted increased mitotic crossing-over frequencies in diploid strains heterozygous for several genetic markers including the ones involved in DNA repair and recombination mechanisms. Genetic markers of several heterozygous strains were individu`ally analyzed disclosing that some markers were especially sensitive to the treatments. These markers should be chosen as bio-indicators in the homozygotization index assay to better detect the recombinogenic/carcinogenic genomic effects of low-dose X-rays.

Characterization of bovine transcripts preferentially expressed in testis and with a putative role in spermatogenesis

Although the number of genes known to be associated with bovine spermatogenesis has increased in the past few years, regulation of this biological process remains poorly understood. Therefore, discovery of new male fertility genetic markers is of great value for assisted selection in commercially important cattle breeds, e.g., Nelore, that have delayed reproductive maturation and low fertility rates. The objective of the present study was to identify sequences associated with spermatogenesis that could be used as fertility markers. With RT-PCR, the following five transcripts preferentially expressed in adult testis were detected: TET(656) detected only in adult testis; TET(868) and TET(515) expressed preferentially in adult testis but also detected in fetal gonads of both sexes; and TET(456) and TET(262). expressed primarily in the testis, but also present in very low amounts in somatic tissues. Based on their homologies and expression profiles, we inferred that they had putative roles in spermatogenesis. Detection of sequences differentially expressed in testis, ovary, or both, was a useful approach for identifying new genes related to bovine spermatogenesis. The data reported here contributed to discovery of gene pathways involved in bovine spermatogenesis, with potential for prediction of fertility. (C) 2011 Elsevier Inc. All rights reserved.

Toxoplasma gondii isolates: Multi locus RFLP-PCR genotyping from human patients in Sao Paulo State, Brazil identified distinct genotypes

This study investigated the genetic characteristics of Toxoplasma gondii samples collected from 62 patients with toxoplasmosis in Sao Paulo State, Brazil. DNA samples were isolated from blood, cerebrospinal fluid and amniotic fluids of 25 patients with cerebral toxoplasmosis and AIDS, two patients with acute toxoplasmosis, 12 patients with ocular toxoplasmosis, six newborns with congenital toxoplasmosis and 17 pregnant women with acute infection. Diagnosis of toxoplasmosis was based in clinical, radiological and laboratory features. Genotyping was performed using multilocus PCR-RFLP genetic markers including SAG1, SAG2, 5`- and 3`-SAG2, alt.SAG2, SAG3, BTUB, GRA6, C22-8, c29-2, L358, PK1 and Apico. Among the 62 clinical samples, 20 (32%) were successfully genotyped at eight or more genetic loci and were grouped to three distinct genotypes. Eighteen samples belonged to ToxoDB Genotype #65 and the other two samples were identified as ToxoDB Genotypes #6 and #71, respectively (http://toxodb.org/toxo/). Patients presenting Genotypes #6 and #71 had severe and atypical cerebral toxoplasmosis, characterized by diffuse encephalitis without extensive brain lesions. These results indicate that T. gondii Genotype #65 may have a high frequency in causing human toxoplasmosis in Sao Paulo State, Brazil. This unusual finding highlights the need to investigate the possible association of parasite genotypes with human toxoplasmosis. (C) 2011 Elsevier Inc. All rights reserved.