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.

Functional studies of a germ-line polymorphism at codon 47 within the p53 gene.

A rare germ-line polymorphism in codon 47 of the p53 gene replaces the wild-type proline (CCG) with a serine (TCG). Restriction analysis of 101 human samples revealed the frequency of the rare allele to be 0% (n = 69) in Caucasians and 4.7% (3/64, n = 32) among African-Americans. To investigate the consequence of this amino acid substitution, a cDNA construct (p53 mut47ser) containing the mutation was introduced into a lung adenocarcinoma cell line (Calu-6) that does not express p53. A growth suppression similar to that obtained after introduction of a wild-type p53 cDNA construct was observed, in contrast to the result obtained by introduction of p53 mut143ala. Furthermore, expression of neither p53 mut47ser nor wild-type p53 was tolerated by growing cells. In transient expression assays, both mut47ser and wild-type p53 activated the expression of a reporter gene linked to a p53 binding sequence (PG13-CAT) and inhibited the expression of the luciferase gene under the control of the Rous sarcoma virus promoter (RSVluc). In the same assay, mut143ala did not activate the expression of PG13-CAT and produced only a slight inhibitory effect on RSVluc. These findings indicate that the p53 variant with a serine at codon 47 should be considered as a rare germ-line polymorphism that does not alter the growth-suppression activity of p53.

Neuropeptide Y at the cellular level – Studies on PreproNPY L7P Polymorphism and the Mitochondrial Form of NPY

Neuropeptide Y (NPY) is a widely expressed neurotransmitter in the central and peripheral nervous systems. Thymidine 1128 to cytocine substitution in the signal sequence of the preproNPY results in a single amino acid change where leucine is changed to proline. This L7P change leads to a conformational change of the signal sequence which can have an effect on the intracellular processing of NPY. The L7P polymorphism was originally associated with higher total and LDL cholesterol levels in obese subjects. It has also been associated with several other physiological and pathophysiological responses such as atherosclerosis and T2 diabetes. However, the changes on the cellular level due to the preproNPY signal sequence L7P polymorphism were not known. The aims of the current thesis were to study the effects of the [p.L7]+[p.L7] and the [p.L7]+[p.P7] genotypes in primary cultured and genotyped human umbilical vein endothelial cells (HUVEC), in neuroblastoma (SK-N-BE(2)) cells and in fibroblast (CHO-K1) cells. Also, the putative effects of the L7P polymorphism on proliferation, apoptosis and LDL and nitric oxide metabolism were investigated. In the course of the studies a fragment of NPY targeted to mitochondria was found. With the putative mitochondrial NPY fragment the aim was to study the translational preferences and the mobility of the protein. The intracellular distribution of NPY between the [p.L7]+[p.L7] and the [p.L7]+[p.P7] genotypes was found to be different. NPY immunoreactivity was prominent in the [p.L7]+[p.P7] cells while the proNPY immunoreactivity was prominent in the [p.L7]+[p.L7] genotype cells. In the proliferation experiments there was a difference in the [p.L7]+[p.L7] genotype cells between early and late passage (aged) cells; the proliferation was raised in the aged cells. NPY increased the growth of the cells with the [p.L7]+[p.P7] genotype. Apoptosis did not seem to differ between the genotypes, but in the aged cells with the [p.L7]+[p.L7] genotype, LDL uptake was found to be elevated. Furthermore, the genotype seemed to have a strong effect on the nitric oxide metabolism. The results indicated that the mobility of NPY protein inside the cells was increased within the P7 containing constructs. The existence of the mitochondria targeted NPY fragment was verified, and translational preferences were proved to be due to the origin of the cells. Cell of neuronal origin preferred the translation of mature NPY (NPY1-36) when compared to the non neuronal cells that translated both, NPY and the mitochondrial fragment of NPY. The mobility of the mitochondrial fragment was found to be minimal. The functionality of the mitochondrial NPY fragment remains to be investigated. L7P polymorphism in the preproNPY causes a series of intracellular changes. These changes may contribute to the state of cellular senescence, vascular tone and lead to endothelial dysfunction and even to increased susceptibility to diseases, like atherosclerosis and T2 diabetes.

