Determining the impact of uncharacterized inversions in the human genome by droplet digital PCR.

Despite the interest in characterizing genomic variation, the presence of large repeats at the breakpoints hinders the analysis of many structural variants. This is especially problematic for inversions, since there is typically no gain or loss of DNA. Here, we tested novel linkage-based droplet digital PCR (ddPCR) assays to study 20 inversions ranging from 3.1 to 742 kb flanked by inverted repeats (IRs) up to 134 kb long. Of those, we validated 13 inversions predicted by different genome-wide techniques. In addition, we obtained new experimental human population information across 95 African, European, and East Asian individuals for 16 inversions, including four already validated variants without high-throughput genotyping methods. Through comparison with previous data, independent replicates and both inversion breakpoints, we demonstrate that the technique is highly accurate and reproducible. Most studied inversions are widespread across continents, and their frequency is negatively correlated with genetic length. Moreover, all except two show clear signs of being recurrent, and we could better define the factors affecting recurrence levels and estimate the inversion rate across the genome. Finally, the generated genotypes have allowed us to check inversion functional effects, validating gene expression differences reported before for two inversions and finding new candidate associations. Therefore, the developed methodology makes it possible to screen these and other complex genomic variants quickly in a large number of samples for the first time, highlighting the importance of direct genotyping to assess their potential consequences and clinical implications.

Oral microbiome homogeneity across diverse human groups from southern Africa: first results from southwestern Angola and Zimbabwe.

BACKGROUND: While the human oral microbiome is known to play an important role in systemic health, its average composition and diversity patterns are still poorly understood. To gain better insights into the general composition of the microbiome on a global scale, the characterization of microbiomes from a broad range of populations, including non-industrialized societies, is needed. Here, we used the portion of non-human reads obtained through an expanded exome capture sequencing approach to characterize the saliva microbiomes of 52 individuals from eight ethnolinguistically diverse southern African populations from Angola (Kuvale, Kwepe, Himba, Tjimba, Kwisi, Twa, !Xun) and Zimbabwe (Tshwa), including foragers, food-producers, and peripatetic groups (low-status communities who provide services to their dominant neighbors). RESULTS: Our results indicate that neither host genetics nor livelihood seem to influence the oral microbiome profile, with Neisseria, Streptococcus, Prevotella, Rothia, and Porphyromonas being the five most frequent genera in southern African groups, in line with what has been shown for other human populations. However, we found that some Tshwa and Twa individuals display an enrichment of pathogenic genera from the Enterobacteriaceae family (i.e. Enterobacter, Citrobacter, Salmonella) of the Proteobacteria phylum, probably reflecting deficient sanitation and poor health conditions associated with social marginalization. CONCLUSIONS: Taken together, our results suggest that socio-economic status, rather than ethnolinguistic affiliation or subsistence mode, is a key factor in shaping the salivary microbial profiles of human populations in southern Africa.

Along the Indian Ocean Coast: Genomic Variation in Mozambique Provides New Insights into the Bantu Expansion.

The Bantu expansion, which started in West Central Africa around 5,000 BP, constitutes a major migratory movement involving the joint spread of peoples and languages across sub-Saharan Africa. Despite the rich linguistic and archaeological evidence available, the genetic relationships between different Bantu-speaking populations and the migratory routes they followed during various phases of the expansion remain poorly understood. Here, we analyze the genetic profiles of southwestern and southeastern Bantu-speaking peoples located at the edges of the Bantu expansion by generating genome-wide data for 200 individuals from 12 Mozambican and 3 Angolan populations using ∼1.9 million autosomal single nucleotide polymorphisms. Incorporating a wide range of available genetic data, our analyses confirm previous results favoring a "late split" between West and East Bantu speakers, following a joint passage through the rainforest. In addition, we find that Bantu speakers from eastern Africa display genetic substructure, with Mozambican populations forming a gradient of relatedness along a North-South cline stretching from the coastal border between Kenya and Tanzania to South Africa. This gradient is further associated with a southward increase in genetic homogeneity, and involved minimum admixture with resident populations. Together, our results provide the first genetic evidence in support of a rapid North-South dispersal of Bantu peoples along the Indian Ocean Coast, as inferred from the distribution and antiquity of Early Iron Age assemblages associated with the Kwale archaeological tradition.

