Genome-wide DNA polymorphism analyses using VariScan

BACKGROUND: DNA sequence polymorphisms analysis can provide valuable information on the evolutionary forces shaping nucleotide variation, and provides an insight into the functional significance of genomic regions. The recent ongoing genome projects will radically improve our capabilities to detect specific genomic regions shaped by natural selection. Current available methods and software, however, are unsatisfactory for such genome-wide analysis. RESULTS: We have developed methods for the analysis of DNA sequence polymorphisms at the genome-wide scale. These methods, which have been tested on a coalescent-simulated and actual data files from mouse and human, have been implemented in the VariScan software package version 2.0. Additionally, we have also incorporated a graphical-user interface. The main features of this software are: i) exhaustive population-genetic analyses including those based on the coalescent theory; ii) analysis adapted to the shallow data generated by the high-throughput genome projects; iii) use of genome annotations to conduct a comprehensive analyses separately for different functional regions; iv) identification of relevant genomic regions by the sliding-window and wavelet-multiresolution approaches; v) visualization of the results integrated with current genome annotations in commonly available genome browsers. CONCLUSION: VariScan is a powerful and flexible suite of software for the analysis of DNA polymorphisms. The current version implements new algorithms, methods, and capabilities, providing an important tool for an exhaustive exploratory analysis of genome-wide DNA polymorphism data.

The expansion of amino-acid repeats is not associated to adaptive evolution in mammalian genes.

BACKGROUND: The expansion of amino acid repeats is determined by a high mutation rate and can be increased or limited by selection. It has been suggested that recent expansions could be associated with the potential of adaptation to new environments. In this work, we quantify the strength of this association, as well as the contribution of potential confounding factors. RESULTS: Mammalian positively selected genes have accumulated more recent amino acid repeats than other mammalian genes. However, we found little support for an accelerated evolutionary rate as the main driver for the expansion of amino acid repeats. The most significant predictors of amino acid repeats are gene function and GC content. There is no correlation with expression level. CONCLUSIONS: Our analyses show that amino acid repeat expansions are causally independent from protein adaptive evolution in mammalian genomes. Relaxed purifying selection or positive selection do not associate with more or more recent amino acid repeats. Their occurrence is slightly favoured by the sequence context but mainly determined by the molecular function of the gene.

Two high throughput technologies to detect segmental aneuploidies identify new Williams-Beuren syndrome patients with atypical deletions.

OBJECTIVE: To develop and compare two new technologies for diagnosing a contiguous gene syndrome, the Williams-Beuren syndrome (WBS). METHODS: The first proposed method, named paralogous sequence quantification (PSQ), is based on the use of paralogous sequences located on different chromosomes and quantification of specific mismatches present at these loci using pyrosequencing technology. The second exploits quantitative real time polymerase chain reaction (QPCR) to assess the relative quantity of an analysed locus. RESULTS: A correct and unambiguous diagnosis was obtained for 100% of the analysed samples with either technique (n = 165 and n = 155, respectively). These methods allowed the identification of two patients with atypical deletions in a cohort of 182 WBS patients. Both patients presented with mild facial anomalies, mild mental retardation with impaired visuospatial cognition, supravalvar aortic stenosis, and normal growth indices. These observations are consistent with the involvement of GTF2IRD1 or GTF2I in some of the WBS facial features. CONCLUSIONS: Both PSQ and QPCR are robust, easy to interpret, and simple to set up. They represent a competitive alternative for the diagnosis of segmental aneuploidies in clinical laboratories. They have advantages over fluorescence in situ hybridisation or microsatellites/SNP genotyping for detecting short segmental aneuploidies as the former is costly and labour intensive while the latter depends on the informativeness of the polymorphisms.

Plastid genome characterisation in Brassica and Brassicaceae using a new set of nine SSRs.

