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/

Mortality from suicide and other external cause injuries in China: a prospective cohort study

<p>Abstract</p> <p>Background</p> <p>Premature death from suicide is a leading cause of death worldwide. However, the pattern and risk factors for suicide and other external cause injuries are not well understood. This study investigates mortality from suicide and other injuries and associated risk factors in China.</p> <p>Methods</p> <p>A prospective cohort study of 169,871 Chinese adults aged 40 years and older was conducted. Mortality due to suicide or other external cause injuries was recorded.</p> <p>Results</p> <p>Mortality from all external causes was 58.7/100,000 (72.3 in men and 44.4 in women): 14.1/100,000 (14.2 in men and 14.2 in women) for suicide and 44.6/100,000 (58.1 in men and 30.2 in women) for other external cause injuries. Transport accidents (17.2/100,000 overall, 23.4 in men and 10.8 in women), accidental poisoning (7.5/100,000 overall, 10.2 in men and 4.8 in women), and accidental falls (5.7/100,000 overall, 6.5 in men and 5.0 in women) were the three leading causes of death from other external cause injuries in China. In the multivariable analysis, male sex (relative risk [RR] 1.56, 95% confidence interval [CI] 1.03-2.38), age 70 years and older (2.27, 1.29-3.98), living in north China (1.68, 1.20-2.36) and rural residence (2.82, 1.76-4.51) were associated with increased mortality from suicide. Male sex (RR 2.50, 95% CI 1.95-3.20), age 60-69 years (1.93, 1.45-2.58) and 70 years and older (3.58, 2.58-4.97), rural residence (2.29, 1.77-2.96), and having no education (1.56, 1.00-2.43) were associated with increased mortality from other external cause injuries, while overweight (0.60, 0.43-0.83) was associated with decreased risk of mortality from other external cause injuries.</p> <p>Conclusions</p> <p>External cause mortality has become a major public health problem in China. Developing an integrated national program for the prevention of mortality due to external cause injuries in China is warranted.</p

Reactions of Cre with Methylphosphonate DNA: Similarities and Contrasts with Flp and Vaccinia Topoisomerase

Chien-Hui Ma is with UT Austin, Aashiq H. Kachroo is with UT Austin, Anna Macieszak is with Polish Academy of Sciences, Tzu-Yang Chen is with UT Austin, Piotr Guga is with Polish Academy of Sciences, Makkuni Jayaram is with UT Austin.Background -- Reactions of vaccinia topoisomerase and the tyrosine site-specific recombinase Flp with methylphosphonate (MeP) substituted DNA substrates, have provided important insights into the electrostatic features of the strand cleavage and strand joining steps catalyzed by them. A conserved arginine residue in the catalytic pentad, Arg-223 in topoisomerase and Arg-308 in Flp, is not essential for stabilizing the MeP transition state. Topoisomerase or its R223A variant promotes cleavage of the MeP bond by the active site nucleophile Tyr-274, followed by the rapid hydrolysis of the MeP-tyrosyl intermediate. Flp(R308A), but not wild type Flp, mediates direct hydrolysis of the activated MeP bond. These findings are consistent with a potential role for phosphate electrostatics and active site electrostatics in protecting DNA relaxation and site-specific recombination, respectively, against abortive hydrolysis. Methodology/Principal Findings -- We have examined the effects of DNA containing MeP substitution in the Flp related Cre recombination system. Neutralizing the negative charge at the scissile position does not render the tyrosyl intermediate formed by Cre susceptible to rapid hydrolysis. Furthermore, combining the active site R292A mutation in Cre (equivalent to the R223A and R308A mutations in topoisomerase and Flp, respectively) with MeP substitution does not lead to direct hydrolysis of the scissile MeP bond in DNA. Whereas Cre follows the topoisomerase paradigm during the strand cleavage step, it follows the Flp paradigm during the strand joining step. Conclusions/Significance -- Collectively, the Cre, Flp and topoisomerase results highlight the contribution of conserved electrostatic complementarity between substrate and active site towards transition state stabilization during site-specific recombination and DNA relaxation. They have potential implications for how transesterification reactions in nucleic acids are protected against undesirable abortive side reactions. Such protective mechanisms are significant, given the very real threat of hydrolytic genome damage or disruption of RNA processing due to the cellular abundance and nucleophilicity of water.This work was supported by the NIH award GM035654 to M. J. Partial support was provided by the Robert F. Welch Foundation (F-1274) and a Faculty Research Award from the University of Texas at Austin. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.Microbiolog