Detailed analysis of inversions predicted between two human genomes: errors, real polymorphisms, and their origin and population distribution.

The growing catalogue of structural variants in humans often overlooks inversions as one of the most difficult types of variation to study, even though they affect phenotypic traits in diverse organisms. Here, we have analysed in detail 90 inversions predicted from the comparison of two independently assembled human genomes: the reference genome (NCBI36/HG18) and HuRef. Surprisingly, we found that two thirds of these predictions (62) represent errors either in assembly comparison or in one of the assemblies, including 27 misassembled regions in HG18. Next, we validated 22 of the remaining 28 potential polymorphic inversions using different PCR techniques and characterized their breakpoints and ancestral state. In addition, we determined experimentally the derived allele frequency in Europeans for 17 inversions (DAF = 0.01-0.80), as well as the distribution in 14 worldwide populations for 12 of them based on the 1000 Genomes Project data. Among the validated inversions, nine have inverted repeats (IRs) at their breakpoints, and two show nucleotide variation patterns consistent with a recurrent origin. Conversely, inversions without IRs have a unique origin and almost all of them show deletions or insertions at the breakpoints in the derived allele mediated by microhomology sequences, which highlights the importance of mechanisms like FoSTeS/MMBIR in the generation of complex rearrangements in the human genome. Finally, we found several inversions located within genes and at least one candidate to be positively selected in Africa. Thus, our study emphasizes the importance of careful analysis and validation of large-scale genomic predictions to extract reliable biological conclusions.

Matriclans shape populations: Insights from the Angolan Namib Desert into the maternal genetic history of southern Africa.

OBJECTIVES: Southern Angola is a poorly studied region, inhabited by populations that have been associated with different migratory movements into southern Africa. Apart from Kx'a-speaking San foragers and Bantu-speaking pastoralists, ethnographic and linguistic studies have suggested the existence of an enigmatic array of pre-Bantu communities, like the Kwepe (formerly Khoe-Kwadi speakers), Twa and Kwisi. Here, we evaluate previous peopling hypotheses by assessing the relationships between different southern Angolan populations, based on newly collected linguistic data and complete mtDNA genomes. MATERIALS AND METHODS: We analyzed 295 complete mtDNA genomes and linguistic data from seven groups from the Namib Desert (Himba, Kuvale, Tjimba, Twa, Kwisi, Kwepe) and Kunene Province (!Xun), placing special emphasis on the evaluation of the genealogical consistency of the matriclanic system that characterizes most of these groups. RESULTS: We found that the maternal genetic structure of all groups from the Namib Desert was strongly shaped by the consistency of their matriclanic system. The tracking of the maternal heritage enhanced population differentiation by genetic drift and is likely to have caused the divergent mtDNA profiles of the Kwepe, Twa, and Kwisi, who probably formed a single population within the spectrum of Bantu genetic variation. Model-based analyses further suggest that the dominant pastoral groups Kuvale and Himba may be grouped into a Bantu proto-population which also included the ancestors of present-day Tjimba and Herero, as well as the Khoe-Kwadi speaking Damara foragers from Namibia. DISCUSSION: The view from southwestern Angola offers a new perspective on the populating history of southern Africa and the Bantu expansions by showing that social stratification and different subsistence patterns are not always indicative of remnant groups, but may reflect Bantu-internal variation and ethnogenesis.

Functional Impact and Evolution of a Novel Human Polymorphic Inversion That Disrupts a Gene and Creates a Fusion Transcript.

