Population Structure, Stratification, and Introgression of Human Structural Variation.

Structural variants contribute substantially to genetic diversity and are important evolutionarily and medically, but they are still understudied. Here we present a comprehensive analysis of structural variation in the Human Genome Diversity panel, a high-coverage dataset of 911 samples from 54 diverse worldwide populations. We identify, in total, 126,018 variants, 78% of which were not identified in previous global sequencing projects. Some reach high frequency and are private to continental groups or even individual populations, including regionally restricted runaway duplications and putatively introgressed variants from archaic hominins. By de novo assembly of 25 genomes using linked-read sequencing, we discover 1,643 breakpoint-resolved unique insertions, in aggregate accounting for 1.9 Mb of sequence absent from the GRCh38 reference. Our results illustrate the limitation of a single human reference and the need for high-quality genomes from diverse populations to fully discover and understand human genetic variation.

Uganda Genome Resource Enables Insights into Population History and Genomic Discovery in Africa.

Genomic studies in African populations provide unique opportunities to understand disease etiology, human diversity, and population history. In the largest study of its kind, comprising genome-wide data from 6,400 individuals and whole-genome sequences from 1,978 individuals from rural Uganda, we find evidence of geographically correlated fine-scale population substructure. Historically, the ancestry of modern Ugandans was best represented by a mixture of ancient East African pastoralists. We demonstrate the value of the largest sequence panel from Africa to date as an imputation resource. Examining 34 cardiometabolic traits, we show systematic differences in trait heritability between European and African populations, probably reflecting the differential impact of genes and environment. In a multi-trait pan-African GWAS of up to 14,126 individuals, we identify novel loci associated with anthropometric, hematological, lipid, and glycemic traits. We find that several functionally important signals are driven by Africa-specific variants, highlighting the value of studying diverse populations across the region.

Grey wolf genomic history reveals a dual ancestry of dogs.

The grey wolf (Canis lupus) was the first species to give rise to a domestic population, and they remained widespread throughout the last Ice Age when many other large mammal species went extinct. Little is known, however, about the history and possible extinction of past wolf populations or when and where the wolf progenitors of the present-day dog lineage (Canis familiaris) lived1-8. Here we analysed 72 ancient wolf genomes spanning the last 100,000 years from Europe, Siberia and North America. We found that wolf populations were highly connected throughout the Late Pleistocene, with levels of differentiation an order of magnitude lower than they are today. This population connectivity allowed us to detect natural selection across the time series, including rapid fixation of mutations in the gene IFT88 40,000-30,000 years ago. We show that dogs are overall more closely related to ancient wolves from eastern Eurasia than to those from western Eurasia, suggesting a domestication process in the east. However, we also found that dogs in the Near East and Africa derive up to half of their ancestry from a distinct population related to modern southwest Eurasian wolves, reflecting either an independent domestication process or admixture from local wolves. None of the analysed ancient wolf genomes is a direct match for either of these dog ancestries, meaning that the exact progenitor populations remain to be located.

Origins of modern human ancestry.

New finds in the palaeoanthropological and genomic records have changed our view of the origins of modern human ancestry. Here we review our current understanding of how the ancestry of modern humans around the globe can be traced into the deep past, and which ancestors it passes through during our journey back in time. We identify three key phases that are surrounded by major questions, and which will be at the frontiers of future research. The most recent phase comprises the worldwide expansion of modern humans between 40 and 60 thousand years ago (ka) and their last known contacts with archaic groups such as Neanderthals and Denisovans. The second phase is associated with a broadly construed African origin of modern human diversity between 60 and 300 ka. The oldest phase comprises the complex separation of modern human ancestors from archaic human groups from 0.3 to 1 million years ago. We argue that no specific point in time can currently be identified at which modern human ancestry was confined to a limited birthplace, and that patterns of the first appearance of anatomical or behavioural traits that are used to define Homo sapiens are consistent with a range of evolutionary histories.

Million-year-old DNA sheds light on the genomic history of mammoths.