Evolution and Phylogenetic Analysis of Full-Length VP3 Genes of Eastern Mediterranean Bluetongue Virus Isolates

Bluetongue virus (BTV) is the ‘type’ species of the genus Orbivirus within the family Reoviridae. The BTV genome is composed of ten linear segments of double-stranded RNA (dsRNA), each of which codes for one of ten distinct viral proteins. Previous phylogenetic comparisons have evaluated variations in genome segment 3 (Seg-3) nucleotide sequence as way to identify the geographical origin (different topotypes) of BTV isolates. The full-length nucleotide sequence of genome Seg-3 was determined for thirty BTV isolates recovered in the eastern Mediterranean region, the Balkans and other geographic areas (Spain, India, Malaysia and Africa). These data were compared, based on molecular variability, positive-selection-analysis and maximum-likelihood phylogenetic reconstructions (using appropriate substitution models) to 24 previously published sequences, revealing their evolutionary relationships. These analyses indicate that negative selection is a major force in the evolution of BTV, restricting nucleotide variability, reducing the evolutionary rate of Seg-3 and potentially of other regions of the BTV genome. Phylogenetic analysis of the BTV-4 strains isolated over a relatively long time interval (1979–2000), in a single geographic area (Greece), showed a low level of nucleotide diversity, indicating that the virus can circulate almost unchanged for many years. These analyses also show that the recent incursions into south-eastern Europe were caused by BTV strains belonging to two different major-lineages: representing an ‘eastern’ (BTV-9, -16 and -1) and a ‘western’ (BTV-4) group/topotype. Epidemiological and phylogenetic analyses indicate that these viruses originated from a geographic area to the east and southeast of Greece (including Cyprus and the Middle East), which appears to represent an important ecological niche for the virus that is likely to represent a continuing source of future BTV incursions into Europe

The Relationship between Population Structure and Aluminum Tolerance in Cultivated Sorghum

Background: Acid soils comprise up to 50% of the world's arable lands and in these areas aluminum (Al) toxicity impairs root growth, strongly limiting crop yield. Food security is thereby compromised in many developing countries located in tropical and subtropical regions worldwide. In sorghum, SbMATE, an Al-activated citrate transporter, underlies the Alt(SB) locus on chromosome 3 and confers Al tolerance via Al-activated root citrate release. Methodology: Population structure was studied in 254 sorghum accessions representative of the diversity present in cultivated sorghums. Al tolerance was assessed as the degree of root growth inhibition in nutrient solution containing Al. A genetic analysis based on markers flanking Alt(SB) and SbMATE expression was undertaken to assess a possible role for Alt(SB) in Al tolerant accessions. In addition, the mode of gene action was estimated concerning the Al tolerance trait. Comparisons between models that include population structure were applied to assess the importance of each subpopulation to Al tolerance. Conclusion/Significance: Six subpopulations were revealed featuring specific racial and geographic origins. Al tolerance was found to be rather rare and present primarily in guinea and to lesser extent in caudatum subpopulations. Alt(SB) was found to play a role in Al tolerance in most of the Al tolerant accessions. A striking variation was observed in the mode of gene action for the Al tolerance trait, which ranged from almost complete recessivity to near complete dominance, with a higher frequency of partially recessive sources of Al tolerance. A possible interpretation of our results concerning the origin and evolution of Al tolerance in cultivated sorghum is discussed. This study demonstrates the importance of deeply exploring the crop diversity reservoir both for a comprehensive view of the dynamics underlying the distribution and function of Al tolerance genes and to design efficient molecular breeding strategies aimed at enhancing Al tolerance.CGIAR[G3007.04]McKnight FoundationFundacao de Amparo a Pesquisa do Estado de Minas Gerais (FAPEMIG)National Council for Scientific and Technological Development (CNPq

