Search

How and when the Americas were populated remains contentious. Using ancient and modern genome-wide data, we found that the ancestors of all present-day Native Americans, including Athabascans and Amerindians, entered the Americas as a single migration wave from Siberia no earlier than 23 thousand years ago (ka) and after no more than an 8000-year isolation period in Beringia. After their arrival to the Americas, ancestral Native Americans diversified into two basal genetic branches around 13 ka, one that is now dispersed across North and South America and the other restricted to North America. Subsequent gene flow resulted in some Native Americans sharing ancestry with present-day East Asians (including Siberians) and, more distantly, Australo-Melanesians. Putative "Paleoamerican" relict populations, including the historical Mexican Pericúes and South American Fuego-Patagonians, are not directly related to modern Australo-Melanesians as suggested by the Paleoamerican Model.

Mountain gorillas are an endangered great ape subspecies and a prominent focus for conservation, yet we know little about their genomic diversity and evolutionary past. We sequenced whole genomes from multiple wild individuals and compared the genomes of all four Gorilla subspecies. We found that the two eastern subspecies have experienced a prolonged population decline over the past 100,000 years, resulting in very low genetic diversity and an increased overall burden of deleterious variation. A further recent decline in the mountain gorilla population has led to extensive inbreeding, such that individuals are typically homozygous at 34% of their sequence, leading to the purging of severely deleterious recessive mutations from the population. We discuss the causes of their decline and the consequences for their future survival.

It is commonly thought that human genetic diversity in non-African populations was shaped primarily by an out-of-Africa dispersal 50-100 thousand yr ago (kya). Here, we present a study of 456 geographically diverse high-coverage Y chromosome sequences, including 299 newly reported samples. Applying ancient DNA calibration, we date the Y-chromosomal most recent common ancestor (MRCA) in Africa at 254 (95% CI 192-307) kya and detect a cluster of major non-African founder haplogroups in a narrow time interval at 47-52 kya, consistent with a rapid initial colonization model of Eurasia and Oceania after the out-of-Africa bottleneck. In contrast to demographic reconstructions based on mtDNA, we infer a second strong bottleneck in Y-chromosome lineages dating to the last 10 ky. We hypothesize that this bottleneck is caused by cultural changes affecting variance of reproductive success among males.

Arctic populations live in an environment characterized by extreme cold and the absence of plant foods for much of the year and are likely to have undergone genetic adaptations to these environmental conditions in the time they have been living there. Genome-wide selection scans based on genotype data from native Siberians have previously highlighted a 3 Mb chromosome 11 region containing 79 protein-coding genes as the strongest candidates for positive selection in Northeast Siberians. However, it was not possible to determine which of the genes might be driving the selection signal. Here, using whole-genome high-coverage sequence data, we identified the most likely causative variant as a nonsynonymous G>A transition (rs80356779; c.1436C>T [p.Pro479Leu] on the reverse strand) in CPT1A, a key regulator of mitochondrial long-chain fatty-acid oxidation. Remarkably, the derived allele is associated with hypoketotic hypoglycemia and high infant mortality yet occurs at high frequency in Canadian and Greenland Inuits and was also found at 68% frequency in our Northeast Siberian sample. We provide evidence of one of the strongest selective sweeps reported in humans; this sweep has driven this variant to high frequency in circum-Arctic populations within the last 6-23 ka despite associated deleterious consequences, possibly as a result of the selective advantage it originally provided to either a high-fat diet or a cold environment.

The Wellcome Trust Case Control Consortium 3 anorexia nervosa genome-wide association scan includes 2907 cases from 15 different populations of European origin genotyped on the Illumina 670K chip. We compared methods for identifying population stratification, and suggest list of markers that may help to counter this problem. It is usual to identify population structure in such studies using only common variants with minor allele frequency (MAF) >5%; we find that this may result in highly informative SNPs being discarded, and suggest that instead all SNPs with MAF >1% may be used. We established informative axes of variation identified via principal component analysis and highlight important features of the genetic structure of diverse European-descent populations, some studied for the first time at this scale. Finally, we investigated the substructure within each of these 15 populations and identified SNPs that help capture hidden stratification. This work can provide information regarding the designing and interpretation of association results in the International Consortia.

Abstract Background Large palindromes (inverted repeats) make up substantial proportions of mammalian sex chromosomes, often contain genes, and have high rates of structural variation arising via ectopic recombination. As a result, they underlie many genomic disorders. Maintenance of the palindromic structure by gene conversion between the arms has been documented, but over longer time periods, palindromes are remarkably labile. Mechanisms of origin and loss of palindromes have, however, received little attention. Results Here, we use fiber-FISH, 10x Genomics Linked-Read sequencing, and breakpoint PCR sequencing to characterize the structural variation of the P8 palindrome on the human Y chromosome, which contains two copies of the VCY (Variable Charge Y) gene. We find a deletion of almost an entire arm of the palindrome, leading to death of the palindrome, a size increase by recruitment of adjacent sequence, and other complex changes including the formation of an entire new palindrome nearby. Together, these changes are found in ~ 1% of men, and we can assign likely molecular mechanisms to these mutational events. As a result, healthy men can have 1–4 copies of VCY. Conclusions Gross changes, especially duplications, in palindrome structure can be relatively frequent and facilitate the evolution of sex chromosomes in humans, and potentially also in other mammalian species.