Temporal genomic data hold great potential for studying evolutionary processes such as speciation. However, sampling across speciation events would, in many cases, require genomic time series that stretch well back into the Early Pleistocene subepoch. Although theoretical models suggest that DNA should survive on this timescale1, the oldest genomic data recovered so far are from a horse specimen dated to 780-560 thousand years ago2. Here we report the recovery of genome-wide data from three mammoth specimens dating to the Early and Middle Pleistocene subepochs, two of which are more than one million years old. We find that two distinct mammoth lineages were present in eastern Siberia during the Early Pleistocene. One of these lineages gave rise to the woolly mammoth and the other represents a previously unrecognized lineage that was ancestral to the first mammoths to colonize North America. Our analyses reveal that the Columbian mammoth of North America traces its ancestry to a Middle Pleistocene hybridization between these two lineages, with roughly equal admixture proportions. Finally, we show that the majority of protein-coding changes associated with cold adaptation in woolly mammoths were already present one million years ago. These findings highlight the potential of deep-time palaeogenomics to expand our understanding of speciation and long-term adaptive evolution.

Genome evolution across 1,011 Saccharomyces cerevisiae isolates.

Large-scale population genomic surveys are essential to explore the phenotypic diversity of natural populations. Here we report the whole-genome sequencing and phenotyping of 1,011 Saccharomyces cerevisiae isolates, which together provide an accurate evolutionary picture of the genomic variants that shape the species-wide phenotypic landscape of this yeast. Genomic analyses support a single 'out-of-China' origin for this species, followed by several independent domestication events. Although domesticated isolates exhibit high variation in ploidy, aneuploidy and genome content, genome evolution in wild isolates is mainly driven by the accumulation of single nucleotide polymorphisms. A common feature is the extensive loss of heterozygosity, which represents an essential source of inter-individual variation in this mainly asexual species. Most of the single nucleotide polymorphisms, including experimentally identified functional polymorphisms, are present at very low frequencies. The largest numbers of variants identified by genome-wide association are copy-number changes, which have a greater phenotypic effect than do single nucleotide polymorphisms. This resource will guide future population genomics and genotype-phenotype studies in this classic model system.

Archaeological Central American maize genomes suggest ancient gene flow from South America.

Maize (Zea mays ssp. mays) domestication began in southwestern Mexico ∼9,000 calendar years before present (cal. BP) and humans dispersed this important grain to South America by at least 7,000 cal. BP as a partial domesticate. South America served as a secondary improvement center where the domestication syndrome became fixed and new lineages emerged in parallel with similar processes in Mesoamerica. Later, Indigenous cultivators carried a second major wave of maize southward from Mesoamerica, but it has been unclear until now whether the deeply divergent maize lineages underwent any subsequent gene flow between these regions. Here we report ancient maize genomes (2,300-1,900 cal. BP) from El Gigante rock shelter, Honduras, that are closely related to ancient and modern maize from South America. Our findings suggest that the second wave of maize brought into South America hybridized with long-established landraces from the first wave, and that some of the resulting newly admixed lineages were then reintroduced to Central America. Direct radiocarbon dates and cob morphological data from the rock shelter suggest that more productive maize varieties developed between 4,300 and 2,500 cal. BP. We hypothesize that the influx of maize from South America into Central America may have been an important source of genetic diversity as maize was becoming a staple grain in Central and Mesoamerica.

Incidental cardiac findings on somatostatin receptor PET/CT: What do they indicate and are they of clinical relevance?

We present the case of a 47-year-old man with a history of recurrent episodes of frontal headache, fever, and chest discomfort as well as longstanding, difficult to treat arterial hypertension. Clinical work-up revealed the unexpected finding of an underlying pheochromocytoma as well as recent "silent" myocardial infarction. Our case highlights the importance of paying attention to incidental cardiac findings on somatostatin receptor positron emission tomography/computed tomography, as routinely performed in patients with clinically suspected neuroendocrine tumors. These incidental cardiac findings cannot only indicate a primary or secondary (metastatic) neuroendocrine tumor, but also areas of myocardial inflammation, as somatostatin receptors cannot only be found on the majority of neuroendocrine tumors, but also among other tissues on the surface of activated macrophages and lymphocytes. The detection of myocardial inflammation is of clinical importance and its underlying etiology should be evaluated to prompt eventual necessary treatment, as it is a potential driving force for cardiac remodeling and poor prognosis.

A genomic history of Aboriginal Australia.