Despite many years of study into inversions, very little is known about their functional consequences, especially in humans. A common hypothesis is that the selective value of inversions stems in part from their effects on nearby genes, although evidence of this in natural populations is almost nonexistent. Here we present a global analysis of a new 415-kb polymorphic inversion that is among the longest ones found in humans and is the first with clear position effects. This inversion is located in chromosome 19 and has been generated by non-homologous end joining between blocks of transposable elements with low identity. PCR genotyping in 541 individuals from eight different human populations allowed the detection of tag SNPs and inversion genotyping in multiple populations worldwide, showing that the inverted allele is mainly found in East Asia with an average frequency of 4.7%. Interestingly, one of the breakpoints disrupts the transcription factor gene ZNF257, causing a significant reduction in the total expression level of this gene in lymphoblastoid cell lines. RNA-Seq analysis of the effects of this expression change in standard homozygotes and inversion heterozygotes revealed distinct expression patterns that were validated by quantitative RT-PCR. Moreover, we have found a new fusion transcript that is generated exclusively from inverted chromosomes around one of the breakpoints. Finally, by the analysis of the associated nucleotide variation, we have estimated that the inversion was generated ~40,000-50,000 years ago and, while a neutral evolution cannot be ruled out, its current frequencies are more consistent with those expected for a deleterious variant, although no significant association with phenotypic traits has been found so far.

Validation and genotyping of multiple human polymorphic inversions mediated by inverted repeats reveals a high degree of recurrence.

In recent years different types of structural variants (SVs) have been discovered in the human genome and their functional impact has become increasingly clear. Inversions, however, are poorly characterized and more difficult to study, especially those mediated by inverted repeats or segmental duplications. Here, we describe the results of a simple and fast inverse PCR (iPCR) protocol for high-throughput genotyping of a wide variety of inversions using a small amount of DNA. In particular, we analyzed 22 inversions predicted in humans ranging from 5.1 kb to 226 kb and mediated by inverted repeat sequences of 1.6-24 kb. First, we validated 17 of the 22 inversions in a panel of nine HapMap individuals from different populations, and we genotyped them in 68 additional individuals of European origin, with correct genetic transmission in ∼ 12 mother-father-child trios. Global inversion minor allele frequency varied between 1% and 49% and inversion genotypes were consistent with Hardy-Weinberg equilibrium. By analyzing the nucleotide variation and the haplotypes in these regions, we found that only four inversions have linked tag-SNPs and that in many cases there are multiple shared SNPs between standard and inverted chromosomes, suggesting an unexpected high degree of inversion recurrence during human evolution. iPCR was also used to check 16 of these inversions in four chimpanzees and two gorillas, and 10 showed both orientations either within or between species, providing additional support for their multiple origin. Finally, we have identified several inversions that include genes in the inverted or breakpoint regions, and at least one disrupts a potential coding gene. Thus, these results represent a significant advance in our understanding of inversion polymorphism in human populations and challenge the common view of a single origin of inversions, with important implications for inversion analysis in SNP-based studies.

Food and pathogen adaptations in the Angolan Namib desert: Tracing the spread of lactase persistence and human African trypanosomiasis resistance into southwestern Africa.

OBJECTIVES: We investigated the frequency distribution and haplotype diversity of human African trypanosomiasis (HAT) resistance and lactase persistence (LP) variants in populations from the Angolan Namib to trace the spread of these genetic adaptations into southwestern Africa. MATERIALS AND METHODS: We resequenced two fragments of the LCT enhancer and the APOL1 gene and genotyped flanking short tandem repeat loci in six groups with different subsistence traditions living in the Angolan Namib, and in a comparative dataset including other populations from Africa and Europe. LP in the Angolan Namib is represented by the -14010*C allele, which is associated with a predominant haplotype that is shared with other southern and eastern African populations. While LP was found to be more frequent in foragers than in pastoralists, the frequencies of the two APOL1 variants associated with HAT-resistance (G1 and G2) did not differ between the two groups. The G1 allele is mostly associated with a single widespread haplotype. The G2 allele is linked to several haplotypes that are molecularly related to haplotypes found in other African Bantu-speaking populations. The putatively archaic G3 variant displayed more intra-allelic diversity in Africa than in Europe. DISCUSSION: The LP adaptation was carried to southern Africa by non-Bantu speaking pastoralists from eastern Africa, but an obvious link between its presence in southern Angola and groups speaking languages of the Khoe-Kwadi family, as previously found in other areas, could not be confirmed. The presence of APOL1 variants G1 and G2 is linked to the Bantu expansions. Our results suggest that the G3 variant was retained in modern humans by incomplete lineage sorting.

