Genetic diversity of Leishmania amazonensis strains isolated in northeastern Brazil as revealed by DNA sequencing, PCR-based analyses and molecular karyotyping

Abstract Background Leishmania (Leishmania) amazonensis infection in man results in a clinical spectrum of disease manifestations ranging from cutaneous to mucosal or visceral involvement. In the present study, we have investigated the genetic variability of 18 L. amazonensis strains isolated in northeastern Brazil from patients with different clinical manifestations of leishmaniasis. Parasite DNA was analyzed by sequencing of the ITS flanking the 5.8 S subunit of the ribosomal RNA genes, by RAPD and SSR-PCR and by PFGE followed by hybridization with gene-specific probes. Results ITS sequencing and PCR-based methods revealed genetic heterogeneity among the L. amazonensis isolates examined and molecular karyotyping also showed variation in the chromosome size of different isolates. Unrooted genetic trees separated strains into different groups. Conclusion These results indicate that L. amazonensis strains isolated from leishmaniasis patients from northeastern Brazil are genetically diverse, however, no correlation between genetic polymorphism and phenotype were found.

Permanent Genetic Resources added to Molecular Ecology Resources Database 1 December 2011-31 January 2012

This article documents the addition of 473 microsatellite marker loci and 71 pairs of single-nucleotide polymorphism (SNP) sequencing primers to the Molecular Ecology Resources Database. Loci were developed for the following species: Barteria fistulosa, Bombus morio, Galaxias platei, Hematodinium perezi, Macrocentrus cingulum Brischke (a.k.a. M.abdominalis Fab., M.grandii Goidanich or M.gifuensis Ashmead), Micropogonias furnieri, Nerita melanotragus, Nilaparvata lugens Stal, Sciaenops ocellatus, Scomber scombrus, Spodoptera frugiperda and Turdus lherminieri. These loci were cross-tested on the following species: Barteria dewevrei, Barteria nigritana, Barteria solida, Cynoscion acoupa, Cynoscion jamaicensis, Cynoscion leiarchus, Cynoscion nebulosus, Cynoscion striatus, Cynoscion virescens, Macrodon ancylodon, Menticirrhus americanus, Nilaparvata muiri and Umbrina canosai. This article also documents the addition of 116 sequencing primer pairs for Dicentrarchus labrax.

Genetic characterization of fig tree mutants with molecular markers

The fig (Ficus carica L.) is a fruit tree of great world importance and, therefore, the genetic improvement becomes an important field of research for better crops, being necessary to gather information on this species, mainly regarding its genetic variability so that appropriate propagation projects and management are made. The improvement programs of fig trees using conventional procedures in order to obtain new cultivars are rare in many countries, such as Brazil, especially due to the little genetic variability and to the difficulties in obtaining plants from gamete fusion once the wasp Blastophaga psenes, responsible for the natural pollinating, is not found in Brazil. In this way, the mutagenic genetic improvement becomes a solution of it. For this reason, in an experiment conducted earlier, fig plants formed by cuttings treated with gamma ray were selected based on their agronomic characteristics of interest. We determined the genetic variability in these fig tree selections, using RAPD and AFLP molecular markers, comparing them to each other and to the Roxo-de-Valinhos, used as the standard. For the reactions of DNA amplification, 140 RAPD primers and 12 primer combinations for AFLP analysis were used. The selections did not differ genetically between themselves and between them and the Roxo-de-Valinhos cultivar. Techniques that can detect polymorphism between treatments, such as DNA sequencing, must be tested. The phenotypic variation of plants may be due to epigenetic variation, necessitating the use of techniques with methylation-sensitive restriction enzymes.

Whole-genome sequencing of the efficient industrial fuel-ethanol fermentative Saccharomyces cerevisiae strain CAT-1

