Associations of HPV-16 Gene Methylation With Oral HPV-16 Persistence Among a Multinational Sample of Men.

Using data from the Human Papillomavirus (HPV) Infection in Men cohort study, we demonstrate HPV-16 methylation associations with persistent oral HPV infection, the obligate precursor to oropharyngeal cancer. Human papillomavirus type 16 persistence was significantly associated with methylation of HPV-16 L2 CpG-4268 (Wilcoxon P = 0.04), and methylation of HPV-16 E2 CpG Pos 4 (Wilcoxon P = 0.04)

Inference of population splits and mixtures from genome-wide allele frequency data

Many aspects of the historical relationships between populations in a species are reflected in genetic data. Inferring these relationships from genetic data, however, remains a challenging task. In this paper, we present a statistical model for inferring the patterns of population splits and mixtures in multiple populations. In this model, the sampled populations in a species are related to their common ancestor through a graph of ancestral populations. Using genome-wide allele frequency data and a Gaussian approximation to genetic drift, we infer the structure of this graph. We applied this method to a set of 55 human populations and a set of 82 dog breeds and wild canids. In both species, we show that a simple bifurcating tree does not fully describe the data; in contrast, we infer many migration events. While some of the migration events that we find have been detected previously, many have not. For example, in the human data we infer that Cambodians trace approximately 16% of their ancestry to a population ancestral to other extant East Asian populations. In the dog data, we infer that both the boxer and basenji trace a considerable fraction of their ancestry (9% and 25%, respectively) to wolves subsequent to domestication, and that East Asian toy breeds (the Shih Tzu and the Pekingese) result from admixture between modern toy breeds and "ancient" Asian breeds. Software implementing the model described here, called TreeMix, is available at http://treemix.googlecode.comComment: 28 pages, 6 figures in main text. Attached supplement is 22 pages, 15 figures. This is an updated version of the preprint available at http://precedings.nature.com/documents/6956/version/

Bayesian estimates of linkage disequilibrium

[Background] The maximum likelihood estimator of D' – a standard measure of linkage disequilibrium – is biased toward disequilibrium, and the bias is particularly evident in small samples and rare haplotypes. [Results] This paper proposes a Bayesian estimation of D' to address this problem. The reduction of the bias is achieved by using a prior distribution on the pair-wise associations between single nucleotide polymorphisms (SNP)s that increases the likelihood of equilibrium with increasing physical distances between pairs of SNPs. We show how to compute the Bayesian estimate using a stochastic estimation based on MCMC methods, and also propose a numerical approximation to the Bayesian estimates that can be used to estimate patterns of LD in large datasets of SNPs. [Conclusion] Our Bayesian estimator of D' corrects the bias toward disequilibrium that affects the maximum likelihood estimator. A consequence of this feature is a more objective view about the extent of linkage disequilibrium in the human genome, and a more realistic number of tagging SNPs to fully exploit the power of genome wide association studies.Research supported by NIH/NHLBI grant R21 HL080463-01, NIH/NIDDK 1R01DK069646-01A1 and the Spanish research program [projects TIN2004-06204-C03-02 and TIN2005-02516]

Empirical vs Bayesian approach for estimating haplotypes from genotypes of unrelated individuals