The population history of Aboriginal Australians remains largely uncharacterized. Here we generate high-coverage genomes for 83 Aboriginal Australians (speakers of Pama-Nyungan languages) and 25 Papuans from the New Guinea Highlands. We find that Papuan and Aboriginal Australian ancestors diversified 25-40 thousand years ago (kya), suggesting pre-Holocene population structure in the ancient continent of Sahul (Australia, New Guinea and Tasmania). However, all of the studied Aboriginal Australians descend from a single founding population that differentiated ~10-32 kya. We infer a population expansion in northeast Australia during the Holocene epoch (past 10,000 years) associated with limited gene flow from this region to the rest of Australia, consistent with the spread of the Pama-Nyungan languages. We estimate that Aboriginal Australians and Papuans diverged from Eurasians 51-72 kya, following a single out-of-Africa dispersal, and subsequently admixed with archaic populations. Finally, we report evidence of selection in Aboriginal Australians potentially associated with living in the desert.

Chad Genetic Diversity Reveals an African History Marked by Multiple Holocene Eurasian Migrations.

Understanding human genetic diversity in Africa is important for interpreting the evolution of all humans, yet vast regions in Africa, such as Chad, remain genetically poorly investigated. Here, we use genotype data from 480 samples from Chad, the Near East, and southern Europe, as well as whole-genome sequencing from 19 of them, to show that many populations today derive their genomes from ancient African-Eurasian admixtures. We found evidence of early Eurasian backflow to Africa in people speaking the unclassified isolate Laal language in southern Chad and estimate from linkage-disequilibrium decay that this occurred 4,750-7,200 years ago. It brought to Africa a Y chromosome lineage (R1b-V88) whose closest relatives are widespread in present-day Eurasia; we estimate from sequence data that the Chad R1b-V88 Y chromosomes coalesced 5,700-7,300 years ago. This migration could thus have originated among Near Eastern farmers during the African Humid Period. We also found that the previously documented Eurasian backflow into Africa, which occurred ∼3,000 years ago and was thought to be mostly limited to East Africa, had a more westward impact affecting populations in northern Chad, such as the Toubou, who have 20%-30% Eurasian ancestry today. We observed a decline in heterozygosity in admixed Africans and found that the Eurasian admixture can bias inferences on their coalescent history and confound genetic signals from adaptation and archaic introgression.

Extensive Recombination of a Yeast Diploid Hybrid through Meiotic Reversion.

In somatic cells, recombination between the homologous chromosomes followed by equational segregation leads to loss of heterozygosity events (LOH), allowing the expression of recessive alleles and the production of novel allele combinations that are potentially beneficial upon Darwinian selection. However, inter-homolog recombination in somatic cells is rare, thus reducing potential genetic variation. Here, we explored the property of S. cerevisiae to enter the meiotic developmental program, induce meiotic Spo11-dependent double-strand breaks genome-wide and return to mitotic growth, a process known as Return To Growth (RTG). Whole genome sequencing of 36 RTG strains derived from the hybrid S288c/SK1 diploid strain demonstrates that the RTGs are bona fide diploids with mosaic recombined genome, derived from either parental origin. Individual RTG genome-wide genotypes are comprised of 5 to 87 homozygous regions due to the loss of heterozygous (LOH) events of various lengths, varying between a few nucleotides up to several hundred kilobases. Furthermore, we show that reiteration of the RTG process shows incremental increases of homozygosity. Phenotype/genotype analysis of the RTG strains for the auxotrophic and arsenate resistance traits validates the potential of this procedure of genome diversification to rapidly map complex traits loci (QTLs) in diploid strains without undergoing sexual reproduction.

Copy number variation arising from gene conversion on the human Y chromosome.

We describe the variation in copy number of a ~ 10 kb region overlapping the long intergenic noncoding RNA (lincRNA) gene, TTTY22, within the IR3 inverted repeat on the short arm of the human Y chromosome, leading to individuals with 0-3 copies of this region in the general population. Variation of this CNV is common, with 266 individuals having 0 copies, 943 (including the reference sequence) having 1, 23 having 2 copies, and two having 3 copies, and was validated by breakpoint PCR, fibre-FISH, and 10× Genomics Chromium linked-read sequencing in subsets of 1234 individuals from the 1000 Genomes Project. Mapping the changes in copy number to the phylogeny of these Y chromosomes previously established by the Project identified at least 20 mutational events, and investigation of flanking paralogous sequence variants showed that the mutations involved flanking sequences in 18 of these, and could extend over > 30 kb of DNA. While either gene conversion or double crossover between misaligned sister chromatids could formally explain the 0-2 copy events, gene conversion is the more likely mechanism, and these events include the longest non-allelic gene conversion reported thus far. Chromosomes with three copies of this CNV have arisen just once in our data set via another mechanism: duplication of 420 kb that places the third copy 230 kb proximal to the existing proximal copy. Our results establish gene conversion as a previously under-appreciated mechanism of generating copy number changes in humans and reveal the exceptionally large size of the conversion events that can occur.

