Current range characteristics of Swiss stone pine (Pinus cembra L.) along the Carpathians revealed by chloroplast SSR markers.

We investigated the diversity pattern of nine Swiss stone pine (Pinus cembra L.) populations along the Carpathian range including the High Tatras, by using six chloroplast DNA microsatellites (cpSSR). Our aim was to detect genetically distinct regions by clustering of populations, and to tackle possible historical colonization routes. Our analysis referred to an investigated geographical range with the two most distant populations situated at about 500 air km. We found that the most diverse populations are situated at the two edges of the investigated part, in the Retezat Mts. (South Carpathians) and the High Tatras, and diversity decreases towards the populations of the Eastern Carpathians. Hierarchical clustering and NMDS revealed that the populations of the South Carpathians with the Tatras form a distinct cluster, significantly separated from those of the Eastern Carpathians. Moreover, based on the most variable chloroplast microsatellites, the four populations of the two range edges are not significantly different. Our results, supported also by palynological and late glacial macrofossil evidences, indicate refugial territories within the Retezat Mts. that conserved rich haplotype composition. From this refugial territory Pinus cembra might have colonized the Eastern Carpathians, and this was accompanied by a gradual decrease in population diversity. Populations of the High Tatras might have had the same role in the colonizing events of the Carpathians, as positive correlation was detected among populations lying from each other at a distance of 280 km, the maximum distance between neighbouring populations.

Molecular analysis of evolutionary relationships of Phaseolus dumosus with the gene pool of common bean

Valuable genetic variation for bean breeding programs is held within the common bean secondary gene pool which consists of Phaseolus albescens, P. coccineus, P. costaricensis, and P. dumosus. However, the use of close relatives for bean improvement is limited due to the lack of knowledge about genetic variation and genetic plasticity of many of these species. Characterisation and analysis of the genetic diversity is necessary among beans' wild relatives; in addition, conflicting phylogenies and relationships need to be understood and a hypothesis of a hybrid origin of P. dumosus needs to be tested. This thesis research was orientated to generate information about the patterns of relationships among the common bean secondary gene pool, with particular focus on the species Phaseolus dumosus. This species displays a set of characteristics of agronomic interest, not only for the direct improvement of common bean but also as a source of valuable genes for adaptation to climate change. Here I undertake the first comprehensive study of the genetic diversity of P. dumosus as ascertained from both nuclear and chloroplast genome markers. A germplasm collection of the ancestral forms of P. dumosus together with wild, landrace and cultivar representatives of all other species of the common bean secondary gene pool, were used to analyse genetic diversity, phylogenetic relationships and structure of P. dumosus. Data on molecular variation was generated from sequences of cpDNA loci accD-psaI spacer, trnT-trnL spacer, trnL intron and rps14-psaB spacer and from the nrDNA the ITS region. A whole genome DArT array was developed and used for the genotyping of P. dumosus and its closes relatives. 4208 polymorphic markers were generated in the DArT array and from those, 742 markers presented a call rate >95% and zero discordance. DArT markers revealed a moderate genetic polymorphism among P. dumosus samples (13% of polymorphic loci), while P. coccineus presented the highest level of polymorphism (88% of polymorphic loci). At the cpDNA one ancestral haplotype was detected among all samples of all species in the secondary genepool. The ITS region of P. dumosus revealed high homogeneity and polymorphism bias to P. coccineus genome. Phylogenetic reconstructions made with Maximum likelihood and Bayesian methods confirmed previously reported discrepancies among the nuclear and chloroplast genomes of P. dumosus. The outline of relationships by hybridization networks displayed a considerable number of interactions within and between species. This research provides compelling evidence that P. dumosus arose from hybridisation between P. vulgaris and P. coccineus and confirms that P. costaricensis has likely been involved in the genesis or backcrossing events (or both) in the history of P. dumosus. The classification of the specie P. persistentus was analysed based on cpDNA and ITS sequences, the results found this species to be highly related to P. vulgaris but not too similar to P. leptostachyus as previously proposed. This research demonstrates that wild types of the secondary genepool carry a significant genetic variation which makes this a valuable genetic resource for common bean improvement. The DArT array generated in this research is a valuable resource for breeding programs since it has the potential to be used in several approaches including genotyping, discovery of novel traits, mapping and marker-trait associations. Efforts should be made to search for potential populations of P. persistentus and to increase the collection of new populations of P. dumosus, P. albescens and P. costaricensis that may provide valuable traits for introgression into common bean and other Phaseolus crops.