BACKGROUND: The completion of the HapMap project has stimulated further development of haplotype-based methodologies for disease associations. A key aspect of such development is the statistical inference of individual diplotypes from unphased genotypes. Several methodologies for inferring haplotypes have been developed, but they have not been evaluated extensively to determine which method not only performs well, but also can be easily incorporated in downstream haplotype-based association analyses. In this paper, we attempt to do so. Our evaluation was carried out by comparing the two leading Bayesian methods, implemented in PHASE and HAPLOTYPER, and the two leading empirical methods, implemented in PL-EM and HPlus. We used these methods to analyze real data, namely the dense genotypes on X-chromosome of 30 European and 30 African trios provided by the International HapMap Project, and simulated genotype data. Our conclusions are based on these analyses. RESULTS: All programs performed very well on X-chromosome data, with an average similarity index of 0.99 and an average prediction rate of 0.99 for both European and African trios. On simulated data with approximation of coalescence, PHASE implementing the Bayesian method based on the coalescence approximation outperformed other programs on small sample sizes. When the sample size increased, other programs performed as well as PHASE. PL-EM and HPlus implementing empirical methods required much less running time than the programs implementing the Bayesian methods. They required only one hundredth or thousandth of the running time required by PHASE, particularly when analyzing large sample sizes and large umber of SNPs. CONCLUSION: For large sample sizes (hundreds or more), which most association studies require, the two empirical methods might be used since they infer the haplotypes as accurately as any Bayesian methods and can be incorporated easily into downstream haplotype-based analyses such as haplotype-association analyses

Genome sequencing of the extinct Eurasian wild aurochs, Bos primigenius, illuminates the phylogeography and evolution of cattle

Background Domestication of the now-extinct wild aurochs, Bos primigenius, gave rise to the two major domestic extant cattle taxa, B. taurus and B. indicus. While previous genetic studies have shed some light on the evolutionary relationships between European aurochs and modern cattle, important questions remain unanswered, including the phylogenetic status of aurochs, whether gene flow from aurochs into early domestic populations occurred, and which genomic regions were subject to selection processes during and after domestication. Here, we address these questions using whole-genome sequencing data generated from an approximately 6,750-year-old British aurochs bone and genome sequence data from 81 additional cattle plus genome-wide single nucleotide polymorphism data from a diverse panel of 1,225 modern animals. Results Phylogenomic analyses place the aurochs as a distinct outgroup to the domestic B. taurus lineage, supporting the predominant Near Eastern origin of European cattle. Conversely, traditional British and Irish breeds share more genetic variants with this aurochs specimen than other European populations, supporting localized gene flow from aurochs into the ancestors of modern British and Irish cattle, perhaps through purposeful restocking by early herders in Britain. Finally, the functions of genes showing evidence for positive selection in B. taurus are enriched for neurobiology, growth, metabolism and immunobiology, suggesting that these biological processes have been important in the domestication of cattle. Conclusions This work provides important new information regarding the origins and functional evolution of modern cattle, revealing that the interface between early European domestic populations and wild aurochs was significantly more complex than previously thought

Haplotype block structure study of the CFTR gene. Most variants are associated with the M470 allele in several European populations

An average of about 1700 CFTR (cystic fibrosis transmembrane conductance regulator) alleles from normal individuals from different European populations were extensively screened for DNA sequence variation. A total of 80 variants were observed: 61 coding SNSs (results already published), 13 noncoding SNSs, three STRs, two short deletions, and one nucleotide insertion. Eight DNA variants were classified as non-CF causing due to their high frequency of occurrence. Through this survey the CFTR has become the most exhaustively studied gene for its coding sequence variability and, though to a lesser extent, for its noncoding sequence variability as well. Interestingly, most variation was associated with the M470 allele, while the V470 allele showed an 'extended haplotype homozygosity' (EHH). These findings make us suggest a role for selection acting either on the M470V itself or through an hitchhiking mechanism involving a second site. The possible ancient origin of the V allele in an 'out of Africa' time frame is discussed

The landscape of Neandertal ancestry in present-day humans