Crop Updates 2007 - Lupins, Pulses and Oilseeds

This session covers forty eight papers from different authors: 2006 REGIONAL ROUNDUP 1. South east agricultural region, Mark Seymour1 and Jacinta Falconer2, 1Department of Agriculture and Food, 2Cooperative Bulk Handling Group 2. Central agricultural region, Ian Pritchard, Department of Agriculture and Food 3. Great Southern and Lakes region, Rodger Beermier, Department of Agriculture and Food 4. Northern agricultural region, Wayne Parker and Martin Harries, Department of Agriculture and Food LUPINS 5. Development of anthracnose resistant and early flowering albus lupins (Lupinus albus L) in Western Australia, Kedar Adhikari and Geoff Thomas, Department of Agriculture and Food 6. New lupins adapted to the south coast, Peter White, Bevan Buirchell and Mike Baker, Department of Agriculture and Food 7. Lupin species and row spacing interactions by environment, Martin Harries, Peter White, Bob French, Jo Walker, Mike Baker and Laurie Maiolo, Department of Agriculture and Food 8. The interaction of lupin species row spacing and soil type, Martin Harries, Bob French, Laurie Maiolo and Jo Walker, Department of Agriculture and Food 9. The effects of row spacing and crop density on competitiveness of lupins with wild radish, Bob French and Laurie Maiolo, Department of Agriculture and Food 10. The effect of time of sowing and radish weed density on lupin yield, Martin Harries and Jo Walker, Department of Agriculture and Food 11. Interaction of time of sowing and weed management in lupins, Martin Harries and Jo Walker, Department of Agriculture and Food 12. Delayed sowing as a strategy to manage annual ryegrass, Bob French and Laurie Maiolo, Department of Agriculture and Food 13. Is delayed sowing a good strategy for weed management in lupins? Bob French, Department of Agriculture and Food 14. Lupins aren’t lupins when it comes to simazine, Peter White and Leigh Smith, Department of Agriculture and Food 15. Seed yield and anthracnose resistance of Tanjil mutants tolerant to metribuzin, Ping Si1, Bevan Buirchell1,2 and Mark Sweetingham1,2, 1Centre for Legumes in Mediterranean Agriculture, Australia; 2Department of Agriculture and Food 16. The effect of herbicides on nodulation in lupins, Lorne Mills1, Harmohinder Dhammu2 and Beng Tan1, 1Curtin University of Technology and 2Department of Agriculture and Food 17. Effect of fertiliser placements and watering regimes on lupin growth and seed yield in the central grain belt of Western Australia, Qifu Ma1, Zed Rengel1, Bill Bowden2, Ross Brennan2, Reg Lunt2 and Tim Hilder2, 1Soil Science & Plant Nutrition UWA, 2Department of Agriculture and Food 18. Development of a forecasting model for Bean Yellow Mosaic Virus in lupins, T. Maling1,2, A. Diggle1, D. Thackray1,2, R.A.C. Jones2, and K.H.M. Siddique1, 1Centre for Legumes in Mediterranean Agriculture, The University of Western Australia; 2Department of Agriculture and Food 19. Manufacturing of lupin tempe,Vijay Jayasena1,4, Leonardus Kardono2,4, Ken Quail3,4 and Ranil Coorey1,4, 1Curtin University of Technology, Perth, Australia, 2Indonesian Institute of Sciences (LIPI), Indonesia, 3BRI Australia Ltd, Sydney, Australia, 4Grain Foods CRC, Sydney, Australia 20. The impact of lupin based ingredients in ice-cream, Hannah Williams, Lee Sheer Yap and Vijay Jayasena, Curtin University of Technology, Perth WA 21. The acceptability of muffins substituted with varying concentrations of lupin flour, Anthony James, Don Elani Jayawardena and Vijay Jayasena, Curtin University of Technology, PerthWA PULSES 22. Chickpea variety evaluation, Kerry Regan1, Rod Hunter1, Tanveer Khan1,2and Jenny Garlinge1, 1Department of Agriculture and Food, 2CLIMA, The University of Western Australia 23. Advanced breeding trials of desi chickpea, Khan, T.N.1, Siddique, K.H.M.3, Clarke, H.2, Turner, N.C.2, MacLeod, W.1, Morgan, S.1, and Harris, A.1, 1Department of Agriculture and Food, 2Centre for Legumes in Mediterranean Agriculture, 3TheUniversity of Western Australia 24. Ascochyta resistance in chickpea lines in Crop Variety Testing (CVT) of 2006, Tanveer Khan1 2, Bill MacLeod1, Alan Harris1, Stuart Morgan1and Kerry Regan1, 1Department of Agriculture and Food, 2CLIMA, The University of Western Australia 25. Yield evaluation of ascochyta blight resistant Kabuli chickpeas, Kerry Regan1and Kadambot Siddique2, 1Department of Agriculture and Food, 2Institute of Agriculture, The University of Western Australia 26. Pulse WA Chickpea Industry Survey 2006, Mark Seymour1, Ian Pritchard1, Wayne Parker1and Alan Meldrum2, 1Department of Agriculture and Food, 2Pulse Australia 27. Genes from the wild as a valuable genetic resource for chickpea improvement, Heather Clarke1, Helen Bowers1and Kadambot Siddique2, 1Centre for Legumes in Mediterranean Agriculture, 2Institute of Agriculture, The University of Western Australia 28. International screening of chickpea for resistance to Botrytis grey mould, B. MacLeod1, Dr T. Khan1, Prof. K.H.M. Siddique2and Dr A. Bakr3, 1Department of Agriculture and Food, 2The University of Western Australia, 3Bangladesh Agricultural Research Institute 29. Balance® in chickpea is safest applied post sowing to a level seed bed, Wayne Parker, Department of Agriculture and Food, 30. Demonstrations of Genesis 510 chickpea, Wayne Parker, Department of Agriculture and Food 31. Field pea 2006, Ian Pritchard, Department of Agriculture and Food 32. Field pea variety evaluation, Kerry Regan1, Rod Hunter1, Tanveer Khan1,2 and Jenny Garlinge1, 1Department of Agriculture and Food, 2CLIMA, The University of Western Australia 33. Breeding highlights of the Australian Field Pea Improvement Program (AFPIP),Kerry Regan1, Tanveer Khan1,2, Phillip Chambers1, Chris Veitch1, Stuart Morgan1 , Alan Harris1and Tony Leonforte3, 1Department of Agriculture and Food, 2CLIMA, The University of Western Australia, 3Department of Primary Industries, Victoria 34. Field pea germplasm enhancement for black spot resistance, Tanveer Khan, Kerry Regan, Stuart Morgan, Alan Harris and Phillip Chambers, Department of Agriculture and Food 35. Validation of Blackspot spore release model and testing moderately resistant field pea line, Mark Seymour, Ian Pritchard, Rodger Beermier, Pam Burgess and Leanne Young, Department of Agriculture and Food 36. Yield losses from sowing field pea seed infected with Pea Seed-borne Mosaic Virus, Brenda Coutts, Donna O’Keefe, Rhonda Pearce, Monica Kehoe and Roger Jones, Department of Agriculture and Food 37. Faba bean in 2006, Mark Seymour, Department of Agriculture and Food 38. Germplasm evaluation – faba bean, Mark Seymour1, Terri Jasper1, Ian Pritchard1, Mike Baker1 and Tim Pope1,2, 1Department of Agriculture and Food, , 2CLIMA, The University of Western Australia 39. Breeding highlights of the Coordinated Improvement Program for Australian Lentils (CIPAL), Kerry Regan1, Chris Veitch1, Phillip Chambers1 and Michael Materne2, 1Department of Agriculture and Food, 2Department of Primary Industries, Victoria 40. Screening pulse lentil germplasm for tolerance to alternate herbicides, Ping Si1, Mike Walsh2 and Mark Sweetingham1,3, 1Centre for Legumes in Mediterranean Agriculture, 2West Australian Herbicide Resistance Initiative, 3Department of Agriculture and Food 41. Genomic synteny in legumes: Application to crop breeding, Phan, H.T.T.1, Ellwood, S.R.1, Hane, J.1, Williams, A.1, Ford, R.2, Thomas, S.3 and Oliver R1, 1Australian Centre of Necrotrophic Plant Pathogens, Murdoch University, 2BioMarka, University of Melbourne, 3NSW Department of Primary Industries 42. Tolerance of lupins, chickpeas and canola to Balanceâ(Isoxaflutole) and Galleryâ (Isoxaben), Leigh Smith and Peter White, Department of Agriculture and Food CANOLA AND OILSEEDS 43. The performance of TT Canola varieties in the National Variety Test (NVT),WA,2006,Katie Robinson, Research Agronomist, Agritech Crop Research 44. Evaluation of Brassica crops for biodiesel in Western Australia, Mohammad Amjad, Graham Walton, Pat Fels and Andy Sutherland, Department of Agriculture and Food 45. Production risk of canola in different rainfall zones in Western Australia, Imma Farré1, Michael Robertson2 and Senthold Asseng3, 1Department of Agriculture and Food, 2CSIRO Sustainable Ecosystems, 3CSIRO Plant Industry 46. Future directions of blackleg management – dynamics of blackleg susceptibility in canola varieties, Ravjit Khangura, Moin Salam and Bill MacLeod, Department of Agriculture and Food 47. Appendix 1: Contributors 48. Appendix 2: List of common acronym