miRNA profiling of small intestinal neuroendocrine tumors defines novel molecular subtypes and identifies miR-375 as a biomarker of patient survival.

The aim of this study was to define the miRNA profile of small intestinal neuroendocrine tumors and to search for novel molecular subgroups and prognostic biomarkers. miRNA profiling was conducted on 42 tumors from 37 patients who underwent surgery for small intestinal neuroendocrine tumors. Unsupervised hierarchical clustering analysis of miRNA profiles identified two groups of tumor metastases, denoted cluster M1 and M2. The smaller cluster M1 was associated with shorter overall survival and contained tumors with higher grade (WHO grade G2/3) and multiple chromosomal gains including gain of chromosome 14. Tumors of cluster M1 had elevated expression of miR-1246 and miR-663a, and reduced levels of miR-488-3p. Pathway analysis predicted Wnt signaling to be the most significantly altered signaling pathway between clusters M1 and M2. Analysis of miRNA expression in relation to tumor proliferation rate showed significant alterations including downregulation of miR-137 and miR-204-5p in tumors with Ki67 index above 3%. Similarly, tumor progression was associated with significant alterations in miRNA expression, e.g. higher expression of miR-95 and miR-210, and lower expression of miR-378a-3p in metastases. Pathway analysis predicted Wnt signaling to be altered during tumor progression, which was supported by decreased nuclear translocation of ß-catenin in metastases. Survival analysis revealed that downregulation of miR-375 was associated with shorter overall survival. We performed in situ hybridization on biopsies from an independent cohort of small intestinal neuroendocrine tumors using tissue microarrays. Expression of miR-375 was found in 578/635 (91%) biopsies and survival analysis confirmed that there was a correlation between downregulation of miR-375 in tumor metastases and shorter patient survival. We conclude that miRNA profiling defines novel molecular subgroups of metastatic small intestinal neuroendocrine tumors and identifies miRNAs associated with tumor proliferation rate and progression. miR-375 is highly expressed in small intestinal neuroendocrine tumors and may be used as a prognostic biomarker.

Y-chromosomal sequences of diverse Indian populations and the ancestry of the Andamanese.

We present 42 new Y-chromosomal sequences from diverse Indian tribal and non-tribal populations, including the Jarawa and Onge from the Andaman Islands, which are analysed within a calibrated Y-chromosomal phylogeny incorporating South Asian (in total 305 individuals) and worldwide (in total 1286 individuals) data from the 1000 Genomes Project. In contrast to the more ancient ancestry in the South than in the North that has been claimed, we detected very similar coalescence times within Northern and Southern non-tribal Indian populations. A closest neighbour analysis in the phylogeny showed that Indian populations have an affinity towards Southern European populations and that the time of divergence from these populations substantially predated the Indo-European migration into India, probably reflecting ancient shared ancestry rather than the Indo-European migration, which had little effect on Indian male lineages. Among the tribal populations, the Birhor (Austro-Asiatic-speaking) and Irula (Dravidian-speaking) are the nearest neighbours of South Asian non-tribal populations, with a common origin in the last few millennia. In contrast, the Riang (Tibeto-Burman-speaking) and Andamanese have their nearest neighbour lineages in East Asia. The Jarawa and Onge shared haplogroup D lineages with each other within the last ~7000 years, but had diverged from Japanese haplogroup D Y-chromosomes ~53000 years ago, most likely by a split from a shared ancestral population. This analysis suggests that Indian populations have complex ancestry which cannot be explained by a single expansion model.

A high-definition view of functional genetic variation from natural yeast genomes.

