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/

Maternal care in rural China: a case study from Anhui province

<p>Abstract</p> <p>Background</p> <p>Studies on prenatal care in China have focused on the timing and frequency of prenatal care and relatively little information can be found on how maternal care has been organized and funded or on the actual content of the visits, especially in the less developed rural areas. This study explored maternal care in a rural county from Anhui province in terms of care organization, provision and utilization.</p> <p>Methods</p> <p>A total of 699 mothers of infants under one year of age were interviewed with structured questionnaires; the county health bureau officials and managers of township hospitals (n = 10) and county level hospitals (n = 2) were interviewed; the process of the maternal care services was observed by the researchers. In addition, statistics from the local government were used.</p> <p>Results</p> <p>The county level hospitals were well staffed and equipped and served as a referral centre for women with a high-risk pregnancy. Township hospitals had, on average, 1.7 midwives serving an average population of 15,000 people. Only 10–20% of the current costs in county level hospitals and township hospitals were funded by the local government, and women paid for delivery care. There was no systematic organized prenatal care and referrals were not mandatory. About half of the women had their first prenatal visit before the 13th gestational week, 36% had fewer than 5 prenatal visits, and about 9% had no prenatal visits. A major reason for not having prenatal care visits was that women considered it unnecessary. Most women (87%) gave birth in public health facilities, and the rest in a private clinic or at home. A total of 8% of births were delivered by caesarean section. Very few women had any postnatal visits. About half of the women received the recommended number of prenatal blood pressure and haemoglobin measurements.</p> <p>Conclusion</p> <p>Delivery care was better provided than both prenatal and postnatal care in the study area. Reliance on user fees gave the hospitals an incentive to put more emphasis on revenue generating activities such as delivery care instead of prenatal and postnatal care.</p

Molecular Evolution of the Primate Antiviral Restriction Factor Tetherin

Background: Tetherin is a recently identified antiviral restriction factor that restricts HIV-1 particle release in the absence of the HIV-1 viral protein U (Vpu). It is reminiscent of APOBEC3G and TRIM5a that also antagonize HIV. APOBEC3G and TRIM5a have been demonstrated to evolve under pervasive positive selection throughout primate evolution, supporting the redqueen hypothesis. Therefore, one naturally presumes that Tetherin also evolves under pervasive positive selection throughout primate evolution and supports the red-queen hypothesis. Here, we performed a detailed evolutionary analysis to address this presumption. Methodology/Principal Findings: Results of non-synonymous and synonymous substitution rates reveal that Tetherin as a whole experiences neutral evolution rather than pervasive positive selection throughout primate evolution, as well as in non-primate mammal evolution. Sliding-window analyses show that the regions of the primate Tetherin that interact with viral proteins are under positive selection or relaxed purifying selection. In particular, the sites identified under positive selection generally focus on these regions, indicating that the main selective pressure acting on the primate Tetherin comes from virus infection. The branch-site model detected positive selection acting on the ancestral branch of the New World Monkey lineage, suggesting an episodic adaptive evolution. The positive selection was also found in duplicated Tetherins in ruminants. Moreover, there is no bias in the alterations of amino acids in the evolution of the primate Tetherin, implyin

Ancestral Components of Admixed Genomes in a Mexican Cohort

For most of the world, human genome structure at a population level is shaped by interplay between ancient geographic isolation and more recent demographic shifts, factors that are captured by the concepts of biogeographic ancestry and admixture, respectively. The ancestry of non-admixed individuals can often be traced to a specific population in a precise region, but current approaches for studying admixed individuals generally yield coarse information in which genome ancestry proportions are identified according to continent of origin. Here we introduce a new analytic strategy for this problem that allows fine-grained characterization of admixed individuals with respect to both geographic and genomic coordinates. Ancestry segments from different continents, identified with a probabilistic model, are used to construct and study “virtual genomes” of admixed individuals. We apply this approach to a cohort of 492 parent–offspring trios from Mexico City. The relative contributions from the three continental-level ancestral populations—Africa, Europe, and America—vary substantially between individuals, and the distribution of haplotype block length suggests an admixing time of 10–15 generations. The European and Indigenous American virtual genomes of each Mexican individual can be traced to precise regions within each continent, and they reveal a gradient of Amerindian ancestry between indigenous people of southwestern Mexico and Mayans of the Yucatan Peninsula. This contrasts sharply with the African roots of African Americans, which have been characterized by a uniform mixing of multiple West African populations. We also use the virtual European and Indigenous American genomes to search for the signatures of selection in the ancestral populations, and we identify previously known targets of selection in other populations, as well as new candidate loci. The ability to infer precise ancestral components of admixed genomes will facilitate studies of disease-related phenotypes and will allow new insight into the adaptive and demographic history of indigenous people