The onset of lactase persistence in Europe

The genomic region containing the lactase (LCT) gene shows one of the strongest signals of positive selection in Europeans, detectable using a range of approaches including haplotype length, linked microsatellite variation and population-differentiation-based tests. Lactase is the enzyme that carries out the digestion of the milk sugar lactose. Its expression decreases at some point after the weaning period is over in most mammals and in around 68% of all living adult humans. However, in some humans, particularly those from populations with a history of dairying, lactase is expressed throughout adulthood. This trait is called lactase persistence (LP), and in people of European ancestry, it is associated with a single mutation (-13910*T). Evidence from the detection of dairy fat residues in potsherds, and allele frequencies in ancient DNA samples suggest that LP arose after dairying practices had developed. However, the reasons why LP may have been advantageous are still debated, and the respective contribution of demography and natural selection remains to be disentangled. This paper discusses various studies, from archaeology to population genetics, that have shed some light on the subject by investigating the evolution of LP in Europe.

The relationship between vegetation and modern pollen assemblages on Mount Paggeo (NE Greece)

In this study, we investigated the relationship between vegetation and modern-pollen rain along the elevational gradient of Mount Paggeo. We apply multivariate data analysis to assess the relationship between vegetation and modern-pollen rain and quantify the representativeness of forest zones. This study represents the first statistical analysis of pollen-vegetation relationship along an elevational gradient in Greece. Hence, this paper improves confidence in interpretation of palynological records from north-eastern Greece and may refine past climate reconstructions for a more accurate comparison of data and modelling. Numerical classification and ordination were performed on pollen data to assess differences among plant communities that beech (Fagus sylvatica) dominates or co-dominates. The results show a strong relationship between altitude, arboreal cover, human impact and variations in pollen and nonpollen palynomorph taxa percentages.

DNA from historical and trophy samples provides insights into white shark population origins and genetic diversity

Characterizing genetic variation by retrospective genotyping of trophy or historical artifacts from endangered species is an important conservation tool. Loss of genetic diversity in top predators such as the white shark Carcharodon carcharias remains an issue, exacerbated in this species by declining, sometimes isolated philopatric populations. We successfully sequenced mitochondrial DNA (mtDNA) D-loop from osteodentine of contemporary South African white shark teeth (from 3 jaws), and from 34 to 129 yr old dried cartilage and skin samples from 1 Pacific Ocean and 5 Mediterranean sharks. Osteodentine-derived sequences from South African fish matched those derived from an individual’s finclips, but were generally of poorer quality than those from skin and cartilage of historical samples. Three haplotypes were identified from historical Mediterranean samples (n = 5); 2 individuals had unique sequences and 3 shared the contemporary Mediterranean haplotype. Placement of previously undescribed mtDNA haplotypes from historical material within both the Mediterranean and Pacific clades fits with the accepted intra-specific phylogeny derived from contemporary material, verifying our approaches. The utility of our methodology is in its provision of additional genetic resources from osteodentine (for species lacking tooth pulp) and cartilage of rare and endangered species held in often uncurated, contemporary and historical dry collections. Such material can usefully supplement estimates of connectivity, population history, and stock viability. We confirm the depauperate haplotype diversity of historical Mediterranean sharks, consistent with founding by a small number of Pacific colonizers. The consequent lack of diversity suggests serious challenges for the maintenance of this top predator and the Mediterranean ecosystem.

DNA from historical and trophy samples provides insights into white shark population origins and genetic diversity

Characterizing genetic variation by retrospective genotyping of trophy or historical artifacts from endangered species is an important conservation tool. Loss of genetic diversity in top predators such as the white shark Carcharodon carcharias remains an issue, exacerbated in this species by declining, sometimes isolated philopatric populations. We successfully sequenced mitochondrial DNA (mtDNA) D-loop from osteodentine of contemporary South African white shark teeth (from 3 jaws), and from 34 to 129 yr old dried cartilage and skin samples from 1 Pacific Ocean and 5 Mediterranean sharks. Osteodentine-derived sequences from South African fish matched those derived from an individual’s finclips, but were generally of poorer quality than those from skin and cartilage of historical samples. Three haplotypes were identified from historical Mediterranean samples (n = 5); 2 individuals had unique sequences and 3 shared the contemporary Mediterranean haplotype. Placement of previously undescribed mtDNA haplotypes from historical material within both the Mediterranean and Pacific clades fits with the accepted intra-specific phylogeny derived from contemporary material, verifying our approaches. The utility of our methodology is in its provision of additional genetic resources from osteodentine (for species lacking tooth pulp) and cartilage of rare and endangered species held in often uncurated, contemporary and historical dry collections. Such material can usefully supplement estimates of connectivity, population history, and stock viability. We confirm the depauperate haplotype diversity of historical Mediterranean sharks, consistent with founding by a small number of Pacific colonizers. The consequent lack of diversity suggests serious challenges for the maintenance of this top predator and the Mediterranean ecosystem.