The question of how genetic variation in a population influences phenotypic variation and evolution is of major importance in modern biology. Yet much is still unknown about the relative functional importance of different forms of genome variation and how they are shaped by evolutionary processes. Here we address these questions by population level sequencing of 42 strains from the budding yeast Saccharomyces cerevisiae and its closest relative S. paradoxus. We find that genome content variation, in the form of presence or absence as well as copy number of genetic material, is higher within S. cerevisiae than within S. paradoxus, despite genetic distances as measured in single-nucleotide polymorphisms being vastly smaller within the former species. This genome content variation, as well as loss-of-function variation in the form of premature stop codons and frameshifting indels, is heavily enriched in the subtelomeres, strongly reinforcing the relevance of these regions to functional evolution. Genes affected by these likely functional forms of variation are enriched for functions mediating interaction with the external environment (sugar transport and metabolism, flocculation, metal transport, and metabolism). Our results and analyses provide a comprehensive view of genomic diversity in budding yeast and expose surprising and pronounced differences between the variation within S. cerevisiae and that within S. paradoxus. We also believe that the sequence data and de novo assemblies will constitute a useful resource for further evolutionary and population genomics studies.

Having older siblings is associated with gut microbiota development during early childhood.

BACKGROUND: Evidence suggests that early life infections, presence of older siblings and furred pets in the household affect the risk of developing allergic diseases through altered microbial exposure. Recently, low gut microbial diversity during infancy has also been linked with later development of allergies. We investigated whether presence of older siblings, furred pets and early life infections affected gut microbial communities at 9 and 18 months of age and whether these differences were associated with the cumulative prevalence of atopic symptoms of eczema and asthmatic bronchitis at 3 years of age. Bacterial compositions and diversity indices were determined in fecal samples collected from 114 infants in the SKOT I cohort at age 9 and 18 months by 16S rRNA gene sequencing. These were compared to the presence of older siblings, furred pets and early life infections and the cumulative prevalence of diagnosed asthmatic bronchitis and self-reported eczema at 3 years of age. RESULTS: The number of older siblings correlated positively with bacterial diversity (p = 0.030), diversity of the phyla Firmicutes (p = 0.013) and Bacteroidetes (p = 0.004) and bacterial richness (p = 0.006) at 18 months. Further, having older siblings was associated with increased relative abundance of several bacterial taxa at both 9 and 18 months of age. Compared to the effect of having siblings, presence of household furred pets and early life infections had less pronounced effects on the gut microbiota. Gut microbiota characteristics were not significantly associated with cumulative occurrence of eczema and asthmatic bronchitis during the first 3 years of life. CONCLUSIONS: Presence of older siblings is associated with increased gut microbial diversity and richness during early childhood, which could contribute to the substantiation of the hygiene hypothesis. However, no associations were found between gut microbiota and atopic symptoms of eczema and asthmatic bronchitis during early childhood and thus further studies are required to elucidate whether sibling-associated gut microbial changes influence development of allergies later in childhood.

Establishment of intestinal microbiota during early life: a longitudinal, explorative study of a large cohort of Danish infants.

Fecal samples were obtained from a cohort of 330 healthy Danish infants at 9, 18, and 36 months after birth, enabling characterization of interbacterial relationships by use of quantitative PCR targeting 31 selected bacterial 16S rRNA gene targets representing different phylogenetic levels. Nutritional parameters and measures of growth and body composition were determined and investigated in relation to the observed development in microbiota composition. We found that significant changes in the gut microbiota occurred, particularly from age 9 to 18 months, when cessation of breastfeeding and introduction of a complementary feeding induce replacement of a microbiota characterized by lactobacilli, bifidobacteria, and Enterobacteriaceae with a microbiota dominated by Clostridium spp. and Bacteroides spp. Classification of samples by a proxy enterotype based on the relative levels of Bacteroides spp. and Prevotella spp. showed that enterotype establishment occurs between 9 and 36 months. Thirty percent of the individuals shifted enterotype between 18 and 36 months. The composition of the microbiota was most pronouncedly influenced by the time of cessation of breastfeeding. From 9 to 18 months, a positive correlation was observed between the increase in body mass index and the increase of the short-chain-fatty-acid-producing clostridia, the Clostridum leptum group, and Eubacterium hallii. Considering previously established positive associations between rapid infant weight gain, early breastfeeding discontinuation, and later-life obesity, the corresponding microbial findings seen here warrant attention.