The Saccharomyces cerevisiae strains widely used for industrial fuel-ethanol production have been developed by selection, but their underlying beneficial genetic polymorphisms remain unknown. Here, we report the draft whole-genome sequence of the S. cerevisiae strain CAT-1, which is a dominant fuel-ethanol fermentative strain from the sugarcane industry in Brazil. Our results indicate that strain CAT-1 is a highly heterozygous diploid yeast strain, and the similar to 12-Mb genome of CAT-1, when compared with the reference S228c genome, contains similar to 36,000 homozygous and similar to 30,000 heterozygous single nucleotide polymorphisms, exhibiting an uneven distribution among chromosomes due to large genomic regions of loss of heterozygosity (LOH). In total, 58 % of the 6,652 predicted protein-coding genes of the CAT-1 genome constitute different alleles when compared with the genes present in the reference S288c genome. The CAT-1 genome contains a reduced number of transposable elements, as well as several gene deletions and duplications, especially at telomeric regions, some correlated with several of the physiological characteristics of this industrial fuel-ethanol strain. Phylogenetic analyses revealed that some genes were likely associated with traits important for bioethanol production. Identifying and characterizing the allelic variations controlling traits relevant to industrial fermentation should provide the basis for a forward genetics approach for developing better fermenting yeast strains.

Predictive models for mutations in mismatch repair genes: implication for genetic counseling in developing countries

Background: Lynch syndrome (LS) is the most common form of inherited predisposition to colorectal cancer (CRC), accounting for 2-5% of all CRC. LS is an autosomal dominant disease characterized by mutations in the mismatch repair genes mutL homolog 1 (MLH1), mutS homolog 2 (MSH2), postmeiotic segregation increased 1 (PMS1), post-meiotic segregation increased 2 (PMS2) and mutS homolog 6 (MSH6). Mutation risk prediction models can be incorporated into clinical practice, facilitating the decision-making process and identifying individuals for molecular investigation. This is extremely important in countries with limited economic resources. This study aims to evaluate sensitivity and specificity of five predictive models for germline mutations in repair genes in a sample of individuals with suspected Lynch syndrome. Methods: Blood samples from 88 patients were analyzed through sequencing MLH1, MSH2 and MSH6 genes. The probability of detecting a mutation was calculated using the PREMM, Barnetson, MMRpro, Wijnen and Myriad models. To evaluate the sensitivity and specificity of the models, receiver operating characteristic curves were constructed. Results: Of the 88 patients included in this analysis, 31 mutations were identified: 16 were found in the MSH2 gene, 15 in the MLH1 gene and no pathogenic mutations were identified in the MSH6 gene. It was observed that the AUC for the PREMM (0.846), Barnetson (0.850), MMRpro (0.821) and Wijnen (0.807) models did not present significant statistical difference. The Myriad model presented lower AUC (0.704) than the four other models evaluated. Considering thresholds of >= 5%, the models sensitivity varied between 1 (Myriad) and 0.87 (Wijnen) and specificity ranged from 0 (Myriad) to 0.38 (Barnetson). Conclusions: The Barnetson, PREMM, MMRpro and Wijnen models present similar AUC. The AUC of the Myriad model is statistically inferior to the four other models.

High-throughput sequencing of black pepper root transcriptome

Background: Black pepper (Piper nigrum L.) is one of the most popular spices in the world. It is used in cooking and the preservation of food and even has medicinal properties. Losses in production from disease are a major limitation in the culture of this crop. The major diseases are root rot and foot rot, which are results of root infection by Fusarium solani and Phytophtora capsici, respectively. Understanding the molecular interaction between the pathogens and the host's root region is important for obtaining resistant cultivars by biotechnological breeding. Genetic and molecular data for this species, though, are limited. In this paper, RNA-Seq technology has been employed, for the first time, to describe the root transcriptome of black pepper. Results: The root transcriptome of black pepper was sequenced by the NGS SOLiD platform and assembled using the multiple-k method. Blast2Go and orthoMCL methods were used to annotate 10338 unigenes. The 4472 predicted proteins showed about 52% homology with the Arabidopsis proteome. Two root proteomes identified 615 proteins, which seem to define the plant's root pattern. Simple-sequence repeats were identified that may be useful in studies of genetic diversity and may have applications in biotechnology and ecology. Conclusions: This dataset of 10338 unigenes is crucially important for the biotechnological breeding of black pepper and the ecogenomics of the Magnoliids, a major group of basal angiosperms.