Mathematical Modeling to Assess the Drivers of the Recent Emergence of Typhoid Fever in Blantyre, Malawi

BACKGROUND: Multiyear epidemics of Salmonella enterica serovar Typhi have been reported from countries across eastern and southern Africa in recent years. In Blantyre, Malawi, a dramatic increase in typhoid fever cases has recently occurred, and may be linked to the emergence of the H58 haplotype. Strains belonging to the H58 haplotype often exhibit multidrug resistance and may have a fitness advantage relative to other Salmonella Typhi strains.

METHODS: To explore hypotheses for the increased number of typhoid fever cases in Blantyre, we fit a mathematical model to culture-confirmed cases of Salmonella enterica infections at Queen Elizabeth Central Hospital, Blantyre. We explored 4 hypotheses: (1) an increase in the basic reproductive number (R0) in response to increasing population density; (2) a decrease in the incidence of cross-immunizing infection with Salmonella Enteritidis; (3) an increase in the duration of infectiousness due to failure to respond to first-line antibiotics; and (4) an increase in the transmission rate following the emergence of the H58 haplotype.

RESULTS: Increasing population density or decreasing cross-immunity could not fully explain the observed pattern of typhoid emergence in Blantyre, whereas models allowing for an increase in the duration of infectiousness and/or the transmission rate of typhoid following the emergence of the H58 haplotype provided a good fit to the data.

CONCLUSIONS: Our results suggest that an increase in the transmissibility of typhoid due to the emergence of drug resistance associated with the H58 haplotype may help to explain recent outbreaks of typhoid in Malawi and similar settings in Africa.

Genomic and archaeological evidence suggest a dual origin of domestic dogs

The geographic and temporal origins of dogs remain controversial. We generated genetic sequences from 59 ancient dogs and a complete (28x) genome of a late Neolithic dog (dated to ~4800 calendar years before the present) from Ireland. Our analyses revealed a deep split separating modern East Asian and Western Eurasian dogs. Surprisingly, the date of this divergence (~14,000 to 6400 years ago) occurs commensurate with, or several millennia after, the first appearance of dogs in Europe and East Asia. Additional analyses of ancient and modern mitochondrial DNA revealed a sharp discontinuity in haplotype frequencies in Europe. Combined, these results suggest that dogs may have been domesticated independently in Eastern and Western Eurasia from distinct wolf populations. East Eurasian dogs were then possibly transported to Europe with people, where they partially replaced European Paleolithic dogs.

Population structure and conservation implications for the loggerhead sea turtle of the Cape Verde Islands

[EN] The Cape Verde Islands harbour the second largest nesting aggregation of the globally endangered loggerhead sea turtle in the Atlantic. To characterize the unknown genetic structure, connectivity, and demographic history of this population, we sequenced a segment of the mitochondrial (mt) DNA control region (380 bp, n = 186) and genotyped 12 microsatellite loci (n = 128) in females nesting at three islands of Cape Verde. No genetic differentiation in either haplotype or allele frequencies was found among the islands (mtDNA FST = 0. 001, P > 0. 02; nDNA FST = 0. 001, P > 0. 126). However, population pairwise comparisons of the mtDNA data revealed significant differences between Cape Verde and all previously sequenced Atlantic and Mediterranean rookeries (FST = 0. 745; P < 0. 000).

Origin and dispersal routes of foreign green and Kemp's ridley turtles in Spanish Atlantic and Mediterranean waters

[EN] The presence of the green and Kemp's ridley turtles is rare at Atlantic and Mediterranean Spanish waters, but the records have increased during the last decades. We reported a new set of records and reviewed all the historical observations of these species. The analysis of a mitochondrial DNA fragment of the newest records provided insights about the origin of the individuals. The Kemp's ridley turtles arrived from the western Atlantic nesting beaches, although the discovering of a new haplotype suggested the existence of an unknown or low sampled nesting area of origin.