Characterization of six rat strains (Rattus norvegicus) by mitochondrial DNA restriction fragment length polymorphism

Restriction fragment length polymorphism (RFLP) was used to examine the extent of mtDNA polymorphism among six strains of rats (Rattus norvegicus) - Wistar, Wistar Munich, Brown Norway, Wistar Kyoto, SHR and SHR-SP. A survey of 26 restriction enzymes has revealed a low level of genetic divergence among strains. The sites of cleavage by EcoRI, NcoI and XmnI were shown to be polymorphic. The use of these three enzymes allows the 6 strains to be classified into 4 haplotypes and identifies specific markers for each one. The percentage of sequence divergence among all pairs of haplotypes ranged from 0.035 to 0.33%, which is the result of a severe population constriction undergone by the strains. These haplotypes are easily demonstrable and therefore RFLP analysis can be employed for genetic monitoring of rats within animal facilities or among different laboratories.

The P48T germline mutation and polymorphism in the CDKN2A gene of patients with melanoma

CDKN2A has been implicated as a melanoma susceptibility gene in some kindreds with a family history of this disease. Mutations in CDKN2A may produce an imbalance between functional p16ink4a and cyclin D causing abnormal cell growth. We searched for germline mutations in this gene in 22 patients with clinical criteria of hereditary cancer (early onset, presence of multiple primary melanoma or 1 or more first- or second-degree relatives affected) by secondary structural content prediction, a mutation scanning method that relies on the propensity for single-strand DNA to take on a three-dimensional structure that is highly sequence dependent, and sequencing the samples with alterations in the electrophoretic mobility. The prevalence of CDKN2A mutation in our study was 4.5% (1/22) and there was a correlation between family history and probability of mutation detection. We found the P48T mutation in 1 patient with 2 melanoma-affected relatives. The patient descends from Italian families and this mutation has been reported previously only in Italian families in two independent studies. This leads us to suggest the presence of a mutational "hotspot" within this gene or a founder mutation. We also detected a high prevalence (59.1%) of polymorphisms, mainly alleles 500 C/G (7/31.8%) or 540 C/T (6/27.3%), in the 3' untranslated region of exon 3. This result reinforces the idea that these rare polymorphic alleles have been significantly associated with the risk of developing melanoma.

Polymorphism of the promoter region and exon 1 of the CTLA4 gene in endemic pemphigus foliaceus (fogo selvagem)

Endemic pemphigus foliaceus (EPF) is an autoimmune bullous skin disease characterized by acantholysis and antibodies against a desmosomal protein, desmoglein 1. Genetic and environmental factors contribute to development of this multifactorial disease. HLA class II and some cytokine gene polymorphisms are the only genetic markers thus far known to be associated with susceptibility to or protection from EPF. The cytotoxic T-lymphocyte antigen-4 gene (CTLA4) encodes a key immunoreceptor molecule that regulates and inhibits T-cell proliferation. It participates in the regulatory process controlling autoreactivity and therefore has been considered a strong candidate gene in autoimmune diseases. In the search for genes that might influence EPF pathogenesis, we analyzed variants of the CTLA4 gene in a sample of 118 patients and 291 controls from a Brazilian population. This is the first study investigating the possible role of polymorphisms of the 2q33 chromosomal region in differential susceptibility to pemphigus foliaceus. Promoter region and exon 1 single nucleotide polymorphisms -318 (C,T) and 49 (A,G) were genotyped using sequence-specific oligonucleotide probes after amplification by the polymerase chain reaction. The allelic and genotypic frequencies did not differ significantly between the patient and the control groups (-318T: 9.8 and 10.9%, 49G: 33.0 and 35.2% were the allelic frequencies in patients and controls, respectively). In addition, no significant difference was found when the patient and control population samples were stratified by the presence of HLA-DRB1 alleles. We conclude that the CTLA4 -318 (C,T) and 49 (A,G) polymorphisms do not play a major role in EPF development.

