Identifying Selected Regions from Heterozygosity and Divergence Using a Light-Coverage Genomic Dataset from Two Human Populations

When a selective sweep occurs in the chromosomal region around a target gene in two populations that have recently separated, it produces three dramatic genomic consequences: 1) decreased multi-locus heterozygosity in the region; 2) elevated or diminished genetic divergence (FST) of multiple polymorphic variants adjacent to the selected locus between the divergent populations, due to the alternative fixation of alleles; and 3) a consequent regional increase in the variance of FST (S2FST) for the same clustered variants, due to the increased alternative fixation of alleles in the loci surrounding the selection target. In the first part of our study, to search for potential targets of directional selection, we developed and validated a resampling-based computational approach; we then scanned an array of 31 different-sized moving windows of SNP variants (5–65 SNPs) across the human genome in a set of European and African American population samples with 183,997 SNP loci after correcting for the recombination rate variation. The analysis revealed 180 regions of recent selection with very strong evidence in either population or both. In the second part of our study, we compared the newly discovered putative regions to those sites previously postulated in the literature, using methods based on inspecting patterns of linkage disequilibrium, population divergence and other methodologies. The newly found regions were cross-validated with those found in nine other studies that have searched for selection signals. Our study was replicated especially well in those regions confirmed by three or more studies. These validated regions were independently verified, using a combination of different methods and different databases in other studies, and should include fewer false positives. The main strength of our analysis method compared to others is that it does not require dense genotyping and therefore can be used with data from population-based genome SNP scans from smaller studies of humans or other species

SUP: an extension to SLINK to allow a larger number of marker loci to be simulated in pedigrees conditional on trait values

BACKGROUND: With the recent advances in high-throughput genotyping technologies that allow for large-scale association mapping of human complex traits, promising statistical designs and methods have been emerging. Efficient simulation software are key elements for the evaluation of the properties of new statistical tests. SLINK is a flexible simulation tool that has been widely used to generate the segregation and recombination processes of markers linked to, and possibly associated with, a trait locus, conditional on trait values in arbitrary pedigrees. In practice, its most serious limitation is the small number of loci that can be simulated, since the complexity of the algorithm scales exponentially with this number. RESULTS: I describe the implementation of a two-step algorithm to be used in conjunction with SLINK to enable the simulation of a large number of marker loci linked to a trait locus and conditional on trait values in families, with the possibility for the loci to be in linkage disequilibrium. SLINK is used in the first step to simulate genotypes at the trait locus conditional on the observed trait values, and also to generate an indicator of the descent path of the simulated alleles. In the second step, marker alleles or haplotypes are generated in the founders, conditional on the trait locus genotypes simulated in the first step. Then the recombination process between the marker loci takes place conditionally on the descent path and on the trait locus genotypes. This two-step implementation is often computationally faster than other software that are designed to generate marker data linked to, and possibly associated with, a trait locus. CONCLUSION: Because the proposed method uses SLINK to simulate the segregation process, it benefits from its flexibility: the trait may be qualitative with the possibility of defining different liability classes (which allows for the simulation of gene-environment interactions or even the simulation of multi-locus effects between unlinked susceptibility regions) or it may be quantitative and normally distributed. In particular, this implementation is the only one available that can generate a large number of marker loci conditional on the set of observed quantitative trait values in pedigrees

PGA: power calculator for case-control genetic association analyses

<p>Abstract</p> <p>Background</p> <p>Statistical power calculations inform the design and interpretation of genetic association studies, but few programs are tailored to case-control studies of single nucleotide polymorphisms (SNPs) in unrelated subjects.</p> <p>Results</p> <p>We have developed the "Power for Genetic Association analyses" (PGA) package which comprises algorithms and graphical user interfaces for sample size and minimum detectable risk calculations using SNP or haplotype effects under different genetic models and study constrains. The software accounts for linkage disequilibrium and statistical multiple comparisons. The results are presented in graphs or tables and can be printed or exported in standard file formats.</p> <p>Conclusion</p> <p>PGA is user friendly software that can facilitate decision making for association studies of candidate genes, fine-mapping studies, and whole-genome scans. Stand-alone executable files and a Matlab toolbox are available for download at: <url>http://dceg.cancer.gov/bb/tools/pga</url></p