We report a new set of nine primer pairs specifically developed for amplification of Brassica plastid SSR markers. The wide utility of these markers is demonstrated for haplotype identification and detection of polymorphism in B. napus, B. nigra, B. oleracea, B. rapa and in related genera Arabidopsis, Camelina, Raphanus and Sinapis. Eleven gene regions (ndhB-rps7 spacer, rbcL-accD spacer, rpl16 intron, rps16 intron, atpB-rbcL spacer, trnE-trnT spacer, trnL intron, trnL-trnF spacer, trnM-atpE spacer, trnR-rpoC2 spacer, ycf3-psaA spacer) were sequenced from a range of Brassica and related genera for SSR detection and primer design. Other sequences were obtained from GenBank/EMBL. Eight out of nine selected SSR loci showed polymorphism when amplified using the new primers and a combined analysis detected variation within and between Brassica species, with the number of alleles detected per locus ranging from 5 (loci MF-6, MF-1) to 11 (locus MF-7). The combined SSR data were used in a neighbour-joining analysis (SMM, D (DM) distances) to group the samples based on the presence and absence of alleles. The analysis was generally able to separate plastid types into taxon-specific groups. Multi-allelic haplotypes were plotted onto the neighbour joining tree. A total number of 28 haplotypes were detected and these differentiated 22 of the 41 accessions screened from all other accessions. None of these haplotypes was shared by more than one species and some were not characteristic of their predicted type. We interpret our results with respect to taxon differentiation, hybridisation and introgression patterns relating to the 'Triangle of U'.

Three-stranded DNA structure; is this the secret of DNA homologous recognition?

A novel type of triple-stranded DNA structure was proposed by several groups to play a crucial role in homologous recognition between single- and double-stranded DNA molecules. In this still putative structure a duplex DNA was proposed to co-ordinate a homologous single strand in its major groove side. In contrast to the well-characterized pyrimidine-purine-pyrimidine triplexes in which the two like strands are antiparallel and which are restricted to poly-pyrimidine-containing stretches, the homology-specific triplexes would have like strands in parallel orientation and would not be restricted to any particular sequence provided that there is a homology between interacting DNA molecules. For many years the stereo-chemical possibility of forming homology-dependent three- or four-stranded DNA structures during the pairing stage of recombination reactions was seriously considered in published papers. However, only recently has there been a marked increase in the number of papers that have directly tested the formation of triple-stranded DNA structures during the actual pairing stage of the recombination reaction. Unfortunately the results of these tests are not totally clear cut; while some laboratories presented experimental evidence consistent with the formation of triplexes, others studying the same or very similar systems offered alternative explanations. The aim of this review is to present the current state of the central question in the mechanism of homologous recombination, namely, what kind of DNA structure is responsible for DNA homologous recognition. Is it a novel triplex structure or just a classical duplex?

Study of chromosome rearrangements associated with the trpE26 mutation of Bacillus subtilis.

Chromosome rearrangements involved in the formation of merodiploid strains in the Bacillus subtilis 168-166 system were explained by postulating the existence of intrachromosomal homology regions. This working hypothesis was tested by analysing sequences and restriction patterns of the, as yet uncharacterized, junctions between chromosome segments undergoing rearrangements in parent, 168 trpC2 and 166 trpE26, as well as in derived merodiploid strains. Identification, at the Ia/Ib chromosome junction of both parent strains, of a 1.3 kb segment nearly identical to a segment of prophage SPbeta established the existence of one of the postulated homology sequences. Inspection of relevant junctions revealed that a set of different homology regions, derived from prophage SPbeta, plays a key role in the formation of so-called trpE30, trpE30+, as well as of new class I merodiploids. Analysis of junctions involved in the transfer of the trpE26 mutation, i.e. simultaneous translocation of chromosome segment C and rotation of the terminal relative to the origin moiety of the chromosome, did not confirm the presence of any sequence suitable for homologous recombination. We propose a model involving simultaneous introduction of four donor DNA molecules, each comprising a different relevant junction, and their pairing with the junction regions of the recipient chromosome. The resolution of this structure, resting on homologous recombination, would confer the donor chromosome structure to the recipient, achieving some kind of 'transstamping'. In addition, a rather regular pattern of inverse and direct short sequence repeats in regions flanking the breaking points could be correlated with the initial, X-ray-induced, rearrangement.

Phylogeography and recolonization of the Swiss Alps by the Valais shrew (Sorex antinorii), inferred with autosomal and sex-specific markers.

Using one male-inherited, one female-inherited and eight biparentally inherited markers, we investigate the population genetic structure of the Valais shrew (Sorex antinorii) in the Swiss Alps. Bayesian analysis on autosomal microsatellites suggests a clear genetic differentiation between two groups of populations. This geographically based structure is consistent with two separate postglacial recolonization routes of the species into Switzerland from Italian refugia after the last Pleistocene glaciations. Sex-specific markers also confirm genetic structuring among western and eastern areas, since very few haplotypes for either Y chromosome or mtDNA genome are shared between the two regions. Overall, these results suggest that two already well-differentiated genetic lineages colonized the Swiss Alps and came into secondary contact in the Rhône Valley. Low level of admixture between the two lineages is likely explained by the mountainous landscape structure of lateral valleys orthogonal to the main Rhône valley.