Exploration génétique de l'hypothyroïdie congénitale par dysgénésie thyroïdienne

L'hypothyroïdie congénitale par dysgénésie thyroïdienne (HCDT, ectopie dans plus de 80 %) a une prévalence de 1 cas sur 4000 naissances vivantes. L’HCDT est la conséquence d'une défaillance de la thyroïde embryonnaire à se différencier, à se maintenir ou à migrer vers sa localisation anatomique (partie antérieure du cou), qui aboutit à une absence totale de la thyroïde (athyréose) ou à une ectopie thyroïdienne (linguale ou sublinguale). Les HCDT sont principalement non-syndromiques (soit 98% des cas sont non-familiale), ont un taux de discordance de 92% chez les jumeaux monozygotes, et ont une prédominance féminine et ethnique (i.e., Caucasienne). La majorité des cas d’HCDT n’a pas de cause connue, mais est associée à un déficit sévère en hormones thyroïdiennes (hypothyroïdie). Des mutations germinales dans les facteurs de transcription liés à la thyroïde (NKX2.1, FOXE1, PAX8, NKX2.5) ont été identifiées dans seulement 3% des patients atteints d’HCDT sporadiques et l’analyse de liaisons exclue ces gènes dans les rares familles multiplex avec HCDT. Nous supposons que le manque de transmission familiale claire d’HCDT peut résulter de la nécessité d’au moins deux « hits » génétiques différents dans des gènes importants pour le développement thyroïdien. Pour répondre au mieux nos questions de recherche, nous avons utilisé deux approches différentes: 1) une approche gène candidat, FOXE1, seul gène impliqué dans l’ectopie dans le modèle murin et 2) une approche en utilisant les techniques de séquençage de nouvelle génération (NGS) afin de trouver des variants génétiques pouvant expliquer cette pathologie au sein d’une cohorte de patients avec HCDT. Pour la première approche, une étude cas-contrôles a été réalisée sur le promoteur de FOXE1. Il a récemment été découvert qu’une région du promoteur de FOXE1 est différentiellement méthylée au niveau de deux dinucléotides CpG consécutifs, définissant une zone cruciale de contrôle de l’expression de FOXE1. L’analyse d’association basée sur les haplotypes a révélé qu’un haplotype (Hap1: ACCCCCCdel1C) est associé avec le HCDT chez les Caucasiens (p = 5x10-03). Une réduction significative de l’activité luciférase est observée pour Hap1 (réduction de 68%, p<0.001) comparé au promoteur WT de FOXE1. Une réduction de 50% de l’expression de FOXE1 dans une lignée de cellules thyroïdienne humaine est suffisante pour réduire significativement la migration cellulaire (réduction de 55%, p<0.05). Un autre haplotype (Hap2: ACCCCCCC) est observé moins fréquemment chez les Afro-Américain comparés aux Caucasiens (p = 1.7x10-03) et Hap2 diminue l’activité luciférase (réduction de 26%, p<0.001). Deux haplotypes distincts sont trouvés fréquemment dans les contrôles Africains (Black-African descents). Le premier haplotype (Hap3: GTCCCAAC) est fréquent (30.2%) chez les contrôles Afro-Américains comparés aux contrôles Caucasiens (6.3%; p = 2.59 x 10-9) tandis que le second haplotype (Hap4: GTCCGCAC) est trouvé exclusivement chez les contrôles Afro-Américains (9.4%) et est absent chez les contrôles Caucasiens (P = 2.59 x 10-6). Pour la deuxième approche, le séquençage de l’exome de l’ADN leucocytaire entre les jumeaux MZ discordants n’a révélé aucune différence. D'où l'intérêt du projet de séquençage de l’ADN et l’ARN de thyroïdes ectopiques et orthotopiques dans lesquelles de l'expression monoallélique aléatoire dans a été observée, ce qui pourrait expliquer comment une mutation monoallélique peut avoir des conséquences pathogéniques. Finalement, le séquençage de l’exome d’une cohorte de 36 cas atteints d’HCDT a permis d’identifier de nouveaux variants probablement pathogéniques dans les gènes récurrents RYR3, SSPO, IKBKE et TNXB. Ces quatre gènes sont impliqués dans l’adhésion focale (jouant un rôle dans la migration cellulaire), suggérant un rôle direct dans les défauts de migration de la thyroïde. Les essais de migration montrent une forte diminution (au moins 60% à 5h) de la migration des cellules thyroïdiennes infectées par shRNA comparés au shCtrl dans 2 de ces gènes. Des zebrafish KO (-/- et +/-) pour ces nouveaux gènes seront réalisés afin d’évaluer leur impact sur l’embryologie de la thyroïde.