No evidence of association of MUC-1 genetic polymorphism with embryo implantation failure

Pregnancy loss can be caused by several factors involved in human reproduction. Although up to 50% of cases remain unexplained, it has been postulated that the major cause of failed pregnancy is an error of embryo implantation. Transmembrane mucin-1 (MUC-1) is a glycoprotein expressed on the endometrial cell surface which acts as a barrier to implantation. The gene that codes for this molecule is composed of a polymorphic tandem repeat of 60 nucleotides. Our objective was to determine if MUC-1 genetic polymorphism is associated with implantation failure in patients with a history of recurrent abortion. The study was conducted on 10 women aged 25 to 35 years with no history of successful pregnancy and with a diagnosis of infertility. The control group consisted of 32 patients aged 25 to 35 years who had delivered at least two full-term live children and who had no history of abortions or fetal losses. MUC-1 amplicons were obtained by PCR and observed on agarose and polyacrylamide gel after electrophoresis. Statistical analysis showed no significant difference in the number of MUC-1 variable number of tandem repeats between these groups (P > 0.05). Our results suggest that there is no effect of the polymorphic MUC-1 sequence on the implantation failure. However, the data do not exclude MUC-1 relevance during embryo implantation. The process is related to several associated factors such as the mechanisms of gene expression in the uterus, specific MUC-1 post-translational modifications and appropriate interactions with other molecules during embryo implantation.

Apolipoprotein CIII polymorphism and triglyceride levels of a Japanese population living in Southern Brazil

Apolipoprotein CIII (apo-CIII) participates in the regulation of triglyceride-rich lipoprotein metabolism. Several polymorphic sites have been detected within and around the apo-CIII gene. Here, we examined the relationship between apo-CIII SstI polymorphism (CC, CG, GG genotypes) and plasma triglyceride (TG) levels in a group of 159 Japanese individuals living in Southern Brazil. The sample was divided into a group of Japanese descendants (N = 51) with high TG (HTG; >200 mg/dL) and a group of Japanese descendants (N = 108) with normal TG (NTG; <200 mg/dL). TG and total cholesterol levels were analyzed by an enzymatic method using the Labtest-Diagnostic kit and high- and low-density lipoproteins by a direct method using the Labtest-Diagnostic kit and DiaSys Diagnostic System International kit, respectively. A 428-bp sequence of apo-CIII gene was amplified using oligonucleotide primers 5' GGT GAC CGA TGG CTT CAG TTC CCT GA 3' and 5' CAG AAG GTG GAT AGA GCG CTG GCC T 3'. The PCR products were digested with a restriction endonuclease SstI. Rare G allele was highly prevalent in our study population (0.416) compared to Caucasians (0.00-0.11). G allele was almost two times more prevalent in the HTG group compared to the NTG group (P < 0.001). The genotype distribution was consistent with the Hardy-Weinberg equilibrium. There was a significant association between rare G allele and HTG in Japanese individuals living in Southern Brazil as indicated by one-way ANOVA, P < 0.05.

Correlation between interleukin-18promoter -607C/A polymorphism and susceptibility to ischemic stroke

