Genomic reconstruction of the history of extant populations of India reveals five distinct ancestral components and a complex structure.

India, occupying the center stage of Paleolithic and Neolithic migrations, has been underrepresented in genome-wide studies of variation. Systematic analysis of genome-wide data, using multiple robust statistical methods, on (i) 367 unrelated individuals drawn from 18 mainland and 2 island (Andaman and Nicobar Islands) populations selected to represent geographic, linguistic, and ethnic diversities, and (ii) individuals from populations represented in the Human Genome Diversity Panel (HGDP), reveal four major ancestries in mainland India. This contrasts with an earlier inference of two ancestries based on limited population sampling. A distinct ancestry of the populations of Andaman archipelago was identified and found to be coancestral to Oceanic populations. Analysis of ancestral haplotype blocks revealed that extant mainland populations (i) admixed widely irrespective of ancestry, although admixtures between populations was not always symmetric, and (ii) this practice was rapidly replaced by endogamy about 70 generations ago, among upper castes and Indo-European speakers predominantly. This estimated time coincides with the historical period of formulation and adoption of sociocultural norms restricting intermarriage in large social strata. A similar replacement observed among tribal populations was temporally less uniform.

Signatures of natural selection are not uniform across genes of innate immune system, but purifying selection is the dominant signature.

We tested the opposing views concerning evolution of genes of the innate immune system that (i) being evolutionary ancient, the system may have been highly optimized by natural selection and therefore should be under purifying selection, and (ii) the system may be plastic and continuing to evolve under balancing selection. We have resequenced 12 important innate-immunity genes (CAMP, DEFA4, DEFA5, DEFA6, DEFB1, MBL2, and TLRs 1, 2, 4, 5, 6, and 9) in healthy volunteers (n = 171) recruited from a region of India with high microbial load. We have compared these data with those of European-Americans (EUR) and African-Americans (AFR). We have found that most of the human haplotypes are many mutational steps away from the ancestral (chimpanzee) haplotypes, indicating that humans may have had to adapt to new pathogens. The haplotype structures in India are significantly different from those of EUR and AFR populations, indicating local adaptation to pathogens. In these genes, there is (i) generally an excess of rare variants, (ii) high, but variable, degrees of extended haplotype homozygosity, (iii) low tolerance to nonsynonymous changes, (iv) essentially one or a few high-frequency haplotypes, with star-like phylogenies of other infrequent haplotypes radiating from the modal haplotypes. Purifying selection is the most parsimonious explanation operating on these innate immunity genes. This genetic surveillance system recognizes motifs in pathogens that are perhaps conserved across a broad range of pathogens. Hence, functional constraints are imposed on mutations that diminish the ablility of these proteins to detect pathogens.

Comparative analyses of genetic risk prediction methods reveal extreme diversity of genetic predisposition to nonalcoholic fatty liver disease (NAFLD) among ethnic populations of India.

Nonalcoholic fatty liver disease (NAFLD) is a distinct pathologic condition characterized by a disease spectrum ranging from simple steatosis to steato-hepatitis, cirrhosis and hepatocellular carcinoma. Prevalence of NAFLD varies in different ethnic groups, ranging from 12% in Chinese to 45% in Hispanics. Among Indian populations, the diversity in prevalence is high, ranging from 9% in rural populations to 32% in urban populations, with geographic differences as well. Here, we wished to find out if this difference is reflected in their genetic makeup. To date, several candidate genes and a few genomewide association studies (GWAS) have been carried out, and many associations between single nucleotide polymorphisms (SNPs) and NAFLD have been observed. In this study, the risk allele frequencies (RAFs) of NAFLD-associated SNPs in 20 Indian ethnic populations (376 individuals) were analysed. We used two different measures for calculating genetic risk scores and compared their performance. The correlation of additive risk scores of NAFLD for three Hapmap populations with their weighted mean prevalence was found to be high (R(2) = 0.93). Later we used this method to compare NAFLD risk among ethnic Indian populations. Based on our observation, the Indian caste populations have high risk scores compared to Caucasians, who are often used as surrogate and similar to Indian caste population in disease gene association studies, and is significantly higher than the Indian tribal populations.

Variant allele frequency enrichment analysis in vitro reveals sonic hedgehog pathway to impede sustained temozolomide response in GBM.