Genome profiling of ERBB2-amplified breast cancers

<p>Abstract</p> <p>Background</p> <p>Around 20% of breast cancers (BC) show <it>ERBB2 </it>gene amplification and overexpression of the ERBB2 tyrosine kinase receptor. They are associated with a poor prognosis but can benefit from targeted therapy. A better knowledge of these BCs, genomically and biologically heterogeneous, may help understand their behavior and design new therapeutic strategies.</p> <p>Methods</p> <p>We defined the high resolution genome and gene expression profiles of 54 <it>ERBB2</it>-amplified BCs using 244K oligonucleotide array-comparative genomic hybridization and whole-genome DNA microarrays. Expression of ERBB2, phosphorylated ERBB2, EGFR, IGF1R and FOXA1 proteins was assessed by immunohistochemistry to evaluate the functional ERBB2 status and identify co-expressions.</p> <p>Results</p> <p>First, we identified the <it>ERBB2</it>-<it>C17orf37</it>-<it>GRB7 </it>genomic segment as the minimal common 17q12-q21 amplicon, and <it>CRKRS </it>and <it>IKZF3 </it>as the most frequent centromeric and telomeric amplicon borders, respectively. Second, GISTIC analysis identified 17 other genome regions affected by copy number aberration (CNA) (amplifications, gains, losses). The expression of 37 genes of these regions was deregulated. Third, two types of heterogeneity were observed in <it>ERBB2</it>-amplified BCs. The genomic profiles of estrogen receptor-postive (ER+) and negative (ER-) <it>ERBB2</it>-amplified BCs were different. The WNT/β-catenin signaling pathway was involved in ER- <it>ERBB2</it>-amplified BCs, and <it>PVT1 </it>and <it>TRPS1 </it>were candidate oncogenes associated with ER+ <it>ERBB2</it>-amplified BCs. The size of the <it>ERBB2 </it>amplicon was different in inflammatory (IBC) and non-inflammatory BCs. <it>ERBB2</it>-amplified IBCs were characterized by the downregulated and upregulated mRNA expression of ten and two genes in proportion to CNA, respectively. IHC results showed (i) a linear relationship between <it>ERBB2 </it>gene amplification and its gene and protein expressions with a good correlation between ERBB2 expression and phosphorylation status; (ii) a potential signaling cross-talk between EGFR or IGF1R and ERBB2, which could influence response of <it>ERBB2</it>-positive BCs to inhibitors. FOXA1 was frequently coexpressed with ERBB2 but its expression did not impact on the outcome of patients with <it>ERBB2</it>-amplified tumors.</p> <p>Conclusion</p> <p>We have shown that ER+ and ER- <it>ERBB2</it>-amplified BCs are different, distinguished <it>ERBB2 </it>amplicons in IBC and non-IBC, and identified genomic features that may be useful in the design of alternative therapeutical strategies.</p