Single nucleotide polymorphisms in the promoter region ofinterleukin-18 (IL-18), an inflammatory cytokine, have been linked to susceptibility to many diseases, including cancer and immune dysfunction. Here, we explored the potential association between theIL-18 -607C/A (rs1946518) promoter region polymorphism and susceptibility to ischemic stroke (IS). This locus was amplified from peripheral blood samples of 386 IS patients (cases) and 364 healthy individuals (controls) by the polymerase chain reaction with sequence-specific primers. Significant differences were observed by the χ2 test in the -607C/A (rs1946518) genotype and allele frequencies between cases and controls (P < 0.05). Furthermore, after excluding for age, gender, smoking status, and hypertension, logistic regression indicated that IS susceptibility of -607C carriers increased 1.6 times (OR = 1.601, 95%CI = 1.148-2.233, P = 0.006) compared to -607A carriers. Additionally, similar increases in IS risk were noted for male patients or patients less than 65 years old. In conclusion,IL-18 -607C/A (rs1946518) promoter polymorphism is associated with IS susceptibility, and the C allele may confer increased IS risk.

DNA polymorphism in Cab locus of tomato induced by tissue culture

Plants were regenerated from callus induced from leaf disc explants of a tomato F, hybrid heterozygous for three marker loci (a), without anthocyanin (aw), and hairless (hl). Regenerants were studied for somaclonal variation at the phenotypic level by scoring for variation in the marker loci, and at the DNA level by probing geomic DNA blots with a chlorophyll a/b binding protein (Cab-3C) cDNA sequence. While no variation was observed at the phenotypic level in over 950 somaclones studied, DNA polymorphism for the Cab locus could be detected in two out of 17 somaclones tested. Tissue culture induced variation at the phenotypic level for specific loci is very low (less than 0.001 for a, awor hl) but DNA sequence changes are induced at much greater frequency (- 0.1 for a multicopy gene family such as Cab).

Plasmodium vivax lineages: geographical distribution, tandem repeat polymorphism, and phylogenetic relationship

Background: Multi-drug resistance and severe/ complicated cases are the emerging phenotypes of vivax malaria, which may deteriorate current anti-malarial control measures. The emergence of these phenotypes could be associated with either of the two Plasmodium vivax lineages. The two lineages had been categorized as Old World and New World, based on geographical sub-division and genetic and phenotypical markers. This study revisited the lineage hypothesis of P. vivax by typing the distribution of lineages among global isolates and evaluated their genetic relatedness using a panel of new mini-satellite markers. Methods: 18S SSU rRNA S-type gene was amplified from 420 Plasmodium vivax field isolates collected from different geographical regions of India, Thailand and Colombia as well as four strains each of P. vivax originating from Nicaragua, Panama, Thailand (Pak Chang), and Vietnam (ONG). A mini-satellite marker panel was then developed to understand the population genetic parameters and tested on a sample subset of both lineages. Results: 18S SSU rRNA S-type gene typing revealed the distribution of both lineages (Old World and New World) in all geographical regions. However, distribution of Plasmodium vivax lineages was highly variable in every geographical region. The lack of geographical sub-division between lineages suggests that both lineages are globally distributed. Ten mini-satellites were scanned from the P. vivax genome sequence; these tandem repeats were located in eight of the chromosomes. Mini-satellites revealed substantial allelic diversity (7-21, AE = 14.6 +/- 2.0) and heterozygosity (He = 0.697-0.924, AE = 0.857 +/- 0.033) per locus. Mini-satellite comparison between the two lineages revealed high but similar pattern of genetic diversity, allele frequency, and high degree of allele sharing. A Neighbour-Joining phylogenetic tree derived from genetic distance data obtained from ten mini-satellites also placed both lineages together in every cluster. Conclusions: The global lineage distribution, lack of genetic distance, similar pattern of genetic diversity, and allele sharing strongly suggested that both lineages are a single species and thus new emerging phenotypes associated with vivax malaria could not be clearly classified as belonging to a particular lineage on basis of their geographical origin.

Detection of somaclonal variation during cocoa somatic embryogenesis characterised using cleaved amplified polymorphic sequence and the new freeware Artbio