Expanding the Clinical and Genetic Spectrum of Human CD40L Deficiency: The Occurrence of Paracoccidioidomycosis and Other Unusual Infections in Brazilian Patients

CD40 ligand (CD40L) deficiency or X-linked hyper-IgM syndrome (X-HIGM) is a well-described primary immunodeficiency in which Pneumocystis jiroveci pneumonia is a common clinical feature. We have identified an unusual high incidence of fungal infections and other not yet described infections in a cohort of 11 X-HIGM patients from nine unrelated Brazilian families. Among these, we describe the first case of paracoccidioidomycosis (PCM) in X-HIGM. The molecular genetic analysis of CD40L was performed by gene sequencing and evaluation of CD40L protein expression. Nine of these 11 patients (82%) had fungal infections. These included fungal species common to CD40L deficiency (P. jiroveci and Candida albicans) as well as Paracoccidioides brasiliensis. One patient presented with PCM at age 11 years and is now doing well at 18 years of age. Additionally, one patient presented with a simultaneous infection with Klebsiella and Acinetobacter, and one with condyloma caused by human papilloma virus. Molecular analysis revealed four previously described CD40L mutations, two novel missense mutations (c.433 T>G and c.476 G>C) resulting in the absence of CD40L protein expression by activated CD4(+) cells and one novel insertion (c.484_485insAA) within the TNFH domain leading to a frame shift and premature stop codon. These observations demonstrated that the susceptibility to fungal infections in X-HIGM extends beyond those typically associated with X-HIGM (P. jiroveci and C. albicans) and that these patients need to be monitored for those pathogens.

Permanent Genetic Resources added to Molecular Ecology Resources Database 1 February 2012-31 March 2012

This article documents the addition of 171 microsatellite marker loci and 27 pairs of single nucleotide polymorphism (SNP) sequencing primers to the Molecular Ecology Resources Database. Loci were developed for the following species: Bombus pauloensis, Cephalorhynchus heavisidii, Cercospora sojina, Harpyhaliaetus coronatus, Hordeum vulgare, Lachnolaimus maximus, Oceanodroma monteiroi, Puccinia striiformis f. sp. tritici, Rhea americana, Salmo salar, Salmo trutta, Schistocephalus solidus, Sousa plumbea and Tursiops aduncus. These loci were cross-tested on the following species: Aquila heliaca, Bulweria bulwerii, Buteo buteo, Buteo swainsoni, Falco rusticolus, Haliaeetus albicilla, Halobaena caerulea, Hieraaetus fasciatus, Oceanodroma castro, Puccinia graminis f. sp. Tritici, Puccinia triticina, Rhea pennata and Schistocephalus pungitii. This article also documents the addition of 27 sequencing primer pairs for Puffinus baroli and Bulweria bulwerii and cross-testing of these loci in Oceanodroma castro, Pelagodroma marina, Pelecanoides georgicus, Pelecanoides urinatrix, Thalassarche chrysostoma and Thalassarche melanophrys.

Evaluation of the genetic diversity of Histoplasma capsulatum var. capsulatum isolates from north-eastern Brazil