Detecting 22q11.2 deletion in Chinese children with conotruncal heart defects and single nucleotide polymorphisms in the haploid TBX1 locus

<p>Abstract</p> <p>Background</p> <p>Conotruncal heart defects (CTDs) are present in 75-85% of patients suffering from the 22q11.2 deletion syndrome. To date, no consistent phenotype has been consistently correlated with the 22q11.2 deletions. Genetic studies have implicated <it>TBX1 </it>as a critical gene in the pathogenesis of the syndrome. The aim of study was to determine the incidence of the 22q11.2 deletion in Chinese patients with CTDs and the possible mechanism for pathogenesis of CTDs.</p> <p>Methods</p> <p>We enrolled 212 patients with CTDs and 139 unrelated healthy controls. Both karyotypic analysis and multiplex ligation-dependent probe amplification were performed for all CTDs patients. Fluorescence <it>in situ </it>hybridization was performed for the patients with genetic deletions and their relatives. The <it>TBX1 </it>gene was sequenced for all patients and healthy controls. The <it>χ</it>²and Fisher's exact test were used in the statistical analysis.</p> <p>Results</p> <p>Thirteen of the 212 patients with CTDs (6.13%) were found to have the 22q11.2 deletion syndrome. Of the 13 cases, 11 presented with a hemizygous interstitial microdeletion from <it>CLTCL1 </it>to <it>LZTR1</it>; one presented with a regional deletion from <it>CLTCL1 </it>to <it>DRCR8</it>; and one presented with a regional deletion from <it>CDC45L </it>to <it>LZTR1</it>. There were eight sequence variants in the haploid <it>TBX1 </it>genes of the del22q11 CTDs patients. The frequency of one single nucleotide polymorphism (SNP) in the del22q11 patients was different from that of the non-del patients (<it>P </it>< 0.05), and the frequencies of two other SNPs were different between the non-del CTDs patients and controls (<it>P </it>< 0.05).</p> <p>Conclusions</p> <p>CTDs, especially pulmonary atresia with ventricular septal defect and tetralogy of Fallot, are the most common disorders associated with the 22q11.2 deletion syndrome. Those patients with both CTDs and 22q11.2 deletion generally have a typical or atypical deletion region within the <it>TBX1 </it>gene. Our results indicate that <it>TBX1 </it>genetic variants may be associated with CTDs.</p

QTLs for oil yield components in an elite oil palm (Elaeis guineensis) cross

Increased modern farming of superior types of the oil palm, Elaeis guineensis Jacq., which has naturally efficient oil biosynthesis, has made it the world’s foremost edible oil crop. Breeding improvement is, however, circumscribed by time and costs associated with the tree’s long reproductive cycle, large size and 10–15 years of field testing. Marker-assisted breeding has considerable potential for improving this crop. Towards this, quantitative trait loci (QTL) linked to oil yield component traits were mapped in a high-yield population. In total, 164 QTLs associated with 21 oil yield component traits were discovered, with cumulative QTL effects increasing in tandem with the number of QTL markers and matching the QT+ alleles for each trait. The QTLs confirmed all traits to be polygenic, with many genes of individual small effects on independent loci, but epistatic interactions are not ruled out. Furthermore, several QTLs maybe pleiotropic as suggested by QTL clustering of inter-related traits on almost all linkage groups. Certain regions of the chromosomes seem richer in the genes affecting a particular yield component trait and likely encompass pleiotropic, epistatic and heterotic effects. A large proportion of the identified additive effects from QTLs may actually arise from genic interactions between loci. Comparisons with previous mapping studies show that most of the QTLs were for similar traits and shared similar marker intervals on the same linkage groups. Practical applications for such QTLs in marker-assisted breeding will require seeking them out in different genetic backgrounds and environments

A Latent Variable Partial Least Squares Path Modeling Approach to Regional Association and Polygenic Effect with Applications to a Human Obesity Study