Identification of field caught Anopheles gambiae s.s. and Anopheles arabiensis by TaqMan single nucleotide polymorphism genotyping

BACKGROUND: Identification of Anopheles gambiae s.s. and Anopheles arabiensis from field-collected Anopheles gambiae s.l. is often necessary in basic and applied research, and in operational control programmes. The currently accepted method involves use of standard polymerase chain reaction amplification of ribosomal DNA (rDNA) from the 3' 28S to 5' intergenic spacer region of the genome, and visual confirmation of amplicons of predicted size on agarose gels, after electrophoresis. This report describes development and evaluation of an automated, quantitative PCR method based upon TaqMan™ single nucleotide polymorphism (SNP) genotyping. METHODS: Standard PCR, and TaqMan SNP genotyping with newly designed primers and fluorophore-labeled probes hybridizing to sequences of complementary rDNA specific for either An. gambiae s.s. or An. arabiensis, were conducted in three experiments involving field-collected An. gambiae s.l. from western Kenya, and defined laboratory strains. DNA extraction was from a single leg, sonicated for five minutes in buffer in wells of 96-well PCR plates. RESULTS: TaqMan SNP genotyping showed a reaction success rate, sensitivity, and species specificity comparable to that of standard PCR. In an extensive field study, only 29 of 3,041 (0.95%) were determined to be hybrids by TaqMan (i.e., having rDNA sequences from both species), however, all but one were An. arabiensis by standard PCR, suggesting an acceptably low (ca. 1%) error rate for TaqMan genotyping in mistakenly identifying species hybrids. CONCLUSION: TaqMan SNP genotyping proved to be a sensitive and rapid method for identification of An. gambiae s.l. and An. arabiensis, with a high success rate, specific results, and congruence with the standard PCR method

Strong signature of natural selection within an FHIT intron implicated in prostate cancer risk