Since the beginning of the HIV epidemic, there has been a significant increase in the number of histoplasmosis cases in Ceara, a state in north-east Brazil. The lack of epidemiological data on the genotypes circulating in the north-east region shows the importance of more detailed studies on the molecular epidemiology of Histoplasma capsulatum var. capsulatum in this region. Different molecular techniques have been used to better characterize the genetic profile of H. capsulatum var. capsulatum strains. The aim of this study was to analyse the genetic diversity of H. capsulatum var. capsulatum isolates in Fortaleza, the capital of Ceara, through the sequencing of the internal transcribed spacer (ITS)1-5.8S-ITS2 region, and establish the molecular profile of these isolates, along with strains from south-east Brazil, by RAPD analysis, featuring the different clusters in those regions. The isolates were grouped into two clusters. Cluster 1 included strains from the south-east and north-east regions with separation of isolates into three distinct subgroups (subgroups 1a, 1 b and 1 c). Cluster 2 included only samples from north-east Brazil. Sequencing of the ITS1 -5.8S-ITS2 region allowed the detection of two major clades, which showed geographical correlation between them and their subgroups. Therefore, it can be concluded that the H. capsulatum var. capsulatum isolates from Ceara have a high degree of genetic polymorphism. The molecular data also confirm that populations of this fungus are composed of different genotypes in Brazil and worldwide.

Genetic polymorphisms modulate the folate metabolism of Brazilian individuals with Down syndrome

Individuals with Down syndrome (DS) carry three copies of the Cystathionine beta-synthase (C beta S) gene. The increase in the dosage of this gene results in an altered profile of metabolites involved in the folate pathway, including reduced homocysteine (Hcy), methionine, S-adenosylhomocysteine (SAH) and S-adenosylmethionine (SAM). Furthermore, previous studies in individuals with DS have shown that genetic variants in genes involved in the folate pathway influence the concentrations of this metabolism's products. The purpose of this study is to investigate whether polymorphisms in genes involved in folate metabolism affect the plasma concentrations of Hcy and methylmalonic acid (MMA) along with the concentration of serum folate in individuals with DS. Twelve genetic polymorphisms were investigated in 90 individuals with DS (median age 1.29 years, range 0.07-30.35 years; 49 male and 41 female). Genotyping for the polymorphisms was performed either by polymerase chain reaction (PCR) based techniques or by direct sequencing. Plasma concentrations of Hcy and MMA were measured by liquid chromatography-tandem mass spectrometry as previously described, and serum folate was quantified using a competitive immunoassay. Our results indicate that the MTHFR C677T, MTR A2756G, TC2 C776G and BHMT G742A polymorphisms along with MMA concentration are predictors of Hcy concentration. They also show that age and Hcy concentration are predictors of MMA concentration. These findings could help to understand how genetic variation impacts folate metabolism and what metabolic consequences these variants have in individuals with trisomy 21.

ALDH1A2 (RALDH2) genetic variation in human congenital heart disease

Abstract Background Signaling by the vitamin A-derived morphogen retinoic acid (RA) is required at multiple steps of cardiac development. Since conversion of retinaldehyde to RA by retinaldehyde dehydrogenase type II (ALDH1A2, a.k.a RALDH2) is critical for cardiac development, we screened patients with congenital heart disease (CHDs) for genetic variation at the ALDH1A2 locus. Methods One-hundred and thirty-three CHD patients were screened for genetic variation at the ALDH1A2 locus through bi-directional sequencing. In addition, six SNPs (rs2704188, rs1441815, rs3784259, rs1530293, rs1899430) at the same locus were studied using a TDT-based association approach in 101 CHD trios. Observed mutations were modeled through molecular mechanics (MM) simulations using the AMBER 9 package, Sander and Pmemd programs. Sequence conservation of observed mutations was evaluated through phylogenetic tree construction from ungapped alignments containing ALDH8 s, ALDH1Ls, ALDH1 s and ALDH2 s. Trees were generated by the Neighbor Joining method. Variations potentially affecting splicing mechanisms were cloned and functional assays were designed to test splicing alterations using the pSPL3 splicing assay. Results We describe in Tetralogy of Fallot (TOF) the mutations Ala151Ser and Ile157Thr that change non-polar to polar residues at exon 4. Exon 4 encodes part of the highly-conserved tetramerization domain, a structural motif required for ALDH oligomerization. Molecular mechanics simulation studies of the two mutations indicate that they hinder tetramerization. We determined that the SNP rs16939660, previously associated with spina bifida and observed in patients with TOF, does not affect splicing. Moreover, association studies performed with classical models and with the transmission disequilibrium test (TDT) design using single marker genotype, or haplotype information do not show differences between cases and controls. Conclusion In summary, our screen indicates that ALDH1A2 genetic variation is present in TOF patients, suggesting a possible causal role for this gene in rare cases of human CHD, but does not support the hypothesis that variation at the ALDH1A2 locus is a significant modifier of the risk for CHD in humans.