Genetic association studies are now routinely used to identify single nucleotide polymorphisms (SNPs) linked with human diseases or traits through single SNP-single trait tests. Here we introduced partial least squares path modeling (PLSPM) for association between single or multiple SNPs and a latent trait that can involve single or multiple correlated measurement(s). Furthermore, the framework naturally provides estimators of polygenic effect by appropriately weighting trait-attributing alleles. We conducted computer simulations to assess the performance via multiple SNPs and human obesity-related traits as measured by body mass index (BMI), waist and hip circumferences. Our results showed that the associate statistics had type I error rates close to nominal level and were powerful for a range of effect and sample sizes. When applied to 12 candidate regions in data (N = 2,417) from the European Prospective Investigation of Cancer (EPIC)-Norfolk study, a region in FTO was found to have stronger association (rs7204609∼rs9939881 at the first intron P = 4.29×10−7) than single SNP analysis (all with P>10−4) and a latent quantitative phenotype was obtained using a subset sample of EPIC-Norfolk (N = 12,559). We believe our method is appropriate for assessment of regional association and polygenic effect on a single or multiple traits

Inoculations with Arbuscular Mycorrhizal Fungi Increase Vegetable Yields and Decrease Phoxim Concentrations in Carrot and Green Onion and Their Soils

Background As one of the most widely used organophosphate insecticides in vegetable production, phoxim (C12H15N2O3PS) is often found as residues in crops and soils and thus poses a potential threat to public health and environment. Arbuscular mycorrhizal (AM) fungi may make a contribution to the decrease of organophosphate residues in crops and/or the degradation in soils, but such effects remain unknown. Methodology/Principal Findings A greenhouse pot experiment studied the influence of AM fungi and phoxim application on the growth of carrot and green onion, and phoxim concentrations in the two vegetables and their soil media. Treatments included three AM fungal inoculations with Glomus intraradices BEG 141, G. mosseae BEG 167, and a nonmycorrhizal control, and four phoxim application rates (0, 200, 400, 800 mg l−1, while 400 mg l−1 rate is the recommended dose in the vegetable production system). Carrot and green onion were grown in a greenhouse for 130 d and 150 d. Phoxim solution (100 ml) was poured into each pot around the roots 14d before plant harvest. Results showed that mycorrhizal colonization was higher than 70%, and phoxim application inhibited AM colonization on carrot but not on green onion. Compared with the nonmycorrhizal controls, both shoot and root fresh weights of these two vegetables were significantly increased by AM inoculations irrespective of phoxim application rates. Phoxim concentrations in shoots, roots and soils were increased with the increase of phoxim application rate, but significantly decreased by the AM inoculations. Soil phosphatase activity was enhanced by both AM inocula, but not affected by phoxim application rate. In general, G. intraradices BEG 141 had more pronounced effects than G. mosseae BEG 167 on the increase of fresh weight production in both carrot and green onion, and the decrease of phoxim concentrations in plants and soils. Conclusions/Significance Our results indicate a promising potential of AM fungi for enhancing vegetable production and reducing organophosphorus pesticide residues in plant tissues and their growth media, as well as for the phytoremediation of organophosphorus pesticide-contaminated soils

Shared midgut binding sites for Cry1A.105, Cry1Aa, Cry1Ab, Cry1Ac and Cry1Fa proteins from Bacillus thuringiensis in two important corn pests, Ostrinia nubilalis and Spodoptera frugiperda

First generation of insect-protected transgenic corn (Bt-corn) was based on the expression of Cry1Ab or Cry1Fa proteins. Currently, the trend is the combination of two or more genes expressing proteins that bind to different targets. In addition to broadening the spectrum of action, this strategy helps to delay the evolution of resistance in exposed insect populations. One of such examples is the combination of Cry1A.105 with Cry1Fa and Cry2Ab to control O. nubilalis and S. frugiperda. Cry1A.105 is a chimeric protein with domains I and II and the C-terminal half of the protein from Cry1Ac, and domain III almost identical to Cry1Fa. The aim of the present study was to determine whether the chimeric Cry1A.105 has shared binding sites either with Cry1A proteins, with Cry1Fa, or with both, in O. nubilalis and in S. frugiperda. Brush-border membrane vesicles (BBMV) from last instar larval midguts were used in competition binding assays with 125I-labeled Cry1A.105, Cry1Ab, and Cry1Fa, and unlabeled Cry1A.105, Cry1Aa, Cry1Ab, Cry1Ac, Cry1Fa, Cry2Ab and Cry2Ae. The results showed that Cry1A.105, Cry1Ab, Cry1Ac and Cry1Fa competed with high affinity for the same binding sites in both insect species. However, Cry2Ab and Cry2Ae did not compete for the binding sites of Cry1 proteins. Therefore, according to our results, the development of cross-resistance among Cry1Ab/Ac, Cry1A.105, and Cry1Fa proteins is possible in these two insect species if the alteration of shared binding sites occurs. Conversely, cross-resistance between these proteins and Cry2A proteins is very unlikely in such case