Analyses of Neandertal genomes have revealed that Neandertals have contributed genetic variants to modern humans1–2. The antiquity of Neandertal gene flow into modern humans means that regions that derive from Neandertals in any one human today are usually less than a hundred kilobases in size. However, Neandertal haplotypes are also distinctive enough that several studies have been able to detect Neandertal ancestry at specific loci1,3–8. Here, we have systematically inferred Neandertal haplotypes in the genomes of 1,004 present-day humans12. Regions that harbor a high frequency of Neandertal alleles in modern humans are enriched for genes affecting keratin filaments suggesting that Neandertal alleles may have helped modern humans adapt to non-African environments. Neandertal alleles also continue to shape human biology, as we identify multiple Neandertal-derived alleles that confer risk for disease. We also identify regions of millions of base pairs that are nearly devoid of Neandertal ancestry and enriched in genes, implying selection to remove genetic material derived from Neandertals. Neandertal ancestry is significantly reduced in genes specifically expressed in testis, and there is an approximately 5-fold reduction of Neandertal ancestry on chromosome X, which is known to harbor a disproportionate fraction of male hybrid sterility genes20–22. These results suggest that part of the reduction in Neandertal ancestry near genes is due to Neandertal alleles that reduced fertility in males when moved to a modern human genetic background

Effect of trabecular bone loss on cortical strain rate during impact in an in vitro model of avian femur

BACKGROUND: Osteoporotic hip fractures occur due to loss of cortical and trabecular bone mass and consequent degradation in whole bone strength. The direct cause of most fractures is a fall, and hence, characterizing the mechanical behavior of a whole osteopenic bone under impact is important. However, very little is known about the mechanical interactions between cortical and trabecular bone during impact, and it is specifically unclear to what extent epiphyseal trabecular bone contributes to impact resistance of whole bones. We hypothesized that trabecular bone serves as a structural support to the cortex during impact, and hence, loss of a critical mass of trabecular bone reduces internal constraining of the cortex, and, thereby, decreases the impact tolerance of the whole bone. METHODS: To test this hypothesis, we conducted cortical strain rate measurements in adult chicken's proximal femora subjected to a Charpy impact test, after removing different trabecular bone core masses to simulate different osteopenic severities. RESULTS: We found that removal of core trabecular bone decreased by ~10-fold the cortical strain rate at the side opposite to impact (p < 0.01), i.e. from 359,815 ± 1799 μm/m per second (mean ± standard error) for an intact (control) specimen down to 35,997 ± 180 μm/m per second where 67% of the total trabecular bone mass (~0.7 grams in adult chicken) were removed. After normalizing the strain rate by the initial weight of bone specimens, a sigmoid relation emerged between normalized strain rate and removed mass of trabecular bone, showing very little effect on the cortex strain rate if below 10% of the trabecular mass is removed, but most of the effect was already apparent for less than 30% trabecular bone loss. An analytical model of the experiments supported this behavior. CONCLUSION: We conclude that in our in vitro avian model, loss of over 10% of core trabecular bone substantially altered the deformation response of whole bone to impact, which supports the above hypothesis and indicates that integrity of trabecular bone is critical for resisting impact loads

Endometrial stromal sarcoma: a population-based analysis

To determine independent prognostic factors for the survival of patients with endometrial stromal sarcoma (ESS), data were abstracted from the Surveillance, Epidemiology, and End Results (SEER) database of the National Cancer Institute from 1988 to 2003. Kaplan–Meier and Cox proportional hazards models were used for analyses. Of 831 women diagnosed with ESS, the median age was 52 years (range: 17–96 years). In total, 59.9% had stage I, 5.1% stage II, 14.9% stage III, and 20.1% had stage IV disease. Overall, 13.0, 36.1, and 34.7% presented with grades 1, 2, and 3, respectively. Patients with stage I–II vs III–IV disease had 5 years DSS of 89.3% vs 50.3% (P<0.001) and those with grades 1, 2, and 3 cancers had survivals of 91.4, 95.4, and 42.1% (P<0.001). In multivariate analysis, older patients, black race, advanced stage, higher grade, lack of primary surgery, and nodal metastasis were independent prognostic factors for poorer survival. In younger women (<50 years) with stage I–II disease, ovarian-sparing procedures did not adversely impact survival (91.9 vs 96.2%; P=0.1). Age, race, primary surgery, stage, and grade are important prognostic factors for ESS. Excellent survival in patients with grade 1 and 2 disease of all stages supports the concept that these tumors are significantly different from grade 3 tumors. Ovarian-sparing surgeries may be considered in younger patients with early-stage disease