Uncovering adaptive evolution in the human lineage.

BACKGROUND: The recent increase in human polymorphism data, together with the availability of genome sequences from several primate species, provides an unprecedented opportunity to investigate how natural selection has shaped human evolution. RESULTS: We compared human branch-specific substitutions with variation data in the current human population to measure the impact of adaptive evolution on human protein coding genes. The use of single nucleotide polymorphisms (SNPs) with high derived allele frequencies (DAFs) minimized the influence of segregating slightly deleterious mutations and improved the estimation of the number of adaptive sites. Using DAF ≥ 60% we showed that the proportion of adaptive substitutions is 0.2% in the complete gene set. However, the percentage rose to 40% when we focused on genes that are specifically accelerated in the human branch with respect to the chimpanzee branch, or on genes that show signatures of adaptive selection at the codon level by the maximum likelihood based branch-site test. In general, neural genes are enriched in positive selection signatures. Genes with multiple lines of evidence of positive selection include taxilin beta, which is involved in motor nerve regeneration and syntabulin, and is required for the formation of new presynaptic boutons. CONCLUSIONS: We combined several methods to detect adaptive evolution in human coding sequences at a genome-wide level. The use of variation data, in addition to sequence divergence information, uncovered previously undetected positive selection signatures in neural genes.

Apolipoprotein E/C1/C4/C2 gene cluster diversity in two native Andean populations: Aymaras and Quechuas.

The APOE/C1/C4/C2 gene cluster presents high relevance in lipid metabolism and, therefore, has important epidemiological implications. Here, we study for the first time the variation patterns of 25 polymorphisms (10 short tandem repeats, STRs, and 15 single nucleotide polymorphismas, SNPs) in two native Andean samples from Bolivia (45 Aymaras and 45 Quechuas) as well as one European sample (n = 41) as external reference. We estimated diversity parameters, linkage disequilibrium patterns, population structure, and possible selective effects. In general, diversity was low and could be partly attributed to selection (probably due to its physiological importance), since the APOE/C1/C4/C2 region was highly conserved compared to the flanking genes in both Bolivians and Europeans. Moreover, the lower gene diversity in Bolivians compared to Europeans for some markers might indicate different demographic histories. Regarding the APOE isoforms, in addition to ɛ3 (94%) and ɛ4 (5%), isoform ɛ2 (1%) was also detected in Bolivians. In relation to previous hypotheses, our results support that genetic drift or founder effects rather than selection for increased cholesterol absorption are the main factors that have shaped the distribution of APOE isoforms observed in South America.

mtDNA and Y-chromosome diversity in Aymaras and Quechuas from Bolivia: different stories and special genetic traits of the Andean Altiplano populations.

Two Bolivian samples belonging to the two main Andean linguistic groups (Aymaras and Quechuas) were studied for mtDNA and Y-chromosome uniparental markers to evaluate sex-specific differences and give new insights into the demographic processes of the Andean region. mtDNA-coding polymorphisms, HVI-HVII control regions, 17 Y-STRs, and three SNPs were typed in two well-defined populations with adequate size samples. The two Bolivian samples showed more genetic differences for the mtDNA than for the Y-chromosome. For the mtDNA, 81% of Aymaras and 61% of Quechuas presented haplogroup B2. Native American Y-chromosomes were found in 97% of Aymaras (89% hg Q1a3a and 11% hg Q1a3*) and 78% of Quechuas (100% hg Q1a3a). Our data revealed high diversity values in the two populations, in agreement with other Andean studies. The comparisons with the available literature for both sets of markers indicated that the central Andean area is relatively homogeneous. For mtDNA, the Aymaras seemed to have been more isolated throughout time, maintaining their genetic characteristics, while the Quechuas have been more permeable to the incorporation of female foreigners and Peruvian influences. On the other hand, male mobility would have been widespread across the Andean region according to the homogeneity found in the area. Particular genetic characteristics presented by both samples support a past common origin of the Altiplano populations in the ancient Aymara territory, with independent, although related histories, with Peruvian (Quechuas) populations.