RDM1, a novel RNA recognition motif (RRM)-containing protein involved in the cell response to cisplatin in vertebrates.

A variety of cellular proteins has the ability to recognize DNA lesions induced by the anti-cancer drug cisplatin, with diverse consequences on their repair and on the therapeutic effectiveness of this drug. We report a novel gene involved in the cell response to cisplatin in vertebrates. The RDM1 gene (for RAD52 Motif 1) was identified while searching databases for sequences showing similarities to RAD52, a protein involved in homologous recombination and DNA double-strand break repair. Ablation of RDM1 in the chicken B cell line DT40 led to a more than 3-fold increase in sensitivity to cisplatin. However, RDM1-/- cells were not hypersensitive to DNA damages caused by ionizing radiation, UV irradiation, or the alkylating agent methylmethane sulfonate. The RDM1 protein displays a nucleic acid binding domain of the RNA recognition motif (RRM) type. By using gel-shift assays and electron microscopy, we show that purified, recombinant chicken RDM1 protein interacts with single-stranded DNA as well as double-stranded DNA, on which it assembles filament-like structures. Notably, RDM1 recognizes DNA distortions induced by cisplatin-DNA adducts in vitro. Finally, human RDM1 transcripts are abundant in the testis, suggesting a possible role during spermatogenesis.

Copy number variations and cognitive phenotypes in unselected populations

IMPORTANCE: The association of copy number variations (CNVs), differing numbers of copies of genetic sequence at locations in the genome, with phenotypes such as intellectual disability has been almost exclusively evaluated using clinically ascertained cohorts. The contribution of these genetic variants to cognitive phenotypes in the general population remains unclear. OBJECTIVE: To investigate the clinical features conferred by CNVs associated with known syndromes in adult carriers without clinical preselection and to assess the genome-wide consequences of rare CNVs (frequency ≤0.05%; size ≥250 kilobase pairs [kb]) on carriers' educational attainment and intellectual disability prevalence in the general population. DESIGN, SETTING, AND PARTICIPANTS: The population biobank of Estonia contains 52,000 participants enrolled from 2002 through 2010. General practitioners examined participants and filled out a questionnaire of health- and lifestyle-related questions, as well as reported diagnoses. Copy number variant analysis was conducted on a random sample of 7877 individuals and genotype-phenotype associations with education and disease traits were evaluated. Our results were replicated on a high-functioning group of 993 Estonians and 3 geographically distinct populations in the United Kingdom, the United States, and Italy. MAIN OUTCOMES AND MEASURES: Phenotypes of genomic disorders in the general population, prevalence of autosomal CNVs, and association of these variants with educational attainment (from less than primary school through scientific degree) and prevalence of intellectual disability. RESULTS: Of the 7877 in the Estonian cohort, we identified 56 carriers of CNVs associated with known syndromes. Their phenotypes, including cognitive and psychiatric problems, epilepsy, neuropathies, obesity, and congenital malformations are similar to those described for carriers of identical rearrangements ascertained in clinical cohorts. A genome-wide evaluation of rare autosomal CNVs (frequency, ≤0.05%; ≥250 kb) identified 831 carriers (10.5%) of the screened general population. Eleven of 216 (5.1%) carriers of a deletion of at least 250 kb (odds ratio [OR], 3.16; 95% CI, 1.51-5.98; P = 1.5e-03) and 6 of 102 (5.9%) carriers of a duplication of at least 1 Mb (OR, 3.67; 95% CI, 1.29-8.54; P = .008) had an intellectual disability compared with 114 of 6819 (1.7%) in the Estonian cohort. The mean education attainment was 3.81 (P = 1.06e-04) among 248 (≥250 kb) deletion carriers and 3.69 (P = 5.024e-05) among 115 duplication carriers (≥1 Mb). Of the deletion carriers, 33.5% did not graduate from high school (OR, 1.48; 95% CI, 1.12-1.95; P = .005) and 39.1% of duplication carriers did not graduate high school (OR, 1.89; 95% CI, 1.27-2.8; P = 1.6e-03). Evidence for an association between rare CNVs and lower educational attainment was supported by analyses of cohorts of adults from Italy and the United States and adolescents from the United Kingdom. CONCLUSIONS AND RELEVANCE: Known pathogenic CNVs in unselected, but assumed to be healthy, adult populations may be associated with unrecognized clinical sequelae. Additionally, individually rare but collectively common intermediate-size CNVs may be negatively associated with educational attainment. Replication of these findings in additional population groups is warranted given the potential implications of this observation for genomics research, clinical care, and public health.