The Evolutionary Dynamics of the Lion Panthera leo Revealed by Host and Viral Population Genomics

The lion Panthera leo is one of the world's most charismatic carnivores and is one of Africa's key predators. Here, we used a large dataset from 357 lions comprehending 1.13 megabases of sequence data and genotypes from 22 microsatellite loci to characterize its recent evolutionary history. Patterns of molecular genetic variation in multiple maternal (mtDNA), paternal (Y-chromosome), and biparental nuclear (nDNA) genetic markers were compared with patterns of sequence and subtype variation of the lion feline immunodeficiency virus (FIVPle), a lentivirus analogous to human immunodeficiency virus (HIV). In spite of the ability of lions to disperse long distances, patterns of lion genetic diversity suggest substantial population subdivision (mtDNA ΦST = 0.92; nDNA FST = 0.18), and reduced gene flow, which, along with large differences in sero-prevalence of six distinct FIVPle subtypes among lion populations, refute the hypothesis that African lions consist of a single panmictic population. Our results suggest that extant lion populations derive from several Pleistocene refugia in East and Southern Africa (∼324,000–169,000 years ago), which expanded during the Late Pleistocene (∼100,000 years ago) into Central and North Africa and into Asia. During the Pleistocene/Holocene transition (∼14,000–7,000 years), another expansion occurred from southern refugia northwards towards East Africa, causing population interbreeding. In particular, lion and FIVPle variation affirms that the large, well-studied lion population occupying the greater Serengeti Ecosystem is derived from three distinct populations that admixed recently