Neoplastic cells of Glioblastoma multiforme (GBM) may or may not show sustained response to temozolomide (TMZ) chemotherapy. We hypothesize that TMZ chemotherapy response in GBM is predetermined in its neoplastic clones via a specific set of mutations that alter relevant pathways. We describe exome-wide enrichment of variant allele frequencies (VAFs) in neurospheres displaying contrasting phenotypes of sustained versus reversible TMZ-responses in vitro. Enrichment of VAFs was found on genes ST5, RP6KA1 and PRKDC in cells showing sustained TMZ-effect whereas on genes FREM2, AASDH and STK36, in cells showing reversible TMZ-effect. Ingenuity pathway analysis (IPA) revealed that these genes alter cell-cycle, G2/M-checkpoint-regulation and NHEJ pathways in sustained TMZ-effect cells whereas the lysine-II&V/phenylalanine degradation and sonic hedgehog (Hh) pathways in reversible TMZ-effect cells. Next, we validated the likely involvement of the Hh-pathway in TMZ-response on additional GBM neurospheres as well as on GBM patients, by extracting RNA-sequencing-based gene expression data from the TCGA-GBM database. Finally, we demonstrated TMZ-sensitization of a TMZ non-responder neurosphere in vitro by treating them with the FDA-approved pharmacological Hh-pathway inhibitor vismodegib. Altogether, our results indicate that the Hh-pathway impedes sustained TMZ-response in GBM and could be a potential therapeutic target to enhance TMZ-response in this malignancy.

Genomic correlates of variability in immune response to an oral cholera vaccine.

Cholera is endemic to many countries. Recent major outbreaks of cholera have prompted World Health Organization to recommend oral cholera vaccination as a public-health strategy. Variation in percentage of seroconversion upon cholera vaccination has been recorded across populations. Vaccine-induced responses are influenced by host genetic differences. We have investigated association between single-nucleotide polymorphic (SNP) loci in and around 296 immunologically relevant genes and total anti-lipopolysaccharide (LPS) antibody response to a killed whole-cell vaccine, comprising LPS from multiple strains of Vibrio cholerae. Titers derived from standard vibriocidal assays were also analyzed to gain further insights on validated SNP associations. Vaccination was administered to 1000 individuals drawn from India. Data on two independent random subsets, each comprising ∼500 vaccinees, were used for discovery of genomic associations and validation, respectively. Significant associations of four SNPs and haplotypes in three genes (MARCO, TNFAIP3 and CXCL12) with AR were discovered and validated, of which two in TNFAIP3 and CXCL12 were also significantly associated with immunity (fourfold increase in vibriocidal titers). CXCL12 is a neutrophil and lymphocyte chemoattractant that is upregulated in response to V. cholerae infection. LPS in the vaccine possibly provides signals that mimic those of the live bacterium. TNFAIP3 promotes intestinal epithelial barrier integrity and provides tight junction protein regulation; possible requirements for adequate response to the vaccine. LPS is a potent activator of innate immune responses and a ligand of MARCO. Variants in this gene have been found to be associated with LPS response, but not with high vibriocidal titer level.

A novel statistical algorithm for enhancing the utility of HapMap data to design genomic association studies in non-HapMap populations.

The HapMap database should be effectively used in designing disease association studies in non-HapMap populations. The efficiency of portability of tagSNPs from HapMap to non-HapMap populations is widely variable. A new algorithm is proposed for selecting SNPs from HapMap for use in non-HapMap populations by simultaneously considering and combining data on allele frequencies and linkage-disequilibrium values in the four HapMap populations. Empirical comparison and validation of the algorithm are provided by using Tagger, available HapMap data and data from an Indian population. The proposed method is shown to be efficient and effective. A software implementing this algorithm is freely available.

Genetic determinants of immune-response to a polysaccharide vaccine for typhoid.

UNLABELLED: Differences in immunological response among vaccine recipients are determined both by their genetic differences and environmental factors. Knowledge of genetic determinants of immunological response to a vaccine can be used to design a vaccine that circumvents immunogenetic restrictions. The currently available vaccine for typhoid is a pure polysaccharide vaccine, immune response to which is T-cell independent. Little is known about whether genetic variation among vaccinees associates with variation in their antibody response to a polysaccharide vaccine. We conducted a study on 1,000 individuals resident in an area at high-risk for typhoid; vaccinated them with the typhoid vaccine, measured their antibody response to the vaccine, assayed >2,000 curated SNPs chosen from 283 genes that are known to participate in immune-response; and analyzed these data using a strategy to (a) minimize the statistical problems associated with testing of multiple hypotheses, and (b) internally cross-validate inferences, using a half-sample design, with little loss of statistical power. The first stage analysis, using the first half-sample, identified 54 SNPs in 43 genes to be significantly associated with immune response. In the second-stage, these inferences were cross-validated using the second half-sample. First-stage results of only 8 SNPs (out of 54) in 7 genes (out of 43) were cross-validated. We tested additional SNPs in these 7 genes, and found 8 more SNPs to be significantly associated. Haplotypes constructed with these SNPs in these 7 genes also showed significant association. These 7 genes are DEFB1, TLR1, IL1RL1, CTLA4, MAPK8, CD86 and IL17D. The overall picture that has emerged from this study is that (a) immune response to polysaccharide antigens is qualitatively different from that to protein antigens, and (b) polymorphisms in genes involved in polysaccharide recognition, signal transduction, inhibition of T-cell proliferation, pro-inflammatory signaling and eventual production of antimicrobial peptides are associated with antibody response to the polysaccharide vaccine for typhoid. ELECTRONIC SUPPLEMENTARY MATERIAL: The online version of this article (doi:10.1007/s11568-010-9134-1) contains supplementary material, which is available to authorized users.