Novel pharmacogenomic markers of irinotecan-induced severe toxicity in metastatic colorectal cancer patients

L’irinotécan est un agent de chimiothérapie largement utilisé pour le traitement de tumeurs solides, particulièrement pour le cancer colorectal métastatique (mCRC). Fréquemment, le traitement par l’irinotécan conduit à la neutropénie et la diarrhée, des effets secondaires sévères qui peuvent limiter la poursuite du traitement et la qualité de vie des patients. Plusieurs études pharmacogénomiques ont évalué les risques associés à la chimiothérapie à base d’irinotécan, en particulier en lien avec le gène UGT1A, alors que peu d’études ont examiné l’impact des gènes codant pour des transporteurs. Par exemple, le marqueur UGT1A1*28 a été associé à une augmentation de 2 fois du risque de neutropénie, mais ce marqueur ne permet pas de prédire la toxicité gastrointestinale ou l’issue clinique. L’objectif de cette étude était de découvrir de nouveaux marqueurs génétiques associés au risque de toxicité induite par l’irinotécan, en utilisant une stratégie d’haplotype/SNP-étiquette permettant de maximiser la couverture des loci génétiques ciblés. Nous avons analysé les associations génétiques des loci UGT1 et sept gènes codants pour des transporteurs ABC impliqués dans la pharmacocinétique de l’irinotécan, soient ABCB1, ABCC1, ABCC2, ABCC5, ABCG1, ABCG2 ainsi que SLCO1B1. Les profils de 167 patients canadiens atteints de mCRC sous traitement FOLFIRI (à base d’irinotécan) ont été examinés et les marqueurs significatifs ont par la suite été validés dans une cohorte indépendante de 250 patients italiens. Nous avons découvert dans la région intergénique en aval du gène UGT1, un nouveau marqueur (rs11563250G) associé à un moindre risque de neutropénie sévère (rapport des cotes (RC)=0.21; p=0.043 chez les canadiens, RC=0.27; p=0.036 chez les italiens, et RC=0.31 p=0.001 pour les deux cohortes combinées). De plus, le RC est demeuré significatif après correction pour multiples comparaisons (p=0.041). Par ailleurs, pour l’haplotype défini par les marqueurs rs11563250G et UGT1A1*1 (rs8175347 TA6), le RC était de 0.17 (p=0.0004). Un test génétique évaluant ces marqueurs permettrait d’identifier les patients susceptibles de bénéficier d’une augmentation de dose d’irinotécan. En revanche, une autre combinaison de marqueurs, ABCC5 rs3749438 et rs10937158 (T–C), a prédit un risque plus faible de diarrhée sévère dans les deux cohortes (RC = 0.43; p=0.001). La coexistence des marqueurs ABCG1 rs225440T et ABCC5 rs2292997A a prédit un risque accru de neutropénie (RC=5.93; p=0.0002), alors qu’une prédiction encore plus significative a été obtenue lorsque ces marqueurs sont combinés au marqueur de risque bien établi UGT1A1*28 rs8175347 (RC=7.68; p<0.0001). Enfin, les porteurs de l’allèle de protection UGT1 rs11563250G en absence d’allèles de risque, ont montré une incidence réduite de neutropénie sévère (8.2% vs. 34.0%; p<0.0001). Nous concluons que ces nouveaux marqueurs génétiques prédictifs pourraient permettre d’améliorer l’évaluation du risque de toxicité et personnaliser le traitement à base d’irinotécan pour les patients atteints du cancer colorectal métastatique.

Non-Celiac Gluten Sensitivity and Autoimmunity. A Case Report