The scarcity and stochastic nature of genetic mutations presents a significant challenge for scientists seeking to characterise de novo mutation frequency at specific loci. Such mutations can be particularly numerous during regeneration of plants from in vitro culture and can undermine the value of germplasm conservation efforts. We used cleaved amplified polymorphic sequence (CAPS) analysis to characterise new mutations amongst a clonal population of cocoa plants regenerated via a somatic embryogenesis protocol used previously for cocoa cryopreservation. Efficacy of the CAPS system for mutation detection was greatly improved after an ‘a priori’ in silico screen of reference target sequences for actual and potential restriction enzyme recognition sites using a new freely available software called Artbio. Artbio surveys known sequences for existing restriction enzyme recognition sites but also identifies all single nucleotide polymorphism (SNP) deviations from such motifs. Using this software, we performed an in silico screen of seven loci for restriction sites and their potential mutant SNP variants that were possible from 21 restriction enzymes. The four most informative locus-enzyme combinations were then used to survey the regenerant populations for de novo mutants. We characterised the pattern of point mutations and, using the outputs of Artbio, calculated the ratio of base substitution in 114 somatic embryo-derived cocoa regenerants originating from two explant genotypes. We found 49 polymorphisms, comprising 26.3% of the samples screened, with an inferred rate of 2.8 × 10−3 substitutions/screened base. This elevated rate is of a similar order of magnitude to previous reports of de novo microsatellite length mutations arising in the crop and suggests caution should be exercised when applying somatic embryogenesis for the conservation of plant germplasm.

Sequence heterogeneity of an mpb70 gene analogue in Mycobacterium kansasii

The protein antigen MPB70 is a major component of culture supernatants of Mycobacterium bovis and Is an active ingredient of bovine PPD used for skin-testing cattle for tuberculosis. we have shown that Mycobacterium kansasii possesses a similar gene that cross-reacts in a PCR test for M. bovis. Single strand conformational polymorphism analysis, and the DNA sequence of the PCR product, shows differences between M. kansasii strains, supporting the suggestion that M. kansasii is not a homogeneous species.

Analysis of DNA polymorphism in ancient barley herbarium material: validation of the KASP SNP genotyping platform

The availability of crop specimens archived in herbaria and old seed collections represent valuable resources for the analysis of plant genetic diversity and crop domestication. The ability to extract ancient DNA (aDNA) from such samples has recently allowed molecular genetic investigations to be undertaken in ancient materials. While analyses of aDNA initially focused on the use of markers which occur in multiple copies such as the internal transcribed spacer region (ITS) within ribosomal DNA and those requiring amplification of short DNA regions of variable length such as simple sequence repeats (SSRs), emphasis is now moving towards the genotyping of single nucleotide polymorphisms (SNPs), traditionally undertaken in aDNA by Sanger sequencing. Here, using a panel of barley aDNA samples previously surveyed by Sanger sequencing for putative causative SNPs within the flowering-time gene PPD-H1, we assess the utility of the Kompetitive Allele Specific PCR (KASP) genotyping platform for aDNA analysis. We find KASP to out-perform Sanger sequencing in the genotyping of aDNA samples (78% versus 61% success, respectively), as well as being robust to contamination. The small template size (≥46 bp) and one-step, closed-tube amplification/genotyping process make this platform ideally suited to the genotypic analysis of aDNA, a process which is often hampered by template DNA degradation and sample cross-contamination. Such attributes, as well as its flexibility of use and relatively low cost, make KASP particularly relevant to the genetic analysis of aDNA samples. Furthermore, KASP provides a common platform for the genotyping and analysis of corresponding SNPs in ancient, landrace and modern plant materials. The extended haplotype analysis of PPD-H1 undertaken here (allelic variation at which is thought to be important for the spread of domestication and local adaptation) provides further resolution to the previously identified geographic cline of flowering-time allele distribution, illustrating how KASP can be used to aid genetic analyses of aDNA from plant species. We further demonstrate the utility of KASP by genotyping ten additional genetic markers diagnostic for morphological traits in barley, shedding light on the phenotypic traits, alleles and allele combinations present in these unviable ancient specimens, as well as their geographic distributions.