Chandana Basu Mallick

Genomic analyses inform on migration events during the peopling of Eurasia

High-coverage whole-genome sequence studies have so far focused on a limited number of geographically restricted populations, or been targeted at specific diseases, such as cancer. Nevertheless, the availability of high-resolution genomic data has led to the development of new methodologies for inferring population history and refuelled the debate on the mutation rate in humans. Here we present the Estonian Biocentre Human Genome Diversity Panel (EGDP), a dataset of 483 high-coverage human genomes from 148 populations worldwide, including 379 new genomes from 125 populations, which we group into diversity and selection sets. We analyse this dataset to refine estimates of continent-wide patterns of heterozygosity, long- and short-distance gene flow, archaic admixture, and changes in effective population size through time as well as for signals of positive or balancing selection. We find a genetic signature in present-day Papuans that suggests that at least 2% of their genome originates from an early and largely extinct expansion of anatomically modern humans (AMHs) out of Africa. Together with evidence from the western Asian fossil record, and admixture between AMHs and Neanderthals predating the main Eurasian expansion, our results contribute to the mounting evidence for the presence of AMHs out of Africa earlier than 75,000 years ago.

Shared and Unique Components of Human Population Structure and Genome-Wide Signals of Positive Selection in South Asia

South Asia harbors one of the highest levels genetic diversity in Eurasia, which could be interpreted as a result of its long-term large effective population size and of admixture during its complex demographic history. In contrast to Pakistani populations, populations of Indian origin have been underrepresented in previous genomic scans of positive selection and population structure. Here we report data for more than 600,000 SNP markers genotyped in 142 samples from 30 ethnic groups in India. Combining our results with other available genome-wide data, we show that Indian populations are characterized by two major ancestry components, one of which is spread at comparable frequency and haplotype diversity in populations of South and West Asia and the Caucasus. The second component is more restricted to South Asia and accounts for more than 50% of the ancestry in Indian populations. Haplotype diversity associated with these South Asian ancestry components is significantly higher than that of the components dominating the West Eurasian ancestry palette. Modeling of the observed haplotype diversities suggests that both Indian ancestry components are older than the purported Indo-Aryan invasion 3,500 YBP. Consistent with the results of pairwise genetic distances among world regions, Indians share more ancestry signals with West than with East Eurasians. However, compared to Pakistani populations, a higher proportion of their genes show regionally specific signals of high haplotype homozygosity. Among such candidates of positive selection in India are MSTN and DOK5, both of which have potential implications in lipid metabolism and the etiology of type 2 diabetes.

Different population histories of the Mundari- and Mon-Khmer-speaking Austro-Asiatic tribes inferred from the mtDNA 9-bp deletion/insertion polymorphism in Indian populations.

Length variation in the human mtDNA intergenic region between the cytochrome oxidase II (COII) and tRNA lysine (tRNAlys) genes has been widely studied in world populations. Specifically, Austronesian populations of the Pacific and Austro-Asiatic populations of southeast Asia most frequently carry the 9-bp deletion in that region implying their shared common ancestry in haplogroup B. Furthermore, multiple independent origins of the 9-bp deletion at the background of other mtDNA haplogroups has been shown in populations of Africa, Europe, Australia, and India. We have analyzed 3293 Indian individuals belonging to 58 populations, representing different caste, tribal, and religious groups, for the length variation in the 9-bp motif. The 9-bp deletion (one copy) and insertion (three copies) alleles were observed in 2.51% (2.15% deletion and 0.36% insertion) of the individuals. The maximum frequency of the deletion (45.8%) was observed in the Nicobarese in association with the haplogroup B5a D-loop motif that is common throughout southeast Asia. The low polymorphism in the D-loop sequence of the Nicobarese B5a samples suggests their recent origin and a founder effect, probably involving migration from southeast Asia. Interestingly, none of the 302 (except one Munda sample, which has 9-bp insertion) from Mundari-speaking Austro-Asiatic populations from the Indian mainland showed the length polymorphism of the 9-bp motif, pointing either to their independent origin from the Mon-Khmeric-speaking Nicobarese or to an extensive admixture with neighboring Indo-European-speaking populations. Consistent with previous reports, the Indo-European and Dravidic populations of India showed low frequency of the 9-bp deletion/insertion. More than 18 independent origins of the deletion or insertion mutation could be inferred in the phylogenetic analysis of the D-loop sequences.