High Connectivity in the Deepwater Snapper Pristipomoides filamentosus (Lutjanidae) across the Indo-Pacific with Isolation of the Hawaiian Archipelago

In the tropical Indo-Pacific, most phylogeographic studies have focused on the shallow-water taxa that inhabit reefs to approximately 30 m depth. Little is known about the large predatory fishes, primarily snappers (subfamily Etelinae) and groupers (subfamily Epinephelinae) that occur at 100–400 m. These long-lived, slow-growing species support fisheries across the Indo-Pacific, yet no comprehensive genetic surveys within this group have been conducted. Here we contribute the first range-wide survey of a deepwater Indo-Pacific snapper, Pristipomoides filamentosus, with special focus on Hawai'i. We applied mtDNA cytochrome b and 11 microsatellite loci to 26 samples (N = 1,222) collected across 17,000 km from Hawai'i to the western Indian Ocean. Results indicate that P. filamentosus is a highly dispersive species with low but significant population structure (mtDNA ΦST = 0.029, microsatellite FST = 0.029) due entirely to the isolation of Hawai'i. No population structure was detected across 14,000 km of the Indo-Pacific from Tonga in the Central Pacific to the Seychelles in the western Indian Ocean, a pattern rarely observed in reef species. Despite a long pelagic phase (60–180 days), interisland dispersal as adults, and extensive gene flow across the Indo-Pacific, P. filamentosus is unable to maintain population connectivity with Hawai'i. Coalescent analyses indicate that P. filamentosus may have colonized Hawai'i 26 K–52 K y ago against prevailing currents, with dispersal away from Hawai'i dominating migration estimates. P. filamentosus harbors low genetic diversity in Hawai'i, a common pattern in marine fishes, and our data indicate a single archipelago-wide stock. However, like the Hawaiian Grouper, Hyporthodus quernus, this snapper had several significant pairwise comparisons (FST) clustered around the middle of the archipelago (St. Rogatien, Brooks Banks, Gardner) indicating that this region may be isolated or (more likely) receives input from Johnston Atoll to the south