Genetic relatedness of Plasmodium falciparum isolates and the origin of allelic diversity at the merozoite surface protein-1 (MSP-1) locus in Brazil and Vietnam

Abstract Background Despite the extensive polymorphism at the merozoite surface protein-1 (MSP-1) locus of Plasmodium falciparum, that encodes a major repetitive malaria vaccine candidate antigen, identical and nearly identical alleles frequently occur in sympatric parasites. Here we used microsatellite haplotyping to estimate the genetic distance between isolates carrying identical and nearly identical MSP-1 alleles. Methods We analyzed 28 isolates from hypoendemic areas in north-western Brazil, collected between 1985 and 1998, and 23 isolates obtained in mesoendemic southern Vietnam in 1996. MSP-1 alleles were characterized by combining PCR typing with allele-specific primers and partial DNA sequencing. The following single-copy microsatellite markers were typed : Polyα, TA42 (only for Brazilian samples), TA81, TA1, TA87, TA109 (only for Brazilian samples), 2490, ARAII, PfG377, PfPK2, and TA60. Results The low pair-wise average genetic distance between microsatellite haplotypes of isolates sharing identical MSP-1 alleles indicates that epidemic propagation of discrete parasite clones originated most identical MSP-1 alleles in parasite populations from Brazil and Vietnam. At least one epidemic clone propagating in Brazil remained relatively unchanged over more than one decade. Moreover, we found no evidence that rearrangements of MSP-1 repeats, putatively created by mitotic recombination events, generated new alleles within clonal lineages of parasites in either country. Conclusion Identical MSP-1 alleles originated from co-ancestry in both populations, whereas nearly identical MSP-1 alleles have probably appeared independently in unrelated parasite lineages.

Genomic analysis of ERVWE2 locus in patients with multiple sclerosis: absence of genetic association but potential role of human endogenous retrovirus type W elements in molecular mimicry with myelin antigen

Human endogenous retroviruses (HERVs) arise from ancient infections of the host germline cells by exogenous retroviruses, constituting 8% of the human genome. Elevated level of envelope transcripts from HERVs-W has been detected in CSF, plasma and brain tissues from patients with Multiple Sclerosis (MS), most of them from Xq22.3, 15q21.3, and 6q21 chromosomes. However, since the locus Xq22.3 (ERVWE2) lack the 5' LTR promoter and the putative protein should be truncated due to a stop codon, we investigated the ERVWE2 genomic loci from 84 individuals, including MS patients with active HERV-W expression detected in PBMC. In addition, an automated search for promoter sequences in 20 kb nearby region of ERVWE2 reference sequence was performed. Several putative binding sites for cellular cofactors and enhancers were found, suggesting that transcription may occur via alternative promoters. However, ERVWE2 DNA sequencing of MS and healthy individuals revealed that all of them harbor a stop codon at site 39, undermining the expression of a full-length protein. Finally, since plaque formation in central nervous system (CNS) of MS patients is attributed to immunological mechanisms triggered by autoimmune attack against myelin, we also investigated the level of similarity between envelope protein and myelin oligodendrocyte glycoprotein (MOG). Comparison of the MOG to the envelope identified five retroviral regions similar to the Ig-like domain of MOG. Interestingly, one of them includes T and B cell epitopes, capable to induce T effector functions and circulating Abs in rats. In sum, although no DNA substitutions that would link ERVWE2 to the MS pathogeny was found, the similarity between the envelope protein to MOG extends the idea that ERVEW2 may be involved on the immunopathogenesis of MS, maybe facilitating the MOG recognizing by the immune system. Although awaiting experimental evidences, the data presented here may expand the scope of the endogenous retroviruses involvement on MS pathogenesis