Herders of Indian and European Cattle Share Their Predominant Allele for Lactase Persistence

Milk consumption and lactose digestion after weaning are exclusively human traits made possible by the continued production of the enzyme lactase in adulthood. Multiple independent mutations in a 100-bp region--part of an enhancer--approximately 14-kb upstream of the LCT gene are associated with this trait in Europeans and pastoralists from Saudi Arabia and Africa. However, a single mutation of purported western Eurasian origin accounts for much of observed lactase persistence outside Africa. Given the high levels of present-day milk consumption in India, together with archaeological and genetic evidence for the independent domestication of cattle in the Indus valley roughly 7,000 years ago, we sought to determine whether lactase persistence has evolved independently in the subcontinent. Here, we present the results of the first comprehensive survey of the LCT enhancer region in south Asia. Having genotyped 2,284 DNA samples from across the Indian subcontinent, we find that the previously described west Eurasian -13910 C>T mutation accounts for nearly all the genetic variation we observed in the 400- to 700-bp LCT regulatory region that we sequenced. Geography is a significant predictor of -13910*T allele frequency, and consistent with other genomic loci, its distribution in India follows a general northwest to southeast declining pattern, although frequencies among certain neighboring populations vary substantially. We confirm that the mutation is identical by descent to the European allele and is associated with the same>1 Mb extended haplotype in both populations.

The light skin allele of SLC24A5 in South Asians and Europeans shares identity by descent

Skin pigmentation is one of the most variable phenotypic traits in humans. A non-synonymous substitution (rs1426654) in the third exon of SLC24A5 accounts for lighter skin in Europeans but not in East Asians. A previous genome-wide association study carried out in a heterogeneous sample of UK immigrants of South Asian descent suggested that this gene also contributes significantly to skin pigmentation variation among South Asians. In the present study, we have quantitatively assessed skin pigmentation for a largely homogeneous cohort of 1228 individuals from the Southern region of the Indian subcontinent. Our data confirm significant association of rs1426654 SNP with skin pigmentation, explaining about 27% of total phenotypic variation in the cohort studied. Our extensive survey of the polymorphism in 1573 individuals from 54 ethnic populations across the Indian subcontinent reveals wide presence of the derived-A allele, although the frequencies vary substantially among populations. We also show that the geospatial pattern of this allele is complex, but most importantly, reflects strong influence of language, geography and demographic history of the populations. Sequencing 11.74 kb of SLC24A5 in 95 individuals worldwide reveals that the rs1426654-A alleles in South Asian and West Eurasian populations are monophyletic and occur on the background of a common haplotype that is characterized by low genetic diversity. We date the coalescence of the light skin associated allele at 22–28 KYA. Both our sequence and genome-wide genotype data confirm that this gene has been a target for positive selection among Europeans. However, the latter also shows additional evidence of selection in populations of the Middle East, Central Asia, Pakistan and North India but not in South India.

Maternal Footprints of Southeast Asians in North India

We have analyzed 7,137 samples from 125 different caste, tribal and religious groups of India and 99 samples from three populations of Nepal for the length variation in the COII/tRNALys region of mtDNA. Samples showing length variation were subjected to detailed phylogenetic analysis based on HVS-I and informative coding region sequence variation. The overall frequencies of the 9-bp deletion and insertion variants in South Asia were 1.9 and 0.6%, respectively. We have also defined a novel deep-rooting haplogroup M43 and identified the rare haplogroup H14 in Indian populations carrying the 9-bp deletion by complete mtDNA sequencing. Moreover, we redefined haplogroup M6 and dissected it into two well-defined subclades. The presence of haplogroups F1 and B5a in Uttar Pradesh suggests minor maternal contribution from Southeast Asia to Northern India. The occurrence of haplogroup F1 in the Nepalese sample implies that Nepal might have served as a bridge for the flow of eastern lineages to India. The presence of R6 in the Nepalese, on the other hand, suggests that the gene flow between India and Nepal has been reciprocal.

