Whole-genome sequencing reveals a complex African population demographic history and signatures of local adaptation.

We conduct high coverage (>30×) whole-genome sequencing of 180 individuals from 12 indigenous African populations. We identify millions of unreported variants, many predicted to be functionally important. We observe that the ancestors of southern African San and central African rainforest hunter-gatherers (RHG) diverged from other populations >200 kya and maintained a large effective population size. We observe evidence for ancient population structure in Africa and for multiple introgression events from "ghost" populations with highly diverged genetic lineages. Although currently geographically isolated, we observe evidence for gene flow between eastern and southern Khoesan-speaking hunter-gatherer populations lasting until ∼12 kya. We identify signatures of local adaptation for traits related to skin color, immune response, height, and metabolic processes. We identify a positively selected variant in the lightly pigmented San that influences pigmentation in vitro by regulating the enhancer activity and gene expression of PDPK1.

Remarkably Low KIR and HLA Diversity in Amerindians Reveals Signatures of Strong Purifying Selection Shaping the Centromeric KIR Region.

The killer-cell immunoglobulin-like receptors (KIR) recognize human leukocyte antigen (HLA) molecules to regulate the cytotoxic and inflammatory responses of natural killer cells. KIR genes are encoded by a rapidly evolving gene family on chromosome 19 and present an unusual variation of presence and absence of genes and high allelic diversity. Although many studies have associated KIR polymorphism with susceptibility to several diseases over the last decades, the high-resolution allele-level haplotypes have only recently started to be described in populations. Here, we use a highly innovative custom next-generation sequencing method that provides a state-of-art characterization of KIR and HLA diversity in 706 individuals from eight unique South American populations: five Amerindian populations from Brazil (three Guarani and two Kaingang); one Amerindian population from Paraguay (Aché); and two urban populations from Southern Brazil (European and Japanese descendants from Curitiba). For the first time, we describe complete high-resolution KIR haplotypes in South American populations, exploring copy number, linkage disequilibrium, and KIR-HLA interactions. We show that all Amerindians analyzed to date exhibit the lowest numbers of KIR-HLA interactions among all described worldwide populations, and that 83-97% of their KIR-HLA interactions rely on a few HLA-C molecules. Using multiple approaches, we found signatures of strong purifying selection on the KIR centromeric region, which codes for the strongest NK cell educator receptors, possibly driven by the limited HLA diversity in these populations. Our study expands the current knowledge of KIR genetic diversity in populations to understand KIR-HLA coevolution and its impact on human health and survival.

Ancestry, diversity, and genetics of health-related traits in African-derived communities (quilombos) from Brazil.

Brazilian quilombos are communities formed by enslaved Africans and their descendants all over the country during slavery and shortly after its abolition. Quilombos harbor a great fraction of the largely unknown genetic diversity of the African diaspora in Brazil. Thus, genetic studies in quilombos have the potential to provide important insights not only into the African roots of the Brazilian population but also into the genetic bases of complex traits and human adaptation to diverse environments. This review summarizes the main results of genetic studies performed on quilombos so far. Here, we analyzed the patterns of African, Amerindian, European, and subcontinental ancestry (within Africa) of quilombos from the five different geographic regions of Brazil. In addition, uniparental markers (from the mtDNA and the Y chromosome) studies are analyzed together to reveal demographic processes and sex-biased admixture that occurred during the formation of these unique populations. Lastly, the prevalence of known malaria-adaptive African mutations and other African-specific variants discovered in quilombos, as well as the genetic bases of health-related traits, are discussed here, together with their implication for the health of populations of African descent.

Effects of MASP2 haplotypes and MASP-2 levels in hepatitis C-infected patients.

Mannan-binding lectin (MBL) and MBL-associated serine protease 2 (MASP-2) are components of the lectin pathway, which activate the complement system after binding to the HCV structural proteins E1 and E2. We haplotyped 11 MASP2 polymorphisms in 103 HCV patients and 205 controls and measured MASP-2 levels in 67 HCV patients and 77 controls to better understand the role of MASP-2 in hepatitis C susceptibility and disease severity according to viral genotype and fibrosis levels. The haplotype block MASP2*ARDP was associated with protection against HCV infection (OR = 0.49, p = .044) and lower MASP-2 levels in controls (p = .021), while haplotype block AGTDVRC was significantly increased in patients (OR = 7.58, p = .003). MASP-2 levels were lower in patients than in controls (p < .001) and in patients with viral genotype 1 or 4 (poor responders to treatment) than genotype 3 (p = .022) and correlated inversely with the levels of alkaline phosphatase, especially in individuals with fibrosis 3 or 4 (R = -.7, p = .005). MASP2 gene polymorphisms modulate basal gene expression, which may influence the quality of complement response against HCV. MASP-2 levels decrease during chronic disease, independently of MASP2 genotypes, most probably due to consumption and attenuation mechanisms of viral origin and by the reduced liver function, the site of MASP-2 production.