Characterization of the Oral Microbiome of Medicated Type-2 Diabetes Patients.

Type 2 diabetes mellitus (T2DM) is a chronic metabolic disease that is becoming a significant global health care problem. Several studies have shown that people with diabetes are more susceptible to oral problems, such as periodontitis and, although the causes are still inconclusive, oral microbiota is considered to play a major role in oral health. This study aimed to characterize the oral microbiome of a sample representing T2DM patients from Portugal and exploit potential associations between some microorganisms and variables like teeth brushing, smoking habits, average blood sugar levels, medication and nutrient intake. By sequencing the hypervariable regions V3-V4 of the 16S rRNA gene in 50 individuals belonging to a group of diabetes patients and a control group, we found a total of 232 taxa, from which only 65% were shared between both groups. No differences were found in terms of alpha and beta diversity between categories. We did not find significant differences in the oral microbiome profiles of control and diabetes patients. Only the class Synergistia and the genus TG5, which are related to periodontitis, were statistically more frequent in the control group. The similar microbiome profiles of medicated diabetics and the control group indicates that the relationship between the T2DM and the oral microbiome might be more related to either the lifestyle/diet rather than diabetes per se. Moreover, this study provides, for the first time, insights into the oral microbiome of a population with a high prevalence of diabetes.

The Genes of Freedom: Genome-Wide Insights into Marronage, Admixture and Ethnogenesis in the Gulf of Guinea.

The forced migration of millions of Africans during the Atlantic Slave Trade led to the emergence of new genetic and linguistic identities, thereby providing a unique opportunity to study the mechanisms giving rise to human biological and cultural variation. Here we focus on the archipelago of São Tomé and Príncipe in the Gulf of Guinea, which hosted one of the earliest plantation societies relying exclusively on slave labor. We analyze the genetic variation in 25 individuals from three communities who speak distinct creole languages (Forros, Principenses and Angolares), using genomic data from expanded exomes in combination with a contextual dataset from Europe and Africa, including newly generated data from 28 Bantu speakers from Angola. Our findings show that while all islanders display mixed contributions from the Gulf of Guinea and Angola, the Angolares are characterized by extreme genetic differentiation and inbreeding, consistent with an admixed maroon isolate. In line with a more prominent Bantu contribution to their creole language, we additionally found that a previously reported high-frequency Y-chromosome haplotype in the Angolares has a likely Angolan origin, suggesting that their genetic, linguistic and social characteristics were influenced by a small group of dominant men who achieved disproportionate reproductive success.

Different evolutionary histories of the coagulation factor VII gene in human populations?

Immoderate blood clotting constitutes a risk factor for cardiovascular disease in modern industrialised societies, but is believed to have conferred a survival advantage, i.e. faster recovery from bleeding, on our ancestors. Here, we investigate the evolutionary history of the Coagulation Factor VII gene (F7) by analysing five cardiovascular-risk-associated mutations from the F7 promoter and nine neutral polymorphisms (six SNPs and three microsatellites) from the flanking region in 16 populations from the broader Mediterranean region, South Saharan Africa and Bolivia (687 individuals in total). Population differentiation and selection tests were performed and linkage disequilibrium patterns were investigated. In all samples, no linkage disequilibrium between adjacent F7 promoter mutations -402 and -401 was observed. No selection signals were detected in any of the samples from the broader Mediterranean region and South Saharan Africa, while some of the data suggested a potential signal of positive selection for the F7 promoter in the Native American samples from Bolivia. In conclusion, our data suggest, although do not prove, different evolutionary histories in the F7 promoter region between Mediterraneans and Amerindians.

Mixed origin of the current Tunisian population from the analysis of Alu and Alu/STR compound systems.