Genotype-Phenotype Study of the Middle Gangetic Plain in India Shows Association of rs2470102 with Skin Pigmentation

Our understanding of the genetics of skin pigmentation has been largely skewed towards populations of European ancestry, imparting less attention to South Asian populations, who behold huge pigmentation diversity. Here, we investigate skin pigmentation variation in a cohort of 1,167 individuals in the Middle Gangetic Plain of the Indian subcontinent. Our data confirm the association of rs1426654 with skin pigmentation among South Asians, consistent with previous studies, and also show association for rs2470102 single nucleotide polymorphism. Our haplotype analyses further help us delineate the haplotype distribution across social categories and skin color. Taken together, our findings suggest that the social structure defined by the caste system in India has a profound influence on the skin pigmentation patterns of the subcontinent. In particular, social category and associated single nucleotide polymorphisms explain about 32% and 6.4%, respectively, of the total phenotypic variance. Phylogeography of the associated single nucleotide polymorphisms studied across 52 diverse populations of the Indian subcontinent shows wide presence of the derived alleles, although their frequencies vary across populations. Our results show that both polymorphisms (rs1426654 and rs2470102) play an important role in the skin pigmentation diversity of South Asians.

The genetic structure of south Asian populations as revealed by 650 000 SNPs.

The analyses of dense marker sets covering the whole genome has revolutionised the field of (human) population genetics. Driven largely by the needs of biomedical research, these new data are helping to unveil our demographic past, exemplified by the study of mtDNA and Y-chromosome variation during the past ~20 years. We have analysed (Illumina 650K SNPs) over 320 new samples from South and Central Asia and the Caucasus, together with the publicly available databases (HGDP panel and our published data set of ~600 Eurasian samples) and illustrated the power of full genome analyses by addressing two specific questions. (i) What is the nature of genetic continuity and discontinuity between South Asia, Middle East and Central Asia? (ii) What are the genetic origins of the Munda speakers of India? We use principal component and structure-like analyses to reveal the structure in the genome wide SNP data. The most striking feature of the genetic structure of South Asian populations is the clear separation of the Indus valley and southern India populations. The genetic component prevalent in the latter region is marginal in the former and absent outside South Asia. By contrast, the component ubiquitous to Indus valley is also present (~30 - 40%) among Indo-European speakers from Ganges valley and Dravidic speakers in southern India. Furthermore, this component can also be found in Central Asia and the Caucasus as well as in Middle East. We explored possibilities to identify the source region for this genetic component. Alternative models put the origins of Munda languages speakers either in South Asia (the Munda speakers sport exclusively autochthonous South Asian mtDNA variants) or in Southeast Asia, where the other Austro Asiatic languages have spread. Y-chromosome variation supports the latter model through sharing of hg O2a in both regions. We show that in addition to the dominant ancestry component being shared between the Indian Dravidic and Munda speakers, up to 30% of Munda speakers retain an ancestry component otherwise prevalent in East Asia. There is no widespread sign of South Asian ancestry component in Southeast Asia. This provides genomic support to the model by which Indian Austro-Asiatic populations derive from dispersal from Southeast/East Asia, followed by an extensive admixture with local Indian populations.

The influences of genes, the environment, and social factors on the evolution of skin color diversity in India

Skin color is a highly visible and variable trait across human populations. It is not yet clear how evolutionary forces interact to generate phenotypic diversity. Here, we sought to unravel through an integrative framework the role played by three factors—demography and migration, sexual selection, and natural selection—in driving skin color diversity in India.