Introduction, objective: To present a case report in which the finding of non-coeliac gluten sensitivity was decisive for the treatment of a complex autoimmune disease. Materials and methods: A 43-year-old woman with polyarthritis, psoriatic features, anti-SSA/Ro and anti-cyclic citrullinated peptide antibodies, with refractory course, was evaluated for gluten sensitivity despite negative serology for coeliac disease. Results: The patient carried the HLA DQ2 haplotype and duodenal biopsy showed lymphocytic enteritis. A gluten-free diet resolved the clinical picture and permitted tapering of immunosuppressive therapy. Conclusion: Non-coeliac gluten sensitivity can be associated with autoimmunity despite the absence of the specific autoantibodies of coeliac disease.

Soybean Stem Fly, Melanagromyza sojae (Diptera: Agromyzidae), in the New World: detection of high genetic diversity from soybean fields in Brazil

Soybean Stem Fly (SSF), Melanagromyza sojae (Zehntner), belongs to the family Agromyzidae and is highly polyphagous, attacking many plant species of the family Fabaceae, including soybean and other beans. SSF is regarded as one of the most important pests in soybean fields of Asia (e.g., China, India), North East Africa (e.g., Egypt), parts of Russia, and South East Asia. Despite reports of Agromyzidae flies infesting soybean fields in Rio Grande do Sul State (Brazil) in 1983 and 2009 and periodic interceptions of SSF since the 1940s by the USA quarantine authorities, SSF has not been officially reported to have successfully established in the North and South Americas. In South America, M. sojae was recently confirmed using morphology and its complete mitochondrial DNA (mtDNA) was characterized. In the present study, we surveyed the genetic diversity of M. sojae, collected directly from soybean host plants, using partial mtDNA cytochrome oxidase I (COI) gene, and provide evidence of multiple (>10) maternal lineages in SSF populations in South America, potentially representing multiple incursion events. However, a single incursion involving multiple-female founders could not be ruled out. We identified a haplotype that was common in the fields of two Brazilian states and the individuals collected from Australia in 2013. The implications of SSF incursions in southern Brazil are discussed in relation to the current soybean agricultural practices, highlighting an urgent need for better understanding of SSF population movements in the New World, which is necessary for developing effective management options for this significant soybean pest. © FUNPEC-RP.

Determinantes genéticos da hanseníase em uma população do Rio Grande do Norte

Background: Leprosy can cause severe disability and disfigurement and is still a major health in different parts of the world. Only a subset of those individuals exposed to the pathogen will go on to develop clinical disease and there is a broad clinical spectrum amongst leprosy patients. The outcome of infection is in part due to host genes that influence control of the initial infection and the host´s immune response to that infection. Aim: Evaluate if polymorphisms type SNP in the 17q118q21 chromosomic region contribute to development of leprosy in Rio Grande do Norte population. Material and methods: A sample composed of 215 leprosy patients and 229 controls drawn from the same population were genotyped by using a Snapshot assay for eight genes (NOS2A, CCL18, CRLF3, CCL23, TNFAIP1, STAT5B, CCR7 and CSF3) located in chromosomic region 17q118q21. The genotype and allele frequency were measured and statistical analysis was performed by chi-square in SPSS version 15 and graph prism pad version 4 software. Results: Ours results indicated that the markers NOS2A8277, NOS2A8rs16949, CCR78rs11574663 and CSF38rs2227322 presented strong association with leprosy and their risk genotype were GG, TT, AA and GG respectively. The risk genotypes for all markers associated to leprosy presented recessive inheritance standard. When we compared the interaction among the markers in different combination we find that the marker NOS2A8277 associated with CCR78rs11574663 presented highest risk probability to development of leprosy. When we evaluated the haplotype of the risk markers it was found a haplotype associated with increase of the protection (CSF38rs22273228CC, CCR78 rs115746638GA, NOS2A8rs169498CT and NOS2A82778GA). The association of the clinical forms paucibacilary and multibacilary with markers showed that to the markers NOS2A8 2778GG, CCR78rs115746638AA and CSF38rs22273228GG there were a strong influence to migration to multibacilary pole and to marker NOS2A8rs169498TT the high proportion was found to the paucibacilary form. Conclusions: Changes in the genes NOS2A, CCR7 and CSF3 can influence the immune response against Mycobacterium leprae. The combination among these polymorphisms alters the risk probability to develop leprosy. The markers type SNP associated to development of the leprosy also are linked to clinical forms and its severity being the polymorphism NOS2A8rs169498TT associated with paucibacilar form and the polymorphisms NOS2A82778GG, CCR78rs115746638AA and CSF38rs22273228GG associated to multibacilar form