High-throughput sequencing of small RNA transcriptomes reveals critical biological features targeted by microRNAs in cell models used for squamous cell cancer research

Abstract Background The implication of post-transcriptional regulation by microRNAs in molecular mechanisms underlying cancer disease is well documented. However, their interference at the cellular level is not fully explored. Functional in vitro studies are fundamental for the comprehension of their role; nevertheless results are highly dependable on the adopted cellular model. Next generation small RNA transcriptomic sequencing data of a tumor cell line and keratinocytes derived from primary culture was generated in order to characterize the microRNA content of these systems, thus helping in their understanding. Both constitute cell models for functional studies of microRNAs in head and neck squamous cell carcinoma (HNSCC), a smoking-related cancer. Known microRNAs were quantified and analyzed in the context of gene regulation. New microRNAs were investigated using similarity and structural search, ab initio classification, and prediction of the location of mature microRNAs within would-be precursor sequences. Results were compared with small RNA transcriptomic sequences from HNSCC samples in order to access the applicability of these cell models for cancer phenotype comprehension and for novel molecule discovery. Results Ten miRNAs represented over 70% of the mature molecules present in each of the cell types. The most expressed molecules were miR-21, miR-24 and miR-205, Accordingly; miR-21 and miR-205 have been previously shown to play a role in epithelial cell biology. Although miR-21 has been implicated in cancer development, and evaluated as a biomarker in HNSCC progression, no significant expression differences were seen between cell types. We demonstrate that differentially expressed mature miRNAs target cell differentiation and apoptosis related biological processes, indicating that they might represent, with acceptable accuracy, the genetic context from which they derive. Most miRNAs identified in the cancer cell line and in keratinocytes were present in tumor samples and cancer-free samples, respectively, with miR-21, miR-24 and miR-205 still among the most prevalent molecules at all instances. Thirteen miRNA-like structures, containing reads identified by the deep sequencing, were predicted from putative miRNA precursor sequences. Strong evidences suggest that one of them could be a new miRNA. This molecule was mostly expressed in the tumor cell line and HNSCC samples indicating a possible biological function in cancer. Conclusions Critical biological features of cells must be fully understood before they can be chosen as models for functional studies. Expression levels of miRNAs relate to cell type and tissue context. This study provides insights on miRNA content of two cell models used for cancer research. Pathways commonly deregulated in HNSCC might be targeted by most expressed and also by differentially expressed miRNAs. Results indicate that the use of cell models for cancer research demands careful assessment of underlying molecular characteristics for proper data interpretation. Additionally, one new miRNA-like molecule with a potential role in cancer was identified in the cell lines and clinical samples.

The humankind genome: from genetic diversity to the origin of human diseases

Genome-wide association studies have failed to establish common variant risk for the majority of common human diseases. The underlying reasons for this failure are explained by recent studies of resequencing and comparison of over 1200 human genomes and 10 000 exomes, together with the delineation of DNA methylation patterns (epigenome) and full characterization of coding and noncoding RNAs (transcriptome) being transcribed. These studies have provided the most comprehensive catalogues of functional elements and genetic variants that are now available for global integrative analysis and experimental validation in prospective cohort studies. With these datasets, researchers will have unparalleled opportunities for the alignment, mining, and testing of hypotheses for the roles of specific genetic variants, including copy number variations, single nucleotide polymorphisms, and indels as the cause of specific phenotypes and diseases. Through the use of next-generation sequencing technologies for genotyping and standardized ontological annotation to systematically analyze the effects of genomic variation on humans and model organism phenotypes, we will be able to find candidate genes and new clues for disease’s etiology and treatment. This article describes essential concepts in genetics and genomic technologies as well as the emerging computational framework to comprehensively search websites and platforms available for the analysis and interpretation of genomic data.