Single Nucleotide Polymorphism in KIR2DL1 Is Associated With HLA-C Expression in Global Populations.

Regulation of NK cell activity is mediated through killer-cell immunoglobulin-like receptors (KIR) ability to recognize human leukocyte antigen (HLA) class I molecules as ligands. Interaction of KIR and HLA is implicated in viral infections, autoimmunity, and reproduction and there is growing evidence of the coevolution of these two independently segregating gene families. By leveraging KIR and HLA-C data from 1000 Genomes consortium we observed that the KIR2DL1 variant rs2304224*T is associated with lower expression of HLA-C in individuals carrying the ligand HLA-C2 (p = 0.0059). Using flow cytometry, we demonstrated that this variant is also associated with higher expression of KIR2DL1 on the NK cell surface (p = 0.0002). Next, we applied next generation sequencing to analyze KIR2DL1 sequence variation in 109 Euro and 75 Japanese descendants. Analyzing the extended haplotype homozygosity, we show signals of positive selection for rs4806553*G and rs687000*G, which are in linkage disequilibrium with rs2304224*T. Our results suggest that lower expression of HLA-C2 ligands might be compensated for higher expression of the receptor KIR2DL1 and bring new insights into the coevolution of KIR and HLA.

African evolutionary history inferred from whole genome sequence data of 44 indigenous African populations.

BACKGROUND: Africa is the origin of modern humans within the past 300 thousand years. To infer the complex demographic history of African populations and adaptation to diverse environments, we sequenced the genomes of 92 individuals from 44 indigenous African populations. RESULTS: Genetic structure analyses indicate that among Africans, genetic ancestry is largely partitioned by geography and language, though we observe mixed ancestry in many individuals, consistent with both short- and long-range migration events followed by admixture. Phylogenetic analysis indicates that the San genetic lineage is basal to all modern human lineages. The San and Niger-Congo, Afroasiatic, and Nilo-Saharan lineages were substantially diverged by 160 kya (thousand years ago). In contrast, the San and Central African rainforest hunter-gatherer (CRHG), Hadza hunter-gatherer, and Sandawe hunter-gatherer lineages were diverged by ~ 120-100 kya. Niger-Congo, Nilo-Saharan, and Afroasiatic lineages diverged more recently by ~ 54-16 kya. Eastern and western CRHG lineages diverged by ~ 50-31 kya, and the western CRHG lineages diverged by ~ 18-12 kya. The San and CRHG populations maintained the largest effective population size compared to other populations prior to 60 kya. Further, we observed signatures of positive selection at genes involved in muscle development, bone synthesis, reproduction, immune function, energy metabolism, and cell signaling, which may contribute to local adaptation of African populations. CONCLUSIONS: We observe high levels of genomic variation between ethnically diverse Africans which is largely correlated with geography and language. Our study indicates ancient population substructure and local adaptation of Africans.

Correction to: African evolutionary history inferred from whole genome sequence data of 44 indigenous African populations.

Following publication of the original article [1], a typographical error in the formula for calculating di in the "Scans for local adaptation" subsection in the Method section, was identified. The correct formula should be.

Loci associated with skin pigmentation identified in African populations.

Despite the wide range of skin pigmentation in humans, little is known about its genetic basis in global populations. Examining ethnically diverse African genomes, we identify variants in or near SLC24A5, MFSD12, DDB1, TMEM138, OCA2, and HERC2 that are significantly associated with skin pigmentation. Genetic evidence indicates that the light pigmentation variant at SLC24A5 was introduced into East Africa by gene flow from non-Africans. At all other loci, variants associated with dark pigmentation in Africans are identical by descent in South Asian and Australo-Melanesian populations. Functional analyses indicate that MFSD12 encodes a lysosomal protein that affects melanogenesis in zebrafish and mice, and that mutations in melanocyte-specific regulatory regions near DDB1/TMEM138 correlate with expression of ultraviolet response genes under selection in Eurasians.

MBL-associated serine proteases (MASPs) and infectious diseases.

The lectin pathway of the complement system has a pivotal role in the defense against infectious organisms. After binding of mannan-binding lectin (MBL), ficolins or collectin 11 to carbohydrates or acetylated residues on pathogen surfaces, dimers of MBL-associated serine proteases 1 and 2 (MASP-1 and MASP-2) activate a proteolytic cascade, which culminates in the formation of the membrane attack complex and pathogen lysis. Alternative splicing of the pre-mRNA encoding MASP-1 results in two other products, MASP-3 and MAp44, which regulate activation of the cascade. A similar mechanism allows the gene encoding MASP-2 to produce the truncated MAp19 protein. Polymorphisms in MASP1 and MASP2 genes are associated with protein serum levels and functional activity. Since the first report of a MASP deficiency in 2003, deficiencies in lectin pathway proteins have been associated with recurrent infections and several polymorphisms were associated with the susceptibility or protection to infectious diseases. In this review, we summarize the findings on the role of MASP polymorphisms and serum levels in bacterial, viral and protozoan infectious diseases.