During successive historical periods, Tunisia has been a crossroads of multiple civilizations and their corresponding key population movements. The aim of this study was to provide genetic information relating to the mixed origin of the Tunisian population, and to analyze its genetic relationship with other North African and Mediterranean populations. A set of 16 Alu and 3 Alu/STR compound systems has been analyzed in 268 autochthonous Tunisians from the north-center and the south of the country. Our two sampled populations showed no significant differentiation from one another in any of the three Alu/STR compound systems, whereas the analysis of the 16 Alu markers revealed a significant genetic differentiation between them. A sub-Saharan component shown by the three Alu/STR combinations is more noticeable in our north-center sample than in that of the south. The presence of two Alu/STR combinations specific to North African ancestral populations also suggests that the ancient Berber component is relatively more substantial in the north and center regions than in the south. Our Tunisian samples cluster together with other Berber samples from Morocco and Algeria, underpinning the genetic similarity among North Africans regardless of their current linguistic status (Berber or Arabic).

Population relationships in the Mediterranean revealed by autosomal genetic data (Alu and Alu/STR compound systems).

The variation of 18 Alu polymorphisms and 3 linked STRs was determined in 1,831 individuals from 15 Mediterranean populations to analyze the relationships between human groups in this geographical region and provide a complementary perspective to information from studies based on uniparental markers. Patterns of population diversity revealed by the two kinds of markers examined were different from one another, likely in relation to their different mutation rates. Therefore, while the Alu biallelic variation underlies general heterogeneity throughout the whole Mediterranean region, the combined use of Alu and STR points to a considerable genetic differentiation between the two Mediterranean shores, presumably strengthened by a considerable sub-Saharan African genetic contribution in North Africa (around 13% calculated from Alu markers). Gene flow analysis confirms the permeability of the Sahara to human passage along with the existence of trans-Mediterranean interchanges. Two specific Alu/STR combinations-CD4 110(-) and DM 107(-)-detected in all North African samples, the Iberian Peninsula, Greece, Turkey, and some Mediterranean islands suggest an ancient genetic background of current Mediterranean peoples.

Autosomal and X chromosome Alu insertions in Bolivian Aymaras and Quechuas: two languages and one genetic pool.

Thirty-two polymorphic Alu insertions (18 autosomal and 14 from the X chromosome) were studied in 192 individuals from two Amerindian populations of the Bolivian Altiplano (Aymara and Quechua speakers: the two main Andean linguistic groups), to provide relevant information about their genetic relationships and demographic processes. The main objective was to determine from genetic data whether the expansion of the Quechua language into Bolivia could be associated with demographic (Inca migration of Quechua-speakers from Peru into Bolivia) or cultural (language imposition by the Inca Empire) processes. Allele frequencies were used to assess the genetic relationships between these two linguistic groups. Our results indicated that the two Bolivian samples showed a high genetic similarity for both sets of markers and were clearly differentiated from the two Peruvian Quechua samples available in the literature. Additionally, our data were compared with the available literature to determine the genetic and linguistic structure, and East-West differentiation in South America. The close genetic relationship between the two Bolivian samples and their differentiation from the Quechua-speakers from Peru suggests that the Quechua language expansion in Bolivia took place without any important demographic contribution. Moreover, no clear geographical or linguistic structure was found for the Alu variation among South Amerindians.

Evolutionary and functional impact of common polymorphic inversions in the human genome.

Inversions are one type of structural variants linked to phenotypic differences and adaptation in multiple organisms. However, there is still very little information about polymorphic inversions in the human genome due to the difficulty of their detection. Here, we develop a new high-throughput genotyping method based on probe hybridization and amplification, and we perform a complete study of 45 common human inversions of 0.1-415 kb. Most inversions promoted by homologous recombination occur recurrently in humans and great apes and they are not tagged by SNPs. Furthermore, there is an enrichment of inversions showing signatures of positive or balancing selection, diverse functional effects, such as gene disruption and gene-expression changes, or association with phenotypic traits. Therefore, our results indicate that the genome is more dynamic than previously thought and that human inversions have important functional and evolutionary consequences, making possible to determine for the first time their contribution to complex traits.