Previously, a candidate gene linkage approach on brother pairs affected with prostate cancer identified a locus of prostate cancer susceptibility at D3S1234 within the fragile histidine triad gene (FHIT), a tumor suppressor that induces apoptosis. Subsequent association tests on 16 SNPs spanning approximately 381 kb surrounding D3S1234 in Americans of European descent revealed significant evidence of association for a single SNP within intron 5 of FHIT. In the current study, resequencing and genotyping within a 28.5 kb region surrounding this SNP further delineated the association with prostate cancer risk to a 15 kb region. Multiple SNPs in sequences under evolutionary constraint within intron 5 of FHIT defined several related haplotypes with an increased risk of prostate cancer in European-Americans. Strong associations were detected for a risk haplotype defined by SNPs 138543, 142413, and 152494 in all cases (Pearson's χ2 = 12.34, df 1, P = 0.00045) and for the homozygous risk haplotype defined by SNPs 144716, 142413, and 148444 in cases that shared 2 alleles identical by descent with their affected brothers (Pearson's χ2 = 11.50, df 1, P = 0.00070). In addition to highly conserved sequences encompassing SNPs 148444 and 152413, population studies revealed strong signatures of natural selection for a 1 kb window covering the SNP 144716 in two human populations, the European American (π = 0.0072, Tajima's D= 3.31, 14 SNPs) and the Japanese (π = 0.0049, Fay & Wu's H = 8.05, 14 SNPs), as well as in chimpanzees (Fay & Wu's H = 8.62, 12 SNPs). These results strongly support the involvement of the FHIT intronic region in an increased risk of prostate cancer. © 2008 Ding et al

Genes encoding critical transcriptional activators for murine neural tube development and human spina bifida: a case-control study

<p>Abstract</p> <p>Background</p> <p>Spina bifida is a malformation of the neural tube and is the most common of neural tube defects (NTDs). The etiology of spina bifida is largely unknown, although it is thought to be multi-factorial, involving multiple interacting genes and environmental factors. Mutations in transcriptional co-activator genes-<it>Cited2</it>, <it>p300</it>, <it>Cbp</it>, <it>Tfap2α</it>, <it>Carm1 </it>and <it>Cart1 </it>result in NTDs in murine models, thus prompt us to investigate whether homologues of these genes are associated with NTDs in humans.</p> <p>Methods</p> <p>Data and biological samples from 297 spina bifida cases and 300 controls were derived from a population-based case-control study conducted in California. 37 SNPs within <it>CITED2</it>, <it>EP300</it>, <it>CREBBP</it>, <it>TFAP2A</it>, <it>CARM1 </it>and <it>ALX1 </it>were genotyped using an ABI SNPlex assay. Odds ratios and 95% confidence intervals were calculated for alleles, genotypes and haplotypes to evaluate the risk for spina bifida.</p> <p>Results</p> <p>Several SNPs showed increased or decreased risk, including <it>CITED2 </it>rs1131431 (OR = 5.32, 1.04~27.30), <it>EP300 </it>rs4820428 (OR = 1.30, 1.01~1.67), <it>EP300 </it>rs4820429 (OR = 0.50, 0.26~0.50, in whites, OR = 0.7, 0.49~0.99 in all subjects), <it>EP300 </it>rs17002284 (OR = 0.43, 0.22~0.84), <it>TFAP2A </it>rs3798691 (OR = 1.78, 1.13~2.87 in Hispanics), <it>CREBBP </it>rs129986 (OR = 0.27, 0.11~0.69), <it>CARM1 </it>rs17616105 (OR = 0.41, 0.22~0.72 in whites). In addition, one haplotype block in <it>EP300 </it>and one in <it>TFAP2A </it>appeared to be associated with increased risk.</p> <p>Conclusions</p> <p>Modest associations were observed in <it>CITED2</it>, <it>EP300</it>, <it>CREBBP</it>, <it>TFAP2A </it>and <it>CARM1 </it>but not <it>ALX1</it>. However, these modest associations were not statistically significant after correction for multiple comparisons. Searching for potential functional variants and rare causal mutations is warranted in these genes.</p

Cell walls of the dimorphic fungal pathogens Sporothrix schenckii and Sporothrix brasiliensis exhibit bilaminate structures and sloughing of extensive and intact layers

This work was supported by the Fundação Carlos Chagas de Amparo à Pesquisa do Estado do Rio de Janeiro (FAPERJ), grants E-26/202.974/2015 and Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq), grants 229755/2013-5, Brazil. LMLB is a senior research fellow of CNPq and Faperj. NG acknowledged support from the Wellcome Trust (Trust (097377, 101873, 200208) and MRC Centre for Medical Mycology (MR/N006364/1). The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.Peer reviewedPublisher PD

Promoter polymorphism -119C/G in MYG1 (C12orf10) gene is related to vitiligo susceptibility and Arg4Gln affects mitochondrial entrance of Myg1

<p>Abstract</p> <p>Background</p> <p><it>MYG1 </it>(<it>Melanocyte proliferating gene 1</it>, also C12orf10 in human) is a ubiquitous nucleo-mitochondrial protein, involved in early developmental processes and in adult stress/illness conditions. We recently showed that <it>MYG1 </it>mRNA expression is elevated in the skin of vitiligo patients. Our aim was to examine nine known polymorphisms in the <it>MYG1 </it>gene, to investigate their functionality, and to study their association with vitiligo susceptibility.</p> <p>Methods</p> <p>Nine single nucleotide polymorphisms (SNPs) in the <it>MYG1 </it>locus were investigated by SNPlex assay and/or sequencing in vitiligo patients (n = 124) and controls (n = 325). <it>MYG1 </it>expression in skin biopsies was detected by quantitative-real time PCR (Q-RT-PCR) and polymorphisms were further analysed using luciferase and YFP reporters in the cell culture.</p> <p>Results</p> <p>Control subjects with -119G promoter allele (rs1465073) exhibited significantly higher <it>MYG1 </it>mRNA levels than controls with -119C allele (<it>P </it>= 0.01). Higher activity of -119G promoter was confirmed by luciferase assay. Single marker association analysis showed that the -119G allele was more frequent in vitiligo patients (47.1%) compared to controls (39.3%, <it>P </it>< 0.05, OR 1.37, 95%CI 1.02-1.85). Analysis based on the stage of progression of the vitiligo revealed that the increased frequency of -119G allele occurred prevalently in the group of patients with active vitiligo (n = 86) compared to the control group (48.2% <it>versus </it>39.3%, <it>P </it>< 0.05; OR 1.44, 95%CI 1.02-2.03). Additionally, we showed that glutamine in the fourth position (in Arg4Gln polymorphism) completely eliminated mitochondrial entrance of YFP-tagged Myg1 protein in cell culture. The analysis of available EST, cDNA and genomic DNA sequences revealed that Myg1 4Gln allele is remarkably present in human populations but is never detected in homozygous state according to the HapMap database.</p> <p>Conclusions</p> <p>Our study demonstrated that both <it>MYG1 </it>promoter polymorphism -119C/G and Arg4Gln polymorphism in the mitochondrial signal of Myg1 have a functional impact on the regulation of the <it>MYG1 </it>gene and promoter polymorphism (-119C/G) is related with suspectibility for actively progressing vitiligo.</p

Novel variants in the PRDX6 Gene and the risk of Acute Lung Injury following major trauma

<p>Abstract</p> <p>Background</p> <p>Peroxiredoxin 6 (<it>PRDX6</it>) is involved in redox regulation of the cell and is thought to be protective against oxidant injury. Little is known about genetic variation within the PRDX6 gene and its association with acute lung injury (ALI). In this study we sequenced the <it>PRDX6 </it>gene to uncover common variants, and tested association with ALI following major trauma.</p> <p>Methods</p> <p>To examine the extent of variation in the <it>PRDX6 </it>gene, we performed direct sequencing of the 5' UTR, exons, introns and the 3' UTR in 25 African American cases and controls and 23 European American cases and controls (selected from a cohort study of major trauma), which uncovered 80 SNPs. <it>In silico </it>modeling was performed using Patrocles and Transcriptional Element Search System (TESS). Thirty seven novel and tagging SNPs were tested for association with ALI compared with ICU at-risk controls who did not develop ALI in a cohort study of 259 African American and 254 European American subjects that had been admitted to the ICU with major trauma.</p> <p>Results</p> <p>Resequencing of critically ill subjects demonstrated 43 novel SNPs not previously reported. Coding regions demonstrated no detectable variation, indicating conservation of the protein. Block haplotype analyses reveal that recombination rates within the gene seem low in both Caucasians and African Americans. Several novel SNPs appeared to have the potential for functional consequence using <it>in silico </it>modeling. Chi²analysis of ALI incidence and genotype showed no significant association between the SNPs in this study and ALI. Haplotype analysis did not reveal any association beyond single SNP analyses.</p> <p>Conclusions</p> <p>This study revealed novel SNPs within the <it>PRDX6 </it>gene and its 5' and 3' flanking regions via direct sequencing. There was no association found between these SNPs and ALI, possibly due to a low sample size, which was limited to detection of relative risks of 1.93 and above. Future studies may focus on the role of <it>PRDX6 </it>genetic variation in other diseases, where oxidative stress is suspected.</p