The link between partial melting, granitization and granulite development in central Ribeira Fold Belt, SE Brazil: New evidence from elemental and Sr-Nd isotopic geochemistry

Elemental and Sr-Nd isotopic data on metatexites, diatexites, orthogneisses and charnockites from the central Ribeira Fold Belt indicate that they are LILE-enriched weakly peraluminous granodiorites. Harker and Th-Hf-La correlation trends suggest that these rocks represent a co-genetic sequence, whereas variations on CaO, MnO, Y and HREE for charnockites can be explained by garnet consumption during granulitic metamorphism. Similar REE patterns and isotopic results of epsilon(565)(Nd) = -5.4 to -7.3 and (87)Sr/(86)Sr(565) = 0.706-0.711 for metatexites, diatexites, orthogneisses and charnockites, as well as similar T(DM) ages between 2.0 and 1.5 Ga are consistent with evolution from a relatively homogeneous and enriched common crustal (metasedimentary) protolith. Results suggest a genetic link between metatexites, diatexites, orthogneisses and charnockites and a two-step process for charnockite development: (a) generation of the hydrated igneous protoliths by anatexis of metasedimentary rocks; (b) continuous high-grade metamorphism that transformed the ""S-type granitoids"" (leucosomes and diatexites) into orthogneisses and, as metamorphism and dehydration progressed, into charnockites. (C) 2011 Elsevier Ltd. All rights reserved.

DNA barcodes discriminate freshwater fishes from the Paraiba do Sul River Basin, São Paulo, Brazil

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Clonal relationship among invasive and non-invasive strains of enteroinvasive Escherichia coli serogroups

The genetic relatedness among 96 invasive Escherichia coli belonging to several serogroups and 13 non-invasive of several serotypes that share the same O antigen was investigated by multilocus enzyme electrophoresis analysis. The invasive strains were isolated in different parts of the world and most of them recovered from dysentery. Twenty-nine electrophoretic types were distinguished and the most invasive strains were found to belong to two major lineages. These results suggested that the invasive ability in these strains has evolved in divergent chromosomal backgrounds, presumably through the horizontal spread of plasmid-borne invasion genes. The maintenance of invasive phenotypes in separate lineages suggests that this ability confers a selective advantage to invasive strains. Copyright (C) 1999 Federation of European Microbiological Societies.

The taxonomic and phylogenetic relationships of Trypanosoma vivax from South America and Africa

The taxonomic and phylogenetic relationships of Trypanosoma vivax are controversial. It is generally suggested that South American, and East and West African isolates could be classified as subspecies or species allied to T. vivax. This is the first phylogenetic study to compare South American isolates (Brazil and Venezuela) with West/East African T. vivax isolates. Phylogeny using ribosomal sequences positioned all T. vivax isolates tightly together on the periphery of the clade containing all Salivarian trypanosomes. The same branching of isolates within T. vivax clade was observed in all inferred phylogenies using different data sets of sequences (SSU, SSU plus 5.8S or whole ITS rDNA). T. vivax from Brazil, Venezuela and West Africa (Nigeria) were closely related corroborating the West African origin of South American T. vivax, whereas a large genetic distance separated these isolates from the East African isolate (Kenya) analysed. Brazilian isolates from cattle asymptomatic or showing distinct pathology were highly homogeneous. This study did not disclose significant polymorphism to separate West African and South American isolates into different species/subspecies and indicate that the complexity of T. vivax in Africa and of the whole subgenus Trypanosoma (Duttonella) might be higher than previously believed. © 2006 Cambridge University Press.

Detection of Paracoccidioides spp. in environmental aerosol samples

Taking into account that paracoccidioidomycosis infection occurs by inhalation of the asexual conidia produced by Paracoccidioides spp. in its saprobic phase, this work presents the collection of aerosol samples as an option for environmental detection of this pathogen, by positioning a cyclonic air sampler at the entrance of armadillo burrows. Methods included direct culture, extinction technique culture and Nested PCR of the rRNA coding sequence, comprising the ITS1-5.8S-ITS2 region. In addition, we evaluated one armadillo (Dasypus novemcinctus) as a positive control for the studied area. Although the pathogen could not be isolated by the culturing strategies, the aerosol sampling associated with molecular detection through Nested PCR proved the best method for discovering Paracoccidioides spp. in the environment. Most of the ITS sequences obtained in this investigation proved to be highly similar with the homologous sequences of Paracoccidioides lutzii from the GenBank database, suggesting that this Paracoccidioides species may not be exclusive to mid-western Brazil as proposed so far. © 2013 ISHAM.