Ancient and recent origins of shared polymorphisms in yeast.

Shared genetic polymorphisms between populations and species can be ascribed to ancestral variation or to more recent gene flow. Here, we mapped shared polymorphisms in Saccharomyces cerevisiae and its sister species Saccharomyces paradoxus, which diverged 4-6 million years ago. We used a dense map of single-nucleotide diagnostic markers (mean distance 15.6 base pairs) in 1,673 sequenced S. cerevisiae isolates to catalogue 3,852 sequence blocks (≥5 consecutive markers) introgressed from S. paradoxus, with most being recent and clade-specific. The highly diverged wild Chinese S. cerevisiae lineages were depleted of introgressed blocks but retained an excess of individual ancestral polymorphisms derived from incomplete lineage sorting, perhaps due to less dramatic population bottlenecks. In the non-Chinese S. cerevisiae lineages, we inferred major hybridization events and detected cases of overlapping introgressed blocks across distinct clades due to either shared histories or convergent evolution. We experimentally engineered, in otherwise isogenic backgrounds, the introgressed PAD1-FDC1 gene pair that independently arose in two S. cerevisiae clades and revealed that it increases resistance against diverse antifungal drugs. Overall, our study retraces the histories of divergence and secondary contacts across S. cerevisiae and S. paradoxus populations and unveils a functional outcome.

High quality de novo sequencing and assembly of the Saccharomyces arboricolus genome.

BACKGROUND: Comparative genomics is a formidable tool to identify functional elements throughout a genome. In the past ten years, studies in the budding yeast Saccharomyces cerevisiae and a set of closely related species have been instrumental in showing the benefit of analyzing patterns of sequence conservation. Increasing the number of closely related genome sequences makes the comparative genomics approach more powerful and accurate. RESULTS: Here, we report the genome sequence and analysis of Saccharomyces arboricolus, a yeast species recently isolated in China, that is closely related to S. cerevisiae. We obtained high quality de novo sequence and assemblies using a combination of next generation sequencing technologies, established the phylogenetic position of this species and considered its phenotypic profile under multiple environmental conditions in the light of its gene content and phylogeny. CONCLUSIONS: We suggest that the genome of S. arboricolus will be useful in future comparative genomics analysis of the Saccharomyces sensu stricto yeasts.

Dietary xylooligosaccharide downregulates IFN-γ and the low-grade inflammatory cytokine IL-1ß systemically in mice.

Dietary carbohydrates improve growth conditions for distinct populations of bacteria that may affect mucosal and systemic immunity. In this study, we fed in a parallel experiment a 10% xylooligosaccharide (XOS)-supplemented diet or a control diet to 2 groups of male C57BL/6NTac mice for 10 wk from weaning. We found that the XOS diet significantly increased Bifidobacterium throughout the intestine compared with control-fed mice, with the highest proportions found in the ileum after XOS feeding (P < 0.001). In the intestinal epithelium, most innate immune-related genes were unaffected by XOS feeding, whereas expression of interleukin 1ß (Il1ß) (P < 0.01) and interferon γ (Ifnγ) (P < 0.05) was significantly less in blood from XOS-fed mice than from control-fed mice. In vitro treatment of blood with propionate significantly decreased Il1ß (P < 0.01), Ifnγ (P < 0.01), and interleukin 18 (Il18) (P < 0.001) expression, supporting our hypothesis that increased production of short-chain fatty acids (SCFAs) in the gut, which are transported across the intestine and into the systemic compartments, results in downregulation of low-grade inflammatory cytokines. The defensin regenerating islet-derived protein 3γ (RegIIIγ) was significantly more highly expressed in the small intestine (P < 0.01) in XOS-fed mice compared with control-fed mice, suggesting only minor contact between bifidobacteria and epithelial cells. In support of this, the SCFA-induced sodium/hydrogen exchanger isoform 3 expression tended to be greater in the XOS group than in the control group (P = 0.06), indicating an indirect SCFA-mediated antiinflammatory effect of XOS. In conclusion, XOS feeding decreases systemic inflammation, and this effect is most likely caused by higher SCFA concentrations as a result of an increased bifidobacterial saccharolytic fermentation in the entire gut and not only in the large intestine.