Male infertility is a prevalent condition, affecting 5–10% of men. So far, few genetic factors have been described as contributors to spermatogenic failure. Here, we report the first re-sequencing study of the Y-chromosomal Azoospermia Factor c (AZFc) region, combined with gene dosage analysis of the multicopy DAZ, BPY2, and CDYgenes and Y-haplogroup determination. In analysing 2324 Estonian men, we uncovered a novel structural variant as a high-penetrance risk factor for male infertility. The Y lineage R1a1-M458, reported at >20% frequency in several European populations, carries a fixed ~1.6 Mb r2/r3 inversion, destabilizing the AZFc region and predisposing to large recurrent microdeletions. Such complex rearrangements were significantly enriched among severe oligozoospermia cases. The carrier vs non-carrier risk for spermatogenic failure was increased 8.6-fold (p=6.0×10−4). This finding contributes to improved molecular diagnostics and clinical management of infertility. Carrier identification at young age will facilitate timely counselling and reproductive decision-making.

The Tibetan and Andean Plateaus and Ethiopian highlands are the largest regions to have long-term high-altitude residents. Such populations are exposed to lower barometric pressures and hence atmospheric partial pressures of oxygen. Such "hypobaric hypoxia" may limit physical functional capacity, reproductive health, and even survival. As such, selection of genetic variants advantageous to hypoxic adaptation is likely to have occurred. Identifying signatures of such selection is likely to help understanding of hypoxic adaptive processes. Here, we seek evidence of such positive selection using five Ethiopian populations, three of which are from high-altitude areas in Ethiopia. As these populations may have been recipients of Eurasian gene flow, we correct for this admixture. Using single-nucleotide polymorphism genotype data from multiple populations, we find the strongest signal of selection in BHLHE41 (also known as DEC2 or SHARP1). Remarkably, a major role of this gene is regulation of the same hypoxia response pathway on which selection has most strikingly been observed in both Tibetan and Andean populations. Because it is also an important player in the circadian rhythm pathway, BHLHE41 might also provide insights into the mechanisms underlying the recognized impacts of hypoxia on the circadian clock. These results support the view that Ethiopian, Andean, and Tibetan populations living at high altitude have adapted to hypoxia differently, with convergent evolution affecting different genes from the same pathway.

The Levant is a region in the Near East with an impressive record of continuous human existence and major cultural developments since the Paleolithic period. Genetic and archeological studies present solid evidence placing the Middle East and the Arabian Peninsula as the first stepping-stone outside Africa. There is, however, little understanding of demographic changes in the Middle East, particularly the Levant, after the first Out-of-Africa expansion and how the Levantine peoples relate genetically to each other and to their neighbors. In this study we analyze more than 500,000 genome-wide SNPs in 1,341 new samples from the Levant and compare them to samples from 48 populations worldwide. Our results show recent genetic stratifications in the Levant are driven by the religious affiliations of the populations within the region. Cultural changes within the last two millennia appear to have facilitated/maintained admixture between culturally similar populations from the Levant, Arabian Peninsula, and Africa. The same cultural changes seem to have resulted in genetic isolation of other groups by limiting admixture with culturally different neighboring populations. Consequently, Levant populations today fall into two main groups: one sharing more genetic characteristics with modern-day Europeans and Central Asians, and the other with closer genetic affinities to other Middle Easterners and Africans. Finally, we identify a putative Levantine ancestral component that diverged from other Middle Easterners ∼23,700-15,500 years ago during the last glacial period, and diverged from Europeans ∼15,900-9,100 years ago between the last glacial warming and the start of the Neolithic.

IntroductionThe COVID-19 pandemic has caused widespread disruptions including to health services. In the early response to the pandemic many countries restricted population movements and some health services were suspended or limited. In late 2020 and early 2021 some countries re-imposed restrictions. Health authorities need to balance the potential harms of additional SARS-CoV-2 transmission due to contacts associated with health services against the benefits of those services, including fewer new HIV infections and deaths. This paper examines these trade-offs for select HIV services.MethodsWe used four HIV simulation models (Goals, HIV Synthesis, Optima HIV and EMOD) to estimate the benefits of continuing HIV services in terms of fewer new HIV infections and deaths. We used three COVID-19 transmission models (Covasim, Cooper/Smith and a simple contact model) to estimate the additional deaths due to SARS-CoV-2 transmission among health workers and clients. We examined four HIV services: voluntary medical male circumcision, HIV diagnostic testing, viral load testing and programs to prevent mother-to-child transmission. We compared COVID-19 deaths in 2020 and 2021 with HIV deaths occurring now and over the next 50 years discounted to present value. The models were applied to countries with a range of HIV and COVID-19 epidemics.ResultsMaintaining these HIV services could lead to additional COVID-19 deaths of 0.002 to 0.15 per 10,000 clients. HIV-related deaths averted are estimated to be much larger, 19-146 discounted deaths per 10,000 clients.DiscussionWhile there is some additional short-term risk of SARS-CoV-2 transmission associated with providing HIV services, the risk of additional COVID-19 deaths is at least 100 times less than the HIV deaths averted by those services. Ministries of Health need to take into account many factors in deciding when and how to offer essential health services during the COVID-19 pandemic. This work shows that the benefits of continuing key HIV services are far larger than the risks of additional SARS-CoV-2 transmission.

The COVID-19 global pandemic has had considerable health impact, including sub-Saharan Africa. In Malawi, a resource-limited setting in Africa, gaining access to data to inform the COVID-19 response is challenging. Information on adherence to physical distancing guidelines and reducing contacts are nonexistent, but critical to understanding and communicating risk, as well as allocating scarce resources. We present a case study which leverages aggregated call detail records into a daily data pipeline which summarize population density and mobility in an easy-to-use dashboard for public health officials and emergency operations. From March to April 2021, we have aggregated 6-billion calls and text messages and continue to process 12 million more daily. These data are summarized into reports which describe, quantify, and locate mass gatherings and travel between subdistricts. These reports are accessible via web dashboards for policymakers within the Ministry of Health and Emergency Operations Center to inform COVID-19 response efforts and resource allocation.

BackgroundExome sequencing, which allows the global analysis of protein coding sequences in the human genome, has become an effective and affordable approach to detecting causative genetic mutations in diseases. Currently, there are several commercial human exome capture platforms; however, the relative performances of these have not been characterized sufficiently to know which is best for a particular study.ResultsWe comprehensively compared three platforms: NimbleGen's Sequence Capture Array and SeqCap EZ, and Agilent's SureSelect. We assessed their performance in a variety of ways, including number of genes covered and capture efficacy. Differences that may impact on the choice of platform were that Agilent SureSelect covered approximately 1,100 more genes, while NimbleGen provided better flanking sequence capture. Although all three platforms achieved similar capture specificity of targeted regions, the NimbleGen platforms showed better uniformity of coverage and greater genotype sensitivity at 30- to 100-fold sequencing depth. All three platforms showed similar power in exome SNP calling, including medically relevant SNPs. Compared with genotyping and whole-genome sequencing data, the three platforms achieved a similar accuracy of genotype assignment and SNP detection. Importantly, all three platforms showed similar levels of reproducibility, GC bias and reference allele bias.ConclusionsWe demonstrate key differences between the three platforms, particularly advantages of solutions over array capture and the importance of a large gene target set.

Plants adapt to heat via thermotolerance pathways in which the activation of protein folding chaperones is essential. In eukaryotes, protein disulfide isomerases (PDIs) facilitate the folding of nascent and misfolded proteins in the secretory pathway by catalyzing the formation and isomerization of disulfide bonds and serving as molecular chaperones. In Arabidopsis, several members of the PDI family are upregulated in response to chemical inducers of the unfolded protein response (UPR), including both members of the non-classical PDI-M subfamily, PDI9 and PDI10. Unlike classical PDIs, which have two catalytic thioredoxin (TRX) domains separated by two non-catalytic TRX-fold domains, PDI-M isoforms are orthologs of mammalian P5/PDIA6 and possess two tandem catalytic domains. Here, PDI9 accumulation was found to be upregulated in pollen in response to heat stress. Histochemical staining of plants harboring the PDI9 and PDI10 promoters fused to the gusA gene indicated they were actively expressed in the anthers of flowers, specifically in the pollen and tapetum. Immunoelectron microscopy revealed that PDI9 localized to the endoplasmic reticulum in root and pollen cells. transfer DNA (T-DNA) insertional mutations in the PDI9 gene disrupted pollen viability and development in plants exposed to heat stress. In particular, the pollen grains of pdi9 mutants exhibited disruptions in the reticulated pattern of the exine and an increased adhesion of pollen grains. Pollen in the pdi10 single mutant did not display similar heat-associated defects, but pdi9 pdi10 double mutants (DMs) completely lost exine reticulation. Interestingly, overexpression of PDI9 partially led to heat-associated defects in the exine. We conclude that PDI9 plays an important role in pollen thermotolerance and exine biogenesis. Its role fits the mechanistic theory of proteostasis in which an ideal balance of PDI isoforms is required in the endoplasmic reticulum (ER) for normal exine formation in plants subjected to heat stress.